ISICSoo: A class for the calculation of ionization cross sections from ECPSSR and PWBA theory

ISICS, originally a C language program for calculating K-, L- and M-shell ionization and X-ray production cross sections from ECPSSR and PWBA theory, has been reengineered into a C++ language class, named ISICSoo. The new software design enables the use of ISICS functionality in other software systems. The code, originally developed for Microsoft Windows operating systems, has been ported to Linux and Mac OS platforms to facilitate its use in a wider scientific environment. The reengineered software also includes some fixes to the original implementation, which ensure more robust computational results and a review of some physics parameters used in the computation. The paper describes the software design and the modifications to the implementation with respect to the previous version; it also documents the test process and provides some indications about the software performance.Program summaryProgram title: ISICSooCatalogue identifier: AEKN_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEKN_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.: 26 420No. of bytes in distributed program, including test data, etc.: 459 548Distribution format: tar.gzProgramming language: C++Computer: 80486 or higher-level PC or MacOperating system: Any OS with gcc compiler version 4.1 (or newer); tested on Scientific Linux 5 (gcc 4.1.2), Mac OS X 10.6.5 (gcc 4.2.1) and Windows XP (MS Visual C++ 2010 Express)Classification: 16.7Nature of problem: Ionization and X-ray production cross section calculations for ion–atom collisions.Solution method: Numerical integration of form factor using a logarithmic transform and Gaussian quadrature, plus exact integration limits.Additional comments: This program is a portable version of the program ADDS_v4_0.Reasons for the new version: Capability of using ISICS physics functionality in other software systems; porting the software to other platforms than Microsoft Windows; improved computational robustness and performance.Summary of revisions: Reengineering into a C++ class; several internal modifications to improve correctness and robustness; updated binding energies tabulations; performance improvements.Running time: The running time depends on the selected atomic shell and the number of polynomials used in the Gaussian quadrature integration. The examples provided only take seconds to run.     

An improved version of ISICS: a program for calculating K-, L- and M-shell cross sections from PWBA and ECPSSR theory using a personal computer

The C program, ISICS [Z. Liu, S.J. Cipolla, Comput. Phys. Comm. 97 (1996) 315-330], which calculates ionization and X-ray production cross-sections using PWBA and ECPSSR theory, has been enhanced to include new options, correct some minor flaws, and to make the program more versatile.

Program summary

Title of program: ISICSCatalog identifier: ADDS_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADDS_v2_0Program available from: CPC Program Library, Queen's University of Belfast, N. IrelandOperating system under which the program has been tested: WINDOWS XPProgram language used: CComputer: 80486 or higher-level PCsNo. of lines in distributed program, including test data, etc.: 5343No. of bytes in distributed program, including test data, etc.: 151 838Distribution format: tar.gzCatalogue identifier of previous version: ADDSJournal reference of previous version: Comput. Phys. Comm. 97 (1996) 315-330Does the new version supersede the previous version: YesNature of the physical problem: Ionization and X-ray production cross-section calculations for ion-atom collisions.Reasons for new version: Increased functionality and new options.Summary of revisions: Option for the united-atom approximation for binding-energy correction; easier inputting of updated atomic parameters; extension of projectile energy down to eV range; accounting for DHS wave function in K-shell ionization; other miscellaneous changes.Method of solution: Numerical integration of form factor using a logarithmic transform and Gaussian quadrature, plus exact integration limits.Restrictions on the complexity of the problem: The consumed CPU time increases with the atomic shell (K, L, M), but execution is still very fast.Typical running time: No change from previous version.Unusual features of the program: No     

TiReX: Replica-exchange molecular dynamics using Tinker

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

Program summary

Program title: TiReXCatalogue identifier: AEEK_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEEK_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 43 385No. of bytes in distributed program, including test data, etc.: 502 262Distribution format: tar.gzProgramming language: Fortran 90/95Computer: Most UNIX machinesOperating system: LinuxHas the code been vectorized or parallelized?: parallelized with MPIClassification: 16.13External routines: TINKER version 4.2 or 5.0, built as a libraryNature of problem: Replica-exchange molecular dynamics.Solution method: Each replica is assigned to a separate process; temperatures are swapped between replicas at regular time intervals.Running time: The sample run may take up to a few minutes.     

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

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

