arxiv:physics/ v3 [physics.optics] 6 Dec 2007

Transcription

1 Fermat s priniple of least time in the presene of uniformly moving boundaries and media Aleksandar Gjurhinovski and Aleksandar Skeparovski Department of Physis, Faulty of Natural Sienes and Mathematis, Sts. Cyril and Methodius University, P. O. Box 16, 1000 Skopje, Maedonia Dated: Deember 06, 007) Eur. J. Phys ) arxiv:physis/070100v3 physis.optis] 6 De 007 The refration of a light ray by a homogeneous, isotropi and non-dispersive transparent material half-spae in uniform retilinear motion is investigated theoretially. The approah is an amalgamation of the original Fermat s priniple and the fat that an isotropi optial medium at rest beomes optially anisotropi in a frame where the medium is moving at a onstant veloity. Two ases of motion are onsidered: a) the material half-spae is moving parallel to the interfae; b) the material half-spae is moving perpendiular to the interfae. In eah ase, a detailed analysis of the obtained refration formula is provided, and in the latter ase, an intriguing bakward refration of light is notied and thoroughly disussed. The results onfirm the validity of Fermat s priniple when the optial media and the boundaries between them are moving at relativisti speeds. PACS numbers: p, i, 0.30.Xx I. INTRODUCTION The variational approah in geometrial optis was established in 1657 when the famous Frenh mathematiian Pierre de Fermat provided a mathematial proof that the straight line is not the fastest way for light to traverse between two optial media 1,, 3]. This simple and new way of looking at the propagation of light was already empirially notied by Hero of Alexandria fifteen enturies before Fermat by observing that the path of the light between two points upon refletion from a plane mirror is the shortest possible path. Fermat s insightful and intriguing idea that the Nature always ats in the shortest ways has influened generations of sientists, like Maupertius, Euler, Lagrange, Hamilton and Feynman, to suessively translate lassial and quantum physis into the language of variational alulus 1, 4, 5, 6]. Today, the pratial value of Fermat s priniple is immense, ranging from appliations in geophysis, oeanography and material physis on one side, to general relativity and astrophysis on the other 7, 8, 9, 10, 11, 1, 13, 14, 15, 16, 17, 18, 19, 0]. The modern formulation of Fermat s priniple states that a light in going between two fixed points traverses an optial path length that is stationary with respet to the variations of that path. In most ases, this formulation redues to the original one that the atual path between the two points is the path whih requires the shortest time 1]. An aurate desription of the proess beyond Fermat s priniple is provided by the modern quantum eletrodynamis, where the light is onsidered to be made out of photons. Aordingly, the photon takes all the available paths in going between two fixed points, eah path at a ertain probability, and the priniple of least Eletroni address: agjurin@iunona.pmf.ukim.edu.mk Eletroni address: skepalek@iunona.pmf.ukim.edu.mk time follows as an approximation of this argument, 3, 4]. Fermat s priniple has been initially established for optial systems whose properties are time-independent, and, as suh, has entered the mainstream of the standard undergraduate physis urriulum through the derivations of the well-known laws of refletion and refration of light at a smooth stationary interfae between two optial materials at rest ]. Reently, Fermat s priniple was proven valid for optial systems with time-dependent parameters, inorporating the situations of light propagation in uniformly moving optial media 5, 6, 7]. It was also shown that Einstein s law of refletion of light from a plane mirror in uniform retilinear motion follows from Fermat s priniple and the onstant speed of light postulate, providing an indiret proof of the validity of Fermat s priniple when the boundaries upon whih the light is refleted are moving at relativisti speeds 8]. Apart from the latter example, it seems that the problems of light propagation involving uniformly moving boundaries and media have never been attaked diretly by using Fermat s priniple of least time. Motivated by the latest extensions to non-stationary ases, in this paper we will use Fermat s priniple to derive the formulas for the law of refration of a light ray in two speifi situations when the boundaries and the media are moving at onstant veloities. We will limit our analysis to optial media that are homogeneous, isotropi and nononduting in their rest frames of referene. Also, we onsider the media to be non-dispersive at rest, in the sense that the orreting terms due to dispersion are so small that an be negleted. We begin our disussion in Se. II by showing that an optial medium, isotropi in its rest frame, exhibits an optial anisotropy with respet to the frame where the medium is in uniform retilinear motion. We derive the expression for the speed of the photon in the moving medium and show that it depends on the angle between the light ray and the veloity of the medium. We de- Typeset by REVTEX

