求助：分子断键形成两个自由基的反应计算

徐丹丹

新手求助：

大家好，我是刚接触量化计算的小白。想求助计算一个分子断键形成两个自由基时候的过渡态，并绘制势能面图。
C6-C8断键形成两个自由基。


以下是我寻找过渡态时的输入文件，请教大家是否正确。万分感谢！！

%nprocshared=36
%mem=100GB
%chk=TS1.chk
# opt=(calcall,ts,noeigen) freq m062x/6-31++g(d,p) temperature=673

Title Card Required

0 1
C                 -4.99697140    1.74769941    0.36747127
C                 -4.11413040    0.58191041   -0.08504173
H                 -6.00186940    1.40617941    0.63679527
H                 -4.56615840    2.25810041    1.23474727

（省略）

chands

谐振近似下设置温度没用，temperature=673去掉。
6-31++G(d,p)可改为6-31+G(d,p)，如果你想得到更好的结果可以后续再用大点的基组对TS优化。另外注意这个反应有一定概率没有过渡态。这时要考虑无过渡态反应（论坛主页搜“无过渡态”“无势垒”）

sobereva

用6-311G**比用6-31++g(d,p)强得多

好好看此文
谈谈弥散函数和“月份”基组
http://sobereva.com/119

这种普通共价键断裂，而且此过程中又没有额外能量降低因素的情况，没有过渡态

写temperature=673没意义。一方面不影响优化结果，另一方面你也不可能直接用Gaussian在这个温度下给出的热力学校正量，这种高度柔性含很低频的体系起码也得用Shermo（http://sobereva.com/552）基于Quasi-RRHO模型才能算得有意义

wjc404

(不好意思误操作重复发帖了，请管理员删除)

wjc404

新手就要做这个问题有些勉为其难。这个问题的难点除了楼上提到的可能没有过渡态以及侧链构象多之外，还有断键过程波函数多参考特征等。可以先听听社长的初级量子化学培训班，然后上手一些简单的计算任务（诸如闭壳层分子的优化和振动分析，简单双原子分子的势能曲线等），一步一步来。