 vasp.5.2.11 18Jan11 complex
 executed on             LinuxIFC date 2015.09.01  13:14:57
 running on   16 nodes
 distr:  one band on    1 nodes,   16 groups


--------------------------------------------------------------------------------------------------------


 INCAR:
 POTCAR:   PAW_GGA C 05Jan2001
 POTCAR:   PAW_GGA H 07Jul1998
 POTCAR:   PAW_GGA C 05Jan2001
   VRHFIN =C: s2p2
   LEXCH  = 91
   EATOM  =   147.4688 eV,   10.8386 Ry

   TITEL  = PAW_GGA C 05Jan2001
   LULTRA =        F    use ultrasoft PP ?
   IUNSCR =        0    unscreen: 0-lin 1-nonlin 2-no
   RPACOR =     .000    partial core radius
   POMASS =   12.011; ZVAL   =    4.000    mass and valenz
   RCORE  =    1.500    outmost cutoff radius
   RWIGS  =    1.630; RWIGS  =     .863    wigner-seitz radius (au A)
   ENMAX  =  400.000; ENMIN  =  300.000 eV
   ICORE  =        2    local potential
   LCOR   =        T    correct aug charges
   LPAW   =        T    paw PP
   EAUG   =  644.873
   DEXC   =     .000
   RMAX   =    2.266    core radius for proj-oper
   RAUG   =    1.300    factor for augmentation sphere
   RDEP   =    1.501    radius for radial grids
   RDEPT  =    1.300    core radius for aug-charge
   QCUT   =   -5.516; QGAM   =   11.033    optimization parameters

   Description
     l     E      TYP  RCUT    TYP  RCUT
     0   .000     23  1.200
     0   .000     23  1.200
     1   .000     23  1.500
     1  2.500     23  1.500
     2   .000      7  1.500
  local pseudopotential read in
  atomic valenz-charges read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
    PAW grid and wavefunctions read in

   number of l-projection  operators is LMAX  =           4
   number of lm-projection operators is LMMAX =           8

 POTCAR:   PAW_GGA H 07Jul1998
   VRHFIN =H: ultrasoft test
   LEXCH  = 91
   EATOM  =    12.5313 eV,     .9210 Ry

   TITEL  = PAW_GGA H 07Jul1998
   LULTRA =        F    use ultrasoft PP ?
   IUNSCR =        0    unscreen: 0-lin 1-nonlin 2-no
   RPACOR =     .000    partial core radius
   POMASS =    1.000; ZVAL   =    1.000    mass and valenz
   RCORE  =    1.100    outmost cutoff radius
   RWIGS  =     .700; RWIGS  =     .370    wigner-seitz radius (au A)
   ENMAX  =  250.000; ENMIN  =  200.000 eV
   RCLOC  =     .701    cutoff for local pot
   LCOR   =        T    correct aug charges
   LPAW   =        T    paw PP
   EAUG   =  400.000
   RMAX   =    2.174    core radius for proj-oper
   RAUG   =    1.200    factor for augmentation sphere
   RDEP   =    1.112    core radius for depl-charge
   QCUT   =   -5.749; QGAM   =   11.498    optimization parameters

   Description
     l     E      TYP  RCUT    TYP  RCUT
     0   .000     23  1.100
     0   .500     23  1.100
     1  -.300     23  1.100
  local pseudopotential read in
  atomic valenz-charges read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           0  read in
    real space projection operators read in
  non local Contribution for L=           1  read in
    real space projection operators read in
    PAW grid and wavefunctions read in

   number of l-projection  operators is LMAX  =           3
   number of lm-projection operators is LMMAX =           5

 PAW_GGA C 05Jan2001                    :
 energy of atom  1       EATOM= -147.4688
 kinetic energy error for atom=    0.0288 (will be added to EATOM!!)
 PAW_GGA H 07Jul1998                    :
 energy of atom  2       EATOM=  -12.5313
 kinetic energy error for atom=    0.0098 (will be added to EATOM!!)


 POSCAR:  C H
  positions in direct lattice
  velocities in cartesian coordinates
 exchange correlation table for  LEXCH =        7
   RHO(1)=    0.500       N(1)  =     2000
   RHO(2)=  100.500       N(2)  =     4000

 VTST: version 2.03d, (02/18/09)

 CHAIN: initializing optimizer

 OPT: Using VASP Dynamics algorithm
 CHAIN: Read ICHAIN            0



