Here we develop simple numerical algorithms for both stationary and non-stationary solutions of the timedependent 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, twodimensional 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 axiallysymmetric 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).

A self-contained discussion of nonrelativistic quantum mechanical potential scattering in two dimensions is presented. The discussion includes, among other topics, partial wave decomposition in coordinate and momentum space, Lippmann–Schwinger integral equations of scattering for the scattering wavefunction and the transition operator, optical theorem, and the unitarity relation for the transition operator. The present definition of the scattering amplitude in terms of the asymptotic wavefunction differs from the usual definition. The present definition has certain advantages, for example, in writing the optical theorem and in studying the analytical properties of the scattering amplitude.

We present C programming language versions of earlier published Fortran programs (Muruganandam and Adhikari (2009) [1]) for calculating both stationary and non-stationary solutions of the time-dependent Gross-Pitaevskii (GP) equation. The GP equation describes the properties of dilute Bose-Einstein condensates at ultra-cold temperatures. C versions of programs use the same algorithms as the Fortran ones, involving real-and imaginary-time propagation based on a split-step Crank-Nicolson method. In a one-spacevariable 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. In addition to these twelve programs, for six algorithms that involve two and three space variables, we have also developed threaded (OpenMP parallelized) programs, which allow numerical simulations to use all available CPU cores on a computer. All 18 programs are optimized and accompanied by makefiles for several popular C compilers. We present typical results for scalability of threaded codes and demonstrate almost linear speedup obtained with the new programs, allowing a decrease in execution times by an order of magnitude on modern multi-core computers. [14][15][16]. This new version represents translation of all programs to the C programming language, which will make it accessible to the wider parts of the corresponding communities. It is well known that numerical simulations of the GP equation in highly experimentally relevant geometries with two or three space variables are computationally very demanding, which presents an obstacle in detailed numerical studies of such systems. For this reason, we have developed threaded (OpenMP parallelized) versions of programs imagtime2d, imagtime3d, imagtimeaxial, realtime2d, realtime3d, realtimeaxial, which are named imagtime2d-th, imagtime3d-th, imagtimeaxial-th, realtime2d-th, realtime3d-th, realtimeaxial-th, respectively. Figure 1 shows the scalability results obtained for OpenMP versions of programs realtime2d and realtime3d. As we can see, the speedup is almost linear, and on a computer with the total of 8 CPU cores we observe in Fig. 1(a) a maximal speedup of around 7, or roughly 90% of the ideal speedup, while on a computer with 12 CPU cores we find in Fig. 1(b) that the maximal speedup is around 9.6, or 80% of the ideal speedup. Such a speedup represents significant improvement in the performance.

We investigate the zero-temperature properties of a superfluid Bose-Fermi mixture by introducing a set of coupled Galilei-invariant nonlinear Schrödinger equations valid from weak coupling to unitarity. The Bose dynamics is described by a Gross-Pitaevskii-type equation including beyond-mean-field corrections possessing the correct weak-coupling and unitarity limits. The dynamics of the two-component Fermi superfluid is described by a density-functional equation including beyond-mean-field terms with correct weak-coupling and unitarity limits. The present set of equations is equivalent to the equations of generalized superfluid hydrodynamics, which take into account also surface effects. The equations describe the mixture properly as the Bose-Bose repulsive ͑positive͒ and Fermi-Fermi attractive ͑negative͒ scattering lengths are varied from zero to infinity in the presence of a Bose-Fermi interaction. The present model is tested numerically as the Bose-Bose and Fermi-Fermi scattering lengths are varied over wide ranges covering the weak-coupling to unitarity transition.

We study certain stationary and time-evolution problems of trapped Bose-Einstein condensates using the numerical solution of the Gross-Pitaevskii equation with both spherical and axial symmetries. We consider time-evolution problems initiated by changing the interatomic scattering length or harmonic trapping potential suddenly in a stationary condensate. These changes introduce oscillations in the condensate which are studied in detail. We use a time iterative split-step method for the solution of the time-dependent Gross-Pitaevskii equation, where all nonlinear and linear nonderivative terms are treated separately from the time propagation with the kinetic energy terms. Even for an arbitrarily strong nonlinear term this leads to extremely accurate and stable results after millions of time iterations of the original equation.

No abstract

It is demonstrated for the three-dimensional three-particle system that a divergence arising from essentially the same singularity structure of the kernel of the scattering integral equation is responsible for both the EAmov and Thomas e6ects. The above divergence implies that the results of a three-particle dynamical calculation will be sensitive to the details of the two-particle interaction.In two-dimensional systems, the above divergence is absent and consequently the three-particle observables become essentially independent of the two-particle interaction model.

Many of the static and dynamic properties of an atomic Bose-Einstein condensate (BEC) are usually studied by solving the mean-field Gross-Pitaevskii (GP) equation, which is a nonlinear partial differential equation for short-range atomic interaction. More recently, BEC of atoms with long-range dipolar atomic interaction are used in theoretical and experimental studies. For dipolar atomic interaction, the GP equation is a partial integro-differential equation, requiring complex algorithm for its numerical solution. Here we present numerical algorithms for both stationary and non-stationary solutions of the full threedimensional (3D) GP equation for a dipolar BEC, including the contact interaction. We also consider the simplified one-(1D) and two-dimensional (2D) GP equations satisfied by cigar-and disk-shaped dipolar BECs. We employ the split-step Crank-Nicolson method with real-and imaginary-time propagations, respectively, for the numerical solution of the GP equation for dynamic and static properties of a dipolar BEC. The atoms are considered to be polarized along the z axis and we consider ten different cases, e.g., stationary and non-stationary solutions of the GP equation for a dipolar BEC in 1D (along x and z axes), 2D (in x-y and x-z planes), and 3D, and we provide working codes in Fortran 90/95 and C for these ten cases (twenty programs in all). We present numerical results for energy, chemical potential, rootmean-square sizes and density of the dipolar BECs and, where available, compare them with results of other authors and of variational and Thomas-Fermi approximations. Program summary Program title: (i) imag1dZ, (ii) imag1dX, (iii) imag2dXY, (iv) imag2dXZ, (v) imag3d, (vi) real1dZ, (vii) real1dX, (viii) real2dXY, (ix) real2dXZ, (x) real3d