Hover over the menu bar to pick a physics animation.
For Snell's law, the critical concept is that the microscopic phase variation at the refractive surface has to be the same on both sides of the interface. Otherwise we would have discontinuities of the sum of the incident, reflected, and refracted wave amplitudes. If we did have a discontinuity of this sum that would mean that we instantaneous change of wave energy there and that is a violation of physical principles. Another way of stating the same principle is that the discontinuities would not satisfy the partial differential equation for wave behavior . To show this phase matching in the most visually definite way, here we will use rectangular wave sectors which is equivalent to saying that the incident wave has a square wave amplitude profile. A square wave is very well approximated by a sum of sinusoidal waves each having multiples of the square wave frequency. Each sector will have a width of half the waveLength, λ . The colors of the sectors will alternate between red for positive amplitude and black for negative amplitude. Derivation of Snell's law for refractive optical media is given in Appendix 1.
First, we'll discuss the incoming wave train. Its propagation angle is θi where θi is measured clockwise with respect to the x axis. It will make a transition to a new wavelength at x=xs where (xs,y,z=-∞→+∞) is the location of the surface of the refractive medium.
It's convenient to start with the tops of the wave train. If the wave sector top starts at (xt0,yt0) then the bottom left of the sector (with the y coordinate increasing downward as is common for computers) will be described by
yt(xt)=yt0+tanθ(xt-xt0)
In particular the point of intersection of yt with xs isyts=yt0+tanθ(xs-xt0)
The bottom left of the wave sector will be
xb0=xt0-wisinθ
yb0=yt0-wicosθ
yb(xb)=yb0+tanθ(xb-xb0)
ybs=yb0+tanθ(xs-xb0)
Now we want to have the wave move to the right so we let the left side of each sector be
xtl=xt0+δscosθ
ytl=yt0+δssinθ)
xtr=xtl+λi2cosθ
ytr=ytl+λi2sinθ
xbl=xtl+wisinθ
ybl=ytl-wicosθ
xtr=xtl+λi2cosθ
ytr=ytl+λi2sinθ
Our next task is to find the intersections of ybl,ybr,ytl,andytr with xs as a function of δs which is parallel to the wave's propagation direction. The general function for y(x) is
y(x)=y0+tanθ(x-x0)
For the top of the left side of the sector we have:
x0=xt0+δscosθ
y0=yt0+δssinθ)
ytl(xs)=yt0+δssinθ-tanθ[xs-(xt0+δscosθ)]
ytr(xs)=ytl(xs)+λi2sinθ
When xt>xs clipping of the first wave sector will begin and the sector will become a six sided polygon where the intersection with xs is yI: The equation for yI will be:
yI=yb-(xs-xb)cotθ
At the substrate surface, the incident medium sectors change shape (morph) into hybrids with two refractive index parts . The hybrids are 6 sided polygons as seen in the animation. The 6 sided polygon represents both the part of the wave sector that is outside the refractive media and the part that is inside the refractive median. Both of these have the same phase relative to the source as is visually shown by their having the same color. This polygon will continue to be 6 sided until the hybrid sector in fully inside the refractive medium after which it will become 4 sided.It is also important to note that the beam width in the refractive media is
wr=wicosθrcosθi
This is a result of the wave's width being perpendicular to its propagation direction. The wave propagation is animated on Canvas 1.
The phase variation of the incident wave along the boundary is
ϕi(y)=nikiysinθi
where ki=2πλi are the wave vector and wavelength of the incident wave and ni is the index of refraction in the incident medium. In the refractive medium we havenrλr=niλi
where nr is the index of refraction of the refractive medium. so kr=2πλr=nrniki Let the angle of propagation of the wave in the refractive medium be θr. Then the phase variation of the refractive wave along the boundary will be:ϕr(y)=krysinθr ϕr(y)=nrnikiysinθr
The phase matching requirement is thatϕr(y)=nrnikiysinθr=kiy=ϕi(y)sinθi
Solving for sinθr we obtain Snell's law:sinθr=ninrsinθi
Snell's law is equally valid for sound waves where the speed ratio of the waves, crci, is the analog of the index ratio of the waves: ninr.