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

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.

Program summary

Title of program:bose.xCatalogue identifier:ADWZ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADWZ_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format:tar.gzComputers:PC's/Linux, Sun Ultra 10/Solaris, HP Alpha/Tru64, IBM/AIXProgramming language used:mostly Fortran 90Number of bytes in distributed program, including test data, etc.:27 313Number of lines in distributed program, including test data, etc.:28 015Card punching code:ASCIINature of physical problem:It is widely believed that the static properties of dilute Bose condensates, as obtained in atomic traps, can be described to a fairly good accuracy by the time-independent Gross-Pitaevskii equation. This program presents an efficient approach of solving this equation.Method of solution:The solutions of the Gross-Pitaevskii equation corresponding to the condensates in atomic traps are expanded as linear combinations of simple-harmonic oscillator eigenfunctions. Thus, the Gross-Pitaevskii equation which is a second-order nonlinear differential equation, is transformed into a matrix eigenvalue problem. Thereby, its solutions are obtained in a self-consistent manner, using methods of computational linear algebra.Unusual features of the program:None     

From data to probability densities without histograms

When one deals with data drawn from continuous variables, a histogram is often inadequate to display their probability density. It deals inefficiently with statistical noise, and binsizes are free parameters. In contrast to that, the empirical cumulative distribution function (obtained after sorting the data) is parameter free. But it is a step function, so that its differentiation does not give a smooth probability density. Based on Fourier series expansion and Kolmogorov tests, we introduce a simple method, which overcomes this problem. Error bars on the estimated probability density are calculated using a jackknife method. We give several examples and provide computer code reproducing them. You may want to look at the corresponding figures 4 to 9 first.

Program summary

Program title: cdf_to_pdCatalogue identifier: AEBC_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBC_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.: 2758No. of bytes in distributed program, including test data, etc.: 18 594Distribution format: tar.gzProgramming language: Fortran 77Computer: Any capable of compiling and executing Fortran codeOperating system: Any capable of compiling and executing Fortran codeClassification: 4.14, 9Nature of problem: When one deals with data drawn from continuous variables, a histogram is often inadequate to display the probability density. It deals inefficiently with statistical noise, and binsizes are free parameters. In contrast to that, the empirical cumulative distribution function (obtained after sorting the data) is parameter free. But it is a step function, so that its differentiation does not give a smooth probability density.Solution method: Based on Fourier series expansion and Kolmogorov tests, we introduce a simple method, which overcomes this problem. Error bars on the estimated probability density are calculated using a jackknife method. Several examples are included in the distribution file.Running time: The test runs provided take only a few seconds to execute.     

TiReX: Replica-exchange molecular dynamics using Tinker

We present a driver program for performing replica-exchange molecular dynamics simulations with the Tinker package. Parallelization is based on the Message Passing Interface, with every replica assigned to a separate process. The algorithm is not communication intensive, which makes the program suitable for running even on loosely coupled cluster systems. Particular attention is paid to the practical aspects of analyzing the program output.

Program summary

Program title: TiReXCatalogue identifier: AEEK_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEEK_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.: 43 385No. of bytes in distributed program, including test data, etc.: 502 262Distribution format: tar.gzProgramming language: Fortran 90/95Computer: Most UNIX machinesOperating system: LinuxHas the code been vectorized or parallelized?: parallelized with MPIClassification: 16.13External routines: TINKER version 4.2 or 5.0, built as a libraryNature of problem: Replica-exchange molecular dynamics.Solution method: Each replica is assigned to a separate process; temperatures are swapped between replicas at regular time intervals.Running time: The sample run may take up to a few minutes.     

Solution of few-body problems with the stochastic variational method II: Two-dimensional systems

A computational approach is presented for efficient solution of two-dimensional few-body problems, such as quantum dots or excitonic complexes, using the stochastic variational method. The computer program can be used to calculate the energies and wave functions of various two-dimensional systems.

Program summary

Program title: svm-2dCatalogue identifier: AEBE_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBE_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.: 5091No. of bytes in distributed program, including test data, etc.: 130 963Distribution format: tar.gzProgramming language: Fortran 90Computer: The program should work on any system with a Fortran 90 compilerOperating system: The program should work on any system with a Fortran 90 compilerClassification: 7.3Nature of problem: Variational calculation of energies and wave functions using Correlated Gaussian basis.Solution method: Two-dimensional few-electron problems are solved by the variational method. The ground state wave function is expanded into Correlated Gaussian basis functions and the parameters of the basis states are optimized by a stochastic selection procedure. Accurate results can be obtained for 2-6 electron systems.Running time: A couple of hours for a typical system.     

