Abstract:This work considers the numerical computation of ground states of rotating Bose-Einstein condensates (BECs) which can exhibit a multiscale lattice of quantized vortices. This problem involves the minimization of an energy functional on a Riemannian manifold. For this we apply the framework of nonlinear conjugate gradient methods in combination with the paradigm of Sobolev gradients to investigate different metrics. Here we build on previous work that proposed to enhance the convergence of regular Riemannian gradients methods by an adaptively changing metric that is based on the current energy. In this work, we extend this approach to the branch of Riemannian conjugate gradient (CG) methods and investigate the arising schemes numerically. Special attention is given to the selection of the momentum parameter in search direction and how this affects the performance of the resulting schemes. As known from similar applications, we find that the choice of the momentum parameter plays a critical role, with certain parameters reducing the number of iterations required to achieve a specified tolerance by a significant factor. Besides the influence of the momentum parameters, we also investigate how the methods with adaptive metric compare to the corresponding realizations with a standard $H^1_0$-metric. As one of our main findings, the results of the numerical experiments show that the Riemannian CG method with the proposed adaptive metric along with a Polak-Ribiére or Hestenes-Stiefel-type momentum parameter show the best performance and highest robustness compared to the other CG methods that were part of our numerical study.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: calculating the ground states in rotating Bose - Einstein condensates (BECs), which can exhibit multi - scale quantum vortex lattices. Specifically, the author focuses on how to minimize an energy functional on a Riemannian manifold to find these ground states.
### Problem Background
1. **Bose - Einstein Condensates (BECs)**
- When a dilute boson gas undergoes a phase transition near absolute zero, it forms a matter state with superfluid characteristics.
- In actual experiments, quantum vortices can be observed in rotating BECs by applying a rotation potential (such as a stirring potential) and a trapping potential (such as a magnetic or optical potential).
2. **Mathematical Description**
- The mathematical description of the ground state is represented by its quantum state \( u: D \to \mathbb{C} \), where \( D \subset \mathbb{R}^d \), and \( d = 2,3 \).
- The density \( |u|^2: D \to \mathbb{R} \) is usually normalized to satisfy the total mass condition \( \int_D |u|^2 \, dx = 1 \).
- The energy functional \( E(u) \) is given by the Gross - Pitaevskii energy functional:
\[
E(u):=\frac{1}{2}\int_D |\nabla u|^2+V|u|^2-\Omega\overline{u}L_3 u+\frac{\kappa}{2}|u|^4 \, dx
\]
where \( V \) is the external trapping potential, \( L_3:=-i(x_1\partial_{x_2}-x_2\partial_{x_1}) \) is the \( x_3 \)-component of the angular momentum, \( \Omega \) is the angular velocity of the stirring potential, and \( \kappa \) is the repulsive parameter of the inter - particle interaction.
### Research Objectives
The goal of the paper is to develop a new numerical method to efficiently calculate the ground states of rotating BECs. Specifically, the author introduces the Riemannian conjugate Sobolev gradient method and explores the influence of different metrics and momentum parameters on the computational efficiency.
### Method Innovation
1. **Riemannian Conjugate Gradient Method**
- Extend the traditional Riemannian gradient method to the conjugate gradient method, combined with the concept of Sobolev gradient.
- Introduce an adaptive metric, that is, dynamically adjust the metric according to the position of the current iteration point, thereby improving the convergence.
2. **Momentum Parameter Selection**
- Explore four common momentum parameter selections: Dai - Yuan, Fletcher - Reeves, Hestenes - Stiefel and Polak - Ribière - type parameters.
- It is found that the Polak - Ribière and Hestenes - Stiefel - type parameters perform best in the experiments, significantly reducing the number of iterations required to reach the specified tolerance.
### Main Findings
- The adaptive metric significantly accelerates the calculation process compared with the standard \( H^1_0 \)-metric.
- The Polak - Ribière and Hestenes - Stiefel - type momentum parameters perform best in all tests, outperforming other types of parameters.
In summary, this paper aims to improve the efficiency and robustness of calculating the ground states of rotating BECs by introducing the new Riemannian conjugate Sobolev gradient method and its variants.