Next Article in Journal
Estimates on Some General Classes of Holomorphic Function Spaces
Next Article in Special Issue
On the Supersymmetry of the Klein–Gordon Oscillator
Previous Article in Journal
Azulene Moiety as Electron Reservoir in Positively Charged Systems; A Short Survey
Previous Article in Special Issue
Power Law Duality in Classical and Quantum Mechanics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

What Is the Size and Shape of a Wave Packet?

by
Larry S. Schulman
Physics Department, Clarkson University, Potsdam, NY 13699-5820, USA
Symmetry 2021, 13(4), 527; https://doi.org/10.3390/sym13040527
Submission received: 17 February 2021 / Revised: 11 March 2021 / Accepted: 16 March 2021 / Published: 24 March 2021
(This article belongs to the Special Issue Symmetries in Quantum Mechanics and Statistical Physics)

Abstract

:
Under pure quantum evolution, for a wave packet that diffuses (like a Gaussian), scattering can cause localization. Other forms of the wave function, spreading more rapidly than a Gaussian, are unlikely to localize.

What is the size and shape of a wave packet? I am talking about a wave packet of a particle (atom or molecule) in a gas. Is it a plane wave that fills the container? Or is it a microscopic (perhaps 100 nm) object? I am not talking about a situation where there is a potential holding the particle in some region, like a hydrogen atom. The only things around are other particles, the same or different. (Please see Appendix A). Strangely, this subject has not attracted a lot of attention.
Our conclusion considers a number of possibilities. Should the eventual wave functions be Gaussian or Gaussian-like (to be later defined) then, yes, there is localization. If not, probably no. However, we give an argument that the wave function most likely does become Gaussian. (Generally speaking, people use Gaussians in descriptions, but this is not always warranted).
This paper has three parts: arguments for a Gaussian, localization in that case and the case(s) that it is not a Gaussian. We emphasize that the first part (eventually) treats non-Gaussians that become Gaussian after a number of collisions.
Arguments for a Gaussian. The wave function for a Gaussian would be
ψ ( r ) = 1 ( π σ 2 ) 3 / 4 exp ( r r 0 ) 2 2 σ 2 + i p 0 ( r + r 0 ) .
In one dimension, a Gaussian in momentum space is
ϕ ( p ) = σ 2 π 2 1 / 4 exp σ 2 ( p p 0 ) 2 2 2 i p x 0 .
(The following paragraph is the treatment in [1]; for more details, see that reference.) We suppose there is a particle of a different mass and that they scatter. The conservation laws (once they are far enough apart that they do not interact) are
P final = P initial and p final = p initial ,
where P and p are collective coordinates: P = p 1 + p 2 and p = ( m 2 p 1 m 1 p 2 ) / M . In these relations, p j and m j ( j = 1 , 2 ) are the momenta and masses of each of the particles, 1 μ = 1 m 1 + 1 m 2 and M = m 1 + m 2 . Furthermore, we let μ i = m i M ( i = 1 , 2 ). It turns out that, for calculation of the spread, all that matters is the real part of the ( 2 × negative) exponent. This changes from σ 1 2 p 1 2 + σ 2 2 p 2 2 = σ 1 2 ( μ 1 P + p ) + σ 2 2 ( μ 2 P p ) , following Equation (3), to
Q ( p 1 , p 2 ) p 1 2 [ σ 1 2 ( μ 1 μ 2 ) 2 + 4 μ 2 2 σ 2 2 ] + p 2 2 [ σ 2 2 ( μ 2 μ 1 ) 2 + 4 μ 1 2 σ 1 2 ] + 4 ( μ 1 μ 2 ) p 1 p 2 ( μ 1 σ 1 2 μ 2 σ 2 2 ) ,
As shown in [1], provided μ 1 and μ 2 are significantly different from one another, within a few scattering events the “cross” term in p 1 p 2 vanishes. This means that the wave functions have ratio of spread such that m σ 2 is constant and there is no cross term—nothing to decohere (and the von Neumann entropy is maximum). In three dimensions the same thing happens, but is more difficult to show [2,3,4]. That, however, is not our main point.
What happens in one dimension if the initial functions are not a Gaussian? We suppose the wave function has the form exp 2 a n x n . This form can fit various other functional forms, e.g., 1 / ( 1 + x 2 ) 2 to within 0.06 using 4 coefficients (i.e., elements of the set { a n } ). (Please see Appendix B) . Following the earlier method (based on Equation (3)), we find that 4th order terms obey
a 4 ( 1 ) = a 4 ( 1 ) ( μ 1 μ 2 ) 4 + a 4 ( 2 ) 2 4 μ 2 4 a 4 ( 2 ) = a 4 ( 2 ) ( μ 1 μ 2 ) 4 + a 4 ( 1 ) 2 4 μ 1 4 ,
with a 4 ( j ) ( j = 1 , 2 ) the coefficients of x 4 for the respective wave functions. In general, if the deviation from a Gaussian begins with a term p k 2 n there is a matrix that takes one from the values of the coefficients multiplying these terms from before scattering to after scattering. That matrix is
μ 1 μ 2 2 n 2 μ 2 2 n 2 μ 1 2 n μ 1 μ 2 2 n .
The eigenvalues of this matrix are (for all n > 1 ) below 1 for m 1 m 2 and m j 0 ( j = 1 , 2 ).
What happens if there is no other-mass particle? I do not know. I would not have thought it should make a difference, but I do not have a proof. In practice, there is almost always some impurity, but it may scatter rarely.
Thus, there is an indication that in one dimension the wave function approaches a Gaussian. In higher dimension—in particular 3—I do not have definitive results. It is still true [1] that for Gaussians the spread approaches a maximum of von Neumann entropy and (if there are two types of particles, # 1 and # 2 then) m 1 σ 1 2 = m 2 σ 2 2 .
For three dimensions, there are many ways that n th power terms can occur; for # 1 one can have anything of the form p 1 x n 1 p 1 y n 2 p 1 z n 3 with n 1 + n 2 + n 3 = n , and similarly for # 2 , leading to ( n + 2 ) ( n + 1 ) coefficients for the two of them. (Please see Appendix C). Moreover, as discussed in [1], the post-scattering values of p 1 and p 2 involve a rotation, R SO(3), not just a flip. Thus, p 1 = ( μ 1 I + μ 2 R ) p 1 + μ 1 ( I R ) p 2 (and 1 2 for p 2 ).
I have examined cubic and quartic components of the logarithm of the wave function. All the indicated operations have been carried out, the cross terms involving momenta of # 1 and of # 2 have been dropped, and a rather complex recursion for the plethora of coefficients evaluated numerically. Sample results are shown in Figure 1 and Figure 2. The “time” represents the number of scattering events, and the ordinate is the sum total of the absolute values of all cubic or quartic coefficients. Remarkably, they all tend to zero.
Evaluations beyond quartic represent a further problem in symbolic manipulation, but based on the higher-power evidence of one dimension together with the cubic and quartic evidence in three dimensions, it is reasonable to make the claim that all higher power coefficients tend to zero under decoherence.
Is this a proof that eventually everything decoheres to a Gaussian? Absolutely not, but it is an indication.
Localization with Gaussian-like behavior. Scattering can localize. This may be surprising, since some may hold that a wave function can only spread. It turns out (as shown in [2]) that scattering can act like a measurement, that scattering alone can localize. Here we extend that result.
I do not deal with the effects of temperature (cf. [5]), nor with off-diagonal elements of the density matrix [6,7]. I am concerned with pure quantum behavior. Nor do the conclusions depend on interpretation, i.e., they are independent of whether one subscribes to the Copenhagen interpretation (in its many variants), Many Worlds, or some other theory.
I will briefly review the results of [2] and then turn to the extension. (There is a slight change in notation: instead of σ 2 , we use 2 Δ 2 for convenience in matching results.) The principal consequence of [2] is that, assuming the wave function is a Gaussian, particles do not spread indefinitely. The proof in [2] is a self-consistency argument. We assume two normalized Gaussians in 3 dimensions of the form
ψ ( r 1 , r 2 , 0 ) = 1 ( 2 π Δ 2 ) 3 / 2 exp ( r 1 + r 0 ) 2 4 Δ 2 + i p 1 ( r 1 + r 0 ) ( r 2 r 0 ) 2 4 Δ 2 + i p 2 ( r 2 r 0 ) .
with parameters r 0 , p 1 , p 2 and Δ . We assume p 1 p 2 , so that at time m | r / p 1 | they scatter. Assuming no interaction, at time t this becomes
ψ ( r 1 , r 2 , t ) = Δ 4 + ( t 2 m ) 2 2 π Δ 2 3 / 2 exp ( r 1 + r 0 p 1 t m ) 2 4 Δ 2 + i t 2 m + i p 1 ( r 1 + r 0 ) i p 1 2 t 2 m × exp ( r 2 r 0 p 2 t m ) 2 4 Δ 2 + i t 2 m + i p 2 ( r 2 r 0 ) i p 2 2 t 2 m .
using center of mass coordinates, the exponent is
( R + P t 2 m ) 2 2 Δ 2 + i t 2 m ( r + 2 r 0 + p t m / 2 ) 2 2 4 Δ 2 + i t m / 2 i R P + ( r 2 r 0 ) p i P 2 4 m + p 2 m t ,
where R = 1 2 r 1 + r 2 , r = r 1 r 2 , and correspondingly for P and p . Now set t m | r / p 1 | (or just afterward). They scatter, but nothing happens to P . In the center of mass, we use the Born approximation. Up to normalization, the new wave function is δ t i V ( r ) Ψ ( R , r , t ) , where we also assumed the interaction was brief, taking place during a short time δ t . With a further condition that V ( r ) = V 0 exp r 2 4 a 2 , for convenience in integrating, we obtain
| Ψ | 2 = exp r 2 2 a 2 ( R + P t 2 m ) 2 Δ 2 + i t 2 m ( r + 2 r 0 + p t m / 2 ) 2 4 Δ 2 + i t m / 2 .
(Normalization cancels and can be ignored.) Using this wave function, one can calculate the spread in r ( Δ r ), in R ( Δ R ) and from them Δ r 1 and Δ r 2 . The latter quantities are then set equal to the original Δ r 1 and Δ r 2 (both equal to Δ ). Using as the time between collisions scattering length velocity (and mass values) we find that indeed the original Δ can be set equal to the final values of spread. To make equations simpler, we define
θ t m = m v ,
where is the scattering length and v a typical velocity. For any given gas, the range of θ is fixed.
The essential steps in the foregoing derivation are the estimation of Δ r , R and Δ R . It is these that can be generalized.
The first step is to formalize V. Instead of an exponential of range a, we take a potential that is cut off at distance a. This restricts the location of r to be within a distance a of R . In other words, Δ r r · r = 3 a 2 . This is an additional assumption and may require a larger “a” than was previously posited.
The expectation of R is simply R = P t 2 m , since it is unaffected by the interaction/scattering.
The spread in R is another matter and is the principal source of uncertainty in our calculation. For a Gaussian-like initial wave function, we can give estimates. From our previous work [1], we have
( Δ R ) 2 = R 2 R 2 = 3 4 Δ 4 + θ 2 8 Δ 2 .
Now we would like to weaken the assumptions. It turns out that the most important property (for our purposes) of the Gaussian is its diffusive behavior. This assumption about the spread in the center of mass coordinate is weakened in a specific way: in place of the 8 that appears in the denominator, we allow smaller values, γ . (In this sense, the wave function is Gaussian-like.) Specifically, 4 < γ 8 and
( Δ R ) 2 = 3 4 Δ 4 + θ 2 γ Δ 2 .
We now use Δ r 1 2 = Δ r 2 2 = Δ R 2 + 1 4 Δ r 2 to arrive at the denominators for the new wave functions of r 1 and r 2 :
( Δ r new ) 2 = ( Δ r k ) 2 = 3 a 2 + 3 4 Δ 4 + θ 2 γ Δ 2 , k = 1.2 .
Setting the old Δ equal to the new one, we arrive at a self-consistency criterion for the spread
3 Δ 2 = 3 a 2 + 3 4 Δ 4 + θ 2 γ Δ 2 .
This is a quadratic equation whose solution is
Δ 2 = a 2 γ 4 ( γ 4 ) ± θ 2 γ 4 + a 2 γ 4 ( γ 4 ) 2 .
Taking = 70 nm, v = 500 m/s, a = 1 nm and m = 29 gm / 6 × 10 23 5 × 10 23 gm and γ = 4.1 gives a value of 4.5 nm, much less than the mean free path. (Please see Appendix D). These are typical parameters for the air and the wave function is localized. Note though that (for air) there is overlap: (number density) 1 / 3 ( 0.025 × 10 27 ) 1 / 3 3 Å. Therefore, although the wave functions occupy common volumes, they mostly do not interact.
Not a Gaussian. The essential feature of non-Gaussian wave functions is that the spread in R grows. I have not examined the “boundary,” that is, the form of, or parameters in, wave functions that eventually become Gaussian and those that do not. I will though examine various wave functions that are non-Gaussian. Consider, for example, an exponential
ψ ( r , 0 ) = λ 3 8 π exp ( λ r ) .
We assume that both scattering particles have this form, with the same value of λ . The spread in the relative coordinate Δ r is still bounded by a (or 3 a ) but the center of mass coordinate grows. The spread for each, before collision, is 1 / λ . The center of mass coordinate for equal mass particles is R = ( r 1 + r 2 ) / 2 , so that the spread for of the center of mass coordinate is also 1 / λ . As a result, the center of mass can also be taken of the form Equation (17). Now apply the free propagator numerically and fit the result to an exponential. Except for particular values of λ the value of the spread has increased. This means there can be no self-consistent solution (as there was for a Gaussian).
The same happens for the wave function taken as a power.
Conclusions. The point of this is not that no wave function spreads. Rather, it places bounds on spreading for certain wave functions and makes it plausible that scattering is sometimes like a measurement, pinning a particle to a small, localized region.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In particular, we do not deal with chains of oscillators or single oscillators, which according to [8,9] (and other literature) become, by decoherence, coherent states, i.e., Gaussian.

Appendix B

The function 1 / ( 1 + x 2 ) 2 is not symmetric about x = 0 ; hence, it requires odd terms in the expansion exp 2 a n x n . This is accomplished with coefficients a 3 and a 5 that are 10 8 or less. The fourth order term is about 1/100 of the quadratic term (which is about 1.46).

Appendix C

This is the number of possibilities for both particles and is calculated by imagining the interval [ 0 , n ] as having partitions at 2 (integer) locations (and multiplying by 2). Just to be clear, the 3-dimensional tendency to have a Gaussian wave function is new material.

Appendix D

The minus sign in Equation (16) is spurious and gives an imaginary value for Δ 2 .

References

  1. Schulman, L.S. Evolution of wave-packet spread under sequential scattering of particles of unequal mass. Phys. Rev. Lett. 2004, 92, 210404. [Google Scholar] [CrossRef] [PubMed]
  2. Gaveau, B.; Schulman, L.S. Reconciling Kinetic and Quantum Theory. Found. Phys. 2020, 50, 55–60. [Google Scholar] [CrossRef] [Green Version]
  3. Schulman, L.S.; Schulman, L.J. Convergence of matrices under random conjugation: Wave packet scattering without kinematic entanglement. J. Phys. A 2006, 39, 1717–1728. [Google Scholar] [CrossRef] [Green Version]
  4. Schulman, L.S.; Schulman, L.J. Wave packet scattering without kinematic entanglement: Convergence of expectation values. IEEE Trans. Nanotech. 2004, 4, 8–13. [Google Scholar] [CrossRef]
  5. Ford, G.W.; O’Connell, R.F. Wave packet spreading: Temperature and squeezing ef- fects with applications to quantum measurement and decoherence. Am. J. Phys. 2002, 70, 319–324. [Google Scholar] [CrossRef] [Green Version]
  6. Joos, E.; Zeh, H. The Emergence of Classical Properties Through Interaction with the Environment. Z. Phys. B 1985, 59, 223–243. [Google Scholar] [CrossRef]
  7. Gaveau, B.; Schulman, L.S. Decoherence, the Density Matrix, the Thermal State and the Classical World. J. Stat. Phys. 2017, 169, 889–901. [Google Scholar] [CrossRef] [Green Version]
  8. Tegmark, M.; Shapiro, H.S. Decoherence produces coherent states: An explicit proof for harmonic chains. Phys. Rev. E 1994, 50, 2538–2547. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Zurek, W.H.; Habib, S.; Paz, J.P. Coherent States via Decoherence. Phys. Rev. Lett. 1993, 70, 1187–1190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Reduction of non-Gaussian terms with unequal mass scattering in three dimensions.Because that the problem has an intrinsic nonlinearity, the initial coefficients were varied. Usually they were taken as random between 0 and 1, but on some occasions they were taken to be 20 times that. In all cases there was convergence to zero. The programs to establish this were combinations of symbolic manipulation and numerical evaluation.
Figure 1. Reduction of non-Gaussian terms with unequal mass scattering in three dimensions.Because that the problem has an intrinsic nonlinearity, the initial coefficients were varied. Usually they were taken as random between 0 and 1, but on some occasions they were taken to be 20 times that. In all cases there was convergence to zero. The programs to establish this were combinations of symbolic manipulation and numerical evaluation.
Symmetry 13 00527 g001
Figure 2. As in Figure 1, except that this looks at quartic terms. As in the other example, all terms shrink to zero.
Figure 2. As in Figure 1, except that this looks at quartic terms. As in the other example, all terms shrink to zero.
Symmetry 13 00527 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schulman, L.S. What Is the Size and Shape of a Wave Packet? Symmetry 2021, 13, 527. https://doi.org/10.3390/sym13040527

AMA Style

Schulman LS. What Is the Size and Shape of a Wave Packet? Symmetry. 2021; 13(4):527. https://doi.org/10.3390/sym13040527

Chicago/Turabian Style

Schulman, Larry S. 2021. "What Is the Size and Shape of a Wave Packet?" Symmetry 13, no. 4: 527. https://doi.org/10.3390/sym13040527

APA Style

Schulman, L. S. (2021). What Is the Size and Shape of a Wave Packet? Symmetry, 13(4), 527. https://doi.org/10.3390/sym13040527

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop