Quantum locality and equilibrium properties in low-temperature parahydrogen: a multiscale simulation study.

Parahydrogen is the spin-zero singlet state of molecular hydrogen, which at low temperature (between 14 and 25 K) is in a fluid state. A classical treatment of the system leads to unphysical freezing, and the inclusion of quantum delocalization of the molecule is then required to obtain a realistic description of its equilibrium properties. In the present work, we employ the classical-quantum adaptive resolution method AdResS to investigate the spatial extension of quantum delocalization effects in the bulk fluid at low temperature. Specifically, we simulate a small, spherical region of the system in full quantum detail: this region is coupled to a bulk of coarse-grained particles with classical, quantum-derived effective interactions obtained from quantum simulations. The two regions are interfaced through open boundaries and in conditions of thermodynamic equilibrium. Structural properties of the fluid, namely, pair distribution functions, are measured for different sizes of the quantum region. The results of this work show that, for the thermodynamic conditions corresponding to the range of temperature between 14 and 25 K, the bead-based, quantum structural properties of low-temperature parahydrogen are deemed local and do not require the support of an explicit quantum bulk.

[1]  Richard Phillips Feynman,et al.  Atomic Theory of the Two-Fluid Model of Liquid Helium , 1954 .

[2]  L Delle Site,et al.  Adaptive resolution simulation of liquid para-hydrogen: testing the robustness of the quantum-classical adaptive coupling. , 2011, Physical chemistry chemical physics : PCCP.

[3]  Matej Praprotnik,et al.  Multiscale simulation of soft matter: from scale bridging to adaptive resolution. , 2008, Annual review of physical chemistry.

[4]  Victor V. Goldman,et al.  The isotropic intermolecular potential for H2 and D2 in the solid and gas phases , 1978 .

[5]  M. Boninsegni.  Quantum statistics and the momentum distribution of liquid parahydrogen , 2009, 0902.2727.

[6]  Matej Praprotnik,et al.  FAST TRACK COMMUNICATION: Fractional dimensions of phase space variables: a tool for varying the degrees of freedom of a system in a multiscale treatment , 2007 .

[7]  L Delle Site,et al.  Classical to path-integral adaptive resolution in molecular simulation: towards a smooth quantum-classical coupling. , 2010, Physical review letters.

[8]  Matej Praprotnik,et al.  Adaptive resolution scheme for efficient hybrid atomistic-mesoscale molecular dynamics simulations of dense liquids. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  K. Gernoth,et al.  Screening of particle exchange in quantum Boltzmann liquids , 2007 .

[10]  Matej Praprotnik,et al.  Adaptive molecular resolution via a continuous change of the phase space dimensionality. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[11]  Dirk Reith,et al.  Deriving effective mesoscale potentials from atomistic simulations , 2002, J. Comput. Chem..

[12]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[13]  Kurt Kremer,et al.  Comparative atomistic and coarse-grained study of water: What do we lose by coarse-graining? , 2009, The European physical journal. E, Soft matter.

[14]  I. Silvera.  The solid molecular hydrogens in the condensed phase: Fundamentals and static properties , 1980 .

[15]  R. Feynman,et al.  Quantum Mechanics and Path Integrals , 1965 .

[16]  Cecilia Clementi,et al.  Communication: On the locality of hydrogen bond networks at hydrophobic interfaces. , 2010, The Journal of chemical physics.

[17]  K. Kremer,et al.  Adaptive resolution molecular-dynamics simulation: changing the degrees of freedom on the fly. , 2005, The Journal of chemical physics.

[18]  L Delle Site,et al.  Comment on "Adaptive multiscale molecular dynamics of macromolecular fluids". , 2011, Physical review letters.

[19]  Alexander Lukyanov,et al.  Versatile Object-Oriented Toolkit for Coarse-Graining Applications. , 2009, Journal of chemical theory and computation.

[20]  S. G. Penoncello,et al.  Fundamental Equations of State for Parahydrogen, Normal Hydrogen, and Orthohydrogen , 2009 .

[21]  L Delle Site,et al.  Adaptive resolution molecular dynamics simulation through coupling to an internal particle reservoir. , 2011, Physical review letters.

[22]  Rüdiger Esser.  Workshop on Molecular Dynamics on Parallel Computers, John von Neumann Institute for Computing (NIC) Research Centre, Jülich, Germany, 8-10 February 1999 , 2000 .

[23]  L.Delle Site,et al.  Some fundamental problems for an energy-conserving adaptive-resolution molecular dynamics scheme. , 2007, 0709.2579.