Next Article in Journal
Quantum-Inspired Evolutionary Approach for the Quadratic Assignment Problem
Next Article in Special Issue
Agent Inaccessibility as a Fundamental Principle in Quantum Mechanics: Objective Unpredictability and Formal Uncomputability
Previous Article in Journal
From Physics to Bioengineering: Microbial Cultivation Process Design and Feeding Rate Control Based on Relative Entropy Using Nuisance Time
Previous Article in Special Issue
Nonlocality in Bell’s Theorem, in Bohm’s Theory, and in Many Interacting Worlds Theorising
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bouncing Oil Droplets, de Broglie’s Quantum Thermostat, and Convergence to Equilibrium

1
Aix Marseille Université, CNRS, Centrale Marseille, Institut Fresnel UMR 7249, 13013 Marseille, France
2
Graduate School of Mathematical Sciences, the University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8914, Japan
3
Centro Brasileiro de Pesquisas Físicas, Rua Dr. Xavier Sigaud 150, 22290-180 Rio de Janeiro, RJ, Brazil
*
Author to whom correspondence should be addressed.
Entropy 2018, 20(10), 780; https://doi.org/10.3390/e20100780
Submission received: 24 August 2018 / Revised: 1 October 2018 / Accepted: 3 October 2018 / Published: 11 October 2018
(This article belongs to the Special Issue Emergent Quantum Mechanics – David Bohm Centennial Perspectives)

Abstract

:
Recently, the properties of bouncing oil droplets, also known as “walkers,” have attracted much attention because they are thought to offer a gateway to a better understanding of quantum behavior. They indeed constitute a macroscopic realization of wave-particle duality, in the sense that their trajectories are guided by a self-generated surrounding wave. The aim of this paper is to try to describe walker phenomenology in terms of de Broglie–Bohm dynamics and of a stochastic version thereof. In particular, we first study how a stochastic modification of the de Broglie pilot-wave theory, à la Nelson, affects the process of relaxation to quantum equilibrium, and we prove an H-theorem for the relaxation to quantum equilibrium under Nelson-type dynamics. We then compare the onset of equilibrium in the stochastic and the de Broglie–Bohm approaches and we propose some simple experiments by which one can test the applicability of our theory to the context of bouncing oil droplets. Finally, we compare our theory to actual observations of walker behavior in a 2D harmonic potential well.

1. Introduction

“Walkers” are realized as oil droplets generated at the surface of a vibrating oil bath. As shown by Couder and Fort [1,2,3], the vibration of the bath prevents the coalescence of the droplets with the surface, allowing them to remain stable for very long times. Moreover, the trajectories of the walkers are guided by an external wave [4,5] that they themselves generate at the surface of the oil bath. From this point of view, walkers are reminiscent of wave-particle duality [2,6], and they seem to offer deep analogies with de Broglie–Bohm particles [7]. Up until now, different aspects of walker dynamics have been studied in a purely classical framework, typically in a hydrodynamical approach [3,5]. For instance, certain models address their deformations due to their bouncing off the surface of the bath, in function of the density and viscosity of the oil and other parameters [5]. Other studies describe the dynamics of the surface waves that the walkers generate during the bouncing process, and how those waves in turn guide their trajectories. In these models, this complex behavior is characterized by a memory time which relates the dynamics of the walker bouncing at time t to its successive bouncing positions in the past [8,9]. The presence of such a memory effect establishes a first difference with quantum mechanics. Normally, in quantum mechanics, it is assumed that all results of any possible future measurements to be performed on a quantum system are encapsulated in its present quantum state [10]: its wave function at the present time t.
Droplets also transcend the most common interpretations of quantum theory which prohibit any description of the system in terms of instantaneous, classical-like trajectories. Droplets and their trajectories are visible with the naked eye at any time and standard interpretations of quantum mechanics do not apply. This is why we believe that it is necessary and worthwhile to adapt realist (causal) formalisms such as de Broglie–Bohm (dBB) dynamics [11,12] or a stochastic version thereof à la Nelson [13] to explore the analogy with quantum systems. This is the main motivation of the present paper.
Another difference between walker trajectories and quantum trajectories is that the quantum description is intrinsically probabilistic and non-classical, while there exist regimes in which the trajectory of the walkers is deterministic and classical (for example, when they bounce exactly in phase with the bath, they can be shown to follow straight lines at constant velocity [14,15,16,17]). However, there also exist regimes in which a Brownian motion is superimposed on their flow lines (e.g., above the Faraday threshold), and other regimes where the trajectories appear to be chaotic [5]. In fact, in several regimes, droplets appear to exhibit ergodic behavior. In practice, ergodicity has been established on the basis of the following observations: if we prepare a walker at the surface of the liquid bath (a corral, for instance), it will progressively explore each part of the surface, following an apparently random motion [4]. If one then visualizes the statistics of the sojourn time of the walker in each of these regions, a striking pattern emerges, bearing more than a simple resemblance to an interference pattern [4,7]. It is this, again remarkable, manifestation of wave-particle duality that first attracted our attention and which lies at the origin of this paper. The onset of quantum equilibrium in the framework of dBB dynamics and in stochastic versions thereof is an important foundational issue in itself, which has motivated numerous studies (see, e.g., [13,18,19,20,21,22,23,24] as well as [25] and references therein). Several authors in the past have indeed tried to explain how the Born rule emerges from individual trajectories, which is a highly non-trivial problem. In the case of dBB dynamics, it is easy to show that in simple situations the relaxation to the Born statistical distribution does not occur at all, but recent studies [26,27,28,29,30,31] show that in sufficiently complex situations (several modes of different energies for instance) the system might exhibit mixing, which explains the onset of quantum equilibrium in such cases. As we shall show in the present paper, in the case of Nelson-type dynamics, the quantum Brownian motion imposed in such a model accelerates the relaxation to Born’s distribution, and in fact ensures that relaxation to the Born rule will almost always occur (as we shall also show). In our view, for the above reasons, de Broglie–Bohm and Nelson-type dynamics are good candidates for explaining how wavelike statistics emerge after averaging a set of apparently chaotic and/or stochastic trajectories.
Briefly summarized, our main goal is to explain the emergence of aforementioned interference patterns in the framework of the dynamical models of de Broglie–Bohm and of a stochastic version thereof which is based on the models of Bohm-Vigier [18] and Bohm-Hiley [19] but which is formally close to Nelson [13]. Both models are introduced in Section 2. Here, it is worth noting that thus far there is no experimental evidence that droplets indeed follow de Broglie–Bohm and/or Nelson trajectories. Our approach therefore differs radically from previous studies on droplets, in the sense that we impose a quantum dynamics by brute force, whereas, until now, the attempt to illustrate how chaos may underlie quantum stochasticity has been a pillar of the research on walkers/droplets. In fact, Nelson’s original goal, in proposing his dynamics, was to derive an effective wave equation from the properties of an underlying Brownian motion, as in classical statistical mechanics where a diffusion equation is derived from microscopic properties of the atoms. There actually exists an impressive number of attempts in that direction, as, e.g., stochastic electro-dynamics [5,32,33]. However, there exists (as far as we know) no way to derive an effective Schrödinger equation from hydrodynamical models of droplets.
By choosing exactly the opposite approach, i.e., by imposing quantum-like dynamics on the droplets, we pursue three goals. The first one is to describe the onset of quantum equilibrium (and ergodicity). A second objective is to formulate precise quantitative predictions regarding this relaxation process, which can possibly be validated by future experiments. A third objective is to show, for the first time, that certain dBB trajectories present a deep structural resemblance with certain trajectories that have been reported in the literature for droplets trapped in a harmonic potential.
A short discussion of the onset of equilibrium in de Broglie–Bohm dynamics and the importance of coarse-graining is given in Section 3. In the case of our stochastic, Nelson-type dynamics, we derive in Section 4 a new H-theorem showing the relaxation to quantum equilibrium, which does not rely on coarse-graining and is valid at all scales. We pay particular attention to the ergodicity of trajectories in the case of our stochastic dynamics (which mix properties of the de Broglie–Bohm dynamics with Brownian motion). We apply these ideas to discuss ergodicity in the case of the stochastic treatment of a particle trapped in a harmonic potential (Section 5) and to describe the dynamics of a droplet trapped in a harmonic potential (Section 6). In this latter section (in Section 6.1), we also propose some simple experiments by which one can test the applicability of a Nelson-type dynamics to the context of bouncing oil droplets, and we briefly discuss the problems caused by the presence of zeros in the interference pattern that is encoded in the statistics of the trajectories. In Section 7, we study a situation during which the attractor of the probability distribution is no longer a static eigenstate of the (static) Hamiltonian, and we compare the onset of equilibrium in the dBB and stochastic formalisms in that special framework. In Section 8, we tackle the dynamics of droplets in a 2D harmonic potential through a simple model where the pilot wave is treated as a dynamical object. This constitutes a preliminary attempt, ultimately aimed at establishing a dynamics that would combine stochastic and/or dBB dynamics with a feedback of the trajectory on the wave, which is a fundamental feature of droplet phenomenology that has never been addressed in the framework of dBB or Nelson dynamics. The last section is devoted to conclusions and open questions. A short overview of the numerical methods used in the paper is given in the Appendix A.

2. dBB and Nelson Dynamics

2.1. The dBB Theory

In the following quick overview of the dBB theory we shall limit ourselves to the case of a single particle. In the dBB theory, particle positions exist at all times and they are merely revealed by position measurements, instead of “originating” with the measurement as the standard interpretation of quantum mechanics would have it. The dynamics is described by a wave function which obeys the Schrödinger equation:
i Ψ ( x , t ) t = 2 2 m Δ Ψ ( x , t ) + V ( x , t ) Ψ ( x , t )
where V ( x , t ) is an external potential and m the mass of the particle, as well as by a position x . In order to reproduce the predictions of standard quantum mechanics, one must have that the positions are distributed according to
P ( x , t ) = | Ψ ( x , t ) | 2
where P ( x , t ) is the distribution of particle positions over an ensemble of trajectories. An ensemble satisfying condition (2) is said to be in quantum equilibrium.
It is also commonly assumed that (2) is satisfied at some initial time. Therefore, in order to be at (quantum) equilibrium for all t, the condition to enforce is
P ( x , t ) t = | Ψ ( x , t ) | 2 t .
As is well known, the probability density | Ψ ( x , t ) | 2 satisfies the continuity equation
| Ψ ( x , t ) | 2 t + · j ( x , t ) = 0
where
j = m Im Ψ * Ψ ,
is the (probability) current describing the flow of the probability due to (1).
The probability density P , on the other hand, will satisfy a continuity equation
P t + · P v = 0
where v is the velocity field for the particle. Therefore, Equation (3) will be satisfied if
v ( x , t ) = j ( x , t ) | Ψ ( x , t ) | 2 .
The expression (7) for the velocity field is of course not the only possible one: any solution of the form
v ( x , t ) = v ( x , t ) + × f ( x , t ) | Ψ ( x , t ) | 2 ,
where f is a scalar function, will also give rise to Equation (3) (see [34] for more details).
Secondly, if one expresses the wave function in terms of its phase S ( x , t ) and modulus R ( x , t ) = | Ψ ( x , t ) | 2 ,
Ψ ( x , t ) = R ( x , t ) e i S ( x , t ) / ,
one finds that
j = | Ψ ( x , t ) | 2 m S
and that the velocity of the particle at time t is given by
d x ( t ) d t = v ( x , t ) = 1 m S ( x , t ) | x = x ( t ) .
Integrating the system given by Equation (11), we recover the dBB trajectory. From the above, it should be clear that the dBB theory is deterministic. Any stochastic element only comes from our lack of knowledge of the initial positions.
In the context of bouncing droplets, we shall view the external wave generated by the droplet as being in one-to-one correspondence with the “pilot wave” Ψ , which guides the position of the dBB particle.

2.2. A Simple Realization of de Broglie’s Quantum Thermostat—Nelson Dynamics

