Search

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.

Shock reverberation compression experiments on dense gaseous deuterium-helium mixtures are carried out to provide thermodynamic parameters relevant to the conditions in planetary interiors. The multishock pressures are determined up to 120 GPa and reshock temperatures to 7400 K. Furthermore, the unique compression path from shock-adiabatic to quasi-isentropic compressions enables a direct estimation of the high-pressure sound velocities in the unexplored range of 50-120 GPa. The equation of state and sound velocity provide particular dual perspectives to validate the theoretical models. Our experimental data are found to agree with several equation of state models widely used in astrophysics within the probed pressure range. The current data improve the experimental constraints on sound velocities in the Jovian insulating-to-metallic transition layer.

The transport properties of warm dense neon (Ne) and krypton (Kr) are studied by combining self-consistent fluid variational theory (SFVT) with linear response theory (LRT). The components are determined using the SFVT, and the transport parameters, including the electrical conductivity, thermal conductivity, and thermopower, are calculated with the LRT. The relevant scattering mechanisms, including electron-ion, electron-electron, and electron-atom scatterings, are taken into account. An effective potential model in combination with the Muffin-tin model is introduced to further improve the description for electron-atom scattering, which not only includes static, exchange, and polarization interactions but also considers the plasma environmental effects. It is found that for electron-atom scattering, the influence of the plasma density is significant at lower scattering energies but the effects are different for electron-Ne and electron-Kr scattering. For electron-Kr scattering, a plasma density-dependent Ramsauer-Townsend minimum is observed. The obtained transport parameters are compared with the available experiments and other simulations. The plasma phase transition of warm dense Kr is revisited from multiple perspectives based on the numerical simulation results for the electrical conductivity and thermopower. These observations may help one to better understand the transport properties of warm dense noble gases and are an important guide for future experimental designs and theoretical developments.

A dynamic mitigation is presented for sausage and kink instability growths of a z-current-driven magnetized plasma column. In this Rapid Communication we found that a wobbling motion of the z-current electron axis induces a phase-controlled perturbation, so that the growths of the sausage and kink instabilities are successfully and remarkably mitigated. In general, plasma instabilities emerge from perturbations, and the perturbation phase is normally unknown. However, if the perturbation phase is known or actively imposed by, for example, a designed driver wobbling behavior, the instability growth would be controlled and mitigated by a superimposition of the perturbations imposed. The results in this Rapid Communication demonstrate that the wobbling z-current electron beam would provide an improvement in the plasma column stability and uniformity.

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.

Extensive quantum molecular dynamics (QMD) simulations are performed to determine the equation of state, sound velocity, and phase diagram of middle-Z krypton in a warm dense regime where the pressure (P) is up to 300 GPa and the temperature is up to 60 kK. The shock wave experimental data are used to validate the present theoretical models. It is found that, within the regime of the current density (ρ) and temperature (T), sound velocity can effectively discriminate differences between different theoretical models, and therefore it is more suitable as a benchmark to verify the practicability of models. The QMD-simulated results of the ionic structures and electronic properties imply the occurrence of two kinds of phase transitions, including transition from a solidlike to fluid state and that from an insulator to conductive fluid in this T-P regime. The calculated electrical conductivities confirm that the metallization transition occurs at about 60 GPa and 17.5 kK along the principal Hugoniot. With the help of simulation results and experimental data, a comprehensive phase diagram for krypton is constructed by using the solid-fluid and insulator-metal fluid phase boundaries, which fills the gap of the experimental work [Proc. Natl. Acad. Sci. USA 112, 7925 (2015)PNASA60027-842410.1073/pnas.1421801112]. These results will provide an instructive basis for the experimental investigations of rare gases over a wide T-P range.

