COSMO solvation model

Last updated
COSMO surface of a pentaacrylate molecule (red = negative, green = positive equilibrium layer). Pentaacrylat.png
COSMO surface of a pentaacrylate molecule (red = negative, green = positive equilibrium layer).
Charge density surface of 4-nitro-benzoicacid. Calculated with COSMO. 4-Nitrobenzoic-acid-elpot-3D-vdW.png
Charge density surface of 4-nitro-benzoicacid. Calculated with COSMO.

COSMO [1] [2] (COnductor-like Screening MOdel) is a calculation method for determining the electrostatic interaction of a molecule with a solvent. COSMO is a dielectric continuum model [1] [3] [4] [5] (a.k.a. continuum solvation model). These models can be used in computational chemistry to model solvation effects. COSMO has become a popular method of these solvation models in recent years. The COSMO formalism is similar to the method proposed earlier by Hoshi et al. [1] [6] The COSMO approach is based - as many other dielectric continuum models - on the surface segmentation of a molecule surface [3] (usually referred to as 'solvent accessible surface' SAS approach).

Contents

Continuum solvation models - such as COSMO - treat each solvent as a continuum with a permittivity . Continuum solvation models approximate the solvent by a dielectric continuum, surrounding the solute molecules outside of a molecular cavity. In most cases it is constructed as an assembly of atom-centered spheres with radii approximately 20% larger than the Van der Waals radius. For the actual calculation the cavity surface is approximated by segments, e.g., hexagons, pentagons, or triangles.

Unlike other continuum solvation models, COSMO derives the polarization charges of the continuum, caused by the polarity of the solute, from a scaled-conductor approximation. If the solvent were an ideal conductor the electric potential on the cavity surface must disappear. If the distribution of the electric charge in the molecule is known, e.g. from quantum chemistry, then it is possible to calculate the charge on the surface segments. For solvents with finite dielectric constant this charge is lower by approximately a factor :

The factor is approximately

where the value of should be set to 0.5 for neutral molecules and to 0.0 for ions, see original derivation. [2] The value of is erroneously set to 0 in the popular C-PCM reimplementation of COSMO in Gaussian.

From the thus determined solvent charges and the known charge distribution of the molecule, the energy of the interaction between the solvent and the solute molecule can be calculated.

The COSMO method can be used for all methods in theoretical chemistry where the charge distribution of a molecule can be determined, for example semiempirical calculations, Hartree–Fock-method calculations or density functional theory (quantum physics) calculations. [1]

Variants and implementations

COSMO has been implemented in a number of quantum chemistry or semi-empirical codes such as ADF, GAMESS-US, Gaussian, MOPAC, NWChem, TURBOMOLE, and Q-Chem. A COSMO version of the polarizable continuum model PCM has also been developed [ citation needed ]. Depending on the implementation, the details of the cavity construction and the used radii, the segments representing the molecule surface and the value for the dielectric scaling function may vary - which at times causes problems regarding the reproducibility of published results.

Comparison with other methods

While models based on the multipole expansion of the charge distribution of a molecule are limited to small, quasi-spherical or ellipsoidal molecules, the COSMO method has the advantage (as many other dielectric continuum models) that it can be applied to large and irregularly formed molecular structures.

In contrast to the polarizable continuum model (PCM), which uses the exact dielectric boundary conditions, the COSMO method uses the approximative scaling function . Though the scaling is an approximation, it turned out to provide a more accurate description of the so-called outlying charge, reducing the corresponding error. A method comparison [7] of COSMO and the integral equation formalism PCM (IEFPCM), which combines the exact dielectric boundary conditions with a reduced outlying charge error, showed that the differences between the methods are small as compared to deviations to experimental solvation data. The errors introduced by treating a solvent as a continuum and thus neglecting effects like hydrogen bonding or reorientation are thus more relevant to reproduce experimental data than the details of the different continuum solvation methods.

See also

Related Research Articles

<span class="mw-page-title-main">Solvation</span> Association of molecules of a solvent with molecules or ions of a solute

