24小时热门版块排行榜    

查看: 3389  |  回复: 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 fatpig8832 at 2011-07-18 17:20:13:
还有一种可能性,这东西说不定是三重态,不知有没有算过,看看和单重态哪个才是ground state

几何优化未完成,计算就结束了。我把最后几行出错的粘贴如下:
Petite list used in FoFCou.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Integral accuracy reduced to 1.0D-05 until final iterations.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UHF) =  -1289.08035487     A.U. after  129 cycles
             Convg  =    0.2150D-04             -V/T =  2.0010
= 0.0000 = 0.0000 = 1.0000 = 2.1029 S= 1.0339
= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation     2.1029,   after     2.0035
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/kss/g09/l502.exe at Sat Jul 23 17:29:37 2011.
33楼2011-07-23 18:10:32
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考研] 272求调剂 +3 材紫有化 2026-02-28 3/150 2026-02-28 22:52 by ms629
[考研] 化工专硕348,一志愿985求调剂 +4 弗格个 2026-02-28 6/300 2026-02-28 22:00 by wang_dand
[考研] 295求调剂 +5 19171856320 2026-02-28 5/250 2026-02-28 21:39 by gaoxiaoniuma
[考博] 26申博 +4 想申博! 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 材料学调剂 +5 提神豆沙包 2026-02-28 5/250 2026-02-28 21:34 by gaoxiaoniuma
[考研] 311求调剂 +8 南迦720 2026-02-28 8/400 2026-02-28 21:30 by gaoxiaoniuma
[考研] 材料类求调剂 +6 wana_kiko 2026-02-28 6/300 2026-02-28 21:20 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 284求调剂 +4 天下熯 2026-02-28 4/200 2026-02-28 21:13 by gaoxiaoniuma
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 7/350 2026-02-28 20:42 by 好好好1233
[教师之家] 版面费该交吗 +15 苹果在哪里 2026-02-22 18/900 2026-02-28 18:20 by mibaomingg
[考研] 285求调剂 +5 满头大汗的学生 2026-02-28 5/250 2026-02-28 18:10 by 材料专硕调剂;
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见