Transient Quantum Beat Oscillations in Extreme-Relativistic Diffraction in Time

Show more

1. Introduction

Quantum beat oscillations are a common subject in Atomic and Molecular Spectroscopy. In Atomic physics, the term quantum beat refers to a superposed oscillatory behavior in the light intensity emitted by some suddenly excited atomic systems in their subsequent decay [1]. In Molecular Spectroscopy, quantum beat spectroscopy is a Doppler-free time domain method based on the creation of molecular coherences with a laser pulse and the measurement of their subsequent time evolution [2].

In this same context is the work of Villavicencio *et al*. [3], where transient phenomena of phase-modulated cutoff wave packets were explored by deriving an exact general solution to Schrödinger’s equation for finite-range potentials involving arbitrary initial quantum states. They show that the dynamical features of the probability density are governed by a virtual two-level system. They also found that for a system with a bound state the interplay between the virtual levels with the latter causes a quantum beat. They also find a regime characterized by a time-diffraction oscillation. This last result came as a big surprise; it happens that in the exact Schrödinger’s solution of the shutter problem [4], where we find diffraction in time as a consequence of the free time-evolution of an initial space-discontinuous beam of particles, the exact analytic expression for the probability density has a structure in which quantum beats do not exist.

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

In this paper we address the following question: are there quantum shutter solutions (not in the Schrödinger’s equation) where quantum beats substitute the usual Fresnel oscillations of diffraction in time? In this paper, we report the positive finding of such result in the diffraction in time of an extreme-relativistic beam of free particles.

The main contribution of this paper is to show the existence of transient quantum beats in the extreme-relativistic quantum shutter problem. After solving the Klein-Gordon equation for the 1D shutter problem, we prove that, at internuclear distances, a relativistic beam of neutral Pi-mesons, has a probability density which oscillates in time in a similar way to the Fresnel oscillations derived in the Schrödinger equation. However, for an extreme-relativistic beam, the Fresnel oscillations turn into quantum damped beats!

2. The Klein-Gordon Shutter Problem

In the relativistic shutter problem for spin-0 particles, we assume for all $t\le 0$ , that we have a discontinuous right-moving plane wave in the left side of a perfectly absorbing shutter, and zero to the right:

