We discuss stochastic phase-space methods within the truncated Wigner approximation and show explicitly that they can be used to solve non-equilibrium dynamics of bosonic atoms in one-dimensional traps. We consider systems both with and without an optical lattice, and address different approximations in the stochastic synthesization of quantum statistical correlations of the initial atomic field. We also present a numerically efficient projection method for analyzing correlation functions of the simulation results, and demonstrate physical examples of non-equilibrium quantum dynamics of solitons and atom number squeezing in optical lattices.We discuss stochastic phase-space methods within the truncated Wigner approximation and show explicitly that they can be used to solve non-equilibrium dynamics of bosonic atoms in one-dimensional traps. We consider systems both with and without an optical lattice, and address different approximations in the stochastic synthesization of quantum statistical correlations of the initial atomic field. We also present a numerically efficient projection method for analyzing correlation functions of the simulation results, and demonstrate physical examples of non-equilibrium quantum dynamics of solitons and atom number squeezing in optical lattices.

Cold atomsMar 18, 201315 citations

Variational truncated Wigner approximation for weakly interacting Bose fields: Dynamics of coupled condensates

Christopher D. Mink, Axel Pelster, Jens Benary +2 more

The truncated Wigner approximation is an established approach that describes the dynamics of weakly interacting Bose gases beyond the mean-field level. Although it allows a quantum field to be expressed by a stochastic c-number field, the simulation of the time evolution is still very demanding for most applications. Here, we develop a numerically inexpensive scheme by approximating the c-number field with a variational ansatz. The dynamics of the ansatz function is described by a tractable set of coupled ordinary stochastic differential equations for the respective variational parameters. We investigate the non-equilibrium dynamics of a three-dimensional Bose gas in a one-dimensional optical lattice with a transverse isotropic harmonic confinement. The accuracy and computational inexpensiveness of our method are demonstrated by comparing its predictions to experimental data.The truncated Wigner approximation is an established approach that describes the dynamics of weakly interacting Bose gases beyond the mean-field level. Although it allows a quantum field to be expressed by a stochastic c-number field, the simulation of the time evolution is still very demanding for most applications. Here, we develop a numerically inexpensive scheme by approximating the c-number field with a variational ansatz. The dynamics of the ansatz function is described by a tractable set of coupled ordinary stochastic differential equations for the respective variational parameters. We investigate the non-equilibrium dynamics of a three-dimensional Bose gas in a one-dimensional optical lattice with a transverse isotropic harmonic confinement. The accuracy and computational inexpensiveness of our method are demonstrated by comparing its predictions to experimental data.

arXiv:2106.05354SciPost PhysicsFeb 2, 202213 citations

Quantum dynamics in splitting a harmonically trapped Bose-Einstein condensate by an optical lattice: Truncated Wigner approximation

Lorenzo Isella, Janne Ruostekoski

We study the splitting of a harmonically trapped atomic Bose-Einstein condensate when we continuously turn up an optical lattice (or a double-well) potential. As the lattice height is increased, quantum fluctuations of atoms are enhanced. The resulting nonequilibrium dynamics of the fragmentation process of the condensate, the loss of the phase coherence of atoms along the lattice, and the reduced atom number fluctuations in individual lattice sites are stochastically studied within the truncated Wigner approximation. We perform a detailed study of the effects of temperature and lattice height on atom dynamics, and investigate the validity of the classical Gross-Pitaevskii equation in optical lattices. We find the atom number squeezing to saturate in deep lattices due to nonadiabaticity in turning up of the lattice potential that is challenging to avoid in experiments when the occupation number of the lattice sites is large, making it difficult to produce strongly number squeezed (or the Mott insulator) states with large filling factors. We also investigate some general numerical properties of the truncated Wigner approximation.We study the splitting of a harmonically trapped atomic Bose-Einstein condensate when we continuously turn up an optical lattice (or a double-well) potential. As the lattice height is increased, quantum fluctuations of atoms are enhanced. The resulting nonequilibrium dynamics of the fragmentation process of the condensate, the loss of the phase coherence of atoms along the lattice, and the reduced atom number fluctuations in individual lattice sites are stochastically studied within the truncated Wigner approximation. We perform a detailed study of the effects of temperature and lattice height on atom dynamics, and investigate the validity of the classical Gross-Pitaevskii equation in optical lattices. We find the atom number squeezing to saturate in deep lattices due to nonadiabaticity in turning up of the lattice potential that is challenging to avoid in experiments when the occupation number of the lattice sites is large, making it difficult to produce strongly number squeezed (or the Mott insulator) states with large filling factors. We also investigate some general numerical properties of the truncated Wigner approximation.

arXiv:cond-mat/0608634Physical Review ADec 22, 200651 citations

SU(3) truncated Wigner approximation for strongly interacting Bose gases

Kazuma Nagao, Yosuke Takasu, Yoshiro Takahashi +1 more

We develop and utilize the SU(3) truncated Wigner approximation (TWA) in order to analyze far-from-equilibrium quantum dynamics of strongly interacting Bose gases in an optical lattice. Specifically, we explicitly represent the corresponding Bose-Hubbard model at an arbitrary filling factor with restricted local Hilbert spaces in terms of SU(3) matrices. Moreover, we introduce a discrete Wigner sampling technique for the SU(3) TWA and examine its performance as well as that of the SU(3) TWA with the Gaussian approximation for the continuous Wigner function. We directly compare outputs of these two approaches with exact computations regarding dynamics of the Bose-Hubbard model at unit filling with a small size and that of a fully connected spin-1 model with a large size. We show that both approaches can quantitatively capture quantum dynamics on a timescale of $\ensuremath{\hbar}/(Jz)$, where $J$ and $z$ denote the hopping energy and the coordination number. We apply the two kinds of SU(3) TWA to dynamical spreading of a two-point correlation function of the Bose-Hubbard model on a square lattice with a large system size, which has been measured in recent experiments. Noticeable deviations between the theories and experiments indicate that proper inclusion of effects of the spatial inhomogeneity, which is not straightforward in our formulation of the SU(3) TWA, may be necessary.We develop and utilize the SU(3) truncated Wigner approximation (TWA) in order to analyze far-from-equilibrium quantum dynamics of strongly interacting Bose gases in an optical lattice. Specifically, we explicitly represent the corresponding Bose-Hubbard model at an arbitrary filling factor with restricted local Hilbert spaces in terms of SU(3) matrices. Moreover, we introduce a discrete Wigner sampling technique for the SU(3) TWA and examine its performance as well as that of the SU(3) TWA with the Gaussian approximation for the continuous Wigner function. We directly compare outputs of these two approaches with exact computations regarding dynamics of the Bose-Hubbard model at unit filling with a small size and that of a fully connected spin-1 model with a large size. We show that both approaches can quantitatively capture quantum dynamics on a timescale of $\ensuremath{\hbar}/(Jz)$, where $J$ and $z$ denote the hopping energy and the coordination number. We apply the two kinds of SU(3) TWA to dynamical spreading of a two-point correlation function of the Bose-Hubbard model on a square lattice with a large system size, which has been measured in recent experiments. Noticeable deviations between the theories and experiments indicate that proper inclusion of effects of the spatial inhomogeneity, which is not straightforward in our formulation of the SU(3) TWA, may be necessary.

arXiv:2008.09900Physical Review ResearchNov 1, 20219 citations

C-Field Methods for Non-Equilibrium Bose Gases

Matthew J. Davis, Tod M. Wright, P. B. Blakie +3 more

We review c-field methods for simulating the non-equilibrium dynamics of degenerate Bose gases beyond the mean-field Gross-Pitaevskii approximation. We describe three separate approaches that utilise similar numerical methods, but have distinct regimes of validity. Systems at finite temperature can be treated with either the closed-system projected Gross-Pitaevskii equation (PGPE), or the open-system stochastic projected Gross-Pitaevskii equation (SPGPE). These are both applicable in quantum degenerate regimes in which thermal fluctuations are significant. At low or zero temperature, the truncated Wigner projected Gross-Pitaevskii equation (TWPGPE) allows for the simulation of systems in which spontaneous collision processes seeded by quantum fluctuations are important. We describe the regimes of validity of each of these methods, and discuss their relationships to one another, and to other simulation techniques for the dynamics of Bose gases. The utility of the SPGPE formalism in modelling non-equilibrium Bose gases is illustrated by its application to the dynamics of spontaneous vortex formation in the growth of a Bose-Einstein condensate.We review c-field methods for simulating the non-equilibrium dynamics of degenerate Bose gases beyond the mean-field Gross-Pitaevskii approximation. We describe three separate approaches that utilise similar numerical methods, but have distinct regimes of validity. Systems at finite temperature can be treated with either the closed-system projected Gross-Pitaevskii equation (PGPE), or the open-system stochastic projected Gross-Pitaevskii equation (SPGPE). These are both applicable in quantum degenerate regimes in which thermal fluctuations are significant. At low or zero temperature, the truncated Wigner projected Gross-Pitaevskii equation (TWPGPE) allows for the simulation of systems in which spontaneous collision processes seeded by quantum fluctuations are important. We describe the regimes of validity of each of these methods, and discuss their relationships to one another, and to other simulation techniques for the dynamics of Bose gases. The utility of the SPGPE formalism in modelling non-equilibrium Bose gases is illustrated by its application to the dynamics of spontaneous vortex formation in the growth of a Bose-Einstein condensate.

arXiv:1206.5470Cold atomsMar 18, 20132 citations

Quantum and Thermal Effects of Dark Solitons in a One-Dimensional Bose Gas

Andrew D. Martin, Janne Ruostekoski

We numerically study the imprinting and dynamics of dark solitons in a bosonic atomic gas in a tightly confined one-dimensional harmonic trap both with and without an optical lattice. Quantum and thermal fluctuations are synthesized within the truncated Wigner approximation in the quasicondensate description. We track the soliton coordinates and calculate position and velocity uncertainties. We find that the phase fluctuations lower the classically predicted soliton speed and seed instabilities. Individual runs show interactions of solitons with sound waves, splitting, and disappearing solitons.We numerically study the imprinting and dynamics of dark solitons in a bosonic atomic gas in a tightly confined one-dimensional harmonic trap both with and without an optical lattice. Quantum and thermal fluctuations are synthesized within the truncated Wigner approximation in the quasicondensate description. We track the soliton coordinates and calculate position and velocity uncertainties. We find that the phase fluctuations lower the classically predicted soliton speed and seed instabilities. Individual runs show interactions of solitons with sound waves, splitting, and disappearing solitons.

arXiv:0909.2621Physical Review LettersMay 12, 201072 citations

J

and

z

denote the hopping energy and the coordination number. We apply the two kinds of SU(3) TWA to dynamical spreading of a two-point correlation function of the Bose-Hubbard model on a square lattice with a large system size, which has been measured in recent experiments. Noticeable deviations between the theories and experiments indicate that proper inclusion of effects of the spatial inhomogeneity, which is not straightforward in our formulation of the SU(3) TWA, may be necessary.