As mentioned in the introduction, the trajectories of walkers are often characterized by a non-negligible stochastic (Brownian) component which sets them apart from the smooth dBB trajectories. From this point of view, it seems worthwhile to try to model walkers dynamics in terms of stochastic generalisations of dBB dynamics.
de Broglie himself, in fact, considered such generalizations of the deterministic dBB dynamics (which he called the “quantum thermostat hypothesis”) to be highly welcome because they might provide a physically sound picture of the hidden dynamics of static quantum states. For instance, if we consider the position of an electron prepared in the ground state of a hydrogen atom, the dBB dynamics predicts that its position will remain frozen at the same place throughout time, which is counterintuitive to say the least. Adding a stochastic component to its velocity could, in principle, explain why averaging the position of the electron over time is characterized by an exponentially decreasing probability density function, in agreement with the Born rule (provided, of course, that ergodicity is present in the problem in exactly the right proportion). A first proposal in this sense was formulated by Bohm and Vigier in 1954 [18] which, later on, was made more precise by Bohm and Hiley [19], but stochastic derivations of Schrödinger’s equation by Nelson [13] (and others [32,33] in the framework of stochastic electrodynamics) can also be considered to provide models of the quantum thermostat. Quoting de Broglie: “...Finally, the particle’s motion is the combination of a regular motion defined by the guidance formula, with a random motion of Brownian character... any particle, even isolated, has to be imagined as in continuous “energetic contact” with a hidden medium, which constitutes a concealed thermostat. This hypothesis was brought forward some fifteen years ago by Bohm and Vigier [18], who named this invisible thermostat the “subquantum medium”... If a hidden sub-quantum medium is assumed, knowledge of its nature would seem desirable...” (In [35] Ch.XI: On the necessary introduction of a random element in the double solution theory. The hidden thermostat and the Brownian motion of the particle in its wave.)
In this paper, we shall consider a particular model of the quantum thermostat in which, as in the Bohm–Vigier model, a single spinless particle suspended in a Madelung fluid moves with the local velocity of the resulting field, given by Equation (11), and is subjected to fluctuations coming from the latter (cf. Figure 1). However, following Nelson, we shall model these fluctuations by means of a particular stochastic process. To be precise: our model is formally the same as Nelson’s in that it relies on the same stochastic process. However, in spirit, it is closer to the Bohm–Hiley model [19] in that we do not assume to be at quantum equilibrium (an assumption which is fundamental to Nelson’s theory, as was already pointed out by Bohm and Hiley [19]; see also [20] for a detailed presentation and a comparison of both approaches).
This stochastic process is defined on a probabilistic space Ω , characterized by a probability distribution P ( x , t ) and obeying an Ito stochastic differential equation of the general form:
d x = 1 m S + γ x = x ( t ) d t + α d W ( t )
where α is the (constant) diffusion coefficient that characterizes the strength of the random part, and d W ( t ) is a Wiener process in three dimensions. The function γ ( x , t ) in (12) is a systematic drift, the so-called osmotic velocity, which we shall fix in the following way.
The conservation equation of the probability distribution (which we denote by P, in order to stress the difference with the probability in the dBB dynamics which is denoted by P ), can be written as the Fokker–Planck equation:
P t = · P m S + γ P + α 2 Δ P .
If we now require that the quantum equilibrium P ( q , t ) = | Ψ ( q , t ) | 2 be a solution of this Fokker–Planck equation, we obtain from Equations (4), (10), and (13) that
· γ | Ψ | 2 α 2 | Ψ | 2 = 0 ,
which is a constraint on the osmotic velocity. The simplest solution of this constraint is [36]
γ ( x , t ) = α 2 | Ψ | 2 | Ψ | 2 .
In the rest of the paper, we choose the osmotic drift velocity to be Equation (15), with α an a priori free parameter, Nelson’s choice for α ( α = / m ) being irrelevant if we apply this formalism to droplets.
In summary, our Nelson dynamics is fully defined by the following Ito equation
d x ( t ) = 1 m S + α 2 | Ψ | 2 | Ψ | 2 x = x ( t ) d t + α d W ( t )
where d W i ( t ) represents a Wiener process with
< d W i ( t ) > = 0 and < d W i ( t ) d W j ( t ) > = 1 2 δ i j δ ( t t ) ,
and by the Fokker–Planck equation for the associated probability distribution P ( x , t )
P t = α 2 Δ P · P m S + α 2 P | Ψ | 2 | Ψ | 2 ,
where Ψ ( x , t ) satisfies the Schrodinger equation:
i Ψ t = 2 2 m Δ Ψ + V Ψ .
At quantum equilibrium, i.e., when P ( x , t ) = | Ψ ( x , t ) | 2 , the diffusion velocity is balanced by the osmotic term and the Bohm velocity is recovered, on average.
We shall now discuss the details of the relaxation towards quantum equilibrium, in the dBB and stochastic formalisms.

3. Relaxation to Quantum Equilibrium in the de Broglie–Bohm Theory

In our presentation of the dBB theory for a single particle, in Section 2.1, we assumed that the particle positions are initially distributed according to Born’s law
P ( x , t i ) = | Ψ ( x , t i ) | 2
over an ensemble. The dynamics then ensure that the same relation will hold for any later time. This is the assumption de Broglie and Bohm made in their original papers [11,12,37]. Although Bohm tried, already in the 1950s (first on his own—see, e.g., [11] (Section 9)—and later with Vigier [18]), to relax this assumption by modifying the dynamics, to many authors working today on the dBB theory it is still an assumption which has to be made (the final objective of de Broglie, Bohm, Vigier, and Nelson–and many other contributors to various realistic hidden variable interpretations in which quantum systems are assumed to be localized in space at any time–was to rationalize wave-like statistics in terms of individual trajectories; the same problem occurs in droplets phenomenology and, according to us, admits no fully satisfying solution yet).
According to Valentini [21,38,39], however, there is no need to assume that the particle positions are initially distributed according to Born’s law or to modify the dynamics. His claim is that an ensemble in which Born’s law is not satisfied (so-called quantum non-equilibrium) will evolve naturally towards quantum equilibrium, provided that the wave function leads to sufficiently complex dynamics. This relaxation process has to take place on a coarse-grained level and can only occur if the initial distributions do not display any fine-grained micro structure.
Let us first explain the need for coarse-graining. Let us introduce the function f = P / | Ψ | 2 , as in [26]. An important implication of (6) is that the function f is conserved along the dBB trajectories:
d f ( x , t ) d t f ( x , t ) t + x ˙ · f ( x , t ) = 0 .
Hence we have that
P ( x , t ) = P ( x i , t i ) | Ψ ( x i , t i ) | 2 | Ψ ( x , t ) | 2
where x i is the initial position of the particle which leads to x , when evolving from t i to t according to the dBB dynamics. If one assumes that P ( x i , t i ) / | Ψ ( x i , t i ) | 2 1 , relaxation to quantum equilibrium is clearly impossible, at least at the microscopic level. However, as argued by Valentini [21], relaxation is possible at the coarse-grained level, provided the initial distribution does not display any fine-grained microstructure.
The operational definition of the coarse-graining is as follows. We divide the domain of interest A Ω into small cubes of equal edge length ϵ (we call them coarse-graining cells, or CG cells for short). These CG cells do not overlap and their union is equal to A. The coarse-grained densities, which we denote by P ¯ ( x , t ) and | Ψ ( x , t ) | 2 ¯ , are then defined as
P ¯ ( x , t ) = 1 ϵ 3 CG cell x d 3 x P ( x , t )
| Ψ ( x , t ) | 2 ¯ = 1 ϵ 3 CG cell x d 3 x | Ψ ( x , t ) | 2
where the domain of integration is the CG cell containing x .
We can now discuss the second assumption: the lack of a fine-grained microstructure in the initial distribution. Let us assume we have a non-equilibrium distribution P ( x , t i ) which relaxes to quantum equilibrium at the coarse-grained level, under the dynamics generated by the wave function Ψ ( x , t ) . As the dBB theory is time-reversal invariant, in the time-reversed situation, under the dynamics generated by Ψ * ( x , t ) , we would have a distribution that moves away from quantum equilibrium. Thus, it would seem that time-reversal invariance contradicts the possibility of relaxation to quantum equilibrium. This conclusion is unwarranted, however: as the initial distribution P ( x , t i ) relaxes to quantum equilibrium, it retains information on the original values of f (which are constant in time) and thereby acquires a fine-grained microstructure, which means that at the final time t f , P ( x , t f ) will differ significantly from P ¯ ( x , t f ) . Therefore, in the time-reversed situation, the initial distribution would exhibit a fine-grained microstructure, which is prohibited under our assumption, thereby breaking the time-reversal invariance.
In order to quantify the difference between the distribution P ( x , t ) and the quantum equilibrium condition | Ψ ( x , t ) | 2 at the coarse-grained level, Valentini [21,38,39] introduced the entropy-like function
H ¯ ( t ) = Ω d 3 x P ¯ ln P ¯ / | Ψ | 2 ¯
where P ¯ and | Ψ | 2 ¯ as in Equations (23) and (24), for which he has shown the (quantum) H-theorem,
H ¯ ( t ) H ¯ ( t i ) ,
under the assumption of no fine-grained microstructure. It should be stressed, however, that this is not necessarily a monotonic decay and therefore does not prove that quantum equilibrium will always be reached. It merely indicates a tendency towards relaxation. The strongest support for the idea of relaxation to quantum equilibrium comes from numerical simulations of the evolution of non-equilibrium distributions for various quantum systems [26,27,28,29,30,31] (see [25] and references therein for a review). The first numerical simulations were performed by Valentini and Westman [26] who showed, in the case of a 2D box, that relaxation quickly takes place for a wave function which is a superposition of the first 16 modes of energy (the superposition being equally weighted). It was also hinted that the nodes of the wave function, with their associated vorticity, play a crucial role in the relaxation process, as purveyors of chaos (or mixing) in the dynamics. This later claim was properly understood in [40]. The dependence of the relaxation timescale on the coarse-graining length ϵ and on the number of energy modes was studied in [28]. In [31], it was shown that quantum systems with a low number of modes are likely to never fully relax, in which case H ¯ reaches a non-zero residue value. However, such a scenario becomes unlikely as the number of modes increases.
According to the quantum non-equilibrium hypothesis, standard quantum mechanics is only one facet of the pilot-wave theory, that of quantum equilibrium, leaving the possibility for new physics: that of quantum non-equilibrium. One should assume of course that during our time we have only had (or can only have) access to systems for which quantum equilibrium has already been reached. But that does not mean that quantum non-equilibrium never existed in the early universe (which could be inferred from the observation of the remnants of the early fractions of seconds of the universe, just after the Big Bang [41]), or that some, yet undetetected, exotic quantum systems cannot still be in quantum non-equilibrium today [42]. This is why droplets are appealing, because if their dynamics do present analogies with dBB dynamics, their study will allow us to observe relaxation to (quantum) equilibrium “in real time” in the lab with the naked eye, which is not possible with quantum systems for which we have no direct access to individual trajectories.

4. An H-Theorem for Nelson Dynamics

