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     

Highly precise solutions of the one-dimensional Schrödinger equation with an arbitrary potential

This work concerns the obtaining of a highly accurate solution of the one-dimensional Schrödinger equation with a practically arbitrary potential. The approach is based on the power series method and it is implemented in the ultra high precision mode. It is shown that such an approach yields not only highly precise values of the energies but also accurate wave functions.     

Applications of critical temperature in minimizing functions of continuous variables with simulated annealing algorithm

The simulated annealing (SA) algorithm has been recognized as a powerful technique for minimizing complicated functions. However, a critical disadvantage of the SA algorithm is its high computational cost. Therefore, it is the goal of this paper to investigate the use of the critical temperature in SA to reduce its computational cost. This paper presents a systematic study of the critical temperature and its applications in the minimization of functions of continuous variables with the SA algorithm. Based on this study, a new algorithm was developed to exploit the unique feature of the critical temperature in SA. The new algorithm combines SA and local search to determine global minimum effectively. Extensive tests on a variety of functions demonstrated that the new algorithm provides comparable performance to well-established SA techniques. Furthermore, the new algorithm also improves the determination of the starting temperature for the SA algorithm. The results obtained in this study are expected to be useful for improving the efficiency of SA algorithms, and for facilitating the development of temperature parallel SA algorithms.     

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.     

Fast LP method for the Schrödinger equation

The LP and CP methods are two versions of the piecewise perturbation methods to solve the Schrödinger equation. On each step the potential function is approximated by a constant (for CP) or by a linear function (for LP) and the deviation of the true potential from this approximation is treated by the perturbation theory.This paper is based on the idea that an LP algorithm can be made faster if expressed in a CP-like form. We obtain a version of order 12 whose two main ingredients are a new set of formulae for the computation of the zeroth-order solution which replaces the use of the Airy functions, and a convenient way of expressing the formulae for the perturbation corrections. Tests on a set of eigenvalue problems with a very big number of eigenvalues show that the proposed algorithm competes very well with a CP version of the same order and is by one order of magnitude faster than the LP algorithms existing in the literature. We also formulate a new technique for the step width adjustment and bring some new elements for a better understanding of the energy dependence of the error for the piecewise perturbation methods.     

A new class of hybrid global optimization algorithms for peptide structure prediction: integrated hybrids

A novel class of hybrid global optimization methods for application to the structure prediction in protein-folding problem is introduced. These optimization methods take the form of a hybrid between a deterministic global optimization algorithm, the αBB, and a stochastically based method, conformational space annealing (CSA), and attempt to combine the beneficial features of these two algorithms. The αBB method as previously extant exhibits consistency, as it guarantees convergence to the global minimum for twice-continuously differentiable constrained nonlinear programming problems, but can benefit from improvements in the computational front. Computational studies for met-enkephalin demonstrate the promise for the proposed hybrid global optimization method.     

Numerical errors resulting from finite-difference approximation in computations of a one-dimensional Schrödinger equation with the Trotter formula

The parallel algorithm for solving time-dependent Schrödinger equations devised by De Raedt and based on the Trotter formula is not only simple but also unconditionally stable, explicit, and local. We consider the numerical errors resulting from the finite-difference approximation of De Raedt's algorithm by comparing an exact solution of a free particle with the approximate solution calculated by using the Trotter formula, which depends on the size of the spatial-temporal lattice.     

Finite difference approach for the two-dimensional Schrödinger equation with application to scission-neutron emission

We present a grid-based procedure to solve the eigenvalue problem for the two-dimensional Schrödinger equation in cylindrical coordinates. The Hamiltonian is discretized by using adapted finite difference approximations of the derivatives and this leads to an algebraic eigenvalue problem with a large (sparse) matrix, which is solved by the method of Arnoldi. By this procedure the single particle eigenstates of nuclear systems with arbitrary deformations can be obtained. As an application we have considered the emission of scission neutrons from fissioning nuclei.     

Solution of the Schrödinger equation by a high order perturbation method based on a linear reference potential

