Simulating Asymmetric Top Impurities in Superfluid Clusters: A para-Water Dopant in para-Hydrogen

Operator Formulation of Feynman Path Centroid Dynamics for Rotations

An operator formulation of centroid molecular dynamics (CMD) for rotational degrees of freedom is presented. The quasi-density operator concept was introduced by Jang and Voth [J. Chem. Phys 111, 2357 (1999)] and is used to obtain a phase-space mapping without the need for discretized path integrals. The approach allows the calculation of approximate Kubo-transformed time correlation functions. The particle on a ring is chosen as an illustrative example. Numerical results demonstrate that the proposed approach leads to accurate results when compared with exact diagonalization calculations for linear operators. At very low temperatures, it is found that rotational CMD yields results that are in exact agreement with the quantum dynamics of a spin-1 system.

An intramolecular vibrationally excited intermolecular potential energy surface and predicted 2OH overtone spectroscopy of H2O-Kr

A new five-dimensional potential energy surface (PES) for H2O-Kr which explicitly includes the intramolecular 2OH overtone state of the H2O monomer is presented. The intermolecular potential energies were evaluated using explicitly correlated coupled cluster theory [CCSD(T)-F12] with a large basis set. Four vibrationally averaged analytical intermolecular PESs for H2O-Kr with H2O molecules in its |00+〉, |02+〉, |02-〉, and |11+〉 states are obtained by fitting to the multi-dimensional Morse/Long-Range potential function form. Each vibrationally averaged PES fitted to 578 points has root-mean-square (RMS) deviations smaller than 0.14 cm-1 and requires only 58 parameters. The combined radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm were employed to calculate the rovibrational energy levels for |00+〉, |02+〉, |02-〉, and |11+〉 states of the H2O-Kr complexes. The calculated |02-〉Πf/e(101) ← |00+〉Σe(000) and |02+〉Πf/e(110) ← |00+〉Σe(101) infrared transitions are in excellent agreement with the experimental values with RMS discrepancies being only 0.007 and 0.016 cm-1, respectively. These analytical PESs can be used to provide reliable theoretical guidance for future infrared overtone spectroscopy of H2O-Kr.

Neural network interaction potentials for para-hydrogen with flexible molecules

The study of molecular impurities in para-hydrogen ( pH2) clusters is key to push forward our understanding of intra- and intermolecular interactions, including their impact on the superfluid response of this bosonic quantum solvent. This includes tagging with only one or very few pH2, the microsolvation regime for intermediate particle numbers, and matrix isolation with many solvent molecules. However, the fundamental coupling between the bosonic pH2 environment and the (ro-)vibrational motion of molecular impurities remains poorly understood. Quantum simulations can, in principle, provide the necessary atomistic insight, but they require very accurate descriptions of the involved interactions. Here, we present a data-driven approach for the generation of impurity⋯ pH2 interaction potentials based on machine learning techniques, which retain the full flexibility of the dopant species. We employ the well-established adiabatic hindered rotor (AHR) averaging technique to include the impact of the nuclear spin statistics on the symmetry-allowed rotational quantum numbers of pH2. Embedding this averaging procedure within the high-dimensional neural network potential (NNP) framework enables the generation of highly accurate AHR-averaged NNPs at coupled cluster accuracy, namely, explicitly correlated coupled cluster single, double, and scaled perturbative triples, CCSD(T*)-F12a/aVTZcp, in an automated manner. We apply this methodology to the water and protonated water molecules as representative cases for quasi-rigid and highly flexible molecules, respectively, and obtain AHR-averaged NNPs that reliably describe the corresponding H2O⋯ pH2 and H3O+⋯ pH2 interactions. Using path integral simulations, we show for the hydronium cation, H3O+, that umbrella-like tunneling inversion has a strong impact on the first and second pH2 microsolvation shells. The automated and data-driven nature of our protocol opens the door to the study of bosonic pH2 quantum solvation for a wide range of embedded impurities.

A path integral ground state approach for asymmetric top rotors with nuclear spin symmetry: Application to water chains

