Dear Molpro developers / users,<div><br></div><div>I need to calculate the numerical values of IRC coordinate for some paths. Although Molpro is able to print them by input request, it would be meaningless when restarting a qsdpath, which is the case, so I have to do it by myself.</div>
<div><br></div><div>I've tried some possible formulas involving the RMS between the transition state and the desired points, but I wasn't successful reproducing some values obtained in single-run complete qsdpaths. </div>
<div><br></div><div>It's not clear to me whether the IRC coordinate printed is:</div><div><br></div><div>- mass-weighted;</div><div>- in atomic units or in user-defined units;</div><div>- calculated through internal or cartesian coordinates;</div>
<div>- averaged over the number of atoms or over the number of coordinates</div><div><br></div><div>I've tried some combinations and looked over the source code, but couldn't get it.</div><div><br></div><div>Any help would be very appreciated.</div>
<div><br></div><div>Best regards,</div><div><br></div><div>-- <br>Gabriel do Nascimento Freitas<br>D.Sc. Student - Graduate Program of Chemistry<br>Molecular Modelling and Theoretical Chemistry Laboratory - Room 412<br>
Chemistry Institute - Universidade Federal do Rio de Janeiro (UFRJ) - Brazil <br><a href="tel:%28%2B5521%298830-9971" value="+552188309971" target="_blank">(+5521)8830-9971</a> / <a href="tel:%28%2B5521%292562-7179" value="+552125627179" target="_blank">(+5521)2562-7179</a></div>