<DIV>
<DIV><FONT size="4">Dear all</FONT></DIV>
<DIV><FONT size="4">&nbsp;&nbsp;&nbsp;&nbsp; In the forces_bp_efield.f90,one part of code is as following:</FONT></DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; !&nbsp; --- Check that we are working with an insulator with no empty bands ---<BR>&nbsp;&nbsp; IF ((degauss &gt; 0.01d0) .OR. (nbnd /= nelec/2)) &amp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 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&nbsp;is half metallic ,can I use this code to relax ? If my system have odd electrons ,how to set it ?</FONT><BR><BR><BR></DIV></DIV>
<DIV><BR>--<BR>
<DIV>Best regards</DIV>
<DIV>&nbsp;</DIV>
<DIV>Q.J.Wang</DIV>
<DIV>&nbsp;</DIV>
<DIV>XiangTan University </DIV></DIV><PRE>&nbsp;</PRE>