MetaDynamics MD for transition
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 3
- Joined: Mon May 20, 2013 9:34 pm
MetaDynamics MD for transition
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!
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!
Last edited by tanlx on Mon May 20, 2013 10:54 pm, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
MetaDynamics MD for transition
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.
Last edited by admin on Tue May 21, 2013 1:08 pm, edited 1 time in total.
-
- Newbie
- Posts: 3
- Joined: Mon May 20, 2013 9:34 pm
MetaDynamics MD for transition
Hi Head admin,
It is very appreciated. Could you show me some sample of the ICONST file?
It is very appreciated. Could you show me some sample of the ICONST file?
Last edited by tanlx on Thu May 23, 2013 6:33 pm, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
MetaDynamics MD for transition
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).
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).
Last edited by admin on Fri May 24, 2013 11:20 am, edited 1 time in total.
-
- Newbie
- Posts: 3
- Joined: Mon May 20, 2013 9:34 pm
MetaDynamics MD for transition
Thanks.
Last edited by tanlx on Mon Jun 03, 2013 7:14 pm, edited 1 time in total.
MetaDynamics MD for transition
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!
Last edited by kmrana on Fri Jul 11, 2014 2:09 am, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
MetaDynamics MD for transition
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.
Last edited by admin on Fri Jul 11, 2014 9:06 am, edited 1 time in total.
MetaDynamics MD for transition
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)
My system contains a total of 102 atoms, however, 64 of those are fixed. (Simulation Box: ~10 A x 10 A x 15 A)
Last edited by kmrana on Sat Jul 12, 2014 12:27 am, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
MetaDynamics MD for transition
Have a look at a simple reaction described in JP-CM 20(2007)064211, eventually contact the author of the paper.
Last edited by admin on Tue Jul 15, 2014 8:25 am, edited 1 time in total.
MetaDynamics MD for transition
Thank you very much for your suggestion and help!
Last edited by kmrana on Sat Jul 19, 2014 6:06 am, edited 1 time in total.