<DIV><STRONG><FONT size="4">Dear all</FONT></STRONG></DIV>
<DIV><STRONG><FONT size="4">     In the forces_bp_efield.f90,one part of code is as following:</FONT></STRONG></DIV>
<DIV>     !  --- Check that we are working with an insulator with no empty bands ---<BR>   IF ((degauss > 0.01d0) .OR. (nbnd /= nelec/2)) &<BR>        WRITE (stdout,*) 'PAY ATTENTION: EL FIELD AND OCCUPATIONS'</DIV>
<DIV><FONT size="4">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 ? </FONT></DIV>
<DIV><FONT size="4">   </FONT><BR>--<BR></DIV>
<DIV>
<DIV>Best regards</DIV>
<DIV> </DIV>
<DIV>Q.J.Wang</DIV>
<DIV> </DIV>
<DIV>XiangTan University </DIV></DIV>