Program summary

Program title: svm-2dCatalogue identifier: AEBE_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBE_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 5091No. of bytes in distributed program, including test data, etc.: 130 963Distribution format: tar.gzProgramming language: Fortran 90Computer: The program should work on any system with a Fortran 90 compilerOperating system: The program should work on any system with a Fortran 90 compilerClassification: 7.3Nature of problem: Variational calculation of energies and wave functions using Correlated Gaussian basis.Solution method: Two-dimensional few-electron problems are solved by the variational method. The ground state wave function is expanded into Correlated Gaussian basis functions and the parameters of the basis states are optimized by a stochastic selection procedure. Accurate results can be obtained for 2-6 electron systems.Running time: A couple of hours for a typical system.     

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

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

Program summary

Program title: FHOLOCatalogue identifier: AEDS_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDS_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 176 336No. of bytes in distributed program, including test data, etc.: 4 294 872Distribution format: tar.gzProgramming language: C++Computer: any X86 micro computerOperating system: Linux (Debian Etch)RAM: 512 MBClassification: 18Nature of problem: To generate a Fourier Hologram in micro computer only by using open source operating system and packages.Running time: Depends on the matrix size. 10 sec for a matrix of size 256×256.     

SMMP v. 3.0—Simulating proteins and protein interactions in Python and Fortran

We describe a revised and updated version of the program package SMMP. SMMP is an open-source FORTRAN package for molecular simulation of proteins within the standard geometry model. It is designed as a simple and inexpensive tool for researchers and students to become familiar with protein simulation techniques. SMMP 3.0 sports a revised API increasing its flexibility, an implementation of the Lund force field, multi-molecule simulations, a parallel implementation of the energy function, Python bindings, and more.

Program summary

Title of program:SMMPCatalogue identifier:ADOJ_v3_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADOJ_v3_0.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.htmlProgramming language used:FORTRAN, PythonNo. of lines in distributed program, including test data, etc.:52 105No. of bytes in distributed program, including test data, etc.:599 150Distribution format:tar.gzComputer:Platform independentOperating system:OS independentRAM:2 MbytesClassification:3Does the new version supersede the previous version?:YesNature of problem:Molecular mechanics computations and Monte Carlo simulation of proteins.Solution method:Utilizes ECEPP2/3, FLEX, and Lund potentials. Includes Monte Carlo simulation algorithms for canonical, as well as for generalized ensembles.Reasons for new version:API changes and increased functionality.Summary of revisions:Added Lund potential; parameters used in subroutines are now passed as arguments; multi-molecule simulations; parallelized energy calculation for ECEPP; Python bindings.Restrictions:The consumed CPU time increases with the size of protein molecule.Running time:Depends on the size of the simulated molecule.     

Visual tool for estimating the fractal dimension of images

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

Program summary

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

A web-deployed interface for performing ab initio molecular dynamics, optimization, and electronic structure in Fireball

Fireball is an ab initio technique for fast local orbital simulations of nanotechnological, solid state, and biological systems. We have implemented a convenient interface for new users and software architects in the platform-independent Java language to access Fireball's unique and powerful capabilities. The graphical user interface can be run directly from a web server or from within a larger framework such as the Computational Science and Engineering Online (CSE-Online) environment or the Distributed Analysis of Neutron Scattering Experiments (DANSE) framework. We demonstrate its use for high-throughput electronic structure calculations and a multi-100 atom quantum molecular dynamics (MD) simulation.

Program summary

Program title: FireballUICatalogue identifier: AECF_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECF_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.: 279 784No. of bytes in distributed program, including test data, etc.: 12 836 145Distribution format: tar.gzProgramming language: JavaComputer: PC and workstationOperating system: The GUI will run under Windows, Mac and Linux. Executables for Mac and Linux are included in the package.RAM: 512 MBWord size: 32 or 64 bitsClassification: 4.14Nature of problem: The set up and running of many simulations (all of the same type), from the command line, is a slow process. But most research quality codes, including the ab initio tight-binding code FIREBALL, are designed to run from the command line. The desire is to have a method for quickly and efficiently setting up and running a host of simulations.Solution method: We have created a graphical user interface for use with the FIREBALL code. Once the user has created the files containing the atomic coordinates for each system that they are going to run a simulation on, the user can set up and start the computations of up to hundreds of simulations.Running time: 3 to 5 minutes on a 2 GHz Pentium IV processor.     

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

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

