ODPEVP: A program for computing eigenvalues and eigenfunctions and their first derivatives with respect to the parameter of the parametric self-adjoined Sturm-Liouville problem

A FORTRAN 77 program is presented for calculating with the given accuracy eigenvalues, eigenfunctions and their first derivatives with respect to the parameter of the parametric self-adjoined Sturm-Liouville problem with the parametric third type boundary conditions on the finite interval. The program calculates also potential matrix elements - integrals of the eigenfunctions multiplied by their first derivatives with respect to the parameter. Eigenvalues and matrix elements computed by the ODPEVP program can be used for solving the bound state and multi-channel scattering problems for a system of the coupled second-order ordinary differential equations with the help of the KANTBP programs [O. Chuluunbaatar, A.A. Gusev, A.G. Abrashkevich, A. Amaya-Tapia, M.S. Kaschiev, S.Y. Larsen, S.I. Vinitsky, Comput. Phys. Commun. 177 (2007) 649-675; O. Chuluunbaatar, A.A. Gusev, S.I. Vinitsky, A.G. Abrashkevich, Comput. Phys. Commun. 179 (2008) 685-693]. As a test desk, the program is applied to the calculation of the potential matrix elements for an integrable 2D-model of three identical particles on a line with pair zero-range potentials, a 3D-model of a hydrogen atom in a homogeneous magnetic field and a hydrogen atom on a three-dimensional sphere.

Program summary

Program title: ODPEVPCatalogue identifier: AEDV_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDV_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC license, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 3001No. of bytes in distributed program, including test data, etc.: 24 195Distribution format: tar.gzProgramming language: FORTRAN 77Computer: Intel Xeon EM64T, Alpha 21264A, AMD Athlon MP, Pentium IV Xeon, Opteron 248, Intel Pentium IVOperating system: OC Linux, Unix AIX 5.3, SunOS 5.8, Solaris, Windows XPRAM: depends on

1.

the number and order of finite elements;

2.

the number of points; and

3.

the number of eigenfunctions required.

Test run requires 4 MBClassification: 2.1, 2.4External routines: GAULEG [3]Nature of problem: The three-dimensional boundary problem for the elliptic partial differential equation with an axial symmetry similar to the Schrödinger equation with the Coulomb and transverse oscillator potentials is reduced to the two-dimensional one. The latter finds wide applications in modeling of photoionization and recombination of oppositively charged particles (positrons, antiprotons) in the magnet-optical trap [4], optical absorption in quantum wells [5], and channeling of likely charged particles in thin doped films [6,7] or neutral atoms and molecules in artificial waveguides or surfaces [8,9]. In the adiabatic approach [10] known in mathematics as Kantorovich method [11] the solution of the two-dimensional elliptic partial differential equation is expanded over basis functions with respect to the fast variable (for example, angular variable) and depended on the slow variable (for example, radial coordinate ) as a parameter. An averaging of the problem by such a basis leads to a system of the second-order ordinary differential equations which contain potential matrix elements and the first-derivative coupling terms (see, e.g., [12,13,14]). The purpose of this paper is to present the finite element method procedure based on the use of high-order accuracy approximations for calculating eigenvalues, eigenfunctions and their first derivatives with respect to the parameter of the parametric self-adjoined Sturm-Liouville problem with the parametric third type boundary conditions on the finite interval. The program developed calculates potential matrix elements - integrals of the eigenfunctions multiplied by their derivatives with respect to the parameter. These matrix elements can be used for solving the bound state and multi-channel scattering problems for a system of the coupled second-order ordinary differential equations with the help of the KANTBP programs [1,2].Solution method: The parametric self-adjoined Sturm-Liouville problem with the parametric third type boundary conditions is solved by the finite element method using high-order accuracy approximations [15]. The generalized algebraic eigenvalue problem AF=EBF with respect to a pair of unknown (E,F) arising after the replacement of the differential problem by the finite-element approximation is solved by the subspace iteration method using the SSPACE program [16]. First derivatives of the eigenfunctions with respect to the parameter which contained in potential matrix elements of the coupled system equations are obtained by solving the inhomogeneous algebraic equations. As a test desk, the program is applied to the calculation of the potential matrix elements for an integrable 2D-model of three identical particles on a line with pair zero-range potentials described in [1,17,18], a 3D-model of a hydrogen atom in a homogeneous magnetic field described in [14,19] and a hydrogen atom on a three-dimensional sphere [20].Restrictions: The computer memory requirements depend on:

1.

the number and order of finite elements;

2.

the number of points; and

3.

the number of eigenfunctions required.

Restrictions due to dimension sizes may be easily alleviated by altering PARAMETER statements (see sections below and listing for details). The user must also supply DOUBLE PRECISION functions POTCCL and POTCC1 for evaluating potential function U(ρ,z) of Eq. (1) and its first derivative with respect to parameter ρ. The user should supply DOUBLE PRECISION functions F1FUNC and F2FUNC that evaluate functions f₁(z) and f₂(z) of Eq. (1). The user must also supply subroutine BOUNCF for evaluating the parametric third type boundary conditions.Running time: The running time depends critically upon:

1.

the number and order of finite elements;

2.

the number of points on interval [zmin,zmax]; and

3.

the number of eigenfunctions required.

The test run which accompanies this paper took 2 s with calculation of matrix potentials on the Intel Pentium IV 2.4 GHz.References:

[1]

O. Chuluunbaatar, A.A. Gusev, A.G. Abrashkevich, A. Amaya-Tapia, M.S. Kaschiev, S.Y. Larsen, S.I. Vinitsky, Comput. Phys. Comm. 177 (2007) 649-675

[2]

O. Chuluunbaatar, A.A. Gusev, S.I. Vinitsky, A.G. Abrashkevich, Comput. Phys. Comm. 179 (2008) 685-693.

[3]

W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical Recipes: The Art of Scientific Computing, Cambridge University Press, Cambridge, 1986.

[4]

O. Chuluunbaatar, A.A. Gusev, S.I. Vinitsky, V.L. Derbov, L.A. Melnikov, V.V. Serov, Phys. Rev. A 77 (2008) 034702-1-4.

[5]

E.M. Kazaryan, A.A. Kostanyan, H.A. Sarkisyan, Physica E 28 (2005) 423-430.

[6]

Yu.N. Demkov, J.D. Meyer, Eur. Phys. J. B 42 (2004) 361-365.

[7]

P.M. Krassovitskiy, N.Zh. Takibaev, Bull. Russian Acad. Sci. Phys. 70 (2006) 815-818.

[8]

V.S. Melezhik, J.I. Kim, P. Schmelcher, Phys. Rev. A 76 (2007) 053611-1-15.

[9]

F.M. Pen'kov, Phys. Rev. A 62 (2000) 044701-1-4.

[10]

M. Born, X. Huang, Dynamical Theory of Crystal Lattices, The Clarendon Press, Oxford, England, 1954.

[11]

L.V. Kantorovich, V.I. Krylov, Approximate Methods of Higher Analysis, Wiley, New York, 1964.

[12]

U. Fano, Colloq. Int. C.N.R.S. 273 (1977) 127;

A.F. Starace, G.L. Webster, Phys. Rev. A 19 (1979) 1629-1640.

[13]

C.V. Clark, K.T. Lu, A.F. Starace, in: H.G. Beyer, H. Kleinpoppen (eds.), Progress in Atomic Spectroscopy, Part C, Plenum, New York, 1984, pp. 247-320.

[14]

O. Chuluunbaatar, A.A. Gusev, V.L. Derbov, M.S. Kaschiev, L.A. Melnikov, V.V. Serov, S.I. Vinitsky, J. Phys. A 40 (2007) 11485-11524.

[15]

A.G. Abrashkevich, D.G. Abrashkevich, M.S. Kaschiev, I.V. Puzynin, Comput. Phys. Comm. 85 (1995) 40-64.

[16]

K.J. Bathe, Finite Element Procedures in Engineering Analysis, Englewood Cliffs, Prentice-Hall, New York, 1982.

[17]

O. Chuluunbaatar, A.A. Gusev, M.S. Kaschiev, V.A. Kaschieva, A. Amaya-Tapia, S.Y. Larsen, S.I. Vinitsky, J. Phys. B 39 (2006) 243-269.

[18]

Yu.A. Kuperin, P.B. Kurasov, Yu.B. Melnikov, S.P. Merkuriev, Ann. Phys. 205 (1991) 330-361.

[19]

O. Chuluunbaatar, A.A. Gusev, V.P. Gerdt, V.A. Rostovtsev, S.I. Vinitsky, A.G. Abrashkevich, M.S. Kaschiev, V.V. Serov, Comput. Phys. Comm. 178 (2008) 301-330.

[20]

A.G. Abrashkevich, M.S. Kaschiev, S.I. Vinitsky, J. Comp. Phys. 163 (2000) 328-348.

     

