24小时热门版块排行榜    

查看: 1254  |  回复: 11
本帖产生 1 个 QC强帖 ,点击这里进行查看

renren123123

金虫 (正式写手)

[求助] 分子优化求助

# opt freq b3lyp/6-31++g(2d,2p) nosymm geom=connectivity

A3

1 1
C                 -0.97218474   -0.04841131    0.00000164
C                  0.74654581    0.02078778    0.00004819
N                  0.00908182    1.08732390    0.00003223
N                 -0.08521263   -1.08127384   -0.00003048
N                 -2.25808564    0.02767855   -0.00010128
H                 -2.82238882   -0.81322348    0.00033844
H                 -2.71249356    0.93088470   -0.00019287

1 3 2.0 4 1.0 5 1.0
2 3 1.0 4 2.0
3
4
5 6 1.0 7 1.0
6
7

这样优化始终优化不出我想要的如opt-model.jpg图结构,不自恰,麻烦哪位高手用高斯03或09帮我找到正确优化的原因,如果有时间,还请优化出来贴上结果,谢谢!

opt-model.JPG
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

gauss98

禁虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
renren123123: 回帖置顶 2012-11-27 23:12:43
cenwanglai: 金币+5, QC强帖+1, 感谢热心回帖~ 2012-11-28 10:21:31
你原来的结构两个碳原子离得太近了。拉远一点形成一个接近正方形就可以了

我给你调整的初始坐标输入文件:
% nproc=8
# opt freq b3lyp/6-31++g(2d,2p) nosymm geom=connectivity

A3

1 1
C                 -1.13450913   -0.01127887   -0.00000220
C                  0.95631432    0.02923363    0.00005361
N                 -0.03836219    1.04174624    0.00003058
N                 -0.08521300   -1.08127400   -0.00003000
N                 -2.42041013    0.06481113   -0.00010520
H                 -2.98471313   -0.77609087    0.00033380
H                 -2.87481813    0.96801713   -0.00019720

1 3 1.0 4 1.5 5 2.0
2 3 2.0 4 1.5
3
4
5 6 1.0 7 1.0
6
7



优化结束后的结果算频率的坐标:
C,0,-1.0994478143,0.0776653719,-0.1739757193
C,0,0.7433080732,-0.0902050705,0.5030218806
N,0,-0.0435391597,0.9003683894,-0.3201383759
N,0,-0.0614862126,-0.8627745918,-0.2995524251
N,0,-2.3634015779,0.0714470718,0.0499234862
H,0,-2.8801374252,-0.8041990431,0.0808993623
H,0,-2.8770072735,0.9428622623,0.1599051813


注意结果是个局部能量极小点,能量很高,比线性构型能量高很多!

help1.JPG

3楼2012-11-27 11:02:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)

引用回帖:
8楼: Originally posted by renren123123 at 2012-11-29 08:43:53
基组不是问题,是应该能计算出来的,而且也不会因基组大而使得计算结果不合理!谢谢你的关注!

上来就上这么大的基组很可能造成不收敛或者多花时间,坚持用这么大的基组的话建议先试个小基组再试大基组,比如先上6-31g甚至3-21g,可以省些时间

早期MP2发展中有一个很经典的过大基组导致错误的例子,基组的作用首先得搞清楚。如果你用6-31+g(d,p)算出来的结构和你现在的基组算出来的结构误差在0.01A一下那你为什么还要用那么强的基组?
秀逗魔导师
9楼2012-11-29 10:58:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

lxhlxh052c

木虫 (著名写手)

小木虫之无虫

楼主,我觉得你想要的优化结果像比目鱼啊~
很有艺术细菌啊
这就是我想要的
2楼2012-11-27 10:47:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renren123123

金虫 (正式写手)

引用回帖:
3楼: Originally posted by gauss98 at 2012-11-27 11:02:53
你原来的结构两个碳原子离得太近了。拉远一点形成一个接近正方形就可以了

我给你调整的初始坐标输入文件:
% nproc=8
# opt freq b3lyp/6-31++g(2d,2p) nosymm geom=connectivity

A3

1 1
C            ...

谢谢楼上gauss98,我按你提供的初始结构优化,还是找不到最终稳定的构型,而且该离子结构完全不在一个平面内,完全变形了,就跟我之前用DFT计算的那样,按照有机化学相关知识,起码这个四元环由于有两个双键应该是在同一平面内的,但末优化(但末收敛)的结构不在同一平面内!麻烦您提供指导!谢谢!感激不尽~
4楼2012-11-27 23:10:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renren123123

金虫 (正式写手)

谢谢楼上!只是我按你上说的方法做,还是一样很难最终收敛,且四元环完全不在同一平面内,且变形严重!为什么会这样?
5楼2012-11-27 23:11:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renren123123

金虫 (正式写手)

我计算了该结构,得到其无虚频的构型,如下图所示,但完全变形,这种结构合理吗?

wei.JPG

6楼2012-11-27 23:36:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)

引用回帖:
4楼: Originally posted by renren123123 at 2012-11-27 23:10:08
谢谢楼上gauss98,我按你提供的初始结构优化,还是找不到最终稳定的构型,而且该离子结构完全不在一个平面内,完全变形了,就跟我之前用DFT计算的那样,按照有机化学相关知识,起码这个四元环由于有两个双键应该是 ...

个人观点,该结构角度位阻太大因而采取了扭曲构型。而且键级不是一个绝对的概念,由于四元环及氮原子自身特性,电子的delocalize体现的很明显,你看看底下自己拿到的构型是不是也是有新的CN双键形成?这就是计算和经典理论不一样的地方,事实上很多经典理论,尤其是经典价态和键级理论在当前研究中都已经模糊化了。
ps:你设定的基组太大,这可能也是难以收敛的原因,先考虑一下需不需要这么强的基组
秀逗魔导师
7楼2012-11-28 22:49:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renren123123

金虫 (正式写手)

基组不是问题,是应该能计算出来的,而且也不会因基组大而使得计算结果不合理!谢谢你的关注!
8楼2012-11-29 08:43:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyin1977

金虫 (小有名气)

支持ls的。我也不太明白为什么LZ要用这么大的基组来优化结构。
10楼2012-11-29 13:22:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 renren123123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +42 AI8RaGaPaSCR 2026-08-07 44/2200 2026-08-14 00:45 by J0NtvlhbsrfV
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +11 医学老男孩 2026-08-13 29/1450 2026-08-14 00:04 by 大快人心啦
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
信息提示
请填处理意见