The optimization of 3 lattice parameters in wien2k is quite some work
and requires large computational effort.
To use just 10 structures is the minimal requirement to determine 10
fit-parameters. Usually, such a solution will suffer from the smallest
numerical noise and not give accurate results.
In any case, your procedure gave you its solution as output (but as
mentioned above, it might suffer from noise and is probably very
inaccurate):
At the bottom you found:
> Parabolic equation of state: info 8
>
> E = x1 + x2(a-x3)^2
>
> + x4(b-x5)^2 + x6(a-x3)(b-x5)
>
> + x7(c-x8)^2 + x9(a-x3)(c-x8) + x10(b-x5)(c-x8)
>
> Fitparameter are
>
> -195722.961991 0.024128 15.665728 0.031781
>
> 15.308083 0.060289 0.069122 11.327643
>
> -5.455925 -5.218681
So this gives your a,b,c as 15.666, 15.308, 11.328
On the other hand it seems that you got the lowest total energies when
a is small and b is large !!! So I doubt that this is an accurate
solution. You need more points (3x3x3 or, of course, even more ...)
On 10/14/2014 09:35 AM, Shafqat Hussain Shah wrote:
Dear All
I have trouble in extracting the optimized lattice parameters (a,b c)
from my optimized structures using option 6 plus MSR1a.
My system of interest has an orthorhombic structure. I want to
simultaneously optimize its volume and internal atomic positions. I have
used option 6 (with number of structures 10 at 1% change) of the
optimizer plus MSR1a with default parameters in case.inM. All the 10
structures are geometry optimized well and forces on atoms are less than
the tolerances. Then I used the command
parabolfit_lapw -t 3 -f BIFOfm -scf *abc*.scf
to fit the data as shown in the end.
My question ; what is the best way to extract the optimized lattice
parameters for my system from this procedure????
I have tried following things but to no avail.
1> I tried to plot the data (a,b,c, E) with a scatter plot in Origin
where E was color coded. The number of data points was not sufficient
(10) and it did not give the exact values of the optimized lattice
parameters.
2> I tried to find the minima of the fitted function with first order
and second order derivatives (Hessian) but it did not give any information.
3> I looked up in the mailing list but could not find the answer to
my problem.
Any suggestion to extract the optimized lattice parameters from the data
will be highly appreciated.
Many thanks in advance.
Dr. Shafqat Hussain Shah
PS: I am using parallel WIEN2k_13 on a multicore machine.
shs28@apollo1:~/WIEN2k/BIFO/BIFOfm> parabolfit_lapw -t 3 -f BIFOfm -scf
*abc*.scf
BIFOfm.ene and BIFOfm.latparam generated
Enter dimension of fit (number of variable lattice parameters, 1-6):
3 fitcase 10 parameter
lowest data point: -195723.042487450 15.0632000000000
15.7881200000000 11.2255400000000
10 -195723.042487450 -195722.955530700
-195722.964214750 -195722.971284580 -195722.954407870
-195722.971888480 -195722.957634800 -195722.968628590
-195722.944729150 -195723.042426910
I INITIAL X(I) D(I)
1 -0.195723D+06 0.100D+01
2 0.100000D+00 0.600D+00
3 0.150632D+02 0.600D+00
4 0.100000D+00 0.600D+00
5 0.157881D+02 0.600D+00
6 0.000000D+00 0.600D+00
7 0.100000D+00 0.600D+00
8 0.112255D+02 0.600D+00
9 0.000000D+00 0.600D+00
10 0.000000D+00 0.600D+00
IT NF F RELDF PRELDF RELDX MODEL STPPAR D*STEP NPRELDF
0 1 0.237D-01
1 3 0.217D-01 0.85D-01 0.83D+00 0.7D-06 G 0.7D-02 0.4D+00 0.10D+01
2 4 0.621D-02 0.71D+00 0.84D+00 0.3D-06 G 0.4D-02 0.2D+00 0.10D+01
3 6 0.333D-02 0.46D+00 0.47D+00 0.7D-07 G 0.1D+00 0.4D-01 0.82D+00
4 8 0.310D-02 0.69D-01 0.14D+00 0.2D-06 G 0.4D-01 0.1D+00 0.67D+00
5 9 0.257D-02 0.17D+00 0.19D+00 0.2D-06 G 0.3D-01 0.9D-01 0.55D+00
6 10 0.234D-02 0.87D-01 0.96D-01 0.2D-06 G 0.2D-01 0.8D-01 0.41D+00
7 11 0.229D-02 0.22D-01 0.19D-01 0.2D-07 G 0.1D-01 0.1D-01 0.35D+00
8 13 0.226D-02 0.16D-01 0.16D-01 0.6D-07 S 0.3D-01 0.3D-01 0.14D+00
9 14 0.225D-02 0.15D-02 0.12D-02 0.8D-08 G 0.9D-01 0.5D-02 0.10D+01
10 16 0.223D-02 0.12D-01 0.12D-01 0.5D-07 G 0.3D-01 0.3D-01 0.23D+00
11 17 0.220D-02 0.93D-02 0.60D-02 0.4D-07 G 0.1D-01 0.2D-01 0.10D+01
12 18 0.219D-02 0.64D-02 0.41D-02 0.2D-07 G 0.1D-01 0.2D-01 0.10D+01
13 20 0.217D-02 0.11D-01 0.12D-01 0.5D-07 G 0.3D-01 0.3D-01 0.28D+00
14 21 0.215D-02 0.77D-02 0.65D-02 0.2D-07 G 0.3D-01 0.1D-01 0.10D+01
15 22 0.212D-02 0.13D-01 0.19D-01 0.8D-07 S 0.2D-01 0.5D-01 0.35D+00
16 23 0.206D-02 0.28D-01 0.31D-01 0.7D-07 S 0.2D-01 0.4D-01 0.00D+00
17 25 0.202D-02 0.23D-01 0.24D-01 0.7D-07 S 0.6D-02 0.5D-01 0.00D+00
18 26 0.197D-02 0.23D-01 0.26D-01 0.4D-07 G 0.5D-02 0.7D-01 0.10D+01
19 27 0.192D-02 0.23D-01 0.55D-01 0.1D-06 G 0.6D-02 0.1D+00 0.10D+01
20 28 0.181D-02 0.56D-01 0.97D-01 0.9D-07 G 0.6D-02 0.1D+00 0.10D+01
21 29 0.168D-02 0.76D-01 0.12D+00 0.9D-07 G 0.6D-02 0.1D+00 0.10D+01
22 30 0.152D-02 0.95D-01 0.16D+00 0.1D-06 G 0.6D-02 0.1D+00 0.10D+01
23 31 0.130D-02 0.14D+00 0.18D+00 0.1D-06 G 0.6D-02 0.1D+00 0.10D+01
24 32 0.111D-02 0.15D+00 0.16D+00 0.1D-06 G 0.6D-02 0.1D+00 0.10D+01
25 33 0.943D-03 0.15D+00 0.15D+00 0.2D-06 G 0.5D-02 0.1D+00 0.10D+01
26 35 0.804D-03 0.15D+00 0.14D+00 0.2D-06 G 0.2D-02 0.1D+00 0.10D+01
27 37 0.747D-03 0.71D-01 0.72D-01 0.7D-07 G 0.2D-01 0.5D-01 0.10D+01
28 39 0.699D-03 0.63D-01 0.64D-01 0.8D-07 G 0.4D-02 0.5D-01 0.10D+01
29 40 0.622D-03 0.11D+00 0.24D+00 0.4D-06 G 0.4D-02 0.2D+00 0.10D+01
30 41 0.430D-03 0.31D+00 0.37D+00 0.5D-06 G 0.3D-02 0.2D+00 0.10D+01
31 42 0.364D-03 0.15D+00 0.45D+00 0.5D-06 S 0.6D-02 0.2D+00 0.62D+00
32 43 0.195D-03 0.46D+00 0.50D+00 0.5D-06 G 0.2D-02 0.2D+00 0.10D+01
33 44 0.117D-03 0.40D+00 0.43D+00 0.5D-06 G 0.2D-02 0.2D+00 0.10D+01
34 45 0.662D-04 0.44D+00 0.51D+00 0.5D-06 G 0.1D-02 0.2D+00 0.10D+01
35 47 0.325D-04 0.51D+00 0.69D+00 0.5D-06 S-G 0.6D-03 0.2D+00 0.10D+01
36 49 0.120D-04 0.63D+00 0.93D+00 0.5D-06 S-G 0.2D-03 0.2D+00 0.10D+01
37 53 0.328D-05 0.73D+00 0.15D+01 0.1D-08 S 0.1D+02 0.1D-02 0.00D+00
38 54 0.238D-05 0.27D+00 0.27D+00 0.2D-08 G 0.6D-01 0.1D-02 0.10D+01
39 58 0.404D-06 0.83D+00 0.97D+00 0.2D-06 G 0.1D-04 0.1D+00 0.10D+01
40 59 0.337D-07 0.92D+00 0.91D+00 0.2D-07 G 0.1D-06 0.6D-02 0.10D+01
41 60 0.318D-07 0.57D-01 0.42D-01 0.5D-08 G 0.6D-07 0.2D-02 0.10D+01
42 61 0.315D-07 0.71D-02 0.53D-02 0.2D-08 G 0.1D-06 0.7D-03 0.10D+01
43 62 0.315D-07 0.87D-03 0.66D-03 0.6D-09 G 0.3D-07 0.2D-03 0.10D+01
44 63 0.315D-07 0.13D-03 0.13D-03 0.3D-09 G 0.0D+00 0.1D-03 0.13D-03
45 67 0.315D-07 0.68D-07 0.13D-06 0.2D-12 G 0.2D+00 0.8D-07 0.10D+01
46 72 0.315D-07-0.44D-07 0.74D-08 0.8D-14 G 0.2D+02 0.3D-08 0.23D-04
***** FALSE CONVERGENCE *****
FUNCTION 0.315053D-07 RELDX 0.794D-14
FUNC. EVALS 72 GRAD. EVALS 460
PRELDF 0.737D-08 NPRELDF 0.227D-04
I FINAL X(I) D(I) G(I)
1 -0.195723D+06 0.100D+01 -0.141D-08
2 0.241284D-01 0.568D+00 0.430D-07
3 0.156657D+02 0.671D+00 -0.947D-08
4 0.317814D-01 0.903D+00 -0.114D-08
5 0.153081D+02 0.646D+00 0.442D-08
6 0.602895D-01 0.508D+00 -0.512D-09
7 0.691222D-01 0.457D+00 -0.595D-10
8 0.113276D+02 0.161D+01 -0.143D-07
9 -0.545592D+01 0.110D+00 -0.489D-09
10 -0.521868D+01 0.535D+00 0.446D-09
Parabolic equation of state: info 8
E = x1 + x2(a-x3)^2
+ x4(b-x5)^2 + x6(a-x3)(b-x5)
+ x7(c-x8)^2 + x9(a-x3)(c-x8) + x10(b-x5)(c-x8)
Fitparameter are
-195722.961991 0.024128 15.665728 0.031781
15.308083 0.060289 0.069122 11.327643
-5.455925 -5.218681
lattic parameters energy de(EOS)
15.063200 15.788120 11.225540 -195723.042487 -0.000001
14.912570 15.947590 11.338930 -195722.955531 -0.000094
15.063200 15.947590 11.338930 -195722.964215 0.000194
15.213830 15.947590 11.338930 -195722.971285 -0.000037
15.063200 15.788120 11.338930 -195722.954408 -0.000100
15.063200 16.107070 11.338930 -195722.971888 -0.000028
15.063200 15.947590 11.225540 -195722.957635 0.000000
15.063200 15.947590 11.452320 -195722.968629 0.000000
14.912570 15.788120 11.338930 -195722.944729 0.000065
14.912570 15.947590 11.225540 -195723.042427 0.000001
Sigma: 0.000084
Optionally create data points from fit function
Enter number of datapoints for your 3 dimensional Energy surface
NI=0 terminates; NI=1 will use 1 specific value in I-th component and
allows to
generate 2D-cuts
0.000u 0.004s 0:00.01 0.0% 0+0k 1952+8io 5pf+0w
_______________________________________________
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
--
P.Blaha
--------------------------------------------------------------------------
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
Phone: +43-1-58801-165300 FAX: +43-1-58801-165982
Email: bl...@theochem.tuwien.ac.at WIEN2k: http://www.wien2k.at
WWW: http://www.imc.tuwien.ac.at/staff/tc_group_e.php
--------------------------------------------------------------------------
_______________________________________________
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