Search NASA⌕ Search

SEARCH · Search NASA

Results for “charged particle transport”

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

Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma

The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The SHMS 11 GeV/c spectrometer in Hall C at Jefferson Lab

The Super High Momentum Spectrometer (SHMS) has been built for Hall C at the Thomas Jefferson National Accelerator Facility (Jefferson Lab). With a momentum capability reaching 11 GeV/ c , the SHMS provides measurements of charged particles produced in electron-scattering experiments using the maximum available beam energy from the upgraded Jefferson Lab accelerator. The SHMS is an ion-optics magnetic spectrometer comprised of a series of new superconducting magnets which transport charged particles through an array of triggering, tracking, and particle-identification detectors that measure momentum, energy, angle and position in order to allow kinematic reconstruction of the events back to their origin at the scattering target. The detector system is protected from background radiation by a sophisticated shielding enclosure. The entire spectrometer is mounted on a rotating support structure which permits measurements to be taken with a large acceptance over laboratory scattering angles from 5.5° to 40°, thus allowing a wide range of low cross-section experiments to be conducted. Finally, these experiments complement and extend the previous Hall C research program to higher energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Multi-fidelity equations of state and transport coefficient datasets for pulsed-power applications

Reliably simulating experiments relevant to the National Nuclear Security Administration (NNSA) requires a detailed description of material properties across a wide range of conditions. Such properties include the equations of state, charged-particle transport coefficients, and optical properties like the opacity. Together, these properties make up the material models used in radiation-magnetohydrodynamic simulations of nuclear fusion experiments. Many of these models do not incorporate uncertainties in the data used to produce them. It is unknown whether these uncertainties significantly impact the interpretation of simulation results and diagnostics. The purpose of this work is to quantify how such uncertainties impact simulations of pulsed-power experiments. We accomplished this task by first assessing discrepancies between approaches used to generate the data. This included bringing together members of the high-energy-density community spanning the three NNSA laboratories and multiple universities. Then, using these data, we developed a general framework that systematically incorporates physical uncertainties within the material models suitable for uncertainty quantification analyses. The framework utilizes machine learning, Bayesian inference, and incorporates multi-fidelity datasets. We demonstrated the framework by quantifying the impact that material model uncertainties have on simulations of pulsed-power experiments underway on Z at Sandia National Laboratories. As a result of this work, we discovered that modest uncertainties in material models (roughly 20%) correspond to significant uncertainties in the outputs from simulations. Our framework has enabled rapid construction of material models through an automated procedure and allows for the generation of material models of interest to the NNSA.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Transport coefficients of warm dense matter from Kohn-Sham density functional theory

We present a comprehensive study of transport coefficients including DC electrical conductivity and related optical properties, electrical contribution to the thermal conductivity, and the shear viscosity via ab initio molecular dynamics and density functional theory calculations on the “priority 1” cases from the “Second Charged-Particle Transport Coefficient Workshop” [Stanek et al., Phys. Plasmas (to be published 2024)]. The purpose of this work is to carefully document the entire workflow used to generate our reported transport coefficients, up to and including our definitions of finite size and statistical convergence, extrapolation techniques, and choice of thermodynamic ensembles. In pursuit of accurate optical properties, we also present a novel, simple, and highly accurate algorithm for evaluating the Kramers–Kronig relations. These heuristics are often not discussed in the literature, and it is hoped that this work will facilitate the reproducibility of our data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A Primer on Nuclear "Recoil" Data

For many years the ENDF-6 format has existed to contain evaluated nuclear data. One primary driver for the format and the data is neutron transport calculations. Evaluated data can be processed by a code like NJOY into either continuous energy form (ACE format) for Monte Carlo codes like MCNP or into multi-group form (NDI tables) for deterministic codes like Partisn. The reaction cross sections are found in the MF 3 section of the ENDF-6 format. If the reaction produces one or more neutrons as outputs, then secondary neutron data must also be given in MF 4,5 or 6 sections of the format. Energy and angular distributions for the output neutrons must be given in one of several available formats. The most general formats are found in MF 6. MF 4 is for angular distributions, MF 5 is for Energy distributions, and MF 6 contains both. In recent years, more interest has developed in the other output particles (i.e., the “recoil” particles) from neutron induced reactions. This has been driven by interest in charged particle transport and in more specialized partial kermas. ( A separate total neutron kerma has been available for a long time.) Partial kermas are a breakdown of the Kinetic Energy Released into the MAterial by output particle or by neutron reaction. The sum of the partial kermas should be equal to the total kerma on a group-wise basis. The ENDF-6 format is general enough for these new data requirements. Ideally, evaluated data would exist for every output particle (including the secondary neutrons) from every neutron-induced reaction. In the MF 6 format section, data for multi-particle outputs may be entered using the LAW =1 option. This option explicitly allows energy and angular output distributions for each particle produced in the reaction. Such information should preserve the balance between partial kermas and the total kerma at the groupwise level as well as the individual particle averaged energies. For simpler 2-body reactions, LAW = 2 is available. This allows the evaluation to specify only enough data to specify the 2-body reaction fully. Such a simplification does not exist for 3 (or more) body breakups. A simplified form of LAW 2, i.e., LAW 3, also exists to generate approximate recoil output distributions in the absence of full data. However, evaluated data does not generally exist at this fine granularity for reactions involving more than 2 output particles ( e.g., the 3-body break-up reaction). When the multi-body detailed output distribution data is not available, LAW=6 may be employed in NJOY. LAW 6 produces approximate output distributions for all the particles produced in the reaction. Just like in the 2-body case, the smaller particles will generally carry off more energy. At any particular energy, the balance between the sum of the partial kermas and the total kerma will not necessarily be preserved. However, the average energy of each output particle from a reaction is preserved.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

