<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 = 'relax' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/home/olh/tmp/' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '/home/olh/pseudo/' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'O2-H7O3-pt111-relax' ,<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 = 'smearing' ,<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 = 'methfessel-paxton' ,<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 = 'david' ,<BR>&nbsp;/<BR>&nbsp;&amp;IONS<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ion_dynamics = 'bfgs' ,<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,<EM>430072</EM>,Hubei Province,China&nbsp;</FONT></DIV>
<DIV><EM><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></EM></DIV>
<DIV><FONT face="Arial" size="3">======================================</FONT></DIV></DIV></DIV>