Program summary

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

JJGEN: A flexible program for generating lists of jj-coupled configuration state functions

The success of large scale relativistic multiconfiguration Dirac-Hartree-Fock calculations for atomic systems rely on judiciously chosen configuration expansions. Dependent on the atomic system as well as on the studied properties, various correlation effects need to be considered. Based on the active set approach, this program allows the user to generate general lists of jj-coupled configuration state functions to be used as input to the grasp2K multiconfiguration Dirac-Hartree-Fock package [P. Jönsson, X. He, C. Froese Fischer, I.P. Grant, Comput. Phys. Comm. (2007), in press].

Program summary

Program title: JJGENCatalogue identifier: ADZG_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZG_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 10 673No. of bytes in distributed program, including test data, etc.: 430 543Distribution format: tar.gzProgramming language: FortranComputer: Intel compatible PCOperating system: Linux, UnixWord size: 32 bitsClassification: 7.3Nature of problem: Generation of lists of jj-coupled configuration state functions to describe different electron correlation effects in many-electron atoms.Solution method: From a set of reference configurations a list of jj-coupled configuration state functions is generated by excitations to an active set of orbitals. Imposing restrictions on the allowed excitations the configuration expansion can be targeted to describe different correlation effects.Restrictions: The complexity of the cases that can be handled is entirely determined by the grasp2K package [P. Jönsson, X. He, C. Froese Fischer, I.P. Grant, Comput. Phys. Comm. (2007), in press] used for the generation of the electronic wave-functions.Running time: CPU time required to execute test cases: few seconds.     

HypExp 2, Expanding hypergeometric functions about half-integer parameters

In this article, we describe a new algorithm for the expansion of hypergeometric functions about half-integer parameters. The implementation of this algorithm for certain classes of hypergeometric functions in the already existing Mathematica package HypExp is described. Examples of applications in Feynman diagrams with up to four loops are given.

New version program summary

Program title:HypExp 2Catalogue identifier:ADXF_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXF_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.:106 401No. of bytes in distributed program, including test data, etc.:2 668 729Distribution format:tar.gzProgramming language:MathematicaComputer:Computers running MathematicaOperating system:Linux, Windows, MacRAM:Depending on the complexity of the problemSupplementary material:Library files which contain the expansion of certain hypergeometric functions around their parameters are availableClassification:4.7, 5Does the new version supersede the previous version?:YesNature of problem:Expansion of hypergeometric functions about parameters that are integer and/or half-integer valued.Solution method:New algorithm implemented in Mathematica.Reasons for new version:Expansion about half-integer parameters.Summary of revisions:Ability to expand about half-integer valued parameters added.Restrictions:The classes of hypergeometric functions with half-integer parameters that can be expanded are listed below.Additional comments:The package uses the package HPL included in the distribution.Running time:Depending on the expansion.     

Real-time Java simulations of multiple interference dielectric filters

An interactive Java applet for real-time simulation and visualization of the transmittance properties of multiple interference dielectric filters is presented. The most commonly used interference filters as well as the state-of-the-art ones are embedded in this platform-independent applet which can serve research and education purposes. The Transmittance applet can be freely downloaded from the site http://cpc.cs.qub.ac.uk.

Program summary

Program title: TransmittanceCatalogue identifier: AEBQ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBQ_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.: 5778No. of bytes in distributed program, including test data, etc.: 90 474Distribution format: tar.gzProgramming language: JavaComputer: Developed on PC-Pentium platformOperating system: Any Java-enabled OS. Applet was tested on Windows ME, XP, Sun Solaris, Mac OSRAM: VariableClassification: 18Nature of problem: Sophisticated wavelength selective multiple interference filters can include some tens or even hundreds of dielectric layers. The spectral response of such a stack is not obvious. On the other hand, there is a strong demand from application designers and students to get a quick insight into the properties of a given filter.Solution method: A Java applet was developed for the computation and the visualization of the transmittance of multilayer interference filters. It is simple to use and the embedded filter library can serve educational purposes. Also, its ability to handle complex structures will be appreciated as a useful research and development tool.Running time: Real-time simulations     

