Analysis of the structural H/D isotope effect with an excess proton/deuteron in light/heavy water solvent using path integral molecular dynamics simulations

We analyzed the difference in the structural H/D isotope effect between an excess proton in light water (H-body) and an excess deuteron in heavy water (D-body), including the nuclear quantum effect, using path integral molecular dynamics simulations. We found that the second peak of the H-body is shorter than that of the D-body in the radial distribution function of O*-O, where O* is the oxygen atom of the H3O+/D3O+ fragment. The main reason for this would be the difference in the ratio of the Zundel structure with the sp3-like configuration, where the Zundel structure in the H-body (14.0%) is greater than that in the D-body (12.0%). We also found rare occurrences of double H3O+/D3O+ configurations, mainly including Zundel-Zundel-like structures such as H7O3+/D7O3+ and H9O4+/D9O4+. The ratios of such configurations appearing in our simulations are 0.89% and 0.20% for the H-body and the D-body, respectively. To elucidate the H/D isotope effect between an excess proton in light water (H-body) and an excess deuteron in heavy water (D-body), we performed path integral molecular dynamics simulations. We found the second peak of the H-body in the O*-O radial distribution function is slightly shorter than that of the D-body, where O* is the oxygen atom of the H3O+/D3O+ fragment.