Skip to main content
eScholarship
Open Access Publications from the University of California
Cover page of JuTrack: A Julia package for auto-differentiable accelerator modeling and particle tracking

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

(2025)

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. 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. Program summary: Program Title: JuTrack CPC Library link to program files: https://doi.org/10.17632/r2g5zkwp7s.1 Developer's repository link: https://github.com/MSU-Beam-Dynamics/JuTrack.jl.git Licensing provisions: MIT Programming language: Julia Nature of problem: Derivatives of the physics parameters calculated in accelerator modeling are critical for sensitivity analysis and optimization of the whole system. Traditional numerical approaches often rely on finite differences for derivative computations, which can prone to numerical inaccuracies. In highly nonlinear accelerator systems, like those encountered in synchrotrons and colliders, accurate sensitivity analysis and optimization require a large number of derivative evaluations. Thus, there is a need for more efficient methods to compute these derivatives accurately, especially when optimizing complex accelerator lattices or studying complicated collective effects, such as space-charge effects, wakefield effects, and beam-beam interaction. Solution method: JuTrack addresses this problem by integrating compiler-level automatic differentiation (AD) into accelerator modeling routines, offering a powerful toolset for rapid derivative computation. Developed in the Julia programming language, JuTrack uses the Enzyme AD package to perform gradient-based analyses with minimal computational overhead. The package provides an efficient way to compute derivatives by directly differentiating through the model code, thus avoiding approximation errors associated with finite difference methods. It is designed to handle complex beam dynamics simulations, including complicated collective effects, such as space-charge effects, wakefield effects, beam-beam interaction, and combination of Truncated Power Series Algebra (TPSA) with AD. It can be applied to lattice optimization and beam dynamics analysis for future accelerators like the Electron-Ion Collider (EIC). Users can easily apply the package to their models, enabling robust optimization and sensitivity analysis in their accelerator studies. Additional comments including restrictions and unusual features: JuTrack is particularly well-suited for scenarios requiring frequent derivative calculations, such as during beam dynamics optimization, sensitivity analysis, and accelerator tuning. Its integration with the Julia programming language provides excellent performance due to Julia's just-in-time (JIT) compilation capabilities. The modular nature of JuTrack and Julia's easy-to-understand syntax allows for future extensions and custom modifications, making it adaptable to a variety of accelerator configurations.

Cover page of Postannealing-induced intermetallic phase formation in NiPt thin films deposited via direct current and high-power impulse magnetron sputtering

Postannealing-induced intermetallic phase formation in NiPt thin films deposited via direct current and high-power impulse magnetron sputtering

(2025)

Intermetallic phases are preferred to reduce the amount of platinum used for catalytic applications as compared to solid solution alloys, due to their stability at elevated temperatures while preserving or even enhancing the catalytic properties. Here, we show a two-step process to form an intermetallic NiPt L10 phase. In this work, NiPt solid solution thin films were fabricated by direct current and high-power impulse magnetron sputtering processes, which allow for precise thickness and chemical composition control. Following deposition, an additional annealing step is used to form the desired intermetallic phase. We show that the required annealing time for intermetallic phase formation is considerably reduced for NiPt thin films with a thickness of 240 nm, as compared to its bulk counterpart.

Cover page of Residual resistance ratio measurement system for Nb3Sn wires extracted from Rutherford cables

Residual resistance ratio measurement system for Nb3Sn wires extracted from Rutherford cables

(2025)

Residual resistance ratio (RRR) of superconducting strands is an important parameter for magnet electrical stability. RRR serves as a measure of the low-temperature electrical conductivity of the copper within a conductor that has a copper stabilization matrix. For Nb3Sn, due to the need of a reaction heat treatment, the technical requirements for high quality measurements of strands extracted from Rutherford cables are particularly demanding. Quality of wire, cabling deformation, heat treatment temperature, heat treatment atmosphere, sample handling, and measurement methods can all affect the RRR. Therefore, as an integral part of the electrical quality control (QC) of Nb3Sn Rutherford cables manufactured at the Lawrence Berkeley National Laboratory, it was prudent that we established a RRR measurement system that can isolate the assessment of cable-fabrication-related impacts from sample preparation and measurement factors. Here we describe a bespoke cryocooler-based measurement system, capable of measuring RRR of over 80 samples in a single cooldown. The samples are mounted on custom-designed printed circuit boards that accommodate the shape of strands extracted from a Rutherford cable without added deformation, which we will show is critical in ensuring that the measurements accurately represent the RRR values of the conductor within the cable. Using this sample mounting solution, we routinely measure the overall RRR of the strand as well as individual intra-strand sections corresponding to both cable edges and cable broad faces with high reproducibility. Such measurements provide valuable information on the variation of RRR along the length of the strands as well as across strand productions and cable runs over time.

