Hi， Do you read the input document of pwscf?<br>it explicitly claims that&nbsp; tot_charge is used only if nelec is unspecified, and otherwise it is ignored !<br><br><div class="gmail_quote">On Dec 24, 2007 10:47 AM, oulihui666 &lt;
<a href="mailto:oulihui666@126.com">oulihui666@126.com</a>&gt; wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div>Dear pwscf users,&nbsp;&nbsp; 
<br>&nbsp;&nbsp;&nbsp; I have asked the similar problem last time, but I do not understand this </div>
<div>problem very explicitly. I want to do a structural relaxation of a charged </div>
<div>system using pwscf. For example, In the process of oxygen reduction reaction, </div>
<div>adsorption of hydration proton(H3O+) on Pt(111) surface,&nbsp; simultaneously, a </div>
<div>extra electron is added, that is to say, O2 + H+ + e- . When I calculate this </div>
<div>system with pwscf, could I add charge and add a electron simultanously? My input </div>
<div>file is displayed as follows, In this system, there is 157 electrons, when I add </div>
<div>a extra electron, it is 158 electron, simultanously, I add +1 charge, when I run </div>
<div>input file, I do not get an error information, is this correct? I need your </div>
<div>help.&nbsp;&nbsp; <br>Thanks in advance&nbsp;&nbsp;&nbsp;&nbsp; <br>Lihui Ou&nbsp;&nbsp; </div>
<div>&nbsp;&amp;CONTROL<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation = &#39;relax&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = &#39;from_scratch&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = &#39;/home/olh/tmp/&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = &#39;/home/olh/pseudo/&#39; ,
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = &#39;O2-H7O3-pt111-relax&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nstep = 200 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tstress = .true. ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tprnfor = .true. ,<br>&nbsp;/<br>&nbsp;&amp;SYSTEM<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 4,
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = 10.486,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(3) = 2.979,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat = 24,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 3,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 30 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 300 ,
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nelec = 158,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tot_charge = 1.000000,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; occupations = &#39;smearing&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.02 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smearing = &#39;methfessel-paxton&#39; ,
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nspin = 2 ,<br>&nbsp;&nbsp; starting_magnetization(1) = 0.2,<br>&nbsp;&nbsp; starting_magnetization(2) = 0.5,<br>&nbsp;&nbsp; starting_magnetization(3) = 1.0,<br>&nbsp;/<br>&nbsp;&amp;ELECTRONS<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.D-6 ,
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.4D0 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; diagonalization = &#39;david&#39; ,<br>&nbsp;/<br>&nbsp;&amp;IONS<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ion_dynamics = &#39;bfgs&#39; ,<br>&nbsp;/<br>ATOMIC_SPECIES<br>&nbsp;&nbsp; Pt&nbsp; 195.08000&nbsp; Pt.pbe-nd-rrkjus.linux.UPF
 <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp; 15.99900&nbsp; O.pbe-rrkjus.UPF <br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 1.00800&nbsp; H.pbe-rrkjus.linux.UPF <br>ATOMIC_POSITIONS angstrom <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.775000000&nbsp;&nbsp;&nbsp; 1.602000000&nbsp;&nbsp; -0.000000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0 <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 3.204000000
&nbsp;&nbsp;&nbsp; 2.265000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0 <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 4.531000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.387000000&nbsp;&nbsp;&nbsp; 4.005000000&nbsp;&nbsp; -0.000000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0 <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.387000000&nbsp;&nbsp;&nbsp; 0.801000000&nbsp;&nbsp;&nbsp; 2.265000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0 
<br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp; -1.387000000&nbsp;&nbsp;&nbsp; 2.403000000&nbsp;&nbsp;&nbsp; 4.531000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp; -1.387000000&nbsp;&nbsp;&nbsp; 4.005000000&nbsp;&nbsp; -0.000000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0 <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4.162000000&nbsp;&nbsp;&nbsp; 0.801000000&nbsp;&nbsp;&nbsp; 2.265000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0 <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.387000000
&nbsp;&nbsp;&nbsp; 2.403000000&nbsp;&nbsp;&nbsp; 4.531000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 1.602000000&nbsp;&nbsp; -0.000000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0 <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.775000000&nbsp;&nbsp;&nbsp; 3.204000000&nbsp;&nbsp;&nbsp; 2.265000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0 <br>&nbsp;&nbsp; Pt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.775000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 
4.531000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.238000000&nbsp;&nbsp;&nbsp; 0.876000000&nbsp;&nbsp;&nbsp; 6.331000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.209000000&nbsp;&nbsp;&nbsp; 1.500000000&nbsp;&nbsp;&nbsp; 6.331000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.443000000&nbsp;&nbsp;&nbsp; 1.185000000&nbsp;&nbsp;&nbsp; 9.230000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 
<br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.605000000&nbsp;&nbsp;&nbsp; 1.362000000&nbsp;&nbsp;&nbsp; 8.820000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.363000000&nbsp;&nbsp;&nbsp; 2.268000000&nbsp;&nbsp;&nbsp; 8.975000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.670000000&nbsp;&nbsp;&nbsp; 1.207000000&nbsp;&nbsp;&nbsp; 7.885000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp; -0.878000000
&nbsp;&nbsp;&nbsp; 3.907000000&nbsp;&nbsp; 11.663000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp; -1.128000000&nbsp;&nbsp;&nbsp; 3.097000000&nbsp;&nbsp; 11.233000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp; -2.045000000&nbsp;&nbsp;&nbsp; 2.917000000&nbsp;&nbsp; 11.406000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.018000000&nbsp;&nbsp;&nbsp; 2.594000000&nbsp;&nbsp; 
11.726000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.652000000&nbsp;&nbsp;&nbsp; 1.744000000&nbsp;&nbsp; 11.513000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.797000000&nbsp;&nbsp;&nbsp; 1.146000000&nbsp;&nbsp; 12.237000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <br>K_POINTS automatic <br>&nbsp; 4 4 1&nbsp;&nbsp; 1 1 1 <br></div>
<div><br>--<br></div>
<div>
<div>
<div><font face="Arial" size="3">======================================</font></div>
<div><font face="Arial" size="3">Lihui&nbsp;Ou</font></div>
<div><font face="Arial" size="3">PH.D Candidate&nbsp;in Electrochemistry&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </font></div>
<div><font face="Arial" size="3">College of Chemistry and Molecular Science&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </font></div>
<div><font face="Arial" size="3">Wuhan University,<i>430072</i>,Hubei Province,China&nbsp;</font></div>
<div><i><font face="Arial" size="3">E-mail:</font><a href="mailto:oulihui666@126.com" target="_blank"><font face="Arial" size="3">oulihui666@126.com</font></a></i></div>
<div><font face="Arial" size="3">======================================</font></div></div></div><br>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org
</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br><br></blockquote></div><br><br clear="all"><br>-- <br>Hai-Ping Lan <br>Department of Electronics ,
<br>Peking University , Bejing, 100871<br><a href="mailto:lanhaiping@gmail.com">lanhaiping@gmail.com</a>, <a href="mailto:hplan@pku.edu.cn">hplan@pku.edu.cn</a>