We derive an approximate analytical expression for the ground state of double-well BEC's, which reproduces highly accurately the numerical solution for the whole parameter regimes of the two-body repulsive interaction strength, the total number of atoms, and the hopping parameter.
The two-site version of the Bose-Hubbard model can be used as a simplified model of N bosons confined to two quantum wells. We study the behavior of low-temperature Bose-Einstein condensates as a function of the strength of the repulsive interactions between the bosons and of the size of the quantum-mechanical tunneling between the quantum wells. The ground-state of weakly interacting bosons tends to be one in which almost all particles are in the same quantum-mechanical state (‘Bose-Einstein Condensate’, or ‘BEC’). It is found that even very small amounts of interwell tunneling lead to a BEC, but that zero tunneling leads to two separate BECs, one in each well (‘fragmentation’). We study systematically the transition between a single BEC comprised of bosons in both wells, and a fragmented condensate consisting of two BECs, one for each well. Our analysis will use matrix computational methods to solve the Schrödinger equation of the two-well boson system, and obtain the energies and wavefunctions as eigenvalues and eigenvectors of the Hamiltonian matrix. We will compare the results obtained with previous analyses, and against analytical approximations developed for the limiting case of a large number of bosons (large N).
We present a method to solve the model describing either a Bose-Einstein condensate (BEC) in a Mott-insulator state or a double-well BEC. We show that all the energy eigenvalues and eigenstates for an arbitrary (small or large) total atom number N can be explicitly expressed analytically in terms of a parameter $\ensuremath{\lambda}$ whose values are determined by the roots of the polynomials of the order of at most 1+int(N/2), with int(x) denoting x's integer part. We also show that \ensuremath{\lambda}'s explicit analytical expressions for N<7 can be readily obtained by a simple MATHEMATICA code. Besides, finding the roots of the polynomials of the order of at most 1+int(N/2) to give explicitly all the energy eigenvalues and eigenstates greatly simplifies the corresponding calculations, particularly when the total atom number N is large.
By developing an approach by which we are able to quickly obtain spectra and eigenstates, we reveal for what is believed to be the first time the two novel phenomena of magic numbers and erratic level crossings in double-well Bose-Einstein condensates of N atoms. For N < or = 27 and values of U/J that are not too small (U is the two-body interaction strength, and J is the hopping parameter), systems with even atoms are shown to be much more stable than those with odd atoms, and hence even integers play a role in such systems as if they were the magic numbers of nuclei. For N > or = 30, erratic level crossings occur as NU >> J.
Dynamics of the repulsive Bose–Einstein condensate (BEC) in a double-well trap is explored within the 3D time-dependent Gross–Pitaevskii equation. The model avoids numerous common approximations (two-mode treatment, time-space factorization, fixed values of the chemical potential and barrier penetrability, etc) and thus may provide a realistic description of BEC dynamics, including both weak-coupling (sub-barrier) and strong-coupling (above-barrier) regimes and their crossover. The strong coupling regime is achieved by increasing the number N of BEC atoms and thus the chemical potential. The evolution with N of Josephson oscillations (JO) and macroscopic quantum self-trapping (MQST) is examined and the crucial impact of the BEC interaction is demonstrated. At weak coupling, the calculations well reproduce the JO/MQST experimental data. At strong coupling with a significant overlap of the left and right BECs, we observe a remarkable persistence of the Josephson-like dynamics: the JO and MQST converge to a high-frequency JO-like mode where both population imbalance and phase difference oscillate around the zero averages.
We present a theoretical study of the ground state of a Bose-Einstein condensate with repulsive interparticle interactions in a double-well potential, using a restricted variational principle. Within such an approach, there is a transition from a single condensate to a fragmented condensate as the strength of the central barrier of the potential is increased. We determine the nature of this transition through an approximate analytic solution as well as a numerical solution of our model, in the regime where the interparticle interactions can be treated perturbatively. The degree of fragmentation of the condensate is characterized by the degrees of first- and second-order spatial coherence across the barrier.