Centroid molecular dynamics approach to the transport properties of liquid para-hydrogen over the wide temperature range

Self-diffusion of liquid deuterium hydride and liquid tritium

We present a quasi-elastic neutron scattering study of liquid deuterium hydride carried out using the Disk Chopper Spectrometer at the National Institute of Standards and Technology. Under saturated vapor pressure, the self-diffusion constant of deuterium hydride obeys an Arrhenius law D=D0exp-EA/kBT, where the prefactor D0 is given by D0=9.5±1.2Å2/ps and the activation energy is given by EA = 58 ± 2 K. We apply the quantum law of corresponding states to the known diffusion constants of the hydrogen isotopologues. From this application, we estimate that D0≈9.1Å2/ps and EA ≈ 75 K in liquid tritium. Young's theory of quantum-mechanical effects in van der Waals fluids [Young, Phys. Rev. A 23, 1498 (1981)] is shown to apply to the diffusion constants of the liquid hydrogens. Our results underscore the importance of nuclear quantum effects in shaping the properties and behavior of the hydrogen isotopologues.

Distinct molecular dynamics dividing liquid-like and gas-like supercritical hydrogens

Understanding how a supercritical fluid is related to normal liquid and gas and separating it into liquid-like and gas-like regions are of fundamental and practical importance. Despite the usefulness of hydrogen storage, molecular dynamics images on supercritical hydrogens exhibiting strong nuclear quantum effects are scarce. Taking advantage of the non-empirical ab initio molecular dynamics method for hydrogen molecules, we found that, while radial distribution functions and diffusion show a monotonic change along the density, van Hove time correlation functions and intramolecular properties such as bond length and vibrational frequency exhibit the anomalous order crossing the Widom line. By demonstrating that the anomalous order stemmed from the largest deviations between liquid-like and gas-like solvations formed around the Widom line, we concluded that this supercritical fluid is a mixture of liquid and gas possessing heterogeneity. The obtained physical insights can be an index to monitor the supercriticality and to identify distinct liquid-like and gas-like supercritical fluids.

Quantum polyamorphism in compressed distinguishable helium-4

We demonstrate that two amorphous solid states can exist in ⁴He consisting of distinguishable Boltzmann atoms under compressed conditions. The isothermal compression of normal or supercritical fluid ⁴He was conducted at 3-25 K using the isobaric-isothermal path integral centroid molecular dynamics simulation. The compression of fluid first produced the low-dispersion amorphous (LDA) state possessing modest extension of atomic necklaces. Further isothermal compression up to the order of 10 kbar to 1 Mbar or an isobaric cooling of LDA induced the transition to the high-dispersion amorphous (HDA) state. The HDA was characterized by long quantum wavelengths of atoms extended over several Angstroms and the promotion of atomic residual diffusion. They were related to the quantum tunneling of atoms bestriding the potential saddle points in this glass. The change in pressure or temperature induced the LDA-HDA transition reversibly with hysteresis, while it resembled the coil-globule transition of classical polymers. The HDA had lower kinetic and higher Gibbs free energies than the LDA at close temperature. The HDA was absent at T ≥ 13 K, while the LDA-HDA transition pressure significantly decreased with lowering temperature. The LDA and HDA correspond to the trapped and tunneling regimes proposed by Markland et al. [J. Chem. Phys. 136, 074511 (2012)], respectively. The same reentrant behavior as they found was observed for the expansion factor of the quantum wavelength as well as for atomic diffusivity.

Nuclear quantum effects in thermal conductivity from centroid molecular dynamics

We show that the centroid molecular dynamics (CMD) method provides a realistic way to calculate the thermal diffusivity a = λ/ρc_V of a quantum mechanical liquid such as para-hydrogen. Once a has been calculated, the thermal conductivity can be obtained from λ = ρc_Va, where ρ is the density of the liquid and c_V is the constant-volume heat capacity. The use of this formula requires an accurate quantum mechanical heat capacity c_V, which can be obtained from a path integral molecular dynamics simulation. The thermal diffusivity can be calculated either from the decay of the equilibrium density fluctuations in the liquid or by using the Green-Kubo relation to calculate the CMD approximation to λ and then dividing this by the corresponding approximation to ρc_V. We show that both approaches give the same results for liquid para-hydrogen and that these results are in good agreement with the experimental measurements of the thermal conductivity over a wide temperature range. In particular, they correctly predict a decrease in the thermal conductivity at low temperatures-an effect that stems from the decrease in the quantum mechanical heat capacity and has eluded previous para-hydrogen simulations. We also show that the method gives equally good agreement with the experimental measurements for the thermal conductivity of normal liquid helium.

Decelerated Liquid Dynamics Induced by Component-Dependent Supercooling in Hydrogen and Deuterium Quantum Mixtures

