X

Search Results

Searching....

11.2 Chemical Solvent Models

11.2.1 Introduction

(July 14, 2022)

Ab initio quantum chemistry makes possible the study of gas-phase molecular properties from first principles. In liquid solution, however, these properties may change significantly, especially in polar solvents. Although it is possible to model solvation effects by including explicit solvent molecules in the quantum-chemical calculation (e.g. a super-molecular cluster calculation, averaged over different configurations of the molecules in the first solvation shell), such calculations are very computationally demanding. Furthermore, cluster calculations typically do not afford accurate solvation energies, owing to the importance of long-range electrostatic interactions. (Hybrid discrete/continuum models, which contain some explicit solvent, can be quite effective, however. 960 Pliego Jr. J. R., Riveros J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2020), 10, pp. e1440.
Link
, 911 Paul S. K., Herbert J. M.
J. Am. Chem. Soc.
(2021), 143, pp. 10189.
Link
) Accurate prediction of solvation free energies is crucial for modeling of chemical reactions but also for relative conformational energies in solution.

Q-Chem contains several different implicit solvent models, which differ greatly in their level of sophistication. These are generally known as self-consistent reaction field (SCRF) models, because the continuum solvent establishes a “reaction field” (additional terms in the solute Hamiltonian) that depends upon the solute electron density, and must therefore be updated self-consistently during the iterative convergence of the wave function. The simplest and oldest of these models that is available in Q-Chem is the multipole expansion method, also known as a Kirkwood-Onsager model, 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
in which the solute molecule is placed inside of a spherical cavity and its electrostatic potential is represented in terms of a single-center multipole expansion. (This should not be confused with the Onsager model, in which a dipole approximation is used, and to avoid confusion the multipolar expansion method is better terminology.) More sophisticated models, which use a molecule-shaped cavity and the full molecular electrostatic potential, include the conductor-like PCM (C-PCM), 67 Barone V., Cossi M.
J. Phys. Chem. A
(1998), 102, pp. 1995.
Link
, 240 Cossi M. et al.
J. Comput. Chem.
(2003), 24, pp. 669.
Link
the conductor-like screening model (COSMO), the “surface and simulation of volume polarization for electrostatics” [SS(V)PE] model, 211 Chipman D. M.
J. Chem. Phys.
(2000), 112, pp. 5558.
Link
and the closely-related “integral equation formalism” (IEF-PCM). 160 Cancès E., Mennucci B., Tomasi J.
J. Chem. Phys.
(1997), 107, pp. 3032.
Link
, 161 Cancès E., Mennucci B.
J. Chem. Phys.
(2001), 114, pp. 4744.
Link
For an overview of all of these methods and their interconnections, see the review by Herbert. 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link

The C-PCM, IEF-PCM, and SS(V)PE are examples of what are called “apparent surface charge” SCRF models, and the term polarizable continuum models (PCMs), as popularized by Tomasi and coworkers, 1192 Tomasi J., Mennucci B., Cammi R.
Chem. Rev.
(2005), 106, pp. 2999.
Link
is now used almost universally to refer to this class of solvation models. Q-Chem employs a Switching/Gaussian or “SwiG” implementation of these PCMs,661, 662, 664, 484, 659 which resolves a long-standing (though little-publicized) problem with standard PCMs, namely, that the boundary-element methods used to discretize the solute/continuum interface may lead to discontinuities in the potential energy surface for the solute molecule. These discontinuities inhibit convergence of geometry optimizations, introduce serious artifacts in vibrational frequency calculations, and make ab initio molecular dynamics calculations virtually impossible. 661 Lange A. W., Herbert J. M.
J. Phys. Chem. Lett.
(2010), 1, pp. 556.
Link
, 662 Lange A. W., Herbert J. M.
J. Chem. Phys.
(2010), 133, pp. 244111.
Link
In contrast, Q-Chem’s SwiG-PCMs afford potential energy surfaces that are rigorously continuous and smooth. Unlike earlier attempts to obtain smooth PCMs, the SwiG approach largely preserves the properties of the underlying integral-equation solvent models, so that solvation energies and molecular surface areas are hardly affected by the smoothing procedure.

