[CPMD-list] free isolated atom energiesB
Axel Kohlmeyer
axel.kohlmeyer at theochem.ruhr-uni-bochum.de
Fri Sep 17 09:33:34 CEST 2004
On Thu, 16 Sep 2004, Eduardo Ariel Menendez P wrote:
hello eduardo,
AM> Hello everybody:
AM>
AM> During last days I have been involved in generation and testing of
AM> pseudopotentials for Cd and Te. After having performed the tests
stupid question: did you already try the equivalent pseudopotentials
from the goedecker hartwigsen library? at least for LDA there is a very
large collection, including Te and two version of Cd potentials (normal
and semi-core). also the way they are defined seems to make them much
more suitable for those elements than the troullier-martins
pseudopotentials in separable form. at least this is my (limited)
experience so far.
AM> available in the pseudopp generation package (Giannozzi's) and convert the
AM> pseudopotentials to CPMD format, I find that the energies of the isolated
AM> atoms do not coincide precisely. The greatest difference betwwen CPMD and
AM> the atomic
AM> package is for the Te 5p level, where I found a difference of nearly
AM> 0.7eV, while 0.01 difference for Te 5s. I have tried increasing the
AM> cutoff and the lattice constant without improvement.
AM> Hence, I decided to
AM> use a database pseupotencial with the most common element on Earth: O.
AM>
AM> Even so I have two problems
AM>
AM> 1) Using a supplied MT pseudopotential (O_MT_LDA.psp)
AM> I have found .
AM> EIGENVALUES(EV) AND OCCUPATION:
AM> 1 -23.9467297NC 2.000 2 -8.7260317NC 2.000
AM> 3 -8.7260316NC 2.000
please note that these values are _not_ converged as indicated
by the NC after the eigenvalue.
AM> CHEMICAL POTENTIAL = -8.7260316211 EV
AM>
AM> The energies are not precise, compared wih available references:
AM>
AM> http://physics.nist.gov/PhysRefData/DFTdata/Tables/08O.html
AM> 2s -0.871362 Ha= -23.70104 eV
AM> 2p -0.338381 Ha = -9.20 eV
AM>
AM> And the output of the atomic pp generation package og Giannozzi
AM> n l nl e(Ry) (eV)
AM> 1 0 1S( 2.00) -37.5178 -510.4594
AM> 2 0 2S( 2.00) -1.7424 -23.7069
AM> 2 1 2P( 4.00) -0.6766 -9.2055
AM>
AM> E(2P)-E(2S)=14.50eV
AM>
AM> CPMD: E(2P)-E(2S)=15.22eV.
AM>
AM>
AM> Morevoer, if I calculate Kohn-Sham empty state energies I find
AM>
AM> EIGENVALUES(EV) AND OCCUPATION:
AM> 1 -23.9412023 2.000 2 -10.3153247 2.000
AM> 3 -8.7206696 1.000 4 -8.7206695 1.000
AM> 5 -0.5820744 0.000
AM> CHEMICAL POTENTIAL = -8.7206708302 EV
AM>
AM> I would expect that the energies of the states 2,3 and 4 be equal, corresponding to the
AM> three 2p levels. However, the state 2 has a differente. It is interesting that the
AM> occupation factors of 3 and 4 are 1, even when I am using LDA and not LSD. I guess it may be
AM> associated to the quasi-degeneracy.
AM> Here is the input file of the last calculation
AM>
AM> &CPMD
AM> RESTART DENSITY
AM> KOHN-SHAM ENERGIES
AM> 2
AM> MAXSTEP
AM> 150
AM> CENTER MOLECULE ON
AM> &END
AM>
AM> &DFT
AM> NEWCODE
AM> FUNCTIONAL LDA
AM> LDA CORRELATION PZ
AM> SLATER
AM> 0.666667
AM> GC-CUTOFF
AM> 0.1E-07
AM> &END
AM>
AM> &SYSTEM
AM> SYMMETRY
AM> 0
AM> ANGSTROM
AM> CELL
AM> 8.0 1.0 1.0 0 0 0
AM> CUTOFF SPHERICAL
AM> 95.0
AM> &END
AM>
AM> &ATOMS
AM> *O_MT_LDA.psp KLEINMAN-BYLANDER
AM> LMAX=D LOC=D #Also tried with LOC=P and LOC=S
AFAIK, it should be LMAX=P LOCAL=P.
AM> 1
AM> 0.000000 0.000000 0.000000
AM> &END
AM>
AM>
AM>
AM>
AM> 2) SECOND PROBLEM
AM> The second test I did was using ultrasoft pp.
USPPs and LANCZOS are not compatible. see the recent
mail by juerg hutter in the mailing list archives.
finally, i am not so sure, if you can compare the
absolute values, since you are using different grids
for the atomic and the cpmd calculation.
regards,
axel kohlmeyer.
AM> For example,
AM> performing optimization of the wavefunction for one O in a box of 8 Angstroms, and then calculating the Khon-Sham energies I found striking diferences:
AM>
AM> With this entry
AM> &CPMD
AM> RESTART DENSITY
AM> KOHN-SHAM ENERGIES
AM> 0
AM> output:
AM> EIGENVALUES(EV) AND OCCUPATION:
AM> 1 -11.2781946 2.000 2 0.7058206 2.000
AM> 3 0.7065389NC 2.000
AM> CHEMICAL POTENTIAL = 0.7065371502 EV
AM> *********************************
AM> With the same RESTART file and this input
AM> &CPMD
AM> RESTART DENSITY
AM> KOHN-SHAM ENERGIES
AM> 2
AM> y get this outpput:
AM> EIGENVALUES(EV) AND OCCUPATION:
AM> 1 -11.2781946 2.000 2 -0.3815193 2.000
AM> 3 -0.2300328 2.000 4 0.7065208 0.000
AM> 5 0.7065390 0.000
AM> CHEMICAL POTENTIAL = -0.2300326449 EV
AM>
AM> First problem: The energies do not seem to be correct. Now the "P energy is terribly wrong
AM>
AM> http://physics.nist.gov/PhysRefData/DFTdata/Tables/08O.html
AM> 2s -0.871362 Ha= -23.70104 eV
AM> 2p -0.338381 Ha = -9.20 eV
AM>
AM> Also, the states 2,3,and 4 have all diferent energies.
AM>
AM>
AM> Follow full input and output files with the USPP.
AM>
AM> 1) Optimization
AM> &CPMD
AM> INITIALIZE WAVEFUNCTION RANDOM
AM> OPTIMIZE WAVEFUNCTION
AM> MAXSTEP
AM> 150
AM> CENTER MOLECULE ON
AM> &END
AM>
AM>
AM> &DFT
AM> NEWCODE
AM> FUNCTIONAL LDA
AM> LDA CORRELATION PZ
AM> SLATER
AM> 0.666667
AM> GC-CUTOFF
AM> 0.1E-07
AM> &END
AM>
AM> &SYSTEM
AM> SYMMETRY
AM> 0
AM> ANGSTROM
AM> CELL
AM> 8.0 1.0 1.0 0 0 0
AM> CUTOFF SPHERICAL
AM> 55.0
AM> &END
AM>
AM> &ATOMS
AM> *O_VDB_LDA.psp FORMATTED
AM> LMAX=P LOC=P
AM> 1
AM> 0.000000 0.000000 0.000000
AM> &END
AM>
AM> 2) KS orbitals. Only occ states
AM>
AM> &CPMD
AM> RESTART DENSITY
AM> KOHN-SHAM ENERGIES
AM> 0
AM> MAXSTEP
AM> 150
AM> CENTER MOLECULE ON
AM> &END
AM>
AM>
AM> &DFT
AM> NEWCODE
AM> FUNCTIONAL LDA
AM> LDA CORRELATION PZ
AM> SLATER
AM> 0.666667
AM> GC-CUTOFF
AM> 0.1E-07
AM> &END
AM>
AM> &SYSTEM
AM> SYMMETRY
AM> 0
AM> ANGSTROM
AM> CELL
AM> 8.0 1.0 1.0 0 0 0
AM> CUTOFF SPHERICAL
AM> 55.0
AM> &END
AM>
AM> &ATOMS
AM> *O_VDB_LDA.psp FORMATTED
AM> LMAX=D LOC=D
AM> 1
AM> 0.000000 0.000000 0.000000
AM> &END
AM> PROGRAM CPMD STARTED AT: Fri Sep 17 00:46:21 2004
AM>
AM>
AM>
AM> ****** ****** **** **** ******
AM> ******* ******* ********** *******
AM> *** ** *** ** **** ** ** ***
AM> ** ** *** ** ** ** ** **
AM> ** ******* ** ** ** **
AM> *** ****** ** ** ** ***
AM> ******* ** ** ** *******
AM> ****** ** ** ** ******
AM>
AM> VERSION 3.9.1
AM>
AM> COPYRIGHT
AM> IBM RESEARCH DIVISION
AM> MPI FESTKOERPERFORSCHUNG STUTTGART
AM>
AM> The CPMD consortium
AM> WWW: http://www.cpmd.org
AM> Mailinglist: cpmd-list at cpmd.org
AM> E-mail: cpmd at cpmd.org
AM>
AM>
AM> *** Aug 13 2004 -- 13:52:08 ***
AM>
AM> THE INPUT FILE IS: o.inp
AM> THIS JOB RUNS ON:
AM> THE CURRENT DIRECTORY IS:
AM> /home/eariel/mycpmd/O
AM> THE TEMPORARY DIRECTORY IS:
AM> /home/eariel/mycpmd/O
AM> THE PROCESS ID IS: 2530
AM> THE JOB WAS SUBMITTED BY: eariel
AM>
AM>
AM> SINGLE POINT DENSITY OPTIMIZATION
AM>
AM>
AM> EXACT DIAGONALIZATION OF KOHN-SHAM MATRIX
AM>
AM> PATH TO THE RESTART FILES: ./
AM> RESTART WITH OLD DENSITY
AM> GRAM-SCHMIDT ORTHOGONALIZATION
AM> MAXIMUM NUMBER OF STEPS: 1 STEPS
AM> MAXIMUM NUMBER OF ITERATIONS FOR SC: 1 STEPS
AM> PRINT INTERMEDIATE RESULTS EVERY 10001 STEPS
AM> STORE INTERMEDIATE RESULTS EVERY 10001 STEPS
AM> STORE ELECTRONIC DENSITY IN RESTART FILE
AM> NUMBER OF DISTINCT RESTART FILES: 1
AM> TEMPERATURE IS CALCULATED ASSUMING EXTENDED BULK BEHAVIOR
AM> TIME STEP FOR ELECTRONS: 5.0000
AM> TIME STEP FOR IONS: 5.0000
AM> CONVERGENCE CRITERIA FOR WAVEFUNCTION OPTIMIZATION: 1.0000E-03
AM> LANCZOS DIAGONALIZATION (KRYLOV SUBSPACE)
AM> MAX. FRIESNER ITERATIONS 200
AM> MAX. KRYLOV SUBSPACE 6
AM> MAX. KRYLOV BLOCK SIZE FIXED LATER (WAITING FOR NSTATE)
AM> MAX. BETA^2 1.0000E-08
AM> ANDERSON MIXING PARAMETER: 2.0000E-01
AM> BROYDEN MIXING PARAMETER [BROYMIX] 1.5000E-01
AM> BROYDEN CUTOFF [ECUTBROY] EQUAL TO THE DENSITY CUTOFF
AM> BROYDEN MIXING STARTS [NFRBROY] AFTER 0 STEPS
AM> BROYDEN MIXING RESET [IBRESET] AFTER 8 STEPS
AM> BROYDEN MIXING W02 1.0000E-02
AM> ALEXANDER MIXING: 0.9000
AM> SPLINE INTERPOLATION IN G-SPACE FOR PSEUDOPOTENTIAL FUNCTIONS
AM> NUMBER OF SPLINE POINTS: 5000
AM>
AM> EXCHANGE CORRELATION FUNCTIONALS
AM> LDA EXCHANGE: SLATER (ALPHA = 0.66667)
AM> LDA CORRELATION: PERDEW & ZUNGER
AM> [J.P. PERDEW AND A ZUNGER, PRB 23 5048 (1981)]
AM>
AM> *** DETSP| THE NEW SIZE OF THE PROGRAM IS 1940/ 22060 kBYTES ***
AM>
AM> >>>>>>>> CENTER OF MASS HAS BEEN MOVED TO CENTER OF BOX <<<<<<<<
AM>
AM>
AM> ***************************** ATOMS ****************************
AM> NR TYPE X(bohr) Y(bohr) Z(bohr) MBL
AM> 1 O 7.558904 7.558904 7.558904 3
AM> ****************************************************************
AM>
AM> NUMBER OF STATES: 3
AM> NUMBER OF ELECTRONS: 6.00000
AM> CHARGE: 0.00000
AM> ELECTRON TEMPERATURE(KELVIN): 0.00000
AM> OCCUPATION
AM> 2.0 2.0 2.0
AM>
AM> LANCZOS DIAGONALIZATION (KRYLOV SUBSPACE)
AM> MAX. KRYLOV BLOCK SIZE 3
AM>
AM> ============================================================
AM> | pseudopotential report: version 7.3.4 date 6-12-2003 |
AM> ------------------------------------------------------------
AM> | oxygen ceperley-alder exchange-corr |
AM> | z = 8.00 zv = 6.00 exfact = 0.00000 |
AM> | etot = -31.51046 |
AM> | index orbital occupation energy |
AM> | 1 200 2.00 -1.74 |
AM> | 2 210 4.00 -0.68 |
AM> | keyps = 3 ifpcor = 0 |
AM> | rinner = 0.80 for L= 1 |
AM> | rinner = 0.80 for L= 2 |
AM> | rinner = 0.80 for L= 3 |
AM> | new generation scheme: |
AM> | nbeta = 4 kkbeta = 507 rcloc = 1.0000 |
AM> | ibeta l epsilon rcut iptype |
AM> | 1 0 -1.74 1.10 2 |
AM> | 2 0 -0.20 1.10 2 |
AM> | 3 1 -0.80 1.10 2 |
AM> | 4 1 0.20 1.10 2 |
AM> | npf = 8 ptryc = 10.000 |
AM> | lloc = 2 eloc = 0.000 |
AM> | ifqopt = 2 nqf = 8 qtryc = 10.000 |
AM> | all electron calculation used schroedinger equation |
AM> | ************logarithmic mesh************ |
AM> ============================================================
AM>
AM> ****************************************************************
AM> * ATOM MASS RAGGIO NLCC PSEUDOPOTENTIAL *
AM> * O 15.9994 1.2000 NO VANDERBILT *
AM> ****************************************************************
AM>
AM>
AM> OPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPEN
AM> NUMBER OF CPUS PER TASK 1
AM> OPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPEN
AM>
AM> *** RGGEN| THE NEW SIZE OF THE PROGRAM IS 8656/ 36592 kBYTES ***
AM>
AM> ************************** SUPERCELL ***************************
AM> THIS IS AN ISOLATED SYSTEM CALCULATION
AM> POISSON EQUATION SOLVER : HOCKNEY
AM> COULOMB SMOOTHING RADIUS : 1.080
AM> SYMMETRY: SIMPLE CUBIC
AM> LATTICE CONSTANT(a.u.): 15.11781
AM> CELL DIMENSION: 15.1178 1.0000 1.0000 0.0000 0.0000 0.0000
AM> VOLUME(OMEGA IN BOHR^3): 3455.14651
AM> LATTICE VECTOR A1(BOHR): 15.1178 0.0000 0.0000
AM> LATTICE VECTOR A2(BOHR): 0.0000 15.1178 0.0000
AM> LATTICE VECTOR A3(BOHR): 0.0000 0.0000 15.1178
AM> RECIP. LAT. VEC. B1(2Pi/BOHR): 0.0661 0.0000 0.0000
AM> RECIP. LAT. VEC. B2(2Pi/BOHR): 0.0000 0.0661 0.0000
AM> RECIP. LAT. VEC. B3(2Pi/BOHR): 0.0000 0.0000 0.0661
AM> REAL SPACE MESH: 72 72 72
AM> WAVEFUNCTION CUTOFF(RYDBERG): 55.00000
AM> DENSITY CUTOFF(RYDBERG): (DUAL= 4.00) 220.00000
AM> NUMBER OF PLANE WAVES FOR WAVEFUNCTION CUTOFF: 11924
AM> NUMBER OF PLANE WAVES FOR DENSITY CUTOFF: 95160
AM> NUMBER OF PLANE WAVES FOR BROYDEN MIXING: 95160
AM> ****************************************************************
AM>
AM> *** RINFORCE| THE NEW SIZE OF THE PROGRAM IS 19336/ 48876 kBYTES ***
AM> *** FFTPRP| THE NEW SIZE OF THE PROGRAM IS 30400/ 50196 kBYTES ***
AM> *** CLUSTER| THE NEW SIZE OF THE PROGRAM IS 34304/ 53980 kBYTES ***
AM>
AM> GENERATE ATOMIC BASIS SET
AM> O SLATER ORBITALS
AM> 2S ALPHA= 2.2458 OCCUPATION= 2.00
AM> 2P ALPHA= 2.2266 OCCUPATION= 4.00
AM>
AM>
AM> INITIALIZATION TIME: 13.06 SECONDS
AM>
AM> *** WFOPTS| THE NEW SIZE OF THE PROGRAM IS 34520/ 60124 kBYTES ***
AM> RV30| WARNING! NO DENSITY INFORMATION
AM>
AM> RESTART INFORMATION READ ON FILE ./RESTART
AM> *** PHFAC| THE NEW SIZE OF THE PROGRAM IS 34780/ 89428 kBYTES ***
AM> *** ATOMWF| THE NEW SIZE OF THE PROGRAM IS 34984/ 90924 kBYTES ***
AM> ATRHO| CHARGE(R-SPACE): 6.000000 (G-SPACE): 6.000000
AM>
AM> ****************************************************************
AM> * ATOMIC COORDINATES *
AM> ****************************************************************
AM> 1 O 7.558904 7.558904 7.558904
AM> ****************************************************************
AM>
AM> !! WARNING FROM DETDOF: WHY IS GLIB LESS THAN UNITY !!
AM> !! WARNING FROM DETDOF: NOW IT IS SET EQUAL TO UNITY !!
AM>
AM> DEGREES OF FREEDOM FOR SYSTEM: 1
AM>
AM> CPU TIME FOR WAVEFUNCTION INITIALIZATION: 3.44 SECONDS
AM> *** RWFOPT| THE NEW SIZE OF THE PROGRAM IS 64920/ 89564 kBYTES ***
AM> EWALD| SUM IN REAL SPACE OVER 1* 1* 1 CELLS
AM>
AM> <<<<<<<<<<<<<<<<<<< LANCZOS DIAGONALIZATION <<<<<<<<<<<<<<<<<<<<
AM> >> TIME FOR INITIAL SUBSPACE DIAGONALIZATION: 0.32
AM> >> CYCLE NCONV B2MAX B2MIN #HPSI TIME
AM> 1 0 1.561E-01 4.184E-02 6.00 1.58
AM> 2 0 1.562E-01 8.824E-03 5.00 1.57
AM> 3 0 6.681E-02 2.637E-03 5.00 1.58
AM> 4 0 5.009E-02 8.258E-04 5.00 1.57
AM> 5 0 2.734E-02 2.903E-04 5.00 1.58
AM> 6 0 2.382E-02 1.033E-04 5.00 1.57
AM> 7 0 1.440E-02 3.971E-05 5.00 1.58
AM> 8 0 1.340E-02 1.528E-05 5.00 1.57
AM> 9 0 8.514E-03 4.825E-06 5.00 1.58
AM> 10 0 8.230E-03 2.850E-06 5.00 1.57
AM> 11 0 5.415E-03 1.034E-06 5.00 1.57
AM> 12 0 5.390E-03 6.196E-07 5.00 1.57
AM> 13 0 3.644E-03 2.334E-07 5.00 1.58
AM> 14 0 3.711E-03 1.416E-07 5.00 1.57
AM> 15 0 2.561E-03 5.407E-08 5.00 1.57
AM> 16 0 2.654E-03 3.308E-08 5.00 1.58
AM> 17 0 1.859E-03 1.274E-08 5.00 1.57
AM> 18 1 1.950E-03 9.159E-09 4.00 1.13
AM> 19 1 1.380E-03 9.159E-09 3.33 0.82
AM> 20 1 1.459E-03 9.159E-09 3.33 0.83
AM> 21 1 1.039E-03 9.159E-09 3.33 0.83
AM> 22 1 1.104E-03 9.159E-09 3.33 0.82
AM> 23 1 7.904E-04 9.159E-09 3.33 0.83
AM> 24 1 8.425E-04 9.159E-09 3.33 0.82
AM> 25 1 6.047E-04 9.159E-09 3.33 0.83
AM> 26 1 6.459E-04 9.159E-09 3.33 0.83
AM> 27 1 4.644E-04 9.159E-09 3.33 0.82
AM> 28 1 4.966E-04 9.159E-09 3.33 0.83
AM> 29 1 3.574E-04 9.159E-09 3.33 0.83
AM> 30 1 3.824E-04 9.159E-09 3.33 0.82
AM> 31 1 2.754E-04 9.159E-09 3.33 0.83
AM> 32 1 2.948E-04 9.159E-09 3.33 0.83
AM> 33 1 2.124E-04 9.159E-09 3.33 0.82
AM> 34 1 2.274E-04 9.159E-09 3.33 0.83
AM> 35 1 1.639E-04 9.159E-09 3.33 0.82
AM> 36 1 1.755E-04 9.159E-09 3.33 0.83
AM> 37 1 1.265E-04 9.159E-09 3.33 0.83
AM> 38 1 1.354E-04 9.159E-09 3.33 0.82
AM> 39 1 9.763E-05 9.159E-09 3.33 0.83
AM> 40 1 1.045E-04 9.159E-09 3.33 0.83
AM> 41 1 7.535E-05 9.159E-09 3.33 0.82
AM> 42 1 8.067E-05 9.159E-09 3.33 0.83
AM> 43 1 5.816E-05 9.159E-09 3.33 0.82
AM> 44 1 6.227E-05 9.159E-09 3.33 0.83
AM> 45 1 4.489E-05 9.159E-09 3.33 0.83
AM> 46 1 4.806E-05 9.159E-09 3.33 0.82
AM> 47 1 3.465E-05 9.159E-09 3.33 0.83
AM> 48 1 3.709E-05 9.159E-09 3.33 0.82
AM> 49 1 2.675E-05 9.159E-09 3.33 0.82
AM> 50 1 2.863E-05 9.159E-09 3.33 0.83
AM> 51 1 2.064E-05 9.159E-09 3.33 0.82
AM> 52 1 2.210E-05 9.159E-09 3.33 0.83
AM> 53 1 1.593E-05 9.159E-09 3.33 0.83
AM> 54 1 1.706E-05 9.159E-09 3.33 0.82
AM> 55 1 1.230E-05 9.159E-09 3.33 0.83
AM> 56 1 1.317E-05 9.159E-09 3.33 0.82
AM> 57 1 9.494E-06 9.159E-09 3.33 0.83
AM> 58 1 1.016E-05 9.159E-09 3.33 0.83
AM> 59 1 7.329E-06 9.159E-09 3.33 0.82
AM> 60 1 7.844E-06 9.159E-09 3.33 0.82
AM> 61 1 5.657E-06 9.159E-09 3.33 0.83
AM> 62 1 6.055E-06 9.159E-09 3.33 0.82
AM> 63 1 4.367E-06 9.159E-09 3.33 0.83
AM> 64 1 4.674E-06 9.159E-09 3.33 0.83
AM> 65 1 3.371E-06 9.159E-09 3.33 0.82
AM> 66 1 3.609E-06 9.159E-09 3.33 0.82
AM> 67 1 2.603E-06 9.159E-09 3.33 0.82
AM> 68 1 2.786E-06 9.159E-09 3.33 0.83
AM> 69 1 2.010E-06 9.159E-09 3.33 0.83
AM> 70 1 2.151E-06 9.159E-09 3.33 0.82
AM> 71 1 1.552E-06 9.159E-09 3.33 0.83
AM> 72 1 1.661E-06 9.159E-09 3.33 0.82
AM> 73 1 1.198E-06 9.159E-09 3.33 0.83
AM> 74 1 1.283E-06 9.159E-09 3.33 0.83
AM> 75 1 9.256E-07 9.159E-09 3.33 0.83
AM> 76 1 9.910E-07 9.159E-09 3.33 0.82
AM> 77 1 7.151E-07 9.159E-09 3.33 0.83
AM> 78 1 7.657E-07 9.159E-09 3.33 0.82
AM> 79 1 5.526E-07 9.159E-09 3.33 0.83
AM> 80 1 5.919E-07 9.159E-09 3.33 0.82
AM> 81 1 4.273E-07 9.159E-09 3.33 0.83
AM> 82 1 4.578E-07 9.159E-09 3.33 0.82
AM> 83 1 3.307E-07 9.159E-09 3.33 0.82
AM> 84 1 3.544E-07 9.159E-09 3.33 0.83
AM> 85 1 2.561E-07 9.159E-09 3.33 0.82
AM> 86 1 2.748E-07 9.159E-09 3.33 0.83
AM> 87 1 1.987E-07 9.159E-09 3.33 0.83
AM> 88 1 2.134E-07 9.159E-09 3.33 0.82
AM> 89 1 1.545E-07 9.159E-09 3.33 0.83
AM> 90 1 1.662E-07 9.159E-09 3.33 0.83
AM> 91 1 1.205E-07 9.159E-09 3.33 0.82
AM> 92 1 1.299E-07 9.159E-09 3.33 0.83
AM> 93 1 9.449E-08 9.159E-09 3.33 0.82
AM> 94 1 1.021E-07 9.159E-09 3.33 0.83
AM> 95 1 7.456E-08 9.159E-09 3.33 0.83
AM> 96 1 8.096E-08 9.159E-09 3.33 0.82
AM> 97 1 5.940E-08 9.159E-09 3.33 0.83
AM> 98 1 6.490E-08 9.159E-09 3.33 0.83
AM> 99 1 4.796E-08 9.159E-09 3.33 0.82
AM> 100 1 5.283E-08 9.159E-09 3.33 0.83
AM> 101 1 3.942E-08 9.159E-09 3.33 0.82
AM> 102 1 4.391E-08 9.159E-09 3.33 0.83
AM> 103 1 3.317E-08 9.159E-09 3.33 0.83
AM> 104 1 3.747E-08 9.159E-09 3.33 0.82
AM> 105 1 2.874E-08 9.159E-09 3.33 0.83
AM> 106 1 3.302E-08 9.159E-09 3.33 0.83
AM> 107 2 2.578E-08 9.159E-09 2.33 0.78
AM> 108 2 3.019E-08 9.159E-09 1.67 0.75
AM> 109 2 2.403E-08 9.159E-09 1.67 0.75
AM> 110 2 2.871E-08 9.159E-09 1.67 0.75
AM> 111 2 2.330E-08 9.159E-09 1.67 0.75
AM> 112 2 2.838E-08 9.159E-09 1.67 0.76
AM> 113 2 2.345E-08 9.159E-09 1.67 0.75
AM> 114 2 2.906E-08 9.159E-09 1.67 0.75
AM> 115 2 2.439E-08 9.159E-09 1.67 0.75
AM> 116 2 3.066E-08 9.159E-09 1.67 0.75
AM> 117 2 2.606E-08 9.159E-09 1.67 0.75
AM> 118 2 3.313E-08 9.159E-09 1.67 0.75
AM> 119 2 2.843E-08 9.159E-09 1.67 0.75
AM> 120 2 3.644E-08 9.159E-09 1.67 0.76
AM> 121 2 3.149E-08 9.159E-09 1.67 0.75
AM> 122 2 4.062E-08 9.159E-09 1.67 0.75
AM> 123 2 3.527E-08 9.159E-09 1.67 0.75
AM> 124 2 4.569E-08 9.159E-09 1.67 0.75
AM> 125 2 3.981E-08 9.159E-09 1.67 0.75
AM> 126 2 5.172E-08 9.159E-09 1.67 0.75
AM> 127 2 4.518E-08 9.159E-09 1.67 0.75
AM> 128 2 5.880E-08 9.159E-09 1.67 0.76
AM> 129 2 5.145E-08 9.159E-09 1.67 0.75
AM> 130 2 6.705E-08 9.159E-09 1.67 0.75
AM> 131 2 5.873E-08 9.159E-09 1.67 0.75
AM> 132 2 7.660E-08 9.159E-09 1.67 0.75
AM> 133 2 6.715E-08 9.159E-09 1.67 0.75
AM> 134 2 8.763E-08 9.159E-09 1.67 0.75
AM> 135 2 7.685E-08 9.159E-09 1.67 0.75
AM> 136 2 1.003E-07 9.159E-09 1.67 0.74
AM> 137 2 8.803E-08 9.159E-09 1.67 0.76
AM> 138 2 1.150E-07 9.159E-09 1.67 0.75
AM> 139 2 1.009E-07 9.159E-09 1.67 0.75
AM> 140 2 1.317E-07 9.159E-09 1.67 0.75
AM> 141 2 1.156E-07 9.159E-09 1.67 0.75
AM> 142 2 1.510E-07 9.159E-09 1.67 0.75
AM> 143 2 1.326E-07 9.159E-09 1.67 0.75
AM> 144 2 1.732E-07 9.159E-09 1.67 0.75
AM> 145 2 1.520E-07 9.159E-09 1.67 0.74
AM> 146 2 1.986E-07 9.159E-09 1.67 0.75
AM> 147 2 1.743E-07 9.159E-09 1.67 0.75
AM> 148 2 2.277E-07 9.159E-09 1.67 0.75
AM> 149 2 1.999E-07 9.159E-09 1.67 0.75
AM> 150 2 2.611E-07 9.159E-09 1.67 0.75
AM> 151 2 2.293E-07 9.159E-09 1.67 0.75
AM> 152 2 2.995E-07 9.159E-09 1.67 0.75
AM> 153 2 2.629E-07 9.159E-09 1.67 0.75
AM> 154 2 3.434E-07 9.159E-09 1.67 0.75
AM> 155 2 3.015E-07 9.159E-09 1.67 0.75
AM> 156 2 3.939E-07 9.159E-09 1.67 0.75
AM> 157 2 3.458E-07 9.159E-09 1.67 0.75
AM> 158 2 4.517E-07 9.159E-09 1.67 0.75
AM> 159 2 3.966E-07 9.159E-09 1.67 0.75
AM> 160 2 5.181E-07 9.159E-09 1.67 0.75
AM> 161 2 4.549E-07 9.159E-09 1.67 0.75
AM> 162 2 5.941E-07 9.159E-09 1.67 0.75
AM> 163 2 5.217E-07 9.159E-09 1.67 0.75
AM> 164 2 6.814E-07 9.159E-09 1.67 0.75
AM> 165 2 5.983E-07 9.159E-09 1.67 0.75
AM> 166 2 7.815E-07 9.159E-09 1.67 0.75
AM> 167 2 6.862E-07 9.159E-09 1.67 0.75
AM> 168 2 8.963E-07 9.159E-09 1.67 0.75
AM> 169 2 7.870E-07 9.159E-09 1.67 0.75
AM> 170 2 1.028E-06 9.159E-09 1.67 0.75
AM> 171 2 9.026E-07 9.159E-09 1.67 0.75
AM> 172 2 1.179E-06 9.159E-09 1.67 0.75
AM> 173 2 1.035E-06 9.159E-09 1.67 0.75
AM> 174 2 1.352E-06 9.159E-09 1.67 0.76
AM> 175 2 1.187E-06 9.159E-09 1.67 0.75
AM> 176 2 1.551E-06 9.159E-09 1.67 0.75
AM> 177 2 1.362E-06 9.159E-09 1.67 0.75
AM> 178 2 1.778E-06 9.159E-09 1.67 0.75
AM> 179 2 1.562E-06 9.159E-09 1.67 0.75
AM> 180 2 2.039E-06 9.159E-09 1.67 0.75
AM> 181 2 1.791E-06 9.159E-09 1.67 0.75
AM> 182 2 2.339E-06 9.159E-09 1.67 0.76
AM> 183 2 2.054E-06 9.159E-09 1.67 0.74
AM> 184 2 2.682E-06 9.159E-09 1.67 0.74
AM> 185 2 2.356E-06 9.159E-09 1.67 0.75
AM> 186 2 3.076E-06 9.159E-09 1.67 0.75
AM> 187 2 2.701E-06 9.159E-09 1.67 0.75
AM> 188 2 3.528E-06 9.159E-09 1.67 0.74
AM> 189 2 3.098E-06 9.159E-09 1.67 0.75
AM> 190 2 4.046E-06 9.159E-09 1.67 0.75
AM> 191 2 3.553E-06 9.159E-09 1.67 0.75
AM> 192 2 4.639E-06 9.159E-09 1.67 0.75
AM> 193 2 4.074E-06 9.159E-09 1.67 0.75
AM> 194 2 5.320E-06 9.159E-09 1.67 0.75
AM> 195 2 4.672E-06 9.159E-09 1.67 0.74
AM> 196 2 6.101E-06 9.159E-09 1.67 0.75
AM> 197 2 5.358E-06 9.159E-09 1.67 0.75
AM> 198 2 6.996E-06 9.159E-09 1.67 0.76
AM> 199 2 6.144E-06 9.159E-09 1.67 0.75
AM> 200 2 8.023E-06 9.159E-09 1.67 0.75
AM>
AM> EIGENVALUES(EV) AND OCCUPATION:
AM> 1 -11.2781946 2.000 2 0.7058206 2.000
AM> 3 0.7065389NC 2.000
AM> CHEMICAL POTENTIAL = 0.7065371502 EV
AM>
AM> TOTAL INTEGRATED ELECTRONIC DENSITY
AM> IN G-SPACE = 6.049092
AM> IN R-SPACE = 6.049092
AM>
AM>
AM> VANDERBILT AUGMENTATION CHARGES (MEAN VALUE PER ATOM)
AM> ATOM TYPE NR. OF ATOMS CHARGE
AM> O 1 0.049
AM>
AM> (B+E2+X-V) TOTAL ENERGY = -16.81175790 A.U.
AM> (B) BAND ENERGY = -0.72510914 A.U.
AM> (E2=I-H-S+R) ELECTROSTATIC ENERGY = -17.33584728 A.U.
AM> (S) ESELF = 11.96826841 A.U.
AM> (R) ESR = 0.00000000 A.U.
AM> (X) EXCHANGE-CORRELATION ENERGY = -3.99099801 A.U.
AM> (V) EXCHANGE-CORRELATION POTEN. = -5.24019667 A.U.
AM>
AM> ==------------------------------------------------------------==
AM> == NFI= 1 ETOT= -16.811758 TCPU= 173.96 ==
AM> == DRHOMAX= 1.407E-01 DETOT= 0.000E+00 THL= 0.000E+00 ==
AM> ==------------------------------------------------------------==
AM>
AM> RESTART INFORMATION WRITTEN ON FILE ./RESTART.1
AM> *** RWFOPT| THE NEW SIZE OF THE PROGRAM IS 75640/ 95644 kBYTES ***
AM>
AM> ****************************************************************
AM> * *
AM> * FINAL RESULTS *
AM> * *
AM> ****************************************************************
AM>
AM> ****************************************************************
AM> * ATOMIC COORDINATES *
AM> ****************************************************************
AM> 1 O 7.558904 7.558904 7.558904
AM> ****************************************************************
AM>
AM>
AM> EIGENVALUES(EV) AND OCCUPATION:
AM> 1 -11.2781946 2.000 2 0.7058206 2.000
AM> 3 0.7065389NC 2.000
AM> CHEMICAL POTENTIAL = 0.7065371502 EV
AM>
AM> ****************************************************************
AM>
AM>
AM> ELECTRONIC GRADIENT:
AM> MAX. COMPONENT = 1.40719E-01 NORM = 6.04874E+00
AM>
AM> TOTAL INTEGRATED ELECTRONIC DENSITY
AM> IN G-SPACE = 6.049092
AM> IN R-SPACE = 6.049092
AM>
AM>
AM> VANDERBILT AUGMENTATION CHARGES (MEAN VALUE PER ATOM)
AM> ATOM TYPE NR. OF ATOMS CHARGE
AM> O 1 0.049
AM>
AM> (B+E2+X-V) TOTAL ENERGY = -16.81175790 A.U.
AM> (B) BAND ENERGY = -0.72510914 A.U.
AM> (E2=I-H-S+R) ELECTROSTATIC ENERGY = -17.33584728 A.U.
AM> (S) ESELF = 11.96826841 A.U.
AM> (R) ESR = 0.00000000 A.U.
AM> (X) EXCHANGE-CORRELATION ENERGY = -3.99099801 A.U.
AM> (V) EXCHANGE-CORRELATION POTEN. = -5.24019667 A.U.
AM>
AM> ****************************************************************
AM>
AM>
AM> ================================================================
AM> BIG MEMORY ALLOCATIONS
AM> SCR 1576657 PSI 1556068
AM> QRL 903600 YLMB 856440
AM> YF 778034 XF 778034
AM> SCR 613400 PME 429264
AM> RHOE 389017 HGPOT 389017
AM> ----------------------------------------------------------------
AM> [PEAK NUMBER 92] PEAK MEMORY 9571930 = 76.6 MBytes
AM> ================================================================
AM>
AM>
AM> ****************************************************************
AM> * *
AM> * TIMING *
AM> * *
AM> ****************************************************************
AM> SUBROUTINE CALLS CPU TIME ELAPSED TIME
AM> FFT-G/S 4372 39.07 38.96
AM> VPSI 1002 37.95 38.13
AM> S_FWFFT 1091 37.65 37.75
AM> S_INVFFT 1095 36.25 38.53
AM> QINIT 1 10.51 10.51
AM> FNONLOC 1002 8.31 8.30
AM> SPSI 1001 5.24 5.17
AM> RNLSM1 1006 3.48 3.46
AM> FRIESNER 1 2.61 2.65
AM> VDBINIT 1 1.54 1.54
AM> RHOV 2 1.08 1.09
AM> HIP 2 1.06 1.08
AM> NEWD 2 1.01 1.01
AM> HPSI 1001 0.80 0.80
AM> ATRHO 1 0.62 0.63
AM> HOCKNEY 1 0.60 0.63
AM> FWFFT 8 0.47 0.46
AM> ----------------------------------------------------------------
AM> TOTAL TIME 188.25 190.71
AM> ****************************************************************
AM>
AM> CPU TIME : 0 HOURS 3 MINUTES 10.52 SECONDS
AM> ELAPSED TIME : 0 HOURS 3 MINUTES 13.71 SECONDS
AM>
AM> PROGRAM CPMD ENDED AT: Fri Sep 17 00:49:35 2004
AM>
AM> 3) KS energies with empty states
AM>
AM> &CPMD
AM> RESTART DENSITY
AM> KOHN-SHAM ENERGIES
AM> 2
AM> MAXSTEP
AM> 150
AM> CENTER MOLECULE ON
AM> &END
AM>
AM>
AM> &DFT
AM> NEWCODE
AM> FUNCTIONAL LDA
AM> LDA CORRELATION PZ
AM> SLATER
AM> 0.666667
AM> GC-CUTOFF
AM> 0.1E-07
AM> &END
AM>
AM> &SYSTEM
AM> SYMMETRY
AM> 0
AM> ANGSTROM
AM> CELL
AM> 8.0 1.0 1.0 0 0 0
AM> CUTOFF SPHERICAL
AM> 55.0
AM> &END
AM>
AM> &ATOMS
AM> *O_VDB_LDA.psp FORMATTED
AM> LMAX=D LOC=D
AM> 1
AM> 0.000000 0.000000 0.000000
AM> &END
AM> PROGRAM CPMD STARTED AT: Fri Sep 17 00:51:40 2004
AM>
AM>
AM>
AM> ****** ****** **** **** ******
AM> ******* ******* ********** *******
AM> *** ** *** ** **** ** ** ***
AM> ** ** *** ** ** ** ** **
AM> ** ******* ** ** ** **
AM> *** ****** ** ** ** ***
AM> ******* ** ** ** *******
AM> ****** ** ** ** ******
AM>
AM> VERSION 3.9.1
AM>
AM> COPYRIGHT
AM> IBM RESEARCH DIVISION
AM> MPI FESTKOERPERFORSCHUNG STUTTGART
AM>
AM> The CPMD consortium
AM> WWW: http://www.cpmd.org
AM> Mailinglist: cpmd-list at cpmd.org
AM> E-mail: cpmd at cpmd.org
AM>
AM>
AM> *** Aug 13 2004 -- 13:52:08 ***
AM>
AM> THE INPUT FILE IS: o.inp-ks
AM> THIS JOB RUNS ON:
AM> THE CURRENT DIRECTORY IS:
AM> /home/eariel/mycpmd/O
AM> THE TEMPORARY DIRECTORY IS:
AM> /home/eariel/mycpmd/O
AM> THE PROCESS ID IS: 2546
AM> THE JOB WAS SUBMITTED BY: eariel
AM>
AM>
AM> SINGLE POINT DENSITY OPTIMIZATION
AM>
AM>
AM> EXACT DIAGONALIZATION OF KOHN-SHAM MATRIX
AM>
AM> PATH TO THE RESTART FILES: ./
AM> RESTART WITH OLD DENSITY
AM> GRAM-SCHMIDT ORTHOGONALIZATION
AM> MAXIMUM NUMBER OF STEPS: 1 STEPS
AM> MAXIMUM NUMBER OF ITERATIONS FOR SC: 1 STEPS
AM> PRINT INTERMEDIATE RESULTS EVERY 10001 STEPS
AM> STORE INTERMEDIATE RESULTS EVERY 10001 STEPS
AM> STORE ELECTRONIC DENSITY IN RESTART FILE
AM> NUMBER OF DISTINCT RESTART FILES: 1
AM> TEMPERATURE IS CALCULATED ASSUMING EXTENDED BULK BEHAVIOR
AM> TIME STEP FOR ELECTRONS: 5.0000
AM> TIME STEP FOR IONS: 5.0000
AM> CONVERGENCE CRITERIA FOR WAVEFUNCTION OPTIMIZATION: 1.0000E-03
AM> LANCZOS DIAGONALIZATION (KRYLOV SUBSPACE)
AM> MAX. FRIESNER ITERATIONS 200
AM> MAX. KRYLOV SUBSPACE 6
AM> MAX. KRYLOV BLOCK SIZE FIXED LATER (WAITING FOR NSTATE)
AM> MAX. BETA^2 1.0000E-08
AM> ANDERSON MIXING PARAMETER: 2.0000E-01
AM> BROYDEN MIXING PARAMETER [BROYMIX] 1.5000E-01
AM> BROYDEN CUTOFF [ECUTBROY] EQUAL TO THE DENSITY CUTOFF
AM> BROYDEN MIXING STARTS [NFRBROY] AFTER 0 STEPS
AM> BROYDEN MIXING RESET [IBRESET] AFTER 8 STEPS
AM> BROYDEN MIXING W02 1.0000E-02
AM> ALEXANDER MIXING: 0.9000
AM> SPLINE INTERPOLATION IN G-SPACE FOR PSEUDOPOTENTIAL FUNCTIONS
AM> NUMBER OF SPLINE POINTS: 5000
AM>
AM> EXCHANGE CORRELATION FUNCTIONALS
AM> LDA EXCHANGE: SLATER (ALPHA = 0.66667)
AM> LDA CORRELATION: PERDEW & ZUNGER
AM> [J.P. PERDEW AND A ZUNGER, PRB 23 5048 (1981)]
AM>
AM> *** DETSP| THE NEW SIZE OF THE PROGRAM IS 1940/ 22060 kBYTES ***
AM>
AM> >>>>>>>> CENTER OF MASS HAS BEEN MOVED TO CENTER OF BOX <<<<<<<<
AM>
AM>
AM> ***************************** ATOMS ****************************
AM> NR TYPE X(bohr) Y(bohr) Z(bohr) MBL
AM> 1 O 7.558904 7.558904 7.558904 3
AM> ****************************************************************
AM>
AM> NUMBER OF STATES: 5
AM> NUMBER OF ELECTRONS: 6.00000
AM> CHARGE: 0.00000
AM> ELECTRON TEMPERATURE(KELVIN): 0.00000
AM> OCCUPATION
AM> 2.0 2.0 2.0 0.0 0.0
AM>
AM> LANCZOS DIAGONALIZATION (KRYLOV SUBSPACE)
AM> MAX. KRYLOV BLOCK SIZE 5
AM>
AM> ============================================================
AM> | pseudopotential report: version 7.3.4 date 6-12-2003 |
AM> ------------------------------------------------------------
AM> | oxygen ceperley-alder exchange-corr |
AM> | z = 8.00 zv = 6.00 exfact = 0.00000 |
AM> | etot = -31.51046 |
AM> | index orbital occupation energy |
AM> | 1 200 2.00 -1.74 |
AM> | 2 210 4.00 -0.68 |
AM> | keyps = 3 ifpcor = 0 |
AM> | rinner = 0.80 for L= 1 |
AM> | rinner = 0.80 for L= 2 |
AM> | rinner = 0.80 for L= 3 |
AM> | new generation scheme: |
AM> | nbeta = 4 kkbeta = 507 rcloc = 1.0000 |
AM> | ibeta l epsilon rcut iptype |
AM> | 1 0 -1.74 1.10 2 |
AM> | 2 0 -0.20 1.10 2 |
AM> | 3 1 -0.80 1.10 2 |
AM> | 4 1 0.20 1.10 2 |
AM> | npf = 8 ptryc = 10.000 |
AM> | lloc = 2 eloc = 0.000 |
AM> | ifqopt = 2 nqf = 8 qtryc = 10.000 |
AM> | all electron calculation used schroedinger equation |
AM> | ************logarithmic mesh************ |
AM> ============================================================
AM>
AM> ****************************************************************
AM> * ATOM MASS RAGGIO NLCC PSEUDOPOTENTIAL *
AM> * O 15.9994 1.2000 NO VANDERBILT *
AM> ****************************************************************
AM>
AM>
AM> OPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPEN
AM> NUMBER OF CPUS PER TASK 1
AM> OPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPEN
AM>
AM> *** RGGEN| THE NEW SIZE OF THE PROGRAM IS 8656/ 36592 kBYTES ***
AM>
AM> ************************** SUPERCELL ***************************
AM> THIS IS AN ISOLATED SYSTEM CALCULATION
AM> POISSON EQUATION SOLVER : HOCKNEY
AM> COULOMB SMOOTHING RADIUS : 1.080
AM> SYMMETRY: SIMPLE CUBIC
AM> LATTICE CONSTANT(a.u.): 15.11781
AM> CELL DIMENSION: 15.1178 1.0000 1.0000 0.0000 0.0000 0.0000
AM> VOLUME(OMEGA IN BOHR^3): 3455.14651
AM> LATTICE VECTOR A1(BOHR): 15.1178 0.0000 0.0000
AM> LATTICE VECTOR A2(BOHR): 0.0000 15.1178 0.0000
AM> LATTICE VECTOR A3(BOHR): 0.0000 0.0000 15.1178
AM> RECIP. LAT. VEC. B1(2Pi/BOHR): 0.0661 0.0000 0.0000
AM> RECIP. LAT. VEC. B2(2Pi/BOHR): 0.0000 0.0661 0.0000
AM> RECIP. LAT. VEC. B3(2Pi/BOHR): 0.0000 0.0000 0.0661
AM> REAL SPACE MESH: 72 72 72
AM> WAVEFUNCTION CUTOFF(RYDBERG): 55.00000
AM> DENSITY CUTOFF(RYDBERG): (DUAL= 4.00) 220.00000
AM> NUMBER OF PLANE WAVES FOR WAVEFUNCTION CUTOFF: 11924
AM> NUMBER OF PLANE WAVES FOR DENSITY CUTOFF: 95160
AM> NUMBER OF PLANE WAVES FOR BROYDEN MIXING: 95160
AM> ****************************************************************
AM>
AM> *** RINFORCE| THE NEW SIZE OF THE PROGRAM IS 19336/ 48876 kBYTES ***
AM> *** FFTPRP| THE NEW SIZE OF THE PROGRAM IS 30400/ 50196 kBYTES ***
AM> *** CLUSTER| THE NEW SIZE OF THE PROGRAM IS 34304/ 53980 kBYTES ***
AM>
AM> GENERATE ATOMIC BASIS SET
AM> O SLATER ORBITALS
AM> 2S ALPHA= 2.2458 OCCUPATION= 2.00
AM> 2P ALPHA= 2.2266 OCCUPATION= 4.00
AM>
AM>
AM> INITIALIZATION TIME: 13.06 SECONDS
AM>
AM> *** WFOPTS| THE NEW SIZE OF THE PROGRAM IS 34520/ 63288 kBYTES ***
AM> RV30! NUMBER OF STATES HAS CHANGED 3 5
AM> RV30| WARNING! NO DENSITY INFORMATION
AM>
AM> RESTART INFORMATION READ ON FILE ./RESTART
AM> *** PHFAC| THE NEW SIZE OF THE PROGRAM IS 34780/ 92592 kBYTES ***
AM> *** ATOMWF| THE NEW SIZE OF THE PROGRAM IS 34984/ 94088 kBYTES ***
AM> ATRHO| WARNING! THE NUMBER OF GENERATED ATOMIC STATE (4)
AM> ATRHO| WARNING! IS LESS THAN THE NUMBER OF STATE (5)
AM> ATRHO| WARNING! 1 STATES ARE EQUAL TO ZERO
AM> ATRHO| USE SECTION BASIS TO SPECIFY ATOMIC BASIS
AM> ATRHO| CHARGE(R-SPACE): 6.000000 (G-SPACE): 6.000000
AM> ATOMWF| WARNING! RANDOMIZATION OF 1 STATES
AM>
AM> ****************************************************************
AM> * ATOMIC COORDINATES *
AM> ****************************************************************
AM> 1 O 7.558904 7.558904 7.558904
AM> ****************************************************************
AM>
AM> !! WARNING FROM DETDOF: WHY IS GLIB LESS THAN UNITY !!
AM> !! WARNING FROM DETDOF: NOW IT IS SET EQUAL TO UNITY !!
AM>
AM> DEGREES OF FREEDOM FOR SYSTEM: 1
AM>
AM> CPU TIME FOR WAVEFUNCTION INITIALIZATION: 3.41 SECONDS
AM> *** RWFOPT| THE NEW SIZE OF THE PROGRAM IS 65308/ 92728 kBYTES ***
AM> EWALD| SUM IN REAL SPACE OVER 1* 1* 1 CELLS
AM>
AM> <<<<<<<<<<<<<<<<<<< LANCZOS DIAGONALIZATION <<<<<<<<<<<<<<<<<<<<
AM> >> TIME FOR INITIAL SUBSPACE DIAGONALIZATION: 0.49
AM> >> CYCLE NCONV B2MAX B2MIN #HPSI TIME
AM> 1 0 1.561E-01 4.181E-02 6.00 2.38
AM> 2 0 1.563E-01 5.488E-03 5.00 2.39
AM> 3 0 6.679E-02 2.159E-03 5.00 2.40
AM> 4 0 5.013E-02 7.902E-04 5.00 2.39
AM> 5 0 2.733E-02 3.050E-04 5.00 2.40
AM> 6 0 2.384E-02 9.732E-05 5.00 2.40
AM> 7 0 1.440E-02 3.847E-05 5.00 2.39
AM> 8 0 1.341E-02 1.197E-05 5.00 2.39
AM> 9 0 8.518E-03 5.429E-06 5.00 2.40
AM> 10 0 8.238E-03 1.632E-06 5.00 2.39
AM> 11 0 5.415E-03 8.208E-07 5.00 2.40
AM> 12 0 5.396E-03 2.484E-07 5.00 2.39
AM> 13 0 3.645E-03 1.290E-07 5.00 2.40
AM> 14 0 3.716E-03 3.949E-08 5.00 2.39
AM> 15 0 2.562E-03 2.070E-08 5.00 2.39
AM> 16 1 2.657E-03 8.884E-09 4.40 1.96
AM> 17 1 1.860E-03 8.884E-09 4.00 1.65
AM> 18 1 1.952E-03 8.884E-09 4.00 1.66
AM> 19 1 1.380E-03 8.884E-09 4.00 1.65
AM> 20 1 1.460E-03 8.884E-09 4.00 1.65
AM> 21 1 1.040E-03 8.884E-09 4.00 1.66
AM> 22 1 1.106E-03 8.884E-09 4.00 1.65
AM> 23 1 7.907E-04 8.884E-09 4.00 1.65
AM> 24 1 8.436E-04 8.884E-09 4.00 1.65
AM> 25 1 6.049E-04 8.884E-09 4.00 1.66
AM> 26 1 6.467E-04 8.884E-09 4.00 1.65
AM> 27 1 4.645E-04 8.884E-09 4.00 1.65
AM> 28 1 4.972E-04 8.884E-09 4.00 1.65
AM> 29 1 3.575E-04 8.884E-09 4.00 1.66
AM> 30 1 3.830E-04 8.884E-09 4.00 1.65
AM> 31 1 2.755E-04 8.884E-09 4.00 1.65
AM> 32 1 2.953E-04 8.884E-09 4.00 1.66
AM> 33 1 2.125E-04 8.884E-09 4.00 1.65
AM> 34 1 2.278E-04 8.884E-09 4.00 1.65
AM> 35 1 1.640E-04 8.884E-09 4.00 1.65
AM> 36 1 1.758E-04 8.884E-09 4.00 1.66
AM> 37 1 1.265E-04 8.884E-09 4.00 1.65
AM> 38 1 1.356E-04 8.884E-09 4.00 1.65
AM> 39 1 9.766E-05 8.884E-09 4.00 1.65
AM> 40 1 1.047E-04 8.884E-09 4.00 1.66
AM> 41 1 7.538E-05 8.884E-09 4.00 1.65
AM> 42 1 8.081E-05 8.884E-09 4.00 1.65
AM> 43 1 5.818E-05 8.884E-09 4.00 1.66
AM> 44 1 6.237E-05 8.884E-09 4.00 1.65
AM> 45 1 4.491E-05 8.884E-09 4.00 1.65
AM> 46 1 4.814E-05 8.884E-09 4.00 1.65
AM> 47 1 3.466E-05 8.884E-09 4.00 1.66
AM> 48 1 3.716E-05 8.884E-09 4.00 1.65
AM> 49 1 2.675E-05 8.884E-09 4.00 1.65
AM> 50 2 2.868E-05 8.884E-09 3.80 1.64
AM> 51 2 2.065E-05 8.884E-09 3.00 1.58
AM> 52 2 2.214E-05 8.884E-09 3.00 1.57
AM> 53 2 1.594E-05 8.884E-09 3.00 1.57
AM> 54 2 1.709E-05 8.884E-09 3.00 1.58
AM> 55 2 1.230E-05 8.884E-09 3.00 1.57
AM> 56 2 1.319E-05 8.884E-09 3.00 1.57
AM> 57 2 9.497E-06 8.884E-09 3.00 1.57
AM> 58 2 1.018E-05 8.884E-09 3.00 1.58
AM> 59 2 7.330E-06 8.884E-09 3.00 1.57
AM> 60 2 7.860E-06 8.884E-09 3.00 1.58
AM> 61 2 5.658E-06 8.884E-09 3.00 1.57
AM> 62 2 6.067E-06 8.884E-09 3.00 1.58
AM> 63 2 4.368E-06 8.884E-09 3.00 1.57
AM> 64 2 4.684E-06 8.884E-09 3.00 1.58
AM> 65 2 3.372E-06 8.884E-09 3.00 1.57
AM> 66 2 3.616E-06 8.884E-09 3.00 1.57
AM> 67 2 2.603E-06 8.884E-09 3.00 1.57
AM> 68 2 2.791E-06 8.884E-09 3.00 1.57
AM> 69 2 2.009E-06 8.884E-09 3.00 1.58
AM> 70 2 2.155E-06 8.884E-09 3.00 1.57
AM> 71 2 1.551E-06 8.884E-09 3.00 1.57
AM> 72 2 1.664E-06 8.884E-09 3.00 1.57
AM> 73 2 1.198E-06 8.884E-09 3.00 1.58
AM> 74 2 1.284E-06 8.884E-09 3.00 1.57
AM> 75 2 9.245E-07 8.884E-09 3.00 1.58
AM> 76 3 9.915E-07 8.884E-09 2.40 1.13
AM> 77 3 7.137E-07 8.884E-09 2.00 0.83
AM> 78 3 7.655E-07 8.884E-09 2.00 0.84
AM> 79 3 5.511E-07 8.884E-09 2.00 0.83
AM> 80 3 5.910E-07 8.884E-09 2.00 0.83
AM> 81 3 4.254E-07 8.884E-09 2.00 0.83
AM> 82 3 4.563E-07 8.884E-09 2.00 0.83
AM> 83 3 3.285E-07 8.884E-09 2.00 0.84
AM> 84 3 3.523E-07 8.884E-09 2.00 0.83
AM> 85 3 2.536E-07 8.884E-09 2.00 0.83
AM> 86 3 2.721E-07 8.884E-09 2.00 0.84
AM> 87 3 1.958E-07 8.884E-09 2.00 0.83
AM> 88 3 2.101E-07 8.884E-09 2.00 0.83
AM> 89 3 1.512E-07 8.884E-09 2.00 0.84
AM> 90 3 1.622E-07 8.884E-09 2.00 0.83
AM> 91 3 1.167E-07 8.884E-09 2.00 0.84
AM> 92 3 1.252E-07 8.884E-09 2.00 0.83
AM> 93 3 9.014E-08 8.884E-09 2.00 0.83
AM> 94 3 9.671E-08 8.884E-09 2.00 0.84
AM> 95 3 6.960E-08 8.884E-09 2.00 0.83
AM> 96 3 7.468E-08 8.884E-09 2.00 0.84
AM> 97 3 5.374E-08 8.884E-09 2.00 0.83
AM> 98 3 5.766E-08 8.884E-09 2.00 0.83
AM> 99 3 4.150E-08 8.884E-09 2.00 0.84
AM> 100 3 4.453E-08 8.884E-09 2.00 0.83
AM> 101 3 3.204E-08 8.884E-09 2.00 0.83
AM> 102 3 3.439E-08 8.884E-09 2.00 0.84
AM> 103 3 2.474E-08 8.884E-09 2.00 0.84
AM> 104 3 2.655E-08 8.884E-09 2.00 0.83
AM> 105 3 1.911E-08 8.884E-09 2.00 0.83
AM> 106 3 2.050E-08 8.884E-09 2.00 0.83
AM> 107 5 9.530E-09 8.884E-09 0.80 0.35
AM>
AM> EIGENVALUES(EV) AND OCCUPATION:
AM> 1 -11.2781946 2.000 2 -0.3815193 2.000
AM> 3 -0.2300328 2.000 4 0.7065208 0.000
AM> 5 0.7065390 0.000
AM> CHEMICAL POTENTIAL = -0.2300326449 EV
AM>
AM> TOTAL INTEGRATED ELECTRONIC DENSITY
AM> IN G-SPACE = 6.034344
AM> IN R-SPACE = 6.034344
AM>
AM>
AM> VANDERBILT AUGMENTATION CHARGES (MEAN VALUE PER ATOM)
AM> ATOM TYPE NR. OF ATOMS CHARGE
AM> O 1 0.034
AM>
AM> (B+E2+X-V) TOTAL ENERGY = -16.96050918 A.U.
AM> (B) BAND ENERGY = -0.87386055 A.U.
AM> (E2=I-H-S+R) ELECTROSTATIC ENERGY = -17.33584728 A.U.
AM> (S) ESELF = 11.96826841 A.U.
AM> (R) ESR = 0.00000000 A.U.
AM> (X) EXCHANGE-CORRELATION ENERGY = -3.99099801 A.U.
AM> (V) EXCHANGE-CORRELATION POTEN. = -5.24019667 A.U.
AM>
AM> ==------------------------------------------------------------==
AM> == NFI= 1 ETOT= -16.960509 TCPU= 162.84 ==
AM> == DRHOMAX= 1.439E-01 DETOT= 0.000E+00 THL= 0.000E+00 ==
AM> ==------------------------------------------------------------==
AM>
AM> RESTART INFORMATION WRITTEN ON FILE ./RESTART.1
AM> *** RWFOPT| THE NEW SIZE OF THE PROGRAM IS 78992/ 98804 kBYTES ***
AM>
AM> ****************************************************************
AM> * *
AM> * FINAL RESULTS *
AM> * *
AM> ****************************************************************
AM>
AM> ****************************************************************
AM> * ATOMIC COORDINATES *
AM> ****************************************************************
AM> 1 O 7.558904 7.558904 7.558904
AM> ****************************************************************
AM>
AM>
AM> EIGENVALUES(EV) AND OCCUPATION:
AM> 1 -11.2781946 2.000 2 -0.3815193 2.000
AM> 3 -0.2300328 2.000 4 0.7065208 0.000
AM> 5 0.7065390 0.000
AM> CHEMICAL POTENTIAL = -0.2300326449 EV
AM>
AM> ****************************************************************
AM>
AM>
AM> ELECTRONIC GRADIENT:
AM> MAX. COMPONENT = 1.43915E-01 NORM = 6.89928E+00
AM>
AM> TOTAL INTEGRATED ELECTRONIC DENSITY
AM> IN G-SPACE = 6.034344
AM> IN R-SPACE = 6.034344
AM>
AM>
AM> VANDERBILT AUGMENTATION CHARGES (MEAN VALUE PER ATOM)
AM> ATOM TYPE NR. OF ATOMS CHARGE
AM> O 1 0.034
AM>
AM> (B+E2+X-V) TOTAL ENERGY = -16.96050918 A.U.
AM> (B) BAND ENERGY = -0.87386055 A.U.
AM> (E2=I-H-S+R) ELECTROSTATIC ENERGY = -17.33584728 A.U.
AM> (S) ESELF = 11.96826841 A.U.
AM> (R) ESR = 0.00000000 A.U.
AM> (X) EXCHANGE-CORRELATION ENERGY = -3.99099801 A.U.
AM> (V) EXCHANGE-CORRELATION POTEN. = -5.24019667 A.U.
AM>
AM> ****************************************************************
AM>
AM>
AM> ================================================================
AM> BIG MEMORY ALLOCATIONS
AM> SCR 1576657 PSI 1556068
AM> QRL 903600 YLMB 856440
AM> YF 778034 XF 778034
AM> PME 715440 SCR 613400
AM> RHOE 389017 HGPOT 389017
AM> ----------------------------------------------------------------
AM> [PEAK NUMBER 92] PEAK MEMORY 9977414 = 79.8 MBytes
AM> ================================================================
AM>
AM>
AM> ****************************************************************
AM> * *
AM> * TIMING *
AM> * *
AM> ****************************************************************
AM> SUBROUTINE CALLS CPU TIME ELAPSED TIME
AM> FFT-G/S 3972 35.52 35.46
AM> VPSI 534 33.82 33.85
AM> S_FWFFT 991 33.41 33.52
AM> S_INVFFT 995 32.97 32.92
AM> FNONLOC 534 11.30 11.34
AM> QINIT 1 10.55 10.55
AM> SPSI 533 5.33 5.38
AM> RNLSM1 538 3.45 3.43
AM> FRIESNER 1 3.23 3.24
AM> VDBINIT 1 1.41 1.41
AM> RHOV 2 1.08 1.07
AM> HPSI 533 1.03 1.01
AM> HIP 2 1.02 1.03
AM> NEWD 2 0.98 0.98
AM> HOCKNEY 1 0.66 0.67
AM> ATRHO 1 0.57 0.59
AM> FWFFT 8 0.45 0.46
AM> OVLAP2 110 0.45 0.42
AM> INVFFT 7 0.40 0.40
AM> ----------------------------------------------------------------
AM> TOTAL TIME 177.63 177.74
AM> ****************************************************************
AM>
AM> CPU TIME : 0 HOURS 2 MINUTES 59.40 SECONDS
AM> ELAPSED TIME : 0 HOURS 3 MINUTES 1.25 SECONDS
AM>
AM> PROGRAM CPMD ENDED AT: Fri Sep 17 00:54:41 2004
AM>
AM> Other options that I have tried:
AM>
AM> &CPMD
AM> OPTIMIZE WAVEFUNCTION
AM> # PCG MINIMIZE
AM> # TIMESTEP
AM> # 5
AM> RESTART ALL
AM> LANCZOS DIAGONALISATION ALL
AM>
AM> &DFT
AM> GC-CUTOFF
AM> 0.1E-05 # Sugeted for ultrasoft PP
AM>
AM> POISSON SOLVER TUCKERMAN
AM>
AM> What is wrong?
AM>
AM> Thank you
AM>
AM> Eduardo A. Menendez Proupin
AM> Department of Physics
AM> Faculty of Science
AM> University of Chile
AM> Las Palmeras 3425
AM> Ñuñoa, Santiago
AM> Chile
AM> Phone: 56+2+678 74 11
AM> http://fisica.ciencias.uchile.cl/~emenendez/
AM>
AM> _______________________________________________
AM> CPMD-list mailing list
AM> CPMD-list at cpmd.org
AM> http://cpmd.org/mailman/listinfo/cpmd-list
AM>
AM>
--
=======================================================================
Dr. Axel Kohlmeyer e-mail: axel.kohlmeyer at rub.de
Lehrstuhl fuer Theoretische Chemie Phone: ++49 (0)234/32-26673
Ruhr-Universitaet Bochum - NC 03/53 Fax: ++49 (0)234/32-14045
D-44780 Bochum http://www.theochem.ruhr-uni-bochum.de/~axel.kohlmeyer/
=======================================================================
More information about the CPMD-list
mailing list