Extension of HPL to complex arguments

In this paper we describe the extension of the Mathematica package HPL to treat harmonic polylogarithms of complex arguments. The harmonic polylogarithms have been introduced by Remiddi and Vermaseren [E. Remiddi, J.A.M. Vermaseren, Int. J. Modern Phys. A 15 (2000) 725, hep-ph/9905237] and have many applications in high energy particle physics.New version program summaryProgram title: HPLCatalogue identifier: ADWX_v2_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADWX_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.: 13 610No. of bytes in distributed program, including test data, etc.: 1 055 706Distribution format: tar.gzProgramming language: Mathematica 7/8.Computer: All computers running Mathematica.Operating system: Operating systems running Mathematica.Supplementary material: Additional “high weight” MinimalSet files available.Classification: 4.7.Catalogue identifier of previous version: ADWX_v1_0Journal reference of previous version: Comput. Phys. Comm. 174 (2006) 222Does the new version supersede the previous version?: YesNature of problem: Computer algebraic treatment of the harmonic polylogarithms which appear in the evaluation of Feynman diagrams.Solution method: Mathematica implementation.Reasons for new version: Added treatment of complex arguments. Details in arXiv:hep-ph/0703052.Summary of revisions: Added treatment of complex arguments. Details in arXiv:hep-ph/0703052.Running time: A few seconds for each function.     

Automatic calculation of supersymmetric renormalization group equations and loop corrections

SARAH is a Mathematica package for studying supersymmetric models. It calculates for a given model the masses, tadpole equations and all vertices at tree-level. This information can be used by SARAH to write model files for CalcHep/CompHep or FeynArts/FormCalc. In addition, the second version of SARAH can derive the renormalization group equations for the gauge couplings, parameters of the superpotential and soft-breaking parameters at one- and two-loop level. Furthermore, it calculates the one-loop self-energies and the one-loop corrections to the tadpoles. SARAH can handle all N=1 SUSY models whose gauge sector is a direct product of SU(N) and U(1) gauge groups. The particle content of the model can be an arbitrary number of chiral superfields transforming as any irreducible representation with respect to the gauge groups. To implement a new model, the user has just to define the gauge sector, the particle, the superpotential and the field rotations to mass eigenstates.

Program summary

Program title: SARAHCatalogue identifier: AEIB_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEIB_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.: 97 577No. of bytes in distributed program, including test data, etc.: 2 009 769Distribution format: tar.gzProgramming language: MathematicaComputer: All systems that Mathematica is available forOperating system: All systems that Mathematica is available forClassification: 11.1, 11.6Nature of problem: A supersymmetric model is usually characterized by the particle content, the gauge sector and the superpotential. It is a time consuming process to obtain all necessary information for phenomenological studies from these basic ingredients.Solution method: SARAH calculates the complete Lagrangian for a given model whose gauge sector can be any direct product of SU(N) gauge groups. The chiral superfields can transform as any, irreducible representation with respect to these gauge groups and it is possible to handle an arbitrary number of symmetry breakings or particle rotations. Also the gauge fixing terms can be specified. Using this information, SARAH derives the mass matrices and Feynman rules at tree-level and generates model files for CalcHep/CompHep and FeynArts/FormCalc. In addition, it can calculate the renormalization group equations at one- and two-loop level and the one-loop corrections to the one- and two-point functions.Unusual features: SARAH just needs the superpotential and gauge sector as input and not the complete Lagrangian. Therefore, the complete implementation of new models is done in some minutes.Running time: Measured CPU time for the evaluation of the MSSM on an Intel Q8200 with 2.33 GHz. Calculating the complete Lagrangian: 12 seconds. Calculating all vertices: 75 seconds. Calculating the one- and two-loop RGEs: 50 seconds. Calculating the one-loop corrections: 7 seconds. Writing a FeynArts file: 1 second. Writing a CalcHep/CompHep file: 6 seconds. Writing the LaTeX output: 1 second.     

SPICE: Simulation Package for Including Flavor in Collider Events

