Ionic transport in high-energy-density matter

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.

Inverse bremsstrahlung absorption rate for super-Gaussian electron distribution functions including plasma screening

We provide analytic expressions for the effective Coulomb logarithm for inverse bremsstrahlung absorption which predict significant corrections to the Langdon effect and overall absorption rate compared to previous estimates. The calculation of the collisional absorption rate of laser energy in a plasma by the inverse bremsstrahlung mechanism usually makes the approximation of a constant Coulomb logarithm. We dispense with this approximation and instead take into account the velocity dependence of the Coulomb logarithm, leading to a more accurate expression for the absorption rate valid in both classical and quantum conditions. In contrast to previous work, the laser intensity enters into the Coulomb logarithm. In most laser-plasma interactions the electron distribution function is super-Gaussian [Langdon, Phys. Rev. Lett. 44, 575 (1980)0031-900710.1103/PhysRevLett.44.575], and we find the absorption rate under these conditions is increased by as much as ≈30% compared to previous estimates at low density. In many cases of interest the correction to Langdon's predicted reduction in absorption is large; for example at Z=6 and T_{e}=400eV the Langdon prediction for the absorption is in error by a factor of ≈2. However, we also account for the additional effect of plasma screening, which predicts a reduction in absorption by a similar amount (up to ≈30%). These two effects compete to determine the overall absorption, which may be increased or decreased, depending on the conditions. The corrections can be incorporated into radiation-hydrodynamics simulation codes by replacing the familiar Coulomb logarithm with an analytic expression which depends on the super-Gaussian order "M" and the screening length.

Buoyant crystals halt the cooling of white dwarf stars

White dwarfs are stellar remnants devoid of a nuclear energy source, gradually cooling over billions of years1,2 and eventually freezing into a solid state from the inside out3,4. Recently, it was discovered that a population of freezing white dwarfs maintains a constant luminosity for a duration comparable with the age of the universe5, signalling the presence of a powerful, yet unknown, energy source that inhibits the cooling. For certain core compositions, the freezing process is predicted to trigger a solid-liquid distillation mechanism, owing to the solid phase being depleted in heavy impurities6-8. The crystals thus formed are buoyant and float up, thereby displacing heavier liquid downward and releasing gravitational energy. Here we show that distillation interrupts the cooling for billions of years and explains all the observational properties of the unusual delayed population. With a steady luminosity surpassing that of some main-sequence stars, these white dwarfs defy their conventional portrayal as dead stars. Our results highlight the existence of peculiar merger remnants9,10 and have profound implications for the use of white dwarfs in dating stellar populations11,12.

Design of the first fusion experiment to achieve target energy gain G>1

In this work we present the design of the first controlled fusion laboratory experiment to reach target gain G>1 N221204 (5 December 2022) [Phys. Rev. Lett. 132, 065102 (2024)10.1103/PhysRevLett.132.065102], performed at the National Ignition Facility, where the fusion energy produced (3.15 MJ) exceeded the amount of laser energy required to drive the target (2.05 MJ). Following the demonstration of ignition according to the Lawson criterion N210808, experiments were impacted by nonideal experimental fielding conditions, such as increased (known) target defects that seeded hydrodynamic instabilities or unintentional low-mode asymmetries from nonuniformities in the target or laser delivery, which led to reduced fusion yields less than 1 MJ. This Letter details design changes, including using an extended higher-energy laser pulse to drive a thicker high-density carbon (also known as diamond) capsule, that led to increased fusion energy output compared to N210808 as well as improved robustness for achieving high fusion energies (greater than 1 MJ) in the presence of significant low-mode asymmetries. For this design, the burnup fraction of the deuterium and tritium (DT) fuel was increased (approximately 4% fuel burnup and a target gain of approximately 1.5 compared to approximately 2% fuel burnup and target gain approximately 0.7 for N210808) as a result of increased total (DT plus capsule) areal density at maximum compression compared to N210808. Radiation-hydrodynamic simulations of this design predicted achieving target gain greater than 1 and also the magnitude of increase in fusion energy produced compared to N210808. The plasma conditions and hotspot power balance (fusion power produced vs input power and power losses) using these simulations are presented. Since the drafting of this manuscript, the results of this paper have been replicated and exceeded (N230729) in this design, together with a higher-quality diamond capsule, setting a new record of approximately 3.88MJ of fusion energy and fusion energy target gain of approximately 1.9.

