Automatic calculation of particle interaction or decay

Last updated

The automatic calculation of particle interaction or decay is part of the computational particle physics branch. It refers to computing tools that help calculating the complex particle interactions as studied in high-energy physics, astroparticle physics and cosmology. The goal of the automation is to handle the full sequence of calculations in an automatic (programmed) way: from the Lagrangian expression describing the physics model up to the cross-sections values and to the event generator software.

Contents

Overview

Particle accelerators or colliders produce collisions (interactions) of particles (like the electron or the proton). The colliding particles form the Initial State. In the collision, particles can be annihilated or/and exchanged producing possibly different sets of particles, the Final States. The Initial and Final States of the interaction relate through the so-called scattering matrix (S-matrix).

For example, at LEP,
e+
+
e

e+
+
e
, or
e+
+
e

μ+
+
μ
are processes where the initial state is an electron and a positron colliding to produce an electron and a positron or two muons of opposite charge: the final states. In these simple cases, no automatic packages are needed and cross-section analytical expressions can be easily derived at least for the lowest approximation: the Born approximation also called the leading order or the tree level (as Feynman diagrams have only trunk and branches, no loops).

But particle physics is now requiring much more complex calculations like at LHC where are protons and is the number of jets of particles initiated by proton constituents (quarks and gluons). The number of subprocesses describing a given process is so large that automatic tools have been developed to mitigate the burden of hand calculations.

Interactions at higher energies open a large spectrum of possible final states and consequently increase the number of processes to compute.

High precision experiments impose the calculation of higher order calculation, namely the inclusion of subprocesses where more than one virtual particle can be created and annihilated during the interaction lapse creating so-called loops which induce much more involved calculations.

Finally new theoretical models like the supersymmetry model (MSSM in its minimal version) predict a flurry of new processes.

The automatic packages, once seen as mere teaching support, have become, this last 10 years an essential component of the data simulation and analysis suite for all experiments. They help constructing event generators and are sometimes viewed as generators of event generators or Meta-generators.

A particle physics model is essentially described by its Lagrangian. To simulate the production of events through event generators, 3 steps have to be taken. The Automatic Calculation project is to create the tools to make those steps as automatic (or programmed) as possible:

I Feynman rules, coupling and mass generation

II Matrix element code generation: Various methods are used to automatically produce the matrix element expression in a computer language (Fortran, C/C++). They use values (i.e. for the masses) or expressions (i.e. for the couplings) produced by step I or model specific libraries constructed by hands (usually heavily relying on Computer algebra languages). When this expression is integrated (usually numerically) over the internal degrees of freedom it will provide the total and differential cross-sections for a given set of initial parameters like the initial state particle energies and polarization.

III Event generator code generation: This code must them be interfaced to other packages to fully provide the actual final state. The various effects or phenomenon that need to be implemeted are:

The interplay or matching of the precise matrix element calculation and the approximations resulting from the simulation of the parton shower gives rise to further complications, either within a given level of precision like at leading order (LO) for the production of n jets or between two levels of precision when tempting to connect matrix element computed at next-to-leading (NLO) (1-loop) or next-to-next-leading order (NNLO) (2-loops) with LO partons shower package.

Several methods have been developed for this matching, including: Subtraction methods.

But the only correct way is to match packages at the same level theoretical accuracy like the NLO matrix element calculation with NLO parton shower packages. This is currently in development.

History

The idea of automation of the calculations in high-energy physics is not new. It dates back to the 1960s when packages such as SCHOONSCHIP and then REDUCE had been developed.

These are symbolic manipulation codes that automatize the algebraic parts of a matrix element evaluation, like traces on Dirac matrices and contraction of Lorentz indices. Such codes have evolved quite a lot with applications not only optimized for high-energy physics like FORM but also more general purpose programs like Mathematica and Maple.

Generation of QED Feynman graphs at any order in the coupling constant was automatized in the late 70's[15]. One of the first major application of these early developments in this field was the calculation of the anomalous magnetic moments of the electron and the muon[16]. The first automatic system incorporating all the steps for the calculation of a cross section, from Feynman graph generation, amplitude generation through a REDUCE source code that produces a FORTRAN code, phase space integration and event generation with BASES/SPRING[17] is GRAND[18]. It was limited to tree-level processes in QED. In the early nineties, a few groups started to develop packages aiming at the automation in the SM[19]. [1] [2] [3] [4] [5] [6] [7] [8] [9] [10]

