Accelerator physics codes

Last updated

A charged particle accelerator is a complex machine that takes elementary charged particles and accelerates them to very high energies. Accelerator physics is a field of physics encompassing all the aspects required to design and operate the equipment and to understand the resulting dynamics of the charged particles. There are software packages associated with each domain. The 1990 edition of the Los Alamos Accelerator Code Group's compendium [1] provides summaries of more than 200 codes. Certain codes are still in use today, although many are obsolete. Another index of existing and historical accelerator simulation codes is located at the CERN CARE/HHH website. [2]

Contents

Single particle dynamics codes

For many applications it is sufficient to track a single particle through the relevant electric and magnetic fields. Old codes no longer maintained by their original authors or home institutions include: BETA, [3] AGS, ALIGN, COMFORT, DESIGN, DIMAD, HARMON, LEGO, LIAR, MAGIC, MARYLIE, PATRICIA, PETROS, RACETRACK, SYNCH, [4] TRANSPORT, TURTLE, and UAL. Some legacy codes are maintained by commercial organizations for academic, industrial and medical accelerator facilities that continue to use those codes. TRANSPORT, TRACE 3-D and TURTLE are among the historic codes that are commercially maintained. [5]

Major maintained codes include:

Single Particle DynamicsSpin TrackingTaylor MapsWeak-Strong Beam-Beam InteractionElectromagnetic Field TrackingHigher Energy Collective EffectsSynchrotron Radiation TrackingWakefieldsExtensibleNotes
Accelerator Toolbox (AT), [6] YesYes [7] NoNoNoYesNoNoYes
ASTRA [8] YesNoNoNoYesYesNoYesNoFor space-charge effects evaluation
BDSIM [9] YesNoNoNoYesNoNoNoYesFor particle-matter interaction studies.
Bmad (contains PTC) [10] YesYesYesYesYesYesYesYesYesReproduces PTC's unique beam line structures. Also simulates X-rays.
COSY INFINITY [11] YesYesYesNoYesNoNoNoNo
DYNAC [12] YesNoNoNoNoNoNoNoNo
Elegant [13] YesNoNoNoYesYesNoYesNo
MAD8 and MAD-X (includes PTC) [14] YesNoYesYesNoNoYesNoNo
MAD-NG [14] YesNoYesYesNoNoYesNoYesVery extensible, embeds LuaJIT
MERLIN++ [15] [16] YesYesNoNoNoNoNoYesYesOther: beam-matter interactions, sliced-macroparticle tracking
OCELOT [17] YesNoNoNoNoYesYesYesYes
OPA [18] YesNoNoNoNoNoNoNoNo
OPAL [19] YesNoYesNoYesYesNoYesYesOpen source, runs on the laptop and on x 10k cores.
PLACET [20] YesNoNoNoNoYesYesYesYesSimulates a LINAC including wakefields.
Propaga [21] YesNoNoNoNoNoNoNoYes
PTC [22] YesYesYesYesYesNoNoNoYes
SAD [23] YesNoNoYesNoYesYesYesNo
SAMM [24] YesYesNoNoNoNoNoNoNo
SixTrack [25] YesNoYesYesNoNoNoNoNoCan run on BOINC
Zgoubi [26] [27] YesYesNoNoYesNoYesNoYesOpen source.

Columns

Spin Tracking
Tracking of a particle's spin.
Taylor Maps
Construction of Taylor series maps to high order that can be used for simulating particle motion and also can be used for such things as extracting single particle resonance strengths.
Weak-Strong Beam-Beam Interaction
Can simulate the beam-beam interaction with the simplification that one beam is essentially fixed in size. See below for a list of strong-strong interaction codes.
Electromagnetic Field Tracking
Can track (ray trace) a particle through arbitrary electromagnetic fields.
Higher Energy Collective effects
The interactions between the particles in the beam can have important effects on the behavior, control and dynamics. Collective effects take different forms from Intrabeam Scattering (IBS) which is a direct particle-particle interaction to wakefields which are mediated by the vacuum chamber wall of the machine the particles are traveling in. In general, the effect of direct particle-particle interactions is less with higher energy particle beams. At very low energies, space charge has a large effect on a particle beam and thus becomes hard to calculate. See below for a list of programs that can handle low energy space charge forces.
Synchrotron radiation tracking
Ability to track the synchrotron radiation (mainly X-rays) produced by the acceleration of charged particles.
Wakefields
The electro-magnetic interaction between the beam and the vacuum chamber wall enclosing the beam are known as wakefields. Wakefields produce forces that affect the trajectory of the particles of the beam and can potentially destabilize the trajectories.
Extensible
Open source and object oriented coding to make it relatively easy to extend the capabilities.

