uu.seUppsala University Publications
Change search
ReferencesLink to record
Permanent link

Direct link
Numerical simulation code for self-gravitating Bose-Einstein condensates
Uppsala University, Disciplinary Domain of Science and Technology, Physics, Department of Physics and Astronomy.
2013 (English)In: Computer Physics Communications, ISSN 0010-4655, E-ISSN 1879-2944, Vol. 184, no 4, 1339-1343 p.Article in journal (Refereed) Published
Abstract [en]

We completed the development of simulation code that is designed to study the behavior of a conjectured dark matter galactic halo that is in the form of a Bose-Einstein Condensate (BEC). The BEC is described by the Gross-Pitaevskii equation, which can be solved numerically using the Crank-Nicholson method. The gravitational potential, in turn, is described by Poisson's equation, that can be solved using the relaxation method. Our code combines these two methods to study the time evolution of a self-gravitating BEC. The inefficiency of the relaxation method is balanced by the fact that in subsequent time iterations, previously computed values of the gravitational field serve as very good initial estimates. The code is robust (as evidenced by its stability on coarse grids) and efficient enough to simulate the evolution of a system over the course of 10(9) years using a finer (100 x 100 x 100) spatial grid, in less than a day of processor time on a contemporary desktop computer. Program summary Program title: bec3p Catalogue identifier: AEOR_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEOR_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 5248 No. of bytes in distributed program, including test data, etc.: 715402 Distribution format: tar.gz Programming language: C++ or FORTRAN. Computer: PCs or workstations. Operating system: Linux or Windows. 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 x 100 x 100 grid can be solved with a time span of 10 Gyr in approx. a day of running time.

Place, publisher, year, edition, pages
2013. Vol. 184, no 4, 1339-1343 p.
Keyword [en]
Gravity, Poisson's equation, Gross-Pitaevskii equation, Dark matter, Galaxy rotation
National Category
Natural Sciences
URN: urn:nbn:se:uu:diva-197959DOI: 10.1016/j.cpc.2012.12.024ISI: 000315974100027OAI: oai:DiVA.org:uu-197959DiVA: diva2:615095
Available from: 2013-04-08 Created: 2013-04-08 Last updated: 2013-04-08Bibliographically approved

Open Access in DiVA

No full text

Other links

Publisher's full text

Search in DiVA

By author/editor
Madarassy, Enikö J. M.
By organisation
Department of Physics and Astronomy
In the same journal
Computer Physics Communications
Natural Sciences

Search outside of DiVA

GoogleGoogle Scholar
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Altmetric score

Total: 238 hits
ReferencesLink to record
Permanent link

Direct link