Master equations for Wigner functions with spontaneous collapse and their relation to thermodynamic irreversibility

Wigner functions, allowing for a reformulation of quantum mechanics in phase space, are of central importance for the study of the quantum-classical transition. A full understanding of the quantum-classical transition, however, also requires an explanation for the absence of macroscopic superpositions to solve the quantum measurement problem. Stochastic reformulations of quantum mechanics based on spontaneous collapses of the wavefunction are a popular approach to this issue. In this article, we derive the dynamic equations for the four most important spontaneous collapse models - Ghirardi-Rimini-Weber (GRW) theory, continuous spontaneous localization (CSL) model, Di\'osi-Penrose model, and dissipative GRW model - in the Wigner framework. The resulting master equations are approximated by Fokker-Planck equations. Moreover, we use the phase-space form of GRW theory to test, via molecular dynamics simulations, David Albert's suggestion that the stochasticity induced by spontaneous collapses is responsible for the emergence of thermodynamic irreversibility. The simulations show that, for initial conditions leading to anti-thermodynamic behavior in the classical case, GRW-type perturbations do not lead to thermodynamic behavior. Consequently, the GRW-based equilibration mechanism proposed by Albert is not observed.

Since they allow to describe quantum mechanics in a way that has strong formal analogies to classical mechanics (Wigner functions are defined on phase space and governed by a dynamic equation that reduces to the classical Liouville equation for the phase-space density in the limit → 0), an important application of Wigner functions is the study of the quantum-classical transition [30,31]. Wigner functions allow to develop a connection between classical and quantum models for liquid crystals [32] and to study mixed quantum-classical dynamics [33][34][35]. However, recovering the Liouville equation is not sufficient for recovering classical mechanics. Quantum systems can be in superpositions, whereas macroscopic * Corresponding author: raphael.wittkowski@uni-muenster.de superpositions (such as a cat being dead and alive at the same time [36]) are never observed. This is commonly attributed to a "collapse of the wavefunction" that takes place when a measurement is performed. However, this approach is often seen as unsatisfactory given that such a collapse is not allowed by the Schrödinger equation, meaning that one has to make the strange assumption that the Schrödinger equation does not apply to measurements. This issue is known as the measurement problem [37].
A popular approach to the measurement problem is to propose a stochastic modification of quantum mechanics in which collapses occur spontaneously [38][39][40][41]. Spontaneous collapse models make empirical predictions that differ slightly from those of "standard" quantum mechanics and can thus be tested experimentally [42,43]. Such experiments have become an important field of research in recent years [44][45][46]. Incorporating such spontaneous collapses in the theory of Wigner functions is desirable for at least four reasons: First, it allows to solve the problem of the classical-quantum transition in the Wigner framework also regarding the measurement problem. Second, if collapse models should be confirmed experimentally, the Wigner formalism should of course allow to describe them in order to remain a powerful alternative to Hilbert-space quantum mechanics. Third, Wigner functions are very useful for the visualization of collapse effects, as has been exploited for both spontaneous [47,48] and environmentinduced [49] collapses. Fourth, modern quantum technology -a field where Wigner functions are of high importance [2] -may be relevant for experimental tests of spontaneous collapse models. For example, it has been discussed whether SQUIDs (superconducting quantum interference devices) can be used for such a test [50,51].
In addition, spontaneous collapse models might provide a solution to another fundamental problem of physics, namely to the aforementioned problem of thermodynamic irreversibility [52,53]. This problem (which actually includes a variety of subproblems [53]) is concerned with the fact that macroscopic thermodynamics has a clear arrow of time associated with the monotonous increase of entropy, whereas the microscopic laws of (classical or quantum) mechanics are invariant under time reversal. David Albert [54][55][56] has suggested that spontaneous collapses of the wavefunction might be responsible for the irreversible approach to equilibrium. Although it has received significant attention in philosophy of physics [57][58][59][60][61][62], rigorous physical tests of this proposal are lacking at present.
In this work, we provide a reformulation of the four most important spontaneous collapse models, the Ghirardi-Rimini-Weber (GRW) theory [38], the continuous spontaneous localization (CSL) model [39,63], the Diósi-Penrose model [64][65][66][67], and the dissipative GRW model [68,69], based on Wigner functions. We then provide an analytical and numerical discussion of Albert's proposal, thereby investigating whether GRW theory might actually contribute to the solution of the problem of irreversibility. The mechanism suggested by Albert (deviations resulting from stochastic perturbations lead from anti-thermodynamic to thermodynamic behavior) is not observed in the simulations.
This article consists of two main parts. Although both contain important results, a reader only interested in one of them may skip the other one. First, we address the reformulation of collapse models in the Wigner framework: In Section II, we provide a brief introduction to Wigner functions. We then apply the Wigner transformation to the GRW theory (Section III), CSL model (Section IV), Diósi-Penrose model (Section V), and dissipative GRW theory (Section VI). Second, we discuss Albert's GRWbased explanation of irreversibility: The connection of GRW collapses to thermodynamic irreversibility is analyzed in Section VII. In Section VIII, decoherence-based alternatives are discussed. A numerical test of the GRWbased explanation of irreversibility is provided in Section IX. We conclude in Section X.

II. WIGNER FUNCTIONS
In "standard" quantum mechanics, a system is described by the density operator 1 (also known as "density matrix" or "statistical operator")ρ whose time evolution is given by the Liouville-von-Neumann equation with time t, imaginary unit i, reduced Planck constant , commutator [·, ·], and HamiltonianĤ. Furthermore, an observable is represented by a self-adjoint operatorÂ. An alternative to the density-operator-based formalism is to use Wigner functions [1]. Here, the state of the system is described by a function W (x, p) (Wigner function) that depends on the phase-space coordinates x and p. (In this work, we only require Wigner functions depending on position x and momentum p. However, the approach is much more general and can also be applied to systems with other degrees of freedom, such as spin [32,70].) The Wigner function is a generalization of a classical phase-space distribution function, from which it differs in that it can also take negative values [71]. For a one-particle system in one dimension, the Wigner function is defined as [72] W Throughout this work, the dependence of the Wigner function W and density operatorρ on time t is not written explicitly to improve readability. For N particles in three dimensions, Eq. (2) generalizes to [73] W ( r 1 , . . . , r N , p 1 , . . . , p N ) Equations (2) and (3) can not only be applied to the density operatorρ, but to an arbitrary operatorÂ. This gives its so-called Weyl symbol A(x, p). 2 The expectation value Â of the observableÂ, which in "standard" quantum mechanics is given by with the quantum-mechanical trace Tr, can then be calculated as [72] Â = dx dp A(x, p)W (x, p) just as in classical mechanics. We now discuss the dynamics of Wigner functions. First, we introduce the star product [74]: When transforming to phase space, one has to replace the product AB of two Hilbert space operatorsÂ andB by their star product A B, which is defined as [75] A (6) Here, the arrows indicate in which direction the derivatives act. The star product (6) allows to introduce the Moyal bracket [74,76] {A(x, p), B(x, p)} which is the Wigner equivalent of the commutator. When neglecting terms of order 2 , the Moyal bracket (7) reduces to the classical Poisson bracket. Combining Eqs.
(1) and (7) and taking into account the correspondence between commutators and Moyal brackets leads to the dynamic equation for the time evolution of Wigner functions where H is the Weyl symbol of the HamiltonianĤ. This Weyl symbol is given by H(x, p) = p 2 /(2m) + U (x) with the particle mass m and potential U [71]. Equation (8) reduces to the Liouville equation of classical mechanics if terms of order are neglected or if the potential is at most a second-order polynomial in x [72]. (Some difficulties of the limit → 0 are discussed in Ref. [71]. In particular, the classical limit cannot be taken in this simple way if the Hamiltonian generates chaotic classical dynamics [77].) An overview over theory and applications of Wigner functions is given in Refs. [2,4,5,73,78].