$\psi \left(x,t\le 0\right)=\{\begin{array}{ll}{\text{e}}^{i\left(kx-\omega t\right)}\hfill & x\le 0\hfill \\ 0\hfill & x>0\hfill \end{array}={\text{e}}^{i\left(kx-\omega t\right)}\theta \left(-x\right),$ (1)

here $\theta \left(x\right)$ denotes the step function, and the angular frequency is given by $\omega =c{\left({k}^{2}+{\mu}^{2}\right)}^{1/2}$ and $\mu \equiv mc/\hslash $ . At $t=0$ , when the shutter is suddenly opened, we have the initial conditions:

$\psi \left(x,0\right)={\text{e}}^{ikx}\theta \left(-x\right),\text{\hspace{1em}}\frac{\partial \psi \left(x,0\right)}{\partial t}=-i\omega {\text{e}}^{ikx}\theta \left(-x\right),$ (2)

which we use to solve the Klein-Gordon (K-G) equation,

$\frac{{\partial}^{2}\psi}{\partial {x}^{2}}-\frac{{\partial}^{2}\psi}{\partial {\left(ct\right)}^{2}}={\mu}^{2}\psi \mathrm{.}$ (3)

The *exact* solution of this initial-value problem is given in the appendix A, where we use dimensionless variables for “position”
$\chi $ and “time”
$\tau $ defined by

$\chi \equiv \mu x=\frac{mc}{\hslash}x,\text{\hspace{1em}}\tau \equiv \mu ct=\frac{m{c}^{2}}{\hslash}t.$ (4)

For arbitrary (
$\tau >0$ ,
$\chi >0$ ), on the right-hand side of the shutter, we have the *exact* K-G wave function solution:

$\begin{array}{c}\frac{2{\text{e}}^{i\Omega \tau}}{\theta \left(\tau -\chi \right)}\psi \left(\chi \mathrm{,}\tau \mathrm{;}\Omega \right)={\text{e}}^{i\Omega \chi}-{\displaystyle {\int}_{\chi}^{\tau}}\text{\hspace{0.05em}}\text{d}u\text{\hspace{0.17em}}{\text{e}}^{i\Omega u}\left\{\chi \frac{{J}_{1}\left[\xi \left(\chi \mathrm{,}u\right)\right]}{\xi \left(\chi \mathrm{,}u\right)}+i\text{\hspace{0.05em}}\kappa {J}_{0}\left[\xi \left(\chi \mathrm{,}u\right)\right]\right\}\\ \equiv {\text{e}}^{i\Omega \chi}+i\text{\hspace{0.05em}}\kappa \left({C}_{0}+i{S}_{0}\right)-\chi \left({C}_{1}+i{S}_{1}\right)\mathrm{.}\end{array}$ (5)

Here, we have defined both dimensionless “wave-number” $\kappa $ and “angular-frequency” $\Omega $ as:

$\xi \left(\chi ,\tau \right)\equiv \sqrt{{\tau}^{2}-{\chi}^{2}},\text{\hspace{1em}}\kappa \equiv \frac{k}{\mu}=\frac{v/c}{\sqrt{1-{\left(v/c\right)}^{2}}},\text{\hspace{1em}}\Omega \equiv \frac{\omega}{\mu c}=\frac{1}{\sqrt{1-{\left(v/c\right)}^{2}}}.$ (6)

We have also defined, the real functions: ( ${C}_{0}$ , ${S}_{0}$ , ${C}_{1}$ , ${S}_{1}$ ) given by:

${C}_{0}\left(\chi \mathrm{,}\tau \mathrm{;}\Omega \right)+i{S}_{0}\left(\chi \mathrm{,}\tau \mathrm{;}\Omega \right)\equiv \theta \left(\tau -\chi \right){\displaystyle {\int}_{\chi}^{\tau}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\text{e}}^{i\Omega u}{J}_{0}\left[\xi \left(\chi \mathrm{,}u\right)\right]\text{d}u\mathrm{,}$ (7a)

${C}_{1}\left(\chi \mathrm{,}\tau \mathrm{;}\Omega \right)+i{S}_{1}\left(\chi \mathrm{,}\tau \mathrm{;}\Omega \right)\equiv \theta \left(\tau -\chi \right){\displaystyle {\int}_{\chi}^{\tau}}\text{\hspace{0.05em}}{\text{e}}^{i\Omega u}\frac{{J}_{1}\left[\xi \left(\chi \mathrm{,}u\right)\right]}{\xi \left(\chi \mathrm{,}u\right)}\text{d}u\mathrm{.}$ (7b)

3. Time-Dependent Density

Given the K-G wave function $\psi \left(x\mathrm{,}t\right)$ , the probability density is calculated by:

$\rho \left(x\mathrm{,}t\right)=\frac{i\hslash}{2m{c}^{2}}{\psi}^{\ast}\frac{\partial \psi}{\partial t}+c\mathrm{.}c\mathrm{.}\text{\hspace{0.05em}}\mathrm{,}$ (8)

which, for dimensionless variables, turn into $\rho \left(\chi \mathrm{,}\tau \right)$ , we have:

$\rho \left(\chi \mathrm{,}\tau \right)=\frac{i}{2}{\psi}^{\ast}\left(\chi \mathrm{,}\tau \right)\frac{\partial \psi \left(\chi \mathrm{,}\tau \right)}{\partial \tau}+c\mathrm{.}c\mathrm{.}$ (9)

Notice that the initial relativistic plane-wave, $\psi \left(\chi \mathrm{,}\tau \le 0\right)=\mathrm{exp}\left[i\left(\kappa \chi -\Omega \tau \right)\right]$ , does not have unit density, in fact we have

$\rho \left(\chi \mathrm{,}\tau \le 0\right)=\Omega \text{\hspace{0.05em}}\text{\hspace{0.05em}}\mathrm{.}$ (10)

This means, as we will show, that the time-evolved density
$\rho \left(\chi >0,\tau >0\right)$ will oscillate in time around the initial value:
$\rho =\Omega $ . Using Equation (5) we get for arbitrary (
$\tau >0$ ,
$\chi >0$ ) the *exact* K-G probability density:

$\begin{array}{l}\frac{8\rho \left(\chi \mathrm{,}\tau \right)}{\Omega \theta \left(\tau -\chi \right)}\\ =\left[{\text{e}}^{i\Omega \chi}-\chi \left({C}_{1}+i{S}_{1}\right)-i\kappa \left({C}_{0}+i{S}_{0}\right)-i\frac{\chi}{\Omega}\frac{{J}_{1}\left(\xi \right)}{\xi}{\text{e}}^{i\Omega \tau}+\frac{\kappa}{\Omega}{J}_{0}\left(\xi \right){\text{e}}^{i\Omega \tau}\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\times \left[{\text{e}}^{-i\Omega \chi}-\chi \left({C}_{1}-i{S}_{1}\right)+i\kappa \left({C}_{0}-i{S}_{0}\right)\right]+c\mathrm{.}c\mathrm{.}\end{array}$ (11)

The presence of the function $\theta \left(\tau -\chi \right)$ implies, as expected, that for times $0\le t\le x/c$ the density vanishes, in full agreement with the relativity theory.

4. Relativistic Diffraction in Time for Pi-Mesons

At this point, we apply our present results to a relativistic beam of neutral Pi-mesons (
${\pi}^{0}$ ) which are 0-spin particles. We know that Pi-mesons are the carriers of nuclear forces and for this, and no other reason, we assume a particle detector fixed at a distance *x* from the shutter of about,
$x=1.7\times {10}^{-15}\text{m}$ , where nuclear forces have a maximum range. All mesons are unstable and decay in various ways. Neutral Pi-mesons have a mean-life given by:
${t}^{0}\approx {10}^{-17}\text{\hspace{0.05em}}\mathrm{sec}$ , which for a mass of,
$m=135\text{\hspace{0.17em}}\text{MeV}/{\text{c}}^{2}$ , we have,
${\tau}^{0}\approx {10}^{+7}$ . Consequently, our theoretical predictions for diffraction in time for neutral Pi-mesons will make sense only if the predicted transient behavior happens in a range of times
$\Delta \tau $ having orders of magnitude less than
${\tau}^{0}$ . Next, we will show that for all assumed conditions of position and energy, this condition is well satisfied.

Let us assume a relativistic beam of monochromatic, $\omega =c{\left({k}^{2}+{\mu}^{2}\right)}^{1/2}$ , $\mu \equiv mc/\hslash $ , neutral Pi-mesons ( $m=135\text{\hspace{0.17em}}\text{MeV}/{\text{c}}^{2}$ ) with “low” velocities, $v/c=0.5$ , and a particle detector at, $x=1.7\times {10}^{-15}\text{\hspace{0.05em}}\text{m}$ . Using the exact result in Equation (12) we show in Figure 1 the plot of the relativistic probability density versus time.

The main features of this plot are: 1) As expected, for times $0\le t\le x/c$ the density vanishes. 2) In spite of the fact that the exact relativistic density given in Equation (12) looks different from the Schrödinger one, however, for $t>x/c$ the damped oscillations we see in Figure 1 look similar to the ones predicted in the optical Fresnel diffraction by a straight edge [10]. We claim that diffraction in time is present in relativistic quantum mechanics! 3) We measure the period of one initial oscillation, $\Delta \tau =41$ , which corresponds to $\Delta t=2\times {10}^{-22}\text{sec}$ ! This range is well below the ${\pi}^{0}$ -meson’s mean-life: ${t}^{0}\approx {10}^{-17}\text{\hspace{0.05em}}\mathrm{sec}$ .

