Molecular dynamics simulations of classical stopping power

Electron-proton relaxation in hot-dense plasmas with a screened quantum statistical potential

Modeling the nonequilibrium process between ions and electrons is of great importance in laboratory fusion ignition, laser-plasma interaction, and astrophysics. For hot and dense plasmas, theoretical descriptions of Coulomb collisions remain complicated due to quantum effect at short distances and screening effect at long distances. In this paper, we propose an analytical screened quantum statistical potential that takes into account both the short-range quantum diffraction effect and the long-range screening effect. By implementing the newly developed potential into the binary scattering framework, the electron-proton temperature relaxation in hot-dense hydrogen plasmas is investigated. In both the classical and quantum limits, analytical expressions for the Coulomb logarithm have been obtained, which are generally embedded in an asymptotic matching formula. Quantitative comparisons with molecular dynamics simulations and recent OMEGA experiments demonstrate that the present modeling is well suited to describe the temperature relaxation process between electrons and ions in hot-dense plasmas.

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.

Simulation study of cone-in-shell target for indirect-drive ion fast ignition concept under the theory of an effective interaction potential

The stopping power of charged particles released by the deuterium-tritium nuclear reactions has been extensively studied in the weakly to moderately coupled plasma regimes. We have modified the conventional effective potential theory (EPT) stopping framework to have a practical connection to investigate the ions energy loss characteristics in fusion plasma. Our modified EPT model differs from the original EPT framework by a coefficient of order [Formula: see text]([Formula: see text] is a velocity-dependent generalization of the Coulomb logarithm). Molecular dynamics simulations agree well with our modified stopping framework. To study the role of related stopping formalisms in ion fast ignition, we simulate the cone-in-shell configuration under laser-accelerated aluminum beam incidence. In ignition/burn phase, the performance of our modified model is in agreement with its original form and the conventional Li-Petrasso (LP) and Brown-Preston-Singleton (BPS) theories. The LP theory indicates the fastest rate in providing ignition/burn condition. Our modified EPT model with a discrepancy of [Formula: see text] 9%, has the most agreement with LP theory, while that of the original EPT (with a discrepancy of [Formula: see text] 47% to LP) and BPS (with a discrepancy of [Formula: see text] 48% to LP) methods maintain the third and fourth contributions in accelerating the ignition time, respectively.

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. 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.

Temperature relaxation in strongly-coupled binary ionic mixtures

New facilities such as the National Ignition Facility and the Linac Coherent Light Source have pushed the frontiers of high energy-density matter. These facilities offer unprecedented opportunities for exploring extreme states of matter, ranging from cryogenic solid-state systems to hot, dense plasmas, with applications to inertial-confinement fusion and astrophysics. However, significant gaps in our understanding of material properties in these rapidly evolving systems still persist. In particular, non-equilibrium transport properties of strongly-coupled Coulomb systems remain an open question. Here, we study ion-ion temperature relaxation in a binary mixture, exploiting a recently-developed dual-species ultracold neutral plasma. We compare measured relaxation rates with atomistic simulations and a range of popular theories. Our work validates the assumptions and capabilities of the simulations and invalidates theoretical models in this regime. This work illustrates an approach for precision determinations of detailed material properties in Coulomb mixtures across a wide range of conditions.

Molecular dynamics investigation of the stopping power of warm dense hydrogen for electrons

A variety of theoretical models have been proposed to calculate the stopping power of charged particles in matter, which is a fundamental issue in many fields. However, the approximation adopted in these theories will be challenged under warm dense matter conditions. Molecular dynamics (MD) simulation is a good way to validate the effectiveness of these models. We investigate the stopping power of warm dense hydrogen for electrons with projectile energies ranging from 400-10000 eV by means of an electron force field (eFF) method, which can effectively avoid the Coulomb catastrophe in conventional MD calculations. It is found that the stopping power of warm dense hydrogen decreases with increasing temperature of the sample at those high projectile velocities. This phenomenon could be explained by the effect of electronic structure dominated by bound electrons, which is further explicated by a modified random phase approximation (RPA) model based on local density approximation proper to inhomogeneous media. Most of the models extensively accepted by the plasma community, e.g., Landau-Spitzer model, Brown-Preston-Singleton model and RPA model, cannot well address the effect caused by bound electrons so that their predictions of stopping power contradict our result. Therefore, the eFF simulations of this paper reveals the important role played by the bound electrons on stopping power in warm dense plasmas.