Predictive scale-bridging simulations through active learning

Throughout computational science, there is a growing need to utilize the continual improvements in raw computational horsepower to achieve greater physical fidelity through scale-bridging over brute-force increases in the number of mesh elements. For instance, quantitative predictions of transport in nanoporous media, critical to hydrocarbon extraction from tight shale formations, are impossible without accounting for molecular-level interactions. Similarly, inertial confinement fusion simulations rely on numerical diffusion to simulate molecular effects such as non-local transport and mixing without truly accounting for molecular interactions. With these two disparate applications in mind, we develop a novel capability which uses an active learning approach to optimize the use of local fine-scale simulations for informing coarse-scale hydrodynamics. Our approach addresses three challenges: forecasting continuum coarse-scale trajectory to speculatively execute new fine-scale molecular dynamics calculations, dynamically updating coarse-scale from fine-scale calculations, and quantifying uncertainty in neural network models.

Measurements of ion-electron energy-transfer cross section in high-energy-density plasmas

We report on measurements of the ion-electron energy-transfer cross section utilizing low-velocity ion stopping in high-energy-density plasmas at the OMEGA laser facility. These measurements utilize a technique that leverages the close relationship between low-velocity ion stopping and ion-electron equilibration. Shock-driven implosions of capsules filled with D^{3}He gas doped with a trace amount of argon are used to generate densities and temperatures in ranges from 1×10^{23} to 2×10^{24} cm^{-3} and from 1.4 to 2.5 keV, respectively. The energy loss of 1-MeV DD tritons and 3.7-MeV D^{3}He alphas that have velocities lower than the average velocity of the thermal electrons is measured. The energy loss of these ions is used to determine the ion-electron energy-transfer cross section, which is found to be in excellent agreement with quantum-mechanical calculations in the first Born approximation. This result provides an experimental constraint on ion-electron energy transfer in high-energy-density plasmas, which impacts the modeling of alpha heating in inertial confinement fusion implosions, magnetic-field advection in stellar atmospheres, and energy balance in supernova shocks.

Lawson Criterion for Ignition Exceeded in an Inertial Fusion Experiment

For more than half a century, researchers around the world have been engaged in attempts to achieve fusion ignition as a proof of principle of various fusion concepts. Following the Lawson criterion, an ignited plasma is one where the fusion heating power is high enough to overcome all the physical processes that cool the fusion plasma, creating a positive thermodynamic feedback loop with rapidly increasing temperature. In inertially confined fusion, ignition is a state where the fusion plasma can begin "burn propagation" into surrounding cold fuel, enabling the possibility of high energy gain. While "scientific breakeven" (i.e., unity target gain) has not yet been achieved (here target gain is 0.72, 1.37 MJ of fusion for 1.92 MJ of laser energy), this Letter reports the first controlled fusion experiment, using laser indirect drive, on the National Ignition Facility to produce capsule gain (here 5.8) and reach ignition by nine different formulations of the Lawson criterion.

Design of an inertial fusion experiment exceeding the Lawson criterion for ignition

We present the design of the first igniting fusion plasma in the laboratory by Lawson's criterion that produced 1.37 MJ of fusion energy, Hybrid-E experiment N210808 (August 8, 2021) [Phys. Rev. Lett. 129, 075001 (2022)10.1103/PhysRevLett.129.075001]. This design uses the indirect drive inertial confinement fusion approach to heat and compress a central "hot spot" of deuterium-tritium (DT) fuel using a surrounding dense DT fuel piston. Ignition occurs when the heating from absorption of α particles created in the fusion process overcomes the loss mechanisms in the system for a duration of time. This letter describes key design changes which enabled a ∼3-6× increase in an ignition figure of merit (generalized Lawson criterion) [Phys. Plasmas 28, 022704 (2021)1070-664X10.1063/5.0035583, Phys. Plasmas 25, 122704 (2018)1070-664X10.1063/1.5049595]) and an eightfold increase in fusion energy output compared to predecessor experiments. We present simulations of the hot-spot conditions for experiment N210808 that show fundamentally different behavior compared to predecessor experiments and simulated metrics that are consistent with N210808 reaching for the first time in the laboratory "ignition."

Random-walk shielding-potential viscosity model for warm dense metals

