Emergent classicality in continuous quantum measurements

We develop a classical theoretical description for nonlinear many-body dynamics that incorporates the back-action of a continuous measurement process. The classical approach is compared with the exact quantum solution in an example with an atomic Bose-Einstein condensate in a double-well potential where the atom numbers in both potential wells are monitored by light scattering. In the classical description the back-action of the measurements appears as diffusion of the relative phase of the condensates on each side of the trap. When the measurements are frequent enough to resolve the system dynamics, the system behaves classically. This happens even deep in the quantum regime, and demonstrates how classical physics emerges from quantum mechanics as a result of measurement back-action.

It is now possible to carry out experiments with systems comprised of a small and controllable number of atoms, and eventually the interface between quantum mechanics and classical physics [1,2,3] must be confronted. Measurement backaction is an essential element of quantum mechanics, and may play a key role also in the transition from quantum mechanics to classical physics [4,5,6,7]. In this paper we introduce a theoretical method to account for measurement back-action in nonlinear dynamics [8,9,10,11,12,13] classically, using a continuously monitored [14] Bose-Einstein condensate (BEC) in a double-well trap as an explicit example. When the continuous measurements are strong enough so that the system dynamics can be resolved, the classical description and a full quantum analysis agree well even deep within what might appear to be the domain of quantum mechanics.
We employ the minimal two-mode model [15] to describe the Josephson oscillations of the atoms between the potential wells [16,17,18,19,20], see figure 1. Given the hopping amplitude J between the wells, the strength of atom-atom interaction inside a well U, and the imbalance in the energies of the atoms between the sides ǫ, the Hamiltonian is Here a and b are the boson operators for the atoms in the two wells, and we set = 1. The atom numbers in the two wells are monitored by coherent scattering of offresonant light; consult figure 1. Counting of photons on two detectors, one for each potential well, constitutes the final measurement. In our quantum description we consider the atoms as an open system, and model the back-action of the measurements J U ε Figure 1. Left: Bose-Einstein condensate in a double-well potential with tunneling amplitude J, atom-atom interaction strength U , and energy difference per atom between the wells ǫ indicated. Right: Basic idea of detection of the atoms in one trap. Since the detection light is far-off resonance, the dominant effect of the atoms is to alter the phase of the light passing through the sample by an amount proportional to atom number. The phase shift is then detected interferometrically. The amplitude of the light passing to the photodetector is proportional to the number of atoms, so that the rate of photon counts is proportional to the square of the atom number. The scheme is heavily simplified. For instance, mode matching of the light to the condensate is not considered.
by adding a linear relaxation term to the equation of motion of the density operator of the atoms,ρ The most general relaxation term Lρ such that a density operator remains valid during the time evolution is of the Lindblad form [21,22], where L i are operators of the double-well system. In the present model [14] we have the indicating that the scattering rates of photons from the two wells are Γ (a † a) 2 and Γ (b † b) 2 , respectively. If the light is predominantly scattered into a narrow cone in the forward direction, Γ also specifies the strength of the measurements. It can be adjusted, for instance, by varying the intensity or the frequency of the probing light.
One could develop a parallel classical description using various phase-space representations such as the Wigner function for the ordinary [23,24,22] or the SU(2) [25,26,27] coherent states, or the number-phase Wigner representation [28]. We choose the standard Wigner function because it is generally used in quantum optics and is not specifically tailored to the present system. We begin by rewriting the master equation (2) equivalently in terms of the Wigner function W (α, α * , β, β * ; t) [24,22], where α and β are complex numbers corresponding to the operators a and b. In the equation of motion of the Wigner function, α, α * and β, β * behave like canonical pairs with Poisson brackets such as {α, α * } = −i. We then change to four real variables N, Φ, Z and ϕ defined by follows for the remaining distribution function W (z, ϕ), with the effective one-particle Hamiltonian Importantly, the quantum mechanical back-action of the measurements of the particle numbers is now accounted for by classical diffusion of the relative phase of the BECs proportional to Γ∂ ϕϕ . This is, perhaps, as expected: The relevant measurement result is the difference of the photon counting rates on the two detectors, and this is proportional to the expectation value of the operator But the first factor in the final form is the conserved total atom number, so the measured operator is actually the difference of atom numbers between the sides of the trap. Given the (in quantum mechanics, qualitative) conjugate relation between number and phase, measurement of atom number difference should indeed increase the uncertainty of the relative phase. We compare the quantum and classical evolutions in the case when the system is first prepared to the ground state at zero temperature in the presence of an energy imbalance, ǫ = 0, and the imbalance is suddenly removed by setting ǫ = 0. The atoms then start oscillating between the two sides of the potential well.
In our quantum analysis the evolution of the density operator according to (2) is solved using quantum trajectory simulations [29,30,31]. For every run of the simulation we have a state vector that undergoes 'quantum jumps' at random times according to a precisely defined probabilistic law. As discussed in, say, Refs. [22,32], here we also assume an "efficient" measurement process in which every scattered photon gets detected. Then a quantum trajectory simulation with its quantum jumps produces a faithful simulation of an individual experimental run, each quantum jump corresponding to a detection of a photon. Even though the master equation (2) is not conditioned on the results of the measurements and as such describes an experiment in which the measurement results are discarded, the results are built into each quantum trajectory.
In a real experiment one would record a stream of photon counts, and the population imbalance z(t) has to be inferred from them in some manner. For instance, if frequent photon counts merge into an essentially continuous photocurrent that varies on a time scale of the back-and-forth oscillations of the atoms, it would be sensible to regard the difference of the photocurrents coming from the two detectors as a measure of the instantaneous population imbalance. However, for the sake of a clean illustration we analyze our quantum simulations differently: Given a quantum trajectory, a time dependent state vector, we use the expectations value ofẑ = (a † a − b † b)/N to represent the value of the experimentally inferred population imbalance z(t). Evidently this works the better, the more numerous the quantum jumps are on the time scale of the variation of z(t).
In the classical analysis we take into account the quantum fluctuations in the initial state in the same way as is customary with the truncated Wigner approximation (TWA) [33,34,35,36]. Thus, suppose equations (5) and (6) were written in terms of the variables Z and ϕ instead of z and ϕ. We expand the corresponding Hamiltonian H(Z, ϕ) to second order in Z and ϕ around the classical lowest-energy state in the presence of the initial energy imbalance ǫ = 0 between the sides. Since the classical variables Z and ϕ are a canonical pair, we may requantize by postulating the commutator [Z, ϕ] = i.
H(Z, ϕ) thereby becomes the quantum Hamiltonian of a harmonic oscillator. The Wigner function corresponding to the ground state of this oscillator is used as the initial distribution of the variables Z and ϕ, and leads to the initital condition for the distribution of the variables z and ϕ We take χ > 0 and choose ǫ in such a way that the classical equilibrium has the phase ϕ = 0 and a predetermined population imbalance z = f . The parameter reflects the size of the ground state wave function of the underlying harmonic oscillator. The evolution (5) with the initial Wigner function W 0 (z, ϕ) is solved by unravelling the dynamics into stochastic trajectories [33,34,37,26]. If there were no phase diffusion in (5), we would have an ensemble of trajectories (z(t), ϕ(t)) evolving under the effective Hamiltonian h, which makes the usual TWA. Given the diffusion, there exists a precise mathematical correspondence between (5) and the stochastic differential equations where dW (t) denotes a random increment of the Wiener process over the interval dt [37]. The evolution of the distribution W (z, ϕ; t) from the initial condition W 0 (z, ϕ) according to (5) is obtained from an ensemble of stochastic trajectories. Specifically, the initial state of each trajectory (z, ϕ) is sampled from the distribution W 0 (z, ϕ), the dynamics of the trajectory (z(t), ϕ(t)) is determined from (10) and (11), and the expectation value of any quantity over the distribution W (z, ϕ; t) is correctly obtained by averaging the value of the quantity over a large number of such trajectories. The stochastic dynamics is integrated numerically using the Milstein algorithm [37].
We are now in a position to state precisely what "classical description" means in this context: unlike quantum mechanics, the analysis conforms to classical logic. The usual mean-field theory [17] is a classical description, as are the TWA and our present approach. The initial Wigner function W 0 is a valid classical probability distribution, even when it approximately incorporates quantum fluctuations of the atomic populations in the initial state. Moreover, subsequent evolution according to (5) that also accounts for measurement back-action keeps the probability distribution function valid at all future times. Correspondingly, we have been able to represent the complete solution in terms of classical stochastic trajectories. Such a trajectory exists objectively even without observations, but otherwise a similar correspondence between simulations and experiments prevails in the classical description as in quantum mechanics: each individual stochastic realization would be a representative outcome of an individual experiment.  Figure 2 shows a classical (blue upper trace) and a quantum rendition (red lower trace) of a single experiment. In the mean-field theory, and for each individual TWA trajectory, one expects strictly periodic oscillations in the classical trace. The deviations from periodicity originate from the phase diffusion. On the other hand, the measurement strength Γ/J = 0.003 is chosen such that each photon detector produces about ten clicks during a typical back-and-forth swing of the atoms, so that the quantum trace is a fair representation of a measurement record. We have deliberately picked these random examples so that qualitatively similar classical and quantum behavior with occasional trapping of the populations [17] in one or the other well is seen.
Nevertheless, even if we are able to model individual experiments both classically and quantum mechanically, these are stochastic stimulations. In order to construct expectation values of experimental observables for quantitative comparisons we have to average over many trajectories. In figure 3 we demonstrate simulation results for population imbalance when averaged over a thousand trajectories. The upper two panels are for N = 100 atoms and the lower two for N = 10, while the differences between the graphs (a)-(d) are in the measurement strength Γ.
Panel (a) for N = 100 atoms plots the quantum and classical expectation values of the population imbalance z(t) with no measurements, Γ = 0. The quantum result then represents the usual quantum average of the operatorẑ = (a † a − b † b)/N, and has nothing to do with the (nonexistent) continuous monitoring. The quantum and classical graphs both show a collapse up to Jt ≃ 10, but only the quantum expectation value has a revival at around Jt = 70. In panel (b) we turn up the measurements strength to Γ/J = 3 × 10 −3 . Both detectors then report about ten photons during a characteristic oscillation time of the atoms between the traps, so that the oscillations of the atoms could be resolved by monitoring the light scattered off the atoms. The quantum expectation value ofẑ should therefore be a passable representation of the population imbalance inferred from photon counts. The collapse part, quantum and classical, is almost unchanged by the measurements, but the quantum revival has disappeared. The quantum and classical analyses now agree very well.
We finally move on to the atoms number N = 10. In panel (c) we show the  population imbalances as a function of time for both the quantum theory and the classical theory without continuous monitoring, Γ = 0. There is little in common with the two traces, and one would conclude that the classical description is totally off. The proximate reason is that with a decreasing atom number the quantum revivals start overlapping and produce irregular-looking behavior that is absent in the classical approach. It appears that we are deep in the domain of quantum mechanics. However, in panel (d) we have continuous measurements with the strength Γ/J = 0.3, which again gives about ten photon counts during a back-and forth oscillation of the atoms; and again, the classical and quantum results agree.
We have surveyed the behavior of the system for parameter values comparable to those in figures 2 and 3. In the absence of measurements, Γ = 0, we have the usual TWA and the corresponding full quantum theory. The TWA relies on a stochastic description and cannot fully incorporate quantum effects. The generic observation is that TWA agrees well with the exact solution for short times when stochastic field is highly occupied, but becomes progressively less accurate a model of quantum dynamics when the average mode population becomes smaller and the simulation time scales longer [33,34,35]. There are also examples, however, especially in lower dimensional systems in which cases the TWA simulations have been successful in providing qualitative descriptions of experimental findings even in the strongly fluctuating quantum limit [38].
In figure 3, in traces (a) with N = 100 and Γ = 0 we find a reasonable agreement between the TWA and quantum solutions up to the quantum revival around Jt = 70, while in figure 3(c) with N = 10 the two cases part ways after one back-and forth swing of the atoms at Jt ≃ 1, a time that is hard to discern in the figure at all. Second, as the measurement strength is increased from zero, the changes in the quantum results start being qualitatively visible at times when there have been, on the average, a few counts on both detectors. The classical results also change with increasing Γ, but the changes are less conspicuous. For example, with N = 100 the maximum amplitude of the quantum revival at Jt ≃ 70 has decreased to half of its value in figure 3(a) for Γ/J = 2×10 −5 , at which point the average time between photon counts at each detector is about 20/J. Let us next consider Γ/J = 3 × 10 −3 , the value used in figure 3(b) for N = 100. Since our definition of the parameter Γ is such that the photon counting rate on the detectors is proportional to the square of the atom number, for N = 10 this measurement strength is still modest, with the time between counts of about 10/J. A plot (not shown) would reveal that by Jt = 100 the oscillations in the quantum result have lost 2/3 of their amplitude, but that the quantum and classical approaches still give quite different results.
Finally, as soon as the atom numbers are monitored vigorously enough so that the classical evolution in the form of the oscillations of the atoms between the two sides of the double-well trap can be resolved in detail in the scattered light, figures 3(b) and (d), the quantum and classical descriptions give materially the same results. In round numbers, the time scale for classical evolution in our examples is 1/J and the scattering rate to both detectors is N 2 Γ/4. The watershed criterion is therefore N 2 Γ/4J ∼ 1, exceeded by a factor of about 10 in both graphs 3(b) and (d). Moreover, and unlike in the usual case with the ordinary TWA, in our amended version of the TWA the quantumclassical agreement seems to persist for a long time, possibly as long as our numerical computations remain reliable.
The emergent classicality is a consequence of measurement back-action suppressing the subtle interference effects associated with quantum revivals that are inherent in the full quantum description [3]. Remarkably, this occurs even deep in the domain of quantum mechanics with only 10 atoms, representing a far more dramatic transition to classicality than, e.g., measurement-induced relative phase between condensates [5,6,7] or nonlinear instabilities [10,12,13,39]. The classical description at this level of accuracy is only possible because we have included classical phase diffusion to account for quantum mechanical measurement back-action. To wit, the only difference between the classical predictions in graphs 3(c) and (d) is the amount of phase diffusion, yet they are clearly different.
The main technical achievement of this paper is that we have introduced a theoretical method for a continuously monitored quantum system that handles both quantum fluctuations and measurement back-action classically. Now, the computational resources required to solve the time evolution from quantum mechanics tend to increase very fast with the system size. We therefore hope to inspire practical algorithms to analyze time evolution in situations where quantum effects are small but not negligible.