Isotopic mixtures of p-H₂ and o-D₂ molecules have been an attractive binary system because they include two kinds of purely isotopic molecules which possess the same electronic potential but the twice different mass inducing differently pronounced nuclear quantum effects (NQEs). Accessing details of structures and dynamics in such quantum mixtures combining complex molecular dynamics with NQEs of different strengths remains a challenging problem. Taking advantage of the nonempirical molecular dynamics method which describes p-H₂ and o-D₂ molecules, we found that the liquid dynamics slows down at a specific mixing ratio, which can be connected to the observed anomalous slowdown of crystallization in the quantum mixtures. We attributed the decelerated dynamics to the component-dependent supercooling of p-H₂ taking place in the mixtures, demonstrating that there is an optimal mixing ratio to hinder crystallization. The obtained physical insights will help in experimentally controlling and achieving unknown quantum mixtures including superfluid.

Isotopic Effects on Intermolecular and Intramolecular Structure and Dynamics in Hydrogen, Deuterium, and Tritium Liquids: Normal Liquid and Weakly and Strongly Cooled Liquids

Differences in properties such as phase-transition temperature and transport coefficients among liquids of different isotopic compositions, hydrogen, deuterium, and tritium, should originate from their differently pronounced nuclear quantum effects (NQEs) rather than from any subtle difference in the electronic interaction potentials. Accurate and efficient determination of structural and dynamical isotopic effects in the quantum liquids still remains as one of the challenging problems in condensed-phase physics. With a recently developed nonempirical real-time molecular dynamics method which describes nonspherical molecules with the NQEs, we computationally realized and investigated dynamical and quantum isotopic effects of not only traditionally studied isotopes, hydrogen, and deuterium but also a lesser known radioisotope, tritium, in broad thermodynamic conditions from normal liquid to weakly and strongly cooled liquids, which have been hindered by rapid crystallization in spite of numerous experimental attempts at supercooling. Reproducing the previously reported experimental isotope dependence on the bond length and vibrational frequencies of hydrogen, deuterium, and tritium liquids, we further demonstrate that distinctive isotope effects appear in their intermolecular and intramolecular structure and dynamics not only at lower temperature but also at higher temperature, which none has so far been able to obtain quantitative results for realistic systems. Rationalization of their physical origins and the obtained physical insights will help future experimental searching and monitoring intermolecular and intramolecular dynamics and structures of these isotopes not only in normal liquid but also in supercooled liquid.

Dynamical Ordering of Hydrogen Molecules Induced by Heat Flux

Achieving a direct nonequilibrium simulation for hydrogen systems has been quite challenging because nuclear quantum effects (NQEs) have to be taken into account. We directly simulated nonequilibrium hydrogen molecules under a temperature gradient with the recently developed nonempirical molecular dynamics method, which describes nonspherical hydrogen molecules with the NQEs. We found dynamical ordering purely induced by heat flux, which should be distinguished from static ordering like orientational alignment, as decelerated translational motions and enhanced intensity of H-H vibrational power spectra despite the little structural ordering. This dynamical ordering, which was enhanced with stronger heat flux while independent of system size, can be regarded as self-solidification of hydrogen molecules for their efficient heat conduction.

A molecular dynamics study of nuclear quantum effect on diffusivity of hydrogen molecule

In this paper, the nuclear quantum effect of the hydrogen molecule on its diffusivity was analyzed using the molecular dynamics (MD) method. The centroid MD (CMD) method was applied to reproduce the time evolution of the molecules. The diffusion coefficient of hydrogen was calculated using the Green-Kubo method over a wide temperature region, and the temperature dependence of the quantum effect of the hydrogen molecule on its diffusivity was addressed. The calculated results were compared with classical MD results based on the principle of corresponding state (PCS). It was confirmed that the difference in the diffusion coefficient calculated in the CMD and classical MD methods was small, and the PCS appears to be satisfied on the temperature dependence of the diffusion coefficient, even though the quantum effect of the hydrogen molecules was taken into account. It was clarified that this result did not suggest that the quantum effect on the diffusivity of the hydrogen molecule was small but that the two changes in the intermolecular interaction of hydrogen due to the quantum effect offset each other. Moreover, it was found that this tendency was related to the temperature dependence of the ratio of the kinetic energy of the quantum fluctuational motion to the classical kinetic energy.

Distinct structural and dynamical difference between supercooled and normal liquids of hydrogen molecules

Supercooled hydrogen liquid as well as superfluid have continued to elude experimental observation due to rapid crystallization. We computationally realized and investigated supercooled hydrogen liquid by a recently developed non-empirical real-time molecular dynamics method, which describes non-spherical hydrogen molecules with the nuclear quantum effects. We demonstrated that the hydrogen supercooled liquid is not a simply cooled liquid but rather exhibits intrinsic structural and dynamical characters including a precursor of tunneling and superfluidity which neither normal hydrogen liquid nor solid possesses. All of the insights provide a milestone for planning experiments of metastable hydrogen systems like glassy and superfluid states and for identifying various unknown hydrogen phases.

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 planet's 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/cm(3), 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.

Inclusion of inversion symmetry in centroid molecular dynamics: a possible avenue to recover quantum coherence

