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

arXiv:1210.1302 · Journal of Mathematical Physics

Functional Wigner representation of quantum dynamics of Bose-Einstein condensate

Bogdan Opanchuk, P. D. Drummond

Published Apr 1, 2013• 48 citations
arXivHTMLPDF

Abstract

We develop a method of simulating the full quantum field dynamics of multi-mode multi-component Bose-Einstein condensates in a trap. We use the truncated Wigner representation to obtain a probabilistic theory that can be sampled. This method produces c-number stochastic equations which may be solved using conventional stochastic methods. The technique is valid for large mode occupation numbers. We give a detailed derivation of methods of functional Wigner representation appropriate for quantum fields. Our approach describes spatial evolution of spinor components and properly accounts for nonlinear losses. Such techniques are applicable to calculating the leading quantum corrections, including effects such as quantum squeezing, entanglement, EPR correlations, and interactions with engineered nonlinear reservoirs. By using a consistent expansion in the inverse density, we are able to explain an inconsistency in the nonlinear loss equations found by earlier authors.

DOI
10.1063/1.4801781

Related papers

Quasiprobability methods in quantum interferometry of ultracold matter

Bogdan Opanchuk

A method of simulating the full quantum field dynamics of multimode, multi-component Bose-Einstein condensates is developed. The truncated Wigner representation is used to obtain a probabilistic theory. The representation is extended to to use a functional calculus and can be applied to quantum fields in the arbitrary number of dimensions. Detailed proofs of the corresponding theorems are given. This approach produces c-number stochastic equations which may be solved using conventional stochastic integration methods, with the limitation that it is only valid for large mode occupation numbers. The equations describe the spatial evolution of spinor components and properly account for nonlinear interactions and losses. The method is further used to model several experiments where quantum effects are significant, and cannot be calculated with simpler approximations. It is shown to describe leading quantum corrections accurately, including effects such as phase noise, quantum squeezing, entanglement, and interactions with engineered nonlinear reservoirs. In addition to this approach, other positive phase-space methods are also analysed for sampling highly nonclassical states. In conclusion, the phase-space representations developed in this thesis provide a relatively fast and accurate method of simulating dynamics of large nonlinear quantum systems, such as trapped Bose-Einstein condensates. The required stochastic equations are easily and scalably integrated in parallel. Algorithms for efficient computation on modern graphical processing unit hardware are presented.

•Jan 1, 2014•2 citations

Quantum noise in three-dimensional BEC interferometry

Bogdan Opanchuk, M. Egorov, Scott E. Hoffmann +2 more

We develop a theory of quantum fluctuations and squeezing in a three-dimensional Bose-Einstein condensate atom interferometer with nonlinear losses. We use stochastic equations in a truncated Wigner representation to treat quantum noise. Our approach includes the multi-mode spatial evolution of spinor components and describes the many-body dynamics of a mesoscopic quantum system.

arXiv:1105.5493•Europhysics Letters (EPL)•Mar 1, 2012•32 citations

Hybrid phase-space simulation method for interacting Bose fields

Scott E. Hoffmann, Joel Frederick Corney, P. D. Drummond

We introduce an approximate phase-space technique to simulate the quantum dynamics of interacting bosons. With the future goal of treating Bose-Einstein condensate systems, the method is designed for systems with a natural separation into highly occupied (condensed) modes and lightly occupied modes. The method self-consistently uses the Wigner representation to treat highly occupied modes and the positive-$P$ representation for lightly occupied modes. In this method, truncation of higher-derivative terms from the Fokker-Planck equation is usually necessary. However, at least in the cases investigated here, the resulting systematic error, over a finite time, vanishes in the limit of large Wigner occupation numbers. We tested the method on a system of two interacting anharmonic oscillators, with high and low occupations, respectively. The hybrid method successfully predicted atomic quadratures to a useful simulation time 60 times longer than that of the positive-$P$ method. The truncated Wigner method also performed well in this test. For the prediction of the correlation in a quantum nondemolition measurement scheme, for this same system, the hybrid method gave excellent agreement with the exact result, while the truncated Wigner method showed a large systematic error.

arXiv:0803.1887•Physical Review A•Jul 18, 2008•22 citations

The truncated Wigner method for Bose-condensed gases: limits of validity and applications

Alice Sinatra, Carlos Lobo, Yvan Castin

