Huayang Sun

Structure, equation of state, diffusion and viscosity of warm dense Fe under the conditions of a giant planet core

Fe exists abundantly in the universe. In particular, the dynamic structures and transport properties of warm dense Fe are crucial for understanding the evolution and structures of giant planets. In this paper, we present the ionic structures, equation of states, diffusion and viscosity of Fe at two typical densities of 33.385 and 45?g?cm?3 in the temperature range of 1?10?eV, giving the data by the first principles calculations using quantum Langevin molecular dynamics. Furthermore, the validation of the Stokes?Einstein (SE) relation in this regime is discussed, showing the importance of choosing the effective atomic diameter. The results remind us of the careful usage of the SE relation under extreme conditions.

Quantum simulation of thermally-driven phase transition and oxygen K-edge x-ray absorption of high-pressure ice

The structure and phase transition of high-pressure ice are of long-standing interest and challenge, and there is still a huge gap between theoretical and experimental understanding. The quantum nature of protons such as delocalization, quantum tunneling and zero-point motion is crucial to the comprehension of the properties of high-pressure ice. Here we investigated the temperature-induced phase transition and oxygen K-edge x-ray absorption spectra of ice VII, VIII and X using ab initio path-integral molecular dynamics simulations. The tremendous difference between experiments and the previous theoretical predictions is closed for the phase diagram of ice below 300 K at pressures up to 110 GPa. Proton tunneling assists the proton-ordered ice VIII to transform into proton-disordered ice VII where only thermal activated proton-transfer cannot occur. The oxygen K edge with its shift is sensitive to the order-disorder transition, and therefore can be applied to diagnose the dynamics of ice structures.

Nuclear quantum dynamics in dense hydrogen

Nuclear dynamics in dense hydrogen, which is determined by the key physics of large-angle scattering or many-body collisions between particles, is crucial for the dynamics of planets evolution and hydrodynamical processes in inertial confinement confusion. Here, using improved ab initio path-integral molecular dynamics simulations, we investigated the nuclear quantum dynamics regarding transport behaviors of dense hydrogen up to the temperatures of 1 eV. With the inclusion of nuclear quantum effects (NQEs), the ionic diffusions are largely higher than the classical treatment by the magnitude from 20% to 146% as the temperature is decreased from 1 eV to 0.3 eV at 10 g/cm3, meanwhile, electrical and thermal conductivities are significantly lowered. In particular, the ionic diffusion is found much larger than that without NQEs even when both the ionic distributions are the same at 1 eV. The significant quantum delocalization of ions introduces remarkably different scattering cross section between protons compared with classical particle treatments, which explains the large difference of transport properties induced by NQEs. The Stokes-Einstein relation, Wiedemann-Franz law, and isotope effects are re-examined, showing different behaviors in nuclear quantum dynamics.Nuclear quantum effects (NQEs) on the structures and transport properties of dense liquid hydrogen at densities of 10-100 g/cm3 and temperatures of 0.1-1 eV are fully assessed using \textit{ab initio} path-integral molecular dynamics simulations. With the inclusion of NQEs, ionic diffusions are strongly enhanced by the magnitude from 100% to 15% with increasing temperature, while electrical conductivities are significantly suppressed. The analyses of ionic structures and zero-point energy show also the importance of NQEs in these regime. The significant quantum delocalization of ions introduces expressively different scattering cross section between protons compared with classical particle treatments, which can explain the large alterability of transport behaviors. Furthermore, the energy, pressure, and isotope effects are also greatly influenced by NQEs. The complex behaviors show that NQEs can not be neglected for dense hydrogen even in the warm dense regime.

First-principles study on equation of states and electronic structures of shock compressed Ar up to warm dense regime

The equation of states (EOS) and electronic structures of argon with temperatures from 0.02 eV to 3 eV and densities from 0.5 g/cm(3) to 5.5 g/cm(3) are calculated using the pair potential and many-body potential molecular dynamics and the density functional theory (DFT) molecular dynamics with van der Waals (vdW) corrections. First-principles molecular dynamics is implemented above 2.0 g/cm(3). For the cases of low densities below 3 g/cm(3), we performed pair potential molecular dynamics in order to obtain the ionic configurations, which are used in density functional theory to calculate the EOS and electronic structures. We checked the validity of different methods at different densities and temperatures, showing their behaviors by comparing EOS. DFT without vdW correction works well above 1 eV and 3.5 g/cm(3). Below 1 eV and 2.0 g/cm(3), it overestimates the pressure apparently and results in incorrect behaviors of the internal energy. With vdW corrections, the semi-empirical force-field correction (DFT-D2) method gives consistent results in the whole density and temperature region, and the vdW density functional (vdW-DF2) method gives good results below 2.5 g/cm(3), but it overestimates the pressure at higher densities. The interactions among the atoms are overestimated by the pair potential above 1 eV, and a temperature dependent scaled pair potential can be used to correct the ionic configurations of the pair potential up to 3 eV. The comparisons between our calculations and the experimental multi-shock compression results show that the Hugoniot line of DFT-D2 and DFT tends to give larger pressure than the results of the self-consistent fluid variational theory, and the difference increases with the density. The electronic energy gap exists for all our cases up to 5.5 g/cm(3) and 1 eV. The effect of vdW interactions on the electronic structures are also discussed.