Kevin Schmidt
Department of Physics and Astronomy
Arizona State University
Tempe, AZ
(_{t},) = x . | (1) |
(_{t},) = x ( x ). | (2) |
(,t) = e^{ it t - it} | (3) |
k_{z} = | (4) |
Our source is a pulse starting far enough away from the material that it can be localized in the region z < z_{1} , and to have only Fourier components with real positive k_{z} so that they propagate toward the interface. In that case the reflected wave components above do not contribute when calculating the amplitudes A_{} and A_{ || } with the result
A_{}(_{t},) = d^{ 2}r_{t}dt^{- i + it} | (5) |
A_{ || }(_{t},) = d^{ 2}r_{t}dt^{- i + it} | (6) |
The algorithm is therefore, given the incident field, calculate the amplitudes A_{} and A_{ || } . For each value of _{t} and that occurs, solve Maxwell's equations for the transmission t and reflection r coefficients for a plane wave incident. Plug all these in to Eqs. 1 and 2 and integrate to get the reflected and transmitted fields. To get the fields in the material, the solutions in the material region need to be included in Eqs. 1 and 2.
Often the solution for a pulse that is a beam of nearly monochromatic radiation at nearly a single angle of incidence is desired. In addition, often the exiting pulses are the same shape as the incoming one just time delayed and/or shifted along the transverse direction. Let's see if we can see how this can come about. If the pulse is strongly peaked around a particular _{t} and , we can expand the r and t functions around these values. As we will see, it is more convenient to write them as expansions in the exponent.
Initially, let's ignore the vector character of the waves. An incoming pulse can be written as
F_{0}(,t) = e^{ i - it}A(_{t},) | (7) |
F_{T}(,t) | = | t(_{t},)e^{ i - it}A(_{t},) | (8) |
F_{T}(,t) | = | e^{ i(t,)}e^{ i - it}A(_{t},) | |
= | e^{ i(t(0),) + i0 (t - t(0)) + i0( - ) + ...}e^{ i - it}A(_{t},) . |
For now we neglect the quadratic terms in the exponent, but if they
are large this will be a poor approximation. We will come back to this
point later. Notice that the lowest order term in the expansion
e^{ i(t(0),)} is an overall constant
that simply multiplies the wave packet. Similarly, the
terms multiplied by
_{t}^{(0)} and
are
also constants. If we take the case where
the first derivatives of are real, we can approximate the
transmitted amplitude by
F_{T}(,t) | e^{ i}t(_{t}^{(0)},)e^{ i i}A(_{t},) | ||
= | e^{ i}t(_{t}^{(0)},)F_{0} | (9) |
Repeating the calculation for the reflected wave gives in the same approximation and writing r(,) = e^{ i(t,)}
F_{R}(,t) = e^{ i}t(_{t}^{(0)},)F_{0} | (10) |
Before looking at these results in more detail lets examine the neglected terms. If we calculate the width of the packet in k_{t} or , the second derivative terms of or will change this. However the magnitude of these terms is of order the second derivative times factors of k_{t}^{(0)} or w^{ (0)} (there are also cross terms from the mixed derivatives). The requirement that these can be neglected is therefore that these are small. That is the change in the first derivatives is small over the range of frequencies and wave vectors contained in the packet. If these terms are included but higher order terms are neglected, and the packet also expanded in this way, we would have a quadratic form in the exponent. Diagonalizing this quadratic form would give terms that would contribute both to a new width of the packets and to a shifted position from the linear terms and from completing the square. The other neglected terms are the real parts of the derivatives of and . Including them would again change the shape of the packet. However, in this case, the real part of the zeroth order term dominates.
To reiterate, we expand the transmission and reflection amplitudes. The key point is that a large overall phase does not change the energy in the wavep packet therefore for the imaginary part of the exponent, we need to keep the gradient terms since they are the lowest order components that change the energy in the wave packet. For the real parts, the constant term dominates the gradient and gives the lowest order component. Therefore, for a wave packet with a sufficiently narrow frequency range, so that only the lowest order physical terms need to be kept, Eqs. 10 and 11 with and given by the imaginary parts of the exponent give the transmitted and reflected packets.
At this point, we can go back to the vector equations. Notice that the polarization factors are real. Therefore, within the same approximation as the expansions above, we can keep just the lowest order terms and use the polarizations with _{0}
Let's take a material has an index of refraction n() which is independent of z except near z_{1} and z_{2} goes smoothly from 1 to n() 1 so that reflections at the boundaries can be ignored. In this case, each component will smoothly change index of refraction through the material, and if z_{2} - z_{1} = L is the thickness, the transmission coefficient will be
t(_{t},) = e^{ iL - iL} . | (11) |
d_{t} = - 0 0 | (12) |
d_{t} = L, | (13) |
Next look at the time delay of the packet through the medium,
= - | (14) |
v_{g} = c 0 | (15) |
Let's look at the Goos-Hänchen effect. Here we have a beam undergoing total internal reflection. Placing the interface at z = 0 , the reflection coefficient has, for a real index of refraction, unit magnitude and is a phase. The reflection coefficient for perpendicular polarization in terms of the incident angle and the critical angle sin is
_{} = | (16) |
= - 2. | (17) |
|_{kt}| = = | (18) |
d_{} = |_{kt}| = . | (19) |
D_{} = , | (20) |