Fast momentum-selective transport of Bose-Einstein condensates via controlled non-adiabatic dynamics in optical lattices

Fast momentum-selective transport of Bose–Einstein condensates via controlled non-adiabatic dynamics in optical lattices investigates a protocol for achieving narrow momentum distributions in ultracold gases using rapid, non-adiabatic manipulation. The key finding is that spectral purity can be maintained even with fast loading and release by synchronizing the ramp durations with internal condensate breathing oscillations, identifying magic times where quasi-monochromatic momentum distributions are achieved.

The gist: Highly monochromatic momentum distributions can be achieved when the loading/release or acceleration duration matches specific values, referred to as magic times, which synchronize with internal breathing oscillations.

Theoretical Framework and Numerical Methods

The study employs the time-dependent Gross–Pitaevskii equation (GPE) to simulate the full transport sequence of a dilute Bose–Einstein condensate of 87Rb atoms in a one-dimensional optical lattice. The dynamics are governed by the GPE, which accounts for both kinetic energy, external potential, and mean-field interactions. The system is treated under tight transverse confinement where the longitudinal dynamics can be accurately described by the effective one-dimensional GPE (Equation 1). Numerical solutions utilize a Fourier pseudo-spectral method for spatial discretization and a second-order split-operator scheme for time evolution. The initial state is obtained via imaginary time propagation, and the lattice depth is ramped up while the harmonic trap is switched off over specified ramps.

Transport Protocol Stages

The complete transport protocol consists of three distinct stages:

1. Loading into the Optical Lattice: This stage involves ramping on a periodic potential while simultaneously switching off the harmonic trap using complementary ramp functions, defined by Equation (2) and (3). The duration of this ramp, denoted as tL, determines the degree of adiabaticity.

2. Lattice Acceleration: Once loaded at time t = tL, the condensate is accelerated using a symmetric trapezoidal acceleration profile. This profile involves a linear increase to a maximum acceleration (aOL), holding it constant for a plateau duration (∆), and then linearly decreasing back to zero over another duration δ. The total acceleration duration is tacc = 2δ + ∆.

3. Release from the Lattice: The final stage involves ramping down the optical lattice over a duration tR = tL, using the same envelope function as in the loading phase, resulting in Equation (8).

Mechanism for Momentum Selectivity

The paper identifies magic times as specific durations—for loading/release or acceleration—that are synchronized with the internal breathing oscillation period of the condensate. This synchronization is a dominant mechanism governing spectral purity under fast loading conditions. The intra-site breathing dynamics are shown to be directly correlated with the final momentum spread; specifically, tracking the condensate’s spatial width demonstrates this correlation. A variational model based on a Gaussian ansatz quantitatively reproduces these dynamics and provides physical insight into the breathing mechanism, where the size oscillations are linked to spectral purity.

Analysis of Momentum Distributions

The final momentum distribution P(k) is computed as the squared modulus of the Fourier transform of the wave function at the end of release. The results show that for short loading durations (tL < 0.3 ms), non-adiabatic excitations manifest as pronounced side peaks at k = 188kL and k = 192kL (i.e., ±2kL from the main peak). As tL increases, these sidebands are progressively suppressed, and the distribution narrows around the target momentum class of k = 190kL. The transition toward adiabatic transport is clearly established for loading durations exceeding 0.4 ms, corresponding to approximately 60 times the characteristic breathing period (π/ωOL). This confirms that sufficient ramp time ensures minimal excitation and optimal spectral purity, with sideband populations suppressed below 1% in the adiabatic limit.

Role of Interactions and Speedup

The analysis incorporates mean-field interactions via the GPE, which are shown to modulate the breathing dynamics. Comparing the full Gross–Pitaevskii equation with a linear Schrödinger equation (setting g1D = 0) reveals that repulsive interactions can modify non-adiabatic excitations; in this specific case, they lead to a final spatial compression and thus a broader momentum distribution and degraded spectral purity. The protocol offers significant speedup factors of 3 to 6 compared to adiabatic protocols while maintaining high transfer fidelities, providing a practical route for coherent transport under stringent timing constraints. The magic-time phenomenon is robust with respect to the specific shape of the acceleration ramp, confirming its reliance on the total acceleration duration rather than the profile itself. Finally, magic loading times are also identified as optimal durations where P0 approaches unity despite non-adiabatic ramps.