Cover page of Enhanced Isomer Population via Direct Irradiation of Solid-Density Targets Using a Compact Laser-Plasma Accelerator

Enhanced Isomer Population via Direct Irradiation of Solid-Density Targets Using a Compact Laser-Plasma Accelerator

(2025)

Excitation of long-lived states in bromine nuclei using a tabletop laser-plasma accelerator providing pulsed (<100  fs) electron beams provided a sensitive probe of γ strength and level densities in the nuclear quasicontinuum and may indicate angular momentum coupling through electron-nuclear interactions. Solid-density active LaBr_{3} targets absorb real and virtual photons up to 35±2.5  MeV and deexcite through γ cascade into different states. A factor of 4.354±0.932 enhancement of the ^{80}Br^{m}/^{80}Br^{g} isomeric ratio was observed following electron irradiation, as compared to bremsstrahlung. Additional angular momentum transfer could possibly occur through nuclear-plasma or electron-nuclear interactions enabled by the ultrashort electron beam. Further investigation of these mechanisms could have far-reaching impact including decreased storage of long-term nuclear waste and an improved understanding of heavy element formation in astrophysical settings.

Cover page of Compact in-vacuum gamma-ray spectrometer for high-repetition rate PW-class laser–matter interaction

Compact in-vacuum gamma-ray spectrometer for high-repetition rate PW-class laser–matter interaction

(2025)

With the advent of high repetition rate laser facilities, novel diagnostic tools compatible with these advanced specifications are required. This paper presents the design of an active gamma-ray spectrometer intended for these high repetition rate experiments, with particular emphasis on functionality within a PW level laser-plasma interaction chamber's extreme conditions. The spectrometer uses stacked scintillators to accommodate a broad range of gamma-ray energies, demonstrating its adaptability for various experimental setups. In addition, it has been engineered to maintain compactness, electromagnetic pulse resistance, and ISO-5 cleanliness requirements while ensuring high sensitivity. The spectrometer has been tested in real conditions inside the PW-class level interaction chamber at the BELLA center, LBNL. The paper further details the calibration process, which utilizes a 60Co radioactive source, and describes the unfolding technique implemented through a stochastic minimization method.

Decoupling the effects of potential energy, kinetic energy, and ion flux on crystallinity of V-Al and V-Al-N thin films in pulsed filtered cathodic arc deposition

(2025)

Deposition using filtered pulsed cathodic arc plasma is known to produce dense and adherent thin films due to energetic ions which carry significant ion kinetic and potential energies. The role of potential energy coming from multiply charged metal ions in film formation remains under-explored, since the enhancement of charge states in cathodic arcs is coupled with an increased flux and kinetic energies of ions. In this work, the influence of ion potential energy on structural properties of thin films is investigated, while keeping the mean ion kinetic energies unchanged. Two material systems are considered: metallic V-Al and compound V-Al-N in non-reactive and reactive deposition, respectively. For V-Al plasma and thin films, the impact of metal ion potential energy is demonstrated. In the V-Al-N case, in addition to metal ions, activated (namely, ionized, dissociated, and excited) nitrogen species are found to be a significant factor for crystalline growth of the metastable cubic phase.

Cover page of Quasiprobabilistic Readout Correction of Midcircuit Measurements for Adaptive Feedback via Measurement Randomized Compiling

Quasiprobabilistic Readout Correction of Midcircuit Measurements for Adaptive Feedback via Measurement Randomized Compiling

(2025)

