NMHDECAY 2.1: An updated program for sparticle masses, Higgs masses, couplings and decay widths in the NMSSM

We describe the improved properties of the NMHDECAY program, that is designed to compute Higgs and sparticle masses and Higgs decay widths in the NMSSM. In the version 2.0, Higgs decays into squarks and sleptons are included, accompanied by a calculation of the squark, gluino and slepton spectrum and tests against constraints from LEP and the Tevatron. Further radiative corrections are included in the Higgs mass calculation. A link to MicrOMEGAs allows to compute the dark matter relic density, and a rough (lowest order) calculation of BR(b→sγ) is performed. Finally, version 2.1 allows to integrate the RGEs for the soft terms up to the GUT scale.

Program summary

Title of program:NMHDECAY_SCAN, NMHDECAY_SLHACatalogue identifier:ADXW_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXW_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions:noneProgramming language used:FortranComputer:Mac, PC, Sun, Dec, AlphaOperating system:Mac OSX, Linux, Unix, WindowsNo. of lines in distributed program, including test data, etc.:20 060No. of bytes in distributed program, including test data, etc.:133 644RAM:2M bytesDistribution format:tar.gzNumber of processors used:1Classification:11.6Journal reference of previous version:JHEP 0502:066, 2005Does the new version supersede the previous version?:YesNature of problem:Computation of the Higgs and sparticle spectrum in the NMSSM and check of theoretical and experimental constraints.Solution method:Mass matrices including up to 2 loop radiative corrections for the Higgs bosons and all sparticles are computed and diagonalized. All Higgs decay widths are computed and branching ratios are compared to experimental bounds. Renormalisation group equations are integrated up to the GUT scale using a modified Runge-Kutta method, in order to check for the absence of a Landau pole. A modified version of MicrOmegas_1.3 can be called in order to compute the relic density of the lightest sparticle.Reasons for the new version:Higgs to sparticle decays added, computation of dark matter relic density added.Summary of revisions:Treatment of RGEs and radiative corrections improved, Higgs to sparticle decays added, new link to MicrOmegas_1.3.Restrictions:noneUnusual features:noneRunning time:<1 s per point in parameter space     

-XSummer- Transcendental functions and symbolic summation in Form

Harmonic sums and their generalizations are extremely useful in the evaluation of higher-order perturbative corrections in quantum field theory. Of particular interest have been the so-called nested sums, where the harmonic sums and their generalizations appear as building blocks, originating for example, from the expansion of generalized hypergeometric functions around integer values of the parameters. In this paper we discuss the implementation of several algorithms to solve these sums by algebraic means, using the computer algebra system Form.

Program summary

Title of program:XSummerCatalogue identifier:ADXQ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXQ_v1_0Program obtainable from:CPC Program Library, Queen's University of Belfast, N. IrelandLicense:GNU Public License and Form LicenseComputers:allOperating system:allProgram language:FormMemory required to execute:Depending on the complexity of the problem, recommended at least 64 MB RAMNo. of lines in distributed program, including test data, etc.:9854No. of bytes in distributed program, including test data, etc.:126 551Distribution format:tar.gzOther programs called:noneExternal files needed:noneNature of the physical problem:Systematic expansion of higher transcendental functions in a small parameter. The expansions arise in the calculation of loop integrals in perturbative quantum field theory.Method of solution:Algebraic manipulations of nested sums.Restrictions on complexity of the problem:Usually limited only by the available disk space.Typical running time:Dependent on the complexity of the problem.     

A program to compute the two-step excitation of mesospheric sodium atoms for the Polychromatic Laser Guide Star Project

A program to compute the two-step excitation of sodium atoms (3S1/2→3P3/2→4D5/2) using the density-matrix formalism is presented. The BEACON program calculates population evolution and the number of photons emitted by fluorescence from the 3P3/2, 4D5/2, 4P3/2, 4S1/2 levels.

Program summary

