Splitting methods for the Gross-Pitaevskii equation on the full space and vortex nucleation

This paper establishes the convergence of Lie-Trotter and Strang splitting schemes for the Gross-Pitaevskii equation with time-dependent potentials and non-zero boundary conditions in Zhidkov spaces, while demonstrating the conservation of generalized mass, near-preservation of energy, and investigating quantum vortex nucleation in relevant experimental settings.

The Gist

Imagine you are trying to predict the weather, but instead of clouds and rain, you are tracking a super-cold, super-fluid substance called a Bose-Einstein Condensate. This is a state of matter where atoms act like a single giant wave. The math that describes how this "quantum fluid" moves is called the Gross-Pitaevskii Equation (GP).

The problem? This equation is incredibly complex. It's like trying to solve a puzzle where the pieces are constantly changing shape, and you have to do it over an infinite field (the "full space") where the fluid never truly stops moving at the edges.

This paper by Quentin Chauleur and Gaspard Kemlin is essentially a guidebook for building a better, more reliable calculator to simulate this fluid. Here is the breakdown of their work using simple analogies:

1. The Problem: The "Two-Step Dance"

The equation governing this fluid has two main forces acting on it at the same time:

• The Wave Force (Diffusion): The fluid wants to spread out and ripple like water in a pond.
• The Push Force (Nonlinearity & Potential): The fluid pushes against itself and reacts to obstacles (like a spoon stirring a cup of coffee).

Solving both forces at the exact same time is mathematically impossible for a computer to do perfectly. So, scientists use a trick called Splitting Methods.

The Analogy: Imagine you are trying to walk a dog that wants to run in a circle (the wave) while you are trying to pull it toward a park bench (the obstacle).

• The Lie-Trotter Method: You take one step forward (pulling), then let the dog run in a circle for a moment, then pull again. It's a bit jerky, like walking with a lurch.
• The Strang Method: You take half a step, let the dog run, take a full step, let the dog run, and take the final half-step. This is smoother and much more accurate.

The authors proved mathematically that these "splitting" steps work perfectly for this specific quantum fluid, even when the fluid is moving toward infinity and hitting obstacles.

2. The "Zhidkov Space": The Infinite Hotel

Usually, math problems are solved in a finite room. But this fluid exists on an infinite plane. To handle this, the authors used a special mathematical "hotel" called Zhidkov spaces.

The Analogy: Imagine a hotel that stretches forever. Most guests (math functions) stay in the lobby (finite space). But this fluid has guests who live in the infinite hallways. The authors proved that their splitting methods work even for these "infinite hallway" guests, ensuring the simulation doesn't break down when the fluid gets far away from the center.

3. The "Vortex Nucleation": The Quantum Whirlpool

The most exciting part of the paper is what happens when you stir this fluid. If you spin a spoon in water, you get a whirlpool. If you spin a "stirring potential" (an invisible force field) in this quantum fluid, you create Quantum Vortices.

The Analogy: Think of the fluid as a perfectly smooth sheet of ice. If you drag a stick across it, it usually just slides. But if you drag it fast enough or spin it, the ice suddenly cracks and forms a tiny, perfect tornado (a vortex).

• The authors used their new "calculator" to simulate dragging an invisible obstacle through this fluid.
• They successfully watched vortices pop into existence (nucleate) out of nowhere, just like in real physics experiments. This proves their method is accurate enough to capture these tiny, dramatic events.

4. Conservation Laws: The "Unbreakable Rules"

In physics, certain things must always stay the same, like the total amount of "stuff" (Mass) or the total energy.

• Mass: The number of particles in the fluid shouldn't magically appear or disappear.
• Energy: The total energy should stay constant (unless you are adding external heat).

The authors showed that their splitting methods are "honest." They don't cheat the math. Even though the computer approximates the steps, it respects these unbreakable rules. The "Mass" stays perfectly preserved, and the "Energy" stays almost perfectly preserved (with only tiny, acceptable errors that get smaller as the simulation gets more detailed).

5. The "Dark Soliton": The Perfect Test