Algorithmic derivation of Dyson-Schwinger equations

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

Program summary

Program title: DoDSECatalogue identifier: AECT_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECT_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 105 874No. of bytes in distributed program, including test data, etc.: 262 446Distribution format: tar.gzProgramming language: Mathematica 6 and higherComputer: all on which Mathematica is availableOperating system: all on which Mathematica is availableClassification: 11.1, 11.4, 11.5, 11.6Nature of problem: Derivation of Dyson-Schwinger equations for a theory with given interactions.Solution method: Implementation of an algorithm for the derivation of Dyson-Schwinger equations.Unusual features: The results can be plotted as Feynman diagrams in Mathematica.Running time: Less than a second to minutes for Dyson-Schwinger equations of higher vertex functions.     

xPerm: fast index canonicalization for tensor computer algebra

We present a very fast implementation of the Butler-Portugal algorithm for index canonicalization with respect to permutation symmetries. It is called xPerm, and has been written as a combination of a Mathematica package and a C subroutine. The latter performs the most demanding parts of the computations and can be linked from any other program or computer algebra system. We demonstrate with tests and timings the effectively polynomial performance of the Butler-Portugal algorithm with respect to the number of indices, though we also show a case in which it is exponential. Our implementation handles generic tensorial expressions with several dozen indices in hundredths of a second, or one hundred indices in a few seconds, clearly outperforming all other current canonicalizers. The code has been already under intensive testing for several years and has been essential in recent investigations in large-scale tensor computer algebra.

Program summary

Program title: xPermCatalogue identifier: AEBH_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBH_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.: 93 582No. of bytes in distributed program, including test data, etc.: 1 537 832Distribution format: tar.gzProgramming language: C and Mathematica (version 5.0 or higher)Computer: Any computer running C and Mathematica (version 5.0 or higher)Operating system: Linux, Unix, Windows XP, MacOSRAM:: 20 MbyteWord size: 64 or 32 bitsClassification: 1.5, 5Nature of problem: Canonicalization of indexed expressions with respect to permutation symmetries.Solution method: The Butler-Portugal algorithm.Restrictions: Multiterm symmetries are not considered.Running time: A few seconds with generic expressions of up to 100 indices. The xPermDoc.nb notebook supplied with the distribution takes approximately one and a half hours to execute in full.     

QUBIT4MATLAB V3.0: A program package for quantum information science and quantum optics for MATLAB

A program package for MATLAB is introduced that helps calculations in quantum information science and quantum optics. It has commands for the following operations: (i) Reordering the qudits of a quantum register, computing the reduced state of a quantum register. (ii) Defining important quantum states easily. (iii) Formatted input and output for quantum states and operators. (iv) Constructing operators acting on given qudits of a quantum register and constructing spin chain Hamiltonians. (v) Partial transposition, matrix realignment and other operations related to the detection of quantum entanglement. (vi) Generating random state vectors, random density matrices and random unitaries.

Program summary

Program title:QUBIT4MATLAB V3.0Catalogue identifier:AEAZ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAZ_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.:5683No. of bytes in distributed program, including test data, etc.: 37 061Distribution format:tar.gzProgramming language:MATLAB 6.5; runs also on OctaveComputer:Any which supports MATLAB 6.5Operating system:Any which supports MATLAB 6.5; e.g., Microsoft Windows XP, LinuxClassification:4.15Nature of problem: Subroutines helping calculations in quantum information science and quantum optics.Solution method: A program package, that is, a set of commands is provided for MATLAB. One can use these commands interactively or they can also be used within a program.Running time:10 seconds-1 minute     

Monte Carlo event generator for black hole production and decay in proton-proton collisions - QBH version 1.02

We describe the Monte Carlo event generator for black hole production and decay in proton-proton collisions - QBH version 1.02. The generator implements a model for quantum black hole production and decay based on the conservation of local gauge symmetries and democratic decays. The code in written entirely in C++ and interfaces to the PYTHIA 8 Monte Carlo code for fragmentation and decays.

Program summary