Quantum measurements are a fundamental component of quantum computing. However, on present-day quantum computers, measurements can be more error prone than quantum gates and are susceptible to nonunital errors as well as nonlocal correlations due to measurement crosstalk. While readout errors can be mitigated in postprocessing, this is inefficient in the number of qubits due to a combinatorially large number of possible states that need to be characterized. In this work, we show that measurement errors can be tailored into a simple stochastic error model using randomized compiling, enabling the efficient mitigation of readout errors via quasiprobability distributions reconstructed from the measurement of a single preparation state in an exponentially large confusion matrix. We demonstrate the scalability and power of this approach by correcting readout errors without matrix inversion on a large number of different preparation states applied to a register of eight superconducting transmon qubits. Moreover, we show that this method can be extended to midcircuit measurements used for active feedback via quasiprobabilistic error cancellation, and we demonstrate the correction of measurement errors on an ancilla qubit used to detect and actively correct bit-flip errors on an entangled memory qubit. Our approach enables the correction of readout errors on large numbers of qubits and offers a strategy for correcting readout errors in adaptive circuits in which the results of midcircuit measurements are used to perform conditional operations on nonlocal qubits in real time.

Cover page of Voltage-Current Behavior of a Superconducting star Wire in a 6-Around-1 Cable Configuration

Voltage-Current Behavior of a Superconducting star Wire in a 6-Around-1 Cable Configuration

(2025)

A 6-around-1 transposed cable using superconducting star® wires can be useful for future circular collider applications. We made three cable samples using single star® wires with a diameter of 1.3 mm. The first two samples, made with a cabling machine, used star® wires consisting of a 0.7 mm diameter Nb-Ti core. The third cable sample was manually wound and used a star® wire made with a 0.81 mm diameter Cu core. The first sample showed severe degradation after the cable was bent to a 75 mm radius. The current-carrying capability of the innermost and outermost rebco tapes in the star® wire degraded by 42% to 98% and the middle rebco tapes remained intact. This was also the case for the second and straight cable sample. After fabrication of the third cable sample, we observed only about 5% reduction in the current along the wire, measured at different locations inside the terminations. The results indicate that the differences in architecture or fabrication of the star® wires could have caused differences in critical current retention after cabling.

Cover page of Resonant Emittance Mixing of Flat Beams in Plasma Accelerators

Resonant Emittance Mixing of Flat Beams in Plasma Accelerators

(2024)

Linear colliders rely on high-quality flat beams to achieve the desired event rate while avoiding potentially deleterious beamstrahlung effects. Here, we show that flat beams in plasma accelerators can be subject to quality degradation due to emittance mixing. This effect occurs when the beam particles' betatron oscillations in a nonlinearly coupled wakefield become resonant in the horizontal and vertical planes. Emittance mixing can lead to a substantial decrease of the luminosity, the main quantity determining the event rate. In some cases, the use of laser drivers or flat particle beam drivers may decrease the fraction of resonant particles and, hence, mitigate emittance deterioration.

Cover page of Matched Guiding and Controlled Injection in Dark-Current-Free, 10-GeV-Class, Channel-Guided Laser-Plasma Accelerators

Matched Guiding and Controlled Injection in Dark-Current-Free, 10-GeV-Class, Channel-Guided Laser-Plasma Accelerators

(2024)

We measure the high-intensity laser propagation throughout meter-scale, channel-guided laser-plasma accelerators by adjusting the length of the plasma channel on a shot-by-shot basis, showing high-quality guiding of 500 TW laser pulses over 30 cm in a hydrogen plasma of density n_{0}≈1×10^{17}  cm^{-3}. We observed transverse energy transport of higher-order modes in the first ≈12  cm of the plasma channel, followed by quasimatched propagation, and the gradual, dark-current-free depletion of laser energy to the wake. We quantify the laser-to-wake transfer efficiency limitations of currently available petawatt-class lasers and demonstrate via simulation how control over the laser mode can significantly improve beam parameters. Using 21.3 J of laser energy, and triggering localized electron injection, we observed electron bunches with single, quasimonoenergetic peaks up to 9.2 GeV with charge extending beyond 10 GeV.