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: since einstein formulated special relativity in 1905 , the basic space - time symmetry has been that of the lorentz group . he established the energy - momentum relation which is valid for slow massive particles , high - speed massive particles , and massless particles . einstein formulated this relation initially for particles without internal space - time structures , but it is widely accepted that it is valid for all particles including particles with spin and/or internal space - time extension .    in 1939 , wigner published his most fundamental paper dealing with internal space - time symmetries of relativistic particles  @xcite . in this paper , wigner introduced the lorentz group to physics . furthermore , by introducing his `` little groups , '' wigner provided the framework for studying the internal space - time symmetries of relativistic particles . the scientific contents of this paper have not yet been fully recognized by the physics community . we are writing this report as a continuation of the work wigner initiated in this history - making paper .    while particle physicists are still struggling to understand internal space - time symmetries of elementary particles , wigner s lorentz group is becoming useful to many other branches of physics . among them are optical sciences , both quantum and classical . in quantum optics , the coherent and squeezed states are representations of the lorentz group  @xcite . recently , the lorentz group is becoming the fundamental language for classical ray optics . it is gratifying to note that optical components , such as lenses , polarizers , interferometers , lasers , and multi - layers can all be formulated in terms of the lorentz group which wigner formulated in his 1939 paper . classical ray optics is of course a very old subject , but we can not do new physics without measurements using optical instruments . indeed , classical ray optics constitutes wigner s frontier in physics .    the word `` group theory '' sounds like abstract mathematics , but it is gratifying to note that wigner s little groups can be formulated in terms of two - by - two matrices , while classical ray optics is largely a physics of two - by - two matrices . the mathematical correspondence is straight - forward .    in order to see this point clearly , let us start with the following classic example . the second - order differential equation @xmath0 is applicable to a driven harmonic oscillator with dissipation . this can also be used for studying an electronic circuits consisting of inductance , resistance , capacitance , and an alternator . thus , it is possible to study the oscillator system using an electronic circuit . likewise , an algebra of two - by - two matrices can serve as the scientific language for several different branches of physics , including special relativity , ray optics , and quantum optics .    rccc + & & & + & massive , slow & covariance & massless , fast + & & & + energy- & & einstein s & + momentum & @xmath1 & @xmath2 & @xmath3 + & & & + internal & @xmath4 & & @xmath5 + space - time & & wigner s & + [ -1 mm ] symmetry & @xmath6 & little group & gauge trans . + & & & + relativistic & & lorentz - squeezed & + extended & quarks & harmonic & partons + [ -1 mm ] particles & & oscillators & +    there are many physical systems which can be formulated in terms of two - by - two matrices . if we restrict that their determinant be one , there is a well established mathematical discipline called the group theory of @xmath7 . this aspect was noted in the study of lorentz transformations . in group theoretical terminology , the group @xmath7 is the universal covering group of the group of the lorentz group . in practical terms , to each two - by - two matrix , there corresponds one four - by - four matrix which performs a lorentz transformation in the four - dimensional minkowskian space . thus , if a physical system can be explained in terms of two - by - two matrices , it can be explained with the language of lorentz transformations . furthermore , the physical system based on two - by - two matrices can serve as an analogue computer for lorentz transformations .    optical filters , polarizers , and interferometers deal with two independent optical rays . they superpose the beams , change the relative phase shift , and change relative amplitudes . the basic language here is called the jones matrix formalism , consisting of the two - by - two matrix representation of the @xmath7 group  @xcite . the four - by - four mueller matrices are derivable from the two - by - two matrices of @xmath7 . for these two - beam systems , the poincar sphere serves as an effective language  @xcite . since the two - beam system is described by the lorentz group , the poincar sphere is necessarily a representation of the lorentz group . we shall use this sphere to study the degree of coherence between the beams . para - axial lens optics can also be formulated in terms of two - by - two matrices , applicable to the two - component vector space consisting of the distance from the optical axis and the slope with respect to the axis . the lens and translation matrices are triangular , but they are basically representations of the @xmath8 group which is the real subgroup of the group @xmath7  @xcite . laser optics is basically multi - lens lens optics . however , the problem here is how to get a simple mathematical expression for the system of a large number of the same lens separateled by equal distance . here again , group theory simplifies calculations  @xcite .    in multi - layer optics , we deal with two optical rays moving in opposite directions . the standard language in this case is the s - matrix formalism  @xcite . this is also a two - by - two matrix formalism . as in the case of laser cavities , the problem is the multiplication of a large number of matrix chains  @xcite . it is shown in this report that the two - by - two representation of the six - parameter lorentz group is the underlying common scientific language for all of the instruments mentioned above . while the abstract group theoretical ideas make two - by - two matrix calculations more systematic and transparent in optics , optical instruments can act as analogue computers for lorentz transformations in special relativity . it is gratifying to note that special relativity and ray optics can be formulated as the physics of two - by - two matrices .    in sec . [ lorentz ] , we explain how the lorentz group can be formulated in terms of four - by - four matrices . it is shown that the group can have six independent parameters . in sec . [ spinor ] , we explain how it is possible to formulate the lorentz group in terms of two - by - two matrices . it is shown that the four - by - four transformation matrices can be constructed from those those of the two - by - two matrices .    in sec . [ wlittle ] , we discuss the historical significance of wigner s 1939 paper  @xcite on the lorentz group and its application to the internal space - time symmetries of relativistic particles . in sec . [ polari ] , we present the basic building blocks for the two - by - two representation of the lorentz group in terms of the matrices commonly seen in ray optics .    in secs . [ polari ] and [ interfer ] , we study the polarizations for the two - beam system . it is shown that both the jones and mueller matrices are representations of the lorentz group . the role of the stokes parameters is discussed in detail . it is shown in sec . [ poinc ] that the poincar sphere is a graphical representation of the poincar group , and serves as a device to describe the degree of coherence .    in secs . [ olens ] , [ mlens ] , [ lcav ] , and [ mlayer ] , we discuss one - lens system , multi - lens system , laser cavities , and multi - layer optics , respectively . in all of these sections , the central scientific language is the lorentz group .    in appendix [ sqosc ] , we expand the content of the third row of table  [ einwig ] . it is noted that the covariant harmonic oscillator formalism can unify the quark model for slow hadrons with the parton model for ultra - relativistic quarks . it is then shown that the same oscillator formalism serves as the basic scientific language for squeezed states of light .    in appendix [ euler ] , it is shown that the lie - group method , in terms of the generators , is not the only method in constructing group representations . for the rotation group and the three - parameter subgroups of the lorentz group , it is simpler to start with the minimum number of starter matrices . for instance , while there are three generators for the rotation group in the lie approach , we can construct the most general form of the rotation matrix from rotations around two directions , as goldstein constructed the euler angles  @xcite .    in appendix  [ conju ] , it is noted that the four - by - four matrices are real , their two - by - two counterparts are complex . however , there is a three - parameter real subgroup called @xmath8 . it is shown that the complex subgroup @xmath9 is equivalent to @xmath8 through conjugate transformation . let us consider the space - time coordinates @xmath10 . then the rotation around the @xmath11 axis is performed by the four - by - four matrix @xmath12 this transformation is generated by @xmath13 likewise , we can write down the generators of rotations @xmath14 and @xmath15 around the @xmath16 and @xmath17 axes respectively . @xmath18 these three generators satisfy the closed set of commutations relations @xmath19 = i \epsilon_{ijk } j_{k } .\ ] ] this set of commutation relations is for the three - dimensional rotation group . the lorentz boost along the @xmath11 axis takes the form @xmath20 which is generated by @xmath21 likewise , we can write generators of boosts @xmath22 and @xmath23 along the @xmath16 and @xmath17 axes respectively , and they take the form @xmath24 these boost generators satisfy the commutation relations @xmath25 = i \epsilon_{ijk } k_{k } , \qquad \left[k_{i } , k_{j } \right ] = -i \epsilon_{ijk } j_{k } .\ ] ]    indeed , the three rotation generators and the three boost generators satisfy the closed set of commutation relations given in eq.([rotcom ] ) and eq.([locom ] ) . these three commutation relations form the starting point of the lorentz group . the generators given in this section are four - by - four matrices , but they are not the only set satisfying the commutation relations . we can construct also six two - by - two matrices satisfying the same set of commutation relations . the group of transformations constructed from these two - by - two matrices is often called @xmath26 or the two - dimensional representation of the lorentz group . throughout the present paper , we used the two - by - two transformation matrices constructed from the generators of the @xmath26 group . in sec .  [ lorentz ] , we started with four - by - four transformation matrices applicable to the four - dimensional space - time . we then ended up with a set of closed commutation relations for the six generators consisting of three rotation and three boost generators . these generators are in the form of four - by - four matrices .    in this section , we shall see first that there is a set of two - by - two matrices satisfying the same set of commutation relations , constituting the two - by - two representations of the lorentz group . the representation so constructed is called the @xmath26 group , or the universal covering group of the lorentz group . the transformation matrices are applicable to two - component @xmath7 spinors . the algebraic property of this two - by - two representation is the same as that of the four - by - four representation . as we shall see in this paper , the spinors and the four - vectors correspond to jones vectors and stokes parameters respectively in polarization optics  @xcite . the question then is whether we can construct the four - vector from the spinors . in the language of polarization optics , the question is whether it is possible to construct the coherency matrix  @xcite from the jones vector .    with this point in mind , let us start from the following form of the pauli spin matrices : @xmath27 these matrices are written in a different convention . here @xmath28 is imaginary , while @xmath29 is imaginary in the traditional notation . also in this convention , we can construct three rotation generators @xmath30 which satisfy the closed set of commutation relations @xmath31 = i \epsilon_{ijk } j_{k } .\ ] ] we can also construct three boost generators @xmath32 which satisfy the commutation relations @xmath33 = -i \epsilon_{ijk } j_{k } .\ ] ] the @xmath34 matrices alone do not form a closed set of commutation relations , and the rotation generators @xmath35 are needed to form a closed set : @xmath36 = i \epsilon_{ijk } k_{k } .\ ] ]    the six matrices @xmath35 and @xmath34 form a closed set of commutation relations , and they are like the generators of the lorentz group applicable to the ( 3 + 1)-dimensional minkowski space . the group generated by the above six matrices is called @xmath26 consisting of all two - by - two complex matrices with unit determinant . let us write the two - by - two transformation matrix as @xmath37 this matrix has four complex elements with eight real parameters . however , the six generators are all traceless and , the determinant of the matrix has to be one . thus , it has six independent parameters .    in order to construct four - vectors , we need two different spinor representations of the lorentz group . let us go to the commutation relations for the generators given in eqs.([comm1 ] ) , ( [ comm2 ] ) and ( [ comm3 ] ) . these commutators are not invariant under the sign change of the rotation generators @xmath35 , but are invariant under the sign change of the squeeze operators @xmath34 . thus , to each spinor representation , there is another representation with the squeeze generators with opposite sign . this allows us to construct another representation with the generators : @xmath38 we call this representation the `` dotted '' representation . if we write the transformation matrix @xmath39 of eq.([alpha ] ) in terms of the generators as @xmath40 then the transformation matrix in the dotted representation becomes @xmath41 in both of the above matrices , hermitian conjugation changes the direction of rotation . however , it does not change the direction of boosts . we can achieve this only by changing @xmath39 to @xmath42 , and we shall call this the `` dot '' conjugation .    likewise , there are two different set of spinors . let us use write @xmath43 for the spinor in the undotted and dotted representations respectively . then the four - vectors are constructed as  @xcite @xmath44 & { } &    u\dot{v } = ( t + z ) , \quad v\dot{u } = -(t - z)\end{aligned}\ ] ] leading to the matrix @xmath45 where @xmath46 and @xmath47 are one if the spin is up , and are zero if the spin is down , while @xmath48 and @xmath49 are zero and one for the spin - up and spin - down cases .    if the two - by - two matrix of eq . [ alpha ] is applicable to the column vector of eq.([uv ] ) , what is the transformation matrix applicable to the row vector @xmath50 from the right - hand side ? it is the transpose of the matrix applicable to the column vector @xmath50 . we can obtain this column vector from @xmath51 by applying to it the matrix @xmath52 this matrix also has the property @xmath53 where the superscript @xmath54 means the transpose of the matrix . the transformation matrix applicable to the column vector of eq.([dotcol ] ) is @xmath42 of eq.([eldot ] ) . thus the matrix applicable to the row vector @xmath50 in eq.([dotmat ] ) is @xmath55 this is precisely the hermitian conjugate of @xmath39 . let us now write the @xmath56 matrix of eq.([dotmat ] ) as @xmath57 where the set of variables @xmath58 is transformed like a four - vector under lorentz transformations . then the lorentz transformation on @xmath56 can be performed as @xmath59 where the transformation matrix @xmath39 is that of eq.([alpha ] ) . as we have seen in this section , the construction of four - vectors from the two - component spinors is not a trivial task , but has been discussed in the literature  @xcite . likewise , it is possible to construct four - by - four lorentz transformation matrices from the two - by - two matrix of eq.([alpha ] )  @xcite .    as we shall see in the present paper , the matrix of the form eq.([dotmat ] ) or eq.([matv ] ) appears in optics as the coherency matrix and the density matrix for two beam systems such as polarization optics and interferometers . in his 1939 paper , wigner introduced his little groups in order to study the internal space - time structures of relativistic particles . the little group is the maximal subgroup of the lorentz group which leaves the momentum of the particle invariant . the little groups , originally developed for studying relativistic symmetries of particles , are now becoming an important scientific language for ray optics as we shall see in this paper . let us first discuss the physical basis of wigner s little groups . if the speed of a particle is much smaller than that of light , energy - momentum relation is @xmath60 . if the speed is close to that of light , the relation is @xmath3 . these two different relations can be combined into one covariant formula @xmath61 . this aspect of einstein s @xmath62 is well known , as indicated in table  [ einwig ] . in addition , particles have internal space - time variables . massive particles have spins while massless particles have their helicities and gauge variables . our first question is whether this aspect of space - time variables can be unified into one covariant concept . the answer to this question is yes . little group does the job , also as indicated table  [ einwig ] . particles can also have space - time extensions . for instance , in the quark model , hadrons are bound states of quarks . however , the hadron appears as a collection of partons when it moves with speed close to the velocity of light . quarks and partons seem to have quite distinct properties . are they different manifestations of a single covariant entity ? this is one of the most pressing issues in high - energy particles physics . the third row of table  [ einwig ] addresses this question . the mathematical framework of this program was developed by eugene wigner in 1939  @xcite . he constructed the maximal subgroups of the lorentz group whose transformations will leave the four - momentum of a given particle invariant . these groups are known as wigner s little groups . thus , the transformations of the little groups change the internal space - time variables of the particle , while leaving its momentum invariant .    in his paper @xcite , wigner shows that , for each massive particle , there is a lorentz frame in which the particle is at rest . then the three - dimensional rotation group leaves its momentum invariant , while changing the direction of its spin . thus , the little group for this particular case is @xmath63 . in other lorentz frames , the little group is a lorentz - boosted @xmath63 group .    for a massless particle , wigner notes that there are no lorentz frames in which the particle is at rest . the best we can do is to align its momentum to a given axis , or the @xmath11 axis . then , its momentum is invariant under rotations around the @xmath11 axis . in addition , wigner found two more transformations which leave the momentum invariant . the physics of these additional degrees was not explained in wigner s original paper  @xcite , but their mathematics has been worked out . the physics of these degrees of freedom was later determined to be that of gauge transformations  @xcite , and its complete understanding was not achieved until 1990  @xcite . the group of lorentz transformations consists of three boosts and three rotations . the rotations therefore constitute a subgroup of the lorentz group . if a massive particle is at rest , its four - momentum is invariant under rotations . thus the little group for a massive particle at rest is the three - dimensional rotation group . then what is affected by the rotation ? the answer to this question is very simple . the particle in general has its spin . the spin orientation is going to be affected by the rotation ! if we use the four - vector coordinate @xmath10 , the four - momentum vector for the particle at rest is @xmath64 , and the three - dimensional rotation group leaves this four - momentum invariant . this little group is generated by the three rotation generators given in eq.([j33 ] ) and eq.([j12 ] ) . they satisfy the commutation relations for the rotation group given eq.([rotcom ] ) . if the rest - particle is boosted along the @xmath11 direction , it up a non - zero momentum component along the same direction . the @xmath63 generators will also be boosted . the boost matrix takes the form of eq.([boostz ] ) , and the boosted generators will be @xmath65 and this boost will not change the commutation relation of eq.([rotcom ] ) .    for a massless particle moving along the @xmath66 direction , wigner observed that the little group generated by the rotation generator around the @xmath11 axis , namely @xmath67 of eq.([j33 ] ) , and two other generators which take the form @xmath68 if we use @xmath34 for the boost generator along the i - th axis , these matrices can be written as @xmath69 with @xmath22 and @xmath23 given in eq.([k12 ] ) . the generators @xmath70 and @xmath71 satisfy the following set of commutation relations . @xmath72 = 0 , \qquad [ j_{3 } , n_{1 } ] = in_{2 } , \qquad [ j_{3 } , n_{2 } ] = -in_{1 } .\ ] ]    in order to understand the mathematical basis of the above commutation relations , let us consider transformations on a two - dimensional plane with the @xmath73 coordinate system . we can then make rotations around the origin and translations along the @xmath16 and @xmath17 directions . if we write these generators as @xmath74 and @xmath75 respectively , they satisfy the commutation relations  @xcite @xmath76 = 0 , \qquad [ l , p_{x } ] = ip_{y } , \qquad [ l , p_{y } ] = -ip_{x } .\ ] ] this is a closed set of commutation relations for the generators of the @xmath77 group , or the two - dimensional euclidean group . if we replace @xmath78 and @xmath71 of eq.([e2lcom ] ) by @xmath79 and @xmath75 , and @xmath67 by @xmath39 , the commutations relations for the generators of the @xmath77-like little group becomes those for the @xmath77-like little group . this is precisely why we say that the little group for massless particles are like @xmath77 . it is not difficult to associate the rotation generator @xmath67 with the helicity degree of freedom of the massless particle . then what physical variable is associated with the @xmath78 and @xmath71 generators ? indeed , wigner was the one who discovered the existence of these generators , but did not give any physical interpretation to these translation - like generators in his original paper  @xcite . for this reason , for many years , only those representations with the zero - eigenvalues of the @xmath80 operators were thought to be physically meaningful representations  @xcite . it was not until 1971 when janner and janssen reported that the transformations generated by these operators are gauge transformations  @xcite . the role of this translation - like transformation has also been studied for spin-1/2 particles , and it was concluded that the polarization of neutrinos is due to gauge invariance  @xcite . the @xmath63-like little group remains @xmath63-like when the particle is lorentz - boosted . then , what happens when the particle speed becomes the speed of light ? the energy - momentum relation @xmath61 become @xmath81 . is there then a limiting case of the @xmath63-like little group ? since those little groups are like the three - dimensional rotation group and the two - dimensional euclidean group respectively , we are first interested in whether @xmath77 can be obtained from @xmath63 . this will then give a clue as to how to obtain the @xmath77-like little group as a limiting case of @xmath63-like little group . with this point in mind , let us look into this geometrical problem . in 1953 , inn and wigner formulated this problem as the contraction of @xmath63 to @xmath77  @xcite . let us see what they did . we always associate the three - dimensional rotation group with a spherical surface . let us consider a circular area of radius one kilometer centered on the north pole of the earth . since the radius of the earth is more than 6,450 times longer , the circular region appears to be flat . thus , within this region , we use the @xmath77 symmetry group . the validity of this approximation depends on the ratio of the two radii . how about then the little groups which are isomorphic to @xmath63 and @xmath77 ? it is reasonable to expect that the @xmath77-like little group can be obtained as a limiting case of the @xmath63-like little group for massless particles . in 1981 , it was observed by bacry and chang  @xcite and by ferrara and savoy  @xcite that this limiting process is the lorentz boost to infinite - momentum frame .    in 1983 , it was noted by han _ et al _ that the large - radius limit in the contraction of @xmath63 to @xmath77 corresponds to the infinite - momentum limit for the case of the @xmath63-like little group to the @xmath77-like little group for massless particles . they showed that transverse rotation generators become the generators of gauge transformations in the limit of infinite momentum  @xcite . let us see how this happens . the @xmath67 operator of eq.([j33 ] ) , which generates rotations around the @xmath11 axis , is not affected by the boost conjugation of the @xmath82 matrix of eq.([boost ] ) . on the other hand , the @xmath14 and @xmath15 matrices become @xmath83 which are given in eq.([n1n2 ] ) . the generators @xmath78 and @xmath71 are the contracted @xmath15 and @xmath14 respectively in the infinite - momentum . in 1987 , kim and wigner studied this problem in more detail and showed that the little group for massless particles is the cylindrical group which is isomorphic to the @xmath77 group  @xcite . this completes the second row in table  [ einwig ] , where wigner s little group unifies the internal space - time symmetries of massive and massless particles . the transverse components of the rotation generators become generators of gauge transformations in the infinite - momentum limit . let us go back to table i given in sec . [ intro ] . as for the third row for relativistic extended particles , the most efficient approach is to construct representations of the little groups using the wave functions which can be lorentz - boosted . this means that we have to construct wave functions which are consistent with all known rules of quantum mechanics and special relativity . it is possible to construct harmonic oscillator wave functions which satisfy these conditions . we can then take the low - speed and high - speed limits of the covariant harmonic oscillator wave functions for the quark model and the parton model respectively . this aspect was extensively discussed in the literature  @xcite , and is beyond the scope of the present report .    however , it is important to note that the covariant harmonic oscillator formalism use for this purpose can serve as the fock space description for the squeezed state of light  @xcite , we give a brief discussion of this aspect in appendix  [ sqosc ] , entitled `` squeezed harmonic oscillators . ''    in this section , we discussed wigner s little groups applicable to internal space - time symmetries of relativistic particles . however , as we shall see in this paper , these little groups play an important role in understanding ray optics . conversely , optical configurations in ray optics can serve as analogue computers for space - time symmetries in particle physics . let us consider two optical beams propagating along the @xmath11 axis . we are then led to the column vector : @xmath84 we can then achieve a phase shift between the beams by applying the two - by - two matrix : @xmath85 if we are interested in mixing up the two beams , we can apply @xmath86 to the column vector . if the amplitudes become changed either by attenuation or by reflection , we can use the matrix @xmath87 for the change . in this paper , we are dealing only with the relative amplitudes , or the ratio of the amplitudes . as we shall see in sec . [ interfer ] , the above two - by - two matrices have their corresponding four - by - four matrices respectively . repeated applications of these matrices lead to the form @xmath88 where the elements are in general complex numbers . the determinant of this matrix is one . thus , the matrix can have six independent parameters . this matrix takes the identical form as the two - by - two matrix given in eq.([alpha2 ] ) . however , the construction process is different . the lie - group generators are used for eq.([alpha2 ] ) , while eq.([alpha2 ] ) is constructed from repeated applications of the transformation matrices . the difference between these two methods is discussed in appendix  [ euler ] .    in either case , this matrix is the most general form of the matrices in the @xmath26 group , which is known to be the universal covering group for the six - parameter lorentz group . this means that , to each two - by - two matrix of @xmath26 , there corresponds one four - by - four matrix of the group of lorentz transformations applicable to the four - dimensional minkowski space  @xcite . it is possible to construct explicitly the four - by - four lorentz transformation matrix from the parameters @xmath89 and @xmath90 . this expression is available in the literature  @xcite , and we consider here only special cases .    the four - by - four representation of the lorentz group can be constructed from the two - by - two representation of the @xmath26 group , which is known as the universal covering group of the six - parameter lorentz group . this aspect is discussed in sec . [ spinor ] . let us go back to eq.([alpha2 ] ) , the @xmath26 group represented by this matrix has many interesting subgroups . if the matrices are to be hermitian , then the subgroup is @xmath91 corresponding to the three - dimensional rotation group . if all the elements are real numbers , the group becomes the three - parameter @xmath8 group . this subgroup is equivalent to @xmath9 which is the primary scientific language for squeezed states of light  @xcite . we can also consider the matrix of eq.([alpha2 ] ) when one of its off - diagonal elements vanishes . then , it takes the form @xmath92 where @xmath93 is a complex number with two real parameters . in 1939 @xcite , wigner observed this form as one of the subgroups of the lorentz group . he observed further that this group is isomorphic to the two - dimensional euclidean group , and that its four - by - four equivalent can explain the internal space - time symmetries of massless particles including the photons .    in ray optics , we often have to deal with this type of triangular matrices , particularly in lens optics and stability problems in laser and multi - layer optics . in the language of mathematics , dealing with this form is called the iwasawa decomposition  @xcite . if the jones matrix contains all the parameters for the polarized light beam , why do we need the mathematics in the four - dimensional space ? the answer to this question is well known . in addition to the basic parameter given by the jones vector , the stokes parameters give the degree of coherence between the two rays . let us go back to the jones spinor of eq.([jones ] ) , and construct the quantities : @xmath94 & { } &   s_{21 } = < \psi_{2}^{*}\psi_{1 } > , \qquad    s_{22 } = < \psi_{2}^{*}\psi_{2 } >   .\end{aligned}\ ] ] then the stokes vector consists of @xmath95 & { } &   s_{2 } = s_{12 } + s_{21 }   , \quad    s_{3 } = -i(s_{11 } + s_{22 } )    .\end{aligned}\ ] ] the four - component vector @xmath96 transforms like the space - time four - vector @xmath10 under lorentz transformations . the mueller matrix is therefore like the lorentz - transformation matrix .    as in the case of special relativity , let us consider the quantity @xmath97 then @xmath98 is like the mass of the particle while the stokes four - vector is like the four - momentum . if @xmath99 , the two - beams are in a pure state . as @xmath98 increases , the system becomes mixed , and the entropy increases . if it reaches the value of @xmath100 , the system becomes completely random . it is gratifying to note that this mechanism can be formulated in terms of the four - momentum in particle physics  @xcite . although we can borrow all the elegant mathematical identities of the two - by - two representations of the lorentz group , this formalism does not allow us to describe the loss of coherence within the interferometer system . in order to study this effect , we have to construct the coherency matrix : @xmath101 under the optical transformations discussed in this section , the coherency matrix is transformed as @xmath102 as in the case of the lorentz transformation given in eq.([ldag ] ) .    using this formalisms based on the lorentz group , we can discuss the group theoretical property of polarization optics in detail  @xcite . however , polarization optics is based on two independent beams propagating in the same direction . since interferometers are also based on the same optical system , we can continue our discussion in the following section on interferometers . typically , one beam is divided into two by a beam splitter . we can write the incoming beam as @xmath103 then , the beam splitter can be written in the form of a rotation matrix  @xcite : @xmath104 which transforms the column vector of eq.([expo1 ] ) into @xmath105\exp{\left\{i(kz - \omega t)\right\ } } \cr -[\sin(\theta/2)]\exp{\left\{i(kz - \omega t)\right\ } } } .\ ] ] the first beam @xmath106 of eq.([expo1 ] ) is now split into @xmath106 and @xmath107 of eq.([expo2 ] ) . the intensity is conserved . if the rotation angle @xmath108 is -@xmath109 , the initial beam is divided into two beams of the same intensity and the same phase  @xcite . these two beams go through two different optical path lengths , resulting in a phase difference . if the phase difference is @xmath110 , the phase shift matrix is @xmath111    when reflected from mirrors , or while going through beam splitters , there are intensity losses for both beams . the rate of loss is not the same for the beams . this results in the attenuation matrix of the form @xmath112 with @xmath113 . this attenuator matrix tells us that the electric fields are attenuated at two different rates . the exponential factor @xmath114 reduces both components at the same rate and does not affect the degree of polarization . the effect of polarization is solely determined by the squeeze matrix @xmath115    in the detector or in the beam synthesizer , the two beams undergo a superposition . this can be achieved by the rotation matrix like the one given in eq.([rot22 ] )  @xcite . for instance , if the angle @xmath108 is @xmath116 , the rotation matrix takes the form @xmath117 if this matrix is applied to the column vector of eq.([expo2 ] ) , the result is @xmath118 the upper and lower components show the interferences with negative and positive signs respectively . we have shown previously  @xcite that the four - by - four transformation matrices applicable to the stokes parameters are like lorentz - transformation matrices applicable to the space - time minkowskian vector @xmath10 . this allows us to study space - time symmetries in terms of the stokes parameters which are applicable to interferometers . let us first see how the rotation matrix of eq.([rot22 ] ) is translated into the four - by - four formalism . in this case , @xmath119 the corresponding four - by - four matrix takes the form  @xcite @xmath120    let us next see how the phase - shift matrix of eq.([shif22 ] ) is translated into this four - dimensional space . for this two - by - two matrix , @xmath121 for these values , the four - by - four transformation matrix takes the form  @xcite @xmath122 for the squeeze matrix of eq.([sq22 ] ) , @xmath123 as a consequence , its four - by - four equivalent is @xmath124 if the above matrices are applied to the four - dimensional minkowskian space of @xmath10 , the above squeeze matrix will perform a lorentz boost along the @xmath11 or @xmath125 axis with @xmath100 as the time variable . the rotation matrix of eq.([rot44 ] ) will perform a rotation around the @xmath17 or @xmath5 axis , while the phase shifter of eq.([shif44 ] ) performs a rotation around the @xmath11 or the @xmath125 axis . matrix multiplications with @xmath126 and @xmath127 lead to the three - parameter group of rotation matrices applicable to the three - dimensional space of @xmath128 . the phase shifter @xmath127 of eq.([shif44 ] ) commutes with the squeeze matrix of eq.([sq44 ] ) , but the rotation matrix @xmath126 does not . this aspect of matrix algebra leads to many interesting mathematical identities which can be tested in laboratories . one of the interesting cases is that we can produce a rotation by performing three squeezes  @xcite . another interesting case is a combination of squeeze and rotation matrices which will lead to a triangular matrix with unit diagonal elements . this aspect is known as the iwasawa decomposition and is discussed in detail in ref .  @xcite . according to the definition of the density matrix  @xcite , the coherency matrix of eq.([cocy ] ) is also the density matrix . in this section , we shall discuss the cohrency matrix as the density matrix .    under the influence of the @xmath129 transformation given in eq.([lt22 ] ) , this coherency matrix is transformed as @xmath130 & { } & \hspace{5ex } = \pmatrix{\alpha & \beta \cr \gamma & \delta } \pmatrix{s_{11 } & s_{12 } \cr s_{21 } & s_{22 } } \pmatrix{\alpha^ { * } & \gamma^ { * } \cr \beta^ { * } & \delta^ { * } } , \end{aligned}\ ] ] where @xmath131 and @xmath129 are the density matrix and the transformation matrix given in eq.([cocy ] ) and eq.([lt22 ] ) respectively . according to the basic property of the lorentz group , these transformations do not change the determinant of the density matrix @xmath131 . transformations which do not change the determinant are called unimodular transformations .    as we shall see in this section , the determinant for pure states is zero , while that of mixed states does not vanish . is there then a transformation matrix which will change this determinant within the lorentz group . the answer is no . this is the basic issue we would like to address in this section . if the phase difference between the two waves remains intact , the system is said to in a pure state , and the density matrix can be brought to the form @xmath132 through the transformation of eq.([trans22 ] ) with a suitable choice of the @xmath129 matrix . for the pure state , the stokes four - vector takes the form @xmath133    in order to study the symmetry properties of the density matrix , let us ask the following question . is there a group of transformation matrices which will leave the above density matrix invariant ? in answering this question , it is more convenient to use the stokes four - vector . the column vector of eq.([pure4 ] ) is invariant under the operation of the phase shifter @xmath127 of eq.([shif44 ] ) . in addition , it is invariant under the following two matrices : @xmath134 & { } & f_{2}(v ) = \pmatrix { 1 + v^{2}/2   & - v^{2}/2 & 0   & v   \cr     v^{2}/2   & 1 - v^{2}/2 & 0 & v    \cr     0   & 0 & 1 & 0 \cr     u & -v & 0 & 1 } .\end{aligned}\ ] ] these mathematical expressions were first discovered by wigner in 1939  @xcite in connection with the internal space - time symmetries of relativistic particles . they went through a stormy history , but it is gratifying to note that they serve a useful purpose for studying interferometers where each matrix corresponds to an operation which can be performed in laboratories . the @xmath135 and @xmath136 matrices commute with each other , and the multiplication of these leads to the form @xmath137 this matrix contains two parameters .    let us go back to the phase - shift matrix of eq.([shif44 ] ) . this matrix also leaves the stokes vector of eq.([pure4 ] ) invariant . if we define the `` little group '' as the maximal subgroup of the lorentz group which leaves a stokes vector invariant , the little group for the stokes vector of eq.([pure4 ] ) consists of the transformation matrices given in eq.([shif44 ] ) and eq.([d44 ] ) .    next , if the phase relation is completely random , and the first and second components have the same amplitude , the density matrix becomes @xmath138 here is the question : is there a two - by - two matrix which will transform the pure - state density matrix of eq.([pure22 ] ) into the impure - state matrix of eq.([imp22 ] ) ? the answer within the system of matrices of the form given in eq.([lt22 ] ) is no , because the determinant of the pure - state density matrix is zero while that of the impure - state matrix is @xmath139 . is there a way to deal with this problem ? this problem was addressed in ref .  @xcite . in this section , we restrict ourselves to the unimodular transformation of eq.([trans22 ] ) which preserves the value of the determinant of the density matrix . the stokes four - vector corresponding to the above density matrix is @xmath140 this vector is invariant under both the rotation matrix of eq.([rot44 ] ) and the phase shift matrix of eq.([shif44 ] ) . repeated applications of these matrices lead to a three - parameter group of rotations applicable to the three - dimensional space of @xmath128 . not all the impure - state density matrices take the form of eq.([imp22 ] ) . in general , if they are brought to a diagonal form , the matrix takes the form @xmath141 and the corresponding stokes four - vector is @xmath142 with @xmath143 the matrix which transforms eq.([imp4 ] ) to eq.([impp4 ] ) is the squeeze matrix of eq.([sq44 ] ) . the question then is whether it is possible to transform the pure state of eq.([pure4 ] ) to the impure state of eq.([impp4 ] ) or to that of eq.([imp4 ] ) .    in order to see the problem in terms of the two - by - two density matrix , let us go back to the pure - state density matrix of eq.([pure22 ] ) . under the rotation of eq.([rot22 ] ) , @xmath144 the pure - state density matrix becomes @xmath145    for the present case of two - by - two density matrices , the trace of the matrix is one for both pure and impure cases . the trace of the @xmath146 is one for the pure state , while it is less than one for impure states . the next question is whether there is a two - by - two matrix which will eliminate the off - diagonal elements of the above expression that will also lead to the expression of eq.([impp22 ] ) . in order to answer this question , let us note that the determinant of the density matrix vanishes for the pure state , while it is non - zero for impure states . the lorentz - like transformations of eq.([trans22 ] ) leave the determinant invariant . thus , it is not possible to transform a pure state into an impure state by means of the transformations from the six - parameter lorentz group . then is it possible to achieve this purpose using two - by - two matrices not belonging to this group . we do not know the answer to this question . we are thus forced to resort to four - by - four matrices applicable to the stokes four - vector . we are interested in a transformation which will change the density matrix of eq.([pure22 ] ) to that of eq.([imp22 ] ) . for this purpose , we can use the stokes four - vector consisting of the four elements of the density matrix . the question then is whether it is possible to find a transformation matrix which will transform the pure - state four - vector of eq.([pure4 ] ) to the impure - state four - vector of eq.([imp4 ] ) . mathematically , it is more convenient to ask whether the inverse of this process is possible : whether it is possible to transform the four - vector of eq.([imp4 ] ) to that of eq.([pure4 ] ) . this is known in mathematics as the contraction of the three - dimensional rotation group into the two - dimensional euclidean group  @xcite . let us apply the squeeze matrix of eq.([sq44 ] ) to the four - vector of eq.([imp4 ] ) . this can be written as @xmath147 after an appropriate normalization , the right - hand side of the above equation becomes like the pure - state vector of eq.([pure4 ] ) in the limit of large @xmath148 , as @xmath149 becomes equal to @xmath150 in the infinite-@xmath148 limit . this transformation is from a mixed state to a pure or almost - pure state . since we are interested in the transformation from the pure state of eq.([pure4 ] ) to the impure state of eq.([imp4 ] ) , we have to consider an inverse of the above equation : @xmath151 however , the above equation does not start with the pure - state four - vector . if we apply the same matrix to the pure state matrix , the result is @xmath152 the resulting four - vector is proportional to the pure - state four - vector and is definitely not an impure - state four - vector . the inverse of the transformation of eq.([sqimp ] ) is not capable of bringing the pure - state vector into an impure - state vector . let us go back to eq.([sqimp ] ) , it is possible to bring an impure - state into a pure state only in the limit of infinite @xmath148 . otherwise , it is not possible . it is definitely not possible if we take into account experimental considerations . the story is different for the little groups . let us start with the rotation matrix of eq.([rot44 ] ) , and apply to this matrix the transformation matrix of eq.([sqimp ] ) . then @xmath153 & { } & \times \pmatrix{\cosh\eta & -\sinh\eta & 0 & 0 \cr    -\sinh\eta & \cosh\eta & 0 & 0 \cr     0 & 0 & 1 & 0 \cr 0 & 0 & 0 & 1 } .\end{aligned}\ ] ] if @xmath148 is zero , the above expression becomes the rotation matrix of eq.([rot44 ] ) . if @xmath148 becomes infinite , it becomes the little - group matrix @xmath154 of eq.([d1d2 ] ) applicable to the pure state of eq.([pure4 ] ) . the details of this calculation for the case of lorentz transformations are given in the 1986 paper by han _ et al . _ we are then led to the question of whether one little - group transformation matrix can be transformed from the other . if we carry out the matrix algebra of eq.([3mats ] ) , the result is @xmath155 where @xmath156 if @xmath157 , the above expression becomes the rotation matrix of eq.([rot44 ] ) . if @xmath158 , it becomes the @xmath135 matrix of eq.([d1d2 ] ) . here we used the parameter @xmath159 instead of @xmath148 . in terms of this parameter , it is possible to make an analytic continuation from the pure state with @xmath158 to an impure state with @xmath160 including @xmath157 .    on the other hand , we should keep in mind that the determinant of the density matrix is zero for the pure state , while it is non - zero for all impure states . for @xmath158 , the determinant vanishes , but it is nonzero and stays the same for all non - zero values of @xmath159 less than one and greater than or equal to zero . the analytic expression of eq.([anal ] ) hides this singular nature of the little group  @xcite . in secs .  [ polari ] , [ interfer ] , and [ little ] , it was noted that the stokes parameters form a four - vector in the minkowskian space . thus , it was possible to discuss the density matrix in terms of the four - vectors . the poincar sphere was originally constructed from polarization optics . therefore , it is also a representation of the lorentz group . the poincar sphere for various polarization states have been thoroughly discussed in a recent book by brosseau  @xcite . what is interesting is that the poincar sphere has two radii . one of them is the maximum radius specified by @xmath100 , and the other radius is the length of the space - like components , namely @xmath161 according to the four - vector property of the stokes parameter , @xmath162 is invariant under lorentz transformations . on the other hand , @xmath99 for a fully coherent state , while it is non - zero for a partially coherent state . the system is totally incoherent if @xmath163 . however , the lorentz group can not handle this decoherence mechanism . we observe here that the pure state with @xmath99 corresponds to the @xmath77-like little group for massless particles while it corresponds to the @xmath63-like little group for non - zero values of @xmath98 . the transition of @xmath63 to @xmath77 is well known as the group contraction . however , the inverse transformation requires further study . in analyzing optical rays in para - axial lens optics , we start with the lens matrix : @xmath164 and the translation matrix @xmath165 assuming that the beam is propagating along the @xmath11 direction .    then the one - lens system consists of @xmath166 if we perform the matrix multiplication , @xmath167 if we assert that the upper - right element be zero , then @xmath168 and the image is focussed , where @xmath169 and @xmath170 are the distance between the lens and object and between the lens and image respectively . they are in general different , but we shall assume for simplicity that they are the same : @xmath171 . we are doing this because this simplicity does not destroy the main point of our discussion , and because the case with two different values has been dealt with in the literature  @xcite . under this assumption , we are left with @xmath172 the diagonal elements of this matrix are dimensionless . in order to make the off - diagonal elements dimensionless , we write this matrix as @xmath173 indeed , the matrix in the middle contains dimensionless elements . the negative sign in front is purely for convenience . we are then led to study the core matrix @xmath174    here , the important point is that the above matrices can be written in terms of transformations of the lorentz group . in the two - by - two matrix representation , the lorentz boost along the @xmath11 direction takes the form @xmath175 and the rotation around the @xmath17 axis can be written as @xmath176 and the boost along the @xmath16 axis takes the form @xmath177 then the core matrix of eq.([core ] ) can be written as @xmath178 or @xmath179 if @xmath180 . if @xmath16 is greater than 2 , the upper - right element of the core is positive and it can take the form @xmath181 or @xmath182    the expressions of eq.([phi11 ] ) and eq.([chi11 ] ) are a lorentz boosted rotation and a lorentz - boosted boost matrix along the @xmath16 direction respectively . these expressions play the key role in understanding wigner s little groups for relativistic particles . let us look at their explicit matrix representations given in eq.([phi22 ] ) and eq.([chi22 ] ) . the transition from eq.([phi22 ] ) to eq.([chi22 ] ) requires the upper right element going through zero . this can only be achieved through @xmath148 going to infinity . if we like to keep the lower - left element finite during this process , the angle @xmath110 and the boost parameter @xmath183 have to approach zero . the process of approaching the vanishing upper - right element is necessarily a singular transformation . this aspect plays the key role in unifying the internal space - time symmetries of massive and massless particles . this is like einstein s @xmath184 becoming @xmath81 in the limit of large momentum .    on the other hand , the core matrix of eq.([core ] ) is an analytic function of the variable @xmath16 . thus , the lens matrix allows a parametrization which allows the transition from massive particle to massless particle analytically . the lens optics indeed serves as the analogue computer for this important transition in particle physics . > from the mathematical point of view , eq.([phi22 ] ) and eq.([chi22 ] ) represent circular and hyperbolic geometries , respectively . the transition from one to the other is not a trivial mathematical procedure . it requires a further investigation . let us go back to the core matrix of eq.([core ] ) . the @xmath16 parameter does not appear to be a parameter of lorentz transformations . however , the matrix can be written in terms of another set of lorentz transformations . this aspect has been discussed in the literature  @xcite . let us consider a co - axial system of an arbitrary number of lens . their focal lengths are not necessarily the same , nor are their separations . we are then led to consider an arbitrary number of the lens matrix given in eq.([lens ] ) and an arbitrary number of translation matrix of eq.([trans ] ) . they are multiplied like @xmath185 where @xmath80 is the number of lenses .    the easiest way to tackle this problem in to use the lie - algebra approach . let us start with the generators of the sp(2 ) group : @xmath186 since the generators are pure imaginary , the transformation matrices are real .    on the other hand , the @xmath39 and @xmath54 matrices of eq.([lens ] ) and eq.([trans ] ) are generated by @xmath187 if we introduce the third matrix @xmath188 all three matrices form a closed set of commutation relations : @xmath189 = ix_{3 } , \qquad \left[x_{1 } , x_{3}\right ] = -ix_{1 } ,   \qquad \left[x_{2 } , x_{3}\right ] = ix_{2 } .\ ] ] thus , these generators also form a closed set of lie algebra generating real two - by - two matrices . what group would this generate ? the answer has to be @xmath8 . the truth is that the three generators given in eq.([shear ] ) can be written as linear combinations of the generators of the @xmath8 group given in eq.([sq11 ] )  @xcite . thus , the @xmath190 matrices given above can also act as the generators of the @xmath8 group , and the lens - system matrix given in eq.([lsystem ] ) is a three - parameter matrix of the form of eq.([alpha2 ] ) with real elements . the resulting real matrix is written as @xmath191 and is called the @xmath192 matrix . the question then is how many lenses are need to give the most general form of the @xmath192 matrix  @xcite . this matrix has three independent parameter .    according to bargmann  @xcite , this three - parameter matrix can be decomposed into @xmath193 which can be written as the product of one symmetric matrix resulting from @xmath194 and one rotation matrix : @xmath195   & -\sin[(\beta + \alpha)/2 ]    \cr \sin[(\beta + \alpha)/2 ]    & \cos[(\beta + \alpha)/2 ] } .\ ] ] the expression in eq.([symm ] ) is a symmetric matrix , while that of eq.([ortho ] ) is an orthogonal one . we can then decompose each of these two matrices into the lens and translation matrices . the net result is that we do not need more than three lenses to describe the lens system consisting of an arbitrary number of lenses . the detailed calculations are given in ref . in a laser cavity , the optical ray makes round trips between two mirrors . one cycle is therefore equivalent to a two - lens system with two identical lenses and the same distance between the lenses . let us rewrite the matrix corresponding to the one - lens system given in eq.([core ] ) . @xmath196 then one complete cycle consists of @xmath197 .    here , the cycle starts from the midpoint between the two mirrors  @xcite , unlike the traditional approach to this problem where the cycle starts from the surface of one of the two mirrors  @xcite . by choosing the midpoint , we can eliminate the auxiliary flat mirror between them needed in the traditional approach . this problem was discussed in detail in ref . @xcite .    for @xmath80 cycles , the expression should be @xmath198 . this calculation , using the above expression , will not lead to a manageable form . however , we can resort to the expressions of eq.([phi11 ] ) and eq.([chi11 ] ) . then one cycle consists of @xmath199 [ b(\eta ) r(\phi/2 ) b(-\eta ) ]     \nonumber \\[2ex ] & { } & \hspace{5 mm } = b(\eta ) r(\phi ) b(-\eta ) , \end{aligned}\ ] ] if the upper - right element is negative . if it is positive , the expression should be @xmath200 [ b(\eta ) x(\chi/2 ) b(-\eta ) ]     \nonumber \\[2ex ] & { } & \hspace{5 mm } = b(\eta ) x(\chi ) b(-\eta ) .\end{aligned}\ ] ]    if these expressions are repeated @xmath80 times , @xmath201 if the upper - right element is negative . if it is positive , the expression should be @xmath202 as @xmath80 becomes large , @xmath203 and @xmath204 become very large , the beam deviates from the laser cavity . thus , we have to restrict ourselves to the case given in eq.([stable ] ) . the core of the expression of eq.([stable ] ) is the rotation matrix @xmath205^{n }      .\ ] ] this means that one complete cycle in the cavity corresponds to the rotation matrix @xmath206 . the rotation continues as the beam continues to repeat the cycle . let us go back to eq.([cycle ] ) . the expression corresponds to a lorentz boosted rotation , or bringing a moving particle to its rest frame , rotate , and boost back to the original momentum . the rotation associated with the momentum - preserving transformation is called the wigner s little - group rotation , which is related to the wigner rotation commonly mentioned in the literature  @xcite . the most efficient way to study multi - layer optics is to use the s - matrix formalism  @xcite . this formalism is also based on two - by - two matrices , and we can write @xmath207 where @xmath208 and @xmath209 are the incoming , reflected and transmitted beams . here we use the matrix notion @xmath192 for the two - by - two @xmath210 matrix . we consider a system of two different optical layers . for convenience , we start from the boundary from @xmath211 to @xmath212 . we can write the boundary matrix as  @xcite @xmath213 taking into account both the transmission and reflection of the beam . the parameter @xmath148 is of course determined by the transmission and reflection coefficient at the surface . this problem was studied in depth by monzon _ et al . _ @xcite , and their research on this subject is still continuing . as the beam goes through the @xmath212 , the beam undergoes the phase shift represented by the matrix @xmath214 when the wave hits the surface of the second medium , the boundary matrix is @xmath215 which is the inverse of the matrix given in eq.([bd11 ] ) . within the second medium , we write the phase - shift matrix as @xmath216 then , when the wave hits the first medium from the second , we have to go back to eq.([bd11 ] ) . thus , the one cycle consists of @xmath217 & { } & \hspace{10 mm } \times \pmatrix{\cosh(\eta/2 ) & -\sinh(\eta/2 ) \cr -\sinh(\eta/2 ) & \cosh(\eta/2 ) } \pmatrix{e^{-i\phi_2/2 } & 0 \cr 0 & e^{i\phi_2/2 } } .\end{aligned}\ ] ]    the above matrices contain complex numbers . however , it is possible to transform simultaneously @xmath218 to @xmath219 and transform @xmath220 to @xmath221 using a conjugate transformation , as is shown in appendix  [ conju ] . it is also possible to transform these expressions back to their original forms . this transformation property has been discussed in detail in ref . @xcite , and also in appendix  [ conju ] of the present paper .    as a consequence , the matrix of eq.([m1 ] ) becomes @xmath222 & { } &   \hspace{10mm}\times \pmatrix { e^{-\eta/2 } & 0 \cr 0 & e^{\eta/2 } } \pmatrix{\cos(\phi_2/2 ) & -\sin(\phi_2/2 ) \cr \sin(\phi_2/2 ) &   \cos(\phi_2/2 ) }   .\end{aligned}\ ] ] in the above expression , the first three matrices are of the same mathematical form as that of the core matrix for the one - lens system given in eq.([phi11 ] ) . the fourth matrix is an additional rotation matrix . this makes the mathematics of repetition more complicated , but this has been done  @xcite .    as a consequence the net result becomes @xmath223 or @xmath224 where the parameters @xmath225 and @xmath226 are to be determined from the input parameters @xmath227 and @xmath228 . detailed calculations are given in ref . @xcite . it is interesting to note that the lorentz group can serve as a computational device also in multi - layer optics . we have seen in this report that the lorentz group provides convenient calculational tools in many branches of ray optics . the reason is that ray optics is largely based on two - by - two matrices . these matrices also constitute the group @xmath26 which serves as the universal covering group of the lorentz group . the optical instruments discussed in this report are the fundamental components in optical circuits . in the world of electronics , electric circuits form the fabric of the system . in the future high - technology world , optical components will hold the key to technological advances . indeed , the lorentz group is the fundamental language for the new world . it is by now well known that the lorentz group is the basic language for quantum optics . coherent and squeezed states are representations of the lorentz group . it is challenging to see how the lorentz nature of the above - mentioned optical components will manifest itself in quantum world . the lorentz group was introduced to physics by einstein and wigner to understand the space - time symmetries of relativistic particles and the covariant world of electromagnetic fields . it is gratifying to note that the lorentz group can serve as the language common both to particle physics and optical sciences . the purpose of this appendix is to expand the third row of table  [ einwig ] using the covariant formalism of harmonic oscillators . we explain first what this formalism does for covariant picture of relativistic extended particles . we then point out that this oscillator formalism can serve as the basic language for two - mode coherent states  @xcite , namely squeezed states of light . the concept of localized probability distribution is the backbone of the present formulation of quantum mechanics . of course , we would like to have more deterministic form of dynamics , and efforts have been and are still being made along this direction . one of the most serious problems with this probabilistic interpretation is whether this concept of probability is consistent with special relativity .    in a given lorentz frame , we know how to do quantum mechanics with a localized probability distribution . how would this distribution look to an observer in a different lorentz frame ? * would the probability distribution appear the same to this observer ? * if different , how is the probability distribution distorted ? * is the total probability conserved ? * the lorentz boost mixes up the spatial coordinate with the time variable . what role does the time - separation variable play in defining the boundary condition for localization and the probability distribution ? we can answer some or all of the above questions only if we construct covariant wave functions , namely wave functions which can be lorentz - boosted . it is easy to construct these wave functions if we know the answer . in the initial development of quantum mechanics , the harmonic oscillator played the pivotal role . thus , it is clear to us that if there is a wave function which can be lorentz - boosted , this has to be the harmonic oscillator wave function . until we construct wave functions which can be lorentz - transformed , we can not say that quantum mechanics is consistent with relativity . indeed , we should examine this problem before attempting to invent more definitive quantum mechanics . since the hadron , in the quark model , is a bound - state of quarks  @xcite , feynman , kislinger , and ravndal , in 1971  @xcite , raised the following question in connection with the quark model for hadrons : the hadronic spectrum can indeed be explained in terms of the degeneracy of three - dimensional harmonic oscillator wave functions in the hadronic rest frame ; however , what happens when the hadron moves ? indeed , et al . _ wrote a lorentz - invariant differential equation whose solutions can become non - relativistic wave functions if the time - separation variable can be ignored . this lorentz - invariant differential equation is a four - dimensional partial differential equation , with many different solutions depending on separation of variables and boundary conditions . there is a set of normalizable solutions which can serve as a representation space for wigner s little group for massive particles  @xcite . we can start with this set of solutions and give physical interpretations , especially to the time - separation variable .    finally , is the covariance of the oscillator wave function consistent with what we observe in the real world . here again , feynman plays the key role . while the quark model can be fit into the oscillator scheme in the hadronic rest frame , feynman in 1969 came up with the idea of partons  @xcite . according to feynman s parton model , the hadron consists of an infinite number of partons when it moves with a velocity close to that of light . quarks and partons are believed to be the same particles , but their properties are totally different . while the quarks inside the hadron interact coherently with external signals , partons interact incoherently . if the partons are lorentz - boosted quarks , does the lorentz boost destroy the coherence ?    before 1964  @xcite , the hydrogen atom was used for illustrating bound states . these days , we use hadrons which are bound states of quarks . let us use the simplest hadron consisting of two quarks bound together with an attractive force , and consider their space - time positions @xmath229 and @xmath230 , and use the variables @xmath231 the four - vector @xmath232 specifies where the hadron is located in space and time , while the variable @xmath16 measures the space - time separation between the quarks . according to einstein , this space - time separation contains a time - like component which actively participates as can be seen from @xmath233 when the hadron is boosted along the @xmath11 direction . in terms of the light - cone variables defined as  @xcite @xmath234 the boost transformation of eq.([boostm ] ) takes the form @xmath235 the @xmath46 variable becomes expanded while the @xmath48 variable becomes contracted . does this time - separation variable exist when the hadron is at rest ? yes , according to einstein . in the present form of quantum mechanics , we pretend not to know anything about this variable . indeed , this variable belongs to feynman s rest of the universe . in this report , we shall see the role of this time - separation variable in the decoherence mechanism . also in the present form of quantum mechanics , there is an uncertainty relation between the time and energy variables . however , there are no known time - like excitations . unlike heisenberg s uncertainty relation applicable to position and momentum , the time and energy separation variables are c - numbers , and we are not allowed to write down the commutation relation between them . how does this space - time asymmetry fit into the world of covariance  @xcite . the answer is that wigner s @xmath63-like little group is not a lorentz - invariant symmetry , but is a covariant symmetry  @xcite . it has been shown that the time - energy uncertainty relation applicable to the time - separation variable fits perfectly into the @xmath63-like symmetry of massive relativistic particles  @xcite . the c - number time - energy uncertainty relation allows us to write down a time distribution function without excitations  @xcite . if we use gaussian forms for both space and time distributions , we can start with the expression @xmath236 for the ground - state wave function . what do feynman _ et al . _ say about this oscillator wave function ?    in their classic 1971 paper @xcite , feynman _ _ start with the following lorentz - invariant differential equation . @xmath237 this partial differential equation has many different solutions depending on the choice of separable variables and boundary conditions . _ insist on lorentz - invariant solutions which are not normalizable . on the other hand , if we insist on normalization , the ground - state wave function takes the form of eq.([ground ] ) . it is then possible to construct a representation of the poincar group from the solutions of the above differential equation  @xcite . if the system is boosted , the wave function becomes @xmath238 this wave function becomes eq.([ground ] ) if @xmath148 becomes zero . the transition from eq.([ground ] ) to eq.([eta ] ) is a squeeze transformation . the wave function of eq.([ground ] ) is distributed within a circular region in the @xmath239 plane , and thus in the @xmath240 plane . on the other hand , the wave function of eq.([eta ] ) is distributed in an elliptic region with the light - cone axes as the major and minor axes respectively . if @xmath148 becomes very large , the wave function becomes concentrated along one of the light - cone axes . indeed , the form given in eq.([eta ] ) is a lorentz - squeezed wave function . this lorentz - squeeze mechanism is illustrated in fig . [ ellipse ] . there are many different solutions of the lorentz invariant differential equation of eq.([osceq ] ) . the solution given in eq.([eta ] ) is not lorentz invariant but is covariant . it is normalizable in the @xmath241 variable , as well as in the space - separation variable @xmath11 . how can we extract probability interpretation from this covariant wave function ? this issue has been discussed thoroughly in the literature  @xcite . another pressing problem in physics is that hadrons , like the proton , can be regarded as quantum bound states of quarks when they move slowly . however , they appear like collections of partons when they move with velocity very close to that of light  @xcite . since the quarks and partons have quite different properties , it is a challenging problem in modern physics to show that they are the same but appear different depending on the observer s lorentz frame . it has been shown that this lorentz - squeezed wave function can explain the quark model and the parton model as two different manifestations of one covariant entity  @xcite . thus , the covariant harmonic oscillator can occupy the third row of table  [ einwig ] . the lorentz - squeezed wave function can be written as @xmath242\right\ } .\ ] ]    as was discussed in the literature for several different purposes  @xcite , this wave function can be expanded as @xmath243 this is an expansion in terms of the oscillations along two different directions . if those directions are @xmath244 and @xmath245 , we can write eq.([expan ] ) as @xmath246 the mathematics of harmonic oscillators can be translated into second - quantized fock space . in the fock space , this expression describes the two - photon coherent state discussed first by yuen in 1976  @xcite . the group @xmath26 consists of two - by - two unimodular matrices whose elements are complex . there are therefore six independent parameters , and thus six generators of the lie algebra . this group is locally isomorphic to the six - parameter lorentz group or @xmath247 applicable to the minkowskian space of three space - like directions and one time - like direction . we can construct representations of the group starting from the lie algebra given in secs . [ lorentz ] and [ spinor ] . we can construct the representation by using a method similar to what goldstein did for the three - dimensional rotation group in terms of the euler angles  @xcite . there are three - generators for the rotation group , but goldstein starts with rotations around the @xmath11 and @xmath16 directions . rotations around the @xmath17 axis and the most general form for the rotation matrix can be constructed from repeated applications of those two starting matrices . let us call this type of approach the `` euler construction . ''    in constructing the lorentz group , we observe first that we need rotations around two different directions . as for the boost , we need boosts along one given direction since boots along other directions can be achieved by rotations . there are three basic advantages of this approach . first , the number of `` starter '' matrices is less than the number of generators . for example , we need only two starters for the three - parameter rotation group . in our case , we started with two matrices for the three - parameter group @xmath8 and also for @xmath9 . second , each starter matrix takes a simple form and has its own physical interpretation . the third advantage can be stated in the following way . repeated applications of the starter matrices will lead to a very complicated expression . however , the complicated expression can be decomposed into the minimum number of starter matrices . for example , this number is three for the three - dimensional rotation group . this number is also three for @xmath91 and @xmath8 . we call this the euler decomposition . the present paper is based on both the euler construction and the euler decomposition .    among the several useful euler decompositions , the iwasawa decomposition plays an important role in the lorentz group . we have seen in this paper what the decomposition does to the two - by - two matrices of @xmath8 , but it has been an interesting subject since iwasawa s first publication on this subject  @xcite . it is beyond the scope of this paper to present a historical review of the subject . however , we would like to point out that there are areas of physics where this important mathematical theorem was totally overlooked .    for instance , in particle theory , wigner s little groups dictate the internal space - time symmetries of massive and massless particles which are locally isomorphic to @xmath63 and @xmath77 respectively  @xcite . the little group is the maximal subgroup of the lorentz group whose transformations do not change the four - momentum of a given particle  @xcite . the @xmath77-like subgroup for massless particles is locally isomorphic to the subgroup of @xmath26 which can be started from one of the matrices in eq.([shear1 ] ) and the diagonal matrix of eq.([ps11 ] ) . thus there was an underlying iwasawa decomposition while the @xmath77-like subgroup was decomposed into rotation and boost matrices  @xcite , but the authors did not know this .    in optics , there are two - by - two matrices with one vanishing off - diagonal element . it was generally known that this has something to do with the iwasawa effect , but simon and mukunda  @xcite and han _ et al . _ @xcite started treating the iwasawa decomposition as the main issue in their papers on polarized light .    in para - axial lens optics , the matrices of the form given in eq.([shear1 ] ) are the starters , and repeated applications of those two starters will lead to the most general form of @xmath8 matrices . it had been a challenging problem since 1985  @xcite to write the most general two - by - two matrix in lens optics in terms the minimum number of those starter matrices . this problem has been solved recently  @xcite , and the central issue in the problem was the iwasawa decomposition . the core matrix of eq.([core ] ) contains the chain of the matrices @xmath248 the lorentz group allows us to simplify this expression under certain conditions .    for this purpose , we transform the above expression into a more convenient form , by taking the conjugate of each of the matrices with @xmath249 then @xmath250 leads to @xmath251 in this way , we have converted @xmath252 of eq.([su11 ] ) into a real matrix , but it is not simple enough .      after multiplication , the matrix of eq.([abcd ] ) will take the form @xmath259 where @xmath260 and @xmath261 are real numbers . if @xmath82 and @xmath131 vanish , this matrix will become diagonal , and the problem will become too simple . if , on the other hand , only one of these two elements become zero , we will achieve a substantial mathematical simplification and will be encouraged to look for physical circumstances which will lead to this simplification . let us summarize . we started in this section with the matrix representation @xmath252 given in eq.([su11 ] ) . this form can be transformed into the @xmath56 matrix of eq.([abcd ] ) through the conjugate transformation @xmath262 where @xmath131 is given in eq.([ccc ] ) . conversely , we can recover the @xmath252 representation by @xmath263 for calculational purposes , the @xmath56 representation is much easier because we are dealing with real numbers . on the other hand , the @xmath252 representation is of the form for the s - matrix we intend to compute . it is gratifying to see that they are equivalent . let us go back to eq.([abcd ] ) and consider the case where the angles @xmath110 and @xmath226 satisfy the following constraints . @xmath264 thus @xmath265 then in terms of @xmath108 , we can reduce the matrix of eq.([abcd ] ) to the form @xmath266 thus the matrix takes a surprisingly simple form if the parameters @xmath108 and @xmath148 satisfy the constraint @xmath267then the matrix becomes @xmath268 this aspect of the lorentz group is known as the iwasawa decomposition  @xcite , and has been discussed in the optics literature  @xcite . matrices of this form are not so strange in optics . in para - axial lens optics , the translation and lens matrices are written as @xmath269 respectively . these matrices have the following interesting mathematical property  @xcite , @xmath270 and @xmath271 we note that the multiplication is commutative , and the parameter becomes additive . these matrices convert multiplication into addition , as logarithmic functions do . s. ferrara and c. savoy , in _ supergravity 1981 _ , s. ferrara and j. g. taylor eds . ( cambridge univ . press , cambridge , 1982 ) , p. 151 . see also p. kwon and m. villasante , j. math . * 29 * , 560 ( 1988 ) ; _ ibid . _ * 30 * , 201 ( 1989 ) . b. c. sanders and a. mann , group 22 , proceedings of the 22nd international colloquium on group theoretical methods in physics , s. p. cornel _ et al_. eds . ( international press , boston , 1999 ) . see pp 474 - 478 . j. j. monzn and l. l. snchez - soto , am . * 64 * , 156 ( 1996 ) , j. j. monzn and l. l. snchez - soto , j. opt . a * 17 * , 1475 ( 2000 ) ; j. j. monzn , t. yonte , l. l. snchez - soto , and j. carinena , j. opt . am . a , * 19 * , 985 ( 2002 ) . r. p. feynman , _ the behavior of hadron collisions at extreme energies _ , in _ high energy collisions _ , proceedings of the third international conference , stony brook , new york , edited by c. n. yang _ et al . _ , pages 237 - 249 ( gordon and breach , new york , 1969 ) .
Abstract:

Answer: it has been almost one hundred years since einstein formulated his special theory of relativity in 1905 . he showed that the basic space - time symmetry is dictated by the lorentz group . it is shown that this group of lorentz transformations is not only applicable to special relativity , but also constitutes the scientific language for optical sciences . it is noted that coherent and squeezed states of light are representations of the lorentz group . the lorentz group is also the basic underlying language for classical ray optics , including polarization optics , interferometers , the poincare sphere , one - lens optics , multi - lens optics , laser cavities , as well multilayer optics . * lorentz group in ray optics *    sibel bakal + department of physics , middle east technical university , + 06531 ankara , turkey    elena georgieva + science systems and applications , inc . , lanham , md 20771 , u.s.a . , and + national aeronautics and space administration , + goddard space flight center , laser and electro - optics branch , + code 554 , greenbelt , maryland 20771 , u.s.a .    y. s. kim + department of physics , university of maryland , + college park , maryland 20742 , u.s.a . +    marilyn e. noz + department of radiology , new york university , new york , new york 10016 , u.s.a .