Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

Posts Tagged ‘transmission’

Polarization angles for normal transmission and reflection

Posted by peeterjoot on January 22, 2014

[Click here for a PDF of this post with nicer formatting]

Question: Polarization angles for normal transmission and reflection ([1] pr 9.14)

For normal incidence, without assuming that the reflected and transmitted waves have the same polarization as the incident wave, prove that this must be so.

Answer

Working with coordinates as illustrated in fig. 1.1, the incident wave can be assumed to have the form

fig 1.1: Normal incidence coordinates

 

\begin{aligned}\tilde{\mathbf{E}}_{\mathrm{I}} = E_{\mathrm{I}} e^{i (k z - \omega t)} \hat{\mathbf{x}}\end{aligned} \hspace{\stretch{1}}(1.0.1a)

\begin{aligned}\tilde{\mathbf{B}}_{\mathrm{I}} = \frac{1}{{v}} \hat{\mathbf{z}} \times \tilde{\mathbf{E}}_{\mathrm{I}} = \frac{1}{{v}} E_{\mathrm{I}} e^{i (k z - \omega t)} \hat{\mathbf{y}}.\end{aligned} \hspace{\stretch{1}}(1.0.1b)

Assuming a polarization \hat{\mathbf{n}} = \cos\theta \hat{\mathbf{x}} + \sin\theta \hat{\mathbf{y}} for the reflected wave, we have

\begin{aligned}\tilde{\mathbf{E}}_{\mathrm{R}} = E_{\mathrm{R}} e^{i (-k z - \omega t)} (\hat{\mathbf{x}} \cos\theta + \hat{\mathbf{y}} \sin\theta)\end{aligned} \hspace{\stretch{1}}(1.0.2a)

\begin{aligned}\tilde{\mathbf{B}}_{\mathrm{R}} = \frac{1}{{v}} (-\hat{\mathbf{z}}) \times \tilde{\mathbf{E}}_{\mathrm{R}} = \frac{1}{{v}} E_{\mathrm{R}} e^{i (-k z - \omega t)} (\hat{\mathbf{x}} \sin\theta - \hat{\mathbf{y}} \cos\theta).\end{aligned} \hspace{\stretch{1}}(1.0.2b)

And finally assuming a polarization \hat{\mathbf{n}} = \cos\phi \hat{\mathbf{x}} + \sin\phi \hat{\mathbf{y}} for the transmitted wave, we have

