Extent of validity of the hydrodynamic description of ions in dense plasmas

Ionization calculations using classical molecular dynamics

By performing an ensemble of molecular dynamics simulations, the model-dependent ionization state is computed for strongly interacting systems self-consistently. This is accomplished through a free energy minimization framework based on the technique of thermodynamic integration. To illustrate the method, two simple models applicable to partially ionized hydrogen plasma are presented in which pair potentials are employed between ions and neutral particles. Within the models, electrons are either bound in the hydrogen ground state or distributed in a uniform charge-neutralizing background. Particular attention is given to the transition between atomic gas and ionized plasma, where the effect of neutral interactions is explored beyond commonly used models in the chemical picture. Furthermore, pressure ionization is observed when short-range repulsion effects are included between neutrals. The developed technique is general, and we discuss the applicability to a variety of molecular dynamics models for partially ionized warm dense matter.

Modeling of warm dense hydrogen via explicit real-time electron dynamics: Dynamic structure factors

We present two methods for computing the dynamic structure factor for warm dense hydrogen without invoking either the Born-Oppenheimer approximation or the Chihara decomposition, by employing a wave-packet description that resolves the electron dynamics during ion evolution. First, a semiclassical method is discussed, which is corrected based on known quantum constraints, and second, a direct computation of the density response function within the molecular dynamics. The wave-packet models are compared to PIMC and DFT-MD for the static and low-frequency behavior. For the high-frequency behavior the models recover the expected behavior in the limits of small and large momentum transfers and show the characteristic flattening of the plasmon dispersion for intermediate momentum transfers due to interactions, in agreement with commonly used models for x-ray Thomson scattering. By modeling the electrons and ions on an equal footing, both the ion and free electron part of the spectrum can now be treated within a single framework where we simultaneously resolve the ion-acoustic and plasmon mode, with a self-consistent description of collisions and screening.

Development of a new quantum trajectory molecular dynamics framework

An extension to the wave packet description of quantum plasmas is presented, where the wave packet can be elongated in arbitrary directions. A generalized Ewald summation is constructed for the wave packet models accounting for long-range Coulomb interactions and fermionic effects are approximated by purpose-built Pauli potentials, self-consistent with the wave packets used. We demonstrate its numerical implementation with good parallel support and close to linear scaling in particle number, used for comparisons with the more common wave packet employing isotropic states. Ground state and thermal properties are compared between the models with differences occurring primarily in the electronic subsystem. Especially, the electrical conductivity of dense hydrogen is investigated where a 15% increase in DC conductivity can be seen in our wave packet model compared with other models. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Collective ion dynamics in Coulomb one-component plasmas within the self-consistent relaxation theory

In this paper, we present the theoretical formalism describing the collective ion dynamics of the nonideal Coulomb classical one-component plasmas on the basis of the self-consistent relaxation theory. The theory is adapted to account for correlations between the frequency relaxation parameters that characterize the three- and four-particle dynamics and the parameters associated with the two-particle dynamics. The dynamic structure factor spectra and dispersion characteristics calculated for a wide range of wave numbers are in agreement with the molecular dynamics simulation data and the results obtained with the theory of the frequency moments. The proposed formalism reproduces all the features inherent to the Coulomb one-component plasmas and requires only knowledge of the coupling parameter and the information about the structure.

Strongly coupled Yukawa trilayer liquid: Structure and dynamics

The equilibrium structure and the dispersion relations of collective excitations in trilayer Yukawa systems in the strongly coupled liquid regime are examined. The equilibrium correlations reveal a variety of structures in the liquid phase, reminiscent of the corresponding structures in the solid phase. At small layer separation substitutional disorder becomes the governing feature. Theoretical dispersion relations are obtained by applying the quasilocalized charge approximation (QLCA) formalism, while numerical data are generated by microcanonical molecular dynamics simulations. The dispersions and polarizations of the collective excitations obtained through both of these methods are compared and discussed in detail. We find that the QLCA method is, in general, very satisfactory, but that there are phenomena not covered by the QLCA. In particular, by analyzing the dynamical longitudinal and transverse current fluctuation spectra we discover the existence of a structure not related to the collective mode spectra. This also provides insight into the long-standing problem of the gap frequency discrepancy, observed in strongly coupled layered systems in earlier studies.

Sound speed in Yukawa one-component plasmas across coupling regimes