2 sribe this relativistially indued optial anisotropy in a more onvenient way, by introduing an effetive refrative index of the moving medium, defined as a ratio between the speed of the photon in vauum and the speed of the photon in the moving medium. The notion that the photon will see the moving medium as a stationary medium having an effetive refrative index different from its refrative index at rest, is inorporated into Fermat s priniple in Se. III and IV to investigate the refration of the photon inident from a vauum half-spae upon a uniformly moving material half-spae. In Se. III we investigate the situation when the material half-spae is moving parallel to the interfae, and in Se. IV we take the material half-spae to move perpendiularly to the interfae. In both ases, we onsider the plane of inidene to be normal to the vauum-material interfae and parallel to the veloity of the medium. A summary of the results and a omparison with other similar treatments are given in Se. V, where the possibility to extend the method to more ompliated onfigurations is also suggested. Alternative derivation of the refration law when the material half-spae is moving perpendiular to the interfae is presented in Appendix A. Finally, in Appendix B, the onditions for the appearane of the bakward refration are derived analytially. II. AN OPTICAL ANISOTROPY OF A TRANSPARENT MEDIUM INDUCED BY ITS UNIFORM MOTION If a flash of light is emitted from a fixed point soure in an empty spae a vauum), the outgoing photons will move radially on a spherial wavefront expanding at a speed of light 9]. With respet to an observer from a frame that moves in a straight line at a onstant veloity, the wavefront of the pulse will remain sphere, expanding at the same onstant speed in all diretions 30]. The result is a orollary of the seond postulate of speial relativity that the speed of light in a vauum is the same in all inertial frames of referene regardless of the motion of the light soure 31]. The situation is hanged if the expansion of the pulse is taking plae in a material medium 3, 33]. If we onsider a homogeneous, isotropi and non-dispersive transparent medium at rest with respet to S -frame, and if we treat the problem two-dimensionally, the equation that desribes the evolution of the pulse is: x + y = /n) t, 1) where n is the refrative index of the material at rest. This is an equation of an expanding irle whose radius at time t is /n)t. By using the Lorentz transformation: x = γx ut), y = y, t = γt ux/ ), ) with γ = 1 u / ) 1/, we an make a transition to S-frame where the medium is moving at a onstant speed FIG. 1: An illustration of the light drag in the S-frame where the medium is moving at a onstant speed u < /n) in the diretion of the positive x-axis. The motion of the medium is a ause for the wavefront of the pulse originating from O to be a dragged ellipse possessing an axial symmetry with respet to the x-axis. Here, the photon emitted from O at an angle ϕ will travel at a speed vϕ). u in the positive diretion of the x-axis. In this way, Eq. 1) is transformed into: x ηt) at) + y = 1, 3) bt) where by η, a and b we denoted: ) 1 1/n η = u 1 u /n, 4) 1 u / ) a = /n) 1 u /n, 5) 1 u / ) 1/ b = /n) 1 u /n. 6) Aording to Eq. 3), the wavefront of the pulse the ray surfae) with respet to an observer in S-frame is an expanding ellipse dragged by the moving medium see Fig. 1). At time t from the emission of the pulse, the produt ηt denotes the dragging parameter, that is, the distane from the origin O of the pulse the origin of the xy-oordinate system) to the enter of the ellipse, and the produts at and bt are the semi-minor and semimajor axes of the ellipse, respetively. Consequently, the parameters a and b are the speeds of expansion of the semi-minor and semi-major axes of the ellipse, respetively, and η is the speed of the enter of the ellipse along

3 3 the x-axis. When η < a, the ellipse will be inompletely dragged, and the origin of the pulse will remain inside the ellipse at all times. This is the ase of a subluminal motion of the medium, when the speed of the medium u is smaller than /n. On the other hand, when η a, the dragging is overwhelming, and the ellipse no longer enloses the origin. In this ase the medium is moving at superluminal speeds /n < u < ), ausing an existene of a Mah one, outside of whih no light signal an ever reah, and inside of whih light is doubled up, with two pulses of light in eah diretion 9, 3, 33]. We may translate Eq. 3) into the polar oordinates by making the substitutions x = vt os ϕ and y = vt sin ϕ: ) 1 u sin ϕ u n os ϕ v u osϕ 1 1n )v + u = 0. 7) n Here, v is the speed of the photon along a light ray direted at an angle ϕ from the veloity of the medium see Fig. 1). Equation 7) is a quadrati in v, and it an be solved for v in terms of the angle ϕ, the speed of the medium u, and the refrative index n of the medium at rest. By rejeting the seond solution of Eq. 7) from the requirement that v = /n when u = 0, we obtain: vϕ) = u1 1/n )osϕ + /n u sin ϕ + os ϕ/n ) ] 1/ 1 u / ] } 1/ 1 u / sin ϕ + os ϕ/n )] 1. 8) From Eq. 8) we onlude that the ray veloity in the moving medium depends on the angle ϕ between the diretion of the ray and the veloity of the medium. In this sense, an optial medium that is optially isotropi in its rest frame will possess an optial anisotropy in the frame in whih it is moving at a onstant veloity. This indued optial anisotropy is of a purely relativisti origin, and it is different from the usual anisotropy in the rystals 34]. The observer in the S-frame may trak the path of a light ray in the moving medium by onsidering it to be a stationary medium that possesses a relativistially indued optial anisotropy haraterized by an effetive refrative index n ef ϕ): n ef ϕ) = vϕ). 9) Realling Eq. 8) for vϕ), the expression for n ef in Eq. 9) an be re-written in terms of u, n and ϕ: n ef ϕ) = u/)1 n )osϕ + n1 u / ] 1/ u/) n 1)os ϕ + 1 n u / ] } 1/ 1 n u / ) 1. 10) FIG. : The least-time derivation of the refration formula when the material half-spae is moving at a onstant speed u in the diretion of the positive x-axis. When n = 1, Eqs. 8) and 10) redue to n ef = 1 and v = for all ϕ, whih desribes an expansion of the pulse in a vauum. Also, if u = 0, then n ef = n and v = /n for all ϕ, whih is the ase when the medium is stationary. If the photon moves parallel to the veloity of the medium, n ef n±u/1 n ) and v /n±u1 1/n ) to the first order in u/, where the sign is taken + or depending on whether ϕ = 0 or ϕ = 180, respetively. This longitudinal light drag has been theoretially predited by Fresnel 1818) on the basis of his ether theory, and experimentally onfirmed by Fizeau 1851) in his elebrated running-water experiment 3]. III. REFRACTION FROM A MATERIAL HALF-SPACE MOVING PARALLEL TO THE INTERFACE A photon traveling in vauum is inident upon the surfae of a semi-infinite optial material whih oupies the positive region of the y-axis see Fig. ). The material half-spae is moving at a onstant speed u in the positive diretion of the x-axis, parallel to the vauum-material interfae. The material in its rest frame has a refrative index n, and is assumed to be homogeneous, isotropi and non-dispersive. We will use Fermat s priniple to derive the formula for the law of refration of the photon from the interfae between the vauum and the moving medium. Let A be the spae-point in the vauum half-spae, belonging to the inident ray of the photon, and B is the spae-point in the moving medium that lies on the refrated ray. By d A and d B we denote the shortest dis-

