24小时热门版块排行榜    

查看: 1521  |  回复: 16

Madlife

铁虫 (初入文坛)

[求助] 分子吸附问题

我用DMOL3做111面吸附CO2的时候,HCP、FCC、桥位、顶位的竖直、水平吸附优化后得到的全是物理吸附,但是我们老师说有化学吸附,我查了相关文献说我这个面如果不加氢活化是没有CO2化学吸附的,可是我们老师又说不用加氢。那我想问,如果要加氢的话,那我之前优化的这些是不是就要全部重新加氢再优化一次,还有就是加氢的话,该怎么操作呢?
另外,我们这个表面是个催化剂表面,并不洁净,有微量的Zr,Zn等元素,排列顺序不明,做不出化学吸附和没考虑这个有关系吗
希望有前辈能指点迷津,谢谢了!
回复此楼

» 收录本帖的淘帖专辑推荐

量子化学和计算化学

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

心之所向,无所畏惧。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur


identation: 金币+1, 谢谢交流 2013-11-20 17:44:41
不知道你怎么定义物理吸附和化学吸附的分界。
你的表面和计算参数似乎也没给出来,就没办法判断是不是在建模或是计算的时候出错的。

我对催化不是特别懂,更详细的信息也许会对后面来回答、诊断的人有帮助。
不过虽然有微量的杂质金属,如果计算足够合理也许能够找到现象的主要矛盾。
不一定挂在论坛,计算问题问题欢迎留言。
3楼2013-11-19 00:59:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gump_813

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
identation: 金币+1, 谢谢交流 2013-11-20 17:44:49
不是很清楚 但是CO2的吸附很弱是肯定的 至于实验上为什么有化学吸附 那就不好说了
你是什么催化剂?
4楼2013-11-19 01:14:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
identation: 金币+2, 谢谢交流 2013-11-20 17:44:58
引用回帖:
7楼: Originally posted by Madlife at 2013-11-19 19:24:11
另外  前辈有用过dmol中的Dynamics么
我问了做实验的学姐,她说吸附是在300+度下进行的...

对于12核心服务器(2.6GHz)你这个体系正常应该大约半分钟一次scf,一般30次scf计算一个构型,即15分钟。动力学模拟我不懂,但怎么着也要取1000个构型,显然不现实。而且有很多东西我还是按照高估的。

smearing=0.005Ha其实已经相当于电子受到温度的作用排布有所变化(smearing=kT,k是Boltzmann常量)。你这个体系其实没那么多复杂能量又接近的构象,优化就能解决问题。先把上述准备工作做好,优化做好,考虑一些杂质、缺陷的问题,以及考虑反应路径、过渡态,这些可能对催化比较重要。

是些个人理解,望大家多指教吧。
不一定挂在论坛,计算问题问题欢迎留言。
9楼2013-11-19 22:59:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

Madlife

铁虫 (初入文坛)

我又看了一遍文献我觉得是我理解错误了,不管那篇文献,做出的全是物理吸附正常吗?
求前辈指点啊,谢谢了!
心之所向,无所畏惧。
2楼2013-11-18 22:07:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Madlife

铁虫 (初入文坛)

引用回帖:
4楼: Originally posted by gump_813 at 2013-11-19 01:14:39
不是很清楚 但是CO2的吸附很弱是肯定的 至于实验上为什么有化学吸附 那就不好说了
你是什么催化剂?

铜,文献说80K时才会产生CO2 2-离子,之后才会发生化学吸附,dmol3是研究0K时的情况吧,我不知道该怎么设置温度了
心之所向,无所畏惧。
5楼2013-11-19 19:04:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Madlife

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by 卡开发发 at 2013-11-19 00:59:12
不知道你怎么定义物理吸附和化学吸附的分界。
你的表面和计算参数似乎也没给出来,就没办法判断是不是在建模或是计算的时候出错的。

我对催化不是特别懂,更详细的信息也许会对后面来回答、诊断的人有帮助。
不 ...

感谢前辈又一次来指导我了
我是根据优化后计算的吸附能来定义的,都小于10KJ/mol。
至于我的计算参数,优化质量采用高质量的默认设置,smearing为0.005Ha,计算都得到了收敛。
分子吸附问题
图1.png

心之所向,无所畏惧。
6楼2013-11-19 19:14:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Madlife

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by 卡开发发 at 2013-11-19 00:59:12
不知道你怎么定义物理吸附和化学吸附的分界。
你的表面和计算参数似乎也没给出来,就没办法判断是不是在建模或是计算的时候出错的。

我对催化不是特别懂,更详细的信息也许会对后面来回答、诊断的人有帮助。
不 ...

另外  前辈有用过dmol中的Dynamics么
我问了做实验的学姐,她说吸附是在300+度下进行的
心之所向,无所畏惧。
7楼2013-11-19 19:24:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
identation: 金币+3, 谢谢交流 2013-11-20 17:45:10
引用回帖:
6楼: Originally posted by Madlife at 2013-11-19 19:14:16
感谢前辈又一次来指导我了
我是根据优化后计算的吸附能来定义的,都小于10KJ/mol。
至于我的计算参数,优化质量采用高质量的默认设置,smearing为0.005Ha,计算都得到了收敛。

图1.png
...

唉,我也算不上什么前辈。

严格来说,对于单个分子的吸附,dmol3处理这个体系应该把晶面取到两倍的orbital cutoff与分子的xy方向长度,这样才会保证临近的分子间的作用对吸附没影响。当然现在很多吸附做的都不是很规范。

明显的问题还有金属层数原则上也要测试的,这个金属原子实在太薄了离真实体系有距离,你可以测试随着金属厚度增加计算结果会有所不同,这个误差有时候不小。

看不出的一些问题你也要注意,对于dmol3,k-point和orbital cutoff也需要测试,有的时候默认的参数满足不了测试的结果。另外计算吸附能的时候,分子的计算不能简单选择为Fine,因为分子的碳氧默认的orbital cutoff和金属的不一样,对于Cu是4.4埃,就要都设定成4.4埃。当然也有其他的方案。 否则做差的能量精度不同。

这个能量有的时候未必可靠。10kJ/mol的吸附能有可能不超过你计算误差的值。除了模型大小和上述的参数,因为不同的原子采用的基组大小有差异,还要考虑BSSE校正。另外,计算方法的选用也很重要,对于Dmol3有些方法有的时候偏差5kcal/mol也不稀奇,如果确定是物理吸附,那么色散力的校正也要考虑考虑。恩。。。Dmol3有时候处理金属体系很不好。

有的时候除了那些杂质元素,有的时候金属本身的缺陷也不容忽视。有些奇特的性质不是完美的结构表现得出来的。

都是些个人观点,希望再多一些大牛指正讨论。有啥问题可以再讨论一下。
不一定挂在论坛,计算问题问题欢迎留言。
8楼2013-11-19 22:43:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Madlife

铁虫 (初入文坛)

引用回帖:
9楼: Originally posted by 卡开发发 at 2013-11-19 22:59:50
对于12核心服务器(2.6GHz)你这个体系正常应该大约半分钟一次scf,一般30次scf计算一个构型,即15分钟。动力学模拟我不懂,但怎么着也要取1000个构型,显然不现实。而且有很多东西我还是按照高估的。

smearing ...

明白了,非常感谢!
心之所向,无所畏惧。
10楼2013-11-19 23:56:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Madlife 的主题更新
信息提示
请填处理意见