III. GRW THEORY
The formalism introduced in Section II is of interest for the Ghirardi-Rimini-Weber (GRW) theory [38], which is a modified form of quantum mechanics. It allows to deal with the famous quantum-mechanical measurement problem [37,[79][80][81][82]: If one performs a measurement on a quantum system that is in a superposition, then, assuming that both system and measurement apparatus obey the (linear, deterministic, and unitary) Schrödinger equation, the apparatus will, at the end, be in a superposition too. The quantum-mechanical time evolution inevitably leads to macroscopic superpositions, as is well known from the thought experiment on Schrödinger's cat [36]. However, this is not what we observe. A measurement on a system in a superposition produces a definite outcome. Therefore, one typically postulates a "collapse of the wavefunction", which is a sudden transition to an eigenstate of the observable one has measured as a consequence of measurement. This postulate, however, has no basis in the deterministic Schrödinger equation. A variety of solutions to this problem have been proposed, including hidden variables ("Bohmian mechanics" [83,84]), superdeterminism [85,86], modal interpretations [87], and the existence of many worlds in which all outcomes are realized ("Everett interpretation", also known as "many-worlds interpretation" [81,88]). It is sometimes claimed that decoherence effects (see Section VIII) provide a solution to the measurement problem. This, however, is wrong, since decoherence alone does not determine which one of two possible measurement outcomes occurs [79,89].
The GRW theory [38] (reviewed in Refs. [40,42,90]) takes a different approach by assuming that the collapse of the wavefunction is not a consequence of measurements, but something that happens spontaneously. For a system with N particles, the deterministic dynamics given by Eq. (1) is replaced by the more general master equation [38] where λ i is the localization rate for the i-th particle. In this work, we assume that λ i = λ with a constant λ for all particles. Since the localization rate is often assumed to depend on the particle mass [91,92], this can be motivated by the assumption that the particles belong to the same species.T i is the operator with the position operator for the i-th particlex i and a constant α determining the width of the collapsed function. It is easily shown that Eq. (9) is trace-preserving [38]. The idea of Eq. (9) is that "spontaneous collapse of the wavefunction" simply means that the particle's wavefunction is, at a certain rate λ, stochastically multiplied by a Gaussian. The parameters are tuned in such a way that these collapses are very rare for small systems (such as quantum particles), but very common for large systems (such as a measurement apparatus). This explains why superpositions are observed in microscopic, but not in macroscopic systems. The GRW theory (as well as other approaches to the measurement problem) is conceptually very important for a particular application of Wigner functions, namely the study of the classical limit. As discussed in Section II, the governing equation (8) for Wigner functions reduces to the classical Liouville equation for → 0. The measurement problem shows that recovering the Liouville equation is not sufficient to recover the observations of classical mechanics: If a cat is in a superposition of "being alive" and "being dead", the Liouville dynamics will not change this (a superposition of a dead and a living cat evolving classically remains a superposition of a dead and a living cat). What we require is a mechanism that explains why such superpositions are not observed in the macroscopic world. The GRW theory is one such mechanism.
A typical problem of spontaneous collapse models is that energy is not conserved. In fact, due to the noise, it increases at a constant rate. Although this effect is extremely small (the temperature of an ideal GRW gas increases by about 10 −15 K per year [38]), it is still an undesirable feature. To solve this problem, dissipative extensions of the GRW theory [68], CSL model [69], and Diósi-Penrose model [93] have been derived. They ensure that the energy remains finite during the time evolution of the system. These approaches are discussed in Section VI. Further extensions of GRW theory (not considered in this work) are models with colored rather than white noise for the stochastic perturbations of the wave-function [94] and relativistic extensions [95].
We now transform the governing equation (9) of the GRW theory to the Wigner formalism. We first demonstrate this derivation for the case of one particle in one spatial dimension to keep the notation compact, and then generalize to the case of N particles in three dimensions which is more relevant for statistical mechanics.
First, we use the fact that [38] x where |x is a position eigenstate with eigenvalue x. The operatorT (ρ) is subject to the same transformation rule (2) as the density operatorρ itself. Using Eq. (11), we find Consequently, the spontaneous multiplication by a Gaussian function of the position mediated by the operatorT (ρ) in the Hilbert space formalism corresponds to a convolution with a Gaussian function of the momentum in the Wigner function formalism.
We can now combine Eqs. (8), (9), and (12) to the GRW master equation for Wigner functions, given by To simplify Eq. (13), we make two approximations. First, we use a Kramers-Moyal expansion [96], which at second order leads to a Fokker-Planck equation [97]. For this purpose, we Taylor expand W (x, p − p ) around p = 0 and find where ! is the factorial and !! is the double factorial. When truncating at order n = 1, we obtain from Eq. (14) the Fokker-Planck equation with the momentum diffusion coefficient D p = λα 2 /4. Equations of the same form (but with a different physical content) have been obtained in the context of decoherence [98]. Here, one typically also has a friction term (see Section VIII), Eq. (15) is then obtained when taking the limit of vanishing friction at fixed diffusion constant [20]. Moreover, equations similar to Eq. (15) have been obtained for the collapse-induced [47] and decoherence-induced [99] dynamics of quantum rotors. Since Eq. (15) is of order 2 , the first term in Eq. (13) (the Moyal bracket) should also be expanded up to this order. This gives We have assumed here that H has the form H(x, p) = p 2 /(2m) + U (x) with a potential U . If U is quadratic, linear, or constant, e.g., if we have a free particle or a harmonic oscillator, the time evolution obtained from Eq. (16) is thus classical. Making this assumption, the Fokker-Planck equation (15) reduces to The same result is obtained if, e.g., due to a weak potential, the last term in Eq. (16) is small compared to the diffusion term in Eq. (17). Now, we consider N particles in three dimensions. A calculation analogous to Eq. (12) gives Hence, the extension of Eq. (13) to N particles in three dimensions is given by with i T i given by Eq. (18). The Fokker-Planck equation (17) generalizes to with D p = 3λα 2 /4. This derivation shows that the Wigner approach is very useful in the study of spontaneous collapse models: Using Hilbert space theory, the derivation of a Fokker-Planck equation for GRW theory (as done in Ref. [38]) requires several pages of calculation and a significant number of auxiliary assumptions. In the Wigner formalism, however, it is obtained almost immediately and in a very natural way. On the other hand, it is notable that Eq. (13) has the form of a convolution in momentum space, whereas in Eq. (9) the density operator is simply multiplied by a Gaussian in position space. For this reason, if we do not operate with a simplified (Fokker-Planck) model, it is easier to work with the Fourier-transformed Wigner functioñ that depends on two position coordinates rather than on a position and a momentum coordinate. A physical reason for why it is easier to work with the Fouriertransformed Wigner function (21) is that Wigner functions typically treat position and momentum on an equal footage, whereas in GRW theory collapses only take place in position space. For the Wigner function (21), Eq. (12) simplifies to The result (22) also directly follows from Eq. (11) together with the fact thatW is related to the density operator in position representation, (2)). The relation of W and ρ via the Fourier transformation [100] also explains why we need a convolution when describing the GRW dynamics in terms of Wigner functions (a Fourier transformation converts a multiplication to a convolution).

