Identification of vortices in quantum fluids: Finite element algorithms and programs.

We present finite-element numerical algorithms for the identification of vortices in quantum fluids described by a macroscopic complex wave function. Their implementation using the free software FreeFem++ is distributed with this paper as a post-processing toolbox that can be used to analyse numerical or experimental data. Applications for Bose-Einstein condensates (BEC) and superfluid helium flows are presented. Programs are tested and validated using either numerical data obtained by solving the Gross-Pitaevskii equation or experimental images of rotating BEC. Vortex positions are computed as topological defects (zeros) of the wave function when numerical data are used. For experimental images, we compute vortex positions as local minima of the atomic density, extracted after a simple image processing. Once vortex centers are identified, we use a fit with a Gaussian to precisely estimate vortex radius. For vortex lattices, the lattice parameter (inter-vortex distance) is also computed. The post-processing toolbox offers a complete description of vortex configurations in superfluids. Tests for two-dimensional (giant vortex in rotating BEC, Abrikosov vortex lattice in experimental BEC) and three-dimensional (vortex rings, Kelvin waves and quantum turbulence fields in superfluid helium) configurations show the robustness of the software. The communication with programs providing the numerical or experimental wave function field is simple and intuitive. The post-processing toolbox can be also applied for the identification of vortices in superconductors. Program Title: FFEM_postproc_data_2D, FFEM_postproc_data_3D and FFEM_postproc_image CPC Library link to program files: https://doi.org/10.17632/bhrtdbkxmy.1 Licensing provisions: Apache 2.0 Programming language: FreeFem++ (v 4.12) free software (www.freefem.org) Nature of problem: The software is scoped to the identification of quantized vortices in 2D or 3D configurations described by a complex wave function. Either simulation data obtained through numerical resolution of the Gross-Pitaevskii equation or 2D experimental images of Bose-Einstein condensates (BEC) can be used as input data. Examples of vortex identification in rotating BECs and superfluid helium are illustrated in the paper. Solution method: We model the complex wave function using P1 (piece-wise linear) Galerkin triangular (in 2D) and tetrahedral (in 3D) finite elements. For simulation data, the zeros of the wave function are directly computed and the circulation of the velocity on triangles sides is assessed. Vortices are obtained when a zero is found in a triangle with non-zero circulation. In 3D, the algorithm is applied on tetrahedron faces and the zero points are linked using mesh connectivity to form vortex lines. For experimental images, vortices are identified as local minima of the atomic density and extracted using the image contrast. Additional comments including restrictions and unusual features: Running time is from seconds to minutes depending on the mesh resolution of the simulation/experimental data. [ABSTRACT FROM AUTHOR]