Dear Forum,<br><br>I calculated the energy of a aluminium surface, but the value that i obtained using this is 'extremely high'  energy, i obtain the same using 4.2.1 and 4.1 version of QE.<br><br>My input file is :<br>

=========================================<br>&CONTROL<br>   calculation           = 'scf'<br>   prefix                = 'al-111-43.5'<br>   tprnfor               = .true.<br>   wfcdir                = './restore/'<br>

   pseudo_dir            = '/home/jperalta/pseudo/'<br>   outdir                = './restore'<br>/<br>&SYSTEM<br>   ibrav                 =  0<br>   celldm(1)             =  1.0<br>   nat                   = 21<br>

   ntyp                  = 1<br>   ecutwfc               = 30<br>   ecutrho               = 500<br>   degauss               = 0.03<br>   occupations           = 'smearing'<br>   smearing              = 'mv'<br>

   nspin                 = 2<br>   starting_magnetization= 0.1<br>   nbnd                  = 300<br>/<br>&ELECTRONS<br>   electron_maxstep      = 300<br>   mixing_mode           = 'local-TF'<br>   mixing_beta           = 0.3<br>

   mixing_ndim           = 10<br>/<br>&IONS<br>!   ion_dynamics          = 'bfgs'<br>/<br>&CELL<br>!   cell_dynamics         = 'bfgs'<br>!   press                 = 0.0D0<br>!   press_conv_thr        = 0.1D0<br>

/<br>ATOMIC_SPECIES<br>Al  26.982 Al.pbe-sp-van.UPF<br>ATOMIC_POSITIONS (angstrom)<br>Al   4.295143   -2.479802   43.5<br>Al   5.726858   -3.306403   41.689899<br>Al   1.431714   -0.826601   41.689899<br>Al   2.863429   -3.306403   41.689899<br>

Al   4.295143   -0.826601   41.689899<br>Al   2.863429   -1.653201   39.351919<br>Al   7.158572   -4.133004   39.351919<br>Al   5.726858   -1.653201   39.351919<br>Al   4.295143   -4.133004   39.351919<br>Al   8.590287   -4.959604   37.013940<br>

Al   4.295143   -2.479802   37.013940<br>Al   5.726858   -4.959604   37.013940<br>Al   7.158572   -2.479802   37.013940<br>Al   5.726858   -3.306403   34.675960<br>Al   1.431714   -0.826601   34.675960<br>Al   2.863429   -3.306403   34.675960<br>

Al   4.295143   -0.826601   34.675960<br>Al   2.863429   -1.653201   32.337980<br>Al   7.158572   -4.133004   32.337980<br>Al   5.726858   -1.653201   32.337980<br>Al   4.295143   -4.133004   32.337980<br>K_POINTS (automatic)<br>

6 6 6 0 0 0<br>CELL_PARAMETERS {cubic}<br>10.822192        0.000000        0.000000<br>5.411096         -9.372293       0.000000<br>0.000000         0.000000        139.892407<br>================================<br><br>i change some electrons parameters in order to improve the convergence. I don't know if the enery problem is caused for some special parameter or some like that. but in a paper of Sanchez et al (Molecular Physics, 2004, vol 102, 1045-1055), the energies are aroung 2 to 20 eV.<br>

<br>And the ouput in my case say :<br><br>    the Fermi energy is    -1.8658 ev<br><br>!    total energy              =   -3357.65873780 Ry  ---------------> 2K eV by atom.<br>     Harris-Foulkes estimate   =   -3357.65873669 Ry<br>

     estimated scf accuracy    <       0.00000057 Ry<br><br>     The total energy is the sum of the following terms:<br><br>     one-electron contribution = -110140.98098627 Ry<br>     hartree contribution      =   54852.05508820 Ry<br>

     xc contribution           =    -409.97773145 Ry<br>     ewald contribution        =   52341.24991020 Ry<br>     smearing contrib. (-TS)   =      -0.00501849 Ry<br><br>     total magnetization       =     0.00 Bohr mag/cell<br>

     absolute magnetization    =     0.00 Bohr mag/cell<br><br>     convergence has been achieved in  34 iterations<br><br><br>Any help i really appreciate it.<br><br>Best Regards<br>Joaquin Peralta<br>Iowa State University.<br clear="all">

<br>-- <br>----------------------------------------------------<br>Group of NanoMaterials<br>----------------------------------------------------<br><a href="http://www.gnm.cl" target="_blank">http://www.gnm.cl</a><br>----------------------------------------------------<br>

Joaquín Andrés Peralta Camposano<br>----------------------------------------------------<br><a href="http://www.lpmd.cl/jperalta" target="_blank">http://www.lpmd.cl/jperalta</a><br><br>In a world without frontiers, <br>who needs Gates and Win.<br>