Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Beam dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Beam results for the new super-conducting booster at Jefferson Lab

The performance results of a new super-conducting booster for the CEBAF injector at Jefferson Lab, could be of interest for other similar electron injectors. A recent addition of this new booster has provided us the ability to achieve a more adiabatic acceleration and therefore an improvement to the beam dynamics and beam brightness. It has also simplified the design and operation of the section of the injector responsible for accelerating the electron beam from a few hundred keV to several MeV (typically 6.7 MeV). The addition of the new booster was part of an upgrade to the CEBAF injector to improve the beam quality for future physics experiments with high sensitivity to beam quality. The booster consists of two cavities: a 2-cell cavity followed by a 7-cell cavity. This combination allows for a wide range of input electron beam energies, from 130 keV to more than 300 keV. In fact, during the last year, the booster was successfully operated with 140,180, and 200 keV input beam energies as the electron gun was being upgraded. This paper describes the new booster, presents beam optics data results from different beam studies, commissioning, and the physics quality beam operation.

Accelerator Physics↗
DOE OSTI2026

Experimental Studies of Nonlinear Integrable Optics with Elliptic Potentials

The stable transport of charged particle beams is the core challenge in designing and operating accelerators. The current paradigm for transverse focusing is based on quadrupole and dipole elements, described as a linear integrable Hamiltonian by Courant and Snyder. The operational limits of high intensity accelerators are often determined by collective instabilities within the beam. Nonlinear elements may be added to suppress certain forms of these collective effects, but restrict the range of stable trajectories. These shortcomings motivate extending to a nonlinear integrable system which can offer the benefits of suppressing collective instabilities without limiting the stable trajectories. Danilov and Nagaitsev proposed a novel nonlinear integrable system with elliptical potentials, which could be implemented as magnetic elements in an accelerator. Such a system has been implemented for practical verification in the integrable optics test accelerator (IOTA), a small storage ring constructed for beam dynamics studies. Electron beam studies in IOTA treat the low-emittance beam as a macroparticle for detailed probing of the expected single particle dynamics. This dissertation details new measurements of the nonlinear integrable system relevant for practical implementation. Turn-by-turn measurements of the kicked beam responses are used for phase space reconstruction and analysis of the predicted nonlinear dynamics. The stability and aperture for various nonlinear configurations are measured with beam losses. Amplitude dependent detuning, a core figure of merit for suppressing instabilities, is measured and compared with high fidelity particle tracking simulations. The synchrotron radiation images of circulating beam allow direct measurements of the topology and lifetime in configurations where nonlinear focusing terms dominate.