A many-body system of charged particles interacting via a pairwise Yukawa potential, the so-called Yukawa one-component plasma (YOCP), is a good approximation for a variety of physical systems. Such systems are completely characterized by two parameters: the screening parameter, κ, and the nominal coupling strength, Γ. It is well known that the collective spectrum of the YOCP is governed by a longitudinal acoustic mode, both in the weakly and strongly coupled regimes. In the long-wavelength limit, the linear term in the dispersion (i.e., ω=sk) defines the sound speed s. We study the evolution of this latter quantity from the weak- through the strong-coupling regimes by analyzing the dynamic structure function S(k,ω) in the low-frequency domain. Depending on the values of Γ and κ and w=s/v_{th} (i.e., the ratio between the phase velocity of the wave and the thermal speed of the particles), we identify five domains in the (κ,Γ) parameter space in which the physical behavior of the YOCP exhibits different features. The competing physical processes are the collective Coulomb-like versus binary-collision-dominated behavior and the individual particle motion versus quasilocalization. Our principal tool of investigation is molecular dynamics (MD) computer simulation from which we obtain S(k,ω). Recent improvements in the simulation technique have allowed us to obtain a large body of high-quality data in the range Γ={0.1-10000} and κ={0.5-5}. The theoretical results based on various models are compared in order to see which one provides the most cogent physical description and the best agreement with MD data in the different domains.

High-frequency response of classical strongly coupled plasmas

The dynamic structure factor (DSF) of the Yukawa system is here obtained with highly converged molecular dynamics (MD) over the entire liquid phase. The data provide a rigorous test of theoretical models of ion-acoustic wave-dispersion relations, the intermediate scattering function, and the high-frequency response. We compare our MD results with seven diverse models, finding good agreement among those that enforce the three basic sum rules for dispersion properties, although one of the models has previously unreported spurious peaks. The MD simulations reveal that at intermediate frequencies ω, the high-frequency response of the DSF follows a power law, going approximately as ω^{-p}, where p>0, and p shows nontrivial dependencies on the wave vector q and the plasma parameters κ and Γ. In contrast, among the seven comparison models, the predicted high-frequency response is found to be independent of {q,κ,Γ}. This high-frequency power suggests a useful fitting form. In addition, these results expose limitations of several models and, moreover, suggest that some approaches are difficult or impossible to extend because of the lack of finite moments. We also find the double-plasmon resonance peak in our MD simulations that none of the theoretical models predicts.

Thomson scattering cross section in a magnetized, high-density plasma

We calculate the Thomson scattering cross section in a nonrelativistic, magnetized, high-density plasma-in a regime where collective excitations can be described by magnetohydrodynamics. We show that, in addition to cyclotron resonances and an elastic peak, the cross section exhibits two pairs of peaks associated with slow and fast magnetosonic waves; by contrast, the cross section arising in pure hydrodynamics possesses just a single pair of Brillouin peaks. Both the position and the width of these magnetosonic-wave peaks depend on the ambient magnetic field and temperature, as well as transport and thermodynamic coefficients, and so can therefore serve as a diagnostic tool for plasma properties that are otherwise challenging to measure.

Dynamical structure factor of strongly coupled ions in a dense quantum plasma

The dynamical structure factor (DSF) of strongly coupled ions in dense plasmas with partially and strongly degenerate electrons is investigated. The main focus is on the impact of electronic correlations (nonideality) on the ionic DSF. The latter is computed by carrying out molecular dynamics (MD) simulations with a screened ion-ion interaction potential. The electronic screening is taken into account by invoking the Singwi-Tosi-Land-Sjölander approximation, and it is compared to the MD simulation data obtained considering the electronic screening in the random phase approximation and using the Yukawa potential. We find that electronic correlations lead to lower values of the ion-acoustic mode frequencies and to an extension of the applicability limit with respect to the wave-number of a hydrodynamic description. Moreover, we show that even in the limit of weak electronic coupling, electronic correlations have a nonnegligible impact on the ionic longitudinal sound speed. Additionally, the applicability of the Yukawa potential with an adjustable screening parameter is discussed, which will be of interest, e.g., for the interpretation of experimental results for the ionic DSF of dense plasmas.

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.

Theory of density fluctuations in strongly radiative plasmas

Derivation of the dynamic structure factor, an important parameter linking experimental and theoretical work in dense plasmas, is possible starting from hydrodynamic equations. Here we obtain, by modifying the governing hydrodynamic equations, a new form of the dynamic structure factor which includes radiative terms. The inclusion of such terms has an effect on the structure factor at high temperatures, which suggests that its effect must be taken into consideration in such regimes.