The paper is devoted to the enhancement of the accuracy of the line-based perturbation method via the introduction of the perturbation corrections. We effectively construct the first and the second order corrections. We also perform the error analysis to predict that the introduction of successive corrections substantially enhances the order of the method from four, for the zeroth order version, to six and ten when the first and the second-order corrections are included. In order to remove the effect of the accuracy loss due to near-cancellation of like-terms when evaluating the perturbation corrections we construct alternative asymptotic formulae using a Maple code. We also propose a procedure for choosing the step size in terms of the preset accuracy and give a number of numerical illustrations.     

Lattice Boltzmann schemes for quantum applications

We review the basic ideas behind the quantum lattice Boltzmann equation (LBE), and present a few thoughts on the possible use of such an equation for simulating quantum many-body problems on both (parallel) electronic and quantum computers.     

Arbitrarily precise numerical solutions of the one-dimensional Schrödinger equation

In this paper, how to overcome the barrier for a finite difference method to obtain the numerical solutions of a one-dimensional Schrödinger equation defined on the infinite integration interval accurate than the computer precision is discussed. Five numerical examples of solutions with the error less than 10⁻⁵⁰ and 10⁻³⁰ for the bound and resonant state, respectively, obtained by the Obrechkoff one-step method implemented in the multi precision mode, which include the harmonic oscillator, the Pöschl-Teller potential, the Morse potential and the Woods-Saxon potential, demonstrate that the finite difference method can yield the eigenvalues of a complex potential with an arbitrarily desired precision within a reasonable efficiency.     

A new kind of discretization scheme for solving a two-dimensional time-independent Schrödinger equation

In this paper we present a new kind of discretization scheme for solving a two-dimensional time-independent Schrödinger equation. The scheme uses a symmetrical multi-point difference formula to represent the partial differentials of the two-dimensional variables, which can improve the accuracy of the numerical solutions to the order of Δx2Nq+2 when a (2Nq+1)-point formula is used for any positive integer Nq with Δx=Δy, while Nq=1 equivalent to the traditional scheme. On the other hand, the new scheme keeps the same form of the traditional matrix equation so that the standard algebraic eigenvalue algorithm with a real, symmetric, large sparse matrix is still applicable. Therefore, for the same dimension, only a little more CPU time than the traditional one should be used for diagonalizing the matrix. The numerical examples of the two-dimensional harmonic oscillator and the two-dimensional Henon-Heiles potential demonstrate that by using the new method, the error in the numerical solutions can be reduced steadily and extensively through the increase of Nq, which is more efficient than the traditional methods through the decrease of the step size.     

Pseudospectral method based on prolate spheroidal wave functions for semiconductor nanodevice simulation

We solve Schrödinger's equation for semiconductor nanodevices by applying prolate spheroidal wave functions of order zero as basis functions in the pseudospectral method. When the functions involved in the problem are bandlimited, the prolate pseudospectral method outperforms the conventional pseudospectral methods based on trigonometric and orthogonal polynomials and related functions, asymptotically achieving similar accuracy using a factor of π/2 less unknowns than the latter. The prolate pseudospectral method also employs a more uniform spatial grid, achieving better resolution near the center of the domain.     

On the numerical treatment of an ordinary differential equation arising in one-dimensional non-Fickian diffusion problems

In a recent study, Chen and Liu [Comput. Phys. Comm. 150 (2003) 31] considered a one-dimensional, linear non-Fickian diffusion problem with a potential field, which, upon application of the Laplace transform, resulted in a second-order linear ordinary differential equation which was solved by means of a control-volume finite difference method that employs exponential shape functions. It is first shown that this formulation does not properly account for the spatial dependence of the drift forces and results in oscillatory solutions near the left boundary when these forces are large. A piecewise linearized method that provides piecewise analytical solutions, is exact in exact arithmetic for constant coefficients, homogeneous, second-order linear ordinary differential equations and results in three-point finite difference equations is then proposed. Numerical simulations indicate that the piecewise linearized method is free from unphysical oscillations and more accurate than that of Chen and Liu, especially for large drift forces. The method is then applied to non-Fickian diffusion problems with non-constant drift forces in order to determine the effects of the potential field on the concentration distribution.     

Optimized evaluation of a large sum of functions using a three-grid approach