5. Damped Beats in Extreme-Relativistic Diffraction in Time

The next obvious question is: how does diffraction in time look like for extreme-

Figure 1. Relativistic diffraction in time of Pi-mesons, $v/c=0.5$ .

relativistic values of *v*/*c*? In this section, we show properties of extreme relativistic diffraction in time which, as far as we know, have never been reported before.

As we gradually increase the relativistic velocities *v*/*c*, the diffraction in time pattern gradually changes from Fresnel into damped beats! To illustrate this, in Figure 2, we have chosen the particular case of an extreme-relativistic beam of neutral Pi-mesons,
$v/c=0.99$ . Assuming the particle detector at a fixed distance,
$x=1.7\times {10}^{-15}\text{\hspace{0.05em}}\text{m}$ , we plot as a function of time the K-G shutter density
$\rho \left(\chi \mathrm{,}\tau \right)$ , given in Equation (12). We clearly see the gradual distorsion of the Fresnel oscillatios.

Finally in Figure 3 for a more energetic beam, $v/c=0.998$ , we clearly see the total transformation of the Fresnel oscillations into a perfect case of quantum beats. Here for the slow, modulating amplitude, we measure its period: $\Delta {\tau}_{slow}=6$ , which means a “slow” oscillation with a modulating period, ${T}_{slow}=3\times {10}^{-23}\mathrm{sec}$ and an angular frequency ${\omega}_{slow}=2\times {10}^{+23}{{\displaystyle \mathrm{sec}}}^{-1}$ . As for the fast oscillations we have $\Delta {\tau}_{fast}=0.4$ , implying a shorter period, ${T}_{fast}=2\times {10}^{-24}\mathrm{sec}$ and a faster angular frequency: ${\omega}_{fast}=3\times {10}^{+24}{{\displaystyle \mathrm{sec}}}^{-1}$ . All periods are shorter than the ${\pi}^{0}$ -meson’s mean-life: ${t}^{0}\approx {10}^{-17}\text{\hspace{0.05em}}\mathrm{sec}$ .

