I am sorry to bother but I am a little bit confused about how to run rate calculations with Py-Chemshell.
I didn’t see any section on rate calculations in the Py-ChemShell manual but I was looking at the Tcl one trying to make it work. Once reactant and ts have been calculated, my doubts are:
How can I let the input file read the qts_reactant.txt, qts_ts.txt, and class.in files? Is there a specific keyword for reading “.txt” and “.in” files?
Is it necessary to read the coordinates of the reactant and TS state with .xyz/.pun files?
Thank you very much in advance for your help
from chemsh import *
rs = Fragment(coords=“reactant_opt.pun”)
ts = Fragment(coords=“ts.pun”)
theory = ORCA(frag=rs, method=“dft”, functional=“blyp”, basis=“def2-svp”)
opt = Opt(theory=theory, frag2=ts, thermal=True, rate=True, inithessian=5)