4 4 tanes between the interfae and the points A and B, respetively, and l is the distane between the orthogonal projetions of the points A and B on the interfae. Aording to Fermat s priniple, the true path between the two fixed points A and B taken by the photon is the one that is traversed in the least time. Of all the points along the interfae, the photon will be refrated at the point O whose x-position minimizes the total time t AB required for the photon to over the path from A to B: Here, t AB = t AO + t OB. 11) t AO = x + d A) 1/ 1) is the time for the photon to travel from A to O at a speed of light, and l x) + d 1/ t OB = B] 13) vϕ) is the orresponding time from O to B. By vϕ) we denote the speed of the photon along the refrated ray at an angle ϕ from the veloity of the medium. If we take into aount that vϕ) = /n ef ϕ), and use Eqs. 1) and 13), we may re-write the total transit time t AB in Eq. 11) as: t AB = x + d A) 1/ + n efϕ) l x) + d ] 1/ B. 14) We substitute n ef from Eq. 10) into Eq. 14), use the relation osϕ = l x l x) + d B ]1/ 15) from Fig., and simplify the result to obtain: x + d 1/ t AB = A) + u/ )1 n ) 1 n u l x) / + n/) 1 u / ) 1/ 1 n u / u/) n 1)l x) + 1 n u / ) l x) + d ] } 1/ B. 16) The priniple of least time requires that dt AB /dx = 0, and thus: 0 = x x + d u/)1 n ) A )1/ 1 n u / n1 u / ) 3/ 1 n u / l x) u/) n 1)l x) l x) + d B ]1/ l x) + d B ] + 1 n u / } 1/. 17) From Fig. we notie the expressions for the angle of inidene α and the angle of refration β of the photon: sin α = sin β = whih we use into Eq. 17) to get: x, x + d A )1/ 18) l x l x) + d, B ]1/ 19) 0 = sin α u/)1 n ) 1 n u / n1 u / ) 3/ 1 n u / sinβ u/) n 1)sin β + 1 n u / ]. 0) 1/ If we employ the identity sinβ = tanβ1 + tan β) 1/ into Eq. 0) and do some algebra, we obtain a quadrati equation in tanβ: n 1 u ) sin α sin α u ) ] tan β 1 n u 1 1 u / ) sin α + u n 1)] = 0. 1) Equation 1) has two solutions in tanβ, and the solution: tan β = n 1)u/) + 1 n u / )sin α ] n 1 u/)sinα ] } 1/ sin α u/) 1 u / ) 1/ ) is the only one that is physially orret, whih we hoose from the requirement that when u = 0 it should redue to the usual Snell s law of refration. Equation ) is the formula for the law of refration of the photon from the vauum-material interfae when the optial material is moving uniformly parallel to its surfae. It gives the angle of refration β of the photon as a funtion of the inident angle α, the speed u of the medium, and the refrative index n of the medium in its rest frame. Two alternative derivations leading to the same refration formula are given in Ref. 3].] The funtional dependene of the angle of refration β on the inident angle α is plotted in Fig. 3 for n = 1.5 and for several different values of the speed u of the medium. We immediately notie that the law of refration in Eq. ) differs onsiderably from the usual Snell s law in the ase of a stationary medium. First of all, when u/ > 0, the angle of refration β attain non-zero values even when the inident angle α is zero. In other words, the photon that enters the moving medium perpendiularly to its surfae will be dragged by the medium in the moving diretion. This is the transverse Fresnel-Fizeau light drag, theoretially predited by Fresnel in 1818 and experimentally deteted by Jones in the early 1970s with an auray up to the first order of u/ 3, 35, 36, 37].

5 5 FIG. 3: The angle of refration β versus the inident angle α for n = 1.5 and for different speeds u of the material halfspae, when the material half-spae is moving parallel to the interfae. The influene of the dragging effet on the refration is also visible for non-zero inident angles, espeially when they take negative values 38]. In the latter ase, if u/ > 0, starting from α = 0 down to a ertain negative inident angle, the angle of refration β is positive see Fig. 3), whih means that both the inident and the refrated ray of the photon will lie on the same side with respet to the interfae normal. Moreover, when the speed u of the medium beomes greater than /n, the angle of refration will be positive for any value of the inident angle α. Equation ) an also be used to alulate the angle of refration of the photon when the medium moves in the opposite diretion to the one in Fig., that is, in the negative diretion of the x-axis. However, the situation when the photon is inident at a positive angle α upon the medium whih moves at a speed u is physially idential to the situation whih orresponds to a negative inident angle α and a speed u. At the end, we emphasize that the absolute value of the angle of refration β in Eq. ) annot beome equal or greater than 90, no matter what the values of α, n and u are. In other words, it is not possible to observe the situation in whih the photon would slide along the surfae of the medium, or to be refleted bak to the vauum half-spae. To prove this property, we will try to find suh a value of α whih will turn the denominator in the right-hand side of Eq. ) into zero. This is the ondition for the angle β to equal ±90. Hene, α should be suh that: n 1 u ) sin α sin α u ) = 0, 3) whih is a quadrati equation in sin α, and has two solutions: sin α) 1, = n ± u/ nu/ ± 1. 4) Having in mind that n > 1 and u <, we have n ± FIG. 4: The least-time derivation of the refration formula when the material half-spae is moving at a onstant speed u in the positive diretion of the y-axis. u/)nu/ ± 1) 1 > 1, whih means that there are no real values for the angle α to satisfy Eq. 3). IV. REFRACTION FROM A MATERIAL HALF-SPACE MOVING PERPENDICULAR TO THE INTERFACE In the following, we onsider the ase in whih the material half-spae in Se. III is moving at a onstant speed u in the positive diretion of the y-axis. The geometry of the problem is shown in Fig. 4. Observe that unlike the ase in Se. III, the vauum-material interfae also is in uniform motion, having the same onstant speed u as the medium in the positive y-diretion. A photon is emitted from point A in the vauum halfspae. At this instant, the vauum-material interfae oinides with the x-axis, and d A is the shortest distane between the interfae and the point A. After being refrated from the moving interfae, the photon reahes the fixed point B in the material half-spae. We denote by d B the shortest distane between the interfae and the point B at time when the photon was emitted from A, and by l the distane between the orthogonal projetions of the points A and B on the interfae. Due to the motion of the vauum-material interfae, the hypothetial points from whih the photon may be refrated at the interfae will not lie on the same horizontal. Aording to Fermat s priniple, the photon will be refrated from a point O that makes AOB the path of

