[Wien] Problem with minimization

Reza Mahani kh.mahani at gmail.com
Wed Sep 12 16:20:26 CEST 2012


Hi Prof Blaha and wien2k users,

   I am trying to optimize internal coordinates of a small semiconductor
surface using mini (the 1st methods suggested in the manual):


   min_lapw -j "runsp_lapw -I -fc 2.0 -p -it -NI -i 100"

I have followed the steps mentioned in the manual but after 5 geometry it
sopped with the message in case.dayfile that force is small. However, when
we checked the force on different atoms (shown below), it shows larger
forces. So we are not sure why the job stopped.


Top of case.inM:

PORT 2.00 0.35      # PORT/NEWT;  tolf, Initial Trust Radius
1.0 1.0 1.0 1.0   #Atom    1 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    2 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    3 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    4 Generated by pairhess


Forces:

:FGL001:   1.ATOM     -24.500500000     0.000000000   -49.905000000 total
forces
:FGL002:   2.ATOM      -4.855000000     0.000000000    -2.205500000 total
forces
:FGL003:   3.ATOM     -15.616000000     0.000000000    -4.224000000 total
forces
:FGL004:   4.ATOM       6.779500000     0.000000000    -6.616500000 total
forces
:FGL005:   5.ATOM       8.279000000     0.000000000    -2.916500000 total
forces
:FGL006:   6.ATOM      33.619000000     0.000000000   -40.479000000 total
forces
:FGL007:   7.ATOM      10.189500000     0.000000000    -1.776500000 total
forces
:FGL008:   8.ATOM      -0.703000000     0.000000000    -2.872000000 total
forces
:FGL009:   9.ATOM       2.128500000     0.000000000    57.396500000 total
forces
:FGL010:  10.ATOM      23.293500000     0.000000000   -34.333000000 total
forces
:FGL011:  11.ATOM      -3.002500000     0.000000000    -1.038000000 total
forces
:FGL012:  12.ATOM     -22.008500000     0.000000000    13.880500000 total
forces
:FGL013:  13.ATOM     -29.674000000     0.000000000    45.145000000 total
forces
:FGL014:  14.ATOM      16.142000000     0.000000000    23.337000000 total
forces


I would appreciate if anyone can tell me if I am doing anything wrong in
setting up the calculations.


Regards
Reza
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20120912/47d9e372/attachment.htm>


More information about the Wien mailing list