Stopping power of hot dense deuterium-tritium plasmas mixed with impurities to charged particles

In this work, we studied the stopping power of deuterium-tritium (DT) plasmas mixed with impurities to the injected charged particles. Based on the Brown-Preston-Singleton model, the analytical expression for the change ratio of stopping power (denoted by η) induced by impurities in DT plasmas is developed, in which both classical short-distance collision part and quantum correction contribution are purely linear response to the impurity concentration ξ_{X}, while the classical long-range collision brings about higher-order nonlinear response to ξ_{X}. Furthermore, the expression for change ratio of deposition depth (denoted by χ) of charged particles induced by impurities in DT plasmas is also derived. As applications, we systemically investigated the energy loss of α particles deposited into a hot dense DT plasma mixed with impurity X(X=C, Si, Ge), where the temperature and density of DT are smaller than 10 keV and 500 g/cm^{3} and the concentration of Xξ_{X} is less than 5%. The numerical results suggest that (i) for the case of C mixed into DT, both change ratios of stopping power and deposition depth of α particles (i.e., η and χ) are linear response to the concentration of C ξ_{C}; (ii) for the case of Si mixed into DT, the second-order nonlinear response of η and χ to ξ_{Si} cannot be ignored when the densities of DT are larger than 200 g/cm^{3}; and (iii) for the case of Ge mixed into DT, the second- and third-order nonlinear response of η and χ to ξ_{Ge} are very remarkable because of the higher ionization degree and heavier atomic mass of Ge. The formulas and findings in this work may be helpful to the research of internal confinement fusion (ICF) related implosion physics and may provide useful theoretical guidance and data for the design of ICF target.

Ion energy-loss characteristics and friction in a free-electron gas at warm dense matter and nonideal dense plasma conditions

We investigate the energy-loss characteristics of an ion in warm dense matter (WDM) and dense plasmas concentrating on the influence of electronic correlations. The basis for our analysis is a recently developed ab initio quantum Monte Carlo- (QMC) based machine learning representation of the static local field correction (LFC) [Dornheim et al., J. Chem. Phys. 151, 194104 (2019)JCPSA60021-960610.1063/1.5123013], which provides an accurate description of the dynamical density response function of the electron gas at the considered parameters. We focus on the polarization-induced stopping power due to free electrons, the friction function, and the straggling rate. In addition, we compute the friction coefficient which constitutes a key quantity for the adequate Langevin dynamics simulation of ions. Considering typical experimental WDM parameters with partially degenerate electrons, we find that the friction coefficient is of the order of γ/ω_{pi}=0.01, where ω_{pi} is the ionic plasma frequency. This analysis is performed by comparing QMC-based data to results from the random-phase approximation (RPA), the Mermin dielectric function, and the Singwi-Tosi-Land-Sjölander (STLS) approximation. It is revealed that the widely used relaxation time approximation (Mermin dielectric function) has severe limitations regarding the description of the energy loss of ions in a correlated partially degenerate electrons gas. Moreover, by comparing QMC-based data with the results obtained using STLS, we find that the ion energy-loss properties are not sensitive to the inaccuracy of the static local field correction (LFC) at large wave numbers, k/k_{F}>2 (with k_{F} being the Fermi wave number), but that a correct description of the static LFC at k/k_{F}≲1.5 is important.

Transport of intense particle beams in large-scale plasmas