Space charge codes

The self interaction (e.g. space charge) of the charged particle beam can cause growth of the beam, such as with bunch lengthening, or intrabeam scattering. Additionally, space charge effects may cause instabilities and associated beam loss. Typically, at relatively low energies (roughly for energies where the relativistic gamma factor is less than 10 or so), the Poisson equation is solved at intervals during the tracking using particle-in-cell algorithms. Space charge effects lessen at higher energies so at higher energies the space charge effects may be modeled using simpler algorithms that are computationally much faster than the algorithms used at lower energies. Codes that handle low energy space charge effects include:

At higher energies, space charge effects include Touschek scattering and coherent synchrotron radiation (CSR). Codes that handle higher energy space charge include:

"Strong-strong" beam-beam effects codes

When two beams collide, the electromagnetic field of one beam will then have strong effects on the other one, called beam-beam effects. So called "weak-strong" simulations model one beam (called the "strong" beam since it affects the other beam) as a fixed distribution (typically a Gaussian distribution) which interacts with the particles of the other "weak" beam. This greatly simplifies the simulation. A full "strong-strong" simulation is more complicated and takes more simulation time. Strong-strong codes include

Impedance computation codes

An important class of collective effects may be summarized in terms of the beams response to an "impedance". An important job is thus the computation of this impedance for the machine. Codes for this computation include

Magnet and other hardware-modeling codes

To control the charged particle beam, appropriate electric and magnetic fields must be created. There are software packages to help in the design and understanding of the magnets, RF cavities, and other elements that create these fields. Codes include

Lattice description and data interchange issues

Given the variety of modeling tasks, there is not one common data format that has developed. For describing the layout of an accelerator and the corresponding elements, one uses a so-called "lattice file". There have been numerous attempts at unifying the lattice file formats used in different codes. One unification attempt is the Accelerator Markup Language, and the Universal Accelerator Parser. [52] Another attempt at a unified approach to accelerator codes is the UAL or Universal Accelerator Library. [53] As of 2023 neither of these formats are maintained.

The file formats used in MAD may be the most common, with translation routines available to convert to an input form needed for a different code. Associated with the Elegant code is a data format called SDDS, with an associated suite of tools. If one uses a Matlab-based code, such as Accelerator Toolbox, one has available all the tools within Matlab.

For the interchange of particle positions and electromagnetic fields, the OpenPMD [54] standard defines a format which can then be implemented with a file format like HDF5.

Codes in applications of particle accelerators

There are many applications of particle accelerators. For example, two important applications are elementary particle physics and synchrotron radiation production. When performing a modeling task for any accelerator operation, the results of charged particle beam dynamics simulations must feed into the associated application. Thus, for a full simulation, one must include the codes in associated applications. For particle physics, the simulation may be continued in a detector with a code such as Geant4.

For a synchrotron radiation facility, for example, the electron beam produces an x-ray beam that then travels down a beamline before reaching the experiment. Thus, the electron beam modeling software must interface with the x-ray optics modelling software such as SRW, [55] Shadow, [56] McXTrace, [57] or Spectra. [58] Bmad [10] can model both X-rays and charged particle beams. The x-rays are used in an experiment which may be modeled and analyzed with various software, such as the DAWN science platform. [59] OCELOT [60] also includes both synchrotron radiation calculation and x-ray propagation models.

Industrial and medical accelerators represent another area of important applications. A 2013 survey estimated that there were about 27,000 industrial accelerators and another 14,000 medical accelerators world wide, [61] and those numbers have continued to increase since that time. [62] Codes used at those facilities vary considerably and often include a mix of traditional codes and custom codes developed for specific applications. The Advanced Orbit Code (AOC) [63] developed at Ion Beam Applications is an example.

See also

Related Research Articles

<span class="mw-page-title-main">Particle physics</span> Study of subatomic particles and forces

Particle physics or high-energy physics is the study of fundamental particles and forces that constitute matter and radiation. The field also studies combinations of elementary particles up to the scale of protons and neutrons, while the study of combination of protons and neutrons is called nuclear physics.

<span class="mw-page-title-main">Cyclotron</span> Type of particle accelerator