6. Conclusions

Two simultaneous properties are needed for the existence of quantum beats. The

Figure 2. Deformed Fresnel oscillations for extreme relativistic Pi-mesons.

Figure 3. Transient quantum beats for extreme relativistic Pi-mesons.

Figure 4. From bottom to top: the trigonometric, the Bessel, and the product.

first is that in the probability density, which is a quadratic expression of the wave function, we have products of two different oscillatory functions. In our case in Equation (12) we have the most important contribution of these products:

$\left[\mathrm{sin}\left(\Omega \tau \right)+\mathrm{cos}\left(\Omega \tau \right)\right]{J}_{0}\left(\sqrt{{\tau}^{2}-{\chi}^{2}}\right)\text{\hspace{1em}}\left(\tau \ge \chi \right)$ (12)

where the agular frequency $\Omega $ is given by $\Omega \equiv \frac{1}{\sqrt{1-{\left(v/c\right)}^{2}}}$ .

The second property is that only for extreme relativistic velocities the two oscillatory functions must have similar angular frequencies. Here, for $v/c=0.998$ , ( $\Omega =15$ ), the Bessel ${J}_{0}$ becomes a slow oscillating function of time and each trigonometric is a fast one. See Figure 4. The products show damped beats! However, for $v/c<0.9$ , the beats are no longer present.

In conclusion, quantum beat oscillations are the consequence, at extreme relativistic velocities, of the interference between the initial incident wave function: ${\text{e}}^{i\Omega \tau}$ , and the Green’s function of the relativistic shutter problem: ${J}_{0}\left({\tau}^{2}-{\chi}^{2}\right)$ . This a pure quantum relativistic phenomenon. We claim that our result is original and has never been reported before.

A. The Klein-Gordon Solution for the 1D Shutter Problem

Let us consider the Klein-Gordon (K-G) equation for $\psi \left(\chi \mathrm{,}\tau \right)$ defined in the infinite range $-\infty <\chi <\infty $ ,

$\frac{{\partial}^{2}\psi}{\partial {\chi}^{2}}-\frac{{\partial}^{2}\psi}{\partial {\tau}^{2}}=\psi .$ (13)

We demand bounded boundary conditions: $\underset{\left|x\right|\to \text{\hspace{0.05em}}\infty}{\mathrm{lim}}\left|\psi \left(\chi ,\tau \right)\right|<\infty $ , and the initial conditions, corresponding to the shutter problem:

$\psi \left(\chi ,0\right)={\text{e}}^{i\kappa \chi}\theta \left(-\chi \right),\text{\hspace{1em}}\frac{\partial \psi \left(\chi ,0\right)}{\partial \tau}=-i\Omega {\text{e}}^{i\kappa \chi}\theta \left(-\chi \right)$ (14)

where $\Omega \equiv \sqrt{1+{\kappa}^{2}}$ .

To solve the K-G equation, we begin taking the Laplace transform of Equation (14). Denoting,

$\varphi \left(\chi \mathrm{,}s\right)\equiv \mathcal{L}\left[\psi \left(\chi \mathrm{,}\tau \right)\right]\equiv {\displaystyle {\int}_{0}^{\infty}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\text{e}}^{-s\tau}\psi \left(\chi \mathrm{,}\tau \right)\text{d}\tau \mathrm{,}$

we find the differential equation

$\frac{{\text{d}}^{2}\varphi \left(\chi \mathrm{,}s\right)}{\text{d}{\chi}^{2}}-\left({s}^{2}+1\right)\varphi \left(\chi \mathrm{,}s\right)=-\left(s-i\Omega \right){\text{e}}^{i\kappa \chi}\theta \left(-\chi \right)$ (15)

which holds in the range $-\infty <\chi <\infty $ .

Due to the presence of the step function $\theta \left(-\chi \right)$ , the origin $\chi =0$ is a singular point where the function and its first derivative must be continuous. This fact suggests to break the infinite range into the left ( $\chi \le 0$ ) and right ( $\chi \ge 0$ ) ranges, having different differential equation for each one. For the left-hand side of the shutter, $\chi \le 0$ , we define ${\varphi}_{<}\left(\chi ,s\right)$ as the solution of the differential equation:

$\frac{{\text{d}}^{2}{\varphi}_{<}}{\text{d}{\chi}^{2}}-\left({s}^{2}+1\right){\varphi}_{<}=-\left(s-i\Omega \right){\text{e}}^{i\kappa \chi},$ (16)

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

$\frac{{\text{d}}^{2}{\varphi}_{>}}{\text{d}{\chi}^{2}}-\left({s}^{2}+1\right){\varphi}_{>}=0.$ (17)

Here both functions ${\varphi}_{<}$ and ${\varphi}_{>}$ must be bounded: ( ${\varphi}_{<}$ at $\chi \to -\infty $ ) and ( ${\varphi}_{>}$ at $\chi \to +\infty $ ). The important boundary condition is that the two functions and their corresponding first derivatives must be continuous at the interface, $\chi =0$ .

Equations (17) and (18) are ordinary second order differential equations and their solution is readily obtained. Taking into account the boundary conditions at $\pm \infty $ we have:

${\varphi}_{<}\left(\chi ,s\right)=A{\text{e}}^{\chi \sqrt{{s}^{2}+1}}+\frac{1}{s+i\Omega}{\text{e}}^{i\kappa \chi}\text{\hspace{1em}}\left(\chi \le 0\right)$ (18)

and

${\varphi}_{>}\left(\chi ,s\right)=B{\text{e}}^{-\chi \sqrt{{s}^{2}+1}}\text{\hspace{1em}}\left(\chi \ge 0\right)$ (19)

The constants *A* and *B* are fixed from the continuity conditions at the interface
$\chi =0$ :

${\varphi}_{<}\left(0,s\right)={\varphi}_{>}\left(0,s\right),\text{\hspace{1em}}\text{d}{\varphi}_{<}\left(0,s\right)/\text{d}\chi =\text{d}{\varphi}_{>}\left(0,s\right)/\text{d}\chi $

We have a set of coupled algebraic equations for the constants *A* and *B*:

$A+\frac{1}{s+i\Omega}=B$ (20)

$A\sqrt{{s}^{2}+1}+\frac{i\kappa}{s+i\Omega}=-B\sqrt{{s}^{2}+1}$ (21)

with solutions given by:

$2A=\frac{-1}{s+i\Omega}\left(1+\frac{i\kappa}{\sqrt{{s}^{2}+1}}\right)$ (22)