SUSY Les Houches Accord 2 I/O made easy

A library for reading and writing data in the SUSY Les Houches Accord 2 format is presented. The implementation is in native Fortran 77. The data are contained in a single array conveniently indexed by preprocessor statements.

Program summary

Program title: SLHA2LibCatalogue identifier: AEDY_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDY_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.: 7550No. of bytes in distributed program, including test data, etc.: 160 123Distribution format: tar.gzProgramming language: FortranComputer: For the build process, a Fortran 77 compiler in a Unixish environment (make, shell) are requiredOperating system: Linux, Mac OS, Windows (Cygwin), Tru64 UnixRAM: The SLHA Record is currently 88 944 bytes longClassification: 4.14, 11.6Nature of problem: Exchange of SUSY parameters and decay information in an ASCII file format.Solution method: The SLHA2Lib provides routines for reading and writing files in the SUSY Les Houches Accord 2 format, a common interchange format for SUSY parameters and decay data.Restrictions: The fixed-sized array that holds the SLHA2 data necessarily limits the amount of decay data that can be stored. This limit can be enlarged by editing and re-running the SLHA2.m program.Unusual features: Data are transported in a single “double complex” array in Fortran, indexed through preprocessor macros. This is about the simplest conceivable container and needs neither dynamic memory allocation nor Fortran extension like structures.Running time: Both reading and writing a SLHA file are typically in the range of a few milliseconds.     

An object-oriented C++ implementation of Davidson method for finding a few selected extreme eigenpairs of a large, sparse, real, symmetric matrix

A C++ class named Davidson is presented for determining a few eigenpairs with lowest or alternatively highest values of a large, real, symmetric matrix. The algorithm described by Stathopoulos and Fischer is used. The exception mechanism is involved to report the errors. The class is written in ANSI C++, so it is fully portable. In addition a console program as well as a program with graphical user interface for Microsoft Windows is attached, which allow one to calculate the lowest eigenstates of time-independent Schrödinger equation for a given binding potential in one, two or three spatial dimensions. The package contains the classes providing often used potential functions (model atom potential, Coulomb potential, square well potential and Kramers-Henneberger well potential) as well as a possibility to use any potential stored in a file (then any dimensionality of the problem is allowed).The described code is the subject of M.Sc. thesis of T.D. prepared under the supervision of J.M.

Program summary

Program title: DavidsonCatalogue identifier: ADZM_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZM_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.: 3 037 055No. of bytes in distributed program, including test data, etc.: 20 002 609Distribution format: tar.gzProgramming language: C++Computer: AllOperating system: AnyRAM: User's parameters dependentWord size: 32 and 64 bitsSupplementary material: Test results for the 2D and 3D cases is availableClassification: 4, 4.8Nature of problem: Finding a few extreme eigenpairs of a real, symmetric, sparse matrix. Examples in quantum optics (interaction of matter with a laser field).Solution method: Davidson algorithmRunning time: The test example included in the distribution package (1D matrix) takes approximately 30 minutes to run. 2D matrix calculations can take hours and 3D, days, to run.     

Motion4D - A library for lightrays and timelike worldlines in the theory of relativity

The Motion4D-library solves the geodesic equation as well as the parallel- and Fermi-Walker-transport in four-dimensional Lorentzian spacetimes numerically. Initial conditions are given with respect to natural local tetrads which are adapted to the symmetries or the coordinates of the spacetime. Beside some already implemented metrics like the Schwarzschild and Kerr metric, the object oriented structure of the library permits to implement other metrics or integrators in a straight forward manner.

Program summary

