******************************************************************************************
                                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.6.5 of MULTIBINIT 
.(MPI version, prepared for a x86_64_linux_gnu11.4 computer) 

.Copyright (C) 1998-2026 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 : Tue  5 May 2026.
- ( at 13h08 )
  
- The starting date is more than 3 years after the initial release
- of this version of ABINIT, namely Jul 2022.
- This version of ABINIT is not supported anymore.
- Action: please, switch to a more recent version of ABINIT.
-  nproc =    1

================================================================================

 Read the information in the reference structure in 
-/home/hexu/projects/atomchain_dev/atomchain/.tmp/batio3_ddb_example/BaTiO3_phonon.ddb
 to initialize the multibinit input

================================================================================

 -outvars_multibinit: echo values of input variables ----------------------

 Flags : 
     ifcflag         1
   prt_model         4
     strcpli        -1
 Fit the coefficients :
         fit_coeff  1
 fit_generateCoeff  1
 fit_initializeDat  0
        fit_cutoff  8.00000000E+00
       fit_droprat  0.00000000E+00
        fit_option  0
         fit_iatom  0
        fit_ncoeff  3
          fit_grid  1  1  1
         ts_option  0
    fit_rangePower  3  3
     fit_dispterms  1
    fit_anhaStrain  0
    fit_SPCoupling  0
      fit_SPC_maxS  1
     fit_max_nbody        999
 Interatomic Force Constants Inputs :
      dipdip         0
      ifcana         0
      ifcout   2000000
      natifc         5
       atifc         1   2   3   4   5
 Description of grid 1 :
        brav         1
       ngqpt         1         1         1
      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)=  7.5589045  0.0000000  0.0000000  G(1)=  0.1322943  0.0000000  0.0000000
 R(2)=  0.0000000  7.5589045  0.0000000  G(2)=  0.0000000  0.1322943  0.0000000
 R(3)=  0.0000000  0.0000000  7.5589045  G(3)=  0.0000000  0.0000000  0.1322943
 Unit cell volume ucvol=  4.3189341E+02 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.0000000  0.0000000  0.0000000  Ba
    2)    0.5000000  0.5000000  0.5000000  Ti
    3)    0.5000000  0.5000000  0.0000000   O
    4)    0.5000000  0.0000000  0.5000000   O
    5)    0.0000000  0.5000000  0.5000000   O

 DDB file with 1 blocks has been read.

================================================================================

 Extraction of the energy of the structure (unit: Hartree)

 Energy  =        0.000000000000E+00

================================================================================

 Extraction of the stress tensor (unit: GPa) and forces (unit: Ha/bohr)

 --- !WARNING:
     The stress tensor of the reference structure is not specify
     The stress tensor will be set to zero
 ---

================================================================================

 Extraction of the clamped elastic tensor (unit:10^2GPa)


 Warning : Elastic Tensor is set to zero (not available in the DDB)

================================================================================

 Calculation of acoustic sum rule


================================================================================

 Calculation of the interatomic forces from DDB

       Homogeneous q point set in the B.Z.
 Grid q points  :        1
  1)   0.00000000E+00  0.00000000E+00  0.00000000E+00

 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 :
  -1.233234E-04 -1.233234E-04 -1.233234E-04  6.667508E-08  6.670936E-08
   6.681207E-08  1.064943E-04  1.064943E-04  1.064943E-04  1.748844E-04
   1.748844E-04  1.748844E-04  2.831090E-04  2.831090E-04  2.831090E-04
 Phonon frequencies in cm-1    :
- -2.706635E+01 -2.706635E+01 -2.706635E+01  1.463349E-02  1.464101E-02
-  1.466355E-02  2.337279E+01  2.337279E+01  2.337279E+01  3.838269E+01
-  3.838269E+01  3.838269E+01  6.213525E+01  6.213525E+01  6.213525E+01

================================================================================

 Calculation of the internal-strain  tensor

 Warning : Internal strain is set to zero (not available in the DDB)



 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

================================================================================

-Reading the training-set file :
-/home/hexu/projects/atomchain_dev/pymultibinit/examples/BaTiO3_training/real_training_run/BaTiO3_multibinit_HIST.nc


================================================================================

  Starting Fit Iterations  
  -----------------------  
  Select in total fit_ncoeff =   3 coefficients
  In   1 iterations
  Over   3 irreducible atoms
  Selecting   1 coefficients per atom in each iteration


--------------------------------------------------------------------------------

 Start Iteration (  1/  1)


--------------------------------------------------------------------------------

 Starting Fit process

