We propose physical interpretations, also valid for temperatures different from zero, for stochastic methods which have been developed recently to describe the evolution of a quantum system interacting with a reservoir. As opposed to the usual reduced density operator approach, which refers to ensemble averages, these methods deal with the dynamics of single realizations, and involve the solution of stochastic Schrödinger equations. These procedures have been shown to be completely equivalent to the master equation approach when ensemble averages are taken over many realizations. We show that these techniques are not only convenient mathematical tools for dissipative systems, but may actually correspond to concrete physical processes, for any temperature of the reservoir. We consider a mode of the electromagnetic field in a cavity interacting with a beam of two- or three-level atoms, the field mode playing the role of a small system and the atomic beam standing for a reservoir at finite temperature, the interaction between them being given by the Jaynes-Cummings model. We show that the evolution of the field states, under continuous monitoring of the state of the atoms which leave the cavity, can be described in terms of either the Monte Carlo wavefunction (quantum jump) method or a stochastic Schrödinger equation, depending on the system configuration. We also show that the Monte Carlo wavefunction approach leads, for finite temperatures, to localization into jumping Fock states, while the diffusion equation method leads to localization into states with a diffusing average photon number, which for sufficiently small temperatures are close approximations to mildly squeezed states. We prove analytically that, in the quantum jump situation, the system evolves in the mean towards a Fock state, even if an infinite number of photon-number amplitudes is present in the initial state.
We propose physical interpretations, also valid for temperatures different from zero, for stochastic methods which have been developed recently to describe the evolution of a quantum system interacting with a reservoir. As opposed to the usual reduced density operator approach, which refers to ensemble averages, these methods deal with the dynamics of single realizations, and involve the solution of stochastic Schrödinger equations. These procedures have been shown to be completely equivalent to the master equation approach when ensemble averages are taken over many realizations. We show that these techniques are not only convenient mathematical tools for dissipative systems, but may actually correspond to concrete physical processes, for any temperature of the reservoir. We consider a mode of the electromagnetic field in a cavity interacting with a beam of two- or three-level atoms, the field mode playing the role of a small system and the atomic beam standing for a reservoir at finite temperature, the interaction between them being given by the Jaynes-Cummings model. We show that the evolution of the field states, under continuous monitoring of the state of the atoms which leave the cavity, can be described in terms of either the Monte Carlo wavefunction (quantum jump) method or a stochastic Schrödinger equation, depending on the system configuration. We also show that the Monte Carlo wavefunction approach leads, for finite temperatures, to localization into jumping Fock states, while the diffusion equation method leads to localization into states with a diffusing average photon number, which for sufficiently small temperatures are close approximations to mildly squeezed states. We prove analytically that, in the quantum jump situation, the system evolves in the mean towards a Fock state, even if an infinite number of photon-number amplitudes is present in the initial state.
A general description of the Quantum Trajectories is presented, both from the Montecarlo and Stochastic Schrodinger equation viewpoint. Also a physical interpretation of these trajectories is given. The example is a mode of the electromagnetic field in a cavity interacting with a beam of two or three level atoms, the field playing the role of the small system and the atomic beam standing for a reservoir at finite temperature. We show that the evolution of the field states, under continuous monitoring of the state of the atom which leave the cavity can be described with the Monte-Carlo method or a diffusion equation, depending on the configuration of the system.
We study the stochastic dynamics of the electromagnetic field in a lossless cavity interacting with a beam of two-level atoms, given that the atomic states are measured after they have crossed the cavity. The atoms first interact at the exit of the cavity with a classical laser field ℰ and then enter into a detector which measures their states. Each measurement disentangles the field and the atoms and changes in a random way the state |ψ(t)〉 of the cavity field. For weak atom-field coupling, the evolution of |ψ(t)〉 when many atoms cross the cavity and the detector is characterized by a succession of quantum jumps occurring at random times, separated by quasi-Hamiltonian evolutions, both of which depend on the laser field ℰ. For ℰ=0, the dynamics is the same as in the Monte Carlo wave function model of Dalibard et al. [Phys. Rev. Lett. 68, 580 (1992)] and Carmichael, An Open System Approach to Quantum Optics, Lecture Notes in Physics Vol. 18 (Springer, Berlin, 1991)]. The density matrix of the quantum field, obtained by averaging the projector |ψ(t)〉〈ψ(t)| over all results of the measurements, is independent of ℰ and follows the master equation of the damped harmonic oscillator at finite temperature. We provide numerical evidence showing that for large ℰ, an arbitrary initial field state |ψ(0)〉 evolves under the monitoring of the atoms and the measurements toward squeezed states |α,re2iφ〉, moving in the α-complex plane but with almost constant squeezing parameters r and φ. The values of r and φ are determined analytically. On the other hand, for ℰ=0, the dynamics transforms the initial state into Fock states |n〉 with fluctuating numbers of photons n, as shown in Kist et al. [J. Opt. B: Quantum Semiclassical Opt. 1, 251 (1999)]. In the last part, we derive the quantum jump dynamics from the linear quantum jump model proposed in Spehner and Bellissard [J. Stat. Phys. 104, 525 (2001)], for arbitrary open quantum systems having a Lindblad-type evolution. A careful derivation of the infinite jump rates limit, where the dynamics can be approximated by a diffusion process of the quantum state, is also presented.
Very recently there have been a number of papers discussing the applicability of stochastic methods as an alternative to the usual treatment of master equations in quantum mechanics.1-5 In this paper we present a wave function approach including a stochastic element to study the evolution of a small system when it is coupled to a reservoir. We show how fluctuations and dissipation originate from quantum jumps occuring randomly during the time evolution of the system, and we present a physical interpretation of the procedure. We demonstrate how this treatment may replace a wide class of master equations, and its use in problems with a large number of states N is discussed. Since the density matrix has N ×N elements, where a wave function is determined by only N variables, there is a computational advantage in using the wave function procedure. This is illustrated by calculations on cooling atoms in a laser field.
We study the time evolution of the quantum field inside a cavity coupled to a beam of two-level atoms of temperature T , given that each atom, after having crossed the cavity, interacts with a classical field and finally with a detector measuring its state. It is found that, if the coupling between the atoms and the quantum field is weak and is not too small, for any given realization of the measurements, an arbitrary initial state of the field localizes after some time into squeezed states. The centre α of the squeezed state moves randomly in time in the complex plane, but the squeezing amplitude r and phase ϕ show very small fluctuations. Their mean values and are independent of the random results of the measurements, of the initial state and of the atom–field coupling constant λ. The time evolution of r and ϕ is determined analytically by deriving and solving the quantum state diffusion equation describing the field dynamics in the limit of small λ, keeping finite. It is shown that increases with T , i.e., the squeezing is enhanced by increasing the temperature of the atomic beam.
The decay of coherence when a quantum system interacts with a much larger environment is usually described by a master equation for the system reduced density matrix and emphasizes the evolution of an entire ensemble. We consider two methods that have been developed recently to simulate the evolution of single realizations. Quantum-state diffusion involves both diffusion, where the individual quantum trajectory fluctuates through a Wiener process deriving from the environment, and localization to a coherent state, an eigenstate of the relevant Lindblad operator describing the coupling of the system to the environment. We demonstrate the localization process for different initial states and utilize the Wigner function to depict this localization in phase space. We concentrate on quantum states that can be expressed as a superposition of appropriate coherent states. For an initial superposition of two coherent states (a Schr\"odinger ``cat''), one of the two components will dominate the evolution. For initial Fock states, which can be described as a continuous superposition of coherent states on a ring, localization takes place when one coherent state is selected from that ring where each component has nearly the same energy as the original Fock state. We also consider the localization from a nonclassical squeezed ground state, which can be expressed as a superposition of coherent states along a line in phase space. The second simulation method considered is the state vector Monte Carlo, or ``quantum jump,'' approach, which relates to the direct counting of decay quanta. In the case of an initial Schr\"odinger ``cat,'' we find that when no quantum is detected the ``cat'' shrinks, but when a quantum is detected, the Schr\"odinger ``cat'' ``jumps'' from one type of ``cat'' to another with different internal phase. For an initial squeezed state we show how quantum jumps lead to individual realizations which are superpositions of two squeezed states.
Dissipation, the irreversible loss of energy and coherence, from a microsystem is the result of coupling to a much larger macrosystem (or reservoir) that is so large that one has no chance of keeping track of all of its degrees of freedom. The microsystem evolution is then described by tracing over the reservoir states, which results in an irreversible decay as excitation leaks out of the initially excited microsystems into the outer reservoir environment. Earlier treatments of this dissipation used density matrices to describe an ensemble of microsystems, either in the Schr\"odinger picture with master equations, or in the Heisenberg picture with Langevin equations. The development of experimental techniques to study single quantum systems (for example, single trapped ions, or cavity-radiation-field modes) has stimulated the construction of theoretical methods to describe individual realizations conditioned on a particular observation record of the decay channel. These methods, variously described as quantum-jump, Monte Carlo wave function, and quantum-trajectory methods, are the subject of this review article. We discuss their derivation, apply them to a number of current problems in quantum optics, and relate them to ensemble descriptions.