24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3551  |  回复: 41
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

jianqin

金虫 (正式写手)

[求助] Gaussian几何优化后的键长不合理

用Gaussian 09做分子几何优化。分子的部分结构见图1,具有两个手性中心(C3和C6)。几何优化分步过程如下:先用chem3D MM2做初步优化,再用PM3优化(包括频率分析),接下来用RHF 6-31G做优化(包括频率分析),计算均成功完成。但是分子中的碳硫键C3-S1(图1),经PM3优化后,键长1.882埃,基本合理,经RHF 6-31G优化后竟然变成3.172埃,显然是不合理的。我的疑问如下:
我的这种分步优化应该是合理的吧?
键长不合理的情况是否与分子中包含手性中心有关?因为C3是手性的,导致gaussian对碳硫键的处理不合理?(另一个手性中心C6周围的键长均在合理范围内)
gaussian中对手性分子的几何优化有什么需要特别注意的吗?
6-31G的基组对于处理硫这样的原子足够吗?
补充:刚看到其他的帖子,认为初始分子的构型是否合理,对找到正确的优化结构很关键?那么我的分步优化产生的构型是否是比较合理的构型呢?


[ Last edited by jianqin on 2011-7-15 at 09:40 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianqin

金虫 (正式写手)

引用回帖:
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
21楼2011-07-18 12:51:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 42 个回答

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
小红豆(金币+2): 谢谢指导 2011-07-15 19:21:15
这分子又不大,起码加个*吧... 6-31g* 即 6-31g(d) ...
2楼2011-07-15 10:20:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianqin

金虫 (正式写手)

引用回帖:
Originally posted by fatpig8832 at 2011-07-15 10:20:09:
这分子又不大,起码加个*吧... 6-31g* 即 6-31g(d) ...

我是准备HF 6-31G优化后,再用B3LYP 6-311+G(d,p)优化
3楼2011-07-15 10:33:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★
小红豆(金币+3): 感谢交流 2011-07-15 19:21:30
引用回帖:
Originally posted by jianqin at 2011-07-15 10:33:52:
我是准备HF 6-31G优化后,再用B3LYP 6-311+G(d,p)优化

可以直接把pm3的结果用你的大基组优化,不一定要先算小基组,要不你可以先试试6-31g*有没有什么不同
4楼2011-07-15 12:00:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211 0703化学 346分求调剂 +24 土豆er? 2026-04-09 26/1300 2026-04-13 11:49 by 电化学及催化
[考研] 化工学硕 285求调剂 +28 Wisjxn 2026-04-07 28/1400 2026-04-13 10:56 by 学zh
[考研] 一志愿南科大生物学297分,求调剂推荐 +10 Y-yyusx 2026-04-06 11/550 2026-04-13 09:30 by mfyue2000
[考研] 材料299专硕求调剂 +11 +21 2026-04-09 11/550 2026-04-13 08:57 by lhj2009
[考研] 材料考研调剂 +28 云木达达 2026-04-11 30/1500 2026-04-13 08:51 by Sammy2
[考研] 材料工程085601,270求调剂 +36 @ASDF1234 2026-04-08 39/1950 2026-04-13 08:37 by Delta2012
[考研] 本科郑州大学,一志愿华东师范大学282求调剂 +24 熊哥xtk 2026-04-07 27/1350 2026-04-13 08:30 by lhj2009
[考研] 290调剂生物0860 +28 哇哈哈,。 2026-04-11 31/1550 2026-04-13 01:16 by 幸免 ..
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen~
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 22408调剂求助 +7 毂12 2026-04-09 9/450 2026-04-11 09:23 by 哦哦123
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 083200 305分 求二轮调剂 不接受跨专业 +9 Claireyyyy 2026-04-09 10/500 2026-04-10 21:21 by Claireyyyy
[考研] 284求调剂 +19 梵@@ 2026-04-06 21/1050 2026-04-10 21:12 by zhouxiaoyu
[考研] 22408 366分,本科211,一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 一志愿中科院105500专业总分315求调剂 +6 lallalh 2026-04-09 7/350 2026-04-09 17:51 by lallalh
[考研] 一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药 +7 超人不会飞@ 2026-04-08 7/350 2026-04-09 09:35 by gong120082
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 专硕085403,291分,有两篇专利,一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
信息提示
请填处理意见