We develop a model, called the "random-walk shielding-potential viscosity model" (RWSP-VM) that introduces the statistics of random-walk ions and the Debye shielding effect to describe the viscosities of warm dense metals. The viscosities of several metals with low to high atomic numbers (Be, Al, Fe, and U) are calculated using the analytical expression of RWSP-VM. Additionally, we simulate the viscosities of Fe and Be by employing the Langevin molecular dynamics (MD) and classical MD, while the MD data for Al and U are obtained from a previous work. The results of the RWSP-VM are in good agreement with the MD results, which validates the proposed model. Furthermore, we compare the RWSP-VM with the one-component plasma model and Yukawa viscosity model and show that the three models yield results in excellent agreement with each other in the regime where the RWSP-VM is applicable. These results indicate that the RWSP-VM is a universal, accurate, and highly efficient model for calculating the viscosity of metals in the warm dense state. The code of the proposed RWSP-VM is provided, and it is envisaged that it will have broad application prospects in numerous fields.

Ultracold neutral plasma expansion in a strong uniform magnetic field

In strongly magnetized neutral plasmas, electron motion is reduced perpendicular to the magnetic field direction. This changes dynamical plasma properties such as temperature equilibration, spatial density evolution, electron pressure, and thermal and electrical conductivity. In this paper we report measurements of free plasma expansion in the presence of a strong magnetic field. We image laser-induced fluorescence from an ultracold neutral Ca^{+} plasma to map the plasma size as a function of time for a range of magnetic field strengths. The asymptotic expansion velocity perpendicular to the magnetic field direction falls rapidly with increasing magnetic field strength. We observe that the initially Gaussian spatial distribution remains Gaussian throughout the expansion in both the parallel and perpendicular directions. We compare these observations with a diffusion model and with a self-similar expansion model and show that neither of these models reproduces the observed behavior over the entire range of magnetic fields used in this study. Modeling the expansion of a magnetized ultracold plasma poses a nontrivial theoretical challenge.

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.

Transport properties of a quasisymmetric binary nitrogen-oxygen mixture in the warm dense regime

Transport properties of mixtures in the warm dense matter (WDM) regime play an important role in natural astrophysics. However, a physical understanding of ionic transport properties in quasisymmetric liquid mixtures has remained elusive. Here, we present extensive ab initio molecular dynamics (AIMD) simulations on the ionic diffusion and viscosity of a quasisymmetric binary nitrogen-oxygen (N-O) mixture in a wide warm dense regime of 8-120 kK and 4.5-8.0 g/cm^{3}. Diffusion and viscosity of N-O mixtures with different compositions are obtained by using the Green-Kubo formula. Unlike asymmetric mixtures, the change of proportions in N-O mixtures slightly affects the viscosity and diffusion in the strong-coupling region. Furthermore, the AIMD results are used to build and verify a global pseudo-ion in jellium (PIJ) model for ionic transport calculations. The PIJ model succeeds in reproducing the transport properties of N-O mixtures where ionization has occurred, and provides a promising alternative approach to obtaining comparable results to AIMD simulations with relatively small computational costs. Our current results highlight the characteristic features of the quasisymmetric binary mixtures and demonstrate the applicability of the PIJ model in the WDM regime.

Multifidelity regression of sparse plasma transport data available in disparate physical regimes

Physical data are typically generated by experiments and computations in limited parameter regimes. When datasets generated using such disparate methods are combined into one dataset, the resulting dataset is typically sparse, with dense "islands" in a potentially high-dimensional parameter space, and predictions must be interpolated among such islands. Using plasma transport data as our example, we propose a multifidelity Gaussian-process regression framework that incorporates physical data from multiple sources at multiple fidelities. The impact of the proposed framework varies from little improvement over simpler approaches to qualitatively changing the prediction with consistently increased confidence in regions lacking high-fidelity data. By varying low- and high-fidelity data sources, we demonstrate an approach for determining when multifidelity Gaussian-process regression adds value over single-fidelity regression and therefore when its additional computational costs are merited. We also examine the case in which the outputs of the low- and high-fidelity models correspond to different physical quantities, one of which may be intrinsically computationally cheaper to compute. We conclude by analyzing strategies for sampling high-fidelity data for use in this framework, and we develop a simple sampling approach for reducing regression error across large gaps in data.

Multiscale simulation of plasma flows using active learning

