Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stopping Power”

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 19 records

Quantum computation of stopping power for inertial fusion target design

Stopping power is the rate at which a material absorbs the kinetic energy of a charged particle passing through it—one of many properties needed over a wide range of thermodynamic conditions in modeling inertial fusion implosions. First-principles stopping calculations are classically challenging because they involve the dynamics of large electronic systems far from equilibrium, with accuracies that are particularly difficult to constrain and assess in the warm-dense conditions preceding ignition. Here, we describe a protocol for using a fault-tolerant quantum computer to calculate stopping power from a first-quantized representation of the electrons and projectile. Our approach builds upon the electronic structure block encodings of Su et al. [ PRX Quant. 2 , 040332 (2021)], adapting and optimizing those algorithms to estimate observables of interest from the non-Born–Oppenheimer dynamics of multiple particle species at finite temperature. We also work out the constant factors associated with an implementation of a high-order Trotter approach to simulating a grid representation of these systems. Ultimately, we report logical qubit requirements and leading-order Toffoli costs for computing the stopping power of various projectile/target combinations relevant to interpreting and designing inertial fusion experiments. We estimate that scientifically interesting and classically intractable stopping power calculations can be quantum simulated with roughly the same number of logical qubits and about one hundred times more Toffoli gates than is required for state-of-the-art quantum simulations of industrially relevant molecules such as FeMoco or P450.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Machine learning and TDDFT software for stopping power computation

(SF-24-012) Stopping power describes the rate that a material slows radiation particles passing through it and is useful in designing many technologies. Few organizations can perform new measurements, which require significant resources and rare equipment, and all others rely on coarse approximations rendered from pre-existing data. Methods for computing stopping power in new materials, such as time-dependent density functional theory (TD-DFT), have only recently (circa-2015) become available but are too computationally costly to use frequently enough to have a pronounced impact. We have created a method that opens a pathway to computing stopping power without any need for experimental data by combining electronic structure computations and machine learning.

Ward, Logan↗

Calculating electronic stopping power in materials from first principles

From the early models of electronic stopping power to the current first principles simulations, the techniques evolved to increase the range of validity and to reduce empiricism. Thanks to a combination of theoretical advances provided by Time Dependent Density Functional Theory and the development of numerical codes, it is currently possible to predict electronic stopping power for realistic materials by performing direct simulations of the electron excitation processes beyond linear response, and including electronic band structure effects. Electronic stopping power is an important quantity used to predict and understand the effects of particle radiation in matter. First principles calculations of electronic stopping power can be applied to any atomistic system, solids, liquids and alloys. This review aims here to help graduate level students and researchers immerse themselves into state-of-the-art techniques to computationally model and calculate electronic stopping power.

36 MATERIALS SCIENCE↗

Ab Initio Studies on the Stopping Power of Warm Dense Matter with Time-Dependent Orbital-Free Density Functional Theory

Here, electronic transport properties of warm dense matter, such as electrical/thermal conductivities and nonadiabatic stopping power, are of particular interest to geophysics, planetary science, astrophysics, and inertial confinement fusion (ICF). One example is the α-particle stopping power of dense deuterium–tritium (DT) plasmas, which must be precisely known for current small-margin ICF target designs to ignite. We have developed a time-dependent orbital-free density functional theory (TD-OF-DFT) method for ab initio investigations of the charged-particle stopping power of warm dense matter. Our current dependent TD-OF-DFT calculations have reproduced the recently well-characterized stopping power experiment in warm dense beryllium. Forα-particle stopping in warm and solid-density DT plasmas, the ab initio TD-OF-DFT simulations show a lower stopping power up to ~25% in comparison with two stopping-power models widely used in the high-energy-density physics community.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electronic stopping power of protons and alpha particles in nickel

The electronic stopping power of nickel for protons and alpha particles at velocities below and around the Fermi velocity has been obtained to high accuracy using time-dependent density functional theory. For the wide range of projectile velocities considered, we observed di erent regimes of electronic stopping due to the alternative participation of s- and d-band electrons. Despite the sharp discontinuity in the electronic density of states near the Fermi energy characteristic of the nickel band structure, we do not nd an anomalous non-linear electronic stopping power limit as a function of velocity. However, we nd a crossover region above v = 0:15 a:u: both for protons and alpha particles, related to the increase in participating host electrons and, in the case of alpha particles, to an increase of the charge state. We compare our calculated results with widely available experimental data and analyze the low velocity limits in the context of Lindhard's linear response theory and previous non-linear density functional calculations. The comparison shows good accord with the lowest velocity experiments available. This may indicate that the adiabatic local density approximation is already a good theory to calculate electronic stopping power in materials at low velocity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sensitivity of an (α,n) neutron source to (α,n) cross sections and stopping powers

