Next Article in Journal
Complex Networks and Interacting Particle Systems
Next Article in Special Issue
Large Coherent States Formed from Disordered k-Regular Random Graphs
Previous Article in Journal
An Entropy-Balanced Orthogonal Learning Bamboo Forest Growth Optimization Algorithm with Quasi-Affine Transformation Evolutionary and Its Application in Capacitated Vehicle Routing Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Predictions for Norm-Conserving Spontaneous Collapse

1
Department of Physics and Astronomy, University of Pittsburgh, 100 Allen Hall, Pittsburgh, PA 15260, USA
2
Department of Mathematics, University of Pittsburgh, 301 Thackeray Hall, Pittsburgh, PA 15260, USA
*
Author to whom correspondence should be addressed.
Entropy 2023, 25(11), 1489; https://doi.org/10.3390/e25111489
Submission received: 25 September 2023 / Revised: 16 October 2023 / Accepted: 17 October 2023 / Published: 27 October 2023

Abstract

:
Previous work has shown that nonlocal collapse in quantum mechanics can be described by a deterministic, non-unitary operator added to the standard Schrödinger equation. In terms of key aspects, this term differs from prior work on spontaneous collapse. In this paper, we discuss the possible predictions of this model that can be tested by experiments. This class of collapse model does not intrinsically imply unique experimental predictions, but it allows for the possibility.

1. Introduction

Spontaneous collapse theories of quantum mechanics have the appeal that they allow the experimental results of quantum mechanics to be accounted for within a mathematical formalism that needs no separate category of “measurement”, which is left off the books. It has also been argued that a consistent approach to gravitation requires some form of spontaneous collapse [1,2,3]. There have been various approaches to spontaneous collapse, including early work by Pearle [4] and Gisin [5,6], the Ghirardi-Rimini-Weber (GRW) model [7], a formalism that puts the GRW model into a relativistic framework [8], and models of continuous spontaneous localization (CSL) [9,10,11,12,13,14]. A common characteristic of these approaches is to argue that the equations of quantum mechanics must be altered to include new, non-unitary terms. This means that they are, at least in principle, experimentally falsifiable, and conversely, there can be predictions for experiments that give evidence in support of them.
The above models involve, in one way or another, a universal source of random fluctuations that interacts with all particles, even those in the vacuum of free space. This randomness is expected to lead to both the washing out of interference and other decoherence effects [15,16,17,18,19] and other ways of detecting noise [20,21,22,23]. So far, experiments have not detected this type of noise directly. In general, these models imply the nonconservation of energy that can be detected as heat because random spatial translations of the particles amount to a heat source [7]. This prediction has recently been tested with null results [24]. An alteration of the GRW model may allow for the dissipation of this extra heat [25].
Recent work [26,27] has presented a new approach inspired by these spontaneous collapse models but without an extrinsic noise term that would wash out interference in a vacuum. This model proposes a specific new term to be added to the Schrödinger equation, which gives the Born rule for measurements but also, as discussed below, conserves energy and the total particle number and conserves the norm of the full many-body wave function. The overall approach is entirely in the context of many-body quantum field theory without committing to a definite ontology of “particles”. All that is needed is to define the localized eigenstates of fermion fields.
Since total energy is conserved in this model, the question remains whether it gives any experimentally testable predictions. In the following, we show that this class of model may allow for novel experimental predictions but need not. In the absence of any unique experimental predictions, this approach would still have the appeal of internal self-consistency and agreement with known experiments.

2. Properties of the Non-Unitary Term

The non-unitary term added to the Hamitonian in the Schrödinger equation in the model derived in Ref. [27] is
V ^ = n i ω R , n N ^ n ( 1 N ^ n ) ( 1 N ^ n ) N ^ n = n i ω R , n N ^ n N ^ n ,
where N ^ n is the number operator acting on a state, n, which is a single-particle state within a general fermionic many-particle state. As shown in Ref. [27], this operator gives a Rabi rotation on the Bloch sphere between two states, namely a fermion that does not exist in state n (state “0”) and a fermion existing in state n (state “1”), and the scalar factor ω R , n is the effective Rabi frequency. This factor is assumed to fluctuate randomly, dependent on local environmental fluctuations, as discussed further below. This model is taken to apply only to fermions, which naturally have two states corresponding to occupation 0 and 1, and does not apply to bosons (One might conjecture an analogous term for bosons but of the opposite nature, of the form α ^ n α ^ n , where α ^ n = a n + a n is the boson wave amplitude, which would favor coherent states over Fock states and give the spontaneous appearance of classical waves. For another approach, see Ref. [28]).
Because this operator gives pure Rabi rotations, it conserves the norm of any many-body state. We can see this in the following calculation. Suppose that the total Hermitian Hamiltonian is given by H = H 0 + V ^ I , where H 0 defines the eigenstates of the system, and V ^ I is the full set of standard Hermitian interactions. We write the time evolution of the norm of the state using the S-matrix expansion for the Hermitian part of the Hamiltonian (see, e.g., Ref. [29], Section 8.1):
d ψ | ψ = ψ t | ψ t ψ 0 | ψ 0 = ψ 0 | e i H 0 t / 1 + i 0 t ( V ^ I ( t ) + V ^ ( t ) ) d t 1 2 0 t ( V ^ I ( t ) + V ^ ( t ) ) d t 0 t ( V ^ I ( t ) + V ^ ( t ) ) d t + × 1 i 0 t ( V ^ I ( t ) + V ^ ( t ) ) d t 1 2 0 t ( V ^ I ( t ) + V ^ ( t ) ) d t 0 t ( V ^ I ( t ) + V ^ ( t ) ) d t + e i H 0 t / | ψ 0 ψ 0 | ψ 0 ,
where for any operator, we write, O ^ , O ^ ( t ) e i H 0 t / O ^ e i H 0 t / .
The set of all terms with just Hermitian V ^ I ( t ) operators must give unitary evolution if the entire wave function is taken into account. However, If the wave function is taken to refer only to a restricted subsystem, then the second-order expansion gives dissipative terms (as shown in Ref. [29], Section 8.1), which corresponds to the decay of the norm of the wave function; if the same second-order terms are used to compute the time-evolution of the density matrix, they give rise to the Lindbladian operator (see, e.g., Ref. [30], Section 19.3).
Moving on to those terms with the products of V ^ I ( t ) and V ^ ( t ) , we note that
V ^ ( t ) = n N ^ n e i H 0 t / N ^ n e i H 0 t / = V ^ ( 0 ) ,
since N ^ n commutes with H 0 . Therefore, integrals over V ^ ( t ) will be proportional to t, and terms with the products of integrals of V ^ I ( t ) and V ^ ( t ) will be proportional to t 2 , and therefore, will vanish for a small t, which is the differential limit taken here.
Finally, we turn to the terms with only a non-Hermitian V ^ ( t ) . The products of integrals of these will be proportional to t 2 or higher, which will vanish according to the same argument given above. We are left with just the linear terms:
ψ t | ψ t ψ 0 | ψ 0 = i t ψ 0 | V ^ V ^ | ψ 0 = n 2 ω R , n t ψ 0 | N ^ n N ^ n | ψ 0 = n 2 ω R , n t N ^ n N ^ n = 0 .
We can, therefore, call the term (1) “quasi-unitary”—although it is not Hermitian, it is strictly norm-conserving, and if ω R , n fluctuates randomly as positive and negative, then the time average of the adjoint of the operator is the same as the time average of the operator.
The form of (1) is similar to the terms in various other proposals in the literature, e.g., Refs. [5,11,12]. It is often assumed, based on the arguments discussed further in Appendix A, that a term like this must give a mechanism for superluminal communication. The general form of the argument is (1) to construct a density matrix for the state of interest and (2) then average this over many instances to get the density matrix for a mixed state that depends on a specific choice of measurement of one degree of freedom of a nonlocally entangled state, (3) act on this density matrix for with the non-Hermitian operator, and (4) from this, show that changes in the measurement can be detected nonlocally by changes in the density matrix.
This argument can only work, however, if the action of the non-Hermitian operator is the same for every instance in the average, which is carried out for the density matrix of a mixed state. If the non-Hermitian operator acts differently every time, namely with random fluctuations that lead to the Born rule, as discussed in Section 3, then its action is exactly the same as a standard quantum measurement, which cannot give superluminal communication. To put it another way, the calculation for the average evolution must take two limits in the proper order, and incorrect results are obtained if they are not taken in the proper order. One must first take the limit d t 0 to obtain the proper differential equation for any single time-evolution for one particular random walk, and then only afterward take the limit of an infinite number of random walks in an ensemble to get the average behavior.
Consider the entangled state of a system in which the fermion is in a superposition of both occupying and not occupying a state, n,
| ψ = α n | ψ 0 | 0 + β n | ψ 1 | 1 ,
where | 0 and | 1 represent the unoccupied and occupied localized fermion state, n, of interest, and | ψ 0 and | ψ 1 represent the associated many-body states of the entire rest of the system (which may be large superpositions of many-body states). The states | ψ 0 and | ψ 1 are orthonormal because they have different total numbers of fermions, and the total fermion number is conserved. For norm-conserving evolution, the coefficients α n and β n must always be normalized according to | α n | 2 + | β n | 2 = 1 . As discussed in Section 3, various physical processes can lead to time-dependent fluctuations of ω R , n , and these fluctuations of ω R , n lead to a random walk of the coefficients α n and β n , which has the end result of either α or β becoming nearly exactly 1 while the other becomes nearly exactly 0, with the probability of each given by the Born rule. This has the same end result as a projection onto one state followed by the normalization of the wave function, as in standard measurement theory. Therefore, this type of random walk cannot produce the results that a standard Born-rule measurement does not produce.
The projection of a state onto one of two localized fermion states is sufficient to give nonlocal correlation. Consider the entangled state
| Ψ = α 1 | 0 , 1 | ψ 1 + α 2 | 1 , 0 | ψ 2 ,
where the 0s and 1s represent the fermion occupation numbers of two states, which can be widely spatially separated, and | ψ 1 and | ψ 2 are the associated environments (again, these may be highly complicated, massive superpositions). Then, according to the model proposed here, fluctuations in the local environment of each of these states leads to a random walk of α n and β n for each electron state n via the term (1). Because the term acts on electrons in each of the atoms, when one electron state experiences an “upward” kick toward occupation N n = 1 , the other experiences a downward kick toward zero occupation. Eventually, this correlated random walk will end with one of the states having an occupation of 1, and the other having an occupation of 0, effectively projecting the full state onto one of the two states in (6). This argument can be generalized to any number of spatially separated states in a superposition.
Conservation rules. While the term (1) is not Hermitian, we can derive several conservation rules for it. First, it conserves the total number of fermions in the system. This can be seen by computing the time evolution of the expectation value of the total number operator N ^ = N ^ n predicted by the Schrödinger equation. By the same arguments as above, we keep only the terms that are linear in t, which gives, for an infinitesimal time lapse, t,
ψ t | N ^ | ψ t ψ 0 | N ^ | ψ 0 = ψ 0 | 1 + i V ^ d t N ^ 1 i V ^ d t | ψ 0 ψ 0 | N ^ | ψ 0 i d t ψ 0 | V ^ N ^ N ^ V ^ | ψ 0 = i N d t ψ 0 | V ^ V ^ | ψ 0 ,
where in the last line, we have used the property of | ψ 0 , meaning that it is a superposition of states all with the same total number of fermions, N, since any standard Hermitian interaction that changes the states of the individual fermions conserves the total fermion number. The change in N ^ is, therefore, 0 when using the results of (4).
This term also conserves the total energy of the system under the same conditions evoked for energy conservation when deriving irreversible behavior under Fermi’s golden rule and the quantum Boltzmann equation (see, e.g., Sections 4.7 and 4.8 of Ref. [29]). In this case, we can calculate the time dependence of the expectation value of H 0 = ω n N ^ n . Analogous to (7) above, we have
ψ t | H 0 | ψ t ψ 0 | H 0 | ψ 0 = i d t ψ 0 | V ^ H 0 H 0 V ^ | ψ 0 = d t n , n ω R , n ω n ψ 0 | N ^ n N ^ n N ^ n + N ^ n N ^ n N ^ n | ψ 0 = d t n , n ω R , n ω n 2 N ^ n N ^ n 2 N ^ n N ^ n .
This is equal to zero if N ^ n N ^ n factors to N ^ n N ^ n , which, as shown in Ref. [31], is equivalent to the strong decoherence limit, which is the same limit used to justify the quantum Boltzmann equation. The same approach can be taken if we add a number-conserving interaction to the Hamiltonian with terms that take the form A a i a j . In this case, we also find energy conservation as long as N ^ n a i a j = N ^ n a i a j , which is also the case when there is strong decoherence.

