My Community

Hi Everyone,
Is there anyone know how to setup the metaDynamics to simulate the alloy transition? It is very appreciated if you could help this.
I found the metaDynamics in the website (http://cms.mpi.univie.ac.at/vasp/vasp/M ... ics_I.html) as following:

For the brief overview of the method, see Sec. 6.62.1 and references therein. The computational setup for metadynamics should contain the following items:
1. set the standard MD-related flags: IBRION=0, TEBEG, POTIM, NSW
2. set MDALGO to 11 or 21 for metadynamics with an Andersen, or Nose-Hoover thermostat, respectively (check description of ANDERSEN_PROB and SMASS for the thermostat-specific setting)
3. define the parameters HILLS_H, HILLS_W, HILLS_BIN, see Sec. 6.62.3.
4. define collective variables in the file ICONST (see Sec. 6.62.4), the input parameter STATUS for collective variables must be set to 5
5. if needed, define bias potential in file PENALTYPOT, see Sec. 6.62.4.

I used the MDALGo = 11, but don't know how to setup the ICONST. Can any body help me? Email: tanliuxi@gmail.com.
Thanks very muck for view!

ICONST defines constraints. These depend on your system and on the collective variable. A set of different collective variables can be used to sample the same reaction coordinate.

Hi Head admin,
It is very appreciated. Could you show me some sample of the ICONST file?

In this example
http://cms.mpi.univie.ac.at/vasp/vasp/P ... penaltypot
ICONST file defines two interatomic distances (between atoms 1 - 5 and 1 - 6) in order to make constraints during the metadynamics run (the distance 1 -5 is not allowed to become larger than 1.6 A). This is achieved using the bias potential defined in PENALTYPOT file (seven lines stand for seven Gaussians located as defined by first two columns, each Gaussian is 1.0 eV high with the halfwidth of 0.2 eV).

Thanks.

In metadynamics calculation with vasp, however, I found there is no clear description how to mention different intermediates along the reaction profile. Should I need to perform separate run for each intermediate or they can be put together in the same run? If you please clarify this I would greatly appreciate. Thank you for your help!

How large system to take and how long simulation to run depends on your computer resources. The ideal case is the simulation of the complete story reactant-intermediate-product.

Is there any example how to define reactant-intermediate-product (in ideal case) in the same input of metadynamics? If you describe it little bit in more details it would be of great help. Thank you very much for your reply!

My system contains a total of 102 atoms, however, 64 of those are fixed. (Simulation Box: ~10 A x 10 A x 15 A)

Have a look at a simple reaction described in JP-CM 20(2007)064211, eventually contact the author of the paper.

Thank you very much for your suggestion and help!