Vibrational Self-Consistent Field (VSCF) and Post-VSCF Method Calculations Combined with the Reference Interaction Site Model Self-Consistent Field Method Coupled with the Constrained Spatial Electron Density Distribution: Applications to NaHCOO in Aqueous Phase

Investigating vibrational behavior in solution is crucial for understanding molecular dynamics within a solvent environment. Notably, the analysis of Raman spectra for molecules in solution is important owing to its ability to unveil intricate solute-solvent interactions. Previous studies have effectively employed frequency calculations utilizing the reference interaction site model self-consistent field method in conjunction with constrained spatial electron density distribution (RISM-SCF-cSED) to understand molecular vibrations in solution, primarily focusing on fundamental vibrational modes. However, the oversight of overtones and combination tones in these studies prompted us to combine the vibrational self-consistent field (VSCF) and vibrational second-order Mo̷ller-Plesset perturbation (VMP2) methods with RISM-SCF-cSED to address these aspects theoretically. Illustrating the efficacy of this integrated approach, we computed the Raman spectra of sodium formate (NaHCOO) in water, revealing the necessity of accounting for molecular anharmonicity in solution vibrational analysis. Our findings underscore the potency of VSCF and VMP2 in conjunction with RISM-SCF-cSED as a robust theoretical framework for such calculations.

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.

Theoretical Study of Raman Intensities of p-Nitroaniline in Different Solvent Conditions by Using a Reference Interaction Site Model Self-Consistent Field Explicitly Including Constrained Spatial Electron Density Distribution

Raman spectroscopy is one of the most powerful tools to understand and characterize the states and structures of systems in several environments. To obtain highly accurate changes in Raman intensities of systems in solution, theoretical treatment, which can deal with not only the states and structures of systems but also the environment around molecules, proves to be significant. Hence, in this study, we developed the calculation of changes in Raman intensities of systems in different solvent conditions by using the reference interaction site model self-consistent field study explicitly including constrained spatial electron density distribution; this model is designed based on elements from both quantum mechanics and statistical mechanics. We showed that our calculation method could reproduce the changes in Raman intensities of p-nitroaniline (pNA) under different solvent conditions, including supercritical water, which has been observed in previous experimental studies. Based on the analysis of the calculation results, we observed that the ratio of the Raman intensity change of pNA in different solvent conditions is strongly correlated with the charge-transfer character of pNA.

Analytical second derivatives of the free energy in solution by the reference interaction site model self-consistent field explicitly including constrained spatial electron density distribution

The application of analytical derivative methods to solution systems is important because several chemical reactions occur in solution. The reference interaction site model (RISM) is one of the solvation theories used to study solution systems and has shown good performance, especially in the polar solvent systems. Although the analytical first derivative based on the RISM coupled with quantum methods (RISM-SCF) has already been derived, the analytical second derivative has not been proposed yet. Therefore, in this study, the analytical second derivative was derived using RISM-SCF explicitly including constrained spatial electron density distribution (RISM-SCF-cSED). The performance of this method was validated with the Hessian calculations of formaldehyde and para-nitroaniline in solution, and the results demonstrated that the method accurately calculated frequency values at a small computational cost.

Accurate Absorption Energy Calculations in Solution Using the Reference Interaction Site Model Self-Consistent Field Including the Constrained Spatial Electron Density Distribution

The solvation effect is an important factor determining the properties of molecules in solution. The reference interaction site model (RISM) is a powerful method to treat the solvation effect with pair-correlation functions, such as a radial distribution function. This study developed a hybrid method between quantum mechanics and RISM using the spatial electron density distributions on each atomic site (RISM-SCF-cSED). Sophisticated quantum mechanical approaches can be used to consider the solvation effect because the computational cost of RISM-SCF-cSED is reasonable. In this study, the absorption energies of 5-(dimethylamino)-2,4-pentadienal in various solutions were calculated using RISM-SCF-cSED. The experimental data were well reproduced with an average errors of ∼0.06 eV, using multi-reference perturbation theory.

