PROGRAM CPMD STARTED AT: Wed Jan 5 18:39:43 2005 ****** ****** **** **** ****** ******* ******* ********** ******* *** ** *** ** **** ** ** *** ** ** *** ** ** ** ** ** ** ******* ** ** ** ** *** ****** ** ** ** *** ******* ** ** ** ******* ****** ** ** ** ****** VERSION 3.9.2 COPYRIGHT IBM RESEARCH DIVISION MPI FESTKOERPERFORSCHUNG STUTTGART The CPMD consortium WWW: http://www.cpmd.org Mailinglist: cpmd-list@cpmd.org E-mail: cpmd@cpmd.org *** Jan 5 2005 -- 12:51:29 *** THE INPUT FILE IS: 3-gly-md-cont.inp THIS JOB RUNS ON: opt13 THE CURRENT DIRECTORY IS: /tmp_mnt/rubberbandman/akohlmey/webpage/myhome/cpmd-tutor/reference/04-gly THE TEMPORARY DIRECTORY IS: /tmp_mnt/rubberbandman/akohlmey/webpage/myhome/cpmd-tutor/reference/04-gly THE PROCESS ID IS: 23540 THE JOB WAS SUBMITTED BY: akohlmey ****************************************************************************** * INFO - INFO - INFO - INFO - INFO - INFO - INFO - INFO - INFO - INFO - INFO * ****************************************************************************** * MD production of a glycine in gas phase. * * no thermostat. * ****************************************************************************** CAR-PARRINELLO MOLECULAR DYNAMICS PATH TO THE RESTART FILES: ./ RESTART WITH OLD ORBITALS RESTART WITH OLD ION POSITIONS RESTART WITH OLD VELOCITIES ITERATIVE ORTHOGONALIZATION MAXIT: 30 EPS: 1.00E-06 MAXIMUM NUMBER OF STEPS: 500 STEPS PRINT INTERMEDIATE RESULTS EVERY 10001 STEPS STORE INTERMEDIATE RESULTS EVERY 10001 STEPS STORE INTERMEDIATE RESULTS EVERY 501 SELF-CONSISTENT STEPS NUMBER OF DISTINCT RESTART FILES: 1 TEMPERATURE IS CALCULATED ASSUMING EXTENDED BULK BEHAVIOR FICTITIOUS ELECTRON MASS: 400.0000 TIME STEP FOR ELECTRONS: 4.0000 TIME STEP FOR IONS: 4.0000 TRAJECTORIES ARE SAVED ON FILE EVERY 10 STEPS TRAJEC.xyz IS SAVED ON FILE EVERY 10 STEPS ELECTRON DYNAMICS: THE TEMPERATURE IS NOT CONTROLLED ION DYNAMICS: THE TEMPERATURE IS NOT CONTROLLED SUBTRACT CENTER OF MASS VELOCITY EVERY 25 STEPS SUBTRACT ROTATION AROUND CENTER OF MASS EVERY 25 STEPS SPLINE INTERPOLATION IN G-SPACE FOR PSEUDOPOTENTIAL FUNCTIONS NUMBER OF SPLINE POINTS: 2000 EXCHANGE CORRELATION FUNCTIONALS LDA EXCHANGE: SLATER (ALPHA = 0.66667) LDA CORRELATION: LEE, YANG & PARR [C.L. LEE, W. YANG, AND R.G. PARR, PRB 37 785 (1988)] GRADIENT CORRECTED FUNCTIONAL DENSITY THRESHOLD: 5.00000E-06 EXCHANGE ENERGY [A.D. BECKE, PHYS. REV. A 38, 3098 (1988)] PARAMETER BETA: 0.004200 CORRELATION ENERGY [LYP: C.L. LEE ET AL. PHYS. REV. B 37, 785 (1988)] *** DETSP| THE NEW SIZE OF THE PROGRAM IS 1716/ 49776 kBYTES *** >>>>>>>> CENTER OF MASS HAS BEEN MOVED TO CENTER OF BOX <<<<<<<< ***************************** ATOMS **************************** NR TYPE X(bohr) Y(bohr) Z(bohr) MBL 1 N 10.987436 13.131875 13.468706 3 2 C 13.438036 12.178637 14.828613 3 3 C 15.268374 14.655943 14.420429 3 4 O 17.266557 14.715263 15.633500 3 5 O 14.338897 16.170299 12.691970 3 6 H 11.969727 14.427276 12.175075 3 7 H 10.193959 14.506482 14.645546 3 8 H 9.725479 11.809643 12.888917 3 9 H 14.214852 10.617428 13.738530 3 10 H 12.991638 11.667640 16.764726 3 **************************************************************** NUMBER OF STATES: 15 NUMBER OF ELECTRONS: 30.00000 CHARGE: 0.00000 ELECTRON TEMPERATURE(KELVIN): 0.00000 OCCUPATION 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 ============================================================ | pseudopotential report: version 7.3.4 date 10-23-2003 | ------------------------------------------------------------ | nitrogen C-A + B88gx + LYPgc exchange-corr | | z = 7.00 zv = 5.00 exfact = 1.00000 | | etot = -19.28646 | | index orbital occupation energy | | 1 200 2.00 -1.35 | | 2 210 3.00 -0.52 | | keyps = 3 ifpcor = 0 | | rinner = 0.80 for L= 1 | | rinner = 0.80 for L= 2 | | rinner = 0.80 for L= 3 | | new generation scheme: | | nbeta = 4 kkbeta = 499 rcloc = 1.0000 | | ibeta l epsilon rcut iptype | | 1 0 -1.50 1.10 2 | | 2 0 0.00 1.10 2 | | 3 1 -0.70 1.10 2 | | 4 1 0.20 1.10 2 | | npf = 8 ptryc = 10.000 | | lloc = 2 eloc = 0.000 | | ifqopt = 3 nqf = 8 qtryc = 10.000 | | all electron calculation used schroedinger equation | | ************logarithmic mesh************ | ============================================================ ============================================================ | pseudopotential report: version 7.3.4 date 10-23-2003 | ------------------------------------------------------------ | carbon C-A + B88gx + LYPgc exchange-corr | | z = 6.00 zv = 4.00 exfact = 1.00000 | | etot = -10.70828 | | index orbital occupation energy | | 1 200 2.00 -1.00 | | 2 210 2.00 -0.38 | | keyps = 3 ifpcor = 0 | | rinner = 0.80 for L= 1 | | rinner = 0.80 for L= 2 | | rinner = 0.80 for L= 3 | | new generation scheme: | | nbeta = 4 kkbeta = 491 rcloc = 1.0000 | | ibeta l epsilon rcut iptype | | 1 0 -1.20 1.10 2 | | 2 0 0.20 1.10 2 | | 3 1 -0.60 1.10 2 | | 4 1 0.20 1.10 2 | | npf = 8 ptryc = 10.000 | | lloc = 2 eloc = 0.000 | | ifqopt = 3 nqf = 8 qtryc = 10.000 | | all electron calculation used schroedinger equation | | ************logarithmic mesh************ | ============================================================ ============================================================ | pseudopotential report: version 7.3.4 date 2- 3-2004 | ------------------------------------------------------------ | oxygen C-A + B88gx + LYPgc exchange-corr | | z = 8.00 zv = 6.00 exfact = 1.00000 | | etot = -31.57898 | | index orbital occupation energy | | 1 200 2.00 -1.75 | | 2 210 4.00 -0.66 | | keyps = 3 ifpcor = 0 | | rinner = 0.80 for L= 1 | | rinner = 0.80 for L= 2 | | rinner = 0.80 for L= 3 | | new generation scheme: | | nbeta = 4 kkbeta = 507 rcloc = 1.0000 | | ibeta l epsilon rcut iptype | | 1 0 -1.75 1.20 2 | | 2 0 -0.20 1.20 2 | | 3 1 -0.80 1.20 2 | | 4 1 0.20 1.20 2 | | npf = 8 ptryc = 10.000 | | lloc = 2 eloc = 0.000 | | ifqopt = 3 nqf = 8 qtryc = 10.000 | | all electron calculation used schroedinger equation | | ************logarithmic mesh************ | ============================================================ ============================================================ | pseudopotential report: version 7.3.4 date 10-23-2003 | ------------------------------------------------------------ | hydrogen C-A + B88gx + LYPgc exchange-corr | | z = 1.00 zv = 1.00 exfact = 1.00000 | | etot = -0.92437 | | index orbital occupation energy | | 1 100 1.00 -0.48 | | keyps = 3 ifpcor = 0 | | rinner = 0.80 for L= 1 | | new generation scheme: | | nbeta = 1 kkbeta = 371 rcloc = 0.8000 | | ibeta l epsilon rcut iptype | | 1 0 -0.48 0.80 2 | | npf = 6 ptryc = 10.000 | | lloc = 1 eloc = 0.000 | | ifqopt = 3 nqf = 8 qtryc = 10.000 | | all electron calculation used schroedinger equation | | ************logarithmic mesh************ | ============================================================ **************************************************************** * ATOM MASS RAGGIO NLCC PSEUDOPOTENTIAL * * N 14.0067 1.2000 NO VANDERBILT * * C 12.0112 1.2000 NO VANDERBILT * * O 15.9994 1.2000 NO VANDERBILT * * H 1.0080 1.2000 NO VANDERBILT * **************************************************************** OPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPEN NUMBER OF CPUS PER TASK 1 OPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPEN *** RGGEN| THE NEW SIZE OF THE PROGRAM IS 19608/ 79524 kBYTES *** ************************** SUPERCELL *************************** THIS IS AN ISOLATED SYSTEM CALCULATION POISSON EQUATION SOLVER : TUCKERMAN & MARTYNA SHORT RANGE POTENTIAL LENGTH * BOX LENGTH 7.000 SYMMETRY: SIMPLE CUBIC LATTICE CONSTANT(a.u.): 28.34589 CELL DIMENSION: 28.3459 1.0000 1.0000 0.0000 0.0000 0.0000 VOLUME(OMEGA IN BOHR^3): 22775.62400 LATTICE VECTOR A1(BOHR): 28.3459 0.0000 0.0000 LATTICE VECTOR A2(BOHR): 0.0000 28.3459 0.0000 LATTICE VECTOR A3(BOHR): 0.0000 0.0000 28.3459 RECIP. LAT. VEC. B1(2Pi/BOHR): 0.0353 0.0000 0.0000 RECIP. LAT. VEC. B2(2Pi/BOHR): 0.0000 0.0353 0.0000 RECIP. LAT. VEC. B3(2Pi/BOHR): 0.0000 0.0000 0.0353 REAL SPACE MESH: 108 108 108 WAVEFUNCTION CUTOFF(RYDBERG): 25.00000 DENSITY CUTOFF(RYDBERG): (DUAL= 5.00) 125.00000 NUMBER OF PLANE WAVES FOR WAVEFUNCTION CUTOFF: 23969 NUMBER OF PLANE WAVES FOR DENSITY CUTOFF: 268677 **************************************************************** *** RINFORCE| THE NEW SIZE OF THE PROGRAM IS 65440/ 134796 kBYTES *** *** FFTPRP| THE NEW SIZE OF THE PROGRAM IS 103560/ 159344 kBYTES *** *** CLUSTER| THE NEW SIZE OF THE PROGRAM IS 107820/ 163544 kBYTES *** GENERATE ATOMIC BASIS SET N SLATER ORBITALS 2S ALPHA= 1.9237 OCCUPATION= 2.00 2P ALPHA= 1.9170 OCCUPATION= 3.00 C SLATER ORBITALS 2S ALPHA= 1.6083 OCCUPATION= 2.00 2P ALPHA= 1.5679 OCCUPATION= 2.00 O SLATER ORBITALS 2S ALPHA= 2.2458 OCCUPATION= 2.00 2P ALPHA= 2.2266 OCCUPATION= 4.00 H SLATER ORBITALS 1S ALPHA= 1.0000 OCCUPATION= 1.00 INITIALIZATION TIME: 11.20 SECONDS *** MDPT| THE NEW SIZE OF THE PROGRAM IS 108296/ 187288 kBYTES *** RESTART INFORMATION READ ON FILE ./RESTART *** PHFAC| THE NEW SIZE OF THE PROGRAM IS 119784/ 318164 kBYTES *** **************************************************************** * ATOMIC COORDINATES * **************************************************************** 1 N 10.892391 13.108611 13.738458 2 C 13.385995 12.215069 14.668606 3 C 15.243338 14.430710 14.431737 4 O 17.171458 14.823837 15.694407 5 O 14.524857 16.220199 12.542552 6 H 12.858927 15.733331 12.121383 7 H 9.600584 13.385277 15.133760 8 H 10.133713 11.696908 12.608707 9 H 13.954923 10.563967 13.517303 10 H 13.358951 11.749606 16.646955 **************************************************************** DEGREES OF FREEDOM FOR SYSTEM: 27 EWALD| SUM IN REAL SPACE OVER 1* 1* 1 CELLS **************************************************************** * ATOMIC COORDINATES * **************************************************************** 1 N 10.892391 13.108611 13.738458 2 C 13.385995 12.215069 14.668606 3 C 15.243338 14.430710 14.431737 4 O 17.171458 14.823837 15.694407 5 O 14.524857 16.220199 12.542552 6 H 12.858927 15.733331 12.121383 7 H 9.600584 13.385277 15.133760 8 H 10.133713 11.696908 12.608707 9 H 13.954923 10.563967 13.517303 10 H 13.358951 11.749606 16.646955 **************************************************************** CPU TIME FOR INITIALIZATION: 16.99 SECONDS ATOM COORDINATES GRADIENTS (-FORCES) 1 N 10.8909 13.1080 13.7388 -1.124E-02 -2.193E-02 -1.845E-02 2 C 13.3859 12.2143 14.6696 -1.308E-03 1.336E-02 -7.439E-03 3 C 15.2428 14.4303 14.4311 2.946E-02 3.339E-02 -7.922E-03 4 O 17.1712 14.8245 15.6937 -2.614E-02 -1.567E-02 -1.428E-02 5 O 14.5269 16.2215 12.5430 4.437E-02 -8.740E-03 2.956E-02 6 H 12.8578 15.7333 12.1210 -4.102E-02 -9.396E-03 -1.260E-02 7 H 9.6010 13.3831 15.1321 -1.313E-03 3.199E-03 7.289E-03 8 H 10.1352 11.6948 12.6085 -5.713E-05 8.985E-03 3.480E-03 9 H 13.9544 10.5608 13.5163 5.301E-03 5.521E-03 5.341E-03 10 H 13.3587 11.7498 16.6457 4.087E-04 -9.293E-03 1.262E-02 TOTAL INTEGRATED ELECTRONIC DENSITY IN G-SPACE = 30.000000 IN R-SPACE = 30.000000 VANDERBILT AUGMENTATION CHARGES (MEAN VALUE PER ATOM) ATOM TYPE NR. OF ATOMS CHARGE N 1 0.304 C 2 0.105 O 2 0.809 H 5 0.032 (K+E1+L+N+X) TOTAL ENERGY = -56.22279580 A.U. (K) KINETIC ENERGY = 28.27123896 A.U. (E1=A-S+R) ELECTROSTATIC ENERGY = -41.18446946 A.U. (S) ESELF = 44.54855464 A.U. (R) ESR = 2.32493600 A.U. (L) LOCAL PSEUDOPOTENTIAL ENERGY = -41.90108558 A.U. (N) N-L PSEUDOPOTENTIAL ENERGY = 13.21713324 A.U. (X) EXCHANGE-CORRELATION ENERGY = -14.62561296 A.U. GRADIENT CORRECTION ENERGY = -0.73486940 A.U. NFI EKINC TEMPP EKS ECLASSIC EHAM DIS TCPU 1 0.00081 275.0 -56.22280 -56.21104 -56.21024 0.406E-05 16.94 2 0.00081 276.5 -56.22287 -56.21105 -56.21024 0.162E-04 17.00 3 0.00082 278.1 -56.22295 -56.21106 -56.21024 0.363E-04 16.99 4 0.00083 279.8 -56.22303 -56.21107 -56.21024 0.643E-04 16.93 5 0.00084 281.5 -56.22311 -56.21108 -56.21024 0.100E-03 16.96 6 0.00085 283.2 -56.22319 -56.21109 -56.21024 0.144E-03 16.94 7 0.00086 284.8 -56.22327 -56.21110 -56.21024 0.196E-03 16.95 8 0.00087 286.2 -56.22334 -56.21111 -56.21024 0.256E-03 16.96 9 0.00088 287.6 -56.22341 -56.21112 -56.21024 0.323E-03 17.07 10 0.00089 288.8 -56.22347 -56.21112 -56.21024 0.399E-03 17.14 11 0.00089 289.8 -56.22352 -56.21113 -56.21024 0.482E-03 16.97 12 0.00090 290.7 -56.22356 -56.21114 -56.21024 0.573E-03 17.21 13 0.00090 291.3 -56.22359 -56.21114 -56.21024 0.672E-03 16.94 14 0.00090 291.7 -56.22361 -56.21114 -56.21024 0.779E-03 16.97 15 0.00091 292.0 -56.22362 -56.21114 -56.21024 0.892E-03 16.94 16 0.00091 292.0 -56.22363 -56.21114 -56.21024 0.101E-02 17.12 17 0.00091 291.8 -56.22362 -56.21114 -56.21024 0.114E-02 16.97 18 0.00091 291.5 -56.22360 -56.21114 -56.21024 0.128E-02 16.92 19 0.00091 290.9 -56.22358 -56.21114 -56.21024 0.142E-02 16.95 20 0.00090 290.2 -56.22355 -56.21114 -56.21024 0.156E-02 16.93 21 0.00090 289.4 -56.22351 -56.21114 -56.21024 0.172E-02 16.97 22 0.00089 288.5 -56.22346 -56.21113 -56.21024 0.188E-02 16.98 23 0.00089 287.5 -56.22341 -56.21113 -56.21024 0.204E-02 17.03 24 0.00088 286.4 -56.22336 -56.21112 -56.21024 0.221E-02 16.94 25 0.00088 285.2 -56.22331 -56.21111 -56.21024 0.238E-02 17.15 26 0.00087 284.1 -56.22325 -56.21111 -56.21024 0.256E-02 17.20 **************************************************************** * * * CPMD RECEIVED THE SOFT EXIT REQUEST * * THE COMMAND WAS ISSUED AT Wed Jan 5 18:47:49 2005 * * * **************************************************************** RESTART INFORMATION WRITTEN ON FILE ./RESTART.1 **************************************************************** * AVERAGED QUANTITIES * **************************************************************** MEAN VALUE +/- RMS DEVIATION [-^2]**(1/2) ELECTRON KINETIC ENERGY 0.878232E-03 0.306799E-04 IONIC TEMPERATURE 286.70 4.96 DENSITY FUNCTIONAL ENERGY -56.223370 0.242782E-03 CLASSICAL ENERGY -56.211114 0.308805E-04 CONSERVED ENERGY -56.210236 0.134870E-05 NOSE ENERGY ELECTRONS 0.000000 0.00000 NOSE ENERGY IONS 0.000000 0.00000 CONSTRAINTS ENERGY 0.000000 0.00000 ION DISPLACEMENT 0.928532E-03 0.798684E-03 CPU TIME 17.0027 **************************************************************** * * * FINAL RESULTS * * * **************************************************************** **************************************************************** * ATOMIC COORDINATES * **************************************************************** 1 N 10.854576 13.089703 13.742106 2 C 13.382372 12.197930 14.690812 3 C 15.237143 14.427513 14.413847 4 O 17.161449 14.836654 15.675553 5 O 14.579599 16.249195 12.556731 6 H 12.771366 15.721666 12.090274 7 H 9.594185 13.342572 15.124766 8 H 10.167153 11.667051 12.612667 9 H 13.952616 10.510369 13.513140 10 H 13.353917 11.728770 16.657480 **************************************************************** **************************************************************** ELECTRONIC GRADIENT: MAX. COMPONENT = 8.81506E-04 NORM = 5.33633E-05 TOTAL INTEGRATED ELECTRONIC DENSITY IN G-SPACE = 30.000000 IN R-SPACE = 30.000000 VANDERBILT AUGMENTATION CHARGES (MEAN VALUE PER ATOM) ATOM TYPE NR. OF ATOMS CHARGE N 1 0.304 C 2 0.104 O 2 0.803 H 5 0.031 (K+E1+L+N+X) TOTAL ENERGY = -56.22325437 A.U. (K) KINETIC ENERGY = 28.17681534 A.U. (E1=A-S+R) ELECTROSTATIC ENERGY = -41.31147167 A.U. (S) ESELF = 44.54855464 A.U. (R) ESR = 2.20838604 A.U. (L) LOCAL PSEUDOPOTENTIAL ENERGY = -41.68923375 A.U. (N) N-L PSEUDOPOTENTIAL ENERGY = 13.19031852 A.U. (X) EXCHANGE-CORRELATION ENERGY = -14.58968281 A.U. GRADIENT CORRECTION ENERGY = -0.73387988 A.U. **************************************************************** ================================================================ BIG MEMORY ALLOCATIONS SCR 7012699 EIGRB 5373540 YF 2590058 PSI 2590058 VR 2590058 XF 2590058 YLMB 2418093 QRAD 1527840 QRL 1454400 RHOE 1295029 ---------------------------------------------------------------- [PEAK NUMBER 107] PEAK MEMORY 34632640 = 277.1 MBytes ================================================================ **************************************************************** * * * TIMING * * * **************************************************************** SUBROUTINE CALLS CPU TIME ELAPSED TIME NEWD 27 147.43 147.43 RHOV 27 56.97 56.95 S_INVFFT 432 37.03 37.01 GCENER 27 30.33 30.35 FWFFT 163 26.94 26.91 INVFFT 135 21.87 21.88 XCENER 27 20.74 20.75 RNLSM2 27 19.41 28.25 S_FWFFT 216 17.36 17.44 FFT-G/S 1296 15.94 15.79 RNLSM1 54 12.28 18.21 VPSI 27 9.94 9.93 RHOOFR 27 9.81 9.90 QINIT 1 8.80 8.81 NLFORCE 27 7.63 7.63 PHASE 298 5.08 5.15 VOFRHOT 27 3.94 3.92 EICALC 27 3.79 3.83 GRADEN 27 3.74 3.74 ROTATE 163 2.16 2.17 VOFRHOB 27 1.85 1.83 PHFAC 27 1.82 1.84 ---------------------------------------------------------------- TOTAL TIME 464.86 479.74 **************************************************************** CPU TIME : 0 HOURS 7 MINUTES 50.48 SECONDS ELAPSED TIME : 0 HOURS 8 MINUTES 7.33 SECONDS PROGRAM CPMD ENDED AT: Wed Jan 5 18:47:51 2005 **************************************************************** * * * CPMD RECEIVED THE SOFT EXIT REQUEST * * THE COMMAND WAS ISSUED AT Wed Jan 5 18:47:51 2005 * * * ****************************************************************