Let us start by introducing an analog of Valentini’s entropy, Equation (25), for the probability distribution P ( x , t ) associated with our Nelson dynamics, as defined by Equations (16)–(19):
H V ( t ) = Ω d 3 x P ln P | Ψ | 2 ,
which is a special instance of a relative entropy known as the Kullback–Leibler divergence [43].
We also define a second non-negative functional
L f ( t ) = Ω d 3 x f ( P | Ψ | 2 )
where
f ( x , t ) = P ( x , t ) | Ψ ( x , t ) | 2 .
Note that we always impose the boundary conditions | Ψ | 2 | Ω = P | Ω = 0 and f | Ω = 1 so as to avoid divergence of these integrals on the boundary of Ω .
It should be noted that the entropy of Equation (27) or the functional of Equation (28) we shall use to quantify the relaxation to quantum equilibrium are very different from the entropies usually considered in the context of classical H-theorems (like e.g., the Boltzmann entropy). One should bear in mind, however, that quantum equilibrium is radically different from classical equilibrium [44] and has no connection whatsoever with relaxation to quantum thermal equilibrium, for the simple reason that the Born distribution of positions reached by an ensemble of trajectories à la Nelson or dBB is not a thermal distribution.
To understand why the functionals in Equations (27) and (28) are non-negative and why they are zero if and only if (quantum) equilibrium is reached (that is to say when f = 1 everywhere in space), it is important to note that the integrands of H V and L f satisfy the inequalities
( P | Ψ | 2 ) P ln P | Ψ | 2 P | Ψ | 2 ( P | Ψ | 2 ) ,
for which any of the possible equalities only hold when P = | Ψ | 2 . (This is immediate from the trivial inequality: x > 0 , ( 1 1 / x ) ln x x 1 .)
Now, since both P ( x , t ) and | Ψ ( x , t ) | 2 are probability distributions, i.e., since we have Ω P d x = Ω | Ψ | 2 d x = 1 , it follows from (30) that whenever H V ( t ) and L f ( t ) are well-defined, they satisfy the following inequalities:
0 H V ( t ) L f ( t ) .
Moreover, for the same reason, L f can be re-expressed as Ω d 3 x f ( P | Ψ | 2 ) ( P | Ψ | 2 ) , the integrand in which is non-negative due to Equation (30). Therefore, L f can only be zero if its integrand is zero, i.e., if P = | Ψ | 2 (if P, | Ψ | 2 , and f are sufficiently smooth, which is something we shall always assume unless otherwise stated). Similarly [21], one also has that H V can only be zero when P = | Ψ | 2 everywhere in Ω .
Let us now prove the relaxation to quantum equilibrium. Substituting P = f | Ψ | 2 in the Fokker–Planck Equation (18), and using Equations (4) and (10), it is easily verified that
| Ψ | 2 f t = α 2 · ( | Ψ | 2 f ) | Ψ | 2 m ( f ) ( S ) .
Rewriting L f as
L f = Ω d 3 x f ( f 1 ) | Ψ | 2 ,
its behavior in time can be calculated using Equations (32), (4), and (10):
d L f d t = Ω d 3 x · | Ψ | 2 m ( f 2 f ) S + α 2 ( 2 f 1 ) · ( | Ψ | 2 f )
= α 2 Ω d 3 x · ( 2 f 1 ) | Ψ | 2 f 2 f 2 | Ψ | 2
= α Ω d 3 x f 2 | Ψ | 2 ,
which is of course strictly negative, for all t, as long as f and | Ψ | 2 are not identically zero. Hence, if | Ψ | 2 is not zero throughout Ω , L f will decrease monotonically for as long as f is not (identically) equal to 1 on Ω , and therefore necessarily converges to 0, a value it can only attain when f 1 or, equivalently, when P | Ψ | 2 . We have thus established a strong H-theorem showing that, in the case of Nelson dynamics, any probability distribution P necessarily converges to | Ψ | 2 , if the latter does not become zero identically. Note that this excludes the case of a free particle for which lim t + | Ψ ( x , t ) | 2 = 0 , for all x, which means that d L f d t tends to zero even when f does not converge to 1.
A result, similar to the above, is also easily established for H V since L f dominates the latter, or alternatively from the formula
d H V d t = α 2 Ω d 3 x f 2 | Ψ | 2 f .
The above results show that (excluding the case of the free particle) Nelson dynamics, naturally, exhibits relaxation towards quantum equilibrium and that it does so for general initial probability distributions (at least, as long as the initial distribution is smooth enough). In this stochastic setting, there is therefore no need for any assumptions on the microstructure of the initial distributions, nor is there any need for the coarse-grained hypothesis when deriving an H-theorem.
Note that these results also show that we have, in fact, convergence of the distribution P to the quantum equilibrium distribution | Ψ | 2 in the L 1 norm. This is a consequence of the so-called Csiszár-Kullback-Pinsker inequality [43]:
L 1 2 H V
where
L 1 = Ω d 3 x | P | Ψ | 2 | .
This generalizes the results by Petroni and Guerra [22,23] obtained in their study of the relaxation towards quantum equilibrium in the framework of the Nelson dynamics of a single particle in a harmonic potential. The L 1 norm is also used by Efthymiopoulos et al. [25] in the context of the dBB theory.
We shall illustrate these results by means of numerical simulations for the case of a ground state for the 1D-harmonic oscillator in Section 5.4, for the case of the 2D-harmonic oscillator in Section 6.2, and in the case of a coherent state in Section 7.1.
A last important remark concerns the influence of possible zeros in the equilibrium distribution Ψ ( x , t ) , which would give rise to singularities in the osmotic velocity terms in the Ito equation, Equation (16), or in the Fokker–Planck equation, Equation (18) (or equivalently in Equation (32)), and might make the functions H V and L f ill-defined. In Section 6.2, we discuss the case of the first excited state of the 1D-harmonic oscillator, for which Ψ ( x , t ) has a node at x = 0 . One could in fact imagine studying higher excited states for which one has a finite number of nodes. In that case, the osmotic velocity of Equation (15) will have simple poles at a finite number of positions in x. At the level of the Ito equation, one would not expect a finite set of poles to cause any particular problems, not only because the probability of hitting a pole exactly in the stochastic evolution is zero but also because the osmotic term tends to move the particle away from the pole very quickly. Similarly, a finite number of simple poles in the convection-diffusion equation, Equation (32), for f only influence the velocity field in the convection term in a finite number of distinct places, and it is to be expected that this would have the effect of actually enhancing the mixing of information in the system.
Moreover, it is also clear that simple nodes in Ψ ( x , t ) only give rise to (a finite number of) logarithmic singularities in the integrand of H V and that the integral in Equation (27) therefore still converges. The H-theorem for H V derived above is thus still valid, and an arbitrary distribution P (sufficiently smooth) will still converge to quantum equilibrium, even in the presence of nodes for Ψ ( x , t ) . The same cannot be said, however, of the function L f , as simple zeros in Ψ ( x , t ) give rise to double poles in the integrand and a possible divergence of the integral of Equation (28). Hence, at the beginning of the evolution, for an arbitrary P, the function L f might take an infinitely large value (the integrand only diverges when | Ψ | 2 P , i.e., when it is positive) but as soon as convergence sets in (which is guaranteed by the H-theorem for H V ), the divergent parts in its integrand will be smoothed out and the function L f will take finite values that converge to zero as time goes on. Of course, when calculating these quantities for the results of numerical simulations, there is always some amount of coarse-graining going on and genuine infinities never occur.

5. Relaxation to Quantum Equilibrium and Nelson Dynamics: The Static Case

In this section, in order to simplify the discussion, we will only consider the case of stationary states Ψ s t ( x ) for the one dimensional Schrödinger equation, i.e., energy levels for which S = E t and which therefore have zero Bohm velocity (11): S S x = 0 .

5.1. Fokker–Planck Operator and a Formal Connection to a Schrödinger Equation

There exists a wide literature [45,46] concerning a particular method for studying the convergence of solutions of the Fokker–Planck equation to a stationary one, which is only sporadically mentioned in the literature devoted to Nelson dynamics [47]. This approach makes it possible to quantify very precisely the speed of convergence to equilibrium, in terms of (negative) eigenvalues of the Fokker–Planck operator. In order to show this, let us rewrite the Fokker–Planck Equation (18) in terms of the Fokker–Planck operator L ^ :
P t = L ^ P = γ x γ ( x ) x + α 2 2 x 2 P
where (15):
γ ( x ) = α | Ψ s t | x | Ψ s t | .
Note that, due to the presence of the first derivative x , the L ^ operator is not Hermitian.
Now, in order to establish the H-theorem, we must prove that in the long-time limit this equation tends to a stationary solution P s t = | Ψ s t | 2 . The key idea here is to transform the Fokker–Planck equation to a simple diffusion equation through the transformation
P ( x , t ) = P s t ( x ) g ( x , t ) ,
under which the r.h.s. of Equation (40) reduces to
L ^ P = P s t ( x ) H ^ s t g ( x , t )
where H ^ s t is now a Hermitian operator:
H ^ s t = α 2 2 x 2 1 2 γ x + γ 2 α .
The function g ( x , t ) thus obeys a “Schrödinger-like” equation (though with imaginary time) with an effective potential ( H ^ s t ) that depends on γ ( x ) :
g ( x , t ) t = H ^ s t g ( x , t ) .
Note that the effective potential is exactly the Bohm-quantum potential defined by
Q Ψ = 2 2 m 1 | Ψ s t | 2 | Ψ s t | x 2 ,
which can be expressed in terms of the osmotic velocity (41) as
Q Ψ m α = 1 2 γ x + γ 2 α .

5.2. Superposition Ansatz

We can now represent the solution of Equation (45) as a superposition of discrete eigenvectors (all orthogonal, as the operator H ^ s t is Hermitian) and impose the superposition ansatz [48]:
g ( x , t ) = k = 0 a k ( t ) g k ( x ) .
Equation (45) is separable and gives rise to the eigenvalue problem:
1 a k ( t ) d a k ( t ) d t = 1 g k ( x ) H ^ s t g k ( x ) = λ k .
As a result, we have
g ( x , t ) = k = 0 a k e λ k t g k ( x )
for a set of constants a k and where all the λ k are real (as H is Hermitian), for eigenfunctions g k ( x ) that satisfy the orthonormality conditions:
d x g k ( x ) g l ( x ) = δ k , l .
Thus, we have the expression
P ( x , t ) = k = 0 a k e λ k t P s t ( x ) g k ( x ) .
By construction, the function P s t ( x ) is an eigenstate of the effective Hamiltonian with energy 0. We shall associate the label λ 0 with this energy level.
In order to have a well defined probability distribution and to avoid any divergence in time, it is clear that all eigenvalues λ k have to be negative, which requires Ψ s t to be the ground state of the effective Hamiltonian H s t : just as in the case of the usual Schrödinger equation, the eigenvalues λ k in Equation (49) are all negative only if Ψ s t ( x ) has no zeros (see also [49], Appendix 2, for an elementary proof that all λ k are indeed positive if Ψ s t ( x ) h as no zeros).
If Ψ s t ( x ) does have zeros, the osmotic velocity will have singularities. In [49] (Appendix 2), we consider what happens in the case when Ψ s t ( x ) is an excited state of the harmonic oscillator and we derive a formal solution in terms of the eigenvalues λ k , which are now not all negative, thus revealing the appearance of instabilities for cases where the above formalism would still be valid.

5.3. One-Dimensional Oscillator and the Evolution of Gaussian Distributions for the Ground State

In [49] (Appendix 2), we discuss the application of the method of the effective Hamiltonian outlined in Section 5.1 to this particular problem, and we derive a Green function for the associated Fokker–Planck equation when Ψ s t is the ground state of the one-dimensional oscillator. This Green function K P ( x , x , t ) is defined through
P ( x , t ) = d x P ( x , 0 ) K P ( x , x , t )
where the kernel K P is given by
K P ( x , x , t ) = a π s i n h ( ω t ) 1 2 e ω n + 1 2 t × e a s i n h ( ω t ) ( x 2 + x 2 ) c o s h ( ω t ) + ( x 2 x 2 ) s i n h ( ω t ) 2 x x .
An important property of the Green function for this case is that, if | Ψ ( x ) | 2 and P ( x , 0 ) are Gaussian, then P ( x , t ) will still be Gaussian (53). Let us define the ground state as
| Ψ s t | 2 | Ψ ( x ) | 2 = 2 a π e 2 a x 2 ,
for which we can then write
P ( x , t ) = 2 b ( t ) π e 2 b ( t ) x x ( t ) 2 .
Injecting Equation (56) in the Fokker–Planck equation, Equation (40), gives a differential equation for x ( t ) ,
d x ( t ) d t = 2 a α x ( t ) ,
which is readily solved as
x ( t ) = x 0 e 2 a α t ,
as well as an equation for b ( t ) ,
1 2 b ( t ) d b ( t ) d t + 2 α b ( t ) a = 0
with solution
b ( t ) = a 1 1 a b 0 e 4 a α t .
From Equations (56) and (60), we can then calculate the width of the non-equilibrium Gaussian as
σ x 2 ( t ) 1 4 b ( t ) = 1 4 a 1 e 4 a α t + a b 0 e 4 a α t = σ e q 2 1 e 4 a α t + σ x 2 ( 0 ) e 4 a α t
where σ e q 2 represents the width 1 / ( 4 a ) of the equilibrium distribution of Equation (55).
Clearly, x = t x e q = 0 with a characteristic relaxation time inversely proportional to the diffusion coefficient α . Moreover,
d σ x ( t ) d t 4 a α σ e q 2 σ x 2 ( 0 ) e 4 a α t ,
which has the same sign as that of the difference ( σ e q σ x ( 0 ) ) . Hence, σ x ( t ) converges monotonically to the equilibrium value σ e q , with a characteristic time inversely proportional to the diffusion coefficient α , as can be seen in Figure 2.

5.4. Ergodicity in the Relaxation to Quantum Equilibrium for the Ground State of the Harmonic Oscillator

