24小时热门版块排行榜    

Znn3bq.jpeg
查看: 893  |  回复: 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的回帖

cysqxr

木虫 (小有名气)

★ ★ ★ ★ ★
majun04(金币+5,VIP+0):谢谢你的帮助!
我说一下自己的看法
单独的用质子与其他物质作用,这个模型不是很合理。通常要做质子与其他的分子作用,需要把质子溶剂化(例如采用水合氢离子H3O+),然后你做计算就不会出现不收敛的问题了。l502的错误意味着你的结构出现了很不合理的情况,像你这种单独解离质子的,把质子离得太远以后,化学上来说就没有电子了,这在量化中很难处理。即使你使用了scf=qc也不会有所改变。看看你的反应体系,采用相应的溶剂来处理你的质子。
另外呢,如果目标分子已经被质子化了,接着再发生其他反应的情况,可以把溶剂分子去掉了做(如果溶剂分子不再参与反应)。

[ Last edited by cysqxr on 2008-11-28 at 11:00 ]
10楼2008-11-28 10:58:40
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 300求调剂 +4 橙a777 2026-04-15 4/200 2026-04-15 20:05 by 学员JpLReM
[考研] 085404 22408 309分求调剂 +9 lzmk 2026-04-14 10/500 2026-04-15 20:02 by 学员JpLReM
[考研] 327求调剂 +24 Xxjc1107. 2026-04-13 26/1300 2026-04-15 19:27 by coolfishwll
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[考研] 279求调剂 +12 张番茄不炒蛋 2026-04-11 12/600 2026-04-14 15:38 by zs92450
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 300分求调剂 (085501机械专硕,本科扬大) +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu~
[基金申请] 2026 WR青拔 +3 冬日阳光CAS 2026-04-09 6/300 2026-04-13 18:40 by liuchb715
[考研] 一志愿双非085400电子信息344 求调剂,对材料和化学方向也感兴趣 +12 无情的小羊 2026-04-09 13/650 2026-04-13 14:17 by 张zhihao
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
[考研] 一志愿沪9,326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
[考研] 招收有机化学、化工,药学,食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
信息提示
请填处理意见