Dirac Quantum Beats in Extreme Relativistic Diffraction in Time

Abstract

We solve the relativistic Dirac equation for the shutter problem. We prove that, at intermediate relativistic energies, the probability density oscillates in time in a similar way to the optical Fresnel oscillations in a straight edge. However, for extreme-relativistic beams, the Fresnel oscillations turn into quantum damped beats.

Share and Cite:

Godoy, S. (2022) Dirac Quantum Beats in Extreme Relativistic Diffraction in Time. Journal of Applied Mathematics and Physics, 10, 1711-1720. doi: 10.4236/jamp.2022.105119.

1. Introduction

Similarities between optics and quantum mechanics have long been recognized. One example of this symmetry was obtained by Moshinsky [1] who addressed the following quantum shutter problem. Consider a monoenergetic beam of free particles, moving parallel to the x-axis. For negative times, the beam is interrupted at $x=0$ by a perfectly absorbing shutter perpendicular to the beam. Suddenly, at time $t=0$, the shutter is opened, allowing for $t>0$ the free time-evolution of the beam of particles. What is the time-dependent density observed at a distance x from the shutter?

The shutter problem implies solving the time-dependent Schrödinger equation with an initial condition given by

$\psi \left(x,0\right)={\text{e}}^{ikx}\theta \left(-x\right),$ (1)

where $\theta \left(x\right)$ denotes the step function. For $t>0$, Moshinsky proved that the free propagation of the beam has a probability density, $\rho \left(x,t\right)={|\psi \left(x,t\right)|}^{2}$, given by:

$\rho \left(x,t\right)=\frac{1}{2}{\left[C\left(\xi \right)+\frac{1}{2}\right]}^{2}+\frac{1}{2}{\left[S\left(\xi \right)+\frac{1}{2}\right]}^{2},$ (2)

here, $C\left(\xi \right)+iS\left(\xi \right)\equiv {\int }_{0}^{\xi }\mathrm{exp}\left(i\pi {u}^{2}/2\right)\text{d}u$, denotes the complex Fresnel function and the argument $\xi$ is given by

$\xi \left(x,t\right)\equiv \sqrt{\frac{m}{\pi \hslash t}}\left(\frac{\hslash k}{m}t-x\right).$ (3)

The right-hand side in Equation (2) looks similar to the mathematical expression for the light intensity in the optical Fresnel diffraction by a straight-edge, Born-Wolf [2]. For a fixed position, $X\equiv px/\hslash =1$, the plot of the probability density $\rho \left(X,T\right)$ as a function of time, $T\equiv Et/\hslash$, is shown in Figure 1.

A good measure of the “width” in time of this diffraction effect, can be obtained from the difference $\Delta t\equiv {t}_{2}-{t}_{1}$ between the first two times at which $\rho$ takes the (classical) mean value. We obtain, for thermal neutrons, $\Delta t\approx {10}^{-9}\text{ }\text{ }\text{s}$. The experimental evidence of this quantum prediction has been confirmed until very recently by Szriftigiser, Guéry-Odelin, Arndt, and Dalibard [3].

These transient oscillations are a pure quantum phenomenon, and similar oscillations arise at the moment of closing and opening gates in nanoscopic circuits [4]. For a review on the subject see [5] [6]. With adequate potentials added to the model, it has been used to study transient dynamics of tunneling matter waves [7] [8] [9] [10], and the transient response to abrupt changes of the interaction potential in semiconductor structures and quantum dots [11] [12]. There is, in summary, a strong motivation for a thorough understanding of transient time oscillation in beams of matter.

The analogy between the quantum shutter problem and optical diffraction raises the question of whether the transient densities for other types of wave equations show this analogy. Using the Dirac equation for the shutter problem at low energies, we found in a previous paper [13], that the algebraic expression for the relativistic quantum density $\rho \left(x,t\right)$ no longer resembles the algebraic Fresnel expression for the optical diffraction. In spite of this, when the exact relativistic density is plotted versus time, the plots show transient oscillations which prove that diffraction in time is present in the relativistic realm.

