[Wien] PORT, atomic position relaxation

Stefaan Cottenier Stefaan.Cottenier at fys.kuleuven.ac.be
Wed Aug 4 23:01:58 CEST 2004


> > Does that mean i wont be able to relax the "y" coordinate
> > in any  of these cases.
> 
> Indeed.

Forgot to add: if you really want to relax also this y=1/4 coordinate, 
then do the following:

* Change it to 0.250001 (consistently, if you have a corresponding 3/4 
position for one of the atoms at that position, then change it to 
0.749999).

* Repeat init_lapw. sgroup and nn will reduce the symmetry of your 
cell, and you will end up with more k-points.

* Converge one scf-cycle, and see what is the force for that y-
coordinate (you will have a force now, because in the lower space 
group this position will not be fixed by symmetry any more). But very 
often, you will see that the corresponding force is zero (within the 
numerical noise).

Stefaan




More information about the Wien mailing list