24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1528  |  回复: 11

schalke

至尊木虫 (著名写手)

[求助] 各位师兄师姐请进,为什么找完过渡态后IRC跑出来的峰顶部有问题?

模拟过渡态后,只有一个虚频,振动方向也对,但是IRC跑出的峰顶是平的,很奇怪,是什么原因?怎么改才是正确的?

输入文件如下:
%mem=8GB
%nprocshared=6
%chk=20120612_Dimer_acid_ts_freq_simplized.chk
#p irc=(maxpoints=10,rcfc) b3lyp/6-31+g(d,p) maxdisk=40GB

20120612_Dimer_acid_ts_irc_simplized

0 1
C                  0.40233500    3.21378000    1.96150600
H                  1.27882900    3.86279900    2.07424400
H                 -0.32329700    3.76141800    1.32885100
C                  2.08893900   -0.31068100   -0.77122100
H                  2.16849500   -0.60881300   -1.82569300
C                  3.23509200   -0.97147600    0.01411800
H                  3.17320900   -0.66520900    1.06593000
H                  3.09715400   -2.06200100    0.00465800
C                  4.62235200   -0.62845200   -0.54533100
H                  4.76005700    0.46287800   -0.53014900
H                  4.67482500   -0.92696400   -1.60287300
C                  5.77270100   -1.29211400    0.22351000
H                  5.71759400   -0.99927900    1.28216000
H                  5.63982200   -2.38379200    0.20308900
C                  7.16299200   -0.94348000   -0.32422400
H                  7.29699800    0.14706400   -0.30113500
H                  7.21798800   -1.23434100   -1.38263800
C                  8.30514200   -1.61338100    0.44852300
H                  8.29797200   -1.31457700    1.50359000
H                  9.28234900   -1.34373400    0.03316900
H                  8.21913400   -2.70607000    0.41444000
C                  0.73921800    1.91133000    1.32206100
H                  2.22621400    0.77479000   -0.75792400
C                  0.25690300    0.69658700    1.77251800
C                  0.70238900   -0.71931100   -0.26267300
C                  0.27353600   -0.52056000    1.07644100
H                 -0.06597500    3.07807400    2.94245300
C                 -0.54391900    1.73160400   -1.10104200
C                 -0.49407100    0.35047800   -1.44024500
C                 -1.78674100   -0.46014200   -1.63761200
H                 -1.54374700   -1.32679100   -2.26625500
H                 -2.46943600    0.15655000   -2.23969700
C                 -2.56147300   -0.97643100   -0.41167900
H                 -2.65381600   -0.19110600    0.34823600
H                 -1.99929300   -1.78831300    0.06378900
C                 -3.95763400   -1.49547500   -0.78372400
H                 -3.86376700   -2.25497000   -1.57423400
H                 -4.54650800   -0.67563900   -1.22031200
C                 -4.72514200   -2.09487100    0.40204300
H                 -4.81499000   -1.33914400    1.19636200
H                 -4.13957400   -2.92012100    0.83325200
C                 -6.12371800   -2.60882400    0.03527700
H                 -6.03428700   -3.36110300   -0.76087100
H                 -6.71169700   -1.78380700   -0.39045800
C                 -6.87941300   -3.21312200    1.22443600
H                 -7.87188700   -3.57056800    0.92886200
H                 -7.01619200   -2.47484700    2.02368300
H                 -6.33280400   -4.06310300    1.65015800
C                 -1.70533700    2.44913300   -0.65386200
H                 -2.56681600    1.86380600   -0.33504700
C                 -1.82616900    3.80196000   -0.62268100
H                 -0.98652000    4.40960000   -0.96381000
C                 -3.07050300    4.54248200   -0.22484700
H                 -3.47401200    5.12732800   -1.06329900
H                 -2.87823800    5.25724000    0.58635200
H                 -3.85424200    3.85553800    0.11160600
H                  1.47452100    1.94070000    0.53045900
H                 -0.32104500    0.71418300    2.69794200
H                 -0.31814600   -1.31824700    1.51877200
H                  0.46032600   -1.72375600   -0.62145300
H                  0.19027900    0.18268200   -2.27161900
H                  0.28564700    2.34742700   -1.44234600
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

schalke

至尊木虫 (著名写手)

人工置顶,求解啊
2楼2012-06-15 09:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

第一次 见 开眼了...我想知道 top的2个结构 是一样的么...把ts再优化一下?
3楼2012-06-15 09:25:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

schalke

至尊木虫 (著名写手)

引用回帖:
3楼: Originally posted by meteoric30 at 2012-06-15 09:25:23
第一次 见 开眼了...我想知道 top的2个结构 是一样的么...把ts再优化一下?

过渡态已经优化完毕,而且过渡态结构对应的是top2个点的中间位置,但并没有在IRC中显示出过渡态的点,但计算结果中含有过渡态结构。。。我是先算出过渡态的频率,然后读取其chk文件做irc,是否和读取chk文件有关?
4楼2012-06-15 09:39:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

引用回帖:
4楼: Originally posted by schalke at 2012-06-15 09:39:25
过渡态已经优化完毕,而且过渡态结构对应的是top2个点的中间位置,但并没有在IRC中显示出过渡态的点,但计算结果中含有过渡态结构。。。我是先算出过渡态的频率,然后读取其chk文件做irc,是否和读取chk文件有关?...

2点的能量和过渡态的能量也不一样吗?我意思 还是再优化下ts 再freq
5楼2012-06-15 17:14:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

schalke

至尊木虫 (著名写手)

引用回帖:
5楼: Originally posted by meteoric30 at 2012-06-15 17:14:48
2点的能量和过渡态的能量也不一样吗?我意思 还是再优化下ts 再freq...

不一样,过渡态的结构能量最高,应该在两点中间。我是做完ts和freq后,用结果做的irc。我现在改用rcfc再计算一遍看看。另外,我算过渡态的时候用了nosymm、noeigentest、maxcyc的指令,但算irc的时候都没有用,会不会影响结果?
6楼2012-06-15 19:05:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

引用回帖:
6楼: Originally posted by schalke at 2012-06-15 19:05:19
不一样,过渡态的结构能量最高,应该在两点中间。我是做完ts和freq后,用结果做的irc。我现在改用rcfc再计算一遍看看。另外,我算过渡态的时候用了nosymm、noeigentest、maxcyc的指令,但算irc的时候都没有用,会不 ...

应该不会~~你有没有重新 算一下freq
7楼2012-06-16 09:18:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

schalke

至尊木虫 (著名写手)

引用回帖:
7楼: Originally posted by meteoric30 at 2012-06-16 09:18:29
应该不会~~你有没有重新 算一下freq...

我先算了初始构型的频率,然后利用chk文件优化过渡态,再计算过渡态的频率,然后利用过渡态chk文件计算irc
8楼2012-06-16 09:20:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

引用回帖:
8楼: Originally posted by schalke at 2012-06-16 09:20:52
我先算了初始构型的频率,然后利用chk文件优化过渡态,再计算过渡态的频率,然后利用过渡态chk文件计算irc...

没问题~~我意思 再优化下过渡态~算freq 再看看irc应该没什么问题了
9楼2012-06-16 09:32:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

schalke

至尊木虫 (著名写手)

引用回帖:
9楼: Originally posted by meteoric30 at 2012-06-16 09:32:57
没问题~~我意思 再优化下过渡态~算freq 再看看irc应该没什么问题了...

我在linux下面算了一次,windows下面算了一次,都是这结果,太诡异。。。
10楼2012-06-16 09:34:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 schalke 的主题更新
信息提示
请填处理意见