We have just shown how Gaussian initial distributions converge towards quantum equilibrium, but one could also ask the same question for non-Gaussian initial distributions. Convergence is guaranteed by the H-theorem of Section 4, but contrary to the Gaussian case, we have no clear measure for the rate of convergence, except for the entropy-like functions H V (27) and L f (28), or the L 1 norm (39), defined in Section 4. The evolution, in time, of these three quantities is shown in Figure 3, for the stochastic trajectories obtained from 20,000 uniformly distributed initial conditions. The relaxation towards quantum equilibrium is clearly visible in all three quantities. As expected, the convergence of H V is extremely fast. Note that, although initially very large, L f quickly matches L 1 , up to numerical fluctuations.
One important question concerning this relaxation process is of course that of possible ergodicity. Since we want to study the ergodic properties of Nelson dynamics in a numerical way, we choose the definition of ergodicity that is, in our approach, the easiest to test. Defining the time average h ^ of a function h on Ω by the limit (if it exists),
h ^ = lim t + 1 t 0 t h x t d t
(where x t represents the position of a particle at time t , as obtained form the Ito stochastic differential equation, Equation (16), for an initial condition x ) we say [50] that the corresponding stochastic process is ergodic if the time average of any bounded function h on Ω is always independent of x . Since for bounded h the time average is also invariant under shifts in time, we can say that we have ergodicity if all time averages of such functions are in fact constants. The main reason for choosing this particular definition is that it is well-suited to empirical testing, since it is of course sufficient to establish constancy of the time averages for all indicator functions χ A of arbitrary (measurable) sets A Ω , for the analogous property to ensue automatically for all bounded functions on Ω . Another reason for choosing this particular definition is that it is also applicable to non-stationary stochastic processes, as in the case of the coherent state of Section 7.
More precisely, we need to verify that
χ ^ A = lim t + 1 t 0 t χ A x t d t
is independent of both t and x , for any measurable A Ω . Remember that one has of course that χ A ( x t ) = χ ϕ t 1 A ( x ) , where ϕ t 1 A = x Ω | x t A .
In the present case, i.e., that of the Nelson dynamics defined by the stationary (ground) state of the 1D harmonic oscillator, it is clear that the distribution | Ψ s t | 2 obtained from the ground state eigenfunction Ψ s t is a stationary solution to the associated Fokker–Planck equation, Equation (18). This distribution provides a natural invariant measure μ on Ω : d μ = | Ψ s t | 2 d x , for which Ω d μ = 1 and
μ ( A ) = A | Ψ s t | 2 d x = μ ( ϕ t 1 A ) , t > 0 , A Ω .
If a stationary stochastic process is ergodic, i.e., if for such a process all χ ^ A are indeed constants, the values of these constants are simply the measures of the subsets A [51]. Therefore, when one needs to decide whether or not a stationary stochastic process is ergodic, it suffices to verify that χ ^ A = μ ( A ) , for any A Ω .
The usual way to check this condition is to consider sampling time averages for a sufficiently refined “binning” of Ω . Starting from a particular initial particle position x , we calculate the trajectory x t that follows from the Ito stochastic equation, Equation (16), for a sufficiently long time t. As was explained for the coarse-graining in Section 3, the configuration space Ω is subdivided into a large number of non-overlapping cells or “bins” A k ( k = 1 , , N b ), each with the same volume Δ x . The trajectory x t ( t [ 0 , t ] ) is then sampled at regular intervals Δ t , yielding N + 1 sample positions x n Δ t ( n = 0 , , N ) , for N = t / Δ t . We then define the sampling function φ N , k
φ N , k = 1 N n = 0 N χ A k ( x n Δ t ) ,
which is a discretization of 1 t 0 t χ A x t d t in Equation (64) and which gives the frequency with which the (sample of the) orbit visited the bin A k . Hence, if in the limit N + , for diminishing bin sizes Δ x and sampling steps Δ t , the normalized distribution obtained from φ N , k / Δ x tends to a constant distribution (and, in particular, does not depend on the initial positions x ), then the stochastic process is ergodic according to the above definition.
Moreover, since in that case χ ^ A k = μ ( A k ) , this normalized distribution must in fact coincide with that for the invariant measure for the stationary process. For example, in the case at hand, if the normalized distribution we obtain is indeed independent of the initial positions, then since μ ( A k ) = | Ψ s t ( x ) | 2 | x = ξ Δ x for some point ξ A k , we must have that for sufficiently large N
φ N , k Δ x μ ( A k ) Δ x = | Ψ s t ( x ) | 2 | x = ξ ,
i.e., the empirical distribution obtained from this sampling time average must coincide with the stationary quantum probability | Ψ s t | 2 . This is exactly what we obtain from our numerical simulations, as can be seen from the histograms depicted in Figure 4. After a certain amount of time, the histograms we obtain indeed converge to the equilibrium distribution, and this for arbitrary initial positions. The convergence clearly improves if we increase the integration time, or if we diminish the spatial size of the bins (while diminishing the sampling time step in order to keep the occupancy rate of each bin high enough). Although purely numerical, we believe this offers conclusive proof for the ergodicity of the Nelson dynamics associated with the ground state of the harmonic oscillator in one dimension.
The same can be said, in fact, for the 2-dimensional oscillator which will be the main topic of Section 6 (and even for the 2D corral as can be seen in [49]). Some results of a simulation of a single trajectory under the Nelson dynamics for the ground state of this system are shown in Figure 5, in which the red dot in the plot on the left-hand side indicates the (final) position of the particle at time t. The probability distribution obtained by sampling the trajectory clearly decreases with the distance to the origin along concentric circles.

6. Nelson Dynamics: A Phenomenological Dynamical Model for Walkers?

6.1. 2D Harmonic Oscillator

Experimentally, it has proven possible to study the dynamics of bouncing droplets under the influence of an effective harmonic potential in two dimensions, thanks to a well-chosen electro-magnetic configuration and magnetic droplets [52]. It is therefore interesting to compare predictions that we, on our side, can make in the framework of Nelson dynamics, with actual experimental observations of droplets dynamics (see [53] for a pioneering work very similar to ours in the case of the double slit experiment). We think that an important comparison to make concerns the convergence to equilibrium.
For example, if the initial distribution of positions projected along a reference axis, say X, fits a mixture of the ground state and the nth Fock state ( n = 1 , 2 ) (Appendix 2, [49]) for the 2D harmonic oscillator (conveniently weighted in order to respect the ineluctable constraint of positivity), our Nelson-like model predicts that the typical time of convergence to equilibrium will scale like the inverse of the eigenvalue of the nth Fock state, i.e., as 1 / n , which constitutes a very precise quantitative prediction. This follows from Equation (52), when P s t ( x ) is the Gaussian ground state of the 1D harmonic oscillator and where the eigenfunctions g k are the energy eigenstates (Fock states) of the harmonic oscillator (this, of course, because of the separability of the Schrödinger equation and of our Nelson dynamics along X and Y in the case of an isotropic 2D oscillator).
A possible way to measure this characteristic time would be to record the projections along X of trajectories that correspond to an equally spaced grid of initial positions, weighted so as to fit a mixture of the ground state with the nth Fock state ( n = 1 , 2 ), and to compare the histogram constructed in this way at different times with theoretical predictions derived from Equation (52).
Another precise quantitative (theoretical) prediction, which is even simpler to verify, is that if we prepare a droplet many times at exactly the same initial position, the position obtained after averaging over all trajectories will (1) decrease exponentially in time and (2) be characterized by a decay time which scales like 1 / a α , by virtue of the discussion in Section 5.3 and particularly Equation (58). It has been suggested that droplet trajectories might be characterized by a quantum-like Zitterbewegung, which can be seen in relativistic quantum dynamics as an intrinsic Brownian motion at the Compton scale [54,55], and various proposals have been formulated in order to express the amplitude and frequency of this Zitterbewegung [5,56] in terms of the parameters characterizing droplet dynamics (these are, e.g., the viscosity of the fluid, the mass of the droplets, the ratio between the amplitude of the vibrations imposed on the bath and the Faraday threshold, the oil temperature [5,57], and so on). Exploring these analogies in depth lies beyond the scope of this paper, but the aforementioned attempts ([56] in particular) pave the way for introducing a Brownian component in the description of droplet trajectories.

6.2. Presence of Zeros in the Interference Pattern

One of our first motivations, when we decided to incorporate a Brownian component in the dBB theory in order to simulate the dynamics of droplets, was the pioneering paper [4] reporting on observations of a walker trapped in a spherical 2D cavity (corral), for which the histogram of positions occupied over time by a single droplet trajectory faithfully reproduces the Bessel function J 0 (this is also related to the Green function of the Helmholtz equation, with a typical length equal to the Faraday wave length of the vibrating bath over which droplets propagate [16]). These observations reveal, in a telling way, the presence of a pilot-wave that guides the dynamics of the particles, and raise the question of ergodicity.
If we try the approach we used for the 2D harmonic oscillator in the case of the corral (effectively replacing the Gaussian ground state of the 2D harmonic oscillator by the zero order Bessel function), we are immediately confronted with problems caused by the presence of zeros in the Bessel function. In particular, the eigenvalues λ k of the Fokker–Planck operator in Equation (49) are not always negative when zeros are present, which of course would menace the stability of the relaxation process. However, as we already indicated in Section 4, although the effect of zeros of the pilot wave in our Nelson dynamics is by no means trivial, there are several observations that indicate that this problem is not crucial.
First of all, as mentioned in Section 4, the Wiener process makes it in principle possible to “jump” over the zeros of the equilibrium distribution. This has actually been confirmed in numerical simulations for the case of the 1D harmonic oscillator, where we imposed that the equilibrium distribution P s t is the square modulus of the first excited (Fock) state ([49], Appendix 2A), with the following amplitude:
P s t = | Ψ s t | 2 = | Ψ 1 ( x , t ) | 2 = 2 a π 1 2 a x 2 e 2 a x 2 .
Indeed, as can be clearly seen from Figure 6, the particle will, from time to time, jump over the zero in the middle (with the same probability from left to right as in the opposite direction), in such a way that finally the trajectory covers the full real axis, while the histogram of positions faithfully reproduces the quantum prediction P s t = | Ψ s t | 2 = | Ψ 1 ( x , t ) | 2 . This indicates that, even in the presence of a zero in the equilibrium distribution, the relaxation process is still ergodic.
The relaxation of a uniform initial distribution to this quantum equilibrium is shown in Figure 7, for the quantities H V , L f , and L 1 .
A second indication that the problem posed by the presence of zeros is not so serious stems in fact from the experimental observations. Indeed, if we study the observations reported in [4] for the case of a corral, it is clear that the minima of the histogram expressing the distribution of positions of the droplet are in fact not zeros. This, undoubtedly, is due to the presence of a non-negligible residual background. Without this background, the droplet would never pass between regions separated by zeros: due to the rotational symmetry of the corral, the zeros form circles centered at the origin and the position histogram obtained from a trajectory would remain confined to a torus comprising the initial position. This, however, is clearly not the case, which thus suggests the following strategy: to simulate Nelson dynamics with a static distribution P s t = | Ψ s t | 2 given by the Bessel function J 0 but supplemented with a constant positive background ϵ :
d x ( t ) = α 2 J 0 ( r ) 2 J 0 ( r ) 2 + ϵ d t + α d W ( t ) .
In this case, the singularities of the Fokker–Planck equation automatically disappear and, despite the fact that we have no analytic expression for the solutions as in the case of the ground state of the harmonic oscillator, we are able to numerically simulate Nelson dynamics without difficulty. The results of these simulations are shown in Figure 8. The osmotic velocity in the Nelson dynamics clearly tends to bring back the particle to regions where | Ψ | 2 has extrema and the resemblance with the plot on the left is striking. The fact that this result again does not depend on the choice of initial condition strongly suggests that the relaxation process to quantum equilibrium is also ergodic in this case.

7. Relaxation to Quantum Equilibrium with dBB and Nelson Dynamics: The Non-Static Case

7.1. Nelson Dynamics and Asymptotic Coherent States