Solvation describes the interaction of a solvent with dissolved molecules. Both ionized and uncharged molecules interact strongly with a solvent, and the strength and nature of this interaction influence many properties of the solute, including solubility, reactivity, and color, as well as influencing the properties of the solvent such as its viscosity and density. If the attractive forces between the solvent and solute particles are greater than the attractive forces holding the solute particles together, the solvent particles pull the solute particles apart and surround them. The surrounded solute particles then move away from the solid solute and out into the solution. Ions are surrounded by a concentric shell of solvent. Solvation is the process of reorganizing solvent and solute molecules into solvation complexes and involves bond formation, hydrogen bonding, and van der Waals forces. Solvation of a solute by water is called hydration.

<span class="mw-page-title-main">Permittivity</span> Measure of the electric polarizability of a dielectric

In electromagnetism, the absolute permittivity, often simply called permittivity and denoted by the Greek letter ε (epsilon), is a measure of the electric polarizability of a dielectric. A material with high permittivity polarizes more in response to an applied electric field than a material with low permittivity, thereby storing more energy in the material. In electrostatics, the permittivity plays an important role in determining the capacitance of a capacitor.

<span class="mw-page-title-main">Molecular dynamics</span> Computer simulations to discover and understand chemical properties

Molecular dynamics (MD) is a computer simulation method for analyzing the physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a fixed period of time, giving a view of the dynamic "evolution" of the system. In the most common version, the trajectories of atoms and molecules are determined by numerically solving Newton's equations of motion for a system of interacting particles, where forces between the particles and their potential energies are often calculated using interatomic potentials or molecular mechanical force fields. The method is applied mostly in chemical physics, materials science, and biophysics.

<span class="mw-page-title-main">Molecular mechanics</span> Use of classical mechanics to model molecular systems

Molecular mechanics uses classical mechanics to model molecular systems. The Born–Oppenheimer approximation is assumed valid and the potential energy of all systems is calculated as a function of the nuclear coordinates using force fields. Molecular mechanics can be used to study molecule systems ranging in size and complexity from small to large biological systems or material assemblies with many thousands to millions of atoms.

<span class="mw-page-title-main">Molecular modelling</span> Discovering chemical properties by physical simulations

Molecular modelling encompasses all methods, theoretical and computational, used to model or mimic the behaviour of molecules. The methods are used in the fields of computational chemistry, drug design, computational biology and materials science to study molecular systems ranging from small chemical systems to large biological molecules and material assemblies. The simplest calculations can be performed by hand, but inevitably computers are required to perform molecular modelling of any reasonably sized system. The common feature of molecular modelling methods is the atomistic level description of the molecular systems. This may include treating atoms as the smallest individual unit, or explicitly modelling protons and neutrons with its quarks, anti-quarks and gluons and electrons with its photons.

Polarizability usually refers to the tendency of matter, when subjected to an electric field, to acquire an electric dipole moment in proportion to that applied field. It is a property of all matter, considering that matter is made up of elementary particles which have an electric charge, namely protons and electrons. When subject to an electric field, the negatively charged electrons and positively charged atomic nuclei are subject to opposite forces and undergo charge separation. Polarizability is responsible for a material's dielectric constant and, at high (optical) frequencies, its refractive index.

Amsterdam Density Functional (ADF) is a program for first-principles electronic structure calculations that makes use of density functional theory (DFT). ADF was first developed in the early seventies by the group of E. J. Baerends from the Vrije Universiteit in Amsterdam, and by the group of T. Ziegler from the University of Calgary. Nowadays many other academic groups are contributing to the software. Software for Chemistry & Materials (SCM), formerly known as Scientific Computing & Modelling is a spin-off company from the Baerends group. SCM has been coordinating the development and distribution of ADF since 1995. Together with the rise in popularity of DFT in the nineties, ADF has become a popular computational chemistry software package used in the industrial and academic research. ADF excels in spectroscopy, transition metals, and heavy elements problems. A periodic structure counterpart of ADF named BAND is available to study bulk crystals, polymers, and surfaces. The Amsterdam Modeling Suite has expanded beyond DFT since 2010, with the semi-empirical MOPAC code, the Quantum_ESPRESSO plane wave code, a density-functional based tight binding (DFTB) module, a reactive force field module ReaxFF, and an implementation of Klamt's COSMO-RS method, which also includes COSMO-SAC, UNIFAC, and QSPR.