IV. CSL MODEL
A drawback of the GRW theory is that it does not preserve the symmetry properties of a system of identical particles [90]. In particular, it does not respect the fact that the wavefunction of a system of identical particles has to be symmetric under particle exchange for bosons and antisymmetric for fermions. This problem was solved by the continuous spontaneous localization (CSL) model introduced in Refs. [39,63]. We here follow Ref. [90]. First, we define the locally averaged density operator whereâ † ( r , s) andâ † ( r , s) are creation and annihilation operators for a particle with spin component s at position r and g is a symmetric positive function typically chosen to be The density operator is assumed to follow the dynamic equation with a constant ξ that is related to the localization rate, see Eq. (29) below, and the anticommutator [·, ·] + . From here on, we ignore the spin degrees of freedom since they are not relevant for the collapse dynamics. See Ref. [90] for a discussion of how Eq. (25) follows from the theory of stochastic processes in Hilbert space.
We now derive the Wigner representation of Eq. (25). In position representation, the last two terms on the righthand side of Eq. (25) read [90] Thus, the Weyl symbol for the operator given by the last two terms on the right-hand side of Eq. (25) is given by Therefore, the CSL master equation for Wigner functions reads Equations (28) and (19) coincide for N = 1 if we set [90] λ = ξ α 4π

V. QUANTUM GRAVITY
We now come back to the measurement problem. An interesting solution in the spirit of dynamical collapse models was proposed by Penrose [64,65] and Diósi [66,67], which is based on the problem of quantum gravity: As is well known, there is at present no generally accepted and experimentally verified unification of quantum mechanics and general relativity. While most approaches to this issue are based on quantizing gravity, an alternative approach consists in "gravitizing quantum mechanics" [64], i.e., in modifying quantum mechanics to incorporate effects of gravity. A spatial quantum superposition generates a superposition of spacetimes, which is supposed to lead to a gravity-induced collapse of the wavefunction [45]. Since this effect will depend on the mass of the system, it will be important for macroscopic, but negligible for microscopic systems; explaining why superpositions are not observed on macroscopic scales. Recent experimental tests [45], however, found no evidence of a gravity-induced collapse of the wave function, putting significant constraints on the possible validity of the Diósi-Penrose model. A simple quantitative model for this effect is given by [45,66,67] with the gravitational constant G, the total mass density operatorM , and the Euclidean norm · . This operator is given by [45] with the mass-density operatorμ i for the i-th particle. The Diósi-Penrose model constitutes a third type of collapse models in addition to GRW theory and CSL model. A difficulty arises from the fact that the standard definitionμ with the mass m i of and the position operatorˆ r i for the i-th particle leads to divergences [45]. To solve this problem, one has to smear out the mass density. This can be done in various ways [93]. Proposals from the literature include [67] with a length scale R 0 and the Heaviside function Θ, as well as [101] µ We are only interested in the non-Hamiltonian term on the right-hand side of Eq. (30). Moreover, we restrict ourselves to a single particle with mass density operatorμ . As shown in the supplementary material of Ref. [45], the non-Hamiltonian term in Eq. (30) can be written as with the functionΓ depending on the Fourier-transformed mass densityμ. The matrix element of Eq. (35) is given by [45] Hence, the Weyl symbol of the expression (35) is given by Using Eq. (38), the Diósi-Penrose master equation for Wigner functions is given by When comparing Eq. (39) with Eq. (13), it is easily seen that the general structure is the same, in particular if we assume that the mass-density distribution is Gaussian. However, the expression that the Wigner function is convoluted with differs due to the presence of the factor 1/( p ) 2 in Eq. (36). Equation (39) is an important result, since it has recently been suggested that Wigner functions could be helpful for experimental tests of quantum gravity [102].

VI. DISSIPATIVE GRW MODEL
Finally, we turn to a more recent extension of the GRW theory, the dissipative GRW model proposed by Smirne et al. [68]. As discussed in Section III, the standard spontaneous collapse models have the problem that the energy increases continuously. This increase is not measurable on typical observational timescales, but leads to a divergence of the temperature for t → ∞. To avoid this problem, the dissipative GRW model includes dissipation to ensure that the system thermalizes to a finite temperature. A similar extension has been developed for the CSL model by Smirne and Bassi [69] and for the Diósi-Penrose model by Bahrami et al. [93].
Since the derivation of the Wigner representation is more involved for this theory, we restrict the discussion to a single particle in one dimension to improve readability (the extension is straightforward). The dissipative GRW model is given by the master equation [68] ∂ which looks exactly like Eq. (9). However, the standard collapse operatorT is replaced by the dissipative collapse operatorT with a localization length scale r c = 1/ √ α and a temperature parameter k. The operator (41) reduces to the standard collapse operator (10) (and the dissipative GRW model (40) to the ordinary GRW model (9)) for k → 0, which corresponds to the high-temperature limit.
In the momentum representation, Eq. (41) gives where |p and |q are momentum eigenstates with eigenvalues p and q, respectively. We now use the fact that the Wigner function (2) can also be written as [32] W (x, p) = 1 2π dp p − Thus, the Weyl symbol ofT d is given by dp dp e − r 2 Hence, the dissipative GRW master equation for Wigner functions reads with T d given by Eq. (44).
Since Eq. (45) with the Weyl symbol (44) is extremely complicated, we now derive a simplified form. First, we note that Eq. (44) can be written as dp dp e − r 2 (Equation (46) is simply the first line of Eq. (44).) Motivated by the fact that the standard (nondissipative) GRW model is recovered for k = 0, we now Taylor expand Eq. (46) to first order in k. This gives The result (47) is already reminiscent of Eq. (12) from the standard GRW model. In analogy to the treatment in Section III, we now Taylor expand W (x, p − p ) around p = p in Eq. (47) and find where we have dropped terms of order k 2 and products of k and 2 in the last step. After inserting Eq. (48) into Eq. (45), we finally obtain the Fokker-Planck equation with damping parameter γ = 2λk, Boltzmann constant k B , and noise temperature T n = 2 /(8mkk B r 2 c ). This is precisely the temperature a dissipative GRW system thermalizes to [103], which confirms the consistency of our approach. If we ignore the fact that Eq. (49) contains a Moyal rather than a Poisson bracket, it is formally identical to the Kramers equation [104] which describes the dynamics of an underdamped Brownian particle in contact with a heat bath T b . In Eq. (49), however, the damping and diffusion terms arise not from the interaction with a solvent, but from the fundamental equations of motion of (GRW-type) quantum mechanics, implying that they are also present in a closed system.