Up to now, we have developed analytic and numerical tools aimed at studying the onset of equilibrium when the asymptotic equilibrium distribution is static. As the H-theorem of Section 4 is also valid for non-stationary processes, one of course expects relaxation to take place even if the asymptotic state is not static, for instance, when it is a Gaussian distribution, the center of which periodically oscillates at the classical frequency ω of the oscillator without deformation (typical for coherent states). In fact, our numerical simulations show not only that equilibrium is reached even in this case, but also that this relaxation is ergodic.
More precisely, we considered a wave function in the coherent state
Ψ ( x , t ) = 2 a π 1 4 e a x x ¯ t 2 + i p ¯ t x + i φ ( t )
where φ is a global phase containing the energy, and x ¯ t ( p ¯ t ) is the mean position (momentum) of a classical oscillator at time t:
x ¯ t = x ¯ 0 c o s ω t and p ¯ t = m x ¯ 0 s i n ω t ,
with ω = 2 a α ( α = / m ). For this ansatz we solved the Ito equation, Equation (16), numerically for a collection of initial conditions.
As can be seen in Figure 9, the trajectories are affected by the stochastic evolution but keep oscillating at the same period because of the deterministic part of the Ito process. Notice, however, that the trajectories seem to be approaching classical trajectories that only differ from each other by a simple shift. This can be explained as follows: at equilibrium (cf. Figure 10), the Brownian motion is balanced by the osmotic velocity and the dBB velocity is recovered “on average.” Now, the center of the Gaussian distribution moves at a classical velocity by virtue of Ehrenfest’s theorem; moreover, in the present case, the dBB velocities can only depend on time and not on space as the envelope of a coherent state moves without deformation. Hence, the dBB trajectories obtained at equilibrium are, in fact, classical trajectories that only differ by a mere shift in space (the magnitude of which, however, may change over time).
Secondly, as can be clearly seen on Figure 10, even for a uniform initial probability distribution, the convergence to the quantum equilibrium is remarkably fast and the converged distribution faithfully follows the oscillating motion of the non-stationary equilibrium distribution. The remarkable speed of the convergence to quantum equilibrium is corroborated by the decay of the functions H V and L f and of the L 1 norm shown in Figure 11.
Moreover, Figure 12 depicts the sampling time average (as defined in Section 5.4) of a single trajectory for this non-stationary stochastic process. The convergence of the sampling distribution to a static distribution Φ ( x ) , described by the integral of | Ψ ( x , t ) | 2 as given by Equation (70), over a period of the oscillation
Φ ( x ) = ω 2 π 0 2 π / ω | Ψ ( x , t ) | 2 d t
is striking. As the asymptotic distribution Φ ( x ) does not depend on the choice of initial condition, we conclude that the relaxation to equilibrium for the non-stationary stochastic process associated with Nelson dynamics for the coherent state (70) is ergodic (in the sense explained in Section 5.4).

7.2. Onset of Equilibrium with a Dynamical Attractor in dBB Dynamics and Nelson Dynamics

If one wants to investigate the onset of equilibrium in dBB dynamics, one obviously has to consider non-static asymptotic distributions since in static cases the dBB dynamics freezes the trajectories (as the phase of the wave function is then position-independent). Even in the case of a coherent state (see Section 7.1), the distribution of dBB positions would never reach equilibrium because it moves as a whole (as the shape of a coherent state remains the same throughout time). In a sense, coherent states behave as solitary waves. Moreover, the absence of zeros in the wave function might explain why mixing does not occur. In Figure 13 we show the result of simulations of dBB trajectories in the case of a 2D harmonic oscillator for a quantum state consisting of a superposition of equally weighted products of states along X and Y, chosen among M energy (Fock) states ([49], Appendix 2A), with randomly chosen initial phases θ n x , n y :
Ψ ( x , y , t ) = 1 M n x = 0 M 1 n y = 0 M 1 e i θ n x , n y i ω n x + n y + 1 t ψ n x x ψ n y y .
We then compared the relaxation process for dBB with the quantum thermostat given by Nelson dynamics for M = 4 2 = 16 energy states. The results are shown in Figure 14 in which H V (for the dBB and for the Nelson dynamics) and L 1 (for both the dBB and Nelson dynamics) are plotted at the (same) coarse-grained level. We started from a uniform distribution of positions; we took α = 0.1 . In both cases, the position distributions P and P converge to | Ψ | 2 . Moreover, we recover an exponential decay for H ¯ V , as already observed in [26], even in the absence of stochastic (Brownian) noise à la Nelson. However, we observe that the convergence to equilibrium occurs faster in the presence of the quantum thermostat.

8. Dynamical Model for Droplets and Double Quantization of the 2-D Harmonic Oscillator

In this section we shall focus on the description of droplets dynamics as described in [9,52], for a magnetized droplet moving in an isotropic 2-D harmonic potential. We shall show that dBB dynamics allows us to reproduce some of the main features of the experimental observations. In [9,52], it is reported that stable structures appear in the droplets dynamics whenever a double quantisation condition is satisfied. The Hamiltonian of the isotropic 2-D harmonic oscillator being invariant under rotations, we may indeed impose a double quantisation constraint, requiring that the energy states of the 2D quantum harmonic oscillator are also eigenstates of the angular momentum. In polar coordinates, these states (which are parameterized by two quantum numbers, the energy number n and the magnetic number m) are expressed as follows [58]:
ψ n , m r , θ , t = a π k ! k + | m | ! e a r 2 2 a r m L k m a r 2 e i ω ( n + 1 ) t + i m θ
where L k m are the generalized Laguerre polynomials and k = n m 2 . Note that these solutions are linear combinations of the product of Fock states in x and y.
A first experimental result reported in [9] is the following: trajectories are chaotic and nearly unpredictable unless the spring constant of the harmonic potential takes quantized values that are strongly reminiscent of energy quantization (under the condition that, during the experiment, the size of the orbits is fixed once and for all). For quantized energies—in our case given by E n = ( n + 1 ) ω , for some “effective” value of to be determined from actual experiments—stable orbits, to which one can attribute yet another quantum number, appear, this time for the angular momentum, which is strongly reminiscent of the magnetic number (the eigenvalue of the orbital momentum, perpendicular to the surface of the vessel, is given by the product of and m). In [9] it is shown, for instance, that for the first excitation (n = 1, m = ± 1 ) droplet orbits are circular or oval, turning clockwise or anti-clockwise depending on the sign of m. At the second energy level (n = 2, m = 2 , 0 , + 2 ), ovals appear again for m = ± 2 and lemniscates for an average value of the angular momentum < m > = 0 . At the fourth energy level (n = 4, m = 4 , 2 , 0 , 2 , 4 ), trefoils appear (for m = ± 2 ).
We simulated dBB trajectories, always considering a superposition of one of the aforementioned doubly quantized eigenstates ψ n , m with the ground state:
Ψ r , θ , t = ξ 0 e i φ 0 ψ 0 , 0 r , θ , t + j = 0 n ξ j + 1 e i φ j + 1 ψ n , n + 2 j r , θ , t
where φ j and ξ j are real numbers with 0 < ξ 0 ξ j 0 . Computing the guidance relation of Equation (11) for a single eigenstate (74), one ends up with a value for S for which the trajectories are circles of radius R around the origin, with tangential velocities proportional to m / R . In particular, the dynamics is frozen when m = 0 .
Mixing the wave function with the ground state, however, generates a periodic (in time) component in the dBB velocity field, which turns circular orbits into ovals when ξ 0 is small enough, and eventually generates more complex structures, such as “rosaces” instead. We also tuned the energy difference between the ground state and the excited states such that two timescales characterize the dynamics. These are the “centrifugal” period, necessary for drawing a full circle around the origin, which varies as m / R 2 , and the “Bohr” period which varies like T / ( n + 1 ) , where T is the classical period of the oscillator. Tuning these parameters, we were able to simulate dBB trajectories very similar to those reported in [9]. For instance, we found circles and ovals (see Figure 15a,b) for ( n , m ) = ( 1 , 1 ) and ( n , m ) = ( 2 , 2 ) . Note that the lemniscate cannot be obtained with a superposition of the ground state and the ( n , m ) = ( 2 , 0 ) state for which dBB velocities are necessarily purely radial, contrary to the suggestion made in [9]. Instead, it should be generated with a superposition of the ground state with ( n , m ) = ( 2 , + 2 ) , ( 2 , 2 ) , and ( 2 , 0 ) in which the weights of the m = + 2 and 2 components are slightly different (see Figure 15c). Figure 16 shows further detail of the evolution along this trajectory. Tuning the energy, we were also able to generate a trefoil and a “rosace” (see Figure 17).
It is worth noting, however, that chaos is omnipresent in the dBB dynamics for this system, in the sense that the trajectories exhibit an extreme sensitivity to the initial conditions, which explains why these dBB orbits mimicking stable droplets orbits are in general unstable. For instance, Figure 18 shows intermittent transitions between an oval trajectory and a lemniscate (as has also been reported in [9]), for a superposition of the ground state with the ( n , m ) = ( 2 , + 2 ) , ( 2 , 2 ) , and ( 2 , 0 ) states. Preliminary results furthermore show that the trajectories are also unstable under Nelson dynamics, i.e., in the presence of “noise,” whenever this noise (parameterized by α in (16)) exceeds a critical value. Note that many experiments involving droplets are characterized by a lack of stability and predictability. For instance, the appearance of interferences similar to those obtained in a double slit experiment (see [49,59] for a description à la Nelson of the double slit experiment) has been attributed to “air currents” in [60]. Therefore, although our approach might not explain every detail of the double quantization reported in [9], it does reproduce many of its essential features, and we believe it would be very interesting to deepen this analogy. For instance, having access to the empirical values of the weights of the ground state, or of the effective values of and of the mass in the case of droplets [56], would allow us to test our model in real detail.
Another experiment, reported in [61], during which both the position of the droplet and the excitation of the bath are monitored, and where a superposition between two distinct modes of the bath is reported, could also provide more insight and might offer some means to test the validity of our model: using exactly the same observation device, but this time in the case where the droplet undergoes a 2-D isotropic potential, would allow one to check whether the modes of the bath are similar to the ( n , m ) quantum modes which we associate with the quantized droplets trajectories.

9. Conclusions and Open Questions

In this paper we studied stochastic, Nelson-like dynamics and dBB dynamics, with the aim of simulating the dynamics of droplets. The stochastic approach has the merit that it explicitly takes into account the influence of noise on the dynamics [59,62]. In contrast to experiments where noise is considered to be a parameter that should be minimized, here, noise is considered to be a relevant parameter for the dynamics (see also [53]). For instance, as we have shown, it plays an essential role in the relaxation towards equilibrium and in the ergodicity of the dynamics. In the dBB approach, on the other hand, the main ingredient is the chaotic nature of the dynamics [25]. Both models thus shed a different light on the dynamics and could possibly fit diverse sets of regimes in droplet dynamics. Note that in the limit where the amplitude of the Brownian motion in our Nelson dynamics tends to zero, the dynamics approaches dBB dynamics very closely. In sufficiently complex situations (e.g., when the mixing process due to the presence of zeros in the wave function becomes effective [26,40]), we expect the relaxation to equilibrium to be accompanied by chaotic rather than stochastic dynamics, as one has in Nelson dynamics (although Nelson dynamics with small but non-zero Brownian motion is hard to distinguish from dBB dynamics, it has the advantage that relaxation is guaranteed to occur, even in the absence of coarse graining and/or mixing).
Ultimately, experiments ought to indicate whether it is relevant, with respect to droplet phenomenology, to formalize the dynamical influence of noise à la Nelson (and/or dBB) as we did in the present paper. We have formulated several proposals in this sense in Section 6.1 and Section 8. As emphasized throughout the paper, however, our models should be seen as a first step in the direction of a dynamical model, which remains to be formulated, combining Nelson’s stochastic dynamics (and/or dBB dynamics) and memory effects. We think that the results of Section 8 show that this is a promising program for future research.
Finally, it is worth recalling some of the problems that arose when first de Broglie and then Bohm and Nelson developed their theories aimed at deriving quantum dynamics (statistics) as an emergent property, i.e., resulting from an underlying “hidden” dynamics.
The most severe problem is undoubtedly non-locality, which was recognized by Bohm [11,12] to be an irreducible feature of dBB dynamics (see also [19,63] for similar conclusions concerning Nelson-type dynamics). Today, under the influence of the work of John Bell [64] and his followers, it is widely recognized that quantum theory is irreducibly non-local, which makes it particularly difficult to mimic using classical models. Note that entanglement and non-locality (as well as decoherence, which is the corollary of entanglement [65]) only appear if we consider more than one particle at a time, which explains why we did not address these fundamental questions in the core of the paper, where a single droplet is described. It would be interesting to enlarge our model such that the presence of the environment can be taken into account. This would require incorporating the description of open quantum systems, for which a generalization of Bohmian dynamics has been developed in the past [66,67,68,69], but obviously this is beyond the scope of the present paper.
Another problem concerns the fact that the pilot wave is a complex function. This poses still unresolved problems in the case of Nelson dynamics because Nelson’s diffusion process does not make it possible [70,71] to fix the phase of the wave function unequivocally (see [72] for an interesting proposal involving a multivalued wave function, also based on Zitterbewegung). In our approach, which is mainly of quantum inspiration, complex wave functions and imaginary phases appear spontaneously, but if we wish to scrutinize the link with the empirically observed modes at the surface of oil baths [9,52,56,61], it will be important to interpret the exact meaning of this complex phase. In the framework of his double solution program [73,74] de Broglie, and others, showed how to derive the Schrödinger equation from a Klein–Gordon equation in the non-relativistic limit. This is only possible provided the real wave bounces at an extremely high frequency (of the order of m c 2 / h ). A similar approach has been proposed in the context of droplets phenomenology in [75], where a complex Schrödinger equation is derived from the Klein–Gordon equation along these lines. Although such (interesting and promising) alternative studies of droplets solve the problem of the appearance of a complex phase in a classical context, it is worth noting that the phenomenological results outlined in Section 8, concerning the quantization of droplet orbits in the case of a harmonic potential [9,52], cannot be explained simply in terms of excited modes of the oil bath, because in these experiments only the droplet undergoes the harmonic potential, the oil bath being electromagnetically neutral. This difficulty actually concerns any classical model in which droplet dynamics is formulated in terms of classical modes of the bath only.
To conclude, in our view, the programs that aim at simulating droplet dynamics with quantum tools or at describing the emergence of quantum dynamics based on droplet dynamics, are still largely incomplete and raise challenging fundamental questions. This Pandora’s box is now open and it will not be closed any time soon, but this is not something to be feared as it offers new and stimulating perspectives for future research in the field.

