Quantum simulation for muoniated and deuterated methyl radicals in implicit water solvent: combined ab initio path integral molecular dynamics and the polarizable continuum model simulation study

We have combined ab initio path integral molecular dynamics (PIMD) simulation and the polarizable continuum model (PCM) method to efficiently incorporate solvent effects into nuclear quantum fluctuation of molecular systems. Our combined ab initio PIMD–PCM simulation was applied to muoniated and deuterated methyl radical immersed in implicit water solvent to gain information on solvent and isotope effects from one simulation run. We found that solvent effects lead to the bond elongation and a decrease in the magnitude of isotropic hyperfine coupling constants. These are consistent with the trends in conventional static calculations and experiments. In addition, the performance of cavity models (universal force field, united atom specified for Kohn–Sham and these hybrid models) and the conservation of the PIMD–PCM Hamiltonian were accessed. We confirmed that solvent effects on nuclear quantum fluctuation are efficiently computed using our combined simulation of quantum solute in implicit solvent.