News
News

\(drf90\)

classical-only part of the Direct (or Discrete) Reaction Field (DRF) approach
Current version: 2012.01 (shared with CC-BY-NC license). See also all python tools and the drf90 pages.

DRF block/namelist for general DRF parameters

to be updated soon
parameter
description
default
\(afact\)
A-factor used in Thole's model
2.1304
\(direct\)
Parameter specifying that the Direct-method should be used for solving the relay equations. This means that the Relay-matrix will NOT be kept in memory, but everytime an element of it is needed, it will be calculated on-the-fly.
-
\(ibem\)
Type of boundary used around system:
  • 0: Do not use boundary
  • 5: Use Boundary Elements for solving Poisson-Boltzmann equations, either by reading Connolly surface from 'csurf'-file, or by constructing a GEPOL93 surface (see also GEPOL-block). This depends on the NOGEPOL parameter.
0
\(idrfout\)
Parameter describing amount of output generated:
  • 1: Normal printing
  • 2: Increased number of lines
  • 3: Maximum output (for debugging purposes)
1
\(ifldqm\)
Include field from QM region (for coupling with Molcas)
  • 0: MASSIETODO
  • 1: MASSIETODO
0
\(iljcomb\)
Combination of vdW parameters taken from LennardJones block: TODOMASSIE
  • 1: use geometric mean for combining radius parameter, \(R_{ij} = \sqrt{ R_{i} · R_{j} }\)
  • 2: use arithmetic mean for combining radius parameter, \(R_{ij} = \dfrac{ R_{i} + R_{j}}{2} \)
  • 3: use summation for combining radius parameter, \(R_{ij} = R_{i} + R_{j} \)
In all cases is the energy combined with geometric mean.
1
\(imdc\)
Redistribute induced dipoles by assigning fractional atomic charges:
  • 0: MASSIETODO
  • 1: MASSIETODO
0
\(iradex\)
Type of atomic radii used:
  • 4: Frecer's charge dependent radii (H-Ar)
  • 5: Miertus's charge dependent radii (H-Ar)

Note:

  • This option will always overrule any radius specified in the EXTERNAL block
  • If no radius is available for an atom with option 4/5, then if a nonzero radius is given in the EXTERNAL block, that one will be used, else a tabulated value will be used (Bondii for some atoms; -99.999 if unknown ⇒ this will make the program stop, and you will have to specify a radius in the EXTERNAL block)
4
\(isodis\)
Isotropic or anisotropic dispersion energy:
  • 0: Use anisotropic dispersion
  • 1: Use isotropic dispersion
0
\(isorep\)
Isotropic or anisotropic repulsion energy:
  • 0: Use anisotropic repulsion
  • 1: Use isotropic repulsion
1
\(itermax\)
Maximum number of iterations for iterative solution of relay equations
100
\(ithole\)
Type of screening function used in Thole's model:
  • 0: No screening function
  • 1: Conical function
  • 2: Exponential function
  • TODOMASSIE
2
\(ithole\_fld\)
Same as ITHOLE (see above), but only for inducing fields
same as choice for \(ithole\)
\(iunit\)
Atomic coordinates unit on input:
  • 0: Atomic units (Bohrs)
  • 1: Angstroms
0
\(ivdwcomb\)
Combination of vdW parameters taken from (UFF, MM3) database:
  • 1: use geometric mean for combining radius parameter, \(R_{ij} = \sqrt{ R_{i} · R_{j} }\)
  • 2: use arithmetic mean for combining radius parameter, \(R_{ij} = \dfrac{ R_{i} + R_{j}}{2} \)
  • 3: use summation for combining radius parameter, \(R_{ij} = R_{i} + R_{j} \)
In all cases is the energy combined with geometric mean.
different for different options of \(i\_ener\_rep\)
DRF block/namelist for general DRF parameters
Parameter Description Default value
AFACT
A-factor used in Thole's model 2.1304
CVGREL
Convergence criterium for solving relay equations 1.0.10-6
DIRECT
Parameter specifying that the Direct-method should be used for solving the relay equations. This means that the Relay-matrix will NOT be kept in memory, but everytime an element of it is needed, it will be calculated on-the-fly. -
IBEM
Type of boundary used around system.
0 : Do not use boundary
5 : Use Boundary Elements for solving Poisson-Boltzmann equations,
either by reading Connolly surface from 'csurf'-file,
or by constructing a GEPOL93 surface (see also GEPOL-block).
This depends on the NOGEPOL parameter.
0
IDRFOUT
Parameter describing amount of output generated
1 : Normal printing
2 : Increased number of lines
3 : Maximum (debug) output
1
IFLDQM
Include field from QM region (for coupling with Molcas) 0
ILJCOMB
Combination of vdW parameters taken from LennardJones block:
1: use geometric mean for combining radius paramater (Rij = sqrt[Ri * Rj])
2: use arithmetic mean for combining radius parameter (Rij = 0.5*[Ri + Rj])
3: use summation for combining radius parameter (Rij = [Ri + Rj])
In all cases is the energy combined with geometric mean.
1
IMDC
Redistribute induced dipoles by assigning fractional atomic charges 0
IRADEX
Type of atomic radii used.
4 : Frecer's charge dependent radii (H-Ar)
5 : Miertus's charge dependent radii (H-Ar)

