Automatic calculation of particle interaction or decay

Automatic calculation of particle interaction or decay

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.

Particle accelerator or colliders produce collisions (interactions) of particle (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 expression 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 p p -> n jets where p are protons and n jets means n (2,3,4,...) 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 particles 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 sometime viewed as "generators of event generators" or "Meta-generators".__TOC__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;* LanHEP is an example of Feynman rules generation.;* Some model needs an additional step to compute, based on some parameters, the mass and coupling of new predicted particles.

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:;* Initial state radiation and beamstrahlung for e+ e- initial states.;* Parton Distribution Functions describing the actual content in terms of gluons and quarks of the p or p-bar initial state particles;* Parton showering describing the way final state quarks or gluons due to the QCD confinement generate additional quark/gluon pairs generating a so-called shower of partons before transforming into hadrons.;* Hadronization describing how the final quark pairs/triplets form the visible and detectable hadrons.;* Underlying event takes care of the way the rest, in term of constituent, of the initial protons also contribute to any given event.

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:
* Substraction 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

" in construction"

The idea of automation of the calculations in high-energy physics is not new. It dates back to the 1960's 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] . [ [http://www.slac.stanford.edu/spires/find/hep?key=5471150 Automatic calculation of Feynman amplitudes.] By T. Kaneko. In *Lyon 1990, New computing techniques in physics research* 555-56.] [ [http://www.slac.stanford.edu/spires/find/hep?key=2574861 Automatic calculation in high-energy physics by Grace/Chanel and CompHEP.] By E.E. Boos, et al., Int.J.Mod.Phys.C5:615-628,1994.] [ [http://www.slac.stanford.edu/spires/find/hep?key=2719657 Automatic calculation of Feynman loop diagram. 1. Generation of simplified form of amplitude.] By Jian-Xiong Wang. Comput.Phys.Commun.77:263-285,1993.] [ [http://www.slac.stanford.edu/spires/find/hep?key=3147525 Automatic calculation of two loop weak corrections to muon anomalous magnetic moment.] By Toshiaki Kaneko, Nobuya Nakazawa.
In *Pisa 1995, New computing techniques in physics research* 173-17. [hep-ph/9505278]
] [ [http://www.slac.stanford.edu/spires/find/hep?key=3360695 Automatic calculation of SUSY particle production.] By Minami-Tateya Collaboration
In *Zvenigorod 1995, High energy physics and quantum field theory* 155-16. [hep-ph/9605414]
] [ [http://www.slac.stanford.edu/spires/find/hep?key=3461025 Automatic calculation of massive two loop self-energies with XLOOPS. ] By Johannes Franzkowski. Nucl.Instrum.Meth.A389:333-338,1997. [hep-ph/9611378] ] [ [http://www.slac.stanford.edu/spires/find/hep?key=4306120 Automatic Feynman diagram calculation with xloops: A Short overview.] By Lars Brucher. [hep-ph/0002028] MZ-TH-98-18 (Mar 1998) 5p.] [ [http://www.slac.stanford.edu/spires/find/hep?key=4515617 Automatic amplitude calculation and event generation for collider physics: GRACE and CompHEP.] By D. Perret-Gallix. In *Moscow 1999, High energy physics and quantum field theory* 270-28.] [ [http://www.slac.stanford.edu/spires/find/hep?key=5631408 Automatic calculations in high energy physics and Grace at one-loop.] By G. Belanger, et al.,
Phys.Rept.430:117-209,2006. [hep-ph/0308080]
] [ [http://www.slac.stanford.edu/spires/find/hep?key=5795915 Automatic one loop calculation of MSSM processes with GRACE.] By J. Fujimoto, T. Ishikawa, M. Jimbo, T. Kon, M. Kuroda. Nucl.Instrum.Meth.A534:246-249,2004. [hep-ph/0402145] ]

Matrix element calculation methods

see Feynman diagram

Helicity amplitude

Feynman amplitudes are written in terms of spinor products of wave functions for massless fermions, andthen evaluated numerically before the amplitudes are squared. Taking into account fermion masses implies that Feynman amplitudes are decomposed into vertexamplitudes 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 [Schwinger-Dyson] [Dyson-Schwinger] equations. The computationalcost 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 availableas well. Additionally, the color and helicity structures are appropriately transformed so the usual summation isreplaced by the Monte Carlo techniques. [HELAC: A Package to compute electroweak helicity amplitudes.Aggeliki Kanaki, Costas G. Papadopoulos (Democritos Nucl. Res. Ctr.) . DEMO-HEP-2000-01, Feb 2000. 14pp.Published in Comput.Phys.Commun.132:306-315,2000.e-Print: hep-ph/0002082] .

Higher order calculations

" in construction " [Automatic calculations in high energy physics and Grace at one-loop.G. Belanger, F. Boudjema (Annecy, LAPP) , J. Fujimoto, T.Ishikawa, T. Kaneko (KEK, Tsukuba) , K. Kato (Kogakuin U.) ,Y. Shimizu (KEK, Tsukuba), LAPTH-982-03, KEK-CP-138, Aug 2003. 83pp.Published in Phys.Rept.430:117-209,2006.e-Print: hep-ph/0308080URL for hep-ph http://arxiv.org/abs/hep-ph/0308080]

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 electron, muon or photons but also partons, the constituent of the proton (quarks and gluons).

More effects must them 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 (there are confined insite the hadrons) must be so to say dressed so that they form the known hadrons or mesons. they is made in two steps: Parton shower and Hadronization.

When the initial state are protons, at high energy it is only their constituant 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" must be also be simulated as they contribute to the complexity of the event: Underlying event

Initial state radiation and beamstrahlung

("to be written")

Parton shower and Hadronization

("to be written")

At leading Order (LO)

("to be written")

At Next-to-Leading order (NLO)

("to be written")

Structure and Fragmentation Functions

("to be written")

Underlying event

("to be written")

Model specific packages

("to be written")

MSSM

Related computational issues

("to be written")

Multi-dimensional integrators

("to be written")

Ultra-High Precision Numerical computation

("to be written")

Existing Packages

Feynman rules generators

* LanHEP
* ...

Tree Level Packages

"PD": Public Domain, "SM": Standard Model, "MSSM": Minimal Supersymmetric Standard Model"HA": Helicity Amplitude"DS": Dyson Schwinger

Higher-order Packages

Additional package for Event generation

References


Wikimedia Foundation. 2010.

Игры ⚽ Поможем сделать НИР

Look at other dictionaries:

  • 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… …   Wikipedia

  • Event generator — Event generators are software libraries that generate simulated high energy particle physics events [ [http://arjournals.annualreviews.org/doi/abs/10.1146/annurev.nucl.55.090704.151505 M. L. Mangano T. J. Stelzer, Ann. Rev. Nucl. Part. Sci. 55,… …   Wikipedia

  • Event generators — are software packages that randomly generate particle physics events produced in particle accelerators, collider experiments, or during the initial phases of the creation of the Universe (Big Bang).Events come in different types called processes… …   Wikipedia

  • cosmos — /koz meuhs, mohs/, n., pl. cosmos, cosmoses for 2, 4. 1. the world or universe regarded as an orderly, harmonious system. 2. a complete, orderly, harmonious system. 3. order; harmony. 4. any composite plant of the genus Cosmos, of tropical… …   Universalium

  • electromagnetic radiation — Physics. radiation consisting of electromagnetic waves, including radio waves, infrared, visible light, ultraviolet, x rays, and gamma rays. [1950 55] * * * Energy propagated through free space or through a material medium in the form of… …   Universalium

  • Europe, history of — Introduction       history of European peoples and cultures from prehistoric times to the present. Europe is a more ambiguous term than most geographic expressions. Its etymology is doubtful, as is the physical extent of the area it designates.… …   Universalium

  • Supernova — This article is about the astronomical event. For other uses, see Supernova (disambiguation). Multiwavelength X ray, infrared, and optical compilation image of Kepler s supernova remnant, SN 1604. A supernova is a stellar explosion that is more… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”