Wieland, John [Michigan U.] (ORCID:000000028971852↗
DOE OSTI2025

Second harmonic Bessel-Gauss beam shaping with elliptic axicon aberrations

The second harmonic (SH) of an axicon generated Bessel-Gauss beam is created through the nonlinear interaction of photons with a crystal, resulting in the energy doubling of the output photons. In this work, we show experimentally that in addition to frequency doubling, the SH of Bessel-Gauss beams under asymmetric aberrations from an elliptic axicon exhibit intriguing beam formation. Particularly, the central region of the SH beam profile is composed of two central spots of various geometries surrounded by nested ellipses; one of which is the configuration of two central gamma dots with similar radius knotted by nested ellipses for a zeroth-order Bessel-Gauss pump. These SH beams consistently maintain their spatial profile throughout propagation, reminiscing the behavior of screw dislocations in wave patterns. Our numerical simulations produce beam dynamics consistent with that of experiments and further implicate the remarkable interweaving of bright spots with dark vortices. This is especially noticeable when the beams dynamically oscillate along the optical axis, resulting in the genesis of spatially polarized beams with a knotted framework. The insights gained from our study establish a novel paradigm for exploring interactions of Bessel-like beams with vortex dynamics. This, in turn, has the potential to spark innovations in optical applications, fostering new methodologies to harness and manipulate complex light structures. Our experimental findings could spur new methods of generating logical states of light and new opportunities for material processing control. Published by the American Physical Society 2025

Wang, Tianhong (ORCID:0000000169104201)↗
DOE OSTI2024

A Novel Two Stage Collimation Unit for Fermilab Booster

A new two-stage collimation unit (2SC) for Fermilab Booster will be installed during 2024 summer shutdown. It is a supplementary collimator for existing single stage Booster collimators. Design details of this 2SC adapted to Booster conditions are described. Results of beam dynamics simulations on collimation efficiency and evaluation of collimator shielding with MARS code on this new system are presented. The analysis on beam dynamics shows that the system will help to collimate up to about 1% of the beam in a confined way and, prompt and residual activation meet Fermilab Radiological Control limits in current as well as during PIP-II era high intensity operation of the Booster.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Extending JuTrack’s capabilities to the FRIB accelerator to enhance online modeling

JuTrack is a Julia-based accelerator modeling and tracking package that utilizes compiler-level automatic differentiation (AD) to enable fast and accurate derivative calculations. While JuTrack provides a solid foundation for beam dynamics simulations, its capabilities must be extended to support the Facility for Rare Isotopes (FRIB) linac. This includes modeling heavy-ion linac accelerator components such as the liquid-lithium charge stripper, which facilitates efficient acceleration by remove electrons from heavy isotopes, and incorporating multi-charge state acceleration tracking, which allows for charge-dependent beam dynamics. These extensions address challenges such as the beam matching and optimization of multi charge state through various accelerating structures and beam-material interaction modeling while maintaining the auto differentiation capability. This work focuses on adapting JuTrack to incorporate these elements, enhancing its online modeling abilities. We present modifications to JuTrack’s framework and demonstrate their performance in FRIB simulations.

Accelerator Physics↗
DOE OSTI2025

On Final Results From the Muon $g\mathrm{-}2$ Experiment at Fermilab

The Muon g-2 Experiment at Fermilab has measured the muon anomalous magnetic moment, a_mu, with unprecedented precision, leveraging a dataset from Runs 1 6 that is 21 times larger than its Brookhaven predecessor. This talk will present the experiment s final result, which serves as a benchmark for testing the Standard Model with high precision. The measurement was performed in a storage ring with a highly uniform magnetic field, where precise beam dynamics understanding is vital to determine the muon anomalous precession frequency. Key to this effort were advanced simulation tools including COSY INFINITY, precise fringe field modeling, and calculations of beam dynamics characteristics like tunes and chromaticity. We will conclude by exploring how the experiment's findings reshape our understanding of particle physics.

Valetov, Eremey [Michigan State U.]↗
DOE OSTI2025

Leveraging prior mean models for faster Bayesian optimization of particle accelerators

Tuning particle accelerators is a challenging and time-consuming task that can be automated and carried out efficiently using suitable optimization algorithms, such as model-based Bayesian optimization techniques. One of the major advantages of Bayesian algorithms is the ability to incorporate prior information about beam physics and historical behavior into the model used to make control decisions. In this work, we examine incorporating prior accelerator physics information into Bayesian optimization algorithms by utilizing fast executing, neural network models trained on simulated or historical datasets as prior mean functions in Gaussian process models. We show that in ideal cases, this technique substantially increases convergence speed to optimal solutions in high-dimensional tuning parameter spaces. Additionally, we demonstrate that even in non-ideal cases, where prior models of beam dynamics do not exactly match experimental conditions, the use of this technique can still enhance convergence speed. Finally, we demonstrate how these methods can be used to improve optimization in practical applications, such as transferring information gained from beam dynamics simulations to online control of the LCLS injector, and transferring knowledge gained from experimental measurements across different operating modes, such as accelerating different ion species at the ATLAS heavy ion accelerator.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Study of Quadrupole Fringe Fields in the Interaction Region of the Hadron Storage Ring of the Electron Ion Collider

Fringe fields in quadrupole magnets are usually neglected in studies of beam dynamics at accelerators. However, the extreme optical parameters present in the final focus of a collider such as the Electron–Ion Collider (EIC) may give rise to effects that should not be overlooked. The calculation of quadrupole fringe fields presented in this study follows the procedure outlined in Ref. [1], specialized to the case of a straight reference orbit (i.e., with no dipole field component). A right-handed Cartesian coordinate system is employed, with the $z$-axis aligned with the quadrupole axis and $x$ and $y$ denoting the horizontal and vertical transverse coordinates, respectively. The magnetic quadrupole field gradient, $G = \partial B_y / \partial x$, transitions from its peak value inside the quadrupole—where it is nearly independent of the longitudinal coordinate $z$—to zero at some distance beyond the magnet edge. Consequently, $G$ is treated as a function of $z$. The region over which this variation occurs is defined as the quadrupole fringe field region. The study begins with the development of a description of the magnetic field in the fringe region using a power-series expansion in the transverse coordinates $x$and $y$, consistent with the longitudinally varying gradient. A model for the $z$-dependence of the gradient is then proposed and adjusted to reproduce magnetic field data obtained from three-dimensional field calculations. To evaluate the impact of the fringe fields on beam dynamics, the corresponding vector potential is derived and incorporated into a Hamiltonian formulation of particle motion. The significance of the fringe fields is quantified by calculating the amplitude-dependent tune shift from the Hamiltonian. Using linear beam optics parameters of the Hadron Storage Ring (HSR) of the EIC, the tune shift due to the fringe fields of all quadrupole magnets in the IR-6 interaction region is evaluated. Finally, the resulting tune shifts are compared with those arising from other nonlinear field components present in the HSR.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Digital twin framework for PIP-II linac: AI-driven multi-scale modeling from ion source to 800 MeV

The PIP-II superconducting linac at Fermilab is designed to deliver multi-megawatt proton beams for neutrino physics and other high-intensity applications. To expedite commissioning and enhance operational reliability, we have developed an EPICS-based data flow framework that seamlessly integrates digital twins (DT) with physical twins (PT). These digital twins comprise high-fidelity beam dynamics models or data-driven surrogate models connected to their physical counterparts through real-time diagnostics and advanced machine-learning algorithms.Central to this framework is Linac_Gen, an accelerated simulation tool that incorporates convolutional neural networks, random forests, and genetic algorithms to provide up to a tenfold speedup in optimizing the accelerator geometry model. An EPICS translator layer ensures interoperability by efficiently mapping lattice parameters across diverse simulation platforms.Our EPICS-based framework supports multiple operational modes—monitoring, passive learning, closed-loop control, and online learning—covering the entire machine lifecycle. By leveraging HPC resources and multi-objective optimization techniques, the digital twin enables adaptive trajectory correction, real-time fault detection, and predictive modeling of beam stability. This comprehensive approach paves the way for robust, high-intensity operation and data-driven accelerator R&D at Fermilab.

Pathak, Abhishek [Fermilab]↗
DOE OSTI2026

Analytic Chromaticity Formulas for the Muon g -2 experiment at fermilab

The Muon g-2 Experiment (E989) at Fermilab measured the muon anomalous magnetic moment aμ with unprecedented precision of 127 ppb using a storage ring. The experiment pushed the limits in terms of accuracy and systematic errors. Achieving the 78-ppb total systematic uncertainty required precise beam dynamics modeling and corrections for effects on the measured muon precession frequency. Here, we derived the first analytic aberration formulas up to the second order for the muon g-2 storage ring’s combined-function electrostatic quadrupoles (superimposed magnetic dipole and electric quadrupole fields) using an order-by-order perturbation method. From these, we obtained the exact chromaticity formulas for three ring models of different granularity and validated them against numerical calculations using COSY INFINITY, achieving analytic-numerical agreement to ⁠$\mathscr{O}(10^{-10})$. This work resolved discrepancies between previous approximate derivations and provided essential beam dynamics results for Runs 4-6 analyses. We also calculated nonlinear chromaticities up to ninth order. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous muon g-2 experiment at Brookhaven National Laboratory, with the final result announced in June 2025.

aberrations↗
DOE OSTI2026

Status of 22 GeV CEBAF Upgrade Using FFA Arcs

Jefferson Lab is exploring an upgrade to extend CEBAF’s energy to ~22 GeV within the existing tunnel by incorporating Fixed-Field Alternating-Gradient (FFA) arcs. In this scheme, the two lowest-energy electromagnetic arcs are removed, the remaining arcs are reassigned to lower energies, and a new pair of high-energy FFA arcs enables six additional recirculations using the existing SRF linacs. The non-scaling FFA lattice employs Halbach-derived permanent magnets with provisions for dipole, quadrupole, and higher-order components, offering large momentum acceptance and reduced operating cost. Implementing this concept requires updated linac optics, modifications to the remaining arcs, and a redesigned electromagnetic switchyard. Ongoing beam-dynamics studies address synchrotron-radiation–induced energy loss and emittance growth, along with strategies for their mitigation.

Ogur, S. [Thomas Jefferson National Accelerator Fa↗
DOE OSTI2026

Dynamic Multipoles in the RCS

The Rapid Cycling Synchrotron (RCS) of the Electron--Ion Collider (EIC) will accelerate polarized electrons up to 18 GeV at a repetition rate of 1 Hz. Dynamic multipoles driven by eddy currents in the vacuum chamber can affect beam dynamics and polarization preservation during the energy ramp. Building on the analytical treatment of S.Y. Lee for a ramped normal-conducting dipole magnet, we use conformal transformations to derive compact expressions for dynamic multipoles in higher-order ramped magnets with circular beam pipes. We further derive analytical frequency-domain transfer functions for both self-response and cross-response multipoles, providing compact dynamical models for their evolution during arbitrary ramp waveforms. The resulting analytical predictions are benchmarked against finite-element simulations, including the effects of practical pole geometry and beam-pipe misalignment. Finally, the expected dynamic multipoles in the RCS magnets are estimated and their implications for beam dynamics and polarization preservation are discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

JuTrack: A Julia package for auto-differentiable accelerator modeling and particle tracking

Efficient accelerator modeling and particle tracking are key for the design and configuration of modern particle accelerators. In this work, we present JuTrack, a nested accelerator modeling package developed in the Julia programming language and enhanced with compiler-level automatic differentiation (AD). With the aid of AD, JuTrack enables rapid derivative calculations in accelerator modeling, facilitating sensitivity analyses and optimization tasks. Here we demonstrate the effectiveness of AD-derived derivatives through several practical applications, including sensitivity analysis of space-charge-induced emittance growth, nonlinear beam dynamics analysis for a synchrotron light source, and lattice parameter tuning of the future Electron-Ion Collider (EIC). Through the incorporation of automatic differentiation, this package opens up new possibilities for accelerator physicists in beam physics studies and accelerator design optimization.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Preliminary Look at the LBEG & MPEG Beam Transmissions between HPSim and Operation Data

This report summarizes recent work on estimating the transmission of LANSCE H- beams from the end of the present DTL to both the PSR stripper foil (LBEG) and WNR target 4 (MPEG). The LBEG beam might be considered a more typical LINAC beam and lends itself to continuous monitoring of transmission through the various stages of the accelerator. For the MPEG, however, the widely-spaced micropulses and lack of the requisite sensitivity current monitors throughout the accelerator make these measurements extremely difficult and more uncertain. Therefore, to assist in estimating the MPEG beam transmission, beam-dynamics simulations using HPSim were employed. These beam transmission estimates from the end of the Drift Tube Linac (DTL) to Target 4 (MPEG) and the PSR stripper foil (LBEG) are vital to determine the charge requirements for LAMP’s front end. In this technote, we demonstrate three major efforts in determining the transmission: (1) Convert the WNR beamline lattice from TRANSPORT to HPSim for use in the simulation; (2) Simulate the optimization process in the Central Control Room (CCR) that brings down the losses by up to a factor of 4000 between the end of the initial physics-based tuneup phase and production beam operation; (3) Analyze operational data to deduce measured transmissions. Table 1 shows the estimated losses with HPSim and operational analysis. Finally, a better measurement and other improvements to refine the results are also proposed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Simulation and measurement of beam-induced heating of ceramic vacuum chambers

In this article, we summarize recent theoretical and experimental studies of the impedance and beam-induced heating of titanium-coated ceramic vacuum chambers used in the NSLS-II injection kickers. The impedance was calculated using the field matching theory assuming planar approximation and compared with the mpedanceake2 code. For the coating thickness of a few microns, we demonstrated that the beam-induced power is dissipated in the titanium coating and that the longitudinally averaged two-dimensional power density is approximated by an analytical expression, thus allowing the use of a simplified model of the power density as input for the code to simulate the temperature distribution with realistic nonuniform thickness of the Ti coating. For a few values of the NSLS-II beam current, we measured the beam-induced heating of two ceramic chambers using thermal sensors installed along the chamber and compared the measurement results with the simulations. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2024

3D small-gain formula allowing strong focusing and harmonic lasing for a ring-based x-ray free electron laser oscillator

We present a detailed derivation of a formula for the small-gain calculation for an x-ray free electron laser oscillator (XFELO) based on a medium-energy (3–4 GeV) storage ring. We found harmonic lasing and strong focusing are essential for this beam energy range. Taking the small-signal low-gain formula developed by Kim and his colleagues, we modified it in such a way that the gain can be calculated without the “no focusing approximation,” and a strong focusing can be applied, as well as harmonic lasing. In this formula, the gain is represented as a product of two factors with one of them depending only on the harmonic number, undulator period, and gap. Using this factor, we show that it is favorable to use harmonic lasing to achieve hard x-ray FEL working in the small-signal low-gain regime with the medium-energy electron beam. Our formula also allows FEL optimization by varying the vertical gradient of the undulator, the vertical dispersion, and the horizontal and vertical focusing, independently. As an example, we applied this formula to study the feasibility of an XFELO option for the National Synchrotron Light Source II (NSLS-II) upgrade. Since a quite high peak current is required for the FEL, collective effects of beam dynamics in medium-energy synchrotrons significantly affect the electron beam parameters. We carried out a multiparameter optimization taking collective effects into account. Note, even though our example is for a ring-based XFELO at 3 to 4 GeV, the formula and, in particular, the approach developed here may be applied to other types of FELs. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2024

Using a particle-in-cell model for accelerator control room applications

Many accelerator control rooms rely on envelope models to simulate beam dynamics because they are fast and accurate at tracking the beam core. Particle-in-Cell (PIC) models, however, can track particles inside and outside the core and, with the improvements of computers, are now fast enough to be used in control rooms. The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory is currently developing a tool to use a Particle-in-Cell model for control room applications. This report covers the progress so far and the future goals of using PyORBIT, a Particle-in-Cell simulation model, in the SNS control room.

Cathey, Brandon↗
DOE OSTI2024

Capture Cavities For The Cw Polarized Positron Source Ce + baf

The initial design of the capture cavities for the continuous wave (CW) polarized positron beams at Jefferson Lab (Ce + BAF) is presented. A chain of standing wave multi-cell copper cavities inside a solenoid tunnel are selected to improve the positron capture efficiency. The cavity design strategy is presented to accommodate constrains from the large phase distribution of the incident beams, RF power and RF heating. A matrix of design parameters’ range are given for future system optimization when the capture cavities are considered together with other sub-systems and beam dynamics. The contents will also be useful for other CW cavity design for beams with large phase space distribution.

43 PARTICLE ACCELERATORS↗