Note:
- This option will always overrule any radius specified in the EXTERNAL block.
- If no radius is available for an atom with option 4/5, then if a nonzero radius is given in the EXTERNAL block, that one will be used, else a tabulated value will be used (Bondii for some atoms; -99.999 if unknown => this will make the program stop, and you will have to specify a radius in the EXTERNAL block).
4
ISODIS
Isotropic or anisotropic dispersion energy
0 : Use anisotropic dispersion
1 : Use isotropic dispersion
0
ISOREP
Isotropic or anisotropic repulsion energy
0 : Use anisotropic repulsion
1 : Use isotropic repulsion
1 (New Default !)
ITERMAX
Maximum number of iterations for iterative solution of relay equations 100
ITHOLE
Type of screening function used in Thole's model
0 : No screening function
1 : Conical function
2 : Exponential function
2
ITHOLE_FLD
Same as ITHOLE (see above), but only for inducing fields same as choice for ITHOLE
IUNIT
Atomic coordinates unit on input
0 : Atomic units (Bohrs)
1 : Angstroms
0
IVDWCOMB
Combination of vdW parameters taken from (UFF, MM3) database:
1: use geometric mean for combining radius paramater (Rij = sqrt[Ri * Rj])
2: use arithmetic mean for combining radius parameter (Rij = 0.5*[Ri + Rj])
3: use summation for combining radius parameter (Rij = [Ri + Rj])
In all cases is the energy combined with geometric mean.
different for different options of i_ener_rep
I_ENER_DIS
Type of dispersion energy:
0: do not calculate dispersion
1: use Slater-Kirkwood type dispersion from atomic polarizability (tensors)
2: use Lennard-Jones potential type, with parameters taken from LennardJones block
4: use LJ potential type, with parameters taken from MM3 database
5: use LJ potential type, with parameters taken from UFF database
6: use LJ potential type, with parameters taken from MM3 database (energy scaled by 0.25)
7: use LJ potential type, with parameters taken from UFF database (energy scaled by 0.25)
1
I_ENER_EST
Type of electrostatic energy:
0: do not calculate electrostatic energy
1: do calculate
1
I_ENER_IND
Type of induction (polarization) energy
0: do not calculate induced dipoles with corresponding induction energy
1: do calculate from (atomic/molecular) polarizability (tensors)
1
I_ENER_REP
Type of short-range repulsion energy
0: do not calculate repulsion energy

Lennard-Jones (rho**[-12], with rho = r / Req) type

1: CHARMM-1983 like, as function of Slater-Kirkwood dispersion, screened for overlapping charge densities
2: Lennard-Jones repulsion, with parameters taken from LennardJones block
3: CHARMM-1983 like (see 1), NOT screened
4: use LJ repulsion, parameters taken from MM3 database
5: use LJ repulsion, parameters taken from UFF database
6: use LJ repulsion, parameters taken from MM3 database (energy parameter scaled by 0.25)
7: use LJ repulsion, parameters taken from UFF database (energy parameter scaled by 0.25)

Exponential (X6) repulsion (B*exp[N(1-rho)]) type

-1: CHARMM-1983 like, as function of Slater-Kirkwood dispersion, screened for overlapping charge densities
-2: X6 repulsion, with parameters taken from LennardJones block
-3: CHARMM-1983 like (see 1), NOT screened
-4: use X6 repulsion, parameters taken from MM3 database
-5: use X6 repulsion, parameters taken from UFF database
-6: use X6 repulsion, parameters taken from MM3 database (energy parameter scaled by 0.25)
-7: use X6 repulsion, parameters taken from UFF database (energy parameter scaled by 0.25)
-3 (New Default !)
NOGEPOL
Use of GEPOL93 for making Boundary Elements (see IBEM) for solving Poisson-equation
0 : Use GEPOL93
1 : Read BEM's from 'csurf' file
0
NSTATE
Number of states to take into account.
If > 1, then the energy of the additional states is included in statistical analysis during molecular simulations.

For each additional state, the force field parameters (charges only for the moment) should be given in a EXTSTATEi block (formatted similar to EXTERNAL block), with i replaced by the number of the additional state (i.e. EXTSTATE2 for i=2, etc.).
1
RFACT
Ad-hoc factor for increasing atomic radii in repulsion. Only to be used for debugging purposes. 1.0
TOTALALPHA
Calculates total "molecular" polarizability of the complete system, if present -
X6ZETA
Normalization factor N for X6 repulsion parameter.
13.772