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

Reduced density matrices and coherence of trapped interacting bosons

Kaspar Sakmann, Alexej I. Streltsov, Ofir E. Alon, L. S. Cederbaum

Published Aug 8, 2008• 150 citations
arXivHTMLPDF

Abstract

The first- and second-order correlation functions of trapped interacting Bose-Einstein condensates are investigated numerically on a many-body level from first principles. Correlations in real space and momentum space are treated. The coherence properties are analyzed. The results are obtained by solving the many-body Schr\"odinger equation. It is shown in an example how many-body effects can be induced by the trap geometry. A generic fragmentation scenario of a condensate is considered. The correlation functions are discussed along a pathway from a single condensate to a fragmented condensate. It is shown that strong correlations can arise from the geometry of the trap, even at weak interaction strengths. The natural orbitals and natural geminals of the system are obtained and discussed. It is shown how the fragmentation of the condensate can be understood in terms of its natural geminals. The many-body results are compared to those of mean-field theory. The best solution within mean-field theory is obtained. The limits in which mean-field theories are valid are determined. In these limits the behavior of the correlation functions is explained within an analytical model.

DOI
10.1103/physreva.78.023615

Related papers

Correlation Effects in a Trapped Bose-Fermi Mixture: Exact Results

Ofir E. Alon, L. S. Cederbaum

Abstract Many-body properties of a fermionic impurity embedded in a Bose-Einstein condensate are analyzed analytically using a solvable model, the harmonic-interaction model for Bose-Fermi mixtures. The one-particle and two-particle densities, reduced density matrices, and correlation functions of the fermions and bosons, both in position and momentum spaces, are prescribed in closed form. The various coherence lengths are analyzed. We show that the first-order coherence lengths in position and momentum spaces are equal whereas the second-order quantities can differ substantially. Illustrative examples where the sole interaction is between the impurity and the condensate are presented. Implications are briefly discussed.

arXiv:2311.12192•Journal of Physics Conference Series•Nov 1, 2024•2 citations

Numerically exact dynamics of the interacting many-body Schrödinger equation for Bose-Einstein condensates : comparison to Bose-Hubbard and Gross-Pitaevskii theory

Kaspar Sakmann

In this thesis, the physics of trapped, interacting Bose-Einstein condensates is analyzed by solving the many-body Schrödinger equation. Particular emphasis is put on coherence, fragmentation and reduced density matrices. First, the ground state of a trapped Bose-Einstein condensate and its correlation functions are obtained. Then the dynamics of a bosonic Josephson junction is investigated by solving the time-dependent many-body Schrödinger equation numerically exactly. These are the first exact results in literature in this context. It is shown that the standard approximations of the field, Gross-Pitaevskii theory and the Bose-Hubbard model fail at weak interaction strength and within their range of expected validity. For stronger interactions the dynamics becomes strongly correlated and a new equilibration phenomenon is discovered. By comparison with exact results it is shown that a symmetry of the Bose-Hubbard model between attractive and repulsive interactions must be considered an artefact of the model. A conceptual innovation of this thesis are time-dependent Wannier functions. Equations of motion for time-dependent Wannier functions are derived from the variational principle. By comparison with exact results it is shown that lattice models can be greatly improved at little computational cost by letting the Wannier functions of a lattice model become time-dependent.

•Jan 1, 2010•1 citation

Bose-Einstein condensation in trapped bosons: A variational Monte Carlo analysis

Jonathan L. DuBois, H. R. Glyde

Several properties of trapped hard-sphere bosons are evaluated using variational Monte Carlo techniques. A trial wave function composed of a renormalized single-particle Gaussian and a hard-sphere Jastrow function for pair correlations is used to study the sensitivity of condensate and noncondensate properties to the hard-sphere radius and the number of particles. Special attention is given to diagonalizing the one-body density matrix and obtaining the corresponding single-particle natural orbitals and their occupation numbers for the system. The condensate wave function and condensate fraction are then obtained from the single-particle orbital with the highest occupation. The effect of interaction on other quantities such as the ground-state energy, the mean radial displacement, and the momentum distribution is calculated as well. Results are compared with mean-field theory in the dilute limit.

arXiv:cond-mat/0008368•Physical Review A•Jan 5, 2001•83 citations

Many-body density and coherence of trapped cold bosons

Camille Lévêque, Fritz Diorico, Jörg Schmiedmayer +1 more

Many-body densities and correlation functions are of paramount importance for understanding quantum many-body physics. Here, we present a method to compute them; our approach is general and based on the action of bosonic or fermionic annihilation field operators on the many-body wavefunction. We analyze $N = 6$ quasi-one-dimensional harmonically-trapped bosons with weak to strong contact interaction strength up to the Tonks-Girardeau limit with infinite repulsion using the MultiConfigurational Time-Dependent Hartree method for indistinguishable particles (MCTDH-X). We compare our MCTDH-X solutions to the analytical ones in the infinite repulsion regime as well as to the so-called correlated pair wavefunction approach and find a good agreement. Since numerical approximations are not bound to the cases where analytical solutions are known, we thus demonstrate a general method to investigate high-order reduced density matrices and correlation functions in systems for which analytical solutions are unknown. We trace the build-up of correlation features in the crossover from weak interactions to the Tonks-Girardeau limit and find that the higher-order correlation functions and densities resemble those in the Tonks-Girardeau limit for way smaller interactions than anticipated from just the one-body density.

arXiv:2006.10755•Jun 18, 2020•8 citations

Correlations in ultracold trapped few-boson systems: Transition from condensation to fermionization

Sascha Zöllner, Hans‐Dieter Meyer, Peter Schmelcher

We study the correlation properties of the ground states of few ultracold bosons, trapped in double wells of varying barrier height in one dimension. Extending previous results on the signature of the transition from a Bose-condensed state via fragmentation to the hard-core limit, we provide a deeper understanding of that transition by relating it to the loss of coherence in the one-body density matrix and to the emerging long-range tail in the momentum spectrum. These are accounted for in detail by discussing the natural orbitals and their occupations. Our discussion is complemented by an analysis of the two-body correlation function.

arXiv:cond-mat/0608639•Physical Review A•Dec 8, 2006•82 citations

Fragmentation and correlations in a rotating Bose–Einstein condensate undergoing breakup

Sunayana Dutta, Axel U. J. Lode, Ofir E. Alon

The theoretical investigation of rotating Bose-Einstein condensates has mainly focused on the emergence of quantum vortex states and the condensed properties of such systems. In the present work, we concentrate on other facets by examining the impact of rotation on the ground state of weakly interacting bosons confined in anharmonic potentials computed both at the mean-field level and particularly at the many-body level of theory. For the many-body computations, we employ the well-established many-body method known as the multiconfigurational time-dependent Hartree method for bosons. We present how various degrees of fragmentation can be generated following the breakup of the ground state densities in anharmonic traps without ramping up a potential barrier for strong rotations. The breakup of the densities is found to be associated with the acquisition of angular momentum in the condensate due to the rotation. In addition to fragmentation, the presence of many-body correlations is examined by computing the variances of the many-particle position and momentum operators. For strong rotations, the many-body variances become smaller than their mean-field counterparts, and one even finds a scenario with opposite anisotropies of the mean-field and many-body variances. Further, it is observed that for higher discrete symmetric systems of order k, namely three-fold and four-fold symmetry, breakup to k sub-clouds and emergence of k-fold fragmentation take place. All in all, we provide a thorough many-body investigation of how and which correlations build up when a trapped Bose-Einstein condensate breaks up under rotation.

