Skip to Main content Skip to Navigation
Preprints, Working Papers, ...

A finite-element toolbox for the stationary Gross-Pitaevskii equation with rotation

Guillaume Vergez 1, 2 Ionut Danaila 1 Sylvain Auliac 2 Frédéric Hecht 2, 3
3 ALPINES - Algorithms and parallel tools for integrated numerical simulations
LJLL - Laboratoire Jacques-Louis Lions, INSMI - Institut National des Sciences Mathématiques et de leurs Interactions, Inria de Paris
Abstract : We present a new numerical system using classical finite elements with mesh adaptivity for computing stationary solutions of the Gross-Pitaevskii equation. The programs are written as a toolbox for FreeFem++ (, a free finite-element software available for all existing operating systems. This offers the advantage to hide all technical issues related to the implementation of the finite element method, allowing to easily implement various numerical algorithms.Two robust and optimised numerical methods were implemented to minimize the Gross-Pitaevskii energy: a steepest descent method based on Sobolev gradients and a minimization algorithm based on the state-of-the-art optimization library Ipopt. For both methods, mesh adaptivity strategies are implemented to reduce the computational time and increase the local spatial accuracy when vortices are present. Different run cases are made available for 2D and 3D configurations of Bose-Einstein condensates in rotation. An optional graphical user interface is also provided, allowing to easily run predefined cases or with user-defined parameter files. We also provide several post-processing tools (like the identification of quantized vortices) that could help in extracting physical features from the simulations. The toolbox is extremely versatile and can be easily adapted to deal with different physical models.
Complete list of metadatas
Contributor : Guillaume Vergez <>
Submitted on : Tuesday, February 23, 2016 - 12:48:40 AM
Last modification on : Friday, March 27, 2020 - 4:04:31 AM
Document(s) archivé(s) le : Sunday, November 13, 2016 - 1:37:08 AM


Files produced by the author(s)


  • HAL Id : hal-01277660, version 1
  • ARXIV : 1602.07071


Guillaume Vergez, Ionut Danaila, Sylvain Auliac, Frédéric Hecht. A finite-element toolbox for the stationary Gross-Pitaevskii equation with rotation. 2016. ⟨hal-01277660⟩



Record views


Files downloads