To prove their calculator works, they tested it on a Dark Soliton.
The Analogy: Imagine a wave on a string. A "Dark Soliton" is a specific, perfect "hole" or dip in the wave that travels without changing shape. It's like a perfect, solitary shadow moving across a bright screen.

• They ran their simulation with this perfect shadow.
• Result: The shadow stayed perfect. The "Lie" method was okay (first-order accurate), but the "Strang" method was excellent (second-order accurate), meaning the error was tiny and disappeared quickly as they made the steps smaller.

Summary

This paper is a mathematical guarantee. It tells physicists and engineers: "You can trust this specific way of calculating the Gross-Pitaevskii equation. It works on infinite fields, it respects the laws of physics (mass and energy), and it is accurate enough to simulate the creation of quantum whirlpools."

They didn't just guess; they built a rigorous mathematical bridge between the messy reality of quantum fluids and the clean, step-by-step logic of computer simulations.

Technical Summary

Here is a detailed technical summary of the paper "Splitting Methods for the Gross-Pitaevskii Equation on the Full Space and Vortex Nucleation" by Quentin Chauleur and Gaspard Kemlin.

1. Problem Statement

The paper addresses the numerical time integration of the Gross-Pitaevskii (GP) equation on the full space $\mathbb{R}^d$ ($d \in \{1, 2, 3\}$) with non-vanishing boundary conditions at infinity. The equation is given by:
$$i\partial_t u = \Delta u + \frac{1}{\varepsilon^2}(1 - |u|^2)u + V(t, x)u$$
subject to $|u(t, x)| \to 1$ as $|x| \to \infty$.

Key Challenges:

• Boundary Conditions: Unlike standard nonlinear Schrödinger equations where $u \to 0$ at infinity, the GP equation models Bose-Einstein condensates where the density approaches a non-zero constant. This requires working in Zhidkov spaces ($X^k$) rather than standard Sobolev spaces $H^k$, as the solution $u$ itself is not in $L^2$, though its gradient and deviation from 1 are.
• Time-Dependent Potentials: The potential $V(t, x)$ acts as a "stirring" mechanism (e.g., moving obstacles), which is crucial for studying vortex nucleation but complicates the analysis of conservation laws and splitting schemes.
• Lack of Convergence Theory: While splitting methods (Lie-Trotter and Strang) are widely used numerically for GP, there was a lack of rigorous theoretical convergence proofs in the functional framework appropriate for non-vanishing boundary conditions, especially with time-dependent potentials.

2. Methodology

The authors develop a rigorous convergence analysis using operator splitting methods (Lie-Trotter and Strang) adapted to the specific geometry of the GP equation.

A. Functional Framework and Decomposition

Instead of working directly with $u$, the authors decompose the solution as $u = \phi + v$, where:

• $\phi$ is a fixed, finite-energy stationary solution (or traveling wave) satisfying the boundary conditions.
• $v \in H^k(\mathbb{R}^d)$ represents the perturbation, which decays at infinity.

This transforms the problem into an equation for $v$:
$$i\partial_t v = -i\Delta v - i\Delta \phi - i(1 - |\phi+v|^2 + V)(\phi+v)$$
This allows the use of standard Sobolev spaces $H^k$ for the variable $v$, while handling the affine nature of the linear flow.

B. Splitting Schemes

The equation is split into two sub-problems:

1. Affine Linear Flow ($A$): $i\partial_t w = \Delta w + \Delta \phi$. The flow is explicitly solvable: $\Phi^t_A(w_0) = e^{-it\Delta}w_0 + e^{-it\Delta}\phi - \phi$.
2. Nonlinear Flow ($B$): $i\partial_t z = (1 - |\phi+z|^2 + V)\phi + z$. Due to the specific structure, this flow preserves $|\phi+z|^2$ and is explicitly solvable.

The authors define:

• Lie-Trotter: $u^{n+1} = \Phi^\tau_B \circ \Phi^\tau_A (u^n)$
• Strang: $u^{n+1} = \Phi^{\tau/2}_A \circ \Phi^\tau_B \circ \Phi^{\tau/2}_A (u^n)$

C. Analytical Tools

