SEARCH · Engineering Papers

Results for “Charged particle motion”

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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Is there an exact magnetic moment for charged particle motion in a time-dependent, homogeneous magnetic field?

The non-perturbative guiding-centre model provides an exact alternative to full-orbit simulations of charged particle dynamics in situations where traditional guiding-centre theory may fail. We demonstrate that the charged particle motion in a homogeneous, time-varying magnetic field is a solvable example of the non-perturbative guiding-centre model. This entails showing that the exact magnetic moment of Qin and Davidson can be constructed to be asymptotic to the adiabatic invariant series of Kruskal. In contrast to the perturbative invariant, the exact invariant contains information about parametric resonances. These resonances destroy the conservation of the usual magnetic moment over very long times. This refutes some previous claims about the all-time invariance of the magnetic moment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

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↗

DOE OSTI2025

Nonperturbative Guiding Center Model for Magnetized Plasmas

Perturbative guiding center theory adequately describes the slow drift motion of charged particles in the strongly magnetized regime characteristic of thermal particle populations in various magnetic fusion devices. However, it breaks down for particles with large-enough energy. Here, we report on a data-driven method for learning a nonperturbative guiding center model from full-orbit particle simulation data. We show the data-driven model significantly outperforms traditional asymptotic theory in magnetization regimes appropriate for fusion-born α particles in stellarators, thus opening the door to nonperturbative guiding center calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Many-body interference in kagome crystals

When electrons in metals act collectively, they enable emergent phenomena and electronic functionalities that transcend the behaviour of individual particles. Coherent collective charge motion has so far been observed primarily in superconductors, in which it arises with the formation of Cooper pairs. Here we report experimental evidence for coherent charge transport in the normal state of the kagome metal CsV 3 Sb 5 , indicative of a distinct collective electronic state. The signature is a set of magnetoresistance oscillations in mesoscopic crystalline pillars under in-plane magnetic fields, with a periodicity determined by the number of magnetic flux quanta h/e threading between adjacent kagome layers—effectively forming an interlayer Aharonov–Bohm interferometer. The cooperative nature of this phenomenon is evidenced by a non-analytic angular dependence characterized by abrupt transitions between discrete oscillation frequencies and its persistence over length scales that exceed the single-particle mean free path. Notably, the oscillation amplitude matches other anomalous electronic responses reported in CsV 3 Sb 5 , pointing to an underlying mechanism that establishes intrinsic coherence. These findings shed new light on the debated nature of correlated order in kagome metals and establish CsV 3 Sb 5 as a platform for realizing long-range coherent charge transport in the absence of superconductivity—opening new directions for coherence in correlated electron systems beyond conventional models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

A Self Consistent 2D Simulation of Coherent Synchrotron Radiation Effects on Beam Dynamics

An increasing interest in high quality and high current electron beams necessitates a thorough understanding and prediction of coherent synchrotron radiation effects. The self-interaction of charged particles in a beam undergoing synchrotron motion is a physically significant process that is all too often computationally intensive with very little analytical results to rely on for the general case. The coherent spectrum of this interaction is of utmost importance to the design of free electron lasers (FELs) and an accurate assessment is imperative for their design. This work presents a novel implementation to the numerical simulation of charged particle beams. The simulation is a self-consistent approach including the self-fields generated by the beam of which coherent synchrotron radiation effects are of primary interest. A particle-in-cell model is used where a planar beam sampled by point particles is deposited on an encompassing grid at each timestep. The electromagnetic fields are calculated on the grid using the retarded potentials according to causality. The electromagnetic forces from the fields are interpolated on each particle which in turn advance in time. The simulation is benchmarked against well-established results for coherent synchrotron radiation effects. In addition, studies are provided that show the convergence of simulation results for increasing resolution. A study into the transverse beam size effects on beam dynamics is performed as well as a proof of concept where the simulation is used by a genetic algorithm to optimize the design parameters of a beam lattice. The results of these studies in tandem verify the efficacy of the simulation for its practical use in accelerator design or the study of synchrotron radiation effects