Other solvent models available in Q-Chem include the “Langevin dipoles” model; 339 Florián J., Warshel A.
J. Phys. Chem. B
(1997), 101, pp. 5583.
Link
, 340 Florián J., Warshel A.
J. Phys. Chem. B
(1999), 103, pp. 10282.
Link
as well as the SM8, 796 Marenich A. V. et al.
J. Chem. Theory Comput.
(2007), 3, pp. 2011.
Link
SM12, 794 Marenich A. V., Cramer C. J., Truhlar D. G.
J. Chem. Theory Comput.
(2013), 9, pp. 609.
Link
and SMD 793 Marenich A. V., Cramer C. J., Truhlar D. G.
J. Phys. Chem. B
(2009), 113, pp. 6378.
Link
models developed at the University of Minnesota. SM8 and SM12 are based upon the generalized Born method for electrostatics, augmented with atomic surface tensions intended to capture nonelectrostatic effects (cavitation, dispersion, exchange repulsion, and changes in solvent structure). These models have been carefully parameterized to reproduce experimental free energies of solvation. 242 Cramer C. J., Truhlar D. G.
Acc. Chem. Res.
(2008), 41, pp. 760.
Link
The SMD model, in which the “D” is for “density”, combines IEF-PCM with similar nonelectrostatic corrections. Statistically speaking, SMD is not any more or less accurate than other SMx models, 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
but has the advantage of being based on rigorous electrostatics derived from the solute’s exact SCF density. The SM12 and SMD models can each be used in arbitrary basis sets. The SM8 model uses generalized Born electrostatics based on “CM4” charges, 594 Kelly C. P., Cramer C. J., Truhlar D. G.
J. Chem. Theory Comput.
(2005), 1, pp. 1133.
Link
, 888 Olson R. M. et al.
J. Chem. Theory Comput.
(2007), 3, pp. 2046–2054.
Link
which are themselves based on Löwdin atomic charges, and as such this model should only be used in the basis sets for which it was parameterized: 6-31G*, 6-31+G*, or 6-31+G**. (Other basis sets, if requested will use the 6-31G* parameters but this is not recommended.) The SM12 model also uses generalized Born electrostatics but substitutes CM4 charges for Hirshfeld charges; the latter are more stable with respect to changes in basis set and therefore SM12 is available in arbitrary basis sets. The trade-off is that an analytic gradient is available for SM8 but not SM12; see Table 11.2. An analytic gradient is available for SMD also.

Model Cavity Non- Basis
Construction Discretization Electrostatic Sets
Terms? Supported
Kirkwood-Onsager spherical point charges no all
Langevin Dipoles atomic spheres dipoles in no all
(user-definable) 3-d space
Poisson Equation atomic spheres grid in no all
Solver (user-definable) 3-d space
C-PCM atomic spheres point charges or user- all
(user-definable) smooth Gaussians specified
SS(V)PE/ atomic spheres point charges or user- all
IEF-PCM (user-definable) smooth Gaussians specified
COSMO predefined point charges none all
atomic spheres
Isodensity SS(V)PE isodensity contour point charges none all
CMIRS isodensity contour point charges automatic all
SM8 predefined automatic 6-31G*
atomic spheres N/Aa 6-31+G*
6-31+G**
SM12 predefined N/Aa automatic all
atomic spheres
SMD predefined point charges automatic all
atomic spheres
aGeneralized Born electrostatic model; does not require cavity construction.
Table 11.1: Summary of implicit solvation models available in Q-Chem, indicating how the solute cavity is constructed and discretized, whether non-electrostatic terms are (or can be) included, and which basis sets are available for use with each model.