Ion-ion dynamic structure factor of warm dense mixtures

The ion-ion dynamic structure factor of warm dense matter is determined using the recently developed pseudoatom molecular dynamics method [Starrett et al., Phys. Rev. E 91, 013104 (2015)]. The method uses density functional theory to determine ion-ion pair interaction potentials that have no free parameters. These potentials are used in classical molecular dynamics simulations. This constitutes a computationally efficient and realistic model of dense plasmas. Comparison with recently published simulations of the ion-ion dynamic structure factor and sound speed of warm dense aluminum finds good to reasonable agreement. Using this method, we make predictions of the ion-ion dynamical structure factor and sound speed of a warm dense mixture-equimolar carbon-hydrogen. This material is commonly used as an ablator in inertial confinement fusion capsules, and our results are amenable to direct experimental measurement.

Fluid approach to evaluate sound velocity in Yukawa systems and complex plasmas

The conventional fluid description of multicomponent plasma, supplemented by an appropriate equation of state for the macroparticle component, is used to evaluate the longitudinal sound velocity of Yukawa fluids. The obtained results are in very good agreement with those obtained earlier employing the quasilocalized charge approximation and molecular dynamics simulations in a rather broad parameter regime. Thus, a simple yet accurate tool to estimate the sound velocity across coupling regimes is proposed, which can be particularly helpful in estimating the dust-acoustic velocity in strongly coupled dusty (complex) plasmas. It is shown that, within the present approach, the sound velocity is completely determined by particle-particle correlations and the neutralizing medium (plasma), apart from providing screening of the Coulomb interaction, has no other effect on the sound propagation. The ratio of the actual sound velocity to its "ideal gas" (weak coupling) scale only weakly depends on the coupling strength in the fluid regime but exhibits a pronounced decrease with the increase of the screening strength. The limitations of the present approach in applications to real complex plasmas are briefly discussed.

Transverse current fluctuations in the Yukawa one-component plasma

Using numerical simulations, we investigate the wave number and frequency dependent transverse current correlation function CT(k,ω) of a single-component fluid with Yukawa interaction potential, also known as the Yukawa one-component plasma. The transverse current correlation function is an important quantity because it contains the microscopic details of the viscoelastic behavior of the fluid. We show that, in the region of densities and temperatures in which shear waves do not propagate, the dynamics of the system are in striking agreement with a simple model of generalized hydrodynamics. As either the density is increased or the temperature decreased, the transverse current correlation function shows additional structure that the simple models fail to capture.

Fluctuating hydrodynamics of multispecies nonreactive mixtures

In this paper we discuss the formulation of the fluctuating Navier-Stokes equations for multispecies, nonreactive fluids. In particular, we establish a form suitable for numerical solution of the resulting stochastic partial differential equations. An accurate and efficient numerical scheme, based on our previous methods for single species and binary mixtures, is presented and tested at equilibrium as well as for a variety of nonequilibrium problems. These include the study of giant nonequilibrium concentration fluctuations in a ternary mixture in the presence of a diffusion barrier, the triggering of a Rayleigh-Taylor instability by diffusion in a four-species mixture, as well as reverse diffusion in a ternary mixture. Good agreement with theory and experiment demonstrates that the formulation is robust and can serve as a useful tool in the study of thermal fluctuations for multispecies fluids.

Orbital-free density-functional theory simulations of the dynamic structure factor of warm dense aluminum

Here, we report orbital-free density-functional theory (OF DFT) molecular dynamics simulations of the dynamic ion structure factor of warm solid density aluminum at T=0.5 eV and T=5 eV. We validate the OF DFT method in the warm dense matter regime through comparison of the static and thermodynamic properties with the more complete Kohn-Sham DFT. This extension of OF DFT to dynamic properties indicates that previously used models based on classical molecular dynamics may be inadequate to capture fully the low frequency dynamics of the response function.

Hydrodynamic theory for ion structure and stopping power in quantum plasmas

