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:0803.1887 · Physical Review A

Hybrid phase-space simulation method for interacting Bose fields

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

Published Jul 18, 2008• 22 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.78.013622

Related papers

Functional Wigner representation of quantum dynamics of Bose-Einstein condensate

Bogdan Opanchuk, P. D. Drummond

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.

arXiv:1210.1302•Journal of Mathematical Physics•Apr 1, 2013•48 citations

Comparative study of dynamical simulation methods for the dissociation of molecular Bose-Einstein condensates

S. L. W. Midgley, Sebastian Wüster, M. K. Olsen +2 more

We describe a pairing mean-field theory related to the Hartree-Fock-Bogoliubov approach, and apply it to the dynamics of dissociation of a molecular Bose-Einstein condensate into correlated bosonic atom pairs. We also perform the same simulation using two stochastic phase-space techniques for quantum dynamics---the positive-$P$ representation method and the truncated Wigner method. By comparing the results of our calculations we are able to assess the relative strengths of these theoretical techniques in describing molecular dissociation in one spatial dimension. An important aspect of our analysis is the inclusion of atom-atom interactions, which can be problematic for the positive-$P$ method. We find that the truncated Wigner method mostly agrees with the positive-$P$ simulations, but can be simulated for significantly longer times. The pairing mean-field theory results diverge from the quantum dynamical methods after relatively short times.

arXiv:0811.2030•Physical Review A•May 20, 2009•13 citations

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

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

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

Non-Gaussian pure states and positive Wigner functions

Joel Frederick Corney, M. K. Olsen

Non-Gaussian correlations in a pure state are inextricably linked with certain nonclassical features, such as a non-positive-definite Wigner function. In a commonly used simulation technique in ultracold atoms and quantum optics, known as the truncated Wigner method, the quantum dynamics is mapped to stochastic trajectories in phase space, governed by a positive approximation to the true Wigner distribution. The question thus arises: How accurate is this approach in predicting truly nonclassical behavior? In this article, we benchmark the ability of the truncated Wigner phase-space method to reproduce the non-Gaussian statistics of the single-mode anharmonic oscillator. We find that the this method can reliably predict departures from Gaussian statistics over a wide range of particle numbers, whereas the positive-$P$ representation, which involves no approximations, is limited by rapidly growing statistical uncertainty. The truncated Wigner function, furthermore, is able to reproduce the non-Gaussian correlations while satisfying the condition for purity.

arXiv:1412.4868•Physical Review A•Feb 19, 2015•25 citations
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 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.
We describe a pairing mean-field theory related to the Hartree-Fock-Bogoliubov approach, and apply it to the dynamics of dissociation of a molecular Bose-Einstein condensate into correlated bosonic atom pairs. We also perform the same simulation using two stochastic phase-space techniques for quantum dynamics---the positive-PPP representation method and the truncated Wigner method. By comparing the results of our calculations we are able to assess the relative strengths of these theoretical techniques in describing molecular dissociation in one spatial dimension. An important aspect of our analysis is the inclusion of atom-atom interactions, which can be problematic for the positive-PPP method. We find that the truncated Wigner method mostly agrees with the positive-PPP simulations, but can be simulated for significantly longer times. The pairing mean-field theory results diverge from the quantum dynamical methods after relatively short times.
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 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 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.
Non-Gaussian correlations in a pure state are inextricably linked with certain nonclassical features, such as a non-positive-definite Wigner function. In a commonly used simulation technique in ultracold atoms and quantum optics, known as the truncated Wigner method, the quantum dynamics is mapped to stochastic trajectories in phase space, governed by a positive approximation to the true Wigner distribution. The question thus arises: How accurate is this approach in predicting truly nonclassical behavior? In this article, we benchmark the ability of the truncated Wigner phase-space method to reproduce the non-Gaussian statistics of the single-mode anharmonic oscillator. We find that the this method can reliably predict departures from Gaussian statistics over a wide range of particle numbers, whereas the positive-PPP representation, which involves no approximations, is limited by rapidly growing statistical uncertainty. The truncated Wigner function, furthermore, is able to reproduce the non-Gaussian correlations while satisfying the condition for purity.