6 6 the shortest transit time see Fig. 4). We denote by t AO the time required for the photon to go from A to O, and by t OB the time from O to B. Observe that when the photon reahes point O, the interfae has already rossed the distane ut AO in the positive y-diretion from its position y = 0 when the photon was emitted from A. The total time t AB for the photon to travel the path AOB is: From Fig. 4 we have: t AB = t AO + t OB. 5) t AO ) = d A + ut AO ) + x, 6) vβ)t OB ] = d B ut AO ) + l x), 7) where vβ) is the speed of the photon in the moving medium in the diretion determined by the angle β. The speed vβ) of the photon follows from Eq. 8) if we simply substitute β for ϕ. In order to simplify the derivation of the refration formula, we re-write Eqs. 6) and 7) as: t AO = 1 da + ut AO ) + x ] 1/, 8) t OB = n efβ) db ut AO ) + l x) ] 1/, 9) with n ef β) being the effetive refrative index of the moving half-spae in this ase: n ef β) = u/)1 n )osβ + n1 u / ] 1/ u/) n 1)os β + 1 n u / ] } 1/ 1 n u / ) 1. 30) Thus, the total transit time t AB in Eq. 5) an be rewritten as: t AB = 1 da + ut AO ) + x ] 1/ + n efβ) db ut AO ) + l x) ] 1/. 31) If we substitute Eq. 30) into Eq. 31), and take into aount the expression: d B ut AO osβ =, 3) d B ut AO ) + l x) 1/ ] from Fig. 4, we obtain an expression for the transit time t AB of the photon as a funtion of the position x of the point of refration O at the moving interfae: t AB = da + ut AO ) + x ] 1/ + n/) 1 u / ) 1/ 1 n u / + u + u1 n )d B ut AO ) n u n 1)d B ut AO ) 1 n u ) 1/ db ut AO) + l x) ]}. 33) We now invoke the priniple of least time that the value of x should be suh that the transit time t AB of the photon is minimum. We take the derivative of t AB in Eq. 33) with respet to x and set the result to zero to obtain: ) 0 = u osα + sin α + u /)n 1) dx 1 n u / ) n 1 u / ) 1/ ) dx 1 n u / u 1 u ) osβ + 1 n u ) ] dx sin β ] 1/ u n 1)os β + 1 n u, 34) where we employed the relations: sin α = osα = sin β = x x + d A + ut AO ) ] 1/, 35) d A + ut AO x + d A + ut AO ) ] 1/, 36) l x l x) + d B ut AO ) ] 1/, 37) from Fig. 4, together with Eq. 3) and the fat that t AO = t AO x). To find dt AO /dx, we take the derivative of Eq. 6) with respet to x and obtain: t AO dx ) = d A + ut AO )u dx ) + x, 38) from whih we get: ) x/t AO ) = dx 1 u/)d A + ut AO )/t AO )]. 39) By reognizing from Fig. 4 that sinα = x/t AO ) and osα = d A + ut AO )/t AO ), we reast Eq. 39) into the form: dx ) = sin α u osα. 40) We substitute Eq. 40) into Eq. 34), apply the trigonometri identities sinβ = tanβ1+tan β) 1/ and osβ = 1 + tan β) 1/, and do some algebra to obtain a quadrati equation in tanβ: n 1 n u / ) 1 u / sin α 1 u/)osα] } tan β + n u/)sinα 1 u/)osα tan β 1 u / )sin α = 0. 41) 1 u/)osα] Equation 41) has two solutions in tanβ, and its only physially orret solution is tanβ = 1 u / )sin α n u/)1 u/)osα]

