[Wien] Regarding lattice constant missing in volume optimization curve data and also in output file of tetragonal lattice

Bhanu Prakash bprakashjune2011 at gmail.com
Wed Dec 28 10:19:12 CET 2022


D*ear Prof. Blaha and Wien2k users,*

I am currently working with Wien2k_21.1 version. I would like to get
some suggestions for doing optimization of lattice parameters for a
tetragonal alloy .The experimental lattice parameter for the alloy is
*a=b=**5**.**426** A ,c=**10.816** A**ng**.*
I have constructed the structure file with the space group I-4 (space
group number:82) and the above mentioned experimental lattice
parameters.
After doing the initialization with 1000 k points energy minimization
has been done by varying the volume -10%,-5%,0%,5% and 10 % with
constant c/a optimization.I obtained a fairly good E Vs c/a curve.
But I am facing some problems.
Now my problems are

1) I did not get the optimized lattice parameters (l*attice constant
after the optimization process)* with varying the c/a ratio with
constant volume.

2) In my case.outputeos file, I did not get a lattice constant value.
I also used 2Doptimizaton process but result is same.
*Here**Also missing the lattice constant value in fourth line of
case.outputeos file.(blank space ##)*
Equation of state: EOS2 (PRB52,8064)        info           2
 a,b,c,d   2611383.290602 -38111329.012348 -59534474.9093151774906515.364957
 V0,B(GPa),BP,E0      1112.1797  12436214.1735         3.3770
 #####################################################################
 Equation of state: Murnaghan                info           7
 E=E0+[B*V/BP*(1/(BP-1)*(V0/V)**BP +1)-B*V0/(BP-1)]/14703.6
 Pressure=B/BP*((V0/V)**BP -1)
 V0,B(GPa),BP,E0      1110.7646       100.0000         5.0000     -25764.476906
         vol       energy         de(EOS2)      de(Murnaghan)  Pressure(GPa)
    1116.6428   -25764.478527     0.001765     0.001726         -0.521
    1116.6429   -25764.478351     0.001993     0.001550         -0.521
    1116.6427   -25764.478004     0.001007     0.001203         -0.521
    1116.6427   -25764.477467     0.000502     0.000666         -0.521
    1116.6428   -25764.475907    -0.000990    -0.000894         -0.521
    1116.6429   -25764.473677    -0.002749    -0.003124         -0.521
    1116.6429   -25764.478544     0.001950     0.001743         -0.521
    1116.6427   -25764.478367     0.001302     0.001566         -0.521
    1116.6428   -25764.478017     0.001390     0.001216         -0.521
    1116.6427   -25764.477497     0.000533     0.000696         -0.521
    1116.6427   -25764.476816    -0.000148     0.000015         -0.521
    1116.6427   -25764.474912    -0.002287    -0.001889         -0.521
    1116.6429   -25764.472327    -0.004267    -0.004474         -0.521
                  Sigma:          0.001923     0.001939

Please give your valuable suggestions. How to get lattice constant
from optimization process.

Thanks & Regards
Bhanu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20221228/7cb6644c/attachment.htm>


More information about the Wien mailing list