arXiv:2211.03176•Scientific Reports•Feb 27, 2023•16 citations
The first- and second-order correlation functions of trapped interacting Bose-Einstein condensates are investigated numerically on a many-body level from first principles. Correlations in real space and momentum space are treated. The coherence properties are analyzed. The results are obtained by solving the many-body Schr"odinger equation. It is shown in an example how many-body effects can be induced by the trap geometry. A generic fragmentation scenario of a condensate is considered. The correlation functions are discussed along a pathway from a single condensate to a fragmented condensate. It is shown that strong correlations can arise from the geometry of the trap, even at weak interaction strengths. The natural orbitals and natural geminals of the system are obtained and discussed. It is shown how the fragmentation of the condensate can be understood in terms of its natural geminals. The many-body results are compared to those of mean-field theory. The best solution within mean-field theory is obtained. The limits in which mean-field theories are valid are determined. In these limits the behavior of the correlation functions is explained within an analytical model.
Abstract Many-body properties of a fermionic impurity embedded in a Bose-Einstein condensate are analyzed analytically using a solvable model, the harmonic-interaction model for Bose-Fermi mixtures. The one-particle and two-particle densities, reduced density matrices, and correlation functions of the fermions and bosons, both in position and momentum spaces, are prescribed in closed form. The various coherence lengths are analyzed. We show that the first-order coherence lengths in position and momentum spaces are equal whereas the second-order quantities can differ substantially. Illustrative examples where the sole interaction is between the impurity and the condensate are presented. Implications are briefly discussed.
In this thesis, the physics of trapped, interacting Bose-Einstein condensates is analyzed by solving the many-body Schrödinger equation. Particular emphasis is put on coherence, fragmentation and reduced density matrices. First, the ground state of a trapped Bose-Einstein condensate and its correlation functions are obtained. Then the dynamics of a bosonic Josephson junction is investigated by solving the time-dependent many-body Schrödinger equation numerically exactly. These are the first exact results in literature in this context. It is shown that the standard approximations of the field, Gross-Pitaevskii theory and the Bose-Hubbard model fail at weak interaction strength and within their range of expected validity. For stronger interactions the dynamics becomes strongly correlated and a new equilibration phenomenon is discovered. By comparison with exact results it is shown that a symmetry of the Bose-Hubbard model between attractive and repulsive interactions must be considered an artefact of the model. A conceptual innovation of this thesis are time-dependent Wannier functions. Equations of motion for time-dependent Wannier functions are derived from the variational principle. By comparison with exact results it is shown that lattice models can be greatly improved at little computational cost by letting the Wannier functions of a lattice model become time-dependent.
Several properties of trapped hard-sphere bosons are evaluated using variational Monte Carlo techniques. A trial wave function composed of a renormalized single-particle Gaussian and a hard-sphere Jastrow function for pair correlations is used to study the sensitivity of condensate and noncondensate properties to the hard-sphere radius and the number of particles. Special attention is given to diagonalizing the one-body density matrix and obtaining the corresponding single-particle natural orbitals and their occupation numbers for the system. The condensate wave function and condensate fraction are then obtained from the single-particle orbital with the highest occupation. The effect of interaction on other quantities such as the ground-state energy, the mean radial displacement, and the momentum distribution is calculated as well. Results are compared with mean-field theory in the dilute limit.
Many-body densities and correlation functions are of paramount importance for understanding quantum many-body physics. Here, we present a method to compute them; our approach is general and based on the action of bosonic or fermionic annihilation field operators on the many-body wavefunction. We analyze N=6N = 6N=6 quasi-one-dimensional harmonically-trapped bosons with weak to strong contact interaction strength up to the Tonks-Girardeau limit with infinite repulsion using the MultiConfigurational Time-Dependent Hartree method for indistinguishable particles (MCTDH-X). We compare our MCTDH-X solutions to the analytical ones in the infinite repulsion regime as well as to the so-called correlated pair wavefunction approach and find a good agreement. Since numerical approximations are not bound to the cases where analytical solutions are known, we thus demonstrate a general method to investigate high-order reduced density matrices and correlation functions in systems for which analytical solutions are unknown. We trace the build-up of correlation features in the crossover from weak interactions to the Tonks-Girardeau limit and find that the higher-order correlation functions and densities resemble those in the Tonks-Girardeau limit for way smaller interactions than anticipated from just the one-body density.
We study the correlation properties of the ground states of few ultracold bosons, trapped in double wells of varying barrier height in one dimension. Extending previous results on the signature of the transition from a Bose-condensed state via fragmentation to the hard-core limit, we provide a deeper understanding of that transition by relating it to the loss of coherence in the one-body density matrix and to the emerging long-range tail in the momentum spectrum. These are accounted for in detail by discussing the natural orbitals and their occupations. Our discussion is complemented by an analysis of the two-body correlation function.
The theoretical investigation of rotating Bose-Einstein condensates has mainly focused on the emergence of quantum vortex states and the condensed properties of such systems. In the present work, we concentrate on other facets by examining the impact of rotation on the ground state of weakly interacting bosons confined in anharmonic potentials computed both at the mean-field level and particularly at the many-body level of theory. For the many-body computations, we employ the well-established many-body method known as the multiconfigurational time-dependent Hartree method for bosons. We present how various degrees of fragmentation can be generated following the breakup of the ground state densities in anharmonic traps without ramping up a potential barrier for strong rotations. The breakup of the densities is found to be associated with the acquisition of angular momentum in the condensate due to the rotation. In addition to fragmentation, the presence of many-body correlations is examined by computing the variances of the many-particle position and momentum operators. For strong rotations, the many-body variances become smaller than their mean-field counterparts, and one even finds a scenario with opposite anisotropies of the mean-field and many-body variances. Further, it is observed that for higher discrete symmetric systems of order k, namely three-fold and four-fold symmetry, breakup to k sub-clouds and emergence of k-fold fragmentation take place. All in all, we provide a thorough many-body investigation of how and which correlations build up when a trapped Bose-Einstein condensate breaks up under rotation.