[Pw_forum] forces_bp_efield.f90

Sun Dec 6 13:30:07 CET 2009

Dear all
     In the forces_bp_efield.f90,one part of code is as following:
     !  --- Check that we are working with an insulator with no empty bands ---
   IF ((degauss > 0.01d0) .OR. (nbnd /= nelec/2)) &
        WRITE (stdout,*) 'PAY ATTENTION: EL FIELD AND OCCUPATIONS'
So I want to know the forces_bp_efield.f90 only fit to insulator ? If my system is half metallic ,can I use this code to relax ? If my system have odd electrons ,how to set it ? 
   
--

Best regards
 
Q.J.Wang
 
XiangTan University 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20091206/59991a94/attachment.html>