<div dir="ltr"><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Debnath Talukdar</b> <span dir="ltr"><<a href="mailto:debnath.talukdar@students.iiserpune.ac.in">debnath.talukdar@students.iiserpune.ac.in</a>></span><br>Date: Tue, Sep 13, 2016 at 4:25 PM<br>Subject: Error in Hybrid Exc calculations<br>To: PWSCF Forum <<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>><br><br><br><div dir="ltr"><div><div>Hi,<br></div>   I was trying to do hybrid calculations in quantum espresso. The calculation stopped by printing the following message<br><br>    Error in routine electrons (1):<br>     dexx is negative!   Check that exxdiv_treatment is appropriate for the system<br><br></div><div>As I know that in hybrid Exc calculation to solve the poisson equation in reciprocal space their is <br></div><div>fourier transform of density from real to q space. That colomb potential in poisson equation is diverge at q=0. To remove this divergencies in coulomb potential in q space is adds a small G to<br></div><div>denomenetor(q+G); known as   gygi-baldereschi method, but it is appropriate  for cubic and quasi cubic supercell (Gygi-Baldereschi PRB 34, 4405 (1986)).<br> <br>At first I have used  gygi-baldereschi method for my system that is hexagonal metallic but I have found above error then I have changed this method from 'gygi-baldereschi' to 'vcut_ws' that is appropriate for anisotropic supercell but I still got that same error message.<br><br></div><div>After that I have changed Hybrid functional from HSE to gau-PBE which used exxdiv_treatment='none' and x_gamma_extrapolation=.false. but I still got that same error <br> message again.<br><br><br></div><div>Can anybody give any suggestion that how can I solve this problem ? <br><br></div><div>I have also doubt about following...<br><br></div><div>1)As we all know that hybrid calculation is very expansive so how can I reduce different   parameter in Hybrid calculation as compared to normal scf calculations.<br> <br></div><div>   Such that in my normal scf with DFT Exc I have used k -mesh 21x21x1, scf threshold 10-8Ry,<br></div><div>   then how can I reduce these parameters in Hybrid calculations; is it arbitrary? <br><br></div><div>2)How can I reduce ecutfock to speed up the calculation. Currently I have used ecutfock=ecutrho (for my system ecutrho=10*ecutwfc), is their is any way to reduce this kinetic energy cutoff for fock (exchange) operator at the cost of loss in accuracy in my calculations.<br><br></div><div>3) In our pwscf input file we give explicitly effect of exchange-correlation functional through card<br></div><div>input_dft, I am doubt that how can I gave both vdW(rev-vdw-b86r) functional with hybrid calculations(where input_dft='HSE').<br><br></div><div>I have read papers but I am unable to solve these issues can anyone will give any other papers that how can I solve this.<br></div><div><br></div><div>Below I have attached input, output files of my calculations.<br>1)PRB 34, 4405 (1986)<br>2)PHYSICAL REVIEW B 80, 115205 ͑2009͒<br>3)PHYSICAL REVIEW B 74, 121102͑R͒ ͑2006͒<br></div><div>4)JCP 118, 8207 (2003)<br>5)THE JOURNAL OF CHEMICAL PHYSICS 134, 024116 (2011)<br><br></div><div><br><br><br></div><div>Thanks<br></div></div>
</div><br></div>