We present a theory for the dynamical ion structure factor (DISF) and ion stopping power in an unmagnetized collisional quantum plasma with degenerate electron fluids and nondegenerate strongly correlated ion fluids. Our theory is based on the fluctuation dissipation theorem and the quantum plasma dielectric constant that is deduced from a linearized viscoelastic quantum hydrodynamical (LVQHD) model. The latter incorporates the essential physics of quantum forces, which are associated with the quantum statistical pressure, electron-exchange, and electron-correlation effects, the quantum electron recoil effect caused by the dispersion of overlapping electron wave functions that control the dynamics of degenerate electron fluids, and the viscoelastic properties of strongly correlated ion fluids. Both degenerate electrons and nondegenerate strongly correlated ions are coupled with each other via the space charge electric force. Thus, our LVQHD theory is valid for a collisional quantum plasma at atomic scales with a wide range of the ion coupling parameter, the plasma composition, and plasma number densities that are relevant for compressed plasmas in laboratories (inertial confinement fusion schemes) and in astrophysical environments (e.g., warm dense matter and the cores of white dwarf stars). It is found that quantum electron effects and viscoelastic properties of strongly correlated ions significantly affect the features of the DISF and the ion stopping power (ISP). Unlike previous theories, which have studied ion correlations in terms of the ion coupling parameter, by neglecting the essential physics of collective effects that are competing among each other, we have here developed a method to evaluate the dependence of the plasma static and dynamical features in terms of individual parameters, like the Wigner-Seitz radius, the ion atomic number, and the ion temperature. It is found that due to the complex nature of charge screening in quantum plasmas, the ion coupling parameter alone cannot be a good measure for determining ion correlation effects in a collisional quantum plasma, and such a characteristic of a dense quantum plasma should be evaluated against each of the plasma parameters involved. The present investigation thus provides testable predictions for the DISF and ISP and is henceforth applicable to a wide range of compressed plasma categories ranging from laboratory to astrophysical warm dense matter.

Dynamic ion structure factor of warm dense matter

The dynamics of the ion structure in warm dense matter is determined by molecular dynamics simulations using an effective ion-ion potential. This potential is obtained from ab initio simulations and has a strong short-range repulsion added to a screened Coulomb potential. Models based on static or dynamic local field corrections are found to be insufficient to describe the data. An extended Mermin approach, a hydrodynamic model, and the method of moments with local constraints are capable of reproducing the numerical results but have rather limited predictive powers as they all need some numerical data as input. The method of moments is found to be the most promising.

Diffusion in ionic mixtures across coupling regimes

Molecular dynamics simulations are used to investigate the diffusion properties of one-component plasmas and binary ionic mixtures from the weakly to the strongly coupled regimes. A physically motivated model for the diffusivities is proposed that reproduces the simulation data and gives insight into the nature of ionic motions and interactions in plasmas across the coupling regimes. The model extends the widely used Chapman-Spitzer theory from the weakly to the moderately coupled regime. In the strongly coupled regime, diffusion is modeled in terms of thermally activated jumps between equilibrium positions separated by an energy barrier. The basic ideas discussed are applicable to the study of other transport coefficients.

Comparative merits of the memory function and dynamic local-field correction of the classical one-component plasma

The complementarity of the liquid and plasma descriptions of the classical one-component plasma is explored by studying wave number and frequency dependent dynamical quantities: the dynamical structure factor (DSF) and the dynamic local field correction (LFC). Accurate molecular dynamics (MD) simulations are used to validate and test models of the DSF and LFC. Our simulations, which span the entire fluid regime (Γ=0.1-175), show that the DSF is very well represented by a simple and well known memory function model of generalized hydrodynamics. On the other hand, the LFC, which we have computed using MD for the first time, is not well described by existing models.

Quantum hydrodynamics of strongly coupled electron fluids

We have extended the classical hydrodynamics formalism to include nonlocal quantum behavior via the phenomenological Bohm potential. We have then solved the quantum hydrodynamics equations to derive an expression for the dynamical structure factor, which describes the density-density correlations in the system. This formalism can be applied to high density strongly coupled electron fluids in the long wavelength domain. We show that at densities above 7×10(25) cm(-3) there are significant differences in the dispersion relation. Future experiments at large laser facilities could provide an experimental test of the theory.

Density fluctuations in the Yukawa one-component plasma: an accurate model for the dynamical structure factor

Using numerical simulations, we investigate the equilibrium dynamics of a single-component fluid with Yukawa interaction potential. We show that, for a wide range of densities and temperatures, the dynamics of the system are in striking agreement with a simple model of generalized hydrodynamics. Since the Yukawa potential can describe the ion-ion interactions in a plasma, our results have significant applicability for both analyzing and interpreting the results of x-ray scattering data from high-power lasers and fourth-generation light sources.