In a mixture of two kinds of identical bosons, there are two types of pairs: identical bosons’ pairs, of either species, and pairs of distinguishable bosons. In the present work, the fragmentation of pairs in a trapped mixture of Bose–Einstein condensates is investigated using a solvable model, the symmetric harmonic-interaction model for mixtures. The natural geminals for pairs made of identical or distinguishable bosons are explicitly contracted by diagonalizing the intra-species and inter-species reduced two-particle density matrices, respectively. Properties of pairs’ fragmentation in the mixture are discussed, the role of the mixture’s center-of-mass and relative center-of-mass coordinates is elucidated, and a generalization to higher-order reduced density matrices is made. As a complementary result, the exact Schmidt decomposition of the wave function of the bosonic mixture is constructed. The entanglement between the two species is governed by the coupling of their individual center-of-mass coordinates, and it does not vanish at the limit of an infinite number of particles where any finite-order intra-species and inter-species reduced density matrix per particle is 100% condensed. Implications are briefly discussed.
In a mixture of two kinds of identical bosons, there are two types of pairs: identical bosons’ pairs, of either species, and pairs of distinguishable bosons. In the present work, the fragmentation of pairs in a trapped mixture of Bose–Einstein condensates is investigated using a solvable model, the symmetric harmonic-interaction model for mixtures. The natural geminals for pairs made of identical or distinguishable bosons are explicitly contracted by diagonalizing the intra-species and inter-species reduced two-particle density matrices, respectively. Properties of pairs’ fragmentation in the mixture are discussed, the role of the mixture’s center-of-mass and relative center-of-mass coordinates is elucidated, and a generalization to higher-order reduced density matrices is made. As a complementary result, the exact Schmidt decomposition of the wave function of the bosonic mixture is constructed. The entanglement between the two species is governed by the coupling of their individual center-of-mass coordinates, and it does not vanish at the limit of an infinite number of particles where any finite-order intra-species and inter-species reduced density matrix per particle is 100% condensed. Implications are briefly discussed.
Abstract A mixture of two kinds of identical bosons, species 1 with N 1 bosons of mass m 1 and species 2 with N 2 bosons of mass m 2 , held in a harmonic potential of frequency ω and interacting by harmonic intra-species and inter-species particle-particle interactions of strengths λ 1 , λ 2 , and λ 12 is discussed. This is an exactly-solvable model of a generic mixture of trapped interacting bosons which allows one to investigate and determine analytically properties of interest. To start, closed form expressions for the frequencies, ground-state energy, and wave-function of the mixture are obtained and briefly analyzed as a function of the masses, numbers of particles, and strengths and signs of interactions. To prove Bose–Einstein condensation of the mixture three steps are needed. First, we integrate the all-particle density matrix, employing a four-parameter matrix-recurrence relations, down to the lowest-order intra-species and inter-species reduced density matrices of the mixture. Second, the coupled Gross–Pitaevskii (mean-field) equations of the mixture are solved analytically. Third, we analyze the mixture’s reduced density matrices in the limit of an infinite number of particles of both species 1 and 2 (when the interaction parameters, i.e. the products of the number of particles times the intra-species and inter-species interaction strengths, are held fixed) and prove that: (i) both species 1 and 2 are 100% condensed; (ii) the inter-species reduced density matrix per particle is separable and given by the product of the intra-species reduced density matrices per particle; and (iii) the mixture’s energy per particle, and reduced density matrices and densities per particle all coincide with the Gross–Pitaevskii quantities. Finally, when the infinite-particle limit is taken with respect to, say, species 1 only (with interaction parameters held fixed) we prove that: (iv) only species 1 is 100% condensed and its reduced density matrix and density per particle, as well as the mixture’s energy per particle, coincide with the Gross–Pitaevskii quantities of species 1 alone; and (v) the inter-species reduced density matrix per particle is nonetheless separable and given by the product of the intra-species reduced density matrices per particle. The results are compared and discussed with respect to the recent work by Bouvrie et al (2014 Eur. Phys. J . D 68 346 ) who found for attractive mixtures vanishing bipartite entanglement between the two species and between a single particle (of either kind) and the remaining particles in the mixture. Implications are briefly discussed.
We consider a multi-species mixture of interacting bosons, N1 bosons of mass m1, N2 bosons of mass m2, and N3 bosons of mass m3, in a harmonic trap with frequency ω. The corresponding intra-species interaction strengths are λ11, λ22, and λ33, and the inter-species interaction strengths are λ12, λ13, and λ23. When the shape of all interactions is harmonic, the system corresponds to the generic multi-species harmonic-interaction model, which is exactly solvable. We start by solving the many-particle Hamiltonian and concisely discussing the ground-state wavefunction and energy in explicit forms as functions of all parameters, the masses, numbers of particles, and the intra-species and inter-species interaction strengths. We then explicitly compute the reduced one-particle density matrices for all the species and diagonalize them, thus generalizing the treatment by the authors earlier. The respective eigenvalues determine the degree of fragmentation of each species. As an application, we focus on phenomena that do not arise in the corresponding single-species or two-species systems. For instance, we consider a mixture of two kinds of bosons in a bath made by a third kind, controlling the fragmentation of the former by coupling to the latter. Another example exploits the possibility of different connectivities (i.e., which species interacts with which species) in the mixture, and demonstrates how the fragmentation of species 3 can be manipulated by the interaction between species 1 and species 2, when species 3 and 1 do not interact with each other. We highlight the properties of fragmentation that only appear in the multi-species mixture. Further applications are briefly discussed.
Abstract Fragmentation of bosons and pairs in a trapped imbalanced bosonic mixture is investigated analytically using an exactly solvable model, the generic harmonic-interaction model for mixtures. Closed-form expressions for the eigenvalues and eigenfunctions of the reduced one-particle and two-particle density matrices as a function of all parameters, the masses, numbers of bosons, and the intraspecies and interspecies interactions, are obtained and analyzed. As an application, we consider a system made of N 1 = 100 non-interacting species 1 bosons embedded in a bath made of N 2 = 106 non-interacting species 2 bosons, and show how fragmentation of the system’s bosons and pairs emerges from the system–bath interaction only. Interestingly, the lighter the bosons comprising the bath are the stronger is the system’s fragmentation. Further applications are briefly discussed.
We investigate a trapped mixture of Bose-Einstein condensates consisting of a multiple number of P species. To be able to do so, an exactly solvable many-body model is called into play. This is the P-species harmonic-interaction model. After presenting the Hamiltonian, the ground-state energy and wavefunction are explicitly calculated. All properties of the mixture's ground state can, in principle, be obtained from the many-particle wavefunction. A scheme to integrate the all-particle density matrix is derived and implemented, leading to closed-form expressions for the reduced one-particle density matrices. Of particular interest is the infinite-particle-number limit, which is obtained when the numbers of bosons are taken to infinity while keeping the interaction parameters fixed. We first prove that at the infinite-particle-number limit all the species are 100% condensed. The mean-field solution of the P-species mixture is also obtained analytically and is used to show that the energy per particle and densities per particle computed at the many-body level of theory boil down to their mean-field counterparts. Despite these, correlations in the mixture exist at the infinite-particle-number limit. To this end, we obtain closed-form expressions for the correlation energy, namely, the difference between the mean-field and many-body energies, and the depletion of the species, i.e., the number of particles residing outside the condensed modes, at the infinite-particle-number limit. The depletion and the correlation energy per species are shown to critically depend on the number of species. Of separate interest is the entanglement between one species of bosons and the other P - 1 species. This quantity is governed by the coupling of the center-of-mass coordinates of the species and is obtained by the respective Schmidt decomposition of the P-species wavefunction. Interestingly, there is an optimal number of species, here P = 3, where the entanglement is maximal. Importantly, the manifestation of this interspecies entanglement in an observable is possible. It is the position-momentum uncertainty product of one species in the presence of the other P - 1 species, which is derived and demonstrated to correlate with the interspecies entanglement. All in all, we show and explain how correlations at the infinite-particle-number limit of a trapped multiple-species bosonic mixture depend on the interactions and how they evolve with the number of species. Generalizations and implications are briefly discussed.
A solvable model of a generic trapped bosonic mixture, N 1 bosons of mass m 1 and N 2 bosons of mass m 2 trapped in an harmonic potential of frequency ω and interacting by harmonic inter-particle interactions of strengths λ 1 , λ 2 , and λ 12 , is discussed. It has recently been shown for the ground state [J. Phys. A 50 , 295002 (2017)] that in the infinite-particle limit, when the interaction parameters λ 1 ( N 1 − 1), λ 2 ( N 2 − 1), λ 12 N 1 , λ 12 N 2 are held fixed, each of the species is 100% condensed and its density per particle as well as the total energy per particle are given by the solution of the coupled Gross-Pitaevskii equations of the mixture. In the present work we investigate properties of the trapped generic mixture at the infinite-particle limit, and find differences between the many-body and mean-field descriptions of the mixture, despite each species being 100%. We compute analytically and analyze, both for the mixture and for each species, the center-of-mass position and momentum variances, their uncertainty product, the angular-momentum variance, as well as the overlap of the exact and Gross-Pitaevskii wavefunctions of the mixture. The results obtained in this work can be considered as a step forward in characterizing how important are many-body effects in a fully condensed trapped bosonic mixture at the infinite-particle limit.
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.