The first-order sensitivities of the (α,n) neutron source rate density in a homogeneous material to isotopic (α,n) cross sections, elemental electronic stopping power data, and general coefficients for the nuclear stopping power are derived for the total neutron source and the energy-dependent neutron source. The physics bases for some of the equations in the SOURCES4C calculation of the stopping power are briefly discussed. The derivatives have been implemented in SOURCES4C. In simple PuO 2 and PuB 2 test problems, the computed derivatives compare extremely well with central differences. As a result, the first-order sensitivities are applied to show that improving the value of one of SOURCES4C’s stopping-power coefficients changes the (α,n) neutron source rate by 1.738 × 10 –6 % in the PuO 2 problem.

(α,n) neutron source↗

Validation of the single-event method for low-energy electron transport via stopping power calculations with $\mathrm{MCNP}$

Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Time-dependent orbital-free density functional theory for electronic stopping power: Comparison to the Mermin-Kohn-Sham theory at high temperatures

Electronic stopping power in warm dense matter can affect energy transport and heating in astrophysical processes and internal confinement fusion. For cold condensed matter systems, stopping power can be modeled from first-principles using real-time time-dependent density functional theory (DFT). However, high temperatures (10's to 100's of eV) may be computationally prohibitive for traditional Mermin-Kohn-Sham DFT. New experimental measurements in the warm dense regime motivates the development of first-principles approaches, which can reach these temperatures. Here, we have developed a time-dependent orbital-free density functional theory, which includes a novel nonadiabatic and temperature-dependent kinetic energy density functional, for the simulation of stopping power at any temperature. The approach is nonlinear with respect to the projectile perturbation, includes all ions and electrons, and does not require a priori determination of screened interaction potentials. Finally, our results compare favorably with Kohn-Sham for temperatures in the WDM regime, especially nearing 100 eV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reproducibility of real-time time-dependent density functional theory calculations of electronic stopping power in warm dense matter

Real-time time-dependent density functional theory (TDDFT) is widely considered to be the most accurate available method for calculating electronic stopping powers from first principles, but there have been relatively few assessments of the consistency of its predictions across different implementations. This problem is particularly acute in the warm dense regime, where computational costs are high and experimental validation is rare and resource intensive. We report a comprehensive cross-verification of stopping power calculations in conditions relevant to inertial confinement fusion conducted using four different TDDFT implementations. We find excellent agreement among both the post-processed stopping powers and relevant time-resolved quantities for alpha particles in warm dense hydrogen. We also analyze sensitivities to a wide range of methodological details, including the exchange-correlation model, pseudopotentials, initial conditions, observable from which the stopping power is extracted, averaging procedures, projectile trajectory, and finite-size effects. We show that among these details, pseudopotentials, trajectory-dependence, and finite-size effects have the strongest influence, and we discuss different strategies for controlling the latter two considerations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of Coulomb coupling on stopping power and a link to macroscopic transport

Molecular dynamics simulations are used to assess the influence of Coulomb coupling on the energy evolution of charged projectiles in the classical one-component plasma. The average projectile kinetic energy is found to decrease linearly with time when ν α /ω p ≲ 10 –2 , where να is the Coulomb collision frequency between the projectile and the medium, and ω p is the plasma frequency. Stopping power is obtained from the slope of this curve. In comparison to the weakly coupled limit, strong Coulomb coupling causes the magnitude of the dimensionless stopping power, ( a/k B T ) dE / dx , to increase, the Bragg peak to shift to several times the plasma thermal speed, and for the stopping power curve to broaden substantially. The rate of change of the total projectile kinetic energy averaged over many independent simulations is shown to consist of two measurable components: a component associated with a one-dimensional friction force and a thermal energy exchange rate. In the limit of a slow and massive projectile, these can be related to the macroscopic transport rates of self-diffusion and temperature relaxation in the plasma. Simulation results are compared with available theoretical models for stopping power, self-diffusion, and temperature relaxation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measuring stopping power in warm dense matter plasmas at OMEGA

A platform has been developed for accurately measuring the stopping power of high energy protons through warm dense matter (WDM) plasmas characterized by x-ray Thomson scattering. Here, in this work, stopping power measurements were successfully made through both WDM Beryllium and Boron plasmas. In the Boron experiments, an increase in stopping was observed over their cold target counter-parts. This increase in stopping was shown to agree well with models that account for the partial ionization of the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Trajectory sampling and finite-size effects in first-principles stopping power calculations

Abstract Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for evaluating and optimizing trajectories in this context. This methodology enables rigorous analysis of the failure modes of various common trajectory choices in crystalline materials. Although randomly selecting trajectories is common practice in stopping power calculations in solids, we show that nearly 30% of random trajectories in an FCC aluminum crystal will not representatively sample the material over the time and length scales feasibly simulated with TDDFT, and unrepresentative choices incur errors of up to 60%. We also show that finite-size effects depend on ion trajectory via “ouroboros” effects beyond the prevailing plasmon-based interpretation, and we propose a cost-reducing scheme to obtain converged results even when expensive core-electron contributions preclude large supercells. This work helps to mitigate poorly controlled approximations in first-principles stopping power calculations, allowing 1–2 order of magnitude cost reductions for obtaining representatively averaged and converged results.