KANTBP 2.0: New version of a program for computing energy levels, reaction matrix and radial wave functions in the coupled-channel hyperspherical adiabatic approach

A FORTRAN 77 program for calculating energy values, reaction matrix and corresponding radial wave functions in a coupled-channel approximation of the hyperspherical adiabatic approach is presented. In this approach, a multi-dimensional Schrödinger equation is reduced to a system of the coupled second-order ordinary differential equations on a finite interval with homogeneous boundary conditions: (i) the Dirichlet, Neumann and third type at the left and right boundary points for continuous spectrum problem, (ii) the Dirichlet and Neumann type conditions at left boundary point and Dirichlet, Neumann and third type at the right boundary point for the discrete spectrum problem. The resulting system of radial equations containing the potential matrix elements and first-derivative coupling terms is solved using high-order accuracy approximations of the finite element method. As a test desk, the program is applied to the calculation of the reaction matrix and radial wave functions for 3D-model of a hydrogen-like atom in a homogeneous magnetic field. This version extends the previous version 1.0 of the KANTBP program [O. Chuluunbaatar, A.A. Gusev, A.G. Abrashkevich, A. Amaya-Tapia, M.S. Kaschiev, S.Y. Larsen, S.I. Vinitsky, Comput. Phys. Commun. 177 (2007) 649-675].

Program summary

Program title: KANTBPCatalogue identifier: ADZH_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZH_v2_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.: 20 403No. of bytes in distributed program, including test data, etc.: 147 563Distribution format: tar.gzProgramming language: FORTRAN 77Computer: Intel Xeon EM64T, Alpha 21264A, AMD Athlon MP, Pentium IV Xeon, Opteron 248, Intel Pentium IVOperating system: OC Linux, Unix AIX 5.3, SunOS 5.8, Solaris, Windows XPRAM: This depends on

1.

the number of differential equations;

2.

the number and order of finite elements;

3.

the number of hyperradial points; and

4.

the number of eigensolutions required.

The test run requires 2 MBClassification: 2.1, 2.4External routines: GAULEG and GAUSSJ [2]Nature of problem: In the hyperspherical adiabatic approach [3-5], a multidimensional Schrödinger equation for a two-electron system [6] or a hydrogen atom in magnetic field [7-9] is reduced by separating radial coordinate ρ from the angular variables to a system of the second-order ordinary differential equations containing the potential matrix elements and first-derivative coupling terms. The purpose of this paper is to present the finite element method procedure based on the use of high-order accuracy approximations for calculating approximate eigensolutions of the continuum spectrum for such systems of coupled differential equations on finite intervals of the radial variable ρ∈[ρmin,ρmax]. This approach can be used in the calculations of effects of electron screening on low-energy fusion cross sections [10-12].Solution method: The boundary problems for the coupled second-order differential equations are solved by the finite element method using high-order accuracy approximations [13]. The generalized algebraic eigenvalue problem AF=EBF with respect to pair unknowns (E,F) arising after the replacement of the differential problem by the finite-element approximation is solved by the subspace iteration method using the SSPACE program [14]. The generalized algebraic eigenvalue problem (A−EB)F=λDF with respect to pair unknowns (λ,F) arising after the corresponding replacement of the scattering boundary problem in open channels at fixed energy value, E, is solved by the LDLT factorization of symmetric matrix and back-substitution methods using the DECOMP and REDBAK programs, respectively [14]. As a test desk, the program is applied to the calculation of the reaction matrix and corresponding radial wave functions for 3D-model of a hydrogen-like atom in a homogeneous magnetic field described in [9] on finite intervals of the radial variable ρ∈[ρmin,ρmax]. For this benchmark model the required analytical expressions for asymptotics of the potential matrix elements and first-derivative coupling terms, and also asymptotics of radial solutions of the boundary problems for coupled differential equations have been produced with help of a MAPLE computer algebra system.Restrictions: The computer memory requirements depend on:

1.

the number of differential equations;

2.

the number and order of finite elements;

3.

the total number of hyperradial points; and

4.

the number of eigensolutions required.

Restrictions due to dimension sizes may be easily alleviated by altering PARAMETER statements (see Section 3 and [1] for details). The user must also supply subroutine POTCAL for evaluating potential matrix elements. The user should also supply subroutines ASYMEV (when solving the eigenvalue problem) or ASYMS0 and ASYMSC (when solving the scattering problem) which evaluate asymptotics of the radial wave functions at left and right boundary points in case of a boundary condition of the third type for the above problems.Running time: The running time depends critically upon:

1.

the number of differential equations;

2.

the number and order of finite elements;

3.

the total number of hyperradial points on interval [ρmin,ρmax]; and

4.

the number of eigensolutions required.

The test run which accompanies this paper took 2 s without calculation of matrix potentials on the Intel Pentium IV 2.4 GHz.References:[1] O. Chuluunbaatar, A.A. Gusev, A.G. Abrashkevich, A. Amaya-Tapia, M.S. Kaschiev, S.Y. Larsen, S.I. Vinitsky, Comput. Phys. Commun. 177 (2007) 649-675; http://cpc.cs.qub.ac.uk/summaries/ADZHv10.html.[2] W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical Recipes: The Art of Scientific Computing, Cambridge University Press, Cambridge, 1986.[3] J. Macek, J. Phys. B 1 (1968) 831-843.[4] U. Fano, Rep. Progr. Phys. 46 (1983) 97-165.[5] C.D. Lin, Adv. Atom. Mol. Phys. 22 (1986) 77-142.[6] A.G. Abrashkevich, D.G. Abrashkevich, M. Shapiro, Comput. Phys. Commun. 90 (1995) 311-339.[7] M.G. Dimova, M.S. Kaschiev, S.I. Vinitsky, J. Phys. B 38 (2005) 2337-2352.[8] O. Chuluunbaatar, A.A. Gusev, V.L. Derbov, M.S. Kaschiev, L.A. Melnikov, V.V. Serov, S.I. Vinitsky, J. Phys. A 40 (2007) 11485-11524.[9] O. Chuluunbaatar, A.A. Gusev, V.P. Gerdt, V.A. Rostovtsev, S.I. Vinitsky, A.G. Abrashkevich, M.S. Kaschiev, V.V. Serov, Comput. Phys. Commun. 178 (2007) 301 330; http://cpc.cs.qub.ac.uk/summaries/AEAAv10.html.[10] H.J. Assenbaum, K. Langanke, C. Rolfs, Z. Phys. A 327 (1987) 461-468.[11] V. Melezhik, Nucl. Phys. A 550 (1992) 223-234.[12] L. Bracci, G. Fiorentini, V.S. Melezhik, G. Mezzorani, P. Pasini, Phys. Lett. A 153 (1991) 456-460.[13] A.G. Abrashkevich, D.G. Abrashkevich, M.S. Kaschiev, I.V. Puzynin, Comput. Phys. Commun. 85 (1995) 40-64.[14] K.J. Bathe, Finite Element Procedures in Engineering Analysis, Englewood Cliffs, Prentice-Hall, New York, 1982.     

POTHMF: A program for computing potential curves and matrix elements of the coupled adiabatic radial equations for a hydrogen-like atom in a homogeneous magnetic field

A FORTRAN 77 program is presented which calculates with the relative machine precision potential curves and matrix elements of the coupled adiabatic radial equations for a hydrogen-like atom in a homogeneous magnetic field. The potential curves are eigenvalues corresponding to the angular oblate spheroidal functions that compose adiabatic basis which depends on the radial variable as a parameter. The matrix elements of radial coupling are integrals in angular variables of the following two types: product of angular functions and the first derivative of angular functions in parameter, and product of the first derivatives of angular functions in parameter, respectively. The program calculates also the angular part of the dipole transition matrix elements (in the length form) expressed as integrals in angular variables involving product of a dipole operator and angular functions. Moreover, the program calculates asymptotic regular and irregular matrix solutions of the coupled adiabatic radial equations at the end of interval in radial variable needed for solving a multi-channel scattering problem by the generalized R-matrix method. Potential curves and radial matrix elements computed by the POTHMF program can be used for solving the bound state and multi-channel scattering problems. As a test desk, the program is applied to the calculation of the energy values, a short-range reaction matrix and corresponding wave functions with the help of the KANTBP program. Benchmark calculations for the known photoionization cross-sections are presented.

Program summary

Program title:POTHMFCatalogue identifier:AEAA_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAA_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.:8123No. of bytes in distributed program, including test data, etc.:131 396Distribution format:tar.gzProgramming language:FORTRAN 77Computer:Intel Xeon EM64T, Alpha 21264A, AMD Athlon MP, Pentium IV Xeon, Opteron 248, Intel Pentium IVOperating system:OC Linux, Unix AIX 5.3, SunOS 5.8, Solaris, Windows XPRAM:Depends on

1.

the number of radial differential equations;

2.

the number and order of finite elements;

3.

the number of radial points.

