Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aeor_v1_0.tar.gz(698 Kbytes)
Manuscript Title: Numerical simulation code for self-gravitating Bose-Einstein condensates
Authors: Eniko J. M. Madarassy, Viktor T. Toth
Program title: bec3p
Catalogue identifier: AEOR_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 184(2013)1339
Programming language: C++ or FORTRAN.
Computer: PCs or workstations.
Operating system: Linux or Windows.
Keywords: Gravity, Poisson's equation, Gross-Pitaevskii equation, Dark matter, Galaxy rotation.
Classification: 1.5.

Nature of problem:
Simulation of a self-gravitating Bose-Einstein condensate by simultaneous solution of the Gross-Pitaevskii and Poisson equations in three dimensions.

Solution method:
The Gross-Pitaevskii equation is solved numerically using the Crank-Nicholson method; Poisson's equation is solved using the relaxation method. The time evolution of the system is governed by the Gross-Pitaevskii equation; the solution of Poisson's equation at each time step is used as an initial estimate for the next time step, which dramatically increases the efficiency of the relaxation method.

Running time:
Depends on the chosen size of the problem. On a typical personal computer, a 100 × 100 × 100 grid can be solved with a time span of 10 Gyr in approx. a day of running time.