7 7 FIG. 5: The angle of refration β versus the inident angle α for n = 1.5 and for different speeds u of the material half-spae moving in the positive diretion of the y-axis. The portions of the urves orresponding to α arosu/) are drawn with dashed lines. At inident angles orresponding to the dashed urves, the refration does not our. } 1/ 1 + n 1 u/)osα] 1 u / )sin α], 4) whih is the formula for the law of refration of the photon in this ase. We rejet the seond solution of Eq. 41) as a physially inorret one, by using the same argument as in the dedution of Eq. ) from Eq. 1). A derivation of Eq. 4) via the relativisti veloity transformation formulas is provided in Appendix A.] The obtained refration formula works both for positive and negative values of the inident angle α, as well as for both moving diretions of the medium. The righthand side of Eq. 4) is an odd funtion of α, and hene the hange in the sign of α does not introdue physially new situation. The path of the photon whih is inident upon the surfae at a negative angle α will be symmetri with respet to the surfae normal, to the path of the photon orresponding to a positive inident angle α. Therefore, in the following, we onsider only positive inident angles 0 α < 90 ). When the medium moves in the positive diretion of the y-axis, u/ > 0 and the right-hand side of Eq. 4) is always positive, whih means that the angle of refration β is between 0 and 90. The dependene of β on the inident angle α for n = 1.5 and for various positive values of u/ is shown in Fig. 5. We see that as u/ varies from 0 to 1, the urves orresponding to different values of u/ are regularly shifted toward the smaller values of β, whih is expetable if one takes into aount the dragging effet of the moving medium. We emphasize that the portions of the urves drawn with dashed lines orrespond to inident angles at whih the photon will never reah the interfae. Namely, for a given speed u of the medium in the positive y-diretion, there exists an interval of values for the inident angle α α max α < 90 ) for whih the y-projetion of the photon s veloity osα) is less than FIG. 6: The angle of refration β as a funtion on the inident angle α for n = 1.5 when the material half-spae is moving perpendiular to the interfae, in the negative diretion of the y-axis. When κ 1 < u/ < 0 the refration is regular, and the angle of refration β of the photon is in the interval 0, 90 ]. For κ < u/ < κ 1, the refration is: a) regular 0 β < 90 ), if 0 α < α ; and b) bakward 90 < β < 180), if α < α < 90. If u/ < κ, the refration is always bakward. or equal to the veloity of the medium. The value of α max follows from: osα max = u/, 43) and it an be shown that it is exatly the inident angle at whih the urve β = βα) attains the maximum value. Hene, for inident angles belonging to the interval α max, 90 ], the refration beomes impossible. Translating to the frame in whih the medium is at rest, and thus taking into aount the aberration phenomenon, we have a photon inident upon the interfae at angles larger than 90. The situation hanges dramatially when u/ < 0. In this ase, the denominator on the right-hand side of Eq. 4) an take zero or even negative values, whih implies that the angle of refration β an beome equal or greater than 90. In Fig. 6 we plotted β as a funtion of α for n = 1.5 and for various negative values of u/. What we see in Fig. 6 an be summarized as follows. When the medium moves in the negative diretion of the y-axis up to some value of u/ denoted as κ 1, the refration seems to be regular. In this range of speeds κ 1 < u/ < 0), as the absolute value of the speed inreases, the urves are shifted toward higher values of β, and β remains smaller than 90 for eah α. When u/ < κ 1, a speifi value for the inident angle α appears for whih the angle of refration of the photon beomes equal to 90. We denote this ritial inident angle by α. For α < α < 90,

8 8 the angle of refration β is greater than 90. By further inreasing the absolute value of u/, the ritial angle α dereases, and for some value of u/, denoted by κ, α beomes 0. For u/ < κ, the angle of refration β is greater than 90 for all values of α. The quantities κ 1 and κ are funtions on the refrative index n of the medium in its rest frame, and as n inreases, the values of κ 1 and κ are beoming loser to eah other. It an be shown see Appendix B) that κ 1 = 1+n ) 1/ and κ = 1/n in our ase, for n = 1.5, κ and κ 0.67). Evidently, the value κ orresponds to the smallest superluminal speed of the medium. The existene of angles of refration greater than 90 in the ase u/ < 0 implies that the photon, instead of moving in the same general diretion as before the refration, moves bakward. At first sight, it might seem that the photon was refleted at the interfae instead of being refrated, but this is not what had happened see Fig. 7). One the photon strikes the vauum-material interfae, it undergoes a refration, enters the medium and ontinues to move into the medium at all times. The urious behavior of the refrated photon is due to the fat that the photon, while penetrating into the medium, is dragged by the overall motion of the medium. This dragging effet is a ause for the y-omponent of the photon s veloity in the moving medium to be negative, and hene, the refrated photon will move in the same general diretion as the medium. We emphasize that the path of the photon that undergoes a bakward refration does not oinide with the path whih it would trae if it were refleted from the interfae see Fig. 7). The path of the bakwardly refrated photon is fully determined by Eq. 4), and the path of the refleted photon obeys the Einstein s formula for refletion from a uniformly moving mirror see Ref. 8, 39, 40]). In the language of modern quantum eletrodynamis, whether the photon will be refleted or refrated at the inerfae is a matter of probability. In this sense, if a beam of photons is inident on the moving interfae at an angle α, and if the onditions for the appearane of the bakward refration are fulfilled, some of the photons will be refleted, and some bakwardly refrated. Therefore, in a more realisti view in whih the refletion is onsidered also, two different beams spreading bakward will be identified. In this way, we have a full onsisteny between the observations with respet to both inertial frames, S and S. Namely, for the observer in S -frame to whih the medium is at rest, both proesses are developing in an ordinary way: the refrated photon moves forward, obeying the Snell s law of refration, while the refleted one moves bakward, aording to the usual law of refletion: inident angle equals the refleted angle. Hene, the observer in S learly distinguishes two rays, and onsequently, two rays must exist with respet to the observer in the S-frame: one refleted, and one bakwardly refrated. V. CONCLUDING REMARKS We have presented a novel approah to the problem of refration of a light ray at an interfae between two homogeneous, isotropi and non-dispersive transparent optial materials in uniform retilinear motion. Our method amalgamates the original Fermat s priniple and the fat that an isotropi optial material at rest beomes optially anisotropi if it is moving at a onstant veloity. The derivation is in the framework of basi optis and speial relativity, requiring some knowledge of alulus at an elementary level. As suh, it may be regarded as an instrutive addendum to the typial introdutory physis ourses. We have analyzed in details the refration at a vauummaterial interfae, when the material half-spae is moving uniformly parallel or uniformly perpendiular to the interfae. In both ases, we onsidered the plane of inidene to be normal to the interfae and parallel to the veloity of the medium. The oinidene of the obtained refration formulas with the ones obtained with the other methods onfirms the validity of Fermat s priniple in the presene of uniformly moving boundaries and media. In addition, when the material half-spae moves perpendiular to the interfae, in the diretion whih is opposite to that of the inident light, we observed a bakward refration of light. To our knowledge, this optial effet has not been notied before in the literature. We emphasize that the refration formulas obtained in this paper refer to refration of a photon, and thus, desribe the refration of a light ray. We note that in a set of papers on the subjet of refletion and transmission of a plane eletromagneti wave by a moving medium, starting with the pioneering papers by Tai and Yeh, the refration formulas for different setups were derived by using the notions from the lassial eletromagneti theory and the Lorentz transformation 41, 4, 43, 44, 45, 46, 47]. Although these publiations also inlude the ases disussed in our Setions III and IV, the resulting refration formulas do not math the refration formulas in Eq. ) and 4) in our paper 48]. This mismath is due to the fat that Tai, Yeh and their suessors onsidered the angles of inidene, refletion and refration of the light wave as the angles between the orresponding propagation vetors and the interfae normal. However, we have shown that an optial isotropi medium in uniform retilinear motion beomes an optial anisotropi medium. Therefore, the propagation vetor in the moving medium will, in general, not oinide with the light ray in the moving medium. In onlusion, the refration formulas by Tai, Yeh and their followers do not refer to refration of a light ray. The approah in this paper an be extended to inlude more sophistiated situations, speifially, by taking into aount the dispersion of the material, or onsidering various three-dimensional ases when the plane of inidene is not neessary normal to the interfae and not parallel to the veloity of the medium. We also suggest a rep-

