<div dir="ltr">Hi Everyone,<div><br></div><div>I am learning the metadynamics by the website: <a href="https://www.cp2k.org/exercises:2015_cecam_tutorial:mtd1">https://www.cp2k.org/exercises:2015_cecam_tutorial:mtd1</a></div><div><br></div><div>However, there is a question on the <a href="https://manual.cp2k.org/trunk/CP2K_INPUT.html" style="font-size: medium; font-family: monospace; text-transform: uppercase;">CP2K_INPUT</a><span style="font-size: medium; color: rgb(0, 0, 0); font-family: monospace; text-transform: uppercase;">/</span><a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL.html" style="font-size: medium; font-family: monospace; text-transform: uppercase;">FORCE_EVAL</a><span style="font-size: medium; color: rgb(0, 0, 0); font-family: monospace; text-transform: uppercase;">/</span><a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/EXTERNAL_POTENTIAL.html" style="font-size: medium; font-family: monospace; text-transform: uppercase;">EXTERNAL_POTENTIAL</a>. In the first task of tutorial, an example on the dynamics of two HNO3 molecules over a graphene sheet are considered. If a constraint is not added, it is observed that the two molecules very fast and expolre a large spece of configurations, being often far from each other and far from graphene. In order to restrict the exploration to region, it is necessary to limit the movement of the two molecules. To this purpose, an external potential defining a spherical potential centered on the center of the system coordinated is added. The potential is given as follows:</div><div><pre class="code" style="padding: 0.7em 1em; font-family: Consolas, "Andale Mono WT", "Andale Mono", "Bitstream Vera Sans Mono", "Nimbus Mono L", Monaco, "Courier New", monospace; font-size: 14px; direction: ltr; background-color: rgb(251, 250, 249); color: rgb(51, 51, 51); box-shadow: rgb(204, 204, 204) 0px 0px 0.5em inset; border-radius: 2px; overflow: auto; overflow-wrap: normal; border-width: 1px; border-style: solid; border-color: rgb(204, 204, 204);"> &EXTERNAL_POTENTIAL
     ATOMS_LIST   61..70
     FUNCTION   0.000000000001*(Z^2)^4
  &END

  &EXTERNAL_POTENTIAL
     ATOMS_LIST   61..70
     FUNCTION   0.0000000000001*(X^2)^4
  &END

  &EXTERNAL_POTENTIAL
     ATOMS_LIST   61..70
     FUNCTION   0.0000000000001*(Y^2)^4
  &END</pre></div><div><br></div><div>Question: How does work on such external_potential? i.e., How does such external_potential constaint the two HNO3 molecules to the vicinity of graphene?</div><div><br></div><div>Ning</div></div>