Dear Eduardo, ecut_ws is a real space cutoff on the wigner Seitz cell for coulomb. We (Giovanni Bussi Amdrea Ferretti and myself) have developped this for SaX (GW calculations and BSE). The issue was that for strongly anysotropic cells the coulomb potential is not behaving any more as the FFT of 1/q**2 cloes to q=0. Usual techniques like Baldereschi-massidda and others "correct " only q=0, and this is not enough. So the ecut_ws is a cutoff that says to the code since how big q+G it has to correct coulomb: the more anysotropic the more you have to have big this ecut_ws, but it will take more time. The method is similar to Martyna-Tuckermans one withe the difference that we just cutoff the fock operator (this is alsways valid as the overlaps goe to 0 when r-r' -> infinity). If you are interested I have a presentation, and we can discuss these questions by skype. Unfortunatelly, we did nerver write the corresponding paper (always to much work to do). But it is planned to do it one day.  <br>
<br><br>hope this can help<br><br>cheers<br><br>Layla<br><div class="gmail_quote">2012/10/29 Eduardo Ariel Menendez Proupin <span dir="ltr"><<a href="mailto:eariel99@gmail.com" target="_blank">eariel99@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">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<span class="HOEnZb"><font color="#888888"><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: <a href="tel:%2B34%2091%20497%206706" value="+34914976706" target="_blank">+34 91 497 6706</a></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>
</font></span><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br>