Transport of particle beams in plasmas is widely employed in fundamental research, industry, and medicine. Due to the high inertia of ion beams, their transport in plasmas is usually assumed to be stable. Here we report the focusing and flapping of intense slab proton beams transporting through large-scale plasmas by using a recently developed kinetic particle-in-cell simulation code. The beam self-focusing effect in the simulation is prominent and agrees well with previous experiments and theories. Moreover, the beam can curve and flap like turbulence as the beam density increases. Simulation and analysis indicate that the self-generated magnetic fields, produced by movement of collisional plasmas, are the dominant driver of such behaviors. By analyzing the spatial growth rate of magnetic energy and energy deposition of injected proton beams, it is found that the focusing and flapping are significantly determined by the injected beam densities and energies. In addition, a remarkable nonlinear beam energy loss is observed. Our research might find application in inertial confinement fusion and also might be of interest to the laboratory astrophysics community.

Experimental Validation of Low-Z Ion-Stopping Formalisms around the Bragg Peak in High-Energy-Density Plasmas

We report on the first accurate validation of low-Z ion-stopping formalisms in the regime ranging from low-velocity ion stopping-through the Bragg peak-to high-velocity ion stopping in well-characterized high-energy-density plasmas. These measurements were executed at electron temperatures and number densities in the range of 1.4-2.8 keV and 4×10^{23}-8×10^{23}  cm^{-3}, respectively. For these conditions, it is experimentally demonstrated that the Brown-Preston-Singleton formalism provides a better description of the ion stopping than other formalisms around the Bragg peak, except for the ion stopping at v_{i}∼0.3v_{th}, where the Brown-Preston-Singleton formalism significantly underpredicts the observation. It is postulated that the inclusion of nuclear-elastic scattering, and possibly coupled modes of the plasma ions, in the modeling of the ion-ion interaction may explain the discrepancy of ∼20% at this velocity, which would have an impact on our understanding of the alpha energy deposition and heating of the fuel ions, and thus reduce the ignition threshold in an ignition experiment.

Anomalous material-dependent transport of focused, laser-driven proton beams

Intense lasers can accelerate protons in sufficient numbers and energy that the resulting beam can heat materials to exotic warm (10 s of eV temperature) states. Here we show with experimental data that a laser-driven proton beam focused onto a target heated it in a localized spot with size strongly dependent upon material and as small as 35 μm radius. Simulations indicate that cold stopping power values cannot model the intense proton beam transport in solid targets well enough to match the large differences observed. In the experiment a 74 J, 670 fs laser drove a focusing proton beam that transported through different thicknesses of solid Mylar, Al, Cu or Au, eventually heating a rear, thin, Au witness layer. The XUV emission seen from the rear of the Au indicated a clear dependence of proton beam transport upon atomic number, Z, of the transport layer: a larger and brighter emission spot was measured after proton transport through the lower Z foils even with equal mass density for supposed equivalent proton stopping range. Beam transport dynamics pertaining to the observed heated spot were investigated numerically with a particle-in-cell (PIC) code. In simulations protons moving through an Al transport layer result in higher Au temperature responsible for higher Au radiant emittance compared to a Cu transport case. The inferred finding that proton stopping varies with temperature in different materials, considerably changing the beam heating profile, can guide applications seeking to controllably heat targets with intense proton beams.

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

Electronic transport properties of warm dense matter, such as electrical or 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 three stopping-power models often used in the high-energy-density physics community.

Observation of a strong-coupling effect on electron-ion collisions in ultracold plasmas