36 MATERIALS SCIENCE↗

Stopping-power enhancement from discrete particle-wake correlations in high-energy-density plasmas

Three-dimensional (3D) simulations of electron beams propagating in high-energy-density plasmas using the quasistatic Particle-in-Cell (PIC) code QuickPIC demonstrate a significant increase in stopping power when beam electrons mutually interact via their wakes. Each beam electron excites a plasma wave wake of wavelength ~ 2 π c / ω p e , where c is the speed of light and ω p e is the background plasma frequency. Here we show that a discrete collection of electrons undergoes a beam-plasma-like instability caused by mutual particle-wake interactions that causes electrons to bunch in the beam, even for beam densities n b for which fluid theory breaks down. This bunching enhances the beam's stopping power, which we call “correlated stopping,” and the effect increases with the “correlation number” N b ≡ n b ( c / ω p e ) 3 . For example, a beam of monoenergetic 9.7 MeV electrons with N b = 1 / 8 , in a cold background plasma with n e = 10 26 cm – 3 (450 g cm – 3 DT), has a stopping power of 2.28 ± 0.04 times the single-electron value, which increases to 1220 ± 5 for N b = 64 . The beam also experiences transverse filamentation, which eventually limits the stopping enhancement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Positron-impact electronic excitations and mass stopping power of H 2

Positron-impact electronic excitation cross sections, mean excitation energies, and mass stopping power of the H 2 molecule have been calculated for energies from 10 eV up to 2 keV using the convergent close-coupling method that utilizes single- and two-center expansions. Overall, results are compared to previous studies. Application of Bragg's rule of stopping power additivity is discussed by comparing results obtained for atomic (H) and molecular (H 2 ) targets for positron impact.

74 ATOMIC AND MOLECULAR PHYSICS↗

Plasma stopping-power measurements reveal transition from non-degenerate to degenerate plasmas

Physically realized electron gas systems usually reside in either the quantum non-degenerate or fully degenerate limit, where the average de Broglie wavelength of the thermal electrons becomes comparable with the interparticle distance between electrons. A few systems, such as young brown dwarfs and the cold dense fuels created in imploded cryogenic capsules at the National Ignition Facility, lie between these two limits and are partially degenerate. The National Ignition Facility has the unique capability of varying the electron quantum degeneracy by adjusting the laser drive used to implode the capsules. This allows experimental studies of the effects of the degeneracy level on plasma transport properties. By measuring rare nuclear reactions in these cold dense fuels, we show that the electron stopping power, which is the rate of energy loss per unit distance travelled by a charged particle, changes with increasing electron density. We observe a quantum-induced shift in the peak of the stopping power using diagnostics that measure above and below this peak. The observed changes in the stopping power are shown to be unique to the transition region between non-degenerate and degenerate plasmas. Furthermore, our results support the screening models applied to partially degenerate astrophysical systems such as young brown dwarfs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗

Electron mass stopping power in H 2

Calculations of electron mass stopping power (SP) of electrons in ${\mathrm{H}}_{2}$ have been performed here using the convergent close-coupling method for incident electron energies up to 2000 eV. Convergence of the calculated SP has been established by increasing the size of the close-coupling expansion from 9 to 491 states. Good agreement was found with the SP measurements of Munoz et al. [Chem. Phys. Lett. 433, 253 (2007)].

74 ATOMIC AND MOLECULAR PHYSICS↗

Mixed stochastic-deterministic time-dependent density functional theory: application to stopping power of warm dense carbon

Warm dense matter (WDM) describes an intermediate phase, between condensed matter and classical plasmas, found in natural and man-made systems. In a laboratory setting, WDM is often created dynamically. It is typically laser or pulse-power generated and can be difficult to characterize experimentally. Measuring the energy loss of high energy ions, caused by a WDM target, is both a promising diagnostic and of fundamental importance to inertial confinement fusion research. However, electron coupling, degeneracy, and quantum effects limit the accuracy of easily calculable kinetic models for stopping power, while high temperatures make the traditional tools of condensed matter, e.g. time-dependent density functional theory (TD-DFT), often intractable. In this study, we have developed a mixed stochastic-deterministic approach to TD-DFT which provides more efficient computation while maintaining the required precision for model discrimination. Recently, this approach showed significant improvement compared to models when compared to experimental energy loss measurements in WDM carbon. Here, we describe this approach and demonstrate its application to warm dense carbon stopping across a range of projectile velocities. We compare direct stopping-power calculation to approaches based on combining homogeneous electron gas response with bound electrons, with parameters extracted from our TD-DFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