Can a Uniform Imaginary Magnetic Field Be Generated for Ultracold Atoms?

A four-level scheme that separates loss from field

The paper asks whether a spatially uniform imaginary magnetic field can be engineered for ultracold atoms, and answers by construction rather than by measurement. The proposal starts from two Raman-coupled internal states, |↑⟩ and |↓⟩, and adds two auxiliary states |a↑⟩ and |a↓⟩. Each Raman-coupled state is coupled to its own auxiliary state with Rabi frequency ωλ, and each auxiliary state decays irreversibly at rate κλ. Because population does not return to the four-level system, the surviving ensemble is described by non-Hermitian evolution.

Eliminating the auxiliary manifold under the conditions κλ ≫ ΩR, |ωλ|, and ℏκλ much larger than the relevant motional energies gives an effective two-level Hamiltonian with a purely imaginary detuning ±iΔ(x)/2, a state-independent loss rate Γ(x), and the Raman coupling ΩR(x)e^{±i2k_R y}. The two loss channels are what make Δ and Γ separately controllable: Γ(x) = |ω↑(x)|²/(2κ↑) + |ω↓(x)|²/(2κ↓) and Δ(x) = |ω↓(x)|²/(2κ↓) − |ω↑(x)|²/(2κ↑).

This separation matters because the paper contrasts it with the two-level situation. With a linear detuning at fixed Raman amplitude, the vector potential is governed by Δ/√(ΩR² − Δ²), so the field is nonuniform and the gap closes at |Δ| = ΩR. Removing atoms from only one of the two Raman-coupled states, as in the earlier experiment cited by the authors, also introduces a scalar loss gradient that affects motion independently of the magnetic field. The four-level design is intended to remove both problems at once.

Uniform field requires a nonzero curl

A uniform imaginary vector potential is not the same object as a uniform imaginary magnetic field. The paper notes that a spatially uniform imaginary vector potential has zero curl, and an imaginary magnetic field instead requires a spatially varying imaginary vector potential with a nonzero curl. The synthetic field is defined as B_s = ∇ × A_{s,s}, with A_{s,s} = iℏ⟨χ_s^L|∇χ_s^R⟩ built from the left and right dressed states of the non-Hermitian atom-light coupling.

Taking the curl gives B_s = is(ℏ/2)∇(sinh ξ) × ∇φ, where tanh ξ = Δ/ΩR and φ is the Raman phase. The field is purely imaginary and has opposite signs in the two dressed branches. It requires nonparallel gradients of the normalized detuning and the Raman phase, so the geometry of the light fields, not just their strength, sets the field.

For the specific configuration φ = −2k_R y and ξ = ξ(x), the field reduces to B_s = −isℏk_R (d/dx)[Δ/√(ΩR² − Δ²)] e_z. A uniform field therefore requires sinh ξ to be affine in x. The authors choose sinh ξ = ηx, Δ(x) = Ω0 ηx, ΩR(x) = Ω0√(1 + η²x²), and Γ(x) = Γ0, which keeps the dressed-state splitting constant at ε_s = sℏΩ0/2 and yields B_s = −isB e_z with B = ℏk_R η. The constant Raman amplitude alternative instead gives B_s = −isℏk_R η/(1 − η²x²)^{3/2} e_z for |ηx| < 1, diverging formally at the closing gap where the single-band adiabatic description ceases to apply.

An inverted oscillator solves the quadratic model

After linearizing the x component of the vector potential to A_x = −iℏη²x/2, compensating the geometric scalar potential with V = −W, and dropping the constant real term in A_y and the constant dressed-state energy, the lower-branch Hamiltonian becomes H = [(p_x − A_x)² + (p_y − iBx)²]/(2m) − iℏΓ0/2. This is quadratic, so the authors can solve it exactly rather than numerically.

The nonunitary substitution ψ(r,t) = e^{−Γ0 t/2} e^{η²x²/4} χ(r,t) turns the problem into an inverted harmonic oscillator with H̃ = (p_x² + p_y²)/(2m) − mω²x²/2 − iωx p_y, where ω = B/m = ℏk_R η/m. Because p_y commutes with H̃, each momentum sector evolves independently, and completing the square produces an imaginary position shift generated by T_p = exp(−p p_x/(ℏmω)). The propagator is the inverted-oscillator kernel, and for the initial Gaussian envelope both integrals remain Gaussian.

The result is a closed set of first moments, a surviving fraction N(t), and a position covariance. The authors report that the terms proportional to σ_x² in the first moments describe initial position uncertainty, while terms proportional to 1/σ_x² describe initial momentum uncertainty, with Var(p_x)_0 = ℏ²/(4σ_x²). A constant loss term −ℏΓ0/2 cancels between the norm-derivative and normalization terms in the non-Hermitian Ehrenfest equation, so it does not affect normalized observables. The derivation of that cancellation requires no Gaussian assumption.

Transport depends on the initial width

The paper's headline dynamical result is that transport in this synthetic field depends on the width of the atomic wavepacket, which the authors present as unique to non-Hermitian systems. For a packet launched along the Raman recoil direction, attenuation produces a transverse density shift that is linear in time, preceding the quadratic onset of transverse momentum. The covariance depends on the initial spreads independently of the first moments.