9 FIG. 7: A spae-time simulation depiting the differene between the refletion and the bakward refration. The material half-spae the shaded region) is moving in the negative y-diretion at a speed 0.7. The refrative index of the material half-spae in its rest frame is n = 1.5. The values of x and y oordinates are given in meters. a) The situation at time t = 0 ns. A photon, loated at point 0., 0.), is inident upon the vauum-material interfae, making an angle of 30 with respet to the interfae normal in the vauum half-spae. The interfae is loated at y = ; b) The path of the photon and the position of the interfae at time t = ns; ) The situation at t = 3 ns. The photon reahes the interfae at point , ); d) The situation at t = 4 ns. The ray br orresponds to the path of the bakwardly refrated photon, making an angle of with respet to the interfae normal in the material half-spae. The speed of the bakwardly refrated photon is 0.151, and its y veloity omponent is The ray r is the path of the photon if it were refleted from the interfae. The angle of refletion is with respet to the interfae normal in the vauum half-spae. The y veloity omponent of the refleted photon is ; e) The situation at t = 5 ns; f) The situation at t = 6 ns. 9

10 10 etition of the derivations in Se. III and IV when the refration ours at an interfae between two uniformly moving non-vauum optial materials having different refrative indies in their rest frames of referene. Aknowledgments The authors are indebted to I. Carusotto for onstrutive omments that substantially improved the disussions in Setions IV and V. We also thank J. Hamel, N. Novkovski, K. Trenčevski and V. Urumov for a ritial reading of the manusript, and V. Jovanov for helping us with the referenes. APPENDIX A: A DERIVATION OF EQ. 4) BY USING THE RELATIVISTIC VELOCITY TRANSFORMATION FORMULAS The refration formula in Eq. 4) an also be obtained by omparing the path of the photon in the S - frame where the medium is stationary, to the orresponding path of the photon in the S-frame where the medium is in uniform retilinear motion. We first desribe the refletion with respet to S -frame. Let the vauum-material interfae oinide with x -axis, with y > 0 being the region of the material half-spae, and y < 0 the region of the vauum half-spae. Let the photon be inident on the interfae from the vauum half-spae in a diretion making a positive angle α with the interfae normal. The photon is refrated in the material half-spae at an angle β with respet to the interfae normal, obeying the Snell s law of refration: sin α = n sin β, A1) where n is the refrative index of the medium at rest. Taking that the inident photon moves at a speed of light in vauum, and the refrated photon at a speed /n, the veloity omponents of the photon in the S -frame are: v x 1 = sin α, v y 1 = osα, v x = /n)sinβ, v y = /n)osβ, A) A3) A4) A5) where index 1 orresponds to the omponents of the inident photon, and index to the omponents of the refrated photon for simpliity, we use v x and v y instead of v x and v y ). With respet to S-frame where the material half-spae is moving at a onstant speed u in the positive y-diretion, the photon is inident on the moving interfae at an angle α at the same onstant speed as in the S -frame, but its speed after being refrated from the interfae is v β whih is a funtion of the angle of refration β. Thus, the orresponding veloity omponents of the photon in S-frame are: v x1 = sinα, v y1 = osα, v x = v β sinβ, v y = v β osβ. A6) A7) A8) A9) Dividing Eq. A8) with Eq. A9), and using the relativisti transformation formulas for the veloity omponents between S and S : we obtain: v x = v x1 u / ) 1/ 1 + v yu/, A10) v y = v y + u 1 + v y A11) u/, tanβ = v x 1 u / ) 1/ v y. A1) + u Combining Eqs. A4) and A5) with Eq. A1), we obtain: v x = /n )sinα, v y = /n)1 sin α /n ) 1/, whih we substitute into Eq. A1) to get: tan β = /n )sin α 1 u / ) 1/ /n)1 sin α /n ) 1/ + u. Using the veloity transformation formula: A13) A14) A15) v x = v x1 u / ) 1/ 1 v y u/ A16) and Eqs. A), A6) and A7), we get: sinα = sinα1 u / ) 1/ 1 u/)osα. A17) We substitute Eq. A17) into Eq. A15) and simplify the result to obtain the law of refration of the photon in Eq. 4). APPENDIX B: CONDITIONS FOR THE APPEARANCE OF THE BACKWARD REFRACTION Beause a bakward refration ours only for u/ < 0, we rewrite Eq. 4) by putting κ instead of u/ κ = u/ ): tan β = 1 κ )sin α κn 1 + κ osα) } 1/ 1 + n 1 + κ osα] 1 κ )sin α]. B1)