VII. IRREVERSIBILITY AND APPROACH TO EQUILIBRIUM
David Albert [54][55][56] has suggested that the GRW theory might be useful in explaining the approach to thermodynamic equilibrium. A typical explanation for the approach to thermodynamic equilibrium runs as follows: Macroscopic states of a system are associated with a large number of possible microscopic states. The equilibrium state is the macrostate with the largest entropy, i.e., the state with the largest number of possible microscopic realizations. Consequently, most of the microstates asso-ciated with an initial macrostate will evolve towards an equilibrium state, and since we do not know the exact microstate of the system, we have to conclude that it is considerably more likely that the system will evolve towards equilibrium rather than away from it. This explanation, Albert argues, is unsatisfactory because it involves an appeal to the limitations of observers (ignorance about the precise initial microstate) in explaining what appears to be an observer-independent phenomenon, namely the approach to equilibrium. GRW theory provides a possible solution by introducing objective probabilities associated with the spontaneous col-lapse of the wavefunction. If the system is initially in an "abnormal" microstate that would lead to an evolution away from equilibrium -which is unlikely, but not impossible -then a GRW collapse will in most cases lead to a "normal" state from which the system equilibrates. According to Albert, this happens since every microscopic neighborhood of the abnormal microstate in phase space contains significantly more normal than abnormal microstates, such that thermodynamic abnormality is unstable under perturbations. Some discussions of Albert's proposal can be found in Refs. [57][58][59]. See Ref. [62] for a review.
Albert's approach belongs to the so-called Boltzmann approach to statistical mechanics [58]. In this framework, one is interested in a single system that has a certain macrostate (characterized by its macroscopic properties). This macrostate depends on the microscopic state of the system. Various microstates can correspond to the same macrostate. The equilibrium state is then the macrostate that the largest number of microstates correspond to. An alternative framework, which is more widely used in practical applications, is the Gibbs approach, where one studies ensembles (hypothetical sets of infinitely many copies of a system with different initial conditions). "Approach to equilibrium" then means that the probability distribution over an ensemble of systems converges towards the equilibrium (in the simplest case: Maxwell-Boltzmann) distribution. A more detailed discussion and comparison of both approaches can be found in Ref. [105].
Albert's proposal is conceptually appealing for at least three reasons [58]: First, it allows to solve two problems -the quantum measurement problem and the problem of thermodynamic irreversibility -at once. Second, it unifies quantum-mechanical probabilities with the probabilities used in statistical mechanics. Third, it allows to explain thermodynamic irreversibility without requiring assumptions about a probability distribution over initial conditions. However, it is only based on a qualitative argument, a quantitative test is lacking at present. Such a test would be required to see how far GRW collapses can actually get us on the road to equilibrium [62].
A simple argument for Albert's proposal can be developed based on Eq. (20). As is well known, the Boltzmann equation [106] provides a useful model for the approach to thermodynamic equilibrium. The one-particle distribution function f ( p), giving the probability that a particle has momentum p, approaches the Maxwell-Boltzmann (equilibrium) distribution due to the presence of collisions. This is a consequence of the fact that, in the derivation of the Boltzmann equation, one assumes that velocities are uncorrelated, i.e., that the two-particle distribution satisfies f 2 ( p 1 , p 2 ) = f ( p 1 )f ( p 2 ) (molecular chaos), before a collision [107]. If one instead assumes that velocities are uncorrelated after a collision, we would have found that the entropy always decreases.
Thus, if Eq. (20) should explain the approach to thermodynamic equilibrium, it should enforce f 2 ( p 1 , p 2 ) = f ( p 1 )f ( p 2 ). There is a good reason to believe that this is the case: The distributions f and f 2 can be obtained by integrating the Wigner function 3 over all positional and over all but one (for f ) or two (for f 2 ) momentum degrees of freedom [21]. Ignoring the reversible Hamiltonian part of the dynamics, we find that the difference between a factorized and a nonfactorized distribution evolves as Obviously, f 2 ( p 1 , p 2 ) = f ( p 1 )f ( p 2 ) is a fixed point of Eq. (51). From the general properties of the diffusion equation (which leads to a homogeneous distribution), we can infer that the fixed point is stable. Consequently, the GRW dynamics indeed tends to destroy correlations in momentum space and thereby leads to molecular chaos. There is also a different aspect that has to be taken into account: When comparing Eq. (20) to the Kramers equation (50) we can note that Eq. (20) has no friction term. However, the friction term in Eq. (50) is known to be essential for the approach to an equilibrium state (corresponding, in the classical case, to the Maxwell-Boltzmann distribution), and is intimately connected to the fluctuations described by the diffusion term in Eq. (50)  This brings us back to the distinction between Boltzmannian and Gibbsian approaches to equilibrium introduced above. If W can be interpreted as a probability distribution (i.e., in the classical case), the Kramers equation (50) describes an approach towards the Maxwell-Boltzmann distribution, which is the equilibrium distribution of Gibbsian statistical mechanics. For this, the friction term plays an essential role. In Eq. (20), in contrast, a friction term is not present, such that we cannot expect that W converges to a Maxwell-Boltzmann distribution. This, however, does not exclude that it leads (with high probability) to a convergence towards the macrostate with the largest phase-space volume, i.e., to Boltzmannian equilibrium.
Up to now, we have only considered the original GRW model by Ghirardi et al. [38]. This is the theory that Albert's original suggestion is based on. 4 The drawbacks of the standard GRW model in explaining equilibration can, however, be avoided by considering the dissipative GRW model introduced in Section VI. As shown there, the Wigner function approximately follows the dynamic equation (49) in this theory. Since Eq. (49) does contain a damping term related to the diffusion term in the standard way (actually, Eq. (49) reduces to Eq. (50) for → 0 if we set T n = T b ), the dissipative GRW model should be perfectly capable of describing the approach to thermodynamic equilibrium. This suggests the following modification of Albert's proposal: The approach to thermodynamic equilibrium might be explained by spontaneous dissipative GRW collapses.
This would have interesting consequences: The Kramers equation (50), of which Eq. (49) is an extension, describes a system coupled to a heat bath at fixed temperature (rather than an isolated system with fixed energy) that should be described in the canonical ensemble. Since spontaneous GRW-type collapses occur also in isolated systems, this would imply that the equilibrium distribution of an isolated system is given by the canonical rather than by the microcanonical distribution [69,93]. This is a very significant difference to "standard" statistical mechanics, where isolated systems are described by a microcanonical distribution, which corresponds to the energy-conserving case. While the different ensembles agree in the thermodynamic limit, they can differ significantly for systems with a small number of particles [109]. A different result would be obtained in an energy-conserving spontaneous collapse model that might be derived in the future. (Note, however, that Albert's original discussion [54,55] is set up in a Boltzmannian framework that is not based on ensembles at all.) In addition, the reduction of the spontaneous heating by dissipation has implications for tests of spontaneous collapse models that rely on this effect, for example those based on measurements of neutron star heating [110].
The idea that the dissipative GRW model (49) is more appropriate as an explanation for thermodynamic irreversibility than the original GRW model is, however, solely based on the fact that it is formally analogous to Eq. (50), which is known to lead to equilibration. Physically, the dissipative GRW model essentially assumes that all particles in the universe are coupled to a universal heat bath with a certain temperature T n . (The nondissipative GRW model is the special case T n → ∞.) Consequently, in the dissipative GRW model, systems will (in the long run) approach a canonical equilibrium state with temperature T n . This is irreversibility, but not of the type known from classical thermodynamics. There, a typical manifestation of equilibration is that if two systems with initially different temperatures T 1 and T 2 are brought in thermal contact, heat will flow from the hotter to the colder system until they have reached equal temperatures [111]. This final temperature will depend on the initial temperatures and other properties of the systems. It is certainly not a prediction of classical thermodynamics that all isolated systems converge to the same equilibrium state characterized by a universal temperature T n . This, of course, does not necessarily imply that the dissipative GRW model is incorrect. It merely implies that it makes a different empirical prediction than standard thermodynamics, namely that there is a universal temperature that all systems, irrespective of their initial temperature and energy, will equilibrate to. This equilibration will take place on a very long timescale, such that this prediction does not contradict the observation of "normal" equilibration in the form of heat exchange on ordinary timescales. However, it was this equilibration that we were seeking an explanation for. In this regard, the dissipative GRW model will play the same role (if any) as the standard one, namely that of erasing correlations by providing dissipation.
These conclusions are not affected in a significant way by working (as Albert does) in a Boltzmannian rather than a Gibbsian framework. As discussed by te Vrugt [53], the distinction between Gibbsian and Boltzmannian statistical mechanics can be linked to the distinction between deterministic and stochastic forms of dynamical density functional theory (DDFT) [33]. DDFT is a theory for the nonequilibrium evolution of the one-body density in classical fluids. Deterministic DDFT describes ensembles, whereas stochastic DDFT a (possibly spatially averaged) single system [112]. We can thus get an idea of the Boltzmannian equilibration predicted by dissipative GRW models by considering the predictions of stochastic DDFT, namely an approach to the minimum state of the free energy functional (which has a different definition compared to deterministic DDFT [33]). Hence, according to dissipative GRW theory, the system will (in equilibrium) fluctuate around the minimum of a free energy functional 5 with temperature T n .
When it comes to solving the fundamental problems of explaining thermodynamic irreversibility, typical arguments against coarse-graining -namely that the equilibrium state of equilibrium statistical mechanics is reproduced only approximately [115] and on certain timescales -also apply to GRW-based explanations. Nevertheless, Albert's approach, if successful, would still explain why anti-thermodynamic behavior is never observed in macroscopic systems. For microscopic systems, it has recently 5 In Gibbsian statistical mechanics, it is well known that equilibrium systems coupled to a heat bath minimize a free energy functional. For the Boltzmannian case, this deserves some comment: As shown by Dean [113], a system of Langevin equations for the position of classical particles can be exactly rewritten as a Langevin equation for the microscopic one-body density. The latter Langevin equation can be written in terms of a (free energy) functional depending on the one-body density. Then, the corresponding Fokker-Planck equation has a stationary solution if one inserts the one-body density that leads to the minimum of this functional. As discussed in Section IX, the classical limit of dissipative GRW theory also leads to Langevin equations, with the difference to Dean's treatment being that we have to take also the momentum density into account. Dean's formalism has been extended to underdamped Langevin equations (for which the momentum density is relevant) by Nakamura and Yoshimori [114].
been demonstrated experimentally that heat can spontaneously flow from a cold to a hot system for initially correlated spins [116]. This is not in contradiction to the GRW-based approach, since GRW collapses are (by construction) extremely rare for small quantum systems and thus allow for such a behavior.

VIII. CONNECTION TO DECOHERENCE
Spontaneous collapse models like the GRW theory have a very close connection to the theory of quantum decoherence [117], which is reviewed in Refs. [77,79,98,118]. In fact, some authors refer to GRW-type collapses as "intrinsic decoherence" [119]. Moreover, there has been a significant amount of research on explaining the emergence of irreversibility based on the interaction of quantum systems with their environment [120], in particular decoherence [77,118,121]. Consequently, it is of interest to relate such approaches to Albert's proposal discussed here. The purpose of this section (which is a slight digression) is thus to present some other approaches from the literature that Albert's theory, on which we focus here, is related to. This allows to see our results in a more general context.
In contrast to the GRW theory, decoherence does by itself neither provide a solution to the measurement problem nor produce stochasticity. However, it is an important ingredient in various interpretations of quantum mechanics [81]. This has been exploited in an explanation of the approach to equilibrium proposed by Hemmo and Shenker [57][58][59]. Like Albert's idea, it is based on explaining probabilities in statistical mechanics based on quantum-mechanical probabilities. However, they take as a starting point not GRW theory, but no-collapse interpretations (such as modal interpretations or the Everett interpretation), in which a collapse of the wavefunction does not occur. Since the quantum dynamics of closed systems is isolated in the absence of collapses, stochasticity then requires external interventions. Thus, the proposal discussed in this section belongs to the tradition of interventionism [115,122,123], which explains the approach to equilibrium based on external perturbations.
We discuss the approach following Ref. [58]. In decoherence theory, one typically assumes that the initial state |Ψ of a quantum system and its environment can be written in the form with the state of the system |ψ , the state of the environment |E , and the tensor product ⊗. Moreover, one assumes that the interaction Hamiltonian describing the interaction of the system with its environment commutes with a certain observable, the so-called pointer variable.
A typical pointer variable is the position. Let {|ψ i } be a basis of the Hilbert space of the quantum system that consists of eigenstates of the pointer variable. The state of the system at time t can then be written as with the relative states |E i and the amplitudes ν i (t).
Due to the coupling with the environment, the relative states will satisfy E i |E j ≈ δ ij with the Kronecker delta δ ij after an extremely short time. Then, the reduced density operator (obtained by tracing over environmental degrees of freedom) of the system takes the form For a harmonic oscillator weakly coupled to an environment in equilibrium, the states |ψ i correspond to narrowly peaked Gaussians, such that the decohered system follows quasi-classical trajectories [124]. This result is, by itself, no solution of the quantum measurement problem since it leads to a superposition of different classical time evolutions. However, in nocollapse interpretations of quantum mechanics, such as modal interpretations [125], it is possible to interpret the reduced state (54) as a probability distribution over the different quasi-classical states {|ψ i }. During the time evolution, there will then be stochastic transitions between such states. These stochastic transitions can then play the same role in the approach to equilibrium as the stochastic GRW collapses in Albert's theory.
What is interesting about this approach is that the Wigner function of a decohering system coupled to a thermal environment can, in the simplest case, be shown to have the form [20,58,121,[126][127][128] which, for → 0, reduces to the Kramers equation (50). Thus, in this approach, the probability distribution also approaches a Maxwell-Boltzmann form in the classical case. Consequently, it is, like the dissipative GRW model, close in form to standard theories of equilibration. However, for these models, the temperature the system equilibrates to is that of the environment (rather than that of a universal heat bath). Determining which of these options is the correct one (if any) is beyond the scope of this work. In particular, it would obviously require an experimental test [44][45][46] of the (dissipative) GRW model. Moreover, it should be noted that there are important conceptual differences between these approaches -if we explain thermodynamic irreversibility by decoherence, then there can be no irreversibility in an isolated system. In contrast, such a behavior is to be expected from the GRW perspective, which would thus (in a sense) imply that the standard models of stochastic dynamics (such as the Fokker-Planck equation or the Langevin equation) are actually more fundamental than the ordinary Hamiltonian dynamics which they are usually thought to be an approximation to.
A further recent approach to the problem of explaining the emergence of statistical mechanics that is interesting in this context was proposed by Drossel [129][130][131][132]. She argues that the practice of statistical mechanics, including (but not limited to) the presence of irreversibility, is incompatible with the universal validity of deterministic unitary quantum mechanics, and has to be explained by the presence of fundamental stochasticity. Her view can be classified as an intermediate position between that of Albert and that of Hemmo and Shenker. Like Albert (and unlike Hemmo and Shenker), she sees the origin of irreversibility in a stochastic quantum dynamics that violates the standard Schrödinger time evolution. However, like Hemmo and Shenker (and unlike Albert), she attributes this stochasticity to interactions of a quantum system with the environment. Finally, Wallace [29] has argued, within an Everettian framework, that the probabilities of statistical mechanics can be understood as arising from quantum-mechanical probabilities based on the fact that the classical limit of quantum mechanics is approximately isomorphic to a theory of classical probability distributions (as is evident from the Wigner function picture).

IX. NUMERICAL EXPERIMENTS
To test Albert's proposal, we have performed simulations of a many-particle quantum system with and without GRW collapses. A first aspect to take into account here is the fact that our dynamics should allow for equilibration. This means that it is not possible to use an ideal gas, where "ideal" means that there are no interactions, not even collisions. Such an ideal gas is not guaranteed to approach equilibrium, since all particles will just move on a straight line without affecting each other. As soon as interactions are present, it is not possible to describe the system using a one-particle Wigner function without a closure for the interaction part of the dynamics. However, such closures often already introduce prior assumptions regarding the approach to equilibrium (such as the time-reversal symmetry breaking associated with the hypothesis of molecular chaos, see Section VII), which we want to avoid here. Thus, we would have to solve the dynamic equation (19) or (20) numerically for a manyparticle system over a significant timescale, which is impossible.
We can simplify the problem by noting that we are mainly interested in the effect of GRW collapses on the approach to equilibrium, and not in quantum effects such as superpositions. Using the Wigner framework, we have been able to derive the Fokker-Planck equation (20) which, although describing a quantum system, is very similar in form to classical Fokker-Planck equations. Thus, we can assume that we work in the classical limit where quantum effects apart from collapses (which, in the Fokker-Planck approximation, correspond to momentum diffusion) are negligible, such that particles have welldefined trajectories that can be described using Langevin equations. 6 The Langevin equations corresponding to the Fokker-Planck equation (20) read where χ i (t) is white noise with the properties with the ensemble average · E , dyadic product ⊗, and identity matrix 1. Thus, we can base our simulations on the Langevin equations (56) and (57). These are numerically much easier to handle than the many-particle Wigner function and thus allow us to model a large system. Moreover, this approach is in the spirit of Albert's discussion, which is based on the idea that the role of quantum collapses is to provide random perturbations, here represented by the noise term in Eq. (57). If we use a dissipative model such as Eq. (49) instead, Eq. (57) has to be replaced by and the diffusion constant D p in Eq. (59) is given by D p = γmk B T n . This then recovers the standard equations of motion for an underdamped Brownian particle, from which the GRW-based model given by Eqs. (56) and (57) can be obtained by taking the limit γ → 0 at fixed D p . We use a smoothed Lennard-Jones potential for the particle interactions (see Appendix A). For our simulations, Eqs. (56) and (57) are nondimensionalized (see Appendix B) and then solved using the Leapfrog algorithm with time step size dt = 0.0025.
To see how our simulations should be set up, we take into account what precisely Albert's suggestion is based on: Although, if the initial state of the system is chosen randomly, it is much more likely that the system evolves towards equilibrium, there are certain ("abnormal") initial conditions for which the system does not evolve towards equilibrium. It is then the (expected) effect of the Normal initial state Abnormal initial state  GRW collapses that the system is brought back onto the right track. To see whether this actually works, we first require such an abnormal initial condition. For this purpose, we consider a many-particle system whose time evolution is governed by Eq. (8) (i.e., by unitary quantum mechanics without spontaneous collapses). Starting from an initial nonequilibrium distribution at t = 0, we study its time evolution by solving the unitary dynamics (in its Newtonian approximation given by Eqs. (56) and (57) with χ = 0) numerically. After a certain time t rev , the system has evolved towards a distribution that is relatively close to an equilibrium state. We then stop and use the distribution at time t = t rev with reversed momenta as the initial condition for a second simulation. By the timereversal invariance of quantum and classical mechanics, the system should then evolve back towards the initial distribution. Since it has thereby moved from a close-to-equilibrium to a far-from-equilibrium state, we have thus found an "abnormal" initial condition (distribution at t = t rev with reversed momenta) with the property that the system spontaneously moves away from equilibrium.
Next, we perform another simulation with the same abnormal initial condition, but this time using GRW theory. As is obvious from Eq. (20), the dynamics with spontaneous collapse is no longer time-reversal invariant (the left-hand side and the first two terms on the right-hand side change signs under (t, p) → (−t, −p), the third term on the right-hand side does not). If spontaneous collapses can enforce the validity of the second law of thermodynamics in the way suggested by Albert, then the system should be observed to evolve towards equilibrium rather than away from it this time.
The results are shown in Fig. 1, visualizing the local temperature measured by the kinetic energy density e kin (see Appendix B, all variables are dimensionless). We consider a dense fluid (the dimensionless 7 number density is = 0.7), which makes the simulations more efficient and therefore allows us to see all relevant effects. For the first simulation, we prepare two two-dimensional systems of the same size and let them equilibrate at different temperatures T 1 = 0.5335 and T 2 = 0.6391. 8 Afterwards, the systems are brought in thermal contact. This reference state, which has a sharp temperature gradient, is shown in Fig. 1a. The total system contains about 8.5 · 10 6 particles. We consider a quadratic domain with boundary length L = 3461.75. 9 If we evolve the system forward in time following standard Hamiltonian dynamics (Eqs. (56) and (57) with χ(t) = 0) up to time t rev = 625, 10 we observe -in agreement with what one would expect from thermodynamics -an equilibration towards a state with a more homogeneous temperature distribution. This relaxed state is shown in Fig. 1b. A more detailed analysis of this classical relaxation process can be found in Ref. [52].
For the second simulation, we use as an initial condition the final state of the first one (shown in Fig. 1b), but with reversed particle momenta. This creates an "abnormal" initial condition that leads to anti-thermodynamic behavior under the Hamiltonian time evolution. The simulations confirm this: if we evolve this state forward in time using Hamilton's equations, we see that the system evolves from a state with a relatively homogeneous temperature (Fig. 1b) towards a state with a sharp temperature gradient shown in Fig. 1c, where the two subsystems have regained their initial temperatures. Consequently, the system has moved away from thermal equilibrium. This is surprising from a thermodynamic perspective, but should be expected from the fact that the underlying Hamiltonian microdynamics is time-reversible. Moreover, effects of this type are known from spin systems [115,116,133].
Finally, Albert's suggestion is tested via the third simulation, which has the same initial condition as the second simulation, but uses the GRW-based model given by Eqs. (56) and (57) with nonvanishing noise χ(t) (see 7 From the dimensionless number density , the dimensional number density can be obtained by multiplication with σ 2 , where σ is a length scale (see Appendix B). 8 We have set the temperatures in such a way that (a) the system behaves as a liquid, and (b) the two temperatures are different but not too far from each other (to avoid strongly nonlinear effects). The values of T 1 and T 2 are then obtained as the long-time averages in the equilibration simulations. These simulations had as an initial condition a homogeneous density with some noise and a Rayleigh velocity magnitude distribution. The standard deviation for the Gaussians used for the initial velocity distributions were 0.75 and 0.85, respectively. 9 We have made the system as large as reasonably possible given the computational constraints. The precise particle number was then N = 2 23 = 8388608, from which the domain length L is obtained as L = N/ ≈ 3461.75. 10 This time was chosen based on the available computing time. It is sufficiently long to achieve a reasonable degree of (but not full) equilibration, which is what we require for the problem at hand.
Appendix B for the noise amplitude). If his suggestion is correct, then the stochastic perturbations induced by the GRW dynamics should lead to "normal" thermodynamic behavior despite the abnormal initial condition. As can be seen from Fig. 1d, which shows the final state of the third simulation, this is not the case: the system moves away from thermal equilibrium towards a state with inhomogeneous temperature even in the presence of stochastic GRW-type perturbations. We have also performed a fourth simulation with the same initial condition as the second and third simulations, but using the dissipative GRW model given by Eqs. (56) and (60). The final state of this simulation is shown in Fig. 1e, revealing that damping also does not lead to a restoration of thermodynamic behavior. Figure 1f, which compares the temperature profiles for the final states of the second (green), third (orange), and fourth (red) simulation to the initial state of the first simulation (blue), shows that the system evolves back to the initial state (with a small difference due to unavoidable numerical errors) regardless of whether or not stochastic GRW-type perturbations are present. Consequently, our results refute Albert's proposal according to which, for anti-thermodynamic initial conditions, thermodynamic behavior can be recovered via stochastic GRW-type perturbations.
The numerical experiments presented here are similar to the ones by Orban and Bellemans [134], where anti-thermodynamic behavior was also achieved by timereversing a molecular dynamics simulation. Orban and Bellemans found that the initial state is not completely recovered due to numerical round-off errors, from which they inferred that anti-thermodynamic behavior is unstable (similar in spirit to Albert's suggestion). Such "numerical irreversibility" has subsequently been investigated also by other authors (some of whom even linked it to quantum effects) [135,136]. Since the effects of the round-off errors are now well understood, we can test what happens if GRW noise is added to the dynamics. Figure 1f shows that the difference between the reference state and the result of the deterministic time-reversed simulation (resulting purely from numerical errors) is larger than the difference between the deterministic and the stochastic simulations, even though the amplitude of the GRW noise is much larger than that of the "numerical noise" (such that our "null result" is not a consequence of numerical errors). The temperature gradient in the final state of the time-reversed simulations is smoother than in the reference state, suggesting that the way back from equilibrium has not been completed. However, this increased smoothness is present for both the deterministic and the stochastic simulations. This implies that the GRW noise has no significant effect on whether or not there is anti-thermodynamic behavior, and that the small differences between the initial state of the first simulation and the final states of the second and third simulations are a consequence of round-off errors.
A similar conclusion can be obtained from Fig. 2, which shows the time dependence of the absolute values of the , and fourth (damped stochastic) simulation. For the second, third, and fourth simulation, time is reversed in this plot to allow for a comparison to the first simulation. The difference between deterministic and stochastic time-reversed dynamics (and between the different forms of stochastic dynamics) is larger for larger wavenumbers, but generally smaller than between the reference and the time-reversed simulations.
Fourier modes of the kinetic energy density |ẽ kin | for the wavenumbers (a) k = 2π/L, (b) k = 8π/L, and (c) k = 28π/L for all three simulations (see Appendix B, all variables are dimensionless). The second, third, and fourth simulations are plotted with reversed time to allow for a comparison to the first one. For all three wavenumbers, the agreement between the two time-reversed simulations is better than their agreement with the first (reference) simulation, implying that the GRW collapses make no significant difference for the macroscopic (thermodynamic) behavior of the system. Therefore, it is unlikely that they are responsible for the emergence of macroscopic irreversibility. However, the collapses do make a certain difference: as can be seen from comparing Figs. 2a and 2b to Fig. 2c, the disagreement between deterministic and stochastic simulation results (i.e., the effect of GRW noise) becomes larger for increasing wavenumber k. Since larger wavenumbers correspond to smaller length scales, this means that there are notable differences between deterministic and GRW dynamics on microscopic scales. However, these differences become less important on the macroscopic (thermodynamic) level. A similar result is found when comparing the damped and the undamped GRW model: there are differences on smaller length scales (larger wavenumbers), but not on the macroscopic level.
At this point, we should discuss three possible objections to our result. First, it can be argued that the relevant GRW parameters (in particular D p ) are so small that observable differences should be expected for macroscopic systems (10 23 particles), but not for systems with only a few particles (Albert [54] explicitly notes that, if his explanation is correct, an absolutely isolated gas with about 10 5 particles would not be expected to show thermodynamic behavior). A macroscopic system would be far too large to be studied with microscopic simulations on reasonable timescales. To compensate for this fact, we have used a very strong noise. Assuming, for the sake of definiteness, that our fluid consists of argon atoms, our simulation parameters correspond to D p ≈ 2.554 · 10 −43 J 2 s/m 2 (see Appendix B), whereas using the parameter values λ = 10 −16 /s and α = 10 14 m −2 used in Ref. [38] (to get an idea for the order of magnitude of D p ) gives (in two dimensions) D p = λα 2 /2 ≈ 5.56 · 10 −71 J 2 s/m 2 . Adler [43] has suggested values of λ that are about 10 orders of magnitude larger, increasing the possible value to D p ≈ 10 −61 J 2 s/m 2 , which is still significantly below our value. If we assume that GRW effects become relevant if λN is of the order 1/s, then for D p ≈ 10 −61 J 2 s/m 2 (corresponding to λ ≈ 10 −6 /s) one would require about 10 6 particles to see GRW effects. Since we have about 10 7 particles and a noise that is about 9 orders of magnitude larger, there is no reason to assume that our "null result" for effects of GRW collapses on thermodynamics is a consequence of the system size. Current bounds for the GRW parameters are given in Ref. [44].
Second, we have tested a quantum-mechanical theory (GRW) using classical molecular dynamics simulations, motivated by the fact that solving the quantummechanical equations for a many-particle system is impossible. This can be justified using the Wigner formalism, which shows that if we ignore all quantum effects except for the spontaneous collapses, GRW theory is equivalent to a classical statistical theory for a system whose dynamics is governed by Langevin equations. Obviously, the result we would have obtained from a full quantum-mechanical calculation would be different. Nevertheless, if a full quantum-mechanical calculation would have led to a different result regarding the problem of irreversibility, then irreversibility would not (solely) be a consequence of GRW collapses, but (also) of the other quantum effects (such as the peculiarities of quantum chaos) that have been ignored in our simulations. Clearly, for a system that is in a macroscopic superposition, the GRW collapses will make an important difference for the dynamics (namely, they will destroy this superposition). However, for explaining why a system of particles that already behave quasi-classically approaches thermodynamic equilibrium, they appear to play no significant role.
Third, we have made the approximation of moving from a master to a Fokker-Planck equation (e.g., from Eq. (13) to Eq. (17)) in the derivation of the Langevin equations (56) and (57). This corresponds to the assumption that the Gaussian function that the Wigner function is convoluted with in Eq. (13) decays very rapidly for larger values of |p |. Since the full width at half maximum of this Gaussian is 2 ln(2)α 2 ≈ 1.756 · 10 −27 kg m/s (which is extremely small), this assumption should also be unproblematic.

X. CONCLUSION
This article has two central results. First, we have derived master and Fokker-Planck equations for the Wigner function based on the GRW theory, the CSL model, the Diósi-Penrose model, and the dissipative GRW model. This provides a self-consistent dynamical theory in phase space that allows to explain the emergence of classicality from quantum mechanics regarding both the emergence of Liouville dynamics and the solution of the quantum measurement problem. Moreover, it makes the advantages of the Wigner framework available to researchers interested in spontaneous collapse models. Second, we have used Langevin equations derived from the GRW Fokker-Planck equation for a numerical test of Albert's suggestion that GRW collapses might be responsible for the approach to thermodynamic equilibrium. The simulations reveal that stochastic GRW-type perturbations do not lead to thermodynamic behavior if it is not already present in the deterministic dynamics. Consequently, our results do not support Albert's idea.

Appendix B: Dimensionless quantities
For the nondimensionalization of our equations, we first write Eqs. (56) and (60) as where A is the noise amplitude and ξ i (t) has zero mean and correlation ξ i (t) ⊗ ξ j (t ) E = 1δ ij δ(t − t ). We then rescale r(t) = σ r nd (t/τ ), t = τ t nd , U ({ r j }) = U nd ({ r j /σ}), ξ i (t) = ξ i,nd (t/τ )/ √ τ , γ = γ nd /τ , and A = mσA nd /τ As discussed in Section IX, it is useful to compare the strength of the noise used in our simulations to the values expected from current estimates of the GRW model parameters. For the purposes of testing Albert's theory, it is good if the noise is "too large", since we can thereby ensure that we are not overlooking any effects due to our "small" system size. An estimate of the order of magnitude is sufficient. The value of D p in SI units is given by In the Lennard-Jones potential, the parameters σ and measure the range and depth of the potential, respectively. To get concrete numbers, we take the example of argon, for which the parameters σ and take the values σ = 3.4 · 10 −10 m and ≈ 1.657 · 10 −21 J, respectively [138]. Moreover, the mass of an argon atom is m ≈ 6.634 · 10 −26 kg [139], such that τ ≈ 2.151 · 10 −12 s. From Eq. (B3) and A nd = 10 −4 (value used in the simulations 11 ), we then find D p ≈ 2.554 · 10 −43 J 2 s/m 2 .
11 To achieve a good description of GRW effects in approximately For the simulation with damping, we have used a noise amplitude A nd = 10 −4 (as in the undamped case) and a noise temperature T n,nd = 0.5863 related to the noise amplitude via the fluctuation-dissipation relation A 2 nd = 2γ nd T n,nd , which determines the dimensionless damping parameter γ nd . For the example of argon, this corresponds to a physical value T n = T n,nd /k B = 70.356 K of the noise temperature. 12 The dimensionless temperature T nd (we here use T 1 = 0.5335 and T 2 = 0.6391) is defined as with the dimensional temperature T , number of spatial dimensions d (we use d = 2 for the simulations), and dimensionless velocity v i,nd of the i-th particle. For argon, the temperatures used here correspond to T 1 /k B = 64.02 K and T 2 /k B = 76.692 K. 13 The local temperature is measured by the local kinetic energy density Moreover, the y-averaged local kinetic energy density shown in Fig. 1d is given by e kin y = 1 L dy nd e kin ( r nd ). For the figures, the local kinetic energy density was calculated via a mass redistribution technique 14 .
classical fluids, we have chosen the noise amplitude in such a way that it is significantly larger than the numerical round-off error, but also significantly smaller than the typical order of magnitude of the force resulting from the particle interactions. 12 To avoid over-or undercooling the system, the dimensionless noise temperature T n,nd was simply chosen as the average of T 1 and T 2 . This specific choice avoids the problem that (as discussed in Section VII) the presence of the noise field will in general lead to a temperature different from the thermodynamically expected one (an effect that is also interesting, but different from the one we wish to investigate here). It is difficult to give an estimate of the real value of Tn (assuming that there is such a real value, i.e., that the dissipative GRW model is correct). Assuming that the noise has a cosmological origin, a value of about 1 K is reasonable [69]. Given that the damping parameter γ is proportional to The Fourier coefficients shown in Fig. 2 are defined as where k = (2π/L) n with n ∈ N 2 . The definition (B6) ensures thatẽ kin ( 0) = T nd for d = 2 (see Eq. (B4)). In Fig. 2, we always set k y = 0, such that k = k x = k .