Title of program: BEACONCatalogue identifier:ADSXProgram Summary URL:http://cpc.cs.qub.ac.uk/cpc/summaries/ADSXProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneOperating systems under which the program has been tested: Win; UnixProgramming language used: FORTRAN 77Memory required to execute with typical data: 1 MwNumber of bits in a word: 32Number of processors used: 1 (a parallel version of this code is also available and can be obtained on request)Number of lines in distributed program, including test data, etc.: 29 287Number of bytes in distributed program, including test data, etc.: 830 331Distribution format: tar.gzCPC Program Library subprograms used: noneNature of physical problem: Resolution of the Bloch equations in the case of the two-step laser excitation of sodium atoms.Method of solution: The program BEACON calculates the evolution of level population versus time using the density-matrix formalism. The number of photons emitted from the 3P3/2, 4D5/2 and 4P3/2 levels is calculated using the branching ratios and the level lifetimes.Restriction on the complexity of the problem: Since the backscatter emission is calculated after the excitation process, excitation with laser pulse duration longer than the 4D5/2 level lifetime cannot be rigorously treated. Particularly, cw laser excitation cannot be calculated with this code.Typical running time:12 h     

BRANECODE: A program for simulations of braneworld dynamics

We describe an algorithm and a C++ implementation that we have written and made available for calculating the fully nonlinear evolution of 5D braneworld models with scalar fields. Bulk fields allow for the stabilization of the extra dimension. However, they complicate the dynamics of the system, so that analytic calculations (performed within an effective 4D theory) are usually only reliable for static bulk configurations or when the evolution of the extra dimension is negligible. In the general case, the nonlinear 5D dynamics can be studied numerically, and the algorithm and code we describe are the first ones of that type designed for this task. The program and its full documentation are available on the Web at http://www.cita.utoronto.ca/~jmartin/BRANECODE/.¹ In this paper we provide a brief overview of what the program does and how to use it.

Program summary

Title of program: BRANECODECatalogue identifier: ADVXProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVXProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneOperating systems under which the program has been tested: LinuxProgramming language used: C++Memory required to execute with typical data: less than 1 MBHas the code been vectorized?: noPeripherals used: noneNo. of lines in distributed program, including test data, etc.: 8277No. of bytes in distributed program, including test data, etc.: 74 939CPC Program Library subprograms used: noneNature of physical problem: Dynamics of two co-dimension one branes in a five-dimensional spacetime with a bulk scalar field and arbitrary potentials. The dynamics is governed by the five dimensional Einstein equations of gravity and the junction conditions at the position of the branes.Method of solution: Leapfrog algorithm to solve system of (1+1)-dimensional partial differential equations; Initial and boundary value problem.Restrictions on the complexity of the problem: Assumption of homogeneity along three spatial dimensions parallel to the branes.Typical running time: Depending on the grid size and length of the time evolution: from ∼1 s to ∼1 h or longer.Unusual features of the program:none     

PDSW: A program for the calculation of photon energy distribution resulting from radioactive elements in seawater

Photons, when emitted from radioactive sources in seawater, are subsequent to multiple scattering mechanisms, namely the photoelectric effect, the Compton scattering and the pair production effect. Thus, the monoenergetic emission of photons in seawater will result in equilibrium in a distribution of photons with different energies. PDSW is a MATLAB program which calculates this distribution and can be found useful for the characterization of measured spectra obtained by gamma detectors such as NaI(Tl). PDSW has been developed as an autonomous MATLAB function in order to make possible to integrate it in other applications. All calculations are performed using a typical value for seawater salinity (3.5%).

Program summary

