Determination of the shear viscosity of the one-component plasma

Classical scattering and fragmentation of clusters of ions in helical confinement

We explore the scattering dynamics of classical Coulomb-interacting clusters of ions confined to a helical geometry. Ion clusters of equally charged particles constrained to a helix can form many-body bound states; i.e., they exhibit stable motion of Coulomb-interacting identical ions. We analyze the scattering and fragmentation behavior of two ion clusters, thereby understanding the rich phenomenology of their dynamics. The scattering dynamics is complex in the sense that it exhibits cascades of decay processes involving strongly varying cluster sizes. These processes are governed by the internal energy flow and the underlying oscillatory many-body potential. We specifically focus on the impact of the collision energy on the dynamics of individual ions during and immediately after the collision of two clusters and on the internal dynamics of ion clusters that are excited during a cluster collision.

Intrinsic bulk viscosity of the one-component plasma

The intrinsic bulk viscosity of the one-component plasma (OCP) is computed and analyzed using equilibrium molecular dynamics simulations and the Green-Kubo formalism. It is found that bulk viscosity exhibits a maximum at Γ≈1, corresponding to the condition that the average kinetic energy of particles equals the potential energy at the average interparticle spacing. The weakly coupled and strongly coupled limits are analyzed and used to construct a model that captures the full range of coupling strengths simulated: Γ≈10^{-2}-10^{2}. Simulations are also run of the Yukawa one-component plasma (YOCP) in order to understand the impact of electron screening. It is found that electron screening leads to a smaller bulk viscosity due to a reduction in the excess heat capacity of the system. Bulk viscosity is shown to be at least an order of magnitude smaller than shear viscosity in both the OCP and YOCP. The generalized frequency-dependent bulk viscosity coefficient is also analyzed. This is found to exhibit a peak near twice the plasma frequency in strongly coupled conditions, which is associated with the oscillatory decay observed in the bulk viscosity autocorrelation function. The generalized shear and bulk viscosity coefficients are found to have a similar magnitude for ω≳2ω_{p} at strongly coupled conditions.

Stokes-Einstein Relation in Different Models of Water

The purpose of this paper is to discuss to which extent a microscopic version of the Stokes-Einstein (SE) relation without the hydrodynamic radius applies to liquid water. We demonstrate that the self-diffusion and shear viscosity data for five popular water models, recently reported by Ando [J. Chem. Phys. 159, 101102 (2023)], are in excellent agreement with the SE relation. The agreement with experimental results is also quite impressive. The limitations on the applicability of the SE relation are briefly discussed.

Quasiuniversal behavior of shear relaxation times in simple fluids

We calculate the shear relaxation times in four important simple monatomic model fluids: Lennard-Jones, Yukawa, soft-sphere, and hard-sphere fluids. It is observed that in properly reduced units, the shear relaxation times exhibit quasiuniversal behavior when the density increases from the gaslike low values to the high-density regime near crystallization. They first decrease with density at low densities, reach minima at moderate densities, and then increase toward the freezing point. The reduced relaxation times at the minima and at the fluid-solid phase transition are all comparable for the various systems investigated, despite more than ten orders of magnitude difference in real systems. Important implications of these results are discussed.

Entropy of strongly coupled Yukawa fluids

The entropy of strongly coupled Yukawa fluids is discussed from several perspectives. First, it is demonstrated that a vibrational paradigm of atomic dynamics in dense fluids can be used to obtain a simple and accurate estimate of the entropy without any adjustable parameters. Second, it is explained why a quasiuniversal value of the excess entropy of simple fluids at the freezing point should be expected, and it is demonstrated that a remaining very weak dependence of the freezing point entropy on the screening parameter in the Yukawa fluid can be described by a simple linear function. Third, a scaling of the excess entropy with the freezing temperature is examined, a modified form of the Rosenfeld-Tarazona scaling is put forward, and some consequences are briefly discussed. Fourth, the location of the Frenkel line on the phase diagram of Yukawa systems is discussed in terms of the excess entropy and compared with some predictions made in the literature. Fifth, the excess entropy scaling of the transport coefficients (self-diffusion, viscosity, and thermal conductivity) is reexamined using the contemporary datasets for the transport properties of Yukawa fluids. The results could be of particular interest in the context of complex (dusty) plasmas, colloidal suspensions, electrolytes, and other related systems with soft pairwise interactions.

Bulk viscosity of the rigid rotor one-component plasma

