In the scf file you find the volume:

       LATTICE                      = B   
:POT  : POTENTIAL OPTION EX_PBE EC_PBE VX_PBE VC_PBE
:LAT  : LATTICE CONSTANTS= 11.31946 11.79152 16.06778    1.571    1.571    1.571
:VOL  : UNIT CELL VOLUME =    1072.31280

and 11.31946 * 11.79152 * 16.06778 / 2 = 1072.312797243878022456

so what is the problem, beside the rounding  ?

By the way, what you are doing is "searching for a local minimum" but not an 
"optimization"!

You have to find the absolute (!) energy minimum in the a, b, c space
In yourprocedure it might depend on which "minimum energy search" you use to 
start, and on the series of seraches.
You can not expect a priori that only one minimum exists in the range
a +- Delta a, b +- Delta b, c +- Delta c 
and yo might end up in a "wrong" local minimum.

Further, your structure has free position parameters, 
therefore, you have to optimize the positions TOGETHER with the lattice 
parameters,
obviously you did not do that, the positions are equal in both files ?

For completeness:
The "Volume optimization" alone has only a physical meaning for cubic systems 
(or spherical as most of the equations of states are for the pressure-volume 
relation inside earth, 
I wonder how they would look like if the earth was a disk, ... or maybe it is 
;-).

Ciao
Gerhard

DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:
"I think the problem, to be quite honest with you,
is that you have never actually known what the question is."

====================================
Dr. Gerhard H. Fecher
Institut of Inorganic and Analytical Chemistry
Johannes Gutenberg - University
55099 Mainz
and
Max Planck Institute for Chemical Physics of Solids
01187 Dresden
________________________________________
Von: Wien [wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von shamik 
chakrabarti [shamik15041...@gmail.com]
Gesendet: Samstag, 13. Juni 2020 07:14
An: A Mailing list for WIEN2k users
Betreff: Re: [Wien] Volume has changed during running optimization of c/a with 
constant volume & b/a

Dear All,
              I am sending the LiNiNbO4_ca_initial.struct file which is 
obtained after volume optimization with a:b:c =constant. Also I am sending  
LiNiNbO4_ba_initial.struct which is obtained after c/a optimization with 
constant volume and b/a. Both the structures are showing same unit cell volume 
as ~ 2144 bohr^3 as I have obtained in my calculation where as the printed 
volume in the w2web page obtained after c/a optimization is  2058.7690 bohr^3 
which is different than the actual optimized volume we have obtained. Please 
comment on this.

with regards,


On Sat, 13 Jun 2020 at 02:03, Tomas Kana 
<k...@seznam.cz<mailto:k...@seznam.cz>> wrote:

Zpět
[https://i.im.cz/private-avatar/6be7aba869b05ecd061bd77ff52f23ec0d7890f91360a35dc19cfad48b41ad7080ea62f5311a154a/96]Dnes
 22:15

Tomas Kana
Komu:
wien@zeus.theochem.tuwien.ac.at<mailto:wien@zeus.theochem.tuwien.ac.at>
Re: [Wien] Volume has changed during running optimization of c/a with constant 
volume & b/a
Dear Shamik, If I understand you correctly then I suggest to visualise better 
your energy-volume curve. Try this (it works for me): 1. Keep c/a and b/a 
experimental ones and run x optimize option 1 and make an energy-volume curve 
by modifying he option save_lapw -d ../$i  $i in the file optimize.job that is 
produced by x optimize. 2. For each separated directory run_lapw and collect 
all the scf files to one one directory and visualise the resulting 
energy-volume curve by  eplot -t vol 3. In each separated directory optimize 
c/a  at constant volume and b/a. Again collect all scf files and visualise both 
"experimental" and "relaxed" energy volume curve in one plot. You will probably 
note and understand the difference between the two mimimums. 4. Do the same 
optimization for b/a at constant volume and c/a and draw the three 
energy-volume curves together. Their minima could again differ. Repeat 3 and 4 
until the resulting e-v curve does not substantially differ from the preceding 
one. Alternatively you may use the package OrthoOpt from wien2k goodies.
_______________________________________________
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at<mailto:Wien@zeus.theochem.tuwien.ac.at>
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html


--
Dr. Shamik Chakrabarti
Research Fellow
Department of Physics
Indian Institute of Technology Patna
Bihta-801103
Patna
Bihar, India
_______________________________________________
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Reply via email to