We develop a coupled-cluster theory for bosonic mixtures of binary species in external traps, providing a promising theoretical approach to demonstrate highly accurately the many-body physics of mixtures of Bose-Einstein condensates. The coupled-cluster wavefunction for the binary species is obtained when an exponential cluster operator eT, where T = T(1) + T(2) + T(12) and T(1) accounts for excitations in species-1, T(2) for excitations in species-2, and T(12) for combined excitations in both species, acts on the ground state configuration prepared by accumulating all bosons in a single orbital for each species. We have explicitly derived the working equations for bosonic mixtures by truncating the cluster operator up to the single and double excitations and using arbitrary sets of orthonormal orbitals for each of the species. Furthermore, the comparatively simplified version of the working equations are formulated using the Fock-like operators. Finally, using an exactly solvable many-body model for bosonic mixtures that exists in the literature allows us to implement and test the performance and accuracy of the coupled-cluster theory for situations with balanced as well as imbalanced boson numbers and for weak to moderately strong intra- and interspecies interaction strengths. The comparison between our computed results using coupled-cluster theory with the respective analytical exact results displays remarkable agreement exhibiting excellent success of the coupled-cluster theory for bosonic mixtures. All in all, the correlation exhaustive coupled-cluster theory shows encouraging results and could be a promising approach in paving the way for high-accuracy modeling of various bosonic mixture systems.
We develop a coupled-cluster theory for bosonic mixtures of binary species in external traps, providing a promising theoretical approach to demonstrate highly accurately the many-body physics of mixtures of Bose-Einstein condensates. The coupled-cluster wavefunction for the binary species is obtained when an exponential cluster operator eT, where T = T(1) + T(2) + T(12) and T(1) accounts for excitations in species-1, T(2) for excitations in species-2, and T(12) for combined excitations in both species, acts on the ground state configuration prepared by accumulating all bosons in a single orbital for each species. We have explicitly derived the working equations for bosonic mixtures by truncating the cluster operator up to the single and double excitations and using arbitrary sets of orthonormal orbitals for each of the species. Furthermore, the comparatively simplified version of the working equations are formulated using the Fock-like operators. Finally, using an exactly solvable many-body model for bosonic mixtures that exists in the literature allows us to implement and test the performance and accuracy of the coupled-cluster theory for situations with balanced as well as imbalanced boson numbers and for weak to moderately strong intra- and interspecies interaction strengths. The comparison between our computed results using coupled-cluster theory with the respective analytical exact results displays remarkable agreement exhibiting excellent success of the coupled-cluster theory for bosonic mixtures. All in all, the correlation exhaustive coupled-cluster theory shows encouraging results and could be a promising approach in paving the way for high-accuracy modeling of various bosonic mixture systems.
A coupled-cluster approach for systems of N bosons in external traps is developed. In the coupled-cluster approach the exact many-body wave function is obtained by applying an exponential operator expT to the ground configuration $\ensuremath{\mid}{\ensuremath{\phi}}_{0}⟩$. The natural ground configuration for bosons is, of course, when all reside in a single orbital. Because of this simple structure of $\ensuremath{\mid}{\ensuremath{\phi}}_{0}⟩$, the appearance of excitation operators $T={\ensuremath{\sum}}_{n=1}^{N}{T}_{n}$ for bosons is much simpler than for fermions. We can treat very large numbers of bosons with coupled-cluster expansions. In a substantial part of this work, we address the issue of size consistency for bosons and enquire whether truncated coupled-cluster expansions are size consistent. We show that, in contrast to the familiar situation for fermions for which coupled-cluster expansions are size consistent, for bosons the answer to this question depends on the choice of ground configuration. Utilizing the natural ground configuration, working equations for the truncated coupled-cluster with T=T1+T2, i.e., coupled-cluster singles doubles are explicitly derived. Finally, an illustrative numerical example for a condensate with up to N=10000 bosons in an harmonic trap is provided and analyzed. The results are highly promising.
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.
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.
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.
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.
This PhD thesis contains new results on the mathematical study of Bose-Einstein condensation and the main part of it is devoted to mixtures of condensates, i.e., systems composed of multiple bosonic species in interaction. We prove the validity of effective ground state theories for mixtures in the Gross-Pitaevskii and mean-field regime. We show that the ground state energy asymptotics, in the large-N limit, is captured by the minimum of a suitable one-body functional. Moreover, we prove that in the ground state all species exhibit Bose-Einstein condensation onto the minimizer of that functional. For mixtures in the mean-field regime, we provide a rigorous justification of Bogoliubov’s theory. This is done by computing the contribution to the ground state energy which is due to excited particles. We also prove a norm approximation for the ground state vector, in the Fock space norm. From the time-dependent viewpoint, we prove for the first time the validity of the effective equations that were previously known due to heuristic physical arguments, and that are confirmed by robust experimental evidence. Our results show that, for mixtures in the Gross-Pitaevskii and mean-field regime, the effective dynamics is governed by a system of non-linear Schrödinger equations, one for each species of the mixture. In the final part of the thesis we present additional results on problems and models related to the study on mixtures. We were able to derive the effective dynamics for spinor- and pseudo-spinor condensates. The equations that we obtain are precisely those of modern experiments with ultra-cold spin bosons. We also show that the mean-field model provide a time-dependent control of condensation that is very accurate for the typical duration times of experiments. A further result is the global well-posedness in the energy space of the singular Hartree equation. Last, we present new remarks on the adaptation of known techniques that one needs in order to prove the derivation of the magnetic Gross-Pitaevskii equation.