Active Shielding Particle Pusher (ASPP): Charged-Particle Tracking Through Electromagnetic Fields
Explore the source record for details and available documents.
D Fry↗
SEARCH · Engineering Papers
Results for “Particle pusher”
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.
New particle pusher with hadronic interactions for modeling multimessenger emission from compact objects
We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Special unitary particle pusher for extreme fields
Not Available
Gordon, D. F. (ORCID:0000000308448988)↗
A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm
A customized finite-difference field solver for the particle-in-cell (PIC) algorithm that provides higher fidelity for wave-particle interactions in intense electromagnetic waves is presented. In many problems of interest, particles with relativistic energies interact with intense electromagnetic fields that have phase velocities near the speed of light. Numerical errors can arise due to (1) dispersion errors in the phase velocity of the wave, (2) the staggering in time between the electric and magnetic fields and between particle velocity and position and (3) errors in the time derivative in the momentum advance. Errors of the first two kinds are analyzed in detail. It is shown that by using field solvers with different -space operators in Faraday’s and Ampere’s law, the dispersion errors and magnetic field time-staggering errors in the particle pusher can be simultaneously removed for electromagnetic waves moving primarily in a specific direction. Here, the new algorithm was implemented into Osiris by using customized higher-order finite-difference operators. Schemes using the proposed solver in combination with different particle pushers are compared through PIC simulation. It is shown that the use of the new algorithm, together with an analytic particle pusher (assuming constant fields over a time step), can lead to accurate modeling of the motion of a single electron in an intense laser field with normalized vector potentials, eA / mc 2 , exceeding for typical cell sizes and time steps.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Accurately simulating nine-dimensional phase space of relativistic particles in strong fields
Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. I. Curvilinear Special Relativistic Module
Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with evolving, GPU-centric computing infrastructures. It achieves hardware-agnostic performance portability across various GPU and CPU architectures using the Kokkos library. Crucially, Entity maintains a high standard for usability, clarity, and customizability, offering a robust and easy-to-use framework for developing new algorithms and grid geometries, which allows extensive control without requiring edits to the core source code. This paper details the core general-coordinate special relativistic module. Entity is the first PIC code designed to solve the Vlasov–Maxwell system in general coordinates, enabling a coordinate-agnostic framework that provides the foundational structure for straightforward extension to arbitrary coordinate geometries. The core methodology achieves numerical stability by solving particle equations of motion in the global orthonormal Cartesian basis, despite using generalized coordinates like Cartesian, axisymmetric spherical, and quasi-spherical grids. Charge conservation is ensured via a specialized current deposition technique using conformal currents. The code exhibits robust scalability and performance portability on major GPU platforms (AMD MI250X, NVIDIA A100, and Intel Max Series), with the 3D particle pusher and the current deposition operating efficiently at about 2 ns per particle per time step. Functionality is validated through a comprehensive suite of standard Cartesian plasma tests and the accurate modeling of relativistic magnetospheres in curvilinear axisymmetric geometries.
Hakobyan, Hayk [Flatiron Institute, New York, NY (↗
Symplectic neural network and its application to charged particle dynamics in electromagnetic fields
Recently, machine learning models have shown many successes in various applications in science and technology. In this work, we focus on the charged particle dynamics, with the development of a class of symplectic neural networks, including a linear version, SympMat, and a nonlinear version, HénonNet. Both are designed to preserve the structure of Hamiltonian systems. We show that they can be used to model relevant Hamiltonian systems of interest in plasma physics and astrophysics, for linear and nonlinear charged particle dynamics, with the potential to bridge multi-scale simulations. These symplectic neural networks are adapted to the applications in plasma simulations and particle-wave interaction with parametric dependence and periodicity, where we have investigated their performance and accuracy. In particular, SympMat is shown to outperform the traditional Boris particle pusher down to the sub-gyroperiod scale in the case of charged particles in uniform magnetic fields. HénonNet successfully predicts the hot electron distribution, which is validated against theoretical results. These results highlight the potential of symplectic neural networks as a trajectory integrator for particle-in-cell simulations or a fast surrogate to replace conventional numerical schemes.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
An implicit particle code with exact energy and charge conservation for studies of dense plasmas in axisymmetric geometries
A collisional particle code based on implicit energy- and charge-conserving methods in axisymmetric geometries is presented. A new particle pusher for axisymmetric systems is introduced that is compatible with exact energy and charge conservation and yields improved accuracy compared to other methods. How to appropriately treat all aspects of the algorithm near the r = 0 axis of symmetry is described in detail. Here, the axisymmetric model is verified by simulating the free expansion of a plasma sphere in 2D cylindrical and 1D spherical geometries. The algorithm's ability to study the dynamic compression of a dense plasma is illustrated by simulating the dynamic Z-pinch in 1D cylindrical geometry.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Exploring exponential time integration for strongly magnetized charged particle motion
Not Available
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Nyström type exponential integrators for strongly magnetized charged particle dynamics
Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.
general physics↗
Hybrid simulation of energetic particles interacting with magnetohydrodynamics using a slow manifold algorithm and GPU acceleration
The hybrid method combining particle-in-cell and magnetohydrodynamics can be used to study the interaction between energetic particles and global plasma modes. In this paper we introduce the M3D-C1-K code, which is developed based on the M3D-C1 finite element code solving the magnetohydrodynamics equations, with a newly developed kinetic module simulating energetic particles. The particle pushing is done using a new algorithm by applying the Boris pusher to the classical Pauli particles to simulate the slow-manifold of particle orbits, with long-term accuracy and fidelity. The particle pushing can be accelerated using GPUs with a significant speedup. The moments of the particles are calculated using the δƒ method, and are coupled into the magnetohydrodynamics simulation through pressure or current coupling schemes. Several linear simulations of magnetohydrodynamics modes driven by energetic particles have been conducted using M3D-C1-K with the δƒ method, including fishbone, toroidal Alfvén eigenmodes and reversed shear Alfvén eigenmodes. Good agreement with previous results from other eigenvalue, kinetic and hybrid codes have been achieved.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Demonstrating imaging plate detector stacks for proton radiography using exploding pusher capsules
A new detector design has been fielded on the Omega-60 laser facility using a D 3 He backlighter capsule to radiograph a magnetized experiment. D 3 He capsules produce a low flux (~ 10 8 total yield in 4π) of protons requiring detectors with close-to single-particle sensitivity to be employed. The new detector stack consisted of both imaging plate (IP) and CR-39 to detect 14.7 MeV protons where, historically, only CR-39 has been used. IP is sensitive to significant contributions in signal from both x-rays and protons so additional filtering has to be added to attenuate the x-rays. Here, the signals and features observed from a single shot are detected by both IP and CR-39, giving confirmation that it is protons creating the signatures on the IP. Measurements of PSL/pixel are used to calculate approximate on-shot proton yields and agree with commonly measured yields. A second IP is used at the rear of the stack to measure, and subsequently remove, a background signal. The spatial resolution at the IP is limited, primarily, by a range of proton energies with varying deflections being recorded at the detector, causing potential ‘blurring’ of features. We find for this setup that the blurring effect is less than the smallest identifiable features observed at the detector, and therefore of minimal impact. A large benefit to this new stack design is that IP can be scanned and processed on much faster timescales than CR-39 allowing for prompt shot feedback. Future designs and modifications to the stack design fielded on this experiment could help improve the contrast of the radiograph on IP, as well as detecting protons produced at alternative energies.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
First measurement of the 10B(α,n)13N reaction in an inertial confinement fusion implosion at the National Ignition Facility: Initial steps toward the development of a radiochemistry mix diagnostic
We report the first measurement of the 10B(α,n)13N reaction in a polar-direct-drive exploding pusher (PDXP) at the National Ignition Facility (NIF). This work is motivated by the need to develop alternative mix diagnostics, radiochemistry being the focus here. The target is composed of a 65/35 at. % deuterium–tritium (DT) fill surrounded by a roughly 30 μm thick beryllium ablator. The inner portion of the beryllium ablator is doped with 10 at. % of 10B. Radiation-hydrodynamics calculations were performed in 1D to optimize both the remaining boron rho-R and the DT neutron yield. A charged-particle transport post-processor has been developed to study α-induced reactions on the ablator material. Results indicate a large 13N production from α-induced reactions on 10B, measurable by the radiochemical analysis of a gaseous samples system at the NIF. The PDXP target N201115-001 was successfully fielded on the NIF, and nitrogen from the 10B(α,n)13N reaction was measured. The 13N production yield, as well as the DT neutron yield, was, however, lower than expected. Some of the reduced yields can be explained by the oblate shape, but the ratios of the various radiochemical signals are not commensurate with expectations based on a simple reduction of the 1D results. Preliminary 2D radiation-hydrodynamics computations are consistent with the experimental measurements, and work is ongoing to extend the radiochemistry analysis into higher dimensions.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Command of active and responsive elastomers by topological defects and patterns
The project resulted in the development of stimuli-responsive elastomer coatings formed by photopolymerized liquid crystal molecules. The molecular orientation of the liquid crystal elastomers is coupled to rubber-like elasticity. The project developed an approach to produce complex patterns of molecular orientation in elastomers by employing photopatterning technique based on plasmonic metamasks that convert unpolarized incoming light into a transmitted light beam with spatially-varying linear polarization. The polarization-modulated light beam aligns the film of azodye molecules, which serves as the substrate for liquid crystal elastomer coating. The alignment of the azodye molecules follows the polarization pattern with a phase shift of the alignment direction by 90 degrees. The azodye substrate imposes orientation of molecules in the adjacent liquid crystal layer which is preserved after the monomers are polymerized into a liquid crystal elastomer. By using differently aligned substrates one could produce various patterns of molecular orientation in the elastomer, which dictate the mechanical properties of the elastomers and their response to external stimuli such as temperature, humidity, and ultraviolet irradiation. In particular, the predesigned molecular orientation produces deterministically defined topography response of the liquid crystal coatings, in which their free surface changes from flat to either locally elevated or locally compressed, depending on whether the pre-inscribed molecular orientation is predesigned with a splay or bend deformation. The research established that the mechanism of this deterministic relation between the dynamic topography and pre-inscribed molecular orientation is rooted in the change of the tensorial order parameter that characterizes the degree of orientational order, and in the contraction/expansion of polymer molecules in response to the order parameter changes. In particular, lowering the order causes contraction of polymer globules. The work demonstrated that the dynamic topography of elastomer coatings with three-dimensional pattern of orientational order includes up and down motion of the material (along the normal to the coating) and also shifts in lateral directions. We demonstrate that the deterministic relationship between the orientational patterns and variations of coatings’ profile is caused by forces that are mathematically equivalent to active forces in the system of “pullers” or “pushers” in active matter, demonstrating universality of the out-of-equilibrium behavior. The project resulted in theoretical models that predictively describe the dynamic response of coatings. The topography of liquid crystal elastomers is sensitive to stimuli such as temperature, ultraviolet irradiation, and humidity. The applicability of the patterned approach to produce dynamic coatings was expanded from nematic elastomers to their smectic counterparts. The project combined a battery of imaging techniques to unveil the dynamic properties of the coatings, ranging from the state-of-the-art dynamic holographic microscopy to fluorescent confocal polarizing microscopy and polscope microscopy. The research uncovered potential applications of dynamic elastomer coatings, such as control placement of colloidal particles, engineering of living tissues with alignment patterns of cells that follow the orientational order of elastomer coatings. The approach developed in the project can be used to design programmable dynamic coatings with functionalities that mimic biological tissues such as skin.
36 MATERIALS SCIENCE↗
Entropy generation from hydrodynamic mixing in inertial confinement fusion indirect-drive targets
The increase in entropy from the physical mixing of two adjacent materials in inertial confinement fusion (ICF) implosions and gas-filled hohlraums is analytically assessed. An idealized model of entropy generation from the mixing of identical ideal-gas particles across a material interface in the presence of pressure and temperature gradients is applied. Physically, mix-driven entropy generation refers to the work done by the gases in expanding into a larger common volume from atomic mixing under the condition of no internal energy change, or work needed to restore the initial unmixed state. Furthermore, the effect of a mix-generated entropy increase is analytically shown to lead to less compression of the composite ICF fluid under adiabatic conditions. The amount of entropy generation is estimated to be ~10 J for a Rayleigh–Taylor-induced micrometer-scale annular mixing layer between the solid deuterium–tritium fuel and (undoped) high-density carbon pusher of an imploding capsule at the National Ignition Facility (NIF). This level of entropy generation is consistent with lower-than-expected fuel compressions measured on the NIF [Hurricane et al., Phys. Plasmas 26, 052704 (2019)]. The degree of entropy increase from mixing of high-Z hohlraum wall material and low-Z, moderate- to high-density gas fills is estimated to lead to ~100 kJ of heat generation for NIF-scale experiments [Moody et al., Phys. Plasmas 21, 056317 (2014)]. This value represents a significant fraction of the inferred missing x-ray drive energy based on observed delays in capsule implosion times compared with mainline simulations [Jones et al., Phys. Plasmas 19, 056315 (2012)].
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A pushered capsule implosion as an alternate approach to the ignition regime for inertial confinement fusion
We report in inertial confinement fusion, the threshold for ignition is a highly dynamic quantity as the sources and sinks of power in the hot spot can vary rapidly. In this article, we consider the ignition condition as a race between heating and disassembly rates and make use of a prior solution to the fusion hot-spot thermodynamics to develop a Lawson-like ignition criteria for pressure × confinement time (p-τ) vs temperature. Low-Z capsule designs reach the temperature for this threshold using as much of the shell as feasible as ablator but then are limited in τ by low stagnated mass. An alternate approach, the pushered single shell (PSS) design [D. D.-M. Ho, S. MacLaren, and Y. Wang, “High-yield implosions via radiation trapping and high rho-R,” paper presented at the 60th Annual Meeting of the APS Division of Plasma Physics, 2018], introduces a dense inner layer of Mo-Be alloy that is smoothly graded outward to pure Be, increasing the confinement time at stagnation and lowering the temperature requirement at the ignition threshold. Here, we describe a PSS ignition design for the National Ignition Facility and use the theory as well as simulations to compare it with the low-Z capsule approach. Additionally, we show how an adjustment to the design is used to anticipate the effects of mixing at the fuel–ablator interface.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