Polarizable continuum model

The polarizable continuum model (PCM) is a commonly used method in computational chemistry to model solvation effects. If it is necessary to consider each solvent molecule as a separate molecule, the computational cost of modeling a solvent-mediated chemical reaction would grow prohibitively high. Modeling the solvent as a polarizable continuum, rather than individual molecules, makes ab initio computation feasible. Two types of PCMs have been popularly used: the dielectric PCM (D-PCM) in which the continuum is polarizable (see dielectrics) and the conductor-like PCM (C-PCM) in which the continuum is conductor-like similar to COSMO Solvation Model.[1][2]

The molecular free energy of solvation is computed as the sum of three terms:

Gsol = Ges + Gdr + Gcav
Ges = electrostatic
Gdr = dispersion-repulsion
Gcav = cavitation[3]

The Charge-transfer effect is also considered as a part of solvation in cases.[1]

The PCM solvation model is available for calculating energies and gradients at the Hartree–Fock and density functional theory (DFT) levels in several quantum chemical computational packages such as Gaussian, GAMESS[3] and JDFTx.

The authors of a 2002 paper observe that PCM has limitations where non-electrostatic effects dominate the solute-solvent interactions. They write in the abstract: "Since only electrostatic solute-solvent interactions are included in the PCM, our results lead to the conclusion that, for the seven molecules studied, in cyclohexane, acetone, methanol, and acetonitrile electrostatic effects are dominant while in carbon tetrachloride, benzene, and chloroform other nonelectrostatic effects are more important."[4]

There is an integral equation formalism (IEF) version of the PCM which is very commonly used.[5]

PCM is also used to model outer solvation layers in multi-layered solvation approach.[6]

  1. ^ a b Jacopo Tomasi, Benedetta Mennucci, and Roberto Cammi (2005). "Quantum Mechanical Continuum Solvation Models." Chem. Rev. 105(8): 2999-3094.[1]
  2. ^ Maurizio Cossi, Nadia Rega, Giovanni Scalmani, Vincenzo Barone (2003). "Energies, structures, and electronic properties of molecules in solution with the C-PCM solvation model." J. Comput. Chem. 24(6): 669-681.[2]
  3. ^ a b Hendrik Zipse (9 February 2004). "The Polarizable Continuum Model (PCM)". Archived from the original on September 28, 2011. Retrieved January 25, 2009.
  4. ^ B. Mennucci et al. "Polarizable Continuum Model (PCM) Calculations of Solvent Effects on Optical Rotations of Chiral Molecules." J. Phys. Chem. A 2002, 106, 6102-6113. Link to full text
  5. ^ Mennucci, B.; Cancès, E.; Tomasi, J. (December 1997). "Evaluation of Solvent Effects in Isotropic and Anisotropic Dielectrics and in Ionic Solutions with a Unified Integral Equation Method: Theoretical Bases, Computational Implementation, and Numerical Applications". The Journal of Physical Chemistry B. 101 (49): 10506–10517. doi:10.1021/jp971959k.
  6. ^ Mark S. Gordon "CLUSTER-BASED APPROACHES TO SOLVATION" Iowa State University, Ames Laboratory.[3] Archived 2012-02-28 at the Wayback Machine