The efficient evaluation of numerous values of functions that vary rapidly only in a small part of the region of interest is presented. An optimized algorithm using a sequence of grids with increasing resolution is developed. The algorithm does not make any assumptions about special properties of the function to be evaluated, e.g., symmetry. An additional speed-up is obtained by exploiting the asymptotic behaviour of the functions to be summed. Two applications from high resolution atmospheric radiative transfer modelling in the infrared and microwave are presented. In a third example asymmetric Rautian line shapes important for high resolution molecular spectroscopy are considered. Computational gains by more than two orders of magnitude with relative errors less than 10⁻³ have been achieved.     

A Mathematica program for the approximate analytical solution to a nonlinear undamped Duffing equation by a new approximate approach

According to Mickens [R.E. Mickens, Comments on a Generalized Galerkin's method for non-linear oscillators, J. Sound Vib. 118 (1987) 563], the general HB (harmonic balance) method is an approximation to the convergent Fourier series representation of the periodic solution of a nonlinear oscillator and not an approximation to an expansion in terms of a small parameter. Consequently, for a nonlinear undamped Duffing equation with a driving force Bcos(ωx), to find a periodic solution when the fundamental frequency is identical to ω, the corresponding Fourier series can be written as

     

Scalable and portable implementation of the fast multipole method on parallel computers

A scalable and portable Fortran code is developed to calculate Coulomb interaction potentials of charged particles on parallel computers, based on the fast multipole method. The code has a unique feature to calculate microscopic stress tensors due to the Coulomb interactions, which is useful in constant-pressure simulations and local stress analyses. The code is applicable to various boundary conditions, including periodic boundary conditions in two and three dimensions, corresponding to slab and bulk systems, respectively. Numerical accuracy of the code is tested through comparison of its results with those obtained by the Ewald summation method and by direct calculations. Scalability tests show the parallel efficiency of 0.98 for 512 million charged particles on 512 IBM SP3 processors. The timing results on IBM SP3 are also compared with those on IBM SP4.     

On time-step bounds in unitary quantum evolution using the Lanczos method

To solve the non-relativistic time dependent Schrödinger equation using the Lanczos method, Park and Light have provided an approximate expression for the time step for a given accuracy. We provide an exact expression for the time step in terms of the eigenvalues and eigenvectors of the resulting tridiagonal matrix. For two test problems, the values of the time step provided by Park and Light differ significantly from the exact values provided by the present method. We also indicate upper and lower bounds for the time step in terms of the maximum and minimum eigenvalues. These bounds indicate the possibility of using a new time step given by the geometric mean of the eigenvalues of the tridiagonal matrix.     

Generalized dynamic programming: A variables freezing method

A variables freezing method is designed for widening the possibilities of numerical application of generalized dynamic programming algorithms. Reducing the size of the memory under certain conditions, the method surmounts the “dimensional curse.” This paper is the continuation of [1–3].     

CANM, a program for numerical solution of a system of nonlinear equations using the continuous analog of Newton's method

A FORTRAN program is presented which solves a system of nonlinear simultaneous equations using the continuous analog of Newton's method (CANM). The user has the option of either to provide a subroutine which calculates the Jacobian matrix or allow the program to calculate it by a forward-difference approximation. Five iterative schemes using different algorithms of determining adaptive step size of the CANM process are implemented in the program.

Program summary

Title of program: CANMCatalogue number: ADSNProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSNProgram available from: CPC Program Library, Queen's University of Belfast, Northern IrelandLicensing provisions: noneComputer for which the program is designed and others on which it has been tested:Computers: IBM RS/6000 Model 320H, SGI Origin2000, SGI Octane, HP 9000/755, Intel Pentium IV PCInstallation: Department of Chemistry, University of Toronto, Toronto, CanadaOperating systems under which the program has been tested: IRIX64 6.1, 6.4 and 6.5, AIX 3.4, HP-UX 9.01, Linux 2.4.7Programming language used: FORTRAN 90Memory required to execute with typical data: depends on the number of nonlinear equations in a system. Test run requires 80 KBNo. of bits in distributed program including test data, etc.: 15283Distribution format: tar gz formatNo. of lines in distributed program, including test data, etc.: 1794Peripherals used: line printer, scratch disc storeExternal subprograms used: DGECO and DGESL [1]Keywords: nonlinear equations, Newton's method, continuous analog of Newton's method, continuous parameter, evolutionary differential equation, Euler's methodNature of physical problem: System of nonlinear simultaneous equations