-
Next Generation Computational Tools for the Modeling and Design of Particle Accelerators at Exascale
Authors:
Axel Huebl,
Remi Lehe,
Chad E. Mitchell,
Ji Qiang,
Robert D. Ryne,
Ryan T. Sandberg,
Jean-Luc Vay
Abstract:
Particle accelerators are among the largest, most complex devices. To meet the challenges of increasing energy, intensity, accuracy, compactness, complexity and efficiency, increasingly sophisticated computational tools are required for their design and optimization. It is key that contemporary software take advantage of the latest advances in computer hardware and scientific software engineering…
▽ More
Particle accelerators are among the largest, most complex devices. To meet the challenges of increasing energy, intensity, accuracy, compactness, complexity and efficiency, increasingly sophisticated computational tools are required for their design and optimization. It is key that contemporary software take advantage of the latest advances in computer hardware and scientific software engineering practices, delivering speed, reproducibility and feature composability for the aforementioned challenges. A new open source software stack is being developed at the heart of the Beam pLasma Accelerator Simulation Toolkit (BLAST) by LBNL and collaborators, providing new particle-in-cell modeling codes capable of exploiting the power of GPUs on Exascale supercomputers. Combined with advanced numerical techniques, such as mesh-refinement, and intrinsic support for machine learning, these codes are primed to provide ultrafast to ultraprecise modeling for future accelerator design and operations.
△ Less
Submitted 9 August, 2022; v1 submitted 3 August, 2022;
originally announced August 2022.
-
Using Kernel-Based Statistical Distance to Study the Dynamics of Charged Particle Beams in Particle-Based Simulation Codes
Authors:
Chad E. Mitchell,
Robert D. Ryne,
Kilean Hwang
Abstract:
Measures of discrepancy between probability distributions (statistical distance) are widely used in the fields of artificial intelligence and machine learning. We describe how certain measures of statistical distance can be implemented as numerical diagnostics for simulations involving charged-particle beams. Related measures of statistical dependence are also described. The resulting diagnostics…
▽ More
Measures of discrepancy between probability distributions (statistical distance) are widely used in the fields of artificial intelligence and machine learning. We describe how certain measures of statistical distance can be implemented as numerical diagnostics for simulations involving charged-particle beams. Related measures of statistical dependence are also described. The resulting diagnostics provide sensitive measures of dynamical processes important for beams in nonlinear or high-intensity systems, which are otherwise difficult to characterize. The focus is on kernel-based methods such as Maximum Mean Discrepancy, which have a well-developed mathematical foundation and reasonable computational complexity. Several benchmark problems and examples involving intense beams are discussed. While the focus is on charged-particle beams, these methods may also be applied to other many-body systems such as plasmas or gravitational systems.
△ Less
Submitted 8 April, 2022;
originally announced April 2022.
-
Simulations of Future Particle Accelerators: Issues and Mitigations
Authors:
D. Sagan,
M. Berz,
N. M. Cook,
Y. Hao,
G. Hoffstaetter,
A. Huebl,
C. -K. Huang,
M. H. Langston,
C. E. Mayes,
C. E. Mitchell,
C. -K. Ng,
J. Qiang,
R. D. Ryne,
A. Scheinker,
E. Stern,
J. -L. Vay,
D. Winklehner,
H. Zhang
Abstract:
The ever increasing demands placed upon machine performance have resulted in the need for more comprehensive particle accelerator modeling. Computer simulations are key to the success of particle accelerators. Many aspects of particle accelerators rely on computer modeling at some point, sometimes requiring complex simulation tools and massively parallel supercomputing. Examples include the modeli…
▽ More
The ever increasing demands placed upon machine performance have resulted in the need for more comprehensive particle accelerator modeling. Computer simulations are key to the success of particle accelerators. Many aspects of particle accelerators rely on computer modeling at some point, sometimes requiring complex simulation tools and massively parallel supercomputing. Examples include the modeling of beams at extreme intensities and densities (toward the quantum degeneracy limit), and with ultra-fine control (down to the level of individual particles). In the future, adaptively tuned models might also be relied upon to provide beam measurements beyond the resolution of existing diagnostics. Much time and effort has been put into creating accelerator software tools, some of which are highly successful. However, there are also shortcomings such as the general inability of existing software to be easily modified to meet changing simulation needs. In this paper possible mitigating strategies are discussed for issues faced by the accelerator community as it endeavors to produce better and more comprehensive modeling tools. This includes lack of coordination between code developers, lack of standards to make codes portable and/or reusable, lack of documentation, among others.
△ Less
Submitted 24 August, 2021;
originally announced August 2021.
-
Extracting Dynamical Frequencies from Invariants of Motion in Finite-Dimensional Nonlinear Integrable Systems
Authors:
Chad E. Mitchell,
Robert D. Ryne,
Kilean Hwang,
Sergei Nagaitsev,
Timofey Zolkin
Abstract:
Integrable dynamical systems play an important role in many areas of science, including accelerator and plasma physics. An integrable dynamical system with $n$ degrees of freedom (DOF) possesses $n$ nontrivial integrals of motion, and can be solved, in principle, by covering the phase space with one or more charts in which the dynamics can be described using action-angle coordinates. To obtain the…
▽ More
Integrable dynamical systems play an important role in many areas of science, including accelerator and plasma physics. An integrable dynamical system with $n$ degrees of freedom (DOF) possesses $n$ nontrivial integrals of motion, and can be solved, in principle, by covering the phase space with one or more charts in which the dynamics can be described using action-angle coordinates. To obtain the frequencies of motion, both the transformation to action-angle coordinates and its inverse must be known in explicit form. However, no general algorithm exists for constructing this transformation explicitly from a set of $n$ known (and generally coupled) integrals of motion. In this paper we describe how one can determine the dynamical frequencies of the motion as functions of these $n$ integrals in the absence of explicitly-known action-angle variables, and we provide several examples.
△ Less
Submitted 10 June, 2021; v1 submitted 4 June, 2021;
originally announced June 2021.
-
A Staged Muon Accelerator Facility For Neutrino and Collider Physics
Authors:
Jean-Pierre Delahaye,
Charles Ankenbrandt,
Stephen Brice,
Alan David Bross,
Dmitri Denisov,
Estia Eichten,
Stephen Holmes,
Ronald Lipton,
David Neuffer,
Mark Alan Palmer,
S. Alex Bogacz,
Patrick Huber,
Daniel M. Kaplan,
Pavel Snopok,
Harold G. Kirk,
Robert B. Palmer,
Robert D. Ryne
Abstract:
Muon-based facilities offer unique potential to provide capabilities at both the Intensity Frontier with Neutrino Factories and the Energy Frontier with Muon Colliders. They rely on a novel technology with challenging parameters, for which the feasibility is currently being evaluated by the Muon Accelerator Program (MAP). A realistic scenario for a complementary series of staged facilities with in…
▽ More
Muon-based facilities offer unique potential to provide capabilities at both the Intensity Frontier with Neutrino Factories and the Energy Frontier with Muon Colliders. They rely on a novel technology with challenging parameters, for which the feasibility is currently being evaluated by the Muon Accelerator Program (MAP). A realistic scenario for a complementary series of staged facilities with increasing complexity and significant physics potential at each stage has been developed. It takes advantage of and leverages the capabilities already planned for Fermilab, especially the strategy for long-term improvement of the accelerator complex being initiated with the Proton Improvement Plan (PIP-II) and the Long Baseline Neutrino Facility (LBNF). Each stage is designed to provide an R&D platform to validate the technologies required for subsequent stages. The rationale and sequence of the staging process and the critical issues to be addressed at each stage, are presented.
△ Less
Submitted 5 February, 2015;
originally announced February 2015.
-
Snowmass 2013 Computing Frontier: Accelerator Science
Authors:
P. Spentzouris,
E. Cormier-Michel,
C. Joshi,
J. Amundson,
W. An,
D. L. Bruhwiler,
J. R. Cary,
B. Cowan,
V. K. Decyk,
E. Esarey,
R. A. Fonseca,
A. Friedman,
C. G. R. Geddes,
D. P. Grote,
I. Kourbanis,
W. P. Leemans,
W. Lu,
W. B. Mori,
C. Ng,
Ji Qiang,
T. Roberts,
R. D. Ryne,
C. B. Schroeder,
L. O. Silva,
F. S. Tsung
, et al. (2 additional authors not shown)
Abstract:
This is the working summary of the Accelerator Science working group of the Computing Frontier of the Snowmass meeting 2013. It summarizes the computing requirements to support accelerator technology in both Energy and Intensity Frontiers.
This is the working summary of the Accelerator Science working group of the Computing Frontier of the Snowmass meeting 2013. It summarizes the computing requirements to support accelerator technology in both Energy and Intensity Frontiers.
△ Less
Submitted 8 October, 2013;
originally announced October 2013.
-
A fast high-order method to calculate wakefield forces in an electron beam
Authors:
Ji Qiang,
Chad Mitchell,
Robert D. Ryne
Abstract:
In this paper we report on a high-order fast method to numerically calculate wakefield forces in an electron beam given a wake function model. This method is based on a Newton-Cotes quadrature rule for integral approximation and an FFT method for discrete summation that results in an $O(Nlog(N))$ computational cost, where $N$ is the number of grid points. Using the Simpson quadrature rule with an…
▽ More
In this paper we report on a high-order fast method to numerically calculate wakefield forces in an electron beam given a wake function model. This method is based on a Newton-Cotes quadrature rule for integral approximation and an FFT method for discrete summation that results in an $O(Nlog(N))$ computational cost, where $N$ is the number of grid points. Using the Simpson quadrature rule with an accuracy of $O(h^4)$, where $h$ is the grid size, we present numerical calculation of the wakefields from a resonator wake function model and from a one-dimensional coherent synchrotron radiation (CSR) wake model. Besides the fast speed and high numerical accuracy, the calculation using the direct line density instead of the first derivative of the line density avoids numerical filtering of the electron density function for computing the CSR wakefield force.
△ Less
Submitted 23 March, 2012;
originally announced March 2012.
-
A Model for One-Dimensional Coherent Synchrotron Radiation including Short-Range Effects
Authors:
Robert D. Ryne,
Bruce Carlsten,
Ji Qiang,
Nikolai Yampolsky
Abstract:
A new model is presented for simulating coherent synchrotron radiation (CSR) in one dimension. The method is based on convolving an integrated Green function (IGF) with the longitudinal charge density. Since it is based on an IGF, the accuracy of this approach is determined by how well one resolves the charge density and not by resolving the single particle wake function. Since short-range wakefie…
▽ More
A new model is presented for simulating coherent synchrotron radiation (CSR) in one dimension. The method is based on convolving an integrated Green function (IGF) with the longitudinal charge density. Since it is based on an IGF, the accuracy of this approach is determined by how well one resolves the charge density and not by resolving the single particle wake function. Since short-range wakefield effects are included analytically, the approach can be much more efficient than ordinary (non-IGF) approaches in situations where the wake function and charge density have disparate spatial scales. Two cases are presented: one derived from the full wake including short-range effects, and one derived from the asymptotic wake. In the latter case the algorithm contains the same physics as others based on the asymptotic approximation, but requires only the line charge density and not its derivative. Examples are presented that illustrate the limitations of the asymptotic-wake approximation, and that illustrate how microbunching can increase the CSR fields by orders of magnitude depending on the microbunching wavelength.
△ Less
Submitted 24 February, 2012; v1 submitted 10 February, 2012;
originally announced February 2012.
-
On FFT-based convolutions and correlations, with application to solving Poisson's equation in an open rectangular pipe
Authors:
Robert D. Ryne
Abstract:
A new method is presented for solving Poisson's equation inside an open-ended rectangular pipe. The method uses Fast Fourier Transforms (FFTs) to perform mixed convolutions and correlations of the charge density with the Green function. Descriptions are provided for algorithms based on the ordinary Green function and for an integrated Green function (IGF). Due to its similarity to the widely used…
▽ More
A new method is presented for solving Poisson's equation inside an open-ended rectangular pipe. The method uses Fast Fourier Transforms (FFTs) to perform mixed convolutions and correlations of the charge density with the Green function. Descriptions are provided for algorithms based on the ordinary Green function and for an integrated Green function (IGF). Due to its similarity to the widely used Hockney algorithm for solving Poisson's equation in free space, this capability can be easily implemented in many existing particle-in-cell beam dynamics codes.
△ Less
Submitted 21 November, 2011;
originally announced November 2011.
-
Exploring Minimal Scenarios to Produce Transversely Bright Electron Beams Using the Eigen-Emittance Concept
Authors:
Leanne D. Duffy,
Kip A. Bishofberger,
Bruce C. Carlsten,
Alex Dragt,
Quinn R. Marksteiner,
Steven J. Russell,
Robert D. Ryne,
Nikolai A. Yampolsky
Abstract:
Next generation hard X-ray free electron lasers require electron beams with low transverse emittance. One proposal to achieve these low emittances is to exploit the eigen-emittance values of the beam. The eigen-emittances are invariant under linear beam transport and equivalent to the emittances in an uncorrelated beam. If a correlated beam with two small eigen-emittances can be produced, removal…
▽ More
Next generation hard X-ray free electron lasers require electron beams with low transverse emittance. One proposal to achieve these low emittances is to exploit the eigen-emittance values of the beam. The eigen-emittances are invariant under linear beam transport and equivalent to the emittances in an uncorrelated beam. If a correlated beam with two small eigen-emittances can be produced, removal of the correlations via appropriate optics will lead to two small emittance values, provided non-linear effects are not too large. We study how such a beam may be produced using minimal linear correlations. We find it is theoretically possible to produce such a beam, however it may be more difficult to realize in practice. We identify linear correlations that may lead to physically realizable emittance schemes and discuss promising future avenues.
△ Less
Submitted 11 July, 2011;
originally announced July 2011.
-
Large-Scale Simulation of Beam Dynamics in High Intensity Ion Linacs Using Parallel Supercomputers
Authors:
Robert D. Ryne,
Ji Qiang
Abstract:
In this paper we present results of using parallel supercomputers to simulate beam dynamics in next-generation high intensity ion linacs. Our approach uses a three-dimensional space charge calculation with six types of boundary conditions. The simulations use a hybrid approach involving transfer maps to treat externally applied fields (including rf cavities) and parallel particle-in-cell techniq…
▽ More
In this paper we present results of using parallel supercomputers to simulate beam dynamics in next-generation high intensity ion linacs. Our approach uses a three-dimensional space charge calculation with six types of boundary conditions. The simulations use a hybrid approach involving transfer maps to treat externally applied fields (including rf cavities) and parallel particle-in-cell techniques to treat the space-charge fields. The large-scale simulation results presented here represent a three order of magnitude improvement in simulation capability, in terms of problem size and speed of execution, compared with typical two-dimensional serial simulations. Specific examples will be presented, including simulation of the spallation neutron source (SNS) linac and the Low Energy Demonstrator Accelerator (LEDA) beam halo experiment.
△ Less
Submitted 30 September, 2000;
originally announced October 2000.
-
Halo Formation in Spheroidal Bunches with Self-Consistent Stationary Distributions
Authors:
A. V. Fedotov,
R. L. Gluckstern,
S. S. Kurennoy,
R. D. Ryne
Abstract:
A new class of self-consistent 6-D phase space stationary distributions is constructed both analytically and numerically. The beam is then mismatched longitudinally and/or transversely, and we explore the beam stability and halo formation for the case of 3-D axisymmetric beam bunches using particle-in-cell simulations. We concentrate on beams with bunch length-to-width ratios varying from 1 to 5…
▽ More
A new class of self-consistent 6-D phase space stationary distributions is constructed both analytically and numerically. The beam is then mismatched longitudinally and/or transversely, and we explore the beam stability and halo formation for the case of 3-D axisymmetric beam bunches using particle-in-cell simulations. We concentrate on beams with bunch length-to-width ratios varying from 1 to 5, which covers the typical range of the APT linac parameters. We find that the longitudinal halo forms first for comparable longitudinal and transverse mismatches. An interesting coupling phenomenon - a longitudinal or transverse halo is observed even for very small mismatches if the mismatch in the other plane is large - is discovered.
△ Less
Submitted 20 July, 1998;
originally announced July 1998.
-
Lyapunov Exponents without Rescaling and Reorthogonalization
Authors:
Govindan Rangarajan,
Salman Habib,
Robert D. Ryne
Abstract:
We present a new method for the computation of Lyapunov exponents utilizing representations of orthogonal matrices applied to decompositions of M or MM_trans where M is the tangent map. This method uses a minimal set of variables, does not require renormalization or reorthogonalization, can be used to efficiently compute partial Lyapunov spectra, and does not break down when the Lyapunov spectru…
▽ More
We present a new method for the computation of Lyapunov exponents utilizing representations of orthogonal matrices applied to decompositions of M or MM_trans where M is the tangent map. This method uses a minimal set of variables, does not require renormalization or reorthogonalization, can be used to efficiently compute partial Lyapunov spectra, and does not break down when the Lyapunov spectrum is degenerate.
△ Less
Submitted 12 March, 1998;
originally announced March 1998.
-
Chaos in Time Dependent Variational Approximations to Quantum Dynamics
Authors:
Fred Cooper,
John Dawson,
Salman Habib,
Robert D. Ryne
Abstract:
Dynamical chaos has recently been shown to exist in the Gaussian approximation in quantum mechanics and in the self-consistent mean field approach to studying the dynamics of quantum fields. In this study, we first show that any variational approximation to the dynamics of a quantum system based on the Dirac action principle leads to a classical Hamiltonian dynamics for the variational parameter…
▽ More
Dynamical chaos has recently been shown to exist in the Gaussian approximation in quantum mechanics and in the self-consistent mean field approach to studying the dynamics of quantum fields. In this study, we first show that any variational approximation to the dynamics of a quantum system based on the Dirac action principle leads to a classical Hamiltonian dynamics for the variational parameters. Since this Hamiltonian is generically nonlinear and nonintegrable, the dynamics thus generated can be chaotic, in distinction to the exact quantum evolution. We then restrict attention to a system of two biquadratically coupled quantum oscillators and study two variational schemes, the leading order large N (four canonical variables) and Hartree (six canonical variables) approximations. The chaos seen in the approximate dynamics is an artifact of the approximations: this is demonstrated by the fact that its onset occurs on the same characteristic time scale as the breakdown of the approximations when compared to numerical solutions of the time-dependent Schrodinger equation.
△ Less
Submitted 16 October, 1996; v1 submitted 9 October, 1996;
originally announced October 1996.
-
Finding matched rms envelopes in rf linacs: A Hamiltonian approach
Authors:
Robert D. Ryne
Abstract:
We present a new method for obtaining matched solutions of the rms envelope equations. In this approach, the envelope equations are first expressed in Hamiltonian form. The Hamiltonian defines a nonlinear mapping, $\cal M$, and for periodic transport systems the fixed points of the one-period map are the matched envelopes. Expanding the Hamiltonian around a fiducial trajectory one obtains a linear…
▽ More
We present a new method for obtaining matched solutions of the rms envelope equations. In this approach, the envelope equations are first expressed in Hamiltonian form. The Hamiltonian defines a nonlinear mapping, $\cal M$, and for periodic transport systems the fixed points of the one-period map are the matched envelopes. Expanding the Hamiltonian around a fiducial trajectory one obtains a linear map, $M$, that describes trajectories (rms envelopes) near the fiducial trajectory. Using $\cal M$ and $M$ we construct a contraction mapping that can be used to obtain the matched envelopes. The algorithm is quadratically convergent. Using the zero-current matched parameters as starting values, the contraction mapping typically converges in a few to several iterations. Since our approach uses numerical integration to obtain all the mappings, it includes the effects of nonidealized, $z$-dependent transverse and longitudinal focusing fields. We present numerical examples including finding a matched beam in a quadrupole channel with rf bunchers.
△ Less
Submitted 11 February, 1995;
originally announced February 1995.
-
Symplectic Computation of Lyapunov Exponents
Authors:
Salman Habib,
Robert D. Ryne
Abstract:
A recently developed method for the calculation of Lyapunov exponents of dynamical systems is described. The method is applicable whenever the linearized dynamics is Hamiltonian. By utilizing the exponential representation of symplectic matrices, this approach avoids the renormalization and reorthogonalization procedures necessary in usual techniques. It is also easily extendible to damped syste…
▽ More
A recently developed method for the calculation of Lyapunov exponents of dynamical systems is described. The method is applicable whenever the linearized dynamics is Hamiltonian. By utilizing the exponential representation of symplectic matrices, this approach avoids the renormalization and reorthogonalization procedures necessary in usual techniques. It is also easily extendible to damped systems. The method is illustrated by considering two examples of physical interest: a model system that describes the beam halo in charged particle beams and the driven van der Pol oscillator.
△ Less
Submitted 1 November, 1994;
originally announced November 1994.
-
Symplectic Calculation of Lyapunov Exponents
Authors:
Salman Habib,
Robert D. Ryne
Abstract:
The Lyapunov exponents of a chaotic system quantify the exponential divergence of initially nearby trajectories. For Hamiltonian systems the exponents are related to the eigenvalues of a symplectic matrix. We make use of this fact to develop a new method for the calculation of Lyapunov exponents of such systems. Our approach avoids the renormalization and reorthogonalization of usual techniques.…
▽ More
The Lyapunov exponents of a chaotic system quantify the exponential divergence of initially nearby trajectories. For Hamiltonian systems the exponents are related to the eigenvalues of a symplectic matrix. We make use of this fact to develop a new method for the calculation of Lyapunov exponents of such systems. Our approach avoids the renormalization and reorthogonalization of usual techniques. It is also easily extendible to damped systems. We apply our method to two examples of physical interest: a model system that describes the beam halo in charged particle beams and the driven van der Pol oscillator.
△ Less
Submitted 25 June, 1994;
originally announced June 1994.