Strong interaction regime of the nonlinear Landau-Zener problem for photo- and magneto-association of cold atoms

We discuss the strong interaction regime of the nonlinear Landau-Zener problem coming up at coherent photo- and magneto-association of ultracold atoms. We apply a variational approach to an exact third-order nonlinear differential equation for the molecular state probability and construct an accurate approximation describing the whole time dynamics of the coupled atom-molecular system. The resultant solution improves the accuracy of the previous approximation by A. Ishkhanyan et al. [J. Phys. A 39, 14887 (2006)]. The obtained results reveal a remarkable observation that in the strong coupling limit the resonance crossing is mostly governed by the nonlinearity while the coherent atom-molecular oscillations coming up soon after the resonance has been crossed are principally of linear nature. This observation is supposed to be general for all the nonlinear quantum systems having the same generic quadratic nonlinearity, due to the basic attributes of the resonance crossing processes in such systems. The constructed approximation turns out to have a larger applicability range (than it was initially expected) covering the whole moderate coupling regime for which the proposed solution accurately describes all the main characteristics of the system’s evolution except the amplitude of the coherent atom-molecule oscillation, which is rather overestimated.

💡 Research Summary

The paper addresses the nonlinear Landau‑Zener (LZ) problem that arises in coherent photo‑association and magnetic Feshbach resonance of ultracold atoms, focusing on the strong‑interaction (large‑λ) regime. Starting from the standard two‑mode mean‑field model, the authors rewrite the coupled complex amplitude equations for the atomic (a₁) and molecular (a₂) fields as a single third‑order nonlinear differential equation for the molecular population probability p(t)=|a₂|². The dimensionless Landau‑Zener parameter λ=U₀²/δ₀ measures the strength of the coupling relative to the sweep rate of the detuning.

In the strong‑coupling limit (λ≫1) the last two terms of the exact equation dominate, suggesting a simplified “limit” equation. The authors improve upon earlier work by adding a small constant A to this limit equation, thereby retaining a term that was previously discarded. By applying a variational approach and integrating the modified equation via a change of variables, they obtain an implicit quartic algebraic relation for p₀(t) (the zero‑order approximation). The integration constant C is fixed by the initial condition p₀(−∞)=0. To determine A, they substitute p₀(t) back into the exact equation, compute the residual, and enforce the cancellation of the divergence at the resonance crossing (t=0). Using a single Newton iteration they find A≈(9/4)^{1/3}·λ^{−1}, which yields a highly accurate description of the dynamics up to the moment when the system passes through resonance.

However, p₀(t) alone fails to capture the coherent atom‑molecule oscillations that appear after the crossing. To remedy this, the authors introduce a correction term u(t) defined by p(t)=p₀(t)+u(t). The exact equation for u(t) is still nonlinear, but assuming u≪p₀ allows them to linearize it, obtaining a driven linear LZ equation with an effective Landau‑Zener parameter λ* and a scaling factor C. They conjecture that u(t) can be represented as C·p_LZ(t;λ*), where p_LZ is the known solution of the linear LZ problem expressed through Kummer confluent hypergeometric functions.

The parameters λ* and C are fixed by minimizing the residual of the full equation after substitution of the ansatz. This leads to the analytic expressions λ*≈−(9/2)·λ and C≈(6/λ)·A. Consequently, the final approximate solution reads

 p(t) ≈ p₀(t) + C·p_LZ(t; λ*),

which combines the nonlinear “crossing” part p₀(t) with a linear‑oscillatory part. Extensive numerical simulations show that this composite solution reproduces the exact dynamics with remarkable accuracy across a wide range of λ, from the strong‑coupling limit down to moderate values (λ≈3–10). The final transition probability follows a power‑law dependence on λ, in contrast to the exponential dependence of the linear LZ formula, confirming the dominant role of nonlinearity during the crossing. The post‑crossing oscillations are essentially linear, as reflected by the effective λ* and the accurate phase and period of the oscillations. The only noticeable discrepancy is a modest overestimation of the oscillation amplitude, which the authors acknowledge.

Physically, the work reveals a universal two‑stage picture for systems with quadratic nonlinearity: (i) the resonance crossing is governed by the nonlinear term, leading to a non‑exponential scaling of the transition probability; (ii) after the crossing, the system behaves like a linear LZ problem with renormalized parameters, producing coherent atom‑molecule oscillations. This separation is argued to be generic for a broad class of bosonic field theories featuring the same quadratic nonlinearity.

In summary, the authors provide a refined analytical framework that improves upon previous approximations, extends the applicability to moderate coupling, and offers clear physical insight into the interplay between nonlinearity and linear dynamics in ultracold atom‑molecule conversion. The methodology—variational treatment of the exact nonlinear equation plus a linear correction—could be adapted to more complex scenarios involving losses, multi‑mode couplings, or external noise, opening avenues for future theoretical and experimental investigations.