Improved RISM-CASSCF Optimization via State-Average Treatment and Damping for Characterizing Excited Molecules in Solution with Multireference Perturbation Theory

Determination of excited states of near-infrared (NIR) bioimaging dyes is a challenging theoretical task because of their energy levels with a small gap and the presence of solvation. In the previous study, we showed that the development of the reference interaction site model coupled with the complete active space second-order perturbation theory, the RISM-CASPT2 method, and its extension with the density matrix renormalization group enabled high accuracy prediction of the photochemical properties of bioimaging-related fluorescent molecules in solution (Shimizu et al., J. Chem. Theory Comput. 2018, 14, 5673-5679). This method, however, has a technical issue in convergence of CASSCF optimization, which was encountered when applying the method to a wider variety of systems; thus, practical applications have been hindered. Here, we present an improved scheme of CASSCF optimization with and without the density matrix renormalization group treatment. Detailed derivations and analysis of the second-order orbital optimization scheme with the inclusion of solvation through RISM revealed the requirement of a correction term to the orbital Hessian matrix. As a practical approach, the state-average RISM-CASPT2 method with damping treatment for solvation is presented for improving the convergence of the calculation under reasonable computational cost. The improved scheme allows for performing accurate and numerically stable theoretical analysis of the bioimaging-related excited state with various types of solvation for a P═O-bridged rhodol derivative, which is recently highlighted as a promising photostable NIR dye molecule.

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.

3D-RISM-MP2 Approach to Hydration Structure of Pt(II) and Pd(II) Complexes: Unusual H-Ahead Mode vs Usual O-Ahead One

Solvation of transition metal complexes with water has been one of the fundamental topics in physical and coordination chemistry. In particular, Pt(II) complexes have recently attracted considerable interest for their relation to anticancer activity in cisplatin and its analogues, yet the interaction of the water molecule and the metal center has been obscured. The challenge from a theoretical perspective remains that both the microscopic solvation effect and the dynamical electron correlation (DEC) effect have to be treated simultaneously in a reasonable manner. In this work we derive the analytical gradient for the three-dimensional reference interaction site model Møller-Plesset second order (3D-RISM-MP2) free energy. On the basis of the three-regions 3D-RISM self-consistent field (SCF) method recently proposed by us, we apply a new layer of the Z-vector method to the CP-RISM equation as well as point-charge approximation to the derivatives with respect to the density matrix elements in the RISM-CPHF equation to remarkably reduce the computational cost. This method is applied to study the interaction of H2O with the d(8) square planar transition metal complexes in aqueous solution, trans-[Pt(II)Cl2(NH3)(glycine)] (1a), [Pt(II)(NH3)4](2+) (1b), [Pt(II)(CN)4](2-) (1c), and their Pd(II) analogues 2a, 2b, and 2c, respectively, to elucidate whether the usual H2O interaction through O atom (O-ahead mode) or unusual one through H atom (H-ahead mode) is stable in these complexes. We find that the interaction energy of the coordinating water and the transition metal complex changes little when switching from gas to aqueous phase, but the solvation free energy differs remarkably between the two interaction modes, thereby affecting the relative stability of the H-ahead and O-ahead modes. Particularly, in contrast to the expectation that the O-ahead mode is preferred due to the presence of positive charges in 1b, the H-ahead mode is also found to be more stable. The O-ahead mode is found to be more stable than the H-ahead one only in 2b. The energy decomposition analysis (EDA) at the 3D-RISM-MP2 level revealed that the O-ahead mode is stabilized by the electrostatic (ES) interaction, whereas the H-ahead one is mainly stabilized by the DEC effect. The ES interaction is also responsible for the difference between the Pd(II) and Pt(II) complexes; because the electrostatic potential is more negative along the z-axis in the Pt(II) complex than in the Pd(II) one, the O-ahead mode prefers the Pd(II) complexes, whereas the H-ahead becomes predominant in the Pt(II) complexes.