Author Contributions

Conceptualization, M.H., R.W., S.C. and T.D.; Formal Analysis, M.H., R.W. and T.D.; Funding Acquisition, R.W., S.C. and T.D.; Investigation, M.H. and R.W.; Methodology, M.H., R.W. and T.D.; Project Administration, R.W. and T.D.; Software, M.H. and S.C.; Supervision, T.D.; Validation, M.H. and R.W.; Visualization, M.H. and R.W.; Writing—Original Draft Preparation, M.H., R.W., S.C. and T.D.; Writing—Review & Editing, R.W. and T.D.

Funding

This research was funded by the John Templeton foundation (grant 60230).

Acknowledgments

The authors gratefully acknowledge funding and support from the John Templeton foundation (grant 60230, Non-Linearity and Quantum Mechanics: Limits of the No-Signaling Condition, 2016–2019) and an FQXi Physics of What Happens grant (Quantum Rogue Waves as Emergent Quantum Events, 2015–2017). RW would also like to thank FQXi for support through a mini-grant (Grant number FQXi-MGA-1819), which enabled him to elucidate crucial aspects of the relaxation to quantum equilibrium. SC also thanks the Foundational Questions Institute (fqxi.org) for its support though a mini-grant (Grant number FQXi-MGA-1705). T.D. also thanks Aurélien Drézet (CNRS-Grenoble) for drawing our attention to the paper of Kyprianidis and for useful comments on the draft of our paper during T.D.’s personal visit to Grenoble (3 May 2018). Sincere thanks to G. Grössing for his kind invitation to contribute to this special issue.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Numerical Simulations

Firstly, we discuss the case of the dBB dynamics. It is assumed that we have an analytical solution of the Schrödinger equation Ψ ( t , x ) . We want to compute the evolution of a given initial non-equilibrium density P ( t i , x ) up to a final time t f and for intermediate time events (we denote all these events by t k , with t 0 = t i and t f = t K ). In particular, we are interested in the coarse-grained non-equilibrium density
P ¯ ( x , t k ) = 1 ϵ 3 CG cell x d 3 x P ( x , t k ) ,
which is defined in Equation (23).
Numerically, we replace that integral by a discrete sum over a finite set of points x l , which are uniformly distributed over the CG cells. In order to obtain the value of each P ( x l , t k ) we use the Liouville relation
P ( x l , t k ) | Ψ ( x l , t k ) | 2 = P ( x i l , t i ) | Ψ ( x i l , t i ) | 2
where x i l is the position of the particle, which, when evolved according to Equation (11) from t i up to t k , gives x l .
In order to obtain x i l for each x l , we consider the time-reversed evolution with wave-function Ψ * ( t , x ) and initial condition x l at time t k . The position x l , if time evolved from t k up to t i according to Equation (11), will give the position x i l . As there is usually no analytical solution of Equation (11), we use a Runge–Kutta (RK) algorithm [76] to obtain a numerical estimate of the position x i l . In order to know if we can trust the result of the RK algorithm, we perform two realizations of the algorithm with different choices of a so-called tolerance parameter (the smaller the value of that tolerance parameter, the more precise the computation), say γ and 10 1 γ . If the distance between the two positions is less than some chosen value δ , the result of the last iteration of the RK algorithm is trusted. Otherwise, we perform another iteration with 10 2 γ and we compare it to the previous realization of the RK algorithm. We repeat the procedure until the constraint on the distance between the two successive results of the RK algorithm is satisfied, or until we reach some minimal value of the tolerance parameter. In that case, the position x l is considered as a bad position and it is discarded from the numerical integration of Equation (23). This method was used in [26].
That is one possible method but we could also adopt a more brute-force method: Randomly generate a set of N initial positions according to P ( t i , x ) and let them evolve according to an Euler algorithm (that is, we divide the time-interval in small time-steps of length Δ t and we increment the position by v ( t ) Δ t at each time step). We record the positions of the N particles for each value of t k , we count the number of particles in each CG cell for each time t k (say n C G ), and we divide n C G by N in order to define P ¯ ( x , t k ) . The first method turns out to be more efficient in the case of the dBB dynamics but it is not applicable in the presence of stochastic terms.
In the case of Nelson dynamics, we used the Euler–Maruyama method for stochastic processes to approximate the solution of the Ito equation, Equation (12). In the same way as Euler’s method, the time T is divided into N small discrete time steps Δ t . For each time t i , we generated a random variable normally distributed Δ W i = Δ t N 0 , 1 . The integration scheme has the form
x i + 1 = x i + v ( x i , i Δ t ) Δ t + α Δ W i .
We invite the reader interested in the details to consult [77]. The remaining question is how to choose the time step Δ t so that one can trust the result of the numerical simulations. One way to do this is the following. We know that the Born distribution remains invariant under Nelson’s dynamics (equivariance). We therefore start with some value for Δ t and decrease it until the result of the numerical simulation confirms this theoretical prediction. We then perform the numerical simulation for the non-equilibrium distribution with the value of Δ t thus obtained.

