Abstract:We develop a simulation method for evaporative cooling of trapped Bose-Einstein condensate at finite temperatures using Zaremba-Nikuni-Griffin (ZNG) formalism. ZNG formalism includes the generalized GP equation and a semiclassical kinetic equation for the thermal cloud, which treats the excitations semiclassically within the Hartree Fock approximation. The generalized GP equation includes the mean field due to the thermal cloud and the source term associated with collisions between the condensate and the thermal cloud. Our method is based on the numerical approach developed by Jackson and Zaremba, which simulates the kinetic equation using test particles. A key point of our method is to mimic the evaporative cooling process by eliminating the test particles with high energy. We show that our method successfully describes condensate growth during evaporative cooling. We also numerically simulate vortex lattice formation during evaporative cooling in the presence of the rotating thermal cloud.
What problem does this paper attempt to address?
This paper aims to solve the problem of how to simulate the evaporative cooling process of Bose - Einstein condensation (BEC) at a finite temperature, especially the formation of vortex lattices in rotating Bose gases. Specifically, the researchers developed a simulation method based on Zaremba - Nikuni - Griffin (ZNG) formalism to describe the growth of condensates and their interactions with the thermal cloud during the evaporative cooling process of trapped Bose gases at a finite temperature.
### Main problems:
1. **Simulation of evaporative cooling**: How to simulate the evaporative cooling process of Bose gases at a finite temperature by numerical methods, especially when the gas is in a rotating state.
2. **Formation of vortex lattices**: How do vortices form in rotating Bose gases during the evaporative cooling process, and how do these vortices affect the dynamic behavior of the system?
### Research background:
Experimentally, atomic Bose - Einstein condensation is achieved by evaporatively cooling trapped Bose atoms. Evaporative cooling techniques include selectively removing high - energy atoms and collision relaxation of the remaining gas. In this finite - temperature dynamic process, the condensate and the thermal cloud coexist. Many experimental phenomena show that the existence of the thermal cloud has a significant impact on the system, such as the formation of vortex lattices and Landau damping.
### Method overview:
To simulate the evaporative cooling process, the authors used the ZNG formalism, which includes the generalized Gross - Pitaevskii (GP) equation and the semiclassical Boltzmann kinetic equation. The generalized GP equation describes the behavior of the condensate, while the Boltzmann equation describes the distribution function of the thermal cloud. Evaporative cooling is simulated by removing test particles with energy higher than the cutoff energy. The key points are:
- **Generalized GP equation**:
\[
i\hbar\frac{\partial\Phi(\mathbf{r},t)}{\partial t}=\left[-\frac{\hbar^{2}\nabla^{2}}{2m}+V_{\text{ext}}(\mathbf{r},t)+gn_{c}(\mathbf{r},t)+2g\tilde{n}(\mathbf{r},t)-iR(\mathbf{r},t)\right]\Phi(\mathbf{r},t)
\]
where \(n_{c} = |\Phi|^{2}\) is the condensate density, and \(\tilde{n}=\int\frac{d\mathbf{p}}{(2\pi\hbar)^{3}}f(\mathbf{p},\mathbf{r},t)\) is the non - condensate density.
- **Semiclassical Boltzmann equation**:
\[
\frac{\partial f(\mathbf{p},\mathbf{r},t)}{\partial t}+\frac{\mathbf{p}}{m}\cdot\nabla f(\mathbf{p},\mathbf{r},t)-\nabla U\cdot\nabla_{\mathbf{p}}f(\mathbf{p},\mathbf{r},t)=C_{12}[f]+C_{22}[f]
\]
where \(U = V_{\text{ext}}+2gn_{c}+2g\tilde{n}\) is the effective potential energy of the non - condensate.
### Simulation results:
- **Condensate growth**: Through evaporative cooling, the condensate fraction \(N_{c}/N\) increases with time and finally reaches a temperature of about \(T\sim0.56T_{c}\).
- **Vortex lattice formation**: In rotating Bose gases, a vortex lattice is formed during the evaporative cooling process, and these vortices gradually decay after the cooling stops.
### Conclusion:
The authors have successfully developed an evaporative - cooling simulation method based on the ZNG formalism, which can effectively describe the condensate growth and vortex lattice formation of Bose gases at a finite temperature. However, since the ZNG formalism is based on...