The inverted oscillator generates quadrature squeezing, while the imaginary-field coupling produces spatial correlations and a rotation of the principal axes of the Gaussian density. The paper derives a threshold: above it, one principal position spread initially narrows while the other broadens. For the circular launch with σ_x = σ_y = σ, writing r = ℏ/(2mωσ²) and ε = η²σ², the short-time covariance expansion has a matrix whose smaller eigenvalue is negative precisely when 4(1 + ε)r² < 1. At equality, λ_−(t) = σ² and λ_+(t) = σ²[1 + ((3 + ε)²/4 + r²)s²/D_η], with θ_+(t) = −arctan(2r/(3 + ε)).

This threshold is the paper's most decision-relevant number for a reader comparing scenarios, because it says the qualitative behavior of the packet is set by the ratio of the initial width to the magnetic length scale, not by the field strength alone. Reversing the field reverses the shear and reflects the principal axes while leaving their widths unchanged. As η → 0 at fixed k_R, both vector-potential components vanish and the covariance reduces to free expansion, v_ii = σ_i² + ℏ²t²/(4m²σ_i²) and v_xy = 0.

Cyclic transport yields an area-dependent population contrast

The paper also treats transport in a moving trap around a closed circuit. For a closed path, ∮ X dX = 0, so the A_x component gives no net geometric contribution. The spatial mode and its normalization return to their initial values, leaving a geometric factor W_C = exp((i/ℏ)∮_C A·dr) = e^{−BA/ℏ}, where the area is defined by the trap path because the density can be displaced from the well.

The surviving populations for the two circulation senses are N_±(T) ≃ e^{−Γ0 T} e^{∓2πBR_xR_y/ℏ}, so ln[N_+(T)/N_−(T)] ≃ −4πBR_xR_y/ℏ. Varying R_y at fixed R_x therefore changes the logarithmic contrast linearly, and adiabatic paths with equal signed area give the same factor. Because the atoms sample a nonzero field, the authors describe this response as a geometric Aharonov–Bohm analogue.

Reversing the transport direction changes the sign of the geometric contribution to the population, and the contrast between clockwise and counterclockwise circuits increases with the enclosed area. This is a predicted signature in the surviving population, not a measured one, and it is the quantity the authors propose as an observable target alongside the density and principal spreads.

Adiabatic projection agrees with two-level numerics

The adiabatic projection is an approximation, and the paper tests it. The authors compare the two-level non-Hermitian Hamiltonian with the adiabatic projection description and report numerical agreement, which they take as confirmation that the adiabatic approach is applicable. The comparison also identifies the interval in which the adiabatic description remains accurate and shows its dependence on both the dressed-state gap and the initial preparation.

The approximation itself is the neglect of transitions between dressed branches, not an expansion in Δ/ΩR. The paper states that the expressions for the geometric vector and scalar potentials involve no expansion in Δ/ΩR. The conditions for eliminating the auxiliary manifold are κλ ≫ ΩR, |ωλ|, and ℏκλ much larger than the relevant motional energies, and the resolvent is approximated as R(E) ≈ −(2i/ℏ)Σλ κλ^{−1}|aλ⟩⟨aλ|.

Two further conditions bound the quadratic model. Nonnegative loss rates require Γ0 ≥ Ω0|η|L over |x| ≤ L, where 2L is the system size, and the linearization A_x ≃ −iℏη²x/2 assumes |ηx| ≪ 1. The uniform-field configuration also requires compensating the geometric scalar potential with V = −W in the populated region. Outside these conditions, the constant-Raman-amplitude alternative shows the field becoming nonuniform and formally diverging at the gap closing, where the single-band adiabatic description breaks down.

The width threshold is the number to watch

The paper is a theoretical proposal, so its results are predictions and analytical derivations rather than measurements. The authors state that spatial beam shaping and state-dependent removal offer a route to observing these effects, while interactions and controlled preparation of moving dressed states remain directions for further study. The numerical data underlying the figures are openly available in Zenodo.

For a reader deciding what to take from this work, the priority is the width threshold 4(1 + ε)r² < 1, with r = ℏ/(2mωσ²) and ε = η²σ², because it separates two qualitatively different transport regimes: below it, one principal position spread initially narrows while the other broadens; at equality, one spread stays at σ². The second priority is the comparison between the four-level design and the two-level alternative. The four-level scheme keeps the dressed-state gap constant and the overall loss coordinate-independent, which is what makes the field uniform; the two-level scheme with linear detuning at fixed Raman amplitude gives a nonuniform field and a gap closing at |Δ| = ΩR. The third is the area-dependent population contrast ln[N_+(T)/N_−(T)] ≃ −4πBR_xR_y/ℏ, which scales linearly with R_y at fixed R_x and reverses sign with the transport direction. These three quantities, not the existence of the scheme alone, are what a follow-up experiment or a reader comparing non-Hermitian transport scenarios would need to check.

Sources used in this answer