3. Numerical and Mathematical Results Showing the Born Rule

Ref. [27] showed that for any single fermion state, n, adding the operator (1) to the Schrödinger equation maps to the motion of a vector on a Bloch sphere, given by U = ( sin θ cos ϕ , sin θ sin ϕ , cos θ ) , representing the superposition (5) with α = e i ϕ sin θ / 2 and β = cos θ / 2 . The motion of this Bloch vector is then governed by the dynamical equations
U 2 t = ω R , n U 3 ( 1 U 3 2 ) 1 / 2 U 3 t = ω R , n U 2 ( 1 U 3 2 ) 1 / 2 .
Here, the factor ( 1 U 3 2 ) 1 / 2 gives two attractors corresponding to N ^ n = 0 and 1; the “collapse” of the wave function occurs when the Bloch vector hits either of these attractors. The dynamics of this Bloch vector model were analyzed numerically in Ref. [26]. When the Rabi factor ω R , n fluctuates randomly in time, the Bloch vector undergoes a random walk between the two attractors.
The t results of many random walks of this type have been shown numerically to give the Born probability rule with high accuracy. Figure 1 shows a comparison of the linear prediction of the Born rule for the probability of collapse to the t results of a numerical calculation using this model, using a value of ω R , n picked randomly for a sequence of time intervals d t , with the probability of the value within each time interval given by the Lorentzian (Cauchy) distribution
P ( ω R , n ) = 1 π γ ω R , n 2 + γ 2 ,
where γ gives the characteristic range of the fluctuations, with γ d t 1 and a cutoff of γ d t 0.5 . This is the expected distribution for a wide range of typical processes, in which the fluctuations have an exponentially decaying correlation function in time, i.e., no long-time correlation (see, e.g., Refs. [32,33]). The Wiener-Khinchtine theorem says that the frequency spectrum of a classical noise source is given by the Fourier transform of the temporal correlation function (see, e.g., Ref. [29], Section 9.5). For an exponentially decaying temporal correlation function C ( t ) e γ | t | , this implies
F ( ω ) = 1 2 π e γ | t | e i ω t d t = 1 2 π 0 e γ / τ e i ω t d t + 1 2 π 0 e t γ e i ω t d t = 1 π γ γ 2 + ω 2 .
Note that the function e γ | t | has a discontinuity in its slope at t = 0 . The Fourier transform of this gives infinitely high-frequency components in the spectral function, which are nonphysical. Realistically, this sharp discontinuity should be replaced by a smoothly curved peak at t = 0 over some small time range. This was accounted for in the numerical model by a high-frequency cutoff in the Lorentzian distribution.
The same linear result was found in numerical simulations using a Gaussian distribution of fluctuations of ω R , n . The Born rule holds for all of the random walks considered numerically on a finite interval in the t limit, which, mathematically, is because these walks are “martingale” walks. Martingales are defined in Appendix B, where it is shown that the numerical simulations will converge to +1 or −1 with probabilities linearly proportional to their starting position.

4. Quasi-Unitary Evolution Is the Same As Weak Measurement