A cyclotron is a type of particle accelerator invented by Ernest Lawrence in 1929–1930 at the University of California, Berkeley, and patented in 1932. A cyclotron accelerates charged particles outwards from the center of a flat cylindrical vacuum chamber along a spiral path. The particles are held to a spiral trajectory by a static magnetic field and accelerated by a rapidly varying electric field. Lawrence was awarded the 1939 Nobel Prize in Physics for this invention.

<span class="mw-page-title-main">Synchrotron radiation</span> Electromagnetic radiation emitted by charged particles accelerated perpendicular to their velocity

Synchrotron radiation is the electromagnetic radiation emitted when relativistically charged particles are subject to an acceleration perpendicular to their velocity. It is produced artificially in some types of particle accelerators or naturally by fast electrons moving through magnetic fields. The radiation produced in this way has a characteristic polarization, and the frequencies generated can range over a large portion of the electromagnetic spectrum.

Accelerator physics is a branch of applied physics, concerned with designing, building and operating particle accelerators. As such, it can be described as the study of motion, manipulation and observation of relativistic charged particle beams and their interaction with accelerator structures by electromagnetic fields.

<span class="mw-page-title-main">ATLAS experiment</span> CERN LHC experiment

ATLAS is the largest general-purpose particle detector experiment at the Large Hadron Collider (LHC), a particle accelerator at CERN in Switzerland. The experiment is designed to take advantage of the unprecedented energy available at the LHC and observe phenomena that involve highly massive particles which were not observable using earlier lower-energy accelerators. ATLAS was one of the two LHC experiments involved in the discovery of the Higgs boson in July 2012. It was also designed to search for evidence of theories of particle physics beyond the Standard Model.

<span class="mw-page-title-main">Beamline</span> Trajectory of a beam of accelerated particles

In accelerator physics, a beamline refers to the trajectory of the beam of particles, including the overall construction of the path segment along a specific path of an accelerator facility. This part is either

<span class="mw-page-title-main">Dipole magnet</span> Simplest type of magnet

A dipole magnet is the simplest type of magnet. It has two poles, one north and one south. Its magnetic field lines form simple closed loops which emerge from the north pole, re-enter at the south pole, then pass through the body of the magnet. The simplest example of a dipole magnet is a bar magnet.

<span class="mw-page-title-main">Large Electron–Positron Collider</span> Particle accelerator at CERN, Switzerland

The Large Electron–Positron Collider (LEP) was one of the largest particle accelerators ever constructed. It was built at CERN, a multi-national centre for research in nuclear and particle physics near Geneva, Switzerland.

In plasma physics, the particle-in-cell (PIC) method refers to a technique used to solve a certain class of partial differential equations. In this method, individual particles in a Lagrangian frame are tracked in continuous phase space, whereas moments of the distribution such as densities and currents are computed simultaneously on Eulerian (stationary) mesh points.

<span class="mw-page-title-main">Synchrotron</span> Type of cyclic particle accelerator

A synchrotron is a particular type of cyclic particle accelerator, descended from the cyclotron, in which the accelerating particle beam travels around a fixed closed-loop path. The magnetic field which bends the particle beam into its closed path increases with time during the accelerating process, being synchronized to the increasing kinetic energy of the particles. The synchrotron is one of the first accelerator concepts to enable the construction of large-scale facilities, since bending, beam focusing and acceleration can be separated into different components. The most powerful modern particle accelerators use versions of the synchrotron design. The largest synchrotron-type accelerator, also the largest particle accelerator in the world, is the 27-kilometre-circumference (17 mi) Large Hadron Collider (LHC) near Geneva, Switzerland, built in 2008 by the European Organization for Nuclear Research (CERN). It can accelerate beams of protons to an energy of 13 tera electronvolts (TeV or 1012 eV).

<span class="mw-page-title-main">Geant4</span> Scientific software for particle physics

Geant4 is a platform for "the simulation of the passage of particles through matter" using Monte Carlo methods. It is the successor of the GEANT series of software toolkits developed by The Geant4 Collaboration, and the first to use object oriented programming. Its development, maintenance and user support are taken care by the international Geant4 Collaboration. Application areas include high energy physics and nuclear experiments, accelerator and space physics studies. The software is used by a number of research projects around the world.

<span class="mw-page-title-main">Gargamelle</span> CERN Bubble chamber particle detector

Gargamelle was a heavy liquid bubble chamber detector in operation at CERN between 1970 and 1979. It was designed to detect neutrinos and antineutrinos, which were produced with a beam from the Proton Synchrotron (PS) between 1970 and 1976, before the detector was moved to the Super Proton Synchrotron (SPS). In 1979 an irreparable crack was discovered in the bubble chamber, and the detector was decommissioned. It is currently part of the "Microcosm" exhibition at CERN, open to the public.