Title of program: PDSWCatalogue identifier: ADWWProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADWWProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandComputer: x86Operating systems: WindowsProgramming language used: MATLABMemory required: 10 MbNumber of bits in a word: 32Number of processors used: 1Vectorized or parallelized?: noNumber of bytes in distributed program, including test data, etc.: 16 378Number of lines in distributed program, including test data, etc.: 3004Distribution format:tar.gzCPC Program Library subprograms used: noneNature of physical problem: Calculation of photon energy distribution in seawater taking into account the photoelectric effect, the Compton scattering and the pair production effect.Method of solution: Analytical calculation of the continuity equation for photon energy distribution in seawater and numerical integration of this equation in equilibrium.Restrictions on the complexity of the program: Very small resolution results in large memory requirements and high execution time.Typical running time: (Maximum energy-minimum resolution) 20 sUnusual features of the program: none     

BoltzTraP. A code for calculating band-structure dependent quantities

A program for calculating the semi-classic transport coefficients is described. It is based on a smoothed Fourier interpolation of the bands. From this analytical representation we calculate the derivatives necessary for the transport distributions. The method is compared to earlier calculations, which in principle should be exact within Boltzmann theory, and a very convincing agreement is found.

Program summary

Title of program:BoltzTraPCatalogue identifier:ADXU_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXU_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions:noneProgramming language used:Fortran 90Computer:The program should work on any system with a F90 compiler. The code has been tested with the Intel Fortran compilerOperating system:Unix/LinuxRAM:bytes up to 2 GB for low symmetry, small unit cell structuresNo. of lines in distributed program, including test data, etc.:1 534 213No. of bytes in distributed program, including test data, etc.:27 473 227Distribution format:tar.gzExternal routines:The LaPack and Blas libraries are neededNature of problem:Analytic expansion of energy-bands. Calculation of semi-classic integrals.Solution method:Smoothed Fourier expansion of bands.Running time:Up to 3 hours for low symmetry, small unit cell structures.     

USFKAD: An expert system for partial differential equations

The computation of the solution, by the separation of variables process, of the Poisson, diffusion, and wave equations in rectangular, cylindrical, or spherical coordinate systems, with Dirichlet, Neumann, or Robin boundary conditions, can be carried out in the time, Laplace, or frequency domains by a decision-tree process, using a library of eigenfunctions. We describe an expert system, USFKAD, that has been constructed for this purpose.

Program summary

Title of program:USFKADCatalogue identifier:ADYN_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYN_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions:noneOperating systems under which the program has been tested: Windows, UNIXProgramming language used:C++, LaTeXNo. of lines in distributed program, including test data, etc.: 11 699No. of bytes in distributed program, including test data, etc.: 537 744Memory required to execute with typical data: 1.3 MegabytesDistribution format: tar.gzNature of mathematical problem: Analytic solution of Poisson, diffusion, and wave equationsMethod of solution: Eigenfunction expansionsRestrictions concerning the complexity of the problem: A few rarely-occurring singular boundary conditions are unavailable, but they can be approximated by regular boundary value problems to arbitrary accuracy.Typical running time:1 secondUnusual features of the program: Solutions are obtained for Poisson, diffusion, or wave PDEs; homogeneous or nonhomogeneous equations and/or boundary conditions; rectangular, cylindrical, or spherical coordinates; time, Laplace, or frequency domains; Dirichlet, Neumann, Robin, singular, periodic, or incoming/outgoing boundary conditions. Output is suitable for pasting into LaTeX documents.     

A numerical algorithm for efficiently obtaining a Feynman parameter representation of one-gluon loop QCD Feynman diagrams for a large number of external gluons

A numerical program is presented which facilitates a computation pertaining to the full set of one-gluon loop diagrams (including ghost loop contributions), with M attached external gluon lines in all possible ways. The feasibility of such a task rests on a suitably defined master formula, which is expressed in terms of a set of Grassmann and a set of Feynman parameters. The program carries out the Grassmann integration and performs the Lorentz trace on the involved functions, expressing the result as a compact sum of parametric integrals. The computation is based on tracing the structure of the final result, thus avoiding all intermediate unnecessary calculations and directly writing the output. Similar terms entering the final result are grouped together. The running time of the program demonstrates its effectiveness, especially for large M.

Program summary