Figure 1. Schrödinger diffraction in time, $X\equiv px/\hslash =1$.

The main contribution of this paper is to show the existence of transient quantum beats in extreme-relativistic diffraction in time. In Section 2, we solve the Dirac equation for the shutter problem and in Section 3 we prove, as expected, that for intermediate energies, at internuclear distances, the plot of the probability density looks similar to the optical Fresnel oscillations derived in the Schrödinger equation. Finally, in Section 4, as we gradually increase the energy until we reach extreme-relativistic energies, the Dirac diffraction changes gradually from a Fresnel pattern into quantum damped beats! For 1/2-spin particles, this has never been reported before. We conclude with the application of the present result to the quantum problem of extreme-relativistic particles suddenly released from an infinite potential well.

2. The Dirac Shutter Problem

For 1/2-spin particles, we assume for all negative times, $t\le 0$, a beam of right-moving free particles in the left side of a perfectly absorbing shutter located at the origin, and none to the right. For the case of propagation into the direction of the z-axis, we want to calculate the spinor wave function, $\psi \left(z,t\right)={\left({\psi }_{1},{\psi }_{2},{\psi }_{3},{\psi }_{4}\right)}^{\text{T}}$, which is the solution of the one-dimensional Dirac equation:

$\left\{\left(\begin{array}{cc}I& 0\\ 0& I\end{array}\right)\frac{1}{i}\frac{\partial }{\partial \left(ct\right)}+\left(\begin{array}{cc}0& {\sigma }_{z}\\ {\sigma }_{z}& 0\end{array}\right)\frac{1}{i}\frac{\partial }{\partial z}+\frac{1}{{\lambda }_{c}}\left(\begin{array}{cc}I& 0\\ 0& -I\end{array}\right)\right\}\psi =0,$ (4)

here ${\sigma }_{z}$ is the Pauli matrix and, ${\lambda }_{c}\equiv \hslash /{m}_{0}c$, the reduced Compton wave-length. Three quantum numbers are needed to classify the Dirac free-particle solutions, namely, the momentum, $p\equiv \hslash k$, positive or negative energies, $E=±\hslash \omega$, where $\omega =c{\left({k}^{2}+{\lambda }_{c}^{-2}\right)}^{1/2}$, and the helicity, ${\Lambda }_{S}=S\cdot p/p$. In our case, we assume that the initial condition corresponds to a right moving plane wave propagating along the $z$ -direction, $k=\left(0,0,k\right)$, a positive energy, $E=+\hslash \omega$, and the spin parallel to the direction of motion, ${S}_{z}=+1/2$. The incident plane wave is then given by:

${\psi }_{k,+E,+1/2}\left(z,t\right)=N{\left(1,0,\frac{k}{{\lambda }_{c}^{-1}+\omega /c},0\right)}^{\text{T}}\text{\hspace{0.17em}}{\text{e}}^{i\left(kz-\omega t\right)}\theta \left(-z\right)\text{ }\left(t\le 0\right),$ (5)

where N is the normalization factor, and $\theta \left(z\right)$ denotes the Heaviside step function defined in Equation (1).

For free particles, the helicity ${\Lambda }_{S}$ is a constant of motion. The initial direction of spin, ${S}_{z}=1/2$, will be conserved. Therefore, the two components of the wave function: ( ${\psi }_{2}$ and ${\psi }_{4}$ ) which are zero at the initial time will remain zero for all positive times: ${\psi }_{2}\left(z,t\right)={\psi }_{4}\left(z,t\right)=0\text{\hspace{0.17em}}\left(t\ge 0\right)$. The two remaining components: $\psi \equiv {\left({\psi }_{1},{\psi }_{3}\right)}^{\text{T}}$ evolve in time as solutions of the equation,