<span class="mw-page-title-main">Force field (chemistry)</span> Concept on molecular modeling

In the context of chemistry and molecular modelling, a force field is a computational method that is used to estimate the forces between atoms within molecules and also between molecules. More precisely, the force field refers to the functional form and parameter sets used to calculate the potential energy of a system of atoms or coarse-grained particles in molecular mechanics, molecular dynamics, or Monte Carlo simulations. The parameters for a chosen energy function may be derived from experiments in physics and chemistry, calculations in quantum mechanics, or both. Force fields are interatomic potentials and utilize the same concept as force fields in classical physics, with the difference that the force field parameters in chemistry describe the energy landscape, from which the acting forces on every particle are derived as a gradient of the potential energy with respect to the particle coordinates.

The Poisson–Boltzmann equation is a useful equation in many settings, whether it be to understand physiological interfaces, polymer science, electron interactions in a semiconductor, or more. It aims to describe the distribution of the electric potential in solution in the direction normal to a charged surface. This distribution is important to determine how the electrostatic interactions will affect the molecules in solution. The Poisson–Boltzmann equation is derived via mean-field assumptions. From the Poisson–Boltzmann equation many other equations have been derived with a number of different assumptions.

Implicit solvation is a method to represent solvent as a continuous medium instead of individual “explicit” solvent molecules, most often used in molecular dynamics simulations and in other applications of molecular mechanics. The method is often applied to estimate free energy of solute-solvent interactions in structural and chemical processes, such as folding or conformational transitions of proteins, DNA, RNA, and polysaccharides, association of biological macromolecules with ligands, or transport of drugs across biological membranes.

<span class="mw-page-title-main">Water model</span>

In computational chemistry, a water model is used to simulate and thermodynamically calculate water clusters, liquid water, and aqueous solutions with explicit solvent. The models are determined from quantum mechanics, molecular mechanics, experimental results, and these combinations. To imitate a specific nature of molecules, many types of models have been developed. In general, these can be classified by the following three points; (i) the number of interaction points called site, (ii) whether the model is rigid or flexible, (iii) whether the model includes polarization effects.

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 and the conductor-like PCM (C-PCM) in which the continuum is conductor-like similar to COSMO Solvation Model.

In chemistry, ion association is a chemical reaction whereby ions of opposite electric charge come together in solution to form a distinct chemical entity. Ion associates are classified, according to the number of ions that associate with each other, as ion pairs, ion triplets, etc. Ion pairs are also classified according to the nature of the interaction as contact, solvent-shared or solvent-separated. The most important factor to determine the extent of ion association is the dielectric constant of the solvent. Ion associates have been characterized by means of vibrational spectroscopy, as introduced by Niels Bjerrum, and dielectric-loss spectroscopy.

Biology Monte Carlo methods (BioMOCA) have been developed at the University of Illinois at Urbana-Champaign to simulate ion transport in an electrolyte environment through ion channels or nano-pores embedded in membranes. It is a 3-D particle-based Monte Carlo simulator for analyzing and studying the ion transport problem in ion channel systems or similar nanopores in wet/biological environments. The system simulated consists of a protein forming an ion channel (or an artificial nanopores like a Carbon Nano Tube, CNT), with a membrane (i.e. lipid bilayer) that separates two ion baths on either side. BioMOCA is based on two methodologies, namely the Boltzmann transport Monte Carlo (BTMC) and particle-particle-particle-mesh (P3M). The first one uses Monte Carlo method to solve the Boltzmann equation, while the later splits the electrostatic forces into short-range and long-range components.

<span class="mw-page-title-main">Electric dipole moment</span> Measure of positive and negative charges

The electric dipole moment is a measure of the separation of positive and negative electrical charges within a system, that is, a measure of the system's overall polarity. The SI unit for electric dipole moment is the coulomb-meter (C⋅m). The debye (D) is another unit of measurement used in atomic physics and chemistry.

The reaction field method is used in molecular simulations to simulate the effect of long range dipole-dipole interactions for simulations with periodic boundary conditions. Around each molecule there is a 'cavity' or sphere within which the Coulomb interactions are treated explicitly. Outside of this cavity the medium is assumed to have a uniform dielectric constant. The molecule induces polarization in this media which in turn creates a reaction field, sometimes called the Onsager reaction field. Although Onsager's name is often attached to the technique, because he considered such a geometry in his theory of the dielectric constant, the method was first introduced by Barker and Watts in 1973.