--------------------------------------------------------------------------------------------------------


 ion  position               nearest neighbor table
   1  0.477  0.509  0.535-   2 1.10   3 1.10   4 1.10   5 1.10
   2  0.448  0.569  0.564-   1 1.10
   3  0.447  0.451  0.567-   1 1.10
   4  0.549  0.511  0.547-   1 1.10
   5  0.465  0.506  0.463-   1 1.10

  LATTYP: Found a simple cubic cell.
 ALAT       =    15.0000000000

  Lattice vectors:

 A1 = (  15.0000000000,   0.0000000000,   0.0000000000)
 A2 = (   0.0000000000,  15.0000000000,   0.0000000000)
 A3 = (   0.0000000000,   0.0000000000,  15.0000000000)
 Subroutine PRICEL returns:
 Original cell was already a primitive cell.


Analysis of symmetry for initial positions (statically):

 Routine SETGRP: Setting up the symmetry group for a
 simple cubic supercell.


 Subroutine GETGRP returns: Found  1 space group operations
 (whereof  1 operations were pure point group operations)
 out of a pool of 48 trial point group operations.


The static configuration has the point symmetry C_1 .

Analysis of symmetry for dynamics (positions and initial velocities):

 Subroutine DYNSYM returns: Found  1 space group operations
 (whereof  1 operations were pure point group operations)
 out of a pool of  1 trial space group operations
 (whereof  1 operations were pure point group operations)
 and found also     1 'primitive' translations


The dynamic configuration has the point symmetry C_1 .

Analysis of constrained symmetry for selective dynamics:

 Subroutine DYNSYM returns: Found  1 space group operations
 (whereof  1 operations were pure point group operations)
 out of a pool of  1 trial space group operations
 (whereof  1 operations were pure point group operations)
 and found also     1 'primitive' translations


The constrained configuration has the point symmetry C_1 .


 KPOINTS: Automatic Mesh

Automatic generation of k-mesh.

 Subroutine IBZKPT returns following result:
 ===========================================

 Found      1 irreducible k-points:

 Following reciprocal coordinates:
            Coordinates               Weight
  0.000000  0.000000  0.000000      1.000000

 Following cartesian coordinates:
            Coordinates               Weight
  0.000000  0.000000  0.000000      1.000000


 Subroutine IBZKPT_HF returns following result:
 ==============================================

 Found      1 k-points in 1st BZ
 the following      1 k-points will be used (e.g. in the exchange kernel)
 Following reciprocal coordinates:   # in IRBZ
  0.000000  0.000000  0.000000       1