$\left(I\frac{\partial }{\partial \left(ct\right)}+{\sigma }_{x}\frac{\partial }{\partial z}+\frac{i}{{\lambda }_{c}}{\sigma }_{z}\right)\psi \left(z,t\right)=0$ (6)

with the initial condition:

$\psi \left(z,0\right)\equiv {\left(\begin{array}{cc}{\psi }_{1}\left(z,0\right)& {\psi }_{3}\left(z,0\right)\end{array}\right)}^{\text{T}}={\left(\begin{array}{cc}N& \gamma \end{array}\right)}^{\text{T}}{\text{e}}^{ikz}\theta \left(-z\right),$ (7)

where $\gamma \equiv Nk/\left({\lambda }_{c}^{-1}+\omega /c\right)$. Using the condition: ${N}^{2}+{\gamma }^{2}=1$, the normalization factor becomes: ${N}^{2}=\left(1+c/{\lambda }_{c}\omega \right)/2$.

For the sake of conciseness, we denote, as usual $\beta \equiv v/c$, and use the Compton wave-length ${\lambda }_{c}$, to define dimensionless variables:

$Z\equiv \frac{z}{{\lambda }_{c}},\text{ }T\equiv \frac{ct}{{\lambda }_{c}},\text{ }\kappa \equiv k{\lambda }_{c}=\frac{p}{{m}_{0}c},\text{ }\Omega \equiv \frac{\omega {\lambda }_{c}}{c}=\frac{E}{{m}_{0}{c}^{2}}$

${N}^{2}=\left(1+\sqrt{1-{\beta }^{2}}\right)/2,\text{ }{\gamma }^{2}=\left(1-\sqrt{1-{\beta }^{2}}\right)/2.$ (8)

Using these variables we derive in Appendix A, for the Dirac shutter, the exact transmitted wave, ${\psi }_{>}\left(Z,T\right)$, valid for ( $Z\ge 0$ ):

$\frac{2{y}_{>}\left(Z,T\right)}{\theta \left(T-Z\right){\text{e}}^{-i\Omega T}}=\left(\begin{array}{c}\gamma +N\\ \gamma +N\end{array}\right){\text{e}}^{+i\Omega Z}+\left(\begin{array}{c}-\gamma \left[i{G}_{0}+{G}_{2}\right]-NZ{G}_{1}\\ +N\left[i{G}_{0}-{G}_{2}\right]-\gamma Z{G}_{1}\end{array}\right)\left(Z\ge 0\right)$ (9)

The functions ${G}_{0}\left(Z,T\right)$, ${G}_{1}\left(Z,T\right)$ and ${G}_{2}\left(Z,T\right)$ are defined in Appendix A. Notice the step function $\theta \left(T-Z\right)$ which shows the correct relativistic prediction that no traveling wave arrives at position z until $ct\ge z$.

3. Intermediate-Energy Diffraction in Time

In this section the case of an incoming monochromatic beam with intermediate relativistic velocities, $\beta =0.4$, and the particle detector located at $Z=20$ ( $z={\lambda }_{c}Z=4.2×{10}^{-15}\text{ }\text{m}$, for neutrons), is investigated. We use the Dirac solution: ${\psi }_{>}\left(z,t\right)$, given in Equation (9), to calculate the probability density $\rho$ :

$\rho \left(Z,T\right)\equiv {|{\psi }_{>}|}^{2}={|{\psi }_{1}|}^{2}+{|{\psi }_{3}|}^{2}\equiv {\rho }_{1}+{\rho }_{3}.$ (10)

For a position, $Z=20$, we show in Figure 2 a plot of the probability density as a function of time. We find, as expected, damped diffraction oscillations which look similar to the ones predicted by the Schrödinger theory (see Figure 1). However, the present relativistic diffraction differs from the Schrödinger one by the presence of a small amplitude oscillation superposed on the main one. Understanding these double oscillations is important. For $\beta =0.4$, the initial condition for $\left({\psi }_{1},{\psi }_{3}\right)$ given in Equation (7) has constant amplitudes $\left(N,\gamma \right)$ which