As discussed in Section 2, the nonlinear operator (1) is similar in form to that suggested by Gisin [5] and others, and it has been argued that the nonlinear terms of this type intrinsically imply the possibility of superluminal communication. As shown in Section 2, however, the operator (1) is strictly norm-conserving, so that the Bloch vector for the full wave function (5) always has an exact unit length. This Bloch vector represents the full density matrix for the system as well as the environment that is entangled with it. When this entanglement with the environment is taken into account, the density matrix for any single random walk is always a pure state of the full many-body wave function, with a Bloch vector of unit length. As discussed in Section 1, in the t limit, a random walk of the Bloch vector implies that both the system and the environment collapse into one of two possibilities in accordance with the Born rule.
In Section 2, we asserted that the average of many random walks using the operator (1) will always give an evolution of the density matrix that is consistent with the standard Born rule of measurements in quantum mechanics. This can be seen explicitly by calculating the prediction of a random walk (“quantum trajectory”) using standard quantum mechanics in the “weak measurement” limit. Weak measurement theory [34,35,36,37] (for reviews, see, e.g., Refs. [38,39,40]) can be viewed as a type of spontaneous collapse, but it does not invoke any explicit non-unitary terms; instead, many weak Copenhagen-type measurements are assumed to occur, without specifying the observer. Recent work [41] has shown that weak measurement theory can be integrated into the decoherence approach of Zurek. Zurek’s approach (e.g., Ref. [42]) shows that the collapse of a density matrix to the diagonal can occur via unitary dynamics (sometimes called “first collapse,”) but cannot reproduce the statistics of individual quantum trajectories, i.e., single measurements (“second collapse”). Therefore, weak measurement theory is completely compatible with standard quantum mechanics but gives no insight into the cause of collapses. In this section, we show that weak measurement gives exactly the same operator as posited in (1).
Consider the following scenario of weak measurement, following the model from Ref. [43]. At time t = 0 , the initial state is
| ψ = | ϕ | X = 0 ,
where | X is the state of an external detector with a center-of-mass value at X (e.g., the position of a needle in a meter), and | ϕ = α | 1 + β | 1 is the internal state, which is in a superposition. The external state, in general, has position uncertainty in a Gaussian
| X = d x e ( x X ) 2 / 2 σ 2 | x ,
where | x 0 is the state of the external detector at exactly x 0 .
At t = 0 , a weak interaction of the system with the detector is turned on briefly, taking the form of
V ^ = g S z P ,
where S z is the standard spin operator, P = i / x is the momentum operator of the external center of mass X, and g is some small number. After a short time d t , the state of the system is
| ψ = | ψ + g d t d x | x α x e x 2 / 2 σ 2 | 1 + β x e x 2 / 2 σ 2 | 1 = | ψ + g d t d x | x α x σ 2 e x 2 / 2 σ 2 | 1 β x σ 2 e x 2 / 2 σ 2 | 1 = d x | x e x 2 / 2 σ 2 α 1 + x g d t σ 2 | 1 + β 1 x g d t σ 2 | 1 d x | x α e ( x g d t ) 2 / 2 σ 2 | 1 + β e ( x + g d t ) 2 / 2 σ 2 | 1 .
The last line is commonly used in discussing this scenario, but we will stick with the third line.
We now do a strong measurement of the external detector to collapse it to a definite value x = x 0 . By defining d f = x g d t / σ 2 1 , the state of the whole system is then
| ψ = 1 | α | 2 ( 1 + d f ) 2 + | β | 2 ( 1 d f ) 2 α ( 1 + d f ) | 1 + β ( 1 d f ) | 1 | x 0 1 1 + 2 | α | 2 d f 2 | β | 2 d f α ( 1 + d f ) | 1 + β ( 1 d f ) | 1 | x 0 1 | α | 2 d f + | β | 2 d f α ( 1 + d f ) | 1 + β ( 1 d f ) | 1 | x 0 .
The change in time is then
d | ψ = | ψ | ψ = g d t x 0 α ( 1 | α | 2 + | β | 2 ) | 1 β ( 1 + | α | 2 | β | 2 ) | 1 | x 0 = g d t x 0 ( | β | 2 | α | 2 ) ( α | 1 + β | 1 ) ( α | 1 + β | 1 | x 0 = g d t x 0 ( S z S z ) | ϕ | x 0 .
We have, thus, obtained a nonlinear operator of exactly the same form as (1) because the S z operator has exactly the same operation as N ^ n acting on two states; the nonlinearity arises from the fact that a standard quantum measurement is an intrinsically non-unitary process.
If σ is large, then x 0 is equally probable to be positive or negative. Therefore, if we apply the above process of weak measurement multiple times, with each time allowing the position of the needle in the detector to gain uncertainty due to normal wavepacket spreading, and then carrying out a measurement of the classical center of mass of the needle each time, we will get a random walk that is exactly equivalent to the results discussed in Section 2.
The significance of this mapping of weak measurement theory to the proposed non-Hermitian operator (1) for spontaneous collapse is that because weak measurement theory is based on standard measurement theory, the random walk with the non-Hermitian operator (1) cannot give any effects that violate known physics; in particular, the term (1) does not intrinsically allow superluminal communication or lead to a violation of energy conservation because weak measurement does not, and the two are mathematically identical. Due to environmental fluctuations, the random walk will give a density matrix corresponding to a non-unit-length Bloch vector only when the results of many trials in an ensemble are averaged; for any single random walk, it will always correspond to a Bloch vector of unit length.
Therefore, we can switch our viewpoint and adopt the operator (1) as a fundamental postulate. By the martingale rule discussed in Appendix B, the result of many sequential weak measurements will obey the Born rule, which is to say, many weak measurements will give a strong measurement, which makes sense since the same information is extracted, whether quickly or slowly. Thus, we can postulate that every strong measurement is actually the end result of many weak measurements of the form (1), which are intrinsic, rather than needing to rig the particular scenario discussed above.
This result depends crucially, however, on the assumption that the steps are random. In the case of weak measurement theory, this was true because of the axiomatic assumption of the Born rule for measurements. In the case of the operator (1), as a new term in the Schrödinger equation, the randomness comes from the randomness of ω R , n ( t ) . This, in turn, depends on what physical processes we believe affect ω R , n ( t ) . If ω R , n ( t ) is nonrandom and deterministically predictable, then the model presented in Section 2 does imply the possibility, at least in principle, of superluminal communication, as we will see below.

5. Possible Sources of Fluctuations and Experimental Implications

Let us now consider three different approaches to the physical source of the fluctuations of ω R , n in this model.
Universal fluctuations from a novel source. The first possibility we consider is that the fluctuations in ω R , n are unrelated to the local environment of the state n, and come from some universal source, such as gravity noise, dark matter, or some other fundamental field. This noise source would presumably still act locally in spacetime but would not be directly related to the presence of regular matter or energy. In this case, ω R , n ( t ) would be fundamentally unpredictable, and no superluminal communication would be possible. This approach has much in common with the GRW and related Penrose-Diósi models—it predicts the decoherence and collapse that are unrelated to the presence of any detector; particle states decohere in a pure vaccum on some length and time scale. However, it does not imply the nonconservation of energy; as is shown in Section 2, energy is conserved under the same conditions of decoherence as Fermi’s golden rule.
Local fluctuations of the environment. A second possibility, proposed in Refs. [26,27], is that ω R , n is directly the result of fluctuations of the local environment. For example, the fluctuation of the local energy density, which gives the fluctuation of the phase precession of state n, could give this. This has the appeal that measurement and collapse are directly related to the decoherence found in detectors and in any macroscopic system with strong decoherence.
In this case, we can posit that the full operator, which gives the nonlinear correction to the Schrödinger equation, is
V ^ ( t ) = i ξ d 3 r t H ( r , t ) Ψ ( r ) Ψ ( r ) Ψ ( r ) Ψ ( r ) ,
where Ψ ( r ) and Ψ ( r ) are spatial field operators (cf. Ref. [29], Section 4.6), H ( r , t ) is the standard energy operator from unitary physics, and ξ is a small parameter with units of time, which is a new physical constant. A fully relativistic version of this term is discussed in Appendix C.
The term (18) is equivalent to the form (1) under the assumption of coarse graining; that is, if we assume that / t H ( r , t ) is slowly varying in space for scales of length that are large compared to the localized states, n. In that case, we can write
V ^ = i ξ R t H ( R , t ) V ( R ) d 3 r Ψ ( r ) Ψ ( r ) Ψ ( r ) Ψ ( r ) ,
where R gives the positions of the coarse grains, and the integral is over the volume inside a coarse grain. We can then write the Fourier series
Ψ ( r ) = n ϕ n * ( r ) a n ,
where ϕ n ( r ) is the wave function, and a n is the creation operator for an eigenstate n, assuming that these eigenstates are localized within a coarse grain. Then, (19) becomes
V ^ = i ξ R t H ( R , t ) n , n V ( R ) d 3 r ϕ n * ( r ) ϕ n ( r ) a n a n a n a n .
When using the orthogonality relation d 3 r ϕ n * ( r ) ϕ n ( r ) = δ n , n then this gives us
V ^ = i ξ R , n t H ( R , t ) a n a n a n a n ,
which is equivalent to (1) under the assumption that the states, n, of interest are localized.
If the local environment has fluctuating energy density, then this will give the type of random walk discussed in Section 2, which reproduces the Born rule. Typical values of fluctuations in atom-based detectors are meV per picosecond, and the numerical results discussed above indicate that the order of 100 steps in a random walk is required to have a “collapse”; therefore, we can estimate ξ is of the order of 10 11 s/eV. This form does not violate energy conservation or produce heat because the fluctuations are driven entirely by real energy fluctuations that exist in the local environment.
With the form (19), however, we have the possibility of the external control of the detection statistics. Suppose that instead of waiting for local environmental fluctuations, we control H / t directly, for example, when using an intense laser pulse to give it a constant value that overcomes any local fluctuations. In that case, at least in principle, we can “jam” the detector locally to always detect a particle that is in violation of the Born rule by producing a rapid increase in the local energy density. Then, any remote particles entangled with the one we have jammed will also experience a violation of the Born rule, which would allow superluminal communication. A person could monitor a local detector for deviations from the Born rule with an agreed-upon flux of entangled particles sent between the two detectors, and a deviation from the Born rule could be registered as one bit of information, e.g., a bit value of 1 = for the Born rule, and a bit value of 0 = for the non-Born rule.
For the estimated value of ξ above and the atomic states of the order of the Bohr radius, jamming the detection by using this method would require a rate of energy density change in the order of 10 16 W/cm 3 . This would correspond to, for example, 10 kW of laser power absorbed in a volume with dimension of 1 micron. This is high but is within the range of modern technology.
Nonjammable environmental fluctuations. Suppose that the above experimental test is carried out and that no change from the Born-rule statistics is found in the entangled pairs. Does this falsify the proposed model? No, because it could be the case that form (1) is still correct but that the connection of ω R , n to the rate of change in the local energy density is incorrect. The form (1) allows for a wide variety of theories for what physical processes lead to the Rabi rotations. For example, the Rabi frequency could be proportional to the second derivative in time, rather than the first derivative, which would still be sensitive to local fluctuations but would require an acceleration of the local energy density increase to give detectable “jamming.” It could also be the case that the Rabi rotations are induced by a term proportional to the product of local fluctuations and a global, universal noise source, as discussed above. In that case, there would still be no decoherence in a vacuum, but the background fluctuations would prevent any jamming of the detection. It is also simply possible that the effect of the local fluctuations is not linear but saturates, e.g., instead of H / t , it is proportional to tan [ H / t ] .
If no jamming is possible and no unique experimental evidence is found, this model is, therefore, not falsified, but it loses the appeal of unique predictions. In that case, the appeal of this spontaneous collapse model would primarily be that it carries less philosophical baggage, with a fairly plebian description of reality, entailing just the evolution of waves with resonances and nonlinearities without such things as parallel universes and the injection of human consciousness as a metaphysical entity while not violating any of the experimental results of standard quantum mechanics.

6. Conclusions

The above considerations show that experimental tests of spontaneous collapse, such as superluminal communication, are not automatically implied by the existence of a nonlinear term in the Schrödinger equation in the form (1), but are, in principle, allowed for some versions of the physical mechanism that gives the fluctuations.
As discussed in Ref. [27], the most natural way to preserve a “narrativity” in spontaneous collapse theories of this type is to posit a preferred reference frame, such as the rest frame of the cosmic microwave background (i.e., the rest frame of the center of mass of the universe). In this case, even in the case of superluminal communication, there will be no grandfather paradoxes because any return communications must occur in the same preferred reference frame. Appendix C gives further discussion on how to implement this type of collapse in a relativistic framework.
As discussed above, if no physical effect such as jamming is seen, it does not automatically falsify the proposed mechanism because there may be intrinsic physical properties that always give a noisy signal. In that case, the appeal of this mechanism is that it puts measurement into the realm of the standard kinematic descriptions of waves. What the present model shows is that a spontaneous collapse model that agrees with all the known experimental results and is as logically coherent as possible.

Author Contributions

Supervision and numerical modeling, D.W.S.; Writing and editing, D.W.S. and D.N.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a guest scientist appointment at the Institute for Quantum Optics and Quantum Information of the Austrian Academy of Sciences, Vienna.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The numerical model used for this study as well as Ref. [26] can be found at https://notebookarchive.org/2022-07-5la04jb.

Acknowledgments

We thank Markus Arndt, Harvey Brown, Časlav Brukner, John Earman, Aritro Mukherjee, and David Wallace for helpful conversations.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Critique of the Historical Argument against a Simple Nonlinear Term

The standard argument that simple, spontaneous collapse models all involve superluminal communication is reviewed here.
Suppose that we have a source that emits a standard two-photon entangled state,
| ψ = α | V | V + β | H | H ,
where | V is a vertically linearly polarized state, and | H is horizontally polarized; the first ket refers to a left-going photon and the second to a right-going one.
The photon on the left then hits a standard detection apparatus oriented to detect vertical and horizontal polarization. It is assumed that this leads to the following density matrix:
ρ ^ 1 = | α | 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 | β | 2 ,
i.e., a mixed state of the projection onto states | V | V , | V | H , | H | V , and | H | H . If we define N ^ L = N ^ H N ^ V for the photon on the left, an operator of the form V ^ = i ω ( N ^ L N ^ L ) will then give
ρ ^ 1 t = ( i / ) ( V ^ ρ ^ 1 ρ ^ 1 V ^ ) = 2 ω ( | β | 2 | α | 2 ) | α | 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 | β | 2 2 ω | α | 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 | β | 2 = 2 ω ( | β | 2 | α | 2 + 1 ) | α | 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( | β | 2 | α | 2 1 ) | β | 2 .
If, for example, | α | 2 = . 5 , this gives
ρ ^ 1 t = ω 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 .
This leads to evolution toward the pure | V | V state.
On the other hand, suppose, at the last second, a person on the left changed the measurement to polarization along ± 45 instead of H and V. The recipe for getting the density matrix of the mixed state, in this case, is to project the state (A1) onto these axes and then create the stochastic mixture. The projector for + 45 acting on the states on the left is
P + 45 = 1 2 ( | H + | V ) ( H | + V | ) .
If we act with this operator on the whole state and normalize, we obtain
| ψ + 45 = 1 2 ( α | H | V + α | V | V + β | H | H + β | V | H ) ,
which occurs with probability ( | α | 2 + | β | 2 ) / 2 = 1 / 2 . The corresponding pure-density matrix is
ρ ^ + 45 = 1 2 | α | 2 α β * | α | 2 α β * α * β | β | 2 α * β | β | 2 | α | 2 α β * | α | 2 α β * α * β | β | 2 α * β | β | 2
The 45 polarization is similar:
ρ ^ 45 = 1 2 | α | 2 α β * | α | 2 α β * α * β | β | 2 α * β | β | 2 | α | 2 α β * | α | 2 α β * α * β | β | 2 α * β | β | 2
The mixed-density matrix that occurs is then
ρ ^ 2 = 1 2 | α | 2 0 0 α β * 0 | β | 2 α * β 0 0 α β * | α | 2 0 α * β 0 0 | β | 2 .
The total probability of getting a measurement of a photon with H-polarization on the right is still the same, namely | β | 2 (the sum of the two diagonal terms corresponding to a ket with a rightmost | H ), but there are now off-diagonal terms.
It is assumed that the non-Hermitian operator V ^ , in this case, involves the number operator N ^ L = N ^ + 45 N ^ 45 for the 45 polarizations. The creation operator for 45 polariztion is a 45 = a V + a H , which implies the number operator N 45 = a 45 a 45 = 1 2 ( a V + a H ) ( a V + a H ) = 1 2 ( N V + N H + a V a H + a H a V ) . Similarly, N 45 = 1 2 ( N V + N H a V a H a H a V ) , and therefore N 45 N 45 = a V a H + a H a V . The time dependence is then given by
ρ ^ 2 t = ( i / ) ( V ^ ρ ^ 2 ρ ^ 2 V ^ ) ,
with V ^ = i ω ( N ^ L N ^ L ) as before, but in this case, the expectation value of N ^ L = N ^ + 45 N ^ 45 is zero since there is an equal likelihood of both polarizations. The time dependence is then given just by the N ^ L operator,
N ^ L = 0 0 1 0 0 0 0 1 1 0 0 0 0 1 0 0 .
This gives
ρ ^ 2 t = ω 0 α β * | α | 2 0 α * β 0 0 | β | 2 | α | 2 0 0 α β * 0 | β | 2 α * β 0 .
By comparing this to (A7), we see that the action is to move the mixed-density matrix toward the density matrix for the pure 45 polarization state. The relative probability of detecting single photons with H and V polarizations is unchanged.
Therefore, the change in the count rates is different for the detector on the right, depending on what measurement setting the person on the left has chosen since in the first case above, with vertical and horizontal detection, the non-Hermitian term was found to give a change of the relative probabilities. The person on the right can, therefore, presumably detect the position of the detector on the left instantaneously, which would imply superluminal communication.
As discussed in the main text, this argument depends crucially on first creating the mixed state and then, afterward, applying the non-Hermitian operator. However, the probabilistic nature of the mixed state assumes the existence of an ensemble of trials in which multiple collapses or projections have occurred. Applying the non-Hermitian operator afterward assumes that it acts the same way, deterministically, for every one of these trials in the ensemble. However, the assumption of our model is that the non-Hermitian operator is the cause of the collapses and that it acts stochastically to give random outcomes in each trial, in accordance with the Born rule. The proper recipe, therefore, is to first form the density matrix of the pure state and then apply the non-Hermitian operator to get a collapse in one particular trial and, last, to average over many of these trials to get the final mixed-density matrix.
One must be careful when carrying this out. The approaches that use Itô calculus (for a review, see [44]) to perform the stochastic average can also implicitly rely on improper assumptions. For example, Gisin and Percival [45] posit the form of the Schrödinger equation (in terms of the operators used here)
| ψ t = | ψ 0 i H | ψ 0 d t 1 2 n N ^ n 2 + N ^ n 2 2 N ^ n N ^ n | ψ 0 d t + n ( N ^ n N ^ n ) | ψ 0 d W n ,
where d W n is a (real-valued) Itô stochastic differential accounting for fluctuations. This can be rewritten as
| ψ t = | ψ 0 i H | ψ 0 d t 1 2 n N ^ n N ^ n 2 | ψ 0 d t + n ( N ^ n N ^ n ) | ψ 0 d W n = | ψ 0 i H | ψ 0 d t 1 2 n V ^ n 2 | ψ 0 d t + n V ^ n | ψ 0 d W n ,
and V ^ n = N ^ n N ^ n . The averaged density matrix can then be constructed as
| ψ t ψ t | = | ψ 0 i H | ψ 0 d t 1 2 n V ^ n 2 | ψ d t + n V ^ n | ψ d W n × ψ 0 | + i ψ 0 | H d t 1 2 n ψ 0 | V ^ n 2 d t + n ψ 0 | V ^ n d W n .
As per Itô calculus, those terms that are linear in d W n vanish (the mean of the random walk remains unchanged) while d W n d W m = δ n m d t (fluctuation of different states are uncorrelated, and the magnitude of the random walk increases as t 1 / 2 ). This gives us
| ψ t ψ t | = | ψ 0 ψ 0 | i [ H , | ψ 0 ψ 0 | ] 1 2 n V ^ n 2 | ψ 0 ψ 0 | d t 1 2 n | ψ 0 ψ 0 | V ^ n 2 d t + n V ^ n | ψ ψ 0 | V ^ n .
The last three terms are exactly the Lindbladian for the operator V ^ n . It can be shown that this term gives the decay of the off-diagonal terms of the density matrix while preserving the diagonal terms, leading to a diagonal density matrix that corresponds to a mixed state. This, in turn, prevents superluminal communication of the type discussed above since the diagonal of the mixed-density matrix is always left unchanged.
This approach of Gisin and Percival requires fine-tuning in that the magnitude of the third term in (A13) must match exactly the amplitude of the noise assumed in the fourth term. Without the third term, the Itô calculus applied to just the fourth term, as taken by Gisin and Percival, will not preserve the diagonal of the density matrix. But we argue here that such a term is not needed when physically realistic fluctuations are accounted for correctly.
Itô calculus intrinsically assumes that the fluctuations are fast compared to all other time scales; in particular, it assumes that fluctuations occur on time scales that are short compared to d t . However, in our model, as discussed in Section 5, the fluctuations are assumed to arise from physical fluctuations in the environment. As such, these have a short but nonzero persistence time, which means that one can always pick a d t that is short enough such that ω R , n can be treated as constant. Therefore, one should first take the d t 0 to derive the time evolution and, only afterward, take the stochastic average.
For the moment, treating ω R , n ( t ) as not time-varying, we can then use standard time-dependent perturbation theory in the interaction representation, as in Section 2 (see, e.g., Ref. [29], Chapter 4). Starting with the Schrödinger equation
i d d t | ψ t = H 0 | ψ t + i n V ^ n ω R , n | ψ t ,
the prescription of this method is to write
| ψ ( t ) = | ψ 0 + 1 i 0 t d t n V ^ n ( t ) | ψ ( t ) ,
where V n ( t ) = e i H 0 t / V ^ n e i H 0 t / and | ψ ( t ) = e i H 0 t / | ψ t . Since N ^ n commutes with the Hamiltonian H 0 , the only time-dependence of V ^ n ( t ) comes from the time dependence of | ψ t , which is used in the term N ^ n , which occurs in V ^ n . The expansion of the time dependence up to the first order is then
| ψ ( t ) 1 + n ω R , n t ψ 0 | N n ^ | ψ 0 N ^ n | ψ 0 .
When substituting this back into (A18) for the second-order expansion, we then have
| ψ ( t ) | ψ 0 + 0 t d t n ω R , n ψ 0 | 1 + m ω R , m t ψ 0 | N m ^ | ψ 0 N ^ m × N ^ n 1 + m ω R , m t ψ 0 | N ^ m | ψ 0 N ^ m | ψ 0 N ^ n × 1 + n ω R , n t ψ 0 | N ^ n | ψ 0 N ^ n | ψ 0 | ψ 0 + n ω R , n t ψ 0 | N ^ n | ψ 0 N ^ n | ψ 0 + t 2 2 n , n ω R , n ω R , n ψ 0 | N ^ n | ψ 0 N ^ n ψ 0 | N ^ n | ψ 0 N ^ n | ψ 0 + t 2 n , m ω R , n ω R , m ψ 0 | ψ 0 | N ^ m | ψ 0 N ^ m N ^ n | ψ 0 .
When taking t 2 as being negligible compared to t, as we did in Section 2, we form the density matrix
| ψ t ψ t | = e i H 0 t / | ψ 0 + t n V ^ n ω R , n | ψ 0 ψ 0 | + t n V ^ n ω R , n ψ 0 | e i H 0 t / .
Again, by keeping just the terms linear in t, we then have
| ψ t ψ t | = | ψ t ψ t | + t n ω R , n | ψ t ψ t | V ^ n + V ^ n | ψ t ψ t | ,
where | ψ t = e i H 0 t / | ψ 0 is the change in the state by normal Hamiltonian evolution. In other words, when taking t = d t as an infinitesimal, the change in the density matrix due to the non-Hermitian term is
d d t | ψ t ψ t | = n ω R , n | ψ t ψ t | V ^ n + V ^ n | ψ t ψ t | .
The action of the operator | ψ ψ | V ^ n + V ^ n | ψ ψ | on the two states | 0 and | 1 , for | ψ = α | 1 + β | 0 is given by the matrix
2 | α | 2 | β | 2 α * β ( 2 | α | 2 1 ) β * α ( 2 | α | 2 1 ) 2 | α | 2 | β | 2 .
As discussed in Ref. [27], this maps to a one-dimensional problem corresponding to a random walk of the polar angle θ on the Bloch sphere, where cos θ = | α | 2 | β | 2 ; the off-diagonal components of the density matrix are not independent of the diagonal components but instead approach zero as the diagonal components approach either | α | = 1 or | β | = 1 .
We can now take the stochastic average for fluctuations in ω R , n . Since the fluctuations of θ are linearly proportional to ω R , n , the Itô calculus prescribes that the average of θ does not change in time. This implies that an average over many random walks will give
| ψ ψ | = | α | 2 | 0 0 | β | 2 ,
in accordance with the Born rule. This is an alternate derivation of the proof given in Appendix B.

Appendix B. Proof That a Martingale Random Walk Satisfies the Born Rule

Let Y 1 , Y 2 , be a sequence of independent and identically distributed random variables with a common symmetric, continuous probability density ϕ that is strictly positive in an open neighborhood of the origin. Examples of such ϕ include the Gaussian and Lorenztian distributions. We also suppose that 0 < L < 1 is a fixed constant, small enough so that
L c 1 / 2 ,
where c > 0 is a cutoff parameter that we will make use of later. For convenience of numerical implementation, we want to ensure that the step sizes are uniformly bounded. We define σ to be the function
σ ( x ) = x if | x | c c if x c + c if x + c .
Finally, we define σ n = σ ( Y n ) . For each fixed n, σ n is, thus, a continuous, symmetric random variable with mean 0 and uniformly bounded range [ c , c ] .
We wish to study the following stochastic finite difference equation (SFDE):
X n + 1 = X n + σ n + 1 L ( 1 X n 2 ) X 0 [ 1 , 1 ] .
This SFDE is a random walk starting at a given X 0 [ 1 , 1 ] . We shall show that this model leads to the Born rule. That is, when given any initial condition X 0 [ 1 , 1 ] , we can assign a probability, p + ( X 0 ) , of collapsing to state 1, and a probability, p ( X 0 ) , of collapsing to state 1 , such that
( Born rule ) p ± ( X 0 ) = 1 ± X 0 2 .
Our strategy is to first prove that we can achieve the Born rule within arbitrary precision ϵ > 0 , and then to take the limit as ϵ tends to 0 to recover the Born rule. The crucial observation we will make is that the stochastic process ( X n ) defined by (A28) is a martingale. The martingale property is that
E [ X n + 1 | X 0 , , X n ] = X n ,
where E denotes the conditional expectation. Intuitively, our best prediction of the next position of a martingale, given the history of the process, is simply the current position (martingales are like stochastic constants). In order to prove that our process is a martingale, we simply compute the conditional expectation.
First, we note that X n + 1 depends only on X n and not X 0 , , X n 1 . Moreover, X n is known to us at time n. Finally, σ n + 1 is independent of X n . According to the properties of conditional expectation, we, therefore, have
E [ X n + 1 | X 0 , , X n ] = E [ X n + 1 | X n ] = E [ X n + σ n + 1 L ( 1 X n 2 ) | X n ] = X n + L ( 1 X n 2 ) E [ σ n + 1 | X n ] = X n + L ( 1 X n 2 ) E [ σ n + 1 ] = X n .
In the last line, we have used E [ σ n + 1 ] = 0 for every n. Therefore, the sequence ( X n ) is a martingale.
Now, we introduce 0 < ϵ < 1 and define a stopping time
N ( ϵ ) = inf { n 0 : | X n | 1 ϵ } .
We state here and prove below that N ( ϵ ) is finite almost surely, which simply means that P [ N ( ϵ ) < ] = 1 . Because of this, we can make use of the optional stopping theorem, ref. [46] which states that if a martingale ( Z n ) is bounded in the sense that there exists some M > 0 , such that | Z n | < M for all n 0 , and if a stopping time τ is finite almost surely, then E [ Z τ ] = E [ Z 0 ] .
We have already established that X n is a martingale, and we are assuming that N ( ϵ ) is finite almost surely. Thus, we only need to find a uniform bound on | X n | to apply the optional stopping theorem. A bound is achieved as follows: for all n, | X n + 1 | | X n | + L c | 1 X n 2 | 1 + L c according to the triangle inequality. Thus, the optional stopping theorem applies, and we can conclude that
E [ X N ( ϵ ) ] = E [ X 0 ] = X 0 .
With (A33) in hand, we proceed as follows: Let p ϵ + ( X 0 ) = P [ X N ( ϵ ) 1 ϵ ] be the probability that collapse to + 1 occurs (before collapse to 1 ) with a tolerance of ϵ , and let p ϵ ( X 0 ) = P [ X N ( ϵ ) 1 + ϵ ] be the probability that collapse to 1 occurs (before collapse to + 1 ) with a tolerance of ϵ . Since we are assuming that P [ N ( ϵ ) < ] = 1 , we know that at least one of the two events ( X N ( ϵ ) 1 ϵ ) or ( X N ( ϵ ) 1 + ϵ ) must occur. Moreover, for sufficiently small ϵ , these two events are complements. Therefore, the probabilities satisfy p ϵ + ( X 0 ) = 1 p ϵ ( X 0 ) .
Let b = b ( ϵ ) = X N ( ϵ ) if the event ( X N ( ϵ ) 1 + ϵ ) occurred, and let a = a ( ϵ ) = X N ( ϵ ) if the event ( X N ( ϵ ) 1 ϵ ) occurred. We have b 1 + ϵ and a ϵ 1 , but we will also need a sufficiently large lower bound on b, as well as a sufficiently small lower bound on a. It turns out that, in fact, 1 b and a 1 . To prove this, consider the function f ( x ) = | x | + 1 2 | 1 x 2 | defined on the interval [ 1 , 1 ] . It is straightforward to find that the maximum value of f on this closed interval is 1. Thus, we can estimate the position of the random walk at the stopping time N = N ( ϵ ) :
| X N | | X N 1 | + L c | 1 X N 1 2 | | X N 1 | + 1 2 | 1 X N 1 2 | by ( A 26 ) = f ( X N 1 ) 1 ,
since, prior to the stopping time, X N 1 must be an element of [ 1 , 1 ] .
Now, we can expand out the expected value in (A33):
X 0 = E [ X N ( ϵ ) ] = p ϵ + ( X 0 ) a + ( 1 p ϵ + ( X 0 ) ) b .
Solving for p ϵ + ( X 0 ) yields
p ϵ + ( X 0 ) = b + X 0 a b .
Since we have shown that 1 b 1 + ϵ and 1 ϵ a 1 , it follows that lim ϵ 0 a = 1 and lim ϵ 0 b = 1 ; hence, lim ϵ 0 p ϵ + ( x 0 ) = 1 + X 0 2 . Likewise, one finds that lim ϵ 0 p ϵ ( x 0 ) = 1 X 0 2 . Therefore, (A29) holds.
Proof that N ( ϵ ) is bounded almost surely. We wish to show that P [ N ( ϵ ) < ] = 1 . Choose δ > 0 so that ϕ > 0 on [ δ / 2 , δ / 2 ] and such that 0 < δ < c . We claim that the stochastic process { X n : n = 0 , 1 , 2 , } defined via (A28) exhibits arbitrarily long sequences of consecutive, equally sized steps in the direction of + 1 . Let K be an arbitrary natural number, and let E n K be the event ( σ n + i δ / 2 for i = 1 , , K ) . Consider the sequence of independent events E 1 K , E K + 1 K , E 2 K + 1 K , . From our assumptions on the probability density ϕ for m = 0 , 1 , 2 , , we have
P [ E m K + 1 K ] = P [ σ n δ / 2 ] K δ / 2 δ ϕ ( x ) d x K = : λ > 0 .
Therefore,
m = 0 P [ E m K + 1 K ] m = 0 λ = .
Since the events E m K + 1 K are independent, the second Borel-Cantelli lemma applies, and P [ E m K + 1 K infinitely often ] = 1 . We can, thus, choose an m so that E m K + 1 K holds almost surely. Consider now the SFDE (after reindexing, we can assume the process starts at time m K )
X n + 1 = X n + σ n + 1 L ( 1 X n 2 ) , n = 1 , , K X 0 [ 1 , 1 ]
Consider the analogous deterministic finite difference equation
y n + 1 = y n + δ 2 L ( 1 y n 2 ) , n 1 y 0 [ 1 , 1 ]
It follows from the consideration of the function f defined above and the fact that δ L / 2 < c L / 2 < c L < 1 / 2 that, for system (A40), | y n | 1 for all n. Moreover, { y n } is monotone-increasing. Thus, the sequence converges upward to a limit, y. Then, y y = δ 2 L ( 1 ( y ) 2 ) , so y = ± 1 . Thus, y = 1 . Hence, we can find a least M < , so that 1 ϵ y M 1 , and { n < : | y n | 1 ϵ } . Now, select K to be large enough so that M < K . For every n = m K , , ( m + 1 ) K , we have σ n δ / 2 ; thus, X n y n . Hence, (almost surely)
N ( ϵ ) = inf { n : | X n | 1 ϵ } inf { n : | y n | 1 ϵ } = M < .

Appendix C. Relativistic Formalism for Nonlocal Projections

The non-unitary term (18) has the appeal that it is a compact addition to the Schrödinger equation that allows numerical methods to be applied. It also has nonlocal effects, acting as a type of projector at equal times over all space.
As discussed in Ref. [27], one can generalize this to other relativistic reference frames by positing that there is a favored, universal reference frame in which these projections occur, which could be the center of mass of the universe, seen in the cosmic microwave background, or some other chosen frame. In this case, we can take the many-body states of the system, defined along equal-time slices in this favored reference frame.
A many-body “state” in the laboratory rest frame at time t 0 can be written as
| ψ = α | 0 , 1 , 0 , 0 , 1 , + β | 1 , 1 , 0 , 1 , 0 , +
where α and β are complex phase factors for a potentially infinite number of different superpositions of Fock states, and the Fock states are generated by the factor ( a ( r , t 0 ) ) N acting at each point r in space, where N can equal 0 or 1. This is somewhat unconventional notation, as the state created by a ( r , t 0 ) is not an eigenstate of the system, but the set of Fock states defined this way is a complete set of orthogonal states that can form the basis of any many-body state. This can be carried out simply by switching the roles of r and k in the standard plane-wave basis of the many-body field theory (see, e.g., Ref. [29], Chapter 4) to give a discrete sum over r -states and a continuum of k -states created by Ψ ( k ) . In this notation, a plane-wave state with momentum k is created as
| k = Ψ k | vac = 1 k max 3 r e i ( k · r ω k t ) a ( r , t ) | vac ,
which is a superposition of Fock states, each with a single particle at one location. Here, k max is a cutoff of momentum, which plays the same role as a system size cutoff V in conventional notation, and the limit is taken of k max . (It has been questioned [47,48] whether the limit V can be taken without causing a breakdown of the conventional many-body theory, as the finite value of V is used to justify discrete k -states. Here we simply note that when taking V as infinite from the start, with k max finite, and, afterward, taking k max , it is formally and mathematically equivalent to the conventional approach of taking k defined over an infinite range with V being finite, and then taking the limit V , since we have just redefined the names of the conjugate variables.)
The state (A41) can be visualized as the horizontal line in Figure A1 with open circles, where each circle represents a local fermionic state in a continuum. It amounts to taking a slice at equal times in the overall hyperspace field. If an operator N ^ ( r , t 0 ) is applied to this state, this will remove any Fock states in the overall superposition that have N ( r , t 0 ) = 0 .
The time-slice on which a many-body “state” is defined is somewhat arbitrary, however. We could also define a “state” as a slice of the global field at equal times in a different, chosen reference frame (e.g., the rest frame of the cosmic microwave background). If that chosen reference frame is moving at speed v relative to the laboratory rest frame, its equal time slice will correspond and have the same form as (A41), but in this case of the local creation operators, it will be Ψ ( r , t 0 r · v / c 2 ) . This corresponds to the tilted line with open circles in Figure A1. If an operator N ^ ( r , t 0 ) is applied to this state, this will remove any Fock states in this overall superposition that have N ( r , t 0 ) = 0 . This implies a backwards-in-time action of a sort, but as discussed in Ref. [27], it does not imply grandfather paradoxes.
Of course, we could just solve all problems in the chosen reference frame, in which case we do not have to worry about backwards-in-time actions. But in many cases, it may be more convenient to solve the normal Hermitian Hamiltonian dynamics in the laboratory frame and then apply the nonlocal, non-Hermitian term separately. In this case, we have a well-defined algorithm for evolving the field in time:
  • For a given many-body “state” along a time slice, first evolve the state using the local, Hermitian interaction Hamiltonian to take the state from time t to time t + d t . This may generate superpositions of different Fock states.
  • At each point of this new “state”, apply the non-Hermitian operator (18). This will have the effect of changing the relative weights of the phase factors α , β of the overall superposition of Fock states.
  • With this new superposition, go on to the next time slice t + d t and start the process over.
This algorithm may seem odd, but it is well-defined, is consistent between different reference frames, and involves no random numbers other than what randomness is generated by actual fluctuations in the local environment via the H term in (18).
Figure A1. Horizontal line with open circles: the standard definition of a many-body state at equal times in the observer’s rest frame. The circles represent creation at multiple locations r via Ψ ( r , t = 0 ) . Tilted line with open circles: the definition of a many-body state at equal times in a universal, favored reference frame, moving at speed v relative to the rest frame of the horizontal axis (with the time axis given by the line with an arrow, t = r / v ). The circles represent creation at multiple locations r via Ψ ( r , t r · v / c 2 ) .
Figure A1. Horizontal line with open circles: the standard definition of a many-body state at equal times in the observer’s rest frame. The circles represent creation at multiple locations r via Ψ ( r , t = 0 ) . Tilted line with open circles: the definition of a many-body state at equal times in a universal, favored reference frame, moving at speed v relative to the rest frame of the horizontal axis (with the time axis given by the line with an arrow, t = r / v ). The circles represent creation at multiple locations r via Ψ ( r , t r · v / c 2 ) .
Entropy 25 01489 g0a1

References

  1. Penrose, R. On the gravitization of quantum mechanics 1: Quantum state reduction. Found. Phys. 2014, 44, 557–575. [Google Scholar] [CrossRef]
  2. Penrose, R. On gravity’s role in quantum state reduction. Gen. Relativ. Gravit. 1996, 28, 581–600. [Google Scholar] [CrossRef]
  3. Diósi, L. On the conjectured gravity-related collapse rate EΔ/ of massive quantum superpositions. AVS Quantum Sci. 2022, 4, 015605. [Google Scholar] [CrossRef]
  4. Pearle, P. Reduction of the state vector by a nonlinear Schrödinger equation. Phys. Rev. D 1976, 13, 857. [Google Scholar] [CrossRef]
  5. Gisin, N. Quantum measurements and stochastic processes. Phys. Rev. Lett. 1984, 52, 1657. [Google Scholar] [CrossRef]
  6. Gisin, N. Weinberg’s non-linear quantum mechanics and supraluminal communications. Phys. Lett. 1990, 143A, 1–2. [Google Scholar] [CrossRef]
  7. Ghirardi, G.C.; Rimini, A.; Weber, T. Unified dynamics for microscopic and macroscopic systems. Phys. Rev. D 1986, 34, 470. [Google Scholar] [CrossRef]
  8. Tumulka, R. A relativistic version of the Ghirardi-Rimini-Weber model. J. Stat. Phys. 2006, 125, 825–840. [Google Scholar] [CrossRef]
  9. Diósi, L. Quantum stochastic processes as models for state vector reduction. J. Phys. A 1988, 21, 2885. [Google Scholar] [CrossRef]
  10. Diósi, L. Models for universal reduction of macroscopic quantum fluctuations. Phys. Rev. A 1989, 40, 1165. [Google Scholar] [CrossRef]
  11. Gisin, N.; Rigo, M. Relevant and irrelevant nonlinear Schrödinger equations. J. Phys. A 1995, 28, 7375. [Google Scholar] [CrossRef]
  12. Smirne, A.; Bassi, A. Dissipative continuous spontaneous localization (CSL) model. Sci. Rep. 2015, 5, 12518. [Google Scholar] [CrossRef] [PubMed]
  13. Bassi, A.; Lochan, K.; Satin, S.; Singh, T.P.; Ulbricht, H. Models of wave-function collapse, underlying theories, and experimental tests. Rev. Mod. Phys. 2013, 85, 471. [Google Scholar] [CrossRef]
  14. Carlesso, M.; Donadi, S.; Ferialdi, L.; Paternostro, M.; Ulbricht, H.; Bassi, A. Present status and future challenges of non-interferometric tests of collapse models. Nat. Phys. 2022, 18, 243–250. [Google Scholar] [CrossRef]
  15. Fein, Y.Y.; Geyer, P.; Zwick, P.; Kiaka, F.; Pedalino, S.; Mayor, M.; Gerlich, S.; Arndt, M. Quantum superposition of molecules beyond 25 kDa. Nat. Phys. 2019, 15, 1242–1245. [Google Scholar] [CrossRef]
  16. Carlesso, M.; Bassi, A.; Falferi, P.; Vinante, A. Experimental bounds on collapse models from gravitational wave detectors. Phys. Rev. D 2016, 94, 124036. [Google Scholar] [CrossRef]
  17. Gasbarri, G.; Belenchia, A.; Carlesso, M.; Donadi, S.; Bassi, A.; Kaltenbaek, R.; Paternostro, M.; Ulbricht, H. Testing the foundation of quantum physics in space via Interferometric and non-interferometric experiments with mesoscopic nanoparticles. Commun. Phys. 2021, 4, 155. [Google Scholar] [CrossRef]
  18. Donadi, S.; Bassi, A.; Ferialdi, L.; Curceanu, C. The effect of spontaneous collapses on neutrino oscillations. Found. Phys. 2013, 43, 1066–1089. [Google Scholar] [CrossRef]
  19. Bahrami, M.; Donadi, S.; Ferialdi, L.; Bassi, A.; Curceanu, C.; Di Domenico, A.; Hiesmayr, B.C. Are collapse models testable with quantum oscillating systems? The case of neutrinos, kaons, chiral molecules. Sci. Rep. 2013, 3, 1952. [Google Scholar] [CrossRef]
  20. Adler, S.L.; Bassi, A.; Carlesso, M.; Vinante, A. Testing continuous spontaneous localization with Fermi liquids. Phys. Rev. D 2019, 99, 103001. [Google Scholar] [CrossRef]
  21. Carlesso, M.; Paternostro, M.; Ulbricht, H.; Vinante, A.; Bassi, A. Non-interferometric test of the continuous spontaneous localization model based on rotational optomechanics. New J. Phys. 2018, 20, 083022. [Google Scholar] [CrossRef]
  22. Vinante, A.; Mezzena, R.; Falferi, P.; Carlesso, M.; Bassi, A. Improved noninterferometric test of collapse models using ultracold cantilevers. Phys. Rev. Lett. 2017, 119, 110401. [Google Scholar] [CrossRef] [PubMed]
  23. Piscicchia, K.; Bassi, A.; Curceanu, C.; Del Grande, R.; Donadi, S.; Hiesmayr, B.C.; Pichler, A. CSL collapse model mapped with spontaneous radiation. Entropy 2017, 19, 319. [Google Scholar] [CrossRef]
  24. Donadi, S.; Piscicchia, K.; Curceanu, C.; Diósi, L.; Laubenstein, M.; Bassi, A. Underground test of gravity-related wave function collapse. Nat. Phys. 2021, 17, 74–78. [Google Scholar] [CrossRef]
  25. Smirne, A.; Vacchini, B.; Bassi, A. Dissipative extension of the Ghirardi-Rimini-Weber model. Phys. Rev. A 2014, 90, 062135. [Google Scholar] [CrossRef]
  26. Snoke, D.W. A model of spontaneous collapse with energy conservation. Found. Phys. 2021, 51, 100. [Google Scholar] [CrossRef]
  27. Snoke, D.W. Mathematical formalism for nonlocal spontaneous collapse in quantum field theory. Found. Phys. 2023, 53, 34. [Google Scholar] [CrossRef]
  28. Mukherjee, A.; Gotur, S.; Aalberts, J.; van den Ende, R.; Mertens, L.; van Wezel, J. Quantum state reduction of general initial states through spontaneous unitarity violation. arXiv 2023, arXiv:2301.03233. [Google Scholar]
  29. Snoke, D.W. Solid State Physics: Essential Concepts, 2nd ed.; Cambridge University Press: Cambridge, UK, 2020. [Google Scholar]
  30. Snoke, D.W. Interpreting Quantum Mechanics: Modern Foundations; Cambridge University Press: Cambridge, UK, in press.
  31. Snoke, D.W.; Liu, G.-Q.; Girvin, S.M. The basis of the second law of thermodynamics in quantum field theory. Ann. Phys. 2012, 327, 1825–1851. [Google Scholar] [CrossRef]
  32. Abragam, A. The Principles of Nuclear Magnetism; Clarendon Press: Oxford, UK, 1961. [Google Scholar]
  33. Chaikin, P.M.; Lubensky, T.C. Principles of Condensed Matter Physics; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  34. Diósi, L. Continuous quantum measurement and Ito formalism. Phys. Lett. A 1988, 129, 419–423. [Google Scholar] [CrossRef]
  35. Aharonov, Y.; Albert, D.Z.; Vaidman, L. How the result of a measurement of a component of the spin of a spin-1/2 particle can turn out to be 100. Phys. Rev. Lett. 1988, 60, 1351. [Google Scholar] [CrossRef] [PubMed]
  36. Duck, I.M.; Stevenson, P.M.; Sudarshan, E.C.G. The sense in which a “weak measurement” of a spin-1/2 particle’s spin component yields a value 100. Phys. Rev. D 1989, 40, 2112. [Google Scholar] [CrossRef] [PubMed]
  37. Knight, J.M.; Vaidman, L. Weak measurement of photon polarization. Phys. Lett. A 1990, 143, 357–361. [Google Scholar] [CrossRef]
  38. Braginsky, V.B.; Khalili, F.Y. Quantum Measurement; Cambridge University Press: Cambridge, UK, 1992. [Google Scholar]
  39. Jacobs, K.; Steck, D.A. A straightforward introduction to continuous quantum measurement. Contemp. Phys. 2006, 47, 279–303. [Google Scholar] [CrossRef]
  40. Clerk, A.; Devoret, M.; Girvin, S.; Marquardt, F.; Schoelkopf, R. Introduction to quantum noise, measurement, and amplification. Rev. Mod. Phys. 2010, 82, 1155. [Google Scholar] [CrossRef]
  41. Dittrich, T.; Rodríguez, O.; Viviescas, C. Simulating spin measurement with a finite heat bath model for the environment. Phys. Rev. A 2022, 106, 042203. [Google Scholar] [CrossRef]
  42. Zurek, W. Probabilities from entanglement, Born’s rule pk=|ψk|2 from envariance. Phys. Rev. A 2005, 71, 052105. [Google Scholar] [CrossRef]
  43. Tamir, B.; Cohen, E. Introduction to weak measurements and weak values. Quanta 2013, 2, 7–17. [Google Scholar] [CrossRef]
  44. Percival, I. Quantum State Diffusion; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  45. Gisin, N.; Percival, I.C. Quantum State Diffusion: From Foundations to Applications. arXiv 1997, arXiv:quant-ph/9701024. [Google Scholar]
  46. Baudoin, F. Diffusion Processes and Stochastic Calculus; European Mathematical Society Press: Zürich, Switzerland, 2014; pp. 18–20. [Google Scholar]
  47. Haag, R. On Quantum Field Theories. Danish Royal Academy of Sciences and Letters: Mathematical-Physical Messages. Dan. Mat. Fys. Medd. 1955, 29, 1–37. [Google Scholar]
  48. Earman, J.; Fraser, D. Haag’s theorem and its implications for the foundations of quantum field theory. Erkenntnis 2006, 63, 305–344. [Google Scholar] [CrossRef]
Figure 1. Comparison of the numerical results for the model with the Born rule (solid line). Circles: probability of end state N ^ n = 1 for γ d t = 0.002 . Squares: γ d t = 0.02 . For each value of γ d t , 8000 trials were run for each data point, and the number of steps in each trial was limited to 4000.
Figure 1. Comparison of the numerical results for the model with the Born rule (solid line). Circles: probability of end state N ^ n = 1 for γ d t = 0.002 . Squares: γ d t = 0.02 . For each value of γ d t , 8000 trials were run for each data point, and the number of steps in each trial was limited to 4000.
Entropy 25 01489 g001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Snoke, D.W.; Maienshein, D.N. Experimental Predictions for Norm-Conserving Spontaneous Collapse. Entropy 2023, 25, 1489. https://doi.org/10.3390/e25111489

AMA Style

Snoke DW, Maienshein DN. Experimental Predictions for Norm-Conserving Spontaneous Collapse. Entropy. 2023; 25(11):1489. https://doi.org/10.3390/e25111489

Chicago/Turabian Style

Snoke, D. W., and D. N. Maienshein. 2023. "Experimental Predictions for Norm-Conserving Spontaneous Collapse" Entropy 25, no. 11: 1489. https://doi.org/10.3390/e25111489

APA Style

Snoke, D. W., & Maienshein, D. N. (2023). Experimental Predictions for Norm-Conserving Spontaneous Collapse. Entropy, 25(11), 1489. https://doi.org/10.3390/e25111489

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