OpenEdge: A collaborative, open-source, multi-purpose direct simulation Monte Carlo for plasma simulation in magnetic fusion environments

OpenEdge is a collaborative, open-source, object-oriented Direct Simulation Monte Carlo (DSMC) code, designed specifically for plasma simulations in magnetic fusion environments. Here, the code features include advanced structures, robust capabilities, and an effective parallelization strategy, all of which significantly enhance performance. It includes specialized modules for managing complex particle interactions, including collisions, ionization/recombination, and reflection/sputtering. Benchmarks and performance analyses have confirmed its efficiency and scalability. Versatile and adaptable, OpenEdge is applied across a broad spectrum of plasma-material interaction studies and charged particle transport in various fusion research settings.

Boundary plasma↗
DOE OSTI2021

Parameter space for magnetization effects in high-energy-density plasmas

Magnetic fields are well known to affect the evolution of fluids via the J × B force, where J is the current density and B is the magnetic field. This force leads to the influence of magnetic fields on hydrodynamics (magnetohydrodynamics). Magnetic fields are often neglected in modeling of high-energy-density plasmas, since J × B is very small compared with the plasma pressure gradients. However, many experiments lie in a separate part of parameter space where the plasma is indirectly affected via magnetization of the heat flux and charged particle transport. This is true even for initially unmagnetized plasmas, since misaligned density and temperature gradients can self-generate magnetic fields. By comparing terms in the induction equation, we go on to estimate the regions of parameter space where these self-generated fields are strong enough to affect the hydrodynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

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↗
DOE OSTI2025

Monte Carlo toolkit for designing and validating step-range-filter spectrometer designs

Here, we present a Monte Carlo toolkit for validating step range filter (SRF) spectrometer designs. Geant4 is used to transport charged particles through the SRF filters to generate synthetic SRF data that include realistic CR-39 effects. Synthetic SRF spectra generated by this method inherently account for instrument response and allow for the quantification of SRF performance before shots. The usefulness of this toolkit is demonstrated through its application to a number of problems. A new broadband SRF for the ∼10 MeV wide 3He3He proton spectrum is validated, and an analysis method for analyzing 3He3He-p SRF data that accounts for instrument response is put forth. In addition, an SRF design for the compact recoil-proton spectrometer (CRS) on the Z-machine is validated. Finally, a new calibration technique for the DD-p SRF is proposed and validated.

Johnson, T. M. (ORCID:0000000193032949)↗
DOE OSTI2026

Utilizing the deuterium-tritium fusion resonance to diagnose thermal runaway in igniting plasmas

For high-efficiency inertial confinement fusion implosions, it is predicted that a burning hot spot will successfully encompass all surrounding fuel and then transition into a thermal runaway where the internal energy increase from fusion occurs on a timescale faster than the expansion of the fuel is able to quench the fusion chain reaction after ignition occurs. Observation of this dynamic phase transition would indicate distinct burn properties and indicate an implosion's robustness. A technique for diagnosing the presence of thermal runaway from measurements of nuclear reaction history is presented. The technique is based on taking the logarithmic derivative of the nuclear reaction history, called the 𝛼 curve, and allowing a mathematical decoupling of the mass, volume, and thermal reactivity in the fusion reaction rate equation. During thermal runaway, where the thermal temperature dominates the burn dynamics, a maximum in the 𝛼 curve is found where there is a maximum in the first derivative of the thermal fusion reactivity, an effect to the deuterium-tritium (DT) fusion cross-section resonance. This provides a distinct signature related to the fundamental nature of the DT fusion nuclear resonance and signifies the transition into the fusion thermal instability. Impacts of charged particle transport on the effect are also assessed and the analytical formulas are compared and found to be in agreement with radiation hydrodynamic codes.

high-energy-density plasmas↗
DOE OSTI2022

Implementation and Results of Extended-MHD Post-processing of xRAGE Simulations

Magnetic fields are known to affect the evolution of conducting fluids via the J× B force, where J is the current density and B is the magnetic field. Typically, J× B in high energy density plasmas (HEDP) systems is considered small compared to the plasma pressure gradients. However, many experiments lie in a part of the parameter space where the plasma is indirectly affected via magnetization of the heat flux and charged particle transport. This is true even for initially unmagnetized plasmas, since misaligned density and temperature gradients can self-generate the so-called Biermann magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Field Study of Nighttime Leakage Currents in Bifacial PV Modules: Correlation with Atmospheric Electric Field Data