Ultracold plasmas (UCPs) provide a well-controlled system for studying multiple aspects in plasma physics that include collisions and strong-coupling effects. By applying a short electric field pulse to an UCP, a plasma electron center-of-mass oscillation can be initiated. For accessible parameter ranges, the damping rate of this oscillation is determined by the electron-ion collision rate. We performed measurements of the oscillation damping rate with such parameters and compared the measured rates to both a molecular dynamics (MD) simulation that includes strong-coupling effects and a Monte Carlo binary collision simulation designed to predict the damping rate including only weak-coupling considerations. We found agreement between the experimentally measured damping rate and the MD result. This agreement did require including the influence of a previously unreported UCP heating mechanism whereby the presence of a dc electric field during ionization increased the electron temperature, but estimations and simulations indicate that such a heating mechanism should be present for our parameters. The measured damping rate at our coldest electron temperature conditions was much faster than the weak-coupling prediction obtained from the Monte Carlo binary collision simulation, which indicates the presence of a significant strong-coupling influence. The density averaged electron strong-coupling parameter Γ measured at our coldest electron temperature conditions was 0.35(8).

Measurement of Preheat Due to Nonlocal Electron Transport in Warm Dense Matter

This Letter presents a novel approach to study electron transport in warm dense matter. It also includes the first x-ray Thomson scattering (XRTS) measurement from low-density CH foams compressed by a strong laser-driven shock at the OMEGA laser facility. The XRTS measurement is combined with velocity interferometry (VISAR) and optical pyrometry (SOP) providing a robust measurement of thermodynamic conditions in the shock. Evidence of significant preheat contributing to elevated temperatures reaching 17.5-35 eV in shocked CH foam is measured by XRTS. These measurements are complemented by abnormally high shock velocities observed by VISAR and early emission seen by SOP. These results are compared to radiation hydrodynamics simulations that include first-principles treatment of nonlocal electron transport in warm dense matter with excellent agreement. Additional simulations confirm that the x-ray contribution to this preheat is negligible.

Microfield dynamics in dense hydrogen plasmas with high-Z impurities

We use large-scale classical molecular dynamics to determine microfield properties for several dense plasma mixtures. By employing quantum statistical potentials (QSPs) to regularize the Coulomb interaction, our simulations follow motions of electrons as well as ions for times long enough to track relaxation phenomena involving both types of particles. Coulomb coupling, relative to temperature, of different pairs of species in the hot, dense matter being simulated ranges from weak to strong. We first study the effect of such coupling differences, along with composition and QSP differences, on the roles of electrons and various mixture components in determining probability distributions of instantaneous, total microfields experienced by the ions. Then, we address two important dynamical questions: (1) How is the quasistatic part of the total field to be extracted from the time-dependent simulation data? (2) Under what conditions does the commonly used approximation of ions with fixed Yukawa-like screening by free electrons accurately describe quasistatic fields? We identify a running, short-time average of the total field at each ion as its slowly evolving, quasistatic part. We consider several ways to specify the averaging interval, and note the influence of ion dynamics in this issue. When all species are weakly coupled, the quasistatic fields have probability distributions agreeing well with those we obtain from simulations of Yukawa-screened ions. However, agreement deteriorates as the coupling between high-Z ions increases well beyond unity, principally because the Yukawa model tends to underestimate the true screening of close high-Z pairs. Examples of this fact are given, and some consequences for the high-field portions of probability distributions are discussed.

Dynamic properties of the energy loss of multi-MeV charged particles traveling in two-component warm dense plasmas

The energy loss of multi-MeV charged particles moving in two-component warm dense plasmas (WDPs) is studied theoretically beyond the random-phase approximation. The short-range correlations between particles are taken into account via dynamic local field corrections (DLFC) in a Mermin dielectric function for two-component plasmas. The mean ionization states are obtained by employing the detailed configuration accounting model. The Yukawa-type effective potential is used to derive the DLFC. Numerically, the DLFC are obtained via self-consistent iterative operations. We find that the DLFC are significant around the maximum of the stopping power. Furthermore, by using the two-component extended Mermin dielectric function model including the DLFC, the energy loss of a proton with an initial energy of ∼15 MeV passing through a WDP of beryllium with an electronic density around the solid value n_{e}≈3×10^{23}cm^{-3} and with temperature around ∼40 eV is estimated numerically. The numerical result is reasonably consistent with the experimental observations [A. B. Zylsta et al., Phys. Rev. Lett. 111, 215002 (2013)PRLTAO0031-900710.1103/PhysRevLett.111.215002]. Our results show that the partial ionization and the dynamic properties should be of importance for the stopping of charged particles moving in the WDP.

