[Wien] (no subject)

Santhy Jaiker k_santhy79 at yahoo.co.in
Mon Jul 27 06:47:33 CEST 2009


Dear Marks,
 
thanks for immediate reply.    i listed below the FGL001 in case.scf_mini and case.scf.
There is a jumb in FGL at last iteration in case.scf.  That makes me doubt about my calculation. 
FGL001 case.scf
:FGL001:   1.ATOM               -11.099558998   -22.199117923     0.000000000 partial forces
:FGL001:   1.ATOM               -11.162490178   -22.324980282     0.000000000 partial forces
:FGL001:   1.ATOM               -11.346664915   -22.693329755     0.000000000 partial forces
:FGL001:   1.ATOM               -12.054496349   -24.108992620     0.000000000 partial forces
:FGL001:   1.ATOM               -13.053889672   -26.107779258     0.000000000 partial forces
:FGL001:   1.ATOM               -12.542934680   -25.085869279     0.000000000 partial forces
:FGL001:   1.ATOM               -11.962697656   -23.925395234     0.000000000 partial forces
:FGL001:   1.ATOM               -12.057383101   -24.114766123     0.000000000 partial forces
:FGL001:   1.ATOM                -0.110273902    -0.220547804     0.000000000 total forces
FGL001 case.scf_mini
:FGL001:   1.ATOM                 2.994138516     5.988277012     0.000000000 t
:FGL001:   1.ATOM                -0.136254664    -0.272509328     0.000000000 t
:FGL001:   1.ATOM                -0.110273902    -0.220547804     0.000000000 t


 
 
 
will you explain, "What converges faster, a heavily damped oscillator or a
slightly underdamped one?" this!

 
Regards
Santhy Jaiker
 
Research Scholar 
Materials Thermodynamics Lab
Department of Metallurgical and Materials Engineering
Indian Institute of Technology Madras
Chennai - 600 036
India
Tel: 91-044-22575763
**********************************************
Whatever the mind can conceive and believe, the mind can achieve. 
Dr. Napoleon Hill
**********************************************


--- On Mon, 27/7/09, Laurence Marks <L-marks at northwestern.edu> wrote:


From: Laurence Marks <L-marks at northwestern.edu>
Subject: Re: [Wien] (no subject)
To: "A Mailing list for WIEN2k users" <wien at zeus.theochem.tuwien.ac.at>
Date: Monday, 27 July, 2009, 8:55 AM


2009/7/26 Santhy Jaiker <k_santhy79 at yahoo.co.in>:
> Dear all,
>
> i have some few doubts in running minimization.
> 1.  i have to come to conclusion about the force from FORxxx or FGLxxx?
> 2. FORxxx or FGLxxx from case.scf_mini or case.scf?

FGL -- they are normally the same

> 3. I got the error "STOP in MINI, FORCES small".  I think, it doesn't mean
> error, it is converged.  am i right?

Read the output file, case.outputM. If you have setup the problem
badly it may not have converged.

> 4. In faq, mentioned that DELTA was too large (force changed sign !!!,), it
> could even diverge.  in my case, FOR001 changed from 5.186 to -0..236.
> DELTA is small, but force
> is changed in sign. how to consider this?

The faq mentioned DELTA only for NEWT, it is ignored for PORT.
There is nothing wrong with forces changing sign, it is to be
expected. What converges faster, a heavily damped oscillator or a
slightly underdamped one?

> 5. There are forces in atom 1 and 3.  After minimization, the position
> coordiante are changed.
>
> atom 1: 0.1697      0.3394        0.25
>
Of course!
>
>
>  I have to check FORxxx a FGLxxx from laves.scf_mini
>
>
>  After running minimization,
> i have to check ENE, GRAD, MIN, WARN, DD from case.outputM and ENE, FORxxx
> from
> case.scf_mini.
>
> In outputM file, ENE is decrease, GRAD is also decreasing and there is no
> warning.
>
> Energy in case.scf_mini
> ENE  : ********** TOTAL ENERGY IN Ry =       -27195.281278
> :ENE  : ********** TOTAL ENERGY IN Ry =       -27195.281764
> :ENE  : ********** TOTAL ENERGY IN Ry =       -27195.281794
> Force in atom 1 obtained from  case.scf_mini
> :FOR001:   1.ATOM          5.186          0.000          0.000
> 5.186 t
> :FOR001:   1.ATOM          0.236          0.000          0.000
> -0.236 t
> :FOR001:   1.ATOM          0.191          0.000          0.000
> -0.191 t
> Force in atom 3 obatined from case.scf_mini
> :FOR003:   3.ATOM          3.408          0.000          0.000
> 3.408 t
> :FOR003:   3.ATOM          2.296          0.000          0.000
> 2.296 t
> :FOR003:   3.ATOM          1.253          0.000          0.000
> 1.253 t
>
> FGL001  from case.scf_mini
> :FGL001:   1.ATOM                 2.994138516     5.988277012
> 0.000000000 t
> :FGL001:   1.ATOM                -0.136254664    -0.272509328
> 0.000000000 t
> :FGL001:   1.ATOM                -0.110273902    -0.220547804
> 0.000000000 t
>
> FGL003  from case.scf_mini
> :FGL003:   3.ATOM                 0.000000000     0.000000000
> 3.408000000 t
> :FGL003:   3.ATOM                 0.000000000     0.000000000
> 2.296000000 t
> :FGL003:   3.ATOM                 0.000000000     0.000000000
> 1.253000000 t
> Santhy Jaiker
>

The minimum energy position has no forces. Think about it.

> Research Scholar
> Materials Thermodynamics Lab
> Department of Metallurgical and Materials Engineering
> Indian Institute of Technology Madras
> Chennai - 600 036
> India
> Tel: 91-044-22575763
> **********************************************
> Whatever the mind can conceive and believe, the mind can achieve.
>
> Dr. Napoleon Hill
>
> **********************************************
>
>
> --- On Mon, 27/7/09, Bahram Abedi <abedi527 at gmail.com> wrote:
>
> From: Bahram Abedi <abedi527 at gmail.com>
> Subject: [Wien] external magnetic field
> To: wien at zeus.theochem.tuwien.ac..at
> Date: Monday, 27 July, 2009, 2:38 AM
>
> Dear wien users,
> I'm trying to check the effect of an externally applied magnetic field
> on some simple metals (following the procedure given in section 7.2 of
> the usergiude). The problem is that I don’t see any changes in the
> calculated properties (dos, magnetic moment, etc.) under the field
> application.
> Can you please tell me if this is okey? If not, what is wrong?
>
> Regards,
> Bahram Abedi.
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>
> ________________________________
> See the Web's breaking stories, chosen by people like you. Check out Yahoo!
> Buzz.
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>
>



-- 
Laurence Marks
Department of Materials Science and Engineering
MSE Rm 2036 Cook Hall
2220 N Campus Drive
Northwestern University
Evanston, IL 60208, USA
Tel: (847) 491-3996 Fax: (847) 491-7820
email: L-marks at northwestern dot edu
Web: www.numis.northwestern.edu
Chair, Commission on Electron Crystallography of IUCR
www.numis.northwestern.edu/
Electron crystallography is the branch of science that uses electron
scattering to study the structure of matter.
_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien



      Looking for local information? Find it on Yahoo! Local http://in.local.yahoo.com/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20090726/dcc2ce1e/attachment.htm>


More information about the Wien mailing list