Bulk viscosity of a plasma consisting of strongly coupled diatomic ions is computed using molecular dynamics simulations. The simulations are based on the rigid rotor one-component plasma, which is introduced as a model system that adds two degrees of molecular rotation to the traditional one-component plasma. It is characterized by two parameters: the Coulomb coupling parameter, Γ, and the bond length parameter, Ω. Results show that the long-range nature of the Coulomb potential can lead to long rotational relaxation times, which in turn yield large values for bulk viscosity. The bulk-to-shear viscosity ratio is found to span from small to large values depending on the values of Γ and Ω. Although bulk viscosity is often neglected in plasma modeling, these results motivate that it can be large in molecular plasmas with rotational degrees of freedom.

Freezing density scaling of transport coefficients in the Weeks-Chandler-Andersen fluid

It is shown that the transport coefficients (self-diffusion, shear viscosity, and thermal conductivity) of the Weeks-Chandler-Andersen (WCA) fluid along isotherms exhibit a freezing density scaling (FDS). The functional form of this FDS is essentially the same or closely related to those in the Lennard-Jones fluid, hard-sphere fluid, and some liquefied noble gases. This proves that this FDS represents a quasi-universal corresponding state principle for simple classical fluids with steep interactions. Some related aspects, such as a Stokes-Einstein relation without a hydrodynamic diameter and gas-to-liquid dynamical crossover, are briefly discussed. Simple fitting formulas for the transport coefficients of the dense WCA fluid are suggested.

System Size Dependence of the Diffusion Coefficients in MD Simulations: A Simple Correction Formula for Pure Dense Fluids

A practical correction formula relating the self-diffusion coefficient of dense liquids from molecular dynamics (MD) simulations with periodic boundary conditions to the self-diffusion coefficient in the thermodynamic limit is discussed. This formula applies to pure dense fluids and has a very simple form D = D0(1 - γN-1/3), where D0 is the self-diffusion coefficient in the thermodynamic limit and N is the number of particles in the simulation. The numerical factor γ is dependent on the geometry of the simulation cell. Remarkably, γ ≃ 1.0 is the most popular cubic geometry. The success of this formula is supported by results from MD simulations, including very recent simulations with a "magic" simulation geometry.

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.

Self-Diffusion in Simple Liquids as a Random Walk Process

It is demonstrated that self-diffusion in dense liquids can be considered a random walk process; its characteristic length and time scales are identified. This represents an alternative to the often assumed hopping mechanism of diffusion in the liquid state. The approach is illustrated using the one-component plasma model.

Excess entropy and Stokes-Einstein relation in simple fluids

The Stokes-Einstein (SE) relation between the self-diffusion and shear viscosity coefficients operates in sufficiently dense liquids not too far from the liquid-solid phase transition. By considering four simple model systems with very different pairwise interaction potentials (Lennard-Jones, Coulomb, Debye-Hückel or screened Coulomb, and the hard sphere limit) we identify where exactly on the respective phase diagrams the SE relation holds. It appears that the reduced excess entropy s_{ex} can be used as a suitable indicator of the validity of the SE relation. In all cases considered the onset of SE relation validity occurs at approximately s_{ex}≲-2. In addition, we demonstrate that the line separating gaslike and liquidlike fluid behaviours on the phase diagram is roughly characterized by s_{ex}≃-1.

Entropy of simple fluids with repulsive interactions near freezing

Among different thermodynamic properties of liquids, the entropy is one of the hardest quantities to estimate. Therefore, the development of models allowing accurate estimations of the entropy for different mechanisms of interatomic interactions represents an important problem. Here, we propose a method for estimating the excess entropy of simple liquids not too far from the liquid-solid phase transition. The method represents a variant of cell theory, which particularly emphasizes relations between liquid state thermodynamics and collective modes properties. The method is applied to calculate the excess entropy of inverse-power-law fluids with ∝r^-n repulsive interactions. The covered range of potential softness is extremely wide, including the very soft Coulomb (n = 1) case, much steeper n = 6 and n = 12 cases, and the opposite hard-sphere interaction limit (n = ∞). An overall reasonably good agreement between the method's outcome and existing "exact" results is documented at sufficiently high fluid densities. Its applicability condition can be conveniently formulated in terms of the excess entropy itself. The method is also applied to the Lennard-Jones potential but demonstrates considerably lower accuracy in this case. Our results should be relevant to a broad range of liquid systems that can be described with isotropic repulsive interactions, including liquid metals, macromolecular systems, globular proteins, and colloidal suspensions.

Prandtl Number in Classical Hard-Sphere and One-Component Plasma Fluids

The Prandtl number is evaluated for the three-dimensional hard-sphere and one-component plasma fluids, from the dilute weakly coupled regime up to a dense strongly coupled regime near the fluid-solid phase transition. In both cases, numerical values of order unity are obtained. The Prandtl number increases on approaching the freezing point, where it reaches a quasi-universal value for simple dielectric fluids of about ≃1.7. Relations to two-dimensional fluids are briefly discussed.