We describe SPICE: Simulation Package for Including Flavor in Collider Events. SPICE takes as input two ingredients: a standard flavor-conserving supersymmetric spectrum and a set of flavor-violating slepton mass parameters, both of which are specified at some high “mediation” scale. SPICE then combines these two ingredients to form a flavor-violating model, determines the resulting low-energy spectrum and branching ratios, and outputs HERWIG and SUSY Les Houches files, which may be used to generate collider events. The flavor-conserving model may be any of the standard supersymmetric models, including minimal supergravity, minimal gauge-mediated supersymmetry breaking, and anomaly-mediated supersymmetry breaking supplemented by a universal scalar mass. The flavor-violating contributions may be specified in a number of ways, from specifying charges of fields under horizontal symmetries to completely specifying all flavor-violating parameters. SPICE is fully documented and publicly available, and is intended to be a user-friendly aid in the study of flavor at the Large Hadron Collider and other future colliders.

Program summary

Program title: SPICECatalogue identifier: AEFL_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFL_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.: 8153No. of bytes in distributed program, including test data, etc.: 67 291Distribution format: tar.gzProgramming language: C++Computer: Personal computerOperating system: Tested on Scientific Linux 4.xClassification: 11.1External routines: SOFTSUSY [1,2] and SUSYHIT [3]Nature of problem: Simulation programs are required to compare theoretical models in particle physics with present and future data at particle colliders. SPICE determines the masses and decay branching ratios of supersymmetric particles in theories with lepton flavor violation. The inputs are the parameters of any of several standard flavor-conserving supersymmetric models, supplemented by flavor-violating parameters determined, for example, by horizontal flavor symmetries. The output are files that may be used for detailed simulation of supersymmetric events at particle colliders.Solution method: Simpson's rule integrator, basic algebraic computation.Additional comments: SPICE interfaces with SOFTSUSY and SUSYHIT to produce the low energy sparticle spectrum. Flavor mixing for sleptons and sneutrinos is fully implemented; flavor mixing for squarks is not included.Running time: <1 minute. Running time is dominated by calculating the possible and relevant three-body flavor-violating decays of sleptons, which is usually 10-15 seconds per slepton.References:

[1]

B.C. Allanach, Comput. Phys. Commun. 143 (2002) 305, arXiv:hep-ph/0104145.

[2]

B.C. Allanach, M.A. Bernhardt, arXiv:0903.1805 [hep-ph].

[3]

A. Djouadi, M.M. Muhlleitner, M. Spira, Acta Phys. Pol. B 38 (2007) 635, arXiv:hep-ph/0609292.

     

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.     

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.     

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

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

Program summary

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

A recursive method to calculate UV-divergent parts at one-loop level in dimensional regularization

