Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids

The significance of solvent effects in electronic structure calculations has long been noted, and various methods have been developed to consider this effect. The reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is a hybrid model that combines the integral equation theory of molecular liquids with quantum chemistry. This method can consider the statistically convergent solvent distribution at a significantly lower cost than molecular dynamics simulations. Because the RISM theory explicitly considers the solvent structure, it performs well for systems where hydrogen bonds are formed between the solute and solvent molecules, which is a challenge for continuum solvent models. Taking advantage of being founded on the variational principle, theoretical developments have been made in calculating various properties and incorporating electron correlation effects. In this review, we organize the theoretical aspects of RISM-SCF-cSED and its distinctions from other hybrid methods involving integral equation theories. Furthermore, we carefully present its progress in terms of theoretical developments and recent applications.

Coupled Cluster Theory Combined with Reference Interaction Site Model Self-Consistent Field Explicitly Including Spatial Electron Density Distribution

Calculating the geometry and energy of a molecule in a solution is one of the most important tasks in chemistry. However, performing an accurate calculation in a solution is still a difficult task because the electronic structure and solute-solvent interactions are required to be accurately evaluated with an efficient computational cost. To overcome this difficulty, we proposed the coupled cluster with single and double excitations and perturbative triple excitations combined with the reference interaction site model (RISM) by employing our fitting approach. Our method correctly reproduced the relative stabilities of 1,2,3-triazole, isonicotinic acid, cytosine, and 6-chloro-2-pyridone in the aqueous phase, whereas the dielectric continuum model provided incorrect results for isonicotinic acid and 6-chloro-2-pyridone. Our method provided accurate results because the RISM captured the local solvation structure, such as hydrogen bonds.