We investigate a noninteracting many-particle bosonic system, placed in an asymmetric double-well potential. We first consider the dynamics of a single particle and determine its time-dependent probabilities to be in the left or the right well of the potential. These probabilities obey the standard Josephson equations, which in their many-particle interpretation also describe a globally coherent system, such as a Bose-Einstein condensate. This system exhibits the widely studied Josephson oscillations of the population imbalance between the wells. In our paper we go beyond the regime of global coherence by developing a formalism based on an effective density matrix. This formalism gives rise to a generalization of Josephson equations, which differ from the standard ones by an additional parameter, that has the meaning of the degree of fragmentation. We first consider the solution of the generalized Josephson equations in the particular case of thermal equilibrium at finite temperatures, and extend our discussion to the nonequilibrium regime afterwards. Our model leads to a constraint on the maximum amplitude of Josephson oscillations for a given temperature and the total number of particles. A detailed analysis of this constraint for typical experimental scenarios is given.
We investigate a noninteracting many-particle bosonic system, placed in an asymmetric double-well potential. We first consider the dynamics of a single particle and determine its time-dependent probabilities to be in the left or the right well of the potential. These probabilities obey the standard Josephson equations, which in their many-particle interpretation also describe a globally coherent system, such as a Bose-Einstein condensate. This system exhibits the widely studied Josephson oscillations of the population imbalance between the wells. In our paper we go beyond the regime of global coherence by developing a formalism based on an effective density matrix. This formalism gives rise to a generalization of Josephson equations, which differ from the standard ones by an additional parameter, that has the meaning of the degree of fragmentation. We first consider the solution of the generalized Josephson equations in the particular case of thermal equilibrium at finite temperatures, and extend our discussion to the nonequilibrium regime afterwards. Our model leads to a constraint on the maximum amplitude of Josephson oscillations for a given temperature and the total number of particles. A detailed analysis of this constraint for typical experimental scenarios is given.
Abstract The out-of-equilibrium quantum dynamics of an interacting Bose gas trapped in a one-dimensional asymmetric double-well potential is studied by solving the many-body Schrödinger equation numerically accurately. We examine how the gradual loss of symmetry of the confining trap affects the macroscopic quantum tunneling dynamics of the system between the two wells. In an asymmetric double well, the two wells are not equivalent anymore, say, the left well is deeper than the right one. Accordingly, we analyze the dynamics by initially preparing the condensate in both the left and the right wells. The dynamics of the system is characterized by the time evolution of a few physical quantities of increasing many-body complexity, namely, the survival probability, depletion and fragmentation, and the many-particle position and momentum variances. In particular, we have examined the frequencies and amplitudes of the oscillations of the survival probabilities, the time scale for the development of fragmentation and its degree, and the growth and oscillatory behavior of the many-particle position and momentum variances. There is an overall suppression of the oscillations of the survival probabilities in an asymmetric double well, except for resonant values of asymmetry for which the one-body ground state energy in the right well matches with one of the one-body excited states in the left well, thereby resulting in resonantly enhanced tunneling from the right well ground state. Overall, depending on whether the condensate is initially prepared in the left or right well, the repulsive inter-atomic interactions affect the survival probabilities differently. For a sufficiently strong repulsive interaction, the system is found to become fragmented. The degree of fragmentation depends both on the asymmetry of the trap and the initial well in which the condensate is prepared in a non-trivial manner. Furthermore, we show that the phenomenon of resonantly enhanced tunneling can be accompanied by a large degree of fragmentation (depletion) for the strong (weak) interaction. The many-particle position and momentum variances follow the density oscillations of the system in the asymmetric double well and bears prominent signatures of the degree of depletion or fragmentation, depending on the strength of the interactions. These quantities further exhibit a fine structure signifying a breathing-mode oscillation. Finally, a universality of fragmentation for systems made of different numbers of particles but the same interaction parameter is also found and its dependence on the asymmetry is investigated. The phenomenon is robust despite the asymmetry of the junction and admits a macroscopically-large fragmented condensate characterized by a diverging many-particle position variance. This is as far as one can get from the dynamics of the density in the junction.
We investigate the finite-temperature properties of a bosonic Josephson junction composed of N interacting atoms confined by a quasi-one-dimensional asymmetric double-well potential, modeled by the two-site Bose-Hubbard Hamiltonian. We compute numerically the spectral decomposition of the statistical ensemble of states, the thermodynamic and entanglement entropies, the population imbalance, the quantum Fisher information, and the coherence visibility. We analyze their dependence on the system parameters, showing in particular how finite temperature and on-site energy asymmetry affect the entanglement and coherence properties of the system. Moreover, starting from a quantum phase model which accurately describes the system over a wide range of interactions, we develop a reliable description of the strong tunneling regime, where thermal averages may be computed analytically using a modified Boltzmann weight involving an effective temperature. We discuss the possibility of applying this effective description to other models in suitable regimes.
In a double-well potential, a Bose-Einstein condensate exhibits Josephson oscillations or self-trapping, depending on its initial preparation and on the ratio of interparticle interaction to interwell tunneling. Here, we elucidate the role of the exchange symmetry for the dynamics with a mixture of two distinguishable species with identical physical properties, that is, which are governed by an isospecific interaction and external potential. In the mean-field limit, the spatial population imbalance of the mixture can be described by the dynamics of a single species in an effective potential with modified properties or, equivalently, with an effective total particle number. The oscillation behavior can be tuned by populating the second species while maintaining the spatial population imbalance and all other parameters constant. In the corresponding many-body approach, the single-species description approximates the full counting statistics well also outside the realm of spin-coherent states. The method is extended to general Bose-Hubbard systems and to their classical mean-field limits, which suggests an effective single-species description of multicomponent Bose gases with weakly an-isospecific interactions.
We unravel the out-of-equilibrium quantum dynamics of a few interacting bosonic clouds in a two-dimensional asymmetric double-well potential at the resonant tunneling scenario. At the single-particle level of resonant tunneling, particles tunnel under the barrier from, typically, the ground-state in the left well to an excited state in the right well, i.e., states of different shapes and properties are coupled when their one-particle energies coincide. In two spatial dimensions, two types of resonant tunneling processes are possible, to which we refer to as longitudinal and transversal resonant tunneling. Longitudinal resonant tunneling implies that the state in the right well is longitudinally-excited with respect to the state in the left well, whereas transversal resonant tunneling implies that the former is transversely-excited with respect to the latter. We show that interaction between bosons makes resonant tunneling phenomena in two spatial dimensions profoundly rich, and analyze these phenomena in terms of the loss of coherence of the junction and development of fragmentation, and coupling between transverse and longitudinal degrees-of-freedom and excitations. To this end, a detailed analysis of the tunneling dynamics is performed by exploring the time evolution of a few physical quantities, namely, the survival probability, occupation numbers of the reduced one-particle density matrix, and the many-particle position, momentum, and angular-momentum variances. To accurately calculate these physical quantities from the time-dependent many-boson wavefunction, we apply a well-established many-body method, the multiconfigurational time-dependent Hartree for bosons (MCTDHB), which incorporates quantum correlations exhaustively. By comparing the survival probabilities and variances at the mean-field and many-body levels of theory and investigating the development of fragmentation, we identify the detailed mechanisms of many-body longitudinal and transversal resonant tunneling in two dimensional asymmetric double-wells. In particular, we find that the position and momentum variances along the transversal direction are almost negligible at the longitudinal resonant tunneling, whereas they are substantial at the transversal resonant tunneling which is caused by the combination of the density and breathing mode oscillations. We show that the width of the interparticle interaction potential does not affect the qualitative physics of resonant tunneling dynamics, both at the mean-field and many-body levels. In general, we characterize the impact of the transversal and longitudinal degrees-of-freedom in the many-boson tunneling dynamics at the resonant tunneling scenarios.
This topical review provides an overview of the key theoretical features of Bose–Einstein condensates (BECs) in cold atomic gases at near zero temperature in the situation where all the bosons occupy at most two single particle states or modes. This situation applies to single-component BECs in double well trap potentials and to two-component BEC in single well trap potentials, such as occur when BEC are used in interferometry experiments. The Hamiltonian is introduced in terms of field operators and mode expansions are restricted to a total of two modes. Spin operators and their eigenstates are introduced as the fundamental basis states for describing the two-mode N boson quantum system. The spin states have a macroscopic angular momentum quantum number of N/2 and the magnetic quantum number k specifies the relative number of bosons in the two modes. The treatment presented involves an extensive use of angular momentum theory, including unitary rotation operators. Important states of the two-mode system such as binomial or coherent states, relative phase eigenstates are discussed. Boson position measurements are specified via quantum correlation functions, and the use of these functions in describing coherence properties, interference patterns and fragmentation effects in BECs is presented. The Bloch vector is defined and related to the quantum correlation functions, with quantum fluctuations of the Bloch vector being treated in terms of the covariance matrix. Applications to important two-mode states are made. Spin squeezing is discussed. Based on applying variational principles, the general dynamical behaviour of the two-mode BEC is determined via generalised Gross–Pitaevskii equations for the modes and matrix mechanics equations for the probability amplitudes of the relative number basis states, the mode and amplitude equations being coupled and self-consistent. The single mode equations are also presented. The Hamiltonian is written in terms of the spin operators and the Josephson Hamiltonian obtained as a simplification in which the dynamical behaviour of the mode functions is ignored – for the one-component case the mode functions are also required to be localised and separate. Coefficients in the Josephson Hamiltonian describe tunneling/intercomponent coupling, asymmetry and collisions and these are defined via integrals involving the mode functions. The Josephson model involves using the Josephson Hamiltonian to give simple predictions of the energy states and dynamical behaviour of the two-mode system, dynamical effects on the mode functions being ignored. The three regimes – Rabi, Josephson and Fock are described, and the energy states obtained for the Fock and Rabi regimes. Dynamical behaviour treatments based on the Josephson model are outlined. In the situation where all bosons are in the same single particle state, semi-classical Bloch equations are derived and their solutions given in terms of elliptic functions. The quantum regime is treated using matrix mechanics equations for the probability amplitudes. Two representative applications of the Josephson model dynamics are treated, with graphs showing the results of numerical work being displayed. The first is in describing Heisenberg limited BEC interferometry for a single-component BEC in a double well, the treatment showing collapses and revivals in the probability distribution for the relative phase. The second treats Ramsey interferometry for a two-component BEC in a single well, the study revealing that oscillations of the Bloch vector collapse and revive, with the Bloch vector's departure from the Bloch sphere during the collapse period revealing that the BEC has fragmented. In both cases collisions cause the dephasing effects that result in the collapse, revival phenomena. The review ends with a brief outline of phase space and other approaches that extend the treatment beyond the two-mode theory, enabling decoherence effects associated with bosons in non-condensate modes to be studied. A summary of the review contents is included. Detailed mathematical derivations are included in several appendices, available as online supplementary material.
The time-dependent non-equilibrium dynamics of a Bose–Einstein condensate (BEC) typically generates incoherent excitations out of the condensate due to the finite frequencies present in the time evolution. We present a detailed derivation of a general non-equilibrium Greenʼs function technique that describes the coupled time evolution of an interacting BEC and its single-particle excitations in a trap, based on an expansion in terms of the exact eigenstates of the trap potential. We analyze the dynamics of a Bose system in a small double-well potential with initially all particles in the condensate. When the trap frequency is larger than the Josephson frequency, , the dynamics changes at a characteristic time, , abruptly from the slow Josephson oscillations of the BEC to fast Rabi oscillations driven by quasiparticle excitations in the trap. For times , the Josephson oscillations are undamped, in agreement with the experiments. We analyze the physical origin of the finite scale as well as its dependence on the trap parameter Δ .