The time dependence of onedimensional quantum mechanical probability densities is presented when the potential in which a particle moves is suddenly changed, called a quench. Quantum quenches are mainly addressed, but a comparison with results for the dynamics in the framework of classical statistical mechanics is useful. Analytical results are presented when the initial and final potentials are harmonic oscillators. When the final potential vanishes, the problem reduces to the broadening of wave packets. A simple introduction to the concept of the Wigner function is presented, which allows a better understanding of the dynamics of general wave packets. It is pointed out how special the broadening of Gaussian wave packets is, the only example usually presented in quantum mechanics textbooks.
A quantum quench is a sudden change in the Hamilton operator of a system. The initial operator $ H \u0302 i$ is changed to the final one $ H \u0302 f$. The study of the dynamics after the quench is an active field of research in quantum manybody systems,^{1} as this presents an interesting example of nonequilibrium physics.
In this paper, we present one simple example of such a quench. Our discussion for a single particle in one dimension complements standard courses as only a few timedependent problems are discussed in most quantum mechanics textbooks.
As the initial state, we take the ground state $  \psi 0 \u27e9$ of the initial Hamiltonian $ H \u0302 i = T \u0302 + V \u0302 i$, where $ T \u0302 = p \u0302 2 / 2 m$ is the operator of the kinetic energy and $ V i ( x )$ is the initial attractive potential. At t = 0, the change to $ H \u0302 f = T \u0302 + V \u0302 f$, is performed. The time dependence can be described in the Schrödinger or the Heisenberg picture. Which is more appropriate depends on H_{f} and the initial state, here taken as the ground state of H_{i}. For the special case of vanishing $ V f ( x )$, the dynamics is just the free wave packet evolution of the ground state $  \psi 0 \u27e9$ of $ H \u0302 i$. The description of quenches in the framework of classical statistical mechanics^{2} is also presented. A comparison with the quantum mechanical results is presented in Secs. III and IV.
A complete analytical description of the dynamics is possible when the initial and final Hamiltonians are of the harmonic oscillator type. We also present results for linear potentials $ V f ( x ) = \u2212 F x$ in Sec. III.
In Sec. IV, we assume that $ V i ( x )$ is a square well potential and $ V f ( x )$ vanishes. We present a simple introduction to the concept of Wigner functions.^{3,4} This allows the description of the broadening of quantum mechanical wave packets similar to the classical case.
The Hamiltonian of a onedimensional harmonic oscillator reads
with
λ being the spring constant. The angular frequency of the oscillator is given by
$ \omega 0 = \lambda / m$.
The quantum mechanical description of the harmonic oscillator is especially simple using the ladder operators. One defines the lowering operator
$ a \u0302$ and its adjoint
$ a \u0302 \u2020$,
which obey the commutation relation
$ [ a \u0302 , a \u0302 \u2020 ] = 1 \u0302$. The position operator
$ x \u0302$ and the momentum operator
$ p \u0302$ read in terms of
$ a \u0302$ and
$ a \u0302 \u2020$,
The ground state
$  0 \u27e9$ is annihilated by
$ a \u0302$, i.e.,
$ a \u0302  0 \u27e9 = 0$ holds.
As $ H \u0302 i$ we use the Hamiltonian in Eq. (6) with $ \lambda \u2192 \lambda i$, and its ground state $  0 \u27e9 i$ as the initial state. Averages of operators $ A \u0302$ in this ground state we denote by $ \u27e8 A \u0302 \u27e9$. Then, $ \u27e8 x \u0302 \u27e9 = 0$ holds as well as $ \u27e8 p \u0302 \u27e9 = 0$.
As the final potential, we use
i.e., a harmonic oscillator with spring constant
λ_{f} shifted to the position
$ a F = F / \lambda f$ and angular frequency
$ \omega f = \lambda f / m$.
The equations of motion and their solutions for the operators
$ x \u0302 ( t )$ and
$ p \u0302 ( t )$ are identical in form to the classical ones for
x(
t) and
p(
t). In the quantum mechanical case, the classical initial conditions
x_{0} and
p_{0} are replaced by
$ x \u0302 = x \u0302 ( 0 )$ and
$ p \u0302 = p \u0302 ( 0 )$. With the definition
$ \Delta A \u0302 \u2261 A \u0302 \u2212 \u27e8 A \u0302 \u27e9$, one obtains
and
In order to calculate
$ \rho q m ( x , t )$, it is useful to express
$ \Delta x \u0302 ( t )$ as a linear combination of the ladder operators corresponding to the harmonic oscillator
$ H \u0302 i$,
where using Eq.
(8),
$ \alpha ( t )$ is given by
Using
$ a \u0302  0 \u27e9 = 0$ and
$ [ a \u0302 , a \u0302 \u2020 ] = 1 \u0302$, the expectation value of
$ ( \Delta x \u0302 ( t ) ) 2$ in the ground state of
$ H \u0302 i$ is given by
$  \alpha ( t )  2$,
In order to calculate
$ \rho q m ( x , t )$, we use the Baker–Hausdorffidentity.
^{5} It states that if
$ [ A \u0302 , [ A \u0302 , B \u0302 ] ] = 0 = [ B \u0302 , [ A \u0302 , B \u0302 ] ]$,
For operators
$ A \u0302$ and
$ B \u0302$ linear in the ladder operators, the requirements are fulfilled. Therefore, the expectation value in Eq.
(3) can easily be calculated. In the integrand in Eq.
(3), we write
and evaluate the expectation value of
$ e i k \Delta x \u0302 ( t )$ using Eqs.
(13) and
(16),
Because of
$ a \u0302  0 \u27e9 i = 0$ which implies
$ i \u27e8 0  a \u0302 \u2020 = 0$, the expectation value in the second equality equals 1. Putting this into the upper part of Eq.
(3), the Gaussian integration can be performed and one obtains
An analogous result is obtained for the momentum probability distribution
with the width determined by
Figure 1 shows the time dependence of $ \u27e8 x \u0302 ( t ) \u27e9$, $ \Delta x ( t ) = \u27e8 ( \Delta x \u0302 ( t ) ) 2 \u27e9$, $ \Delta p ( t ) = \u27e8 ( \Delta p \u0302 ( t ) ) 2 \u27e9$, and the uncertainty product. The initial wave packet makes an oscillatory motion around the center of the harmonic potential in H_{f}. The width $ \Delta x ( t )$ and $ \Delta p ( t )$ oscillate with double frequency and the uncertainty product $ \Delta x ( t ) \Delta p ( t )$ equals $ \u210f / 2$ at positions $ 0 , a F$, and $ 2 a F$ and is larger at the intermediate positions.
Fig. 1.
Time dependence of $ \u27e8 x \u0302 ( t ) \u27e9 / a F$ (solid line), $ \Delta x ( t ) / \Delta x ( 0 )$ (dashed line), $ \Delta p ( t ) / \Delta p ( 0 )$ (doubledashed dotted line), and the uncertainty product $ 2 \Delta x ( t ) \Delta p ( t ) / \u210f$ (dasheddotted line) as a function of $ \omega f t$ for $ \omega f / \omega i = 1 / 2$.
Fig. 1.
Time dependence of $ \u27e8 x \u0302 ( t ) \u27e9 / a F$ (solid line), $ \Delta x ( t ) / \Delta x ( 0 )$ (dashed line), $ \Delta p ( t ) / \Delta p ( 0 )$ (doubledashed dotted line), and the uncertainty product $ 2 \Delta x ( t ) \Delta p ( t ) / \u210f$ (dasheddotted line) as a function of $ \omega f t$ for $ \omega f / \omega i = 1 / 2$.
Close modal
The behavior of $ \rho q m ( x , t )$ in Eq. (19) is well known. It is usually obtained by calculating $ \u27e8 x  0 ( t ) \u27e9 i$ introducing the concepts of coherent and squeezed states^{5,6} as a typical quantum mechanical behavior. The solution presented here without introducing these states is simpler.
Surprisingly, Eq.
(19) can also be obtained purely classically. To show this, we consider the classical quench dynamics switching from
H_{i} to
H_{f} when
$ \rho c l ( x 0 , p 0 ) = \rho 0 ( x 0 ) \rho \u0303 0 ( p 0 )$ and
$ \rho 0 ( x 0 )$ and
$ \rho \u0303 0 ( p 0 )$ are both Gaussians. Then, the integration in Eq.
(5) can be performed analytically by first calculating
by Gaussian integrations. The remaining
$ k \u2212$ integration is also Gaussian and one obtains
with
For the case of the canonical ensemble discussed after Eq.
(5), one has
$ \u27e8 p 0 2 \u27e9 / \u27e8 x 0 2 \u27e9 = m 2 \omega i 2$ and the time dependence of
$ \u27e8 ( \Delta x ( t ) ) 2 \u27e9 c l$, apart from the prefactor is the same as in Eq.
(15). For the special temperature choice,
$ k B T = \u210f \omega i / 2$, also the prefactor is the same. This corresponds to the initial condition
where
$ \psi 0 ( x 0 ) = \u27e8 x 0  \psi 0 \u27e9$ is the real wave function for the Gaussian initial state
$  \psi 0 \u27e9 =  0 \u27e9 i$ considered in this section and
$ \psi \u0303 0 ( p 0 ) = \u27e8 p 0  \psi 0 \u27e9$ is the corresponding real Gaussian momentum amplitude.
It is left as an exercise to repeat the calculation of $ \rho q m ( x , t )$ for $  1 \u27e9 i = a \u0302 \u2020  0 \u27e9 i$ and $  2 \u27e9 i = ( a \u0302 \u2020 ) 2  0 \u27e9 i / 2$ as the initial state by expanding the exponential functions in the expectation value in the second equality in Eq. (18). The results for the special case $ H \u0302 f = T \u0302$ can be found in Ref. 7. There the free time evolution after the quench for the initial state $ c p \u0302 2  0 \u27e9 i$ is also discussed. It is a special linear combination of $  0 \u27e9 i$ and $  2 \u27e9 i$. Results for the general linear combination of these two states are presented in Ref. 8.
We next discuss the case
$ \lambda f = 0$ and a finite value of
F, i.e., the linear potential
$ V f ( x ) = \u2212 F x$. Performing the limit
$ \omega f \u2192 0$ in Eq.
(10) leads to
independent of the value of
F which only appears in the expression for
$ \u27e8 x \u0302 ( t ) \u27e9$ in Eq.
(10), given by
The result for
$  \u27e8 x  \psi ( t ) \u27e9  2$ for the linear potential is given by Eq.
(19) with the results of Eqs.
(26) and
(27) inserted. The time dependence of the broadening is identical to the free particle case.
This holds for linear potentials for arbitrary initial states
$  \psi ( 0 ) \u27e9$. The solution of the Heisenberg equation of motion is given by
Putting this into the second line of Eq.
(3) yields
This is in accord with the (weak) equivalence principle,
^{9,10} which states that all laws of a freely falling particle are the same as in an unaccelerated reference frame. The presented proof for the broadening of wave packets using the Heisenberg picture by directly addressing the measurable probability density
$  \u27e8 x  \psi ( t ) \u27e9  2$ is much simpler than using the Schrödinger picture and calculating
$ \u27e8 x  \psi ( t ) \u27e9$ first.
^{9,10}
In quantum mechanics textbooks, the factor multiplying t^{2} for the free case F = 0 in Eq. (26) is usually expressed differently as the broadening of the Gaussian wave packet is treated before the harmonic oscillator. Instead of $ \omega i 2$, the factor is written ( $ \u210f / ( 2 m \u27e8 x \u0302 2 \u27e9 ) )$^{2}.
For the case $ \lambda f = 0$ and F = 0, the free broadening of a Gaussian wave packet using operator manipulations was presented in this journal recently.^{11}
In Sec. III, we used the Heisenberg picture to calculate the time dependence of the probability density ρ_{qm} after the quench. Here, we use the Schrödinger picture to obtain $ \rho q m ( x , t )$ for the quench in which $ V f ( x ) \u2261 0$.
We begin with the usual approach to obtain $ \rho q m ( x , t )$ by calculating $ \u27e8 x  \psi ( t ) \u27e9$ and taking its absolute square. In the second part of this section, we discuss the additional insight one can obtain by using Eq. (2) instead.
In the first approach, one uses the eigenstates of
$ H \u0302 f$. For the case
$ H \u0302 f = T \u0302$, those are given by the momentum states
$  p \u27e9$. Inserting the unit operator expressed in terms of the momentum states and using
$ \u27e8 x  p \u27e9 = e ipx / \u210f / 2 \pi \u210f$ yields
with
$ \u03f5 p = p 2 / 2 m$ and the momentum representation
of the arbitrary initial state
$  \psi \u27e9$.
If the wave function $ \psi ( x ) = \u27e8 x  \psi \u27e9$ is Gaussian, the same holds for $ \psi \u0303 ( p )$ and the integration in Eq. (30) can be performed analytically. This is presented in almost all quantum mechanics textbooks. For generic $ \psi ( x )$, the integration in Eq. (30) has to be performed numerically.
In this section, we take as the initial state
$  \psi 0 \u27e9$ the ground state of an attractive square well potential
$ V i ( x ) = V i \Theta ( a \u2212  x  )$, with Θ being the step function and
$ V i < 0$. Its wave function is given by
with
$ k 0 = 2 m ( E 0 \u2212 V i ) / \u210f$ and the ground state energy
E_{0}.
In order to have continuous first derivatives at
$ x = \xb1 a$ has to hold with
$ tan ( k 0 a ) > 0$. The normalization constant is given by
$ c 0 = [ 1 / \kappa + a ( 1 + \kappa 2 / k 0 2 ) + \kappa / k 0 2 ] \u2212 1 / 2$. The dependence of
E_{0} on
V_{i} is irrelevant in the following.
There are two interesting limiting cases:

The limit $ a \u2192 0 , V i \u2192 \u2212 \u221e$ with arbitrary $ \kappa > 0$, i.e., an attractive delta potential with the ground state wave function $ \psi 0 ( x ) = \kappa e \u2212 \kappa  x $.

The infinitely deep potential: For $ a k 0 \u2192 \pi / 2$, the ratio $ \kappa / k 0$ tends to infinity and one obtains $ \psi 0 ( x ) = cos ( k 0 x ) \Theta ( a \u2212  x  ) / a$.
The momentum representation of
$  \psi 0 \u27e9$ can be calculated analytically. Using Eqs.
(31) and
(32), one obtains with
$ p = \u210f k$ Inserting this result into the integral in Eq.
(30), it can be calculated numerically to obtain
$ \u27e8 x  \psi ( t ) \u27e9$. Its absolute value squared is shown in Fig.
2 for
$ k 0 a = \pi / 2$, i.e., the ground state of the infinitely deep well, for four different times. It shows that the probability density to find the particle at the origin is larger for
$ t = 0.14 t 0$ than at the initial time
t = 0, where
$ t 0 = m a 2 / \u210f$. For
$ t = 0.07 t 0$, it has a minimum at the origin.
Fig. 2.
Probability density ρ_{qm} times a as a function of x/a for four different times: T = 0 (solid line), $ t = 0.07 t 0$ (dasheddotted line), $ t = 0.14 t 0$ (dashed line), and $ t = 0.28 t 0$ (doubledashed dotted line). Note that $ a \rho q m ( 0 , 0.14 )$ is larger than $ a \rho q m ( 0 , 0 )$.
Fig. 2.
Probability density ρ_{qm} times a as a function of x/a for four different times: T = 0 (solid line), $ t = 0.07 t 0$ (dasheddotted line), $ t = 0.14 t 0$ (dashed line), and $ t = 0.28 t 0$ (doubledashed dotted line). Note that $ a \rho q m ( 0 , 0.14 )$ is larger than $ a \rho q m ( 0 , 0 )$.
Close modal
To elucidate this surprising effect in more detail, $ a \rho q m$ is shown for x = 0 as a function of $ t / t 0$ for three different values of $ k 0 a$ in Fig. 3. The short time oscillatory behavior is more pronounced when the well is deeper.
Fig. 3.
Probability density ρ_{qm} times a for x = 0 as a function of $ t / t 0$ for three different values of $ k 0 a$: $ k 0 a = \pi / 2$, which corresponds to the infinitely deep well (solid line), $ k 0 a = \pi / 2.5$ (dashed line), and $ k 0 a = \pi / 3$ (dasheddotted line).
Fig. 3.
Probability density ρ_{qm} times a for x = 0 as a function of $ t / t 0$ for three different values of $ k 0 a$: $ k 0 a = \pi / 2$, which corresponds to the infinitely deep well (solid line), $ k 0 a = \pi / 2.5$ (dashed line), and $ k 0 a = \pi / 3$ (dasheddotted line).
Close modal
The fact that the probability to find the particle at the origin for a
t > 0 is larger than in the initial state is a purely quantum mechanical effect. It is easy to see that this cannot happen in the classical case when the probability to find the particle at
t = 0 has its maximum at
$ x 0 = 0$. For a free particle, the trajectory is given by
$ x x 0 p 0 ( t ) = x 0 + p 0 t / m$ and using Eq.
(4), one obtains
In order to be as close as possible to the quantum mechanical case for a general initial state
$  \psi \u27e9$, the phase space density
$ \rho 0 ( x 0 , p 0 )$ should yield the quantum mechanical space probability
$  \u27e8 x 0  \psi \u27e9  2$ by integration over
p_{0} and the momentum probability
$  \u27e8 p 0  \psi \u27e9  2$ by integration over
x_{0}. By choosing
this is obviously fulfilled. We call this the “classical approximation” in quotation marks as the dynamics using Eq.
(35) is classical, but this initial condition involves quantum mechanical probability densities.
For $  \psi \u27e9$, the ground state of the square well potential $  \u27e8 x 0  \psi \u27e9  2$ has its maximum at $ x 0 = 0$. For t > 0, one has $ \rho 0 ( \u2212 p 0 t / m , p 0 ) < \rho 0 ( 0 , p 0 )$, which using Eq. (35) implies $ \rho c l ( 0 , t ) < \rho c l ( 0 , 0 )$.
It turns out that the “classical approximation” works rather well for larger times. This is shown in Fig. 4 for $ k 0 a = \pi / 2$.
Fig. 4.
Probability density ρ_{qm} times a for $ k 0 a = \pi / 2$ as a function of $ t / t 0$ for three different values of x/a (solid lines), compared to the “classical approximation” described in the text (dashed lines).
Fig. 4.
Probability density ρ_{qm} times a for $ k 0 a = \pi / 2$ as a function of $ t / t 0$ for three different values of x/a (solid lines), compared to the “classical approximation” described in the text (dashed lines).
Close modal
For a better understanding of this, we show that the “classical approximation” gives the exact result for the time dependence of the quantum mechanical width of the wave packet. With the initial condition in Eq.
(36), one obtains
$ \u27e8 ( x ( 0 ) ) n \u27e9 c l = \u27e8 ( x \u0302 ( 0 ) ) n \u27e9 q m$ and
$ \u27e8 ( p ( 0 ) ) n \u27e9 c l = \u27e8 ( p \u0302 ( 0 ) ) n \u27e9 q m$. As we consider initial states
$  \psi \u27e9$ with even wave functions
$ \u27e8 x ( 0 ) \u27e9 c l = 0$ and
$ \u27e8 ( x ( 0 ) p ( 0 ) \u27e9 c l = 0$ holds. This leads to
To calculate the expectation value
$ \u27e8 x \u0302 p \u0302 \u27e9 + \u27e8 p \u0302 x \u0302 \u27e9$ for a square integrable wave function
$ \psi ( x )$, one can use that the integral over the derivative of
$ x  \psi ( x )  2$ vanishes as
$ x  \psi ( x )  2 \u2192 0 for x \u2192 \xb1 \u221e$. For
$ \psi ( x ) = f ( x ) e i \varphi ( x )$ with
f and
$\varphi $ real functions, one obtains (exercise)
Therefore, the sum vanishes for real
$ \psi ( x )$. This leads to
In order to discuss the long time behavior of
$ \rho c l ( x , t )$, one can alternatively perform the
p_{0} integration in Eq.
(4) first or change the integration variable in Eq.
(35), both leading to
If
$  \psi ( x 0 )  2$ in Eq.
(36) decays fast enough like for the square well potential Eq.
(32), the term
$ m x 0 / t$ in the integrand can be neglected in the long time limit leading to
for finite
x/
t.
In Fig. 5, we show numerical results for the quantum mechanical long time behavior using Eqs. (30) and (34). The function $ ( t / t 0 ) a \rho q m$ is shown for various values of $ t / t 0$ as a function of $ ( x / a ) / ( t / t 0 )$. It quickly approaches $ \u210f  \psi \u0303 ( ( x / a ) / ( t / t 0 ) )  2 / a$, in the agreement with Eq. (41).
Fig. 5.
Long time behavior of $ ( t / t 0 ) a \rho q m$ for the infinitely deep square well as a function of $ ( x / a ) / ( t / t 0 )$ for different values of $ t / t 0$. It quickly approaches $ \u210f  \psi \u0303 ( ( x / a ) / ( t / t 0 ) )  2 / a$ (long dashed curve).
Fig. 5.
Long time behavior of $ ( t / t 0 ) a \rho q m$ for the infinitely deep square well as a function of $ ( x / a ) / ( t / t 0 )$ for different values of $ t / t 0$. It quickly approaches $ \u210f  \psi \u0303 ( ( x / a ) / ( t / t 0 ) )  2 / a$ (long dashed curve).
Close modal
Before we present an explanation why the quantum mechanical and the “classical” long time behavior agree we switch to Eq. (2) to calculate $ \rho q m ( x , t )$ without computing $ \u27e8 x  \psi ( t ) \u27e9$ first.
For the later treatment of the free particle time dependence, it turns out to be useful to factorize
$ e i k x \u0302 = e i k x \u0302 / 2 e i k x \u0302 / 2$ and insert the unit operator expressed in terms of the momentum states in between. Using
$ e i k x \u0302 / 2  p \u27e9 =  p + \u210f k / 2 \u27e9$ as can be proven by multiplying both sides by
$ \u27e8 x $, this leads to
Putting this into Eq.
(2) and changing the integration variable
$ k = p \u0303 / \u210f$, one obtains
with
This is a real function as seen by changing the integration variable
$ p \u0303 \u2192 \u2212 p \u0303$. Integration over
x gives a factor
$ 2 \pi \delta ( p \u0303 / \u210f )$ leading to
It is left as an exercise to show that after inserting unit operators in terms of position states in the integral in Eq.
(44), the function
$ W ( x , p , t )$ can also be written as
In 1932, Wigner
^{3} introduced the quantum phasespace distribution
$ W ( x , p , t )$ for arbitrary time dependent states
$  \psi ( t ) \u27e9$, which has similarities to the classical phasespace distribution,
^{2} but allows to obtain exact quantum mechanical results, as shown above. The higher dimensional generalization of the Wigner function is widely used as a tool in various areas, e.g., quantum optics.
^{12,13} However, it is introduced only in a few quantum mechanics textbooks, e.g., Refs.
5 and
6. A “pedestrian” introduction has been published in this journal.
^{4}
For $ W ( x , p , 0 )$ to be a probability density, i.e., it is positive everywhere, $ \psi ( x )$ has to be the exponential of a quadratic polynomial.^{14} A wellknown example is a realvalued Gaussian. The integration in Eq. (46) can then be easily performed and one obtains $ W ( x , p , 0 ) = \psi ( x ) 2 \psi \u0303 ( p ) 2$.
The form of the Wigner function in Eq.
(44) is well suited to determine its time dependence for free particles discussed in this section. The momentum states are the eigenstates for
$ H \u0302 f = T \u0302$. With
$ \u03f5 p + p \u0303 / 2 \u2212 \u03f5 p \u2212 p \u0303 / 2 = p \u0303 p / m$, one obtains
If this is inserted into Eq.
(43), it leads to
This result for
$ \rho q m ( x , t )$ has the form as in the classical case in Eq.
(35), with
$ W ( x , p , 0 )$ instead of
$ \rho 0 ( x , p )$. Therefore, the arguments leading to Eq.
(41) can be generalized to obtain
if
$ W ( x , p , 0 )$ has a welllocalized
xdependence, which is the case for the ground state of a deep square well potential.
The long time behavior of the broadening of wave packets has been discussed previously in two papers in this journal.^{7,15} In Ref. 7, the result Eq. (49) is presented using the free particle propagator, shortly discussed in the Appendix. The statement in Ref. 15 that all wave packets become “approximately Gaussian” after a long enough time is critically discussed.
The Wigner function for the ground state of the infinitely deep square well can be obtained by simple integration using Eq.
(46),
^{16} It is shown in Fig.
6 as a function of momentum for two different values of
x. One can see that
W is negative in some intervals. This is different for
$  \psi ( 0 )  2  \psi \u0303 ( p )  2$ shown as the dashed curve.
Fig. 6.
Wigner function $ \u210f W ( x , \u210f k , 0 )$ as a function of ak for x = 0 and $ x = a / 2$. The dashed curve is the positive function $ \u210f  \psi ( 0 )  2  \psi \u0303 ( \u210f k )  2$.
Fig. 6.
Wigner function $ \u210f W ( x , \u210f k , 0 )$ as a function of ak for x = 0 and $ x = a / 2$. The dashed curve is the positive function $ \u210f  \psi ( 0 )  2  \psi \u0303 ( \u210f k )  2$.
Close modal
The oscillatory behavior of $ \rho q m ( 0 , t )$ shown in Fig. 3 can be understood without performing a numerical integration by plotting the integrand in Eq. (48) $ W ( \u2212 \u210f k t / m , \u210f k , 0 ) = W ( k a 2 t / t 0 , \u210f k , 0 )$ (symmetric in k). In Fig. 7, this is shown for the times of the deepest minimum ( $ t / t 0 \u2248 0.071$) and highest maximum ( $ t / t 0 \u2248 0.128$) of the full line in Fig. 3. One clearly sees that the integral over the full curve yields a result larger than that of the dotted curve.
Fig. 7.
Integrand $ W ( \u2212 \u210f k t / m , \u210f k , 0 )$ for the calculation of $ \rho q m ( 0 , t )$ using Eq. (48) for $ t = 0.14 t 0$ (full curve) and $ t = 0.07 t 0$ (dashed curve).
Fig. 7.
Integrand $ W ( \u2212 \u210f k t / m , \u210f k , 0 )$ for the calculation of $ \rho q m ( 0 , t )$ using Eq. (48) for $ t = 0.14 t 0$ (full curve) and $ t = 0.07 t 0$ (dashed curve).
Close modal
What we called the “classical approximation” can be given a different interpretation. Comparing Eqs.
(35) and
(48), it can be viewed as the following approximation for the Wigner function:
i.e., a factorized probability distribution. As discussed above, this holds exactly only for a Gaussian wave packet.
It should be mentioned that various other wave packets have been studied which are not ground state wave functions of an initial potential, e.g., a rectangular initial wave function,^{7,17} which at short times leads to a strongly oscillatory behavior. In Ref. 17, this is discussed also using the Wigner function.
The broadening of a free Gaussian wave packet is one of the few timedependent problems treated in courses on quantum mechanics. The wave function $ \u27e8 x  \psi ( t ) \u27e9$ can be calculated analytically to obtain the probability density $  \u27e8 x  \psi ( t ) \u27e9  2$. It is little known that the time dependence of other wave packets can differ considerably from the smooth Gaussian broadening.
It is important to include the treatment of timedependent Hamiltonians when teaching quantum mechanics. A simple case is the sudden change from $ H \u0302 i$ to $ H \u0302 f$ treated in this paper.
It turned out that directly addressing $  \u27e8 x  \psi ( t ) \u27e9  2$ using Eq. (3) allowed a simple way to obtain results for harmonic oscillator systems. It is little known that the oscillatory behavior in Fig. 1 can also happen in the framework of classical statistical mechanics. For the free time evolution after the quench directly addressing $  \u27e8 x  \psi ( t ) \u27e9  2$ provides additional insight, as the introduction of the concept of the Wigner function is straightforward. This function allows an understanding of the oscillatory behavior of the broadening of the initial ground state of an infinitely deep well without performing a numerical integration. This is an argument for presenting the concept of the Wigner function in quantum mechanics courses.
© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Published open access through an agreement with
University of Göttingen Faculty of Physics