A method is introduced to calculate the UV-divergent parts at one-loop level in dimensional regularization. The method is based on the recursion, and the basic integrals are just the scaleless integrals after the recursive reduction, which involve no other momentum scales except the loop momentum itself. The method can be easily implemented in any symbolic computer language, and a implementation in Mathematica is ready to use.Program summaryProgram title: UVPartCatalogue identifier: AELY_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AELY_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 361No. of bytes in distributed program, including test data, etc.: 412 084Distribution format: tar.gzProgramming language: MathematicaComputer: Any computer where the Mathematica is running.Operating system: Any capable of running Mathematica.Classification: 11.1External routines: FeynCalc (http://www.feyncalc.org/), FeynArts (http://www.feynarts.de/)Nature of problem: To get the UV-divergent part of any one-loop expression.Solution method: UVPart is a Mathematica package where the recursive method has been implemented.Running time: In general it is below one second.     

S@M, a mathematica implementation of the spinor-helicity formalism

In this paper we present the package S@M (Spinors@Mathematica) which implements the spinor-helicity formalism in Mathematica. The package allows the use of complex-spinor algebra along with the multi-purpose features of Mathematica. The package defines the spinor objects with their basic properties along with functions to manipulate them. It also offers the possibility of evaluating the spinorial objects numerically at every computational step. The package is therefore well suited to be used in the context of on-shell technology, in particular for the evaluation of scattering amplitudes at tree- and loop-level.

Program summary

Program title: S@MCatalogue identifier: AEBF_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBF_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.: 14 404No. of bytes in distributed program, including test data, etc.: 77 536Distribution format: tar.gzProgramming language: MathematicaComputer: All computers running MathematicaOperating system: Any system running MathematicaClassification: 4.4, 5, 11.1Nature of problem: Implementation of the spinor-helicity formalismSolution method: Mathematica implementationRunning time: The notebooks provided with the package take only a few seconds to run.     

Including R-parity violation in the numerical computation of the spectrum of the minimal supersymmetric standard model: SOFTSUSY3.0

Current publicly available computer programs calculate the spectrum and couplings of the minimal supersymmetric standard model under the assumption of R-parity conservation. Here, we describe an extension to the SOFTSUSY program which includes R-parity violating effects. The user provides a theoretical boundary condition upon the high-scale supersymmetry breaking R-parity violating couplings. Successful radiative electroweak symmetry breaking, electroweak and CKM matrix data are used as weak-scale boundary conditions. The renormalisation group equations are solved numerically between the weak scale and a high energy scale using a nested iterative algorithm. This paper serves as a manual to the R-parity violating mode of the program, detailing the approximations and conventions used.

Program summary

Program title:SOFTSUSY v3.0Catalogue identifier: ADPM_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADPM_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.: 75 927No. of bytes in distributed program, including test data, etc.: 570 916Distribution format: tar.gzProgramming language: C++, FortranComputer: Personal computerOperating system: Tested on Linux 4.xWord size: 32 bitsClassification: 11.6Catalogue identifier of previous version: ADPM_v1_0Journal reference of previous version: Comput. Phys. Comm. 143 (2002) 305Does the new version supersede the previous version?: YesNature of problem: Calculating supersymmetric particle spectrum and mixing parameters in the R-parity violating minimal supersymmetric standard model. The solution to the renormalisation group equations must be consistent with a high-scale boundary condition on supersymmetry breaking parameters and Rp parameters, as well as a weak-scale boundary condition on gauge couplings, Yukawa couplings and the Higgs potential parameters.Solution method: Nested iterative algorithmReasons for new version: This is an extension to the SOFTSUSY program which includes R-parity violating effects. The user provides a theoretical boundary condition upon the high-scale supersymmetry breaking R-parity violating couplings. Successful radiative electroweak symmetry breaking, electroweak and CKM matrix data are used as weak-scale boundary conditions. The renormalisation group equations are solved numerically between the weak scale and a high energy scale using a nested iterative algorithm. The paper serves as a manual to the R-parity violating mode of the program, detailing the approximations and conventions used.Restrictions:SOFTSUSY3.0 will provide a solution only in the perturbative regime and it assumes that all couplings of the MSSM are real (i.e. CP-conserving). The iterative SOFTSUSY algorithm will not converge if parameters are too close to a boundary of successful electroweak symmetry breaking, but a warning flag will alert the user to this fact.Running time: A few seconds per parameter point.     

QCDMAPT: Program package for Analytic approach to QCD

A program package, which facilitates computations in the framework of Analytic approach to QCD, is developed and described in detail. The package includes both the calculated explicit expressions for relevant spectral functions up to the four-loop level and the subroutines for necessary integrals.

Program summary

Program title: QCDMAPTCatalogue identifier: AEGP_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGP_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 2579No. of bytes in distributed program, including test data, etc.: 180 052Distribution format: tar.gzProgramming language: Maple 9 and higherComputer: Any which supports Maple 9Operating system: Any which supports Maple 9Classification: 11.1, 11.5, 11.6Nature of problem: Subroutines helping computations within Analytic approach to QCD.Solution method: A program package for Maple is provided. It includes both the explicit expressions for relevant spectral functions and the subroutines for basic integrals used in the framework of Analytic approach to QCD.Running time: Template program running time is about a minute (depends on CPU).     

LanHEP—a package for the automatic generation of Feynman rules in field theory. Version 3.0

The LanHEP program version 3.0 for Feynman rules generation from the Lagrangian is described. It reads the Lagrangian written in a compact form, close to the one used in publications. It means that Lagrangian terms can be written with summation over indices of broken symmetries and using special symbols for complicated expressions, such as covariant derivative and strength tensor for gauge fields. Supersymmetric theories can be described using the superpotential formalism and the 2-component fermion notation. The output is Feynman rules in terms of physical fields and independent parameters in the form of CompHEP model files, which allows one to start calculations of processes in the new physical model. Alternatively, Feynman rules can be generated in FeynArts format or as LaTeX table. One-loop counterterms can be generated in FeynArts format.

Program summary

Program title: LanHEPCatalogue identifier: ADZV_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECH_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.: 83 041No. of bytes in distributed program, including test data, etc.: 1 090 931Distribution format: tar.gzProgramming language: CComputer: PCOperating system: LinuxRAM: 2 MB (SM), 12 MB (MSSM), 120 MB (MSSM with counterterms)Classification: 4.4Nature of problem: Deriving Feynman rules from the LagrangianSolution method: The program reads the Lagrangian written in a compact form, close to the one used in publications. It means that Lagrangian terms can be written with summation over indices of broken symmetries and using special symbols for complicated expressions, such as covariant derivative and strength tensor for gauge fields. Tools for checking the correctness of the model, and for simplifying the output expressions are provided. The output is Feynman rules in terms of physical fields and independent parameters in the form of CompHEP model files, which allows one to start calculations of processes in the new physical model. Alternatively, Feynman rules can be generated in FeynArts format or as a LaTeX table.Running time: 1 sec (SM), 8 sec (MSSM), 8 min (MSSM with counterterms)     

Strongdeco: Expansion of analytical,strongly correlated quantum states into a many-body basis

We provide a Mathematica code for decomposing strongly correlated quantum states described by a first-quantized, analytical wave function into many-body Fock states. Within them, the single-particle occupations refer to the subset of Fock–Darwin functions with no nodes. Such states, commonly appearing in two-dimensional systems subjected to gauge fields, were first discussed in the context of quantum Hall physics and are nowadays very relevant in the field of ultracold quantum gases. As important examples, we explicitly apply our decomposition scheme to the prominent Laughlin and Pfaffian states. This allows for easily calculating the overlap between arbitrary states with these highly correlated test states, and thus provides a useful tool to classify correlated quantum systems. Furthermore, we can directly read off the angular momentum distribution of a state from its decomposition. Finally we make use of our code to calculate the normalization factors for Laughlin?s famous quasi-particle/quasi-hole excitations, from which we gain insight into the intriguing fractional behavior of these excitations.Program summaryProgram title: StrongdecoCatalogue identifier: AELA_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AELA_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.: 5475No. of bytes in distributed program, including test data, etc.: 31 071Distribution format: tar.gzProgramming language: MathematicaComputer: Any computer on which Mathematica can be installedOperating system: Linux, Windows, MacClassification: 2.9Nature of problem: Analysis of strongly correlated quantum states.Solution method: The program makes use of the tools developed in Mathematica to deal with multivariate polynomials to decompose analytical strongly correlated states of bosons and fermions into a standard many-body basis. Operations with polynomials, determinants and permanents are the basic tools.Running time: The distributed notebook takes a couple of minutes to run.     

AutoDipole - Automated generation of dipole subtraction terms -

We present an automated generation of the subtraction terms for next-to-leading order QCD calculations in the Catani-Seymour dipole formalism. For a given scattering process with n external particles our Mathematica package generates all dipole terms, allowing for both massless and massive dipoles. The numerical evaluation of the subtraction terms proceeds with MadGraph, which provides Fortran code for the necessary scattering amplitudes. Checks of the numerical stability are discussed.

Program summary

Program title: AutoDipoleCatalogue identifier: AEGO_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGO_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.: 138 042No. of bytes in distributed program, including test data, etc.: 1 117 665Distribution format: tar.gzProgramming language: Mathematica and FortranComputer: Computers running Mathematica (version 7.0)Operating system: The package should work on every Linux system supported by Mathematica. Detailed tests have been performed on Scientific Linux as supported by DESY and CERN and on openSUSE and Debian.RAM: Depending on the complexity of the problem, recommended at least 128 MB RAMClassification: 11.5External routines: MadGraph (including HELAS library) available under http://madgraph.hep.uiuc.edu/ or http://madgraph.phys.ucl.ac.be/ or http://madgraph.roma2.infn.it/. A copy of the tar file, MG_ME_SA_V4.4.30, is included in the AutoDipole distribution package.Nature of problem: Computation of next-to-leading order QCD corrections to scattering cross sections, regularization of real emission contributions.Solution method: Catani-Seymour subtraction method for massless and massive partons [1,2]; Numerical evaluation of subtracted matrix elements interfaced to MadGraph [3-5] (stand-alone version) using helicity amplitudes and the HELAS library [6,7] (contained in MadGraph).Restrictions: Limitations of MadGraph are inherited.Running time: Dependent on the complexity of the problem with typical run times of the order of minutes.References:

[1]

S. Catani, M.H. Seymour, Nuclear Phys. B 485 (1997) 291, hep-ph/9605323.

[2]

S. Catani, et al., Nuclear Phys. B 627 (2002) 189, hep-ph/0201036.

[3]

T. Stelzer, W.F. Long, Comput. Phys. Comm. 81 (1994) 357, hep-ph/9401258.

[4]

F. Maltoni, T. Stelzer, JHEP 0302 (2003) 027, hep-ph/0208156.

[5]

J. Alwall, et al., JHEP 0709 (2007) 028, arXiv:0706.2334 [hep-ph].

[6]

K. Hagiwara, H. Murayama, I. Watanabe, Nuclear Phys. B 367 (1991) 257.

[7]

H. Murayama, I. Watanabe, K. Hagiwara, KEK-91-11.

     

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.     

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.     

Symbolic computation of the Hartree-Fock energy from a chiral EFT three-nucleon interaction at NLO

We present the first of a two-part Mathematica notebook collection that implements a symbolic approach for the application of the density matrix expansion (DME) to the Hartree-Fock (HF) energy from a chiral effective field theory (EFT) three-nucleon interaction at N²LO. The final output from the notebooks is a Skyrme-like energy density functional that provides a quasi-local approximation to the non-local HF energy. In this paper, we discuss the derivation of the HF energy and its simplification in terms of the scalar/vector-isoscalar/isovector parts of the one-body density matrix. Furthermore, a set of steps is described and illustrated on how to extend the approach to other three-nucleon interactions.

Program summary

Program title: SymbHFNNNCatalogue identifier: AEGC_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGC_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.: 96 666No. of bytes in distributed program, including test data, etc.: 378 083Distribution format: tar.gzProgramming language: Mathematica 7.1Computer: Any computer running Mathematica 6.0 and later versionsOperating system: Windows Xp, Linux/UnixRAM: 256 MbClassification: 5, 17.16, 17.22Nature of problem: The calculation of the HF energy from the chiral EFT three-nucleon interaction at N²LO involves tremendous spin-isospin algebra. The problem is compounded by the need to eventually obtain a quasi-local approximation to the HF energy, which requires the HF energy to be expressed in terms of scalar/vector-isoscalar/isovector parts of the one-body density matrix. The Mathematica notebooks discussed in this paper solve the latter issue.Solution method: The HF energy from the chiral EFT three-nucleon interaction at N²LO is cast into a form suitable for an automatic simplification of the spin-isospin traces. Several Mathematica functions and symbolic manipulation techniques are used to obtain the result in terms of the scalar/vector-isoscalar/isovector parts of the one-body density matrix.Running time: Several hours     

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

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

Program summary

Program title: MilneCatalogue identifier: AEGS_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGS_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 701No. of bytes in distributed program, including test data, etc.: 6845Distribution format: tar.gzProgramming language: Fortran 77Computer: PC under Linux or WindowsOperating system: Ubuntu 8.04 (Kernel version 2.6.24-16-generic), Windows-XPClassification: 4.11, 21.1, 21.2Nature of problem: The X functions are integral expressions. The convergence of these regular and Cauchy Principal Value integrals are impaired by the singularities of the integrand in the complex plane. The DE quadrature scheme tackles these singularities in a robust manner compared to the standard Gauss quadrature.Running time: The test included in the distribution takes a few seconds to run.     

Generating and using truly random quantum states in Mathematica

The problem of generating random quantum states is of a great interest from the quantum information theory point of view. In this paper we present a package for Mathematica computing system harnessing a specific piece of hardware, namely Quantis quantum random number generator (QRNG), for investigating statistical properties of quantum states. The described package implements a number of functions for generating random states, which use Quantis QRNG as a source of randomness. It also provides procedures which can be used in simulations not related directly to quantum information processing.Program summaryProgram title: TRQSCatalogue identifier: AEKA_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEKA_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.: 7924No. of bytes in distributed program, including test data, etc.: 88 651Distribution format: tar.gzProgramming language: Mathematica, CComputer: Requires a Quantis quantum random number generator (QRNG, http://www.idquantique.com/true-random-number-generator/products-overview.html) and supporting a recent version of MathematicaOperating system: Any platform supporting Mathematica; tested with GNU/Linux (32 and 64 bit)RAM: Case dependentClassification: 4.15Nature of problem: Generation of random density matrices.Solution method: Use of a physical quantum random number generator.Running time: Generating 100 random numbers takes about 1 second, generating 1000 random density matrices takes more than a minute.