<div dir="ltr"><div>Hello,</div><div><br></div><div>I am trying to calculate bond dissociation energy (BDE) for the following reaction   </div><div>N-N-H(2+) --> N-N(2+) + H</div><div>using both CCSD(T) and full valence MRCI methods. Because H-atom has only one electron I am using SCF energy for it. BDE I got is around 205 kcal/mol  which looks too large for me (for example, BDE of the following reaction N-N-H(+) --> N-N + H(+) is 116.9 kcal/mol). You can see inputs below. What I am doing wrong?</div><div><br></div><div>Thank you in advance  David</div><div><br></div><div>N-N(2+) MRCI input  </div><div><br></div><div><div>***,NN++</div><div>memory,350,m</div><div>gprint,orbital,civector;</div><div>gthresh,energy=1.d-12</div><div> r = 1.13662754 ang;</div><div>geometry={N1              !z-matrix geometry input</div><div>          N2,N1,r}</div><div>basis=aug-cc-pVTZ;</div><div>{hf;accu,16;</div><div>wf,12,1,0;</div><div>print,orbprint=32}</div><div>{multi;</div><div>maxiter,39;</div><div>wf,12,1,0}</div><div>{mrci,maxiti=200;</div><div>wf,12,1,0}</div><div>{optg,gradient=1.d-5,fourpoint}</div><div>pop</div><div><br></div><div>N-N(2+) CCSD(T) input</div><div><div>***,NN++</div><div>memory,350,m</div><div>gprint,orbital,civector;</div><div>gthresh,energy=1.d-12</div><div> r = 1.08644919 ang;</div><div>geometry={N1              !z-matrix geometry input</div><div>          N2,N1,r}</div><div>basis=aug-cc-pVTZ;</div><div>{hf;accu,16;</div><div>wf,12,1,0;</div><div>print,orbprint=32}</div><div>{uccsd(t),maxit=100;</div><div>wf,12,1,0}</div><div>{optg,gradient=1.d-5,fourpoint}</div><div>pop</div></div></div><div><br></div><div>N-N-H(2+) MRCI input </div><div><br></div><div><div>***,HNN++</div><div>memory,350,m</div><div>gprint,orbital,civector;</div><div>gthresh,energy=1.d-12</div><div> r = 1.09810143 ang;</div><div> r1 = 1.14511161 ang;</div><div>geometry={N1              !z-matrix geometry input</div><div>          N2,N1,r;</div><div>          H1,N1,r1,N2,180.0}</div><div>basis=aug-cc-pVTZ;</div><div>{hf;accu,16;</div><div>wf,13,1,1;</div><div>print,orbprint=32}</div><div>{multi;</div><div>maxiter,39;</div><div>wf,13,1,1}</div><div>{mrci,maxiti=200;</div><div>wf,13,1,1}</div><div>{optg,gradient=1.d-5,fourpoint}</div><div>pop</div></div><div><br></div><div><br></div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div>____________________________________________________________________________<br>Dr. David Danovich,  The Lise Meitner Minerva Center for Computational Quantum Chemistry,</div>
<div>Institute of Chemistry, The Hebrew University, Edmond J. Safra Campus - Givat Ram,</div>
<div>91904 Jerusalem, Israel <a href="http://yfaat.ch.huji.ac.il/david.html" target="_blank">http://yfaat.ch.huji.ac.il/david.html</a>, <a href="mailto:david.danovich@huji.ac.il" target="_blank">david.danovich@huji.ac.il</a> </div>
<div>FAX:(+972)-2-6584033, Phone:(+972)-2-6586934(w),  Mobile:(+972)-544-768669</div></div></div>
</div>