[求助] 【已经解决，编译出的问题】DAY3_water_inp1计算出错

&CPMD
  OPTIMIZE WAVEFUNCTION
  STRUCTURE BONDS ANGLES
  RHOOUT
&END

&DFT
  FUNCTIONAL PBE
&END

&SYSTEM
  SYMMETRY
  1
  CELL
  15.0000  1.0  1.0  0  0  0
  CUTOFF
  70.0
&END

&ATOMS
*O_MT_PBE.psp  KLEINMAN-BYLANDER
   LMAX=P
   1
     0.0000 -0.1267  0.0000
*H_MT_PBE.psp
   LMAX=S
   2
     0.0000  1.0058  1.4524
     0.0000  1.0058 -1.4524
&END

PROGRAM CPMD STARTED AT: Sat Aug 31 12:00:25 2013  
SETCNST| USING: CODATA 2006 UNITS


               ******  ******    ****  ****  ******   
              *******  *******   **********  *******  
             ***       **   ***  ** **** **  **   ***
             **        **   ***  **  **  **  **    **
             **        *******   **      **  **    **
             ***       ******    **      **  **   ***
              *******  **        **      **  *******  
               ******  **        **      **  ******   

                          VERSION 3.13.2

                            COPYRIGHT
                      IBM RESEARCH DIVISION
                MPI FESTKOERPERFORSCHUNG STUTTGART

                       The CPMD consortium
                  Home Page: http://www.cpmd.org
               Mailing List: cpmd-list@cpmd.org
                     E-mail: cpmd@cpmd.org


                  ***  Jun  5 2013 -- 09:59:13  ***

THE INPUT FILE IS:                                          inp1
THIS JOB RUNS ON:                                         gg0110
THE CURRENT DIRECTORY IS:
                  /home-gg/users/nscc410_CWL/cpmd_work/DAY3/water
THE TEMPORARY DIRECTORY IS:
                  /home-gg/users/nscc410_CWL/cpmd_work/DAY3/water
THE PROCESS ID IS:                                         12657
THE JOB WAS SUBMITTED BY:                             nscc410_CW


SINGLE POINT DENSITY OPTIMIZATION

PATH TO THE RESTART FILES:                                    ./
GRAM-SCHMIDT ORTHOGONALIZATION
MAXIMUM NUMBER OF STEPS:                             10000 STEPS
MAXIMUM NUMBER OF ITERATIONS FOR SC:                 10000 STEPS
PRINT INTERMEDIATE RESULTS EVERY                     10001 STEPS
STORE INTERMEDIATE RESULTS EVERY                     10001 STEPS
NUMBER OF DISTINCT RESTART FILES:                              1
TEMPERATURE IS CALCULATED ASSUMING EXTENDED BULK BEHAVIOR
STORE ELECTRON DENSITY AT THE END OF THE RUN
FICTITIOUS ELECTRON MASS:                               400.0000
TIME STEP FOR ELECTRONS:                                  5.0000
TIME STEP FOR IONS:                                       5.0000
CONVERGENCE CRITERIA FOR WAVEFUNCTION OPTIMIZATION:   1.0000E-05
WAVEFUNCTION OPTIMIZATION BY PRECONDITIONED DIIS
THRESHOLD FOR THE WF-HESSIAN IS                           0.5000
MAXIMUM NUMBER OF VECTORS RETAINED FOR DIIS:                  10
STEPS UNTIL DIIS RESET ON POOR PROGRESS:                      10
FULL ELECTRONIC GRADIENT IS USED
SPLINE INTERPOLATION IN G-SPACE FOR PSEUDOPOTENTIAL FUNCTIONS
    NUMBER OF SPLINE POINTS:                                 5000

EXCHANGE CORRELATION FUNCTIONALS
    LDA EXCHANGE:                                            NONE
    LDA XC THROUGH PADE APPROXIMATION
    S.GOEDECKER, J.HUTTER, M.TETER PRB 54 1703 (1996)
    GRADIENT CORRECTED FUNCTIONAL
    DENSITY THRESHOLD:                                1.00000E-08
    EXCHANGE ENERGY
       [PBE: J.P. PERDEW ET AL. PRL 77, 3865 (1996)]
    CORRELATION ENERGY
       [PBE: J.P. PERDEW ET AL. PRL 77, 3865 (1996)]

***     DETSP| SIZE OF THE PROGRAM IS   12588/ 319884 kBYTES ***

***************************** ATOMS ****************************
   NR   TYPE        X(bohr)        Y(bohr)        Z(bohr)     MBL
    1      O       0.000000      -0.126700       0.000000       3
    2      H       0.000000       1.005800       1.452400       3
    3      H       0.000000       1.005800      -1.452400       3
****************************************************************

