24小时热门版块排行榜    

查看: 1237  |  回复: 9
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

majun04

木虫 (著名写手)

majun

[交流] 【求助】计算过渡态,优化初始反应物遇到困难请帮忙!

想计算一个质子氢参与的化学反应。第一个定点应该是质子氢和反应物放在一起优化结构。因为是能量图中的初始定点,要求氢离子和反应分子没有相互作用,所以放置位置设置为3埃,发现优化以后,氢离子和我的分子中带有负电荷的原子成键了(这重情况在计算带有电荷的离子解离开的时候也发生过,依然没有解决)。调节距离后(6埃),优化结构依然错误,输出文件提示如下:
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB+HF-LYP) =  -132.764483473     A.U. after  129 cycles
             Convg  =    0.1330D-01             -V/T =  2.0229
             S**2   =   0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/majun/gaussian03/g03/l502.exe at Wed Nov 19 00:45:40 2008.
Job cpu time:  0 days  0 hours  2 minutes  2.6 seconds.
File lengths (MBytes):  RWF=     11 Int=      0 D2E=      0 Chk=      4 Scr=      1

因为是第一次计算带有电荷的原子参与的反应过渡态,所以很困惑。请高手帮忙指点一下。
还有就是在计算反应过程中,对A+B-->C+D的反应,感觉初始能量计算时候要优化A+B的时候(把两个分子放到一个输入文件中),比较困难。其中一个部分是带有电荷的(比如说金属原子),电荷要固定一下吗(也就是指定电荷都在这个金属原子上)?
初始的能量是不是也可以用A和B的单独的能量计算加和?我要是想A+B放在一起,为后面用QST2粗略的找过渡态方便。
回复此楼
乐观、执着、勤奋
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

引用回帖:
Originally posted by dolphin618 at 2008-11-19 03:47:
3埃,感觉距离小了点,其实10埃我觉得可以,优化的时候加上scf=qc试试

依然不行,错误信息是一样的。请帮忙想想还有没有别的办法?
乐观、执着、勤奋
6楼2008-11-19 16:02:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

dolphin618

金虫 (正式写手)

★ ★ ★ ★
majun04(金币+2,VIP+0):谢谢你,我加上试试。
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
3埃,感觉距离小了点,其实10埃我觉得可以,优化的时候加上scf=qc试试
生命本就是不可名说的缘法,可又有谁愿做菩提下微笑拈起的莲花?
2楼2008-11-19 03:47:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cysqxr

木虫 (小有名气)

★ ★ ★
majun04(金币+1,VIP+0):谢谢你!
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
从你的描述看应该是有关于吸附的问题(我猜的不对就算了啊)具体的吸附距离是多少要根据实际情况判断  另外你说的电荷 我知道的范围内好像不可以单独固定把电荷加在某个原子上 如果加的话 就是整体加
3楼2008-11-19 08:24:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dongyuhui

银虫 (正式写手)

★ ★ ★
majun04(金币+1,VIP+0):谢谢你,我加大距离看看怎么样。
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
同意楼上的说法,楼主可以试一下把距离变大,我以前改了距离就好了
快乐的心
4楼2008-11-19 08:51:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-13 6/300 2026-09-17 23:28 by cuPZTDXS3VOd
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-13 9/450 2026-09-17 23:05 by cuPZTDXS3VOd
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-17 16:03 by 0ngsZVSCTS7z
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 14:27 by ObVzqQQh1rrE
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 40ms4Wlo4umh 2026-09-14 5/250 2026-09-17 14:05 by ObVzqQQh1rrE
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-17 13:29 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 yaCh4X3wd045 2026-09-14 4/200 2026-09-17 13:17 by ObVzqQQh1rrE
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +8 3n8v2C8RimXI 2026-09-13 8/400 2026-09-17 11:27 by UAIjTRMuUFoj
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-17 05:50 by k06BKNblrGRH
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-17 05:14 by k06BKNblrGRH
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-17 03:14 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-17 03:04 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-13 6/300 2026-09-17 02:14 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 02:23 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见