<p>Hi Dear user,</p>
<p>I have tried to test the k-points and ecutwfc convergence for<br>a Mg(002) slab which contains 7 layers and the middle 5 layers<br>have been fixed.</p>
<p>PROBLEM I:<br>I have tested (1x1x1, 2x2x1, 4x4x1, 5x5x1, 7x7x1, 8x8x1, 10x10x1) M-P<br>k-points grids for one input file, other parameters are the same(see below).</p>
<p>For 2x2x1, 5x5x1 and 8x8x1, I have got<br>"The maximum number of steps has been reached.<br> End of BFGS Geometry Optimization" <br>However, for others, say, (1x1x1, 4x4x1, 7x7x1, 10x10x1)I have got <br>"bfgs converged in xx scf cycles and xx bfgs steps
<br> End of BFGS Geometry Optimization"</p>
<p>For the latter cases, the convergences have completed successfully, <br>but for the former cases, they seemed not finished. I wonder, if just<br>for k-points test, whether it is worth to restart the jobs to make the<br>
optimization converged? Or, choosing new K sampling, such as (3x3x1, <br>6x6x1, 9x9x1, etc.)?</p>
<p>PROBLEM II<br>when testing ecutwfc convergence, I have chosen 10, 14, 18, 22, 26, 30 Ry<br>as the value of the ecutwfc. It works well for ecutwfc=18, 26 and 30. When<br>ecutwfc=14 or 22, the following errors occured:<br>
...<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>     from bfgs : error #         1<br>     bfgs history already reset at previous step<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</p>
<p>     stopping ...</p>
<p>I have checked the mailing list archives, like:<br><a href="http://www.democritos.it/pipermail/pw_forum/2006-November/005331.html">http://www.democritos.it/pipermail/pw_forum/2006-November/005331.html</a><br>suggested that restarting the jobs to make the geometry optimization converged.
<br>However, I wonder why the optimizaition is failed for certain lower values of <br>ecutwfc which seems would be more easily to be completed. </p>
<p>In PROBLEM I, ecutwfc=24Ry; In PROBLEM II, K_POINTS are 5 5 1 0 0 0</p>
<p>the following is the input file:<br>&CONTROL<br>calculation  = 'relax',<br>prefix       = 'Mgrelax',<br>pseudo_dir   = '/....../pseudo/',<br>outdir       = '/....../tmp/',<br>tstress=.true.
<br>tprnfor=.true.<br>/<br>&SYSTEM<br>ibrav       = 14,<br>a           = 3.149307D0,<br>b           = 3.149307D0,<br>c           = 32.624496,<br>cosab       = 0.5D0,<br>cosac       = 0.D0,<br>cosbc       = 0.D0,<br>nat         = 7,
<br>ntyp        = 1,<br>ecutwfc     = xx.D0,<br>occupations = "smearing",<br>smearing    = "methfessel-paxton",<br>degauss     = 0.05D0,<br>/<br>&ELECTRONS<br>/<br>&IONS<br>ion_dynamics='bfgs',
<br>pot_extrapolation = "second_order",<br>wfc_extrapolation = "second_order",<br>/<br>ATOMIC_SPECIES<br>Mg  24.305  Mg.pw91-np-van.UPF<br>ATOMIC_POSITIONS {angstrom}<br>Mg   2.099538   2.099538   0.000000
<br>Mg   1.049769   1.049769   2.718708   0 0 0<br>Mg   2.099538   2.099538   5.437416   0 0 0<br>Mg   1.049769   1.049769   8.156124   0 0 0<br>Mg   2.099538   2.099538  10.874832   0 0 0<br>Mg   1.049769   1.049769  13.59354
    0 0 0<br>Mg   2.099538   2.099538  16.312248<br>K_POINTS {automatic}<br>nk1 nk2 nk3 0 0 0<br></p>