Plasma flows encountered in high-energy-density experiments display features that differ from those of equilibrium systems. Nonequilibrium approaches such as kinetic theory (KT) capture many, if not all, of these phenomena. However, KT requires closure information, which can be computed from microscale simulations and communicated to KT. We present a concurrent heterogeneous multiscale approach that couples molecular dynamics (MD) with KT in the limit of near-equilibrium flows. To reduce the cost of gathering information from MD, we use active learning to train neural networks on MD data obtained by randomly sampling a small subset of the parameter space. We apply this method to a plasma interfacial mixing problem relevant to warm dense matter, showing considerable computational gains when compared with the full kinetic-MD approach. We find that our approach enables the probing of Coulomb coupling physics across a broad range of temperatures and densities that are inaccessible with current theoretical models.

Static and dynamic properties of multi-ionic plasma mixtures

Complex plasma mixtures with three or more components are often encountered in astrophysics or in inertial confinement fusion (ICF) experiments. For mixtures containing species with large differences in atomic number Z, the modeling needs to consider at the same time the kinetic theory for low-Z elements combined with the theory of strongly coupled plasma for high-Z elements, as well as all the intermediate situations that can appear in multicomponent systems. For such cases, we study the pair distribution functions, self-diffusions, mutual diffusion, and viscosity for ternary mixtures at extreme conditions. These quantities can be produced from first principles using orbital free molecular dynamics at the computational expense of very intensive simulations to reach good statistics. Utilizing the first-principles results as reference data, we assess the merit of a global analytic model for transport coefficients, "pseudo-ions in jellium" (PIJ), based on an isoelectronic assumption (iso-n_{e}). With a multicomponent hypernetted-chain integral equation, we verify the quality of the iso-n_{e} prescription for describing the static structure of the mixtures. This semianalytical modeling compares well with the simulation results and allows one to consider plasma mixtures not accessible to simulations. Applications are given for the mix of materials in ICF experiments. A reduction of a multicomponent mixture to an effective binary mixture is also established in the hydrodynamic limit and compared with PIJ estimations for ICF relevant mixtures.

Sculpted ultracold neutral plasmas

Ultracold neutral plasma (UNP) experiments allow for careful control of plasma properties across Coulomb coupling regimes. Here, we examine how UNPs can be used to study heterogeneous, nonequilibrium phenomena, including nonlinear waves, transport, hydrodynamics, kinetics, stopping power, and instabilities. Through a series of molecular dynamics simulations, we have explored UNPs formed with spatially modulated ionizing radiation. We have developed a computational model for such sculpted UNPs that includes an ionic screened Coulomb interaction with a spatiotemporal screening length, and Langevin-based spatial ion-electron and ion-neutral collisions. We have also developed a hydrodynamics model and have extracted its field quantities (density, flow velocity, and temperature) from the molecular dynamics simulation data, allowing us to investigate kinetics by examining moment ratios and phase-space dynamics; we find that it is possible to create UNPs that vary from nearly perfect fluids (Euler limit) to highly kinetic plasmas. We have examined plasmas in three geometries: a solid rod, a hollow rod, and a gapped slab; we have studied basic properties of these plasmas, including the spatial Coulomb coupling parameter. By varying the initial conditions, we find that we can design experimental plasmas that would allow the exploration of a wide range of phenomena, including shock and blast waves, stopping power, two-stream instabilities, and much more. Using an evaporative cooling geometry, our results suggest that much larger Coulomb couplings can be achieved, possibly in excess of 10.

Direct numerical simulation and Reynolds-averaged Navier-Stokes modeling of the sudden viscous dissipation for multicomponent turbulence

Simulations of a turbulent multicomponent fluid mixture undergoing isotropic deformations are carried out to investigate the sudden viscous dissipation. This dissipative mechanism was originally demonstrated using simulations of an incompressible single-component fluid [S. Davidovits and N. J. Fisch, Phys. Rev. Lett. 116, 105004 (2016)10.1103/PhysRevLett.116.105004]. By accounting for the convective and diffusive transfer of various species, the current work aims to increase the physical fidelity of previous simulations and their relevance to inertial confinement fusion applications. Direct numerical simulations of the compressed fluid show that the sudden viscous dissipation of turbulent kinetic energy is unchanged from the single-component scenario. More importantly, the simulations demonstrate that the mass fraction variance and covariance for the various species also exhibit a sudden viscous decay. Reynolds-averaged Navier-Stokes simulations were carried out using the k-l model to assess its ability to reproduce the sudden viscous dissipation. Results show that the standard k-l formulation does not capture the sudden decay of turbulent kinetic energy, mass-fraction variance, and mass-fraction covariance for simulations with various compression and expansion rates, or different exponents for the power-law model of viscosity. A new formulation of the k-l model that is based on previous improvements to the k-ε family of models is proposed, which leads to consistently good agreement with the direct numerical simulations for all the isotropic deformations under consideration.