Viscosity of the magnetized strongly coupled one-component plasma

The viscosity tensor of the magnetized one-component plasma, consisting of five independent shear viscosity coefficients, a bulk viscosity coefficient, and a cross coefficient, is computed using equilibrium molecular dynamics simulations and the Green-Kubo relations. A broad range of Coulomb coupling and magnetization strength conditions are studied. Magnetization is found to strongly influence the shear viscosity coefficients when the gyrofrequency exceeds the Coulomb collision frequency. Three regimes are identified as the Coulomb coupling strength and magnetization strength are varied. The Green-Kubo relations are used to separate kinetic and potential energy contributions to each viscosity coefficient, showing how each contribution depends upon the magnetization strength. The shear viscosity coefficient associated with the component of the pressure tensor parallel to the magnetic field, and the two coefficients associated with the component perpendicular to the magnetic field, are all found to merge to a common value at strong Coulomb coupling.

Vibrational model of thermal conduction for fluids with soft interactions

A vibrational model of heat transfer in simple liquids with soft pairwise interatomic interactions is discussed. A general expression is derived, which involves an averaging over the liquid collective mode excitation spectrum. The model is applied to quantify heat transfer in a dense Lennard-Jones liquid and a strongly coupled one-component plasma. Remarkable agreement with the available numerical results is documented. A similar picture does not apply to the momentum transfer and shear viscosity of liquids.

Temperature relaxation in binary hard-sphere mixture system: Molecular dynamics and kinetic theory study

Computational schemes to describe the temperature relaxation in the binary hard-sphere mixture system are given on the basis of molecular dynamics (MD) simulation and renormalized kinetic theory. Event-driven MD simulations are carried out for three model systems in which the initial temperatures and the ratios of diameter and mass of two components are different to study the temporal evolution of each component temperature in nanoscale molecular conditions mimicking those in living cells. On the other hand, the temperature changes of the two components are also described in terms of a mean-field kinetic theory with the correlation functions calculated in the Percus-Yevick approximation. The calculated results by both the computational approaches have shown fair agreement with each other, whereas slight deviations have been found in the temporal range of femto- to picoseconds when the initial temperatures of the two components are significantly different, such as 300 K vs 1000 K. This discrepancy can be ascribed to the fast intra-component temperature relaxation assumed in the kinetic theory, and its violation in the MD simulations can be evaluated in terms of the Kullback-Leibler divergence between the equilibrated Maxwell-Boltzmann distribution at each temperature and the actual non-equilibrium velocity distribution realized in the MD. Thus, the present analysis provides a quantitative basis for addressing the temperature inhomogeneities experimentally observed in nanoscale crowding conditions.

Benchmarking the effective one-component plasma model for warm dense neon and krypton within quantum molecular dynamics simulation

The effective one-component plasma (EOCP) model has provided an efficient approach to obtaining many important thermophysical parameters of hot dense matter [J. Clérouin, et al., Phys. Rev. Lett. 116, 115003 (2016)PRLTAO0031-900710.1103/PhysRevLett.116.115003]. In this paper, we perform extensive quantum molecular dynamics (QMD) simulations to determine the equations of state, ionic structures, and ionic transport properties of neon and krypton within the warm dense matter (WDM) regime where the density (ρ) is up to 12 g/cm^{3} and the temperature (T) is up to 100 kK. The simulated data are then used as a benchmark to explicitly evaluate the EOCP and Yukawa models. It is found that, within present ρ-T regime, the EOCP model can excellently reproduce the diffusion and viscosity coefficients of neon and krypton due to the fact that this model defines a system which nearly reproduces the actual physical states of WDM. Therefore, the EOCP model may be a promising alternative approach to reasonably predicting the transport behaviors of matter in WDM regime at lower QMD computational cost. The evaluation of Yukawa model shows that the consideration of the energy level broadening effect in the average atom model is necessary. Finally, with the help of EOCP model, the Stokes-Einstein relationships about neon and krypton are discussed, and fruitful plasma parameters as well as a practical ρ-T-dependent formula of the effective coupling parameter are obtained. These results not only provide valuable information for future theoretical and experimental studies on dense neon and krypton but also reveal the applicability of the EOCP model and the limitation of the Yukawa model in WDM regime and further support the continuing search for a unified description of ionic transport in dense plasma.

Testing thermal conductivity models with equilibrium molecular dynamics simulations of the one-component plasma

