|
在学习polyrate时发现计算速率需要MEP的相关信息(hessian矩阵和梯度等),为此,我针对我所研究的体系进行了反应路径计算(#p IRC(cartesian,calcall,maxpoints=52,stepsize=2,verytight) M062X/6-311++G** empiricaldispersion=gd3 int=superfine),随后取-1<s<1范围内的各结构进行频率分析(#p freq M062X/6-311++G** empiricaldispersion=gd3 nosymm int=superfine),结果发现部分点存在除沿反应路径的虚频以外的小虚频。再查看polyrate中CH3+H2→CH4+H的例子发现,其反应路径计算中也存在小虚频。故我想请教各位老师,这种小虚频是否需要去掉以及怎么去掉。
polyrate中的部分反应路径信息:
Classical and adiabatic energies (kcal/mol), effective CD-SC reduced
mass (mu) (a.u.), and frequencies (cm**-1) vs. s
s(angstrom) VMEP Va^G mu^CD-SC frequencies (cm**-1)
0.106 21.1453 52.9224 1044.95 3787 3787 3569 1824 1824 1721 1663
1663 1584 403 403
0.132 19.6131 51.2482 881.12 3786 3786 3571 1840 1840 1798 1677
1676 1608 274 273
0.159 18.0370 49.1829 685.63 3785 3785 3574 1945 1855 1854 1683
1683 1623 29i 30i
0.185 16.5311 48.0715 439.34 3785 3785 3578 2179 1866 1867 1684
1684 1633 249i 249i
0.212 15.1765 47.2412 280.42 3786 3786 3585 2519 1875 1875 1682
1682 1641 337i 337i
|
|