Atomic bound state and scattering properties of effective momentum-dependent potentials

Effective classical dynamics provide a potentially powerful avenue for modeling large-scale dynamical quantum systems. We have examined the accuracy of a Hamiltonian-based approach that employs effective momentum-dependent potentials (MDPs) within a molecular-dynamics framework through studies of atomic ground states, excited states, ionization energies, and scattering properties of continuum states. Working exclusively with the Kirschbaum-Wilets (KW) formulation with empirical MDPs [C. L. Kirschbaum and L. Wilets, Phys. Rev. A 21, 834 (1980)0556-279110.1103/PhysRevA.21.834], optimization leads to very accurate ground-state energies for several elements (e.g., N, F, Ne, Al, S, Ar, and Ca) relative to Hartree-Fock values. The KW MDP parameters obtained are found to be correlated, thereby revealing some degree of transferability in the empirically determined parameters. We have studied excited-state orbits of electron-ion pair to analyze the consequences of the MDP on the classical Coulomb catastrophe. From the optimized ground-state energies, we find that the experimental first- and second-ionization energies are fairly well predicted. Finally, electron-ion scattering was examined by comparing the predicted momentum transfer cross section to a semiclassical phase-shift calculation; optimizing the MDP parameters for the scattering process yielded rather poor results, suggesting a limitation of the use of the KW MDPs for plasmas.

Multipole expansion in plasmas: Effective interaction potentials between compound particles

In this paper, the multipole expansion method is used to determine effective interaction potentials between particles in both classical dusty plasma and dense quantum plasma. In particular, formulas for interactions of dipole-dipole and charge-dipole pairs in a classical nondegenerate plasma as well as in degenerate quantum and semiclassical plasmas were derived. The potentials describe interactions between atoms, atoms and charged particles, dust particles in the complex plasma, atoms and electrons in the degenerate plasma, and metals. Correctness of the results obtained from the multipole expansion is confirmed by their agreement with the results based on other methods of statistical physics and dielectric response function. It is shown that the method of multipole expansion can be used to derive effective interaction potentials of compound particles, if the effect of the medium on the potential of individual particles comprising compound particles is known.

Ionic transport in high-energy-density matter

Ionic transport coefficients for dense plasmas have been numerically computed using an effective Boltzmann approach. We have developed a simplified effective potential approach that yields accurate fits for all of the relevant cross sections and collision integrals. Our results have been validated with molecular-dynamics simulations for self-diffusion, interdiffusion, viscosity, and thermal conductivity. Molecular dynamics has also been used to examine the underlying assumptions of the Boltzmann approach through a categorization of behaviors of the velocity autocorrelation function in the Yukawa phase diagram. Using a velocity-dependent screening model, we examine the role of dynamical screening in transport. Implications of these results for Coulomb logarithm approaches are discussed.

Measurements of Ion Stopping Around the Bragg Peak in High-Energy-Density Plasmas

For the first time, quantitative measurements of ion stopping at energies around the Bragg peak (or peak ion stopping, which occurs at an ion velocity comparable to the average thermal electron velocity), and its dependence on electron temperature (T(e)) and electron number density (n(e)) in the range of 0.5-4.0 keV and 3×10(22) to 3×10(23) cm(-3) have been conducted, respectively. It is experimentally demonstrated that the position and amplitude of the Bragg peak varies strongly with T(e) with n(e). The importance of including quantum diffraction is also demonstrated in the stopping-power modeling of high-energy-density plasmas.

Effective potentials of interactions and thermodynamic properties of a nonideal two-temperature dense plasma