Leakage currents measured on PV modules in the field originate from a potential difference between the modules' frame and the cells. They can be a relative indicator of Potential-Induced Degradation (PID) severity, especially when comparing the same module design in a different environment. As modules are not operating at night, no leakage current should be observed but our team has reported several events of nighttime leakage currents on bifacial PV modules. These events have been firstly observed during a thunderstorm that are characterized by strong atmospheric electrical field values. This lead us to believe that nighttime leakage currents could originate from the atmospheric electric charges. In this paper, we correlate nighttime leakage currents measured on bifacial PV modules with field mill data to identify the origin of nighttime leakage currents. Our results show that so far, no leakage currents at night occur when the atmospheric electric field is between 0 and 150–200 V/m (standard value for fair weather). As soon as the atmospheric electric field is out of this range, leakage currents are observed with or without rain involved. This suggests a transport of charged particles from the atmosphere to the modules' frame. A combination of heavy rain with strong atmospheric electric field results into high nighttime leakage currents with a magnitude up to 8 times higher than what observed during the day with -1500V applied. This is explained by an easier transport of the charged particles through the water droplets. Based on these results, leakage currents observed during the day might not be only due to the inherent potential difference between the frame and the cells depending on the atmospheric electric field activity. We believe that it should be taken into account in PID studies.

14 SOLAR ENERGY↗
DOE OSTI2025

Polaron Delocalization and Transport in Doped Graphene Nanoribbon Thin Films

Graphene nanoribbons (GNRs) are quantum-confined π-conjugated monolayer semiconductors with attractive properties for optoelectronic applications. However, the ground- and excited-state properties of charge carriers in GNRs are still poorly understood, particularly with regards to the coupling between charges and the GNR lattice and the degree to which this coupling impacts local and macroscopic charge transport. To address this issue, we systematically correlate carrier density-dependent charge transport with spectroscopic modulations in chemically doped thin films of armchair graphene nanoribbons (9-aGNRs). This study combines Fourier transform infrared (FTIR) and ultraviolet−visible−near-infrared (UV−vis−NIR) spectroscopy with both local and macroscopic conductivity measurements to arrive at a full and self-consistent picture of transport in doped GNR thin films. Using three different molecular p-type dopants (i.e., oxidants), we demonstrate that hole polarons are the dominant quasi-particle determining charge transport in GNRs and that the degree of polaron delocalization depends sensitively on the dopant and the hole density. For all three dopants, the local conductivity probed by microwave spectroscopy substantially exceeds the long-range conductivity obtained by four-point probe measurements. Interestingly, the dopant size substantially influences charge transport at high hole densities. We ascribe this effect to different propensities for forming bipolarons with lower mobilities than polarons. Comparison of GNR transport and spectral properties to other prototypical π-conjugated semiconductors (e.g., semiconducting polymers or carbon nanotubes) benchmark the charge transport properties of GNR thin films for optoelectronic devices and applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Transport of High-energy Charged Particles through Spatially Intermittent Turbulent Magnetic Fields

Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

High-order diamond differencing schemes for the Boltzmann Fokker-Planck equation in 3D Cartesian geometries

The Boltzmann Fokker-Planck, an approximate form of the linear Boltzmann equation is commonly used to treat efficiently the transport of charged particles in matter. This paper introduces the application of high-order diamond differencing schemes (HODD), specifically the DD1 and DD2 schemes which are 4- and 6-order accurate respectively, to handle the spatial discretization of that equation in 3D Cartesian geometries. The energy deposition solutions for the coupled transport of electrons and photons presented in this work shows that HODD, compared to classical DD scheme, provides correction to the oscillations and a reduced propensity to yield negative fluxes. They are useful tools to minimize local error, notably in regions with abrupt variations of the flux solution. They also can be used to reduced execution time by decreasing the needed number of voxels to obtain a fixed accuracy. On the tested benchmarks, the DD1 scheme is 87%- 92%-91% more accurate than the classical DD scheme for total, mean per-voxels and maximum deviation of energy deposition values respectively. For comparison, a calculation with 8 times more voxels, requiring roughly 2.5 times more time to execute, is 92%-90%-77% more accurate. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Resonant Slow Extraction Simulation using Bmad

Simulations of slow extraction and transport of charged particle beams from synchrotrons requires careful modeling and experimentation. In this paper, we outline how the Bmad modeling software was adapted and developed to run third-integer resonant extraction simulations of beams at Brookhaven National Laboratory’s Booster synchrotron. Further, we show experimental comparisons of the slowly extracted beams transferred to the NASA Space Radiation Laboratory (NSRL) transport line. In this process, beam passes through a stripping foil element at the extraction point, which, along with stripping remaining electrons from the beam ions, acts as a scatterer to modify the phase space, moving to a more Gaussian-like distribution. This modification to the beam helps generate a uniform beam at the beam line’s target location; which is necessary for the variety of experiments performed at the facility. During this work, beam energy loss and multiple scattering by foil routines were built in conjunction with the Bmad code developers and are now integrated into the software.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

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