Figure 2. Dirac transient oscillations with $Z=20$ and $v/c=0.4$.

are quite different in magnitude: ${N}^{2}\approx 1-{\beta }^{2}/4=0.96$ and ${\gamma }^{2}\approx {\beta }^{2}/2=0.04$. This explains the plots shown in Figure 2, where the amplitude oscillations of ${\rho }_{3}$ is very small in comparison with those of ${\rho }_{1}$. For low energies the relativistic density ${\rho }_{1}$ looks like the Schrödinger diffraction oscillations and ${\rho }_{3}$ is just a small oscillatory perturbation. In the plot of Figure 2 we measure the half-period of one initial oscillation, $\Delta T=74.5$ ( $\Delta t={\lambda }_{c}\Delta T/c=1.0×{10}^{-22}$ sec. for neutrons).

Extreme-Relativistic Diffraction in Time

The next question is: for 1/2-spin particles, how does diffraction in time look like for an extreme-relativistic value of $\beta$ ? In this section we show the main contribution of this paper, namely that as we gradually increase the velocity $\beta$, the Dirac diffraction changes gradually from a Fresnel pattern into quantum damped beats!

To illustrate this result we take the case of an extreme-relativistic velocity: $\beta =0.998$, and the particle detector at a fixed distance: $Z=1$, ( $z\equiv Z{\lambda }_{c}=2.1×{10}^{-16}\text{ }\text{ }\text{m}$, for neutrons). For this case we plot in Figure 3, as a function of time T, the exact Dirac density $\rho \left(Z,T\right)$ given in Equation (9).

The origin of these quantum beats is easy to understand. For an extreme-relativistic value of $\beta$ we have: First, two Dirac components $\left({\psi }_{1},{\psi }_{3}\right)$ which have similar amplitudes: ${N}^{2}=0.53$ and ${\gamma }^{2}=0.47$, (in fact, if $\beta \to 1$, then: ${N}^{2}\to {\gamma }^{2}=0.5$ ). Second, looking at the plot on Figure 3, we measure the half-periods of oscillations for ${\rho }_{1}$ and ${\rho }_{3}$, we get: $\Delta {T}_{1}=0.42$ and $\Delta {T}_{3}=0.38$, ( $\Delta {t}_{1}={\lambda }_{c}\Delta {T}_{1}/c=3×{10}^{-25}\mathrm{sec}$, for neutrons). These corresponds to angular frequencies: ${\Omega }_{1}=\pi /\Delta {T}_{1}=7.4$ and ${\Omega }_{3}=8.2={\Omega }_{1}+0.8$. So, we have two density oscillations with similar amplitudes and similar angular frequencies. When we add them together we get quantum beats! This result is well known using trigonometric functions: $\mathrm{sin}\left({A}_{1}\right)+\mathrm{sin}\left({A}_{1}+\epsilon \right)=2\mathrm{cos}\left(\epsilon /2\right)\mathrm{sin}\left({A}_{1}+\epsilon /2\right)$. Therefore, for the extreme-relativistic, $\beta =0.998$, the total probability density may be interpreted as the product of two oscillations, one with slow frequency:

Figure 3. Quantum beats for $Z=1$ and extreme-relativistic $v/c=0.998$.

${\Omega }_{slow}=\left({\Omega }_{3}-{\Omega }_{1}\right)/2=0.4$ and the other with a fast one: ${\Omega }_{fast}=\left({\Omega }_{3}+{\Omega }_{1}\right)/2=7.4+0.4$. The slow frequency oscillation modulates the amplitude of the fast one. Density quantum beats is a pure relativistic quantum result!

4. Conclusions

