Search NASA⌕ Search

Engineering topics

Correa, A. A.

Publications and source records attributed to Correa, A. A..

Review of the first charged-particle transport coefficient comparison workshop

We present the results of the first Charged-Particle Transport Coefficient Code Comparison Workshop, which was held in Albuquerque, NM October 4–6, 2016. In this first workshop, scientists from eight institutions and four countries gathered to compare calculations of transport coefficients including thermal and electrical conduction, electron–ion coupling, inter-ion diffusion, ion viscosity, and charged particle stopping powers. In this paper, we give general background on Coulomb coupling and computational expense, review where some transport coefficients appear in hydrodynamic equations, and present the submitted data. Large variations are found when either the relevant Coulomb coupling parameter is large or computational expense causes difficulties. Understanding the general accuracy and uncertainty associated with such transport coefficients is important for quantifying errors in hydrodynamic simulations of inertial confinement fusion and high-energy density experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reduction of electron-phonon coupling in warm dense iron

The electron-ion relaxation dynamics in warm dense iron were investigated by time-resolved x-ray absorption near-edge spectroscopy (XANES). A novel analysis, combining ab initio density functional theory (DFT) and two temperature model (TTM) simulations, was developed to calculate the x-ray absorption spectra as a function of delay time. Here we present experimental evidence of changes at the XAS L3 edge of iron that are consistent with the reduction of the electron-phonon coupling under warm dense matter conditions. The experimental results are in agreement with the model that takes both the electron ($T_e$) and the ion temperature ($T_i$) dependence of the thermophysical properties into consideration, while models where either constant electron-phonon coupling factor (G) or only $T_e$- dependent G are used do not agree with the observed relaxation dynamics of iron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Role of electrons in collision cascades in solids. I. Dissipative model

In this article, we present a detailed model for the nonadiabatic coupling between ions and electrons in energetic ion-solid interactions over a wide range of energies in concentrated solid-solution fcc alloys of the 3d transition metals Ni, Co, Fe, and Cr. The model is based on general statistical mechanical principles and results in a stochastic modification of the classical nuclei motion which is parameterized by the first-principles calculation of a dissipation function produced by explicit time-dependent electronic evolution. This model provides a full picture of an entire collision process, from the ballistic to the thermal phases of a cascade, giving a detailed description of the energy exchange between ions and electrons till their final thermalization, removing in this way some ad hoc assumptions used in the state-of-the-art atomistic two-temperature models. Here, this work is separated in two papers; in the present Part I, we report on the ab initio methodology used to translate stopping power and electron-phonon interaction into a parameterized dissipation function; Part II, to be published, addresses the nonadiabatic ion dynamics using our modified Langevin dynamics [Tamm et al. Phys. Rev. Lett. 120, 185501 (2018)] applying the dissipation functions developed here to specific collision cascade events.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Role of electrons in collision cascades in solids. II. Molecular dynamics

In this work, we present a model for the role of electrons in collision cascades in solids in which the coupling between ions and electrons is calculated using first-principles models and introduced into the classical ion equations of motion using our modified Langevin dynamics [A. Tamm et al., Phys. Rev. Lett. 120, 185501 (2018)]. This model gives a full picture of the entire collision process, from the ballistic to the thermal phases of a cascade, giving a detailed representation of the energy exchange between ions and electrons until their final thermalization, removing in this way some ad hoc assumptions used in the state-of-the-art two-temperature model. This work is separated into two papers: Part I [M. Caro et al., Phys. Rev. B 99, 174301 (2019)] reports on the ab initio methodology used to translate stopping power into the parametrized dissipation function. Finally, part II applies the nonadiabatic ion dynamics using the dissipation functions developed in Part I to specific collision cascade events.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ab initio modeling of nonequilibrium electron-ion dynamics of iron in the warm dense matter regime

Here, the spatiotemporal electron and ion relaxation dynamics of iron induced by femtosecond laser pulses was studied using a one-dimensional two-temperature model (1D-TTM) where electron and ion temperature-dependent thermophysical parameters such as specific heat ( C ), electron-phonon coupling ( G ), and thermal conductivity ( K ) were calculated with ab initio density-functional-theory (DFT) simulations. Based on the simulated time evolutions of electron and ion temperature distributions [ T e ( x,t ) and T i ( x,t )], the time evolution of x-ray absorption near-edge spectroscopy (XANES) was calculated and compared with experimental results reported by Fernandez-Pañella et al., where the slope of XANES spectrum at the onset of absorption (s) was used due to its excellent sensitivity to the electron temperature. Our results indicate that the ion temperature dependence on G and C , which is largely neglected in the past studies, is very important for studying the nonequilibrium electron-ion relaxation dynamics of iron in warm dense matter (WDM) conditions. It is also shown that the 1/s behavior becomes very sensitive to the thermal gradient profile, in other words, to the values of K in a TTM simulation, for target thickness of about two to four times the mean free path of conduction electrons. Our approach based on 1D-TTM and XANES simulations can be used to determine the optimal combination of target geometry and laser fluence for a given target material, which will enable us to tightly constrain the thermophysical parameters under electron-ion nonequilibrium WDM conditions.