NUMBER OF STATES:                                              4
NUMBER OF ELECTRONS:                                     8.00000
CHARGE:                                                  0.00000
ELECTRON TEMPERATURE(KELVIN):                            0.00000
OCCUPATION
  2.0  2.0  2.0  2.0

    ============================================================  
    |    Pseudopotential Report                                |  
    ------------------------------------------------------------  
    |  Atomic Symbol                   :   O                   |  
    |  Atomic Number                   :   8                   |  
    |  Number of core states           :   1                   |  
    |  Number of valence states        :   2                   |  
    |  Exchange-Correlation Functional :                       |  
    |     Slater exchange :  0.6667                            |  
    |     LDA correlation : Ceperley-Alder                     |  
    |     Exchange GC     : Perdew-Burke-Ernzerhof             |  
    |     Correlation GC  : Perdew-Burke-Ernzerhof             |  
    |  Electron Configuration :   N   L  Occupation            |  
    |                             1   S    2.0000              |  
    |                             2   S    2.0000              |  
    |                             2   P    4.0000              |  
    |  Full Potential Total Energy =      -74.943856           |  
    |  Trouiller-Martins normconserving PP                     |  
    |     n    l        rc       energy                        |  
    |     2    S    1.1200     -0.87882                        |  
    |     2    P    1.1200     -0.33215                        |  
    |     3    D    1.1200     -0.33215                        |  
    |  Number of Mesh Points :   631                           |  
    |  Pseudoatom Total Energy =  -15.770807                   |  
    ============================================================  

    ============================================================  
    |    Pseudopotential Report     Fri Feb 20 19:18:54 1998   |  
    ------------------------------------------------------------  
    |  Atomic Symbol                   :   H                   |  
    |  Atomic Number                   :   1                   |  
    |  Number of core states           :   0                   |  
    |  Number of valence states        :   1                   |  
    |  Exchange-Correlation Functional :                       |  
    |     Slater exchange :   .6667                            |  
    |     LDA correlation : Ceperley-Alder                     |  
    |     Exchange GC     : Perdew-Burke-Ernzerhof             |  
    |     Correlation GC  : Perdew-Burke-Ernzerhof             |  
    |  Electron Configuration :   N   L  Occupation            |  
    |                             1   S     .7000              |  
    |  Full Potential Total Energy     -.474341                |  
    |  Trouiller-Martins normconserving PP                     |  
    |     n    l        rc       energy                        |  
    |     1    S     .5000      -.37203                        |  
    |     2    P     .3828      -.37203                        |  
    |  Number of Mesh Points :   511                           |  
    |  Pseudoatom Total Energy    -.367794                     |  
    ============================================================  

****************************************************************
*   ATOM       MASS   RAGGIO NLCC              PSEUDOPOTENTIAL *
*      O    15.9994   1.2000  NO     KLEINMAN      S  NONLOCAL *
*                                                  P     LOCAL *
*      H     1.0080   1.2000  NO                   S     LOCAL *
****************************************************************


PARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARA
  NCPU     NGW     NHG  PLANES  GXRAYS  HXRAYS ORBITALS Z-PLANES
     0   16686  133590      80    1245    5005       4       1
                G=0 COMPONENT ON PROCESSOR :     0
PARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARAPARA

***    LOADPA| SIZE OF THE PROGRAM IS   19488/ 327236 kBYTES ***

OPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPEN
NUMBER OF CPUS PER TASK                                        1
OPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPENMPOPEN

***     RGGEN| SIZE OF THE PROGRAM IS   22104/ 329844 kBYTES ***

************************** SUPERCELL ***************************
SYMMETRY:                                           SIMPLE CUBIC
LATTICE CONSTANT(a.u.):                                 15.00000
CELL DIMENSION:  15.0000  1.0000  1.0000  0.0000  0.0000  0.0000
VOLUME(OMEGA IN BOHR^3):                              3375.00000
LATTICE VECTOR A1(BOHR):           15.0000     0.0000     0.0000
LATTICE VECTOR A2(BOHR):            0.0000    15.0000     0.0000
LATTICE VECTOR A3(BOHR):            0.0000     0.0000    15.0000
RECIP. LAT. VEC. B1(2Pi/BOHR):      0.0667     0.0000     0.0000
RECIP. LAT. VEC. B2(2Pi/BOHR):      0.0000     0.0667     0.0000
RECIP. LAT. VEC. B3(2Pi/BOHR):      0.0000     0.0000     0.0667
REAL SPACE MESH:                    80           80           80
WAVEFUNCTION CUTOFF(RYDBERG):                           70.00000
DENSITY CUTOFF(RYDBERG):          (DUAL= 4.00)         280.00000
NUMBER OF PLANE WAVES FOR WAVEFUNCTION CUTOFF:             16686
NUMBER OF PLANE WAVES FOR DENSITY CUTOFF:                 133590
****************************************************************

***  RINFORCE| SIZE OF THE PROGRAM IS   28092/ 336616 kBYTES ***
***    FFTPRP| SIZE OF THE PROGRAM IS   44724/ 353224 kBYTES ***

GENERATE ATOMIC BASIS SET
      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:                                0.38 SECONDS

***    WFOPTS| SIZE OF THE PROGRAM IS   48488/ 375928 kBYTES ***
***     PHFAC| SIZE OF THE PROGRAM IS   48644/ 412016 kBYTES ***
***    ATOMWF| SIZE OF THE PROGRAM IS   50476/ 414624 kBYTES ***


PROGRAM STOPS IN SUBROUTINE SCRPTR| RHOG |SCRATCH SPACE NOT FOUND                                                                        [PROC=   0]