Nuclear quantum effects in phase transition between Ice VII and Ice X.

The theoretical modeling of high-pressure ice remains challenging owing to the complexity in accurately reflecting its properties attributable to nuclear quantum effects. To explore the nuclear quantum effects of the phase transition between Ice VII and Ice X, we introduce an approach based on ab initio path-integral molecular dynamics. The results indicate that quantum effects facilitate the phase transition, with the observed isotope effects consistent with the experimental outcomes. We demonstrate that quantum effects manifest differently across ice phases: In Ice VII, quantum effects reduce the pressure through the centralization of protons. In contrast, in Ice X, quantum effects increase the pressure owing to the increased kinetic energy of zero-point vibration.

Location of the Shared Proton in Proton-Bound Dimer Compound of Hydrogen Sulfate and Formate: Path Integral Molecular Dynamics Study.

The structure of the proton-bound dimer compound of hydrogen sulfate and formate has been studied by considering nuclear quantum effects (NQEs) using the path integral molecular dynamics method. This study unveiled the location of the shared proton and answered the following question: "Is the shared proton localized on either an anion or located around the center of two anions?" We have elucidated that the shared proton is distributed in the region beyond the transition state due to the NQEs, even though the shared proton did not completely overcome the transition state for the proton shuttle.

A path integral molecular dynamics study on the muoniated xanthene-thione molecule.

A positive Mu is a useful tool for investigating the spin density of radical species. The theoretical estimation of its behavior in a molecule requires the inclusion of a quantum effect due to the small mass of muonium. Herein, we performed ab initio a path integral molecular dynamics (PIMD) simulation, which accurately included a multi-dimensional quantum effect, for muoniated 9H-xanthene-9-thione (µXT). Our results showed that the quantum effect significantly increased the hyperfine coupling constant (HFCC) value of µXT, which qualitatively improved the calculated HFCC value, compared to the experimental one. In the PIMD simulation, the bond length between muonium and sulfur in µXT is longer than that between hydrogen and sulfur in a hydrogenated 9H-xanthene-9-thione (HXT), leading to a spin density transfer from XT (9H-xanthene-9-thione) to muonium due to neutral dissociations. Additionally, we found that the S-Mu bond in µXT prefers a structure perpendicular to the molecular plane, where the interaction between Mu and the singly occupied molecular orbital of µXT is the strongest. These structural changes resulted in a larger HFCC value in the PIMD simulation of µXT.

Nuclear quantum and H/D isotope effects on intramolecular hydrogen bond in curcumin.

Curcumin and its derivatives possess intramolecular low-barrier hydrogen bonds for intramolecular proton transfer. The π-delocalization in the OCCCO framework of the hydrogen bond in these compounds is reorganized concomitantly with the proton transfer. To characterize the hydrogen bond and π-delocalization, we performed path integral molecular dynamics simulations, revealing that although the proton migration and reorganization of the π-delocalized structure showed a positive correlation, the correlation was weak.

Direct Elucidation of the Vibrationally Averaged Structure of Benzene: A Path Integral Molecular Dynamics Study.

Path integral molecular dynamics (PIMD) simulations for C6H6, C6D6, and C6T6 have been carried out to directly estimate the distribution of projected C-H(D,T) bond lengths onto the principal axis plane. The average values of raw C-H(D,T) bond lengths obtained from PIMD simulations are in the order of ⟨RC-H⟩ > ⟨RC-D⟩ > ⟨RC-T⟩ due to the anharmonicity of the potential energy curve. However, the projected C-H(D,T) bond lengths are almost the same as those reported by Hirano et al. [J. Mol. Struct. 2021, 1243, 130537]. Our PIMD simulations directly and strongly support the explanation by Hirano et al. for the experimental observations that almost the same projected C-H(D) bond lengths are found for C6H6 and C6D6. The PIMD simulations also predicted the same projected bond lengths for C6T6 as those of C6H(D)6. In addition to the previous local mode analysis, the present PIMD simulations predicted, for benzene isotopologues, that the vibrationally averaged structure is planar but non-flat.

A path integral molecular dynamics study on the NH4+ rotation in NH4+⋯XH2 (X = Be or Mg) dihydrogen bond systems.

The nuclear quantum effects (NQEs) in dihydrogen bond (DHB) complexes, i.e., NH4+⋯BeH2 and NH4+⋯MgH2, have been investigated using multicomponent quantum mechanics (MC_QM) calculations and path integral molecular dynamics (PIMD) simulation. The MC_QM method considers the NQEs, whereas PIMD considers both the NQEs and the thermal effects. The linear C3v structure is maintained in the optimized structures obtained by the static MP2 and MC_MP2 calculations, whereas the average structures obtained by the PIMD simulation are nonlinear. The strong DHB interaction in NH4+⋯MgH2 suppresses the fluctuation in the Hδ+NMg and Hδ-MgN angles, and hence, the NH4+ rotation did not occur in the simulation of NH4+⋯MgH2. The analysis of the radius of gyration revealed that the nuclear quantum fluctuation in the perpendicular direction is suppressed by the formation of the DHB complex, whereas that in the parallel direction is slightly enhanced in both the Hδ+ and Hδ- nuclei.

Low-Barrier Hydrogen Bond in Fujikurin A-D: A Computational Study.

The compounds Fujikurin A, B, and D, recently isolated from Fusarium fujikuroi, possess intramolecular low-barrier hydrogen bonds (LBHBs), which are hydrogen bonds with a very low-energy barrier for proton transfer. The isolated compounds have a hydrogen-bonded proton that appears to rapidly switch between two equilibrium states via a transition state (TS). To understand the characteristics of these intramolecular LBHBs in detail, we performed path integral molecular dynamics (PIMD) simulations, which can consider nuclear quantum effects (NQEs) under a finite temperature. The PIMD simulations predicted that the NQE completely washed out the energy barrier for the proton transfer reaction. Consequently, a single-well shape emerged in the results, along with the effective free-energy potential surface for the hydrogen-bonded proton distribution. Thus, we conclude that the hydrogen-bonded proton in Fujikurin does not in fact transfer between two equilibrium structures but widely delocalizes around the global minimum structure involving the TS region.

Carbon Monoxide Hydrogenation on Ice Surfaces.

We have performed density functional calculations to investigate the carbon monoxide hydrogenation reaction (H+CO→HCO), which is important in interstellar clouds. We found that the activation energy of the reaction on amorphous ice is lower than that on crystalline ice. In the course of this study, we demonstrated that it is roughly possible to use the excitation energy of the reactant molecule (CO) in place of the activation energy. This relationship holds also for small water clusters at the CCSD level of calculation and the two-layer-level ONIOM (CCSD : X3LYP) calculation. Generally, since it is computationally demanding to estimate activation energies of chemical reactions in a circumstance of many water molecules, this relationship enables one to determine the activation energy of this reaction on ice surfaces from the knowledge of the excitation energy of CO only. Incorporating quantum-tunneling effects, we discuss the reaction rate on ice surfaces. Our estimate that the reaction rate on amorphous ice is almost twice as large as that on crystalline ice is qualitatively consistent with the experimental evidence reported by Hidaka et al. [Chem. Phys. Lett., 2008, 456, 36.].