[Wien] minimization of positions

Hem Chandra Kandpal kandpal at uni-mainz.de
Sat Dec 24 12:43:25 CET 2005


Dear Wien Users, 

I am trying to optimize the positions of a ternary system.
Below I am giving full detail about that system
system is 
MnCoSn ( a = 4.258 Ang, c = 5.36) (Space group 186)
Mn (0 0 1/4)
Co (1/3 2/3 0.44=z1)
Ge (1/3 2/3 0.06=z2)

To optimize the  z parameters of Co (z1) and Ge (z2) atoms.
I am using PORT option as given in the manual.

Let me tell you one thing, my system positions are far away from the
solution, because this system is only reported as Ni2In  structure type
(space group 194) and I would like to minimize it in LiGaGe structure 
type (space group 186). In space group 186 each atom will have their 
own positions as given above. These positions I took from other ternary system. 

I follow this procedure 
1. runsp_lapw -ec 0.0001 -i 40 -r 15
2. runsp_lapw -cc 0.001 -i 40 -r 15
3. runsp_lapw -fc 1 -i 40 -r 15

then I Have edited case.inM file like this

PORT 2.0           # PORT/NEWT;  tolf (geom.opt. stops when F<;tolf) 
0.0 0.0 1.0 3.0    # Mn
0.0 0.0 1.0 1.0    # Co #PORT: delta1,2,3=precond., 4:Bondorder 
0.0 0.0 1.0 1.0    # Sn

when I start with ---> min -NI -i 4 -s 1 -j 'runsp_lapw -I -i 40 -fc 1.0 '  

It is running and my forces are very far from force criterion and at the 
end in my dayfile it writes line like this -->

##########################dayfile
:ENERGY convergence:  0 0 .0000525000000000
:CHARGE convergence:  0 0.0000 .0002335

>   stop
>   mini        (01:07:58) 0.070u 0.015s 0:00.08 100.0% 0+0k 0+0io 0pf+0w
 :WARNING, CURVATURE CONDITION FAILED
 :WARNING, S.Y WAS -0.304898817526764
 :WARNING: Step size reduced due to overlapping spheres -- check RMT
 :WARNING: Step size reduced due to overlapping spheres -- check RMT
############################################################

my forces are looking like --->

###########################################
:ENE  : ********** TOTAL ENERGY IN Ry =       -34925.629006

       TOTAL FORCE IN mRy/a.u. = |F|     Fx             Fy             Fz
:FOR001:   1.ATOM        114.688     0.000     0.000   114.688
:FOR002:   2.ATOM        992.576     0.000     0.000   -992.576
:FOR003:   3.ATOM       1776.551     0.000     0.000   1776.551

#############################################

Can anyone please explain me, what is the process to do such kind of 
minimization, am I doing the right way? and what to do if I see 
with iterations, forces are increasing on each atom and vice versa.

What is the criterion for choosing RMT values for such complicated 
(I mean without knowing exact positions) systems.



############################

In the begining, I thought, I am doing some thing wrong with 
my calculations. then I tried to check these steps (as given above for 
minimization) on a well reported system 

example LiGaGe  (a = 4.175, c = 6.783 Ang) 
	Li (0 0 1/4)
        Ga (1/3 2/3 0.06)
        Ge (1/3 273 0.44)

and it coverges without any error, in my dayfile 

#############################################
>   stop
>   mini        (16:09:10) 0.029u 0.009s 0:00.04 50.0%  0+0k 0+0io 0pf+0w

>   stop forces small, minimization stopped

>   relaxation steps summarized in LiGaGe_opt.scf_mini

>   relaxed structure and results in LiGaGe_opt.struct, scf, clmsum
##############################################


and in my third optimized scf file :

############################################# 
:ENE  : ********** TOTAL ENERGY IN Ry =       -16202.425763

       TOTAL FORCE IN mRy/a.u. = |F|     Fx             Fy             Fz
:FOR001:   1.ATOM          1.242          0.000          0.000    1.242
:FOR002:   2.ATOM          2.019          0.000          0.000    -2.019
:FOR003:   3.ATOM          0.428          0.000          0.000     0.428

############################################


If it is working for such stabilized system it should work for 
my unrelaxed system. I think somewhere I missing some points.

What would be the good option for such system either PORT 
or NEWT. Can I also edit DELTA values (range???) in my 
case.inM file with PORT option like Prof. Blaha expalined
for NEWT in FAQ of Wien. 


I will be grateful for your kind reply and help

Thanking you in advance 
hem kandpal (mainz)


More information about the Wien mailing list