Hello,<br>         I assumed surface relaxation for the Ag(110) surface would be an easy matter. As it happens, I find non-zero forces parallel to the surface and therefore, upon relaxation, I get a non-neglectable parallel displacement of the ions. I am working with a 2x2sqrt(2) unit cell. I attach my <a href="http://file.in">file.in</a> and the forces I get. Note that the forces along one direction (x) are much larger that along the other one (y). Note as well the large values considered for ecutwfc, ecutrho, kpoints and the small degauss and conv_thr. Any ideas on how to correct this behaviour? Am I missing something? Thanks in advance and sorry if this is a stupid question. Small changes in the input parameters can lower the forces I get but I find nothing systematic. <br>
<br>&control<br>   prefix        = 't_G0',<br>   outdir        = ' ',     <br>   pseudo_dir    = '/home/qe/pseudo',    <br>   calculation   = 'scf',<br>   tprnfor       = .true.,<br>   lkpoint_dir   = .false.,<br>
   verbosity     = 'high',<br>/<br>&system<br>   ibrav       = 8,<br>   celldm(1)   = 11.1298607358762567,<br>   celldm(2)   = 1.4142135623730949,<br>   celldm(3)   = 2.25,<br>   nat         = 16,<br>   ntyp        = 1,<br>
   ecutwfc     = 45,<br>   ecutrho     = 450.00,<br>   occupations = 'smearing',   <br>   smearing    = 'mv',         <br>   degauss     = 0.002,      <br>/<br>&electrons<br>   conv_thr = 1.d-12<br>   diagonalization = 'cg'<br>
/<br>&ions<br>/<br>ATOMIC_SPECIES<br>   Ag   107.8682  Ag.pbe-d-rrkjus.UPF  <br>ATOMIC_POSITIONS  crystal  <br>   Ag   0.0000000000000000   0.0000000000000000   0.3333333333333333   <br>   Ag   0.5000000000000000   0.0000000000000000   0.3333333333333333   <br>
   Ag   0.0000000000000000   0.5000000000000000   0.3333333333333333   <br>   Ag   0.5000000000000000   0.5000000000000000   0.3333333333333333   <br>   Ag   0.2500000000000000   0.2500000000000000   0.2222222222222222   <br>
   Ag   0.7500000000000000   0.2500000000000000   0.2222222222222222   <br>   Ag   0.2500000000000000   0.7500000000000000   0.2222222222222222   <br>   Ag   0.7500000000000000   0.7500000000000000   0.2222222222222222   <br>
   Ag   0.0000000000000000   0.0000000000000000   0.1111111111111111   <br>   Ag   0.5000000000000000   0.0000000000000000   0.1111111111111111   <br>   Ag   0.0000000000000000   0.5000000000000000   0.1111111111111111   <br>
   Ag   0.5000000000000000   0.5000000000000000   0.1111111111111111   <br>   Ag   0.2500000000000000   0.2500000000000000   0.0000000000000000   <br>   Ag   0.7500000000000000   0.2500000000000000   0.0000000000000000   <br>
   Ag   0.2500000000000000   0.7500000000000000   0.0000000000000000   <br>   Ag   0.7500000000000000   0.7500000000000000   0.0000000000000000   <br>K_POINTS automatic<br>   33 23 1  0 0 0<br><br>     atom   1 type  1   force =     0.00000000    0.00000000   -0.01175276<br>
     atom   2 type  1   force =     0.00000000    0.00000000   -0.01175282<br>     atom   3 type  1   force =     0.00000000    0.00000000   -0.01175318<br>     atom   4 type  1   force =     0.00000000    0.00000000   -0.01175325<br>
     atom   5 type  1   force =    -0.00004221    0.00000033    0.01048323<br>     atom   6 type  1   force =     0.00004221    0.00000033    0.01048323<br>     atom   7 type  1   force =    -0.00004221   -0.00000033    0.01048323<br>
     atom   8 type  1   force =     0.00004221   -0.00000033    0.01048323<br>     atom   9 type  1   force =     0.00000000    0.00000000   -0.01048309<br>     atom  10 type  1   force =     0.00000000    0.00000000   -0.01048306<br>
     atom  11 type  1   force =     0.00000000    0.00000000   -0.01048272<br>     atom  12 type  1   force =     0.00000000    0.00000000   -0.01048269<br>     atom  13 type  1   force =    -0.00004148    0.00000003    0.01175266<br>
     atom  14 type  1   force =     0.00004148    0.00000003    0.01175266<br>     atom  15 type  1   force =    -0.00004148   -0.00000003    0.01175266<br>     atom  16 type  1   force =     0.00004148   -0.00000003    0.01175266<br>
<br>Dr. G. A. Bocan<br>Surface Physics Group, <br>Centro Atómico Bariloche, <br>