A path integral ground state (PIGS) approach for the simulation of asymmetric top rotors is presented. The method is based on Monte Carlo sampling of angular degrees of freedom. A symmetry-adapted rotational density matrix is used to account for nuclear spin statistics. To illustrate the method, ground-state properties of collections of para-water molecules confined to a one-dimensional lattice are computed. Those include energetic and structural observables. An advantage of the PIGS method is that expectation values can be obtained directly since the square of the wavefunction is sampled during a simulation. To benchmark the method, ground state energies and orientational distributions are computed using exact diagonalization for a single para-water molecule in an external field using a finite basis of symmetric top eigenfunctions. Benchmark results are also provided for N = 2 para-water molecules pinned to lattice sites at various distances to sample the crossover from hydrogen bonding to the dipole-dipole interaction regime. Excellent agreement between the PIGS results and the finite basis set calculations is observed. A thorough analysis of the convergence in terms of the imaginary time propagation length and systematic Trotter error is performed. The PIGS approach is then applied to a chain of N = 11 water molecules, and an equation of state is constructed in terms of the intermolecular separation. Ordering effects are also studied, and a transition between hydrogen bonding to dipole-dipole alignment is observed. The method is scalable and can also be applied in higher dimensions.

Equation of state and first principles prediction of the vibrational matrix shift of solid parahydrogen

We generate the equation of state (EOS) of solid parahydrogen (para-H₂) using a path-integral Monte Carlo (PIMC) simulation based on a highly accurate first-principles adiabatic hindered rotor potential energy curve for the para-H₂ dimer. The EOS curves for the fcc and hcp structures of solid para-H₂ near the equilibrium density show that the hcp structure is the more stable of the two, in agreement with experiment. To accurately reproduce the structural and energy properties of solid para-H₂, we eliminated by extrapolation the systematic errors associated with the choice of simulation parameters used in the PIMC calculation. We also investigate the temperature dependence of the EOS curves, and the invariance of the equilibrium density with temperature is satisfyingly reproduced. The pressure as a function of density and the compressibility as a function of pressure are both calculated using the obtained EOS and are compared with previous simulation results and experiments. We also report the first ever a priori prediction of a vibrational matrix shift from first-principles two-body potential functions, and its result for the equilibrium state agrees well with experiment.

Microscopic molecular superfluid response: theory and simulations

Since its discovery in 1938, superfluidity has been the subject of much investigation because it provides a unique example of a macroscopic manifestation of quantum mechanics. About 60 years later, scientists successfully observed this phenomenon in the microscopic world though the spectroscopic Andronikashvili experiment in helium nano-droplets. This reduction of scale suggests that not only helium but also para-H2 (pH2) can be a candidate for superfluidity. This expectation is based on the fact that the smaller number of neighbours and surface effects of a finite-size cluster may hinder solidification and promote a liquid-like phase. The first prediction of superfluidity in pH2 clusters was reported in 1991 based on quantum Monte Carlo simulations. The possible superfluidity of pH2 was later indirectly observed in a spectroscopic Andronikashvili experiment in 2000. Since then, a growing number of studies have appeared, and theoretical simulations have been playing a special role because they help guide and interpret experiments. In this review, we go over the theoretical studies of pH2 superfluid clusters since the experiment of 2000. We provide a historical perspective and introduce the basic theoretical formalism along with key experimental advances. We then present illustrative results of the theoretical studies and comment on the possible future developments in the field. We include sufficient theoretical details such that the review can serve as a guide for newcomers to the field.

Probing the Superfluid Response of para-Hydrogen with a Sulfur Dioxide Dopant

We recently presented the first attempt at using an asymmetric top molecule (para-water) to probe the superfluidity of nanoclusters (of para-hydrogen) [ Zeng , T. ; Li , H. ; Roy , P.-N. J. Phys. Chem. Lett. 2013 , 4 , 18 - 22 ]. Unfortunately, para-water could not be used to probe the para-hydrogen superfluid response. We now report a theoretical simulation of sulfur dioxide rotating in para-hydrogen clusters and show that this asymmetric top can serve as a genuine probe of superfluidity. With this probe, we predict that as few as four para-hydrogen molecules are enough to form a superfluid cluster, the smallest superfluid system to date. We also propose the concept of "exchange superfluid fraction" as a more precise measurement. New superfluid scenarios brought about by an asymmetric top dopant and potential experimental measurements are discussed.