[Paper Review] Faraday patterns in Bose-Einstein condensates. Amplitude equation for rolls in the parametrically driven, damped Gross-Pitaevskii equation

Germán J. de Valcárcel|arXiv (Cornell University)|Apr 18, 2002

Cold Atom Physics and Bose-Einstein Condensates3 citations

TL;DR

This paper derives a complex Ginzburg-Landau-type amplitude equation for roll patterns in parametrically driven, damped Bose-Einstein condensates using the Gross-Pitaevskii equation. It shows that spatially modulated roll states emerge via parametric resonance, and the amplitude dynamics are governed by a Landau equation with broken phase symmetry, revealing the onset and stability of Faraday-like patterns in ultracold quantum fluids.

ABSTRACT

The parametrically driven, damped Gross-Pitaevskii equation, which models Bose-Einstein condensates in which the interatomic s-wave scattering length is modulated in time, is shown to support spatially modulated states in the form of rolls. A Landau equation with broken phase symmetry is derived, which governs the dynamics of the roll amplitude.

Motivation & Objective

• To understand the emergence of spatially modulated patterns, specifically rolls, in Bose-Einstein condensates with periodically modulated s-wave scattering length.
• To analyze the linear stability of the homogeneous BEC ground state under parametric driving to identify conditions for spontaneous symmetry breaking.
• To derive a reduced amplitude equation that governs the slow dynamics of the roll pattern envelope in the weakly nonlinear regime.
• To characterize the role of damping and parametric driving in stabilizing or destabilizing roll structures in dilute quantum gases.

Proposed method

• Uses the parametrically driven, damped Gross-Pitaevskii equation in 1D or 2D geometry to model time-modulated BECs with tunable interactions.
• Performs a linear stability analysis around the homogeneous solution, leading to a damped Mathieu equation for perturbations with wavenumber k.
• Applies multiple-scale asymptotic expansion to separate fast time scales (driving frequency) and slow time scales (envelope evolution).
• Derives the amplitude equation by eliminating secular terms via solvability conditions, resulting in a complex Ginzburg-Landau-type equation for the roll envelope.
• Uses Floquet theory and perturbation methods in the limit of weak damping and near-resonance conditions to obtain analytical expressions.
• Transforms the derived amplitude equation into a standard form—Landau equation with broken phase symmetry—by rescaling time and amplitude variables.

Experimental results

Research questions

• RQ1 Under what conditions does parametric driving of a Bose-Einstein condensate lead to spontaneous formation of roll patterns?
• RQ2 How does damping affect the stability and growth rate of spatially modulated structures in the driven BEC?
• RQ3 What is the effective amplitude equation that governs the slow evolution of roll patterns in the weakly nonlinear regime?
• RQ4 How does the broken phase symmetry in the amplitude equation influence the pattern selection and dynamics of the system?
• RQ5 What is the role of the nonlinear dispersion relation and resonance tongues in determining the onset of pattern formation?

Key findings

• Roll patterns emerge in the parametrically driven, damped Gross-Pitaevskii equation due to parametric resonance, with instability occurring in resonance tongues centered at Ω(k) = nω.
• The amplitude dynamics of the rolls are governed by a Landau equation with broken phase symmetry: dR/dt = −[γ(1+k²)+i(ω−Ω)]R + i(αk²/Ω)R̄ − i(3+5k²)/Ω |R|²R.
• The critical wavenumber for roll formation is k_n = √(√(1+n²ω²)−1), corresponding to the n-th resonance tongue in the parametrically driven system.
• The solvability condition derived from the multiple-scale expansion leads to a consistent amplitude equation that captures the nonlinear saturation of pattern growth.
• The derived amplitude equation is valid in the weak damping and near-resonance regime, and it reproduces the neutral stability curve from linear analysis when R=0.
• The roll solution is expressed in terms of the complex amplitude R(t), with the full wavefunction ψ(x,t) reconstructed via ψ = exp[−i(α/ω)sin(2ωt)](1 + w cos(kx)), where w depends on R and its complex conjugate.