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

Physical Review A

Density-matrix equations and photon recoil for multistate atoms

Juha Javanainen

Published Nov 1, 1991• 59 citations
Publisher

Abstract

We derive a Fokker-Planck equation (FPE) for the center-of-mass Wigner function of an atom with an arbitrary level structure, moving in a given three-dimensional light field. This FPE is obtained by eliminating adiabatically the internal degrees of freedom from the expansion of the equation of motion of the atomic density operator into a power series of photon momentum. The force and the diffusion tensor in the FPE may be computed by solving a set of differential equations closely related to the time-evolution equations of the density matrix for the internal state of the atom, as encountered in the conventional recoilless spectroscopy. By studying spontaneous, induced, and free-evolution terms in the equations of motion of the density operator for various level configurations, we devise algorithms that enable a computer to build all the required equations automatically for an arbitrary atomic level scheme.

DOI
10.1103/physreva.44.5857

Related papers

Laser cooling theory of multistate atoms

Juha Javanainen

In the limit of asymptotically small photon momentum, a Fokker–Planck equation for the center-of-mass Wigner function suffices to describe laser cooling of an atom. We have devised a procedure to compute the light pressure force and the diffusion tensor for an atom with an arbitrary level structure, traveling with an arbitrary velocity in any given three-dimensional light field. The scheme is based on solving a set of differential equations closely related to the density matrix equations of motion encountered in conventional recoilless laser spectroscopy. Computer programs have been written that automatically construct and solve the required equations from broad specifications of the problem, such as energies and angular momenta of the atomic levels. A Langevin equation stochastic simulation has been implemented to solve the Fokker-Planck equation of laser cooling using the computed force and diffusion. Comparisons of optical-molasses experiments and simulations carried out for a model sodium atom with a J = 2→J = 3 transition show qualitative similarities such as sub-Doppler limit temperatures, but the quantitative agreement of theory and experiment is not satisfactory.

Annual Meeting Optical Society of America•Jan 1, 1990•0 citations

Atom-field interactions: Density-matrix equations including quantization of the center-of-mass motion

B. Dubetsky, P. R. Berman

Master equations are derived for the time evolution of density-matrix elements characterizing atoms interacting with classical radiation fields that drive transitions between electronic-state manifolds of levels. The atomic center-of-mass motion is quantized and recoil accompanying absorption, stimulated emission, and spontaneous emission is included. An adiabatic elimination of the optical coherences is used to obtain quantum rate equations for ground- and excited-state density-matrix elements. Subsequent adiabatic elimination of excited-state density-matrix elements results in rate equations for the slow evolution of ground-state density-matrix elements. It is shown that recoil during spontaneous emission mixes spatial and internal atomic degrees of freedom. \textcopyright{} 1996 The American Physical Society.

Physical Review A•Jan 1, 1996•10 citations

Light-induced motion of trapped ions. III. Expansion around the recoilless solution

Juha Javanainen

For pt.II see ibid., vol.14, p.4191 (1981). The author derives a Fokker-Planck equation for the Wigner function describing the centre-of-mass motion of a harmonically trapped ion in a travelling wave of light, using an expansion in the recoil velocity around the exact time-dependent solution of the recoilless case. The practical calculations are reduced to the scheme presented in the earlier papers of this series. However, the final energy of the cooling differs from his earlier results, and agrees with more recent calculations employing for the centre of mass the representation in the harmonic oscillator states.

Journal of Physics B Atomic and Molecular Physics•Apr 28, 1985•18 citations

The N-atom laser below saturation: a density matrix approach without large-N scaling

Michael Fleischhauer, S Sandhager

The authors present a method of deriving a reduced density operator equation for a quantized field interacting with an arbitrary number of atoms that does not neglect atom-atom correlations. As a particular example they analyse a three-level laser operating well below saturation and solve the density matrix equations in steady state. For weak pumping the system is equivalent to the Scully-Lamb model (1966) and the corresponding equations are recovered. Using the saturation photon number as a (large) scaling parameter one can derive a Fokker-Planck equation for the Glauber-P-distribution. The stationary solution of the reduced density matrix equation is compared with the solution of the Fokker-Planck equation and-for the case of a single-atom laser-with a numerical solution of the total density matrix equation.

Quantum and Semiclassical Optics Journal of the European Optical Society Part B•Jun 1, 1995•1 citation

Fokker–Planck Equations

Bryan J. Dalton, John Jeffers, Stephen Mark Barnett

Abstract This chapter deals with deriving Fokker–Planck equations (FPEs) that govern the behaviour of phase space distribution functions (normalised and unnormalised) for boson and fermion systems due to dynamical or thermal evolution. Specific forms are obtained for master equation dynamics, relating the drift vector and diffusion matrix elements in the FPE to coefficients specifying terms in the Hamiltonian and relaxation operators. The approach used is to derive correspondence rules, which map changes in the density operator when multiplied from left or right by mode annihilation or creation operators (as in evolution equations) onto corresponding changes in distribution functions, these being differentiations and or multiplications by phase space variables (left or right for fermions). Derivations are carried out via changes to characteristic functions or to canonical forms of the density operator. The correct successive application of correspondence rules for fermions is treated, and alternative forms obtained for bosons to correspondence rules and the FPE.

•Nov 13, 2014•10 citations

Adiabatic expansion for the single-mode laser

Fritz Haake, Maciej Lewenstein

