Processing math: 100%

Hover over the menu bar to pick a physics animation.

Snell's Law of Refraction

Introduction

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.

Geometry

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 is

yts=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)

where wi is the width of the incident beam. In particular, the sector bottom left will intersect xs at

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θ)

where δs is the displacement of the left side from its starting location. Then the top right side of the sector becomes:

xtr=xtl+λi2cosθ
ytr=ytl+λi2sinθ

The analogous values for (xtl,ytl) become:

xbl=xtl+wisinθ
ybl=ytl-wicosθ

and for (xtr,ytr)

xtr=xtl+λi2cosθ
ytr=ytl+λi2sinθ

Displacement of the sector by δs

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θ)

so, noting that ytl(x) is a function of x and different from y_(tl) our expression for ytl(x) becomes:

ytl(xs)=yt0+δssinθ-tanθ[xs-(xt0+δscosθ)]

For the top right side of the sector we have:

ytr(xs)=ytl(xs)+λi2sinθ

each of these equations can be solved for y(xs) and often the values will be outside the wave envelope.

Effects on Wave Sectors when They Encounter the Refractive Surface

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.

Appendix 1: Snell's Law due to Microscopic Phase Matching At Surface

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 have

nrλ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.