Peculiar behavior is found in a Si cluster ion, moving with a speed of ∼0.22c (c: speed of light) in a solid Al plasma in the context of cluster-ion beam driven inertial confinement fusion: The Si ion, moving closely behind the forward-moving Si ion for a distance of several Å in the cluster, feels the wake field generated by the forward Si. The acceleration interaction force on the rear Si ion by the forward-moving ion may balance with the deceleration backward force in the longitudinal-moving direction. The forward-moving Si ion would be decelerated, as an isolated ion is decelerated without correlation. However, the deceleration of the rear Si ion, moving behind closely, would be reduced significantly. The rear Si ion may catch up and overtake the forward-moving Si ion in the cluster during the Si cluster-ion interaction with the high-density Al plasma. This peculiar behavior appears when the ions are aligned well longitudinally. The wake field is confined around the Si ion in the forward and transverse directions for a distance smaller than the Si cluster interionic distance l_{c}. However, the tail of the wake field extends beyond l_{c} due to the Si ion high speed of ∼0.22c. Therefore, the peculiar behavior shown above appears only for the ions in one cluster aligned well longitudinally.

The magnetic quadrupole structure formation during the interaction of two ultrashort high power laser pulses with a collisionless plasma is demonstrated with 2.5-dimensional particle-in-cell simulations. The subsequent expansion of the quadrupole is accompanied by magnetic-field annihilation in the ultrarelativistic regime, when the magnetic field cannot be sustained by the plasma current. This results in a dominant contribution of the displacement current exciting a strong large scale electric field. This field leads to the conversion of magnetic energy into kinetic energy of accelerated electrons inside the thin current sheet.

A consistent theoretical model that can be applied in a wide range of densities and temperatures is necessary for understanding the variation of a material's properties during compression and heating. Taking argon as an example, we show that the combination of self-consistent fluid variational theory and linear response theory is a promising route for studying warm dense matter. Following this route, the compositions, equations of state, and transport properties of argon plasma are calculated in a wide range of densities (0.001-20 g/cm(3)) and temperatures (5-100 kK). The obtained equations of state and electrical conductivities are found in good agreement with available experimental data. The plasma phase transition of argon is observed at temperatures below 30 kK and density about 2-6g/cm(3). The minimum density for the metallization of argon is found to be about 5.8 g/cm(3), occurring at 30-40 kK. The effects of many-particle correlations and dynamic screening on the electrical conductivity are also discussed through the effective potentials.

The current partially ionized plasmas models for xenon show substantial differences since the description of pressure and thermal ionization region becomes a formidable task, prompting the need for an improved understanding of dense xenon plasmas behavior at above 100 GPa. We performed double-shock compression experiments on dense xenon to determine accurately the Hugoniot up to 172 GPa using a time-resolved optical radiation method. The planar strong shock wave was produced using a flyer plate impactor accelerated up to ∼6 km/s with a two-stage light-gas gun. The time-resolved optical radiation histories were acquired by using a multiwavelength channel optical transience radiance pyrometer. Shock velocity was measured and mass velocity was determined by the impedance-matching methods. The experimental equation of state of dense xenon plasmas are compared with the self-consistent fluid variational calculations of dense xenon in the region of partial ionization over a wide range of pressures and temperatures.

Experimental equations of state on generation of nonideal xenon plasma by intense shock wave compression was presented in the ranges of pressure of 2-16 GPa and temperature of 31-50 kK, and the xenon plasma with the nonideal coupling parameter Γ range from 0.6-2.1 was generated. The shock wave was produced using the flyer plate impact and accelerated up to ∼6 km/s with a two-stage light gas gun. Gaseous specimens were shocked from two initial pressures of 0.80 and 4.72 MPa at room temperature. Time-resolved spectral radiation histories were recorded by using a multiwavelength channel pyrometer. The transient spectra with the wavelength range of 460-700 nm were recorded by using a spectrometer to evaluate the shock temperature. Shock velocity was measured and particle velocity was determined by the impedance matching methods. The equations of state of xenon plasma and ionization degree have been discussed in terms of the self-consistent fluid variational theory.

The ionization degree and equation of state of dense xenon plasma were calculated by using self-consistent fluid variational theory for temperature of 4-30kK and density of 0.01-8.5gcm;{3} . The dense fluid xenon will be ionized at high pressures and temperatures. The ionization energy of xenon will be lowered due to the interactions among all particles of Xe, Xe+ , Xe2+ , and e . The ionization degree is obtained from nonideal ionization equilibrium, taking into account the correlative contributions to the chemical potential which is determined self-consistently by the free energy function. The composition of xenon has been calculated with given densities and temperatures in the region of partial ionization. The calculated results show a pressure softening regime at the onset of ionization. Comparison is performed with available shock-wave experiments and other theoretical calculations.