Setting ICONST file with Centre of Mass as Collective Variable

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
pranav_roy
Newbie
Newbie
Posts: 1
Joined: Sun Apr 23, 2023 10:26 pm

Setting ICONST file with Centre of Mass as Collective Variable

#1 Post by pranav_roy » Sun Dec 22, 2024 3:48 am

Dear VASP Team,

I was trying to set up an ICONST file for a calculation where I'm trying to set up my collective variable as the distance between a C atom, and the Centre of Mass of 3 atoms (say H atoms). I know that we define primitive coordinate using R flag, and complex coordinates using S flag, however rest of the mentioned flags are for cartesian coordinates and bond centers as such. Is there any way to implement this in ICONST file for a metadynamics simulation in VASP?

Thank you for the help!


christopher_sheldon1
Global Moderator
Global Moderator
Posts: 74
Joined: Mon Mar 25, 2024 1:36 pm

Re: Setting ICONST file with Centre of Mass as Collective Variable

#2 Post by christopher_sheldon1 » Tue Jan 07, 2025 9:46 am

Dear Pranav,

Thank you for your question. There is no specific way to define the centre of mass using ICONST.

However, ICONST can be used to define more than two coordinates, so you could define three different primitive coordinates,e.g. the distance between the C atoms and each H atom, then monitor these:

Code: Select all

R 3 4 0
R 3 5 0
R 3 6 0
S 1 0 0 7
S 0 1 0 7
S 0 0 1 7

Using your own script, you could calculate the distance to the centre of mass of these atoms. E.g. in my output it gives:

Code: Select all

  Monit_coord
   mc> S          1.07090
   mc> S          1.12301
   mc> S          1.05014

Does this answer your question?

Best wishes,

Chris


Post Reply