Macroscopic Quantum Self-Trapping (MQST)

• MQST is a nonlinear dynamical regime in macroscopic quantum systems, such as Bose–Einstein condensates, characterized by a persistent population imbalance due to interparticle interactions and coherent tunneling.
• The phenomenon is modeled using two-mode approximations of the Gross–Pitaevskii equation, revealing a bifurcation from Josephson oscillations to self-trapped states when interaction strengths exceed critical thresholds, with quantum fluctuations further modulating the dynamics.
• Experimental implementations in atomic Josephson junctions, photonic molecules, and exciton-polariton systems validate the theoretical predictions, offering actionable insights into the control of nonlinear quantum processes and phase dynamics.

Macroscopic Quantum Self-Trapping (MQST) designates a class of nonlinear dynamical regimes in which a macroscopic quantum many–body system, such as a Bose–Einstein condensate (BEC), becomes dynamically localized—typically as an anomalous population imbalance persisting on long timescales—due to the interplay of interparticle interactions and coherent tunneling. Most commonly manifest in double-well BECs, Josephson junctions, and related multi-mode systems, MQST represents a breakdown of symmetric (Josephson) oscillations and can be generalized to a range of platforms, including exciton-polariton gases, supersolids, and even real-space or momentum-space lattices. Unlike single-particle localization, MQST is the macroscopic consequence of nonlinearity and constructive interference of many-body effects.

1. Theoretical Frameworks and Fundamental Equations

The prototypical scenario for MQST is the bosonic Josephson junction—two weakly coupled BECs in a symmetric double-well potential. The system is accurately described within a two-mode (dimer) approximation, reducing the time-dependent Gross–Pitaevskii equation (GPE) to coupled amplitude equations for the left- and right-localized wavefunctions. Defining macroscopic occupation amplitudes $u(t)$ and $v(t)$, the key dynamical variables are the population imbalance $Z=(N_1-N_2)/N$ and relative phase $\theta=\phi_2-\phi_1$.

Dimensionless two-mode equations incorporating both mean-field and quantum fluctuation corrections (e.g., LHY) take the form

$$\frac{dZ}{d\tau} = -\sqrt{1-Z^2} \sin\theta,\quad \frac{d\theta}{d\tau} = \frac{Z}{\sqrt{1-Z^2}}\cos\theta + \Lambda\,Z - \frac{\varepsilon}{2^{3/2}}\frac{1}{(1+Z)^{3/2}-(1-Z)^{3/2}},$$

where $\Lambda$ quantifies the mean-field nonlinearity, and $\varepsilon$ encodes beyond mean-field (Lee–Huang–Yang) quantum fluctuations (Abdullaev et al., 2023).

These equations are Hamiltonian: $$H(Z, \theta) = -\sqrt{1-Z^2}\cos\theta + \frac{\Lambda}{2}Z^2 + \frac{\varepsilon}{5\sqrt{2}} \left[(1+Z)^{5/2} + (1-Z)^{5/2}\right].$$ Generalizations for multispecies, dipolar, supersolid, or momentum-space systems modify the effective parameters and coupling structure but preserve the core MQST phenomenology (Donelli et al., 28 Jan 2025, Schimelfenig et al., 15 Oct 2025, Yi et al., 2015).

2. Critical Conditions and Dynamical Regimes

A hallmark of MQST is a dynamical bifurcation separating Josephson oscillations from self-trapped states, governed by the initial state's energy relative to the system's separatrix. The criterion for self-trapping is: $H(Z_0, \theta_0) > H_{\text{sep}},$ where $H_{\text{sep}}$ is the separatrix energy (e.g., $v(t)$0 with LHY correction). Solving for the critical parameter yields

$v(t)$1

For a pure contact interaction and $v(t)$2, classical MQST onset is at $v(t)$3; inclusion of quantum fluctuations ($v(t)$4) can lower the threshold or induce MQST even in the absence of a mean-field term (Abdullaev et al., 2023).

MQST is characterized by persistent nonzero time-averaged population imbalance; the phase variable evolves unboundedly ("running phase"). In the phase space $v(t)$5, the trajectory encircles a fixed point at $v(t)$6 rather than $v(t)$7 (Abbarchi et al., 2012, Bardin et al., 26 Feb 2026).

In strongly interacting or more complex systems (supersolids, excitonic condensates, or coupled tubes), MQST can additionally arise from long-range interactions, finite-range effects, or correlation blockade mechanisms, extending the critical condition beyond simple tunneling-versus-interaction arguments (Donelli et al., 28 Jan 2025, Li et al., 2013, Yi et al., 2015).

3. Quantum Fluctuations, Finite-Size Effects, and Exact Quantum Dynamics

While mean-field theory predicts sharp transitions and indefinite self-trapping, exact quantum treatments demonstrate qualitative deviations at finite particle number. For the two-site Bose–Hubbard model, the exact many-body dynamics show that for any finite $v(t)$8, strict MQST is replaced by a quasi-self-trapped regime: the population imbalance eventually crosses zero on a timescale $v(t)$9, where $Z=(N_1-N_2)/N$0 is the smallest energy gap in the spectrum. For interaction strengths above threshold and large $Z=(N_1-N_2)/N$1, $Z=(N_1-N_2)/N$2 scales exponentially with $Z=(N_1-N_2)/N$3, realizing “quasi-MQST” on experimentally relevant timescales (Bardin et al., 26 Feb 2026).