36 MATERIALS SCIENCE↗

On the local density dependence of electronic stopping of ions in solids

We use here time-dependent density functional theory to calculate the electronic stopping S e in binary Ni-Ni atomic collisions and for a Ni projectile along channeling directions in a Ni crystal. Our results show that when S e is reported as a function of the ground state target electronic density ρ 0 the stopping is not a single-valued function of the local density, as assumed in formalisms that date back to the origins of quantum mechanics, but shows loops, suggesting that it is inaccurate to model stopping as a dissipative force of the type F = β ( ρ 0 ) v , as it is customarily done in non-adiabatic molecular dynamics simulations of ion-solid interactions. We compare our results with S e in a uniform electron gas where the above definition for the force holds, and conclude on the validity of using jellium as a crude approximation for more realistic inhomogeneous electron gases.

74 ATOMIC AND MOLECULAR PHYSICS↗

Measuring the shock impedance mismatch between high-density carbon and deuterium at the National Ignition Facility

Fine-grained diamond, or high-density carbon (HDC), is being used as an ablator for inertial confinement fusion (ICF) research at the National Ignition Facility (NIF). Accurate equation of state (EOS) knowledge over a wide range of phase space is critical in the design and analysis of integrated ICF experiments. Here in this paper, we report shock and release measurements of the shock impedance mismatch between HDC and liquid deuterium conducted during shock-timing experiments having a first shock in the ablator ranging between 8 and 14 Mbar. Using ultrafast Doppler imaging velocimetry to track the leading shock front, we characterize the shock velocity discontinuity upon the arrival of the shock at the HDC/liquid deuterium interface. Comparing the experimental data with tabular EOS models used to simulate integrated ICF experiments indicates the need for an improved multiphase EOS model for HDC in order to achieve a significant increase in neutron yield in indirect-driven ICF implosions with HDC ablators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stopping power beyond the adiabatic approximation

Energetic ions traveling in solids deposit energy in a variety of ways, being nuclear and electronic stopping the two avenues in which dissipation is usually treated. This separation between electrons and ions relies on the adiabatic approximation in which ions interact via forces derived from the instantaneous electronic ground state. In a more detailed view, in which non-adiabatic effects are explicitly considered, electronic excitations alter the atomic bonding, which translates into changes in the interatomic forces. In this work, we use time dependent density functional theory and forces derived from the equations of Ehrenfest dynamics that depend instantaneously on the time-dependent electronic density. With them we analyze how the inter-ionic forces are affected by electronic excitations in a model of a Ni projectile interacting with a Ni target, a metallic system with strong electronic stopping and shallow core level states. We find that the electronic excitations induce substantial modifications to the inter-ionic forces, which translate into nuclear stopping power well above the adiabatic prediction. Particularly, we observe that most of the alteration of the adiabatic potential in early times comes from the ionization of the core levels of the target ions, not readily screened by the valence electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electron-phonon interaction within classical molecular dynamics

Here, we present a model for nonadiabatic classical molecular dynamics simulations that captures with high accuracy the wave-vector q dependence of the phonon lifetimes, in agreement with quantum mechanics calculations. It is based on a local view of the e-ph interaction where individual atom dynamics couples to electrons via a damping term that is obtained as the low-velocity limit of the stopping power of a moving ion in a host. The model is parameter free, as its components are derived from ab initio -type calculations, is readily extended to the case of alloys, and is adequate for large-scale molecular dynamics computer simulations. We also show how this model removes some oversimplifications of the traditional ionic damped dynamics commonly used to describe situations beyond the Born-Oppenheimer approximation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electron elevator: Excitations across the band gap via a dynamical gap state

We use time-dependent density functional theory to study self-irradiated Si. We calculate the electronic stopping power of Si in Si by evaluating the energy transferred to the electrons per unit path length by an ion of kinetic energy from 1 eV to 100 keV moving through the host. Electronic stopping is found to be significant below the threshold velocity normally identified with transitions across the band gap. A structured crossover at low velocity exists in place of a hard threshold. Lastly, an analysis of the time dependence of the transition rates using coupled linear rate equations enables one of the excitation mechanisms to be clearly identified: a defect state induced in the gap by the moving ion acts like an elevator and carries electrons across the band gap.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