Program title: QBHCatalogue identifier: AEGU_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGU_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 10 048No. of bytes in distributed program, including test data, etc.: 118 420Distribution format: tar.gzProgramming language: C++Computer: x86Operating system: Scientific Linux, Mac OS XRAM: 1 GBClassification: 11.6External routines: PYTHIA 8130 (http://home.thep.lu.se/~torbjorn/pythiaaux/present.html) and LHAPDF (http://projects.hepforge.org/lhapdf/)Nature of problem: Simulate black hole production and decay in proton-proton collision.Solution method: Monte Carlo simulation using importance sampling.Running time: Eight events per second.     

A Breit–Pauli distorted wave implementation for autostructure

We describe the Breit–Pauli distorted wave (BPDW) approach for the electron-impact excitation of atomic ions that we have implemented within the autostructure code.

Program summary

Program title:autostructureCatalogue identifier: AEIV_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEIV_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.: 130 987No. of bytes in distributed program, including test data, etc.: 1 031 584Distribution format: tar.gzProgramming language: Fortran 77/95Computer: GeneralOperating system: UnixHas the code been vectorized or parallelized?: Yes, a parallel version, with MPI directives, is included in the distribution.RAM: From several kbytes to several GbytesClassification: 2, 2.4Nature of problem: Collision strengths for the electron-impact excitation of atomic ions are calculated using a Breit–Pauli distorted wave approach with the optional inclusion of two-body non-fine-structure and fine-structure interactions.Solution method: General multi-configuration Breit–Pauli atomic structure. A jK-coupling partial wave expansion of the collision problem. Slater state angular algebra. Various model potential non-relativistic or kappa-averaged relativistic radial orbital solutions — the continuum distorted wave orbitals are not required to be orthogonal to the bound.Additional comments: Documentation is provided in the distribution file along with the test-case.Running time: From a few seconds to a few hours.     

LIBERI: Library for numerical evaluation of electron-repulsion integrals

We provide a C library, called LIBERI, for numerical evaluation of four-center electron repulsion integrals, based on successive reduction of integral dimension by using Fourier transforms. LIBERI enables us to compute the integrals for numerically defined basis functions within ¹⁰−5 Hartree accuracy as well as their derivatives with respect to the atomic nuclear positions. Damping of the Coulomb interaction can also be imposed to take account of screening effect.

Program summary

Program title: LIBERICatalogue identifier: AEGG_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGG_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.: 44 091No. of bytes in distributed program, including test data, etc.: 1 692 085Distribution format: tar.gzProgramming language: CComputer: allOperating system: any Unix-like systemRAM: 5-10 MbClassification: 7.4External routines: Lapack (http://www.netlib.org/lapack/), Blas (http://www.netlib.org/blas/), FFTW3 (http://www.fftw.org/)Nature of problem: Numerical evaluation of four-center electron-repulsion integrals.Solution method: Four-center electron-repulsion integrals are computed for given basis function set, based on successive reduction of integral dimension using Fourier transform.Running time: 0.5 sec for the demo program supplied with the package.     

The Invar tensor package

The Invar package is introduced, a fast manipulator of generic scalar polynomial expressions formed from the Riemann tensor of a four-dimensional metric-compatible connection. The package can maximally simplify any polynomial containing tensor products of up to seven Riemann tensors within seconds. It has been implemented both in Mathematica and Maple algebraic systems.

Program summary

Program title:Invar Tensor PackageCatalogue identifier:ADZK_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZK_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.: 136 240No. of bytes in distributed program, including test data, etc.:2 711 923Distribution format:tar.gzProgramming language:Mathematica and MapleComputer:Any computer running Mathematica versions 5.0 to 5.2 or Maple versions 9 and 10Operating system:Linux, Unix, Windows XPRAM:30 MbWord size:64 or 32 bitsClassification:5External routines:The Mathematica version requires the xTensor and xPerm packages. These are freely available at http://metric.iem.csic.es/Martin-Garcia/xActNature of problem:Manipulation and simplification of tensor expressions. Special attention on simplifying scalar polynomial expressions formed from the Riemann tensor on a four-dimensional metric-compatible manifold.Solution method:Algorithms of computational group theory to simplify expressions with tensors that obey permutation symmetries. Tables of syzygies of the scalar invariants of the Riemann tensor.Restrictions:The present versions do not fully address the problem of reducing differential invariants or monomials of the Riemann tensor with free indices.Running time:Less than a second to fully reduce a monomial of the Riemann tensor of degree 7 in terms of independent invariants.