11 11 The photon will be refrated at β = 90 when α = α, whih is the ase when the denominator in the right-hand side of Eq. B1) equals zero: n 1 + κ osα ] 1 κ )1 os α )] 1/ = = κn 1 + κ osα ). B) The last equation is an irrational equation in osα, and its solution is equivalent to the solution of the system: 1 κ + κ n 1 κ n )] os α +κn 1 κ n ) osα +1 κ n 1 κ n ) B3) = 0 κn 1 + κ osα ) 0 whih simplifies to: osα ) 1, = κn 1 κ n ) ± 1 κ ) 1 n 1 κ n ) ] 1/ } 1 κ + κ n 1 κ n ) ] 1 osα 1/κ B4) The negative values for osα imply α > 90. Sine these solutions are physially unaeptable, we rejet the solution with the sign in the first equation of the system B4), beause this solution orresponds to osα < 0. What is left to be further analyzed is the expression: osα = κn 1 κ n ) + 1 κ ) 1 n 1 κ n ) ] } 1/ 1 κ + κ n 1 κ n ) ] 1. B5) The right-hand side in Eq. B5) is a real number if the expression under the square-root is non-negative, i.e. 1 n 1 κ n ) 0. Hene κ n 1) 1/ /n = κ min. Speifially, when κ = κ min, Eq. B5) redues to osα = n 1) 1/ /n. Sine 0 α 90, the right-hand side in Eq. B5) should satisfy the inequality: 0 κn 1 κ n ) + 1 κ ) 1 n 1 κ n ) ] 1/ } 1 κ + κ n 1 κ n ) ] 1 1. B6) For a fixed n, one an determine the interval of values for κ, κ κ 1, κ ], suh that for eah κ within this interval, Eq. B5) has a physially aeptable solution for α. This is the ritial inident angle at whih the bakward refration atually appears, and for α < α the refration is regular forward), while for α > α the refration is bakward. The ends of the interval κ 1, κ ] orrespond to the ritial inident angles α = 90 and α = 0, respetively, and an be obtained from the following equations: 1 κ 1 ) 1 n 1 κ 1n ) ] 1/ = κ 1 n 1 κ 1n ), B7) 1 n 1 κ n ) ] 1/ = 1 κ 1 n 1 κ n )] 1 κ. B8) The last two equations are irrational equations in κ 1,. To avoid umbersome alulations, we substitute x = n 1 κ 1n ) in Eq. B7) and x = 1 n 1 κ n ) in Eq. B8) to obtain: 1 x ) 1/ = κ 1x 1 κ, 1 B9) x 1/ = 1 κ x. 1 κ B10) Equation B9) is equivalent to the following system: κ 1 x + ) 1 κ 1 x 1 + κ 1 = 0 x B11) 0 while Eq. B10) orresponds to: κ x ) 1 + κ x + 1 = 0 x B1) 1/κ The quadrati equation in Eq. B11) has two solutions in x : x 1 = 1 κ 1 and x = 1 κ 1 )/κ 1. The solution x is not onsistent with the inequality in Eq. B11), sine x < 0 we have taken into aount that 0 < κ < 1). Thus, the only solution of the system B11) is x 1 = 1 κ 1, leading to the equation n 1 κ 1 n ) = 1 κ 1, whih gives: κ 1 = 1 + n ) 1/. B13) Similarly, by solving the quadrati equation in Eq. B1) one obtains two solutions in x : x 1 = 1 and x = 1/κ. Sine the solution x does not satisfy the inequality in Eq. B1), we rejet it and onsider only the solution x 1 = 1, whih leads to 1 n 1 κ n ) = 1, and hene: κ = 1/n. B14) Equations B13) and B14) determine the length of the interval κ 1, κ ] and its dependene on n. We an now summarize the results. For κ < 1 + n ) 1/, the refration is regular forward) for any inident angle α. When 1 + n ) 1/ κ 1/n, there exists a ritial inident angle α, suh that for α in the interval 0 < α < α the refration is regular forward), while for α < α < 90 the refration is bakward. The angle α an be alulated from Eq. B5). It depends on κ, and tends to zero as κ approahes 1/n. When κ = 1/n, α = 0, and the right-hand side of Eq. B1) beomes an undetermined expression of type 0/0, whih an be resolved