Corrections to mean-field predictions due to number conservation can be computed using atomic coherent state (ACS) variational methods, introducing $Z=(N_1-N_2)/N$4 corrections to the critical interaction and Josephson frequency: $Z=(N_1-N_2)/N$5 with numerically validated accuracy for $Z=(N_1-N_2)/N$6 (Wimberger et al., 2021).

4. Extensions to Multicomponent, Dipolar, and Nonlinear Systems

MQST emerges robustly across platforms once nonlinear interactions induce sufficient bifurcation in the effective Hamiltonian landscape:

• Multicomponent/Binary BECs: Coupled pendulum-type equations describe self-trapping of individual species, with analytic thresholds depending on intra- and inter-species couplings. Out-of-phase and in-phase self-trapping regimes are predicted and confirmed by full coupled GPE numerics (0905.2080).
• Supersolids: Josephson junction physics transposes from spatial weak links to self-induced links created by density modulations. MQST in supersolids follows a generalized two-mode Hamiltonian with modified tunneling and on-site energies, incorporating dipolar, contact, and LHY corrections. The critical imbalance for self-trapping reflects this complex energetic structure (Donelli et al., 28 Jan 2025).
• Momentum-Space Double Wells: For spin-orbit-coupled BECs in optical lattices, MQST manifests as persistent population imbalance between momentum eigenstates, with a dynamical phase transition characterized by diverging oscillation periods and a non-analytic order parameter (Schimelfenig et al., 15 Oct 2025).
• Driven and Periodically Modulated Systems: In dynamically modulated wells, MQST can occur in Floquet-engineered two-mode models, even showing re-entrant tunneling at large nonlinearity, absent in static spatial double wells (Wüster et al., 2011).

5. Experimental Realizations and Detection Protocols

MQST has been realized and measured in diverse settings:

• Atomic Josephson Junctions: MQST is detected via persistent nonzero atom number difference between spatial wells. Feshbach resonances permit tuning of nonlinear interaction strength to access or suppress MQST. Monitoring of phase dynamics and oscillation frequencies allows characterization of the dynamical regime (Nesterenko et al., 2012, Saha et al., 2019).
• Exciton-Polaritons and Photonic Molecules: MQST manifests as localization of condensate emission in one micropillar of a photonic molecule at high density. The finite lifetime of polaritons leads to a dynamical untrapping effect, unique to dissipative condensates (Abbarchi et al., 2012).
• Supersolid and Dipolar Gases: Self-induced weak links in dipolar supersolids host Josephson oscillations and MQST. Protocols involve monitoring of time-dependent density modulations and phase evolution using in-situ imaging and phase-sensitive detection (Donelli et al., 28 Jan 2025).
• Momentum-Space MQST: Imaging of momentum populations following detuning ramps and quench protocols enables mapping of the dynamical phase diagram and the critical MQST threshold, including observation of nonanalytic behavior at the transition (Schimelfenig et al., 15 Oct 2025).

6. Beyond Mean-Field Effects: Dissipation, Quantum Correlations, and Metastable States

MQST regimes are sensitive to quantum and thermal fluctuations, dissipation, and external noise:

• Dissipative Systems: Coupling to baths introduces phase diffusion and reduces the coherence factor $Z=(N_1-N_2)/N$7. The phase-diffusion coefficient and transition from Josephson oscillation to MQST can be quantitatively characterized, with strong dissipation destroying coherent self-trapping (Saha et al., 2022).
• Quantum Correlation Induced Self-Trapping: In low-dimensional systems (e.g., coupled 1D tubes), quantum correlations can block tunneling even without density gradients, a regime not captured by mean-field theory or truncated Wigner approximations. Many-body simulations (t-DMRG) reveal that MQST can arise solely from the energy cost of reshaping correlated states, leading to persistent population imbalance (Li et al., 2013).
• Excitonic Systems with Coupled Order Parameters: Self-trapping can seed local phase transformation, as in excitonic condensates where coupling to charge-transfer and lattice dimerization fields feeds back to stabilize domains beyond the lifetime of the condensate (Yi et al., 2015).

7. Physical Interpretations, Generalizations, and Phenomenological Implications

MQST embodies a classical symmetry-breaking bifurcation—a dynamical analog of spontaneous symmetry breaking in mean-field Hamiltonian systems. The phenomenon's universality extends from ultracold atoms in double wells to rotating toroidal condensates (mapping angular momentum states to “wells”), dynamically modulated potentials, supersolids with self-induced links, and strongly correlated low-dimensional gases (Baharian et al., 2012, Pilipchuk, 2012, Yi et al., 2015). The precise dynamical features—critical threshold, persistent imbalance, running phase, and spectral signatures—are quantitatively accessible, and analytic solutions in terms of elementary or elliptic functions are possible in key limits.

MQST impacts phase-sensitive interferometry, dynamical phase transition studies, band structure engineering, and the fundamental understanding of quantum nonlinear matter-wave dynamics. Mapping between platforms—spatial, momentum, internal state, or correlation-induced—demonstrates the broad utility and reach of MQST as a paradigm of collective nonlinear quantum dynamics in complex systems (Schimelfenig et al., 15 Oct 2025, Yi et al., 2015).