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

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.

Phase Separation in Cold Para-H_{2} D_{2} Clusters

Low temperature phase separation in mixtures of ^{3}He and ^{4}He isotopes is a unique property of quantum fluids. Hydrogen has long been considered as another potential quantum liquid and has been predicted to be superfluid at T≤1  K, well below freezing temperature of ≈14  K. Phase separation has also been predicted in mixtures of para-H_{2} and D_{2} at temperatures ≤3  K. To defer the freezing, we produced clusters containing para-H_{2} and D_{2} at an estimated temperature of ≈2  K whose state was studied by vibrational Raman spectroscopy. The results indicate that the clusters are liquid and show the phase separation of the isotopes. The phase separation is further corroborated by the quantum molecular dynamics simulation.

Anomalously supercooled H2-D2 mixtures flowing inside a carbon nano tube

H₂ and D₂ molecules condensed in a carbon nano tube (CNT) and their nonequilibrium flow through nano pores offer a key test to reveal mass molecular transport and separation of purely isotopic molecules that possess the same electronic potential but a two-times difference in mass inducing differently enhanced nuclear quantum effects (NQEs) such as nuclear delocalization and zero-point energy. Taking advantage of the non-empirical quantum molecular dynamics method developed for condensed H₂-D₂ molecules that can describe various kinds of condensed phases and thermodynamic states including uneven density and a shear flow, we investigated condensed isotopic H₂-D₂ mixtures flowing inside nanoscale adsorbable CNTs. We found that, in any mixture, the more delocalized H₂ molecules are more supercooled than the less delocalized D₂ molecules in a two-dimensional liquid film adsorbed around the CNT well, and that the stronger supercooling of the H₂ molecules than the D₂ molecules in an equilibrium state becomes more enhanced under the nonequilibrium flow due to the isotope-dependent flow-induced condensation, demonstrating the anomalous condensed-phase quantum sieving under the nonequilibrium flow and its dependence on the mixing ratio and temperature. The differently enhanced NQEs of the purely isotopic molecules essentially influence the condensed adsorption and their flows occurring in the nanoscale CNT, which should be distinguished from a dilute gas adsorption. The predicted properties and obtained physical insights in this paper will help in experimentally controlling condensed H₂-D₂ mixtures, and open a new strategy and innovative design of nanoporous materials for adsorptive separation of condensed-phase mixtures under a nonequilibrium flow not of a dilute gas mixture in an equilibrium state.

Real-time hydrogen molecular dynamics satisfying the nuclear spin statistics of a quantum rotor

Apparent presence of the nuclear-spin species of a hydrogen molecule, para-hydrogen and ortho-hydrogen, associated with the quantum rotation is a manifestation of the nuclear quantum nature of hydrogen, governing not only molecular structures but also physical and chemical properties of hydrogen molecules. It has been a great challenge to observe and calculate real-time dynamics of such molecularized fermions. Here, we developed the non-empirical quantum molecular dynamics method that enables real-time molecular dynamics simulations of hydrogen molecules satisfying the nuclear spin statistics of the quantum rotor. While reproducing the species-dependent quantum rotational energy, population ratio, specific heat, and H-H bond length and frequency, we found that their translational, orientational and vibrational dynamics becomes accelerated with the higher rotational excitation, concluding that the nuclear quantum rotation stemmed from the nuclear spin statistics can induce various kinds of dynamics and reactions intrinsic to each hydrogen species.

Flow-Induced Autonomic Ordering of Hydrogen Molecules under a Non-Equilibrium Flow

A non-equilibrium molecular flow through a carbon nanotube (CNT) serves as a key system for revealing molecular transport and establishing nanofluidics. It has been challenging to simulate a non-equilibrium flow of hydrogen molecules exhibiting strong nuclear quantumness. Taking advantage of the quantum molecular dynamics method that can calculate real-time trajectories of hydrogen molecules even under a non-equilibrium flow, we found that the non-equilibrium flow makes hydrogen molecules more condensed and accelerates their adsorption near a CNT surface, letting the molecules flow more smoothly by propagating velocity momenta more efficiently along the CNT axis and by suppressing transverse molecular dynamics on the CNT cross section. Such flow-induced autonomic ordering indicates the importance of monitoring and investigating dynamics and adsorption of hydrogen molecules under a non-equilibrium circumstance as well as in a quiet equilibrium state, opening a new strategy for efficient hydrogen liquefaction and storage.

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.

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.