No Arabic abstract
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.
We present an accurate computational study of the electronic structure and lattice dynamics of solid molecular hydrogen at high pressure. The band-gap energies of the $C2/c$, $Pc$, and $P6_3/m$ structures at pressures of 250, 300, and 350 GPa are calculated using the diffusion quantum Monte Carlo (DMC) method. The atomic configurations are obtained from ab-initio path-integral molecular dynamics (PIMD) simulations at 300 K and 300 GPa to investigate the impact of zero-point energy and temperature-induced motion of the protons including anharmonic effects. We find that finite temperature and nuclear quantum effects reduce the band-gaps substantially, leading to metallization of the $C2/c$ and $Pc$ phases via band overlap; the effect on the band-gap of the $P6_3/m$ structure is less pronounced. Our combined DMC-PIMD simulations predict that there are no excitonic or quasiparticle energy gaps for the $C2/c$ and $Pc$ phases at 300 GPa and 300 K. Our results also indicate a strong correlation between the band-gap energy and vibron modes. This strong coupling induces a band-gap reduction of more than 2.46 eV in high-pressure solid molecular hydrogen. Comparing our DMC-PIMD with experimental results available, we conclude that none of the structures proposed is a good candidate for phases III and IV of solid hydrogen.
The dynamic electron-ion collisions play an important rolein determining the static and transport properties of warmdense matter (WDM). Electron force field (eFF) method is applied to study the ionic transport properties of warm densehydrogen. Compared with the results from quantum moleculardynamics and orbital-free molecular dynamics, the ionicdiffusions are largely reduced by involving the dynamic collisions of electrons and ions. This physics is verified by quantum Langevin molecular dynamics (QLMD) simulations, which includes electron-ion collisions induced friction(EI-CIF) into the dynamic equation of ions. Based on these new results, we proposed a model including the correctionof collisions induced friction of the ionic diffusion. The CIF model has been verified to be valid at a wide range ofdensity and temperature. We also compare the results with the Yukawa one component plasma (YOCP) model andEffective OCP (EOCP) model. We proposed to calculate the self-diffusion coefficients using the EOCP model modifiedby the CIF model to introduce the dynamic electron-ion collisions effect.
The laws of quantum mechanics are often tested against the behaviour of the lightest element in the periodic table, hydrogen. One of the most striking properties of molecular hydrogen is the coupling between molecular rotational properties and nuclear spin orientations, giving rise to the spin isomers ortho- and para-hydrogen. At high pressure, as intermolecular interactions increase significantly, the free rotation of H2 molecules is increasingly hindered, and consequently a modification of the coupling between molecular rotational properties and the nuclear spin system can be anticipated. To date, high-pressure experimental methods have not been able to observe nuclear spin states at pressures approaching 100 GPa and consequently the effect of high pressure on the nuclear spin statistics could not be directly measured. Here, we present in-situ high-pressure nuclear magnetic resonance data on molecular hydrogen in its hexagonal phase I up to 123 GPa at room temperature. While our measurements confirm the presence of I=1 ortho-hydrogen at low pressures, above 70 GPa, where inter- and intramolecular distances become comparable, we observe a crossover in the nuclear spin statistics from a spin-1 quadrupolar to a spin-1/2 dipolar system, evidencing the loss of spin isomer distinction. These observations represent a unique case of a nuclear spin crossover phenomenon in quantum solids.
We study the thermophysical properties of warm dense hydrogen using quantum molecular dynamics simulations. New results are presented for the pair distribution functions, the equation of state, the Hugoniot curve, and the reflectivity. We compare with available experimental data and predictions of the chemical picture. Especially, we discuss the nonmetal-to-metal transition which occurs at about 40 GPa in the dense fluid.
Hydrogen bond symmetrisations in H-bonded systems triggered by pressure induced nuclear quantum effects (NQEs) is a long-known concept1 but experimental evidences in high-pressure ices have remained elusive with conventional methods2,3. Theoretical works predicted quantum-mechanical tunneling of protons within water ices to occur at pressures above 30 GPa and the H-bond symmetrisation transition above 60 GPa4. Here, we used 1H-NMR on high-pressure ice up to 90 GPa, and demonstrate that NQEs govern the behavior of the hydrogen bonded protons in ice VII already at significantly lower pressures than previously expected. A pronounced tunneling mode was found to be present up to the highest pressures of 90 GPa, well into the stability field of ice X, where NQEs are not anticipated in a fully symmetrized H-bond network. We found two distinct transitions in the NMR shift data at about 20 GPa and 75 GPa attributed to the step-wise symmetrization of the H-bond (HB), with high-barrier H-Bonds (HBHB) to low-barrier H-bonds (LBHB) and LBHB to symmetric H-bonds (SHB) respectively. These transitions could have major implication on the physical properties of high-pressure ices and planetary interior models. NQEs observed in this chemically simple system over a wide pressure range could prove to be useful in designing a new generation of electronic devices exploiting protonic tunneling.