Matrix element calculation methods

Helicity amplitude

Feynman amplitudes are written in terms of spinor products of wave functions for massless fermions, and then evaluated numerically before the amplitudes are squared. Taking into account fermion masses implies that Feynman amplitudes are decomposed into vertex amplitudes by splitting the internal lines into wave function of fermions and polarization vectors of gauge bosons.

All helicity configuration can be computed independently.

Helicity amplitude squared

The method is similar to the previous one, but the numerical calculation is performed after squaring the Feynman Amplitude. The final expression is shorter and therefore faster to compute, but independent helicity information are not anymore available.

Dyson-Schwinger recursive equations

The scattering amplitude is evaluated recursively through a set of Dyson-Schwinger equations. The computational cost of this algorithm grows asymptotically as 3n, where n is the number of particles involved in the process, compared to n! in the traditional Feynman graphs approach. Unitary gauge is used and mass effects are available as well. Additionally, the color and helicity structures are appropriately transformed so the usual summation is replaced by the Monte Carlo techniques. [11]

Higher order calculations

[12]

Additional package for Event generation

The integration of the "matrix element" over the multidimensional internal parameters phase space provides the total and differential cross-sections. Each point of this phase space is associated to an event probability. This is used to randomly generate events closely mimicking experimental data. This is called event generation, the first step in the complete chain of event simulation. The initial and final state particles can be elementary particles like electrons, muons, or photons but also partons (protons and neutrons).

More effects must then be implemented to reproduce real life events as those detected at the colliders.

The initial electron or positron may undergo radiation before they actually interact: initial state radiation and beamstrahlung.

The bare partons that do not exist in nature (they are confined inside the hadrons) must be so to say dressed so that they form the known hadrons or mesons. They are made in two steps: parton shower and hadronization.

When the initial state particles are protons at high energy, it is only their constituents which interact. Therefore, the specific parton that will experience the "hard interaction" has to be selected. Structure functions must therefore be implemented. The other parton may interact "softly", and must also be simulated as they contribute to the complexity of the event: the underlying event.

Initial state radiation and beamstrahlung

Parton shower and Hadronization

At leading Order (LO)

At Next-to-Leading order (NLO)

Structure and Fragmentation Functions

The fragmentation function (F.F.) is a probability distribution function. It is used to find the density function of fragmented mesons in hadron -hadron collision.

The structure function, like the fragmentation function, is also a probability density function. It is analogous to the structure factor in solid-state physics.

Underlying event

Model specific packages

SM

MSSM

Automatic software packages can be useful in exploring a number of Beyond the Standard Model (BSM) theories, such as the Minimal Supersymmetric Standard Model (MSSM), to predict and understand possible particle interactions in future physics experiments.

Several computation issues need to be considered for automatic calculations. For example one scenario is the fact that special functions often need to be calculated in these software packages, both/either algebraically and/or numerically. For algebraic calculations, symbolic packages e.g. Maple, Mathematica often need to consider abstract, mathematical structures in subatomic particle collisions and emissions.

Multi-dimensional integrators

Ultra-High Precision Numerical computation

Existing Packages

Feynman rules generators

Tree Level Packages

NameModelMax FSTested FSShort descriptionPublicationMethodOutputStatus
MadGraph5 Any Model1/2->n2->8complete, massive, helicity, color, decay chain what is MG5 HA (automatic generation)Output PD
Grace SM/MSSM2->n2->6complete,massive,helicity,color Manual v2.0 HAOutput PD
CompHEP ModelMax FSTested FSShort descriptionPublicationmethodOutputStatus
CalcHEP ModelMax FSTested FSShort description Publication MethodOutputStatus
SherpaSM/MSSM2->n2->8massive publication HA/DSOutput PD
GenEvaModelMax FSTested FSShort descriptionPublicationMethodOutputStatus
HELAC ModelMax FSTested FSShort descriptionPublicationMethodOutputStatus
NameModelMax FSTested FSShort descriptionPublicationMethodOutputStatus