We carry out a systematic adiabatic elimination of the atomic degrees of freedom from the quantum-mechanical master equation for the single-mode laser. We represent the reduced density operator of the field mode by various quasiprobabilities and construct the respective equations of motion. The generators of infinitesimal time translations are obtained as series in powers of two parameters, the smallness of which defines the adiabatic limit. Our "adiabatic" expansion treats that part of the atom-field interaction as a zeroth-order effect which describes the action of the field on the atoms, while the reaction of the atoms is treated perturbatively. The adiabatic equilibrium thus assigned to the atomic variables at all times is a conditional one, contingent on the current state of the field mode. As a result, saturation effects in the atoms are fully accounted for in low orders of our expansion. In second order, especially, we obtain a Fokker-Planck equation for the Wigner function of the field mode which is valid for arbitrary pump strengths below, near, and above threshold. We compare our results with those of previous theories.

Physical Review A•Feb 1, 1983•35 citations
We derive a Fokker-Planck equation (FPE) for the center-of-mass Wigner function of an atom with an arbitrary level structure, moving in a given three-dimensional light field. This FPE is obtained by eliminating adiabatically the internal degrees of freedom from the expansion of the equation of motion of the atomic density operator into a power series of photon momentum. The force and the diffusion tensor in the FPE may be computed by solving a set of differential equations closely related to the time-evolution equations of the density matrix for the internal state of the atom, as encountered in the conventional recoilless spectroscopy. By studying spontaneous, induced, and free-evolution terms in the equations of motion of the density operator for various level configurations, we devise algorithms that enable a computer to build all the required equations automatically for an arbitrary atomic level scheme.
In the limit of asymptotically small photon momentum, a Fokker–Planck equation for the center-of-mass Wigner function suffices to describe laser cooling of an atom. We have devised a procedure to compute the light pressure force and the diffusion tensor for an atom with an arbitrary level structure, traveling with an arbitrary velocity in any given three-dimensional light field. The scheme is based on solving a set of differential equations closely related to the density matrix equations of motion encountered in conventional recoilless laser spectroscopy. Computer programs have been written that automatically construct and solve the required equations from broad specifications of the problem, such as energies and angular momenta of the atomic levels. A Langevin equation stochastic simulation has been implemented to solve the Fokker-Planck equation of laser cooling using the computed force and diffusion. Comparisons of optical-molasses experiments and simulations carried out for a model sodium atom with a J = 2→J = 3 transition show qualitative similarities such as sub-Doppler limit temperatures, but the quantitative agreement of theory and experiment is not satisfactory.
Master equations are derived for the time evolution of density-matrix elements characterizing atoms interacting with classical radiation fields that drive transitions between electronic-state manifolds of levels. The atomic center-of-mass motion is quantized and recoil accompanying absorption, stimulated emission, and spontaneous emission is included. An adiabatic elimination of the optical coherences is used to obtain quantum rate equations for ground- and excited-state density-matrix elements. Subsequent adiabatic elimination of excited-state density-matrix elements results in rate equations for the slow evolution of ground-state density-matrix elements. It is shown that recoil during spontaneous emission mixes spatial and internal atomic degrees of freedom. \textcopyright{} 1996 The American Physical Society.
For pt.II see ibid., vol.14, p.4191 (1981). The author derives a Fokker-Planck equation for the Wigner function describing the centre-of-mass motion of a harmonically trapped ion in a travelling wave of light, using an expansion in the recoil velocity around the exact time-dependent solution of the recoilless case. The practical calculations are reduced to the scheme presented in the earlier papers of this series. However, the final energy of the cooling differs from his earlier results, and agrees with more recent calculations employing for the centre of mass the representation in the harmonic oscillator states.
The authors present a method of deriving a reduced density operator equation for a quantized field interacting with an arbitrary number of atoms that does not neglect atom-atom correlations. As a particular example they analyse a three-level laser operating well below saturation and solve the density matrix equations in steady state. For weak pumping the system is equivalent to the Scully-Lamb model (1966) and the corresponding equations are recovered. Using the saturation photon number as a (large) scaling parameter one can derive a Fokker-Planck equation for the Glauber-P-distribution. The stationary solution of the reduced density matrix equation is compared with the solution of the Fokker-Planck equation and-for the case of a single-atom laser-with a numerical solution of the total density matrix equation.
Abstract This chapter deals with deriving Fokker–Planck equations (FPEs) that govern the behaviour of phase space distribution functions (normalised and unnormalised) for boson and fermion systems due to dynamical or thermal evolution. Specific forms are obtained for master equation dynamics, relating the drift vector and diffusion matrix elements in the FPE to coefficients specifying terms in the Hamiltonian and relaxation operators. The approach used is to derive correspondence rules, which map changes in the density operator when multiplied from left or right by mode annihilation or creation operators (as in evolution equations) onto corresponding changes in distribution functions, these being differentiations and or multiplications by phase space variables (left or right for fermions). Derivations are carried out via changes to characteristic functions or to canonical forms of the density operator. The correct successive application of correspondence rules for fermions is treated, and alternative forms obtained for bosons to correspondence rules and the FPE.
We carry out a systematic adiabatic elimination of the atomic degrees of freedom from the quantum-mechanical master equation for the single-mode laser. We represent the reduced density operator of the field mode by various quasiprobabilities and construct the respective equations of motion. The generators of infinitesimal time translations are obtained as series in powers of two parameters, the smallness of which defines the adiabatic limit. Our "adiabatic" expansion treats that part of the atom-field interaction as a zeroth-order effect which describes the action of the field on the atoms, while the reaction of the atoms is treated perturbatively. The adiabatic equilibrium thus assigned to the atomic variables at all times is a conditional one, contingent on the current state of the field mode. As a result, saturation effects in the atoms are fully accounted for in low orders of our expansion. In second order, especially, we obtain a Fokker-Planck equation for the Wigner function of the field mode which is valid for arbitrary pump strengths below, near, and above threshold. We compare our results with those of previous theories.