Exact Dirac diffraction in time solution, what for? We claim that the exact Dirac diffraction in time solution can be used as a basic building block to solve a more interesting problem: Consider a relativistic particle which, for all negative times, was inside an infinite potential well with walls at $z=\left(0,a\right)$. Suddenly, at time $t=0$, the particle is released, allowing for $t>0$ the free time-evolution of the particle. For $t>0$, what is the time-dependent density observed at a distance $z>a$ or $z<0$ ?

To solve this problem, for 1/2-spin particles, we need the time-dependent solution of a free-particle Dirac equation, with the following initial condition:

$\psi \left(z,0\right)={\left(\begin{array}{cc}N& \gamma \end{array}\right)}^{\text{T}}\left(\frac{{\text{e}}^{ikz}-{\text{e}}^{-ikz}}{2i}\right)\left[\theta \left(a-z\right)-\theta \left(-z\right)\right].$ (11)

Here k is a quantized, relativistic, wave-number. It is evident that the time-evolution of this initial condition is the superposition of four similar solutions for the Dirac’s diffraction in time solutions.

At first glance we expect, at low and high energies, a sort of Fraunhofer and Fresnel diffraction density patterns which resemble the optical diffraction of light by a narrow slit. The important fact is that, at extreme-relativistic energies, the expected density wave pattern will be a pure relativistic quantum result: outgoing density quantum beats at right and left directions! As far as we know, this has never been reported before. How useful is this model? That remains to be seen. We will get the exact analytic solution of this problem in the near future.

As for the verification of this theoretical prediction, we make clear the difficulty of doing so. For neutrons at internuclear distances, $z~{10}^{-16}\text{ }\text{ }\text{m}$, and velocities $v/c=0.998$, our Dirac shutter solution predicts beats oscillations with periods of the order $\Delta t~{10}^{-25}\mathrm{sec}$ ! As far as we know, these times order of magnitude is out of reach of the present technology.

Appendix

A1. The Dirac Solution

Using dimensionless variables: $Z\equiv z/{\lambda }_{c}$, $T\equiv ct/{\lambda }_{c}$, $\kappa \equiv k{\lambda }_{c}$ and $\Omega \equiv \omega {\lambda }_{c}/c$, the Dirac Equation (6) and the initial condition, Equation (7) become:

$\left(I\frac{\partial }{\partial T}+{\sigma }_{x}\frac{\partial }{\partial Z}+i\text{ }{\sigma }_{z}\right)\psi \left(Z,T\right)=0,$ (12)

$\psi \left(Z,0\right)\equiv {\left(\begin{array}{cc}{\psi }_{1}\left(Z,0\right)& {\psi }_{3}\left(Z,0\right)\end{array}\right)}^{\text{T}}={\left(\begin{array}{cc}N& \gamma \end{array}\right)}^{\text{T}}{\text{e}}^{i\kappa Z}\theta \left(-Z\right).$ (13)

Taking the Laplace transform in Equation (12), and denoting,

$\varphi \left(Z,s\right)\equiv \mathcal{L}\left[\psi \left(Z,T\right)\right]\equiv {\int }_{0}^{\infty }\text{ }\text{ }{\text{e}}^{-sT}\psi \left(Z,T\right)\text{d}T,$

Equation (12) becomes:

