Computer Physics Communications > 2009 > 180 > 10 > 1888-1912
Fortran programs for the time-dependent Gross–Pitaevskii equation in a fully anisotropic trap
Fortran programs for the time-dependent Gross–Pitaevskii equation in a fully anisotropic trap
Source
Abstract
Here we develop simple numerical algorithms for both stationary and non-stationary solutions of the time-dependent Gross–Pitaevskii (GP) equation describing the properties of Bose–Einstein condensates at ultra low temperatures. In particular, we consider algorithms involving real- and imaginary-time propagation based on a split-step Crank–Nicolson method. In a one-space-variable form of the GP equation we consider the one-dimensional, two-dimensional circularly-symmetric, and the three-dimensional spherically-symmetric harmonic-oscillator traps. In the two-space-variable form we consider the GP equation in two-dimensional anisotropic and three-dimensional axially-symmetric traps. The fully-anisotropic three-dimensional GP equation is also considered. Numerical results for the chemical potential and root-mean-square size of stationary states are reported using imaginary-time propagation programs for all the cases and compared with previously obtained results. Also presented are numerical results of non-stationary oscillation for different trap symmetries using real-time propagation programs. A set of convenient working codes developed in Fortran 77 are also provided for all these cases (twelve programs in all). In the case of two or three space variables, Fortran 90/95 versions provide some simplification over the Fortran 77 programs, and these programs are also included (six programs in all). Program title: (i) imagetime1d, (ii) imagetime2d, (iii) imagetime3d, (iv) imagetimecir, (v) imagetimesph, (vi) imagetimeaxial, (vii) realtime1d, (viii) realtime2d, (ix) realtime3d, (x) realtimecir, (xi) realtimesph, (xii) realtimeaxialCatalogue identifier: AEDU_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEDU_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 122 907No. of bytes in distributed program, including test data, etc.: 609 662Distribution format: tar.gzProgramming language: FORTRAN 77 and Fortran 90/95Computer: PCOperating system: Linux, UnixRAM: 1 GByte (i, iv, v), 2 GByte (ii, vi, vii, x, xi), 4 GByte (iii, viii, xii), 8 GByte (ix)Classification: 2.9, 4.3, 4.12Nature of problem: These programs are designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in one-, two- or three-space dimensions with a harmonic, circularly-symmetric, spherically-symmetric, axially-symmetric or anisotropic trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Solution method: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation, in either imaginary or real time, over small time steps. The method yields the solution of stationary and/or non-stationary problems.Additional comments: This package consists of 12 programs, see “Program title”, above. FORTRAN77 versions are provided for each of the 12 and, in addition, Fortran 90/95 versions are included for ii, iii, vi, viii, ix, xii. For the particular purpose of each program please see the below.Running time: Minutes on a medium PC (i, iv, v, vii, x, xi), a few hours on a medium PC (ii, vi, viii, xii), days on a medium PC (iii, ix). Title of program: imagtime1d.FTitle of electronic file: imagtime1d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 1 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in one-space dimension with a harmonic trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems. Title of program: imagtimecir.FTitle of electronic file: imagtimecir.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 1 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in two-space dimensions with a circularly-symmetric trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems. Title of program: imagtimesph.FTitle of electronic file: imagtimesph.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 1 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in three-space dimensions with a spherically-symmetric trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems. Title of program: realtime1d.FTitle of electronic file: realtime1d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 2 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in one-space dimension with a harmonic trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems. Title of program: realtimecir.FTitle of electronic file: realtimecir.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 2 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in two-space dimensions with a circularly-symmetric trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems. Title of program: realtimesph.FTitle of electronic file: realtimesph.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 2 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in three-space dimensions with a spherically-symmetric trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems. Title of programs: imagtimeaxial.F and imagtimeaxial.f90Title of electronic file: imagtimeaxial.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 2 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time: Few hours on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in three-space dimensions with an axially-symmetric trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems. Title of program: imagtime2d.F and imagtime2d.f90Title of electronic file: imagtime2d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 2 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time: Few hours on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in two-space dimensions with an anisotropic trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems. Title of program: realtimeaxial.F and realtimeaxial.f90Title of electronic file: realtimeaxial.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 4 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time Hours on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in three-space dimensions with an axially-symmetric trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems. Title of program: realtime2d.F and realtime2d.f90Title of electronic file: realtime2d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 4 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time: Hours on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in two-space dimensions with an anisotropic trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems. Title of program: imagtime3d.F and imagtime3d.f90Title of electronic file: imagtime3d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 4 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time: Few days on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in three-space dimensions with an anisotropic trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems. Title of program: realtime3d.F and realtime3d.f90Title of electronic file: realtime3d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum Ram Memory: 8 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time: Days on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross–Pitaevskii nonlinear partial differential equation in three-space dimensions with an anisotropic trap. The Gross–Pitaevskii equation describes the properties of a dilute trapped Bose–Einstein condensate.Method of solution: The time-dependent Gross–Pitaevskii equation is solved by the split-step Crank–Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems.
Identifiers
journal ISSN : | 0010-4655 |
DOI | 10.1016/j.cpc.2009.04.015 |
Authors
P. Muruganandam
- Instituto de Física Teórica, UNESP – São Paulo State University, Barra Funda, 01140-070 São Paulo, São Paulo, Brazil
- School of Physics, Bharathidasan University, Palkalaiperur Campus, Tiruchirappalli – 620024, Tamil Nadu, India
S.K. Adhikari
- Instituto de Física Teórica, UNESP – São Paulo State University, Barra Funda, 01140-070 São Paulo, São Paulo, Brazil
Keywords
Bose–Einstein condensate Gross–Pitaevskii equation Split-step Crank–Nicolson scheme Real- and imaginary-time propagation Fortran program Partial differential equation
Bose–Einstein condensate Gross–Pitaevskii equation Split-step Crank–Nicolson scheme Real- and imaginary-time propagation Fortran program Partial differential equation
Additional information
Publication languages:
English
Data set: Elsevier
Publisher
Elsevier Science
Fields of science
article