Next: Output files
Up: Metadynamics
Previous: The Implemented Types of
Contents
Index
- ANALYSIS: A standard MD run is performed, where the equations of motions are not
affected by the hills-potential or the coupling potential.
The selected CV are monitored and the values are reported in the
output file, after every 10 MD steps.
This option is useful in order to observe the behavior of the selected CV in equilibrium conditions.
With this option only two output files are written: istvar_mtd, and enevar_mtd (see section
9.11.6). The former file contains the values of the
and their averages in time.
- METASTEPNUM: The maximum number of MTD steps is read from the
next line,
(default: 100).
- META_RESTART:
To restart a metadynamic's run from where the former run has stopprd, one can use this keyword, and write in the
following line, the number of meta-steps completed already
(default: 0).
Beware that for restarting the MTD in this way, the output
files of the previous run are to be available in the run's directory and must contain a number
of lines at least equal to the
.
From this files the previous history of the MTD is read and the MTD is initialized accordingly.
Otherwise, it is possible to restart from the restart file of the MTD,
.
This is an unformatted file, which is written whenever the standard CPMD restart file,
, is also written.
It contains the number of meta-steps already performed, the number of CV used in the previous run and the information
about the position and the size of the hills which have been already located.
To restart the MTD from this file, the same keyword is used, and the keyword
is added in the same line,
providing that the unformatted restart file is available in the run's directory.
In this second case the number of performed meta-step is not read from the input file but from the restart file
.
Obviously, when a run is restarted, the same number and the same kinds of CV must be used. However masses, force
constants, scaling factors, and the width and height of the hills can be changed.
- MINSTEPNUM INTERMETA: The minimum number of MD steps between two
MTD steps (in general, the MTD step is characterized by the positioning
of a new hill in the CV-space) is read from the next line,
(default: 100). This is a lower bound because,
before the construction of a new hill, the displacement in the
CV-space is checked, and the new step is accepted only if the
calculated displacement is above a given tolerance.
- MOVEMENT CHECK: The tolerance for the acceptance of a
new MTD step is read from the next line (default: 0.001D0)
- CHECK DELAY: The number of MD steps to be run, before a
new check of the displacement is done, is read from the following line (default: 20).
- MAXSTEPNUM INTERMETA: The maximum number of MD steps that can be run,
before a new MTD step is accepted anyway, is read from the following line (default: 300).
- METASTORE [NO TRAJECTORY]:
In the following line, three integer numbers are given,
which indicate respectively how often (in terms of MTD steps)
the
and the
files are over-written (default: 50),
how often the trajectory files are appended (default: 1),
and how often a quench of the electronic wavefunctions onto the
BO surface is performed (default: 10000).
With the additional flag NO TRAJECTORY, the trajectories are
still written according to the settings as indicated by the
TRAJECTORY keyword in the &CPMD section.
The selection of the files (e.g. turning on TRAJEC.xyz via the XYZ flag
and turning TRAJECOTORY off via a negative value of NTRAJ in combination
with the SAMPLE flag) is always honored.
- MAXKINEN: From the following line, the maximum electronic
kinetic energy is given, above which a quench of the
electronic wavefunctions onto the BO surface is performed
anyway (by default no quench is done whatever is the kinetic energy).
- LAGRANGE TEMPERATURE: The temperature
used to initialize
the velocities of the
CV is read from the next line.
By default
is chosen equal to the temperature of the ions,
the units are Kelvin. Notice that this keyword causes the
initialization of the Lagrangian formulation of MTD.
- LAGRANGE TEMPCONTROL: The control of the
is activated,
and the rescaling velocities algorithm is used.
The average temperature and the permitted range of variation are
read from the next line. By default
is not controlled.
Notice that this keyword causes the initialization of the extended degrees of freedom of the Lagrangian formulation of MTD.
- LAGRANGE LANGEVIN: Performs Langevin dynamics for the Lagrangian
formulation of MTD. In the next line the Temperature (Kelvin) and the
friction
(a.u.) are read. The Langevin equation in its standard
form writes (
is a CV):
 |
|
|
(272) |
where
is the CV mass,
a generic force field,
the friction
coefficient,
the temperature and
is a white noise. The integration
algorithm is a Velocity-Verlet which can be written as:
the
are independent Gaussian random numbers with mean zero and
variance one.
- HILLS: With this keyword it is defined the shape of
.
If
is read in the same line, no hill potential is used.
The default hills' shape is the Gaussian-like one described above.
If
is read from the same line, the second exponential term in equation 263
is not applied. i.e., a normal Gaussian function rather than a Gaussian tube formalism is used.
If, instead,
is read in the same line as
,
Lorentzian functions are used in place of Gaussians. If, instead,
is read in the same line as
, the rational function
described in the previous section is used; in this case, if
is read,
the exponents
and
, and the boosting factor
are also
read immediately after (defaults:
,
,
).
If, instead,
is read on the same line as
,
the shifted Gaussians are used, where the tails after a given
cutoff are set equal to zero; in this case, if
is read,
the cutoff
and the boosting factor
are also read
immediately after (defaults:
,
).
In all this cases, if the symbol
is read in the same
line as
, the perpendicular width
and the height
are read immediately after (defaults:
,
).
- TUNING HHEIGHT: With this keyword the height of the hills is
tuned according to the curvature of the underlying potential.
If the symbol
is read in the same line as
,
the lower and upper bounds of
are read immediately after
(defaults:
,
).
- HILLVOLUME: With this keyword the volume of the hills,
, is kept constant
during the MTD run, i.e. when the height changes due to the
tuning (see the previous keyword), the width is changed accordingly.
This option is can be used only if the tuning of the hills' height is active.
- CVSPACE BOUNDARIES: With this keyword the confinement of the
CV-space is required, in the direction of a selected group of CV.
The number of dimensions of the CV-space, in which the confinement
is applied, is read from the next line,
(default: 0).
The following
lines describe the type of confinement.
For each line the following parameters are required: index of the
CV (as from the list given in the definition of the CV),
the strength of the confinement potential
in Hartree, two real numbers,
and
, that determine from which value of the CV the
confining potential is active. Finally, if the keyword
is read in the same line, the real number, which is read
immediately after, determines how smoothly the confining
potential is switched on (default
).
The confining potential can be used with the following CV:
DIST:
and it becomes active only if
.
DIFFER:
and it becomes active only if
or
.
Coordination numbers: if
,
, if
,
.
- RESTRAIN VOLUME: With this keyword a confining potential is
applied to the volume variations. The option can be used only in
combination with the NPE/NPT MD. From the next line, the following parameters are
required:
,
, and
. The factors
and
multiplied by the initial volume of the cell, give, respectively,
the lower and upper bounds for the volume, whereas
gives the
strength of the confining potential.
- MULTI NUM: his keyword should be used when a multiple set of
MTD runs are performed simultaneously on the same system. Here the
number of separated sets CV for each subsystem has to be given in
the following line,
.
This means that the first
CV, in the list of those
defined in the
section, will belong to
the first subset, the next
to the second, and so on.
This option is implemented only together with the extended
Lagrangian formulation.
- MONITOR: This keyword requires that an additional monitoring of the values
of the CV is performed along the MD trajectory. This means that the values are written
on an output file every WCV_FREQ MD steps, even if no hill is added at that step.
The frequency for updating the file is read from the following line, the file name is
cvmdck_mtd, and it is not created if this option is not activated.
Next: Output files
Up: Metadynamics
Previous: The Implemented Types of
Contents
Index
Costas Bekas
2008-09-04