Program title:DILOG2Program identifier:ADXN_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXN_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandProgramming language:FORTRAN 90Computer(s) for which the program has been designed:Personal ComputerOperating system(s) for which the program has been designed: Windows 98, XP, LINUXNumber of processors used:oneNo. of lines in distributed program, including test data, etc.:2000No. of bytes in distributed program, including test data, etc.:16 249Distribution format:tar.gzExternal routines/libraries used:noneCPC Program Library subprograms used:noneNature of problem:The computation of one gluon/ghost loop diagrams in QCD with many external gluon lines is a time consuming task, practically beyond reasonable reach of analytic procedures. We apply recently proposed master formulas towards the computation of such diagrams with an arbitrary number (M) of external gluon lines, achieving a final result which reduces the problem to one involving integrals over the standard set, for given M, of Feynman parameters.Solution method:The structure of the master expressions is analyzed from a numerical computation point of view. Using the properties of Grassmann variables we identify all the different forms of terms that appear in the final result. Each form is called “structure”. We calculate theoretically the number of terms belonging to every “structure”. We carry out the calculation organizing the whole procedure into separate calculations of the terms belonging to every “structure”. Terms which do not contribute to the final result are thereby avoided. The final result, extending to large values of M, is also presented with terms belonging to the same “structure” grouped together.Restrictions:M is coded as a 2-digit integer. Overflow in the dimension of used array is expected to appear for M?20 in a processor that uses 4-bytes integers or for M?34 in a processor with 8-bytes integers.Running time:Depends on M, see enclosed figures.     

FracMAP: A user-interactive package for performing simulation and orientation-specific morphology analysis of fractal-like solid nano-agglomerates

Computer simulation techniques have found extensive use in establishing empirical relationships between three-dimensional (3d) and two-dimensional (2d) projected properties of particles produced by the process of growth through the agglomeration of smaller particles (monomers). In this paper, we describe a package, FracMAP, that has been written to simulate 3d quasi-fractal agglomerates and create their 2d pixelated projection images by restricting them to stable orientations as commonly encountered for quasi-fractal agglomerates collected on filter media for electron microscopy. Resulting 2d images are analyzed for their projected morphological properties.

Program summary

Program title: FracMAPCatalogue identifier: AEDD_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDD_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.: 4722No. of bytes in distributed program, including test data, etc.: 27 229Distribution format: tar.gzProgramming language: C++Computer: PCOperating system: Windows, LinuxRAM: 2.0 MegabytesClassification: 7.7Nature of problem: Solving for a suitable fractal agglomerate construction under constraints of typical morphological parameters.Solution method: Monte Carlo approximation.Restrictions: Problem complexity is not representative of run-time, since Monte Carlo iterations are of a constant complexity.Additional comments: The distribution file contains two versions of the FracMAP code, one for Windows and one for Linux.Running time: 1 hour for a fractal agglomerate of size 25 on a single processor.     

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.     

HypExp, a Mathematica package for expanding hypergeometric functions around integer-valued parameters

We present the Mathematica package HypExp which allows to expand hypergeometric functions around integer parameters to arbitrary order. At this, we apply two methods, the first one being based on an integral representation, the second one on the nested sums approach. The expansion works for both symbolic argument z and unit argument. We also implemented new classes of integrals that appear in the first method and that are, in part, yet unknown to Mathematica.

Program summary

Title of program:HypExpCatalogue identifier:ADXF_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXF_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicence:noneComputers:Computers running Mathematica under Linux or WindowsOperating system:Linux, WindowsProgram language:MathematicaNo. of bytes in distributed program, including test data, etc.:739 410No. of lines in distributed program, including test data, etc.:89 747Distribution format:tar.gzOther package needed:the package HPL, included in the distributionExternal file required:noneNature of the physical problem:Expansion of hypergeometric functions around integer-valued parameters. These are needed in the context of dimensional regularization for loop and phase space integrals.Method of solution:Algebraic manipulation of nested sums and integral representation.Restrictions on complexity of the problem:Limited by the memory availableTypical running time:Strongly depending on the problem and the availability of libraries.     