We study the truncated Wigner method applied to a weakly interacting spinless Bose-condensed gas which is perturbed away from thermal equilibrium by a time-dependent external potential. The principle of the method is to generate an ensemble of classical fields ψ( r ) which samples the Wigner quasi-distribution function of the initial thermal equilibrium density operator of the gas, and then to evolve each classical field with the Gross–Pitaevskii equation. In the first part of the paper we improve the sampling technique over our previous work (Sinatra et al 2000 J. Mod. Opt. 47 2629–44) and we test its accuracy against the exactly solvable model of the ideal Bose gas. In the second part of the paper we investigate the conditions of validity of the truncated Wigner method. For short evolution times it is known that the time-dependent Bogoliubov approximation is valid for almost pure condensates. The requirement that the truncated Wigner method reproduces the Bogoliubov prediction leads to the constraint that the number of field modes in the Wigner simulation must be smaller than the number of particles in the gas. For longer evolution times the nonlinear dynamics of the noncondensed modes of the field plays an important role. To demonstrate this we analyse the case of a three-dimensional spatially homogeneous Bose-condensed gas and we test the ability of the truncated Wigner method to correctly reproduce the Beliaev–Landau damping of an excitation of the condensate. We have identified the mechanism which limits the validity of the truncated Wigner method: the initial ensemble of classical fields, driven by the time-dependent Gross–Pitaevskii equation, thermalizes to a classical field distribution at a temperature T class which is larger than the initial temperature T of the quantum gas. When T class significantly exceeds T a spurious damping is observed in the Wigner simulation. This leads to the second validity condition for the truncated Wigner method, T class − T ≪ T , which requires that the maximum energy max of the Bogoliubov modes in the simulation does not exceed a few k B T .

arXiv:cond-mat/0201217•Journal of Physics B Atomic Molecular and Optical Physics•Aug 21, 2002•352 citations

Dynamical quantum noise in trapped Bose-Einstein condensates

Michael J. Steel, M. K. Olsen, Lev I. Plimak +5 more

We introduce the study of dynamical quantum noise in Bose-Einstein condensates through numerical simulation of stochastic partial differential equations obtained using phase-space representations. We derive evolution equations for a single trapped condensate in both the positive-$P$ and Wigner representations and perform simulations to compare the predictions of the two methods. The positive-$P$ approach is found to be highly susceptible to the stability problems that have been observed in other strongly nonlinear, weakly damped systems. Using the Wigner representation, we examine the evolution of several quantities of interest using from a variety of choices of initial state for the condensate and compare results to those for single-mode models.

arXiv:cond-mat/9807349•Physical Review A•Dec 1, 1998•433 citations

Quantum dynamics in ultracold atomic physics

Qiongyi He, M. D. Reid, Bogdan Opanchuk +3 more

We review recent developments in the theory of quantum dynamics in ultracold atomic physics, including exact techniques and methods based on phase-space mappings that are applicable when the complexity becomes exponentially large. Phase-space representations include the truncated Wigner, positive- P and general Gaussian operator representations which can treat both bosons and fermions. These phase-space methods include both traditional approaches using a phase-space of classical dimension, and more recent methods that use a non-classical phase-space of increased dimensionality. Examples used include quantum Einstein-Podolsky-Rosen (EPR) entanglement of a four-mode BEC, time-reversal tests of dephasing in single-mode traps, BEC quantum collisions with up to 106 modes and 105 interacting particles, quantum interferometry in a multi-mode trap with nonlinear absorption, and the theory of quantum entropy in phase-space. We also treat the approach of variational optimization of the sampling error, giving an elementary example of a nonlinear oscillator.

