The Coherence Field Interaction Chamber

ADI Crank–Nicolson Propagator (all scenarios)

All five simulations use the same 2D ADI Crank–Nicolson scheme with split \(x/y\) Laplacian sweeps, second-order accuracy in both space and time, and Dirichlet boundary conditions. The nonlinear GPE coupling \(g|\psi|^2\) is handled via a frozen-density phase-rotation half-step applied before and after the linear ADI sweeps.

Absorber: cosine-taper imaginary potential \(\Gamma(r) = \Gamma_0\bigl[(r/L - f_0)/(1-f_0)\bigr]^4\) for \(r > f_0 L\), with \(f_0 = 0.65\), \(\Gamma_0 = 10\) (IX1–2, IX4–5) or \(f_0 = 0.70\), \(\Gamma_0 = 8\) (IX3). The wave function is renormalised after each step.

Reproduce all animations:

source .venv/bin/activate
python3 chapters/gen_interaction_animations.py          # all five
python3 chapters/gen_interaction_animations.py IX1 IX3  # subset
python3 chapters/gen_interaction_animations.py --list   # show parameters