计算化学公社

标题: 关于是否需要消除MEP上各结构小虚频的疑问 [打印本页]

作者
Author: 400__LUX 时间: 2025-7-3 10:49
标题: 关于是否需要消除MEP上各结构小虚频的疑问
在学习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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)