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
|
\(idrfout\)
|
Parameter describing amount of output generated:
|
1
|
\(ifldqm\)
|
Include field from QM region (for coupling with Molcas)
|
0
|
\(iljcomb\)
|
Combination of vdW parameters taken from LennardJones block: TODOMASSIE
|
1
|
\(imdc\)
|
Redistribute induced dipoles by assigning fractional atomic charges:
|
0
|
\(iradex\)
|
Type of atomic radii used:
Note:
|
4
|
\(isodis\)
|
Isotropic or anisotropic dispersion energy:
|
0
|
\(isorep\)
|
Isotropic or anisotropic repulsion energy:
|
1
|
\(itermax\)
|
Maximum number of iterations for iterative solution of relay equations
|
100
|
\(ithole\)
|
Type of screening function used in Thole's model:
|
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
|
\(ivdwcomb\)
|
Combination of vdW parameters taken from (UFF, MM3) database:
|
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 | |||||||||
IDRFOUT |
Parameter describing amount of output generated
|
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.
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 | |||||||||
ISOREP |
Isotropic or anisotropic repulsion energy
|
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
|
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 | |||||||||
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 | |||||||||
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 |