Solvable Model of a Generic Driven Mixture of Trapped Bose-Einstein Condensates and Properties of a Many-Boson Floquet State at the Limit of an Infinite Number of Particles

Coupled-cluster theory for trapped bosonic mixtures

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.

Fragmentation of Identical and Distinguishable Bosons’ Pairs and Natural Geminals of a Trapped Bosonic Mixture

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.