abinit/tests/tutomultibinit/Refs/tmulti1_1.abo

407 lines
20 KiB
Plaintext

******************************************************************************************
Welcome to MULTIBINIT,
a software platform designed for the construction and use of second-principles models
for lattice, spin and electron degrees of freedom.
.Version 10.1.4.5 of MULTIBINIT
.(MPI version, prepared for a x86_64_linux_gnu13.2 computer)
.Copyright (C) 1998-2025 ABINIT group .
MULTIBINIT comes with ABSOLUTELY NO WARRANTY.
It is free software, and you are welcome to redistribute it
under certain conditions (GNU General Public License,
see ~abinit/COPYING or http://www.gnu.org/copyleft/gpl.txt).
MULTIBINIT is a software project of the University of Liege
(PHYTHEMA & NANOMAT groups), in collaboration with other partners.
-----------------------------------------------------------------------------------------
MULTIBINIT - LATTICE MODELS
Project initiated and coordinated by Philippe GHOSEZ and his group at ULiege
(Philippe.Ghosez@uliege.be).
Main contributors: Alexandre MARTIN, Jordan BIEDER, Michael Marcus SCHMITT,
Louis BASTOGNE, Xu HE, Alireza SASANI, Huazhang ZHANG, Subhadeep BANDYOPADHYAY,
Philippe GHOSEZ.
Technical support: Xu HE (X.He@uliege.be)
*****************************************************************************************
.Starting date : Fri 13 Sep 2024.
- ( at 19h02 )
- The starting date is more than 2 years after the initial release
- of this version of ABINIT, namely Jul 2022.
- Note that the use beyond 3 years after the release will not be supported.
- Action: please, switch to a more recent version of ABINIT.
- nproc = 1
================================================================================
Read the information in the reference structure in
-/home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/tutomultibinit/Input/tmulti1_DDB
to initialize the multibinit input
================================================================================
-outvars_multibinit: echo values of input variables ----------------------
Flags :
ifcflag 1
prt_model 1
strcpli -1
Miscellaneous information :
asr 2
Interatomic Force Constants Inputs :
dipdip 0
ifcana 0
ifcout 2000000
natifc 20
atifc 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
16 17 18 19 20
Description of grid 1 :
brav 1
ngqpt 2 2 2
nqshft 1
q1shft
0.00000000E+00 0.00000000E+00 0.00000000E+00
First list of wavevector (reduced coord.) :
nph1l 1
qph1l
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.000E+00
================================================================================
Read the DDB information of the reference system and perform some checks
==== Info on the Cryst% object ====
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
R(1)= 9.9768860 0.0000000 0.0000000 G(1)= 0.1002317 0.0000000 0.0000000
R(2)= 0.0000000 10.1890339 0.0000000 G(2)= 0.0000000 0.0981447 0.0000000
R(3)= 0.0000000 0.0000000 14.2034780 G(3)= 0.0000000 0.0000000 0.0704053
Unit cell volume ucvol= 1.4438521E+03 bohr^3
Angles (23,13,12)= 9.00000000E+01 9.00000000E+01 9.00000000E+01 degrees
Time-reversal symmetry is present
Reduced atomic positions [iatom, xred, symbol]:
1) -0.0102018 0.0466924 0.2500000 Ca
2) 0.0102018 0.9533076 0.7500000 Ca
3) 0.5102018 0.5466924 0.2500000 Ca
4) 0.4897982 0.4533076 0.7500000 Ca
5) 0.0000000 0.5000000 0.0000000 Ti
6) 0.0000000 0.5000000 0.5000000 Ti
7) 0.5000000 -0.0000000 0.0000000 Ti
8) 0.5000000 -0.0000000 0.5000000 Ti
9) 0.0815525 0.4788799 0.2500000 O
10) 0.9184475 0.5211201 0.7500000 O
11) 0.5815525 0.0211201 0.7500000 O
12) 0.4184475 0.9788799 0.2500000 O
13) 0.7072278 0.2919715 0.0429245 O
14) 0.2927722 0.7080285 0.9570755 O
15) 0.2927722 0.7080285 0.5429245 O
16) 0.7072278 0.2919715 0.4570755 O
17) 0.7927722 0.7919715 0.0429245 O
18) 0.2072278 0.2080285 0.9570755 O
19) 0.2072278 0.2080285 0.5429245 O
20) 0.7927722 0.7919715 0.4570755 O
DDB file with 10 blocks has been read.
================================================================================
Extraction of the energy of the structure (unit: Hartree)
Energy = -5.785439537588E+02
================================================================================
Dielectric Tensor and Effective Charges
anaddb : Zero the imaginary part of the Dynamical Matrix at Gamma,
and impose the ASR on the effective charges
The violation of the charge neutrality conditions
by the effective charges is as follows :
atom electric field
displacement direction
1 1 -2.403821 0.000000
1 2 -0.000000 0.000000
1 3 0.000000 0.000000
2 1 0.000000 0.000000
2 2 -1.866807 0.000000
2 3 0.000000 0.000000
3 1 -0.000000 0.000000
3 2 0.000000 0.000000
3 3 -1.422018 0.000000
Effective charge tensors after
imposition of the charge neutrality (if requested by user),
and eventual restriction to some part :
atom displacement
1 1 1.453433E+00 8.057855E-02 -1.665606E-17
1 2 3.625397E-02 1.481414E+00 -2.775558E-18
1 3 5.507747E-18 -1.387779E-18 1.570687E+00
2 1 1.453433E+00 8.057855E-02 -1.665606E-17
2 2 3.625397E-02 1.481414E+00 -2.775558E-18
2 3 5.507747E-18 -1.387779E-18 1.570687E+00
3 1 1.453433E+00 -8.057855E-02 -1.665606E-17
3 2 -3.625397E-02 1.481414E+00 -2.775558E-18
3 3 5.507747E-18 -1.387779E-18 1.570687E+00
4 1 1.453433E+00 -8.057855E-02 -1.665606E-17
4 2 -3.625397E-02 1.481414E+00 -2.775558E-18
4 3 5.507747E-18 -1.387779E-18 1.570687E+00
5 1 7.204839E+00 5.113710E-01 7.438256E-01
5 2 -6.876662E-01 7.319820E+00 -2.248431E-01
5 3 -9.178499E-01 -6.062353E-02 7.272712E+00
6 1 7.204839E+00 5.113710E-01 -7.438256E-01
6 2 -6.876662E-01 7.319820E+00 2.248431E-01
6 3 9.178499E-01 6.062353E-02 7.272712E+00
7 1 7.204839E+00 -5.113710E-01 -7.438256E-01
7 2 6.876662E-01 7.319820E+00 -2.248431E-01
7 3 9.178499E-01 -6.062353E-02 7.272712E+00
8 1 7.204839E+00 -5.113710E-01 7.438256E-01
8 2 6.876662E-01 7.319820E+00 2.248431E-01
8 3 -9.178499E-01 6.062353E-02 7.272712E+00
9 1 -1.939073E+00 -1.649446E-01 -1.665606E-17
9 2 2.208078E-02 -1.754548E+00 -2.775558E-18
9 3 5.507747E-18 -1.387779E-18 -5.118020E+00
10 1 -1.939073E+00 -1.649446E-01 -1.665606E-17
10 2 2.208078E-02 -1.754548E+00 -2.775558E-18
10 3 5.507747E-18 -1.387779E-18 -5.118020E+00
11 1 -1.939073E+00 1.649446E-01 -1.665606E-17
11 2 -2.208078E-02 -1.754548E+00 -2.775558E-18
11 3 5.507747E-18 -1.387779E-18 -5.118020E+00
12 1 -1.939073E+00 1.649446E-01 -1.665606E-17
12 2 -2.208078E-02 -1.754548E+00 -2.775558E-18
12 3 5.507747E-18 -1.387779E-18 -5.118020E+00
13 1 -3.359599E+00 -1.623543E+00 -2.028591E-04
13 2 -1.574037E+00 -3.523343E+00 -9.997223E-02
13 3 2.370312E-03 -7.357064E-02 -1.862689E+00
14 1 -3.359599E+00 -1.623543E+00 -2.028591E-04
14 2 -1.574037E+00 -3.523343E+00 -9.997223E-02
14 3 2.370312E-03 -7.357064E-02 -1.862689E+00
15 1 -3.359599E+00 -1.623543E+00 2.028591E-04
15 2 -1.574037E+00 -3.523343E+00 9.997223E-02
15 3 -2.370312E-03 7.357064E-02 -1.862689E+00
16 1 -3.359599E+00 -1.623543E+00 2.028591E-04
16 2 -1.574037E+00 -3.523343E+00 9.997223E-02
16 3 -2.370312E-03 7.357064E-02 -1.862689E+00
17 1 -3.359599E+00 1.623543E+00 2.028591E-04
17 2 1.574037E+00 -3.523343E+00 -9.997223E-02
17 3 -2.370312E-03 -7.357064E-02 -1.862689E+00
18 1 -3.359599E+00 1.623543E+00 2.028591E-04
18 2 1.574037E+00 -3.523343E+00 -9.997223E-02
18 3 -2.370312E-03 -7.357064E-02 -1.862689E+00
19 1 -3.359599E+00 1.623543E+00 -2.028591E-04
19 2 1.574037E+00 -3.523343E+00 9.997223E-02
19 3 2.370312E-03 7.357064E-02 -1.862689E+00
20 1 -3.359599E+00 1.623543E+00 -2.028591E-04
20 2 1.574037E+00 -3.523343E+00 9.997223E-02
20 3 2.370312E-03 7.357064E-02 -1.862689E+00
Now, the imaginary part of the dynamical matrix is zeroed
================================================================================
Extraction of the stress tensor (unit: GPa) and forces (unit: Ha/bohr)
Cartesian components of forces (hartree/bohr)
1 -0.30800054E-05 -0.10255015E-04 -0.31554436E-29
2 0.30800054E-05 0.10255015E-04 -0.15777218E-29
3 0.30800054E-05 -0.10255015E-04 0.00000000E+00
4 -0.30800054E-05 0.10255015E-04 0.00000000E+00
5 0.00000000E+00 0.00000000E+00 0.00000000E+00
6 0.00000000E+00 0.00000000E+00 0.00000000E+00
7 0.00000000E+00 0.00000000E+00 0.00000000E+00
8 0.00000000E+00 0.00000000E+00 0.00000000E+00
9 -0.19666635E-03 -0.33774783E-03 0.00000000E+00
10 0.19666635E-03 0.33774783E-03 0.00000000E+00
11 -0.19666635E-03 0.33774783E-03 0.00000000E+00
12 0.19666635E-03 -0.33774783E-03 0.00000000E+00
13 0.24085511E-03 -0.36542027E-03 0.13299857E-03
14 -0.24085511E-03 0.36542027E-03 -0.13299857E-03
15 -0.24085511E-03 0.36542027E-03 0.13299857E-03
16 0.24085511E-03 -0.36542027E-03 -0.13299857E-03
17 -0.24085511E-03 -0.36542027E-03 0.13299857E-03
18 0.24085511E-03 0.36542027E-03 -0.13299857E-03
19 0.24085511E-03 0.36542027E-03 0.13299857E-03
20 -0.24085511E-03 -0.36542027E-03 -0.13299857E-03
Cartesian components of stress tensor (hartree/bohr^3)
sigma(1 1)= -1.27676002E-07 sigma(3 2)= 0.00000000E+00
sigma(2 2)= -5.42027786E-09 sigma(3 1)= 0.00000000E+00
sigma(3 3)= -3.33970359E-08 sigma(2 1)= 0.00000000E+00
================================================================================
Extraction of the clamped elastic tensor (unit:10^2GPa)
3.9133823 1.3361329 1.1936863 -0.0000000 -0.0000000 -0.0000000
1.3361172 3.8023382 1.1335897 -0.0000000 -0.0000000 -0.0000000
1.1936770 1.1335796 4.0167474 -0.0000000 -0.0000000 -0.0000000
-0.0000000 -0.0000000 0.0000000 1.2884614 -0.0000000 0.0000000
-0.0000000 0.0000000 -0.0000000 -0.0000000 1.2158450 0.0000000
-0.0000000 -0.0000000 -0.0000000 0.0000000 -0.0000000 1.4230204
================================================================================
Calculation of acoustic sum rule
================================================================================
Calculation of the interatomic forces from DDB
Homogeneous q point set in the B.Z.
Grid q points : 8
1) 0.00000000E+00 0.00000000E+00 0.00000000E+00
2) 5.00000000E-01 0.00000000E+00 0.00000000E+00
3) 0.00000000E+00 5.00000000E-01 0.00000000E+00
4) 5.00000000E-01 5.00000000E-01 0.00000000E+00
5) 0.00000000E+00 0.00000000E+00 5.00000000E-01
6) 5.00000000E-01 0.00000000E+00 5.00000000E-01
7) 0.00000000E+00 5.00000000E-01 5.00000000E-01
8) 5.00000000E-01 5.00000000E-01 5.00000000E-01
The interatomic forces have been obtained
================================================================================
Calculation of dynamical matrix for each ph1l points
Phonon at Gamma, with non-analyticity in the
direction (cartesian coordinates) 0.00000 0.00000 0.00000
Phonon energies in Hartree :
-2.294305E-06 -1.589880E-06 -1.457514E-06 4.403333E-04 4.476135E-04
5.009874E-04 5.081990E-04 6.094765E-04 6.228453E-04 6.841952E-04
7.149261E-04 7.563444E-04 7.823985E-04 8.187584E-04 8.503830E-04
8.756807E-04 9.115740E-04 9.303034E-04 9.451129E-04 9.477895E-04
1.008750E-03 1.073909E-03 1.087361E-03 1.094627E-03 1.194544E-03
1.215171E-03 1.243183E-03 1.299532E-03 1.326571E-03 1.393394E-03
1.420781E-03 1.496467E-03 1.550377E-03 1.629425E-03 1.632909E-03
1.672657E-03 1.693239E-03 1.740117E-03 1.970287E-03 1.976703E-03
2.016032E-03 2.076543E-03 2.094483E-03 2.116629E-03 2.126912E-03
2.233949E-03 2.246967E-03 2.248255E-03 2.250899E-03 2.332700E-03
2.346309E-03 2.378494E-03 2.382681E-03 2.422239E-03 2.441826E-03
2.465514E-03 2.565346E-03 3.445785E-03 3.566167E-03 3.726860E-03
Phonon frequencies in cm-1 :
- -5.035417E-01 -3.489383E-01 -3.198875E-01 9.664199E+01 9.823980E+01
- 1.099540E+02 1.115368E+02 1.337646E+02 1.366987E+02 1.501635E+02
- 1.569081E+02 1.659984E+02 1.717166E+02 1.796967E+02 1.866375E+02
- 1.921897E+02 2.000674E+02 2.041780E+02 2.074283E+02 2.080157E+02
- 2.213950E+02 2.356958E+02 2.386481E+02 2.402428E+02 2.621722E+02
- 2.666991E+02 2.728472E+02 2.852143E+02 2.911487E+02 3.058145E+02
- 3.118255E+02 3.284365E+02 3.402684E+02 3.576175E+02 3.583822E+02
- 3.671058E+02 3.716230E+02 3.819116E+02 4.324281E+02 4.338362E+02
- 4.424679E+02 4.557485E+02 4.596859E+02 4.645464E+02 4.668033E+02
- 4.902951E+02 4.931522E+02 4.934350E+02 4.940152E+02 5.119685E+02
- 5.149552E+02 5.220192E+02 5.229380E+02 5.316199E+02 5.359190E+02
- 5.411179E+02 5.630283E+02 7.562623E+02 7.826833E+02 8.179513E+02
================================================================================
Calculation of the internal-strain tensor
Force-response internal strain tensor(Unit:Hartree/bohr)
Atom dir strainxx strainyy strainzz strainyz strainxz strainxy
1 x 0.0941002 -0.0408713 -0.0533039 0.0000000 0.0000000 -0.0609607
1 y -0.0828675 0.1403061 -0.0248158 0.0000000 0.0000000 -0.0169793
1 z 0.0000000 0.0000000 0.0000000 -0.0368215 -0.0711707 -0.0000000
2 x -0.0941002 0.0408713 0.0533039 -0.0000000 -0.0000000 0.0609607
2 y 0.0828675 -0.1403061 0.0248158 -0.0000000 -0.0000000 0.0169793
2 z 0.0000000 0.0000000 0.0000000 0.0368215 0.0711707 -0.0000000
3 x -0.0941002 0.0408713 0.0533039 -0.0000000 -0.0000000 -0.0609607
3 y -0.0828675 0.1403061 -0.0248158 0.0000000 0.0000000 0.0169793
3 z 0.0000000 0.0000000 0.0000000 -0.0368215 0.0711707 0.0000000
4 x 0.0941002 -0.0408713 -0.0533039 0.0000000 0.0000000 0.0609607
4 y 0.0828675 -0.1403061 0.0248158 -0.0000000 -0.0000000 -0.0169793
4 z 0.0000000 0.0000000 0.0000000 0.0368215 -0.0711707 0.0000000
5 x -0.0000000 0.0000000 0.0000000 0.0000000 -0.0000000 0.0000000
5 y 0.0000000 0.0000000 0.0000000 -0.0000000 0.0000000 -0.0000000
5 z -0.0000000 -0.0000000 -0.0000000 0.0000000 0.0000000 0.0000000
6 x -0.0000000 0.0000000 0.0000000 -0.0000000 -0.0000000 0.0000000
6 y 0.0000000 0.0000000 0.0000000 -0.0000000 -0.0000000 -0.0000000
6 z -0.0000000 -0.0000000 -0.0000000 0.0000000 -0.0000000 -0.0000000
7 x -0.0000000 -0.0000000 -0.0000000 -0.0000000 -0.0000000 0.0000000
7 y 0.0000000 0.0000000 0.0000000 -0.0000000 -0.0000000 -0.0000000
7 z -0.0000000 -0.0000000 -0.0000000 0.0000000 -0.0000000 -0.0000000
8 x 0.0000000 -0.0000000 -0.0000000 0.0000000 -0.0000000 0.0000000
8 y 0.0000000 0.0000000 0.0000000 -0.0000000 -0.0000000 -0.0000000
8 z -0.0000000 -0.0000000 -0.0000000 0.0000000 -0.0000000 0.0000000
9 x 0.0165308 -0.0318013 -0.1284522 -0.0000000 -0.0000000 0.0027962
9 y 0.0090276 -0.0682327 0.0203506 0.0000000 0.0000000 -0.0401272
9 z -0.0000000 -0.0000000 -0.0000000 0.0337670 -0.0716536 -0.0000000
10 x -0.0165308 0.0318013 0.1284522 0.0000000 0.0000000 -0.0027962
10 y -0.0090276 0.0682327 -0.0203506 -0.0000000 -0.0000000 0.0401272
10 z -0.0000000 -0.0000000 -0.0000000 -0.0337670 0.0716536 -0.0000000
11 x 0.0165308 -0.0318013 -0.1284522 -0.0000000 -0.0000000 -0.0027962
11 y -0.0090276 0.0682327 -0.0203506 -0.0000000 -0.0000000 -0.0401272
11 z -0.0000000 -0.0000000 -0.0000000 -0.0337670 -0.0716536 0.0000000
12 x -0.0165308 0.0318013 0.1284522 0.0000000 0.0000000 0.0027962
12 y 0.0090276 -0.0682327 0.0203506 0.0000000 0.0000000 0.0401272
12 z -0.0000000 -0.0000000 -0.0000000 0.0337670 0.0716536 0.0000000
13 x 0.1036986 0.0028863 -0.0013285 -0.0294432 -0.0159880 0.0323406
13 y 0.0023604 -0.0644718 0.0033325 -0.0085421 -0.0319129 -0.0243703
13 z -0.0364388 -0.0301746 -0.0242418 0.0104449 0.0047313 -0.0546885
14 x -0.1036986 -0.0028863 0.0013285 0.0294432 0.0159880 -0.0323406
14 y -0.0023604 0.0644718 -0.0033325 0.0085421 0.0319129 0.0243703
14 z 0.0364388 0.0301746 0.0242418 -0.0104449 -0.0047313 0.0546885
15 x -0.1036986 -0.0028863 0.0013285 -0.0294432 -0.0159880 -0.0323406
15 y -0.0023604 0.0644718 -0.0033325 -0.0085421 -0.0319129 0.0243703
15 z -0.0364388 -0.0301746 -0.0242418 -0.0104449 -0.0047313 -0.0546885
16 x 0.1036986 0.0028863 -0.0013285 0.0294432 0.0159880 0.0323406
16 y 0.0023604 -0.0644718 0.0033325 0.0085421 0.0319129 -0.0243703
16 z 0.0364388 0.0301746 0.0242418 0.0104449 0.0047313 0.0546885
17 x -0.1036986 -0.0028863 0.0013285 0.0294432 -0.0159880 0.0323406
17 y 0.0023604 -0.0644718 0.0033325 -0.0085421 0.0319129 0.0243703
17 z -0.0364388 -0.0301746 -0.0242418 0.0104449 -0.0047313 0.0546885
18 x 0.1036986 0.0028863 -0.0013285 -0.0294432 0.0159880 -0.0323406
18 y -0.0023604 0.0644718 -0.0033325 0.0085421 -0.0319129 -0.0243703
18 z 0.0364388 0.0301746 0.0242418 -0.0104449 0.0047313 -0.0546885
19 x 0.1036986 0.0028863 -0.0013285 0.0294432 -0.0159880 -0.0323406
19 y -0.0023604 0.0644718 -0.0033325 -0.0085421 0.0319129 -0.0243703
19 z -0.0364388 -0.0301746 -0.0242418 -0.0104449 0.0047313 0.0546885
20 x -0.1036986 -0.0028863 0.0013285 -0.0294432 0.0159880 0.0323406
20 y 0.0023604 -0.0644718 0.0033325 0.0085421 -0.0319129 0.0243703
20 z 0.0364388 0.0301746 0.0242418 0.0104449 -0.0047313 -0.0546885
Bound for ifc SR:
x=[ -1 1], y=[ -1 1] and z=[ -1 1]
================================================================================
Impose acoustic sum rule on total ifc
================================================================================
================================================================================
There is no file for the coefficients from polynomial fitting
================================================================================
Generation of the xml file for the model in tmulti1_1_model.xml
================================================================================
-
- Proc. 0 individual time (sec): cpu= 0.3 wall= 0.3
================================================================================
+Total cpu time 0.295 and wall time 0.329 sec
multibinit : the run completed succesfully.