Program title: Motion4D-libraryCatalogue identifier: AEEX_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEEX_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.: 150 425No. of bytes in distributed program, including test data, etc.: 5 139 407Distribution format: tar.gzProgramming language: C++Computer: All platforms with a C++ compilerOperating system: Linux, Unix, WindowsRAM: 39 MBytesClassification: 1.5External routines: Gnu Scientific Library (GSL) (http://www.gnu.org/software/gsl/)Nature of problem: Solve geodesic equation, parallel and Fermi-Walker transport in four-dimensional Lorentzian spacetimes.Solution method: Integration of ordinary differential equationsRunning time: The test runs provided with the distribution require only a few seconds to run.     

PArthENoPE: Public algorithm evaluating the nucleosynthesis of primordial elements

We describe a program for computing the abundances of light elements produced during Big Bang Nucleosynthesis which is publicly available at http://parthenope.na.infn.it/. Starting from nuclear statistical equilibrium conditions the program solves the set of coupled ordinary differential equations, follows the departure from chemical equilibrium of nuclear species, and determines their asymptotic abundances as function of several input cosmological parameters as the baryon density, the number of effective neutrino, the value of cosmological constant and the neutrino chemical potential. The program requires commercial NAG library routines.

Program summary

Program title: PArthENoPECatalogue identifier: AEAV_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAV_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.: 10 033No. of bytes in distributed program, including test data, etc.: 46 002Distribution format: tar.gzProgramming language: Fortran 77Computer: PC-compatible running Fortran on Unix, MS Windows or LinuxOperating system: Windows 2000, Windows XP, LinuxClassification: 1.2, 1.9, 17.8External routines: NAG LibrariesNature of problem: Computation of yields of light elements synthesized in the primordial universe.Solution method: BDF method for the integration of the ODEs, implemented in a NAG routine.Running time: 90 sec with default parameters on a Dual Xeon Processor 2.4 GHz with 2 GB RAM.     

golem95: A numerical program to calculate one-loop tensor integrals with up to six external legs

We present a program for the numerical evaluation of form factors entering the calculation of one-loop amplitudes with up to six external legs. The program is written in Fortran95 and performs the reduction to a certain set of basis integrals numerically, using a formalism where inverse Gram determinants can be avoided. It can be used to calculate one-loop amplitudes with massless internal particles in a fast and numerically stable way.

Program summary

Program title: golem95_v1.0Catalogue identifier: AEEO_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEEO_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.: 50 105No. of bytes in distributed program, including test data, etc.: 241 657Distribution format: tar.gzProgramming language: Fortran95Computer: Any computer with a Fortran95 compilerOperating system: Linux, UnixRAM: RAM used per form factor is insignificant, even for a rank six six-point form factorClassification: 4.4, 11.1External routines: Perl programming language (http://www.perl.com/)Nature of problem: Evaluation of one-loop multi-leg tensor integrals occurring in the calculation of next-to-leading order corrections to scattering amplitudes in elementary particle physics.Solution method: Tensor integrals are represented in terms of form factors and a set of basic building blocks (“basis integrals”). The reduction to the basis integrals is performed numerically, thus avoiding the generation of large algebraic expressions.Restrictions: The current version contains basis integrals for massless internal particles only. Basis integrals for massive internal particles will be included in a future version.Running time: Depends on the nature of the problem. A rank 6 six-point form factor at a randomly chosen kinematic point takes 0.13 seconds on an Intel Core 2 Q9450 2.66 GHz processor, without any optimisation. With compiler optimisation flag -O3 the same point takes 0.09 seconds. Timings for lower point form factors are: All form factors for five-point functions from rank 0 to rank 4: 0.04 s. All form factors for rank 5 five-point functions: 0.05 s. All form factors for four-point functions from rank 0 to rank 4: 0.01 s.     

A Fortran program for calculating electron or hole mobility in disordered semiconductors from first-principles

A Fortran program is developed to calculate charge carrier (electron or hole) mobility in disordered semiconductors from first-principles. The method is based on non-adiabatic ab initio molecular dynamics and static master equation, treating dynamic and static disorder on the same footing. We have applied the method to calculate the hole mobility in disordered poly(3-hexylthiophene) conjugated polymers as a function of temperature and electric field and obtained excellent agreements with experimental results. The program could be used to explore structure–mobility relation in disordered semiconducting polymers/organic semiconductors and aid rational design of these materials.

Program summary

Program title: FPMuCatalogue identifier: AEJV_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEJV_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.: 788 580No. of bytes in distributed program, including test data, etc.: 8 433 024Distribution format: tar.gzProgramming language: Fortran 90Computer: Any architecture with a Fortran 90 compilerOperating system: Linux, WindowsRAM: Proportional to the system size, in our example, 1.2 GBClassification: 7.9Nature of problem: Determine carrier mobility from first-principles in disordered semiconductors as a function of temperature, electric field and carrier concentration.Solution method: Iteratively solve master equation with carrier state energy and transition rates determined from first-principles.Restrictions: Mobility for disordered semiconductors where the carrier wave-functions are localized and the carrier transport is due to phonon-assisted hopping mechanism.Running time: Depending on the system size (about an hour for the example here).     

Program package for multicanonical simulations of U(1) lattice gauge theory

We document our Fortran 77 code for multicanonical simulations of 4D U(1) lattice gauge theory in the neighborhood of its phase transition. This includes programs and routines for canonical simulations using biased Metropolis heatbath updating and overrelaxation, determination of multicanonical weights via a Wang-Landau recursion, and multicanonical simulations with fixed weights supplemented by overrelaxation sweeps. Measurements are performed for the action, Polyakov loops and some of their structure factors. Many features of the code transcend the particular application and are expected to be useful for other lattice gauge theory models as well as for systems in statistical physics.

Program summary

Program title: STMC_U1MUCACatalogue identifier: AEET_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEET_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.: 18 376No. of bytes in distributed program, including test data, etc.: 205 183Distribution format: tar.gzProgramming language: Fortran 77Computer: Any capable of compiling and executing Fortran codeOperating system: Any capable of compiling and executing Fortran codeClassification: 11.5Nature of problem: Efficient Markov chain Monte Carlo simulation of U(1) lattice gauge theory close to its phase transition. Measurements and analysis of the action per plaquette, the specific heat, Polyakov loops and their structure factors.Solution method: Multicanonical simulations with an initial Wang-Landau recursion to determine suitable weight factors. Reweighting to physical values using logarithmic coding and calculating jackknife error bars.Running time: The prepared tests runs took up to 74 minutes to execute on a 2 GHz PC.     

A completely open source based computing system for computer generation of Fourier holograms

Computer generated holograms are usually generated using commercial software like MATLAB, MATHCAD, Mathematica, etc. This work is an approach in doing the same using freely distributed open source packages and Operating System. A Fourier hologram is generated using this method and tested for simulated and optical reconstruction. The reconstructed images are in good agreement with the objects chosen. The significance of using such a system is also discussed.

Program summary

Program title: FHOLOCatalogue identifier: AEDS_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDS_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.: 176 336No. of bytes in distributed program, including test data, etc.: 4 294 872Distribution format: tar.gzProgramming language: C++Computer: any X86 micro computerOperating system: Linux (Debian Etch)RAM: 512 MBClassification: 18Nature of problem: To generate a Fourier Hologram in micro computer only by using open source operating system and packages.Running time: Depends on the matrix size. 10 sec for a matrix of size 256×256.     

Visual tool for estimating the fractal dimension of images

This work presents a new Visual Basic 6.0 application for estimating the fractal dimension of images, based on an optimized version of the box-counting algorithm. Following the attempt to separate the real information from “noise”, we considered also the family of all band-pass filters with the same band-width (specified as parameter). The fractal dimension can be thus represented as a function of the pixel color code. The program was used for the study of paintings cracks, as an additional tool which can help the critic to decide if an artistic work is original or not.

Program summary

Program title: Fractal Analysis v01Catalogue identifier: AEEG_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEEG_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.: 29 690No. of bytes in distributed program, including test data, etc.: 4 967 319Distribution format: tar.gzProgramming language: MS Visual Basic 6.0Computer: PCOperating system: MS Windows 98 or laterRAM: 30MClassification: 14Nature of problem: Estimating the fractal dimension of images.Solution method: Optimized implementation of the box-counting algorithm. Use of a band-pass filter for separating the real information from “noise”. User friendly graphical interface.Restrictions: Although various file-types can be used, the application was mainly conceived for the 8-bit grayscale, windows bitmap file format.Running time: In a first approximation, the algorithm is linear.     

Milne, a routine for the numerical solution of Milne's problem

The routine Milne provides accurate numerical values for the classical Milne's problem of neutron transport for the planar one speed and isotropic scattering case. The solution is based on the Case eigen-function formalism. The relevant X functions are evaluated accurately by the Double Exponential quadrature. The calculated quantities are the extrapolation distance and the scalar and the angular fluxes. Also, the H function needed in astrophysical calculations is evaluated as a byproduct.

Program summary

Program title: MilneCatalogue identifier: AEGS_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGS_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.: 701No. of bytes in distributed program, including test data, etc.: 6845Distribution format: tar.gzProgramming language: Fortran 77Computer: PC under Linux or WindowsOperating system: Ubuntu 8.04 (Kernel version 2.6.24-16-generic), Windows-XPClassification: 4.11, 21.1, 21.2Nature of problem: The X functions are integral expressions. The convergence of these regular and Cauchy Principal Value integrals are impaired by the singularities of the integrand in the complex plane. The DE quadrature scheme tackles these singularities in a robust manner compared to the standard Gauss quadrature.Running time: The test included in the distribution takes a few seconds to run.     

Fortran programs for the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

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 summary

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).

Program summary (1)

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.

Program summary (2)

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.

Program summary (3)

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.

Program summary (4)

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.

Program summary (5)

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.

Program summary (6)

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.

Program summary (7)

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.

Program summary (8)

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.

Program summary (9)

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.

Program summary (10)

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.

Program summary (11)

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.

Program summary (12)

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.     

Algorithmic derivation of Dyson-Schwinger equations

We present an algorithm for the derivation of Dyson-Schwinger equations of general theories that is suitable for an implementation within a symbolic programming language. Moreover, we introduce the Mathematica package DoDSE¹ which provides such an implementation. It derives the Dyson-Schwinger equations graphically once the interactions of the theory are specified. A few examples for the application of both the algorithm and the DoDSE package are provided.

Program summary

Program title: DoDSECatalogue identifier: AECT_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECT_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.: 105 874No. of bytes in distributed program, including test data, etc.: 262 446Distribution format: tar.gzProgramming language: Mathematica 6 and higherComputer: all on which Mathematica is availableOperating system: all on which Mathematica is availableClassification: 11.1, 11.4, 11.5, 11.6Nature of problem: Derivation of Dyson-Schwinger equations for a theory with given interactions.Solution method: Implementation of an algorithm for the derivation of Dyson-Schwinger equations.Unusual features: The results can be plotted as Feynman diagrams in Mathematica.Running time: Less than a second to minutes for Dyson-Schwinger equations of higher vertex functions.     

LATTICEEASY: A program for lattice simulations of scalar fields in an expanding universe

We describe a C++ program that we have written and made available for calculating the evolution of interacting scalar fields in an expanding universe. The program is particularly useful for the study of reheating and thermalization after inflation. The program and its full documentation are available on the Web at http://www.science.smith.edu/departments/Physics/fstaff/gfelder/latticeeasy/. In this paper we provide a brief overview of what the program does and what it is useful for.

Program summary

Program title: LATTICEEASYCatalog identifier: AEAW_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAW_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.: 2579No. of bytes in distributed program, including test data, etc.: 34 521Distribution format: tar.gzProgramming language: C++Computer: AnyOperating system: AnyRAM: Typically 4 MB to 800 MBClassification: 1.9Nature of problem: After inflation the universe consisted of interacting fields in a high energy, nonthermal state [1]. The evolution of these fields can not be described with standard approximation techniques such as linearization, kinetic theory, or Hartree expansion, and must thus be simulated numerically. Fortunately, the fields rapidly acquire large occupation numbers over a range of frequencies, so their evolution can be accurately modeled with classical field theory [2]. The specific fields and interactions relevant at these high energies are not known, so different models must be tested phenomenologically.Solution method: LATTICEEASY solves the equations of motion for interacting scalar fields in an expanding universe. The user describes a particular theory by entering the field potential and its derivatives in a “model file” and the program then uses a staggered leapfrog method to evolve the field equations and Friedmann equation for the fields and the expansion of the universe.Restrictions: In its current form LATTICEEASY only includes scalar fields and does not include metric perturbations.Running time: The running time can range from minutes to weeks.References: [1] A.D. Linde, Particle Physics and Inflationary Cosmology, Harwood, Chur, Switzerland, 1990. [2] S. Khlebnikov, I. Tkachev, Phys. Rev. Lett. 77 (1996) 219, hep-ph 9603378.     

Object-oriented design patterns in Fortran 90/95: mazev1, mazev2 and mazev3

This paper discusses the concept, application, and usefulness of software design patterns for scientific programming in Fortran 90/95. An example from the discipline of object-oriented design patterns, that of a game based on navigation through a maze, is used to describe how some important patterns can be implemented in Fortran 90/95 and how the progressive introduction of design patterns can usefully restructure Fortran software as it evolves. This example is complemented by a discussion of how design patterns have been used in a real-life simulation of Particle-in-Cell plasma physics. The following patterns are mentioned in this paper: Factory, Strategy, Template, Abstract Factory and Facade.

Program summary

Program title: mazev1, mazev2, mazev3Catalogue identifier: AEAI_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAI_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.: 1958No. of bytes in distributed program, including test data, etc.: 17 100Distribution format: tar.gzProgramming language: Fortran 95Computer: PC/MacOperating system: Unix/Linux/Mac (FreeBSD)/Windows (Cygwin)RAM: These are interactive programs with small (KB) memory requirementsClassification: 6.5, 20Nature of problem: A sequence of programs which demonstrate the use of object oriented design patterns for the restructuring of Fortran 90/95 software. The programs implement a simple maze game similar to that described in [1].Solution method: Restructuring uses versions of the Template, Strategy and Factory design patterns.Running time: Interactive.References:

[1] 

E. Gamma, R. Helm, R. Johnson, J. Vlissides, Design Patterns: Elements of Reusable Object Oriented Software, Addison-Wesley, 1995, ISBN 0201633612.