Hi, All,<br>I wonder how to do broken symmetry calculations for a hydrogen molecule with a separation larger than 3.0 a.u. using pw.x code. As I understand,&nbsp; the in the region larger then 3.0 a.u. separation, the restricted wavefountion method, for instance the RHF, would fail, and UHF method or correlation method is needed. I don&#39;t know how to break the symmetry of alpha spin and beta spin using pw.x code. I tried to do that by spin-polarized calculations using the following input. The energy is 0.05 a.u.lower than the sipn-unpolarized method at he hydrogen separation of 2.4333672 angstom. But at the sepration of 2.0100254, the spin-polarized method gave much energy then  sipn-unpolarized. Could you please check my input file? Mabye I have smething wrong.<br>
<br>Another question is about the binding energy of hydrgen molecule. I calculated the total engergy of H2 at 0.74 angstrom bond distanced and the energy H atom, and comapred the energy difference, I got a binding energy much higher than the experimental value of 436 kj/mol. I guess I did not calculate the energy of Hydrogen atom properly. Could you give me any suggestion how to do that? <br>
Thank you.<br><br>Qian<br><br><br><b>INPUT: Spin-polarized calculation</b><br>&amp;control<br>/<br>&amp;system<br>&nbsp;&nbsp;&nbsp; ibrav&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp; = 1<br>&nbsp;&nbsp;&nbsp; celldm(1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 30.0<br>&nbsp;&nbsp;&nbsp; nat&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; &nbsp; = 2<br>&nbsp;&nbsp;&nbsp; ntyp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; = 1<br>
&nbsp; &nbsp; nosym=.true.,<br>&nbsp; &nbsp; starting_magnetization(1)=1.0<br>&nbsp;&nbsp;&nbsp; nspin =2 <br>&nbsp;&nbsp;&nbsp; ecutwfc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 60.0<br>&nbsp;&nbsp;&nbsp; occupations = &#39;smearing&#39;<br>&nbsp; &nbsp; degauss = 0.01<br>&nbsp;&nbsp;&nbsp; smearing = &#39;gaussian&#39;<br>/<br>&amp;electrons<br>
&nbsp;&nbsp;&nbsp; conv_thr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0d-8<br>/<br>ATOMIC_SPECIES<br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp; 1.00&nbsp;&nbsp;&nbsp; H.vbc.UPF<br>ATOMIC_POSITIONS {angstrom}<br>&nbsp;&nbsp; H&nbsp;&nbsp; 0.0&nbsp;&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp; 0.0&nbsp; &nbsp; 0.0<br>&nbsp;&nbsp; H&nbsp;&nbsp; 2.4333672 &nbsp; 0.0&nbsp;&nbsp; 0.0<br>K_POINTS&nbsp; AUTOMATIC<br>1&nbsp;&nbsp; 1&nbsp;&nbsp; 1&nbsp;&nbsp; 0&nbsp;&nbsp; 0&nbsp;&nbsp; 0<br>
<br><br>