We investigate the fully quantum evolution of the population imbalance in a perfectly symmetric Bose-Josephson junction modeled by a two-mode Bose-Hubbard Hamiltonian, focusing on the validity of macroscopic quantum self-trapping beyond the mean-field theory. We show that for any finite number of particles, the exact quantum dynamics leads to the breakdown of macroscopic quantum self-trapping after a finite time, regardless of the initial state. Using the symmetries of the Bose-Hubbard Hamiltonian, we provide a mathematical demonstration of this result and analyze the spectral properties governing the dynamics. We identify a branching behavior in the eigenvalue differences and a nontrivial structure of the population-imbalance amplitudes. These features allow us to distinguish two clearly different dynamical regimes and to elucidate the mechanism leading to the emergence of a quasi–macroscopic quantum self-trapping regime for large particle numbers. These findings bridge the gap between mean-field predictions and exact quantum dynamics and provide insight into the emergence of classical nonlinear behavior from finite quantum many-body systems.
We investigate the fully quantum evolution of the population imbalance in a perfectly symmetric Bose-Josephson junction modeled by a two-mode Bose-Hubbard Hamiltonian, focusing on the validity of macroscopic quantum self-trapping beyond the mean-field theory. We show that for any finite number of particles, the exact quantum dynamics leads to the breakdown of macroscopic quantum self-trapping after a finite time, regardless of the initial state. Using the symmetries of the Bose-Hubbard Hamiltonian, we provide a mathematical demonstration of this result and analyze the spectral properties governing the dynamics. We identify a branching behavior in the eigenvalue differences and a nontrivial structure of the population-imbalance amplitudes. These features allow us to distinguish two clearly different dynamical regimes and to elucidate the mechanism leading to the emergence of a quasi–macroscopic quantum self-trapping regime for large particle numbers. These findings bridge the gap between mean-field predictions and exact quantum dynamics and provide insight into the emergence of classical nonlinear behavior from finite quantum many-body systems.
We investigate finite-size quantum effects in the dynamics of N bosonic particles which are tunneling between two sites adopting the two-site Bose-Hubbard model. By using time-dependent atomic coherent states (ACSs) we extend the standard mean-field equations of this bosonic Josephson junction, which are based on time-dependent Glauber coherent states. In this way we find 1/N corrections to familiar mean-field (MF) results: the frequency of macroscopic oscillation between the two sites, the critical parameter for the dynamical macroscopic quantum self-trapping (MQST), and the attractive critical interaction strength for the spontaneous symmetry breaking (SSB) of the ground state. To validate our analytical results we perform numerical simulations of the quantum dynamics. In the case of Josephson oscillations around a balanced configuration we also find that, for a few atoms, the numerical results are in good agreement with the predictions of time-dependent ACS variational approach, provided that the time evolution is not too long. Also, the numerical results of SSB are better reproduced by the ACS approach than with the MF approach. Instead, the onset of MQST is correctly reproduced by ACS theory only in the large N regime and, for this phenomenon, the 1/N correction to the MF formula is not reliable.
Understanding how coherent quantum dynamics give way to correlation-dominated behavior in low-dimensional systems remains a central challenge in quantum many-body physics. Here, we investigate a one-dimensional Bose-Josephson junction confined in a box trap using the multiconfigurational time-dependent Hartree method for bosons (MCTDHB). By varying the interaction strength and initial population imbalance, we identify distinct dynamical regimes governed by the competition between coherence and correlation-induced fragmentation. Weak interactions support coherent Josephson oscillations, whereas increasing imbalance leads to damping. At intermediate interaction strength, varying only the initial imbalance induces a crossover from nearly pure coherent oscillations to many-body dephasing with collapse-and-revival dynamics, and ultimately to equilibration accompanied by strong fragmentation and the saturation of many-body observables. In the strongly interacting regime, the system enters a dynamical freezing regime characterized by pronounced fragmentation, well-separated particle-resolved density peaks, and strongly suppressed tunneling. A systematic comparison with the Bose-Hubbard model reveals excellent agreement in the weakly interacting regime, while progressively larger deviations emerge as higher-orbital occupations beyond the two-mode approximation become significant. These results provide a unified picture of the emergence and competition of coherence, many-body dephasing, equilibration, and dynamical freezing, while delineating the regime of validity of the Bose-Hubbard description.
We study atomic Josephson junctions (AJJs) with one and two bosonic species confined by a double-well potential. Proceeding from the second quantized Hamiltonian, we show that it is possible to describe the zero-temperature AJJ microscopic dynamics by means of extended Bose–Hubbard (EBH) models, which include usually neglected nonlinear terms. Within the mean-field approximation, the Heisenberg equations derived from such two-mode models provide a description of AJJ macroscopic dynamics in terms of ordinary differential equations (ODEs). We discuss the possibility of distinguishing the Rabi, Josephson and Fock regimes in terms of the macroscopic parameters which appear in the EBH Hamiltonians, and then in the ODEs. We compare the predictions for the relative populations of the Bose gas atoms in the two wells obtained from the numerical solutions of the two-mode ODEs, with those deriving from the direct numerical integration of the Gross–Pitaevskii equations (GPEs). Our investigations show that the nonlinear terms of the ODEs are crucial to achieve a good agreement between the ODE and GPE approaches, and in particular to give quantitative predictions of the self-trapping regime.
We study atomic Josephson junctions (AJJs) with one and two bosonic species confined by a double-well potential. Proceeding from the second quantized Hamiltonian, we show that it is possible to describe the zero-temperature AJJ microscopic dynamics by means of extended Bose–Hubbard (EBH) models, which include usually neglected nonlinear terms. Within the mean-field approximation, the Heisenberg equations derived from such two-mode models provide a description of AJJ macroscopic dynamics in terms of ordinary differential equations (ODEs). We discuss the possibility of distinguishing the Rabi, Josephson and Fock regimes in terms of the macroscopic parameters which appear in the EBH Hamiltonians, and then in the ODEs. We compare the predictions for the relative populations of the Bose gas atoms in the two wells obtained from the numerical solutions of the two-mode ODEs, with those deriving from the direct numerical integration of the Gross–Pitaevskii equations (GPEs). Our investigations show that the nonlinear terms of the ODEs are crucial to achieve a good agreement between the ODE and GPE approaches, and in particular to give quantitative predictions of the self-trapping regime.
In this paper we consider a bosonic Josephson junction described by a two-mode Bose-Hubbard model, and we thoroughly analyze a quantum phase transition occurring in the system in the limit of infinite bosonic population. We discuss the relation between this quantum phase transition and the dynamical bifurcation occurring in the spectrum of the discrete self-trapping equations describing the system at the semiclassical level. In particular, we identify five regimes depending on the strength of the effective interaction among bosons, and study the finite-size effects arising from the finiteness of the bosonic population. We devote special attention to the critical regime which reduces to the dynamical bifurcation point in the thermodynamic limit of infinite bosonic population. Specifically, we highlight an anomalous scaling in the population imbalance between the two wells of the trapping potential, as well as in two quantities borrowed from quantum information theory, i.e., the entropy of entanglement and the ground-state fidelity. Our analysis is not limited to the zero-temperature case, but considers thermal effects as well.
Abstract Typical treatments of superconducting or superfluid Josephson junctions rely on mean-field or two-mode models; we explore many-body dynamics of an isolated, ultracold, Bose-gas long Josephson junction using time-evolving block decimation simulations. We demonstrate that with increasing repulsive interaction strength, localized dynamics emerge that influence macroscopic condensate behavior and can lead to the formation of solitons that directly oppose the symmetry of the junction. Initial state population and phase yield insight into dynamic tunneling regimes of a quasi one-dimensional double well potential, from Josephson oscillations to macroscopic-self-trapping. Population imbalance simulations reveal substantial deviation of many-body dynamics from mean-field Gross–Pitaevskii predictions, particularly as the barrier height and interaction strength increase. In addition, the sudden approximation supports localized particle–hole formation after a diabatic quench, and correlation measures unveil a new dynamic regime: the Fock flashlight.