Record Details

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

DSpace at IIT Bombay

View Archive Info
 
 
Field Value
 
Title A basis-set based Fortran program to solve the Gross-Pitaevskii equation for dilute Bose gases in harmonic and anharmonic traps
 
Creator TIWARI, RAKESH PRABHAT
SHUKLA, ALOK
 
Subject harmonic analysis
magnetic properties
parameter estimation
thermal effects
oscillators (electronic)
 
Description 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 Schrödinger 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. Iterations can be initiated using either the simple-harmonic-oscillator ground state solution, or the Thomas–Fermi (TF) solution. It is found that for condensates containing up to a few hundred atoms, both approaches lead to rapid convergence. However, in the strong interaction limit of condensates containing thousands of atoms, it is the SD approach coupled with the TF starting orbitals, which leads to quick convergence. 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.
 
Publisher Elsevier
 
Date 2009-04-03T09:05:19Z
2011-11-25T20:39:01Z
2011-12-26T13:08:41Z
2011-12-27T05:33:58Z
2009-04-03T09:05:19Z
2011-11-25T20:39:01Z
2011-12-26T13:08:41Z
2011-12-27T05:33:58Z
2006
 
Type Article
 
Identifier Computer Physics Communications 174(12), 966-982
http://dx.doi.org/10.1016/j.cpc.2005.10.014
http://hdl.handle.net/10054/1125
http://dspace.library.iitb.ac.in/xmlui/handle/10054/1125
 
Language en