High Power Particle Beams for Radiation Science
Review of how we use particle beams to support our radiation effects science mission
SEARCH · Search NASA
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.
Review of how we use particle beams to support our radiation effects science mission
Explore the source record for details and available documents.
The analytical expressions for the electromagnetic potential generated from a focusing charged particle beam are indispensable in various beam physics problems. In this article, we review the theory in detail and point out the necessary assumptions made in the derivation.
The Zettawatt-Equivalent Ultrashort pulse laser System (ZEUS) is presently operational at the Gerard Mourou Center for Ultrafast Optical Science (CUOS) at the University of Michigan. ZEUS is a significant upgrade of the previous high power laser systems at CUOS and consists of two beamlines thatoperate in perfect synchronization. The 500 TW beamline became operational in 2023, 2 PW operation started in 2025 and full 3 PW power levels will be available in 2026. It is presently the highest power laser system in the US. In this grant the high field science group at CUOS has leveraged this unique high power laser facility to investigate laser wake field acceleration (LWFA) in ultra-high power laser plasma interactions and have shown how this can scale for future electron–positron colliders at high energy. The dual beam experimental configuration enables flexibility for many frontier experiments in laser-driven acceleration research, in particular, enabling extended channelling/acceleration experiments, positron generation/acceleration experiments and proof-of-principle transverse pumping “dephasingless” electron acceleration experiment and theory. LWFA may be able to miniaturize particle accelerators for high energy physics and also enable new sources of ultrafast, extreme brightness and precise x-rays for a wide variety of applications. In laser wake field acceleration, an electron bunch “surfs” on the electron plasma wave (the “wake field”) generated by the ponderomotive force of an intense laser. The plasma wave has a strong longitudinal electric field that stays in phase with the relativistic driver. A relativistic charged particle may, therefore, remain in phase with the accelerating field over long distances and gain ultra-relativistic energies. The accelerating electric field strength that the plasma wave can support can be many orders of magnitude higher than that of conventional accelerators, which makes laser wakefield acceleration an exciting prospect as an advanced accelerator concept. In this research project we have investigated the scaling of this mechanism to laser powers of 2 PW and have measured the x-ray emission and radio frequency emission resulting from the acceleration process. We have also performed theoretical investigation of mechanisms to scale laser driven accelerators to much higher energy using dephasingless acceleration processes.
Gallium Oxide (Ga 2 O 3 ), particularly in its β-phase, is attracting lots of interest for high-power and high-voltage electronics due to its wide bandgap, high breakdown field, and thermal stability. This study investigates the radiation response of Ga 2 O 3 Metal-Oxide-Semiconductor Field-Effect Transistors (MOSFETs) under Total Ionizing Dose (TID) and Displacement Damage (DD) conditions, which are critical for applications in radiation environments. Utilizing a dual-beam Focused Ion Beam and Scanning Electron Microscope setup, microscale analysis of radiation effects on individual devices is performed. The findings reveal distinct behaviors under TID and DD + TID conditions, with TID leading to threshold voltage shifts due to trapped charges, while DD results in decreased drive current attributed to increased carrier scattering from lattice defects. Notably, it is demonstrated that the TID effect can be mitigated through dynamic threshold voltage adjustments and that the predicted TID from ions calculated by Monte Carlo simulations overestimates actual TID due to unaccounted charge yield effects. In conclusion, this research enhances the understanding of Ga2O3 MOSFETs' performance in harsh radiation environments, providing insights for the design of robust electronic devices for space and nuclear applications.
Imaging methods driven by probes, electrons, and ions have played a dominant role in modern science and engineering. Opportunities for machine vision and AI that focus on consumer problems like driving and feature recognition, are now presenting themselves for automating aspects of the scientific processes. This proposal aims to enable and drive discovery in ultra-low energy implantation by taking advantage of faster processing, flexible control and detection methods, and architecture-agnostic workflows that will result in higher efficiency and shorter scientific development cycles. Custom microscope control, collection and analysis hardware will provide a framework for conducting novel in situ experiments revealing unprecedented insight into surface dynamics at the nanoscale. Ion implantation is a key capability for the semiconductor industry. As devices shrink, novel materials enter the manufacturing line, and quantum technologies transition to being more mainstream. Traditional implantation methods fall short in terms of energy, ion species, and positional precision. Here we demonstrate 1 keV focused ion beam Au implantation into Si and validate the results via atom probe tomography. We show the Au implant depth at 1 keV is 0.8 nm and that identical results for low energy ion implants can be achieved by either lowering the column voltage, or decelerating ions using bias – while maintaining a sub-micron beam focus. We compare our experimental results to static calculations using SRIM and dynamic calculations using binary collision approximation codes TRIDYN and IMSIL. A large discrepancy between the static and dynamic simulation is found that is due to lattice enrichment with high stopping power Au and surface sputtering. Additionally, we demonstrate how model details are particularly important to the simulation of these low-energy heavy-ion implantations. Finally, we discuss how our results pave a way to much lower implantation energies, while maintaining high spatial resolution.
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.
A ‘Phase Space’ is a mathematical representation of all possible physical states of matter of a physical volume or system. Particle beams at Fermilab exist in a six-dimensional phase space: three position dimensions, and three momentum dimensions. To ‘reconstruct’ that space means to take measurement data from virtual detectors along those beams and map out where particles are and how they’re moving, using computational methods. In this experiment, we simulated a particle beam (Fig 1) and used the position distribution of each particle at each detector to map out the magnitude of the beam’s initial momentum.
Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.
Not Available
Millicharged particles appear in several extensions of the standard model, but have not yet been detected. These hypothetical particles could be produced by an intense proton beam striking a fixed target. We use data collected in 2020 by the SENSEI experiment in the MINOS cavern at the Fermi National Accelerator Laboratory to search for ultrarelativistic millicharged particles produced in collisions of protons in the NuMI beam with a fixed graphite target. The absence of any ionization events with 3 to 6 electrons in the SENSEI data allow us to place world-leading constraints on millicharged particles for masses between 30 to 380 MeV. This work also demonstrates the potential of utilizing low-threshold detectors to investigate new particles in beam-dump experiments, and motivates a future experiment designed specifically for this purpose. Published by the American Physical Society 2024
Beam-intercepting devices such as beam windows and particle-production targets are critical components of accelerator target facilities for High Energy Physics (HEP) experiments. The high-power, pulsed structure of the particle beams used for these experiments leads to thermal shock and high-cycle fatigue in addition to radiation damage resulting from the accumulated particle fluence. This can lead to degradation of the target system s mechanical and thermal properties; considerably reducing their lifetimes and presenting substantial challenges to reliable operation of multi-MW class facilities. Recently several major accelerator facilities have been forced to operate at reduced power levels due to target survivability concerns. Furthermore, at Fermilab it is planned to increase the neutrino production beam power up to 2.4 MW in coming years. Therefore, timely R&D on the irradiated behavior of target system materials is critical to efficient operation of accelerator facilities and full utilization of recent accelerator power upgrades for HEP research. This talk will begin with an overview of high-power targetry, and the significant challenges presented by beam power increases expected for future HEP experiments. We will then cover several past materials irradiation studies that have been completed by the High-Power Targetry R&D group at Fermilab and its collaborators on common accelerator and target materials such as graphite, beryllium, titanium, and tungsten. Finally, we will conclude with a discussion of two novel materials investigations under way within the group; high-entropy alloys for beam window applications, and electrospun nanofibers to serve as particle production targets.
The interaction of charged particle beams with vacuum chamber components gives rise to electromagnetic wakefields, whose frequency-domain representation is known as beam coupling impedance. Geometric impedance arising from discontinuities and transitions in the vacuum chamber is the focus of this study. Minimizing this impedance is essential to mitigate adverse collective effects in modern storage rings operating with high-intensity particle beams. Accurate and reliable impedance simulations is a key factor of the vacuum chamber design. This paper presents the results of a convergence study of two widely used electromagnetic solvers, GdfidL and ECHO3D, applied to key vacuum-chamber components of the National Synchrotron Light Source II (NSLS-II) storage ring. Detailed comparisons are performed for several geometries, including flange absorbers, RF bellows, button-type beam position monitors, and an in-vacuum undulator (IVU). The results show notable differences in convergence and computational efficiency between the two codes. While GdfidL provides highly resolved results and serves as a common benchmark tool, ECHO3D yields consistent results with coarser meshes, significantly reducing simulation time and memory demands. Simulations with a full-geometry IVU model demonstrate that simplified taper-transition models can miss important impedance contributions. In conclusion, these findings provide practical guidelines for efficient and accurate impedance modeling to optimize design of vacuum chamber components for accelerators.
Abstract Advanced accelerator-based light sources such as free electron lasers (FEL) accelerate highly relativistic electron beams to generate incredibly short (10s of femtoseconds) coherent flashes of light for dynamic imaging, whose brightness exceeds that of traditional synchrotron-based light sources by orders of magnitude. FEL operation requires precise control of the shape and energy of the extremely short electron bunches whose characteristics directly translate into the properties of the produced light. Control of short intense beams is difficult due to beam characteristics drifting with time and complex collective effects such as space charge and coherent synchrotron radiation. Detailed diagnostics of beam properties are therefore essential for precise beam control. Such measurements typically rely on a destructive approach based on a combination of a transverse deflecting resonant cavity followed by a dipole magnet in order to measure a beam’s 2D time vs energy longitudinal phase-space distribution. In this paper, we develop a non-invasive virtual diagnostic of an electron beam’s longitudinal phase space at megapixel resolution (1024 × 1024) based on a generative conditional diffusion model. We demonstrate the model’s generative ability on experimental data from the European X-ray FEL.
A phase space is a mathematical representation of all possible physical states of a system. Particle beams at Fermilab exist within a six-dimensional (6D) phase space defined by three positional components, (x, y, z) and three momentum components, (px, py, pz). To reconstruct this space implies taking measurement data from detectors and mapping out particle behavior using computational methods. The beam detectors, however, are only able to detect spatial distribution among the events of the beam, therefore being limited to positional data. Also, due to the vast number of events in a particle beam, it is extremely difficult to analyze and differentiate every single one’s behavior. However, with Machine Learning (ML), which can distinguish between patterns and map out particle behavior more efficiently. We first used the particle beam software, G4beamline, to simulate a 10,000-event muon beam, adjusting parameters such as initial momentum magnitude (p¬0) and virtual detector position. Using ten virtual detectors, we analyzed p0 values such that minimum 9,990 events were analyzed by every detector. We then input the data from these beam simulations to a C++ program, that randomly selects 100 events, and creates a 2D histogram based on spatial distribution, detector position, and event intensity. This process is repeated 100 times to create 100 histograms per p0 value. These images were then input to a modified ResNet18 Convolutional Neural Network (CNN) for training, and to predict p0 from some unseen set of histograms. The model was accurate when trained on momentum increments of 5 MeV/c and provided with denser training samples around highly variable test values. These results displayed machine learning being able to accurately predict p0 from being trained on different particle behaviors.
The Proton Improvement Project-2 centers on building a new superconducting linear particle accelerator (Linac) at Fermilab. At the heart of the accelerator is the reference line, a critical system that defines the ideal path for the particle beam as it passes through magnets, RF cavities, and other beamline elements. Temperature stability is crucial for the reliable operation of RF components, such as mixers and filters. Fluctuations affect key performance parameters like conversion loss, isolation, and linearity. To mitigate any drift caused by ambient temperature changes, a heating plate assembly is utilized to maintain key components at a controlled temperature of 40°C. The system utilizes an aluminum 36”x36”x0.5” heat plate powered by a MOSFET-based control circuit, delivering approximately 460 W of thermal energy through a resistor array. Real-time temperature feedback is provided by a PT100 Resistance Temperature Detector (RTD), which interfaces with a Proportional–Integral–Derivative (PID) control algorithm to maintain closed-loop temperature regulation. The control signal actively modulates the gate voltage of an N channel MOSFET, dynamically adjusting power delivery in response to deviations from the temperature setpoint. Simulations and LTspice models validate the functionality and responsiveness of the circuit under varying conditions. The prototype has successfully demonstrated stable thermal control, paving the way for integration into the PIP-II infrastructure. The final design will feature an expanded resistor array, as well as communication with a PLC for continuous data acquisition and diagnostics. This work directly supports Fermilab’s broader mission by contributing to the stability and reliability of core accelerator systems, enhancing the precision of particle beam delivery for future physics experiments.
This article presents a study of single 𝜋 + production in neutrino neutral-current interactions (NC1𝜋 + ) using the FGD1 hydrocarbon target of the ND280 detector of the T2K experiment. We report the largest sample of such events selected by any experiment, providing the first new data for this channel in over four decades and the first using a sub-GeV neutrino flux. The signal selection strategy and its performance are detailed together with validations of a robust cross section extraction methodology. The measured flux-averaged integrated cross-section is 𝜎 = (6.07 ± 1.22) × 10 −41 cm 2 /nucleon, 1.3𝜎 above the NEUT v5.4.0 expectation.
We present the optimization and performance characterization of the second-generation Fermilab Constant Fraction Discriminator ASIC (FCFD), designed for the readout of AC-coupled low-gain avalanche detector (LGAD) strip-sensors. The FCFD is explicitly engineered to be insensitive to signal-amplitude variations, thereby removing the need for time-walk correction that is required in other LGAD time-stamping readout ASICs. This updated version, referred to as FCFD1.1, incorporates several enhancements over the first iteration to address key challenges in AC-coupled LGAD front-end design. We outline the primary readout-ASIC design considerations for these applications, describe the methodology used to evaluate critical sensor and system parameters, and summarize the additionally implemented features. Performance measurements using injected charge signals and minimum-ionizing particles in test-beams demonstrate a time resolution of approximately 40 ps and a position resolution of roughly \SI{15}{\micro\meter} when tested with beam particles uniformly over the active area of the sensor.