$\frac{\text{d}\varphi }{\text{d}Z}+\left({\sigma }_{x}s+{\sigma }_{y}\right)\varphi =\psi \left(Z,0\right),\text{ }\left(-\infty (14)

This is a linear matrix differential equation in the Z variable, and the solution is readily obtained. Due to the presence of the step function $\theta \left(-Z\right)$, the origin $Z=0$ is a singular point, where we demand that the function $\varphi \left(Z,s\right)$ must be continuous. This fact suggest breaking the infinite range $\left(-\infty into the left ( $Z<0$ ) and right ( $Z>0$ ) ranges.

For the left-side of the shutter, $Z<0$, we define the function ${\varphi }_{<}\left(Z,s\right)$ as the solution of the differential equation:

$\frac{\text{d}{\varphi }_{<}}{\text{d}Z}+\left({\sigma }_{x}s+{\sigma }_{y}\right){\varphi }_{<}={\left(\begin{array}{cc}N& \gamma \end{array}\right)}^{\text{T}}{\text{e}}^{i\kappa Z},\text{ }\left(Z<0\right)$ (15)

and for the right-side, $Z>0$, we define ${\varphi }_{>}\left(Z,s\right)$ as the solution of

$\frac{\text{d}{\varphi }_{>}}{\text{d}Z}+\left({\sigma }_{x}s+{\sigma }_{y}\right){\varphi }_{>}=0,\text{ }\left(Z>0\right).$ (16)

Both functions ${\varphi }_{<}$ and ${\varphi }_{>}$ must be bounded, ( ${\varphi }_{<}$ at $Z=-\infty$ ) and ( ${\varphi }_{>}$ at $Z=+\infty$ ), and be continuous at the interface, $Z=0$.

The Hermitian matrix: ${\sigma }_{x}s+{\sigma }_{y}$ has eigenvalues: ${\lambda }_{1}=\sqrt{{s}^{2}+1}=-{\lambda }_{2}$, with associated orthonormal eigenvectors:

${u}_{1}=\frac{1}{\sqrt{2}}{\left(\begin{array}{cc}1& \sqrt{\frac{s+i}{s-i}}\end{array}\right)}^{\text{T}},\text{ }{u}_{2}=\frac{1}{\sqrt{2}}{\left(\begin{array}{cc}1& -\sqrt{\frac{s+i}{s-i}}\end{array}\right)}^{\text{T}}.$ (17)

Taking into account the boundary conditions, we have the general solutions for both differential equations:

${\varphi }_{>}\left(Z,s\right)=\frac{A}{\sqrt{2}}{\left(\begin{array}{cc}1& \sqrt{\frac{s+i}{s-i}}\end{array}\right)}^{\text{T}}{\text{e}}^{-Z\sqrt{{s}^{2}+1}}\text{ }\left(Z>0\right),$ (18)

${\varphi }_{<}\left(Z,s\right)=\frac{B}{\sqrt{2}}{\left(\begin{array}{cc}1& -\sqrt{\frac{s+i}{s-i}}\end{array}\right)}^{\text{T}}{\text{e}}^{+Z\sqrt{{s}^{2}+1}}+{\left(\begin{array}{cc}N& \gamma \end{array}\right)}^{\text{T}}\frac{{\text{e}}^{i\kappa Z}}{s+i\Omega }\text{ }\left(Z<0\right).$ (19)

The constants A and B are fixed by the continuity condition: ${\varphi }_{<}\left(0,s\right)={\varphi }_{>}\left(0,s\right)$, we have:

$A=\frac{1}{\sqrt{2}}\frac{1}{s+i\Omega }\left[\gamma \sqrt{\frac{s-i}{s+i}}+N\right],\text{ }B=\frac{1}{\sqrt{2}}\frac{1}{s+i\Omega }\left[\gamma \sqrt{\frac{s-i}{s+i}}-N\right]$ (20)

Substituting (A, and B) into Equations (18) and (19) we get the solution for Dirac shutter problem in the $\left(Z,s\right)$ space.

$\text{ }{\psi }_{>}\left(Z,s\right)=\frac{{\text{e}}^{-Z\sqrt{{s}^{2}+1}}}{2\left(s+i\Omega \right)}\left(\begin{array}{c}\gamma \sqrt{\frac{s-i}{s+i}}+N\\ \gamma +N\sqrt{\frac{s+i}{s-i}}\end{array}\right)\text{ }\left(Z>0\right)$

${\psi }_{<}\left(Z,s\right)=\frac{{\text{e}}^{+Z\sqrt{{s}^{2}+1}}}{2\left(s+i\Omega \right)}\left(\begin{array}{c}\gamma \sqrt{\frac{s-i}{s+i}}-N\\ -\gamma +N\sqrt{\frac{s+i}{s-i}}\end{array}\right)+\left(\begin{array}{c}N\\ \gamma \end{array}\right)\frac{{\text{e}}^{i\kappa Z}}{s+i\Omega }\text{ }\left(Z<0\right)$ (21)

Notice the simple pole, $s=-i\Omega$, and the branch points, $s=±i$, all of them located in the imaginary axis. By the Nyquist stability criterion, the time dependent solution $\psi \left(Z,T\right)$ will be an oscillatory bounded solution.

Using the convolution theorem and Laplace Transforms Tables [14] we find the following results, valid only for $Z\ge 0$.

${\mathcal{L}}^{-1}\left[\frac{{\text{e}}^{-Z\sqrt{{s}^{2}+1}}}{s+i\Omega }\right]=\theta \left(T-Z\right)\text{ }{\text{e}}^{-i\Omega T}\left[{\text{e}}^{+i\Omega Z}-Z{G}_{1}\right]\text{ }\left(Z\ge 0\right),$ (22)

${\mathcal{L}}^{-1}\left[\sqrt{\frac{s-i}{s+i}}\frac{{\text{e}}^{-Z\sqrt{{s}^{2}+1}}}{s+i\Omega }\right]=\theta \left(T-Z\right){\text{e}}^{-i\Omega T}\left[-i{G}_{0}-{G}_{2}+{\text{e}}^{+i\Omega Z}\right]\text{ }\left(Z\ge 0\right).$ (23)

Here, to simplify the notation, we have denoted the integral-defined, complex functions, valid only for ( $0 ):

${G}_{0}\left(Z,T;\Omega \right)\equiv {\int }_{Z}^{T}\text{ }\text{ }\text{d}u\text{ }\text{ }{\text{e}}^{i\Omega u}{J}_{0}\left(\sqrt{{u}^{2}-{Z}^{2}}\right),$ (24)

${G}_{1}\left(Z,T;\Omega \right)\equiv {\int }_{Z}^{T}\text{ }\text{ }\text{d}u\text{ }\text{ }{\text{e}}^{i\Omega u}\frac{{J}_{1}\left(\sqrt{{u}^{2}-{Z}^{2}}\right)}{\sqrt{{u}^{2}-{Z}^{2}}},$ (25)

${G}_{2}\left(Z,T;\Omega \right)\equiv {\int }_{Z}^{T}\text{ }\text{ }\text{d}u\text{ }\text{ }{\text{e}}^{i\Omega u}u\frac{{J}_{1}\left(\sqrt{{u}^{2}-{Z}^{2}}\right)}{\sqrt{{u}^{2}-{Z}^{2}}}$ (26)

where ${J}_{0}\left(u\right)$ and ${J}_{1}\left(u\right)$ are Bessel functions of first kind and order 0 and 1, respectively.

Therefore, for $Z\ge 0$, the exact transmitted wave is given by:

$\frac{2{\psi }_{>}\left(Z,T\right)}{\theta \left(Z-T\right){\text{e}}^{-i\Omega T}}=\left(\begin{array}{c}N+\gamma \\ N+\gamma \end{array}\right){\text{e}}^{i\Omega Z}-\left(\begin{array}{c}\gamma \left[i{G}_{0}+{G}_{2}\right]+NZ{G}_{1}\\ \gamma Z\text{ }{G}_{1}-N\left[i{G}_{0}-{G}_{2}\right]\end{array}\right)$ (27)

In the same way we obtain, as expected, for $Z\le 0$ the incident and reflected wave:

$\begin{array}{c}\frac{2{\psi }_{<}\left(Z,T\right)}{\theta \left(Z+T\right){\text{e}}^{-i\Omega T}}=2\left(\begin{array}{c}N\\ \gamma \end{array}\right){\text{e}}^{i\kappa Z}+\left(\begin{array}{c}\gamma -N\\ N-\gamma \end{array}\right){\text{e}}^{i\Omega |Z|}\\ \text{\hspace{0.17em}}\text{ }\text{ }-\left(\begin{array}{c}\gamma \left[i{G}_{0}\left(|Z|\right)+{G}_{2}\left(|Z|\right)\right]+NZ{G}_{1}\left(|Z|\right)\\ \gamma Z{G}_{1}\left(|Z|\right)-N\left[i{G}_{0}\left(|Z|\right)-{G}_{2}\left(|Z|\right)\right]\end{array}\right)\end{array}$ (28)

We see that the 1D shutter problem is, in fact, a particular time-dependent scattering problem, and here we have the exact relativistic Dirac solution for spin-1/2 particles.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

 [1] Moshinsky, M. (1951) Boundary Conditions and Time-Dependent States. Physical Review, 84, 525. https://doi.org/10.1103/PhysRev.84.525 [2] Born, M. and Wolf, E. (1965) Principles of Optics. Pergamon Press, Oxford, 192-195. [3] Szriftigiser, P., Guery-Odelin, D., Arndt, M. and Dalibard, J. (1996) Atomic Wave Diffraction and Interference Using Temporal Slits. Physical Review Letters, 77, 4. https://doi.org/10.1103/PhysRevLett.77.4 [4] Schneble, D., Hasuo, M., Anker, T., Pfau, T. and Mlynek, J. (2003) Integrated Atom-Optical Circuit with Continuous-Wave Operation. Journal of the Optical Society of America B, 20, 648-651. https://doi.org/10.1364/JOSAB.20.000648 [5] Garca-Calderón, G., Rubio, G. and Villavicencio, J. (1999) Low-Energy Relativistic Effects and Nonlocality in Time-Dependent Tunneling. Physical Review A, 59, 1758. https://doi.org/10.1103/PhysRevA.59.1758 [6] Delgado, F., Muga, J.G., Ruschhaput, A., Garca-Calderón, G. and Villavicencio, J. (2003) Tunneling Dynamics in Relativistic and Nonrelativistic Wave Equations. Physical Review A, 68, Article ID: 032101. https://doi.org/10.1103/PhysRevA.68.032101 [7] Brouard, M. and Muga, J.G. (1996) Transient and Asymptotic Effects in Tunneling. Physical Review A, 54, 3055. https://doi.org/10.1103/PhysRevA.54.3055 [8] Garca-Calderón, G. and Rubio, G. (1997) Transient Effects and Delay Time in the Dynamics of Resonant Tunneling. Physical Review A, 55, 3361.https://doi.org/10.1103/PhysRevA.55.3361 [9] Garca-Calderón, G. and Villavicencio, J., Delgado, F. and Muga, J.G. (2002) Time Scale of Forerunners in Quantum Tunneling. Physical Review A, 66, Article ID: 042119. https://doi.org/10.1103/PhysRevA.66.042119 [10] Stevens, K.W.H. (1980) A Note on Quantum Mechanical Tunneling. European Journal of Physics, 1, 98. https://doi.org/10.1088/0143-0807/1/2/006 [11] Delgado, F., Cruz, H. and Muga, J.G. (2002) The Transient Response of a Quantum Wave to an Instantaneous Potential Step Switching. Journal of Physics A: Mathematical and General, 35, Article ID: 10377. https://doi.org/10.1088/0305-4470/35/48/311 [12] Delgado, F., Muga, J.G. and Garca-Calderón, G. (2005) Resonant Tunneling Transients and Decay for a One-Dimensional Double Barrier Potential. Journal of Applied Physics 97, Article ID: 013705. https://doi.org/10.1063/1.1826215 [13] Godoy, S. and Villa, K. (2016) A Basis for Causal Scattering Waves, Relativistic Diffraction in Time Functions. Journal of Modern Physics, 7, 1181-1191. https://doi.org/10.4236/jmp.2016.710107 [14] Abramowitz, M. and Stegun, I.A. (1965) Handbook of Mathematical Functions. Dover Publications, NY, 1019-1030.