Abstract:This work is concerned with the construction and analysis of structure-preserving Galerkin methods for computing the dynamics of rotating Bose-Einstein condensate (BEC) based on the Gross-Pitaevskii equation with angular momentum rotation. Due to the presence of the rotation term, constructing finite element methods (FEMs) that preserve both mass and energy remains an unresolved issue, particularly in the context of nonconforming FEMs. Furthermore, in comparison to existing works, we provide a comprehensive convergence analysis, offering a thorough demonstration of the methods' optimal and high-order convergence properties. Finally, extensive numerical results are presented to check the theoretical analysis of the structure-preserving numerical method for rotating BEC, and the quantized vortex lattice's behavior is scrutinized through a series of numerical tests.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to construct and analyze structure - preserving finite element methods (FEMs) in the dynamic calculations of rotating Bose - Einstein condensates (BECs), especially based on the Gross - Pitaevskii equation (GPE) with angular momentum rotation terms. Specifically:
1. **Mass and energy conservation**: Due to the existence of the rotation term, constructing a finite element method that can preserve both mass and energy is an unsolved problem, especially in the context of nonconforming FEMs.
2. **Convergence analysis**: Compared with existing works, this paper provides a more comprehensive convergence analysis, demonstrating the optimal and high - order convergence properties of these methods. In particular, the paper solves the problem of unconditional error estimation, that is, without the coupling conditions of time step and spatial grid width, it proves the optimal \(L^2\) and \(H^1\) estimates of the numerical solution, and has a high - precision convergence rate in the \(H^1\) norm.
3. **Numerical verification**: Through a series of numerical experiments, the correctness of the theoretical analysis is verified, and the behavior of the quantized vortex lattice is studied.
### Formula presentation
The form of the Gross - Pitaevskii equation (GPE) is:
\[
i\partial_t u(x,t)=\left[-\frac{1}{2}\nabla^2 + V(x)-\Omega L_z+\beta|u(x,t)|^2\right]u(x,t),\quad x\in U\subset\mathbb{R}^d,\quad t > 0,
\]
where the boundary conditions and initial conditions are respectively:
\[
u(x,t) = 0,\quad x\in\Gamma=\partial U,\quad t\geq0,
\]
\[
u(x,0)=u_0(x),\quad x\in U.
\]
The system preserves the total mass and total energy conservation:
\[
N(u(\cdot,t)):=\int_U|u(x,t)|^2dx\equiv N(u(\cdot,0)) = N(u_0),\quad t\geq0,
\]
\[
E(u(\cdot,t)):=\int_U\left[\frac{1}{2}|\nabla u(x,t)|^2+V(x)|u(x,t)|^2+\frac{\beta}{2}|u(x,t)|^4-\Omega\overline{u}(x,t)L_z u(x,t)\right]dx\equiv E(u(\cdot,0)) = E(u_0),\quad t\geq0.
\]
### Main contributions
1. **Innovative stabilization term**: A new stabilization term is introduced, enabling nonconforming FEMs to preserve conservation in the sense of discrete mass and energy without affecting the convergence rate.
2. **Unconditional error estimation**: Unconditional error estimation is established, applicable to any time step and spatial grid width, including the optimal \(L^2\) and \(H^1\) estimates and a high - precision \(H^1\) convergence rate.
3. **Handling of attractive interactions**: It can handle attractive interactions (i.e., the case where \(\beta < 0\)), which is less involved in the existing literature.
4. **Handling of complex models**: For the first time, a high - order spatial convergence rate for the Gross - Pitaevskii equation containing angular momentum rotation terms is achieved.
In summary, this paper solves the key problems in the dynamic calculations of rotating Bose - Einstein condensates by constructing and analyzing structure - preserving finite element methods, and provides strong theoretical and numerical support.