Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: over the last more than three decades ferroelectric liquid crystals ( flcs ) have attracted considerable attention as promising chiral liquid crystal materials for applications in fast switching display devices ( a detailed description of flcs can be found , e.g. , in monographs  @xcite ) . equilibrium orientational structures in flcs are represented by helical twisting patterns where flc molecules align on average along a local unit director @xmath2 where @xmath3 is the smectic tilt angle ; @xmath4 is the twisting axis normal to the smectic layers and @xmath5 is the @xmath6-director . the flc director   lies on the smectic cone depicted in fig . [ subfig : director ] with the _ smectic tilt angle _ @xmath3 and rotates in a helical fashion about a uniform twisting axis @xmath4 forming the flc helix with the _ helix pitch _ , @xmath7 . this rotation is described by the azimuthal angle around the cone @xmath8 that specifies orientation of the @xmath6-director in the plane perpendicular to @xmath4 and depends on the dimensionless coordinate along the twisting axis @xmath9 where @xmath10 is the helix twist wave number . the important case of a uniform lying flc helix in the slab geometry with the smectic layers normal to the substrates and @xmath11 where @xmath12 is the electric field applied across the cell , is illustrated in fig . [ fig : geom ] . this is the geometry of surface stabilized flcs ( ssflcs ) pioneered by clark and lagerwall in  @xcite . they studied electro - optic response of flc cells confined between two parallel plates subject to homogeneous boundary conditions and made thin enough to suppress the bulk flc helix . it was found that such cells exhibit high - speed , bistable electro - optical switching between orientational states stabilized by surface interactions . the response of flcs to an applied electric field @xmath12 is characterized by fast switching times due to linear coupling between the field and the _ spontaneous ferroelectric polarization _ @xmath13 where @xmath14 is the _ polarization unit vector_. there is also a threshold voltage necessary for switching to occur and the process of bistable switching is typically accompanied by a hysteresis . figure  [ subfig : cell ] also describes the geometry of deformed helix flcs ( dhflcs ) as it was introduced in  @xcite . this case will be of our primary concern .    in dhflc cells , the flc helix is characterized by a short submicron _ helix pitch _ , @xmath15  @xmath16 m , and a relatively large _ tilt angle _ , @xmath17  deg . by contrast to ssflc cells , where the surface induced unwinding of the bulk helix requires the helix pitch of a flc mixture to be greater than the cell thickness , a dhflc helix pitch is 5 - 10 times smaller than the thickness . this allows the helix to be retained within the cell boundaries . electro - optical response of dhflc cells exhibits a number of peculiarities that make them useful for lc devices such as high speed spatial light modulators  @xcite , colour - sequential liquid crystal display cells  @xcite and optic fiber sensors  @xcite . the effects caused by electric - field - induced distortions of the helical structure underline the mode of operation of such cells . in a typical experimental setup , these effects are probed by performing measurements of the transmittance of normally incident linearly polarized light through a cell placed between crossed polarizers . a more general case of oblique incidence has not received a fair amount of attention . theoretically , a powerful tool to deal with this case is the transfer matrix method which has been widely used in studies of both quantum mechanical and optical wave fields  @xcite . in this work we apply the method for systematic treatment of the technologically important case of dhflcs with subwavelength pitch also known as the _ short - pitch dhflcs_.    recently , the transfer matrix approach to polarization gratings was employed to define the effective dielectric tensor of short - pitch dhflcs  @xcite that gives the principal values and orientation of the optical axes as a function of the applied electric field . biaxial anisotropy and rotation of the in - plane optical axes produced by the electric field can be interpreted as the _ orientational kerr effect _ @xcite . it can be expected that the electric field dependence of the effective dielectric tensor will also manifest itself as electric - field - induced transformations of the polarization - resolved angular ( conoscopic ) patterns in the observation plane after the dhflc cells illuminated by convergent light beam . these patterns are represented by the two - dimensional ( 2d ) fields of polarization ellipses describing the polarization structure behind the conoscopic images  @xcite .    as it was originally recognized by nye  @xcite , the key elements characterizing geometry of such stokes parameter fields are the _ polarization singularities _ that play the fundamentally important role of structurally stable topological defects ( a recent review can be found in ref . in particular , the polarization singularities such as the _ c points _ ( the points where the light wave is circular polarized ) and the _ l lines _ ( the curves along which the polarization is linear ) frequently emerge as the characteristic feature of certain polarization state distributions . for nematic and cholescteric ( chiral nematic ) liquid crystals , the singularity structure of the polarization - resolved angular patterns is generally found to be sensitive to both the director configuration and the polarization characteristics of incident light  @xcite .    in this study , we consider the polarization - resolved angular patterns of dhflc cells as the stokes parameter fields giving detailed information on the incidence angles dependence of the polarization state of light transmitted through the cells . in particular , we explore how the polarization singularities transform under the action of the electric field . our analysis will utilize the transfer matrix approach in combination with the results for the effective dielectric tensor of biaxial flcs evaluated using an improved technique of averaging over distorted helical structures . we also emphasize the role of phase singularities of a different kind and discuss the electro - optic behavior of dhflcs near the exceptional point where the condition of zero - field isotropy is fulfilled . the layout of the paper is as follows . in sec . [ sec : theory ] we introduce our notations and describe the transfer matrix formalism rendered into the @xmath18 matrix form suitable for our purposes . this formalism is employed to deduce a number of the unitarity and symmetry relations with emphasis on the planar anisotropic structures that represent dhflc cells and posses two optical axes lying in the plane of substrates . in sec . [ sec : optics - dhf ] we evaluate the effective dielectric tensor of dhflc cells , discuss the orientational kerr effect and show that electro - optic response of dhflc cells is enhanced near the exceptional point determined by the condition of zero - field isotropy . geometry of the polarization - resolved angular patterns emerging after dhflc cells is considered in sec . [ sec : conoscopy ] . after providing necessary details on our computational approach and the polarization singularities , we present the numerical results describing how the singularity structure of polarization ellipse fields transforms under the action of the electric field . finally , in sec . [ sec : discussion ] we draw the results together and make some concluding remarks . details on some technical results are relegated to appendixes  [ sec : op - evol][sec : flc - spiral ] . in order to describe both the electro - optical properties and the polarization - resolved angular patterns of deformed helix ferroelectric liquid crystal layers with subwavelength pitch we adapt a general theoretical approach which can be regarded as a modified version of the well - known transfer matrix method  @xcite and was previously applied to study the polarization - resolved conoscopic patterns of nematic liquid crystal cells  @xcite . this approach has also been extended to the case of polarization gratings and used to deduce the general expression for the effective dielectric tensor of dhflc cells  @xcite .    in this section , we present the transfer matrix approach as the starting point of our theoretical considerations , with emphasis on its general structure and the symmetry relations . the analytical results for uniformly anisotropic planar structures representing homogenized dhflc cells are given in appendix  [ sec : planar ] . we deal with a harmonic electromagnetic field characterized by the free - space wave number @xmath19 , where @xmath20 is the frequency ( time - dependent factor is @xmath21 ) , and consider the slab geometry shown in fig . [ fig : geom_4rays ] . in this geometry , an optically anisotropic layer of thickness @xmath22 is sandwiched between the bounding surfaces ( substrates ) : @xmath23 and @xmath24 ( the @xmath25 axis is normal to the substrates ) and is characterized by the dielectric tensor @xmath26 and the magnetic permittivity @xmath16    further , we restrict ourselves to the case of stratified media and assume that the electromagnetic fields can be taken in the following factorized form @xmath27 where the vector @xmath28 represents the lateral component of the wave vector . then we write down the representation for the electric and magnetic fields , @xmath12 and @xmath29 , @xmath30 where the components directed along the normal to the bounding surface ( the @xmath25 axis ) are separated from the tangential ( lateral ) ones . in this representation , the vectors @xmath31 and @xmath32 are parallel to the substrates and give the lateral components of the electromagnetic field .    substituting the relations   into the maxwell equations and eliminating the @xmath25 components of the electric and magnetic fields gives equations for the tangential components of the electromagnetic field that can be written in the following @xmath18 matrix form  @xcite : @xmath33 where @xmath34 is the _ differential propagation matrix _ and its @xmath35 block matrices @xmath36 are given by    [ eq : mij ] @xmath37    general solution of the system   @xmath38 can be conveniently expressed in terms of the _ evolution operator _ which is also known as the _ propagator _ and is defined as the matrix solution of the initial value problem    [ eq : evol_problem ] @xmath39    where @xmath40 is the @xmath41 identity matrix . basic properties of the evolution operator are reviewed in appendix  [ sec : op - evol ] .      in the ambient medium with @xmath42 and @xmath43 , the general solution   can be expressed in terms of plane waves propagating along the wave vectors with the tangential component  . for such waves , the result is given by @xmath44 where @xmath45 is the eigenvector matrix for the ambient medium given by @xmath46 @xmath47 are the pauli matrices @xmath48    from eq .  , the vector amplitudes @xmath49 and @xmath50 correspond to the forward and backward eigenwaves with @xmath51 and @xmath52 , respectively . figure  [ fig : geom_4rays ] shows that , in the half space @xmath53 before the entrance face of the layer @xmath23 , these eigenwaves describe the _ incoming and outgoing waves _ @xmath54 whereas , in the half space @xmath55 after the exit face of the layer , these waves are given by @xmath56 in this geometry , there are two plane waves , @xmath57 and @xmath58 , incident on the bounding surfaces of the anisotropic layer , @xmath23 and @xmath24 , respectively . then the standard linear input - output relations @xmath59 linking the vector amplitudes of transmitted and reflected waves , @xmath60 and @xmath61 with the amplitude of the incident wave , @xmath62 through the transmission and reflection matrices , @xmath63 and @xmath64 , assume the following generalized form : @xmath65 where @xmath66 is the matrix     the so - called _ scattering matrix _   that relates the outgoing and incoming waves ; @xmath67 ( @xmath68 ) is the transmission ( reflection ) matrix for the case when the incident wave is incoming from the half space @xmath53 bounded by the entrance face , whereas the mirror symmetric case where the incident wave is impinging onto the exit face of the sample is described by the transmission ( reflection ) matrix @xmath69 ( @xmath70 ) . so , we have    @xmath71    it is our task now to relate these matrices and the evolution operator given by eq .  . to this end , we use the boundary conditions requiring the tangential components of the electric and magnetic fields to be continuous at the boundary surfaces : @xmath72 and @xmath73 , and apply the relation   to the anisotropic layer of the thickness @xmath22 to yield the following result @xmath74      on substituting eqs .   into eq . we have @xmath75 where the matrix @xmath76 linking the electric field vector amplitudes of the waves in the half spaces @xmath77 and @xmath78 bounded by the faces of the layer will be referred to as the _ transfer ( linking ) matrix_. the expression for the transfer matrix is as follows @xmath79 where @xmath80 is the rotated operator of evolution . this operator is the solution of the initial value problem   with @xmath81 replaced with @xmath82 .    from eqs .   and , the block structure of the transfer matrix can be expressed in terms of the transmission and reflection matrices as follows @xmath83 similarly , for inverse of the transfer matrix , @xmath84 we have @xmath85      in appendix  [ sec : op - evol ] , it is shown that , for non - absorbing media with symmetric dielectric tensor , @xmath86 , the operator of evolution satisfies the unitarity relation  . by using eq .   in combination with the algebraic identity    [ eq : vm_identity ] @xmath87    where @xmath88 , for the eigenvector matrix given in eq . , we can deduce the unitarity relation for the transfer matrix   @xmath89    the unitarity relation   for non - absorbing layers can now be used to derive the energy conservation laws    [ eq : energy_laws ] @xmath90}^{\dagger } } { { { \mathbf{t}}}^{t}}_{\pm}+ { { [ { { { \mathbf{r}}}^{t}}_{\mp}]}^{\dagger } } { { { \mathbf{r}}}^{t}}_{\mp } = { \mathbf{i}}_2 , \ ] ]    where a dagger and the superscript @xmath91 will denote hermitian conjugation and matrix transposition , respectively , along with the relations for the block matrices    [ eq : w_ij_unit ] @xmath92}^{\dagger } } , \\ &   \label{eq : w_12_unit }    { \mathbf{w}}_{12}=-{\mathbf{t}}_{+}^{-1 } { \mathbf{r}}_{-}= { { [ { \mathbf{t}}_{-}^{-1 } { \mathbf{r}}_{+}]}^{\dagger } } , \\ &   \label{eq : w_21_unit }    { \mathbf{w}}_{21}={\mathbf{r}}_{+ } { \mathbf{t}}_{+}^{-1}= -{{[{\mathbf{r}}_{- } { \mathbf{t}}_{-}^{-1}]}^{\dagger}}.\end{aligned}\ ] ]    note that eqs .   and   can be conveniently rewritten in the following form    [ eq : rel - w_12_21 ] @xmath93}^{\dagger } } { { { \mathbf{r}}}^{\dagger}}_{+},\end{aligned}\ ] ]    so that multiplying these identities and using the energy conservation law   gives the relations  .    in the translation invariant case of uniform anisotropy , the matrix @xmath34 is independent of @xmath94 and the operator of evolution is given by @xmath95 then , the unitarity condition  @xcite @xmath96 can be combined with eq .   to yield the additional symmetry relations for @xmath97 @xmath98 where an asterisk will indicate complex conjugation , that give the following algebraic identities for the transmission and reflection matrices : @xmath99 it can be readily seen that the relation for the transposed matrices   can be derived by substituting eq .   into the conservation law  .    for the important special case of uniformly anisotropic planar structures with @xmath100 , the algebraic structure of the transfer matrix is described in appendix  [ sec : planar ] . equation   shows that the symmetry relations   remain valid even if the dielectric constants are complex - valued and the medium is absorbing . since identities   are derived from eqs .   and   without recourse to the unitarity relations , they also hold for lossy materials . similar remark applies to the expression for inverse of the transfer matrix  . from eq . , it follows that the relation between the transmission ( reflection ) matrix , @xmath67 ( @xmath68 ) , and its mirror symmetric counterpart @xmath69 ( @xmath70 ) can be further simplified and is given by @xmath101 from eqs .   and  , we have the relation for the transposed reflection matrices @xmath102 whereas the transmission matrices are symmetric . we now pass on to the electro - optical properties of dhflc cells and extend the results of ref . @xcite to the case of biaxial ferroelectric liquid crystals with subwavelength pitch . in addition , the theoretical treatment will be significantly improved by using an alternative fully consistent procedure to perform averaging over distorted flc helix that goes around the limitations of the first - order approximation . we consider a flc film of thickness @xmath22 with the @xmath25 axis which , as is indicated in fig . [ fig : geom ] , is normal to the bounding surfaces : @xmath23 and @xmath24 , and introduce the effective dielectric tensor , @xmath103 , describing a homogenized dhflc helical structure . for a biaxial flc , the components of the dielectric tensor , @xmath104 , are given by @xmath105 where @xmath106 , @xmath107 is the kronecker symbol ; @xmath108 ( @xmath109 ) is the @xmath110th component of the flc director ( unit polarization vector ) given by eq .   ( eq .  ) ; @xmath111 are the _ anisotropy parameters _ and @xmath112 ( @xmath113 ) is the _ anisotropy ( biaxiality ) ratio_. note that , in the case of uniaxial anisotropy with @xmath114 , the principal values of the dielectric tensor are : @xmath115 and @xmath116 , where @xmath117 ( @xmath118 ) is the ordinary ( extraordinary ) refractive index and the magnetic tensor of flc is assumed to be isotropic with the magnetic permittivity @xmath16 . as in sec . [ sec : theory ] ( see fig .  [ fig : geom_4rays ] ) , the medium surrounding the layer is optically isotropic and is characterized by the dielectric constant @xmath119 , the magnetic permittivity @xmath120 and the refractive index @xmath121 .    at @xmath122 , the ideal flc helix @xmath123 where @xmath124 is the free twist wave number and @xmath125 is the equilibrium helical pitch , is defined through the azimuthal angle around the smectic cone @xmath8 ( see fig . [ fig : geom ] and eq .  ) and represents the undistorted structure . for sufficiently small electric fields , the standard perturbative technique applied to the euler - lagrange equation gives the first - order expression  @xcite for the azimuthal angle of a weakly distorted helical structure @xmath126 where @xmath127 is the electric field parameter linearly proportional to the ratio of the applied and critical electric fields : @xmath128 , and @xmath129 .    according to ref . @xcite , normally incident light feels effective in - plane anisotropy described by the averaged tensor , @xmath130 : @xmath131 where @xmath132 , and the effective dielectric tensor @xmath133 can be expressed in terms of the averages @xmath134^{-1}}\rangle } } , \label{eq : eta } \\ &      \beta_{z\alpha}= { \ensuremath{\langle{\epsilon_{z\alpha}/\epsilon_{zz}}\rangle } } = \left\langle \frac{u_1 d_z d_\alpha+u_2 p_z p_\alpha}{1+u_1 d_z^2+u_2 p_z^2 } \right\rangle , \label{eq : beta }    \end{aligned}\ ] ] as follows @xmath135    general formulas  - give the zero - order approximation for homogeneous models describing the optical properties of short pitch dhflcs  @xcite . assuming that the pitch - to - wavelength ratio @xmath136 is sufficiently small , these formulas can now be used to derive the effective dielectric tensor of homogenized short - pitch dhflc cell for both vertically and planar aligned flc helix . the results for vertically aligned dhflc cells were recently published in ref . @xcite and we concentrate on the geometry of planar aligned dhflc helix shown in fig . [ fig : geom ] . for this geometry , the parameters needed to compute the averages @xmath130 ( see eq .  ) , @xmath137 ( see eq .  ) and @xmath138 ( see eq .  ) are given by @xmath139    formulas   can now be inserted into eqs .   to yield the explicit expressions for the elements of the dielectric tensor  :    [ eq : epsilon - ij - planar ] @xmath140      the simplest averaging procedure previously used in refs . @xcite involves substituting the formula for a weakly distorted flc helix   into eqs .   and performing integrals over @xmath141 . this procedure thus heavily relies on the first - order approximation where the director distortions are described by the term linearly proportional to the electric field ( the second term on the right hand side of eq .  ) . quantitatively , the difficulty with this approach is that the linear approximation may not be suffice for accurate computing of the second - order contributions to the diagonal elements of the dielectric tensor  . in this approximation , the second - order corrections describing the helix distortions that involve the change of the helix pitch have been neglected .    in order to circumvent the problem , in this paper , we apply an alternative approach that allows to go beyond the first - order approximation without recourse to explicit formulas for the azimuthal angle . this method is detailed in appendix  [ sec : flc - spiral ] . the analytical results   substituted into eqs .   give the effective dielectric tensor in the following form : @xmath142 the _ zero - field dielectric constants _ , @xmath143 and @xmath144 , that enter the tensor   are defined in eqs .   and  , respectively , and can be conveniently rewritten as follows    [ eq : epsilon_ph ] @xmath145    a similar result for the _ coupling coefficients _ @xmath146 , @xmath147 and @xmath148 ( see eq .  ) reads    [ eq : coupling - coeffs ] @xmath149    note that , following ref . @xcite , we have used the relation   to introduce the electric field parameter @xmath150 where @xmath151 is the dielectric susceptibility of the goldstone mode  @xcite .    the above dielectric tensor is characterized by the three generally different principal values ( eigenvalues ) and the corresponding optical axes ( eigenvectors ) as follows @xmath152 ^ 2+[\gamma_{xy } \alpha_{e}]^2}\end{aligned}\ ] ] where @xmath153 the in - plane optical axes are given by @xmath154.\end{aligned}\ ] ] from eq . , it is clear that , similar to the case of uniaxial flcs studied in ref . @xcite , the zero - field dielectric tensor is uniaxially anisotropic with the optical axis directed along the twisting axis @xmath155 . the applied electric field changes the principal values ( see eqs .   and  ) so that the electric - field - induced anisotropy is generally biaxial . in addition , the in - plane principal optical axes are rotated about the vector of electric field , @xmath156 , by the angle @xmath157 given in eq .  . in the low electric field region , the electrically induced part of the principal values is typically dominated by the kerr - like nonlinear terms proportional to @xmath158 , whereas the electric field dependence of the angle @xmath157 is approximately linear : @xmath159 . this effect is caused by the electrically induced distortions of the helical structure and bears some resemblance to the electro - optic kerr effect . following refs . @xcite , it will be referred to as the _ orientational kerr effect_.    it should be emphasized that this effect differs from the well - known kerr effect which is a quadratic electro - optic effect related to the electrically induced birefringence in optically isotropic ( and transparent ) materials and which is mainly caused by the electric - field - induced orientation of polar molecules  @xcite . by contrast , in our case , similar to polymer stabilized blue phase liquid crystals  @xcite , we deal with the effective dielectric tensor of a nanostructured chiral smectic liquid crystal . this tensor   is defined through averaging over the flc orientational structure . typically , in experiments dealing with the electro - optic response of dhflc cells , the transmittance of normally incident light passing through crossed polarizers is measured as a function of the applied electric field . for normal incidence , the transmission and reflection matrices can be easily obtained from the results given in appendix  [ sec : planar ] by substituting eq .   into eqs .  -  . when the incident wave is linearly polarized along the @xmath160 axis ( the helix axis ) , the transmittance coefficient @xmath161 where @xmath162 is the thickness parameter , describes the intensity of the light passing through crossed polarizers . note that , under certain conditions such as @xmath163 , @xmath164 and the transmittance   can be approximated by simpler formula @xmath165 where @xmath166 is the difference in optical path of the ordinary and extraordinary waves known as the _ phase retardation_.    in ref . @xcite , the relation   was used to fit the experimental data using the theory based on the linear approximation for the helix distortions ( see eq .  ) . these results are reproduced in figure  [ fig : expt_th ] along with the theoretical curve computed using the modified averaging technique . from fig . [ fig : expt_th ] , it is seen that , in the range of relatively high voltages , the averaging method described in appendix  [ sec : flc - spiral ] improves agreement between the theory and the experiment , whereas , at small voltages , the difference between the fitting curves is negligibly small . given the anisotropy and biaxiality ratios , @xmath167 and @xmath168 , the zero - field dielectric constants   and the coupling coefficients   are determined by the smectic tilt angle , @xmath3 . figure  [ fig : gamma_theta ] shows how the coupling coefficients depend on @xmath3 for both uniaxially and biaxially anisotropic flcs . as it can be seen in fig . [ subfig : gamma_uni ] , in the case of conventional flcs with @xmath169 , all the coefficients are positive and the difference of the coupling constants @xmath170 that define the electrically induced part of @xmath171 ( see eq .  ) is negative at @xmath172 .    from eqs .   and , it follows that , at @xmath169 , the principal values of dielectric constants @xmath173 and @xmath174 are decreasing functions of the electric field parameter @xmath175 so that anisotropy of the effective dielectric tensor   is weakly biaxial . in addition , for non - negative @xmath175 and @xmath176 , the azimuthal angle of in - plane optical axis , @xmath177 , given in eq . increases with @xmath175 from zero to @xmath178 . figure  [ subfig : gamma_bia ] demonstrates that this is no longer the case for biaxial flcs . it is seen that , at @xmath179 , the coupling coefficient @xmath147 and the difference @xmath170 both change in sign when the tilt angle @xmath3 is sufficiently small . at such angles , the dielectric constant @xmath173 increases with @xmath175 and the electric field induced anisotropy of dhflc cell becomes strongly biaxial . when @xmath170 and @xmath175 are positive , electric field dependence of the azimuthal angle @xmath177 is non - monotonic and the angle decays to zero in the range of high voltages where @xmath180 .      at @xmath122 , the zero - field anisotropy is uniaxial and is described by the dielectric constants , @xmath143 and @xmath181 , given in eq .  . [ fig : eph_theta ] , these constants are plotted against the tilt angle . it is shown that , at small tilt angles , the anisotropy @xmath182 is positive . it decreases with @xmath3 and the zero - field state becomes isotropic when , at certain critical angle @xmath183 , the condition of zero - field isotropy @xmath184 is fulfilled and @xmath185 . so , the angle @xmath186 can be referred to as the _ isotropization angle_. in what follows we discuss peculiarities of the electro - optic response in the vicinity of the _ isotropization point _ where @xmath171 is proportional to @xmath158 ( see eq .  ) and the kerr - like regime breaks down . mathematically , the isotropization point represents a square root branch - point singularity of the eigenvalues   of the dielectric tensor which is known as the _ exceptional point _  @xcite . in the electric field dependence of the in - plane dielectric constants , @xmath187 and @xmath174 , this singularity reveals itself as a cusp where the derivatives of @xmath188 with respect to @xmath175 are discontinuous . more precisely , we have @xmath189 as is illustrated in fig . [ subfig : epsilon_c_uni ] , the cusp is related to the effect of reconnection of different branches representing solutions of an algebraic equation .    since the azimuthal angle @xmath177 is undetermined at @xmath183 @xmath190 the isotropization point also represents a phase singularity . the electric field dependence of @xmath177 is thus discontinuous and the relation @xmath191 describes its jumplike behaviour at @xmath122 . this behaviour is demonstrated in fig . [ fig : psi_alp ] .    we can now use eq .   and write down the condition of zero - field isotropy   in the following explicit form : @xmath192 the case of a uniaxially anisotropic flc with @xmath169 can be treated analytically . in this case , it is not difficult to check that @xmath193 gives the special solution of eq . that does not depend on the tilt angle and corresponds to an isotropic material with @xmath194 . another solution is given by the relation @xmath195 linking the isotropization angle @xmath186 and the anisotropy parameter , @xmath196 . in fig . [ fig : theta_r1 ] , this solution is represented by the solid line curve . the isotropization angle is shown to be a slowly decreasing function of the anisotropy ratio @xmath167 . from eq . , it starts from the maximal value of @xmath186 is @xmath197 and decays approaching @xmath198 . when the biaxiality ratio @xmath168 differs from unity , the solution of the isotropy condition   can only be written in the parametrized form as follows @xmath199 where @xmath200 the @xmath186 versus @xmath167 curves computed from the representation   are shown in fig . [ fig : theta_r1 ] . it can be seen that , by contrast to the case of uniaxial anisotropy [ @xmath169 ] , for biaxial flcs with @xmath202 , each curve has two branches separated by a gap . the isotropization angle vanishes , @xmath203 , at one of the endpoints of the gap , @xmath204 , whereas the angle @xmath186 equals @xmath178 at the other endpoint . thus , dependence of the isotropization angle on the anisotropy ratio @xmath167 being smooth and continuous for uniaxially anisotropic flcs is found to be splitted into two branches when the flc anisotropy is biaxial . from fig . [ fig : theta_r1 ] , one of the branches with @xmath205 is associated with the endpoint at @xmath204 and lies below of the solid line curve representing conventional flcs . for this branch , the angle @xmath186 decreases with the biaxiality ratio reaching zero at @xmath206 . we conclude this section with the remark on the electro - optic response of dhflc cells in the vicinity of the isotropization point . one of the important factors governing the electric field dependence of the transmittance   is the phase retardation @xmath207 ( see eq .  ) proportional the effective birefringence @xmath208 . the electrically dependent part of this birefringence is plotted as a function of the electric field parameter in fig . [ fig : dn_alp ] . it shown that , at @xmath183 , the kerr - like regime breaks down and the birefringence is dominated by the terms linearly dependent on the electric field . such a pockels - like behaviour manifests itself in the perfectly harmonic dependence of the transmittance , @xmath209 , on the electric field parameter depicted in fig . [ subfig : transm_c_uni ] .    figure  [ subfig : transm_c-2_uni ] illustrates the effect of small deviations from the isotropization angle . though the curve presented in fig . [ subfig : transm_c-2_uni ] and the ones for flc-576a ( see fig . [ fig : expt_th ] ) are quite similar in shape , it is clear that sensitivity to the electric field and the magnitude of transmission peaks are both considerably enhanced near the isotropization point . such behaviour comes as no surprise and derives from the above discussed fact that this point plays the role of a singularity ( an exceptional point ) ( see eqs .   ) . we can now combine the general relations deduced in sec .  [ sec : theory ] ( and in appendix  [ sec : planar ] ) using the transfer matrix method with the results of sec . [ sec : optics - dhf ] to study the _ polarization - resolved angular ( conoscopic ) patterns _ describing the polarization structure behind the conoscopic images of short - pitch dhflc cells that are characterized by the effective dielectric tensor  . this polarization structure is represented by a two - dimensional ( 2d ) distribution of polarization ellipses and results from the interference of eigenmodes excited in the dhflc cells by the plane waves with varying direction of incidence . geometrically , the important elements of the 2d stokes parameter fields are the polarization singularities such as @xmath1 points ( the points where the light wave is circular polarized ) and @xmath0 lines ( the curves along which the polarization is linear ) . in this section the focus of our attention will be on the singularity structure of the polarization - resolved angular patterns emerging after the dhflc cells . our starting point is the computational method used to evaluate the patterns as the polarization ellipse fields . we shall use the electric field vector amplitudes of incident , reflected and transmitted waves conveniently rewritten in the circular basis @xmath210 and the incidence angles , @xmath211 and @xmath212 , related to the lateral component of the wave vector   as follows @xmath213 where @xmath211 ( @xmath212 ) is the _ polar ( azimuthal ) angle of incidence_. dependence of the polarization properties of the waves transmitted through the dhflc cell on the incidence angles , @xmath211 and @xmath212 , will be of our primary concern . the transmission matrix describing conoscopic patterns on the transverse plane of projection is given by  @xcite @xmath214 where @xmath215 and @xmath216 are the polar coordinates in the observation plane ( @xmath217 and @xmath218 are the cartesian coordinates ) and @xmath219 is the aperture dependent scale factor .    the transmission matrix of dhflc cells , @xmath220 , can be computed from general formulas given in appendix  [ sec : planar ] ( see eq .  ) . for this matrix , the parameters @xmath221 that enter the expression for the dielectric tensor of planar structures   should be replaced with the characteristics @xmath222 of the effective dielectric tensor   given in eqs .   . in what follows the incident light is assumed to be linearly polarized @xmath223 where @xmath224 is the polarization azimuth of the incident wave , and the state of polarization of the transmitted wave @xmath225 is defined by the polarization ellipse characteristics . the orientation of the polarization ellipse is specified by the azimuthal angle of polarization ( _ polarization azimuth _ ) @xmath226}^{\ast } } e_{-}^{({\mathrm{trm}})},\end{aligned}\ ] ] where @xmath227 is the @xmath110th component of the stokes vector , and its eccentricity is described by the signed ellipticity parameter @xmath228 that will be referred to as the _ ellipticity_. the ellipse is considered to be right handed ( left handed ) if its helicity is positive ( negative ) , so that @xmath229 ( @xmath230 ) . +    from eq .  , the incidence angles and the points in the observation plane are in one - to - one correspondence . so , computing the polarization azimuth , @xmath231 , the ellipticity , @xmath232 , at each point of the projection plane yields the 2d field of polarization ellipses which is called the polarization - resolved angular ( conoscopic ) pattern . the point where @xmath233 and thus the transmitted wave is circularly polarized with @xmath234 will be referred to as the _ c@xmath235 point_.    this is an example of the _ polarization singularity _ that can be viewed as the _ phase singularities _ of the complex scalar field @xmath236 where the phase @xmath237 ( see eq .  ) become indeterminate . such singularities are characterized by the _ winding number _ which is the signed number of rotations of the two - component field @xmath238 around the circuit surrounding the singularity  @xcite . the winding number also known as the _ signed strength of the dislocation _ is generically @xmath239 . since the polarization azimuth   is defined modulo @xmath240 and @xmath241 , the dislocation strength is twice the index of the corresponding @xmath242 point , @xmath243 . for generic @xmath1 points , @xmath244 and the topological index can be computed as the closed - loop contour integral of the phase @xmath245 modulo @xmath246 @xmath247 where @xmath0 is the closed path around the singularity . in addition to the handedness and the index , the @xmath1 points are classified according to the number of streamlines , which are polarization lines whose tangent gives the polarization azimuth , terminating on the singularity . this is the so - called _ line classification _ that was initially studied in the context of umbilic points  @xcite . mathematically , the straight streamlines that terminate on the singularity are of particular importance as they play the role of separatrices , separating regions of streamlines with differently signed curvature . as is illustrated in fig . [ fig : singul - types ] , for generic @xmath1 points , the number of the straight lines , @xmath248 , may either be 1 or 3 . this number is @xmath249 provided the index equals @xmath250 , @xmath251 , and such @xmath1 points are called _ stars_. at @xmath252 , there are two characteristic patterns of polarization ellipses around a @xmath1 point : ( a )  _ lemon _ with @xmath253 and ( b )  _ monstar _ with @xmath254  @xcite . different quantitative criteria to distinguish between the @xmath1 points of the lemon and the monstar types were deduced in refs . @xcite . from these criteria it can be inferred that a lemon becomes a monstar as it approaches a star and @xmath1 point annihilation occurs only between stars and monstars  @xcite . the case of linearly polarized wave with @xmath255 provides another example of the polarization singularity where the handedness is undefined . the curves along which the polarization is linear are called the _ l lines_.    from eq . the @xmath1 points are nodal points of the scalar complex function @xmath256 which can be found as intersection points of the stokes parameter nodal lines @xmath257 and @xmath258 . similarly , equation   implies that nodes of the stokes parameter field @xmath259 provide the @xmath0 lines where @xmath260 and @xmath255 . now we present the theoretical results for the polarization - resolved patterns of the dhflc cells . these patterns are computed for the cell of thickness @xmath261 @xmath16 m   filled with the flc mixture flc-576a which was studied in ref . @xcite and described at the end of sec . [ subsec : kerr - effect ] . our first remark is that the angular dependence of the elements of the transmission matrix   is determined by the angle difference @xmath262 which is the angle between the in - plane optical axis @xmath263 ( see eq .  ) and the lateral wave vector @xmath264 ( see eq .  ) . then the vector amplitudes    [ eq : tilde_e_phi ] @xmath265    where @xmath266 is the angle between the optical axis @xmath263 and the polarization plane of the linearly polarized incident wave ( see eq .  ) , describing the incident and transmitted waves with polarization ellipses rotated by the angle @xmath267 are related by the transformed transmission matrix @xmath268\,{\mathbf{t}}_{c}(\rho,-\tilde{\phi})\ , \exp[i\tilde{\phi}{\boldsymbol{\sigma}}_3].\end{aligned}\ ] ]    from relation   it follows that , given the angle @xmath266 , the sole effect of changing the azimuthal angle of the optical axis : @xmath269 is the rotation of the polarization ellipse field by the angle @xmath270 . in dhflc cells , this effect manifests itself as the electric field induced rotation and can be clearly seen in fig . [ fig : pol - f-20 ] that shows the patterns emerging after the dhflc cell calculated at @xmath271  deg for two values of the electric filed parameter : @xmath272 ( see fig . [ subfig : pol-00 - 20 ] ) and @xmath273 ( see fig . [ subfig : pol-03 - 20 ] ) . figure  [ fig : pol - f-20 ] also illustrates the case of angular patterns that do not contain @xmath1 points . the geometry of such patterns is completely characterized by the @xmath0 lines . interestingly , at @xmath274  deg , it turned out that the @xmath259 nodal lines can be evaluated using the simplified equation @xmath275 where @xmath207 is the phase retardation expressed in terms of the eigenvalues of the matrix   for uniformly anisotropic planar layers ( see eqs .   and  ) , that thus gives a sufficiently accurate approximation for @xmath0 lines . the angle @xmath266 can be regarded as the governing parameter whose magnitude determines the formation of @xmath1 points . given the aperture , the latter occurs only if the magnitude of @xmath266 exceeds its critical value . the case where the angle @xmath266 is close to the critical value is illustrated in fig . [ subfig : pol-00 - 39_5 ] . it can be seen that , at @xmath276  deg , the singularity structure of the polarization ellipse fields becomes complicated and is characterized by the presence of symmetrically arranged star - monstar pairs of @xmath1 points . in addition to the above discussed electric - field - induced rotation , the electric field is shown to facilitate the formation of @xmath1 points . clearly , the field induced biaxial anisotropy is responsible for this effect .    in fig . [ fig : pol - f-41 ] , we show how the singularity structure of the polarization resolved patterns develops when the angle @xmath266 further increases . this structure can be described as symmetrically arranged chains of star - monstar pairs of @xmath1 points . as is indicated in fig . [ subfig : pol-03 - 41 ] , these tightly packed chains of @xmath1 points are generally located in the vicinity of the lines @xmath277 that give a high accuracy approximation for the @xmath278 nodal lines where @xmath258 and , similar to eq .  , are determined by the phase retardation @xmath279 . note that , since @xmath280}^{\ast } } e_{p}^{({\mathrm{trm}})}$ ] , applicability of approximate formulas   and   implies that the phase difference between the components of the transmitted waves , @xmath281 and @xmath282 , is close to the phase retardation : @xmath283}^{\ast } } e_{p}^{({\mathrm{trm}})}\}\approx \delta$ ] . in sec . [ subsec : kerr - effect ] , a similar approximation has been used to derive the expression for the transmittance given by eq .  . as it was mentioned in the previous section , the loci of @xmath1 points on the projection plane are determined by intersections of the @xmath278 and @xmath284 nodal lines . in our case , the star - monstar pairs are produced as a result of small - scale oscillations of the nodal line around the smooth curve described by eq .  . experimentally , it is a challenging task to resolve accurately the chains of @xmath1 points resulting from such ripplelike oscillations in polarimetry measurements . the 2d polarization - resolved patterns are centrally symmetric being invariant under inversion through the origin : @xmath285 . the reason is that optical properties of planar anisotropic structures are unchanged under a 180-degree rotation about the normal to the cell ( the @xmath25 axis ) . more specifically , we have the central symmetry relation @xmath286 which is an immediate consequence of the fact that the matrix @xmath34 given in eqs .   and   remains intact when the azimuthal angle of the in - plane optic axis is changed by @xmath240 . another symmetry relation @xmath287 describes the transformation of the transmission matrix   under reflection in the mirror symmetry axis directed along @xmath263 : @xmath288 . this relation immediately follows from eq . deduced in appendix  [ sec : planar ] . by using formula  , it is not difficult to show that the polarization ellipse field @xmath289 transforms into its mirror symmetric counterpart @xmath290 when the polarization azimuth @xmath266 changes its sign : @xmath291 . in this paper , we have performed transfer matrix analysis of polarization - resolved angular patterns emerging after electrically controlled short - pitch dhflc cells . our formulation of the transfer matrix method , which is a suitably modified version of the approach developed in refs . @xcite , involves the following steps : ( a )  derivation of the system of equations for the tangential components of the wave field in the @xmath18 matrix form ( see eq .  ) ; ( b )  introducing the evolution operator ( propagator )   and the scattering matrix  ; ( c )  defining the transfer matrix   through the propagator and , finally , ( d )  deducing formulas   that link the transfer and scattering matrices . description of this method is augmented by discussion of a variety of unitarity and symmetry relations ( see sec . [ subsec : symmetry ] and appendix  [ sec : planar ] ) , with an emphasis on the special case of anisotropic planar structures representing homogenized dhflc cells . interestingly , the relations given in eqs .  ,   and   are shown to be essentially independent of the assumption of lossless materials and thus can be used when the medium is absorbing .    in general , we found that , owing to its mathematical structure , the transfer matrix approach provides the framework particularly useful for in - depth analysis of symmetry related properties ( recent examples of such analysis can be found , e.g. , in refs . @xcite ) . similarly , one of the important results of a rigorous analysis performed within such a framework in ref . @xcite is the expression for the effective dielectric tensor   describing the electro - optical properties of uniform lying flc helical structures with subwavelengh pitch .    in sec . [ sec : optics - dhf ] we have extended theoretical considerations of ref . @xcite to the case of biaxial flcs and have applied an alternative technique of averaging over distorted helix to evaluate the dielectric tensor . this technique is presented in appendix  [ sec : flc - spiral ] and gets around the difficulties of the method that relies on the well - known first - order expression for a weakly distorted helix  . the modified averaging procedure allows high - order corrections to the dielectric tensor to be accurately estimated and improves agreement between the theory and the experimental data in the high - field region ( see fig . [ fig : expt_th ] ) . the resulting electric field dependence of the effective dielectric tensor   is linear ( quadratic ) for non - diagonal ( diagonal ) elements with the coupling coefficients given by eq .  . these coupling coefficients along with the zero - field dielectric constants   determine how the applied electric field @xmath292 changes the principal values of the effective dielectric tensor ( see eqs .   and  ) and the azimuthal angle of optical axis  . generally , at @xmath122 , the dhflc cell is uniaxially anisotropic with the dielectric constants  , @xmath144 and @xmath143 , and the optical axis directed along the helix axis . then there are two most important effects induced by the electric field : ( a )  producing biaxial anisotropy by changing the eigenvalues of the dielectric tensor ; and ( b )  rotation of in - plane optical axes by the field dependent angle @xmath157 defined in eq  . at sufficiently low electric field @xmath292 and non - vanishing zero - field anisotropy , the kerr - like regime takes place so that the principal values depend on the electric field quadratically whereas the optical axis angle @xmath157 is approximately proportional to @xmath292 . this is the orientational kerr effect previously studied in refs . @xcite for different geometries .    our results on dependence of the coupling coefficients and the zero - field dielectric constants on the smectic tilt angle @xmath3 described in sec . [ subsec : tilt ] indicate a number of differences between uniaxial and biaxial flcs . what is more important , they show that the zero - field anisotropy may vanish at certain value of @xmath3 which might be called the isotropization angle : @xmath183 ( see fig . [ fig : eph_theta ] ) .    in sec . [ subsec : isotropy ] , the isotropization point determined by the condition of zero - field isotropy   is found to represent a singularity known as the exceptional point  @xcite . for analytic continuation of the dielectric tensor   in the complex @xmath175 plane , the exceptional points occur at the zeros of the square root in eq . where @xmath293 ^ 2+[\gamma_{xy}\alpha_e]^2=0 $ ] . in general , there are two pairs of complex conjugate values of electric field parameter representing four exceptional ( branch ) points . when the difference @xmath294 vanishes , the two branch points coalesce on the real axis at the origin .    in the case of conventional uniaxial flcs , the analytic solution of the condition of zero - field isotropy can be obtained in the closed form and is given by simple formula   where the isotropization angle , @xmath186 , is found to be a decreasing function of the anisotropy parameter @xmath295 . for biaxial flcs with @xmath296 , the solution can only be written in the parametrized form  . as it can be seen in fig . [ fig : theta_r1 ] , the corresponding @xmath186-@xmath167 curves are splitted into two branches separated by the gap . these results significantly differ from the relation @xmath297 that can be easily obtained  @xcite for the dielectric tensor   averaged over the flc helix : @xmath298{\mathbf{i}}_3 + [ u_1\cos^2\theta - ( u_1\sin^2\theta+u_2)/2]{\hat{\mathbf{x}}}\otimes{\hat{\mathbf{x}}}$ ] . the difference stems from the fact that , in our approach , the effective dielectric tensor @xmath299 is defined through the averaged differential propagation matrix @xmath300 and thus is not equal to the averaged dielectric tensor  : @xmath301    at the exceptional point , the kerr - like regime breaks down and the electric field dependence of the birefringence becomes linear ( see fig . [ fig : dn_alp ] ) . this might be called the pockels - like regime which is characterized by the harmonic electric field dependence of the transmittance of light passing through crossed polarizers ( see fig . [ subfig : transm_c_uni ] ) . the curve shown in fig . [ subfig : transm_c-2_uni ] illustrates the electro - optical response of a dhflc cell near the exceptional point . it is seen that sensitivity to the electric field and the magnitude of the transmittance at peaks are both considerably enhanced as compared to the case studied in ref . @xcite ( see also fig .  [ fig : expt_th ] ) . we now try to put these results in a more general physical context . in quantum physics , the exceptional points are known to produce a variety of interesting phenomena including level repulsion and crossing , bifurcation , chaos and quantum phase transitions  @xcite . for optical wave fields , a recent example is unidirectional propagation ( reflection ) of light at the exceptional points in parity - time ( @xmath302 ) symmetric periodic structures and metamaterials that has been the subject of intense studies  @xcite .    to the best of our knowledge , the role of exceptional points in optics of liquid crystal systems has yet to be recognized . the main problem with conventional uniaxial flcs is that the isotropy condition   requires large values of the smectic tilt angle that are typically well above @xmath303  deg . though there are no fundamental limitations preventing preparation of flc mixtures with large tilt angles , this task still remains a challenge to deal with in the future . biaxial flcs , where the isotropization tilt angle can be sufficiently small when @xmath304 is close to @xmath305 , also present a promising alternative approach for future work .    in sec . [ sec : conoscopy ] , in order to gain further insight into the electro - optical properties of the dhflc cells , we have combined the transfer matrix approach and the results for the effective dielectric tensor to explore the polarization - resolved angular patterns which are the polarization ellipse fields representing the polarization structure of conoscopic images of dhflc cells . in the observation plane , such 2d patterns encode information on how the polarization state of transmitted light is changed with the incidence angles and exhibit singularities of a different kind , the polarization singularities such as @xmath0 lines ( lines of linear polarization ) and @xmath1 points ( points of circular polarization ) . note that , similar to the above discussed exceptional point at which the angle @xmath157 becomes undetermined , @xmath1 points can be regarded as phase singularities ( optical phase singularities are reviewed in ref . @xcite ) .    since the differential propagation matrix of planar structures is invariant under rotation of in - plane optical axes by @xmath240 , the patterns are centrally symmetric [ see figs .  [ fig : pol - f-20]-  [ fig : pol - f-41 ] ] . it was shown that , at fixed the angle @xmath306 between the optical axis @xmath263 and the polarization plane of incident wave , the sole effect of the electric - field - induced rotation by the angle @xmath157 is rotation of the polarization ellipse field as a whole by the same angle . the symmetry axis of the @xmath259 nodal lines ( @xmath0 lines ) is found to be directed along @xmath263 . when the @xmath266 is not too small , they can be approximated by solving eq .   and thus are mainly determined by the phase retardation @xmath207 . similar remark applies to the @xmath278 lines and approximate formula  . it turned out that this is the angle @xmath266 that plays the role of the parameter governing formation of @xmath1 points . when the magnitude of @xmath266 exceeds its critical value which , in our case , is close to @xmath307  deg , @xmath1 points emerge as symmetrically arranged and densely packed chains of star - monstar pairs ( see figs . [ fig : pol - f-39_5]-[fig : pol - f-41 ] ) . so , in dhflc cells , rotation of polarization ellipse fields and formation of @xmath1 points are two most important effects describing electrically induced transformations of the polarization - resolved angular patterns . these predictions can be verified experimentally . this work is now in progress . this work is supported by the hkust grant itp/039/12np . we begin with the relation @xmath308 known as the _ composition law_. this result derives from the fact that the operator @xmath309 is the solution of the system   that satisfies the initial condition   with @xmath310 replaced by @xmath311 .    from the composition law it immediately follow that the inverse of the evolution operator is given by @xmath312 and can be found by solving the initial value problem @xmath313    for non - absorbing media with symmetric dielectric tensor , the matrix @xmath34 is real - valued , @xmath314 , and meets the following symmetry identities  @xcite : @xmath315 where an asterisk and the superscript @xmath91 indicate complex conjugation and matrix transposition , respectively . in this case , the evolution operator and its inverse are related as follows : @xmath316 where a dagger will denote hermitian conjugation . by using the relations  , it is not difficult to verify that the operator on the right hand side of eq . is the solution of the cauchy problem  . in this section we present the results for anisotropic planar structures characterized by the dielectric tensor of the following form : @xmath317 where the optical axes @xmath318 lie in the plane of substrates ( the @xmath160-@xmath319 plane ) . the operator of evolution can be expressed in terms of the eigenvalue and eigenvector matrices , @xmath320 and @xmath321 , as follows @xmath322 for the dielectric tensor  , @xmath323 and @xmath324 ( see eq .  ) . assuming that @xmath325 , we have @xmath326 where @xmath327 and @xmath328 are the principal refractive indices ; @xmath329 is the parameter of in - plane anisotropy . for the case where the diagonal block - matrices , @xmath330 and @xmath331 , vanish , it is not difficult to show that the eigenvector and eigenvalue matrices can be taken in the following form : @xmath332 in addition , the eigenvectors satisfy the orthogonality conditions ( a proof can be found , e.g. , in appendix a of ref . @xcite ) that , for the eigenvector matrix of the form  , can be written as follows @xmath333    upon substituting eqs .  -   into eq . , some rather straightforward algebraic manipulations give the transfer matrix @xmath334 { \mathbf{n}}^{-1 } , \\ &    \label{eq : a_pm - pln }    { \mathbf{a}}_{\pm}={\mathbf{e}}_{{\mathrm{m } } } { \mathbf{h}}\pm{\mathbf{h}}_{{\mathrm{m } } } { \mathbf{e}},\end{aligned}\ ] ] where @xmath335 . from eq . , the block @xmath35 matrices of @xmath336 are given by    [ eq : tw_ij - pln ] @xmath337}^{t}}= \notag \\ & = { \mathbf{a}}_{- } { \mathbf{w}}_{- } { \mathbf{a}}_{+}^{t}-{\mathbf{a}}_{+ } { \mathbf{w}}_{+ } { { { \mathbf{a}}}^{t}}_{-}.   \end{aligned}\ ] ]    finally , we can combine eq . and eq .   with eq .   to derive the expressions for the transmission and reflection matrices @xmath338 describing the case where the incident wave is impinging onto the entrance face of the layer , @xmath23 .    as it can be seen from formulas  , the symmetry relations   are satisfied even if the dielectric constants @xmath339 , @xmath340 and @xmath173 are complex - valued . so , the applicability range of identities   includes lossy ( absorbing ) anisotropic materials described by the dielectric tensor of the form given in eq .  . interestingly , inverse of the transfer matrix , @xmath341 , can be obtained from formula   by changing sign of the thickness parameter @xmath342 : @xmath343 . in formulas  , this transformation interchanges the matrices @xmath344 and @xmath345 , so that @xmath346 and @xmath347 . so , from eq .  , the block matrices of @xmath341 are given by    [ eq : w_inv_pln ] @xmath348    from the other hand , eqs .   and give the transfer matrix and its inverse , respectively , expressed in terms of the transmission and reflection matrices , @xmath349 and @xmath350 . these expressions can now be substituted into eq .   to yield the relations @xmath351 linking the transmission ( reflection ) matrix , @xmath352 ( @xmath353 ) , and its mirror symmetric counterpart @xmath69 ( @xmath70 ) .    in conclusion of this section , we consider how the transmission and reflection matrices transform under the reflection in the @xmath354 plane when the azimuthal angle @xmath355 changes its sign : @xmath356 . from eqs .   and , we have @xmath357 by using eq . it is not difficult to deduce a similar relation for the transfer matrix @xmath358 that can be combined with eq .   to yield the result for the transmission and reflection matrices in the final form : @xmath359 an important point is that , similar to identities  , the assumption of lossless ( non - absorbing ) medium is not required to derive the symmetry relations   and eq .  .      for the case of uniaxially anisotropic structure with @xmath360 , it is not difficult to find the expressions for the eigenvalues that enter the eigenvalue matrix   @xmath361 where @xmath362 ( @xmath363 ) is the refractive index for ordinary ( extraordinary ) waves and @xmath329 is the anisotropy parameter . similarly , after computing the eigenvectors , we obtain the eigenvector matrix in the following form : @xmath364 & m_yq_o \\ m_y & -m_xq_o \end{pmatrix } , \quad { \mathbf{h}}= \begin{pmatrix }                m_x q_e & m_y n_o^2 \\               m_y q_e & -m_x [ n_o^2-q_p^2 ] \end{pmatrix } , \\ & \label{eq : n_oe - pln - uni } n_{e}=\dfrac{2q_e\mu}{n_o^2}(n_o^2-q_p^2m_x^2 ) , \quad n_{o}=2q_o\mu ( n_o^2-q_p^2m_x^2).\end{aligned}\ ] ] equations  -   can now be substituted into the general expression for the transfer matrix defined by formulas  -   so as to obtain the transmission and reflection matrices  .      for the general case of biaxial anisotropy , the expressions for the eigenvalues are more complicated than those for uniaxially anisotropic layer ( see eq .  ) . these can be written in the following form : @xmath365 ^ 2 - 4\det{\tilde{{\mathbf{m } } } }   } , \end{aligned}\ ] ] where the matrix @xmath366 is given by    [ eq : tilde - m - pln ] @xmath367 \cdot{\mathbf{rt}}(\psi)= \begin{pmatrix }    \tilde{m}_{11 } & \tilde{m}_{12}\\ \tilde{m}_{21 } & \tilde{m}_{22 } \end{pmatrix } , \\ & \label{eq : tilde - m-11 }        \tilde{m}_{11}=n_{e}^2-(u_a+u_z[1+u_a])\,q_p^2 m_x^2 , \quad        \tilde{m}_{12}=(u_a+u_z[1+u_a])\ , q_p^2 m_x m_y , \\ & \label{eq : tilde - m-22 }       \tilde{m}_{22}=n_{o}^2-u_z q_p^2 m_y^2 , \quad        \tilde{m}_{21}=u_z q_p^2 m_x m_y,\quad u_z=(\epsilon_{\perp}-\epsilon_z)/\epsilon_z.\end{aligned}\ ] ]    it can be readily checked that the result for uniaxial anisotropy   is recovered from eq .   as the limiting case where the parameter of out - of - plane anisotropy @xmath368 is negligible and @xmath360 .    for the eigenvector and normalization matrices , @xmath321 and @xmath369 , given in eq .   and eq .  , respectively , the results are @xmath370 these relations along with formulas   give the transfer matrix for biaxially anisotropic films with two in - plane optical axes .    before closing this section we briefly comment on the important special case of normal incidence that occurs at @xmath371 . in this case , the matrices @xmath372 defined in eq .   can be written in the factorized form @xmath373 where @xmath374 is the matrix describing rotation about the @xmath25 axis by the angle @xmath216 . substituting eq .   into eq . gives the block matrices @xmath375 expressed as a function of the director azimuthal angle @xmath355 .    the result for the transmission and reflection matrices @xmath376 where the diagonal matrices @xmath377 and @xmath378 describe the case in which the director   lies in the incidence plane , immediately follows from the relations   and  . in sec .  [ subsec : diel - tensor - dhf ] , the effective dielectric tensor   of a deformed helix flc cell is expressed in terms of the averages given in eqs .   and  . in this apppendix , we describe how to perform averaging over the helix pitch without recourse to explicit formulas for the azimuthal angle @xmath8 the flc director  . we assume that the azimuthal angle is a function of @xmath160 , so that the free energy density can be written in the following form : @xmath379 where @xmath124 is the free twist wave number . then the free energy functional per unit volume can be written as the free energy density averaged over the helix pitch @xmath380/v={\ensuremath{\langle{f}\rangle}}_{x}\equiv \frac{1}{p}\int_0^p f { \mathrm{d}}x.\end{aligned}\ ] ] the first integral of the stationary point ( euler - lagrange ) equation @xmath381 is given by @xmath382 ^ 2 + u_{e}\cos\phi = e.\end{aligned}\ ] ] assuming that @xmath383 and @xmath384 is non - negative , equation   can be recast into the differential form @xmath385 integrating eq .   over the period yields the relation for the helix wave number @xmath386 where @xmath387 . this relation gives the helix pitch , @xmath7 , expressed in terms of the dimensionless parameter @xmath388 . we can now use equation   to rewrite eq .   in the following form @xmath389^{-1/2}.\end{aligned}\ ] ] an important consequence of this equation is the relation @xmath390 that allows to perform averaging over the helix pitch by computing integrals over the azimuthal angle @xmath8 . in particular , with the help of eqs .   and  , it is not difficult to deduce the following expression for the free energy  : @xmath391)= \notag \\ u_k\left [ 1-m^{-1}\bigl\ { 1 + 2 ( \sqrt{m}-{\ensuremath{\langle{r^{-1}}\rangle}}_{\phi})/{\ensuremath{\langle{r}\rangle}}_{\phi } \bigr\ } \right].\end{aligned}\ ] ]    in the low voltage regime , the parameter @xmath392 is small and the left hand side of eq . can be expanded into the power series in @xmath392 . the expansion up to the second order terms is given by @xmath393 and can be used to average the @xmath25 component of the polarization vector @xmath394 defined in eq .  . the result reads @xmath395 where @xmath151 is the dielectric susceptibility of the goldstone mode  @xcite and @xmath396 . similarly , the averages that enter the formulas for the elements of the effective dielectric tensor   can be expressed in terms of the electric field parameter @xmath175 as follows :    [ eq : averages - planar ] @xmath397^{-1/2}(1 + 3 v\gamma_v^2\alpha_e^2 ) , \\ & { \ensuremath{\langle{v_{zz}^{-1}\cos^2\phi}\rangle}}_{\phi } \approx \gamma_v(1 + 3 [ 1+v]^{1/2}\gamma_v\alpha_e^2 ) , \\ & { \ensuremath{\langle{v_{zz}^{-1}\cos\phi}\rangle}}_{\phi } \approx 2 \gamma_v\alpha_{e } , \quad \gamma_v= [ \sqrt{1+v}+1]^{-1 } , \\    & { \ensuremath{\langle{v_{zz}^{-1}\sin\phi}\rangle}}_{\phi}= { \ensuremath{\langle{v_{zz}^{-1}\sin\phi\cos\phi}\rangle}}_{\phi}= 0.\end{aligned}\ ] ]    substituting the relations   into eqs . give the effective dielectric tensor   which is expressed in terms of the zero - field dielectric constants @xmath398 , \\ & \label{eq : epsilon_p - v } \epsilon_p = \epsilon_2\sqrt{1+v } , \quad \gamma_v=[\sqrt{1+v}+1]^{-1},\end{aligned}\ ] ] and the coupling coefficients @xmath399 , \\ & \gamma_{yy } =3 \epsilon_2 v \gamma_v^2\sqrt{1+v } , \quad \gamma_{xy } = 2 ( \epsilon_{1}-\epsilon_{\perp})\gamma_v\cos\theta\sin\theta.\end{aligned}\ ] ]    45ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]  + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty @noop _ ( ,  ,  ) p.   @noop _ _ ,  the liquid crystals book series ( ,  ,  )  p.   @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * , ( ) @noop _ _  ( , ,  )  p. @noop _ _ ,   ed . ( ,  ,  )  p.   @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop _ _  ( , ,  ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop _ _  ( ,  ,  )  p. @noop * * ,   ( ) @noop * * ,   ( ) @noop * * , ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop _ _ ,   ed . ,  , vol . ( ,  , )  p.   @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) link:\doibase 10.1007/978 - 94 - 007 - 1530 - 1 [ _ _ ] ,  ed . ( ,  ,  )  p. link:\doibase 10.1109/tap.2012.2220316 [ * * ,   ( ) ] @noop * * , ( ) @noop * * ,   ( ) link:\doibase 10.1103/physrevlett.105.053902 [ * * , ( ) ] link:\doibase 10.1103/physreva.85.064103 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.106.213901 [ * * ,   ( ) ] http://dx.doi.org/10.1038/nmat3576 [ * * ,   ( ) ] link:\doibase    10.1103/physreva.89.065801 [ * * ,   ( ) ]
Abstract:

Answer: in order to explore electric - field - induced transformations of polarization singularities in the polarization - resolved angular ( conoscopic ) patterns emerging after deformed helix ferroelectric liquid crystal ( dhflc ) cells with subwavelength helix pitch , we combine the transfer matrix formalism with the results for the effective dielectric tensor of biaxial flcs evaluated using an improved technique of averaging over distorted helical structures . within the framework of the transfer matrix method , we deduce a number of symmetry relations and show that the symmetry axis of @xmath0 lines ( curves of linear polarization ) is directed along the major in - plane optical axis which rotates under the action of the electric field . when the angle between this axis and the polarization plane of incident linearly polarized light is above its critical value , the @xmath1 points ( points of circular polarization ) appear in the form of symmetrically arranged chains of densely packed star - monstar pairs . we also emphasize the role of phase singularities of a different kind and discuss the enhanced electro - optic response of dhflcs near the exceptional point where the condition of zero - field isotropy is fulfilled .