Inversion symmetry is included in the operator formulation of the centroid molecular dynamics (CMD). This work involves the development of a symmetry-adapted CMD (SA-CMD), here particularly for symmetrization and antisymmetrization projections. A symmetry-adapted quasidensity operator, as defined by Blinov and Roy [J. Chem. Phys. 115, 7822 (2001)], is employed to obtain the centroid representation of quantum mechanical operators. Numerical examples are given for a single particle confined to one-dimensional symmetric quartic and symmetric double-well potentials. Two SA-CMD simulations are performed separately for both projections, and centroid position autocorrelation functions are obtained. For each projection, the quality of the approximation as well as the accuracy are similar to those of regular CMD. It is shown that individual trajectories from two separate SA-CMD simulations can be properly combined to recover trajectories for Boltzmann statistics. Position autocorrelation functions are compared to the exact quantum mechanical ones. This explicit account of inversion symmetry provides a qualitative improvement on the conventional CMD approach and allows the recovery of some quantum coherence.

Symmetrized correlation function for liquidpara-hydrogen using complex-time pair-product propagators

We present a simple and efficient method for calculating symmetrized time correlation functions of neat quantum fluids. Using the pair-product approximation to each complex-time quantum mechanical propagator, symmetrized correlation functions are written in terms of a double integral for each degree of freedom with a purely positive integrand. At moderate temperatures and densities, where the pair-product approximation to the Boltzmann operator is sufficiently accurate, the method leads to quantitative results for the early time part of the correlation function. The method is tested extensively on liquid para-hydrogen at 25 K and used to obtain accurate quantum mechanical results for the initial 0.2 ps segment of the symmetrized velocity autocorrelation function of this system, as well as the incoherent dynamic structure factor at certain momentum transfer values.

Fast centroid molecular dynamics: a force-matching approach for the predetermination of the effective centroid forces

A fast centroid molecular dynamics (CMD) methodology is proposed in which the effective centroid forces are predetermined through a force-matching algorithm applied to a standard path integral molecular dynamics simulation. The resulting method greatly reduces the computational cost of generating centroid trajectories, thus extending the applicability of CMD. The method is applied to the study of liquid para-hydrogen at two state points and liquid ortho-deuterium at one state point. The static and dynamical results are compared to those obtained from full adiabatic CMD simulations and found to be in excellent agreement for all three systems; the transport properties are also compared to experiment and found to have a similar level of agreement.

A comparative study of imaginary time path integral based methods for quantum dynamics

The recently introduced approximate many-body quantum simulation method, ring polymer molecular dynamics (RPMD), is compared to the centroid molecular dynamics method (CMD). Comparisons of simulation results for liquid para-hydrogen at two state points and liquid ortho-deuterium at one state point are presented. The calculated quantum correlation functions for the two methods are shown to be in good agreement with one another for a large portion of the time spectrum. However, as the quantum mechanical nature of the system increases, RPMD is less accurate in predicting the kinetic energy of the system than is CMD. A simplified and highly efficient algorithm is proposed which largely corrects this deficiency.

QUANTUM MODE-COUPLING THEORY: Formulation and Applications to Normal and Supercooled Quantum Liquids

▪ Abstract  We review our recent efforts to formulate and study a mode-coupling approach to real-time dynamic fluctuations in quantum liquids. Comparison is made between the theory and recent neutron scattering experiments performed on liquid ortho-deuterium and para-hydrogen. We discuss extensions of the theory to supercooled and glassy states where quantum fluctuations compete with thermal fluctuations. Experimental scenarios for quantum glassy liquids are briefly discussed.

Effective potential analytic continuation approach for real time quantum correlation functions involving nonlinear operators

We apply the effective potential analytic continuation (EPAC) method to the calculation of real time quantum correlation functions involving operators nonlinear in the position operator q. For a harmonic system the EPAC method provides the exact correlation function at all temperature ranges, while the other quantum dynamics methods, the centroid molecular dynamics and the ring polymer molecular dynamics, become worse at lower temperature. For an asymmetric anharmonic system, the EPAC correlation function is in very good agreement with the exact one at t = 0. When the time increases from zero, the EPAC method gives good coincidence with the exact result at lower temperature. Finally, we propose a simplified version of the EPAC method to reduce the computational cost required for the calculation of the standard effective potential.

A centroid molecular dynamics study of liquidpara-hydrogen andortho-deuterium

Centroid molecular dynamics (CMD) is applied to the study of collective and single-particle dynamics in liquid para-hydrogen at two state points and liquid ortho-deuterium at one state point. The CMD results are compared with the results of classical molecular dynamics, quantum mode coupling theory, a maximum entropy analytic continuation approach, pair-product forward- backward semiclassical dynamics, and available experimental results. The self-diffusion constants are in excellent agreement with the experimental measurements for all systems studied. Furthermore, it is shown that the method is able to adequately describe both the single-particle and collective dynamics of quantum liquids.