Status: PD: Public Domain,
Model: SM: Standard Model, MSSM: Minimal Supersymmetric Standard Model
Method: HA: Helicity Amplitude, DS: Dyson Schwinger
Output: ME: Matrix Element, CS: Cross-Sections, PEG: Parton level Event Generation, FEG: Full particle level Event Generation

Higher-order Packages

NameModelOrder testedMax FSTested FSShort descriptionPublicationMethodStatus
Grace L-1SM/MSSM1-loop2->n2->4complete,massive,helicity,colorNAMethodNA
NameOrderModelMax FSTested FSShort descriptionPublicationMethodStatus

Additional package for Event generation

Related Research Articles

<span class="mw-page-title-main">Quark</span> Elementary particle, main constituent of matter

A quark is a type of elementary particle and a fundamental constituent of matter. Quarks combine to form composite particles called hadrons, the most stable of which are protons and neutrons, the components of atomic nuclei. All commonly observable matter is composed of up quarks, down quarks and electrons. Owing to a phenomenon known as color confinement, quarks are never found in isolation; they can be found only within hadrons, which include baryons and mesons, or in quark–gluon plasmas. For this reason, much of what is known about quarks has been drawn from observations of hadrons.

<span class="mw-page-title-main">Quantum chromodynamics</span> Theory of the strong nuclear interactions

In theoretical physics, quantum chromodynamics (QCD) is the theory of the strong interaction between quarks mediated by gluons. Quarks are fundamental particles that make up composite hadrons such as the proton, neutron and pion. QCD is a type of quantum field theory called a non-abelian gauge theory, with symmetry group SU(3). The QCD analog of electric charge is a property called color. Gluons are the force carriers of the theory, just as photons are for the electromagnetic force in quantum electrodynamics. The theory is an important part of the Standard Model of particle physics. A large body of experimental evidence for QCD has been gathered over the years.

The up quark or u quark is the lightest of all quarks, a type of elementary particle, and a significant constituent of matter. It, along with the down quark, forms the neutrons and protons of atomic nuclei. It is part of the first generation of matter, has an electric charge of +2/3 e and a bare mass of 2.2+0.5
−0.4
 MeV/c2
. Like all quarks, the up quark is an elementary fermion with spin 1/2, and experiences all four fundamental interactions: gravitation, electromagnetism, weak interactions, and strong interactions. The antiparticle of the up quark is the up antiquark, which differs from it only in that some of its properties, such as charge have equal magnitude but opposite sign.

The down quark is a type of elementary particle, and a major constituent of matter. The down quark is the second-lightest of all quarks, and combines with other quarks to form composite particles called hadrons. Down quarks are most commonly found in atomic nuclei, where it combines with up quarks to form protons and neutrons. The proton is made of one down quark with two up quarks, and the neutron is made up of two down quarks with one up quark. Because they are found in every single known atom, down quarks are present in all everyday matter that we interact with.

<span class="mw-page-title-main">Minimal Supersymmetric Standard Model</span> Simplest supersymmetric extension to the Standard Model

The Minimal Supersymmetric Standard Model (MSSM) is an extension to the Standard Model that realizes supersymmetry. MSSM is the minimal supersymmetrical model as it considers only "the [minimum] number of new particle states and new interactions consistent with "Reality". Supersymmetry pairs bosons with fermions, so every Standard Model particle has a superpartner yet undiscovered. If discovered, such superparticles could be candidates for dark matter, and could provide evidence for grand unification or the viability of string theory. The failure to find evidence for MSSM using the Large Hadron Collider has strengthened an inclination to abandon it.

In quantum physics, Regge theory is the study of the analytic properties of scattering as a function of angular momentum, where the angular momentum is not restricted to be an integer multiple of ħ but is allowed to take any complex value. The nonrelativistic theory was developed by Tullio Regge in 1959.

In particle physics, flavour or flavor refers to the species of an elementary particle. The Standard Model counts six flavours of quarks and six flavours of leptons. They are conventionally parameterized with flavour quantum numbers that are assigned to all subatomic particles. They can also be described by some of the family symmetries proposed for the quark-lepton generations.

