10
10
reinit: -----------------------------------------------------------------------
11
11
reinit: System Label: md_pr
12
12
reinit: -----------------------------------------------------------------------
13
Siesta Version: trunk-554
14
Architecture : x86_64-linux-gcc
15
Compiler flags: mpif90 -m64 -fPIC -O3 -ftree-vectorize -fexpensive-optimizations -funroll-loops -fprefetch-loop-arrays -fgraphite -fipa-icf -fipa-pure-const -march=native -fschedule-fusion -fselective-scheduling -fipa-sra -fipa-cp -fno-second-underscore -I/opt/zlib/1.2.8/gnu-6.1.0/include -I/opt/hdf5/1.8.16/gnu-6.1.0/include -I/opt/pnetcdf/1.7.0/gnu-6.1.0/include -I/opt/netcdf/4.4.0/gnu-6.1.0/include -I/opt/openmpi/1.10.2/gnu-6.1.0/include -I/opt/mumps/5.0.1/gnu-6.1.0/include
16
PP flags : -DSIESTA__FLOOK -DSIESTA__METIS -DMPI -DFC_HAVE_FLUSH -DFC_HAVE_ABORT -DCDF -DGRID_DP -DMPI_TIMING -DTRANSIESTA_TIMING -DTBTRANS_TIMING -DSIESTA__MUMPS -DNCDF_PARALLEL -DNCDF -DNCDF_4 -DSCALAPACK_DEBUG
17
Libraries : -lzmumps -lmumps_common -lpord -lparmetis -lmetis -L/opt/scalapack/204/gnu-6.1.0/lib -Wl,-rpath=/opt/scalapack/204/gnu-6.1.0/lib -lscalapack -L/opt/openblas/0.2.17/gnu-6.1.0/lib -Wl,-rpath=/opt/openblas/0.2.17/gnu-6.1.0/lib -lopenblas -L/home/nicpa/phd/esl/flook -lflookall -ldl -lnetcdff -lnetcdf -lpnetcdf -lhdf5_hl -lhdf5 -lz -m64 -fPIC -O3 -ftree-vectorize -fexpensive-optimizations -funroll-loops -fprefetch-loop-arrays -fgraphite -fipa-icf -fipa-pure-const -march=native -fschedule-fusion -fselective-scheduling -fipa-sra -fipa-cp -fno-second-underscore -I/opt/zlib/1.2.8/gnu-6.1.0/include -I/opt/hdf5/1.8.16/gnu-6.1.0/include -I/opt/pnetcdf/1.7.0/gnu-6.1.0/include -I/opt/netcdf/4.4.0/gnu-6.1.0/include -I/opt/openmpi/1.10.2/gnu-6.1.0/include -I/opt/mumps/5.0.1/gnu-6.1.0/include
13
Siesta Version: siesta-4.1--731
14
Architecture : x86_64-linux-n-62-18-14
15
Compiler flags: mpifort -m64 -fPIC -O3 -march=native -ftree-vectorize -fexpensive-optimizatioons -funroll-loops -fprefetch-loop-arrays -fno-second-underscore -flto -fuse-linker-plugin
16
PP flags : -I/zdata/groups/common/nicpa/2016-oct/generic/build-tools/1.0/include -I/zdatta/groups/common/nicpa/2016-oct/generic/gmp/6.1.0/include -I/zdata/groups/common/nicpa/2016-oct/generic/mpfr/3.1.4/include -I/zdata/groups/common/nicpa/2016-oct/generic/mpc/1.0.3/include -I/zdata/groups/common/nicpa/2016-oct/generic/isl/0.16.1/include -I/zdata/groups/common/nicpa/2016-oct/generic/gcc/6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/generic/numactl/2.0.11/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/zlib/1.2.8/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/libxml2/2.9.4/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/hwloc/1.11.4/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/openmpi/2.0.1/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/hdf5/1.8.17/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/pnetcdf/1.7.0/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/netcdf/4.4.1/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/flook/0.7.0/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/parmetis/4.0.3/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/scalapack/204/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/openblas/0.2.19/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/mumps/5.0.2/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/fftw/3.3.5/gnu-6.2.0/include -DMPI -DFC_HAVE_FLUSH -DFC_HAVE_ABORT -DCDF -DCDF4 -DSIESTA__METIS -DSIESTA__MUMPS -DTS_NOCHECKS -DNCDF -DNCDF_4 -DNCDF_PARALLEL -DSIESTA__FLOOK
17
Libraries : -lnetcdff -lnetcdf -lpnetcdf -lhdf5_hl -lhdf5 -lz -lzmumps -lmumps_common -llpord -lparmetis -lmetis -lscalapack -lopenblas -lmetis -flto -fuse-linker-plugin -L/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/flook/0.7.0/gnu-6.2.0/lib -Wl,-rpath=/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/flook/0.7.0/gnu-6.2.0/lib -lflookall -ldl
21
21
NetCDF-4 MPI-IO support
22
22
METIS ordering support
24
* Running on 4 nodes in parallel
25
>> Start of run: 18-AUG-2016 8:43:55
24
* Running on 8 nodes in parallel
25
>> Start of run: 2-JUL-2017 11:41:33
27
27
initatom: Reading input for the pseudopotentials and atomic orbitals ----------
28
28
Species number: 1 Atomic number: 12 Label: Mg
525
522
Folding of H and S is implicitly performed
526
523
<dSpData1D:S at geom step 1
527
524
<sparsity:sparsity for geom step 1
528
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
529
<dData1D:(new from dSpData1D) n=306, refcount: 1>
525
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
526
<dData1D:(new from dSpData1D) n=204, refcount: 1>
532
529
Initializing Density Matrix...
533
DM after filling with atomic data:
534
<dSpData2D:(DM initialized from atoms)
530
DM filled with atomic data:
531
<dSpData2D:DM initialized from atoms
535
532
<sparsity:sparsity for geom step 1
536
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
537
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
533
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
534
<dData2D:DM n=204 m=1, refcount: 1>
539
537
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.58964258
538
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
540
539
New grid distribution: 1
546
549
InitMesh: MESH = 24 x 24 x 24 = 13824
547
550
InitMesh: (bp) = 12 x 12 x 12 = 1728
548
551
InitMesh: Mesh cutoff (required, used) = 50.000 63.482 Ry
549
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
552
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
550
553
New grid distribution: 2
555
562
New grid distribution: 3
560
571
Setting up quadratic distribution...
561
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
562
PhiOnMesh: Number of (b)points on node 0 = 432
563
PhiOnMesh: nlist on node 0 = 9991
564
rdiag-debug: Node=1, lwork= 5744>= lworkq= 5742, liwork= 204>= liworkq= 204
565
rdiag-debug: Node=2, lwork= 5744>= lworkq= 5742, liwork= 204>= liworkq= 204
566
rdiag-debug: Node=3, lwork= 5744>= lworkq= 5742, liwork= 204>= liworkq= 204
567
rdiag-debug: Node=0, lwork= 5744>= lworkq= 5742, liwork= 204>= liworkq= 204
572
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
573
PhiOnMesh: Number of (b)points on node 0 = 252
574
PhiOnMesh: nlist on node 0 = 5554
569
576
stepf: Fermi-Dirac step function
571
578
siesta: Program's energy decomposition (eV):
572
siesta: Ebs = -718.671412
579
siesta: Ebs = -718.671408
573
580
siesta: Eions = 5228.674091
574
581
siesta: Ena = 1092.429444
575
582
siesta: Ekin = 2069.792750
576
583
siesta: Enl = -254.036133
577
584
siesta: Eso = 0.000000
578
585
siesta: Eldau = 0.000000
579
siesta: DEna = 15.961140
586
siesta: DEna = 15.961139
580
587
siesta: DUscf = 0.158223
581
588
siesta: DUext = 0.000000
582
589
siesta: Exc = -679.402766
585
592
siesta: Emeta = 0.000000
586
593
siesta: Emolmec = 0.000000
587
594
siesta: Ekinion = 0.000000
588
siesta: Eharris = -2966.596232
595
siesta: Eharris = -2966.596231
589
596
siesta: Etot = -2983.771434
590
597
siesta: FreeEng = -2983.771434
592
siesta: Eharris(eV) = -2966.596232
595
siesta: Eharris(eV) = -2966.596232
597
timer: Routine,Calls,Time,% = IterSCF 1 0.020 2.49
599
siesta: Eharris(eV) = -2966.596231
602
siesta: Eharris(eV) = -2966.596231
604
timer: Routine,Calls,Time,% = IterSCF 1 0.015 2.53
605
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
606
Geom step, scf iteration, dmax: 1 1 1.795009
599
608
Using DM_out to compute the final energy and forces
609
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
601
611
siesta: E_KS - E_eggbox = -2945.6640
603
613
siesta: Atomic forces (eV/Ang):
604
614
1 0.000000 -0.000000 0.000000
605
615
2 -0.000000 0.000000 0.000000
606
3 -0.000000 0.000000 0.000000
607
4 -0.000000 -0.000000 -0.000000
608
5 2.998821 -6.707247 0.000000
609
6 -0.000000 4.024012 -6.147146
610
7 -2.998821 2.683234 6.147146
611
8 -2.998821 6.707247 0.000000
612
9 -0.000000 -4.024012 6.147146
613
10 2.998821 -2.683234 -6.147146
616
3 -0.000001 0.000001 -0.000001
617
4 0.000001 -0.000001 0.000001
618
5 2.998821 -6.707248 -0.000000
619
6 0.000000 4.024012 -6.147148
620
7 -2.998820 2.683234 6.147148
621
8 -2.998821 6.707248 0.000000
622
9 -0.000000 -4.024012 6.147148
623
10 2.998820 -2.683234 -6.147148
614
624
----------------------------------------
615
625
Tot -0.000000 -0.000000 -0.000000
616
626
----------------------------------------
618
Res 3.285729 sqrt( Sum f_i^2 / 3N )
628
Res 3.285730 sqrt( Sum f_i^2 / 3N )
619
629
----------------------------------------
620
Max 6.707247 constrained
630
Max 6.707248 constrained
622
632
Stress-tensor-Voigt (kbar): 720.98 465.34 428.84 142.85 41.81 93.51
623
633
(Free)E + p*V (eV/cell) -2976.9707
668
678
Folding of H and S is implicitly performed
669
679
<dSpData1D:S at geom step 2
670
680
<sparsity:sparsity for geom step 2
671
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
672
<dData1D:(new from dSpData1D) n=306, refcount: 1>
681
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
682
<dData1D:(new from dSpData1D) n=204, refcount: 1>
675
685
Initializing Density Matrix...
676
DM after filling with atomic data:
677
<dSpData2D:(DM initialized from atoms)
686
DM filled with atomic data:
687
<dSpData2D:DM initialized from atoms
678
688
<sparsity:sparsity for geom step 2
679
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
680
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
689
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
690
<dData2D:DM n=204 m=1, refcount: 1>
682
693
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.58990828
694
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
683
695
New grid distribution: 1
689
705
InitMesh: MESH = 24 x 24 x 24 = 13824
690
706
InitMesh: (bp) = 12 x 12 x 12 = 1728
691
707
InitMesh: Mesh cutoff (required, used) = 50.000 63.484 Ry
692
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
708
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
693
709
New grid distribution: 2
698
718
New grid distribution: 3
703
727
Setting up quadratic distribution...
704
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
705
PhiOnMesh: Number of (b)points on node 0 = 432
706
PhiOnMesh: nlist on node 0 = 9991
708
siesta: Eharris(eV) = -2966.604740
711
siesta: Eharris(eV) = -2966.604740
728
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
729
PhiOnMesh: Number of (b)points on node 0 = 252
730
PhiOnMesh: nlist on node 0 = 5554
732
siesta: Eharris(eV) = -2966.604744
735
siesta: Eharris(eV) = -2966.604744
737
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
738
Geom step, scf iteration, dmax: 2 1 1.795020
714
740
Using DM_out to compute the final energy and forces
741
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
716
743
siesta: Atomic forces (eV/Ang):
717
744
1 0.000000 -0.000000 0.000000
718
745
2 -0.000000 0.000000 0.000000
719
3 -0.000000 0.000000 0.000000
720
4 -0.000000 -0.000000 -0.000000
721
5 2.995559 -6.699949 0.000000
722
6 -0.000000 4.019634 -6.140459
723
7 -2.995559 2.680315 6.140459
746
3 0.000000 -0.000001 0.000000
747
4 -0.000000 0.000001 -0.000000
748
5 2.995559 -6.699949 -0.000000
749
6 0.000000 4.019634 -6.140457
750
7 -2.995559 2.680314 6.140458
724
751
8 -2.995559 6.699949 0.000000
725
9 -0.000000 -4.019634 6.140459
726
10 2.995559 -2.680315 -6.140459
752
9 -0.000000 -4.019634 6.140457
753
10 2.995559 -2.680314 -6.140458
727
754
----------------------------------------
728
755
Tot -0.000000 -0.000000 -0.000000
729
756
----------------------------------------
731
Res 3.282155 sqrt( Sum f_i^2 / 3N )
758
Res 3.282154 sqrt( Sum f_i^2 / 3N )
732
759
----------------------------------------
733
760
Max 6.699949 constrained
781
808
Folding of H and S is implicitly performed
782
809
<dSpData1D:S at geom step 3
783
810
<sparsity:sparsity for geom step 3
784
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
785
<dData1D:(new from dSpData1D) n=306, refcount: 1>
811
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
812
<dData1D:(new from dSpData1D) n=204, refcount: 1>
788
815
Initializing Density Matrix...
789
DM after filling with atomic data:
790
<dSpData2D:(DM initialized from atoms)
816
DM filled with atomic data:
817
<dSpData2D:DM initialized from atoms
791
818
<sparsity:sparsity for geom step 3
792
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
793
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
819
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
820
<dData2D:DM n=204 m=1, refcount: 1>
795
823
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.59043989
824
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
796
825
New grid distribution: 1
802
835
InitMesh: MESH = 24 x 24 x 24 = 13824
803
836
InitMesh: (bp) = 12 x 12 x 12 = 1728
804
837
InitMesh: Mesh cutoff (required, used) = 50.000 63.487 Ry
805
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
838
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
806
839
New grid distribution: 2
811
848
New grid distribution: 3
816
857
Setting up quadratic distribution...
817
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
818
PhiOnMesh: Number of (b)points on node 0 = 432
819
PhiOnMesh: nlist on node 0 = 10003
821
siesta: Eharris(eV) = -2966.630710
824
siesta: Eharris(eV) = -2966.630710
858
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
859
PhiOnMesh: Number of (b)points on node 0 = 252
860
PhiOnMesh: nlist on node 0 = 5566
862
siesta: Eharris(eV) = -2966.630707
865
siesta: Eharris(eV) = -2966.630707
867
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
868
Geom step, scf iteration, dmax: 3 1 1.795026
827
870
Using DM_out to compute the final energy and forces
871
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
829
873
siesta: Atomic forces (eV/Ang):
830
874
1 0.000000 -0.000000 0.000000
831
875
2 -0.000000 0.000000 0.000000
832
3 -0.000000 0.000000 0.000000
833
4 -0.000000 -0.000000 -0.000000
834
5 2.984218 -6.674585 0.000000
835
6 -0.000000 4.004417 -6.117212
836
7 -2.984218 2.670168 6.117212
837
8 -2.984218 6.674585 0.000000
838
9 -0.000000 -4.004417 6.117212
839
10 2.984218 -2.670168 -6.117212
876
3 0.000000 -0.000001 -0.000000
877
4 -0.000000 0.000001 0.000000
878
5 2.984218 -6.674584 0.000000
879
6 -0.000001 4.004418 -6.117212
880
7 -2.984218 2.670167 6.117212
881
8 -2.984218 6.674584 -0.000000
882
9 0.000001 -4.004418 6.117212
883
10 2.984218 -2.670167 -6.117212
840
884
----------------------------------------
841
885
Tot -0.000000 -0.000000 -0.000000
842
886
----------------------------------------
844
888
Res 3.269729 sqrt( Sum f_i^2 / 3N )
845
889
----------------------------------------
846
Max 6.674585 constrained
890
Max 6.674584 constrained
848
892
Stress-tensor-Voigt (kbar): 721.17 466.44 430.07 142.34 41.66 93.18
849
893
(Free)E + p*V (eV/cell) -2977.0239
894
938
Folding of H and S is implicitly performed
895
939
<dSpData1D:S at geom step 4
896
940
<sparsity:sparsity for geom step 4
897
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
898
<dData1D:(new from dSpData1D) n=306, refcount: 1>
941
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
942
<dData1D:(new from dSpData1D) n=204, refcount: 1>
901
945
Initializing Density Matrix...
902
DM after filling with atomic data:
903
<dSpData2D:(DM initialized from atoms)
946
DM filled with atomic data:
947
<dSpData2D:DM initialized from atoms
904
948
<sparsity:sparsity for geom step 4
905
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
906
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
949
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
950
<dData2D:DM n=204 m=1, refcount: 1>
908
953
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.59123844
954
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
909
955
New grid distribution: 1
915
965
InitMesh: MESH = 24 x 24 x 24 = 13824
916
966
InitMesh: (bp) = 12 x 12 x 12 = 1728
917
967
InitMesh: Mesh cutoff (required, used) = 50.000 63.492 Ry
918
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
968
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
919
969
New grid distribution: 2
924
978
New grid distribution: 3
929
987
Setting up quadratic distribution...
930
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
931
PhiOnMesh: Number of (b)points on node 0 = 432
932
PhiOnMesh: nlist on node 0 = 10003
934
siesta: Eharris(eV) = -2966.673907
937
siesta: Eharris(eV) = -2966.673907
988
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
989
PhiOnMesh: Number of (b)points on node 0 = 252
990
PhiOnMesh: nlist on node 0 = 5566
992
siesta: Eharris(eV) = -2966.673907
995
siesta: Eharris(eV) = -2966.673907
997
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
998
Geom step, scf iteration, dmax: 4 1 1.795029
940
1000
Using DM_out to compute the final energy and forces
1001
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
942
1003
siesta: Atomic forces (eV/Ang):
943
1004
1 0.000000 -0.000000 0.000000
944
1005
2 -0.000000 0.000000 0.000000
945
3 -0.000000 -0.000000 0.000000
946
4 -0.000000 -0.000000 -0.000000
947
5 2.964988 -6.631575 0.000000
948
6 -0.000000 3.978613 -6.077794
949
7 -2.964988 2.652962 6.077794
950
8 -2.964988 6.631575 0.000000
951
9 -0.000000 -3.978613 6.077794
952
10 2.964988 -2.652962 -6.077794
1006
3 -0.000000 0.000000 0.000001
1007
4 0.000000 -0.000000 -0.000001
1008
5 2.964990 -6.631576 -0.000001
1009
6 0.000000 3.978613 -6.077795
1010
7 -2.964990 2.652962 6.077794
1011
8 -2.964990 6.631576 0.000001
1012
9 -0.000000 -3.978613 6.077795
1013
10 2.964990 -2.652962 -6.077794
953
1014
----------------------------------------
954
1015
Tot -0.000000 -0.000000 -0.000000
955
1016
----------------------------------------
957
Res 3.248659 sqrt( Sum f_i^2 / 3N )
1018
Res 3.248660 sqrt( Sum f_i^2 / 3N )
958
1019
----------------------------------------
959
Max 6.631575 constrained
1020
Max 6.631576 constrained
961
1022
Stress-tensor-Voigt (kbar): 722.05 468.69 432.51 141.58 41.44 92.68
962
1023
(Free)E + p*V (eV/cell) -2977.1415
1007
1068
Folding of H and S is implicitly performed
1008
1069
<dSpData1D:S at geom step 5
1009
1070
<sparsity:sparsity for geom step 5
1010
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
1011
<dData1D:(new from dSpData1D) n=306, refcount: 1>
1071
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
1072
<dData1D:(new from dSpData1D) n=204, refcount: 1>
1014
1075
Initializing Density Matrix...
1015
DM after filling with atomic data:
1016
<dSpData2D:(DM initialized from atoms)
1076
DM filled with atomic data:
1077
<dSpData2D:DM initialized from atoms
1017
1078
<sparsity:sparsity for geom step 5
1018
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
1019
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
1079
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
1080
<dData2D:DM n=204 m=1, refcount: 1>
1021
1083
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.59230592
1084
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1022
1085
New grid distribution: 1
1093
8 1: 12 7: 12 10: 12
1028
1095
InitMesh: MESH = 24 x 24 x 24 = 13824
1029
1096
InitMesh: (bp) = 12 x 12 x 12 = 1728
1030
1097
InitMesh: Mesh cutoff (required, used) = 50.000 63.499 Ry
1031
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1098
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
1032
1099
New grid distribution: 2
1037
1108
New grid distribution: 3
1042
1117
Setting up quadratic distribution...
1043
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1044
PhiOnMesh: Number of (b)points on node 0 = 432
1045
PhiOnMesh: nlist on node 0 = 10004
1047
siesta: Eharris(eV) = -2966.734047
1050
siesta: Eharris(eV) = -2966.734047
1118
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
1119
PhiOnMesh: Number of (b)points on node 0 = 252
1120
PhiOnMesh: nlist on node 0 = 5567
1122
siesta: Eharris(eV) = -2966.734043
1125
siesta: Eharris(eV) = -2966.734043
1127
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
1128
Geom step, scf iteration, dmax: 5 1 1.795030
1053
1130
Using DM_out to compute the final energy and forces
1131
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1055
1133
siesta: Atomic forces (eV/Ang):
1056
1134
1 0.000000 -0.000000 0.000000
1057
1135
2 -0.000000 0.000000 0.000000
1058
3 -0.000000 -0.000000 0.000000
1059
4 -0.000000 -0.000000 -0.000000
1136
3 -0.000000 0.000000 -0.000001
1137
4 0.000000 -0.000000 0.000001
1060
1138
5 2.938079 -6.571388 0.000000
1061
6 -0.000000 3.942504 -6.022633
1139
6 0.000000 3.942504 -6.022632
1062
1140
7 -2.938079 2.628884 6.022633
1063
8 -2.938079 6.571388 0.000000
1064
9 -0.000000 -3.942504 6.022633
1141
8 -2.938079 6.571388 -0.000000
1142
9 -0.000000 -3.942504 6.022632
1065
1143
10 2.938079 -2.628884 -6.022633
1066
1144
----------------------------------------
1067
Tot -0.000000 0.000000 -0.000000
1145
Tot -0.000000 -0.000000 -0.000000
1068
1146
----------------------------------------
1070
1148
Res 3.219175 sqrt( Sum f_i^2 / 3N )
1120
1198
Folding of H and S is implicitly performed
1121
1199
<dSpData1D:S at geom step 6
1122
1200
<sparsity:sparsity for geom step 6
1123
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
1124
<dData1D:(new from dSpData1D) n=306, refcount: 1>
1201
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
1202
<dData1D:(new from dSpData1D) n=204, refcount: 1>
1127
1205
Initializing Density Matrix...
1128
DM after filling with atomic data:
1129
<dSpData2D:(DM initialized from atoms)
1206
DM filled with atomic data:
1207
<dSpData2D:DM initialized from atoms
1130
1208
<sparsity:sparsity for geom step 6
1131
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
1132
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
1209
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
1210
<dData2D:DM n=204 m=1, refcount: 1>
1134
1213
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.59364513
1214
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1135
1215
New grid distribution: 1
1223
8 1: 12 7: 12 10: 12
1141
1225
InitMesh: MESH = 24 x 24 x 24 = 13824
1142
1226
InitMesh: (bp) = 12 x 12 x 12 = 1728
1143
1227
InitMesh: Mesh cutoff (required, used) = 50.000 63.507 Ry
1144
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1228
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
1145
1229
New grid distribution: 2
1150
1238
New grid distribution: 3
1155
1247
Setting up quadratic distribution...
1156
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1157
PhiOnMesh: Number of (b)points on node 0 = 432
1158
PhiOnMesh: nlist on node 0 = 9985
1160
siesta: Eharris(eV) = -2966.810258
1163
siesta: Eharris(eV) = -2966.810258
1248
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
1249
PhiOnMesh: Number of (b)points on node 0 = 252
1250
PhiOnMesh: nlist on node 0 = 5567
1252
siesta: Eharris(eV) = -2966.810256
1255
siesta: Eharris(eV) = -2966.810256
1257
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
1258
Geom step, scf iteration, dmax: 6 1 1.795033
1166
1260
Using DM_out to compute the final energy and forces
1261
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1168
1263
siesta: Atomic forces (eV/Ang):
1169
1264
1 0.000000 -0.000000 0.000000
1170
1265
2 -0.000000 0.000000 0.000000
1171
3 -0.000000 -0.000000 0.000000
1172
4 -0.000000 -0.000000 -0.000000
1173
5 2.893815 -6.472386 0.000000
1174
6 -0.000000 3.883107 -5.931898
1175
7 -2.893815 2.589278 5.931898
1176
8 -2.893815 6.472386 0.000000
1177
9 -0.000000 -3.883107 5.931898
1178
10 2.893815 -2.589278 -5.931898
1266
3 0.000000 0.000000 -0.000000
1267
4 -0.000000 -0.000000 0.000000
1268
5 2.893815 -6.472386 -0.000001
1269
6 -0.000000 3.883107 -5.931897
1270
7 -2.893815 2.589280 5.931898
1271
8 -2.893815 6.472386 0.000001
1272
9 0.000000 -3.883107 5.931897
1273
10 2.893815 -2.589280 -5.931898
1179
1274
----------------------------------------
1180
Tot -0.000000 0.000000 -0.000000
1275
Tot -0.000000 -0.000000 -0.000000
1181
1276
----------------------------------------
1183
1278
Res 3.170676 sqrt( Sum f_i^2 / 3N )
1191
1286
siesta: Stress tensor (static) (eV/Ang**3):
1192
1287
0.452342 0.085570 0.056015
1193
0.085570 0.299213 0.025045
1194
0.056015 0.025045 0.277349
1288
0.085570 0.299213 0.025044
1289
0.056015 0.025044 0.277349
1196
siesta: Pressure (static): -549.50118518 kBar
1291
siesta: Pressure (static): -549.50117807 kBar
1198
1293
siesta: Stress tensor (total) (eV/Ang**3):
1199
1294
0.452016 0.086081 0.056350
1200
1295
0.086081 0.297971 0.025194
1201
1296
0.056350 0.025194 0.275977
1203
siesta: Pressure (total): -547.93083700 kBar
1298
siesta: Pressure (total): -547.93082988 kBar
1205
1300
siesta: E_kin PR = 642.244 K
1233
1328
Folding of H and S is implicitly performed
1234
1329
<dSpData1D:S at geom step 7
1235
1330
<sparsity:sparsity for geom step 7
1236
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
1237
<dData1D:(new from dSpData1D) n=306, refcount: 1>
1331
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
1332
<dData1D:(new from dSpData1D) n=204, refcount: 1>
1240
1335
Initializing Density Matrix...
1241
DM after filling with atomic data:
1242
<dSpData2D:(DM initialized from atoms)
1336
DM filled with atomic data:
1337
<dSpData2D:DM initialized from atoms
1243
1338
<sparsity:sparsity for geom step 7
1244
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
1245
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
1339
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
1340
<dData2D:DM n=204 m=1, refcount: 1>
1247
1343
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.59526075
1344
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1248
1345
New grid distribution: 1
1353
8 1: 12 7: 12 10: 12
1254
1355
InitMesh: MESH = 24 x 24 x 24 = 13824
1255
1356
InitMesh: (bp) = 12 x 12 x 12 = 1728
1256
1357
InitMesh: Mesh cutoff (required, used) = 50.000 63.517 Ry
1257
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1358
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
1258
1359
New grid distribution: 2
1263
1368
New grid distribution: 3
1268
1377
Setting up quadratic distribution...
1269
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1270
PhiOnMesh: Number of (b)points on node 0 = 432
1271
PhiOnMesh: nlist on node 0 = 9985
1273
siesta: Eharris(eV) = -2966.901939
1276
siesta: Eharris(eV) = -2966.901939
1378
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
1379
PhiOnMesh: Number of (b)points on node 0 = 252
1380
PhiOnMesh: nlist on node 0 = 5567
1382
siesta: Eharris(eV) = -2966.901942
1385
siesta: Eharris(eV) = -2966.901942
1387
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
1388
Geom step, scf iteration, dmax: 7 1 1.795038
1279
1390
Using DM_out to compute the final energy and forces
1391
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1281
1393
siesta: Atomic forces (eV/Ang):
1282
1394
1 0.000000 -0.000000 0.000000
1283
1395
2 -0.000000 0.000000 0.000000
1284
1396
3 -0.000000 -0.000000 0.000000
1285
4 -0.000000 -0.000000 -0.000000
1286
5 2.844475 -6.362032 0.000000
1287
6 -0.000000 3.816901 -5.830760
1288
7 -2.844475 2.545131 5.830760
1289
8 -2.844475 6.362032 0.000000
1290
9 -0.000000 -3.816901 5.830760
1291
10 2.844475 -2.545131 -5.830760
1292
----------------------------------------
1293
Tot -0.000000 0.000000 -0.000000
1294
----------------------------------------
1397
4 0.000000 0.000000 -0.000000
1398
5 2.844476 -6.362033 0.000000
1399
6 -0.000000 3.816902 -5.830761
1400
7 -2.844476 2.545132 5.830760
1401
8 -2.844476 6.362033 -0.000000
1402
9 0.000000 -3.816902 5.830761
1403
10 2.844476 -2.545132 -5.830760
1404
----------------------------------------
1405
Tot -0.000000 -0.000000 -0.000000
1406
----------------------------------------
1296
1408
Res 3.116617 sqrt( Sum f_i^2 / 3N )
1297
1409
----------------------------------------
1298
Max 6.362032 constrained
1410
Max 6.362033 constrained
1300
1412
Stress-tensor-Voigt (kbar): 727.94 486.02 451.48 135.19 39.57 88.50
1301
1413
(Free)E + p*V (eV/cell) -2978.0057
1302
1414
Target enthalpy (eV/cell) -2940.0005
1304
1416
siesta: Stress tensor (static) (eV/Ang**3):
1305
0.454846 0.083583 0.054715
1417
0.454847 0.083583 0.054715
1306
1418
0.083583 0.305273 0.024463
1307
1419
0.054715 0.024463 0.283917
1309
siesta: Pressure (static): -557.58267845 kBar
1421
siesta: Pressure (static): -557.58271300 kBar
1311
1423
siesta: Stress tensor (total) (eV/Ang**3):
1312
1424
0.454340 0.084376 0.055234
1313
1425
0.084376 0.303347 0.024695
1314
0.055234 0.024695 0.281788
1426
0.055234 0.024695 0.281789
1316
siesta: Pressure (total): -555.14728748 kBar
1428
siesta: Pressure (total): -555.14732203 kBar
1318
1430
siesta: E_kin PR = 903.330 K
1346
1458
Folding of H and S is implicitly performed
1347
1459
<dSpData1D:S at geom step 8
1348
1460
<sparsity:sparsity for geom step 8
1349
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
1350
<dData1D:(new from dSpData1D) n=306, refcount: 1>
1461
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
1462
<dData1D:(new from dSpData1D) n=204, refcount: 1>
1353
1465
Initializing Density Matrix...
1354
DM after filling with atomic data:
1355
<dSpData2D:(DM initialized from atoms)
1466
DM filled with atomic data:
1467
<dSpData2D:DM initialized from atoms
1356
1468
<sparsity:sparsity for geom step 8
1357
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
1358
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
1469
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
1470
<dData2D:DM n=204 m=1, refcount: 1>
1360
1473
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.59715941
1474
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1361
1475
New grid distribution: 1
1483
8 1: 12 7: 12 10: 12
1367
1485
InitMesh: MESH = 24 x 24 x 24 = 13824
1368
1486
InitMesh: (bp) = 12 x 12 x 12 = 1728
1369
1487
InitMesh: Mesh cutoff (required, used) = 50.000 63.529 Ry
1370
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1488
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
1371
1489
New grid distribution: 2
1376
1498
New grid distribution: 3
1381
1507
Setting up quadratic distribution...
1382
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1383
PhiOnMesh: Number of (b)points on node 0 = 432
1384
PhiOnMesh: nlist on node 0 = 9998
1386
siesta: Eharris(eV) = -2967.007750
1389
siesta: Eharris(eV) = -2967.007750
1508
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
1509
PhiOnMesh: Number of (b)points on node 0 = 252
1510
PhiOnMesh: nlist on node 0 = 5577
1512
siesta: Eharris(eV) = -2967.007753
1515
siesta: Eharris(eV) = -2967.007753
1517
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
1518
Geom step, scf iteration, dmax: 8 1 1.795049
1392
1520
Using DM_out to compute the final energy and forces
1521
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1394
1523
siesta: Atomic forces (eV/Ang):
1395
1524
1 0.000000 -0.000000 0.000000
1396
1525
2 -0.000000 0.000000 0.000000
1397
3 -0.000000 -0.000000 0.000000
1526
3 0.000000 0.000000 0.000000
1398
1527
4 -0.000000 -0.000000 -0.000000
1399
1528
5 2.738140 -6.124198 0.000000
1400
6 -0.000000 3.674213 -5.612787
1529
6 0.000000 3.674213 -5.612788
1401
1530
7 -2.738140 2.449986 5.612787
1402
8 -2.738140 6.124198 0.000000
1403
9 -0.000000 -3.674213 5.612787
1531
8 -2.738140 6.124198 -0.000000
1532
9 -0.000000 -3.674213 5.612788
1404
1533
10 2.738140 -2.449986 -5.612787
1405
1534
----------------------------------------
1406
Tot -0.000000 0.000000 -0.000000
1535
Tot -0.000000 -0.000000 -0.000000
1407
1536
----------------------------------------
1409
1538
Res 3.000107 sqrt( Sum f_i^2 / 3N )
1459
1588
Folding of H and S is implicitly performed
1460
1589
<dSpData1D:S at geom step 9
1461
1590
<sparsity:sparsity for geom step 9
1462
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
1463
<dData1D:(new from dSpData1D) n=306, refcount: 1>
1591
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
1592
<dData1D:(new from dSpData1D) n=204, refcount: 1>
1466
1595
Initializing Density Matrix...
1467
DM after filling with atomic data:
1468
<dSpData2D:(DM initialized from atoms)
1596
DM filled with atomic data:
1597
<dSpData2D:DM initialized from atoms
1469
1598
<sparsity:sparsity for geom step 9
1470
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
1471
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
1599
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
1600
<dData2D:DM n=204 m=1, refcount: 1>
1473
1603
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.59935459
1604
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1474
1605
New grid distribution: 1
1613
8 1: 12 7: 12 10: 12
1480
1615
InitMesh: MESH = 24 x 24 x 24 = 13824
1481
1616
InitMesh: (bp) = 12 x 12 x 12 = 1728
1482
1617
InitMesh: Mesh cutoff (required, used) = 50.000 63.542 Ry
1483
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1618
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
1484
1619
New grid distribution: 2
1489
1628
New grid distribution: 3
1494
1637
Setting up quadratic distribution...
1495
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1496
PhiOnMesh: Number of (b)points on node 0 = 432
1497
PhiOnMesh: nlist on node 0 = 10012
1499
siesta: Eharris(eV) = -2967.125475
1502
siesta: Eharris(eV) = -2967.125475
1638
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
1639
PhiOnMesh: Number of (b)points on node 0 = 252
1640
PhiOnMesh: nlist on node 0 = 5593
1642
siesta: Eharris(eV) = -2967.125475
1645
siesta: Eharris(eV) = -2967.125475
1647
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
1648
Geom step, scf iteration, dmax: 9 1 1.795081
1505
1650
Using DM_out to compute the final energy and forces
1651
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1507
1653
siesta: Atomic forces (eV/Ang):
1508
1654
1 0.000000 -0.000000 0.000000
1509
1655
2 -0.000000 0.000000 0.000000
1510
3 -0.000000 -0.000000 -0.000000
1511
4 -0.000000 -0.000000 -0.000000
1512
5 2.679555 -5.993167 0.000000
1513
6 -0.000000 3.595600 -5.492698
1514
7 -2.679555 2.397567 5.492698
1515
8 -2.679555 5.993167 0.000000
1516
9 -0.000000 -3.595600 5.492698
1517
10 2.679555 -2.397567 -5.492698
1518
----------------------------------------
1519
Tot -0.000000 0.000000 -0.000000
1520
----------------------------------------
1656
3 -0.000000 0.000000 0.000000
1657
4 0.000000 -0.000000 -0.000000
1658
5 2.679555 -5.993166 0.000000
1659
6 -0.000000 3.595600 -5.492697
1660
7 -2.679554 2.397567 5.492697
1661
8 -2.679555 5.993166 -0.000000
1662
9 0.000000 -3.595600 5.492697
1663
10 2.679554 -2.397567 -5.492697
1664
----------------------------------------
1665
Tot -0.000000 -0.000000 -0.000000
1666
----------------------------------------
1522
1668
Res 2.935918 sqrt( Sum f_i^2 / 3N )
1523
1669
----------------------------------------
1524
Max 5.993167 constrained
1670
Max 5.993166 constrained
1526
1672
Stress-tensor-Voigt (kbar): 737.20 512.03 479.89 125.82 36.83 82.37
1527
1673
(Free)E + p*V (eV/cell) -2979.2695
1572
1718
Folding of H and S is implicitly performed
1573
1719
<dSpData1D:S at geom step 10
1574
1720
<sparsity:sparsity for geom step 10
1575
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
1576
<dData1D:(new from dSpData1D) n=306, refcount: 1>
1721
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
1722
<dData1D:(new from dSpData1D) n=204, refcount: 1>
1579
1725
Initializing Density Matrix...
1580
DM after filling with atomic data:
1581
<dSpData2D:(DM initialized from atoms)
1726
DM filled with atomic data:
1727
<dSpData2D:DM initialized from atoms
1582
1728
<sparsity:sparsity for geom step 10
1583
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
1584
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
1729
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
1730
<dData2D:DM n=204 m=1, refcount: 1>
1586
1733
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.60185267
1734
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1587
1735
New grid distribution: 1
1743
8 1: 12 7: 12 10: 12
1593
1745
InitMesh: MESH = 24 x 24 x 24 = 13824
1594
1746
InitMesh: (bp) = 12 x 12 x 12 = 1728
1595
1747
InitMesh: Mesh cutoff (required, used) = 50.000 63.557 Ry
1596
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1748
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
1597
1749
New grid distribution: 2
1602
1758
New grid distribution: 3
1607
1767
Setting up quadratic distribution...
1608
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1609
PhiOnMesh: Number of (b)points on node 0 = 432
1610
PhiOnMesh: nlist on node 0 = 10024
1612
siesta: Eharris(eV) = -2967.255031
1615
siesta: Eharris(eV) = -2967.255031
1768
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
1769
PhiOnMesh: Number of (b)points on node 0 = 252
1770
PhiOnMesh: nlist on node 0 = 5605
1772
siesta: Eharris(eV) = -2967.255029
1775
siesta: Eharris(eV) = -2967.255029
1777
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
1778
Geom step, scf iteration, dmax: 10 1 1.795128
1618
1780
Using DM_out to compute the final energy and forces
1781
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1620
1783
siesta: Atomic forces (eV/Ang):
1621
1784
1 0.000000 -0.000000 0.000000
1622
1785
2 -0.000000 0.000000 0.000000
1623
3 -0.000000 -0.000000 -0.000000
1624
4 -0.000000 -0.000000 -0.000000
1625
5 2.604368 -5.825002 0.000000
1626
6 -0.000000 3.494710 -5.338575
1627
7 -2.604368 2.330292 5.338575
1628
8 -2.604368 5.825002 0.000000
1629
9 -0.000000 -3.494710 5.338575
1630
10 2.604368 -2.330292 -5.338575
1786
3 0.000000 0.000000 -0.000001
1787
4 -0.000000 -0.000000 0.000001
1788
5 2.604369 -5.825002 -0.000001
1789
6 0.000001 3.494710 -5.338576
1790
7 -2.604368 2.330293 5.338576
1791
8 -2.604369 5.825002 0.000001
1792
9 -0.000001 -3.494710 5.338576
1793
10 2.604368 -2.330293 -5.338576
1631
1794
----------------------------------------
1632
1795
Tot -0.000000 -0.000000 -0.000000
1633
1796
----------------------------------------
1685
1848
Folding of H and S is implicitly performed
1686
1849
<dSpData1D:S at geom step 11
1687
1850
<sparsity:sparsity for geom step 11
1688
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
1689
<dData1D:(new from dSpData1D) n=306, refcount: 1>
1851
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
1852
<dData1D:(new from dSpData1D) n=204, refcount: 1>
1692
1855
Initializing Density Matrix...
1693
DM after filling with atomic data:
1694
<dSpData2D:(DM initialized from atoms)
1856
DM filled with atomic data:
1857
<dSpData2D:DM initialized from atoms
1695
1858
<sparsity:sparsity for geom step 11
1696
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
1697
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
1859
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
1860
<dData2D:DM n=204 m=1, refcount: 1>
1699
1863
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.60466106
1864
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1700
1865
New grid distribution: 1
1873
8 1: 12 7: 12 10: 12
1706
1875
InitMesh: MESH = 24 x 24 x 24 = 13824
1707
1876
InitMesh: (bp) = 12 x 12 x 12 = 1728
1708
1877
InitMesh: Mesh cutoff (required, used) = 50.000 63.574 Ry
1709
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1878
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
1710
1879
New grid distribution: 2
1715
1888
New grid distribution: 3
1720
1897
Setting up quadratic distribution...
1721
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1722
PhiOnMesh: Number of (b)points on node 0 = 432
1723
PhiOnMesh: nlist on node 0 = 10045
1725
siesta: Eharris(eV) = -2967.394667
1728
siesta: Eharris(eV) = -2967.394667
1898
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
1899
PhiOnMesh: Number of (b)points on node 0 = 252
1900
PhiOnMesh: nlist on node 0 = 5617
1902
siesta: Eharris(eV) = -2967.394665
1905
siesta: Eharris(eV) = -2967.394665
1907
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
1908
Geom step, scf iteration, dmax: 11 1 1.795193
1731
1910
Using DM_out to compute the final energy and forces
1911
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1733
1913
siesta: Atomic forces (eV/Ang):
1734
1914
1 0.000000 -0.000000 0.000000
1735
1915
2 -0.000000 0.000000 0.000000
1736
3 -0.000000 -0.000000 -0.000000
1737
4 -0.000000 -0.000000 -0.000000
1738
5 2.521180 -5.638941 0.000000
1739
6 -0.000000 3.383082 -5.168052
1740
7 -2.521180 2.255859 5.168052
1741
8 -2.521180 5.638941 0.000000
1742
9 -0.000000 -3.383082 5.168052
1743
10 2.521180 -2.255859 -5.168052
1916
3 -0.000000 -0.000000 0.000001
1917
4 0.000000 0.000000 -0.000001
1918
5 2.521180 -5.638941 0.000001
1919
6 0.000000 3.383083 -5.168053
1920
7 -2.521180 2.255858 5.168052
1921
8 -2.521180 5.638941 -0.000001
1922
9 -0.000000 -3.383083 5.168053
1923
10 2.521180 -2.255858 -5.168052
1744
1924
----------------------------------------
1745
1925
Tot -0.000000 -0.000000 -0.000000
1746
1926
----------------------------------------
1748
Res 2.762390 sqrt( Sum f_i^2 / 3N )
1928
Res 2.762391 sqrt( Sum f_i^2 / 3N )
1749
1929
----------------------------------------
1750
1930
Max 5.638941 constrained
1756
1936
siesta: Stress tensor (static) (eV/Ang**3):
1757
1937
0.465634 0.072617 0.047536
1758
0.072617 0.335684 0.021253
1759
0.047536 0.021253 0.317129
1938
0.072617 0.335684 0.021254
1939
0.047536 0.021254 0.317129
1761
siesta: Pressure (static): -597.32297797 kBar
1941
siesta: Pressure (static): -597.32292967 kBar
1763
1943
siesta: Stress tensor (total) (eV/Ang**3):
1764
1944
0.464087 0.075038 0.049121
1765
1945
0.075038 0.329804 0.021962
1766
1946
0.049121 0.021962 0.310631
1768
siesta: Pressure (total): -589.88642104 kBar
1948
siesta: Pressure (total): -589.88637259 kBar
1770
1950
siesta: E_kin PR = 2489.020 K
1780
1960
0.50000000 0.50000000 0.50000000 1 2 Mg
1781
1961
0.25000000 0.25000000 0.25000000 2 3 C
1782
1962
-0.25000000 -0.25000000 -0.25000000 2 4 C
1783
0.53242408 -0.03242408 0.25000000 3 5 O
1963
0.53242409 -0.03242408 0.25000000 3 5 O
1784
1964
0.25000000 0.53242408 -0.03242408 3 6 O
1785
-0.03242408 0.25000000 0.53242408 3 7 O
1786
-0.53242408 0.03242408 -0.25000000 3 8 O
1965
-0.03242409 0.25000000 0.53242408 3 7 O
1966
-0.53242409 0.03242408 -0.25000000 3 8 O
1787
1967
-0.25000000 -0.53242408 0.03242408 3 9 O
1788
0.03242408 -0.25000000 -0.53242408 3 10 O
1968
0.03242409 -0.25000000 -0.53242408 3 10 O
1790
1970
outcell: Unit cell vectors (Ang):
1791
1971
5.669739 0.005004 0.003276
1798
1978
Folding of H and S is implicitly performed
1799
1979
<dSpData1D:S at geom step 12
1800
1980
<sparsity:sparsity for geom step 12
1801
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
1802
<dData1D:(new from dSpData1D) n=306, refcount: 1>
1981
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
1982
<dData1D:(new from dSpData1D) n=204, refcount: 1>
1805
1985
Initializing Density Matrix...
1806
DM after filling with atomic data:
1807
<dSpData2D:(DM initialized from atoms)
1986
DM filled with atomic data:
1987
<dSpData2D:DM initialized from atoms
1808
1988
<sparsity:sparsity for geom step 12
1809
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
1810
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
1989
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
1990
<dData2D:DM n=204 m=1, refcount: 1>
1812
1993
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.60778792
1994
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1813
1995
New grid distribution: 1
2003
8 1: 12 7: 12 10: 12
1819
2005
InitMesh: MESH = 24 x 24 x 24 = 13824
1820
2006
InitMesh: (bp) = 12 x 12 x 12 = 1728
1821
2007
InitMesh: Mesh cutoff (required, used) = 50.000 63.592 Ry
1822
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2008
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
1823
2009
New grid distribution: 2
1828
2018
New grid distribution: 3
1833
2027
Setting up quadratic distribution...
1834
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1835
PhiOnMesh: Number of (b)points on node 0 = 432
1836
PhiOnMesh: nlist on node 0 = 10048
1838
siesta: Eharris(eV) = -2967.541044
1841
siesta: Eharris(eV) = -2967.541044
2028
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
2029
PhiOnMesh: Number of (b)points on node 0 = 252
2030
PhiOnMesh: nlist on node 0 = 5620
2032
siesta: Eharris(eV) = -2967.541045
2035
siesta: Eharris(eV) = -2967.541045
2037
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
2038
Geom step, scf iteration, dmax: 12 1 1.795278
1844
2040
Using DM_out to compute the final energy and forces
2041
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1846
2043
siesta: Atomic forces (eV/Ang):
1847
2044
1 0.000000 -0.000000 0.000000
1848
2045
2 -0.000000 0.000000 0.000000
1849
3 -0.000000 -0.000000 -0.000000
1850
4 -0.000000 -0.000000 -0.000000
1851
5 2.349544 -5.255055 0.000000
1852
6 -0.000000 3.152770 -4.816223
1853
7 -2.349544 2.102285 4.816223
2046
3 0.000000 -0.000001 0.000002
2047
4 -0.000000 0.000001 -0.000002
2048
5 2.349544 -5.255055 -0.000000
2049
6 -0.000000 3.152770 -4.816225
2050
7 -2.349545 2.102286 4.816224
1854
2051
8 -2.349544 5.255055 0.000000
1855
9 -0.000000 -3.152770 4.816223
1856
10 2.349544 -2.102285 -4.816223
2052
9 0.000000 -3.152770 4.816225
2053
10 2.349545 -2.102286 -4.816224
1857
2054
----------------------------------------
1858
2055
Tot -0.000000 -0.000000 -0.000000
1859
2056
----------------------------------------
1911
2108
Folding of H and S is implicitly performed
1912
2109
<dSpData1D:S at geom step 13
1913
2110
<sparsity:sparsity for geom step 13
1914
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
1915
<dData1D:(new from dSpData1D) n=306, refcount: 1>
2111
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
2112
<dData1D:(new from dSpData1D) n=204, refcount: 1>
1918
2115
Initializing Density Matrix...
1919
DM after filling with atomic data:
1920
<dSpData2D:(DM initialized from atoms)
2116
DM filled with atomic data:
2117
<dSpData2D:DM initialized from atoms
1921
2118
<sparsity:sparsity for geom step 13
1922
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
1923
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
2119
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
2120
<dData2D:DM n=204 m=1, refcount: 1>
1925
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.61125667
2123
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.61125666
2124
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1926
2125
New grid distribution: 1
2133
8 1: 12 7: 12 10: 12
1932
2135
InitMesh: MESH = 24 x 24 x 24 = 13824
1933
2136
InitMesh: (bp) = 12 x 12 x 12 = 1728
1934
2137
InitMesh: Mesh cutoff (required, used) = 50.000 63.612 Ry
1935
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2138
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
1936
2139
New grid distribution: 2
1941
2148
New grid distribution: 3
1946
2157
Setting up quadratic distribution...
1947
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
1948
PhiOnMesh: Number of (b)points on node 0 = 432
1949
PhiOnMesh: nlist on node 0 = 10060
1951
siesta: Eharris(eV) = -2967.692957
1954
siesta: Eharris(eV) = -2967.692957
2158
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
2159
PhiOnMesh: Number of (b)points on node 0 = 252
2160
PhiOnMesh: nlist on node 0 = 5626
2162
siesta: Eharris(eV) = -2967.692957
2165
siesta: Eharris(eV) = -2967.692957
2167
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
2168
Geom step, scf iteration, dmax: 13 1 1.795373
1957
2170
Using DM_out to compute the final energy and forces
2171
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
1959
2173
siesta: Atomic forces (eV/Ang):
1960
2174
1 0.000000 -0.000000 0.000000
1961
2175
2 -0.000000 0.000000 0.000000
1962
3 -0.000000 -0.000000 -0.000000
1963
4 -0.000000 -0.000000 -0.000000
1964
5 2.241020 -5.012326 0.000000
1965
6 -0.000000 3.007145 -4.593763
1966
7 -2.241020 2.005181 4.593763
1967
8 -2.241020 5.012326 0.000000
1968
9 -0.000000 -3.007145 4.593763
1969
10 2.241020 -2.005181 -4.593763
2176
3 0.000001 -0.000001 -0.000000
2177
4 -0.000001 0.000001 0.000000
2178
5 2.241019 -5.012325 0.000000
2179
6 -0.000000 3.007146 -4.593764
2180
7 -2.241020 2.005179 4.593764
2181
8 -2.241019 5.012325 -0.000000
2182
9 0.000000 -3.007146 4.593764
2183
10 2.241020 -2.005179 -4.593764
1970
2184
----------------------------------------
1971
2185
Tot -0.000000 -0.000000 -0.000000
1972
2186
----------------------------------------
1974
2188
Res 2.455426 sqrt( Sum f_i^2 / 3N )
1975
2189
----------------------------------------
1976
Max 5.012326 constrained
2190
Max 5.012325 constrained
1978
2192
Stress-tensor-Voigt (kbar): 760.49 571.67 544.72 105.51 30.88 69.07
1979
2193
(Free)E + p*V (eV/cell) -2982.1138
2024
2238
Folding of H and S is implicitly performed
2025
2239
<dSpData1D:S at geom step 14
2026
2240
<sparsity:sparsity for geom step 14
2027
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
2028
<dData1D:(new from dSpData1D) n=306, refcount: 1>
2241
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
2242
<dData1D:(new from dSpData1D) n=204, refcount: 1>
2031
2245
Initializing Density Matrix...
2032
DM after filling with atomic data:
2033
<dSpData2D:(DM initialized from atoms)
2246
DM filled with atomic data:
2247
<dSpData2D:DM initialized from atoms
2034
2248
<sparsity:sparsity for geom step 14
2035
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
2036
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
2249
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
2250
<dData2D:DM n=204 m=1, refcount: 1>
2038
2253
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.61507784
2254
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2039
2255
New grid distribution: 1
2263
8 1: 12 7: 12 10: 12
2045
2265
InitMesh: MESH = 24 x 24 x 24 = 13824
2046
2266
InitMesh: (bp) = 12 x 12 x 12 = 1728
2047
2267
InitMesh: Mesh cutoff (required, used) = 50.000 63.633 Ry
2048
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2268
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
2049
2269
New grid distribution: 2
2054
2278
New grid distribution: 3
2059
2287
Setting up quadratic distribution...
2060
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2061
PhiOnMesh: Number of (b)points on node 0 = 432
2062
PhiOnMesh: nlist on node 0 = 10086
2064
siesta: Eharris(eV) = -2967.849049
2067
siesta: Eharris(eV) = -2967.849049
2288
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
2289
PhiOnMesh: Number of (b)points on node 0 = 252
2290
PhiOnMesh: nlist on node 0 = 5645
2292
siesta: Eharris(eV) = -2967.849053
2295
siesta: Eharris(eV) = -2967.849053
2297
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
2298
Geom step, scf iteration, dmax: 14 1 1.795506
2070
2300
Using DM_out to compute the final energy and forces
2301
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2072
2303
siesta: Atomic forces (eV/Ang):
2073
2304
1 0.000000 -0.000000 0.000000
2074
2305
2 -0.000000 0.000000 0.000000
2075
3 -0.000000 -0.000000 -0.000000
2076
4 -0.000000 -0.000000 -0.000000
2077
5 2.120882 -4.743624 0.000000
2078
6 -0.000000 2.845937 -4.347500
2079
7 -2.120882 1.897687 4.347500
2080
8 -2.120882 4.743624 0.000000
2081
9 -0.000000 -2.845937 4.347500
2082
10 2.120882 -1.897687 -4.347500
2306
3 0.000001 0.000000 -0.000000
2307
4 -0.000001 -0.000000 0.000000
2308
5 2.120883 -4.743624 -0.000001
2309
6 -0.000001 2.845937 -4.347500
2310
7 -2.120883 1.897687 4.347500
2311
8 -2.120883 4.743624 0.000001
2312
9 0.000001 -2.845937 4.347500
2313
10 2.120883 -1.897687 -4.347500
2083
2314
----------------------------------------
2084
Tot -0.000000 0.000000 -0.000000
2315
Tot -0.000000 -0.000000 -0.000000
2085
2316
----------------------------------------
2087
2318
Res 2.323795 sqrt( Sum f_i^2 / 3N )
2137
2368
Folding of H and S is implicitly performed
2138
2369
<dSpData1D:S at geom step 15
2139
2370
<sparsity:sparsity for geom step 15
2140
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
2141
<dData1D:(new from dSpData1D) n=306, refcount: 1>
2371
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
2372
<dData1D:(new from dSpData1D) n=204, refcount: 1>
2144
2375
Initializing Density Matrix...
2145
DM after filling with atomic data:
2146
<dSpData2D:(DM initialized from atoms)
2376
DM filled with atomic data:
2377
<dSpData2D:DM initialized from atoms
2147
2378
<sparsity:sparsity for geom step 15
2148
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
2149
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
2379
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
2380
<dData2D:DM n=204 m=1, refcount: 1>
2151
2383
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.61926305
2384
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2152
2385
New grid distribution: 1
2393
8 1: 12 7: 12 10: 12
2158
2395
InitMesh: MESH = 24 x 24 x 24 = 13824
2159
2396
InitMesh: (bp) = 12 x 12 x 12 = 1728
2160
2397
InitMesh: Mesh cutoff (required, used) = 50.000 63.657 Ry
2161
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2398
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
2162
2399
New grid distribution: 2
2167
2408
New grid distribution: 3
2172
2417
Setting up quadratic distribution...
2173
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2174
PhiOnMesh: Number of (b)points on node 0 = 432
2175
PhiOnMesh: nlist on node 0 = 10094
2177
siesta: Eharris(eV) = -2968.007264
2180
siesta: Eharris(eV) = -2968.007264
2418
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
2419
PhiOnMesh: Number of (b)points on node 0 = 252
2420
PhiOnMesh: nlist on node 0 = 5644
2422
siesta: Eharris(eV) = -2968.007266
2425
siesta: Eharris(eV) = -2968.007266
2427
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
2428
Geom step, scf iteration, dmax: 15 1 1.795681
2183
2430
Using DM_out to compute the final energy and forces
2431
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2185
2433
siesta: Atomic forces (eV/Ang):
2186
2434
1 0.000000 -0.000000 0.000000
2187
2435
2 -0.000000 0.000000 0.000000
2188
3 -0.000000 -0.000000 -0.000000
2189
4 -0.000000 -0.000000 -0.000000
2436
3 0.000000 0.000001 -0.000001
2437
4 -0.000000 -0.000001 0.000001
2190
2438
5 1.999624 -4.472414 0.000000
2191
6 -0.000000 2.683225 -4.098937
2192
7 -1.999624 1.789189 4.098937
2193
8 -1.999624 4.472414 0.000000
2439
6 0.000000 2.683225 -4.098937
2440
7 -1.999625 1.789190 4.098938
2441
8 -1.999624 4.472414 -0.000000
2194
2442
9 -0.000000 -2.683225 4.098937
2195
10 1.999624 -1.789189 -4.098937
2443
10 1.999625 -1.789190 -4.098938
2196
2444
----------------------------------------
2197
Tot -0.000000 -0.000000 -0.000000
2445
Tot -0.000000 0.000000 -0.000000
2198
2446
----------------------------------------
2200
2448
Res 2.190935 sqrt( Sum f_i^2 / 3N )
2250
2498
Folding of H and S is implicitly performed
2251
2499
<dSpData1D:S at geom step 16
2252
2500
<sparsity:sparsity for geom step 16
2253
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
2254
<dData1D:(new from dSpData1D) n=306, refcount: 1>
2501
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
2502
<dData1D:(new from dSpData1D) n=204, refcount: 1>
2257
2505
Initializing Density Matrix...
2258
DM after filling with atomic data:
2259
<dSpData2D:(DM initialized from atoms)
2506
DM filled with atomic data:
2507
<dSpData2D:DM initialized from atoms
2260
2508
<sparsity:sparsity for geom step 16
2261
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
2262
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
2509
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
2510
<dData2D:DM n=204 m=1, refcount: 1>
2264
2513
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.62382277
2514
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2265
2515
New grid distribution: 1
2523
8 1: 12 7: 12 10: 12
2271
2525
InitMesh: MESH = 24 x 24 x 24 = 13824
2272
2526
InitMesh: (bp) = 12 x 12 x 12 = 1728
2273
2527
InitMesh: Mesh cutoff (required, used) = 50.000 63.682 Ry
2274
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2528
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
2275
2529
New grid distribution: 2
2280
2538
New grid distribution: 3
2285
2547
Setting up quadratic distribution...
2286
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2287
PhiOnMesh: Number of (b)points on node 0 = 432
2288
PhiOnMesh: nlist on node 0 = 10095
2290
siesta: Eharris(eV) = -2968.165848
2293
siesta: Eharris(eV) = -2968.165848
2548
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
2549
PhiOnMesh: Number of (b)points on node 0 = 252
2550
PhiOnMesh: nlist on node 0 = 5649
2552
siesta: Eharris(eV) = -2968.165850
2555
siesta: Eharris(eV) = -2968.165850
2557
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
2558
Geom step, scf iteration, dmax: 16 1 1.795899
2296
2560
Using DM_out to compute the final energy and forces
2561
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2298
2563
siesta: Atomic forces (eV/Ang):
2299
2564
1 0.000000 -0.000000 0.000000
2300
2565
2 -0.000000 0.000000 0.000000
2301
3 -0.000000 -0.000000 -0.000000
2302
4 -0.000000 -0.000000 -0.000000
2303
5 1.866363 -4.174358 0.000000
2304
6 -0.000000 2.504406 -3.825771
2305
7 -1.866363 1.669952 3.825771
2306
8 -1.866363 4.174358 0.000000
2307
9 -0.000000 -2.504406 3.825771
2308
10 1.866363 -1.669952 -3.825771
2566
3 -0.000001 -0.000000 0.000000
2567
4 0.000001 0.000000 -0.000000
2568
5 1.866362 -4.174358 -0.000000
2569
6 0.000000 2.504407 -3.825771
2570
7 -1.866361 1.669952 3.825771
2571
8 -1.866362 4.174358 0.000000
2572
9 -0.000000 -2.504407 3.825771
2573
10 1.866361 -1.669952 -3.825771
2309
2574
----------------------------------------
2310
Tot -0.000000 0.000000 -0.000000
2575
Tot -0.000000 -0.000000 -0.000000
2311
2576
----------------------------------------
2313
2578
Res 2.044924 sqrt( Sum f_i^2 / 3N )
2321
2586
siesta: Stress tensor (static) (eV/Ang**3):
2322
2587
0.485873 0.055139 0.036095
2323
0.055139 0.387199 0.016138
2588
0.055139 0.387200 0.016138
2324
2589
0.036095 0.016138 0.373111
2326
siesta: Pressure (static): -665.54241574 kBar
2591
siesta: Pressure (static): -665.54245335 kBar
2328
2593
siesta: Stress tensor (total) (eV/Ang**3):
2329
2594
0.482599 0.060266 0.039451
2330
2595
0.060266 0.374751 0.017639
2331
2596
0.039451 0.017639 0.359353
2333
siesta: Pressure (total): -649.79814944 kBar
2598
siesta: Pressure (total): -649.79818621 kBar
2335
2600
siesta: E_kin PR = 5982.393 K
2337
siesta: Temp_ion = 1300.399 K
2602
siesta: Temp_ion = 1300.400 K
2339
2604
====================================
2340
2605
Begin MD step = 17
2363
2628
Folding of H and S is implicitly performed
2364
2629
<dSpData1D:S at geom step 17
2365
2630
<sparsity:sparsity for geom step 17
2366
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
2367
<dData1D:(new from dSpData1D) n=306, refcount: 1>
2631
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
2632
<dData1D:(new from dSpData1D) n=204, refcount: 1>
2370
2635
Initializing Density Matrix...
2371
DM after filling with atomic data:
2372
<dSpData2D:(DM initialized from atoms)
2636
DM filled with atomic data:
2637
<dSpData2D:DM initialized from atoms
2373
2638
<sparsity:sparsity for geom step 17
2374
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
2375
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
2639
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
2640
<dData2D:DM n=204 m=1, refcount: 1>
2377
2643
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.62876809
2644
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2378
2645
New grid distribution: 1
2653
8 1: 12 7: 12 10: 12
2384
2655
InitMesh: MESH = 24 x 24 x 24 = 13824
2385
2656
InitMesh: (bp) = 12 x 12 x 12 = 1728
2386
2657
InitMesh: Mesh cutoff (required, used) = 50.000 63.709 Ry
2387
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2658
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
2388
2659
New grid distribution: 2
2393
2668
New grid distribution: 3
2398
2677
Setting up quadratic distribution...
2399
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2400
PhiOnMesh: Number of (b)points on node 0 = 432
2401
PhiOnMesh: nlist on node 0 = 10109
2403
siesta: Eharris(eV) = -2968.322741
2406
siesta: Eharris(eV) = -2968.322741
2678
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
2679
PhiOnMesh: Number of (b)points on node 0 = 252
2680
PhiOnMesh: nlist on node 0 = 5660
2682
siesta: Eharris(eV) = -2968.322738
2685
siesta: Eharris(eV) = -2968.322738
2687
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
2688
Geom step, scf iteration, dmax: 17 1 1.796159
2409
2690
Using DM_out to compute the final energy and forces
2691
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2411
2693
siesta: Atomic forces (eV/Ang):
2412
2694
1 0.000000 -0.000000 0.000000
2413
2695
2 -0.000000 0.000000 0.000000
2414
3 -0.000000 -0.000000 -0.000000
2415
4 -0.000000 -0.000000 -0.000000
2416
5 1.727130 -3.862947 0.000000
2417
6 -0.000000 2.317575 -3.540365
2418
7 -1.727130 1.545372 3.540365
2419
8 -1.727130 3.862947 0.000000
2420
9 -0.000000 -2.317575 3.540365
2421
10 1.727130 -1.545372 -3.540365
2422
----------------------------------------
2423
Tot -0.000000 0.000000 -0.000000
2424
----------------------------------------
2696
3 0.000001 0.000000 -0.000001
2697
4 -0.000001 -0.000000 0.000001
2698
5 1.727129 -3.862946 -0.000001
2699
6 -0.000000 2.317574 -3.540364
2700
7 -1.727130 1.545371 3.540365
2701
8 -1.727129 3.862946 0.000001
2702
9 0.000000 -2.317574 3.540364
2703
10 1.727130 -1.545371 -3.540365
2704
----------------------------------------
2705
Tot -0.000000 -0.000000 -0.000000
2706
----------------------------------------
2426
2708
Res 1.892371 sqrt( Sum f_i^2 / 3N )
2427
2709
----------------------------------------
2428
Max 3.862947 constrained
2710
Max 3.862946 constrained
2430
2712
Stress-tensor-Voigt (kbar): 778.13 610.11 586.12 93.89 27.48 61.46
2431
2713
(Free)E + p*V (eV/cell) -2983.7776
2476
2758
Folding of H and S is implicitly performed
2477
2759
<dSpData1D:S at geom step 18
2478
2760
<sparsity:sparsity for geom step 18
2479
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
2480
<dData1D:(new from dSpData1D) n=306, refcount: 1>
2761
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
2762
<dData1D:(new from dSpData1D) n=204, refcount: 1>
2483
2765
Initializing Density Matrix...
2484
DM after filling with atomic data:
2485
<dSpData2D:(DM initialized from atoms)
2766
DM filled with atomic data:
2767
<dSpData2D:DM initialized from atoms
2486
2768
<sparsity:sparsity for geom step 18
2487
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
2488
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
2769
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
2770
<dData2D:DM n=204 m=1, refcount: 1>
2490
2773
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.63411128
2774
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2491
2775
New grid distribution: 1
2783
8 1: 12 7: 12 10: 12
2497
2785
InitMesh: MESH = 24 x 24 x 24 = 13824
2498
2786
InitMesh: (bp) = 12 x 12 x 12 = 1728
2499
2787
InitMesh: Mesh cutoff (required, used) = 50.000 63.737 Ry
2500
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2788
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
2501
2789
New grid distribution: 2
2506
2798
New grid distribution: 3
2511
2807
Setting up quadratic distribution...
2512
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2513
PhiOnMesh: Number of (b)points on node 0 = 432
2514
PhiOnMesh: nlist on node 0 = 10112
2516
siesta: Eharris(eV) = -2968.475869
2519
siesta: Eharris(eV) = -2968.475869
2808
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
2809
PhiOnMesh: Number of (b)points on node 0 = 252
2810
PhiOnMesh: nlist on node 0 = 5666
2812
siesta: Eharris(eV) = -2968.475869
2815
siesta: Eharris(eV) = -2968.475869
2817
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
2818
Geom step, scf iteration, dmax: 18 1 1.796449
2522
2820
Using DM_out to compute the final energy and forces
2821
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2524
2823
siesta: Atomic forces (eV/Ang):
2525
2824
1 0.000000 -0.000000 0.000000
2526
2 -0.000000 0.000000 0.000000
2527
3 -0.000000 -0.000000 -0.000000
2528
4 -0.000000 -0.000000 -0.000000
2529
5 1.564995 -3.500312 0.000000
2530
6 -0.000000 2.100012 -3.208012
2531
7 -1.564995 1.400300 3.208012
2532
8 -1.564995 3.500312 0.000000
2533
9 -0.000000 -2.100012 3.208012
2534
10 1.564995 -1.400300 -3.208012
2825
2 0.000000 0.000000 0.000000
2826
3 -0.000000 -0.000000 0.000000
2827
4 0.000000 0.000000 -0.000000
2828
5 1.564994 -3.500311 -0.000000
2829
6 -0.000001 2.100012 -3.208011
2830
7 -1.564994 1.400299 3.208012
2831
8 -1.564994 3.500311 0.000000
2832
9 0.000001 -2.100012 3.208011
2833
10 1.564994 -1.400299 -3.208012
2535
2834
----------------------------------------
2536
2835
Tot -0.000000 -0.000000 -0.000000
2537
2836
----------------------------------------
2539
2838
Res 1.714724 sqrt( Sum f_i^2 / 3N )
2540
2839
----------------------------------------
2541
Max 3.500312 constrained
2840
Max 3.500311 constrained
2543
2842
Stress-tensor-Voigt (kbar): 785.31 625.18 602.31 89.49 26.19 58.58
2544
2843
(Free)E + p*V (eV/cell) -2984.4152
2589
2888
Folding of H and S is implicitly performed
2590
2889
<dSpData1D:S at geom step 19
2591
2890
<sparsity:sparsity for geom step 19
2592
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
2593
<dData1D:(new from dSpData1D) n=306, refcount: 1>
2891
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
2892
<dData1D:(new from dSpData1D) n=204, refcount: 1>
2596
2895
Initializing Density Matrix...
2597
DM after filling with atomic data:
2598
<dSpData2D:(DM initialized from atoms)
2896
DM filled with atomic data:
2897
<dSpData2D:DM initialized from atoms
2599
2898
<sparsity:sparsity for geom step 19
2600
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
2601
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
2899
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
2900
<dData2D:DM n=204 m=1, refcount: 1>
2603
2903
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.63986865
2904
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2604
2905
New grid distribution: 1
2913
8 1: 12 7: 12 10: 12
2610
2915
InitMesh: MESH = 24 x 24 x 24 = 13824
2611
2916
InitMesh: (bp) = 12 x 12 x 12 = 1728
2612
2917
InitMesh: Mesh cutoff (required, used) = 50.000 63.767 Ry
2613
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2918
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
2614
2919
New grid distribution: 2
2619
2928
New grid distribution: 3
2624
2937
Setting up quadratic distribution...
2625
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2626
PhiOnMesh: Number of (b)points on node 0 = 432
2627
PhiOnMesh: nlist on node 0 = 10141
2629
siesta: Eharris(eV) = -2968.623561
2632
siesta: Eharris(eV) = -2968.623561
2938
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
2939
PhiOnMesh: Number of (b)points on node 0 = 252
2940
PhiOnMesh: nlist on node 0 = 5689
2942
siesta: Eharris(eV) = -2968.623562
2945
siesta: Eharris(eV) = -2968.623562
2947
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
2948
Geom step, scf iteration, dmax: 19 1 1.796715
2635
2950
Using DM_out to compute the final energy and forces
2951
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2637
2953
siesta: Atomic forces (eV/Ang):
2638
2954
1 0.000000 -0.000000 0.000000
2639
2 -0.000000 0.000000 0.000000
2640
3 -0.000000 -0.000000 -0.000000
2641
4 -0.000000 -0.000000 -0.000000
2955
2 0.000000 0.000000 0.000000
2956
3 -0.000000 0.000001 0.000001
2957
4 0.000000 -0.000001 -0.000001
2642
2958
5 1.407123 -3.147211 0.000000
2643
6 -0.000000 1.888169 -2.884398
2644
7 -1.407123 1.259042 2.884398
2645
8 -1.407123 3.147211 0.000000
2646
9 -0.000000 -1.888169 2.884398
2647
10 1.407123 -1.259042 -2.884398
2959
6 -0.000000 1.888169 -2.884397
2960
7 -1.407123 1.259041 2.884397
2961
8 -1.407123 3.147211 -0.000000
2962
9 0.000000 -1.888169 2.884397
2963
10 1.407123 -1.259041 -2.884397
2648
2964
----------------------------------------
2649
Tot -0.000000 -0.000000 -0.000000
2965
Tot -0.000000 -0.000000 0.000000
2650
2966
----------------------------------------
2652
2968
Res 1.541748 sqrt( Sum f_i^2 / 3N )
2702
3018
Folding of H and S is implicitly performed
2703
3019
<dSpData1D:S at geom step 20
2704
3020
<sparsity:sparsity for geom step 20
2705
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
2706
<dData1D:(new from dSpData1D) n=306, refcount: 1>
3021
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
3022
<dData1D:(new from dSpData1D) n=204, refcount: 1>
2709
3025
Initializing Density Matrix...
2710
DM after filling with atomic data:
2711
<dSpData2D:(DM initialized from atoms)
3026
DM filled with atomic data:
3027
<dSpData2D:DM initialized from atoms
2712
3028
<sparsity:sparsity for geom step 20
2713
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
2714
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
3029
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
3030
<dData2D:DM n=204 m=1, refcount: 1>
2716
3033
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.64605798
3034
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2717
3035
New grid distribution: 1
3043
8 1: 12 7: 12 10: 12
2723
3045
InitMesh: MESH = 24 x 24 x 24 = 13824
2724
3046
InitMesh: (bp) = 12 x 12 x 12 = 1728
2725
3047
InitMesh: Mesh cutoff (required, used) = 50.000 63.799 Ry
2726
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
3048
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
2727
3049
New grid distribution: 2
2732
3058
New grid distribution: 3
2737
3067
Setting up quadratic distribution...
2738
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2739
PhiOnMesh: Number of (b)points on node 0 = 432
2740
PhiOnMesh: nlist on node 0 = 10168
2742
siesta: Eharris(eV) = -2968.763427
2745
siesta: Eharris(eV) = -2968.763427
3068
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
3069
PhiOnMesh: Number of (b)points on node 0 = 252
3070
PhiOnMesh: nlist on node 0 = 5697
3072
siesta: Eharris(eV) = -2968.763428
3075
siesta: Eharris(eV) = -2968.763428
3077
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
3078
Geom step, scf iteration, dmax: 20 1 1.797006
2748
3080
Using DM_out to compute the final energy and forces
3081
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2750
3083
siesta: Atomic forces (eV/Ang):
2751
3084
1 0.000000 -0.000000 0.000000
2752
2 -0.000000 0.000000 0.000000
2753
3 -0.000000 -0.000000 -0.000000
2754
4 -0.000000 -0.000000 -0.000000
2755
5 1.210650 -2.707774 0.000000
2756
6 -0.000000 1.624529 -2.481657
3085
2 0.000000 0.000000 0.000000
3086
3 -0.000000 -0.000000 -0.000001
3087
4 0.000000 0.000000 0.000001
3088
5 1.210651 -2.707775 0.000001
3089
6 -0.000000 1.624530 -2.481657
2757
3090
7 -1.210650 1.083245 2.481657
2758
8 -1.210650 2.707774 0.000000
2759
9 -0.000000 -1.624529 2.481657
3091
8 -1.210651 2.707775 -0.000001
3092
9 0.000000 -1.624530 2.481657
2760
3093
10 1.210650 -1.083245 -2.481657
2761
3094
----------------------------------------
2762
Tot -0.000000 -0.000000 -0.000000
3095
Tot -0.000000 -0.000000 0.000000
2763
3096
----------------------------------------
2765
3098
Res 1.326478 sqrt( Sum f_i^2 / 3N )
2766
3099
----------------------------------------
2767
Max 2.707774 constrained
3100
Max 2.707775 constrained
2769
Stress-tensor-Voigt (kbar): 801.74 662.11 642.18 78.03 22.84 51.08
3102
Stress-tensor-Voigt (kbar): 801.74 662.12 642.18 78.03 22.84 51.08
2770
3103
(Free)E + p*V (eV/cell) -2985.9298
2771
3104
Target enthalpy (eV/cell) -2940.3648
2815
3148
Folding of H and S is implicitly performed
2816
3149
<dSpData1D:S at geom step 21
2817
3150
<sparsity:sparsity for geom step 21
2818
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
2819
<dData1D:(new from dSpData1D) n=306, refcount: 1>
3151
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
3152
<dData1D:(new from dSpData1D) n=204, refcount: 1>
2822
3155
Initializing Density Matrix...
2823
DM after filling with atomic data:
2824
<dSpData2D:(DM initialized from atoms)
3156
DM filled with atomic data:
3157
<dSpData2D:DM initialized from atoms
2825
3158
<sparsity:sparsity for geom step 21
2826
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
2827
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
3159
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
3160
<dData2D:DM n=204 m=1, refcount: 1>
2829
3163
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.65270095
3164
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2830
3165
New grid distribution: 1
3173
8 1: 12 7: 12 10: 12
2836
3175
InitMesh: MESH = 24 x 24 x 24 = 13824
2837
3176
InitMesh: (bp) = 12 x 12 x 12 = 1728
2838
3177
InitMesh: Mesh cutoff (required, used) = 50.000 63.833 Ry
2839
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
3178
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
2840
3179
New grid distribution: 2
2845
3188
New grid distribution: 3
2850
3197
Setting up quadratic distribution...
2851
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2852
PhiOnMesh: Number of (b)points on node 0 = 432
2853
PhiOnMesh: nlist on node 0 = 10188
2855
siesta: Eharris(eV) = -2968.893522
2858
siesta: Eharris(eV) = -2968.893522
3198
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
3199
PhiOnMesh: Number of (b)points on node 0 = 252
3200
PhiOnMesh: nlist on node 0 = 5712
3202
siesta: Eharris(eV) = -2968.893523
3205
siesta: Eharris(eV) = -2968.893523
3207
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
3208
Geom step, scf iteration, dmax: 21 1 1.797320
2861
3210
Using DM_out to compute the final energy and forces
3211
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2863
3213
siesta: Atomic forces (eV/Ang):
2864
3214
1 0.000000 -0.000000 0.000000
2865
2 -0.000000 0.000000 0.000000
2866
3 -0.000000 -0.000000 -0.000000
2867
4 -0.000000 -0.000000 -0.000000
2868
5 1.037056 -2.319508 0.000000
2869
6 -0.000000 1.391589 -2.125814
2870
7 -1.037056 0.927919 2.125814
2871
8 -1.037056 2.319508 0.000000
2872
9 0.000000 -1.391589 2.125814
2873
10 1.037056 -0.927919 -2.125814
2874
----------------------------------------
2875
Tot -0.000000 -0.000000 -0.000000
2876
----------------------------------------
3215
2 0.000000 -0.000000 0.000000
3216
3 0.000000 0.000000 -0.000000
3217
4 -0.000000 -0.000000 0.000000
3218
5 1.037055 -2.319509 0.000000
3219
6 0.000000 1.391590 -2.125814
3220
7 -1.037056 0.927920 2.125814
3221
8 -1.037055 2.319509 -0.000000
3222
9 -0.000000 -1.391590 2.125814
3223
10 1.037056 -0.927920 -2.125814
3224
----------------------------------------
3225
Tot -0.000000 -0.000000 0.000000
3226
----------------------------------------
2878
3228
Res 1.136275 sqrt( Sum f_i^2 / 3N )
2879
3229
----------------------------------------
2880
Max 2.319508 constrained
3230
Max 2.319509 constrained
2882
3232
Stress-tensor-Voigt (kbar): 809.58 678.10 659.32 73.47 21.50 48.10
2883
3233
(Free)E + p*V (eV/cell) -2986.5395
2928
3278
Folding of H and S is implicitly performed
2929
3279
<dSpData1D:S at geom step 22
2930
3280
<sparsity:sparsity for geom step 22
2931
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
2932
<dData1D:(new from dSpData1D) n=306, refcount: 1>
3281
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
3282
<dData1D:(new from dSpData1D) n=204, refcount: 1>
2935
3285
Initializing Density Matrix...
2936
DM after filling with atomic data:
2937
<dSpData2D:(DM initialized from atoms)
3286
DM filled with atomic data:
3287
<dSpData2D:DM initialized from atoms
2938
3288
<sparsity:sparsity for geom step 22
2939
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
2940
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
3289
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
3290
<dData2D:DM n=204 m=1, refcount: 1>
2942
3293
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.65981711
3294
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2943
3295
New grid distribution: 1
3303
8 1: 12 7: 12 10: 12
2949
3305
InitMesh: MESH = 24 x 24 x 24 = 13824
2950
3306
InitMesh: (bp) = 12 x 12 x 12 = 1728
2951
3307
InitMesh: Mesh cutoff (required, used) = 50.000 63.868 Ry
2952
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
3308
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
2953
3309
New grid distribution: 2
2958
3318
New grid distribution: 3
2963
3327
Setting up quadratic distribution...
2964
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
2965
PhiOnMesh: Number of (b)points on node 0 = 432
2966
PhiOnMesh: nlist on node 0 = 10201
2968
siesta: Eharris(eV) = -2969.013599
2971
siesta: Eharris(eV) = -2969.013599
3328
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
3329
PhiOnMesh: Number of (b)points on node 0 = 252
3330
PhiOnMesh: nlist on node 0 = 5728
3332
siesta: Eharris(eV) = -2969.013596
3335
siesta: Eharris(eV) = -2969.013596
3337
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
3338
Geom step, scf iteration, dmax: 22 1 1.797676
2974
3340
Using DM_out to compute the final energy and forces
3341
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
2976
3343
siesta: Atomic forces (eV/Ang):
2977
3344
1 0.000000 -0.000000 0.000000
2978
2 -0.000000 0.000000 0.000000
2979
3 -0.000000 -0.000000 -0.000000
2980
4 -0.000000 -0.000000 -0.000000
2981
5 0.872567 -1.951608 0.000000
2982
6 -0.000000 1.170867 -1.788636
2983
7 -0.872567 0.780741 1.788636
3345
2 0.000000 0.000000 0.000000
3346
3 -0.000000 -0.000001 0.000001
3347
4 0.000000 0.000001 -0.000001
3348
5 0.872567 -1.951608 -0.000000
3349
6 0.000001 1.170867 -1.788636
3350
7 -0.872567 0.780741 1.788635
2984
3351
8 -0.872567 1.951608 0.000000
2985
9 -0.000000 -1.170867 1.788636
2986
10 0.872567 -0.780741 -1.788636
3352
9 -0.000001 -1.170867 1.788636
3353
10 0.872567 -0.780741 -1.788635
2987
3354
----------------------------------------
2988
Tot -0.000000 -0.000000 -0.000000
3355
Tot -0.000000 -0.000000 0.000000
2989
3356
----------------------------------------
2991
3358
Res 0.956049 sqrt( Sum f_i^2 / 3N )
3041
3408
Folding of H and S is implicitly performed
3042
3409
<dSpData1D:S at geom step 23
3043
3410
<sparsity:sparsity for geom step 23
3044
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
3045
<dData1D:(new from dSpData1D) n=306, refcount: 1>
3411
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
3412
<dData1D:(new from dSpData1D) n=204, refcount: 1>
3048
3415
Initializing Density Matrix...
3049
DM after filling with atomic data:
3050
<dSpData2D:(DM initialized from atoms)
3416
DM filled with atomic data:
3417
<dSpData2D:DM initialized from atoms
3051
3418
<sparsity:sparsity for geom step 23
3052
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
3053
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
3419
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
3420
<dData2D:DM n=204 m=1, refcount: 1>
3055
3423
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.66742475
3424
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3056
3425
New grid distribution: 1
3433
8 1: 12 7: 12 10: 12
3062
3435
InitMesh: MESH = 24 x 24 x 24 = 13824
3063
3436
InitMesh: (bp) = 12 x 12 x 12 = 1728
3064
3437
InitMesh: Mesh cutoff (required, used) = 50.000 63.905 Ry
3065
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
3438
ExtMesh (bp) on 0 = 56 x 50 x 47 = 131600
3066
3439
New grid distribution: 2
3071
3448
New grid distribution: 3
3076
3457
Setting up quadratic distribution...
3077
ExtMesh (bp) on 0 = 56 x 50 x 50 = 140000
3078
PhiOnMesh: Number of (b)points on node 0 = 432
3079
PhiOnMesh: nlist on node 0 = 10247
3081
siesta: Eharris(eV) = -2969.124156
3084
siesta: Eharris(eV) = -2969.124156
3458
ExtMesh (bp) on 0 = 51 x 50 x 50 = 127500
3459
PhiOnMesh: Number of (b)points on node 0 = 252
3460
PhiOnMesh: nlist on node 0 = 5750
3462
siesta: Eharris(eV) = -2969.124159
3465
siesta: Eharris(eV) = -2969.124159
3467
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
3468
Geom step, scf iteration, dmax: 23 1 1.798071
3087
3470
Using DM_out to compute the final energy and forces
3471
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3089
3473
siesta: Atomic forces (eV/Ang):
3090
1 0.000000 -0.000000 0.000000
3091
2 -0.000000 0.000000 0.000000
3092
3 -0.000000 -0.000000 -0.000000
3093
4 -0.000000 -0.000000 -0.000000
3094
5 0.715692 -1.600736 0.000000
3095
6 -0.000000 0.960362 -1.467064
3474
1 -0.000000 -0.000000 0.000000
3475
2 0.000000 -0.000000 0.000000
3476
3 -0.000000 -0.000000 0.000000
3477
4 0.000000 0.000000 -0.000000
3478
5 0.715692 -1.600736 -0.000000
3479
6 0.000000 0.960362 -1.467065
3096
3480
7 -0.715692 0.640375 1.467064
3097
3481
8 -0.715692 1.600736 0.000000
3098
9 -0.000000 -0.960362 1.467064
3482
9 -0.000000 -0.960362 1.467065
3099
3483
10 0.715692 -0.640375 -1.467064
3100
3484
----------------------------------------
3101
Tot -0.000000 -0.000000 -0.000000
3485
Tot -0.000000 -0.000000 0.000000
3102
3486
----------------------------------------
3104
3488
Res 0.784165 sqrt( Sum f_i^2 / 3N )
3154
3538
Folding of H and S is implicitly performed
3155
3539
<dSpData1D:S at geom step 24
3156
3540
<sparsity:sparsity for geom step 24
3157
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
3158
<dData1D:(new from dSpData1D) n=306, refcount: 1>
3541
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
3542
<dData1D:(new from dSpData1D) n=204, refcount: 1>
3161
3545
Initializing Density Matrix...
3162
DM after filling with atomic data:
3163
<dSpData2D:(DM initialized from atoms)
3546
DM filled with atomic data:
3547
<dSpData2D:DM initialized from atoms
3164
3548
<sparsity:sparsity for geom step 24
3165
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
3166
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
3549
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
3550
<dData2D:DM n=204 m=1, refcount: 1>
3168
3553
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.67554600
3554
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3169
3555
New grid distribution: 1
3563
8 1: 12 7: 12 10: 12
3175
3565
InitMesh: MESH = 24 x 24 x 24 = 13824
3176
3566
InitMesh: (bp) = 12 x 12 x 12 = 1728
3177
3567
InitMesh: Mesh cutoff (required, used) = 50.000 63.944 Ry
3178
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3568
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
3179
3569
New grid distribution: 2
3184
3578
New grid distribution: 3
3189
3587
Setting up quadratic distribution...
3190
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3191
PhiOnMesh: Number of (b)points on node 0 = 432
3192
PhiOnMesh: nlist on node 0 = 10271
3194
siesta: Eharris(eV) = -2969.224857
3197
siesta: Eharris(eV) = -2969.224857
3588
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
3589
PhiOnMesh: Number of (b)points on node 0 = 216
3590
PhiOnMesh: nlist on node 0 = 4897
3592
siesta: Eharris(eV) = -2969.224859
3595
siesta: Eharris(eV) = -2969.224859
3597
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
3598
Geom step, scf iteration, dmax: 24 1 1.798506
3200
3600
Using DM_out to compute the final energy and forces
3601
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3202
3603
siesta: Atomic forces (eV/Ang):
3203
3604
1 0.000000 -0.000000 0.000000
3204
3605
2 -0.000000 0.000000 0.000000
3205
3 -0.000000 -0.000000 -0.000000
3606
3 0.000000 0.000000 0.000000
3206
3607
4 -0.000000 -0.000000 -0.000000
3207
5 0.544764 -1.218434 0.000000
3208
6 -0.000000 0.731000 -1.116687
3209
7 -0.544764 0.487435 1.116687
3210
8 -0.544764 1.218434 0.000000
3211
9 -0.000000 -0.731000 1.116687
3212
10 0.544764 -0.487435 -1.116687
3608
5 0.544764 -1.218434 0.000001
3609
6 -0.000001 0.730999 -1.116688
3610
7 -0.544763 0.487435 1.116687
3611
8 -0.544764 1.218434 -0.000001
3612
9 0.000001 -0.730999 1.116688
3613
10 0.544763 -0.487435 -1.116687
3213
3614
----------------------------------------
3214
Tot -0.000000 -0.000000 -0.000000
3615
Tot -0.000000 -0.000000 0.000000
3215
3616
----------------------------------------
3217
Res 0.596884 sqrt( Sum f_i^2 / 3N )
3618
Res 0.596883 sqrt( Sum f_i^2 / 3N )
3218
3619
----------------------------------------
3219
3620
Max 1.218434 constrained
3249
3650
0.50000000 0.50000000 0.50000000 1 2 Mg
3250
3651
0.25000000 0.25000000 0.25000000 2 3 C
3251
3652
-0.25000000 -0.25000000 -0.25000000 2 4 C
3252
0.54839186 -0.04839186 0.25000000 3 5 O
3253
0.25000000 0.54839186 -0.04839186 3 6 O
3254
-0.04839186 0.25000000 0.54839186 3 7 O
3255
-0.54839186 0.04839186 -0.25000000 3 8 O
3256
-0.25000000 -0.54839186 0.04839186 3 9 O
3257
0.04839186 -0.25000000 -0.54839186 3 10 O
3653
0.54839187 -0.04839186 0.25000000 3 5 O
3654
0.25000000 0.54839187 -0.04839186 3 6 O
3655
-0.04839187 0.25000000 0.54839186 3 7 O
3656
-0.54839187 0.04839186 -0.25000000 3 8 O
3657
-0.25000000 -0.54839187 0.04839186 3 9 O
3658
0.04839187 -0.25000000 -0.54839186 3 10 O
3259
3660
outcell: Unit cell vectors (Ang):
3260
3661
5.639103 0.025647 0.016789
3267
3668
Folding of H and S is implicitly performed
3268
3669
<dSpData1D:S at geom step 25
3269
3670
<sparsity:sparsity for geom step 25
3270
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
3271
<dData1D:(new from dSpData1D) n=306, refcount: 1>
3671
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
3672
<dData1D:(new from dSpData1D) n=204, refcount: 1>
3274
3675
Initializing Density Matrix...
3275
DM after filling with atomic data:
3276
<dSpData2D:(DM initialized from atoms)
3676
DM filled with atomic data:
3677
<dSpData2D:DM initialized from atoms
3277
3678
<sparsity:sparsity for geom step 25
3278
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
3279
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
3679
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
3680
<dData2D:DM n=204 m=1, refcount: 1>
3281
3683
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.68420591
3684
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3282
3685
New grid distribution: 1
3693
8 1: 12 7: 12 10: 12
3288
3695
InitMesh: MESH = 24 x 24 x 24 = 13824
3289
3696
InitMesh: (bp) = 12 x 12 x 12 = 1728
3290
3697
InitMesh: Mesh cutoff (required, used) = 50.000 63.985 Ry
3291
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3698
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
3292
3699
New grid distribution: 2
3297
3708
New grid distribution: 3
3302
3717
Setting up quadratic distribution...
3303
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3304
PhiOnMesh: Number of (b)points on node 0 = 432
3305
PhiOnMesh: nlist on node 0 = 10297
3307
siesta: Eharris(eV) = -2969.315981
3310
siesta: Eharris(eV) = -2969.315981
3718
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
3719
PhiOnMesh: Number of (b)points on node 0 = 216
3720
PhiOnMesh: nlist on node 0 = 4911
3722
siesta: Eharris(eV) = -2969.315980
3725
siesta: Eharris(eV) = -2969.315980
3727
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
3728
Geom step, scf iteration, dmax: 25 1 1.799014
3313
3730
Using DM_out to compute the final energy and forces
3731
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3315
3733
siesta: Atomic forces (eV/Ang):
3316
3734
1 0.000000 -0.000000 0.000000
3317
2 -0.000000 0.000000 0.000000
3318
3 -0.000000 -0.000000 -0.000000
3319
4 -0.000000 -0.000000 -0.000000
3320
5 0.408490 -0.913641 0.000000
3321
6 -0.000000 0.548139 -0.837346
3735
2 0.000000 -0.000000 0.000000
3736
3 -0.000000 0.000001 -0.000000
3737
4 0.000000 -0.000001 0.000000
3738
5 0.408490 -0.913639 -0.000000
3739
6 0.000001 0.548137 -0.837347
3322
3740
7 -0.408490 0.365502 0.837346
3323
8 -0.408490 0.913641 0.000000
3324
9 -0.000000 -0.548139 0.837346
3741
8 -0.408490 0.913639 0.000000
3742
9 -0.000001 -0.548137 0.837347
3325
3743
10 0.408490 -0.365502 -0.837346
3326
3744
----------------------------------------
3327
Tot -0.000000 -0.000000 -0.000000
3745
Tot 0.000000 -0.000000 -0.000000
3328
3746
----------------------------------------
3330
3748
Res 0.447572 sqrt( Sum f_i^2 / 3N )
3331
3749
----------------------------------------
3332
Max 0.913641 constrained
3750
Max 0.913639 constrained
3334
3752
Stress-tensor-Voigt (kbar): 839.74 742.85 729.01 54.14 15.85 35.44
3335
3753
(Free)E + p*V (eV/cell) -2988.6990
3380
3798
Folding of H and S is implicitly performed
3381
3799
<dSpData1D:S at geom step 26
3382
3800
<sparsity:sparsity for geom step 26
3383
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
3384
<dData1D:(new from dSpData1D) n=306, refcount: 1>
3801
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
3802
<dData1D:(new from dSpData1D) n=204, refcount: 1>
3387
3805
Initializing Density Matrix...
3388
DM after filling with atomic data:
3389
<dSpData2D:(DM initialized from atoms)
3806
DM filled with atomic data:
3807
<dSpData2D:DM initialized from atoms
3390
3808
<sparsity:sparsity for geom step 26
3391
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
3392
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
3809
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
3810
<dData2D:DM n=204 m=1, refcount: 1>
3394
3813
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.69342822
3814
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3395
3815
New grid distribution: 1
3823
8 1: 12 7: 12 10: 12
3401
3825
InitMesh: MESH = 24 x 24 x 24 = 13824
3402
3826
InitMesh: (bp) = 12 x 12 x 12 = 1728
3403
3827
InitMesh: Mesh cutoff (required, used) = 50.000 64.027 Ry
3404
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3828
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
3405
3829
New grid distribution: 2
3410
3838
New grid distribution: 3
3415
3847
Setting up quadratic distribution...
3416
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3417
PhiOnMesh: Number of (b)points on node 0 = 432
3418
PhiOnMesh: nlist on node 0 = 10342
3420
siesta: Eharris(eV) = -2969.398631
3423
siesta: Eharris(eV) = -2969.398631
3848
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
3849
PhiOnMesh: Number of (b)points on node 0 = 216
3850
PhiOnMesh: nlist on node 0 = 4938
3852
siesta: Eharris(eV) = -2969.398629
3855
siesta: Eharris(eV) = -2969.398629
3857
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
3858
Geom step, scf iteration, dmax: 26 1 1.799602
3426
3860
Using DM_out to compute the final energy and forces
3861
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3428
3863
siesta: Atomic forces (eV/Ang):
3429
3864
1 0.000000 -0.000000 0.000000
3430
2 -0.000000 0.000000 0.000000
3431
3 -0.000000 -0.000000 -0.000000
3432
4 -0.000000 -0.000000 -0.000000
3433
5 0.262386 -0.586860 0.000000
3434
6 -0.000000 0.352087 -0.537853
3435
7 -0.262386 0.234773 0.537853
3436
8 -0.262386 0.586860 0.000000
3437
9 -0.000000 -0.352087 0.537853
3438
10 0.262386 -0.234773 -0.537853
3865
2 0.000000 0.000000 0.000000
3866
3 -0.000000 0.000000 -0.000000
3867
4 0.000000 -0.000000 0.000000
3868
5 0.262385 -0.586860 0.000001
3869
6 0.000001 0.352087 -0.537853
3870
7 -0.262386 0.234772 0.537852
3871
8 -0.262385 0.586860 -0.000001
3872
9 -0.000001 -0.352087 0.537853
3873
10 0.262386 -0.234772 -0.537852
3439
3874
----------------------------------------
3440
3875
Tot -0.000000 -0.000000 -0.000000
3441
3876
----------------------------------------
3443
Res 0.287490 sqrt( Sum f_i^2 / 3N )
3878
Res 0.287489 sqrt( Sum f_i^2 / 3N )
3444
3879
----------------------------------------
3445
3880
Max 0.586860 constrained
3493
3928
Folding of H and S is implicitly performed
3494
3929
<dSpData1D:S at geom step 27
3495
3930
<sparsity:sparsity for geom step 27
3496
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
3497
<dData1D:(new from dSpData1D) n=306, refcount: 1>
3931
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
3932
<dData1D:(new from dSpData1D) n=204, refcount: 1>
3500
3935
Initializing Density Matrix...
3501
DM after filling with atomic data:
3502
<dSpData2D:(DM initialized from atoms)
3936
DM filled with atomic data:
3937
<dSpData2D:DM initialized from atoms
3503
3938
<sparsity:sparsity for geom step 27
3504
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
3505
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
3939
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
3940
<dData2D:DM n=204 m=1, refcount: 1>
3507
3943
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.70324236
3944
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3508
3945
New grid distribution: 1
3953
8 1: 12 7: 12 10: 12
3514
3955
InitMesh: MESH = 24 x 24 x 24 = 13824
3515
3956
InitMesh: (bp) = 12 x 12 x 12 = 1728
3516
3957
InitMesh: Mesh cutoff (required, used) = 50.000 64.072 Ry
3517
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3958
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
3518
3959
New grid distribution: 2
3523
3968
New grid distribution: 3
3528
3977
Setting up quadratic distribution...
3529
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3530
PhiOnMesh: Number of (b)points on node 0 = 432
3531
PhiOnMesh: nlist on node 0 = 10372
3533
siesta: Eharris(eV) = -2969.473973
3536
siesta: Eharris(eV) = -2969.473973
3978
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
3979
PhiOnMesh: Number of (b)points on node 0 = 216
3980
PhiOnMesh: nlist on node 0 = 4959
3982
siesta: Eharris(eV) = -2969.473974
3985
siesta: Eharris(eV) = -2969.473974
3987
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
3988
Geom step, scf iteration, dmax: 27 1 1.800299
3539
3990
Using DM_out to compute the final energy and forces
3991
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3541
3993
siesta: Atomic forces (eV/Ang):
3542
3994
1 0.000000 -0.000000 0.000000
3543
2 -0.000000 0.000000 0.000000
3544
3 -0.000000 -0.000000 -0.000000
3545
4 -0.000000 -0.000000 -0.000000
3546
5 0.153372 -0.343036 0.000000
3547
6 -0.000000 0.205804 -0.314390
3548
7 -0.153372 0.137231 0.314390
3549
8 -0.153372 0.343036 0.000000
3550
9 -0.000000 -0.205804 0.314390
3551
10 0.153372 -0.137231 -0.314390
3995
2 0.000000 0.000000 0.000000
3996
3 0.000000 0.000000 0.000001
3997
4 -0.000000 -0.000000 -0.000001
3998
5 0.153371 -0.343036 -0.000000
3999
6 -0.000000 0.205803 -0.314389
4000
7 -0.153371 0.137232 0.314389
4001
8 -0.153371 0.343036 0.000000
4002
9 0.000000 -0.205803 0.314389
4003
10 0.153371 -0.137232 -0.314389
3552
4004
----------------------------------------
3553
4005
Tot -0.000000 -0.000000 -0.000000
3554
4006
----------------------------------------
3556
Res 0.168046 sqrt( Sum f_i^2 / 3N )
4008
Res 0.168045 sqrt( Sum f_i^2 / 3N )
3557
4009
----------------------------------------
3558
4010
Max 0.343036 constrained
3606
4058
Folding of H and S is implicitly performed
3607
4059
<dSpData1D:S at geom step 28
3608
4060
<sparsity:sparsity for geom step 28
3609
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
3610
<dData1D:(new from dSpData1D) n=306, refcount: 1>
4061
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
4062
<dData1D:(new from dSpData1D) n=204, refcount: 1>
3613
4065
Initializing Density Matrix...
3614
DM after filling with atomic data:
3615
<dSpData2D:(DM initialized from atoms)
4066
DM filled with atomic data:
4067
<dSpData2D:DM initialized from atoms
3616
4068
<sparsity:sparsity for geom step 28
3617
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
3618
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
4069
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
4070
<dData2D:DM n=204 m=1, refcount: 1>
3620
4073
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.71367330
4074
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3621
4075
New grid distribution: 1
4083
8 1: 12 7: 12 10: 12
3627
4085
InitMesh: MESH = 24 x 24 x 24 = 13824
3628
4086
InitMesh: (bp) = 12 x 12 x 12 = 1728
3629
4087
InitMesh: Mesh cutoff (required, used) = 50.000 64.118 Ry
3630
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4088
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
3631
4089
New grid distribution: 2
3636
4098
New grid distribution: 3
3641
4107
Setting up quadratic distribution...
3642
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3643
PhiOnMesh: Number of (b)points on node 0 = 432
3644
PhiOnMesh: nlist on node 0 = 10419
3646
siesta: Eharris(eV) = -2969.544047
3649
siesta: Eharris(eV) = -2969.544047
4108
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
4109
PhiOnMesh: Number of (b)points on node 0 = 216
4110
PhiOnMesh: nlist on node 0 = 4991
4112
siesta: Eharris(eV) = -2969.544043
4115
siesta: Eharris(eV) = -2969.544043
4117
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
4118
Geom step, scf iteration, dmax: 28 1 1.801110
3652
4120
Using DM_out to compute the final energy and forces
4121
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3654
4123
siesta: Atomic forces (eV/Ang):
3655
1 0.000000 -0.000000 0.000000
3656
2 -0.000000 0.000000 0.000000
3657
3 -0.000000 -0.000000 -0.000000
3658
4 -0.000000 -0.000000 -0.000000
3659
5 0.057364 -0.128301 0.000000
3660
6 -0.000000 0.076974 -0.117587
3661
7 -0.057364 0.051327 0.117587
3662
8 -0.057364 0.128301 0.000000
4124
1 -0.000000 0.000000 -0.000000
4125
2 0.000000 -0.000000 -0.000000
4126
3 0.000000 0.000000 -0.000001
4127
4 -0.000000 -0.000000 0.000001
4128
5 0.057363 -0.128301 0.000000
4129
6 0.000000 0.076974 -0.117587
4130
7 -0.057363 0.051327 0.117587
4131
8 -0.057363 0.128301 -0.000000
3663
4132
9 -0.000000 -0.076974 0.117587
3664
10 0.057364 -0.051327 -0.117587
4133
10 0.057363 -0.051327 -0.117587
3665
4134
----------------------------------------
3666
4135
Tot -0.000000 -0.000000 -0.000000
3667
4136
----------------------------------------
3719
4188
Folding of H and S is implicitly performed
3720
4189
<dSpData1D:S at geom step 29
3721
4190
<sparsity:sparsity for geom step 29
3722
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
3723
<dData1D:(new from dSpData1D) n=306, refcount: 1>
4191
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
4192
<dData1D:(new from dSpData1D) n=204, refcount: 1>
3726
4195
Initializing Density Matrix...
3727
DM after filling with atomic data:
3728
<dSpData2D:(DM initialized from atoms)
4196
DM filled with atomic data:
4197
<dSpData2D:DM initialized from atoms
3729
4198
<sparsity:sparsity for geom step 29
3730
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
3731
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
4199
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
4200
<dData2D:DM n=204 m=1, refcount: 1>
3733
4203
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.72474934
4204
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3734
4205
New grid distribution: 1
4213
8 1: 12 7: 12 10: 12
3740
4215
InitMesh: MESH = 24 x 24 x 24 = 13824
3741
4216
InitMesh: (bp) = 12 x 12 x 12 = 1728
3742
4217
InitMesh: Mesh cutoff (required, used) = 50.000 64.165 Ry
3743
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4218
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
3744
4219
New grid distribution: 2
3749
4228
New grid distribution: 3
3754
4237
Setting up quadratic distribution...
3755
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3756
PhiOnMesh: Number of (b)points on node 0 = 432
3757
PhiOnMesh: nlist on node 0 = 10440
3759
siesta: Eharris(eV) = -2969.609964
3762
siesta: Eharris(eV) = -2969.609964
4238
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
4239
PhiOnMesh: Number of (b)points on node 0 = 216
4240
PhiOnMesh: nlist on node 0 = 5021
4242
siesta: Eharris(eV) = -2969.609967
4245
siesta: Eharris(eV) = -2969.609967
4247
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
4248
Geom step, scf iteration, dmax: 29 1 1.802021
3765
4250
Using DM_out to compute the final energy and forces
4251
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3767
4253
siesta: Atomic forces (eV/Ang):
3768
1 0.000000 -0.000000 0.000000
3769
2 -0.000000 0.000000 0.000000
3770
3 -0.000000 -0.000000 -0.000000
3771
4 -0.000000 -0.000000 -0.000000
4254
1 -0.000000 0.000000 -0.000000
4255
2 0.000000 -0.000000 0.000000
4256
3 -0.000000 0.000000 -0.000000
4257
4 0.000000 -0.000000 0.000000
3772
4258
5 -0.022971 0.051378 0.000000
3773
6 -0.000000 -0.030824 0.047088
3774
7 0.022971 -0.020554 -0.047088
3775
8 0.022971 -0.051378 0.000000
3776
9 -0.000000 0.030824 -0.047088
3777
10 -0.022971 0.020554 0.047088
4259
6 0.000000 -0.030824 0.047087
4260
7 0.022971 -0.020553 -0.047087
4261
8 0.022971 -0.051378 -0.000000
4262
9 -0.000000 0.030824 -0.047087
4263
10 -0.022971 0.020553 0.047087
3778
4264
----------------------------------------
3779
4265
Tot -0.000000 -0.000000 -0.000000
3780
4266
----------------------------------------
3832
4318
Folding of H and S is implicitly performed
3833
4319
<dSpData1D:S at geom step 30
3834
4320
<sparsity:sparsity for geom step 30
3835
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
3836
<dData1D:(new from dSpData1D) n=306, refcount: 1>
4321
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
4322
<dData1D:(new from dSpData1D) n=204, refcount: 1>
3839
4325
Initializing Density Matrix...
3840
DM after filling with atomic data:
3841
<dSpData2D:(DM initialized from atoms)
4326
DM filled with atomic data:
4327
<dSpData2D:DM initialized from atoms
3842
4328
<sparsity:sparsity for geom step 30
3843
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
3844
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
4329
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
4330
<dData2D:DM n=204 m=1, refcount: 1>
3846
4333
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.73649927
4334
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3847
4335
New grid distribution: 1
4343
8 1: 12 7: 12 10: 12
3853
4345
InitMesh: MESH = 24 x 24 x 24 = 13824
3854
4346
InitMesh: (bp) = 12 x 12 x 12 = 1728
3855
4347
InitMesh: Mesh cutoff (required, used) = 50.000 64.215 Ry
3856
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4348
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
3857
4349
New grid distribution: 2
3862
4358
New grid distribution: 3
3867
4367
Setting up quadratic distribution...
3868
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3869
PhiOnMesh: Number of (b)points on node 0 = 432
3870
PhiOnMesh: nlist on node 0 = 10447
3872
siesta: Eharris(eV) = -2969.674009
3875
siesta: Eharris(eV) = -2969.674009
4368
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
4369
PhiOnMesh: Number of (b)points on node 0 = 216
4370
PhiOnMesh: nlist on node 0 = 5037
4372
siesta: Eharris(eV) = -2969.674013
4375
siesta: Eharris(eV) = -2969.674013
4377
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
4378
Geom step, scf iteration, dmax: 30 1 1.803047
3878
4380
Using DM_out to compute the final energy and forces
4381
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3880
4383
siesta: Atomic forces (eV/Ang):
3881
1 0.000000 -0.000000 0.000000
4384
1 0.000000 0.000000 0.000000
3882
4385
2 -0.000000 0.000000 0.000000
3883
3 0.000000 -0.000000 -0.000000
3884
4 -0.000000 -0.000000 -0.000000
3885
5 -0.089244 0.199605 0.000000
3886
6 -0.000000 -0.119753 0.182937
3887
7 0.089244 -0.079852 -0.182937
3888
8 0.089244 -0.199605 0.000000
3889
9 -0.000000 0.119753 -0.182937
3890
10 -0.089244 0.079852 0.182937
4386
3 0.000000 -0.000000 0.000000
4387
4 -0.000000 0.000000 -0.000000
4388
5 -0.089244 0.199605 0.000001
4389
6 0.000001 -0.119753 0.182937
4390
7 0.089245 -0.079852 -0.182937
4391
8 0.089244 -0.199605 -0.000001
4392
9 -0.000001 0.119753 -0.182937
4393
10 -0.089245 0.079852 0.182937
3891
4394
----------------------------------------
3892
Tot -0.000000 -0.000000 -0.000000
4395
Tot 0.000000 -0.000000 -0.000000
3893
4396
----------------------------------------
3895
4398
Res 0.097782 sqrt( Sum f_i^2 / 3N )
3945
4448
Folding of H and S is implicitly performed
3946
4449
<dSpData1D:S at geom step 31
3947
4450
<sparsity:sparsity for geom step 31
3948
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
3949
<dData1D:(new from dSpData1D) n=306, refcount: 1>
4451
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
4452
<dData1D:(new from dSpData1D) n=204, refcount: 1>
3952
4455
Initializing Density Matrix...
3953
DM after filling with atomic data:
3954
<dSpData2D:(DM initialized from atoms)
4456
DM filled with atomic data:
4457
<dSpData2D:DM initialized from atoms
3955
4458
<sparsity:sparsity for geom step 31
3956
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
3957
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
4459
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
4460
<dData2D:DM n=204 m=1, refcount: 1>
3959
4463
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.74895443
4464
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3960
4465
New grid distribution: 1
4473
8 1: 12 7: 12 10: 12
3966
4475
InitMesh: MESH = 24 x 24 x 24 = 13824
3967
4476
InitMesh: (bp) = 12 x 12 x 12 = 1728
3968
4477
InitMesh: Mesh cutoff (required, used) = 50.000 64.266 Ry
3969
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4478
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
3970
4479
New grid distribution: 2
3975
4488
New grid distribution: 3
3980
4497
Setting up quadratic distribution...
3981
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
3982
PhiOnMesh: Number of (b)points on node 0 = 432
3983
PhiOnMesh: nlist on node 0 = 10512
3985
siesta: Eharris(eV) = -2969.738085
3988
siesta: Eharris(eV) = -2969.738085
4498
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
4499
PhiOnMesh: Number of (b)points on node 0 = 216
4500
PhiOnMesh: nlist on node 0 = 5084
4502
siesta: Eharris(eV) = -2969.738081
4505
siesta: Eharris(eV) = -2969.738081
4507
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
4508
Geom step, scf iteration, dmax: 31 1 1.804175
3991
4510
Using DM_out to compute the final energy and forces
4511
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
3993
4513
siesta: Atomic forces (eV/Ang):
3994
4514
1 0.000000 -0.000000 0.000000
3995
2 -0.000000 0.000000 0.000000
3996
3 0.000000 -0.000000 -0.000000
3997
4 -0.000000 -0.000000 -0.000000
3998
5 -0.138557 0.309901 0.000000
3999
6 -0.000000 -0.185925 0.284022
4000
7 0.138557 -0.123976 -0.284022
4001
8 0.138557 -0.309901 0.000000
4002
9 -0.000000 0.185925 -0.284022
4003
10 -0.138557 0.123976 0.284022
4515
2 -0.000000 0.000000 -0.000000
4516
3 -0.000000 -0.000001 0.000001
4517
4 0.000000 0.000001 -0.000001
4518
5 -0.138557 0.309901 0.000001
4519
6 -0.000000 -0.185925 0.284021
4520
7 0.138559 -0.123976 -0.284022
4521
8 0.138557 -0.309901 -0.000001
4522
9 0.000000 0.185925 -0.284021
4523
10 -0.138559 0.123976 0.284022
4004
4524
----------------------------------------
4005
Tot -0.000000 -0.000000 -0.000000
4525
Tot 0.000000 -0.000000 -0.000000
4006
4526
----------------------------------------
4008
Res 0.151813 sqrt( Sum f_i^2 / 3N )
4528
Res 0.151814 sqrt( Sum f_i^2 / 3N )
4009
4529
----------------------------------------
4010
4530
Max 0.309901 constrained
4058
4578
Folding of H and S is implicitly performed
4059
4579
<dSpData1D:S at geom step 32
4060
4580
<sparsity:sparsity for geom step 32
4061
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
4062
<dData1D:(new from dSpData1D) n=306, refcount: 1>
4581
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
4582
<dData1D:(new from dSpData1D) n=204, refcount: 1>
4065
4585
Initializing Density Matrix...
4066
DM after filling with atomic data:
4067
<dSpData2D:(DM initialized from atoms)
4586
DM filled with atomic data:
4587
<dSpData2D:DM initialized from atoms
4068
4588
<sparsity:sparsity for geom step 32
4069
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
4070
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
4589
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
4590
<dData2D:DM n=204 m=1, refcount: 1>
4072
4593
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.76214823
4594
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4073
4595
New grid distribution: 1
4603
8 1: 12 7: 12 10: 12
4079
4605
InitMesh: MESH = 24 x 24 x 24 = 13824
4080
4606
InitMesh: (bp) = 12 x 12 x 12 = 1728
4081
4607
InitMesh: Mesh cutoff (required, used) = 50.000 64.320 Ry
4082
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4608
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
4083
4609
New grid distribution: 2
4088
4618
New grid distribution: 3
4093
4627
Setting up quadratic distribution...
4094
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4095
PhiOnMesh: Number of (b)points on node 0 = 432
4096
PhiOnMesh: nlist on node 0 = 10534
4098
siesta: Eharris(eV) = -2969.803142
4101
siesta: Eharris(eV) = -2969.803142
4628
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
4629
PhiOnMesh: Number of (b)points on node 0 = 216
4630
PhiOnMesh: nlist on node 0 = 5106
4632
siesta: Eharris(eV) = -2969.803140
4635
siesta: Eharris(eV) = -2969.803140
4637
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
4638
Geom step, scf iteration, dmax: 32 1 1.805438
4104
4640
Using DM_out to compute the final energy and forces
4641
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4106
4643
siesta: Atomic forces (eV/Ang):
4107
4644
1 0.000000 -0.000000 0.000000
4108
4645
2 -0.000000 0.000000 0.000000
4109
3 0.000000 -0.000000 -0.000000
4110
4 -0.000000 -0.000000 -0.000000
4111
5 -0.170742 0.381886 0.000000
4112
6 -0.000000 -0.229112 0.349996
4646
3 -0.000000 -0.000000 0.000001
4647
4 0.000000 0.000000 -0.000001
4648
5 -0.170741 0.381886 -0.000002
4649
6 -0.000000 -0.229111 0.349997
4113
4650
7 0.170742 -0.152773 -0.349996
4114
8 0.170742 -0.381886 0.000000
4115
9 -0.000000 0.229112 -0.349996
4651
8 0.170741 -0.381886 0.000002
4652
9 0.000000 0.229111 -0.349997
4116
4653
10 -0.170742 0.152773 0.349996
4117
4654
----------------------------------------
4118
Tot -0.000000 -0.000000 -0.000000
4655
Tot 0.000000 -0.000000 -0.000000
4119
4656
----------------------------------------
4121
4658
Res 0.187077 sqrt( Sum f_i^2 / 3N )
4131
4668
0.003327 0.563087 0.000974
4132
4669
0.002178 0.000974 0.562237
4134
siesta: Pressure (static): -904.90023346 kBar
4671
siesta: Pressure (static): -904.90032507 kBar
4136
4673
siesta: Stress tensor (total) (eV/Ang**3):
4137
4674
0.562621 0.013380 0.008758
4138
4675
0.013380 0.538678 0.003916
4139
0.008758 0.003916 0.535259
4141
siesta: Pressure (total): -874.02801154 kBar
4143
siesta: E_kin PR = 34329.488 K
4145
siesta: Temp_ion = 2186.786 K
4676
0.008758 0.003916 0.535260
4678
siesta: Pressure (total): -874.02810472 kBar
4680
siesta: E_kin PR = 34329.487 K
4682
siesta: Temp_ion = 2186.785 K
4147
4684
====================================
4148
4685
Begin MD step = 33
4171
4708
Folding of H and S is implicitly performed
4172
4709
<dSpData1D:S at geom step 33
4173
4710
<sparsity:sparsity for geom step 33
4174
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
4175
<dData1D:(new from dSpData1D) n=306, refcount: 1>
4711
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
4712
<dData1D:(new from dSpData1D) n=204, refcount: 1>
4178
4715
Initializing Density Matrix...
4179
DM after filling with atomic data:
4180
<dSpData2D:(DM initialized from atoms)
4716
DM filled with atomic data:
4717
<dSpData2D:DM initialized from atoms
4181
4718
<sparsity:sparsity for geom step 33
4182
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
4183
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
4719
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
4720
<dData2D:DM n=204 m=1, refcount: 1>
4185
4723
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.77611636
4724
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4186
4725
New grid distribution: 1
4733
8 1: 12 7: 12 10: 12
4192
4735
InitMesh: MESH = 24 x 24 x 24 = 13824
4193
4736
InitMesh: (bp) = 12 x 12 x 12 = 1728
4194
4737
InitMesh: Mesh cutoff (required, used) = 50.000 64.375 Ry
4195
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4738
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
4196
4739
New grid distribution: 2
4201
4748
New grid distribution: 3
4206
4757
Setting up quadratic distribution...
4207
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4208
PhiOnMesh: Number of (b)points on node 0 = 432
4209
PhiOnMesh: nlist on node 0 = 10562
4211
siesta: Eharris(eV) = -2969.871479
4214
siesta: Eharris(eV) = -2969.871479
4758
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
4759
PhiOnMesh: Number of (b)points on node 0 = 216
4760
PhiOnMesh: nlist on node 0 = 5128
4762
siesta: Eharris(eV) = -2969.871479
4765
siesta: Eharris(eV) = -2969.871479
4767
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
4768
Geom step, scf iteration, dmax: 33 1 1.806812
4217
4770
Using DM_out to compute the final energy and forces
4771
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4219
4773
siesta: Atomic forces (eV/Ang):
4220
4774
1 0.000000 -0.000000 0.000000
4221
4775
2 -0.000000 0.000000 0.000000
4222
3 0.000000 -0.000000 -0.000000
4223
4 -0.000000 -0.000000 -0.000000
4224
5 -0.221493 0.495398 0.000000
4225
6 -0.000000 -0.297214 0.454029
4226
7 0.221493 -0.198184 -0.454029
4227
8 0.221493 -0.495398 -0.000000
4228
9 -0.000000 0.297214 -0.454029
4229
10 -0.221493 0.198184 0.454029
4230
----------------------------------------
4231
Tot -0.000000 -0.000000 -0.000000
4232
----------------------------------------
4776
3 0.000000 0.000000 -0.000001
4777
4 -0.000000 -0.000000 0.000001
4778
5 -0.221494 0.495399 0.000002
4779
6 0.000000 -0.297213 0.454029
4780
7 0.221494 -0.198183 -0.454029
4781
8 0.221494 -0.495399 -0.000002
4782
9 -0.000000 0.297213 -0.454029
4783
10 -0.221494 0.198183 0.454029
4784
----------------------------------------
4785
Tot 0.000000 -0.000000 -0.000000
4786
----------------------------------------
4234
4788
Res 0.242684 sqrt( Sum f_i^2 / 3N )
4235
4789
----------------------------------------
4236
Max 0.495398 constrained
4790
Max 0.495399 constrained
4238
4792
Stress-tensor-Voigt (kbar): 914.11 888.55 884.90 14.28 4.18 9.35
4239
4793
(Free)E + p*V (eV/cell) -2992.3336
4284
4838
Folding of H and S is implicitly performed
4285
4839
<dSpData1D:S at geom step 34
4286
4840
<sparsity:sparsity for geom step 34
4287
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
4288
<dData1D:(new from dSpData1D) n=306, refcount: 1>
4841
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
4842
<dData1D:(new from dSpData1D) n=204, refcount: 1>
4291
4845
Initializing Density Matrix...
4292
DM after filling with atomic data:
4293
<dSpData2D:(DM initialized from atoms)
4846
DM filled with atomic data:
4847
<dSpData2D:DM initialized from atoms
4294
4848
<sparsity:sparsity for geom step 34
4295
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
4296
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
4849
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
4850
<dData2D:DM n=204 m=1, refcount: 1>
4298
4853
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.79090171
4854
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4299
4855
New grid distribution: 1
4863
8 1: 12 7: 12 10: 12
4305
4865
InitMesh: MESH = 24 x 24 x 24 = 13824
4306
4866
InitMesh: (bp) = 12 x 12 x 12 = 1728
4307
4867
InitMesh: Mesh cutoff (required, used) = 50.000 64.432 Ry
4308
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4868
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
4309
4869
New grid distribution: 2
4314
4878
New grid distribution: 3
4319
4887
Setting up quadratic distribution...
4320
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4321
PhiOnMesh: Number of (b)points on node 0 = 432
4322
PhiOnMesh: nlist on node 0 = 10605
4324
siesta: Eharris(eV) = -2969.944549
4327
siesta: Eharris(eV) = -2969.944549
4888
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
4889
PhiOnMesh: Number of (b)points on node 0 = 216
4890
PhiOnMesh: nlist on node 0 = 5147
4892
siesta: Eharris(eV) = -2969.944550
4895
siesta: Eharris(eV) = -2969.944550
4897
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
4898
Geom step, scf iteration, dmax: 34 1 1.808208
4330
4900
Using DM_out to compute the final energy and forces
4901
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4332
4903
siesta: Atomic forces (eV/Ang):
4333
1 0.000000 -0.000000 0.000000
4904
1 -0.000000 -0.000000 0.000000
4334
4905
2 -0.000000 0.000000 0.000000
4335
3 0.000000 -0.000000 -0.000000
4336
4 -0.000000 -0.000000 -0.000000
4337
5 -0.234366 0.524189 -0.000000
4338
6 -0.000000 -0.314487 0.480416
4339
7 0.234366 -0.209702 -0.480416
4340
8 0.234366 -0.524189 -0.000000
4341
9 -0.000000 0.314487 -0.480416
4342
10 -0.234366 0.209702 0.480416
4906
3 -0.000000 -0.000001 0.000001
4907
4 0.000000 0.000001 -0.000001
4908
5 -0.234365 0.524189 0.000000
4909
6 -0.000000 -0.314486 0.480416
4910
7 0.234364 -0.209701 -0.480416
4911
8 0.234365 -0.524189 -0.000000
4912
9 0.000000 0.314486 -0.480416
4913
10 -0.234364 0.209701 0.480416
4343
4914
----------------------------------------
4344
Tot -0.000000 -0.000000 -0.000000
4915
Tot 0.000000 -0.000000 -0.000000
4345
4916
----------------------------------------
4347
4918
Res 0.256788 sqrt( Sum f_i^2 / 3N )
4397
4968
Folding of H and S is implicitly performed
4398
4969
<dSpData1D:S at geom step 35
4399
4970
<sparsity:sparsity for geom step 35
4400
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
4401
<dData1D:(new from dSpData1D) n=306, refcount: 1>
4971
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
4972
<dData1D:(new from dSpData1D) n=204, refcount: 1>
4404
4975
Initializing Density Matrix...
4405
DM after filling with atomic data:
4406
<dSpData2D:(DM initialized from atoms)
4976
DM filled with atomic data:
4977
<dSpData2D:DM initialized from atoms
4407
4978
<sparsity:sparsity for geom step 35
4408
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
4409
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
4979
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
4980
<dData2D:DM n=204 m=1, refcount: 1>
4411
4983
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.80654555
4984
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4412
4985
New grid distribution: 1
4993
8 1: 12 7: 12 10: 12
4418
4995
InitMesh: MESH = 24 x 24 x 24 = 13824
4419
4996
InitMesh: (bp) = 12 x 12 x 12 = 1728
4420
4997
InitMesh: Mesh cutoff (required, used) = 50.000 64.490 Ry
4421
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4998
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
4422
4999
New grid distribution: 2
4427
5008
New grid distribution: 3
4432
5017
Setting up quadratic distribution...
4433
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4434
PhiOnMesh: Number of (b)points on node 0 = 432
4435
PhiOnMesh: nlist on node 0 = 10668
4437
siesta: Eharris(eV) = -2970.021453
4440
siesta: Eharris(eV) = -2970.021453
5018
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
5019
PhiOnMesh: Number of (b)points on node 0 = 216
5020
PhiOnMesh: nlist on node 0 = 5180
5022
siesta: Eharris(eV) = -2970.021450
5025
siesta: Eharris(eV) = -2970.021450
5027
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
5028
Geom step, scf iteration, dmax: 35 1 1.809743
4443
5030
Using DM_out to compute the final energy and forces
5031
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4445
5033
siesta: Atomic forces (eV/Ang):
4446
5034
1 0.000000 -0.000000 0.000000
4447
2 -0.000000 0.000000 0.000000
4448
3 0.000000 -0.000000 -0.000000
5035
2 0.000000 -0.000000 0.000000
5036
3 0.000000 0.000000 0.000000
4449
5037
4 -0.000000 -0.000000 -0.000000
4450
5 -0.256128 0.572863 -0.000000
4451
6 -0.000000 -0.343689 0.525025
4452
7 0.256128 -0.229174 -0.525025
4453
8 0.256128 -0.572863 -0.000000
4454
9 -0.000000 0.343689 -0.525025
4455
10 -0.256128 0.229174 0.525025
5038
5 -0.256128 0.572864 0.000000
5039
6 -0.000001 -0.343690 0.525025
5040
7 0.256128 -0.229174 -0.525026
5041
8 0.256128 -0.572864 -0.000000
5042
9 0.000001 0.343690 -0.525025
5043
10 -0.256128 0.229174 0.525026
4456
5044
----------------------------------------
4457
5045
Tot -0.000000 -0.000000 -0.000000
4458
5046
----------------------------------------
4460
5048
Res 0.280633 sqrt( Sum f_i^2 / 3N )
4461
5049
----------------------------------------
4462
Max 0.572863 constrained
5050
Max 0.572864 constrained
4464
5052
Stress-tensor-Voigt (kbar): 935.53 922.87 921.06 7.07 2.07 4.63
4465
5053
(Free)E + p*V (eV/cell) -2992.8760
4466
5054
Target enthalpy (eV/cell) -2937.8845
4468
5056
siesta: Stress tensor (static) (eV/Ang**3):
4469
0.590359 -0.005692 -0.003726
5057
0.590358 -0.005692 -0.003726
4470
5058
-0.005692 0.600545 -0.001666
4471
5059
-0.003726 -0.001666 0.602000
4473
siesta: Pressure (static): -957.52641833 kBar
5061
siesta: Pressure (static): -957.52638408 kBar
4475
5063
siesta: Stress tensor (total) (eV/Ang**3):
4476
5064
0.583903 0.004415 0.002890
4477
5065
0.004415 0.576002 0.001292
4478
5066
0.002890 0.001292 0.574873
4480
siesta: Pressure (total): -926.48355654 kBar
5068
siesta: Pressure (total): -926.48352393 kBar
4482
5070
siesta: E_kin PR = 43893.767 K
4510
5098
Folding of H and S is implicitly performed
4511
5099
<dSpData1D:S at geom step 36
4512
5100
<sparsity:sparsity for geom step 36
4513
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
4514
<dData1D:(new from dSpData1D) n=306, refcount: 1>
5101
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
5102
<dData1D:(new from dSpData1D) n=204, refcount: 1>
4517
5105
Initializing Density Matrix...
4518
DM after filling with atomic data:
4519
<dSpData2D:(DM initialized from atoms)
5106
DM filled with atomic data:
5107
<dSpData2D:DM initialized from atoms
4520
5108
<sparsity:sparsity for geom step 36
4521
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
4522
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
5109
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
5110
<dData2D:DM n=204 m=1, refcount: 1>
4524
5113
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.82308667
5114
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4525
5115
New grid distribution: 1
5123
8 1: 12 7: 12 10: 12
4531
5125
InitMesh: MESH = 24 x 24 x 24 = 13824
4532
5126
InitMesh: (bp) = 12 x 12 x 12 = 1728
4533
5127
InitMesh: Mesh cutoff (required, used) = 50.000 64.551 Ry
4534
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5128
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
4535
5129
New grid distribution: 2
4540
5138
New grid distribution: 3
4545
5147
Setting up quadratic distribution...
4546
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4547
PhiOnMesh: Number of (b)points on node 0 = 432
4548
PhiOnMesh: nlist on node 0 = 10735
4550
siesta: Eharris(eV) = -2970.101652
4553
siesta: Eharris(eV) = -2970.101652
5148
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
5149
PhiOnMesh: Number of (b)points on node 0 = 216
5150
PhiOnMesh: nlist on node 0 = 5211
5152
siesta: Eharris(eV) = -2970.101651
5155
siesta: Eharris(eV) = -2970.101651
5157
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
5158
Geom step, scf iteration, dmax: 36 1 1.811385
4556
5160
Using DM_out to compute the final energy and forces
5161
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4558
5163
siesta: Atomic forces (eV/Ang):
4559
1 0.000000 -0.000000 0.000000
4560
2 -0.000000 0.000000 0.000000
4561
3 0.000000 -0.000000 -0.000000
4562
4 -0.000000 -0.000000 -0.000000
4563
5 -0.267947 0.599298 -0.000000
4564
6 -0.000000 -0.359549 0.549253
4565
7 0.267947 -0.239749 -0.549253
4566
8 0.267947 -0.599298 -0.000000
4567
9 -0.000000 0.359549 -0.549253
4568
10 -0.267947 0.239749 0.549253
4569
----------------------------------------
4570
Tot -0.000000 -0.000000 -0.000000
4571
----------------------------------------
4573
Res 0.293583 sqrt( Sum f_i^2 / 3N )
4574
----------------------------------------
4575
Max 0.599298 constrained
5164
1 -0.000000 -0.000000 -0.000000
5165
2 -0.000000 -0.000000 0.000000
5166
3 -0.000001 0.000000 -0.000000
5167
4 0.000001 -0.000000 0.000000
5168
5 -0.267947 0.599297 0.000001
5169
6 -0.000001 -0.359548 0.549252
5170
7 0.267947 -0.239749 -0.549252
5171
8 0.267947 -0.599297 -0.000001
5172
9 0.000001 0.359548 -0.549252
5173
10 -0.267947 0.239749 0.549252
5174
----------------------------------------
5175
Tot 0.000000 0.000000 -0.000000
5176
----------------------------------------
5178
Res 0.293582 sqrt( Sum f_i^2 / 3N )
5179
----------------------------------------
5180
Max 0.599297 constrained
4577
5182
Stress-tensor-Voigt (kbar): 947.76 937.99 936.59 5.46 1.60 3.57
4578
5183
(Free)E + p*V (eV/cell) -2993.0400
4623
5228
Folding of H and S is implicitly performed
4624
5229
<dSpData1D:S at geom step 37
4625
5230
<sparsity:sparsity for geom step 37
4626
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
4627
<dData1D:(new from dSpData1D) n=306, refcount: 1>
5231
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
5232
<dData1D:(new from dSpData1D) n=204, refcount: 1>
4630
5235
Initializing Density Matrix...
4631
DM after filling with atomic data:
4632
<dSpData2D:(DM initialized from atoms)
5236
DM filled with atomic data:
5237
<dSpData2D:DM initialized from atoms
4633
5238
<sparsity:sparsity for geom step 37
4634
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
4635
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
5239
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
5240
<dData2D:DM n=204 m=1, refcount: 1>
4637
5243
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.84056646
5244
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4638
5245
New grid distribution: 1
5253
8 1: 12 7: 12 10: 12
4644
5255
InitMesh: MESH = 24 x 24 x 24 = 13824
4645
5256
InitMesh: (bp) = 12 x 12 x 12 = 1728
4646
5257
InitMesh: Mesh cutoff (required, used) = 50.000 64.614 Ry
4647
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5258
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
4648
5259
New grid distribution: 2
4653
5268
New grid distribution: 3
4658
5277
Setting up quadratic distribution...
4659
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4660
PhiOnMesh: Number of (b)points on node 0 = 432
4661
PhiOnMesh: nlist on node 0 = 10833
4663
siesta: Eharris(eV) = -2970.186321
4666
siesta: Eharris(eV) = -2970.186321
5278
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
5279
PhiOnMesh: Number of (b)points on node 0 = 216
5280
PhiOnMesh: nlist on node 0 = 5252
5282
siesta: Eharris(eV) = -2970.186321
5285
siesta: Eharris(eV) = -2970.186321
5287
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
5288
Geom step, scf iteration, dmax: 37 1 1.813110
4669
5290
Using DM_out to compute the final energy and forces
5291
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4671
5293
siesta: Atomic forces (eV/Ang):
4672
1 0.000000 -0.000000 0.000000
4673
2 -0.000000 0.000000 0.000000
4674
3 0.000000 -0.000000 -0.000000
4675
4 -0.000000 -0.000000 -0.000000
4676
5 -0.282595 0.632060 -0.000000
4677
6 -0.000000 -0.379205 0.579279
4678
7 0.282595 -0.252856 -0.579279
4679
8 0.282595 -0.632060 -0.000000
4680
9 -0.000000 0.379205 -0.579279
4681
10 -0.282595 0.252856 0.579279
4682
----------------------------------------
4683
Tot -0.000000 -0.000000 -0.000000
4684
----------------------------------------
5294
1 -0.000000 0.000000 0.000000
5295
2 -0.000000 -0.000000 -0.000000
5296
3 0.000001 -0.000000 -0.000000
5297
4 -0.000001 0.000000 0.000000
5298
5 -0.282596 0.632061 -0.000001
5299
6 0.000000 -0.379204 0.579279
5300
7 0.282595 -0.252856 -0.579280
5301
8 0.282596 -0.632061 0.000001
5302
9 -0.000000 0.379204 -0.579279
5303
10 -0.282595 0.252856 0.579280
5304
----------------------------------------
5305
Tot 0.000000 -0.000000 -0.000000
5306
----------------------------------------
4686
5308
Res 0.309632 sqrt( Sum f_i^2 / 3N )
4687
5309
----------------------------------------
4688
Max 0.632060 constrained
5310
Max 0.632061 constrained
4690
5312
Stress-tensor-Voigt (kbar): 959.27 955.29 954.72 2.22 0.65 1.45
4691
5313
(Free)E + p*V (eV/cell) -2993.2423
4736
5358
Folding of H and S is implicitly performed
4737
5359
<dSpData1D:S at geom step 38
4738
5360
<sparsity:sparsity for geom step 38
4739
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
4740
<dData1D:(new from dSpData1D) n=306, refcount: 1>
5361
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
5362
<dData1D:(new from dSpData1D) n=204, refcount: 1>
4743
5365
Initializing Density Matrix...
4744
DM after filling with atomic data:
4745
<dSpData2D:(DM initialized from atoms)
5366
DM filled with atomic data:
5367
<dSpData2D:DM initialized from atoms
4746
5368
<sparsity:sparsity for geom step 38
4747
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
4748
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
5369
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
5370
<dData2D:DM n=204 m=1, refcount: 1>
4750
5373
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.85903053
5374
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4751
5375
New grid distribution: 1
5383
8 1: 12 7: 12 10: 12
4757
5385
InitMesh: MESH = 24 x 24 x 24 = 13824
4758
5386
InitMesh: (bp) = 12 x 12 x 12 = 1728
4759
5387
InitMesh: Mesh cutoff (required, used) = 50.000 64.678 Ry
4760
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5388
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
4761
5389
New grid distribution: 2
4766
5398
New grid distribution: 3
4771
5407
Setting up quadratic distribution...
4772
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4773
PhiOnMesh: Number of (b)points on node 0 = 432
4774
PhiOnMesh: nlist on node 0 = 10909
4776
siesta: Eharris(eV) = -2970.272687
4779
siesta: Eharris(eV) = -2970.272687
5408
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
5409
PhiOnMesh: Number of (b)points on node 0 = 216
5410
PhiOnMesh: nlist on node 0 = 5301
5412
siesta: Eharris(eV) = -2970.272685
5415
siesta: Eharris(eV) = -2970.272685
5417
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
5418
Geom step, scf iteration, dmax: 38 1 1.814940
4782
5420
Using DM_out to compute the final energy and forces
5421
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4784
5423
siesta: Atomic forces (eV/Ang):
4785
1 0.000000 -0.000000 0.000000
4786
2 -0.000000 0.000000 0.000000
4787
3 0.000000 -0.000000 -0.000000
4788
4 0.000000 -0.000000 -0.000000
4789
5 -0.309008 0.691135 -0.000000
4790
6 -0.000000 -0.414647 0.633421
5424
1 -0.000000 -0.000000 -0.000000
5425
2 -0.000000 -0.000000 0.000000
5426
3 0.000001 -0.000001 0.000000
5427
4 -0.000001 0.000001 -0.000000
5428
5 -0.309009 0.691135 0.000000
5429
6 -0.000001 -0.414646 0.633421
4791
5430
7 0.309008 -0.276489 -0.633421
4792
8 0.309008 -0.691135 -0.000000
4793
9 -0.000000 0.414647 -0.633421
5431
8 0.309009 -0.691135 -0.000000
5432
9 0.000001 0.414646 -0.633421
4794
5433
10 -0.309008 0.276489 0.633421
4795
5434
----------------------------------------
4796
Tot -0.000000 -0.000000 -0.000000
5435
Tot 0.000000 0.000000 0.000000
4797
5436
----------------------------------------
4799
5438
Res 0.338572 sqrt( Sum f_i^2 / 3N )
4849
5488
Folding of H and S is implicitly performed
4850
5489
<dSpData1D:S at geom step 39
4851
5490
<sparsity:sparsity for geom step 39
4852
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
4853
<dData1D:(new from dSpData1D) n=306, refcount: 1>
5491
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
5492
<dData1D:(new from dSpData1D) n=204, refcount: 1>
4856
5495
Initializing Density Matrix...
4857
DM after filling with atomic data:
4858
<dSpData2D:(DM initialized from atoms)
5496
DM filled with atomic data:
5497
<dSpData2D:DM initialized from atoms
4859
5498
<sparsity:sparsity for geom step 39
4860
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
4861
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
5499
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
5500
<dData2D:DM n=204 m=1, refcount: 1>
4863
5503
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.87851805
5504
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4864
5505
New grid distribution: 1
5513
8 1: 12 7: 12 10: 12
4870
5515
InitMesh: MESH = 24 x 24 x 24 = 13824
4871
5516
InitMesh: (bp) = 12 x 12 x 12 = 1728
4872
5517
InitMesh: Mesh cutoff (required, used) = 50.000 64.745 Ry
4873
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5518
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
4874
5519
New grid distribution: 2
4879
5528
New grid distribution: 3
4884
5537
Setting up quadratic distribution...
4885
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4886
PhiOnMesh: Number of (b)points on node 0 = 432
4887
PhiOnMesh: nlist on node 0 = 10951
4889
siesta: Eharris(eV) = -2970.359177
4892
siesta: Eharris(eV) = -2970.359177
5538
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
5539
PhiOnMesh: Number of (b)points on node 0 = 216
5540
PhiOnMesh: nlist on node 0 = 5331
5542
siesta: Eharris(eV) = -2970.359177
5545
siesta: Eharris(eV) = -2970.359177
5547
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
5548
Geom step, scf iteration, dmax: 39 1 1.816910
4895
5550
Using DM_out to compute the final energy and forces
5551
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4897
5553
siesta: Atomic forces (eV/Ang):
4898
1 0.000000 -0.000000 0.000000
4899
2 -0.000000 0.000000 0.000000
4900
3 0.000000 -0.000000 -0.000000
4901
4 0.000000 -0.000000 -0.000000
4902
5 -0.344531 0.770587 -0.000000
4903
6 -0.000000 -0.462314 0.706238
4904
7 0.344531 -0.308273 -0.706238
5554
1 -0.000000 -0.000000 0.000000
5555
2 -0.000000 0.000000 -0.000000
5556
3 0.000000 0.000000 0.000000
5557
4 -0.000000 -0.000000 -0.000000
5558
5 -0.344531 0.770587 0.000000
5559
6 0.000001 -0.462314 0.706237
5560
7 0.344531 -0.308274 -0.706238
4905
5561
8 0.344531 -0.770587 -0.000000
4906
9 -0.000000 0.462314 -0.706238
4907
10 -0.344531 0.308273 0.706238
5562
9 -0.000001 0.462314 -0.706237
5563
10 -0.344531 0.308274 0.706238
4908
5564
----------------------------------------
4909
Tot -0.000000 -0.000000 -0.000000
5565
Tot 0.000000 0.000000 0.000000
4910
5566
----------------------------------------
4912
5568
Res 0.377493 sqrt( Sum f_i^2 / 3N )
4962
5618
Folding of H and S is implicitly performed
4963
5619
<dSpData1D:S at geom step 40
4964
5620
<sparsity:sparsity for geom step 40
4965
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
4966
<dData1D:(new from dSpData1D) n=306, refcount: 1>
5621
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
5622
<dData1D:(new from dSpData1D) n=204, refcount: 1>
4969
5625
Initializing Density Matrix...
4970
DM after filling with atomic data:
4971
<dSpData2D:(DM initialized from atoms)
5626
DM filled with atomic data:
5627
<dSpData2D:DM initialized from atoms
4972
5628
<sparsity:sparsity for geom step 40
4973
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
4974
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
5629
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
5630
<dData2D:DM n=204 m=1, refcount: 1>
4976
5633
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.89907685
5634
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
4977
5635
New grid distribution: 1
5643
8 1: 12 7: 12 10: 12
4983
5645
InitMesh: MESH = 24 x 24 x 24 = 13824
4984
5646
InitMesh: (bp) = 12 x 12 x 12 = 1728
4985
5647
InitMesh: Mesh cutoff (required, used) = 50.000 64.813 Ry
4986
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5648
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
4987
5649
New grid distribution: 2
4992
5658
New grid distribution: 3
4997
5667
Setting up quadratic distribution...
4998
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
4999
PhiOnMesh: Number of (b)points on node 0 = 432
5000
PhiOnMesh: nlist on node 0 = 11014
5002
siesta: Eharris(eV) = -2970.445862
5005
siesta: Eharris(eV) = -2970.445862
5668
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
5669
PhiOnMesh: Number of (b)points on node 0 = 216
5670
PhiOnMesh: nlist on node 0 = 5358
5672
siesta: Eharris(eV) = -2970.445863
5675
siesta: Eharris(eV) = -2970.445863
5677
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
5678
Geom step, scf iteration, dmax: 40 1 1.818983
5008
5680
Using DM_out to compute the final energy and forces
5681
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
5010
5683
siesta: Atomic forces (eV/Ang):
5011
1 0.000000 -0.000000 0.000000
5684
1 -0.000000 0.000000 -0.000000
5012
5685
2 -0.000000 0.000000 0.000000
5013
3 0.000000 -0.000000 -0.000000
5014
4 0.000000 -0.000000 -0.000000
5015
5 -0.387952 0.867705 -0.000000
5016
6 -0.000000 -0.520580 0.795246
5686
3 -0.000000 0.000000 -0.000000
5687
4 0.000000 -0.000000 0.000000
5688
5 -0.387952 0.867705 -0.000001
5689
6 -0.000001 -0.520579 0.795246
5017
5690
7 0.387952 -0.347125 -0.795246
5018
8 0.387952 -0.867705 -0.000000
5019
9 -0.000000 0.520580 -0.795246
5691
8 0.387952 -0.867705 0.000001
5692
9 0.000001 0.520579 -0.795246
5020
5693
10 -0.387952 0.347125 0.795246
5021
5694
----------------------------------------
5022
Tot -0.000000 -0.000000 -0.000000
5695
Tot 0.000000 0.000000 0.000000
5023
5696
----------------------------------------
5025
5698
Res 0.425069 sqrt( Sum f_i^2 / 3N )
5035
5708
-0.010056 0.638196 -0.002943
5036
5709
-0.006583 -0.002943 0.640766
5038
siesta: Pressure (static): -1014.27595912 kBar
5711
siesta: Pressure (static): -1014.27596406 kBar
5040
5713
siesta: Stress tensor (total) (eV/Ang**3):
5041
5714
0.613744 0.000054 0.000036
5042
5715
0.000054 0.613647 0.000016
5043
5716
0.000036 0.000016 0.613633
5045
siesta: Pressure (total): -983.22542515 kBar
5047
siesta: E_kin PR = 63450.672 K
5049
siesta: Temp_ion = 2090.748 K
5718
siesta: Pressure (total): -983.22543168 kBar
5720
siesta: E_kin PR = 63450.671 K
5722
siesta: Temp_ion = 2090.747 K
5051
5724
====================================
5052
5725
Begin MD step = 41
5075
5748
Folding of H and S is implicitly performed
5076
5749
<dSpData1D:S at geom step 41
5077
5750
<sparsity:sparsity for geom step 41
5078
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
5079
<dData1D:(new from dSpData1D) n=306, refcount: 1>
5751
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
5752
<dData1D:(new from dSpData1D) n=204, refcount: 1>
5082
5755
Initializing Density Matrix...
5083
DM after filling with atomic data:
5084
<dSpData2D:(DM initialized from atoms)
5756
DM filled with atomic data:
5757
<dSpData2D:DM initialized from atoms
5085
5758
<sparsity:sparsity for geom step 41
5086
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
5087
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
5759
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
5760
<dData2D:DM n=204 m=1, refcount: 1>
5089
5763
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.92075217
5764
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
5090
5765
New grid distribution: 1
5773
8 1: 12 7: 12 10: 12
5096
5775
InitMesh: MESH = 24 x 24 x 24 = 13824
5097
5776
InitMesh: (bp) = 12 x 12 x 12 = 1728
5098
5777
InitMesh: Mesh cutoff (required, used) = 50.000 64.883 Ry
5099
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5778
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
5100
5779
New grid distribution: 2
5105
5788
New grid distribution: 3
5110
5797
Setting up quadratic distribution...
5111
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5112
PhiOnMesh: Number of (b)points on node 0 = 432
5113
PhiOnMesh: nlist on node 0 = 11066
5115
siesta: Eharris(eV) = -2970.533301
5118
siesta: Eharris(eV) = -2970.533301
5798
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
5799
PhiOnMesh: Number of (b)points on node 0 = 216
5800
PhiOnMesh: nlist on node 0 = 5386
5802
siesta: Eharris(eV) = -2970.533301
5805
siesta: Eharris(eV) = -2970.533301
5807
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
5808
Geom step, scf iteration, dmax: 41 1 1.821157
5121
5810
Using DM_out to compute the final energy and forces
5811
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 85
5123
5813
siesta: Atomic forces (eV/Ang):
5124
1 0.000000 -0.000000 0.000000
5125
2 -0.000000 0.000000 0.000000
5126
3 0.000000 -0.000000 -0.000000
5127
4 0.000000 0.000000 -0.000000
5128
5 -0.420637 0.940809 -0.000000
5129
6 -0.000000 -0.564439 0.862245
5130
7 0.420637 -0.376371 -0.862245
5131
8 0.420637 -0.940809 -0.000000
5132
9 -0.000000 0.564439 -0.862245
5133
10 -0.420637 0.376371 0.862245
5814
1 0.000000 0.000000 0.000000
5815
2 0.000000 -0.000000 0.000000
5816
3 0.000001 -0.000000 0.000000
5817
4 -0.000001 0.000000 -0.000000
5818
5 -0.420638 0.940809 -0.000000
5819
6 -0.000000 -0.564438 0.862246
5820
7 0.420637 -0.376370 -0.862246
5821
8 0.420638 -0.940809 0.000000
5822
9 0.000000 0.564438 -0.862246
5823
10 -0.420637 0.376370 0.862246
5134
5824
----------------------------------------
5135
Tot -0.000000 -0.000000 -0.000000
5825
Tot 0.000000 0.000000 0.000000
5136
5826
----------------------------------------
5138
5828
Res 0.460881 sqrt( Sum f_i^2 / 3N )
5188
5878
Folding of H and S is implicitly performed
5189
5879
<dSpData1D:S at geom step 42
5190
5880
<sparsity:sparsity for geom step 42
5191
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
5192
<dData1D:(new from dSpData1D) n=306, refcount: 1>
5881
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
5882
<dData1D:(new from dSpData1D) n=204, refcount: 1>
5195
5885
Initializing Density Matrix...
5196
DM after filling with atomic data:
5197
<dSpData2D:(DM initialized from atoms)
5886
DM filled with atomic data:
5887
<dSpData2D:DM initialized from atoms
5198
5888
<sparsity:sparsity for geom step 42
5199
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
5200
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
5889
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
5890
<dData2D:DM n=204 m=1, refcount: 1>
5202
5893
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.94360035
5894
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5203
5895
New grid distribution: 1
5903
8 1: 12 7: 12 10: 12
5209
5905
InitMesh: MESH = 24 x 24 x 24 = 13824
5210
5906
InitMesh: (bp) = 12 x 12 x 12 = 1728
5211
5907
InitMesh: Mesh cutoff (required, used) = 50.000 64.956 Ry
5212
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5908
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
5213
5909
New grid distribution: 2
5218
5918
New grid distribution: 3
5223
5927
Setting up quadratic distribution...
5224
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5225
PhiOnMesh: Number of (b)points on node 0 = 432
5226
PhiOnMesh: nlist on node 0 = 11119
5228
siesta: Eharris(eV) = -2970.622780
5231
siesta: Eharris(eV) = -2970.622780
5928
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
5929
PhiOnMesh: Number of (b)points on node 0 = 216
5930
PhiOnMesh: nlist on node 0 = 5436
5932
siesta: Eharris(eV) = -2970.622779
5935
siesta: Eharris(eV) = -2970.622779
5937
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
5938
Geom step, scf iteration, dmax: 42 1 1.823434
5234
5940
Using DM_out to compute the final energy and forces
5941
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5236
5943
siesta: Atomic forces (eV/Ang):
5237
1 0.000000 -0.000000 0.000000
5944
1 -0.000000 0.000000 0.000000
5238
5945
2 -0.000000 0.000000 0.000000
5239
3 0.000000 -0.000000 -0.000000
5240
4 0.000000 0.000000 -0.000000
5241
5 -0.473087 1.058119 -0.000000
5242
6 -0.000000 -0.634818 0.969759
5243
7 0.473087 -0.423300 -0.969759
5244
8 0.473087 -1.058119 -0.000000
5245
9 -0.000000 0.634818 -0.969759
5246
10 -0.473087 0.423300 0.969759
5247
----------------------------------------
5248
Tot -0.000000 -0.000000 -0.000000
5249
----------------------------------------
5251
Res 0.518349 sqrt( Sum f_i^2 / 3N )
5252
----------------------------------------
5253
Max 1.058119 constrained
5946
3 0.000000 0.000001 -0.000000
5947
4 -0.000000 -0.000001 0.000000
5948
5 -0.473087 1.058117 -0.000000
5949
6 -0.000000 -0.634818 0.969757
5950
7 0.473087 -0.423301 -0.969758
5951
8 0.473087 -1.058117 0.000000
5952
9 0.000000 0.634818 -0.969757
5953
10 -0.473087 0.423301 0.969758
5954
----------------------------------------
5955
Tot 0.000000 0.000000 0.000000
5956
----------------------------------------
5958
Res 0.518348 sqrt( Sum f_i^2 / 3N )
5959
----------------------------------------
5960
Max 1.058117 constrained
5255
5962
Stress-tensor-Voigt (kbar): 1004.78 1017.35 1019.15 -7.02 -2.06 -4.60
5256
5963
(Free)E + p*V (eV/cell) -2992.6085
5301
6008
Folding of H and S is implicitly performed
5302
6009
<dSpData1D:S at geom step 43
5303
6010
<sparsity:sparsity for geom step 43
5304
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
5305
<dData1D:(new from dSpData1D) n=306, refcount: 1>
6011
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
6012
<dData1D:(new from dSpData1D) n=204, refcount: 1>
5308
6015
Initializing Density Matrix...
5309
DM after filling with atomic data:
5310
<dSpData2D:(DM initialized from atoms)
6016
DM filled with atomic data:
6017
<dSpData2D:DM initialized from atoms
5311
6018
<sparsity:sparsity for geom step 43
5312
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
5313
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
6019
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
6020
<dData2D:DM n=204 m=1, refcount: 1>
5315
6023
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.96767783
6024
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5316
6025
New grid distribution: 1
6033
8 1: 12 7: 12 10: 12
5322
6035
InitMesh: MESH = 24 x 24 x 24 = 13824
5323
6036
InitMesh: (bp) = 12 x 12 x 12 = 1728
5324
6037
InitMesh: Mesh cutoff (required, used) = 50.000 65.030 Ry
5325
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
6038
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
5326
6039
New grid distribution: 2
5331
6048
New grid distribution: 3
5336
6057
Setting up quadratic distribution...
5337
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5338
PhiOnMesh: Number of (b)points on node 0 = 432
5339
PhiOnMesh: nlist on node 0 = 11205
5341
siesta: Eharris(eV) = -2970.708152
5344
siesta: Eharris(eV) = -2970.708152
6058
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
6059
PhiOnMesh: Number of (b)points on node 0 = 216
6060
PhiOnMesh: nlist on node 0 = 5468
6062
siesta: Eharris(eV) = -2970.708148
6065
siesta: Eharris(eV) = -2970.708148
6067
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
6068
Geom step, scf iteration, dmax: 43 1 1.825817
5347
6070
Using DM_out to compute the final energy and forces
6071
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5349
6073
siesta: Atomic forces (eV/Ang):
5350
1 0.000000 -0.000000 0.000000
5351
2 -0.000000 0.000000 0.000000
5352
3 0.000000 -0.000000 -0.000000
5353
4 0.000000 0.000000 -0.000000
5354
5 -0.556469 1.244613 -0.000000
5355
6 -0.000000 -0.746705 1.140679
5356
7 0.556469 -0.497907 -1.140679
5357
8 0.556469 -1.244613 -0.000000
5358
9 -0.000000 0.746705 -1.140679
5359
10 -0.556469 0.497907 1.140679
5360
----------------------------------------
5361
Tot -0.000000 -0.000000 -0.000000
5362
----------------------------------------
6074
1 -0.000000 0.000000 0.000000
6075
2 -0.000000 0.000000 -0.000000
6076
3 -0.000001 -0.000000 0.000000
6077
4 0.000001 0.000000 -0.000000
6078
5 -0.556469 1.244614 -0.000000
6079
6 -0.000001 -0.746706 1.140679
6080
7 0.556469 -0.497908 -1.140680
6081
8 0.556469 -1.244614 0.000000
6082
9 0.000001 0.746706 -1.140679
6083
10 -0.556469 0.497908 1.140680
6084
----------------------------------------
6085
Tot 0.000000 0.000000 0.000000
6086
----------------------------------------
5364
6088
Res 0.609708 sqrt( Sum f_i^2 / 3N )
5365
6089
----------------------------------------
5366
Max 1.244613 constrained
6090
Max 1.244614 constrained
5368
6092
Stress-tensor-Voigt (kbar): 1011.94 1023.62 1025.28 -6.52 -1.91 -4.27
5369
6093
(Free)E + p*V (eV/cell) -2992.1260
5414
6138
Folding of H and S is implicitly performed
5415
6139
<dSpData1D:S at geom step 44
5416
6140
<sparsity:sparsity for geom step 44
5417
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
5418
<dData1D:(new from dSpData1D) n=306, refcount: 1>
6141
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
6142
<dData1D:(new from dSpData1D) n=204, refcount: 1>
5421
6145
Initializing Density Matrix...
5422
DM after filling with atomic data:
5423
<dSpData2D:(DM initialized from atoms)
6146
DM filled with atomic data:
6147
<dSpData2D:DM initialized from atoms
5424
6148
<sparsity:sparsity for geom step 44
5425
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
5426
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
6149
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
6150
<dData2D:DM n=204 m=1, refcount: 1>
5428
6153
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 48.99303621
6154
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5429
6155
New grid distribution: 1
6163
8 1: 12 7: 12 10: 12
5435
6165
InitMesh: MESH = 24 x 24 x 24 = 13824
5436
6166
InitMesh: (bp) = 12 x 12 x 12 = 1728
5437
6167
InitMesh: Mesh cutoff (required, used) = 50.000 65.106 Ry
5438
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
6168
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
5439
6169
New grid distribution: 2
5444
6178
New grid distribution: 3
5449
6187
Setting up quadratic distribution...
5450
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5451
PhiOnMesh: Number of (b)points on node 0 = 432
5452
PhiOnMesh: nlist on node 0 = 11300
5454
siesta: Eharris(eV) = -2970.789167
5457
siesta: Eharris(eV) = -2970.789167
6188
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
6189
PhiOnMesh: Number of (b)points on node 0 = 216
6190
PhiOnMesh: nlist on node 0 = 5518
6192
siesta: Eharris(eV) = -2970.789168
6195
siesta: Eharris(eV) = -2970.789168
6197
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
6198
Geom step, scf iteration, dmax: 44 1 1.828259
5460
6200
Using DM_out to compute the final energy and forces
6201
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5462
6203
siesta: Atomic forces (eV/Ang):
5463
1 0.000000 -0.000000 0.000000
6204
1 -0.000000 0.000000 -0.000000
5464
6205
2 -0.000000 0.000000 0.000000
5465
3 0.000000 -0.000000 -0.000000
5466
4 0.000000 0.000000 0.000000
5467
5 -0.619094 1.384682 -0.000000
5468
6 -0.000000 -0.830740 1.269052
5469
7 0.619094 -0.553942 -1.269052
5470
8 0.619094 -1.384682 -0.000000
5471
9 -0.000000 0.830740 -1.269052
5472
10 -0.619094 0.553942 1.269052
6206
3 -0.000001 -0.000000 0.000000
6207
4 0.000001 0.000000 -0.000000
6208
5 -0.619095 1.384682 -0.000000
6209
6 0.000001 -0.830740 1.269051
6210
7 0.619094 -0.553941 -1.269053
6211
8 0.619095 -1.384682 0.000000
6212
9 -0.000001 0.830740 -1.269051
6213
10 -0.619094 0.553941 1.269053
5473
6214
----------------------------------------
5474
Tot -0.000000 -0.000000 -0.000000
6215
Tot 0.000000 0.000000 0.000000
5475
6216
----------------------------------------
5477
6218
Res 0.678325 sqrt( Sum f_i^2 / 3N )
5483
6224
Target enthalpy (eV/cell) -2934.6946
5485
6226
siesta: Stress tensor (static) (eV/Ang**3):
5486
0.644757 -0.015413 -0.010089
6227
0.644756 -0.015413 -0.010089
5487
6228
-0.015413 0.672338 -0.004511
5488
6229
-0.010089 -0.004511 0.676276
5490
siesta: Pressure (static): -1064.58877928 kBar
6231
siesta: Pressure (static): -1064.58873805 kBar
5492
6233
siesta: Stress tensor (total) (eV/Ang**3):
5493
6234
0.638405 -0.005468 -0.003580
5494
6235
-0.005468 0.648191 -0.001600
5495
6236
-0.003580 -0.001600 0.649588
5497
siesta: Pressure (total): -1034.04761951 kBar
6238
siesta: Pressure (total): -1034.04757936 kBar
5499
siesta: E_kin PR = 82442.569 K
6240
siesta: E_kin PR = 82442.568 K
5501
6242
siesta: Temp_ion = 1983.493 K
5527
6268
Folding of H and S is implicitly performed
5528
6269
<dSpData1D:S at geom step 45
5529
6270
<sparsity:sparsity for geom step 45
5530
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
5531
<dData1D:(new from dSpData1D) n=306, refcount: 1>
6271
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
6272
<dData1D:(new from dSpData1D) n=204, refcount: 1>
5534
6275
Initializing Density Matrix...
5535
DM after filling with atomic data:
5536
<dSpData2D:(DM initialized from atoms)
6276
DM filled with atomic data:
6277
<dSpData2D:DM initialized from atoms
5537
6278
<sparsity:sparsity for geom step 45
5538
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
5539
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
6279
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
6280
<dData2D:DM n=204 m=1, refcount: 1>
5541
6283
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 49.01973671
6284
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5542
6285
New grid distribution: 1
6293
8 1: 12 7: 12 10: 12
5548
6295
InitMesh: MESH = 24 x 24 x 24 = 13824
5549
6296
InitMesh: (bp) = 12 x 12 x 12 = 1728
5550
6297
InitMesh: Mesh cutoff (required, used) = 50.000 65.184 Ry
5551
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
6298
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
5552
6299
New grid distribution: 2
5557
6308
New grid distribution: 3
5562
6317
Setting up quadratic distribution...
5563
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5564
PhiOnMesh: Number of (b)points on node 0 = 432
5565
PhiOnMesh: nlist on node 0 = 11420
5567
siesta: Eharris(eV) = -2970.867119
5570
siesta: Eharris(eV) = -2970.867119
6318
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
6319
PhiOnMesh: Number of (b)points on node 0 = 216
6320
PhiOnMesh: nlist on node 0 = 5572
6322
siesta: Eharris(eV) = -2970.867119
6325
siesta: Eharris(eV) = -2970.867119
6327
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
6328
Geom step, scf iteration, dmax: 45 1 1.830712
5573
6330
Using DM_out to compute the final energy and forces
6331
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5575
6333
siesta: Atomic forces (eV/Ang):
5576
1 0.000000 0.000000 0.000000
5577
2 -0.000000 0.000000 0.000000
5578
3 0.000000 -0.000000 -0.000000
6334
1 0.000000 -0.000000 -0.000000
6335
2 0.000000 -0.000000 0.000000
6336
3 -0.000000 -0.000000 0.000000
5579
6337
4 0.000000 0.000000 -0.000000
5580
5 -0.671184 1.501189 -0.000000
5581
6 -0.000000 -0.900638 1.375829
5582
7 0.671184 -0.600551 -1.375829
5583
8 0.671184 -1.501189 -0.000000
5584
9 -0.000000 0.900638 -1.375829
5585
10 -0.671184 0.600551 1.375829
5586
----------------------------------------
5587
Tot -0.000000 -0.000000 -0.000000
5588
----------------------------------------
5590
Res 0.735399 sqrt( Sum f_i^2 / 3N )
5591
----------------------------------------
5592
Max 1.501189 constrained
6338
5 -0.671184 1.501188 -0.000001
6339
6 0.000000 -0.900639 1.375829
6340
7 0.671184 -0.600550 -1.375829
6341
8 0.671184 -1.501188 0.000001
6342
9 -0.000000 0.900639 -1.375829
6343
10 -0.671184 0.600550 1.375829
6344
----------------------------------------
6345
Tot 0.000000 0.000000 -0.000000
6346
----------------------------------------
6348
Res 0.735398 sqrt( Sum f_i^2 / 3N )
6349
----------------------------------------
6350
Max 1.501188 constrained
5594
6352
Stress-tensor-Voigt (kbar): 1031.18 1045.45 1047.49 -7.98 -2.33 -5.22
5595
6353
(Free)E + p*V (eV/cell) -2991.4655
5622
6380
0.50000000 0.50000000 0.50000000 1 2 Mg
5623
6381
0.25000000 0.25000000 0.25000000 2 3 C
5624
6382
-0.25000000 -0.25000000 -0.25000000 2 4 C
5625
0.57974821 -0.07974821 0.25000000 3 5 O
6383
0.57974820 -0.07974820 0.25000000 3 5 O
5626
6384
0.25000000 0.57974821 -0.07974821 3 6 O
5627
-0.07974821 0.25000000 0.57974821 3 7 O
5628
-0.57974821 0.07974821 -0.25000000 3 8 O
6385
-0.07974820 0.25000000 0.57974820 3 7 O
6386
-0.57974820 0.07974820 -0.25000000 3 8 O
5629
6387
-0.25000000 -0.57974821 0.07974821 3 9 O
5630
0.07974821 -0.25000000 -0.57974821 3 10 O
6388
0.07974820 -0.25000000 -0.57974820 3 10 O
5632
6390
outcell: Unit cell vectors (Ang):
5633
6391
5.530139 0.109015 0.071363
5640
6398
Folding of H and S is implicitly performed
5641
6399
<dSpData1D:S at geom step 46
5642
6400
<sparsity:sparsity for geom step 46
5643
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
5644
<dData1D:(new from dSpData1D) n=306, refcount: 1>
6401
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
6402
<dData1D:(new from dSpData1D) n=204, refcount: 1>
5647
6405
Initializing Density Matrix...
5648
DM after filling with atomic data:
5649
<dSpData2D:(DM initialized from atoms)
6406
DM filled with atomic data:
6407
<dSpData2D:DM initialized from atoms
5650
6408
<sparsity:sparsity for geom step 46
5651
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
5652
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
6409
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
6410
<dData2D:DM n=204 m=1, refcount: 1>
5654
6413
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 49.04783669
6414
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5655
6415
New grid distribution: 1
6423
8 1: 12 7: 12 10: 12
5661
6425
InitMesh: MESH = 24 x 24 x 24 = 13824
5662
6426
InitMesh: (bp) = 12 x 12 x 12 = 1728
5663
6427
InitMesh: Mesh cutoff (required, used) = 50.000 65.264 Ry
5664
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
6428
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
5665
6429
New grid distribution: 2
5670
6438
New grid distribution: 3
5675
6447
Setting up quadratic distribution...
5676
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5677
PhiOnMesh: Number of (b)points on node 0 = 432
5678
PhiOnMesh: nlist on node 0 = 11536
5680
siesta: Eharris(eV) = -2970.941360
5683
siesta: Eharris(eV) = -2970.941360
6448
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
6449
PhiOnMesh: Number of (b)points on node 0 = 216
6450
PhiOnMesh: nlist on node 0 = 5610
6452
siesta: Eharris(eV) = -2970.941359
6455
siesta: Eharris(eV) = -2970.941359
6457
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
6458
Geom step, scf iteration, dmax: 46 1 1.833231
5686
6460
Using DM_out to compute the final energy and forces
6461
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5688
6463
siesta: Atomic forces (eV/Ang):
5689
1 0.000000 0.000000 0.000000
6464
1 -0.000000 -0.000000 0.000000
5690
6465
2 -0.000000 0.000000 0.000000
5691
3 0.000000 -0.000000 -0.000000
5692
4 0.000000 0.000000 -0.000000
5693
5 -0.722108 1.615086 -0.000000
5694
6 -0.000000 -0.968971 1.480215
5695
7 0.722108 -0.646115 -1.480215
6466
3 0.000000 0.000000 -0.000000
6467
4 -0.000000 -0.000000 0.000000
6468
5 -0.722108 1.615086 0.000000
6469
6 0.000000 -0.968971 1.480216
6470
7 0.722107 -0.646114 -1.480215
5696
6471
8 0.722108 -1.615086 -0.000000
5697
9 -0.000000 0.968971 -1.480215
5698
10 -0.722108 0.646115 1.480215
6472
9 -0.000000 0.968971 -1.480216
6473
10 -0.722107 0.646114 1.480215
5699
6474
----------------------------------------
5700
Tot -0.000000 -0.000000 -0.000000
6475
Tot 0.000000 0.000000 0.000000
5701
6476
----------------------------------------
5703
6478
Res 0.791194 sqrt( Sum f_i^2 / 3N )
5753
6528
Folding of H and S is implicitly performed
5754
6529
<dSpData1D:S at geom step 47
5755
6530
<sparsity:sparsity for geom step 47
5756
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
5757
<dData1D:(new from dSpData1D) n=306, refcount: 1>
6531
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
6532
<dData1D:(new from dSpData1D) n=204, refcount: 1>
5760
6535
Initializing Density Matrix...
5761
DM after filling with atomic data:
5762
<dSpData2D:(DM initialized from atoms)
6536
DM filled with atomic data:
6537
<dSpData2D:DM initialized from atoms
5763
6538
<sparsity:sparsity for geom step 47
5764
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
5765
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
6539
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
6540
<dData2D:DM n=204 m=1, refcount: 1>
5767
6543
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 49.07739911
6544
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5768
6545
New grid distribution: 1
6553
8 1: 12 7: 12 10: 12
5774
6555
InitMesh: MESH = 24 x 24 x 24 = 13824
5775
6556
InitMesh: (bp) = 12 x 12 x 12 = 1728
5776
6557
InitMesh: Mesh cutoff (required, used) = 50.000 65.347 Ry
5777
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
6558
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
5778
6559
New grid distribution: 2
5783
6568
New grid distribution: 3
5788
6577
Setting up quadratic distribution...
5789
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5790
PhiOnMesh: Number of (b)points on node 0 = 432
5791
PhiOnMesh: nlist on node 0 = 11662
5793
siesta: Eharris(eV) = -2971.011535
5796
siesta: Eharris(eV) = -2971.011535
6578
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
6579
PhiOnMesh: Number of (b)points on node 0 = 216
6580
PhiOnMesh: nlist on node 0 = 5661
6582
siesta: Eharris(eV) = -2971.011535
6585
siesta: Eharris(eV) = -2971.011535
6587
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
6588
Geom step, scf iteration, dmax: 47 1 1.835886
5799
6590
Using DM_out to compute the final energy and forces
6591
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5801
6593
siesta: Atomic forces (eV/Ang):
5802
1 0.000000 0.000000 -0.000000
6594
1 -0.000000 0.000000 -0.000000
5803
6595
2 -0.000000 0.000000 0.000000
5804
3 0.000000 -0.000000 -0.000000
5805
4 0.000000 0.000000 -0.000000
5806
5 -0.788851 1.764366 -0.000000
5807
6 -0.000000 -1.058532 1.617030
5808
7 0.788851 -0.705835 -1.617030
5809
8 0.788851 -1.764366 -0.000000
5810
9 -0.000000 1.058532 -1.617030
5811
10 -0.788851 0.705835 1.617030
5812
----------------------------------------
5813
Tot -0.000000 -0.000000 -0.000000
5814
----------------------------------------
6596
3 0.000000 0.000000 0.000001
6597
4 -0.000000 -0.000000 -0.000001
6598
5 -0.788852 1.764367 -0.000001
6599
6 0.000000 -1.058532 1.617029
6600
7 0.788851 -0.705834 -1.617030
6601
8 0.788852 -1.764367 0.000001
6602
9 -0.000000 1.058532 -1.617029
6603
10 -0.788851 0.705834 1.617030
6604
----------------------------------------
6605
Tot 0.000000 -0.000000 0.000000
6606
----------------------------------------
5816
6608
Res 0.864323 sqrt( Sum f_i^2 / 3N )
5817
6609
----------------------------------------
5818
Max 1.764366 constrained
6610
Max 1.764367 constrained
5820
6612
Stress-tensor-Voigt (kbar): 1047.15 1065.58 1068.21 -10.29 -3.01 -6.74
5821
6613
(Free)E + p*V (eV/cell) -2990.5586
5826
6618
-0.016033 0.688403 -0.004693
5827
6619
-0.010495 -0.004693 0.692500
5829
siesta: Pressure (static): -1089.82005652 kBar
6621
siesta: Pressure (static): -1089.82001093 kBar
5831
6623
siesta: Stress tensor (total) (eV/Ang**3):
5832
6624
0.653575 -0.006425 -0.004206
5833
6625
-0.006425 0.665074 -0.001881
5834
6626
-0.004206 -0.001881 0.666716
5836
siesta: Pressure (total): -1060.31347765 kBar
5838
siesta: E_kin PR = 98646.328 K
5840
siesta: Temp_ion = 1858.661 K
6628
siesta: Pressure (total): -1060.31343317 kBar
6630
siesta: E_kin PR = 98646.327 K
6632
siesta: Temp_ion = 1858.660 K
5842
6634
====================================
5843
6635
Begin MD step = 48
5866
6658
Folding of H and S is implicitly performed
5867
6659
<dSpData1D:S at geom step 48
5868
6660
<sparsity:sparsity for geom step 48
5869
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
5870
<dData1D:(new from dSpData1D) n=306, refcount: 1>
6661
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
6662
<dData1D:(new from dSpData1D) n=204, refcount: 1>
5873
6665
Initializing Density Matrix...
5874
DM after filling with atomic data:
5875
<dSpData2D:(DM initialized from atoms)
6666
DM filled with atomic data:
6667
<dSpData2D:DM initialized from atoms
5876
6668
<sparsity:sparsity for geom step 48
5877
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
5878
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
6669
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
6670
<dData2D:DM n=204 m=1, refcount: 1>
5880
6673
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 49.10849019
6674
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5881
6675
New grid distribution: 1
6683
8 1: 12 7: 12 10: 12
5887
6685
InitMesh: MESH = 24 x 24 x 24 = 13824
5888
6686
InitMesh: (bp) = 12 x 12 x 12 = 1728
5889
6687
InitMesh: Mesh cutoff (required, used) = 50.000 65.431 Ry
5890
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
6688
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
5891
6689
New grid distribution: 2
5896
6698
New grid distribution: 3
5901
6707
Setting up quadratic distribution...
5902
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
5903
PhiOnMesh: Number of (b)points on node 0 = 432
5904
PhiOnMesh: nlist on node 0 = 11736
5906
siesta: Eharris(eV) = -2971.072742
5909
siesta: Eharris(eV) = -2971.072742
6708
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
6709
PhiOnMesh: Number of (b)points on node 0 = 216
6710
PhiOnMesh: nlist on node 0 = 5708
6712
siesta: Eharris(eV) = -2971.072741
6715
siesta: Eharris(eV) = -2971.072741
6717
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
6718
Geom step, scf iteration, dmax: 48 1 1.838629
5912
6720
Using DM_out to compute the final energy and forces
6721
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5914
6723
siesta: Atomic forces (eV/Ang):
5915
1 0.000000 0.000000 -0.000000
5916
2 -0.000000 0.000000 -0.000000
5917
3 0.000000 0.000000 -0.000000
5918
4 0.000000 0.000000 -0.000000
5919
5 -0.892698 1.996632 -0.000000
5920
6 -0.000000 -1.197879 1.829900
5921
7 0.892698 -0.798753 -1.829900
5922
8 0.892698 -1.996632 -0.000000
5923
9 -0.000000 1.197879 -1.829900
5924
10 -0.892698 0.798753 1.829900
5925
----------------------------------------
5926
Tot -0.000000 -0.000000 -0.000000
5927
----------------------------------------
6724
1 -0.000000 0.000000 -0.000000
6725
2 -0.000000 0.000000 0.000000
6726
3 -0.000000 0.000000 -0.000000
6727
4 0.000000 -0.000000 0.000000
6728
5 -0.892697 1.996631 -0.000000
6729
6 -0.000000 -1.197880 1.829900
6730
7 0.892698 -0.798752 -1.829900
6731
8 0.892697 -1.996631 0.000000
6732
9 0.000000 1.197880 -1.829900
6733
10 -0.892698 0.798752 1.829900
6734
----------------------------------------
6735
Tot 0.000000 -0.000000 0.000000
6736
----------------------------------------
5929
6738
Res 0.978105 sqrt( Sum f_i^2 / 3N )
5930
6739
----------------------------------------
5931
Max 1.996632 constrained
6740
Max 1.996631 constrained
5933
6742
Stress-tensor-Voigt (kbar): 1053.13 1065.09 1066.80 -6.68 -1.96 -4.38
5934
6743
(Free)E + p*V (eV/cell) -2989.6491
5979
6788
Folding of H and S is implicitly performed
5980
6789
<dSpData1D:S at geom step 49
5981
6790
<sparsity:sparsity for geom step 49
5982
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
5983
<dData1D:(new from dSpData1D) n=306, refcount: 1>
6791
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
6792
<dData1D:(new from dSpData1D) n=204, refcount: 1>
5986
6795
Initializing Density Matrix...
5987
DM after filling with atomic data:
5988
<dSpData2D:(DM initialized from atoms)
6796
DM filled with atomic data:
6797
<dSpData2D:DM initialized from atoms
5989
6798
<sparsity:sparsity for geom step 49
5990
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
5991
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
6799
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
6800
<dData2D:DM n=204 m=1, refcount: 1>
5993
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 49.14116907
6803
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 49.14116906
6804
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
5994
6805
New grid distribution: 1
6813
8 1: 12 7: 12 10: 12
6000
6815
InitMesh: MESH = 24 x 24 x 24 = 13824
6001
6816
InitMesh: (bp) = 12 x 12 x 12 = 1728
6002
6817
InitMesh: Mesh cutoff (required, used) = 50.000 65.517 Ry
6003
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
6818
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
6004
6819
New grid distribution: 2
6009
6828
New grid distribution: 3
6014
6837
Setting up quadratic distribution...
6015
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
6016
PhiOnMesh: Number of (b)points on node 0 = 432
6017
PhiOnMesh: nlist on node 0 = 11797
6019
siesta: Eharris(eV) = -2971.126261
6022
siesta: Eharris(eV) = -2971.126261
6838
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
6839
PhiOnMesh: Number of (b)points on node 0 = 216
6840
PhiOnMesh: nlist on node 0 = 5727
6842
siesta: Eharris(eV) = -2971.126261
6845
siesta: Eharris(eV) = -2971.126261
6847
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
6848
Geom step, scf iteration, dmax: 49 1 1.841396
6025
6850
Using DM_out to compute the final energy and forces
6851
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
6027
6853
siesta: Atomic forces (eV/Ang):
6028
6854
1 0.000000 0.000000 -0.000000
6029
2 -0.000000 0.000000 -0.000000
6030
3 0.000000 0.000000 -0.000000
6031
4 0.000000 0.000000 -0.000000
6032
5 -0.959685 2.146457 -0.000000
6033
6 -0.000000 -1.287767 1.967213
6034
7 0.959685 -0.858690 -1.967213
6035
8 0.959685 -2.146457 -0.000000
6036
9 -0.000000 1.287767 -1.967213
6037
10 -0.959685 0.858690 1.967213
6855
2 -0.000000 0.000000 0.000000
6856
3 0.000001 -0.000000 -0.000001
6857
4 -0.000001 0.000000 0.000001
6858
5 -0.959684 2.146457 -0.000000
6859
6 -0.000000 -1.287768 1.967213
6860
7 0.959684 -0.858690 -1.967214
6861
8 0.959684 -2.146457 0.000000
6862
9 0.000000 1.287768 -1.967213
6863
10 -0.959684 0.858690 1.967214
6038
6864
----------------------------------------
6039
Tot -0.000000 -0.000000 -0.000000
6865
Tot 0.000000 -0.000000 0.000000
6040
6866
----------------------------------------
6042
6868
Res 1.051501 sqrt( Sum f_i^2 / 3N )
6092
6918
Folding of H and S is implicitly performed
6093
6919
<dSpData1D:S at geom step 50
6094
6920
<sparsity:sparsity for geom step 50
6095
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 7>
6096
<dData1D:(new from dSpData1D) n=306, refcount: 1>
6921
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 7>
6922
<dData1D:(new from dSpData1D) n=204, refcount: 1>
6099
6925
Initializing Density Matrix...
6100
DM after filling with atomic data:
6101
<dSpData2D:(DM initialized from atoms)
6926
DM filled with atomic data:
6927
<dSpData2D:DM initialized from atoms
6102
6928
<sparsity:sparsity for geom step 50
6103
nrows_g=34 nrows=9 sparsity=.2647 nnzs=306, refcount: 8>
6104
<dData2D:(DMatomic) n=306 m=1, refcount: 1>
6929
nrows_g=34 nrows=6 sparsity=.1765 nnzs=204, refcount: 8>
6930
<dData2D:DM n=204 m=1, refcount: 1>
6106
6933
Note: For starting DM, Qtot, Tr[D*S] = 48.00000000 49.17551009
6934
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
6107
6935
New grid distribution: 1
6943
8 1: 12 7: 12 10: 12
6113
6945
InitMesh: MESH = 24 x 24 x 24 = 13824
6114
6946
InitMesh: (bp) = 12 x 12 x 12 = 1728
6115
6947
InitMesh: Mesh cutoff (required, used) = 50.000 65.606 Ry
6116
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
6948
ExtMesh (bp) on 0 = 60 x 54 x 51 = 165240
6117
6949
New grid distribution: 2
6122
6958
New grid distribution: 3
6127
6967
Setting up quadratic distribution...
6128
ExtMesh (bp) on 0 = 60 x 54 x 54 = 174960
6129
PhiOnMesh: Number of (b)points on node 0 = 432
6130
PhiOnMesh: nlist on node 0 = 11913
6132
siesta: Eharris(eV) = -2971.183253
6135
siesta: Eharris(eV) = -2971.183253
6968
ExtMesh (bp) on 0 = 54 x 54 x 54 = 157464
6969
PhiOnMesh: Number of (b)points on node 0 = 216
6970
PhiOnMesh: nlist on node 0 = 5813
6972
siesta: Eharris(eV) = -2971.183252
6975
siesta: Eharris(eV) = -2971.183252
6977
SCF_NOT_CONV: SCF did not converge in maximum number of steps.
6978
Geom step, scf iteration, dmax: 50 1 1.844185
6138
6980
Using DM_out to compute the final energy and forces
6981
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps: 10 91
6140
6983
siesta: Atomic forces (eV/Ang):
6141
1 0.000000 0.000000 -0.000000
6142
2 -0.000000 0.000000 -0.000000
6143
3 0.000000 0.000000 -0.000000
6144
4 0.000000 0.000000 -0.000000
6145
5 -1.042792 2.332337 -0.000000
6146
6 -0.000000 -1.399286 2.137572
6147
7 1.042792 -0.933052 -2.137572
6148
8 1.042792 -2.332337 -0.000000
6984
1 0.000000 -0.000000 0.000000
6985
2 0.000000 0.000000 0.000000
6986
3 0.000000 0.000001 -0.000000
6987
4 -0.000000 -0.000001 0.000000
6988
5 -1.042792 2.332337 -0.000001
6989
6 0.000000 -1.399286 2.137572
6990
7 1.042792 -0.933052 -2.137573
6991
8 1.042792 -2.332337 0.000001
6149
6992
9 -0.000000 1.399286 -2.137572
6150
10 -1.042792 0.933052 2.137572
6993
10 -1.042792 0.933052 2.137573
6151
6994
----------------------------------------
6152
Tot -0.000000 -0.000000 -0.000000
6995
Tot 0.000000 -0.000000 0.000000
6153
6996
----------------------------------------
6155
6998
Res 1.142560 sqrt( Sum f_i^2 / 3N )
6165
7008
-0.023378 0.720304 -0.006842
6166
7009
-0.015303 -0.006842 0.726277
6168
siesta: Pressure (static): -1134.91455924 kBar
7011
siesta: Pressure (static): -1134.91453076 kBar
6170
7013
siesta: Stress tensor (total) (eV/Ang**3):
6171
7014
0.672678 -0.014311 -0.009368
6172
7015
-0.014311 0.698287 -0.004188
6173
7016
-0.009368 -0.004188 0.701944
6175
siesta: Pressure (total): -1107.06757282 kBar
6177
siesta: E_kin PR = 116487.371 K
6179
siesta: Temp_ion = 1689.674 K
7018
siesta: Pressure (total): -1107.06754536 kBar
7020
siesta: E_kin PR = 116487.369 K
7022
siesta: Temp_ion = 1689.673 K
6181
7024
outcell: Unit cell vectors (Ang):
6182
7025
5.499877 0.133574 0.087440
6188
7031
outcell: Cell volume (Ang**3) : 76.9031
6190
7033
siesta: Program's energy decomposition (eV):
6191
siesta: Ebs = -736.099978
7034
siesta: Ebs = -736.099975
6192
7035
siesta: Eions = 5228.674091
6193
7036
siesta: Ena = 1092.841594
6194
siesta: Ekin = 2411.129938
7037
siesta: Ekin = 2411.129934
6195
7038
siesta: Enl = -399.387042
6196
7039
siesta: Eso = 0.000000
6197
7040
siesta: Eldau = 0.000000
6198
siesta: DEna = -121.413048
6199
siesta: DUscf = 36.081406
7041
siesta: DEna = -121.413046
7042
siesta: DUscf = 36.081404
6200
7043
siesta: DUext = 0.000000
6201
siesta: Exc = -727.238789
7044
siesta: Exc = -727.238788
6202
7045
siesta: eta*DQ = 0.000000
6203
7046
siesta: Emadel = 0.000000
6204
7047
siesta: Emeta = 0.000000
6205
7048
siesta: Emolmec = 0.000000
6206
7049
siesta: Ekinion = 1.965589
6207
siesta: Eharris = -2969.217664
6208
siesta: Etot = -2934.694444
6209
siesta: FreeEng = -2934.694444
7050
siesta: Eharris = -2969.217663
7051
siesta: Etot = -2934.694445
7052
siesta: FreeEng = -2934.694445
6211
7054
siesta: Final energy (eV):
6212
siesta: Band Struct. = -736.099978
6213
siesta: Kinetic = 2411.129938
6214
siesta: Hartree = 894.739006
7055
siesta: Band Struct. = -736.099975
7056
siesta: Kinetic = 2411.129934
7057
siesta: Hartree = 894.739002
6215
7058
siesta: Eldau = 0.000000
6216
7059
siesta: Eso = 0.000000
6217
7060
siesta: Ext. field = 0.000000
6218
siesta: Exch.-corr. = -727.238789
6219
siesta: Ion-electron = -3750.353251
6220
siesta: Ion-ion = -1764.936936
7061
siesta: Exch.-corr. = -727.238788
7062
siesta: Ion-electron = -3750.353242
7063
siesta: Ion-ion = -1764.936940
6221
7064
siesta: Ekinion = 1.965589
6222
siesta: Total = -2934.694444
7065
siesta: Total = -2934.694445
6223
7066
siesta: Fermi = -5.089759
6225
7068
siesta: Atomic forces (eV/Ang):
6226
siesta: 1 0.000000 0.000000 -0.000000
6227
siesta: 2 -0.000000 0.000000 -0.000000
6228
siesta: 3 0.000000 0.000000 -0.000000
6229
siesta: 4 0.000000 0.000000 -0.000000
6230
siesta: 5 -1.042792 2.332337 -0.000000
6231
siesta: 6 -0.000000 -1.399286 2.137572
6232
siesta: 7 1.042792 -0.933052 -2.137572
6233
siesta: 8 1.042792 -2.332337 -0.000000
7069
siesta: 1 0.000000 -0.000000 0.000000
7070
siesta: 2 0.000000 0.000000 0.000000
7071
siesta: 3 0.000000 0.000001 -0.000000
7072
siesta: 4 -0.000000 -0.000001 0.000000
7073
siesta: 5 -1.042792 2.332337 -0.000001
7074
siesta: 6 0.000000 -1.399286 2.137572
7075
siesta: 7 1.042792 -0.933052 -2.137573
7076
siesta: 8 1.042792 -2.332337 0.000001
6234
7077
siesta: 9 -0.000000 1.399286 -2.137572
6235
siesta: 10 -1.042792 0.933052 2.137572
7078
siesta: 10 -1.042792 0.933052 2.137573
6236
7079
siesta: ----------------------------------------
6237
siesta: Tot -0.000000 -0.000000 -0.000000
7080
siesta: Tot 0.000000 -0.000000 0.000000
6239
7082
siesta: Stress tensor (static) (eV/Ang**3):
6240
7083
siesta: 0.678469 -0.023378 -0.015303