Gaussian几何优化后的键长不合理 - 量子化学 - Gaussian - 第3页小木虫论坛-学术科研互动平台

版块导航: 正在加载中...

客户端APP下载

论文辅导

申博辅导

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

jianqin

金虫 (正式写手)

应助: 7 (幼儿园)
金币: 3728.4
红花: 2
帖子: 627
在线: 394小时
虫号: 863555
注册: 2009-10-05
专业: 质谱分析

引用回帖:

Originally posted by zhangmt at 2011-07-18 12:17:53:
别隐藏H啊，你直接贴gjf文件吧。看起来方便些。。这样看不方便。

%mem=2GB
%nproc=4
%chk=298.chk
# opt freq rb3lyp/6-31+g(d,p) geom=connectivity iop(1/8=10)

298 B3LYP 6-31+G(d p)

1 1
N                2.22760000 1.72230000 1.05790000
C                1.83990000 1.04100000 -0.16230000
C                2.03220000 -0.30580000 0.00070000
C                2.69640000 -0.51510000 1.35790000
C                2.71530000 0.86840000 1.94170000
C                1.37260000 1.83320000 -1.32620000
O                0.65420000 1.50310000 -2.24220000
O                1.85080000 3.10450000 -1.31520000
C                2.01170000 -1.52080000 2.26630000
S                1.78010000 -1.63190000 -1.09060000
N                -1.05170000 -1.17570000 -1.58920000
C                -0.03270000 -2.13170000 -1.17050000
C                -0.54220000 -2.78790000 0.11500000
C                -1.93870000 -2.22800000 0.38070000
C                -2.10220000 -1.01740000 -0.54760000
C                -2.02880000 0.29220000 0.24680000
N                -3.17820000 1.08830000 0.40970000
O                -0.96860000 0.71730000 0.71280000
C                -4.52800000 0.52100000 0.19490000
C                -3.13680000 2.24610000 1.33200000
H                2.12800000 2.69930000 1.18130000
H                3.75900000 -0.84250000 1.18900000
H                3.07550000 1.11010000 2.95210000
H                1.55940000 3.58050000 -2.08770000
H                2.23790000 -1.32120000 3.32170000
H                2.35480000 -2.54110000 2.04090000
H                0.91560000 -1.51540000 2.13980000
H                -0.73190000 -0.33980000 -2.02400000
H                0.11560000 -2.86170000 -2.01020000
H                -0.55840000 -3.89000000 0.01830000
H                0.13400000 -2.57400000 0.97260000
H                -2.07050000 -1.95970000 1.44550000
H                -2.71660000 -2.98640000 0.16780000
H                -3.07320000 -1.06790000 -1.11020000
H                -4.51470000 -0.22770000 -0.61370000
H                -5.21720000 1.32420000 -0.10090000
H                -4.92180000 0.04440000 1.10420000
H                -4.05570000 2.83880000 1.22480000
H                -3.04440000 1.93280000 2.38180000
H                -2.28440000 2.89360000 1.08750000

1 2 1.0 5 2.0 21 1.0
2 3 2.0 6 1.0
3 4 1.0 10 1.0
4 5 1.0 9 1.0 22 1.0
5 23 1.0
6 7 2.0 8 1.0
7
8 24 1.0
9 25 1.0 26 1.0 27 1.0
10 12 1.0
11 12 1.0 15 1.0 28 1.0
12 13 1.0 29 1.0
13 14 1.0 30 1.0 31 1.0
14 15 1.0 32 1.0 33 1.0
15 16 1.0 34 1.0
16 17 1.0 18 2.0
17 19 1.0 20 1.0
18
19 35 1.0 36 1.0 37 1.0
20 38 1.0 39 1.0 40 1.0
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40

赞一下(1人)

回复此楼

21楼2011-07-18 12:51:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jianqin

金虫 (正式写手)

应助: 7 (幼儿园)
金币: 3728.4
红花: 2
帖子: 627
在线: 394小时
虫号: 863555
注册: 2009-10-05
专业: 质谱分析

引用回帖:

Originally posted by coolrainbow at 2011-07-18 12:02:55:
含有S这种原子的体系必须加极化才行
比如6-311g(d)等等
如果非要在win下算也没关系，优化个三四天不算长

用opt freq DFT B3LYP 6-31+g(d,p)做了优化，gjf文件见楼上。
离子结构式（含手性中心）如图
PM3优化后的结构见下图，碳硫键键长在合理范围内。
我自己的分析是：由于氮上带正电荷，导致氮正离子旁边的羧基取向与环共平面，与吡咯环边上的羰基之间有斥力作用，导致碳硫键被拉长了，我再用不带电荷的分子再做一下优化，试一试。

赞一下

回复此楼

22楼2011-07-18 13:54:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