<span class="mw-page-title-main">Drell–Yan process</span> Process in high-energy hadron–hadron scattering

The Drell–Yan process occurs in high energy hadron–hadron scattering. It takes place when a quark of one hadron and an antiquark of another hadron annihilate, creating a virtual photon or Z boson which then decays into a pair of oppositely-charged leptons. Importantly, the energy of the colliding quark-antiquark pair can be almost entirely transformed into the mass of new particles. This process was first suggested by Sidney Drell and Tung-Mow Yan in 1970 to describe the production of lepton–antilepton pairs in high-energy hadron collisions. Experimentally, this process was first observed by J.H. Christenson et al. in proton–uranium collisions at the Alternating Gradient Synchrotron.

Event generators are software libraries that generate simulated high-energy particle physics events. They randomly generate events as those produced in particle accelerators, collider experiments or the early universe. Events come in different types called processes as discussed in the Automatic calculation of particle interaction or decay article.

In particle physics, the parton model is a model of hadrons, such as protons and neutrons, proposed by Richard Feynman. It is useful for interpreting the cascades of radiation produced from quantum chromodynamics (QCD) processes and interactions in high-energy particle collisions.

The Les Houches Accords are agreements between particle physicists to standardize the interface between the matrix element programs and the event generators used to calculate different quantities. The original accord was initially formed in 2001, at a conference in Les Houches, in the French Alps, before it was subsequently expanded.

<span class="mw-page-title-main">CompHEP</span> High energy physics computation software package

CompHEP is a software package for automatic computations in high energy physics from Lagrangians to collision events or particle decays.

S-matrix theory was a proposal for replacing local quantum field theory as the basic principle of elementary particle physics.

Computational particle physics refers to the methods and computing tools developed in and used by particle physics research. Like computational chemistry or computational biology, it is, for particle physics both a specific branch and an interdisciplinary field relying on computer science, theoretical and experimental particle physics and mathematics. The main fields of computational particle physics are: lattice field theory, automatic calculation of particle interaction or decay and event generators.

<span class="mw-page-title-main">Light front quantization</span> Technique in computational quantum field theory

The light-front quantization of quantum field theories provides a useful alternative to ordinary equal-time quantization. In particular, it can lead to a relativistic description of bound systems in terms of quantum-mechanical wave functions. The quantization is based on the choice of light-front coordinates, where plays the role of time and the corresponding spatial coordinate is . Here, is the ordinary time, is one Cartesian coordinate, and is the speed of light. The other two Cartesian coordinates, and , are untouched and often called transverse or perpendicular, denoted by symbols of the type . The choice of the frame of reference where the time and -axis are defined can be left unspecified in an exactly soluble relativistic theory, but in practical calculations some choices may be more suitable than others.

<span class="mw-page-title-main">Quark–gluon plasma</span> Phase of quantum chromodynamics (QCD)

Quark–gluon plasma is an interacting localized assembly of quarks and gluons at thermal and chemical (abundance) equilibrium. The word plasma signals that free color charges are allowed. In a 1987 summary, Léon van Hove pointed out the equivalence of the three terms: quark gluon plasma, quark matter and a new state of matter. Since the temperature is above the Hagedorn temperature—and thus above the scale of light u,d-quark mass—the pressure exhibits the relativistic Stefan-Boltzmann format governed by temperature to the fourth power and many practically massless quark and gluon constituents. It can be said that QGP emerges to be the new phase of strongly interacting matter which manifests its physical properties in terms of nearly free dynamics of practically massless gluons and quarks. Both quarks and gluons must be present in conditions near chemical (yield) equilibrium with their colour charge open for a new state of matter to be referred to as QGP.

In quantum field theory, initial and final state radiation refers to certain kinds of radiative emissions that are not due to particle annihilation. It is important in experimental and theoretical studies of interactions at particle colliders.

<span class="mw-page-title-main">Light-front quantization applications</span> Quantization procedure in quantum field theory