In this article a dense nonideal, nonisothermal plasma is considered. New effective screened interaction potentials taking into account quantum-mechanical diffraction and symmetry effects have been obtained. The effective potential of the ion-ion interaction in plasmas with a strongly coupled ion subsystem and semiclassical electron subsystem is presented. Based on the obtained effective potentials the analytical expressions for internal energy and the pressure of the considered plasma were obtained.

Measurement of charged-particle stopping in warm dense plasma

We measured the stopping of energetic protons in an isochorically heated solid-density Be plasma with an electron temperature of ∼32  eV, corresponding to moderately coupled [(e^{2}/a)/(k_{B}T_{e}+E_{F})∼0.3] and moderately degenerate [k_{B}T_{e}/E_{F}∼2] "warm-dense matter" (WDM) conditions. We present the first high-accuracy measurements of charged-particle energy loss through dense plasma, which shows an increased loss relative to cold matter, consistent with a reduced mean ionization potential. The data agree with stopping models based on an ad hoc treatment of free and bound electrons, as well as the average-atom local-density approximation; this work is the first test of these theories in WDM plasma.

Electric microfields in dense carbon-hydrogen plasmas

Classical molecular dynamics is used to investigate stationary and time-dependent properties of microfields in hot, solid density, electron-ion plasmas. Even at the high temperatures considered here, such simulations require quantum statistical potentials (QSPs) to mimic the essential effects of diffraction and exchange symmetry for electrons. Fortunately, key results relevant to microfield distributions are found to be insensitive to different, plausible QSP choices. Atomic processes in plasmas will depend on the time average of the microfields. It is not clear, a priori, what the time duration of this average should be. The question of how best to extract the quasistatic (low-frequency) microfield from a classical molecular dynamics simulation is explored in some detail, and the time-averaging approach we adopt involves both plasma and atomic time scale constraints. One of the major findings described in the paper is that for a large time interval, the time-averaged microfield does not significantly change. Our discussion of this suite of large simulations for plasma mixtures focuses on understanding various features and trends revealed by data for C-H plasmas having carbon fractions ranging from 0.01 to 1, and different temperatures well above TFermi.

Reproducible radiation-damage processes in proteins irradiated by intense x-ray pulses

X-ray free-electron lasers have enabled femtosecond protein nanocrystallography, a novel method to determine the structure of proteins. It allows time-resolved imaging of nanocrystals that are too small for conventional crystallography. The short pulse duration helps in overcoming the detrimental effects of radiation damage because x rays are scattered before the sample has been significantly altered. It has been suggested that, fortuitously, the diffraction process self-terminates abruptly once radiation damage destroys the crystalline order. Our calculations show that high-intensity x-ray pulses indeed trigger a cascade of damage processes in ferredoxin crystals, a particular metalloprotein of interest. However, we found that the damage process is initially not completely random. Correlations exist among the protein monomers, so that Bragg diffraction still occurs in the damaged crystals, despite significant atomic displacements. Our results show that the damage process is reproducible to a certain degree, which is potentially beneficial for the orientation step in single-molecule imaging.

Ion potential in warm dense matter: wake effects due to streaming degenerate electrons

The effective dynamically screened potential of a classical ion in a stationary flowing quantum plasma at finite temperature is investigated. This is a key quantity for thermodynamics and transport of dense plasmas in the warm-dense-matter regime. This potential has been studied before within hydrodynamic approaches or based on the zero temperature Lindhard dielectric function. Here we extend the kinetic analysis by including the effects of finite temperature and of collisions based on the Mermin dielectric function. The resulting ion potential exhibits an oscillatory structure with attractive minima (wakes) and, thus, strongly deviates from the static Yukawa potential of equilibrium plasmas. This potential is analyzed in detail for high-density plasmas with values of the Brueckner parameter in the range 0.1≤r(s)≤1 for a broad range of plasma temperature and electron streaming velocity. It is shown that wake effects become weaker with increasing temperature of the electrons. Finally, we obtain the minimal electron streaming velocity for which attraction between ions occurs. This velocity turns out to be less than the electron Fermi velocity. Our results allow for reliable predictions of the strength of wake effects in nonequilibrium quantum plasmas with fast streaming electrons showing that these effects are crucial for transport under warm-dense-matter conditions, in particular for laser-matter interaction, electron-ion temperature equilibration, and stopping power.