A charged particle beam is a spatially localized group of electrically charged particles that have approximately the same position, kinetic energy, and direction. The kinetic energies of the particles are much larger than the energies of particles at ambient temperature. The high energy and directionality of charged particle beams make them useful for many applications in particle physics.

<span class="mw-page-title-main">AWAKE</span>

The AWAKE facility at CERN is a proof-of-principle experiment, which investigates wakefield plasma acceleration using a proton bunch as a driver, a world-wide first. It aims to accelerate a low-energy witness bunch of electrons from 15 to 20 MeV to several GeV over a short distance by creating a high acceleration gradient of several GV/m. Particle accelerators currently in use, like CERN's LHC, use standard or superconductive RF-cavities for acceleration, but they are limited to an acceleration gradient in the order of 100 MV/m.

<span class="mw-page-title-main">Bevatron</span> Particle accelerator at Lawrence Berkeley National Laboratory

The Bevatron was a particle accelerator — specifically, a weak-focusing proton synchrotron — at Lawrence Berkeley National Laboratory, U.S., which began operating in 1954. The antiproton was discovered there in 1955, resulting in the 1959 Nobel Prize in physics for Emilio Segrè and Owen Chamberlain. It accelerated protons into a fixed target, and was named for its ability to impart energies of billions of eV.

The EGS computer code system is a general purpose package for the Monte Carlo simulation of the coupled transport of electrons and photons in an arbitrary geometry for particles with energies from a few keV up to several hundreds of GeV. It originated at SLAC but National Research Council of Canada and KEK have been involved in its development since the early 80s.

Dassault Systèmes Simulia Corp. is a computer-aided engineering (CAE) vendor. Formerly known as Abaqus Inc. and previously Hibbitt, Karlsson & Sorensen, Inc., (HKS), the company was founded in 1978 by David Hibbitt, Bengt Karlsson and Paul Sorensen, and has its headquarters in Providence, Rhode Island.

<span class="mw-page-title-main">Particle accelerator</span> Research apparatus for particle physics

A particle accelerator is a machine that uses electromagnetic fields to propel charged particles to very high speeds and energies, and to contain them in well-defined beams.

<span class="mw-page-title-main">Storage ring</span> Type of particle accelerator

A storage ring is a type of circular particle accelerator in which a continuous or pulsed particle beam may be kept circulating typically for many hours. Storage of a particular particle depends upon the mass, momentum and usually the charge of the particle to be stored. Storage rings most commonly store electrons, positrons, or protons.

<span class="mw-page-title-main">NA63 experiment</span>

The NA63 experiment aims to study the radiation process in strong electromagnetic fields. Located at CERN, in the North Area. It is a fixed-target experiment which uses the H4 secondary electron beams from the SPS, which are directed onto different targets. Those are made from a variety of elements, ranging from the relatively light carbon and silicon, through the heavier iron and tin to tungsten, gold and lead and are either amorphous or mono-crystals.