Ion friction at small values of the Coulomb logarithm

Transport properties of high-energy-density plasmas are influenced by the ion collision rate. Traditionally, this rate involves the Coulomb logarithm, lnΛ. Typical values of lnΛ are ≈10-20 in kinetic theories where transport properties are dominated by weak-scattering events caused by long-range forces. The validity of these theories breaks down for strongly coupled plasmas, when lnΛ is of order one. We present measurements and simulations of collision data in strongly coupled plasmas when lnΛ is small. Experiments are carried out in the first dual-species ultracold neutral plasma (UNP), using Ca^{+} and Yb^{+} ions. We find strong collisional coupling between the different ion species in the bulk of the plasma. We simulate the plasma using a two-species fluid code that includes Coulomb logarithms derived from either a screened Coulomb potential or a the potential of mean force. We find generally good agreement between the experimental measurements and the simulations. With some improvements, the mixed Ca^{+} and Yb^{+} dual-species UNP will be a promising platform for testing theoretical expressions for lnΛ and collision cross-sections from kinetic theories through measurements of energy relaxation, stopping power, two-stream instabilities, and the evolution of sculpted distribution functions in an idealized environment in which the initial temperatures, densities, and charge states are accurately known.

Excess pressure and electric fields in nonideal plasma hydrodynamics

Nonideal plasmas have nontrivial space and time correlations, which simultaneously impact both the excess thermodynamic quantities as well as the collision processes. However, hydrodynamics models for designing and interpreting nonideal plasma experiments, such as inertial-confinement fusion experiments, typically neglect electrodynamics, although some models include electric fields indirectly through a generalized Fick's law. However, because most transport models are not computed self-consistently with the equation of state, there is double counting of the forces in the excess thermodynamic quantities and the collision terms. Here we employ the statistical mechanical hydrodynamic theory of Irving and Kirkwood [J. Chem. Phys. 18, 817 (1950)JCPSA60021-960610.1063/1.1747782] to examine inhomogeneous, nonideal plasmas that contain electric fields. We show that it is not possible to simultaneously separate terms that correspond to electric fields and excess pressure; rather, these quantities arise from the same interparticle Coulomb forces. Moreover, new terms associated with nonlocality appear in the presence of strong inhomogeneities.

Extremely Low Electron-ion Temperature Relaxation Rates in Warm Dense Hydrogen: Interplay between Quantum Electrons and Coupled Ions

Theoretical and computational modeling of nonequilibrium processes in warm dense matter represents a significant challenge. The electron-ion relaxation process in warm dense hydrogen is investigated here by nonequilibrium molecular dynamics using the constrained electron force field (CEFF) method. CEFF evolves wave packets that incorporate dynamic quantum diffraction that obviates the Coulomb catastrophe. Predictions from this model reveal temperature relaxation times as much as three times longer than prior molecular dynamics results based on quantum statistical potentials. Through analyses of energy distributions and mean free paths, this result can be traced to delocalization. Finally, an improved GMS [Gericke, Murillo, and Schlanges, Phys. Rev. E 78, 025401 (2008)PRESCM1539-375510.1103/PhysRevE.78.025401] model is proposed, in which the Coulomb logarithms are in good agreement with CEFF results.

Laser cooling of ions in a neutral plasma

Laser cooling of a neutral plasma is a challenging task because of the high temperatures typically associated with the plasma state. By using an ultracold neutral plasma created by photoionization of an ultracold atomic gas, we avoid this obstacle and demonstrate laser cooling of ions in a neutral plasma. After 135 microseconds of cooling, we observed a reduction in ion temperature by up to a factor of four, with the temperature reaching as low as 50(4) millikelvin. This pushes laboratory studies of neutral plasmas deeper into the strongly coupled regime, beyond the limits of validity of current kinetic theories for calculating transport properties. The same optical forces also retard the plasma expansion, opening avenues for neutral-plasma confinement and manipulation.

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