Test run requires 4 MBClassification:2.5External routines:POTHMF uses some Lapack routines, copies of which are included in the distribution (see README file for details).Nature of problem:In the multi-channel adiabatic approach the Schrödinger equation for a hydrogen-like atom in a homogeneous magnetic field of strength γ (γ=B/B₀, B₀≅2.35×¹⁰⁵ T is a dimensionless parameter which determines the field strength B) is reduced by separating the radial coordinate, r, from the angular variables, (θ,φ), and using a basis of the angular oblate spheroidal functions [3] to a system of second-order ordinary differential equations which contain first-derivative coupling terms [4]. The purpose of this program is to calculate potential curves and matrix elements of radial coupling needed for calculating the low-lying bound and scattering states of hydrogen-like atoms in a homogeneous magnetic field of strength 0<γ?1000 within the adiabatic approach [5]. The program evaluates also asymptotic regular and irregular matrix radial solutions of the multi-channel scattering problem needed to extract from the R-matrix a required symmetric shortrange open-channel reaction matrix K [6] independent from matching point [7]. In addition, the program computes the dipole transition matrix elements in the length form between the basis functions that are needed for calculating the dipole transitions between the low-lying bound and scattering states and photoionization cross sections [8].Solution method:The angular oblate spheroidal eigenvalue problem depending on the radial variable is solved using a series expansion in the Legendre polynomials [3]. The resulting tridiagonal symmetric algebraic eigenvalue problem for the evaluation of selected eigenvalues, i.e. the potential curves, is solved by the LDL^T factorization using the DSTEVR program [2]. Derivatives of the eigenfunctions with respect to the radial variable which are contained in matrix elements of the coupled radial equations are obtained by solving the inhomogeneous algebraic equations. The corresponding algebraic problem is solved by using the LDL^T factorization with the help of the DPTTRS program [2]. Asymptotics of the matrix elements at large values of radial variable are computed using a series expansion in the associated Laguerre polynomials [9]. The corresponding matching points between the numeric and asymptotic solutions are found automatically. These asymptotics are used for the evaluation of the asymptotic regular and irregular matrix radial solutions of the multi-channel scattering problem [7]. As a test desk, the program is applied to the calculation of the energy values of the ground and excited bound states and reaction matrix of multi-channel scattering problem for a hydrogen atom in a homogeneous magnetic field using the KANTBP program [10].Restrictions:The computer memory requirements depend on:

1.

the number of radial differential equations;

2.

the number and order of finite elements;

3.

the total number of radial points.

Restrictions due to dimension sizes can be changed by resetting a small number of PARAMETER statements before recompiling (see Introduction and listing for details).Running time:The running time depends critically upon:

1.

the number of radial differential equations;

2.

the number and order of finite elements;

3.

the total number of radial points on interval [rmin,rmax].

The test run which accompanies this paper took 7 s required for calculating of potential curves, radial matrix elements, and dipole transition matrix elements on a finite-element grid on interval [rmin=0, rmax=100] used for solving discrete and continuous spectrum problems and obtaining asymptotic regular and irregular matrix radial solutions at rmax=100 for continuous spectrum problem on the Intel Pentium IV 2.4 GHz. The number of radial differential equations was equal to 6. The accompanying test run using the KANTBP program took 2 s for solving discrete and continuous spectrum problems using the above calculated potential curves, matrix elements and asymptotic regular and irregular matrix radial solutions. Note, that in the accompanied benchmark calculations of the photoionization cross-sections from the bound states of a hydrogen atom in a homogeneous magnetic field to continuum we have used interval [rmin=0, rmax=1000] for continuous spectrum problem. The total number of radial differential equations was varied from 10 to 18.References:

[1]

W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical Recipes: The Art of Scientific Computing, Cambridge University Press, Cambridge, 1986.

[2]

http://www.netlib.org/lapack/.

[3]

M. Abramovits, I.A. Stegun, Handbook of Mathematical Functions, Dover, New York, 1965.

[4]

U. Fano, Colloq. Int. C.N.R.S. 273 (1977) 127; A.F. Starace, G.L. Webster, Phys. Rev. A 19 (1979) 1629-1640; C.V. Clark, K.T. Lu, A.F. Starace, in: H.G. Beyer, H. Kleinpoppen (Eds.), Progress in Atomic Spectroscopy, Part C, Plenum, New York, 1984, pp. 247-320; U. Fano, A.R.P. Rau, Atomic Collisions and Spectra, Academic Press, Florida, 1986.

[5]

M.G. Dimova, M.S. Kaschiev, S.I. Vinitsky, J. Phys. B 38 (2005) 2337-2352; O. Chuluunbaatar, A.A. Gusev, V.L. Derbov, M.S. Kaschiev, V.V. Serov, T.V. Tupikova, S.I. Vinitsky, Proc. SPIE 6537 (2007) 653706-1-18.

[6]

M.J. Seaton, Rep. Prog. Phys. 46 (1983) 167-257.

[7]

M. Gailitis, J. Phys. B 9 (1976) 843-854; J. Macek, Phys. Rev. A 30 (1984) 1277-1278; S.I. Vinitsky, V.P. Gerdt, A.A. Gusev, M.S. Kaschiev, V.A. Rostovtsev, V.N. Samoylov, T.V. Tupikova, O. Chuluunbaatar, Programming and Computer Software 33 (2007) 105-116.

[8]

H. Friedrich, Theoretical Atomic Physics, Springer, New York, 1991.

[9]

R.J. Damburg, R.Kh. Propin, J. Phys. B 1 (1968) 681-691; J.D. Power, Phil. Trans. Roy. Soc. London A 274 (1973) 663-702.

[10]

O. Chuluunbaatar, A.A. Gusev, A.G. Abrashkevich, A. Amaya-Tapia, M.S. Kaschiev, S.Y. Larsen, S.I. Vinitsky, Comput. Phys. Comm. 177 (2007) 649-675.

     

A new version of Scilab software package for the study of dynamical systems

This work presents a new version of a software package for the study of chaotic flows, maps and fractals [1]. The codes were written using Scilab, a software package for numerical computations providing a powerful open computing environment for engineering and scientific applications. It was found that Scilab provides various functions for ordinary differential equation solving, Fast Fourier Transform, autocorrelation, and excellent 2D and 3D graphical capabilities. The chaotic behaviors of the nonlinear dynamics systems were analyzed using phase-space maps, autocorrelation functions, power spectra, Lyapunov exponents and Kolmogorov-Sinai entropy. Various well-known examples are implemented, with the capability of the users inserting their own ODE or iterative equations.

New version program summary

Program title: Chaos v2.0Catalogue identifier: AEAP_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAP_v2_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.: 1275No. of bytes in distributed program, including test data, etc.: 7135Distribution format: tar.gzProgramming language: Scilab 5.1.1. Scilab 5.1.1 should be installed before running the program. Information about the installation can be found at http://wiki.scilab.org/howto/install/windows.Computer: PC-compatible running Scilab on MS Windows or LinuxOperating system: Windows XP, LinuxRAM: below 150 MegabytesClassification: 6.2Catalogue identifier of previous version: AEAP_v1_0Journal reference of previous version: Comput. Phys. Comm. 178 (2008) 788Does the new version supersede the previous version?: YesNature of problem: Any physical model containing linear or nonlinear ordinary differential equations (ODE).Solution method:

1.

Numerical solving of ordinary differential equations for the study of chaotic flows. The chaotic behavior of the nonlinear dynamical system is analyzed using Poincare sections, phase-space maps, autocorrelation functions, power spectra, Lyapunov exponents and Kolmogorov-Sinai entropies.

2.

Numerical solving of iterative equations for the study of maps and fractals.

Reasons for new version: The program has been updated to use the new version 5.1.1 of Scilab with new graphical capabilities [2]. Moreover, new use cases have been added which make the handling of the program easier and more efficient.Summary of revisions:

1.

A new use case concerning coupled predator-prey models has been added [3].

2.

