24小时热门版块排行榜    

查看: 2507  |  回复: 12

cmyy1988

金虫 (小有名气)

[交流] 【求助】重金求解:如何用castep算小分子裂解的过渡态 已有5人参与

如题,我要用castep算羟基裂解为氢原子与氧原子的过渡态,先在一个空晶胞里放置OH进行优化作为反应物,产物的处理上,我是在相同大小的晶胞里放置氢原子与氧原子,二者距离为2.5A,但是不管这个距离怎么设置,计算总是失败,我把圈数设成1000圈依然失败,显示距离为0.785A,以下为每次都出现的错误提示:
Warning: electronic minimisation did not converge when finding ground state.
  ****************************************************************************

Writing model to 3D_Atomistic__10_.check                                                                                                                                                                                                                                         
Error in geom_get_forces - electronic_minimisation of current_cell failed
Error in geom_get_forces - electronic_minimisation of current_cell failed
Error in geom_get_forces - electronic_minimisation of current_cell failed
Error in geom_get_forces - electronic_minimisation of current_cell failed
Error in geom_get_forces - electronic_minimisation of current_cell failed
Error in geom_get_forces - electronic_minimisation of current_cell failed
Error in geom_get_forces - electronic_minimisation of current_cell failed
Error in geom_get_forces - electronic_minimisation of current_cell failed
MPI Application rank 3 exited before MPI_Finalize() with status 1
请问大伙,该怎么设置产物氢原子与氧原子的构型,怎么去优化,金币是小,不胜感激是真,谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yb2006

铁杆木虫 (正式写手)

★
cmyy1988(金币+1):谢谢参与
cmyy1988(金币+5): 2010-05-12 15:14:36
不知道你这个羟基是在什么条件下裂解的,是在催化剂的催化下呢还是热裂解?
2楼2010-05-11 13:26:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)

★
cmyy1988(金币+1):谢谢参与
cmyy1988(金币+5): 2010-05-12 15:14:50
模型不合理
3楼2010-05-11 20:52:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

franch

版主 (著名写手)

★
cmyy1988(金币+1):谢谢参与
cmyy1988(金币+5): 2010-05-12 15:15:03
引用回帖:
Originally posted by cmyy1988 at 2010-05-11 09:13:25:
如题,我要用castep算羟基裂解为氢原子与氧原子的过渡态,先在一个空晶胞里放置OH进行优化作为反应物,产物的处理上,我是在相同大小的晶胞里放置氢原子与氧原子,二者距离为2.5A,但是不管这个距离怎么设置,计算 ...

空晶胞里放置OH进行优化作为反应物,,,


这好像不行吧,, 应该放到你产生作用的晶胞上面
遇弱则弱
4楼2010-05-12 14:40:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cmyy1988

金虫 (小有名气)

引用回帖:
Originally posted by franch at 2010-05-12 14:40:12:






空晶胞里放置OH进行优化作为反应物,,,


这好像不行吧,, 应该放到你产生作用的晶胞上面

是这样,我是做了羟基在金属表面上的裂解,得到了一个能垒,老师说为了比较这个金属的催化性能,必须要知道没有催化剂下羟基裂解的能垒大小,所以我得知道在没有催化剂下羟基裂解的能量,希望您能指点一二
5楼2010-05-12 15:17:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cmyy1988

金虫 (小有名气)

引用回帖:
Originally posted by franch at 2010-05-12 14:40:12:






空晶胞里放置OH进行优化作为反应物,,,


这好像不行吧,, 应该放到你产生作用的晶胞上面

是这样,我是做了羟基在金属表面上的裂解,得到了一个能垒,老师说为了比较这个金属的催化性能,必须要知道没有催化剂下羟基裂解的能垒大小,所以我得知道在没有催化剂下羟基裂解的能量,希望您能指点一二
6楼2010-05-12 15:17:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cmyy1988

金虫 (小有名气)

引用回帖:
Originally posted by jghe at 2010-05-11 20:52:10:
模型不合理

是这样,我是做了羟基在金属表面上的裂解,得到了一个能垒,老师说为了比较这个金属的催化性能,必须要知道没有催化剂下羟基裂解的能垒大小,所以我得知道在没有催化剂下羟基裂解的能量,因此我在相同的计算方法下设置了那么一个空晶胞,想分别优化作为反应物的羟基和作为产物的氧原子与氢原子,在找过渡态。希望您能指点一二
7楼2010-05-12 15:19:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yb2006

铁杆木虫 (正式写手)

★ ★
zzy870720z(金币+2):感谢交流 2010-05-13 08:35:22
cmyy1988(金币+5): 2010-05-13 09:25:56
我觉得没有必要去找没有催化剂条件下的能垒,因为催化剂的作用就是降低反应的能垒的,有了催化剂这个能垒肯定会比较低。。。如果一定要求的话,嘿嘿,所谓的能垒不就是氢和氧成键的键能吗,你把H+O能量减去OH能量不就行了?还有,我还怀疑没有催化剂情况下这反应有没有过渡态呢,呵呵。顺带问一下,你算出来在催化剂上的能垒是多少?
8楼2010-05-12 15:38:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

franch

版主 (著名写手)

引用回帖:
Originally posted by yb2006 at 2010-05-12 15:38:12:
我觉得没有必要去找没有催化剂条件下的能垒,因为催化剂的作用就是降低反应的能垒的,有了催化剂这个能垒肯定会比较低。。。如果一定要求的话,嘿嘿,所谓的能垒不就是氢和氧成键的键能吗,你把H+O能量减去OH能量 ...

正解
遇弱则弱
9楼2010-05-12 21:46:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cmyy1988

金虫 (小有名气)

[quote]Originally posted by yb2006 at 2010-05-12 15:38:12:
我觉得没有必要去找没有催化剂条件下的能垒,因为催化剂的作用就是降低反应的能垒的,有了催化剂这个能垒肯定会比较低。。。如果一定要求的话,嘿嘿,所谓的能垒不就是氢和氧成键的键能吗,你把H+O能量减去OH能量 ... [/quo
用H和O的能量和减去OH的能量,这应该是一个反应函变,是一条曲线头和尾的差距,而并不是头到过渡态的差距
10楼2010-05-13 11:03:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cmyy1988 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
请填处理意见