References

  1. Couder, Y.; Fort, E. Single-Particle Diffraction and Interference at a Macroscopic Scale. Phys. Rev. Lett. 2006, 97, 15410. [Google Scholar] [CrossRef] [PubMed]
  2. Couder, Y.; Boudaoud, A.; Protière, S.; Fort, E. Walking droplets, a form of wave-particle duality at macroscopic scale? Europhys. News 2010, 41, 14–18. [Google Scholar] [CrossRef] [Green Version]
  3. Couder, Y.; Protière, S.; Fort, E.; Boudaoud, A. Dynamical phenomena: Walking and orbiting droplets. Nature 2005, 437, 208. [Google Scholar] [CrossRef] [PubMed]
  4. Harris, D.M.; Moukhtar, J.; Fort, E.; Couder, Y.; Bush, J.W. Wavelike statistics from pilot-wave dynamics in a circular corral. Phys. Rev. E 2013, 88, 011001. [Google Scholar] [CrossRef] [PubMed]
  5. Bush, J.W.M. Pilot-wave hydrodynamics. Annu. Rev. Fluid Mech. 2015, 47, 269–292. [Google Scholar] [CrossRef]
  6. Couder, Y.; Fort, E. Probabilities and trajectories in a classical wave-particle duality. J. Phys. Conf. Ser. 2012, 361, 012001. [Google Scholar] [CrossRef] [Green Version]
  7. Bush, J.W.M. The new wave of pilot-wave theory. Phys. Today 2015, 68, 47–53. [Google Scholar] [CrossRef] [Green Version]
  8. Eddi, A.; Sultan, E.; Moukhtar, J.; Fort, E.; Rossi, M.; Couder, Y. Information stored in Faraday waves: The origin of a path memory. J. Fluid Mech. 2011, 674, 433–463. [Google Scholar] [CrossRef]
  9. Perrard, S.; Labousse, M.; Miskin, M.; Fort, E.; Couder, Y. Self-organization into quantized eigenstates of a classical wave-driven particle. Nat. Commun. 2014, 5, 3219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Durt, T. Do dice remember? Int. J. Theor. Phys. 1999, 38, 457–473. [Google Scholar] [CrossRef]
  11. Bohm, D. A Suggested Interpretation of the Quantum Theory in Terms of “Hidden” Variables. I. Phys. Rev. 1952, 85, 166–179. [Google Scholar] [CrossRef]
  12. Bohm, D. A Suggested Interpretation of the Quantum Theory in Terms of “Hidden” Variables. II. Phys. Rev. 1952, 85, 180–193. [Google Scholar] [CrossRef]
  13. Nelson, E. Dynamical Theories of Brownian Motion; Princeton University Press: Princeton, NJ, USA, 1967; Volume 131, pp. 2381–2396. [Google Scholar]
  14. Labousse, M. Etude D’une Dynamique à Mémoire de Chemin: une Expérimentation Théorique. Ph.D. Thesis, Université Pierre et Marie Curie UPMC Paris VI, Paris, France, 2014. (In French). [Google Scholar]
  15. Fort, E.; Eddi, A.; Boudaoud, A.; Moukhtar, J.; Couder, Y. Path-memory induced quantization of classical orbits. Proc. Natl. Acad. Sci. USA 2010, 107, 17515–17520. [Google Scholar] [CrossRef] [Green Version]
  16. Dubertrand, R.; Hubert, M.; Schlagheck, P.; Vandewalle, N.; Bastin, T.; Martin, J. Scattering theory of walking droplets in the presence of obstacles. New J. Phys. 2016, 18, 113037. [Google Scholar] [CrossRef] [Green Version]
  17. Tadrist, L.; Shim, J.B.; Gilet, T.; Schlagheck, P. Faraday instability and subthreshold Faraday waves: surface waves emitted by walkers. arXiv, 2017; arXiv:1711.06791. [Google Scholar] [CrossRef]
  18. Bohm, D.; Vigier, J.P. Model of the causal interpretation of quantum theory in terms of a fluid with irregular fluctuations. Phys. Rev. 1954, 96, 208–216. [Google Scholar] [CrossRef]
  19. Bohm, D.; Hiley, B. Non-locality and locality in the stochastic interpretation of quantum mechanics. Phys. Rep. 1989, 172, 93–122. [Google Scholar] [CrossRef]
  20. Kyprianidis, P. The Principles of a Stochastic Formulation of Quantum Theory. Found. Phys. 1992, 22, 1449–1483. [Google Scholar] [CrossRef]
  21. Valentini, A. Signal locality, uncertainty and the subquantum H-theorem. I. Phys. Lett. A 1991, 156, 5–11. [Google Scholar] [CrossRef]
  22. Petroni, N.C. Asymptotic behaviour of densities for Nelson processes. In Quantum Communications and Measurement; Springer: New York, NY, USA, 1995; pp. 43–52. [Google Scholar]
  23. Petroni, N.C.; Guerra, F. Quantum Mechanical States as Attractors for Nelson Processes. Found. Phys. 1995, 25, 297–315. [Google Scholar] [CrossRef]
  24. Guerra, F. Introduction to Nelson Stochastic Mechanics as a Model for Quantum Mechanics. In The Foundations of Quantum Mechanics—Historical Analysis and Open Questions; Garola, C., Rossi, A., Eds.; Springer: Dordrecht, The Netherlands, 1995; pp. 339–355. [Google Scholar]
  25. Efthymiopoulos, C.; Contopoulos, G.; Tzemos, A.C. Chaos in de Broglie–Bohm quantum mechanics and the dynamics of quantum relaxation. Ann. Fond. Louis Broglie 2017, 42, 133–160. [Google Scholar]
  26. Valentini, A.; Westman, H. Dynamical origin of quantum probabilities. Proc. R. Soc. A 2005, 461, 253–272. [Google Scholar] [CrossRef] [Green Version]
  27. Colin, S.; Struyve, W. Quantum non-equilibrium and relaxation to quantum equilibrium for a class of de Broglie–Bohm-type theories. New J. Phys. 2010, 12, 043008. [Google Scholar] [CrossRef]
  28. Towler, M.; Russell, N.J.; Valentini, A. Time scales for dynamical relaxation to the Born rule. Proc. R. Soc. A 2011, 468, 990–1013. [Google Scholar] [CrossRef] [Green Version]
  29. Colin, S. Relaxation to quantum equilibrium for Dirac fermions in the de Broglie–Bohm pilot-wave theory. Proc. R. Soc. A 2012, 468, 1116–1135. [Google Scholar] [CrossRef]
  30. Contopoulos, G.; Delis, N.; Efthymiopoulos, C. Order in de Broglie–Bohm quantum mechanics. J. Phys. A Math. Theor. 2012, 45, 165301. [Google Scholar] [CrossRef] [Green Version]
  31. Abraham, E.; Colin, S.; Valentini, A. Long-time relaxation in the pilot-wave theory. J. Phys. A Math. Theor. 2014, 47, 395306. [Google Scholar] [CrossRef]
  32. De la Peña, L.; Cetto, A.M. The Quantum Dice: An Introduction to Stochastic Electrodynamics; Springer: Dordrecht, The Netherlands, 2013. [Google Scholar]
  33. De la Peña, L.; Cetto, A.M.; Valdés Hernández, A. The Emerging Quantum; Springer: Basel, Switzerland, 2015. [Google Scholar]
  34. Deotto, E.; Ghirardi, G.C. Bohmian mechanics revisited. Found. Phys. 1998, 28, 1–30. [Google Scholar] [CrossRef]
  35. De Broglie, L. Interpretation of quantum mechanics by the double solution theory. Ann. Fond. Louis Broglie 1987, 12, 399–421. [Google Scholar]
  36. Bacciagaluppi, G. Nelsonian mechanics revisited. Found. Phys. Lett. 1999, 12, 1–16. [Google Scholar] [CrossRef] [Green Version]
  37. de Broglie, L. La mécanique ondulatoire et la structure atomique de la matière et du rayonnement. J. Phys. Radium 1927, 8, 225–241. (In French) [Google Scholar] [CrossRef]
  38. Valentini, A. On the Pilot-Wave Theory of Classical, Quantum and Subquantum Physics. Ph.D Thesis, Scuola Internazionale Superiore di Studi Avanzati (SISSA), Trieste, Italy, 1992. [Google Scholar]
  39. Valentini, A. Signal locality, uncertainty and the subquantum H-theorem. II. Phys. Lett. A 1991, 158, 1–8. [Google Scholar] [CrossRef]
  40. Efthymiopoulos, C.; Kalapotharakos, C.; Contopoulos, G. Origin of chaos near critical points of quantum flow. Phys. Rev. E 2009, 79, 036203. [Google Scholar] [CrossRef] [PubMed]
  41. Valentini, A. Inflationary cosmology as a probe of primordial quantum mechanics. Phys. Rev. D 2010, 82, 063513. [Google Scholar] [CrossRef]
  42. Underwood, N.G.; Valentini, A. Quantum field theory of relic nonequilibrium systems. Phys. Rev. D 2015, 92, 063531. [Google Scholar] [CrossRef]
  43. Jüngel, A. Entropy Methods for Partial Differential Equations; Springer: New York, NY, USA, 2016. [Google Scholar]
  44. Bricmont, J. Bayes, Boltzmann and Bohm: Probabilities in Physics. In Chances in Physics; Bricmont, J., Ghirardi, G., Dürr, D., Petruccione, F., Galavotti, M.C., Zanghi, N., Eds.; Springer: Berlin/Heidelberg, Germany, 2001; pp. 3–21. [Google Scholar]
  45. Gardiner, C.W. Handbook of Stochastic Processes; Springer: New York, NY, USA, 1985. [Google Scholar]
  46. Risken, H. Fokker–Planck Equation. In The Fokker–Planck Equation; Risken, H., Frank, T., Eds.; Springer: Berlin/Heidelberg, Germany, 1996; pp. 63–95. [Google Scholar]
  47. Petroni, N.C.; De Martino, S.; De Siena, S. Exact solutions of Fokker–Planck equations associated to quantum wave functions. Phys. Lett. A 1998, 245, 1–10. [Google Scholar] [CrossRef]
  48. Brics, M.; Kaupuzs, J.; Mahnke, R. How to solve Fokker–Planck equation treating mixed eigenvalue spectrum? Cond. Matter Phys. 2013, 16, 13002. [Google Scholar] [CrossRef]
  49. Hatifi, M.; Willox, R.; Colin, S.; Durt, T. Bouncing oil droplets, de Broglie’s quantum thermostat and convergence to equilibrium. arXiv, 2018; arXiv:quant-ph/1807.00569v3. [Google Scholar]
  50. Gray, R.M. Probability, Random Processes, and Ergodic Properties; Springer: New York, NY, USA, 2009. [Google Scholar]
  51. Arnold, V.; Avez, A. Problèmes Ergodiques de la Méchanique Classique; Gauthier-Villars: Paris, France, 1967. (In French) [Google Scholar]
  52. Labousse, M.; Oza, A.U.; Perrard, S.; Bush, J.W. Pilot-wave dynamics in a harmonic potential: Quantization and stability of circular orbits. Phys. Rev. E 2016, 93, 033122. [Google Scholar] [CrossRef] [PubMed]
  53. Grössing, G. Sub-quantum thermodynamics as a basis of emergent quantum mechanics. Entropy 2010, 12, 1975–2044. [Google Scholar] [CrossRef]
  54. Hestenes, D. The Zitterbewegung Interpretation of Quantum Mechanics. Founds. Phys. 1990, 20, 1213–1232. [Google Scholar] [CrossRef]
  55. Colin, S.; Wiseman, H.M. The zig-zag road to reality. J. Phys. A Math. Theor. 2011, 44, 345304. [Google Scholar] [CrossRef]
  56. Gilet, T. Quantumlike statistics of deterministic wave-particle interactions in a circular cavity. Phys. Rev. E 2016, 93, 042202. [Google Scholar] [CrossRef] [PubMed]
  57. Bechhoefer, J.; Ego, V.; Manneville, S.; Johnson, B. An experimental study of the onset of parametrically pumped surface waves in viscous fluids. J. Fluid Mech. 1995, 288, 325–350. [Google Scholar] [CrossRef]
  58. Cohen-Tannoudji, C.; Diu, B.; Laloe, F. Quantum Mechanics; Wiley: New York, NY, USA, 1977. [Google Scholar]
  59. Hatifi, M.; Lopez-Fortin, C.; Durt, T. de Broglie’s double solution: Limitations of the self-gravity approach. Ann. Fond. Louis Broglie 2018, 43, 63–90. [Google Scholar]
  60. Pucci, G.; Harris, D.M.; Faria, L.M.; Bush, J.W.M. Walking droplets interacting with single and double slits. J. Fluid Mech. 2018, 835, 1136–1156. [Google Scholar] [CrossRef]
  61. Sáenz, P.J.; Cristea-Platon, T.; Bush, J. Statistical projection effects in a hydrodynamic pilot-wave system. Nature 2017, 14, 3. [Google Scholar]
  62. Durt, T. L. de Broglie’s double solution and gravitation. Ann. Fond. Louis Broglie 2017, 42, 73–102. [Google Scholar]
  63. Nelson, E. Review of stochastic mechanics. J. Phys. Conf. Ser. 2012, 361, 012011. [Google Scholar] [CrossRef] [Green Version]
  64. Bell, J. On the Einstein–Podolsky–Rosen paradox. Physics 1964, 1, 195–200. [Google Scholar] [CrossRef]
  65. Durt, T. Characterisation of an entanglement-free evolution. arXiv, 2001; arXiv:quant-ph/0109112. [Google Scholar]
  66. Kostin, M. On the Schrödinger-Langevin Equation. J. Chem. Phys. 1972, 57, 3589–3591. [Google Scholar] [CrossRef]
  67. Nassar, A. Fluid formulation of a generalised Schrodinger-Langevin equation. J. Phys. A Math. General 1985, 18, L509. [Google Scholar] [CrossRef]
  68. Olavo, L.; Lapas, L.; Figueiredo, A. Foundations of quantum mechanics: The Langevin equations for QM. Ann. Phys. 2012, 327, 1391–1407. [Google Scholar] [CrossRef]
  69. Nassar, A.B.; Miret-Artés, S. Bohmian Mechanics, Open Quantum Systems and Continuous Measurements; Springer: Basel, Switzerland, 2017. [Google Scholar]
  70. Bacciagaluppi, G. A Conceptual Introduction to Nelson’s Mechanics. In Endophysics, Time, Quantum and the Subjective; Buccheri, R., Saniga, M., Elitzur, A., Eds.; World Scientific: Singapore, 2005; pp. 367–388. [Google Scholar]
  71. Wallstrom, T. Inequivalence between the Schrödinger equation and the Madelung hydrodynamic equations. Phys. Rev. A 1994, 49, 1613–1617. [Google Scholar] [CrossRef] [PubMed]
  72. Derakhshani, M. A Suggested Answer to Wallstrom’s Criticism (I): Zitterbewegung Stochastic Mechanics. arXiv, 2015; arXiv:quant-ph/1510.06391. [Google Scholar]
  73. Bacciagaluppi, G.; Valentini, A. Quantum Theory at the Crossroads: Reconsidering the 1927 Solvay Conference; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  74. Colin, S.; Durt, T.; Willox, R. de Broglie’s double solution program: 90 years later. Ann. Fond. Louis Broglie 2017, 42, 19–71. [Google Scholar]
  75. Borghesi, C. Equivalent quantum equations with effective gravity in a system inspired by bouncing droplets experiments. arXiv, 2017; arXiv:1706.05640. [Google Scholar]
  76. Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P. Numerical Recipes 3rd Edition: The Art of Scientific Computing, 3rd ed.; Cambridge University Press: New York, NY, USA, 2007. [Google Scholar]
  77. Higham, D.J. An algorithmic introduction to numerical simulation of stochastic differential equations. SIAM Rev. 2001, 43, 525–546. [Google Scholar] [CrossRef]