Table 11.1 summarizes the implicit solvent models that are available in Q-Chem. Solvent models are invoked via the SOLVENT_METHOD keyword, as shown below. Additional details about each particular solvent model can be found in the sections that follow. In general, these methods are available for any SCF level of electronic structure theory, with the aforementioned caveat about basis sets for SM8. Post-Hartree–Fock calculations (such as MP2 or EOM-CC) can be performed by first running an SCF + PCM job, in which case the correlated wave function will employ MOs and Hartree-Fock energy levels that are polarized by the solvent. This represents a “zeroth-order” inclusion of the solvent effects at the correlated level of theory, but is perfectly adequate for many applications as higher-order corrections are usually small. 817 Mewes J.-M., Herbert J. M., Dreuw A.
Phys. Chem. Chem. Phys.
(2017), 19, pp. 1644.
Link
, 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
Table 11.2 also summarizes the availability of analytical energy gradients for implicit solvent models. (Finite-difference gradients and Hessians are requested automatically for calculations where the requisite analytic derivatives are not available.)

Note:  The format for specifying implicit solvent models changed significantly starting in Q-Chem version 4.2.1. This change was made in an attempt to simply and unify the input notation for a large number of different models.

Energy Derivatives Solvent Model
C-PCM SS(V)PE & CMIRS COSMO SM8 SM12 SMD PEqS Langevin
  IEF-PCM Dipoles
SCF energy gradient yesa yesa no yes yes no yes no yes
SCF energy Hessian yesa no no yesc no no no no no
CIS/TDDFT energy gradient yesa no — unsupported —
CIS/TDDFT energy Hessian yesa no — unsupported —
MP2 & double-hybridb — unsupported —
Coupled cluster methodsb — unsupported —
aGradients available for van der Waals cavities and solvent-accessible surface (SAS) only
bGradients are not supported but single-point calculations can be performed using solvent-polarized MOs
cHessians of COSMO with the outlying charge correction (SOLVENT_METHOD = COSMO) are not supported.
Table 11.2: Summary of analytic energy gradient and Hessian capabilities with implicit solvent models.

SOLVENT_METHOD

SOLVENT_METHOD
       Sets the preferred solvent method.
TYPE:
       STRING
DEFAULT:
       0
OPTIONS:
       0 Do not use a solvation model. KIRKWOOD Use the Kirkwood-Onsager model (Section 11.2.2). PCM Use an apparent surface charge, polarizable continuum model (Section 11.2.3). ISOSVP Use the isodensity implementation of the SS(V)PE model (Section 11.2.6). COSMO Use COSMO (Section 11.2.8). SM8 Use version 8 of the Cramer-Truhlar SMx model (Section 11.2.9.1). SM12 Use version 12 of the SMx model (Section 11.2.9.2). SMD Use SMD (Section 11.2.9.3). CHEM_SOL Use the Langevin Dipoles model (Section 11.2.10). PEQS Use the Poisson Equation Solver (Section 11.2.11).
RECOMMENDATION:
       Consult the literature (e.g., Ref.  490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
). PCM is a collective name for a family of models and additional input options may be required in this case, in order to fully specify the model; see Section 11.2.3. Several versions of SM12 are available as well, as discussed in Section 11.2.9.2.

Before going into detail about each of these models, a few potential points of confusion warrant mention, with regards to nomenclature. First, “PCM” refers to a family of models that includes C-PCM, COSMO, SS(V)PE, and IEF-PCM. The latter two models are formally equivalent at the level of integral equations, 161 Cancès E., Mennucci B.
J. Chem. Phys.
(2001), 114, pp. 4744.
Link
, 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
but exhibit some some differences in their numerical implementation. 664 Lange A. W., Herbert J. M.
Chem. Phys. Lett.
(2011), 509, pp. 77.
Link
, 1330 You Z.-Q. et al.
J. Chem. Phys.
(2015), 143, pp. 204104.
Link
, 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
One or the other of these models can be selected by additional job control variables in a $pcm input section, as described in Section 11.2.3. COSMO is very similar to C-PCM but includes a correction for that part of the solute’s electron density that penetrates beyond the cavity (the so-called “outlying charge”), 607 Klamt A., Jonas V.
J. Chem. Phys.
(1996), 105, pp. 9972.
Link
although later work cast doubt on the theoretical justification for this and other ad hoc charge renormalization procedures. 210 Chipman D. M.
J. Chem. Phys.
(1999), 110, pp. 8012.
Link
[The IEF-PCM and SS(V)PE methods already contain an implicit correction for outlying charge, as does the C-PCM method that is derived as an approximation to these models, 210 Chipman D. M.
J. Chem. Phys.
(1999), 110, pp. 8012.
Link
, 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
although this was not recognized at the time that COSMO was formulated. See Ref.  490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
for a historical discussion of these developments.] In any case, COSMO is described in Section 11.2.8.

