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
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
In one dimension, a Gaussian in momentum space is
(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
where
P and
p are collective coordinates:
and
. In these relations,
and
(
) are the momenta and masses of each of the particles,
and
. Furthermore, we let
(
). It turns out that, for calculation of the spread, all that matters is the real part of the (
negative) exponent. This changes from
, following Equation (
3), to
As shown in [
1], provided
and
are significantly different from one another, within a few scattering events the “cross” term in
vanishes. This means that the wave functions have ratio of spread such that
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
. This form can fit various other functional forms, e.g.,
to within 0.06 using 4 coefficients (i.e., elements of the set
). (Please see
Appendix B) . Following the earlier method (based on Equation (
3)), we find that 4th order terms obey
with
(
) the coefficients of
for the respective wave functions. In general, if the deviation from a Gaussian begins with a term
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
The eigenvalues of this matrix are (for all
) below 1 for
and
(
).
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,
and
then)
.
For three dimensions, there are many ways that
power terms can occur; for
one can have anything of the form
with
, and similarly for
, leading to
coefficients for the two of them. (Please see
Appendix C). Moreover, as discussed in [
1], the post-scattering values of
and
involve a rotation,
SO(3), not just a flip. Thus,
(and
for
).
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
and of
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
, we use
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
with parameters
,
,
and
. We assume
, so that at time
they scatter. Assuming no interaction, at time
t this becomes
using center of mass coordinates, the exponent is
where
,
, and correspondingly for
and
. Now set
(or just afterward). They scatter, but nothing happens to
. In the center of mass, we use the Born approximation. Up to normalization, the new wave function is
, where we also assumed the interaction was brief, taking place during a short time
. With a further condition that
, for convenience in integrating, we obtain
(Normalization cancels and can be ignored.) Using this wave function, one can calculate the spread in
(
), in
(
) and from them
and
. The latter quantities are then set equal to the original
and
(both equal to
). Using as the time between collisions
(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
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 , and . 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 to be within a distance a of . In other words, . This is an additional assumption and may require a larger “a” than was previously posited.
The expectation of is simply , since it is unaffected by the interaction/scattering.
The
spread in
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
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,
and
We now use
to arrive at the denominators for the new wave functions of
and
:
Setting the old
equal to the new one, we arrive at a self-consistency criterion for the spread
This is a quadratic equation whose solution is
Taking
nm,
m/s,
nm and
gm
gm and
gives a value of
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)
Å. 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
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
We assume that both scattering particles have this form, with the same value of
. The spread in the relative coordinate
is still bounded by
a (or
) but the center of mass coordinate grows. The spread for each, before collision, is
. The center of mass coordinate for equal mass particles is
, so that the spread for of the center of mass coordinate is also
. 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.