Projected Sobolev Gradient Flows for Computing Ground State of Ultracold Dipolar Fermi Gas Based on Density Functional Theory
摘要
The ground state for the ultracold Fermi gas with dipole-dipole interaction is a functional minimization problem based on density functional theory (DFT). We extend the recent work on Sobolev gradient flows for the Gross-Pitaevskii eigenvalue problem [SIAM J. Numer. Anal. 58(3), 1744–1772 (2020)], and present continuous projected Sobolev gradient flows for computing the DFT-based ground state solution of ultracold dipolar Fermi gas. We prove that the gradient flows have the properties of orthonormality preserving and energy diminishing, which are desirable for the computation of the ground state solution. Many numerical techniques for partial differential equations can be used to discretize the time-dependent projected Sobolev gradient flows, which may be an advantage of the method. We propose an efficient and accurate numerical scheme – semi-implicit Euler method in time and Fourier spectral method in space for discretizing these projected Sobolev gradient flows and use them to find the ground state of the ultracold dipolar Fermi gas numerically. Extensive numerical examples in three dimensions for the ground state are reported to demonstrate the power of the numerical methods.