GeM software package for computation of symmetries and conservation laws of differential equations

We present a recently developed Maple-based “GeM” software package for automated symmetry and conservation law analysis of systems of partial and ordinary differential equations (DE). The package contains a collection of powerful easy-to-use routines for mathematicians and applied researchers. A standard program that employs “GeM” routines for symmetry, adjoint symmetry or conservation law analysis of any given DE system occupies several lines of Maple code, and produces output in the canonical form. Classification of symmetries and conservation laws with respect to constitutive functions and parameters present in the given DE system is implemented. The “GeM” package is being successfully used in ongoing research. Run examples include classical and new results.

Program summary

Title of program: GeMCatalogue identifier: ADYK_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYK_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputers: PC-compatible running Maple on MS Windows or Linux; SUN systems running Maple for Unix on OS SolarisOperating systems under which the program has been tested: Windows 2000, Windows XP, Linux, SolarisProgramming language used: Maple 9.5Memory required to execute with typical data: below 100 MegabytesNo. of lines in distributed program, including test data, etc.: 4939No. of bytes in distributed program, including test data, etc.: 166 906Distribution format: tar.gzNature of physical problem: Any physical model containing linear or nonlinear partial or ordinary differential equations.Method of solution: Symbolic computation of Lie, higher and approximate symmetries by Lie's algorithm. Symbolic computation of conservation laws and adjoint symmetries by using multipliers and Euler operator properties. High performance is achieved by using an efficient representation of the system under consideration and resulting symmetry/conservation law determining equations: all dependent variables and derivatives are represented as symbols rather than functions or expressions.Restrictions on the complexity of the problem: The GeM module routines are normally able to handle ODE/PDE systems of high orders (up to order seven and possibly higher), depending on the nature of the problem. Classification of symmetries/conservation laws with respect to one or more arbitrary constitutive functions of one or two arguments is normally accomplished successfully.Typical running time: 1-20 seconds for problems that do not involve classification; 5-1000 seconds for problems that involve classification, depending on complexity.     

Debris disk radiative transfer simulation tool (DDS)

A WWW interface for the simulation of spectral energy distributions of optically thin dust configurations with an embedded radiative source is presented. The density distribution, radiative source, and dust parameters can be selected either from an internal database or defined by the user. This tool is optimized for studying circumstellar debris disks where large grains (agrain ?1 μm) are expected to determine the far-infrared through millimeter dust reemission spectral energy distribution. The tool is available at http://aida28.mpia-hd.mpg.de/~swolf/dds.

Program summary

