An alternative representation for the vector potential of the nonlinear magnetic insert for the Integrable Optics Test Accelerator (IOTA), first described in Sec. V.A. of \cite{Danilov}, is determined from first principles using standard complex variable methods. In particular, it is shown that the coupled system consisting of the 2D Laplace equation and the Bertrand-Darboux equation is equivalent to a single ordinary differential equation in the complex plane, and a simple solution is constructed. The results are consistent with \cite{Danilov}, and this concise representation provides computational advantages for particle tracking through the nonlinear insert by avoiding numerical errors caused by small denominators that appear when evaluating transverse derivatives of the vector potential near the midplane. A similar representation is provided for the spatial dependence of the two invariants of motion.

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.

A numerical method is described for producing stationary solutions of the Vlasov-Poisson system describing a relativistic charged-particle beam in a constant focusing accelerator channel, confined transversely by a general (linear or nonlinear) focusing potential. The method utilizes a variant of the spectral Galerkin algorithm to solve a nonlinear partial differential equation (PDE) in two degrees of freedom for the beam space charge potential in equilibrium. Numerical convergence with an increasing number of computed spectral modes is investigated for several benchmark problems. Preservation of the stationary phase space density is verified using a strongly nonlinear focusing channel based on the Integrable Optics Test Accelerator at Fermi National Accelerator Laboratory.

The Integrable Optics Test Accelerator (IOTA) is a novel storage ring at Fermi National Accelerator Laboratory designed (in part) to investigate the dynamics of beams in the presence of highly nonlinear transverse focusing fields that generate integrable single-particle motion with a large spread in the intrinsic betatron tunes. We describe how contemporary geometrical methods from the theory of integrable Hamiltonian systems may be used to locate all critical separatrixlike structures in the 4D transverse phase space, and to construct a complete analysis of the dynamical bifurcations of the system. Application of these techniques results in a global picture of the nominal on-energy transverse dynamics, revealing a rich diversity of accessible dynamical behavior. Similar techniques may be applied to future facilities that exploit the concept of nonlinear integrable optics.

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 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.

Laser-driven ion acceleration provides ultrashort, high-charge, low-emittance beams, which are desirable for a wide range of high-impact applications. Yet after decades of research, a significant increase in maximum ion energy is still needed. This paper introduces a quality-preserving staging concept for ultraintense ion bunches that is seamlessly applicable from the nonrelativistic plasma source to the relativistic regime. Full three-dimensional particle-in-cell simulations prove robustness and capture of a high-charge proton bunch, suitable for readily available and near-term laser facilities.