References

  1. ↑ Computer Codes for Particle Accelerator Design and Analysis: A Compendium, Second Edition, Helen Stokes Deaven and Kwok Chi Dominic Chen, Los Alamos National Laboratory report number LA-UR-90-1766, 290 pages (1990).
  2. ↑ the CERN CARE/HHH website Archived December 13, 2012, at the Wayback Machine
  3. ↑ "user's guide" (PDF). Archived from the original (PDF) on 2018-07-04. Retrieved 2013-11-16.
  4. ↑ libtracy at sourceforge.net
  5. ↑ AccelSoft Inc. website
  6. ↑ ATcollab website
  7. ↑ See https://github.com/carmignani/festa
  8. 1 2 ASTRA Homepage
  9. ↑ BDSIM Homepage
  10. 1 2 Bmad home page
  11. ↑ "COSY". Archived from the original on 2018-07-07. Retrieved 2018-06-14.
  12. ↑ "DYNAC".
  13. ↑ ELEGANT, a Flexible SDDS Compliant Code for Accelerator Simulation software
  14. 1 2 "MAD - Methodical Accelerator Design". mad@cern.ch. Retrieved 2020-09-09.
  15. ↑ Appleby, Robert; Barlow, Roger J.; Bungau, Adriana; Fallon, James; Kruecker, Dirk; Molson, James; Rafique, Haroon; Rowan, Scott; Serluca, Maurizio; Sjøbæk, Kyrre Ness; Toader, Adina; Tygier, Sam; Walker, Nick; Wolski, Andy (2019-03-03). "Github Merlin-Collaboration/Merlin". doi:10.5281/zenodo.2598428.{{cite journal}}: Cite journal requires |journal= (help)
  16. ↑ Appleby, Robert; Barlow, Roger J.; Bungau, Adriana; Fallon, James; Kruecker, Dirk; Molson, James; Rafique, Haroon; Rowan, Scott; Serluca, Maurizio; Sjøbæk, Kyrre Ness; Toader, Adina; Tygier, Sam; Walker, Nick; Wolski, Andy (2019). "Merlin++". doi:10.5281/zenodo.2598428.{{cite journal}}: Cite journal requires |journal= (help)
  17. ↑ OCELOT collaboration on GitHub
  18. ↑ OPA website
  19. ↑ "Home · Wiki · OPAL / SRC".
  20. ↑ PLACET manual
  21. ↑ Propaga GitHub repository
  22. ↑ "GitHub - jceepf/fpp_book". GitHub . 2019-02-06.
  23. ↑ SAD home page at kek.jp
  24. ↑ SAMM, another Matlab based tracking code, at liv.ac.uk
  25. ↑ SixTrack home page at cern.ch
  26. ↑ Zgoubi home page at sourceforge.net
  27. ↑ Zgoubi Users' Guide
  28. ↑ PIC solver at cst.com
  29. ↑ General Particle Tracer (GPT) from Pulsar Physics
  30. ↑ "IMPACT homepage at Berkeley Lab". Archived from the original on 2015-04-16. Retrieved 2015-04-09.
  31. ↑ ImpactX: an s-based beam dynamics code including space charge effects from Berkeley Lab
  32. ↑ THE MULTIPARTICLE TRACKING CODES SBTRACK AND MBTRACK. R. Nagaoka, PAC '09 paper here
  33. ↑ ORBIT home page at ornl.gov
  34. ↑ PyORBIT Collaboration
  35. ↑ OPAL homepage
  36. ↑ PyHEADTAIL wiki
  37. ↑ Synergia home page at fnal.gov
  38. ↑ TraceWin at CEA Saclay
  39. ↑ TRANFT user's manual, BNL--77074-2006-IR http://www.osti.gov/scitech/biblio/896444
  40. 1 2 3 VSim at Tech-X
  41. ↑ Warp wiki
  42. ↑ "GUINEA-PIG Twiki". twiki.cern.ch. Archived from the original on 2022-01-20. Retrieved 2020-07-03.
  43. ↑ "BeamBeam3D GitHub Repo". "J. Qiang, M. Furman, and R. Ryne, "A parallel particle-in-cell model for beam–beam interaction in high energy ring colliders"". J. Comp. Phys. 2004. doi:10.1016/j.jcp.2004.01.008.
  44. ↑ ABCI home page at kek.jp
  45. 1 2 ACE3P at slac.stanford.gov
  46. ↑ CST Archived 2018-07-29 at the Wayback Machine , Computer Simulation Technology at cst.com
  47. ↑ GdfidL, Gitter drueber, fertig ist die Laube at gdfidl.de
  48. ↑ T. Weiland, DESY
  49. ↑ COMSOL home page at comsol.com
  50. ↑ CST Electromagnetic Studio at cst.com [ permanent dead link ]
  51. ↑ "OPERA at magnet-design-software.com". Archived from the original on 2013-12-24. Retrieved 2013-11-15.
  52. ↑ Description of AML and UAP at cornell.edu
  53. ↑ See references by N. Malitsky and Talman such as this manual from 2002.
  54. ↑ OpenPMD GitHub repository.
  55. ↑ SRW home page at esrf.eu
  56. ↑ Shadow home page at esrf.eu
  57. ↑ McXTrace home page at mcxtrace.org
  58. ↑ "Spectra home page at riken.go.jp". Archived from the original on 2013-08-27. Retrieved 2013-11-15.
  59. ↑ DAWN science platform website
  60. ↑ "An Introduction to Ocelot". GitHub . 16 December 2021.
  61. ↑ R. W. Hamm and M. E. Hamm, Industrial Accelerators
  62. ↑ session on accelerator business opportunities at IPAC-17
  63. ↑ AOC, A Beam Dynamics Code for Medical and Industrial Accelerators
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.