--------------------------------------------------------------------------------------------------------




 Dimension of arrays:
   k-points           NKPTS =      1   k-points in BZ     NKDIM =      1   number of bands    NBANDS=     16
   number of dos      NEDOS =    301   number of ions     NIONS =      5
   non local maximal  LDIM  =      4   non local SUM 2l+1 LMDIM =      8
   total plane-waves  NPLWV = 512000
   max r-space proj   IRMAX =      1   max aug-charges    IRDMAX=   1990
   dimension x,y,z NGX =    80 NGY =   80 NGZ =   80
   dimension x,y,z NGXF=   160 NGYF=  160 NGZF=  160
   support grid    NGXF=   160 NGYF=  160 NGZF=  160
   ions per type =               1   4
 NGX,Y,Z   is equivalent  to a cutoff of   8.87,  8.87,  8.87 a.u.
 NGXF,Y,Z  is equivalent  to a cutoff of  17.73, 17.73, 17.73 a.u.


 I would recommend the setting:
   dimension x,y,z NGX =    73 NGY =   73 NGZ =   73
 SYSTEM =  Pd111
 POSCAR =   C H

 Startparameter for this run:
   NWRITE =      2    write-flag & timer
   PREC   = normal    normal or accurate (medium, high low for compatibility)
   ISTART =      0    job   : 0-new  1-cont  2-samecut
   ICHARG =      2    charge: 1-file 2-atom 10-const
   ISPIN  =      1    spin polarized calculation?
   LNONCOLLINEAR =      F non collinear calculations
   LSORBIT =      F    spin-orbit coupling
   INIWAV =      1    electr: 0-lowe 1-rand  2-diag
   LASPH  =      F    aspherical Exc in radial PAW
   METAGGA=      F    non-selfconsistent MetaGGA calc.

 Electronic Relaxation 1
   ENCUT  =  400.0 eV  29.40 Ry    5.42 a.u.  24.46 24.46 24.46*2*pi/ulx,y,z
   ENINI  =  400.0     initial cutoff
   ENAUG  =  644.9 eV  augmentation charge cutoff
   NELM   =     60;   NELMIN=  2; NELMDL=  0     # of ELM steps
   EDIFF  = 0.1E-05   stopping-criterion for ELM
   LREAL  =      F    real-space projection
   LCOMPAT=      F    compatible to vasp.4.4
   GGA_COMPAT  = T    GGA compatible to vasp.4.4-vasp.4.6
   LMAXPAW     = -100 max onsite density
   LMAXMIX     =    2 max onsite mixed and CHGCAR
   VOSKOWN=      1    Vosko Wilk Nusair interpolation
   ROPT   =    0.00000   0.00000
 Ionic relaxation
   EDIFFG = -.1E-01   stopping-criterion for IOM
   NSW    =     97    number of steps for IOM
   NBLOCK =      1;   KBLOCK =      1    inner block; outer block
   IBRION =      5    ionic relax: 0-MD 1-quasi-New 2-CG
   NFREE  =      2    steps in history (QN), initial steepest desc. (CG)
   ISIF   =      0    stress and relaxation
   IWAVPR =     11    prediction:  0-non 1-charg 2-wave 3-comb
   ISYM   =      2    0-nonsym 1-usesym 2-fastsym
   LCORR  =      T    Harris-Foulkes like correction to forces

   POTIM  = 0.0150    time-step for ionic-motion
   TEIN   =    0.0    initial temperature
   TEBEG  =    0.0;   TEEND  =   0.0 temperature during run
   SMASS  =  -3.00    Nose mass-parameter (am)
   estimated Nose-frequenzy (Omega)   =  0.10E-29 period in steps =****** mass=  -0.514E-26a.u.
   NPACO  =    256;   APACO  = 16.0  distance and # of slots for P.C.
   PSTRESS=    0.0 pullay stress

  Mass of Ions in am
   POMASS =  12.01  1.00
  Ionic Valenz
   ZVAL   =   4.00  1.00
  Atomic Wigner-Seitz radii
   RWIGS  =  -1.00 -1.00
  virtual crystal weights
   VCA    =   1.00  1.00
   NELECT =       8.0000    total number of electrons
   NUPDOWN=      -1.0000    fix difference up-down

 DOS related values:
   EMIN   =  10.00;   EMAX   =-10.00  energy-range for DOS
   EFERMI =   0.00
   ISMEAR =     1;   SIGMA  =   0.20  broadening in eV -4-tet -1-fermi 0-gaus

 Electronic relaxation 2 (details)
   IALGO  =     48    algorithm
   LDIAG  =      T    sub-space diagonalisation (order eigenvalues)
   LSUBROT=      T    optimize rotation matrix (better conditioning)
   TURBO    =      0    0=normal 1=particle mesh
   IRESTART =      0    0=no restart 2=restart with 2 vectors
   NREBOOT  =      0    no. of reboots
   NMIN     =      0    reboot dimension
   EREF     =   0.00    reference energy to select bands
   IMIX   =      4    mixing-type and parameters
     AMIX     =   0.40;   BMIX     =  1.00
     AMIX_MAG =   1.60;   BMIX_MAG =  1.00
     AMIN     =   0.10
     WC   =   100.;   INIMIX=   1;  MIXPRE=   1

 Intra band minimization:
   WEIMIN = 0.0010     energy-eigenvalue tresh-hold
   EBREAK =  0.31E-07  absolut break condition
   DEPER  =   0.30     relativ break condition

   TIME   =   0.40     timestep for ELM

  volume/ion in A,a.u.               =     675.00      4555.12
  Fermi-wavevector in a.u.,A,eV,Ry     =   0.218280  0.412489  0.648264  0.047646
  Thomas-Fermi vector in A             =   0.996232

 Write flags
   LWAVE  =      T    write WAVECAR
   LCHARG =      F    write CHGCAR
   LVTOT  =      F    write LOCPOT, total local potential
   LVHAR  =      F    write LOCPOT, Hartree potential only
   LELF   =      F    write electronic localiz. function (ELF)
   LORBIT =      0    0 simple, 1 ext, 2 COOP (PROOUT)


 Dipole corrections
   LMONO  =      F    monopole corrections only (constant potential shift)
   LDIPOL =      F    correct potential (dipole corrections)
   IDIPOL =      0    1-x, 2-y, 3-z, 4-all directions
   EPSILON=  1.0000000 bulk dielectric constant

 Exchange correlation treatment:
   GGA     =    91    GGA type
   LEXCH   =     7    internal setting for exchange type
   VOSKOWN=      1    Vosko Wilk Nusair interpolation
   LHFCALC =     F    Hartree Fock is set to
   LHFONE  =     F    Hartree Fock one center treatment
   AEXX    =    0.0000 exact exchange contribution

 Linear response parameters
   LEPSILON=     F    determine dielectric tensor
   LRPA    =     F    only Hartree local field effects (RPA)
   LNABLA  =     F    use nabla operator in PAW spheres
   LVEL    =     F    velocity operator in full k-point grid
   LINTERFAST=   F  fast interpolation
   KINTER  =     0    interpolate to denser k-point grid
   CSHIFT  =0.1000    complex shift for real part using Kramers Kronig
   OMEGAMAX=  -1.0    maximum frequency
   RTIME   =    0.100 relaxation time in fs

 Orbital magnetization related:
   ORBITALMAG=     F  switch on orbital magnetization
   LCHIMAG   =     F  perturbation theory with respect to B field



