Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Annual Meeting Optical Society of America

Monte Carlo simulation of the quantum master equation for vacuum, thermal, and squeezed reservoirs

P. Zoller

Published Jan 1, 1992• 0 citations
Publisher

Abstract

We discuss Monte Carlo simulation of the quantum master equation in terms of system wave functions undergoing a sequence of quantum jumps based on the delay function for the occurrence of the “next jump”.1-3 The formalism can be applied to calculate the density operator and correlation functions (spectra) in quantum systems with a large number of degrees of freedom. We present the formalism and applications to (i) a multilevel atom undergoing spontaneous emission including mechanical light effects (in particular, laser cooling in the quantum limit), (ii) coupling of an atomic system to a squeezed vacuum, and (iii) cavity QED problems with coupling to thermal and squeezed baths.

DOI
10.1364/oam.1992.mff2

Related papers

Monte Carlo simulation of master equations in quantum optics for vacuum, thermal, and squeezed reservoirs

Ralph Dum, A. S. Parkins, P. Zoller +1 more

Wave-function simulation of the master equation in terms of quantum jumps is illustrated for vacuum, thermal, and squeezed reservoirs. We discuss simulation techniques for (i) atomic density matrices, and resonance fluoresence and weak-field absorption spectra of atoms, (ii) decay of a two-level system in a squeezed vacuum, and (iii) a strongly coupled atom-cavity system driven by thermal light.

Physical Review A•Oct 1, 1992•245 citations

Monte Carlo wave-function method in quantum optics

Klaus Mølmer, Yvan Castin, Jean Dalibard

We present a wave-function approach to the study of the evolution of a small system when it is coupled to a large reservoir. Fluctuations and dissipation originate in this approach from quantum jumps that occur randomly during the time evolution of the system. This approach can be applied to a wide class of relaxation operators in the Markovian regime, and it is equivalent to the standard master-equation approach. For systems with a number of states N much larger than unity this Monte Carlo wave-function approach can be less expensive in terms of calculation time than the master-equation treatment. Indeed, a wave function involves only N components, whereas a density matrix is described by N2 terms. We evaluate the gain in computing time that may be expected from such a formalism, and we discuss its applicability to several examples, with particular emphasis on a quantum description of laser cooling.

Journal of the Optical Society of America B•Mar 1, 1993•1,238 citations

Monte Carlo simulation of the atomic master equation for spontaneous emission

Ralph Dum, P. Zoller, Helmut Ritsch

A Monte Carlo simulation of the atomic master equation for spontaneous emission in terms of atomic wave functions is developed. Realizations of the time evolution of atomic wave functions are constructed that correspond to an ensemble of atoms driven by laser light undergoing a sequence of spontaneous emissions. The atomic decay times are drawn according to the photon count distribution of the driven atom. Each quantum jump of the atomic electron projects the atomic wave function to the ground state of the atom. Our theory is based on a stochastic interpretation and generalization of Mollow's pure-state analysis of resonant light scattering, and the Srinivas-Davies theory of continuous measurements in photodetection. An extension of the theory to include mechanical light effects and a generalization to atomic systems with Zeeman substructure are given. We illustrate the method by simulating the solutions of the optical Bloch equations for two-level systems, and laser cooling of a two-level atom in an ion trap where the center-of-mass motion of the atom is described quantum mechanically.

Physical Review A•Apr 1, 1992•656 citations

A Monte Carlo wave function method in quantum optics

Klaus Mølmer, Kirstine Berg‐Sørensen, Yvan Castin +1 more

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.

Annual Meeting Optical Society of America•Jan 1, 1992•44 citations

Stochastic Schrödinger equations in cavity QED: physical interpretation and localization

Tarso B. Ledur Kist, Miguel Orszag, Todd A. Brun +1 more

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.

arXiv:quant-ph/9805027•Journal of Optics B Quantum and Semiclassical Optics•Jan 1, 1999•37 citations

The quantum-jump approach to dissipative dynamics in quantum optics

Martin Bodo Plenio, P. L. Knight

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.

arXiv:quant-ph/9702007•Reviews of Modern Physics•Jan 1, 1998•1,514 citations
We discuss Monte Carlo simulation of the quantum master equation in terms of system wave functions undergoing a sequence of quantum jumps based on the delay function for the occurrence of the “next jump”.1-3 The formalism can be applied to calculate the density operator and correlation functions (spectra) in quantum systems with a large number of degrees of freedom. We present the formalism and applications to (i) a multilevel atom undergoing spontaneous emission including mechanical light effects (in particular, laser cooling in the quantum limit), (ii) coupling of an atomic system to a squeezed vacuum, and (iii) cavity QED problems with coupling to thermal and squeezed baths.
We present a wave-function approach to the study of the evolution of a small system when it is coupled to a large reservoir. Fluctuations and dissipation originate in this approach from quantum jumps that occur randomly during the time evolution of the system. This approach can be applied to a wide class of relaxation operators in the Markovian regime, and it is equivalent to the standard master-equation approach. For systems with a number of states N much larger than unity this Monte Carlo wave-function approach can be less expensive in terms of calculation time than the master-equation treatment. Indeed, a wave function involves only N components, whereas a density matrix is described by N2 terms. We evaluate the gain in computing time that may be expected from such a formalism, and we discuss its applicability to several examples, with particular emphasis on a quantum description of laser cooling.
A Monte Carlo simulation of the atomic master equation for spontaneous emission in terms of atomic wave functions is developed. Realizations of the time evolution of atomic wave functions are constructed that correspond to an ensemble of atoms driven by laser light undergoing a sequence of spontaneous emissions. The atomic decay times are drawn according to the photon count distribution of the driven atom. Each quantum jump of the atomic electron projects the atomic wave function to the ground state of the atom. Our theory is based on a stochastic interpretation and generalization of Mollow's pure-state analysis of resonant light scattering, and the Srinivas-Davies theory of continuous measurements in photodetection. An extension of the theory to include mechanical light effects and a generalization to atomic systems with Zeeman substructure are given. We illustrate the method by simulating the solutions of the optical Bloch equations for two-level systems, and laser cooling of a two-level atom in an ion trap where the center-of-mass motion of the atom is described quantum mechanically.
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 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.
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.