Hi Layla and everybody,<br> <br>Let us forget the nscf calculation for the moment, as my scf calculation was probably not fine. <br>I have a strongly anisotropic unit cell, and large enough not to dare to attempt to build a quasi-cubi one, that would have about 360 atoms. For anysotropic cell,  INPUT_PW instructs to use <br>
<br>                      exxdiv_treatment = 'vcut_ws',<br>for this to work one must use <br>                          x_gamma_extrapolation = .false.,   (b.t.w., it is not mentioned at INPUT_PW.txt )<br>then I have to set the value <br>
                          ecutvcut = X ,  ! which is 0.7 Ry in the Si example. <br> so, what is the method followed by using vcut_ws ,  and how to set ecutvcut ? Are low or high values of ecutwfc more or less accurate or more or less expensive? Is it the method given by Sorouri, Folkes and Hine in JCP124, 064105 (2006)?<br>
<br clear="all">Thanks<br>-- <br><div><br></div>
<div><br></div>Eduardo Menendez Proupin<div>Departamento de Química Fisica Aplicada<br>Facultad de Ciencias<br>Universidad Autónoma de Madrid<br>28049 Madrid, Spain<br></div><div>Phone: +34 91 497 6706</div><div><br></div>
<div>On leave from: Departamento de Fisica, Facultad de Ciencias, Universidad de Chile URL: <a href="http://fisica.ciencias.uchile.cl/%7Eemenendez" target="_blank">http://fisica.ciencias.uchile.cl/~emenendez</a><div><br></div>
</div><br>