24小时热门版块排行榜

返回列表

fjxmwgh

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 174.8
红花: 1
帖子: 34
在线: 31.6小时
虫号: 3318721
注册: 2014-07-11
专业: 无机纳米化学

[求助] Gaussian 【势能面】【新手】为什么我的过渡态势能面扫描没有鞍点？

这个过渡态我扫描过freq，只有一个虚频；然后IRC扫出来能量作图是一个馒头峰的形状（据说这样就表示这个的确是过渡态？）
接着开始扫描，我用scan和opt(modredundant)都扫过，然后用matlab查看势能面。scan扫出来直接一个斜坡连最低点都没有，opt(modredundant)扫出来是一个最低点，都没有马鞍形T T。(我的坐标都是直接从CHEM3D里面拷出来的，然后按照标的原子序数调好顺序)，没有Gaussian View。
以下是modredundant的
我是把扫描的过渡态的两个变量稍微调小了一点然后扫描的modredundant，都覆盖了理论值
#B3LYP/6-31G opt(modredundant)

T1

0 1
C 2.868 -0.138 0.265
C 1.611 -1.065 0.402
C 0.576 -0.286 -0.394
N 0.824 1.143 0.074
C 2.314 1.299 -0.026
C -0.962 -0.372 -0.352
O -1.602 -1.422 -0.182
O -1.420 0.864 -0.574
O 3.666 -0.467 -0.906
H 3.489 -0.162 1.169
H 1.308 -1.173 1.451
H 1.806 -2.061 0.001
H 0.873 -0.295 -1.451
H 0.551 1.205 1.062
H 2.587 1.571 -1.047
H 2.666 2.061 0.669
H 4.176 -1.289 -0.763
H -0.159 0.968 -0.419
O -4.176 0.199 -0.052
H -3.824 -0.712 -0.022
H -3.409 0.751 -0.314

17 4 S 10 0.02
17 8 S 10 0.02

以下是scan的
#B3LYP/6-31G scan

t1

0 1
C
C 1 1.568
C 2 1.521 1 101.354
N 3 1.524 2 102.041 1 -43.175
C 1 1.567 2 105.936 3 23.912
C 3 1.541 2 129.534 4 101.961
O 6 1.242 3 124.474 2 -36.040
O 1 1.455 2 112.118 5 104.600
H 1 1.097 2 111.589 5 111.921
H 2 1.097 1 111.310 3 111.277
H 2 1.092 1 111.366 3 113.355
H 3 1.099 2 108.393 4 105.066
H 4 1.027 3 107.961 5 108.833
H 5 1.091 1 108.256 4 109.719
H 5 1.090 1 114.095 4 110.541
H 4 R1 3 91.51 6 -33.21
O 16 R2 3 106.605 7 128.886
H 8 0.978 1 111.009 2 -74.833
O 17 2.883 6 94.853 3 -173.051
H 19 0.978 1 66.082 2 -10.243
H 19 0.981 20 104.520 1 17.603

R1 1.1000 10 0.02
R2 1.2000 10 0.02

1-10就要交作业啦T T 大神们求助啊！！已经扫了一周失败了
我现在也不知道怎么弄，就不停地换变量尝试，键长键角，二面角什么的都试过了(之前主要尝试scan，接下来想试modredundant，主要modredundant太慢了)，scan的结果总是一个斜坡，连最低点都没有……

贴图是modredundant在matlab里面扫出来的东东(matlab程序老师编的)，其实那个modredundant看上去还蛮像马鞍形的是不是我把势能面扫大一点就能出来形状 - -

untitled.jpg

回复此楼