Diagnostic Observables

The intra-site spatial width ∆x(t) is used as a sensitive realspace indicator of spectral purity under non-adiabatic conditions.

Based on the provided scientific paper, here are specific ways an AI system could be improved and what those improvements would enable:

) The core improvement lies in developing a Breathing-Mode Predictive Transport Model (BMPTM). This model integrates non-adiabatic dynamics, intra-site breathing modes, and spectral purity into a single framework.

) This BMPTM can perform the following specific tasks:

1. It can predict the optimal loading/acceleration/release durations (magic times) required to achieve near-perfect momentum selectivity (e.g., >98% population in a single momentum class).

2. It can dynamically determine the necessary acceleration profiles for a given target momentum transfer, accounting for lattice depth and atomic interactions.

3. It can act as a spectral purity diagnostic by calculating the expected final momentum distribution based on real-space observables (like intra-site width, ∆x(t)) during the transport sequence.

) The improved AI system can enable the following capabilities:

1. Predicting Optimal Quantum Control Sequences for Quantum Sensors:

2. Designing High-Fidelity Matter-Wave Sources for Interferometry:

3. Optimizing Coherent Transport Protocols in Time-Constrained Environments:

4. Characterizing Non-Equilibrium Dynamics in Cold Atom Systems:

) Specific Enhancements to AI Capabilities (What the System Can Do):

) Detailed Functionality Enabled by the Improved AI System (Specific Actions):

1. When presented with target momentum transfer requirements and system parameters (lattice depth, atom number), the AI will output a precise sequence of loading time and acceleration duration (magic times) required to ensure spectral purity exceeding a user-defined threshold (e.g., 98% population in the target state).

2. The AI can simulate and suggest specific trapezoidal acceleration profiles that minimize unwanted sideband excitations (at ±2hk¯ L) based on real-time monitoring of the predicted intra-site breathing dynamics, effectively tuning the acceleration ramp in real-time during simulation or experimental execution.

3. The system will use measured or simulated real-space data (like the time evolution of ∆x(t)) to diagnose whether a transport protocol is achieving its intended momentum selectivity, flagging deviations that indicate non-adiabatic excitation or interaction effects (e.g., distinguishing between the effect of repulsive vs. non-interacting atoms).

4. The AI can rapidly estimate the required total sequence duration and speedup factor by comparing the predicted performance against known adiabatic protocols, providing an immediate quantitative assessment of protocol efficiency for time-sensitive quantum sensing applications.

We present a detailed numerical study of a protocol for momentum-selective transport of a Bose-Einstein condensate (BEC) in a one-dimensional optical lattice, achieving narrow momentum distributions through controlled non-adiabatic dynamics. The protocol consists of non-adiabatic loading into the lattice, coherent acceleration using a symmetric trapezoidal acceleration profile, and non-adiabatic release into free space. Using the time-dependent Gross-Pitaevskii equation, we simulate the full sequence and analyze the role of non-adiabatic excitations on the final momentum distribution. We identify the intra-site breathing dynamics as the dominant mechanism governing spectral purity under fast loading conditions. By tracking the condensate's spatial width during the evolution, we demonstrate a direct correlation with the final momentum spread. A variational model based on a Gaussian ansatz quantitatively reproduces the observed dynamics and provides physical insight into the breathing mechanism. Our results reveal the existence of "magic" times, i.e., specific loading or acceleration durations synchronized with the breathing oscillation period, where quasi-monochromatic momentum distributions can be achieved even with loading times as short as 100 microseconds. In the tight-binding regime, this approach offers speedup factors of 3 to 6 compared to adiabatic protocols while maintaining high transfer fidelities, providing a practical route to coherent transport for quantum sensors operating under stringent timing constraints.