Quick question on selective dynamics

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


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
richard_kang
Newbie
Newbie
Posts: 6
Joined: Sun Feb 21, 2021 2:01 am

Quick question on selective dynamics

#1 Post by richard_kang » Mon Feb 14, 2022 1:28 am

Hello,

I think I found a related question which was uploaded long time ago relevant to this but needed a further clarification.
I am using selective dynamics tag for ionic relaxation, but I am specifically fixing only the Z-coordinate of a certain atom. (i.e. T T F)
I'm doing a quite coarse optimization, so using EDIFFG=-0.05 for the convergence criterion.

I was wondering about how VASP treats the force on this atom.
So for example for a finished electronic step which gives this on the above atom:

25166 POSITION TOTAL-FORCE (eV/Angst)
25167 -----------------------------------------------------------------------------------
25168 -0.00039 0.00029 0.00000 0.002144 0.001395 0.036577


Because I specified, TTF, would projection of the force to the xy-plane still be considered for convergence criterion?
(i.e. evaluate if 0.002144^2 + 0.001395^2 = 0.0367 is smaller than 0.05)
Or would the force be just ignored for convergence criterion whatsoever per the F?

Thank you very much!
Richard

ferenc_karsai
Global Moderator
Global Moderator
Posts: 460
Joined: Mon Nov 04, 2019 12:44 pm

Re: Quick question on selective dynamics

#2 Post by ferenc_karsai » Tue Feb 15, 2022 9:31 am

In the code all atoms contribute to the force threshold, but the forces on the atoms with .FALSE. are set to zero. This is way they effectively don't contribute to the threshold.

richard_kang
Newbie
Newbie
Posts: 6
Joined: Sun Feb 21, 2021 2:01 am

Re: Quick question on selective dynamics

#3 Post by richard_kang » Wed Feb 16, 2022 1:32 am

Ah okay, thanks!
But since the atom is labeled TTF, it will still be able to move around in the direction of first two lattice vectors right?
Thanks for the reply.

ferenc_karsai
Global Moderator
Global Moderator
Posts: 460
Joined: Mon Nov 04, 2019 12:44 pm

Re: Quick question on selective dynamics

#4 Post by ferenc_karsai » Wed Feb 16, 2022 7:44 am

Yes of course, the forces will fully act on the first two dimensions.

Post Reply