[CPMD-list] MAXIT
Sergei Lisenkov
proffess at yandex.ru
Fri Dec 12 14:57:09 CET 2003
Dear CPMD users,
I have met an error in MD run. From manual I don't know how to solve it:
....
(K+E1+L+N+X) TOTAL ENERGY = -305.18262043 A.U.
(K) KINETIC ENERGY = 307.83507678 A.U.
(E1=A-S+R) ELECTROSTATIC ENERGY = -303.98985038 A.U.
(S) ESELF = 319.15382432 A.U.
(R) ESR = 13.10960226 A.U.
(L) LOCAL PSEUDOPOTENTIAL ENERGY = -242.09751838 A.U.
(N) N-L PSEUDOPOTENTIAL ENERGY = 38.65432287 A.U.
(X) EXCHANGE-CORRELATION ENERGY = -105.58465132 A.U.
GRADIENT CORRECTION ENERGY = -1.09383941 A.U.
NFI EKINC TEMPP EKS ECLASSIC EHAM DIS TCPU
1343.03311 23665.6 -305.18262 -298.66287 44.37024 0.002 26.09
FILE ./TRAJECTORY EXISTS, NEW DATA WILL BE APPENDED
FILE ./TRAJEC.xyz EXISTS, NEW DATA WILL BE APPENDED
FILE ./MOVIE EXISTS, NEW DATA WILL BE APPENDED
RORTHO:DIFGAM= ??????????????????????? ITER = 31 > 30
MAXIMUM NUMBER OF ITERATIONS EXCEEDED
PROGRAM STOPS IN SUBROUTINE RORTOG| [PROC= 0]
And another thing. This calculations was done with Troul.Martins pseudo. I tried also a Vanderbil potential, and I got:
(K+E1+L+N+X) TOTAL ENERGY = 274.90586009 A.U.
(K) KINETIC ENERGY = 222.98588130 A.U.
(E1=A-S+R) ELECTROSTATIC ENERGY = 173.95634078 A.U.
(S) ESELF = 319.15382432 A.U.
(R) ESR = 13.05706950 A.U.
(L) LOCAL PSEUDOPOTENTIAL ENERGY = -77.79274991 A.U.
(N) N-L PSEUDOPOTENTIAL ENERGY = 23.86876994 A.U.
(X) EXCHANGE-CORRELATION ENERGY = -68.11238202 A.U.
GRADIENT CORRECTION ENERGY = -0.22392827 A.U.
NFI EKINC TEMPP EKS ECLASSIC EHAM DIS TCPU
1137.98105 1171.3 274.90586 275.22854 413.20959 0.000 27.45
FILE ./TRAJECTORY EXISTS, NEW DATA WILL BE APPENDED
FILE ./MOVIE EXISTS, NEW DATA WILL BE APPENDED
2498.01012 3789.1 -164.32603 -163.28215 334.72797 0.000 27.33
3449.64236 3870.6 -137.66681 -136.60049 313.04187 0.000 27.27
4297.53068 5889.6 96.62376 98.24631 395.77698 0.000 27.38
5400.23621 9141.1 -63.55281 -61.03448 339.20173 0.001 27.28
6436.88253 11604.6 62.38252 65.57953 502.46206 0.002 27.35
7487.16340 13737.2 8.22484 12.00937 499.17276 0.003 27.22
8482.57376 14948.3 1.25990 5.37808 487.95184 0.004 27.60
9437.18224 17654.0 54.20901 59.07260 496.25484 0.006 27.25
10440.60420 21855.3 39.80961 45.83064 486.43484 0.008 27.31
11457.87573 25703.0 18.48448 25.56552 483.44125 0.011 27.37
12445.60822 29257.3 32.02978 40.09001 485.69823 0.015 27.42
13433.76948 33646.8 44.99931 54.26884 488.03832 0.020 27.31
14443.40019 38383.0 31.58492 42.15926 485.55944 0.026 27.33
...
You can see that process is disconverged. Why?
Thanks,
Sergey
More information about the CPMD-list
mailing list