Dielectric function of dense plasmas, their stopping power, and sum rules

Mathematical, particularly, asymptotic properties of the random-phase approximation, Mermin approximation, and extended Mermin-type approximation of the coupled plasma dielectric function are analyzed within the method of moments. These models are generalized for two-component plasmas. Some drawbacks and advantages of the above models are pointed out. The two-component plasma stopping power is shown to be enhanced with respect to that of the electron fluid.

Kinetic theory molecular dynamics and hot dense matter: theoretical foundations

Electrons are weakly coupled in hot, dense matter that is created in high-energy-density experiments. They are also mildly quantum mechanical and the ions associated with them are classical and may be strongly coupled. In addition, the dynamical evolution of plasmas under these hot, dense matter conditions involve a variety of transport and energy exchange processes. Quantum kinetic theory is an ideal tool for treating the electrons but it is not adequate for treating the ions. Molecular dynamics is perfectly suited to describe the classical, strongly coupled ions but not the electrons. We develop a method that combines a Wigner kinetic treatment of the electrons with classical molecular dynamics for the ions. We refer to this hybrid method as "kinetic theory molecular dynamics," or KTMD. The purpose of this paper is to derive KTMD from first principles and place it on a firm theoretical foundation. The framework that KTMD provides for simulating plasmas in the hot, dense regime is particularly useful since current computational methods are generally limited by their inability to treat the dynamical quantum evolution of the electronic component. Using the N-body von Neumann equation for the electron-proton plasma, three variations of KTMD are obtained. Each variant is determined by the physical state of the plasma (e.g., collisional versus collisionless). The first variant of KTMD yields a closed set of equations consisting of a mean-field quantum kinetic equation for the electron one-particle distribution function coupled to a classical Liouville equation for the protons. The latter equation includes both proton-proton Coulombic interactions and an effective electron-proton interaction that involves the convolution of the electron density with the electron-proton Coulomb potential. The mean-field approach is then extended to incorporate equilibrium electron-proton correlations through the Singwi-Tosi-Land-Sjolander (STLS) ansatz. This is the second variant of KTMD. The STLS contribution produces an effective electron-proton interaction that involves the electron-proton structure factor, thereby extending the usual mean-field theory to correlated but near equilibrium systems. Finally, a third variant of KTMD is derived. It includes dynamical electrons and their correlations coupled to a MD description for the ions. A set of coupled equations for the one-particle electron Wigner function and the electron-electron and electron-proton correlation functions are coupled to a classical Liouville equation for the protons. This latter variation has both time and momentum dependent correlations.

Diffusivity in asymmetric Yukawa ionic mixtures in dense plasmas

In this paper we present molecular dynamics (MD) calculations of the interdiffusion coefficient for asymmetric mixed plasma for thermodynamic conditions relevant to astrophysical and inertial confinement fusion plasmas. Specifically, we consider mixtures of deuterium and argon at temperatures of 100-500 eV and a number density ∼10(25) ions/cm(3). The motion of 30,000-120,000 ions is simulated in which the ions interact via the Yukawa (screened Coulomb) potential. The electric field of the electrons is included in this effective interaction; the electrons are not simulated explicitly. The species diffusivity is then calculated using the Green-Kubo approach using an integral of the interdiffusion current autocorrelation function, a quantity calculated in the equilibrium MD simulations. Our MD simulation results show that a widely used expression relating the interdiffusion coefficient with the concentration-weighted sum of self-diffusion coefficients overestimates the interdiffusion coefficient. We argue that this effect due to cross-correlation terms in velocities is characteristic of asymmetric mixed plasmas. Comparison of the MD results with predictions of kinetic theories also shows a discrepancy with MD giving effectively a larger Coulomb logarithm.