[Pw_forum] calculating equilibrium lattice constant

Lorenzo Paulatto paulatto at sissa.it
Mon Oct 26 15:28:29 CET 2009


In data 26 ottobre 2009 alle ore 15:00:25, Hande Ustunel  
<hande at newton.physics.metu.edu.tr> ha scritto:
> If it is a shell script that you are running you might try something like

I don't want to be picky, but there was a missing "$", this is a correct  
script:

    for file in out.*
    do
      E=`grep ! $file | tail -1 | awk '{print $5}'`
      echo $E >> energies-file
    done

If you only want to find the energy minimum you can just do a vc-relax  
calculation, which will do all the work for you. It will not give you the  
energy/volume curve though, and it may relax more than you expect (e.g.  
moving some parameter that you wanted to keep fixed).

regards

-- 
Lorenzo Paulatto
SISSA  &  DEMOCRITOS (Trieste)
phone: +39 040 3787 511
skype: paulatz
www:   http://people.sissa.it/~paulatto/

     *** save italian brains ***
  http://saveitalianbrains.wordpress.com/


More information about the Pw_forum mailing list