--------------------------------------------------------------------------------

 The coefficients for the fit around atom  1: Ba, will be generated

 1731 coefficients generated 


 Goal function values at the begining of the fit process (eV^2/A^2):
   Energy          :   1.2505720487946192E-05
   Forces+Stresses :   4.8013172522383104E-04
   Forces          :   4.8013090841234036E-04
   Stresses        :   8.1681149064418615E-10

 N   Selecting          MSDE              MSDFS              MSDF              MSDS
    Coefficient      (eV^2/A^2)        (eV^2/A^2)        (eV^2/A^2)        (eV^2/A^2)
 1     93          1.1148226179E-05  4.5261271669E-04  4.5260202712E-04  1.0689569560E-08

 Fitted coefficients at the end of the fit process:
 93 =>   1.4069158084E-02 (Ba_x-O1_x)^1(Ba_y-O3_y)^1(Ba_y-O2_y[-1 0 -1])^1

 Goal function values at the end of the fit process (eV^2/A^2):
   Energy          :   1.1148226178581105E-05
   Forces+Stresses :   4.5261271669135743E-04
   Forces          :   4.5260202712179729E-04
   Stresses        :   1.0689569560144761E-08


--------------------------------------------------------------------------------

 Starting Fit process

--------------------------------------------------------------------------------

 The coefficients present in the effective potential will be used for the fit

 The coefficients for the fit around atom  2: Ti, will be generated

 738 coefficients generated 


 Goal function values at the begining of the fit process (eV^2/A^2):
   Energy          :   1.2505720487946192E-05
   Forces+Stresses :   4.8013172522383104E-04
   Forces          :   4.8013090841234036E-04
   Stresses        :   8.1681149064418615E-10

 N   Selecting          MSDE              MSDFS              MSDF              MSDS
    Coefficient      (eV^2/A^2)        (eV^2/A^2)        (eV^2/A^2)        (eV^2/A^2)
 2     194         9.7468720578E-06  4.1953096036E-04  4.1950125711E-04  2.9703245658E-08

 Fitted coefficients at the end of the fit process:
 1 =>   1.3522411855E-02 (Ba_x-O1_x)^1(Ba_y-O3_y)^1(Ba_y-O2_y[-1 0 -1])^1
 194 =>   1.7316101467E-02 (Ti_x-O1_x)^1(Ti_x-Ba_x[1 1 1])^1(Ti_z-Ba_z[0 1 1])^1

 Goal function values at the end of the fit process (eV^2/A^2):
   Energy          :   9.7468720577631271E-06
   Forces+Stresses :   4.1953096035761384E-04
   Forces          :   4.1950125711195593E-04
   Stresses        :   2.9703245657946742E-08


--------------------------------------------------------------------------------

 Starting Fit process

--------------------------------------------------------------------------------

 The coefficients present in the effective potential will be used for the fit

 The coefficients for the fit around atom  3: O1, will be generated

 1387 coefficients generated 


 Goal function values at the begining of the fit process (eV^2/A^2):
   Energy          :   1.2505720487946192E-05
   Forces+Stresses :   4.8013172522383104E-04
   Forces          :   4.8013090841234036E-04
   Stresses        :   8.1681149064418615E-10

 N   Selecting          MSDE              MSDFS              MSDF              MSDS
    Coefficient      (eV^2/A^2)        (eV^2/A^2)        (eV^2/A^2)        (eV^2/A^2)
 3     1110        8.8941912966E-06  3.9344607095E-04  3.9340185195E-04  4.4218997661E-08

 Fitted coefficients at the end of the fit process:
 1 =>   1.2619831978E-02 (Ba_x-O1_x)^1(Ba_y-O3_y)^1(Ba_y-O2_y[-1 0 -1])^1
 2 =>   1.7069657781E-02 (Ti_x-O1_x)^1(Ti_x-Ba_x[1 1 1])^1(Ti_z-Ba_z[0 1 1])^1
 1110 =>  -2.6980868933E-02 (O1_x-O2_x)^1(O1_y-O2_y[0 1 0])^1(O1_x-O1_x[0 1 0])^1

 Goal function values at the end of the fit process (eV^2/A^2):
   Energy          :   8.8941912966286722E-06
   Forces+Stresses :   3.9344607094673800E-04
   Forces          :   3.9340185194907675E-04
   Stresses        :   4.4218997661284441E-08


================================================================================

 Generation of the xml file for the fitted polynomial in BaTiO3_fit_coeffs.xml

================================================================================

-
- Proc.   0 individual time (sec): cpu=          9.6  wall=          9.6

================================================================================

+Total cpu time      9.594  and wall time      9.597 sec

 multibinit : the run completed successfully.
