Solvent Effects on the Symmetric and Asymmetric SN2 Reactions in the Acetonitrile Solution: A Reaction Density Functional Theory Study

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.

Solvent Effects on the Menshutkin Reaction

The Menshutkin reaction is a methyl transfer reaction relevant in fields ranging from biochemistry to chemical synthesis. In the present work, the energetics and solvent distributions for NH₃+MeCl and Pyr+MeBr reactions were investigated in explicit solvent (water, methanol, acetonitrile, benzene, cyclohexane) by means of reactive molecular dynamics simulations. For polar solvents (water, methanol, and acetonitrile) and benzene, strong to moderate catalytic effects for both reactions were found, whereas apolar and bulky cyclohexane interacts weakly with the solute and does not show pronounced barrier reduction. The calculated barrier heights for the Pyr+MeBr reaction in acetonitrile and cyclohexane are 23.2 and 28.1 kcal/mol compared with experimentally measured barriers of 22.5 and 27.6 kcal/mol, respectively. The solvent distributions change considerably between reactant and TS but comparatively little between TS and product conformations of the solute. As the system approaches the transition state, correlated solvent motions occur which destabilize the solvent-solvent interactions. This is required for the system to surmount the barrier. Finally, it is found that the average solvent-solvent interaction energies in the reactant, TS, and product state geometries are correlated with changes in the solvent structure around the solute.

Tackling Solvent Effects by Coupling Electronic and Molecular Density Functional Theory

Solvation effects can have a tremendous influence on chemical reactions. However, precise quantum chemistry calculations are most often done either in vacuum neglecting the role of the solvent or using continuum solvent model ignoring its molecular nature. We propose a new method coupling a quantum description of the solute using electronic density functional theory with a classical grand-canonical treatment of the solvent using molecular density functional theory. Unlike a previous work, both densities are minimized self-consistently, accounting for mutual polarization of the molecular solvent and the solute. The electrostatic interaction is accounted using the full electron density of the solute rather than fitted point charges. The introduced methodology represents a good compromise between the two main strategies to tackle solvation effects in quantum calculation. It is computationally more effective than a direct quantum mechanics/molecular mechanics coupling, requiring the exploration of many solvent configurations. Compared to continuum methods, it retains the full molecular-level description of the solvent. We validate this new framework onto two usual benchmark systems: a water solvated in water and the symmetrical nucleophilic substitution between chloromethane and chloride in water. The prediction for the free energy profiles are not yet fully quantitative compared to experimental data, but the most important features are qualitatively recovered. The method provides a detailed molecular picture of the evolution of the solvent structure along the reaction pathway.