Title of program:Debris Disk Radiative Transfer Simulator (DDS)Catalogue identifier:ADVVProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVVProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions:noneComputers:PC with Intel(R) XEON(TM) 2.80 GHz processorOperating systems or monitors under which the program has been tested:SUSE Linux 9.1Programming language used:Fortran 90 (for the main program; furthermore Perl, CGI and HTML)Memory required to execute with typical data:10⁸ wordsNo. of bits in a word:8No. of lines in distributed program, including test data, etc.:44 636No. of bytes in distributed program, including test data, etc.: 4 806 280Distribution format:tar.gzNature of the physical problem:Simulation of scattered light and thermal reemission in arbitrary optically dust distributions with spherical, homogeneous grains where the dust parameters (optical properties, sublimation temperature, grain size) and SED of the illuminating/heating radiative source can be arbitrarily defined (example application: [S. Wolf, L.A. Hillenbrand, Astrophys. J. 596 (2003) 603]). The program is optimized for studying circumstellar debris disks where large grains (i.e. with large size parameters) are expected to determine the far-infrared through millimeter dust reemission spectral energy distribution.Method of solution:Calculation of the dust temperature distribution and dust reemission and scattering spectrum in the optically thin limit.Restrictions on the complexity of the problem:(1) The approach to calculate dust temperatures and dust reemission spectra is only valid in the optically thin regime. The validity of this constraint is verified for each model during the runtime of the code. (2) The relative abundances of different grains can be arbitrarily chosen, but must be constant outside the dust sublimation region, i.e. the shape of the (arbitrary) radial dust density distribution outside the dust sublimation region is the same for all grain sizes and chemistries. (3) The size of upload files (such as the dust density distribution, optical data of the dust grains, stellar spectral energy distribution, etc.) is limited (see http://aida28.mpia-hd.mpg.de/~swolf/dds/ for current file size limits). However, the resulting limitation to the complexity of possible model definitions is marginal only.Typical running time:3 sec-30 min (depending on the complexity of the model).Unusual features of the program:The program as provided through the CPC Program Library is equipped with an HTML user interface. It is installed and available at http://aida28.mpia-hd.mpg.de/~swolf/dds.     

NearFar: A computer program for nearside-farside decomposition of heavy-ion elastic scattering amplitude

The NearFar program is a package for carrying out an interactive nearside-farside decomposition of heavy-ion elastic scattering amplitude. The program is implemented in Java to perform numerical operations on the nearside and farside angular distributions. It contains a graphical display interface for the numerical results. A test run has been applied to the elastic scattering at Elab=1503 MeV.

Program summary

Title of program: NearFarCatalogue identifier: ADYP_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYP_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputers: designed for any machine capable of running Java, developed on PC-Pentium-4Operating systems under which the program has been tested: Microsoft Windows XP (Home Edition)Program language used: JavaNumber of bits in a word: 64Memory required to execute with typical data: case dependentNo. of lines in distributed program, including test data, etc.: 3484Number of bytes distributed program, including test data, etc.: 142 051Distribution format: tar.gzOther software required: A Java runtime interpreter, or the Java Development Kit, version 5.0Nature of physical problem: Interactive nearside-farside decomposition of heavy-ion elastic scattering amplitude.Method of solution: The user must supply a external data file or PPSM parameters which calculates theoretical values of the quantities to be decomposed.Typical running time: Problem dependent. In a test run, it is about 35 s on a 2.40 GHz Intel P4-processor machine.     

HPL, a Mathematica implementation of the harmonic polylogarithms

In this paper, we present an implementation of the harmonic polylogarithm of Remiddi and Vermaseren [E. Remiddi, J.A.M. Vermaseren, Int. J. Modern Phys. A 15 (2000) 725, hep-ph/9905237] for Mathematica. It contains an implementation of the product algebra, the derivative properties, series expansion and numerical evaluation. The analytic continuation has been treated carefully, allowing the user to keep the control over the definition of the sign of the imaginary parts. Many options enables the user to adapt the behavior of the package to his specific problem.

Program summary

Program title: HPLCatalogue identifier:ADWXProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADWXProgram obtained from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions:noneProgramming language: MathematicaNo. of lines in distributed program, including test data, etc.:13 310No. of bytes in distributed program, including test data, etc.: 1 990 584Distribution format: tar.gzComputer:all computers running MathematicaOperating systems:operating systems running MathematicaNature of problem: Computer algebraic treatment of the harmonic polylogarithms which appear in the evaluation of Feynman diagramsSolution method: Mathematica implementation     

PETAG01: A program for the direct simulation of a pellet target

We describe a numerical model of an internal pellet target to study the beam dynamics in storage rings, where the nuclear experiments with such type of target are planned. In this model the Monte Carlo algorithm is applied to evaluate the particle coordinates and momentum deviation depending on time and parameters of the target. One has to mention that due to statistical character of the pellet distribution in the target the analytical techniques are not applicable. This is also true for the particle distribution in the stored beam, which is influenced by various effects (such as a cooling process, intra-beam scattering, betatron oscillation, space charge effect). In this case only the Monte Carlo technique to model energy straggling in combination with the pellet distribution in the target should be considered.

Program summary

Program title: PETAG01Catalogue identifier: ADZV_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZV_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.: 1068No. of bytes in distributed program, including test data, etc.: 11 314Distribution format: tar.gzProgramming language: Fortran 77, C/C++Computer: Platform independentOperating system: MS Windows 95/2000/XP, Linux (Unix)RAM: 128 MBClassification: 11.10Nature of problem: Particle beam dynamics with use of the pellet target.Solution method: Monte Carlo with analytical approximation.Running time: dozens of seconds     

Computing zeros of analytic functions in the complex plane without using derivatives

We present a package in Fortran 90 which solves f(z)=0, where z∈W⊂C without requiring the evaluation of derivatives, f^′(z). W is bounded by a simple closed curve and f(z) must be holomorphic within W.We have developed and tested the package to support our work in the modeling of high frequency and optical wave guiding and resonant structures. The respective eigenvalue problems are particularly challenging because they require the high precision computation of all multiple complex roots of f(z) confined to the specified finite domain. Generally f(z), despite being holomorphic, does not have explicit analytical form thereby inhibiting evaluation of its derivatives.

Program summary

Title of program:EZEROCatalogue identifier:ADXY_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXY_v1_0Program obtainable from:CPC Program Library, Queen's University of Belfast, N. IrelandComputer:IBM compatible desktop PCOperating system:Fedora Core 2 Linux (with 2.6.5 kernel)Programming languages used:Fortran 90No. of bits in a word:32No. of processors used:oneHas the code been vectorized:noNo. of lines in distributed program, including test data, etc.:21045Number of bytes in distributed program including test data, etc.:223 756Distribution format:tar.gzPeripherals used:noneMethod of solution:Our package uses the principle of the argument to count the number of zeros encompassed by a contour and then computes estimates for the zeros. Refined results for each zero are obtained by application of the derivative-free Halley method with or without Aitken acceleration, as the user wishes.     

Vbfnlo: A parton level Monte Carlo for processes with electroweak bosons

Vbfnlo is a fully flexible parton level Monte Carlo program for the simulation of vector boson fusion, double and triple vector boson production in hadronic collisions at next-to-leading order in the strong coupling constant. Vbfnlo includes Higgs and vector boson decays with full spin correlations and all off-shell effects. In addition, Vbfnlo implements CP-even and CP-odd Higgs boson via gluon fusion, associated with two jets, at the leading-order one-loop level with the full top- and bottom-quark mass dependence in a generic two-Higgs-doublet model.A variety of effects arising from beyond the Standard Model physics are implemented for selected processes. This includes anomalous couplings of Higgs and vector bosons and a Warped Higgsless extra dimension model. The program offers the possibility to generate Les Houches Accord event files for all processes available at leading order.

Program summary

Program title:VbfnloCatalogue identifier: AEDO_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDO_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GPL version 2No. of lines in distributed program, including test data, etc.: 339 218No. of bytes in distributed program, including test data, etc.: 2 620 847Distribution format: tar.gzProgramming language: Fortran, parts in C++Computer: AllOperating system: Linux, should also work on other systemsClassification: 11.1, 11.2External routines: Optionally Les Houches Accord PDF Interface library and the GNU Scientific libraryNature of problem: To resolve the large scale dependence inherent in leading order calculations and to quantify the cross section error induced by uncertainties in the determination of parton distribution functions, it is necessary to include NLO corrections. Moreover, whenever stringent cuts are required on decay products and/or identified jets the question arises whether the scale dependence and a k-factor, defined as the ratio of NLO to LO cross section, determined for the inclusive production cross sections are valid for the search region one is interested in.Solution method: The problem is best addressed by implementing the one-loop QCD corrections in a fully flexible NLO parton-level Monte Carlo program, where arbitrary cuts can be specified as well as various scale choices. In addition, any currently available parton distribution function set can be used through the LHAPDF library.Running time: Depending on the process studied. Usually from minutes to hours.     

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.