• Lie Derivatives Framework: The authors adapt the abstract Lie derivative framework (typically used for autonomous systems) to non-autonomous, infinite-dimensional affine systems. They treat the time-dependent potential by extending the phase space to include time as a variable ($\tilde{v} = (v, t)$), converting the non-autonomous system into an autonomous one.
• Commutator Estimates: Convergence rates are derived by analyzing the commutators $[A, B]$ and $[A, [A, B]]$ in the context of the specific nonlinearities and the affine term $\Delta \phi$.
• Stability Analysis: Rigorous stability estimates are proven for the flows in $H^k$ spaces, utilizing the algebra properties of Sobolev spaces in dimensions $d \le 3$.

3. Key Contributions

1. First Rigorous Convergence Proofs: The paper provides the first theoretical convergence estimates for Lie-Trotter (first-order) and Strang (second-order) splitting schemes for the GP equation on the full space with non-vanishing boundary conditions and time-dependent potentials.
2. Adaptation to Affine Flows: The authors successfully extend the standard splitting analysis (which assumes linear flows) to affine flows arising from the decomposition $u = \phi + v$. This required careful handling of the non-homogeneous term $\Delta \phi$.
3. Non-Autonomous Extension: The framework is generalized to handle time-dependent potentials $V(t, x)$ by embedding the system into an extended phase space, proving that the standard convergence rates hold even when the potential varies in time.
4. Conservation Properties:
   • Generalized Mass: Proven to be exactly preserved by the splitting schemes.
   • Energy: Proven to be "near-conserved" (quasi-preserved) with an error scaling consistent with the order of the scheme ($O(\tau)$ for Lie, $O(\tau^2)$ for Strang) when $\partial_t V = 0$.
5. Vortex Nucleation Analysis: The paper investigates the physical phenomenon of quantum vortex nucleation in 2D, simulating condensates flowing past moving and rotating obstacles.

4. Main Results

Theoretical Convergence

Let $u(t)$ be the exact solution and $u^n$ the numerical approximation at time $t_n = n\tau$. Under appropriate regularity assumptions on the initial data and potential:

• Lie-Trotter Scheme:
  $$\| u(t_n) - u^n_L \|_{X^2} \le C \tau$$
• Strang Scheme:
  $$\| u(t_n) - u^n_S \|_{X^2} \le C \tau^2$$
  The constants $C$ depend on the dimension, time horizon, and norms of the initial data, potential, and the background solution $\phi$.

Conservation Laws

• Mass: The generalized mass $M(u) = \lim_{R\to\infty} \int (1-|u|^2)\eta_R$ is exactly preserved by both schemes.
• Energy: For time-independent potentials, the Ginzburg-Landau energy is nearly preserved. The error in energy scales as $O(\tau)$ for Lie and $O(\tau^2)$ for Strang, matching the global error order.

Numerical Validation

• 1D Dark Soliton: Simulations confirm the theoretical convergence rates. Interestingly, for the specific case of a dark soliton initial condition, the energy error exhibits super-convergence (higher than expected order) due to symmetry cancellations, though this disappears with generic perturbations.
• 2D Vortex Nucleation:
  • Case 1 (Moving Obstacle): A Gaussian obstacle moving linearly through the condensate generates vortex-antivortex pairs in a conical wake.
  • Case 2 (Rotating Obstacle): A rotating Gaussian obstacle leads to periodic nucleation of vortex pairs, both inside and outside the rotation radius, mimicking experimental setups for quantum turbulence.

5. Significance

This work bridges a critical gap between the numerical practice and theoretical analysis of the Gross-Pitaevskii equation.

• Mathematical Rigor: It validates the use of efficient splitting methods for problems where standard $L^2$-based analysis fails due to non-zero boundary conditions.
• Physical Relevance: By proving convergence for time-dependent potentials, the paper provides a solid theoretical foundation for simulating quantum turbulence and vortex dynamics in Bose-Einstein condensates, which are central to modern condensed matter physics.
• Methodological Innovation: The technique of transforming non-autonomous affine systems into autonomous ones for error analysis offers a template for analyzing other splitting schemes in complex physical settings.