Equilibrium molecular dynamics simulations are used to calculate the thermal conductivity of the one-component plasma via the Green-Kubo formalism over a broad range of Coulomb coupling strength, 0.1≤Γ≤180. These simulations address previous discrepancies between computations using equilibrium versus nonequilibrium methods. Analysis of heat flux autocorrelation functions show that very long (6×10^{5}ω_{p}^{-1}) time series are needed to reduce the noise level to allow ≲2% accuracy. The simulations provide accurate data for Γ≲1. This enables a test of the traditional Landau-Spitzer theory, which is found to agree with the simulations for Γ≲0.3. It also enables tests of theories to address moderate and strong Coulomb coupling. Two are found to provide accurate extensions to the moderate coupling regime of Γ≲10, but none are accurate in the Γ≳10 regime where potential energy transport and coupling between mass flow and stress dominate thermal conduction.

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.

Overestimation of Viscosity by the Green-Kubo Method in a Dusty Plasma Experiment

The Green-Kubo (GK) method is widely used in simulations of strongly coupled plasmas to obtain the viscosity coefficient. However, the method's applicability, which is often taken for granted, has not been tested experimentally. We report an experimental test using a two-dimensional strongly coupled dusty plasma. We find that the GK viscosity is ≈60% larger than the result of a benchmark hydrodynamic method, obtained in the same experiment with the same conditions except for the presence of shear.

Ionic Transport Coefficients of Dense Plasmas without Molecular Dynamics

We present a theoretical model that allows a fast and accurate evaluation of ionic transport properties of realistic plasmas spanning from warm and dense to hot and dilute conditions, including mixtures. This is achieved by combining a recent kinetic theory based on effective interaction potentials with a model for the equilibrium radial density distribution based on an average atom model and the integral equations theory of fluids. The model should find broad use in applications where nonideal plasma conditions are traversed, including inertial confinement fusion, compact astrophysical objects, solar and extrasolar planets, and numerous present-day high energy density laboratory experiments.

Ionic and electronic transport properties in dense plasmas by orbital-free density functional theory

We validate the application of our recent orbital-free density functional theory (DFT) approach [Phys. Rev. Lett. 113, 155006 (2014);] for the calculation of ionic and electronic transport properties of dense plasmas. To this end, we calculate the self-diffusion coefficient, the viscosity coefficient, the electrical and thermal conductivities, and the reflectivity coefficient of hydrogen and aluminum plasmas. Very good agreement is found with orbital-based Kohn-Sham DFT calculations at lower temperatures. Because the computational costs of the method do not increase with temperature, we can produce results at much higher temperatures than is accessible by the Kohn-Sham method. Our results for warm dense aluminum at solid density are inconsistent with the recent experimental results reported by Sperling et al. [Phys. Rev. Lett. 115, 115001 (2015)].

Shear viscosity for dense plasmas by equilibrium molecular dynamics in asymmetric Yukawa ionic mixtures

We present molecular dynamics (MD) calculations of shear viscosity 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 of 10^{25} ions/cc. 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. Shear viscosity is calculated using the Green-Kubo approach with an integral of the shear stress autocorrelation function, a quantity calculated in the equilibrium MD simulations. We systematically study different mixtures through a series of simulations with increasing fraction of the minority high-Z element (Ar) in the D-Ar plasma mixture. In the more weakly coupled plasmas, at 500 eV and low Ar fractions, results from MD compare very well with Chapman-Enskog kinetic results. In the more strongly coupled plasmas, the kinetic theory does not agree well with the MD results. We develop a simple model that interpolates between classical kinetic theories at weak coupling and the Murillo Yukawa viscosity model at higher coupling. This hybrid kinetics-MD viscosity model agrees well with the MD results over the conditions simulated, ranging from moderately weakly coupled to moderately strongly coupled asymmetric plasma mixtures.

Modified Enskog kinetic theory for strongly coupled plasmas

Concepts underlying the Enskog kinetic theory of hard-spheres are applied to include short-range correlation effects in a model for transport coefficients of strongly coupled plasmas. The approach is based on an extension of the effective potential transport theory [S. D. Baalrud and J. Daligault, Phys. Rev. Lett. 110, 235001 (2013)] to include an exclusion radius surrounding individual charged particles that is associated with Coulomb repulsion. This is obtained by analogy with the finite size of hard spheres in Enskog's theory. Predictions for the self-diffusion and shear viscosity coefficients of the one-component plasma are tested against molecular dynamics simulations. The theory is found to accurately capture the kinetic contributions to the transport coefficients, but not the potential contributions that arise at very strong coupling (Γ≳30). Considerations related to a first-principles generalization of Enskog's kinetic equation to continuous potentials are also discussed.