Figure 1. A particle suspended in a Madelung fluid and subject to local fluctuations.
Figure 1. A particle suspended in a Madelung fluid and subject to local fluctuations.
Entropy 20 00780 g001
Figure 2. Simulations of 10,000 trajectories (calculated from the Ito equation, Equation (16), for the ground state (Equation (55)) of the 1D harmonic oscillator), whose initial positions are normally distributed, for 5 different choices of distribution width (for a = 0.5 and α = 1 ). We observe, in each case, convergence to the equilibrium of Equation (55) as predicted by the theory.
Figure 2. Simulations of 10,000 trajectories (calculated from the Ito equation, Equation (16), for the ground state (Equation (55)) of the 1D harmonic oscillator), whose initial positions are normally distributed, for 5 different choices of distribution width (for a = 0.5 and α = 1 ). We observe, in each case, convergence to the equilibrium of Equation (55) as predicted by the theory.
Entropy 20 00780 g002
Figure 3. Time evolution of H V , L f and the L 1 norm, for a uniform initial probability distribution, calculated from the Ito equation, Equation (16), for the ground state of the 1D harmonic oscillator. Relaxation towards the distribution of the ground state | Ψ s t | 2 of Equation (55) is clearly visible. The simulation was performed for α = 1 , a = 0.5 , and Δ t = 0.01 , for 20,000 uniformly distributed initial conditions.
Figure 3. Time evolution of H V , L f and the L 1 norm, for a uniform initial probability distribution, calculated from the Ito equation, Equation (16), for the ground state of the 1D harmonic oscillator. Relaxation towards the distribution of the ground state | Ψ s t | 2 of Equation (55) is clearly visible. The simulation was performed for α = 1 , a = 0.5 , and Δ t = 0.01 , for 20,000 uniformly distributed initial conditions.
Entropy 20 00780 g003
Figure 4. Histograms of the positions of a single particle, subject to Nelson dynamics for the ground state of the 1D harmonic oscillator. The full (red) curve corresponds to the quantum probability | Ψ s t | 2 . Here, a = 0.5 , α = 1 , and the total simulation time (t = 10,000) is sampled with Δ t = 0.01 . ( a ) The initial particle position is x 0 = 2.5 , and the number of bins N b = 100 (each with spatial size Δ x = 0.0635 ); ( b ) Same as (a) but with N b = 50 and Δ x = 0.1270 ; ( c ) Same as (a) but with t = 200 ; ( d ) Same as (a) but for x 0 = 0.85 .
Figure 4. Histograms of the positions of a single particle, subject to Nelson dynamics for the ground state of the 1D harmonic oscillator. The full (red) curve corresponds to the quantum probability | Ψ s t | 2 . Here, a = 0.5 , α = 1 , and the total simulation time (t = 10,000) is sampled with Δ t = 0.01 . ( a ) The initial particle position is x 0 = 2.5 , and the number of bins N b = 100 (each with spatial size Δ x = 0.0635 ); ( b ) Same as (a) but with N b = 50 and Δ x = 0.1270 ; ( c ) Same as (a) but with t = 200 ; ( d ) Same as (a) but for x 0 = 0.85 .
Entropy 20 00780 g004
Figure 5. (a) A point-particle (the dot near the center) subject to the osmotic velocity field 2 a α ( x ( t ) , y ( t ) ) , due to the ground state of the 2D harmonic oscillator at time t; (b) Color plot of the velocities along a trajectory for the evolution under Nelson dynamics, for the ground state of the 2D harmonic oscillator. The simulation (for a = 0.5 and α = 1 ) started from the initial position ( 2 , 1 ) and was sampled up to t = 1000 with step Δ t = 0.01 .
Figure 5. (a) A point-particle (the dot near the center) subject to the osmotic velocity field 2 a α ( x ( t ) , y ( t ) ) , due to the ground state of the 2D harmonic oscillator at time t; (b) Color plot of the velocities along a trajectory for the evolution under Nelson dynamics, for the ground state of the 2D harmonic oscillator. The simulation (for a = 0.5 and α = 1 ) started from the initial position ( 2 , 1 ) and was sampled up to t = 1000 with step Δ t = 0.01 .
Entropy 20 00780 g005
Figure 6. Histogram of the positions in x of a single particle, in the case of the first Fock state given by Equation (68). The full curve (red) corresponds to the quantum probability | Ψ 1 | 2 . Here, a = 0.5 and α = 1 . The total simulation time t is t = 1000 , and the sampling time step is Δ t = 0.01 . The initial position is x i = 1 , and the number of bins N b = 75 , each with width Δ x = 0.08 .
Figure 6. Histogram of the positions in x of a single particle, in the case of the first Fock state given by Equation (68). The full curve (red) corresponds to the quantum probability | Ψ 1 | 2 . Here, a = 0.5 and α = 1 . The total simulation time t is t = 1000 , and the sampling time step is Δ t = 0.01 . The initial position is x i = 1 , and the number of bins N b = 75 , each with width Δ x = 0.08 .
Entropy 20 00780 g006
Figure 7. Evolution in time of H V (27), L f (28), and of the L 1 (39) norm, for a uniform initial probability distribution, showing the relaxation towards the distribution of the first excited state | Ψ 1 | 2 (68). The simulation was performed for α = 1 , a = 0.5 , and Δ t = 0.01 and from 20,000 uniformly distributed initial conditions.
Figure 7. Evolution in time of H V (27), L f (28), and of the L 1 (39) norm, for a uniform initial probability distribution, showing the relaxation towards the distribution of the first excited state | Ψ 1 | 2 (68). The simulation was performed for α = 1 , a = 0.5 , and Δ t = 0.01 and from 20,000 uniformly distributed initial conditions.
Entropy 20 00780 g007
Figure 8. (a) The quantum probability associated to the Bessel function of the first kind J 0 ; (b) Color plot of the velocities reached along the trajectory for an evolution corresponding to (69). The initial position was ( 1 , 1 ) , the simulation time t = 5000 , and the sampling time step Δ t = 0.005 . We chose α = 0.1 and ϵ = 0.2 , and the size of the domain is L = 2 . On the boundary we impose a harmonic field force of the form 2 a α r .
Figure 8. (a) The quantum probability associated to the Bessel function of the first kind J 0 ; (b) Color plot of the velocities reached along the trajectory for an evolution corresponding to (69). The initial position was ( 1 , 1 ) , the simulation time t = 5000 , and the sampling time step Δ t = 0.005 . We chose α = 0.1 and ϵ = 0.2 , and the size of the domain is L = 2 . On the boundary we impose a harmonic field force of the form 2 a α r .
Entropy 20 00780 g008
Figure 9. Numerical solutions of the Ito stochastic differential equation, Equation (16), corresponding to the coherent state of Equation (70), for three different initial conditions. We used x ¯ 0 = 1 , a = 0.5 , and α = 1 and expressed the results in natural units.
Figure 9. Numerical solutions of the Ito stochastic differential equation, Equation (16), corresponding to the coherent state of Equation (70), for three different initial conditions. We used x ¯ 0 = 1 , a = 0.5 , and α = 1 and expressed the results in natural units.
Entropy 20 00780 g009
Figure 10. The time evolution of a non-equilibrium ensemble, illustrated with position histograms at six different times ((a): t = 0 , (b): t = 1.2 , (c): t = 2.4 , (d): t = 3.6 , (e): t = 4.8 , (f): t = 6 ). The continuous curve is the squared modulus | Ψ | 2 for the coherent state of Equation (70). As can be seen in (df), once equilibrium is reached, the distribution clings to the coherent state and follows its oscillation faithfully. The center of the wave packet moves between 2 and 2 with a period 2 π . We started from a uniform distribution of initial conditions and chose a = 0.5 , α = 1 , and x 0 = 2 . The sampling time step is Δ t = 0.01 , and the number of bins is N b = 50 , each with width Δ x = 0.0461 .
Figure 10. The time evolution of a non-equilibrium ensemble, illustrated with position histograms at six different times ((a): t = 0 , (b): t = 1.2 , (c): t = 2.4 , (d): t = 3.6 , (e): t = 4.8 , (f): t = 6 ). The continuous curve is the squared modulus | Ψ | 2 for the coherent state of Equation (70). As can be seen in (df), once equilibrium is reached, the distribution clings to the coherent state and follows its oscillation faithfully. The center of the wave packet moves between 2 and 2 with a period 2 π . We started from a uniform distribution of initial conditions and chose a = 0.5 , α = 1 , and x 0 = 2 . The sampling time step is Δ t = 0.01 , and the number of bins is N b = 50 , each with width Δ x = 0.0461 .
Entropy 20 00780 g010
Figure 11. Time evolution of H V (27), L f (28), and L 1 (39), for a uniform initial probability distribution, showing the relaxation towards the distribution | Ψ | 2 of the coherent state of Equation (70). The simulation was performed for α = 1 , a = 0.5 , and Δ t = 0.01 and from 20,000 uniformly distributed initial conditions.
Figure 11. Time evolution of H V (27), L f (28), and L 1 (39), for a uniform initial probability distribution, showing the relaxation towards the distribution | Ψ | 2 of the coherent state of Equation (70). The simulation was performed for α = 1 , a = 0.5 , and Δ t = 0.01 and from 20,000 uniformly distributed initial conditions.
Entropy 20 00780 g011
Figure 12. Histogram of the positions for a unique trajectory satisfying the Ito equation, Equation (16), for Equation (70). The full curve corresponds to the integration of | Ψ | 2 over one period. The center of the wave packet moves between 2 and 2 with a period 2 π . Here, a = 0.5 and α = 1 . Total simulation time t is t = 30 , 000 and the sampling time step is Δ t = 0.01 . The initial position is x i = 1 , and the number of bins N b = 100 , each with width Δ x = 0.1 .
Figure 12. Histogram of the positions for a unique trajectory satisfying the Ito equation, Equation (16), for Equation (70). The full curve corresponds to the integration of | Ψ | 2 over one period. The center of the wave packet moves between 2 and 2 with a period 2 π . Here, a = 0.5 and α = 1 . Total simulation time t is t = 30 , 000 and the sampling time step is Δ t = 0.01 . The initial position is x i = 1 , and the number of bins N b = 100 , each with width Δ x = 0.1 .
Entropy 20 00780 g012
Figure 13. Plots showing three possible de Broglie–Bohm (dBB) trajectories for a single point particle in the case of (73) with M = 2 2 = 4 . Each plot (ac) is associated with different initial random phases ( θ n x , n y with n x ( n y ) taking the values 0, 1) and different initial positions.
Figure 13. Plots showing three possible de Broglie–Bohm (dBB) trajectories for a single point particle in the case of (73) with M = 2 2 = 4 . Each plot (ac) is associated with different initial random phases ( θ n x , n y with n x ( n y ) taking the values 0, 1) and different initial positions.
Entropy 20 00780 g013
Figure 14. Plots of the evolution in time of the coarse-grained H-functions H V (a) and L 1 (b) for the Nelson and dBB dynamics. The full line corresponds to the dBB dynamics and the dashed line corresponds to the quantum thermostat. We started from 10,000 initial positions uniformly distributed in a box of size 10 × 10; we chose a = 0.5 , α = 0.1 , and M = 4 2 = 16 energy states.
Figure 14. Plots of the evolution in time of the coarse-grained H-functions H V (a) and L 1 (b) for the Nelson and dBB dynamics. The full line corresponds to the dBB dynamics and the dashed line corresponds to the quantum thermostat. We started from 10,000 initial positions uniformly distributed in a box of size 10 × 10; we chose a = 0.5 , α = 0.1 , and M = 4 2 = 16 energy states.
Entropy 20 00780 g014
Figure 15. dBB trajectories obtained for a single point particle in a superposition of eigenstates (75). Each plot is associated with a different combination ( n , m ) , as indicated. In the (a,b) graphs, we imposed a = 1 and, respectively, ω = 1 , ξ 0 ξ 2 = 0.05 and ω = 0.5 , ξ 0 ξ 3 = 0.05 ; for (c), we imposed a = 3 , ω = 0.5 , ξ 0 ξ 3 = 0.0708 , ξ 0 ξ 2 = 0.0456 , and ξ 0 ξ 1 = 0.0773 .
Figure 15. dBB trajectories obtained for a single point particle in a superposition of eigenstates (75). Each plot is associated with a different combination ( n , m ) , as indicated. In the (a,b) graphs, we imposed a = 1 and, respectively, ω = 1 , ξ 0 ξ 2 = 0.05 and ω = 0.5 , ξ 0 ξ 3 = 0.05 ; for (c), we imposed a = 3 , ω = 0.5 , ξ 0 ξ 3 = 0.0708 , ξ 0 ξ 2 = 0.0456 , and ξ 0 ξ 1 = 0.0773 .
Entropy 20 00780 g015
Figure 16. Plots of three quantities associated to the lemniscate in Figure 15c. ( a ) shows the L z -component of the angular momentum, and the polar plots ( b , c ) show the probability density | ψ | 2 (b) and the θ -component of the probability current (5) along the trajectory (c).
Figure 16. Plots of three quantities associated to the lemniscate in Figure 15c. ( a ) shows the L z -component of the angular momentum, and the polar plots ( b , c ) show the probability density | ψ | 2 (b) and the θ -component of the probability current (5) along the trajectory (c).
Entropy 20 00780 g016
Figure 17. dBB trajectories obtained for a single point particle in a superposition of eigenstates (75). Plots (a,b) correspond to ω = 0 , 7 and ω = 1 , respectively. Case (c) is obtained after multiplying the amplitude of the ( n , m ) = ( 4 , 2 ) state by a complex phase ( e ( 0.3 i ) ). We took a = 1 in all cases.
Figure 17. dBB trajectories obtained for a single point particle in a superposition of eigenstates (75). Plots (a,b) correspond to ω = 0 , 7 and ω = 1 , respectively. Case (c) is obtained after multiplying the amplitude of the ( n , m ) = ( 4 , 2 ) state by a complex phase ( e ( 0.3 i ) ). We took a = 1 in all cases.
Entropy 20 00780 g017
Figure 18. dBB trajectories obtained for a single point particle in a superposition of eigenstates (75) (with m = 2 , 0 , 2 ) showing intermittent transitions between two types of trajectories. The relevant parameter values are ω = 0.2 , a = 1 , ξ 0 ξ 3 = 0.0342 , ξ 0 ξ 2 = 0.2547 , and ξ 0 ξ 1 = 0.0505 .
Figure 18. dBB trajectories obtained for a single point particle in a superposition of eigenstates (75) (with m = 2 , 0 , 2 ) showing intermittent transitions between two types of trajectories. The relevant parameter values are ω = 0.2 , a = 1 , ξ 0 ξ 3 = 0.0342 , ξ 0 ξ 2 = 0.2547 , and ξ 0 ξ 1 = 0.0505 .
Entropy 20 00780 g018

Share and Cite

MDPI and ACS Style

Hatifi, M.; Willox, R.; Colin, S.; Durt, T. Bouncing Oil Droplets, de Broglie’s Quantum Thermostat, and Convergence to Equilibrium. Entropy 2018, 20, 780. https://doi.org/10.3390/e20100780

AMA Style

Hatifi M, Willox R, Colin S, Durt T. Bouncing Oil Droplets, de Broglie’s Quantum Thermostat, and Convergence to Equilibrium. Entropy. 2018; 20(10):780. https://doi.org/10.3390/e20100780

Chicago/Turabian Style

Hatifi, Mohamed, Ralph Willox, Samuel Colin, and Thomas Durt. 2018. "Bouncing Oil Droplets, de Broglie’s Quantum Thermostat, and Convergence to Equilibrium" Entropy 20, no. 10: 780. https://doi.org/10.3390/e20100780

APA Style

Hatifi, M., Willox, R., Colin, S., & Durt, T. (2018). Bouncing Oil Droplets, de Broglie’s Quantum Thermostat, and Convergence to Equilibrium. Entropy, 20(10), 780. https://doi.org/10.3390/e20100780

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