Ab initio path integral Monte Carlo simulations for water trimer with electron correlation effects

Structures of liquid and aqueous water isotopologues at ambient temperature from ab initio path integral simulations

The heavy hydrogen isotopes D and T are found in trace amounts in water, and when their concentration increases they can play an intricate role in modulating the physical properties of the liquid. We present an analysis of the microscopic structures of ambient light water (H₂O(l)), heavy water (D₂O(l)), T₂O(l), HDO(aq) and HTO(aq) studied by ab initio path integral molecular dynamics (PIMD). Unlike previous ab initio PIMD investigations of H₂O(l) and D₂O(l) [Chen et al., Phys. Rev. Lett., 2003, 91, 215503] [Machida et al., J. Chem. Phys., 2017, 148, 102324] we find that D₂O(l) is more structured than H₂O(l), as is predicted by the experiment. The agreement between the experiment and our simulation for H₂O(l) and D₂O(l) allows us to accurately predict the intra- and intermolecular structures of T₂O(l) HDO(aq) and HTO(aq). T₂O(l) is found to have a similar intermolecular structure to that of D₂O(l), while the intramolecular structure is more compact, giving rise to a smaller dipole moment than those of H₂O(l) and D₂O(l). For the mixed isotope species, HDO(aq) and HTO(aq), we find smaller dipole moments and fewer hydrogen bonds when compared with the pure species H₂O and D₂O. We can attribute this effect to the relative compactness of the mixed isotope species, which results in a lower dipole moment than that of the pure species.

Ab initio study of nuclear quantum effects on sub- and supercritical water

The structures of water in the ambient, subcritical, and supercritical conditions at various densities were studied systematically by ab initio path integral molecular dynamics simulations. It was found that the nuclear quantum effects (NQEs) have a significant impact on the structure of hydrogen bonds in close contact, not only in the ambient condition but also in the sub- and supercritical conditions. The NQEs on the structure beyond the hydrogen bond contact are important in ambient water, but not much for water in the sub- and supercritical conditions. The NQEs are furthermore important for determining the number of hydrogen bonds in the ambient conditions, and this role is, however, diminished in the sub- and supercritical conditions. The NQEs do, nevertheless, show their importance in determining the intramolecular structure of water and the close contact structures of the hydrogen bonds, even at sub- and supercritical conditions. Using the RPBE-D3 functional, the computed radial distribution functions for ambient water are in excellent agreement with experimental data, upgrading our previous results using the BLYP-D2 functional [Machida et al., J. Chem. Phys. 148, 102324 (2018)]. The computed radial distribution functions for water in the sub- and supercritical conditions were carefully compared with experiment. In particular, we found that the first peak in hydrogen pair distribution functions matches only when the NQEs are taken into account.

Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems

Recent developments in the fragment molecular orbital (FMO) method for theoretical formulation, implementation, and application to nano and biomolecular systems are reviewed. The FMO method has enabled ab initio quantum-mechanical calculations for large molecular systems such as protein-ligand complexes at a reasonable computational cost in a parallelized way. There have been a wealth of application outcomes from the FMO method in the fields of biochemistry, medicinal chemistry and nanotechnology, in which the electron correlation effects play vital roles. With the aid of the advances in high-performance computing, the FMO method promises larger, faster, and more accurate simulations of biomolecular and related systems, including the descriptions of dynamical behaviors in solvent environments. The current status and future prospects of the FMO scheme are addressed in these contexts.