PENALTYPOT: Difference between revisions

From VASP Wiki
No edit summary
 
(3 intermediate revisions by 2 users not shown)
Line 1: Line 1:
At the beginning of each metadynamics simulation, VASP attempts to read the file {{TAG|PENALTYPOT}} containing the bias potential acting on the geometric parameters with STATUS=5 defined in the {{TAG|ICONST}} file. In analogy to the time-dependent bias potential generated in metadynamics, the bias potential is defined as a superposition of Gaussian hills. Each line in the {{TAG|PENALTYPOT}} file represents one (multidimensional) Gaussian:
File that defines the bias potential for a [[:Category:Biased molecular dynamics|biased molecular-dynamics]] run.


<math>x_1 x_2 ... x_m  h w </math>
At the beginning of each metadynamics simulation, VASP attempts to read the file {{TAG|PENALTYPOT}} containing the bias potential acting on the geometric parameters with <code>STATUS=5</code> defined in the {{TAG|ICONST}} file. In analogy to the time-dependent bias potential generated in metadynamics, the bias potential is defined as a superposition of Gaussian hills. Each line in the {{TAG|PENALTYPOT}} file represents one (multidimensional) Gaussian:


where <math>x_1</math> to <math>x_m</math> stand for the position in the space of active coordinates, and <math>h</math> and <math>w</math> are the height and width of the Gaussian, respectively (note that both positive and negative values are allowed for the parameter <math>h</math>). For example, if two active coordinates with {{TAG|STATUS}}=5 are defined in the {{TAG|ICONST}} file:
::<math>x_1 x_2 ... x_m  h w </math>,
 
where <math>x_1</math> to <math>x_m</math> stand for the position in the space of active coordinates, and <math>h</math> and <math>w</math> are the height and width of the Gaussian, respectively. Note that both positive and negative values are allowed for the parameter <math>h</math>.  
 
For example, if two active coordinates with <code>STATUS=5</code> are defined in the {{TAG|ICONST}} file:
  R 1 5 5
  R 1 5 5
  R 1 6 5  
  R 1 6 5  
Line 16: Line 20:


== Related tags and articles ==
== Related tags and articles ==
{{TAG|MDALGO}},{{TAG|ICONST}}
{{TAG|MDALGO}}, {{TAG|ICONST}}


----
----


[[The_VASP_Manual|Contents]]


[[Category:Files]][[Category:Input files]][[Category:Biased molecular dynamics]]
[[Category:Files]][[Category:Input files]][[Category:Molecular dynamics]][[Category:Biased molecular dynamics]]

Latest revision as of 10:52, 18 October 2023

File that defines the bias potential for a biased molecular-dynamics run.

At the beginning of each metadynamics simulation, VASP attempts to read the file PENALTYPOT containing the bias potential acting on the geometric parameters with STATUS=5 defined in the ICONST file. In analogy to the time-dependent bias potential generated in metadynamics, the bias potential is defined as a superposition of Gaussian hills. Each line in the PENALTYPOT file represents one (multidimensional) Gaussian:

,

where to stand for the position in the space of active coordinates, and and are the height and width of the Gaussian, respectively. Note that both positive and negative values are allowed for the parameter .

For example, if two active coordinates with STATUS=5 are defined in the ICONST file:

R 1 5 5
R 1 6 5 

then each line in the PENALTYPOT must contain four items. The bias potential is defined in the following lines:

1.6 0.8 1.0 0.2
1.6 1.0 1.0 0.2
1.6 1.2 1.0 0.2
1.6 1.4 1.0 0.2
1.6 1.6 1.0 0.2
1.6 1.8 1.0 0.2
1.6 2.0 1.0 0.2

Related tags and articles

MDALGO, ICONST