COSMO: Conductor like Screening Model and the Solvation-key¶

You can study chemistry in solution, as contrasted to the gas phase, with the implementation in BAND of the Conductor like Screening Model (COSMO) of solvation [1].

In the COSMO model all solvents are roughly the same, and approximated by an enveloping metal sheet. One explicit dependency on the solvent is that the solvation energy is scaled by

and this depends on the dielectric constant of the solvent, and an empirical factor \(\chi\). The other is that the shape of the surface is influenced by the Rad parameter, see below.

The solvent information is specified in the solvent key of the solvation block. The simplest option is to use one of the pre-defined solvents:

Solvation
   Enabled Yes/No
   Solvent
      Name [...]
   End
End

Solvation

Type:

Block

Description:

Options for the COSMO (Conductor like Screening Model) solvation model.

Enabled

Type:

Bool

Default value:

No

GUI name:

Include COSMO solvation

Description:

Use the Conductor like Screening Model (COSMO) to include solvent effects.

Solvent

Type:

Block

Description:

Solvent details

Name

Type:

Multiple Choice

Default value:

Water

Options:

[AceticAcid, Acetone, Acetonitrile, Ammonia, Aniline, Benzene, BenzylAlcohol, Bromoform, Butanol, isoButanol, tertButanol, CarbonDisulfide, CarbonTetrachloride, Chloroform, Cyclohexane, Cyclohexanone, Dichlorobenzene, DiethylEther, Dioxane, DMFA, DMSO, Ethanol, EthylAcetate, Dichloroethane, EthyleneGlycol, Formamide, FormicAcid, Glycerol, HexamethylPhosphoramide, Hexane, Hydrazine, Methanol, MethylEthylKetone, Dichloromethane, Methylformamide, Methypyrrolidinone, Nitrobenzene, Nitrogen, Nitromethane, PhosphorylChloride, IsoPropanol, Pyridine, Sulfolane, Tetrahydrofuran, Toluene, Triethylamine, TrifluoroaceticAcid, Water]

GUI name:

Solvent

Description:

Name of a pre-defined solvent. A solvent is characterized by the dielectric constant (Eps) and the solvent radius (Rad).

This is the list of possible solvents and their corresponding Eps and Rad values:

Several other options can be defined in the Solvation block:

Solvation
   CVec [EXACT | FITPOT]
   Charge
      Conv float
      Corr Yes/No
      Iter integer
      Method [CONJ | INVER]
   End
   Enabled Yes/No
   Radii # Non-standard block. See details.
      ...
   End
   SCF [VAR | PERT | NONE]
   Solvent
      Del float
      Emp float
      Eps float
      Name [...]
      Rad float
   End
   Surf [Delley | Wsurf | Asurf | Esurf | Klamt]
End

Solvation

Type:

Block

Description:

Options for the COSMO (Conductor like Screening Model) solvation model.

CVec

Type:

Multiple Choice

Default value:

EXACT

Options:

[EXACT, FITPOT]

GUI name:

Calculate Coulomb interaction

Description:

Choose how to calculate the Coulomb interaction matrix between the molecule and the point charges on the surface: - EXACT: use exact density, and integrate against the potential of the point charges. This may have inaccuracies when integration points are close to the point charges. - FITPOT: evaluate the molecular potential at the positions of the point charges, and multiply with these charges.

Charge

Type:

Block

Description:

Select the algorithm to determine the charges.

Conv

Type:

Float

Default value:

1e-08

Description:

Charge convergence threshold in iterative COSMO solution.

Corr

Type:

Bool

Default value:

Yes

GUI name:

Correct for outlying charge

Description:

Correct for outlying charge.

Iter

Type:

Integer

Default value:

1000

Description:

Maximum number of iterations to solve COSMO equations.

Method

Type:

Multiple Choice

Default value:

CONJ

Options:

[CONJ, INVER]

GUI name:

Charge determination method

Description:

INVER: matrix inversion, CONJ: biconjugate gradient method. The CONJ method is guaranteed to converge with small memory requirements and is normally the preferred method.

Enabled

Type:

Bool

Default value:

No

GUI name:

Include COSMO solvation

Description:

Use the Conductor like Screening Model (COSMO) to include solvent effects.

Radii

Type:

Non-standard block

Description:

The values are the radii of the atomic spheres. If not specified the default values are those by Allinger. Format: ‘AtomType value’. e.g.: ‘H 0.7’

SCF

Type:

Multiple Choice

Default value:

VAR

Options:

[VAR, PERT, NONE]

GUI name:

Handle charges

Description:

Determine the point charges either Variational (VAR) or after the SCF as a Perturbation (PERT).

Solvent

Type:

Block

Description:

Solvent details

Del

Type:

Float

Description:

Del is the value of Klamt’s delta_sol parameter, only relevant in case of Klamt surface.

Emp

Type:

Float

Description:

Emp is the empirical scaling factor x for the energy scaling.

Eps

Type:

Float

Description:

User-defined dielectric constant of the solvent (overrides the Eps value of the solvent defined in ‘Name’)

Name

Type:

Multiple Choice

Default value:

Water

Options:

[AceticAcid, Acetone, Acetonitrile, Ammonia, Aniline, Benzene, BenzylAlcohol, Bromoform, Butanol, isoButanol, tertButanol, CarbonDisulfide, CarbonTetrachloride, Chloroform, Cyclohexane, Cyclohexanone, Dichlorobenzene, DiethylEther, Dioxane, DMFA, DMSO, Ethanol, EthylAcetate, Dichloroethane, EthyleneGlycol, Formamide, FormicAcid, Glycerol, HexamethylPhosphoramide, Hexane, Hydrazine, Methanol, MethylEthylKetone, Dichloromethane, Methylformamide, Methypyrrolidinone, Nitrobenzene, Nitrogen, Nitromethane, PhosphorylChloride, IsoPropanol, Pyridine, Sulfolane, Tetrahydrofuran, Toluene, Triethylamine, TrifluoroaceticAcid, Water]

GUI name:

Solvent

Description:

Name of a pre-defined solvent. A solvent is characterized by the dielectric constant (Eps) and the solvent radius (Rad).

Rad

Type:

Float

Unit:

Angstrom

Description:

User-defined radius of the solvent molecule (overrides the Rad value of the solvent defined in ‘Name’).

Surf

Type:

Multiple Choice

Default value:

Delley

Options:

[Delley, Wsurf, Asurf, Esurf, Klamt]

GUI name:

Surface type

Description:

Within the COSMO model the molecule is contained in a molecule shaped cavity. Select one of the following surfaces to define the cavity: - Wsurf: Van der Waals surface - Asurf: solvent accessible surface - Esurf: solvent excluding surface - Klamt: Klamt surface - Delley: Delley surface.

Additional keys for periodic systems¶

For the simulation of periodic structures ICW solvation, you may specify the following options:

PeriodicSolvation
   RemovePointsWithNegativeZ Yes/No
   NStar integer
End

PeriodicSolvation

Type:

Block

Description:

Additional options for simulations of periodic structures with solvation.

RemovePointsWithNegativeZ

Type:

Bool

Default value:

No

GUI name:

Only above slab

Description:

Whether the COSMO surface is constructed on both sides of a surface. If one is only interested in the solvation effect on the upper side of a surface (in the Z direction), then this option should be set to ‘True’

NStar

Type:

Integer

Default value:

4

Description:

This option, expecting an integer number (>2), handles the accuracy for the construction of the COMSO surface. The larger the given number the more accurate the construction.

General remarks: The accuracy of the result and the calculation time is influenced by the screening radius SCREENING%RMADEL (see Screening block). If the calculation does take too long, defining a smaller radius does help. But: too small radii, especially smaller than the lattice constants, will give unphysical results.

References