Duffin, Dallan [Old Dominion Univ., Norfolk, VA (U↗

DOE OSTI2025

Ponderomotive barriers in rotating mirror devices using static fields

Particularly for aneutronic fusion schemes, it is advantageous to manipulate the fuel species differently from one another and to expel ash promptly. The ponderomotive effect can be used to selectively manipulate particles. It is commonly a result of particle–wave interactions and has a complex dependence on the particle charge and mass, enabling species selectivity. If the plasma is rotating, e.g., due to E x B motion, the ponderomotive effect can be generated using static (i.e., time-independent) perturbations to the electric and magnetic fields, which can be significantly cheaper to produce than time-dependent waves. We propose that this feature can be particularly useful in rotating mirror machines where mirror confinement can be enhanced by rotation, both through centrifugal confinement and additionally through a ponderomotive interaction with a static azimuthal perturbation. We identify specific static perturbations that generate a ponderomotive barrier and other perturbations that can generate either a repulsive barrier or an attractive ponderomotive well, which can be used to attract particles of a certain species while repelling another. We identify the regimes in which the ponderomotive potential can enhance net plasma confinement and the regime in which plasma confinement is not enhanced. The viability of each of these effects is found to be dependent on the specifics of the rotation profile and the resultant dispersion relation in the rotating plasma.

Aneutronic fusion↗

DOE OSTI2026

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 (↗

DOE OSTI2026

Verification of nonperturbative guiding center theory in symmetric fields

We verify a recently-developed nonperturbative guiding center formalism to charged particle dynamics in fields with two-parameter continuous symmetry groups. This entails finding exact constants of motion, valid in the nonperturbative regime, that agree with Kruskal’s adiabatic invariant series to all orders in the perturbative regime, when the field scale length is large compared with a typical gyroradius. We demonstrate that the nonperturbative guiding center model makes exact predictions in these cases, even though it eliminates the cyclotron timescale, thereby establishing a theoretical baseline for performance of the nonperturbative formalism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Electrode strain dynamics in layered intercalation battery cathodes

Rechargeable batteries using electrodes based on intercalation chemistry exhibit notable cyclability, yet their performance still suffers from chemomechanical degradation. In this study, by combining a suite of operando microscopy methods, we explored electrode strain evolution and observed intricate particle cluster rearrangement under electrochemical stimuli. We show that early-stage strain accumulation in intercalation cathodes occurs during the period of interparticle charge transfer and redox reactions stemming from asynchronous coupling and decoupling between chemical (de)intercalation and physical grain motion. This interplay drives heterogeneous redox activity, localized charge equilibration, and multiscale strain cascades that propagate through an asynchronous network of chemical-mechanical interactions. Together, these findings reveal how collective particle dynamics and hierarchical strain transmission dictate electrode deformation and degradation in intercalation cathodes.

25 ENERGY STORAGE↗

DOE OSTI2025

Measurement of Local Electric Field Turbulence in Magnetically Confined Plasmas (Final Report)

This final report describes the development of high-bandwidth plasma fluctuation diagnostics, encompassing the design and implementation of Beam Emission Spectroscopy systems on several fusion experiments such as TFTR, DIII-D and NSTX, the advancement of Alpha particle and Ultrafast Charge-Exchange Recombination Spectroscopy, and the creation of advanced spectrometers for localized studies of electric and magnetic field fluctuations based on the Motional Stark Effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Time-resolved X-ray imaging of the current filamentation instability in solid-density plasmas

The streaming of energetic charged particles can magnetize astrophysical and laboratory plasmas via the current filamentation instability. Despite its importance, the experimental characterization of this instability has remained a challenge. Here, we report an experiment combining a high-intensity optical laser with a high-brightness X-ray free electron laser that successfully images the instability in solid-density plasmas with 200nm spatial and 50fs temporal resolution. We characterize the development of µm-scale filamentary structures and their evolution over tens of picoseconds through a non-linear merging process. The measured plasma density modulations and long merging time reveal the critical importance of space-charge effects and ion motion on this electron-driven instability. Supporting theoretical analysis and kinetic simulations help distinguish the relative role of space-charge and resistive effects. Our findings indicate that magnetic fields on the order of 10megagauss are produced, with important implications for transport and radiation emission of energetic particles in plasmas.

Schoenwaelder, Christopher↗

DOE OSTI2026

Time-resolved X-ray imaging of the current filamentation instability in solid-density plasmas

The streaming of energetic charged particles can magnetize astrophysical and laboratory plasmas via the current filamentation instability. Despite its importance, the experimental characterization of this instability has remained a challenge. Here, we report an experiment combining a high-intensity optical laser with a high-brightness X-ray free electron laser that successfully images the instability in solid-density plasmas with 200 nm spatial and 50 fs temporal resolution. We characterize the development of μm-scale filamentary structures and their evolution over tens of picoseconds through a non-linear merging process. The measured plasma density modulations and long merging time reveal the critical importance of space-charge effects and ion motion on this electron-driven instability. Supporting theoretical analysis and kinetic simulations help distinguish the relative role of space-charge and resistive effects. Our findings indicate that magnetic fields on the order of 10 megagauss are produced, with important implications for transport and radiation emission of energetic particles in plasmas.

Schoenwaelder, Christopher [SLAC National Accelera↗

DOE OSTI2026

Search for the chiral magnetic effect through beam energy dependence of charge separation using event shape selection

High-energy, heavy-ion collisions can create local domains of chirality-imbalanced quarks, reflecting the topological features of quantum chromodynamics. The chiral magnetic effect (CME) predicts an electric charge separation of quarks in such topological domains along the magnetic field ($\vec{B}$) generated by the passing of two high-Z nuclei. Here, we use a correlation observable Δ⁢𝛾$^{112}$ between charged meson pairs to detect the CME-induced charge separation and a novel event shape selection (ESS) method to mitigate the background effects related to elliptic flow (𝑣 2 ). The ESS method classifies events based on the emission pattern of final-state particles and determines Δ⁢𝛾$^{112}_{ESS}$ from the zero-flow limit. We reconstruct the $\vec{B}$ field direction from the spectator nucleons, which minimizes backgrounds unrelated to the collective motion of the system. In this work, we report the measurements of Δ⁢𝛾$^{112}$ and a background indicator Δ⁢𝛾$^{132}$ in Au+Au collisions from the Brookhaven National Laboratory Relativistic Heavy Ion Collider (RHIC) Beam Energy Scan phase II and at the top RHIC energy. After background suppression, Δ⁢𝛾$^{132}_{ESS}$ aligns with zero, and Δ⁢𝛾$^{112}_{ESS}$ is reduced to no more than 20% of Δ⁢𝛾$^{112}$. We observe a finite residual charge separation with 2.5⁢𝜎, 3⁢𝜎, and 3.2⁢𝜎 significance in the 20–50% centrality range of Au + Au collisions at 11.5, 14.6, and 19.6 GeV. The results at 17.3 and 27 GeV also show positive values but with a lower significance of 1.3⁢𝜎 and 1.1⁢𝜎, respectively. The corresponding Δ⁢𝛾$^{112}_{ESS}$ values at 7.7, 9.2, and 200 GeV are consistent with zero within uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Studies of laser stimulated photodetachment from nanoparticles for particle charge measurements

Determining nanoparticle charge is more challenging than that for microparticles due to change in the particle size during the synthesis substantial plasma property variations, and difficulties in visualizing individual particles, rendering conventional microparticle charge diagnostics ineffective in dusty plasma. In this work, we utilized laser-stimulated photodetachment (LSPD) to deduce the mean charge of nanoparticles. Nanoparticles were grown in an Ar/C 2 H 2 mixture using a capacitively coupled RF discharge and the LSPD induced changes in the electron current monitored by a cylindrical Langmuir probe. LSPD signals were obtained and analyzed across different dust growth phases. The prolonged decay of electron current pulses was attributed to the presence of residual negative ions, caused by the effective electrostatic trapping of these ions and the potential post—LSPD re-formation of new ones. The particle charge was estimated by combining the laser-stimulated photodetachment signal from the probe with the dust density obtained from laser-light extinction using the measured nanoparticle size distribution. For a nanoparticles size range of approximately 100–250 nm and mean diameter of $d$ p ∼154.37 nm, the effective mean charge was estimated to be $\langle$$Q$$\rangle$ d $≈$ 37 elementary charge units. The measured charge values are lower than those predicted by orbital motion limited theory, which may be attributed to significant electron depletion in the nanodusty plasma. LSPD results in Ar/C 2 H 2 nano-dusty plasma confirm the applicability of this method for estimating individual nanoparticle charges. However, it has also been demonstrated that electron detachment from residual background negative ions can influence the detachment current decay and must be carefully considered.

dust particle charge↗

DOE OSTI2025

New model for the ion collection by cylindrical probes over a wide range of collisionality

Langmuir probes remain one of the most important diagnostic tools for plasma processing applications. Modern probe analysis usually relies on the electron current part of the Langmuir probe characteristic using the Druyvesteyn method. However, for electronegative plasmas or for discharges containing dust the analysis of the ion current attracted by the probe can be desirable to determine the ion density. But, even at low pressures of a few Pa, the ion current is affected by collisions due to the large cross section for charge exchange. Available theories for collisional or collision-enhanced ion currents onto probes are complex and not well validated. Thus, in this contribution, we compare available collisional probe theories for the ion current to results of particle-in-cell (PIC) simulations. To this end, the probe surrounded by a semi-infinite plasma is simulated using a modified version of the open-source code EDIPIC. A dataset of currents for different neutral gas pressures is obtained and compared to the different theories from the literature. Based on these results, we propose a simpler and more intuitive model for the ion current collected by the probe, based on the model of Gatti and Kortshagen (Phys. Rev. E 78, 046402, 2008), developed for the charging of dust particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Dynamics of McMillan mappings II. axially symmetric map

Here, in this article, we investigate the transverse dynamics of a single particle in a model integrable accelerator lattice, based on a McMillan axially-symmetric electron lens. Although the McMillan e-lens has been considered as a device potentially capable of mitigating collective space charge forces, some of its fundamental properties have not been described yet. The main goal of our work is to close this gap and understand the limitations and potentials of this device. It is worth mentioning that the McMillan axially symmetric map provides the first-order approximations of dynamics for a general linear lattice plus an arbitrary thin lens with motion separable in polar coordinates. Therefore, advancements in its understanding should give us a better picture of more generic and not necessarily integrable round beams. In the first part of the article, we classify all possible regimes with stable trajectories and find the canonical action-angle variables. This provides an evaluation of the dynamical aperture, Poincaré rotation numbers as functions of amplitudes, and thus determines the spread in nonlinear tunes. Also, we provide a parameterization of invariant curves, allowing for the immediate determination of the map image forward and backward in time. The second part investigates the particle dynamics as a function of system parameters. We show that there are three fundamentally different configurations of the accelerator optics causing different regimes of nonlinear oscillations. Each regime is considered in great detail, including the limiting cases of large and small amplitudes. In addition, we analyze the dynamics in Cartesian coordinates and provide a description of observable variables and corresponding spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Impact of solvation on the electronic resonances in uracil

Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Skyrmion-skyrmionium phase separation and laning transitions via spin-orbit torque currents

Many driven binary systems can exhibit laning transitions when the two species have different mobilities, such as colloidal particles with opposite charges in electric fields. Another example is pedestrian or active-matter systems, where particles moving in opposite directions form a phase-separated state that enhances the overall mobility. In this work, we use atomistic simulations to demonstrate that mixtures of skyrmions and skyrmioniums also exhibit pattern formation and laning transitions. Skyrmions move more slowly and at a finite skyrmion Hall angle compared to skyrmioniums, which move faster and without a skyrmion Hall effect. At low drives, the system forms a partially jammed phase where the skyrmioniums are dragged by the surrounding skyrmions, resulting in a finite angle of motion for the skyrmioniums. At higher drives, the system transitions into a laned state, but unlike colloidal systems, the lanes in the skyrmion-skyrmionium mixture are tilted relative to the driving direction due to the intrinsic skyrmion Hall angle. In the laned state, the skyrmionium angle of motion is reversed when it aligns with the tilted lane structure. At even higher drives, the skyrmioniums collapse into skyrmions. Below a critical skyrmion density, both textures can move independently with few collisions, but above this density, the laning state disappears entirely, and the system transitions to a skyrmion-only state. We map out the velocity and Hall responses of the different textures and identify three distinct phases: partially jammed, laned, and skyrmion-only moving crystal states. In conclusion, we compare our results to recent observations of tilted laning phases in pedestrian flows, where chiral symmetry breaking in the particle interactions leads to similar behavior.

36 MATERIALS SCIENCE↗