
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">

<HTML>

<HEAD>

<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">

<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7638.1">

<TITLE>Charge is wrong ...in the output file </TITLE>

</HEAD>

<BODY>

<!-- Converted from text/plain format -->

<BR>


<P><FONT SIZE=2>Hi QEs users,<BR>

<BR>

I got an error in my output file as follows<BR>

<BR>

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<BR>

&nbsp;&nbsp;&nbsp;&nbsp; from electrons : error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<BR>

&nbsp;&nbsp;&nbsp;&nbsp; charge is wrong<BR>

&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<BR>

However from my previous survey in PW_forum, I have found that this problem may be due to not doing spin polarized calculation with an odd number of electrons, what does it mean?<BR>

<BR>

My input file is<BR>

<BR>

&amp;CONTROL<BR>

&nbsp;&nbsp;&nbsp;&nbsp; calculation ='scf'<BR>

&nbsp;&nbsp;&nbsp; restart_mode ='from_scratch'<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/sfihome/dsharma/Pscf/test'<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '/sfiwork/dsharma/QE_Pseudos/'<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'Pscf1'<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tstress = .true.<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tprnfor = .true.<BR>

&nbsp;&nbsp; etot_conv_thr = 1.D-4<BR>

&nbsp;&nbsp; forc_conv_thr = 1.D-3<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nstep = 600<BR>

/<BR>

&nbsp;&amp;SYSTEM<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 0<BR>

&nbsp;&nbsp;&nbsp; celldm(1) =20.40066<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat = 94<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 3<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 30<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 120<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nbnd = 220<BR>

/<BR>

&nbsp;&amp;ELECTRONS<BR>

&nbsp;&nbsp;&nbsp;&nbsp; diagonalization ='cg'<BR>

&nbsp;&nbsp;&nbsp; mixing_mode = 'plain'<BR>

&nbsp;&nbsp;&nbsp; conv_thr = 1.0d-6<BR>

&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.7<BR>

/<BR>

&nbsp;CELL_PARAMETERS cubic<BR>

<BR>

&nbsp;&nbsp;&nbsp;&nbsp; 1.000000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 0.000000000<BR>

&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 1.000000000&nbsp;&nbsp;&nbsp; 0.000000000<BR>

&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 5.000000000<BR>

ATOMIC_SPECIES<BR>

&nbsp;&nbsp; Si&nbsp;&nbsp; 28.08600&nbsp; Si.pz-vbc.UPF<BR>

&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 1.00800&nbsp; H.pz-vbc.UPF<BR>

&nbsp;&nbsp;&nbsp; O&nbsp; 15.9994&nbsp;&nbsp;&nbsp; O.pz-rrkjus.UPF<BR>

&nbsp;ATOMIC_POSITIONS angstroms<BR>

O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7.0781&nbsp;&nbsp;&nbsp; 6.9497&nbsp;&nbsp; 27.1761<BR>

O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4.3952&nbsp;&nbsp;&nbsp; 6.9992&nbsp;&nbsp; 27.2391<BR>

H&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7.4599&nbsp;&nbsp;&nbsp; 6.0212&nbsp;&nbsp; 27.1358<BR>

H&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5.3921&nbsp;&nbsp;&nbsp; 7.0051&nbsp;&nbsp; 27.1349<BR>

H&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.5615&nbsp;&nbsp;&nbsp; 4.4933&nbsp;&nbsp; 27.1304<BR>

Si&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.4277&nbsp;&nbsp;&nbsp; 0.5123&nbsp;&nbsp; 26.6794<BR>

Si&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7.2634&nbsp;&nbsp;&nbsp; 3.3235&nbsp;&nbsp; 26.6750<BR>

Si&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.4470&nbsp;&nbsp;&nbsp; 0.4413&nbsp;&nbsp; 26.6627<BR>

........<BR>

<BR>

K_POINTS automatic<BR>

&nbsp; 2 2 1&nbsp;&nbsp; 0 0 0<BR>

<BR>

&nbsp;&nbsp;&nbsp;&nbsp; Can anyboy please help?<BR>

DS.<BR>

</FONT>

</P>


</BODY>

</HTML>