12 1 by using L Hospital s rule. When κ > 1/n, the denominator in Eq. B1) is negative for any α. On the other hand, the numerator in Eq. B1) is always non-negative. Therefore, for κ > 1/n, β > 90 for eah α. Hene, when the medium is moving in the negative diretion of the y-axis, there exists a ritial speed for the medium u = κ = /n, orresponding to the treshold for the superluminal speeds. Beyond this ritial speed, the photon will always be bakwardly refrated. 1] R. Dugas, A History of Mehanis Dover, New York, 1988), pp ] E. Heht, Optis Addison-Wesley,Reading, 1990), nd ed., pp ] M. Born and E. Wolf, Priniples of Optis Cambridge U.P., Cambridge, 1999), 7th ed., pp. xxv xxxiii, ] R. P. Feynman, Rev. Mod. Phys. 0, ). 5] R. P. Feynman, Quantum Eletrodynamis W. A. Benjamin, New York, 196). 6] R. P. Feynman and A. R. Hibbs, Quantum mehanis and path integrals MGraw-Hill, New York, 1965). 7] V. Meshbey, E. Ragoza, D. Kosloff, U. Egozi and D. Wexler, Pure Appl. Geophys. 159, ). 8] G. Nolet ed.), Seismi tomography with appliations in global seismology and exploration geophysis Reidel, Dordreht, 1987). 9] J. G. Berryman, Phys. Rev. Lett. 6, ). 10] N. S. Shashidhar and G. V. Anand, J. Sound. Vib. 186, ). 11] A. Cho, Siene 31, ). 1] U. Leonhardt, Siene 31, ). 13] U. Leonhardt, New J. Phys. 8, ). 14] F. Giannoni, A. Masiello and P. Piione, J. Math. Phys. 43, ). 15] V. Faraoni, Astrophys. J. 398, ); O. Wuknitz and U. Sperhake, Phys. Rev. D 69, ). 16] I. Kovner, Astrophys. J. 351, ); G. Temple, Pro. Roy. So. London A 168, ); V. Perlik, Gen. Rel. Grav. 38, ). 17] U. Leonhardt and P. Piwniki, Contemp. Phys. 41, ). 18] U. Leonhardt and P. Piwniki, Phys. Rev. A 60, ). 19] U. Leonhardt and P. Piwniki, J. Mod. Opt. 48, ). 0] U. Leonhardt and T. G. Philbin, New J. Phys. 8, ). 1] We note that the original statement of Fermat s priniple is not quite orret, sine there are some ases in whih the transit time is a saddle point. See, e.g., Ref. ]. ] R. P. Feynman, R. B. Leighton and M. Sands, The Feynman Letures on Physis Addison-Wesley, Reading, 1989), Vol. 1, Chap. 6. 3] R. P. Feynman, QED, The Strange Theory of Light and Matter Prineton U.P., Prineton, 1985). 4] E. F. Taylor, S. Vokos, J. M. O Meara and N. S. Thornber, Comput. Phys. 1, ). 5] A. G. Voronovih and O. A. Godin, Phys. Rev. Lett. 91, ). 6] O. A. Godin and A. G. Voronovih, Pro. R. So. Lond. A 460, ). 7] O. A. Godin and A. G. Voronovih, in Proeedings of the International Seminar Days on Diffration, St. Petersburg, Russia, June 8 July 1, 005 St. Petersburg Univ., 005), pp ] A. Gjurhinovski, Am. J. Phys. 7, ). 9] G. P. Sastry and T. R. Ravuri, Am. J. Phys. 58, ). 30] Although the speed of eah photon on the wavefront of the pulse will be frame-invariant, the frequeny, the intensity and the angular distribution of the photons will depend on the referene frame from where the observation is made. See, e.g., Ref. 9]. 31] The invariane of the shape and the speed of expansion of the pulse also follows from the Lorentz-invariane of the equation for the wavefront of the pulse. 3] A. Gjurhinovski, Am. J. Phys. 7, ). 33] H. M. Shey, Am. J. Phys. 37, ). 34] B. M. Bolotovskii and S. N. Stolyarov, Usp. Fiz. Nauk 114, ). 35] R. V. Jones, J. Phys. A 4, L1 1971). 36] R. V. Jones, Pro. R. So. Lond. A 38, ). 37] For a transverse light drag in the ase of a dispersive media, see I. Carusotto, M. Artoni, G. C. La Roa and F. Bassani, Phys. Rev. A 68, ), and the referenes therein. 38] The values for the angles in a given half-spae are taken positive/negative if they are in the lokwise/ounterlokwise diretion from the interfae normal at that partiular half-spae. 39] A. Einstein, Ann. Phys. Leipzig) 17, ). 40] A. Gjurhinovski, Am. J. Phys. 7, ). 41] C. T. Tai, Refletion and refration of a plane eletromagneti wave at the boundary of a semi-infinite moving medium, paper presented at URSI Spring Meeting, Washington D. C. 1965). 4] C. T. Tai, Two boundary-value problems involving moving media, paper presented at URSI Spring Meeting, Washington D. C. 1965). 43] C. Yeh, J. Appl. Phys. 36, ). 44] C. Yeh, J. Appl. Phys. 37, ). 45] V. P. Pyati, J. Appl. Phys. 38, ). 46] T. Shiozawa, K. Hazama and N. Kumagai, J. Appl. Phys. 38, ). 47] Y. Huang, J. Appl. Phys. 76, ). 48] For example, in the ase when the material half-spae is moving parallel to the interfae, Yeh s refration formula predits that when the inident angle is zero, the refrated angle also will be zero for any value for the speed of the medium see Ref. 43] above). This seems ontrary to the observed transverse Fresnel-Fizeau light drag, whih suggests non-zero angles of refration even when the inident angle is zero. Atually, from an eletrodynamis viewpoint, Yeh onsidered the phase veloity of light. In this sense, beause of momentum onservation along the interfae, an inident beam at normal inidene must refrat into a refrated beam whose phase veloity is also normal to the inerfae. On the other hand, the

13 13 formulas in Eq. ) and 4) above are desribing the refration of a light ray. The light ray orresponds to the group veloity, and orretly desribes the path followed by light if the experiment is performed with a penil of laser light. See, e.g., Ref. 34] and 37].