$2B=\frac{1}{s+i\Omega}\left(1-\frac{i\kappa}{\sqrt{{s}^{2}+1}}\right)$ (23)

Substituting Equation (23) into Equation (19) and Equation (24) into Equation (20) we get, in space (
$\chi \mathrm{,}s$ ), the* exact* solution for relativistic diffraction in time.

For ( $\chi \le 0$ ) we have the solution:

${\varphi}_{<}\left(\chi ,s\right)=\frac{-1}{2}\frac{1}{s+i\Omega}\left(1+\frac{i\kappa}{\sqrt{{s}^{2}+1}}\right){\text{e}}^{\chi \sqrt{{s}^{2}+1}}+\frac{1}{s+i\Omega}{\text{e}}^{i\kappa \chi},$ (24)

and for ( $\chi \ge 0$ ) we have:

${\varphi}_{>}\left(\chi ,s\right)=\frac{1}{2}\frac{1}{s+i\Omega}\left(1-\frac{i\kappa}{\sqrt{{s}^{2}+1}}\right){\text{e}}^{-\chi \sqrt{{s}^{2}+1}}.$ (25)

To get the time solution in the variable $\tau $ , we must invert the Laplace transforms. We readily find in Laplace Transforms Tables [11] the following results.

${\mathcal{L}}^{-1}\left[\frac{1}{s+i\Omega}\right]={\text{e}}^{-i\Omega \tau}$ (26)

${\mathcal{L}}^{-1}\left[{\text{e}}^{-\left|\chi \right|\sqrt{{s}^{2}+1}}\right]=\delta \left(\tau -\left|\chi \right|\right)-\left|\chi \right|\frac{{J}_{1}\left(\sqrt{{\tau}^{2}-{\chi}^{2}}\right)}{\sqrt{{\tau}^{2}-{\chi}^{2}}}\theta \left(\tau -\left|\chi \right|\right)$ (27)

${\mathcal{L}}^{-1}\left[\frac{{\text{e}}^{-\left|\chi \right|\sqrt{{s}^{2}+1}}}{\sqrt{{s}^{2}+1}}\right]={J}_{0}\left(\sqrt{{\tau}^{2}-{\chi}^{2}}\right)\theta \left(\tau -\left|\chi \right|\right)$ (28)

where ${J}_{0}\left(x\right)$ and ${J}_{1}\left(x\right)$ denote Bessel functions. Next, using the convolution theorem, and after some simplifications we have:

${\mathcal{L}}^{-1}\left[\frac{{\text{e}}^{-\left|\chi \right|\sqrt{{s}^{2}+1}}}{\left(s+i\Omega \right)\sqrt{{s}^{2}+1}}\right]=\theta \left(\tau -\left|\chi \right|\right){\displaystyle {\int}_{\left|\chi \right|}^{\tau}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{d}u\text{\hspace{0.17em}}{\text{e}}^{-i\Omega \left(\tau -u\right)}{J}_{0}\left(\sqrt{{u}^{2}-{\chi}^{2}}\right)$ (29)

and

${\mathcal{L}}^{-1}\left[\frac{{\text{e}}^{-\left|\chi \right|\sqrt{{s}^{2}+1}}}{s+i\Omega}\right]=\theta \left(\tau -\left|\chi \right|\right)\left[{\text{e}}^{-i\Omega \left(\tau -\left|\chi \right|\right)}-\left|\chi \right|{\displaystyle {\int}_{\left|\chi \right|}^{\tau}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{d}u\text{\hspace{0.17em}}{\text{e}}^{-i\Omega \left(\tau -u\right)}\frac{{J}_{1}\left(\sqrt{{u}^{2}-{\chi}^{2}}\right)}{\sqrt{{u}^{2}-{\chi}^{2}}}\right]$ (30)

With the help of Equations (30) and (31) we have the final solutions:

For $\chi \le 0$ we get, as expected, the incident and reflected wave:

$\begin{array}{l}{\psi}_{<}\left(\chi ,\tau \right)={\text{e}}^{i\left(\kappa \chi -\Omega \tau \right)}-\frac{1}{2}\theta \left(\tau +\chi \right)[{\text{e}}^{-i\Omega \left(\chi +\tau \right)}+\chi {\displaystyle {\int}_{-\chi}^{\tau}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{d}u\text{\hspace{0.17em}}{\text{e}}^{i\Omega \left(u-\tau \right)}\frac{{J}_{1}\left(\sqrt{{u}^{2}-{\chi}^{2}}\right)}{\sqrt{{u}^{2}-{\chi}^{2}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\begin{array}{c}\stackrel{{}^{}}{\text{\hspace{0.05em}}}\\ \text{\hspace{0.05em}}\end{array}+i\kappa {\displaystyle {\int}_{-\chi}^{\tau}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{d}u\text{\hspace{0.17em}}{\text{e}}^{i\Omega \left(u-\tau \right)}{J}_{0}\left(\sqrt{{u}^{2}-{\chi}^{2}}\right)],\end{array}$ (31)

and for $\chi \ge 0$ we get the transmitted wave:

$\begin{array}{l}{\psi}_{>}\left(\chi \mathrm{,}\tau \right)=\frac{1}{2}\theta \left(\tau -\chi \right)[{\text{e}}^{i\Omega \left(\chi -\tau \right)}-\chi {\displaystyle {\int}_{\chi}^{\tau}}\text{\hspace{0.05em}}\text{d}u\text{\hspace{0.17em}}{\text{e}}^{-i\Omega \left(\tau -u\right)}\frac{{J}_{1}\left(\sqrt{{u}^{2}-{\chi}^{2}}\right)}{\sqrt{{u}^{2}-{\chi}^{2}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\begin{array}{c}\stackrel{{}^{}}{\text{\hspace{0.05em}}}\\ \text{\hspace{0.05em}}\end{array}-i\kappa {\displaystyle {\int}_{\chi}^{\tau}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{d}u\text{\hspace{0.17em}}{\text{e}}^{-i\Omega \left(\tau -u\right)}{J}_{0}\left(\sqrt{{u}^{2}-{\chi}^{2}}\right)]\mathrm{,}\end{array}$ (32)

We see that in fact, the 1D shutter problem is nothing but a particular time-dependent scattering problem, and here we have the exact relativistic solution.

References

[1] Hegerfeldt, G.C. and Plenio, M.B. (1994) Quantum Optics: Journal of the European Optical Society Part B, 6, 15.

https://doi.org/10.1088/0954-8998/6/1/003

[2] Carter, R.E. and Huber, J.R. (2000) Chemical Society Reviews, 29, 305-314.

https://doi.org/10.1039/a900724e

[3] Villavicencio, J. and Hernández-Maldonado, A. (1996) Physical Review A, 54, 3055.

https://doi.org/10.1103/PhysRevA.101.042109

[4] Moshinsky, M. (1951) Physical Review, 84, 525.

https://doi.org/10.1103/PhysRev.84.525

[5] Schneble, D., Hasuo, M., Anker, T., Pfau, T. and Mlynek, J. (2003) Journal of the Optical Society of America B, 20, 648-651.

https://doi.org/10.1364/JOSAB.20.000648

[6] Brouard, M. and Muga, J.G. (1996) Physical Review A, 54, 3055.

https://doi.org/10.1103/PhysRevA.54.3055

[7] Delgado, F., Cruz, H. and Muga, J.G. (2002) Journal of Physics A: Mathematical and General, 35, 1037.

https://doi.org/10.1088/0305-4470/35/48/311

[8] Garca-Calderón, G., Rubio, G. and Villavicencio, J. (1999) Physical Review A, 59, 1758.

https://doi.org/10.1103/PhysRevA.59.1758

[9] Delgado, F., Muga, F., Ruschhaput, A., Garca-Calderón, G. and Villavicencio, J. (2003) Physical Review A, 68, Article ID: 032101.

https://doi.org/10.1103/PhysRevA.68.032101

[10] Born, M. and Wolf, E. (1965) Principles of Optics. Pergamon Press, Oxford.

[11] Abramowitz, M. and Stegun, I.A. (1965) Handbook of Mathematical Functions. Dover Publications, New York.