Communication: Quantum molecular dynamics simulation of liquid para-hydrogen by nuclear and electron wave packet approach

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.

Computation of static quantum triplet structure factors of liquid para-hydrogen

The instantaneous and centroid triplet structure factors, S ( 3 ) ( k 1 , k 2 ) , of liquid (one-center) para-hydrogen are computed on the crystallization line for temperatures T/K ≤ 33. The focus is on salient equilateral and isosceles features, and the methods utilized are path integral Monte Carlo (PIMC) simulations and Ornstein-Zernike (OZ) integral equations, which involve Jackson-Feenberg convolution (JF3) and other distinct closures. Long path integral simulation runs are carried out in the canonical ensemble, so as to obtain sufficiently accurate direct PI triplet results. Conclusions are drawn regarding general triplet structure features and the role and usefulness of the OZ closures. The equilateral features are studied in more detail, and one finds that (a) PIMC results point to the existence of regularity in the centroid main peak amplitudes; (b) some of the studied closures give qualitative descriptions for wave numbers below k ≈ 1 Å^-1, but they all fail to describe the main peak amplitude regions (1.75 < k/Å^-1 < 2.5); and (c) JF3 plays the role of a limit closure that is valid for increasing wave numbers (k ≥ 2.6 Å^-1). In addition, representative isosceles PI features turn out to be reasonably bounded (within Δk = 0.1 Å^-1) by those of some closures.

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 corpuscular picture of electrons in chemical bond

We introduce a theory of chemical bond with a corpuscular picture of electrons. It employs a minimal set of localized electron wave packets with "floating and breathing" degrees of freedom and the spin-coupling of non-orthogonal valence-bond theory. Its accuracy for describing potential energy curves of chemical bonds in ground and excited states of spin singlet and triplet is examined.

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.