Analytic expressions for electron-ion temperature equilibration rates from the Lenard-Balescu equation

In this work, we elucidate the mathematical structure of the integral that arises when computing the electron-ion temperature equilibration time for a homogeneous weakly coupled plasma from the Lenard-Balescu equation. With some minor approximations, we derive an analytic formula, requiring no input Coulomb logarithm, for the equilibration rate that is valid for moderate electron-ion temperature ratios and arbitrary electron degeneracy. For large temperature ratios, we derive the necessary correction to account for the coupled-mode effect, which can be evaluated very efficiently using ordinary Gaussian quadrature.

Interfacial mixing in high-energy-density matter with a multiphysics kinetic model

We have extended a recently developed multispecies, multitemperature Bhatnagar-Gross-Krook model [Haack et al., J. Stat. Phys. 168, 822 (2017)JSTPBS0022-471510.1007/s10955-017-1824-9], to include multiphysics capabilities that enable modeling of a wider range of physical conditions. In terms of geometry, we have extended from the spatially homogeneous setting to one spatial dimension. In terms of the physics, we have included an atomic ionization model, accurate collision physics across coupling regimes, self-consistent electric fields, and degeneracy in the electronic screening. We apply the model to a warm dense matter scenario in which the ablator-fuel interface of an inertial confinement fusion target is heated, but for larger length and time scales and for much higher temperatures than can be simulated using molecular dynamics. Relative to molecular dynamics, the kinetic model greatly extends the temperature regime and the spatiotemporal scales over which we are able to model. In our numerical results we observe hydrogen from the ablator material jetting into the fuel during the early stages of the implosion and compare the relative size of various diffusion components (Fickean diffusion, electrodiffusion, and barodiffusion) that drive this process. We also examine kinetic effects, such as anisotropic distributions and velocity separation, in order to determine when this problem can be described with a hydrodynamic model.

A viscous quantum hydrodynamics model based on dynamic density functional theory

Dynamic density functional theory (DDFT) is emerging as a useful theoretical technique for modeling the dynamics of correlated systems. We extend DDFT to quantum systems for application to dense plasmas through a quantum hydrodynamics (QHD) approach. The DDFT-based QHD approach includes correlations in the the equation of state self-consistently, satisfies sum rules and includes irreversibility arising from collisions. While QHD can be used generally to model non-equilibrium, heterogeneous plasmas, we employ the DDFT-QHD framework to generate a model for the electronic dynamic structure factor, which offers an avenue for measuring hydrodynamic properties, such as transport coefficients via x-ray Thomson scattering.

Correlation and transport properties for mixtures at constant pressure and temperature

Transport properties of mixtures of elements in the dense plasma regime play an important role in natural astrophysical and experimental systems, e.g., inertial confinement fusion. We present a series of orbital-free molecular dynamics simulations on dense plasma mixtures with comparison to a global pseudo ion in jellium model. Hydrogen is mixed with elements of increasingly high atomic number (lithium, carbon, aluminum, copper, and silver) at a fixed temperature of 100 eV and constant pressure set by pure hydrogen at 2g/cm^{3}, namely, 370 Mbars. We compute ionic transport coefficients, such as self-diffusion, mutual diffusion, and viscosity for various concentrations. Small concentrations of the heavy atoms significantly change the density of the plasma and decrease the transport coefficients. The structure of the mixture evidences a strong Coulomb coupling between heavy ions and the appearance of a broad correlation peak at short distances between hydrogen atoms. The concept of an effective one component plasma is used to quantify the overcorrelation of the light element induced by the admixture of a heavy element.

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.

Transport properties of an asymmetric mixture in the dense plasma regime

We study how concentration changes ionic transport properties along isobars-isotherms for a mixture of hydrogen and silver, representative of turbulent layers relevant to inertial confinement fusion and astrophysics. Hydrogen will typically be fully ionized while silver will be only partially ionized but can have a large effective charge. This will lead to very different physical conditions for the H and Ag. Large first principles orbital free molecular dynamics simulations are performed and the resulting transport properties are analyzed. Comparisons are made with transport theory in the kinetic regime and in the coupled regime. The addition of a small amount of heavy element in a light material has a dramatic effect on viscosity and diffusion of the mixture. This effect is explained through kinetic theory as a manifestation of a crossover between classical diffusion and Lorentz diffusion.