The light-front quantization of quantum field theories provides a useful alternative to ordinary equal-time quantization. In particular, it can lead to a relativistic description of bound systems in terms of quantum-mechanical wave functions. The quantization is based on the choice of light-front coordinates, where plays the role of time and the corresponding spatial coordinate is . Here, is the ordinary time, is a Cartesian coordinate, and is the speed of light. The other two Cartesian coordinates, and , are untouched and often called transverse or perpendicular, denoted by symbols of the type . The choice of the frame of reference where the time and -axis are defined can be left unspecified in an exactly soluble relativistic theory, but in practical calculations some choices may be more suitable than others. The basic formalism is discussed elsewhere.

<span class="mw-page-title-main">Zvi Bern</span> American theoretical particle physicist

Zvi Bern is an American theoretical particle physicist. He is a professor at University of California, Los Angeles (UCLA).

In quantum field theory, a sum rule is a relation between a static quantity and an integral over a dynamical quantity. Therefore, they have a form such as:

References

  1. Kaneko, T. (1990). "Automatic calculation of Feynman amplitudes". New computing techniques in physics research. p. 555. Archived from the original on 2012-12-11.
  2. Boos, E.E; et al. (1994). "Automatic calculation in high-energy physics by Grace/Chanel and CompHEP". International Journal of Modern Physics C . 5 (4): 615. Bibcode:1994IJMPC...5..615B. doi:10.1142/S0129183194000787.
  3. Wang, J.-X. (1993). "Automatic calculation of Feynman loop-diagrams I. Generation of a simplified form of the amplitude". Computer Physics Communications . 77 (2): 263. Bibcode:1993CoPhC..77..263W. doi:10.1016/0010-4655(93)90010-A.
  4. Kaneko, T.; Nakazawa, N. (1995). "Automatic calculation of two loop weak corrections to muon anomalous magnetic moment". New computing techniques in physics research. p. 173. arXiv: hep-ph/9505278 . Bibcode:1995hep.ph....5278K. Archived from the original on 2012-12-10.
  5. Jimbo, M.; (Minami-Tateya Collaboration); et al. (1995). "Automatic calculation of SUSY particle production". High energy physics and quantum field theory. p. 155. arXiv: hep-ph/9605414 . Bibcode:1996hep.ph....5414J.
  6. Franzkowski, J. (1997). "Automatic calculation of massive two-loop self-energies with XLOOPS". Nuclear Instruments and Methods in Physics Research A . 389 (1–2): 339–342. arXiv: hep-ph/9611378 . Bibcode:1997NIMPA.389..333F. doi:10.1016/S0168-9002(97)00121-6.
  7. Brucher, L. (2000). "Automatic Feynman diagram calculation with xloops: A Short overview". arXiv: hep-ph/0002028 .
  8. Perret-Gallix, D. (1999). "Automatic amplitude calculation and event generation for collider physics: GRACE and CompHEP". High energy physics and quantum field theory. p. 270. Archived from the original on 2012-12-11.
  9. Belanger, G.; et al. (2006). "Automatic calculations in high energy physics and GRACE at one-loop". Physics Reports . 430 (3): 117–209. arXiv: hep-ph/0308080 . Bibcode:2006PhR...430..117B. doi:10.1016/j.physrep.2006.02.001. S2CID   7049291.
  10. Fujimoto, J.; et al. (2004). "Automatic one-loop calculation of MSSM processes with GRACE". Nuclear Instruments and Methods in Physics Research A . 534 (1–2): 246. arXiv: hep-ph/0402145 . Bibcode:2004NIMPA.534..246F. doi:10.1016/j.nima.2004.07.095. S2CID   7717301.
  11. Kanaki, A.; Papadopoulos, C.G. (2000). "HELAC: A Package to compute electroweak helicity amplitudes". Computer Physics Communications . 132 (3): 306–315. arXiv: hep-ph/0002082 . Bibcode:2000CoPhC.132..306K. doi:10.1016/S0010-4655(00)00151-X. S2CID   14533093.
  12. Belanger, G.; et al. (2006). "Automatic calculations in high energy physics and Grace at one-loop". Physics Reports . 430 (3): 117–209. arXiv: hep-ph/0308080 . Bibcode:2006PhR...430..117B. doi:10.1016/j.physrep.2006.02.001. S2CID   7049291.