Two implementations of the SS(V)PE model are also available. The PCM implementation (which is requested by setting SOLVENT_METHOD = PCM in conjunction with appropriate job-control variables in the $pcm input section) uses a solute cavity constructed from atom-centered spheres, in keeping with other PCMs. On the other hand, setting SOLVENT_METHOD = ISOSVP requests an SS(V)PE calculation in which the solute cavity is defined by an isocontour of the solute’s own electron density. 208 Chipman D. M., Dupuis M.
Theor. Chem. Acc.
(2002), 107, pp. 90.
Link
, 204 Chen F., Chipman D. M.
J. Chem. Phys.
(2003), 119, pp. 10289.
Link
, 213 Chipman D. M.
J. Chem. Phys.
(2006), 124, pp. 224111.
Link
This is an appealing, one-parameter cavity construction that avoid many of the problems with cusps in “van der Waals” cavity surfaces that are constructed from atom-centered spheres, and the isodensity implementation of SS(V)PE forms the basis of a physics-based continuum solvation model called CMIRS that is described in Section 11.2.7. 967 Pomogaeva A., Chipman D. M.
J. Chem. Theory Comput.
(2011), 7, pp. 3952.
Link
, 968 Pomogaeva A., Chipman D. M.
J. Phys. Chem. A
(2013), 117, pp. 5812.
Link
, 969 Pomogaeva A., Chipman D. M.
J. Chem. Theory Comput.
(2014), 10, pp. 211.
Link
, 970 Pomogaeva A., Chipman D. M.
J. Phys. Chem. A
(2015), 119, pp. 5173.
Link
, 1327 You Z.-Q., Herbert J. M.
J. Chem. Theory Comput.
(2016), 12, pp. 4338.
Link
CMIRS is competitive in accuracy (for solvation free energies) with the best-available SMx solvation models, 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
despite using far fewer fitting parameters. However, analytic energy gradients are not available for the isodensity cavity construction and therefore not available for CMIRS, whereas these gradients are available when the cavity surface is constructed from atom-centered spheres.

Regarding the accuracy of these models for solvation free energies (ΔG298), the SMx models generally achieve sub-kcal/mol accuracy for neutral molecules, based on comparison to a large database of experimental values, although average errors for ions are more like 4 kcal/mol. 242 Cramer C. J., Truhlar D. G.
Acc. Chem. Res.
(2008), 41, pp. 760.
Link
, 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
(Note that the SM12 and SMD models generally do not improve too much upon SM8 in any statistical sense, 793 Marenich A. V., Cramer C. J., Truhlar D. G.
J. Phys. Chem. B
(2009), 113, pp. 6378.
Link
, 794 Marenich A. V., Cramer C. J., Truhlar D. G.
J. Chem. Theory Comput.
(2013), 9, pp. 609.
Link
but do extend these models to arbitrary basis sets whereas SM8 is limited to a few small basis sets.) To achieve accuracy comparable to the SMx models within the PCM class of solvent models, nonelectrostatic terms must be added. 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
Among the various PCMs described above, the only one that constitutes a “black box” model for solvation energies is CMIRS, which slightly outperforms the SMx models in a statistical sense, 1327 You Z.-Q., Herbert J. M.
J. Chem. Theory Comput.
(2016), 12, pp. 4338.
Link
, 490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
although it is only available for a few solvents.

The following sections provide more details regarding theory and job control for the various implicit solvent models that are available in Q-Chem. Ref.  490 Herbert J. M.
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
(2021), 11, pp. e1519.
Link
contains both formal comparisons amongst these models as well as a side-by-side comparison of the accuracy of solvation free energies.