The Born equation can be used for estimating the electrostatic component of Gibbs free energy of solvation of an ion. It is an electrostatic model that treats the solvent as a continuous dielectric medium.

COSMO-RS is a quantum chemistry based equilibrium thermodynamics method with the purpose of predicting chemical potentials µ in liquids. It processes the screening charge density σ on the surface of molecules to calculate the chemical potential µ of each species in solution. Perhaps in dilute solution a constant potential must be considered. As an initial step a quantum chemical COSMO calculation for all molecules is performed and the results are stored in a database. In a separate step COSMO-RS uses the stored COSMO results to calculate the chemical potential of the molecules in a liquid solvent or mixture. The resulting chemical potentials are the basis for other thermodynamic equilibrium properties such as activity coefficients, solubility, partition coefficients, vapor pressure and free energy of solvation. The method was developed to provide a general prediction method with no need for system specific adjustment.

DMol3 is a commercial software package which uses density functional theory with a numerical radial function basis set to calculate the electronic properties of molecules, clusters, surfaces and crystalline solid materials from first principles. DMol3 can either use gas phase boundary conditions or 3D periodic boundary conditions for solids or simulations of lower-dimensional periodicity. It has also pioneered the use of the conductor-like screening model COSMO Solvation Model for quantum simulations of solvated molecules and recently of wetted surfaces. DMol3 permits geometry optimisation and saddle point search with and without geometry constraints, as well as calculation of a variety of derived properties of the electronic configuration. DMol3 development started in the early eighties with B. Delley then associated with A.J. Freeman and D.E. Ellis at Northwestern University. In 1989 DMol3 appeared as DMol, the first commercial density functional package for industrial use by Biosym Technologies now Accelrys. Delley's 1990 publication was cited more than 3000 times.

In computational chemistry, a solvent model is a computational method that accounts for the behavior of solvated condensed phases. Solvent models enable simulations and thermodynamic calculations applicable to reactions and processes which take place in solution. These include biological, chemical and environmental processes. Such calculations can lead to new predictions about the physical processes occurring by improved understanding.

References

  1. 1 2 3 4 A., Klamt; G., Schüürmann (1993). "COSMO: a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient". J. Chem. Soc. Perkin Trans.2. 2 (5): 799–805. doi:10.1039/P29930000799.
  2. 1 2 Klamt, Andreas (2005). From Quantum Chemistry to Fluid Phase Thermodynamics and Drug Design. Boston, MA, USA: Elsevier. ISBN   9780444519948.
  3. 1 2 Herbert, John M. (2021-03-23). "Dielectric continuum methods for quantum chemistry". WIREs Computational Molecular Science. 11 (4). arXiv: 2203.06846 . doi:10.1002/wcms.1519. ISSN   1759-0876. S2CID   233629977.
  4. Cramer, Christopher J. (2004). Essentials of computational chemistry: theories and models (2nd ed.). Chichester, West Sussex, England: Wiley. ISBN   0-470-09182-7. OCLC   55887497.
  5. Frank, Jensen (2017). Introduction to computational chemistry. John Wiley & Sons. ISBN   978-1-118-82599-0. OCLC   989360916.
  6. Hoshi, Hajime; Sakurai, Minoru; Inoue, Yoshio; Chûjô, Riichirô (1987-07-15). "Medium effects on the molecular electronic structure. I. The formulation of a theory for the estimation of a molecular electronic structure surrounded by an anisotropic medium". The Journal of Chemical Physics. 87 (2): 1107–1115. Bibcode:1987JChPh..87.1107H. doi:10.1063/1.453343. ISSN   0021-9606.
  7. Klamt, A.; Moya, C.; Palomar, J. (2015). "A Comprehensive Comparison of the IEFPCM and SS(V)PE Continuum Solvation Methods with the COSMO Approach". Journal of Chemical Theory and Computation. 11 (9): 4220–4225. doi:10.1021/acs.jctc.5b00601. PMID   26575917.