############################ A LOT OF LINES WERE REMOVED HERE ############################


 POSITION                                       TOTAL-FORCE (eV/Angst)
 -----------------------------------------------------------------------------------
      7.15782      7.63752      8.02990         0.005432     -0.010490     -0.003315
      6.71400      8.53969      8.46467        -0.000780     -0.002283      0.004837
      6.69848      6.76009      8.49860        -0.002744      0.008204     -0.000537
      8.23960      7.66298      8.20237         0.000601      0.005181      0.002440
      6.98190      7.58555      6.94952        -0.002509     -0.000612     -0.003425
 -----------------------------------------------------------------------------------
    total drift:                                0.000896     -0.000486     -0.000729


--------------------------------------------------------------------------------------------------------



  FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV)
  ---------------------------------------------------
  free  energy   TOTEN  =       -24.098768 eV

  energy  without entropy=      -24.123642  energy(sigma->0) =      -24.107059


############################ A LOT OF LINES WERE REMOVED HERE ############################


 POSITION                                       TOTAL-FORCE (eV/Angst)
 -----------------------------------------------------------------------------------
      7.15782      7.63752      8.02990        -0.126084      0.088256     -0.400029
      6.71400      8.53969      8.46467         0.019886     -0.033257     -0.026713
      6.69848      6.76009      8.49860         0.014606     -0.029373     -0.020520
      8.23960      7.66298      8.20237         0.023705     -0.023411      0.002259
      6.98190      7.58555      6.93452         0.067888     -0.002216      0.445002
 -----------------------------------------------------------------------------------
    total drift:                                0.000869     -0.000546     -0.000813


--------------------------------------------------------------------------------------------------------



  FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV)
  ---------------------------------------------------
  free  energy   TOTEN  =       -24.095383 eV

  energy  without entropy=      -24.119696  energy(sigma->0) =      -24.103487

 d Force = 0.8542991E-03[-0.134E-01, 0.151E-01]  d Energy = 0.4182863E-03 0.436E-03
 d Force = 0.1740500E+01[ 0.170E+01, 0.178E+01]  d Ewald  = 0.1739999E+01 0.501E-03


--------------------------------------------------------------------------------------------------------


    POTLOK:  cpu time    0.37: real time    0.38


--------------------------------------------------------------------------------------------------------


  FORCES: max atom, RMS     0.450156    0.202610
  FORCE total and by dimension    0.453049    0.445002
 Finite differences progress:
  Degree of freedom:   9/  9
  Displacement:        2/  2
  Total:              18/ 18

 SECOND DERIVATIVES (NOT SYMMETRIZED)
 ------------------------------------
               3X          3Y          3Z          4X          4Y          4Z          5X          5Y          5Z
  3X   -19.248262   -4.329092   -0.060057    2.836454   13.726885   -8.959895    7.398112  -10.773197   -1.116170
  3Y    -3.332995  -25.578420   11.838015    2.361044   -4.411883    5.537489   -1.875964    2.571285    2.954898
  3Z    -0.727966   11.815416  -13.605773   -1.594963    8.082568   -1.653277   -1.452319   -7.332186    0.917624
  4X     2.600373    1.876728   -1.187104  -31.624347    0.198924   -4.608783   -0.622354    2.693614   -1.829436
  4Y    12.957773   -3.262637    6.449113    0.130296  -14.320054    6.103052   -3.997455   10.526720    1.606990
  4Z    -9.520088    5.226909   -1.546106   -4.837118    6.890380  -11.022359    8.088659   -4.817398   -0.404036
  5X     7.517545   -1.769121   -1.354656   -0.357433   -4.075842    8.204166  -10.617809    6.726545   -4.921521
  5Y    -9.994642    2.859550   -6.636530    1.770754   10.585980   -4.792977    5.840417  -14.881737   -0.555992
  5Z    -1.531617    2.896356    0.926072   -1.800645    2.287990   -0.541162   -4.825002   -0.166897  -31.565272


############################ A LOT OF LINES WERE REMOVED HERE ############################
