ICONST

From VASP Wiki
Revision as of 09:54, 9 March 2022 by Tbucko (talk | contribs)

Geometric parameters that are monitored or controlled in molecular dynamics (e.g. constrained or affected by the action of a bias potential) are defined in the file ICONST. Two kinds of geometric parameters can be defined: primitive (such as bond lengths or angles), and complex (e.g. linear combinations of primitive coordinates). Each coordinate is defined in a separate line, the complex coordinates must be defined after primitive ones. In order to define a primitive coordinate, the following syntax is used: 

FLAG item(1) ... item(N) STATUS


where FLAG is a character used to define a type of primitive coordinate:

  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): R - interatomic distance between atoms item(1) and item(2).
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): A - angle defined by atoms item(1), item(2) and item(3) (with the atom item(2) being the apex).
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): T - torsional angle defined by atoms item(1), item(2), item(3) and item(4).
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): M - distance between atom item(1) and the center of bond between atoms item(2) and item(3).
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): B - distance between the center of bond between atoms item(1) and item(2) and the center of bond between atoms item(3) and item(4)
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): {\displaystyle P} - ratio of length of the bond between atoms item(1) and item(2) and the length of the bond between atoms item(3) and item(4)
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): W - function Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): {\displaystyle \frac{1-\left(R/c\right)^M}{1-\left(R/c\right)^{N}}} with Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): R being the bond length between the atoms item(1) and item(2), Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): c is the reference bond length specified as item(3), and the exponents Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): M and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): N are defined as item(4) and item(5), respectively.
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): X,Y,Z - fractional (direct) coordinates linked with the lattice vectors Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): a , Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): b , and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): c .
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): {\displaystyle cX, cY, cZ} - Cartesian coordinates Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): x , Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): y , and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): z .
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): {\displaystyle LR} - length of lattice vector item(1)
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): {\displaystyle LA} - angle between lattice vectors item(1) and item(2)
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): {\displaystyle LV} - cell volume (no item(i) is defined in this case)

In most cases, item(i) is an integer specifying either the position of the atom or the lattice vector in the POSCAR file (obviously, two atoms are needed to define a bond length, three atoms are required for a bonding angle, etc...). In the special case of the coordinate Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): W , also some additional parameters - the reference bond length (generally a floating point number) and exponents used in the definition (integers) - are defined via item(i), see the corresponding description above. STATUS is an integer distinguishing between the constraint (STATUS=0), the coordinate defining the collective variable in metadynamics (STATUS=5), the coordinate affected by Fermi-type step (STATUS=4) or harmonic potential (STATUS=8) or the monitored (but otherwise unaffected) coordinate (STATUS=7). The following example shows the ICONST file specifying two constraints - bond lengths between the atoms 1 and 5, and between the atoms 1 and 6: 

R 1 5 0
R 1 6 0

The complex coordinates are functions defined in the space of primitive coordinates. All complex coordinates must be defined after the last primitive coordinate. Assuming that Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): M primitive coordinates (Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): q_{i} ) were specified on the first Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): M lines of the ICONST file, the following syntax is used to define a complex coordinate: \\ 

FLAG  c_1  c_2  ...  c_M  STATUS

where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): c_{i} is a coefficient for the primitive coordinate defined at the line Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): i . The number of coefficients must be the same as the number of primitive coordinates. The following types of complex coordinates are supported:

  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): S - linear combination of primitive coordinates, i.e., Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): \left(\xi =\sum _{{i=1}}^{{M}}c_{i}\,q_{i}\right)
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): C - norm of vector of primitive coordinates, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): \left(\xi ={\sqrt {\sum _{{i=1}}^{{M}}\,(c_{i}\,q_{i})^{2}}}\right)
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): D - coordination number Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): \left(\xi =\sum _{{i=1}}^{{M}}{\frac {1-\left(q_{{i}}/c_{{i}}\right)^{9}}{1-\left(q_{{i}}/c_{{i}}\right)^{{14}}}}\right)

As in the case of primitive coordinates, STATUS allows to distinguish between the geometric constraint, the coordinate affected by bias potentials and the monitored coordinate. Whenever complex coordinates are defined, the primitives are used only as a basis for their definition. Consider, for instance, the ICONST file with the following lines: 

R 1 6 0
R 1 5 0
S 1 -1 0

The first two lines define two primitive coordinates - bonds between the atoms Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): 1 and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): 6 , and between the atoms Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): 1 and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): 5 . The complex coordinate - difference between the two bond lengths - is defined on the third line. Consequently, the two primitive coordinates are not constrained in the simulation (despite STATUS=0), the only controlled parameter is the complex coordinate. Clearly, in order to fix the first bond length and the complex coordinate at the same time, the ICONST file should be modified as follows: 

R 1 6 0
R 1 5 0
S 1 -1 0
S 1 0 0

In the context of MD in NpT ensemble, ICONST can be used to impose restrictions on the volume and/or shape of the simulation cell. The following examples cover some of the most important scenarios: 1.) Simulation with a constant cell volume LV 0 2.) Variable lengths of lattice vectors but fixed lattice angles LA 1 2 0 LA 1 3 0 LA 2 3 0 3.) A cubic cell with variable volume but fixed shape LA 1 2 0 LA 1 3 0 LA 2 3 0 LR 1 0 LR 2 0 LR 3 0 S 1 0 0 0 0 0 0 S 0 1 0 0 0 0 0 S 0 0 1 0 0 0 0 S 0 0 0 1 -1 0 0 S 0 0 0 1 0 -1 0 S 0 0 0 0 1 -1 0 4.) An orthorhombic cell variable volume but fixed shape. Here, in order to fix the ratios between lengths of lattice vectors (Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): {\displaystyle a_1:a_2:a_3} ) we define constraints of the form Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): {\displaystyle c_i*a_i + c_j*a_j = 0} . For instance, if the shape of cell is such that Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): {\displaystyle a_1:a_2:a_3=1:1.5:2} , the following ICONST can be used: LA 1 2 0 LA 1 3 0 LA 2 3 0 LR 1 0 LR 2 0 LR 3 0 S 1 0 0 0 0 0 0 S 0 1 0 0 0 0 0 S 0 0 1 0 0 0 0 S 0 0 0 1.5 -1.0 0.0 0 S 0 0 0 2.0 0 -1.0 0 S 0 0 0 0 3.0 -2.0 0

Related Tags and Sections

MDALGO

Retrieved from "https://www.vasp.at/wiki/index.php?title=ICONST&oldid=14507"