Three new use cases concerning fractals (Sierpinsky gasket, Barnsley's Fern and Tree) have been added [3].

3.

The graphical user interface (GUI) of the program has been reconstructed to include the new use cases.

4.

The program has been updated to use Scilab 5.1.1 with the new graphical capabilities.

Additional comments: The program package contains 12 subprograms.

•

interface.sce - the graphical user interface (GUI) that permits the choice of a routine as follows

•

1.sci - Lorenz dynamical system

•

2.sci - Chua dynamical system

•

3.sci - Rosler dynamical system

•

4.sci - Henon map

•

5.sci - Lyapunov exponents for Lorenz dynamical system

•

6.sci - Lyapunov exponent for the logistic map

•

7.sci - Shannon entropy for the logistic map

•

8.sci - Coupled predator-prey model

•

1f.sci - Sierpinsky gasket

•

2f.sci - Barnsley's Fern

•

3f.sci - Barnsley's Tree

Running time: 10 to 20 seconds for problems that do not involve Lyapunov exponents calculation; 60 to 1000 seconds for problems that involve high orders ODE, Lyapunov exponents calculation and fractals.References:

[1]

C.C. Bordeianu, C. Besliu, Al. Jipa, D. Felea, I. V. Grossu, Comput. Phys. Comm. 178 (2008) 788.

[2]

S. Campbell, J.P. Chancelier, R. Nikoukhah, Modeling and Simulation in Scilab/Scicos, Springer, 2006.

[3]

R.H. Landau, M.J. Paez, C.C. Bordeianu, A Survey of Computational Physics, Introductory Computational Science, Princeton University Press, 2008.

     

JADAMILU: a software code for computing selected eigenvalues of large sparse symmetric matrices

A new software code for computing selected eigenvalues and associated eigenvectors of a real symmetric matrix is described. The eigenvalues are either the smallest or those closest to some specified target, which may be in the interior of the spectrum. The underlying algorithm combines the Jacobi-Davidson method with efficient multilevel incomplete LU (ILU) preconditioning. Key features are modest memory requirements and robust convergence to accurate solutions. Parameters needed for incomplete LU preconditioning are automatically computed and may be updated at run time depending on the convergence pattern. The software is easy to use by non-experts and its top level routines are written in FORTRAN 77. Its potentialities are demonstrated on a few applications taken from computational physics.

Program summary

Program title: JADAMILUCatalogue identifier: ADZT_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZT_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.: 101 359No. of bytes in distributed program, including test data, etc.: 7 493 144Distribution format: tar.gzProgramming language: Fortran 77Computer: Intel or AMD with g77 and pgf; Intel EM64T or Itanium with ifort; AMD Opteron with g77, pgf and ifort; Power (IBM) with xlf90.Operating system: Linux, AIXRAM: problem dependentWord size: real:8; integer: 4 or 8, according to user's choiceClassification: 4.8Nature of problem: Any physical problem requiring the computation of a few eigenvalues of a symmetric matrix.Solution method: Jacobi-Davidson combined with multilevel ILU preconditioning.Additional comments: We supply binaries rather than source code because JADAMILU uses the following external packages:

•

MC64. This software is copyrighted software and not freely available. COPYRIGHT (c) 1999 Council for the Central Laboratory of the Research Councils.

•

AMD. Copyright (c) 2004-2006 by Timothy A. Davis, Patrick R. Amestoy, and Iain S. Duff. All Rights Reserved. Source code is distributed by the authors under the GNU LGPL licence.

•

BLAS. The reference BLAS is a freely-available software package. It is available from netlib via anonymous ftp and the World Wide Web.

•

LAPACK. The complete LAPACK package or individual routines from LAPACK are freely available on netlib and can be obtained via the World Wide Web or anonymous ftp.

•

For maximal benefit to the community, we added the sources we are proprietary of to the tar.gz file submitted for inclusion in the CPC library. However, as explained in the README file, users willing to compile the code instead of using binaries should first obtain the sources for the external packages mentioned above (email and/or web addresses are provided).

Running time: Problem dependent; the test examples provided with the code only take a few seconds to run; timing results for large scale problems are given in Section 5.     

Multithreaded transactions in scientific computing. The Growth06_v2 program

Writing a concurrent program can be more difficult than writing a sequential program. Programmer needs to think about synchronization, race conditions and shared variables. Transactions help reduce the inconvenience of using threads. A transaction is an abstraction, which allows programmers to group a sequence of actions on the program into a logical, higher-level computation unit. This paper presents a new version of the GROWTHGr and GROWTH06 programs.

New version program summary

Program title: GROWTH06_v2Catalogue identifier: ADVL_v2_1Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVL_v2_1.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.: 65 255No. of bytes in distributed program, including test data, etc.: 865 985Distribution format: tar.gzProgramming language: Object PascalComputer: Pentium-based PCOperating system: Windows 9x, XP, NT, VistaRAM: more than 1 MBClassification: 4.3, 7.2, 6.2, 8, 14Catalogue identifier of previous version: ADVL_v2_0Journal reference of previous version: Comput. Phys. Comm. 175 (2006) 678Does the new version supersede the previous version?: YesNature of problem: The programs compute the RHEED intensities during the growth of thin epitaxial structures prepared using the molecular beam epitaxy (MBE). The computations are based on the use of kinematical diffraction theory.Solution method: Epitaxial growth of thin films is modelled by a set of non-linear differential equations [1]. The Runge-Kutta method with adaptive stepsize control was used for solving initial value problem for non-linear differential equations [2].Reasons for new version: According to the users' suggestions functionality of the program has been improved. Moreover, new use cases have been added which make the handling of the program easier and more efficient than the previous ones [3].Summary of revisions:

1.

The design pattern (See Fig. 2 of Ref. [3]) has been modified according to the scheme shown on Fig. 1.

2.

A graphical user interface (GUI) for the program has been reconstructed. Fig. 2 presents a hybrid diagram of a GUI that shows how onscreen objects connect to use cases.

3.

The program has been compiled with English/USA regional and language options.

Note: The figures mentioned above are contained in the program distribution file.Unusual features: The program is distributed in the form of source project GROWTH06_v2.dpr with associated files, and should be compiled using Borland Delphi compilers versions 6 or latter (including Borland Developer Studio 2006 and Code Gear compilers for Delphi).Additional comments: Two figures are included in the program distribution file. These are captioned

1.

Static classes model for Transaction design pattern.

2.

A model of a window that shows how onscreen objects connect to use cases.

Running time: The typical running time is machine and user-parameters dependent.References: [1] A. Daniluk, Comput. Phys. Comm. 170 (2005) 265.[2] W.H. Press, B.P. Flannery, S.A. Teukolsky, W.T. Vetterling, Numerical Recipes in Pascal: The Art of Scientific Computing, first ed., Cambridge University Press, 1989.[3] M. Brzuszek, A. Daniluk, Comput. Phys. Comm. 175 (2006) 678.     

LCG MCDB—a knowledgebase of Monte-Carlo simulated events

In this paper we report on LCG Monte-Carlo Data Base (MCDB) and software which has been developed to operate MCDB. The main purpose of the LCG MCDB project is to provide a storage and documentation system for sophisticated event samples simulated for the LHC Collaborations by experts. In many cases, the modern Monte-Carlo simulation of physical processes requires expert knowledge in Monte-Carlo generators or significant amount of CPU time to produce the events. MCDB is a knowledgebase mainly dedicated to accumulate simulated events of this type. The main motivation behind LCG MCDB is to make the sophisticated MC event samples available for various physical groups. All the data from MCDB is accessible in several convenient ways. LCG MCDB is being developed within the CERN LCG Application Area Simulation project.

Program summary

Program title: LCG Monte-Carlo Data BaseCatalogue identifier: ADZX_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZX_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU General Public LicenceNo. of lines in distributed program, including test data, etc.: 30 129No. of bytes in distributed program, including test data, etc.: 216 943Distribution format: tar.gzProgramming language: PerlComputer: CPU: Intel Pentium 4, RAM: 1 Gb, HDD: 100 GbOperating system: Scientific Linux CERN 3/4RAM: 1 073 741 824 bytes (1 Gb)Classification: 9External routines:

•

perl >= 5.8.5;

•

Perl modules

○

DBD-mysql >= 2.9004,

○

File::Basename,

○

GD::SecurityImage,

○

GD::SecurityImage::AC,

○

Linux::Statistics,

○

XML::LibXML > 1.6,

○

XML::SAX,

○

XML::NamespaceSupport;

•

Apache HTTP Server >= 2.0.59;

•

mod auth external >= 2.2.9;

•

edg-utils-system RPM package;

•

gd >= 2.0.28;

•

rpm package CASTOR-client >= 2.1.2-4;

•

arc-server (optional)

Nature of problem: Often, different groups of experimentalists prepare similar samples of particle collision events or turn to the same group of authors of Monte-Carlo (MC) generators to prepare the events. For example, the same MC samples of Standard Model (SM) processes can be employed for the investigations either in the SM analyses (as a signal) or in searches for new phenomena in Beyond Standard Model analyses (as a background). If the samples are made available publicly and equipped with corresponding and comprehensive documentation, it can speed up cross checks of the samples themselves and physical models applied. Some event samples require a lot of computing resources for preparation. So, a central storage of the samples prevents possible waste of researcher time and computing resources, which can be used to prepare the same events many times.Solution method: Creation of a special knowledgebase (MCDB) designed to keep event samples for the LHC experimental and phenomenological community. The knowledgebase is realized as a separate web-server (http://mcdb.cern.ch). All event samples are kept on types at CERN. Documentation describing the events is the main contents of MCDB. Users can browse the knowledgebase, read and comment articles (documentation), and download event samples. Authors can upload new event samples, create new articles, and edit own articles.Restrictions: The software is adopted to solve the problems, described in the article and there are no any additional restrictions.Unusual features: The software provides a framework to store and document large files with flexible authentication and authorization system. Different external storages with large capacity can be used to keep the files. The WEB Content Management System provides all of the necessary interfaces for the authors of the files, end-users and administrators.Running time: Real time operations.References:[1] The main LCG MCDB server, http://mcdb.cern.ch/.[2] P. Bartalini, L. Dudko, A. Kryukov, I.V. Selyuzhenkov, A. Sherstnev, A. Vologdin, LCG Monte-Carlo data base, hep-ph/0404241.[3] J.P. Baud, B. Couturier, C. Curran, J.D. Durand, E. Knezo, S. Occhetti, O. Barring, CASTOR: status and evolution, cs.oh/0305047.     

GeodesicViewer - A tool for exploring geodesics in the theory of relativity

The GeodesicViewer realizes exocentric two- and three-dimensional illustrations of lightlike and timelike geodesics in the general theory of relativity. By means of an intuitive graphical user interface, all parameters of a spacetime as well as the initial conditions of the geodesics can be modified interactively. This makes the GeodesicViewer a useful instrument for the exploration of geodesics in four-dimensional Lorentzian spacetimes.

Program summary

Program title: GeodesicViewerCatalogue identifier: AEFP_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFP_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.: 168 868No. of bytes in distributed program, including test data, etc.: 6 076 202Distribution format: tar.gzProgramming language: C++, Qt, Qwt, OpenGLComputer: All platforms with a C++ compiler, Qt, Qwt, OpenGLOperating system: Linux, Mac OS XRAM: 24 MbytesClassification: 1.5External routines:

•

Gnu Scientific Library (GSL) (http://www.gnu.org/software/gsl/)

•

Motion4D (included in the package). The Motion4D library can also be downloaded from CPC. Catalogue identifier: AEEX

•

Qt (http://qt.nokia.com/downloads)

•

Qwt (http://qwt.sourceforge.net/)

•

OpenGL (http://www.opengl.org/)

Nature of problem: Illustrate geodesics in four-dimensional Lorentzian spacetimes.Solution method: Integration of ordinary differential equations. 3D-Rendering via OpenGL.Running time: Interactive. The examples given take milliseconds.     

micrOMEGAs 2.0.7: a program to calculate the relic density of dark matter in a generic model

micrOMEGAs2.0.7 is a code which calculates the relic density of a stable massive particle in an arbitrary model. The underlying assumption is that there is a conservation law like R-parity in supersymmetry which guarantees the stability of the lightest odd particle. The new physics model must be incorporated in the notation of CalcHEP, a package for the automatic generation of squared matrix elements. Once this is done, all annihilation and coannihilation channels are included automatically in any model. Cross-sections at v=0, relevant for indirect detection of dark matter, are also computed automatically. The package includes three sample models: the minimal supersymmetric standard model (MSSM), the MSSM with complex phases and the NMSSM. Extension to other models, including non supersymmetric models, is described.

Program summary

Title of program:micrOMEGAs2.0.7Catalogue identifier:ADQR_v2_1Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADQR_v2_1.htmlProgram obtainable from: CPC Program Library, Queen's University of 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.:216 529No. of bytes in distributed program, including test data, etc.:1 848 816Distribution format:tar.gzProgramming language used:C and FortranComputer:PC, Alpha, Mac, SunOperating system:UNIX (Linux, OSF1, SunOS, Darwin, Cygwin)RAM:17 MB depending on the number of processes requiredClassification:1.9, 11.6Catalogue identifier of previous version:ADQR_v2_0Journal version of previous version:Comput. Phys. Comm. 176 (2007) 367Does the new version supersede the previous version?:YesNature of problem:Calculation of the relic density of the lightest stable particle in a generic new model of particle physics.Solution method:In numerically solving the evolution equation for the density of dark matter, relativistic formulae for the thermal average are used. All tree-level processes for annihilation and coannihilation of new particles in the model are included. The cross-sections for all processes are calculated exactly with CalcHEP after definition of a model file. Higher-order QCD corrections to Higgs couplings to quark pairs are included.Reasons for new version:The main changes in this new version consist, on the one hand, in improvements of the user interface and treatment of error codes when using spectrum calculators in the MSSM and, on the other hand, on a completely revised code for the calculation of the relic density in the NMSSM based on the code NMSSMTools1.0.2 for the computation of the spectrum.Summary of revisions:

•

The version of CalcHEP was updated to CalcHEP 2.4.

•

The procedure for shared library generation has been improved. Now the libraries are recalculated each time the model is modified.

•

The default value for the top quark mass has been set to 171.4 GeV.

Changes specific to the MSSM model.

•

The deltaMb correction is now included in the B,t,H-vertex and is always included for other Higgs vertices.

•

In case of a fatal error in an RGE program, micrOMEGAs now continues operation while issuing a warning that the given point is not valid. This is important when running scans over parameter space. However this means that the standard ˆC command that could be used to cancel a job now only cancels the RGE program. To cancel a job, use “kill -9 -N” where N is the micrOMEGAs process id, all child processes launched by micrOMEGAs will be killed at once.

•

Following the last SLHA2 release, we use key=26 item of EXTPAR block for the pole mass of the CP-odd Higgs so that micrOMEGAs can now use SoftSUSY for spectrum calculation with EWSB input. The Isajet interface was corrected too, so the user has to recompile the isajet_slha executable. For SuSpect we still support an old “wrong” interface where key=24 is used for the mass of the CP-odd Higgs.

•

In the non-universal SUGRA model, we set the value of M₀ (M1/2,A₀) to the value of the largest subset of equal parameters among scalar masses (gaugino masses, trilinear couplings). In the previous version these parameters were set arbitrarily to be equal to MH2, MG2 and At respectively. The spectrum calculators need an input value for M₀,M1/2 and A₀ for initialisation purposes.

•

We have removed bugs in micrOMEGAs-Isajet interface in case of non-universal SUGRA.

•

$(FFLAGS) is added to compilation instruction of suspect.exe. It was omitted in version 2.0.

•

The treatment of errors in reading of the LesHouches accord file is improved. Now, if the SPINFO block is absent in the SLHA output it is considered as a fatal error.

•

Instructions for calculation of Δρ, μ(g−2), Br(b→sγ) and Br(Bs→μ⁺μ⁻) constraints are included in EWSB sample main programs omg.c/omg.cpp/omg.F.

•

We have corrected the name of the library for neutralino-neutralino annihilation in our sample files MSSM/cs br.*.

Changes specific to the NMSSM model.

•

The NMSSM has been completely revised. Now it is based on NMSSMTools_1.0.2.

•

The deltaMb corrections in the NMSSM are included in the Higgs potential.

CP violation model.

•

We have included in our package the MSSM with CP violation. Our implementation was described in Phys. Rev. D 73 (2006) 115007. It is based on the CPSUPERH package published in Comput. Phys. Comm. 156 (2004) 283.

Unusual features:Depending on the parameters of the model, the program generates additional new code, compiles it and loads it dynamically.Running time:0.2 seconds     

GDF v2.0, an enhanced version of GDF

An improved version of the function estimation program GDF is presented. The main enhancements of the new version include: multi-output function estimation, capability of defining custom functions in the grammar and selection of the error function. The new version has been evaluated on a series of classification and regression datasets, that are widely used for the evaluation of such methods. It is compared to two known neural networks and outperforms them in 5 (out of 10) datasets.

Program summary

Title of program: GDF v2.0Catalogue identifier: ADXC_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXC_v2_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.: 98 147No. of bytes in distributed program, including test data, etc.: 2 040 684Distribution format: tar.gzProgramming language: GNU C++Computer: The program is designed to be portable in all systems running the GNU C++ compilerOperating system: Linux, Solaris, FreeBSDRAM: 200000 bytesClassification: 4.9Does the new version supersede the previous version?: YesNature of problem: The technique of function estimation tries to discover from a series of input data a functional form that best describes them. This can be performed with the use of parametric models, whose parameters can adapt according to the input data.Solution method: Functional forms are being created by genetic programming which are approximations for the symbolic regression problem.Reasons for new version: The GDF package was extended in order to be more flexible and user customizable than the old package. The user can extend the package by defining his own error functions and he can extend the grammar of the package by adding new functions to the function repertoire. Also, the new version can perform function estimation of multi-output functions and it can be used for classification problems.Summary of revisions: The following features have been added to the package GDF:

•

Multi-output function approximation. The package can now approximate any function . This feature gives also to the package the capability of performing classification and not only regression.

•

User defined function can be added to the repertoire of the grammar, extending the regression capabilities of the package. This feature is limited to 3 functions, but easily this number can be increased.

•

Capability of selecting the error function. The package offers now to the user apart from the mean square error other error functions such as: mean absolute square error, maximum square error. Also, user defined error functions can be added to the set of error functions.

•

More verbose output. The main program displays more information to the user as well as the default values for the parameters. Also, the package gives to the user the capability to define an output file, where the output of the gdf program for the testing set will be stored after the termination of the process.

Additional comments: A technical report describing the revisions, experiments and test runs is packaged with the source code.Running time: Depending on the train data.     

Maple procedures for the coupling of angular momenta. An up-date of the Racah module

An up-date of the Racah module is presented, adopted to Maple 11 and 12, which supports both, algebraic manipulations of expressions from Racah's algebra as well as numerical computations of many functions and symbols from the theory of angular momentum. The functions that are known to the program include the Wigner rotation matrices and n-j symbols, Clebsch-Gordan and Gaunt coefficients, spherical harmonics of various kinds as well as several others.

Program summary

Program title:RacahCatalogue identifier: ADFV_v10_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADFV_v10_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.: 30 436No. of bytes in distributed program, including test data, etc.: 544 866Distribution format: tar.gzProgramming language: Maple 11 and 12Computer: All computers with a license for the computer algebra package Maple [1]Operating system: Suse Linux 10.2+ and Ubuntu 8.10Classification: 4.1, 5Catalogue identifier of previous version: ADFV_v9_0Journal reference of previous version: Comput. Phys. Comm. 174 (2006) 616Does the new version supersede the previous version?: YesNature of problem: The theories of angular momentum and spherical tensor operators, sometimes known also as Racah's algebra, provide a powerful calculus for studying spin networks and (quantum) many-particle systems. For an efficient use of these theories, however, one requires not only a reliable handling of a large number of algebraic transformations and rules but, more often than not, also a fast access to their standard quantities, such as the Wigner n-j symbols, Clebsch-Gordan coefficients, spherical harmonics of various kinds, the rotation matrices, and many others.Solution method: A set of Maple procedures has been developed and maintained during the last decade which supports both, algebraic manipulations as well as fast computations of the standard expressions and symbols from the theory of angular momentum [2,3]. These procedures are based on a sizeable set of group-theoretical (and often rather sophisticated) relations which has been discussed and proven in the literature; see the monograph by Varshalovich et al. [4] for a comprehensive compilation. In particular the algebraic manipulation of complex (Racah) expressions may result in considerable simplifications, thus reducing the ‘numerical costs’, and often help obtain further insight into the behaviour of physical systems.Reasons for new version: A revision of the Racah module became necessary for mainly three reasons: (i) Since the last extension of the Racah procedures [5], which was developed within the framework of Maple 8, several updates of Maple were distributed by the vendors (currently Maple 13) and required a number of adaptations to the source code; (ii) the increasing size and program structure of the Racah module made it advisible to separate the (procedures for the treatment of the) atomic shell model from the manipulation and computation of Racah expressions. Therefore, the computation of angular coefficients for different coupling schemes, (grand) coefficients of fractional parentage as well as the matrix elements (of various irreducible tensors from the shell model) is to be maintained from now on independently within the Jucys module; (iii) a number of bugs and inconsistencies have been reported to us and corrected in the present version.Summary of revisions: In more detail, the following changes have been made:

1.

Since recent versions of Maple now support the automatic type checking of all incoming arguments and the definition of user-defined types; we have adapted most of the code to take advantage of these features, and especially those commands that are accessible by the user.

2.

In the computation of the Wigner n-j symbols and Clebsch-Gordan coefficients, we now return a ‘0’ in all cases in which the triangular rules are not fulfilled, for example, if δ(a,b,c)=0 for or . This change in the program saves the user making these tests on the quantum numbers explicitly everytime (in the summation over more complex expressions) that such a symbol or coefficient is invoked. The program still terminates with an error message if the (half-integer and integer) angular momentum quantum numbers appear in an inproper combination.

3.

While a recursive generation of the Wigner 3-j and 6-j symbols [6] may reduce the costs of some computations (and has thus been utilized in the past), it also makes the program rather sophisticated, especially if an algebraic evaluation or computations with a high number of Digits need to be supported by the same generic commands. The following procedures are therefore no longer supported by the Racah module:Racah_compute_w3j_jrange(), Racah_compute_w3j_mrange(),Racah_compute_w3j_recursive(), Racah_compute_w6j_range(), andRacah_compute_w6j_recursive().On most PCs, a sequential computation of all requested symbols is carried out within the same time basically.

4.

Because the module Jucys has grown to a size of about 35,000 lines of code and data, it appears helpful and necessary to maintain it independently. The procedures from the Jucys modules were designed to facilitate the computation of matrix elements of the unit tensors, the coefficients of fractional parentage (of various types) as well as transformation matrices between different coupling schemes [7] and are, thus, independent of the Racah module (although they typically require that the Racah code is available). The Jucys module is no longer distributed together with the present code.

5.

Apart from the Wigner n-j symbols (see above), some minor bugs have been reported and corrected in Racah_expand() and Racah_set().

6.

To facilitate the test of the installation and as a first tutorial on the module, we now provide the Maple worksheet Racah-tests-2009-maple12.mw in the Racah2009 root directory. This worksheet contains the examples and test cases from the previous versions. For the test of the installation, it is recommended that a ‘copy’ of this worksheet is saved and compared to the results from the re-run. It can be used also as a helpful source to define new examples in interactive work with the Racah module.

The Racah module is distributed in a tar file ADFV_v10_0.tar.gz from which the RACAH2009 root directory is (re-)generated by the command tar -zxvf ADFV_v10_0.tar.gz. This directory contains the source code libraries (tested for Maple 11 and 12), a Read.me for the installation of the program, the worksheet Racah-tests-2009-maple12.mw as well as the document Racah-commands-2009.pdf. This .pdf document serves as a Short Reference Manual and provides the definition of all the data structures of the Racah program together with an alphabetic list of all user relevant (and exported) commands. Although emphasis was placed on preserving the compatibility of the program with earlier releases of Maple, this cannot always be guaranteed due to changes in the Maple syntax. The Racah2009 root also contains an example of a .mapleinit file that can be modified and incorporated into the user's home directory to make the Racah module accessible like any other module of Maple. As mentioned above, the worksheet Racah-tests-2009-maple12.mw, help test the installation and may serve as a first tutorial.Restrictions: The (Racah) program is based on the concept of Racah expressions [cf. Fig. 1 in Ref. [4]] which, in principle, may contain any number of Wigner n-j symbols (n?9), Clebsch-Gordan coefficients, spherical harmonics and/or rotation matrices. In practise, of course, the required time and the success of an evaluation procedure depends on the complexity of the expressions and on the storage available, sometimes also on Maple's internal garbage treatment. In some cases, it is advisable to attempt first a simplification of the magnetic quantum numbers for a given expression before the summation over further 6-j and 9-j symbols should be taken into account. For all other quantities (that are compiled in Ref. [8], Tables 1 and 2, and explained in more detail in the Short Reference Manual, Racah-commands-2009.pdf), we currently just facilitate fast numerical computations by exploiting, as far as possible, Maple's hardware floating-point model. The program also supports simplifications on the Wigner rotation matrices. In integrals over the rotation matrices, products of up to three Wigner D-functions or reduced matrices (with the same angular arguments) are recognized; for the integration over a solid angle, however, the domain of integration must be specified explicitly for the Euler angles α and γ in order to force Maple to generate a constant of integration. In the course of the evaluation of Racah expressions, it is, in practice, often difficult to check internally whether all substructures of an expression are defined properly. Therefore, the user must ensure that all angular momenta (if given explicitly) must finally evaluate to integer and half-integer values and that they satisfy proper coupling conditions.Unusual features: The Racah program is designed for interactive use and for providing a quick and algebraic evaluation of (complex) expressions from Racah's algebra. In the evaluation, it exploits a large set of sum rules which are known from Racah's algebra and which may include (multiple) summations over dummy indices; see Varshalovich et al. [5] for a more detailed account of the theory. One strength of the program is that it recognizes automatically the symmetries of the symbols and functions, and that it applies also (some of) the graphical rules due to Yutsis and coworkers [9]. As before, the result of the evaluation process will be provided as Racah expressions, if a further simplification could be achieved, and may hence be used for further derivations and calculations within the given framework. In dealing with recoupling coefficients, these coefficients can be entered simply as a string of angular momenta (variables), separated by commas, and very similar to how they appear in mathematical texts. This is a crucial advantage of the program, compared with previous developments, for which the angular momenta and coupling coefficients had often to be given in a very detailed format. A Short Reference Manual to all procedures of the Racah program is provided by this distribution; it also contains the worksheet Racah-tests-2009-maple12.mw that contains the examples from all previous versions and may help test the installation. This worksheet can serve as a first tutorial to the Racah procedures. In the past, the Racah program has been utilized extensively in a number of applications including angular and polarization studies of heavy ions [10], angular distributions and correlation functions following photon-induced excitation processes [11], entanglement studies [12], in application of point-group symmetries and several others.Running time: The worksheet supplied with the distribution takes about 1 minute to run.References:

[1] Maple is a registered trademark of Waterloo Maple Inc.

[2] S. Fritzsche, Comp. Phys. Commun. 103 (1997) 51.

[3] S. Fritzsche, S. Varga, D. Geschke, B. Fricke, Comp. Phys. Commun. 111 (1998) 167;

T. Ingho, S. Fritzsche, B. Fricke, Comp. Phys. Commun. 139 (2001) 297;

S. Fritzsche, T. Ingho, T. Bastug, M. Tomaselli, Comp. Phys. Commun. 139 (2001) 314.

[4] D.A. Varshalovich, A.N. Moskalev, V.K. Khersonskii, Quantum Theory of Angular Momentum, World Scientific, Singapore a.o., 1988.

[5] J. Pagaran, S. Fritzsche, G. Gaigalas, Comp. Phys. Commun. 174 (2006) 616.

[6] K. Schulten, R.G. Gordon, Comp. Phys. Commun. 11 (1976) 269.

[7] G. Gaigalas, S. Fritzsche, B. Fricke, Comp. Phys. Commun. 135 (2001) 219;

G. Gaigalas, S. Fritzsche, Comp. Phys. Commun. 149 (2002) 39;

G. Gaigalas, O. Scharf, S. Fritzsche, Comp. Phys. Commun. 166 (2005) 141.

[8] S. Fritzsche, T. Ingho, M. Tomaselli, Comp. Phys. Commun. 153 (2003) 424.

[9] A.P. Yutsis, I.B. Levinson, V.V. Vanagas, The Theory of Angular Momentum, Israel Program for Scientific Translation, Jerusalem, 1962.

[10] S. Fritzsche, P. Indelicato, T. Stöhlker, J. Phys. B 38 (2005) S707.

[11] M. Kitajima, M. Okamoto, M. Hoshino, et al., J. Phys. B 35 (2002) 3327;

N.M. Kabachnik, S. Fritzsche, A.N. Grum-Grzhimailo, et al., Phys. Reports 451 (2007) 155;

S. Fritzsche, A.N. Grum-Grzhimailo, E.V. Gryzlova, N.M. Kabachnik, J. Phys. B 41 (2008) 165601;

T. Radtke, et al., Phys. Rev. A 77 (2008) 022507.

[12] T. Radtke, S. Fritzsche, Comp. Phys. Commun. 175 (2006) 145.

     

Model-Driven Development for scientific computing. An upgrade of the RHEEDGr program

Model-Driven Engineering (MDE) is the software engineering discipline, which considers models as the most important element for software development, and for the maintenance and evolution of software, through model transformation. Model-Driven Architecture (MDA) is the approach for software development under the Model-Driven Engineering framework. This paper surveys the core MDA technology that was used to upgrade of the RHEEDGR program to C++0x language standards.

New version program summary

Program title: RHEEDGR-09Catalogue identifier: ADUY_v3_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADUY_v3_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.: 21 263No. of bytes in distributed program, including test data, etc.: 1 266 982Distribution format: tar.gzProgramming language: Code Gear C++ BuilderComputer: Intel Core Duo-based PCOperating system: Windows XP, Vista, 7RAM: more than 1 MBClassification: 4.3, 7.2, 6.2, 8, 14Does the new version supersede the previous version?: YesNature of problem: Reflection High-Energy Electron Diffraction (RHEED) is a very useful technique for studying growth and surface analysis of thin epitaxial structures prepared by the Molecular Beam Epitaxy (MBE). The RHEED technique can reveal, almost instantaneously, changes either in the coverage of the sample surface by adsorbates or in the surface structure of a thin film.Solution method: The calculations are based on the use of a dynamical diffraction theory in which the electrons are taken to be diffracted by a potential, which is periodic in the dimension perpendicular to the surface.Reasons for new version: Responding to the user feedback the graphical version of the RHEED program has been upgraded to C++0x language standards. Also, functionality and documentation of the program have been improved.Summary of revisions:

1.

Model-Driven Architecture (MDA) is the approach defined by the Object Management Group (OMG) for software development under the Model-Driven Engineering framework [1]. The MDA approach shifts the focus of software development from writing code to building models. By adapting a model-centric approach, the MDA approach hopes to automate the generation of system implementation artifacts directly from the model. The following three models are the core of the MDA: (i) the Computation Independent Model (CIM), which is focused on basic requirements of the system, (ii) the Platform Independent Model (PIM), which is used by software architects and designers, and is focused on the operational capabilities of a system outside the context of a specific platform, and (iii) the Platform Specific Model (PSM), which is used by software developers and programmers, and includes details relating to the system for a specific platform. Basic requirements for the calculation of the RHEED intensity rocking curves in the one-beam condition have been described in Ref. [2]. Fig. 1 shows the PIM for the present version of the program. Fig. 2 presents the PSM for the program.

2.

The TGraph2D.bpk package has been recompiled to Graph2D0x.bpl and upgraded according to C++0x language standards. Fig. 3 shows the PSM of the Graph2D component, which is manifested by the Graph2D0x.bpl package presently. This diagram is a graphic presentation of the static view, which shows a collection of declarative model elements and their relationships. Installation instructions of the Graph2D0x package can be found in the new distribution.

3.

The program requires the user to provide the appropriate parameters for the crystal structure under investigation. These parameters are loaded from the parameters.ini file at run-time. Instructions for the preparation of the .ini files can be found in the new distribution.

4.

The program enables carrying out one-dimensional dynamical calculations for the fcc lattice, with a two-atoms basis and fcc lattice, with one atom basis but yet the zeroth Fourier component of the scattering potential in the TRHEED1D::crystPotUg() function can be modified according to users' specific application requirements.

5.

A graphical user interface (GUI) for the program has been reconstructed.

6.

The program has been compiled with English/USA regional and language options.

Unusual features: The program is distributed in the form of main projects RHEEDGr_09.cbproj and Graph2D0x.cbproj with associated files, and should be compiled using Code Gear C++ Builder 2009 compilers.Running time: The typical running time is machine and user-parameters dependent.References:

[1]

OMG, Model Driven Architecture Guide Version 1.0.1, 2003, http://www.omg.org/cgi-bin/doc?omg/03-06-01.

[2]

A. Daniluk, Comput. Phys. Comm. 166 (2005) 123.

     

FERM3D: A finite element R-matrix electron molecule scattering code

FERM3D is a three-dimensional finite element program, for the elastic scattering of a low energy electron from a general polyatomic molecule, which is converted to a potential scattering problem. The code is based on tricubic polynomials in spherical coordinates. The electron-molecule interaction is treated as a sum of three terms: electrostatic, exchange, and polarization. The electrostatic term can be extracted directly from ab initio codes (GAUSSIAN 98 in the work described here), while the exchange term is approximated using a local density functional. A local polarization potential based on density functional theory [C. Lee, W. Yang, R.G. Parr, Phys. Rev. B 37 (1988) 785] describes the long range attraction to the molecular target induced by the scattering electron. Photoionization calculations are also possible and illustrated in the present work. The generality and simplicity of the approach is important in extending electron-scattering calculations to more complex targets than it is possible with other methods.

Program summary

Title of program:FERM3DCatalogue identifier:ADYL_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYL_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandComputer for which the program is designed and others on which it has been tested:Intel Xeon, AMD Opteron 64 bit, Compaq AlphaOperating systems or monitors under which the program has been tested:HP Tru64 Unix v5.1, Red Hat Linux Enterprise 3Programming language used:Fortran 90Memory required to execute with typical data:900 MB (neutral CO₂), 2.3 GB (ionic CO₂), 1.4 GB (benzene)No. of bits in a word:32No. of processors used:1Has the code been vectorized?:NoNo. of lines in distributed program, including test data, etc.:58 383No. of bytes in distributed program, including test data, etc.:561 653Distribution format:tar.gzip fileCPC Program library subprograms used:ADDA, ACDPNature of physical problem:Scattering of an electron from a complex polyatomic molecular target.Method of solution:Solution of a partial differential equation using a finite element basis, and direct sparse linear solvers.Restrictions on the complexity of the problem:Memory constraints.Typical running time:2 hours.Unusual features of the program:

•

very extensive use of memory,

•

requires installation of Lapack, Blas, a direct sparse solver library (SuperLU, freely available, or Pardiso, which requires a license, but is free of charge for academic use), and optionally the Cernlib and Arpack libraries, freely available,

•

requires input from quantum chemistry programs (Gaussian, Molpro or PC Gamess).

     

yambo: An ab initio tool for excited state calculations

yambo is an ab initio code for calculating quasiparticle energies and optical properties of electronic systems within the framework of many-body perturbation theory and time-dependent density functional theory. Quasiparticle energies are calculated within the GW approximation for the self-energy. Optical properties are evaluated either by solving the Bethe-Salpeter equation or by using the adiabatic local density approximation. yambo is a plane-wave code that, although particularly suited for calculations of periodic bulk systems, has been applied to a large variety of physical systems. yambo relies on efficient numerical techniques devised to treat systems with reduced dimensionality, or with a large number of degrees of freedom. The code has a user-friendly command-line based interface, flexible I/O procedures and is interfaced to several publicly available density functional ground-state codes.

Program summary

Program title:yamboCatalogue identifier: AEDH_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDH_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU General Public Licence v2.0No. of lines in distributed program, including test data, etc.: 149 265No. of bytes in distributed program, including test data, etc.: 2 848 169Distribution format: tar.gzProgramming language: Fortran 95, CComputer: any computer architecture, running any flavor of UNIXOperating system: GNU/Linux, AIX, Irix, OS/XHas the code been vectorised or parallelized?: YesRAM: 10-1000 MbytesClassification: 7.3, 4.4, 7.2External routines:

•

BLAS (http://www.netlib.org/blas/)

•

LAPACK (http://www.netlib.org/lapack/)

•

MPI (http://www-unix.mcs.anl.gov/mpi/) is optional.

•

BLACS (http://www.netlib.org/scalapack/) is optional.

•

SCALAPACK (http://www.netlib.org/scalapack/) is optional.

•

FFTW (http://www.fftw.org/) is optional.

•

netCDF (http://www.unidata.ucar.edu/software/netcdf/) is optional.

Nature of problem: Calculation of excited state properties (quasiparticles, excitons, plasmons) from first principles.Solution method: Many body perturbation theory (Dyson equation, Bethe Salpeter equation) and time-dependent density functional theory. Quasiparticle approximation. Plasmon-pole model for the dielectric screening. Plane wave basis set with norm conserving pseudopotentials.Unusual features: During execution, yambo supplies estimates of the elapsed and remaining time for completion of each runlevel. Very friendly shell-based user-interface.Additional comments:yambo was known as “SELF” prior to GPL release. It belongs to the suite of codes maintained and used by the European Theoretical Spectroscopy Facility (ETSF) [1].Running time: The typical yambo running time can range from a few minutes to some days depending on the chosen level of approximation, and on the property and physical system under study.References:[1] The European Theoretical Spectroscopy Facility, http://www.etsf.eu.     

Ten commandments for good default expression simplification

This article provides goals for the design and improvement of default computer algebra expression simplification. These goals can also help users recognize and partially circumvent some limitations of their current computer algebra systems. Although motivated by computer algebra, many of the goals are also applicable to manual simplification, indicating what transformations are necessary and sufficient for good simplification when no particular canonical result form is required.After motivating the ten goals, the article then explains how the Altran partially factored form for rational expressions was extended for Derive and for the computer algebra in Texas Instruments products to help fulfill these goals. In contrast to the distributed Altran representation, this recursive partially factored semi-fraction form:

•

does not unnecessarily force common denominators,

•

discovers and preserves significantly more factors,

•

can represent general expressions, and

•

can produce an entire spectrum from fully factored over a common denominator through complete multivariate partial fractions, including a dense subset of all intermediate forms.

     

A Maple package for verifying ultradiscrete soliton solutions

We present a computer algebra program for verifying soliton solutions of ultradiscrete equations in which both dependent and independent variables take discrete values. The package is applicable to equations and solutions that include the max function. The program is implemented using Maple software.

Program summary

Program title: UltdeCatalogue identifier: AEDB_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDB_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.: 3171No. of bytes in distributed program, including test data, etc.: 13 633Distribution format: tar.gzProgramming language: Maple 10Computer: PC/AT compatible machineOperating system: Windows 2000, Windows XPRAM: Depends on the problem; minimum about 1 GBWord size: 32 bitsClassification: 5Nature of problem: The existence of multi-soliton solutions strongly suggest the integrability of nonlinear evolution equations. However enormous calculation is required to verify multi-soliton solutions of ultradiscrete equations. The use of computer algebra can be helpful in such calculations.Solution method: Simplification by using the properties of max-plus algebra.Restrictions: The program can only handle single ultradiscrete equations.Running time: Depends on the complexity of the equation and solution. For the examples included in the distribution the run times are as follows. (Core 2 Duo 3 GHz, Windows XP)

•

Example 1: 2725 sec.

•

Example 2: 33 sec.

•

Example 3: 1 sec.

     

FFT-split-operator code for solving the Dirac equation in 2+1 dimensions

The main part of the code presented in this work represents an implementation of the split-operator method [J.A. Fleck, J.R. Morris, M.D. Feit, Appl. Phys. 10 (1976) 129-160; R. Heather, Comput. Phys. Comm. 63 (1991) 446] for calculating the time-evolution of Dirac wave functions. It allows to study the dynamics of electronic Dirac wave packets under the influence of any number of laser pulses and its interaction with any number of charged ion potentials. The initial wave function can be either a free Gaussian wave packet or an arbitrary discretized spinor function that is loaded from a file provided by the user. The latter option includes Dirac bound state wave functions. The code itself contains the necessary tools for constructing such wave functions for a single-electron ion. With the help of self-adaptive numerical grids, we are able to study the electron dynamics for various problems in 2+1 dimensions at high spatial and temporal resolutions that are otherwise unachievable.Along with the position and momentum space probability density distributions, various physical observables, such as the expectation values of position and momentum, can be recorded in a time-dependent way. The electromagnetic spectrum that is emitted by the evolving particle can also be calculated with this code. Finally, for planning and comparison purposes, both the time-evolution and the emission spectrum can also be treated in an entirely classical relativistic way.Besides the implementation of the above-mentioned algorithms, the program also contains a large C++ class library to model the geometric algebra representation of spinors that we use for representing the Dirac wave function. This is why the code is called “Dirac++”.

Program summary

Program title: Dirac++ or (abbreviated) d++Catalogue identifier: AEAS_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAS_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.: 474 937No. of bytes in distributed program, including test data, etc.: 4 128 347Distribution format: tar.gzProgramming language: C++Computer: Any, but SMP systems are preferredOperating system: Linux and MacOS X are actively supported by the current version. Earlier versions were also tested successfully on IRIX and AIXNumber of processors used: Generally unlimited, but best scaling with 2-4 processors for typical problemsRAM: 160 Megabytes minimum for the examples given hereClassification: 2.7External routines: FFTW Library [3,4], Gnu Scientific Library [5], bzip2, bunzip2Nature of problem: The relativistic time evolution of wave functions according to the Dirac equation is a challenging numerical task. Especially for an electron in the presence of high intensity laser beams and/or highly charged ions, this type of problem is of considerable interest to atomic physicists.Solution method: The code employs the split-operator method [1,2], combined with fast Fourier transforms (FFT) for calculating any occurring spatial derivatives, to solve the given problem. An autocorrelation spectral method [6] is provided to generate a bound state for use as the initial wave function of further dynamical studies.Restrictions: The code in its current form is restricted to problems in two spatial dimensions. Otherwise it is only limited by CPU time and memory that one can afford to spend on a particular problem.Unusual features: The code features dynamically adapting position and momentum space grids to keep execution time and memory requirements as small as possible. It employs an object-oriented approach, and it relies on a Clifford algebra class library to represent the mathematical objects of the Dirac formalism which we employ. Besides that it includes a feature (typically called “checkpointing”) which allows the resumption of an interrupted calculation.Additional comments: Along with the program's source code, we provide several sample configuration files, a pre-calculated bound state wave function, and template files for the analysis of the results with both MatLab and Igor Pro.Running time: Running time ranges from a few minutes for simple tests up to several days, even weeks for real-world physical problems that require very large grids or very small time steps.References:

[1]

J.A. Fleck, J.R. Morris, M.D. Feit, Time-dependent propagation of high energy laser beams through the atmosphere, Appl. Phys. 10 (1976) 129-160.

[2]

R. Heather, An asymptotic wavefunction splitting procedure for propagating spatially extended wavefunctions: Application to intense field photodissociation of H⁺₂, Comput. Phys. Comm. 63 (1991) 446.

[3]

M. Frigo, S.G. Johnson, FFTW: An adaptive software architecture for the FFT, in: Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing, vol. 3, IEEE, 1998, pp. 1381-1384.

[4]

M. Frigo, S.G. Johnson, The design and implementation of FFTW3, in: Proceedings of the IEEE, vol. 93, IEEE, 2005, pp. 216-231. URL: http://www.fftw.org/.

[5]

M. Galassi, J. Davies, J. Theiler, B. Gough, G. Jungman, M. Booth, F. Rossi, GNU Scientific Library Reference Manual, second ed., Network Theory Limited, 2006. URL: http://www.gnu.org/software/gsl/.

[6]

M.D. Feit, J.A. Fleck, A. Steiger, Solution of the Schrödinger equation by a spectral method, J. Comput. Phys. 47 (1982) 412-433.