\begin{aligned}\tilde{\mathbf{E}}_{\mathrm{T}} = E_{\mathrm{T}} e^{i (k' z - \omega t)} (\hat{\mathbf{x}} \cos\phi + \hat{\mathbf{y}} \sin\phi)\end{aligned} \hspace{\stretch{1}}(1.0.3a)

\begin{aligned}\tilde{\mathbf{B}}_{\mathrm{T}} = \frac{1}{{v}} \hat{\mathbf{z}} \times \tilde{\mathbf{E}}_{\mathrm{T}} = \frac{1}{{v'}} E_{\mathrm{T}} e^{i (k' z - \omega t)} (-\hat{\mathbf{x}} \sin\phi + \hat{\mathbf{y}} \cos\phi).\end{aligned} \hspace{\stretch{1}}(1.0.3b)

With no components of any of the \tilde{\mathbf{E}} or \tilde{\mathbf{B}} waves in the \hat{\mathbf{z}} directions the boundary value conditions at z = 0 require the equality of the \hat{\mathbf{x}} and \hat{\mathbf{y}} components of

\begin{aligned}\left( \tilde{\mathbf{E}}_{\mathrm{I}} + \tilde{\mathbf{E}}_{\mathrm{R}} \right)_{x,y} = \left( \tilde{\mathbf{E}}_{\mathrm{T}} \right)_{x,y}\end{aligned} \hspace{\stretch{1}}(1.0.4a)

\begin{aligned} \left( \frac{1}{\mu} \left( \tilde{\mathbf{B}}_{\mathrm{I}} + \tilde{\mathbf{B}}_{\mathrm{R}} \right) \right)_{x,y} = \left( \frac{1}{\mu'} \tilde{\mathbf{B}}_{\mathrm{T}} \right)_{x,y}.\end{aligned} \hspace{\stretch{1}}(1.0.4b)

With \beta = \mu v/\mu' v', those components are

\begin{aligned}E_{\mathrm{I}} + E_{\mathrm{R}} \cos\theta = E_{\mathrm{T}} \cos\phi \end{aligned} \hspace{\stretch{1}}(1.0.5a)

\begin{aligned}E_{\mathrm{R}} \sin\theta = E_{\mathrm{T}} \sin\phi\end{aligned} \hspace{\stretch{1}}(1.0.5b)

\begin{aligned}E_{\mathrm{R}} \sin\theta = - \beta E_{\mathrm{T}} \sin\phi\end{aligned} \hspace{\stretch{1}}(1.0.5c)

\begin{aligned}E_{\mathrm{I}} - E_{\mathrm{R}} \cos\theta = \beta E_{\mathrm{T}} \cos\phi\end{aligned} \hspace{\stretch{1}}(1.0.5d)

Equality of eq. 1.0.5b, and eq. 1.0.5c require

\begin{aligned}- \beta E_{\mathrm{T}} \sin\phi = E_{\mathrm{T}} \sin\phi,\end{aligned} \hspace{\stretch{1}}(1.0.6)

or (\theta, \phi) \in \{(0, 0), (0, \pi), (\pi, 0), (\pi, \pi)\}. It turns out that all of these solutions correspond to the same physical waves. Let’s look at each in turn

  • (\theta, \phi) = (0, 0). The system eq. 1.0.5.5 is reduced to

    \begin{aligned}\begin{aligned}E_{\mathrm{I}} + E_{\mathrm{R}} &= E_{\mathrm{T}} \\ E_{\mathrm{I}} - E_{\mathrm{R}} &= \beta E_{\mathrm{T}},\end{aligned}\end{aligned} \hspace{\stretch{1}}(1.0.7)

    with solution

    \begin{aligned}\begin{aligned}\frac{E_{\mathrm{T}}}{E_{\mathrm{I}}} &= \frac{2}{1 + \beta} \\ \frac{E_{\mathrm{R}}}{E_{\mathrm{I}}} &= \frac{1 - \beta}{1 + \beta}.\end{aligned}\end{aligned} \hspace{\stretch{1}}(1.0.8)

  • (\theta, \phi) = (\pi, \pi). The system eq. 1.0.5.5 is reduced to

    \begin{aligned}\begin{aligned}E_{\mathrm{I}} - E_{\mathrm{R}} &= -E_{\mathrm{T}} \\ E_{\mathrm{I}} + E_{\mathrm{R}} &= -\beta E_{\mathrm{T}},\end{aligned}\end{aligned} \hspace{\stretch{1}}(1.0.9)

    with solution

    \begin{aligned}\begin{aligned}\frac{E_{\mathrm{T}}}{E_{\mathrm{I}}} &= -\frac{2}{1 + \beta} \\ \frac{E_{\mathrm{R}}}{E_{\mathrm{I}}} &= -\frac{1 - \beta}{1 + \beta}.\end{aligned}\end{aligned} \hspace{\stretch{1}}(1.0.10)

    Effectively the sign for the magnitude of the transmitted and reflected phasors is toggled, but the polarization vectors are also negated, with \hat{\mathbf{n}} = -\hat{\mathbf{x}}, and \hat{\mathbf{n}}' = -\hat{\mathbf{x}}. The resulting \tilde{\mathbf{E}}_{\mathrm{R}} and \tilde{\mathbf{E}}_{\mathrm{T}} are unchanged relative to those of the (0,0) solution above.

  • (\theta, \phi) = (0, \pi). The system eq. 1.0.5.5 is reduced to

    \begin{aligned}\begin{aligned}E_{\mathrm{I}} + E_{\mathrm{R}} &= -E_{\mathrm{T}} \\ E_{\mathrm{I}} - E_{\mathrm{R}} &= -\beta E_{\mathrm{T}},\end{aligned}\end{aligned} \hspace{\stretch{1}}(1.0.11)

    with solution

    \begin{aligned}\begin{aligned}\frac{E_{\mathrm{T}}}{E_{\mathrm{I}}} &= -\frac{2}{1 + \beta} \\ \frac{E_{\mathrm{R}}}{E_{\mathrm{I}}} &= \frac{1 - \beta}{1 + \beta}.\end{aligned}\end{aligned} \hspace{\stretch{1}}(1.0.12)

    Effectively the sign for the magnitude of the transmitted phasor is toggled. The polarization vectors in this case are \hat{\mathbf{n}} = \hat{\mathbf{x}}, and \hat{\mathbf{n}}' = -\hat{\mathbf{x}}, so the transmitted phasor magnitude change of sign does not change \tilde{\mathbf{E}}_{\mathrm{T}} relative to that of the (0,0) solution above.

  • (\theta, \phi) = (\pi, 0). The system eq. 1.0.5.5 is reduced to

    \begin{aligned}\begin{aligned}E_{\mathrm{I}} - E_{\mathrm{R}} &= E_{\mathrm{T}} \\ E_{\mathrm{I}} + E_{\mathrm{R}} &= \beta E_{\mathrm{T}},\end{aligned}\end{aligned} \hspace{\stretch{1}}(1.0.13)

    with solution

    \begin{aligned}\begin{aligned}\frac{E_{\mathrm{T}}}{E_{\mathrm{I}}} &= \frac{2}{1 + \beta} \\ \frac{E_{\mathrm{R}}}{E_{\mathrm{I}}} &= -\frac{1 - \beta}{1 + \beta}.\end{aligned}\end{aligned} \hspace{\stretch{1}}(1.0.14)

    This time, the sign for the magnitude of the reflected phasor is toggled. The polarization vectors in this case are \hat{\mathbf{n}} = -\hat{\mathbf{x}}, and \hat{\mathbf{n}}' = \hat{\mathbf{x}}. In this final variation the reflected phasor magnitude change of sign does not change \tilde{\mathbf{E}}_{\mathrm{R}} relative to that of the (0,0) solution.

We see that there is only one solution for the polarization angle of the transmitted and reflected waves relative to the incident wave. Although we fixed the incident polarization with \mathbf{E} along \hat{\mathbf{x}}, the polarization of the incident wave is maintained regardless of TE or TM labeling in this example, since our system is symmetric with respect to rotation.

References

[1] D.J. Griffith. Introduction to Electrodynamics. Prentice-Hall, 1981.

Advertisements

Posted in Math and Physics Learning. | Tagged: , , , , , | 1 Comment »

PHY456H1F: Quantum Mechanics II. Lecture 21 (Taught by Prof J.E. Sipe). Scattering theory

Posted by peeterjoot on November 24, 2011

[Click here for a PDF of this post with nicer formatting and figures if the post had any (especially if my latex to wordpress script has left FORMULA DOES NOT PARSE errors.)]

Disclaimer.

Peeter’s lecture notes from class. May not be entirely coherent.

Scattering theory.

READING: section 19, section 20 of the text [1].

Here’s (\ref{fig:qmTwoL21:qmTwoL21Fig1}) a simple classical picture of a two particle scattering collision

\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{qmTwoL21Fig1}
\caption{classical collision of particles.}
\end{figure}

We will focus on point particle elastic collisions (no energy lost in the collision). With particles of mass m_1 and m_2 we write for the total and reduced mass respectively

\begin{aligned}M = m_1 + m_2\end{aligned} \hspace{\stretch{1}}(2.1)

\begin{aligned}\frac{1}{{\mu}} = \frac{1}{{m_1}} + \frac{1}{{m_2}},\end{aligned} \hspace{\stretch{1}}(2.2)

so that interaction due to a potential V(\mathbf{r}_1 - \mathbf{r}_2) that depends on the difference in position \mathbf{r} = \mathbf{r}_1 - \mathbf{r} has, in the center of mass frame, the Hamiltonian

\begin{aligned}H = \frac{\mathbf{p}^2}{2 \mu} + V(\mathbf{r})\end{aligned} \hspace{\stretch{1}}(2.3)

In the classical picture we would investigate the scattering radius r_0 associated with the impact parameter \rho as depicted in figure (\ref{fig:qmTwoL21:qmTwoL21Fig2})

\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{qmTwoL21Fig2}
\caption{Classical scattering radius and impact parameter.}
\end{figure}

1D QM scattering. No potential wave packet time evolution.

Now lets move to the QM picture where we assume that we have a particle that can be represented as a wave packet as in figure (\ref{fig:qmTwoL21:qmTwoL21Fig3})
\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{qmTwoL21Fig3}
\caption{Wave packet for a particle wavefunction \Re(\psi(x,0))}
\end{figure}

First without any potential V(x) = 0, lets consider the evolution. Our position and momentum space representations are related by

\begin{aligned}\int {\left\lvert{\psi(x, t)}\right\rvert}^2 dx = 1 = \int {\left\lvert{\psi(p, t)}\right\rvert}^2 dp,\end{aligned} \hspace{\stretch{1}}(2.4)

and by Fourier transform

\begin{aligned}\psi(x, t) = \int \frac{dp}{\sqrt{2 \pi \hbar}} \overline{\psi}(p, t) e^{i p x/\hbar}.\end{aligned} \hspace{\stretch{1}}(2.5)

Schr\”{o}dinger’s equation takes the form

\begin{aligned}i \hbar \frac{\partial {\psi(x,t)}}{\partial {t}} = - \frac{\hbar^2}{2 \mu} \frac{\partial^2 {{\psi(x, t)}}}{\partial {{x}}^2},\end{aligned} \hspace{\stretch{1}}(2.6)

or more simply in momentum space

\begin{aligned}i \hbar \frac{\partial {\overline{\psi}(p,t)}}{\partial {t}} = \frac{p^2}{2 \mu} \frac{\partial^2 {{\overline{\psi}(p, t)}}}{\partial {{x}}^2}.\end{aligned} \hspace{\stretch{1}}(2.7)

Rearranging to integrate we have

\begin{aligned}\frac{\partial {\overline{\psi}}}{\partial {t}} = -\frac{i p^2}{2 \mu \hbar} \overline{\psi},\end{aligned} \hspace{\stretch{1}}(2.8)

and integrating

\begin{aligned}\ln \overline{\psi} = -\frac{i p^2 t}{2 \mu \hbar} + \ln C,\end{aligned} \hspace{\stretch{1}}(2.9)

or

\begin{aligned}\overline{\psi} = C e^{-\frac{i p^2 t}{2 \mu \hbar}} = \overline{\psi}(p, 0) e^{-\frac{i p^2 t}{2 \mu \hbar}}.\end{aligned} \hspace{\stretch{1}}(2.10)

Time evolution in momentum space for the free particle changes only the phase of the wavefunction, the momentum probability density of that particle.

Fourier transforming, we find our position space wavefunction to be

\begin{aligned}\psi(x, t) = \int \frac{dp}{\sqrt{2 \pi \hbar}} \overline{\psi}(p, 0) e^{i p x/\hbar} e^{-i p^2 t/2 \mu \hbar}.\end{aligned} \hspace{\stretch{1}}(2.11)

To clean things up, write

\begin{aligned}p = \hbar k,\end{aligned} \hspace{\stretch{1}}(2.12)

for

\begin{aligned}\psi(x, t) = \int \frac{dk}{\sqrt{2 \pi}} a(k, 0) ) e^{i k x} e^{-i \hbar k^2 t/2 \mu},\end{aligned} \hspace{\stretch{1}}(2.13)

where

\begin{aligned}a(k, 0) = \sqrt{\hbar} \overline{\psi}(p, 0).\end{aligned} \hspace{\stretch{1}}(2.14)

Putting

\begin{aligned}a(k, t) = a(k, 0) e^{ -i \hbar k^2/2 \mu},\end{aligned} \hspace{\stretch{1}}(2.15)

we have

\begin{aligned}\psi(x, t) = \int \frac{dk}{\sqrt{2 \pi}} a(k, t) ) e^{i k x} \end{aligned} \hspace{\stretch{1}}(2.16)

Observe that we have

\begin{aligned}\int dk {\left\lvert{ a(k, t)}\right\rvert}^2 = \int dp {\left\lvert{ \overline{\psi}(p, t)}\right\rvert}^2 = 1.\end{aligned} \hspace{\stretch{1}}(2.17)

A Gaussian wave packet

Suppose that we have, as depicted in figure (\ref{fig:qmTwoL21:qmTwoL21Fig4})
\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{qmTwoL21Fig4}
\caption{Gaussian wave packet.}
\end{figure}

a Gaussian wave packet of the form

\begin{aligned}\psi(x, 0) = \frac{ (\pi \Delta^2)^{1/4}} e^{i k_0 x} e^{- x^2/2 \Delta^2}.\end{aligned} \hspace{\stretch{1}}(2.18)

This is actually a minimum uncertainty packet with

\begin{aligned}\Delta x &= \frac{\Delta}{\sqrt{2}} \\ \Delta p &= \frac{\hbar}{\Delta \sqrt{2}}.\end{aligned} \hspace{\stretch{1}}(2.19)

Taking Fourier transforms we have

\begin{aligned}a(k, 0) &= \left(\frac{\Delta^2}{\pi}\right)^{1/4} e^{-(k - k_0)^2 \Delta^2/2} \\ a(k, t) &= \left(\frac{\Delta^2}{\pi}\right)^{1/4} e^{-(k - k_0)^2 \Delta^2/2} e^{ -i \hbar k^2 t/ 2\mu} \equiv \alpha(k, t)\end{aligned} \hspace{\stretch{1}}(2.21)

For t > 0 our wave packet will start moving and spreading as in figure (\ref{fig:qmTwoL21:qmTwoL21Fig5})
\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{qmTwoL21Fig5}
\caption{moving spreading Gaussian packet.}
\end{figure}

With a potential.

Now “switch on” a potential, still assuming a wave packet representation for the particle. With a positive (repulsive) potential as in figure (\ref{fig:qmTwoL21:qmTwoL21Fig6}), at a time long before the interaction of the wave packet with the potential we can visualize the packet as heading towards the barrier.

\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{qmTwoL21Fig6}
\caption{QM wave packet prior to interaction with repulsive potential.}
\end{figure}

After some time long after the interaction, classically for this sort of potential where the particle kinetic energy is less than the barrier “height”, we would have total reflection. In the QM case, we’ve seen before that we will have a reflected and a transmitted portion of the wave packet as depicted in figure (\ref{fig:qmTwoL21:qmTwoL21Fig7})
\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{qmTwoL21Fig7}
\caption{QM wave packet long after interaction with repulsive potential.}
\end{figure}

Even if the particle kinetic energy is greater than the barrier height, as in figure (\ref{fig:qmTwoL21:qmTwoL21Fig8}), we can still have a reflected component.
\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{qmTwoL21Fig8}
\caption{Kinetic energy greater than potential energy.}
\end{figure}

This is even true for a negative potential as depicted in figure (\ref{fig:qmTwoL21:qmTwoL21Fig9})!

\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{qmTwoL21Fig9}
\caption{qmTwoL21Fig9}
\end{figure}

Consider the probability for the particle to be found anywhere long after the interaction, summing over the transmitted and reflected wave functions, we have

\begin{aligned}1 &= \int {\left\lvert{\psi_r + \psi_t}\right\rvert}^2 \\ &= \int {\left\lvert{\psi_r}\right\rvert}^2  + \int {\left\lvert{\psi_t}\right\rvert}^2 + 2 \Re \int \psi_r^{*} \psi_t\end{aligned}

Observe that long after the interaction the cross terms in the probabilities will vanish because they are non-overlapping, leaving just the probably densities for the transmitted and reflected probably densities independently.

We define

\begin{aligned}T &= \int {\left\lvert{\psi_t(x, t)}\right\rvert}^2 dx \\ R &= \int {\left\lvert{\psi_r(x, t)}\right\rvert}^2 dx.\end{aligned} \hspace{\stretch{1}}(2.23)

The objective of most of our scattering problems will be the calculation of these probabilities and the comparisons of their ratios.

Question. Can we have more than one wave packet reflect off. Yes, we could have multiple wave packets for both the reflected and the transmitted portions. For example, if the potential has some internal structure there could be internal reflections before anything emerges on either side and things could get quite messy.

Considering the time independent case temporarily.

We are going to work through something that is going to seem at first to be completely unrelated. We will (eventually) see that this can be applied to this problem, so a bit of patience will be required.

We will be using the time independent Schr\”{o}dinger equation

\begin{aligned}- \frac{\hbar^2}{2 \mu} \psi_k''(x) = V(x) \psi_k(x) = E \psi_k(x),\end{aligned} \hspace{\stretch{1}}(3.25)

where we have added a subscript k to our wave function with the intention (later) of allowing this to vary. For “future use” we define for k > 0

\begin{aligned}E = \frac{\hbar^2 k^2}{2 \mu}.\end{aligned} \hspace{\stretch{1}}(3.26)

Consider a potential as in figure (\ref{fig:qmTwoL21:qmTwoL21Fig10}), where V(x) = 0 for x > x_2 and x < x_1.

\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{qmTwoL21Fig10}
\caption{potential zero outside of a specific region.}
\end{figure}

We won't have bound states here (repulsive potential). There will be many possible solutions, but we want to look for a solution that is of the form

\begin{aligned}\psi_k(x) = C e^{i k x}, \qquad x > x_2\end{aligned} \hspace{\stretch{1}}(3.27)

Suppose x = x_3 > x_2, we have

\begin{aligned}\psi_k(x_3) = C e^{i k x_3}\end{aligned} \hspace{\stretch{1}}(3.28)

\begin{aligned}{\left.{{\frac{d\psi_k}{dx}}}\right\vert}_{{x = x_3}} = i k C e^{i k x_3} \equiv \phi_k(x_3)\end{aligned} \hspace{\stretch{1}}(3.29)

\begin{aligned}{\left.{{\frac{d^2\psi_k}{dx^2}}}\right\vert}_{{x = x_3}} = -k^2 C e^{i k x_3} \end{aligned} \hspace{\stretch{1}}(3.30)

Defining

\begin{aligned}\phi_k(x) = \frac{d\psi_k}{dx},\end{aligned} \hspace{\stretch{1}}(3.31)

we write Schr\”{o}dinger’s equation as a pair of coupled first order equations

\begin{aligned}\frac{d\psi_k}{dx} &= \phi_k(x) \\ -\frac{\hbar^2}{2 \mu} \frac{d\phi_k(x)}{dx} = - V(x) \psi_k(x) + \frac{\hbar^2 k^2}{2\mu} \psi_k(x).\end{aligned} \hspace{\stretch{1}}(3.32)

At this x = x_3 specifically, we “know” both \phi_k(x_3) and \psi_k(x_3) and have

\begin{aligned}{\left.{{\frac{d\psi_k}{dx}}}\right\vert}_{{x_3}} &= \phi_k(x) \\ -\frac{\hbar^2}{2 \mu} {\left.{{\frac{d\phi_k(x)}{dx}}}\right\vert}_{{x_3}} = - V(x_3) \psi_k(x_3) + \frac{\hbar^2 k^2}{2\mu} \psi_k(x_3),\end{aligned} \hspace{\stretch{1}}(3.34)

This allows us to find both

\begin{aligned}{dx}}}\right\vert}_{{x_3}} \\ {dx}}}\right\vert}_{{x_3}} \end{aligned} \hspace{\stretch{1}}(3.36)

then proceed to numerically calculate \phi_k(x) and \psi_k(x) at neighboring points x = x_3 + \epsilon. Essentially, this allows us to numerically integrate backwards from x_3 to find the wave function at previous points for any sort of potential.

References

[1] BR Desai. Quantum mechanics with basic field theory. Cambridge University Press, 2009.

Posted in Math and Physics Learning. | Tagged: , , , , , , , , , , | Leave a Comment »