A basis-set based fortran program to solve the Gross-Pitaevskii equation for dilute Bose gases in harmonic and anharmonic traps

DOI10.1016/J.CPC.2005.10.014MaRDI QIDQ709846zbMATHOpenAlexFDO

Authors Rakesh Prabhat Tiwari, Alok Kumar Shukla

Publication date 18 October 2010

Published in Computer Physics Communications (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/cond-mat/0603732

zbMATH Keywords

Bose-Einstein condensationGross-Pitaevskii equationnumerical solutionsanharmonic potential

Mathematics Subject Classification ID

Closed and approximate solutions to the Schrödinger, Dirac, Klein-Gordon and other equations of quantum mechanics (81Q05) Software, source code, etc. for problems pertaining to quantum theory (81-04)

Abstract: Inhomogeneous boson systems, such as the dilute gases of integral spin atoms in low-temperature magnetic traps, are believed to be well described by the Gross-Pitaevskii equation (GPE). GPE is a nonlinear Schroedinger equation which describes the order parameter of such systems at the mean field level. In the present work, we describe a Fortran 90 computer program developed by us, which solves the GPE using a basis set expansion technique. In this technique, the condensate wave function (order parameter) is expanded in terms of the solutions of the simple-harmonic oscillator (SHO) characterizing the atomic trap. Additionally, the same approach is also used to solve the problems in which the trap is weakly anharmonic, and the anharmonic potential can be expressed as a polynomial in the position operators x, y, and z. The resulting eigenvalue problem is solved iteratively using either the self-consistent-field (SCF) approach, or the imaginary time steepest-descent (SD) approach. Our results for harmonic traps are also compared with those published by other authors using different numerical approaches, and excellent agreement is obtained. GPE is also solved for a few anharmonic potentials, and the influence of anharmonicity on the condensate is discussed. Additionally, the notion of Shannon entropy for the condensate wave function is defined and studied as a function of the number of particles in the trap. It is demonstrated numerically that the entropy increases with the particle number in a monotonic way.

Recommendations

• Fortran programs for the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap
• C programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap
• Fortran and C programs for the time-dependent dipolar Gross-Pitaevskii equation in an anisotropic trap
• Numerical solution for the Gross-Pitaevskii equation
• Ground state of the time-independent Gross-Pitaevskii equation

Cites work

• A Mathematical Theory of Communication
• Ground-state solution of Bose--Einstein condensate by directly minimizing the energy functional
• Matrix eigensystem routines - EISPACK guide. 2nd ed
• Some Integrals Involving Hermite Polynomials
• Structure of a quantized vortex in boson systems

Cited in

(9)

• A finite-element toolbox for the stationary Gross-Pitaevskii equation with rotation
• Fortran programs for the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap
• Modeling and computation of Bose-Einstein condensates: stationary states, nucleation, dynamics, stochasticity
• Computing wave functions of nonlinear Schrödinger equations: a time-independent approach
• bose.x
• Fortran and C programs for the time-dependent dipolar Gross-Pitaevskii equation in an anisotropic trap
• GPELab, a Matlab toolbox to solve Gross-Pitaevskii equations. I: Computation of stationary solutions
• A minimisation approach for computing the ground state of Gross-Pitaevskii systems
• GSGPEs: a MATLAB code for computing the ground state of systems of Gross-Pitaevskii equations

Describes a project that uses

Uses Software

• EISPACK

This page was built for publication: A basis-set based fortran program to solve the Gross-Pitaevskii equation for dilute Bose gases in harmonic and anharmonic traps

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q709846)