arXiv:1112.0380•Frontiers of Physics•Jan 21, 2012•16 citations
We develop a method of simulating the full quantum field dynamics of multi-mode multi-component Bose-Einstein condensates in a trap. We use the truncated Wigner representation to obtain a probabilistic theory that can be sampled. This method produces c-number stochastic equations which may be solved using conventional stochastic methods. The technique is valid for large mode occupation numbers. We give a detailed derivation of methods of functional Wigner representation appropriate for quantum fields. Our approach describes spatial evolution of spinor components and properly accounts for nonlinear losses. Such techniques are applicable to calculating the leading quantum corrections, including effects such as quantum squeezing, entanglement, EPR correlations, and interactions with engineered nonlinear reservoirs. By using a consistent expansion in the inverse density, we are able to explain an inconsistency in the nonlinear loss equations found by earlier authors.
A method of simulating the full quantum field dynamics of multimode, multi-component Bose-Einstein condensates is developed. The truncated Wigner representation is used to obtain a probabilistic theory. The representation is extended to to use a functional calculus and can be applied to quantum fields in the arbitrary number of dimensions. Detailed proofs of the corresponding theorems are given. This approach produces c-number stochastic equations which may be solved using conventional stochastic integration methods, with the limitation that it is only valid for large mode occupation numbers. The equations describe the spatial evolution of spinor components and properly account for nonlinear interactions and losses. The method is further used to model several experiments where quantum effects are significant, and cannot be calculated with simpler approximations. It is shown to describe leading quantum corrections accurately, including effects such as phase noise, quantum squeezing, entanglement, and interactions with engineered nonlinear reservoirs. In addition to this approach, other positive phase-space methods are also analysed for sampling highly nonclassical states. In conclusion, the phase-space representations developed in this thesis provide a relatively fast and accurate method of simulating dynamics of large nonlinear quantum systems, such as trapped Bose-Einstein condensates. The required stochastic equations are easily and scalably integrated in parallel. Algorithms for efficient computation on modern graphical processing unit hardware are presented.
We introduce an approximate phase-space technique to simulate the quantum dynamics of interacting bosons. With the future goal of treating Bose-Einstein condensate systems, the method is designed for systems with a natural separation into highly occupied (condensed) modes and lightly occupied modes. The method self-consistently uses the Wigner representation to treat highly occupied modes and the positive-PPP representation for lightly occupied modes. In this method, truncation of higher-derivative terms from the Fokker-Planck equation is usually necessary. However, at least in the cases investigated here, the resulting systematic error, over a finite time, vanishes in the limit of large Wigner occupation numbers. We tested the method on a system of two interacting anharmonic oscillators, with high and low occupations, respectively. The hybrid method successfully predicted atomic quadratures to a useful simulation time 60 times longer than that of the positive-PPP method. The truncated Wigner method also performed well in this test. For the prediction of the correlation in a quantum nondemolition measurement scheme, for this same system, the hybrid method gave excellent agreement with the exact result, while the truncated Wigner method showed a large systematic error.
We study the truncated Wigner method applied to a weakly interacting spinless Bose-condensed gas which is perturbed away from thermal equilibrium by a time-dependent external potential. The principle of the method is to generate an ensemble of classical fields ψ( r ) which samples the Wigner quasi-distribution function of the initial thermal equilibrium density operator of the gas, and then to evolve each classical field with the Gross–Pitaevskii equation. In the first part of the paper we improve the sampling technique over our previous work (Sinatra et al 2000 J. Mod. Opt. 47 2629–44) and we test its accuracy against the exactly solvable model of the ideal Bose gas. In the second part of the paper we investigate the conditions of validity of the truncated Wigner method. For short evolution times it is known that the time-dependent Bogoliubov approximation is valid for almost pure condensates. The requirement that the truncated Wigner method reproduces the Bogoliubov prediction leads to the constraint that the number of field modes in the Wigner simulation must be smaller than the number of particles in the gas. For longer evolution times the nonlinear dynamics of the noncondensed modes of the field plays an important role. To demonstrate this we analyse the case of a three-dimensional spatially homogeneous Bose-condensed gas and we test the ability of the truncated Wigner method to correctly reproduce the Beliaev–Landau damping of an excitation of the condensate. We have identified the mechanism which limits the validity of the truncated Wigner method: the initial ensemble of classical fields, driven by the time-dependent Gross–Pitaevskii equation, thermalizes to a classical field distribution at a temperature T class which is larger than the initial temperature T of the quantum gas. When T class significantly exceeds T a spurious damping is observed in the Wigner simulation. This leads to the second validity condition for the truncated Wigner method, T class − T ≪ T , which requires that the maximum energy max of the Bogoliubov modes in the simulation does not exceed a few k B T .
We introduce the study of dynamical quantum noise in Bose-Einstein condensates through numerical simulation of stochastic partial differential equations obtained using phase-space representations. We derive evolution equations for a single trapped condensate in both the positive-PPP and Wigner representations and perform simulations to compare the predictions of the two methods. The positive-PPP approach is found to be highly susceptible to the stability problems that have been observed in other strongly nonlinear, weakly damped systems. Using the Wigner representation, we examine the evolution of several quantities of interest using from a variety of choices of initial state for the condensate and compare results to those for single-mode models.
We review recent developments in the theory of quantum dynamics in ultracold atomic physics, including exact techniques and methods based on phase-space mappings that are applicable when the complexity becomes exponentially large. Phase-space representations include the truncated Wigner, positive- P and general Gaussian operator representations which can treat both bosons and fermions. These phase-space methods include both traditional approaches using a phase-space of classical dimension, and more recent methods that use a non-classical phase-space of increased dimensionality. Examples used include quantum Einstein-Podolsky-Rosen (EPR) entanglement of a four-mode BEC, time-reversal tests of dephasing in single-mode traps, BEC quantum collisions with up to 106 modes and 105 interacting particles, quantum interferometry in a multi-mode trap with nonlinear absorption, and the theory of quantum entropy in phase-space. We also treat the approach of variational optimization of the sampling error, giving an elementary example of a nonlinear oscillator.