7035
7522
<dData1D:(new from dSpData1D) n=5616, refcount: 1>
7037
7524
new_DM -- step: 63
7038
Initializing Density Matrix...
7039
DM filled with atomic data:
7040
<dSpData2D:DM initialized from atoms
7525
Re-using DM from previous geometries...
7526
Number of DMs in history: 1
7527
DM extrapolation coefficients:
7529
New DM after history re-use:
7530
<dSpData2D:SpM extrapolated using coords
7041
7531
<sparsity:sparsity for geom step 63
7042
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 8>
7043
<dData2D:DM n=5616 m=1, refcount: 1>
7046
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7047
New grid distribution: 1
7050
3 1: 24 1: 12 13: 18
7051
4 1: 24 1: 12 19: 24
7053
6 1: 24 13: 24 7: 12
7054
7 1: 24 13: 24 13: 18
7055
8 1: 24 13: 24 19: 24
7057
InitMesh: MESH = 48 x 48 x 48 = 110592
7058
InitMesh: (bp) = 24 x 24 x 24 = 13824
7059
InitMesh: Mesh cutoff (required, used) = 40.000 53.991 Ry
7060
ExtMesh (bp) on 0 = 48 x 36 x 30 = 51840
7061
New grid distribution: 2
7062
1 1: 12 13: 24 13: 24
7063
2 13: 24 1: 12 1: 12
7064
3 13: 24 1: 12 13: 24
7065
4 13: 24 13: 24 1: 12
7066
5 1: 12 1: 12 13: 24
7068
7 13: 24 13: 24 13: 24
7069
8 1: 12 13: 24 1: 12
7070
New grid distribution: 3
7071
1 1: 12 13: 24 13: 24
7072
2 13: 24 1: 12 1: 12
7073
3 13: 24 1: 12 13: 24
7074
4 13: 24 13: 24 1: 12
7075
5 1: 12 1: 12 13: 24
7077
7 13: 24 13: 24 13: 24
7078
8 1: 12 13: 24 1: 12
7079
Setting up quadratic distribution...
7080
ExtMesh (bp) on 0 = 36 x 36 x 36 = 46656
7081
PhiOnMesh: Number of (b)points on node 0 = 1728
7082
PhiOnMesh: nlist on node 0 = 32411
7084
siesta: Eharris(eV) = -6832.124735
7087
siesta: Eharris(eV) = -6832.124735
7089
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
7090
Geom step, scf iteration, dmax: 63 1 0.863540
7092
Using DM_out to compute the final energy and forces
7093
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7532
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 9>
7533
<dData2D:(temp array for extrapolation) n=5616 m=1, refcount: 1>
7535
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7536
New grid distribution: 1
7539
3 1: 24 1: 12 13: 18
7540
4 1: 24 1: 12 19: 24
7542
6 1: 24 13: 24 7: 12
7543
7 1: 24 13: 24 13: 18
7544
8 1: 24 13: 24 19: 24
7546
InitMesh: MESH = 48 x 48 x 48 = 110592
7547
InitMesh: (bp) = 24 x 24 x 24 = 13824
7548
InitMesh: Mesh cutoff (required, used) = 40.000 53.991 Ry
7549
ExtMesh (bp) on 0 = 48 x 36 x 30 = 51840
7550
New grid distribution: 2
7551
1 1: 12 13: 24 13: 24
7552
2 13: 24 1: 12 1: 12
7553
3 13: 24 1: 12 13: 24
7554
4 13: 24 13: 24 1: 12
7555
5 1: 12 1: 12 13: 24
7557
7 13: 24 13: 24 13: 24
7558
8 1: 12 13: 24 1: 12
7559
New grid distribution: 3
7560
1 1: 12 13: 24 13: 24
7561
2 13: 24 1: 12 1: 12
7562
3 13: 24 1: 12 13: 24
7563
4 13: 24 13: 24 1: 12
7564
5 1: 12 1: 12 13: 24
7566
7 13: 24 13: 24 13: 24
7567
8 1: 12 13: 24 1: 12
7568
Setting up quadratic distribution...
7569
ExtMesh (bp) on 0 = 36 x 36 x 36 = 46656
7570
PhiOnMesh: Number of (b)points on node 0 = 1728
7571
PhiOnMesh: nlist on node 0 = 32391
7573
iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) dHmax(eV)
7574
scf: 1 -6830.881918 -6830.881662 -6830.916320 0.000089 -1.331521 0.001762
7575
scf: 2 -6830.881640 -6830.881657 -6830.916314 0.000139 -1.331070 0.003497
7576
scf: 3 -6830.881667 -6830.881664 -6830.916322 0.000092 -1.331297 0.000107
7578
SCF Convergence by DM+H criterion
7579
max |DM_out - DM_in| : 0.0000922096
7580
max |H_out - H_in| (eV) : 0.0001066803
7581
SCF cycle converged after 3 iterations
7583
Using DM_out to compute the final energy and forces
7584
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7586
siesta: E_KS(eV) = -6830.8817
7588
siesta: Atomic forces (eV/Ang):
7589
----------------------------------------
7590
Tot 0.000000 -0.000000 -0.000000
7591
----------------------------------------
7593
Res 0.005886 sqrt( Sum f_i^2 / 3N )
7594
----------------------------------------
7595
Max 0.015040 constrained
7597
Stress-tensor-Voigt (kbar): -51.24 -51.24 -51.24 -3.20 -3.20 -3.20
7598
(Free)E + p*V (eV/cell) -6789.9518
7599
Target enthalpy (eV/cell) -6830.9163
7601
siesta: Stress tensor (static) (eV/Ang**3):
7602
-0.031983 -0.001997 -0.001997
7603
-0.001997 -0.031983 -0.001997
7604
-0.001997 -0.001997 -0.031983
7606
siesta: Pressure (static): 51.24289348 kBar
7608
siesta: Stress tensor (total) (eV/Ang**3):
7609
-0.031983 -0.001997 -0.001997
7610
-0.001997 -0.031983 -0.001997
7611
-0.001997 -0.001997 -0.031983
7613
siesta: Pressure (total): 51.24292760 kBar
7614
FIRE: istep, dt, firenpos, firealf, magv
7615
FIRE: 63 3.6289 2 0.1000 0.0007
7617
siesta: Temp_ion = 0.016 K
7619
====================================
7620
Begin MD quenched step = 64
7621
====================================
7623
outcell: Unit cell vectors (Ang):
7624
10.860000 0.000000 0.000000
7625
0.000000 10.860000 0.000000
7626
0.000000 0.000000 10.860000
7628
outcell: Cell vector modules (Ang) : 10.860000 10.860000 10.860000
7629
outcell: Cell angles (23,13,12) (deg): 90.0000 90.0000 90.0000
7630
outcell: Cell volume (Ang**3) : 1280.8241
7631
<dSpData1D:S at geom step 64
7632
<sparsity:sparsity for geom step 64
7633
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 7>
7634
<dData1D:(new from dSpData1D) n=5616, refcount: 1>
7637
Re-using DM from previous geometries...
7638
Number of DMs in history: 1
7639
DM extrapolation coefficients:
7641
New DM after history re-use:
7642
<dSpData2D:SpM extrapolated using coords
7643
<sparsity:sparsity for geom step 64
7644
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 9>
7645
<dData2D:(temp array for extrapolation) n=5616 m=1, refcount: 1>
7647
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7648
New grid distribution: 1
7651
3 1: 24 1: 12 13: 18
7652
4 1: 24 1: 12 19: 24
7654
6 1: 24 13: 24 7: 12
7655
7 1: 24 13: 24 13: 18
7656
8 1: 24 13: 24 19: 24
7658
InitMesh: MESH = 48 x 48 x 48 = 110592
7659
InitMesh: (bp) = 24 x 24 x 24 = 13824
7660
InitMesh: Mesh cutoff (required, used) = 40.000 53.991 Ry
7661
ExtMesh (bp) on 0 = 48 x 36 x 30 = 51840
7662
New grid distribution: 2
7663
1 1: 12 13: 24 13: 24
7664
2 13: 24 1: 12 1: 12
7665
3 13: 24 1: 12 13: 24
7666
4 13: 24 13: 24 1: 12
7667
5 1: 12 1: 12 13: 24
7669
7 13: 24 13: 24 13: 24
7670
8 1: 12 13: 24 1: 12
7671
New grid distribution: 3
7672
1 1: 12 13: 24 13: 24
7673
2 13: 24 1: 12 1: 12
7674
3 13: 24 1: 12 13: 24
7675
4 13: 24 13: 24 1: 12
7676
5 1: 12 1: 12 13: 24
7678
7 13: 24 13: 24 13: 24
7679
8 1: 12 13: 24 1: 12
7680
Setting up quadratic distribution...
7681
ExtMesh (bp) on 0 = 36 x 36 x 36 = 46656
7682
PhiOnMesh: Number of (b)points on node 0 = 1728
7683
PhiOnMesh: nlist on node 0 = 32391
7685
iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) dHmax(eV)
7686
scf: 1 -6830.882137 -6830.881752 -6830.916410 0.000117 -1.331373 0.002892
7687
scf: 2 -6830.881713 -6830.881743 -6830.916400 0.000188 -1.330517 0.004759
7688
scf: 3 -6830.881761 -6830.881756 -6830.916414 0.000124 -1.330993 0.000154
7690
SCF Convergence by DM+H criterion
7691
max |DM_out - DM_in| : 0.0001237182
7692
max |H_out - H_in| (eV) : 0.0001544132
7693
SCF cycle converged after 3 iterations
7695
Using DM_out to compute the final energy and forces
7696
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7698
siesta: E_KS(eV) = -6830.8818
7700
siesta: Atomic forces (eV/Ang):
7701
----------------------------------------
7702
Tot -0.000000 -0.000000 -0.000000
7703
----------------------------------------
7705
Res 0.003979 sqrt( Sum f_i^2 / 3N )
7706
----------------------------------------
7707
Max 0.010437 constrained
7709
Stress-tensor-Voigt (kbar): -51.12 -51.12 -51.12 -3.05 -3.05 -3.05
7710
(Free)E + p*V (eV/cell) -6790.0476
7711
Target enthalpy (eV/cell) -6830.9164
7713
siesta: Stress tensor (static) (eV/Ang**3):
7714
-0.031908 -0.001906 -0.001906
7715
-0.001906 -0.031908 -0.001906
7716
-0.001906 -0.001906 -0.031908
7718
siesta: Pressure (static): 51.12307775 kBar
7720
siesta: Stress tensor (total) (eV/Ang**3):
7721
-0.031908 -0.001906 -0.001906
7722
-0.001906 -0.031908 -0.001906
7723
-0.001906 -0.001906 -0.031908
7725
siesta: Pressure (total): 51.12318837 kBar
7726
FIRE: istep, dt, firenpos, firealf, magv
7727
FIRE: 64 3.6289 3 0.1000 0.0009
7729
siesta: Temp_ion = 0.025 K
7731
====================================
7732
Begin MD quenched step = 65
7733
====================================
7735
outcell: Unit cell vectors (Ang):
7736
10.860000 0.000000 0.000000
7737
0.000000 10.860000 0.000000
7738
0.000000 0.000000 10.860000
7740
outcell: Cell vector modules (Ang) : 10.860000 10.860000 10.860000
7741
outcell: Cell angles (23,13,12) (deg): 90.0000 90.0000 90.0000
7742
outcell: Cell volume (Ang**3) : 1280.8241
7743
<dSpData1D:S at geom step 65
7744
<sparsity:sparsity for geom step 65
7745
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 7>
7746
<dData1D:(new from dSpData1D) n=5616, refcount: 1>
7749
Re-using DM from previous geometries...
7750
Number of DMs in history: 1
7751
DM extrapolation coefficients:
7753
New DM after history re-use:
7754
<dSpData2D:SpM extrapolated using coords
7755
<sparsity:sparsity for geom step 65
7756
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 9>
7757
<dData2D:(temp array for extrapolation) n=5616 m=1, refcount: 1>
7759
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7760
New grid distribution: 1
7763
3 1: 24 1: 12 13: 18
7764
4 1: 24 1: 12 19: 24
7766
6 1: 24 13: 24 7: 12
7767
7 1: 24 13: 24 13: 18
7768
8 1: 24 13: 24 19: 24
7770
InitMesh: MESH = 48 x 48 x 48 = 110592
7771
InitMesh: (bp) = 24 x 24 x 24 = 13824
7772
InitMesh: Mesh cutoff (required, used) = 40.000 53.991 Ry
7773
ExtMesh (bp) on 0 = 48 x 36 x 30 = 51840
7774
New grid distribution: 2
7775
1 1: 12 13: 24 13: 24
7776
2 13: 24 1: 12 1: 12
7777
3 13: 24 1: 12 13: 24
7778
4 13: 24 13: 24 1: 12
7779
5 1: 12 1: 12 13: 24
7781
7 13: 24 13: 24 13: 24
7782
8 1: 12 13: 24 1: 12
7783
New grid distribution: 3
7784
1 1: 12 13: 24 13: 24
7785
2 13: 24 1: 12 1: 12
7786
3 13: 24 1: 12 13: 24
7787
4 13: 24 13: 24 1: 12
7788
5 1: 12 1: 12 13: 24
7790
7 13: 24 13: 24 13: 24
7791
8 1: 12 13: 24 1: 12
7792
Setting up quadratic distribution...
7793
ExtMesh (bp) on 0 = 36 x 36 x 36 = 46656
7794
PhiOnMesh: Number of (b)points on node 0 = 1728
7795
PhiOnMesh: nlist on node 0 = 32383
7797
iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) dHmax(eV)
7798
scf: 1 -6830.882302 -6830.881791 -6830.916448 0.000148 -1.331291 0.004502
7799
scf: 2 -6830.881742 -6830.881779 -6830.916436 0.000253 -1.329607 0.009533
7800
scf: 3 -6830.881800 -6830.881795 -6830.916452 0.000165 -1.330739 0.000352
7802
SCF Convergence by DM+H criterion
7803
max |DM_out - DM_in| : 0.0001649630
7804
max |H_out - H_in| (eV) : 0.0003519680
7805
SCF cycle converged after 3 iterations
7807
Using DM_out to compute the final energy and forces
7808
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7810
siesta: E_KS(eV) = -6830.8818
7812
siesta: Atomic forces (eV/Ang):
7813
----------------------------------------
7814
Tot 0.000000 -0.000000 -0.000000
7815
----------------------------------------
7817
Res 0.005511 sqrt( Sum f_i^2 / 3N )
7818
----------------------------------------
7819
Max 0.015514 constrained
7821
Stress-tensor-Voigt (kbar): -50.88 -50.88 -50.88 -2.98 -2.98 -2.98
7822
(Free)E + p*V (eV/cell) -6790.2440
7823
Target enthalpy (eV/cell) -6830.9165
7825
siesta: Stress tensor (static) (eV/Ang**3):
7826
-0.031755 -0.001857 -0.001857
7827
-0.001857 -0.031755 -0.001857
7828
-0.001857 -0.001857 -0.031755
7830
siesta: Pressure (static): 50.87739184 kBar
7832
siesta: Stress tensor (total) (eV/Ang**3):
7833
-0.031755 -0.001857 -0.001857
7834
-0.001857 -0.031755 -0.001857
7835
-0.001857 -0.001857 -0.031755
7837
siesta: Pressure (total): 50.87756589 kBar
7838
FIRE: istep, dt, firenpos, firealf, magv
7839
FIRE: 65 3.6289 4 0.1000 0.0009
7841
siesta: Temp_ion = 0.026 K
7843
====================================
7844
Begin MD quenched step = 66
7845
====================================
7847
outcell: Unit cell vectors (Ang):
7848
10.860000 0.000000 0.000000
7849
0.000000 10.860000 0.000000
7850
0.000000 0.000000 10.860000
7852
outcell: Cell vector modules (Ang) : 10.860000 10.860000 10.860000
7853
outcell: Cell angles (23,13,12) (deg): 90.0000 90.0000 90.0000
7854
outcell: Cell volume (Ang**3) : 1280.8241
7855
<dSpData1D:S at geom step 66
7856
<sparsity:sparsity for geom step 66
7857
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 7>
7858
<dData1D:(new from dSpData1D) n=5616, refcount: 1>
7861
Re-using DM from previous geometries...
7862
Number of DMs in history: 1
7863
DM extrapolation coefficients:
7865
New DM after history re-use:
7866
<dSpData2D:SpM extrapolated using coords
7867
<sparsity:sparsity for geom step 66
7868
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 9>
7869
<dData2D:(temp array for extrapolation) n=5616 m=1, refcount: 1>
7871
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7872
New grid distribution: 1
7875
3 1: 24 1: 12 13: 18
7876
4 1: 24 1: 12 19: 24
7878
6 1: 24 13: 24 7: 12
7879
7 1: 24 13: 24 13: 18
7880
8 1: 24 13: 24 19: 24
7882
InitMesh: MESH = 48 x 48 x 48 = 110592
7883
InitMesh: (bp) = 24 x 24 x 24 = 13824
7884
InitMesh: Mesh cutoff (required, used) = 40.000 53.991 Ry
7885
ExtMesh (bp) on 0 = 48 x 36 x 30 = 51840
7886
New grid distribution: 2
7887
1 1: 12 13: 24 13: 24
7888
2 13: 24 1: 12 1: 12
7889
3 13: 24 1: 12 13: 24
7890
4 13: 24 13: 24 1: 12
7891
5 1: 12 1: 12 13: 24
7893
7 13: 24 13: 24 13: 24
7894
8 1: 12 13: 24 1: 12
7895
New grid distribution: 3
7896
1 1: 12 13: 24 13: 24
7897
2 13: 24 1: 12 1: 12
7898
3 13: 24 1: 12 13: 24
7899
4 13: 24 13: 24 1: 12
7900
5 1: 12 1: 12 13: 24
7902
7 13: 24 13: 24 13: 24
7903
8 1: 12 13: 24 1: 12
7904
Setting up quadratic distribution...
7905
ExtMesh (bp) on 0 = 36 x 36 x 36 = 46656
7906
PhiOnMesh: Number of (b)points on node 0 = 1728
7907
PhiOnMesh: nlist on node 0 = 32383
7909
iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) dHmax(eV)
7910
scf: 1 -6830.882477 -6830.881838 -6830.916495 0.000159 -1.331294 0.006008
7911
scf: 2 -6830.881749 -6830.881814 -6830.916471 0.000311 -1.328581 0.015193
7912
scf: 3 -6830.881851 -6830.881844 -6830.916501 0.000222 -1.330537 0.000513
7914
SCF Convergence by DM+H criterion
7915
max |DM_out - DM_in| : 0.0002215329
7916
max |H_out - H_in| (eV) : 0.0005132950
7917
SCF cycle converged after 3 iterations
7919
Using DM_out to compute the final energy and forces
7920
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7922
siesta: E_KS(eV) = -6830.8818
7095
7924
siesta: Atomic forces (eV/Ang):
7096
7925
----------------------------------------
7097
7926
Tot 0.000000 -0.000000 0.000000
7098
7927
----------------------------------------
7100
Res 0.004009 sqrt( Sum f_i^2 / 3N )
7101
----------------------------------------
7102
Max 0.009790 constrained
7104
Stress-tensor-Voigt (kbar): -49.70 -49.70 -49.70 -3.00 -3.00 -3.00
7105
(Free)E + p*V (eV/cell) -6791.1421
7106
Target enthalpy (eV/cell) -6830.8713
7108
siesta: Stress tensor (static) (eV/Ang**3):
7109
-0.031018 -0.001872 -0.001872
7110
-0.001872 -0.031018 -0.001872
7111
-0.001872 -0.001872 -0.031018
7113
siesta: Pressure (static): 49.69761431 kBar
7115
siesta: Stress tensor (total) (eV/Ang**3):
7116
-0.031018 -0.001872 -0.001872
7117
-0.001872 -0.031018 -0.001872
7118
-0.001872 -0.001872 -0.031018
7120
siesta: Pressure (total): 49.69764862 kBar
7929
Res 0.005177 sqrt( Sum f_i^2 / 3N )
7930
----------------------------------------
7931
Max 0.015353 constrained
7933
Stress-tensor-Voigt (kbar): -50.87 -50.87 -50.87 -2.98 -2.98 -2.98
7934
(Free)E + p*V (eV/cell) -6790.2526
7935
Target enthalpy (eV/cell) -6830.9165
7937
siesta: Stress tensor (static) (eV/Ang**3):
7938
-0.031748 -0.001858 -0.001858
7939
-0.001858 -0.031748 -0.001858
7940
-0.001858 -0.001858 -0.031748
7942
siesta: Pressure (static): 50.86660649 kBar
7944
siesta: Stress tensor (total) (eV/Ang**3):
7945
-0.031748 -0.001858 -0.001858
7946
-0.001858 -0.031748 -0.001858
7947
-0.001858 -0.001858 -0.031748
7949
siesta: Pressure (total): 50.86678674 kBar
7950
FIRE: istep, dt, firenpos, firealf, magv
7951
FIRE: 66 3.6289 5 0.1000 0.0010
7953
siesta: Temp_ion = 0.029 K
7955
====================================
7956
Begin MD quenched step = 67
7957
====================================
7959
outcell: Unit cell vectors (Ang):
7960
10.860000 0.000000 0.000000
7961
0.000000 10.860000 0.000000
7962
0.000000 0.000000 10.860000
7964
outcell: Cell vector modules (Ang) : 10.860000 10.860000 10.860000
7965
outcell: Cell angles (23,13,12) (deg): 90.0000 90.0000 90.0000
7966
outcell: Cell volume (Ang**3) : 1280.8241
7967
<dSpData1D:S at geom step 67
7968
<sparsity:sparsity for geom step 67
7969
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 7>
7970
<dData1D:(new from dSpData1D) n=5616, refcount: 1>
7973
Re-using DM from previous geometries...
7974
Number of DMs in history: 1
7975
DM extrapolation coefficients:
7977
New DM after history re-use:
7978
<dSpData2D:SpM extrapolated using coords
7979
<sparsity:sparsity for geom step 67
7980
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 9>
7981
<dData2D:(temp array for extrapolation) n=5616 m=1, refcount: 1>
7983
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
7984
New grid distribution: 1
7987
3 1: 24 1: 12 13: 18
7988
4 1: 24 1: 12 19: 24
7990
6 1: 24 13: 24 7: 12
7991
7 1: 24 13: 24 13: 18
7992
8 1: 24 13: 24 19: 24
7994
InitMesh: MESH = 48 x 48 x 48 = 110592
7995
InitMesh: (bp) = 24 x 24 x 24 = 13824
7996
InitMesh: Mesh cutoff (required, used) = 40.000 53.991 Ry
7997
ExtMesh (bp) on 0 = 48 x 36 x 30 = 51840
7998
New grid distribution: 2
7999
1 1: 12 13: 24 13: 24
8000
2 13: 24 1: 12 1: 12
8001
3 13: 24 1: 12 13: 24
8002
4 13: 24 13: 24 1: 12
8003
5 1: 12 1: 12 13: 24
8005
7 13: 24 13: 24 13: 24
8006
8 1: 12 13: 24 1: 12
8007
New grid distribution: 3
8008
1 1: 12 13: 24 13: 24
8009
2 13: 24 1: 12 1: 12
8010
3 13: 24 1: 12 13: 24
8011
4 13: 24 13: 24 1: 12
8012
5 1: 12 1: 12 13: 24
8014
7 13: 24 13: 24 13: 24
8015
8 1: 12 13: 24 1: 12
8016
Setting up quadratic distribution...
8017
ExtMesh (bp) on 0 = 36 x 36 x 36 = 46656
8018
PhiOnMesh: Number of (b)points on node 0 = 1728
8019
PhiOnMesh: nlist on node 0 = 32391
8021
iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) dHmax(eV)
8022
scf: 1 -6830.882605 -6830.881902 -6830.916559 0.000184 -1.331340 0.007290
8023
scf: 2 -6830.881736 -6830.881856 -6830.916513 0.000352 -1.327781 0.020667
8024
scf: 3 -6830.881923 -6830.881910 -6830.916567 0.000265 -1.330431 0.000494
8026
SCF Convergence by DM+H criterion
8027
max |DM_out - DM_in| : 0.0002646385
8028
max |H_out - H_in| (eV) : 0.0004939356
8029
SCF cycle converged after 3 iterations
8031
Using DM_out to compute the final energy and forces
8032
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
8034
siesta: E_KS(eV) = -6830.8819
8036
siesta: Atomic forces (eV/Ang):
8037
----------------------------------------
8038
Tot -0.000000 -0.000000 0.000000
8039
----------------------------------------
8041
Res 0.004022 sqrt( Sum f_i^2 / 3N )
8042
----------------------------------------
8043
Max 0.012163 constrained
8045
Stress-tensor-Voigt (kbar): -50.87 -50.87 -50.87 -3.00 -3.00 -3.00
8046
(Free)E + p*V (eV/cell) -6790.2523
8047
Target enthalpy (eV/cell) -6830.9166
8049
siesta: Stress tensor (static) (eV/Ang**3):
8050
-0.031748 -0.001873 -0.001873
8051
-0.001873 -0.031748 -0.001873
8052
-0.001873 -0.001873 -0.031748
8054
siesta: Pressure (static): 50.86714726 kBar
8056
siesta: Stress tensor (total) (eV/Ang**3):
8057
-0.031749 -0.001873 -0.001873
8058
-0.001873 -0.031749 -0.001873
8059
-0.001873 -0.001873 -0.031749
8061
siesta: Pressure (total): 50.86734413 kBar
8062
FIRE: istep, dt, firenpos, firealf, magv
8063
FIRE: 67 3.6289 6 0.1000 0.0010
8065
siesta: Temp_ion = 0.034 K
8067
====================================
8068
Begin MD quenched step = 68
8069
====================================
8071
outcell: Unit cell vectors (Ang):
8072
10.860000 0.000000 0.000000
8073
0.000000 10.860000 0.000000
8074
0.000000 0.000000 10.860000
8076
outcell: Cell vector modules (Ang) : 10.860000 10.860000 10.860000
8077
outcell: Cell angles (23,13,12) (deg): 90.0000 90.0000 90.0000
8078
outcell: Cell volume (Ang**3) : 1280.8241
8079
<dSpData1D:S at geom step 68
8080
<sparsity:sparsity for geom step 68
8081
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 7>
8082
<dData1D:(new from dSpData1D) n=5616, refcount: 1>
8085
Re-using DM from previous geometries...
8086
Number of DMs in history: 1
8087
DM extrapolation coefficients:
8089
New DM after history re-use:
8090
<dSpData2D:SpM extrapolated using coords
8091
<sparsity:sparsity for geom step 68
8092
nrows_g=257 nrows=48 sparsity=.0850 nnzs=5616, refcount: 9>
8093
<dData2D:(temp array for extrapolation) n=5616 m=1, refcount: 1>
8095
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
8096
New grid distribution: 1
8099
3 1: 24 1: 12 13: 18
8100
4 1: 24 1: 12 19: 24
8102
6 1: 24 13: 24 7: 12
8103
7 1: 24 13: 24 13: 18
8104
8 1: 24 13: 24 19: 24
8106
InitMesh: MESH = 48 x 48 x 48 = 110592
8107
InitMesh: (bp) = 24 x 24 x 24 = 13824
8108
InitMesh: Mesh cutoff (required, used) = 40.000 53.991 Ry
8109
ExtMesh (bp) on 0 = 48 x 36 x 30 = 51840
8110
New grid distribution: 2
8111
1 1: 12 13: 24 13: 24
8112
2 13: 24 1: 12 1: 12
8113
3 13: 24 1: 12 13: 24
8114
4 13: 24 13: 24 1: 12
8115
5 1: 12 1: 12 13: 24
8117
7 13: 24 13: 24 13: 24
8118
8 1: 12 13: 24 1: 12
8119
New grid distribution: 3
8120
1 1: 12 13: 24 13: 24
8121
2 13: 24 1: 12 1: 12
8122
3 13: 24 1: 12 13: 24
8123
4 13: 24 13: 24 1: 12
8124
5 1: 12 1: 12 13: 24
8126
7 13: 24 13: 24 13: 24
8127
8 1: 12 13: 24 1: 12
8128
Setting up quadratic distribution...
8129
ExtMesh (bp) on 0 = 36 x 36 x 36 = 46656
8130
PhiOnMesh: Number of (b)points on node 0 = 1728
8131
PhiOnMesh: nlist on node 0 = 32391
8133
iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) dHmax(eV)
8134
scf: 1 -6830.882650 -6830.881955 -6830.916612 0.000237 -1.331489 0.008513
8135
scf: 2 -6830.881671 -6830.881874 -6830.916532 0.000393 -1.327143 0.026227
8136
scf: 3 -6830.881987 -6830.881967 -6830.916624 0.000307 -1.330449 0.000422
8138
SCF Convergence by DM+H criterion
8139
max |DM_out - DM_in| : 0.0003074656
8140
max |H_out - H_in| (eV) : 0.0004222841
8141
SCF cycle converged after 3 iterations
8143
Using DM_out to compute the final energy and forces
8144
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 34 24
8146
siesta: E_KS(eV) = -6830.8820
8148
siesta: Atomic forces (eV/Ang):
8149
----------------------------------------
8150
Tot -0.000000 -0.000000 -0.000000
8151
----------------------------------------
8153
Res 0.002678 sqrt( Sum f_i^2 / 3N )
8154
----------------------------------------
8155
Max 0.006870 constrained
8157
Stress-tensor-Voigt (kbar): -50.88 -50.88 -50.88 -3.04 -3.04 -3.04
8158
(Free)E + p*V (eV/cell) -6790.2461
8159
Target enthalpy (eV/cell) -6830.9166
8161
siesta: Stress tensor (static) (eV/Ang**3):
8162
-0.031753 -0.001895 -0.001895
8163
-0.001895 -0.031753 -0.001895
8164
-0.001895 -0.001895 -0.031753
8166
siesta: Pressure (static): 50.87492559 kBar
8168
siesta: Stress tensor (total) (eV/Ang**3):
8169
-0.031753 -0.001895 -0.001895
8170
-0.001895 -0.031753 -0.001895
8171
-0.001895 -0.001895 -0.031753
8173
siesta: Pressure (total): 50.87515925 kBar
7122
8175
outcoor: Relaxed atomic coordinates (scaled):
7123
-0.00000485 -0.00000485 -0.00000485 1 1 Si
7124
0.25000485 0.25000485 0.25000485 1 2 Si
7125
-0.00174822 0.49947078 0.49947078 1 3 Si
7126
0.24914622 0.74911820 0.74911820 1 4 Si
7127
0.49947078 -0.00174822 0.49947078 1 5 Si
7128
0.74911820 0.24914622 0.74911820 1 6 Si
7129
0.49947078 0.49947078 -0.00174822 1 7 Si
7130
0.74911820 0.74911820 0.24914622 1 8 Si
7131
0.99965940 0.00169873 0.00169873 1 9 Si
7132
1.25034060 0.24830127 0.24830127 1 10 Si
7133
0.99847216 0.49523600 0.49523600 1 11 Si
7134
1.25256849 0.73904643 0.73904643 1 12 Si
7135
1.49897585 0.00003695 0.50070940 1 13 Si
7136
1.74929060 0.24996305 0.75102415 1 14 Si
7137
1.49897585 0.50070940 0.00003695 1 15 Si
7138
1.74929060 0.75102415 0.24996305 1 16 Si
7139
0.00169873 0.99965940 0.00169872 1 17 Si
7140
0.24830127 1.25034060 0.24830127 1 18 Si
7141
0.00003695 1.49897585 0.50070940 1 19 Si
7142
0.24996305 1.74929060 0.75102415 1 20 Si
7143
0.49523600 0.99847216 0.49523600 1 21 Si
7144
0.73904643 1.25256849 0.73904643 1 22 Si
7145
0.50070940 1.49897585 0.00003695 1 23 Si
7146
0.75102415 1.74929060 0.24996305 1 24 Si
7147
0.00169873 0.00169873 0.99965940 1 25 Si
7148
0.24830127 0.24830127 1.25034060 1 26 Si
7149
0.00003695 0.50070940 1.49897585 1 27 Si
7150
0.24996305 0.75102415 1.74929060 1 28 Si
7151
0.50070940 0.00003695 1.49897585 1 29 Si
7152
0.75102415 0.24996305 1.74929060 1 30 Si
7153
0.49523600 0.49523600 0.99847216 1 31 Si
7154
0.73904643 0.73904643 1.25256849 1 32 Si
7155
0.99984001 0.99984001 0.00085257 1 33 Si
7156
1.25015999 1.25015999 0.24914743 1 34 Si
7157
1.00146700 1.49930501 0.49957947 1 35 Si
7158
1.24853300 1.75042053 0.75069500 1 36 Si
7159
1.49930501 1.00146700 0.49957947 1 37 Si
7160
1.75042053 1.24853300 0.75069500 1 38 Si
7161
1.50088180 1.50088180 0.00085378 1 39 Si
7162
1.75052922 1.75052922 0.25174822 1 40 Si
7163
0.99984001 0.00085257 0.99984001 1 41 Si
7164
1.25015999 0.24914743 1.25015999 1 42 Si
7165
1.00146700 0.49957947 1.49930501 1 43 Si
7166
1.24853300 0.75069500 1.75042053 1 44 Si
7167
1.50088180 0.00085378 1.50088180 1 45 Si
7168
1.75052922 0.25174822 1.75052922 1 46 Si
7169
1.49930501 0.49957947 1.00146700 1 47 Si
7170
1.75042053 0.75069500 1.24853300 1 48 Si
7171
0.00085257 0.99984001 0.99984001 1 49 Si
7172
0.24914743 1.25015999 1.25015999 1 50 Si
7173
0.00085378 1.50088180 1.50088180 1 51 Si
7174
0.25174822 1.75052922 1.75052922 1 52 Si
7175
0.49957947 1.00146700 1.49930501 1 53 Si
7176
0.75069500 1.24853300 1.75042053 1 54 Si
7177
0.49957947 1.49930501 1.00146700 1 55 Si
7178
0.75069500 1.75042053 1.24853300 1 56 Si
7179
0.94954240 0.94954240 0.94954240 1 57 Si
7180
1.30045760 1.30045760 1.30045760 1 58 Si
7181
0.99743151 1.51095357 1.51095357 1 59 Si
7182
1.25152784 1.75476400 1.75476400 1 60 Si
7183
1.51095357 0.99743151 1.51095357 1 61 Si
7184
1.75476400 1.25152784 1.75476400 1 62 Si
7185
1.51095357 1.51095357 0.99743151 1 63 Si
7186
1.75476400 1.75476400 1.25152784 1 64 Si
8176
0.00003134 0.00003134 0.00003134 1 1 Si
8177
0.24996866 0.24996866 0.24996866 1 2 Si
8178
-0.00170142 0.49955598 0.49955598 1 3 Si
8179
0.24964404 0.74890222 0.74890222 1 4 Si
8180
0.49955598 -0.00170142 0.49955598 1 5 Si
8181
0.74890222 0.24964404 0.74890222 1 6 Si
8182
0.49955598 0.49955598 -0.00170142 1 7 Si
8183
0.74890222 0.74890222 0.24964404 1 8 Si
8184
0.99988295 0.00178259 0.00178259 1 9 Si
8185
1.25011705 0.24821741 0.24821741 1 10 Si
8186
0.99885097 0.49480518 0.49480518 1 11 Si
8187
1.25182688 0.73889755 0.73889755 1 12 Si
8188
1.49892747 0.00027352 0.50047594 1 13 Si
8189
1.74952406 0.24972648 0.75107253 1 14 Si
8190
1.49892747 0.50047594 0.00027352 1 15 Si
8191
1.74952406 0.75107253 0.24972648 1 16 Si
8192
0.00178259 0.99988295 0.00178259 1 17 Si
8193
0.24821741 1.25011705 0.24821741 1 18 Si
8194
0.00027352 1.49892747 0.50047594 1 19 Si
8195
0.24972648 1.74952406 0.75107253 1 20 Si
8196
0.49480518 0.99885096 0.49480518 1 21 Si
8197
0.73889755 1.25182688 0.73889755 1 22 Si
8198
0.50047594 1.49892747 0.00027352 1 23 Si
8199
0.75107253 1.74952406 0.24972648 1 24 Si
8200
0.00178259 0.00178259 0.99988295 1 25 Si
8201
0.24821741 0.24821741 1.25011705 1 26 Si
8202
0.00027352 0.50047594 1.49892747 1 27 Si
8203
0.24972648 0.75107253 1.74952406 1 28 Si
8204
0.50047594 0.00027352 1.49892747 1 29 Si
8205
0.75107253 0.24972648 1.74952406 1 30 Si
8206
0.49480518 0.49480518 0.99885097 1 31 Si
8207
0.73889755 0.73889755 1.25182688 1 32 Si
8208
0.99960595 0.99960595 0.00111246 1 33 Si
8209
1.25039405 1.25039405 0.24888754 1 34 Si
8210
1.00127997 1.49919798 0.49955145 1 35 Si
8211
1.24872003 1.75044854 0.75080203 1 36 Si
8212
1.49919797 1.00127997 0.49955146 1 37 Si
8213
1.75044854 1.24872003 0.75080203 1 38 Si
8214
1.50109778 1.50109778 0.00035596 1 39 Si
8215
1.75044402 1.75044402 0.25170142 1 40 Si
8216
0.99960595 0.00111246 0.99960595 1 41 Si
8217
1.25039405 0.24888754 1.25039405 1 42 Si
8218
1.00127997 0.49955146 1.49919797 1 43 Si
8219
1.24872003 0.75080203 1.75044855 1 44 Si
8220
1.50109778 0.00035596 1.50109778 1 45 Si
8221
1.75044402 0.25170142 1.75044402 1 46 Si
8222
1.49919797 0.49955146 1.00127997 1 47 Si
8223
1.75044855 0.75080203 1.24872003 1 48 Si
8224
0.00111246 0.99960595 0.99960595 1 49 Si
8225
0.24888754 1.25039405 1.25039405 1 50 Si
8226
0.00035596 1.50109778 1.50109778 1 51 Si
8227
0.25170142 1.75044402 1.75044402 1 52 Si
8228
0.49955146 1.00127997 1.49919797 1 53 Si
8229
0.75080203 1.24872003 1.75044854 1 54 Si
8230
0.49955145 1.49919798 1.00127997 1 55 Si
8231
0.75080203 1.75044854 1.24872003 1 56 Si
8232
0.95050167 0.95050168 0.95050167 1 57 Si
8233
1.29949832 1.29949832 1.29949832 1 58 Si
8234
0.99817312 1.51110245 1.51110245 1 59 Si
8235
1.25114903 1.75519482 1.75519482 1 60 Si
8236
1.51110245 0.99817312 1.51110245 1 61 Si
8237
1.75519482 1.25114904 1.75519482 1 62 Si
8238
1.51110245 1.51110245 0.99817312 1 63 Si
8239
1.75519482 1.75519482 1.25114903 1 64 Si
7187
8240
1.12500000 1.12500000 1.12500000 2 65 H
7189
8242
siesta: Program's energy decomposition (eV):
7190
siesta: Ebs = -1985.324036
8243
siesta: Ebs = -2019.264087
7191
8244
siesta: Eions = 12205.375272
7192
siesta: Ena = 3580.587855
7193
siesta: Ekin = 2732.770470
7194
siesta: Enl = 958.324401
8245
siesta: Ena = 3580.580177
8246
siesta: Ekin = 2733.209032
8247
siesta: Enl = 960.447518
7195
8248
siesta: Eso = 0.000000
7196
8249
siesta: Eldau = 0.000000
7197
siesta: DEna = 197.450969
7198
siesta: DUscf = 5.558051
8250
siesta: DEna = 195.273286
8251
siesta: DUscf = 5.414802
7199
8252
siesta: DUext = 0.000000
7200
siesta: Exc = -2100.153129
8253
siesta: Enegf = 0.000000
8254
siesta: Exc = -2100.431510
7201
8255
siesta: eta*DQ = 0.000000
7202
8256
siesta: Emadel = 0.000000
7203
8257
siesta: Emeta = 0.000000
7204
8258
siesta: Emolmec = 0.000000
7205
siesta: Ekinion = 0.000041
7206
siesta: Eharris = -6832.124694
7207
siesta: Etot = -6830.836613
7208
siesta: FreeEng = -6830.871270
8259
siesta: Ekinion = 0.000282
8260
siesta: Eharris = -6830.881705
8261
siesta: Etot = -6830.881684
8262
siesta: FreeEng = -6830.916342
7210
8264
siesta: Final energy (eV):
7211
siesta: Band Struct. = -1985.324036
7212
siesta: Kinetic = 2732.770470
7213
siesta: Hartree = 493.162455
8265
siesta: Band Struct. = -2019.264087
8266
siesta: Kinetic = 2733.209032
8267
siesta: Hartree = 493.999666
7214
8268
siesta: Eldau = 0.000000
7215
8269
siesta: Eso = 0.000000
7216
8270
siesta: Ext. field = 0.000000
7217
siesta: Exch.-corr. = -2100.153129
7218
siesta: Ion-electron = -3000.435874
7219
siesta: Ion-ion = -4956.180577
7220
siesta: Ekinion = 0.000041
7221
siesta: Total = -6830.836613
7222
siesta: Fermi = -1.286537
8271
siesta: Enegf = 0.000000
8272
siesta: Exch.-corr. = -2100.431510
8273
siesta: Ion-electron = -3001.657749
8274
siesta: Ion-ion = -4956.001406
8275
siesta: Ekinion = 0.000282
8276
siesta: Total = -6830.881684
8277
siesta: Fermi = -1.330449
7224
8279
siesta: Stress tensor (static) (eV/Ang**3):
7225
siesta: -0.031018 -0.001872 -0.001872
7226
siesta: -0.001872 -0.031018 -0.001872
7227
siesta: -0.001872 -0.001872 -0.031018
8280
siesta: -0.031753 -0.001895 -0.001895
8281
siesta: -0.001895 -0.031753 -0.001895
8282
siesta: -0.001895 -0.001895 -0.031753
7229
8284
siesta: Cell volume = 1280.824056 Ang**3
7231
8286
siesta: Pressure (static):
7232
8287
siesta: Solid Molecule Units
7233
siesta: 0.00033783 0.00033874 Ry/Bohr**3
7234
siesta: 0.03101848 0.03110228 eV/Ang**3
7235
siesta: 49.69761431 49.83188278 kBar
7236
(Free)E+ p_basis*V_orbitals = -6818.317332
7237
(Free)Eharris+ p_basis*V_orbitals = -6819.605413
7238
>> End of run: 2-JUL-2017 12:01:12
8288
siesta: 0.00034583 0.00034474 Ry/Bohr**3
8289
siesta: 0.03175329 0.03165245 eV/Ang**3
8290
siesta: 50.87492559 50.71336063 kBar
8291
(Free)E+ p_basis*V_orbitals = -6818.362644
8292
(Free)Eharris+ p_basis*V_orbitals = -6818.362665
8293
>> End of run: 5-NOV-2018 11:14:04