24小时热门版块排行榜    

查看: 3466  |  回复: 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的回帖

fatpig8832

铁杆木虫 (著名写手)

★ ★
zhou2009(金币+2): 2011-08-06 15:19:21
引用回帖:
Originally posted by jianqin at 2011-07-25 08:42:08:
谢谢!虽然我不知道这个E(uhf)=-1289.08035487代表的物理意义,但是这个值应该比单重态的小。
单重态,rhf 6-31g能量结果:
Zero-point correction=                           0.358864 (Hartree/Particle) ...

这个E(UHF)是指electronic energy,每优化一步,就会给出一个E(UHF)。对大部分寻找minimum的优化来说,每一步的E(UHF)总体趋势上是一个降低的过程,

在你的单重态输出中,electronic energy应该是
=-1288.405788-0.358864
=-1288.385237-0.379415
=...
=...
=-1288.764652

你可以再搜一下三重态的输出文件的其他E(UHF),如果都在-1289.08035487的附近,都比单重态能量低的话,那么基本可以断定,它的ground state是三重态(也可以试试五重,我想机会不大)。然后,你现在做的单重态的这些东西就没有太大意义了...

» 本帖已获得的红花(最新10朵)

38楼2011-07-25 13:15:03
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工科 320总分 求调剂 +7 梨花珞晚风 2026-03-17 7/350 2026-03-22 10:03 by 2026paper
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 324求调剂 +5 lucky呀呀呀鸭 2026-03-20 5/250 2026-03-20 22:30 by 促天成
[考研] 一志愿华中农业071010,总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 281求调剂(0805) +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
信息提示
请填处理意见