| 查看: 3186 | 回复: 17 | ||||
[求助]
求助：用CASTEP计算O2在Au(111)面的吸附！ 已有1人参与
|
||||
|
各位：大家好！ 如题，本人用CASTEP计算O2在Au(111)面的吸附能。 O2分子及Au(111)面已经会建了。谁能帮忙介绍一下，下面的操作流程吗？ 疑点三： 1. MS CASTEP help文件中关于CO吸附在Pd(110)面上的例子是在已经知道了C-Pd之间形成化学键（已知其键长）的基础上做的，而现在我不知道我的O2分子与Au的相互作用，这样怎么计算吸附能？ 2. 怎么把O2放到Au(111)面上，直接copy进去吗？关键是这样放进去后怎么摆放O2分子，离Au(111)面多高，还是随意都行？ 3. 把O2分子放到Au(111)面后，是直接进行几何优化吗？不用进行其他的运算吗？（比如动力学，因为用forcite模块进行吸附能计算的时候，弛豫过程是要进行动力学运算的，不懂CASTEP模块怎么样）。 先谢谢各位了！ |
回复此楼» 猜你喜欢
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有10人回复
-
基金系统什么内容也没有
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3567.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
|
随便怎么画，所谓第一原理当然是不包含这些成键参数的，要不然就不叫第一原理了。另一个问题其实你给得信息很笼统，计算失败的原因有很多，比如SCF难收敛或者构型优化难收敛，这样的情况最好是把输出文件贴过来，大家才容易帮忙诊断。姑且按照以下因素先做排查： 1、首先所有的动能截断和k点最好都要经过合适的测试，如果参数不合理算出来的结果基本上意义不大，MS默认的参数并不是100%就可靠的。 2、放的比较远一般是物理吸附。如果是物理吸附的话，O2可能会维持气态状态的三重态，三重态计算必须开自旋极化，或者有些程序叫非限制自旋，你可以在modify-electronic configuration去给这两个原子各+1的spin up，这样在SCF过程中就不至于掉到其他的自旋态了；如果要做解离吸附，有可能就得多试试了。如果你放的实在太远就偏离实际平衡结构太多，程序需要相当长的优化步才能得到合适的结构，这样的做法当然也就不合适了。 3、另外还可能是这样的，物理吸附一般色散作用可能占不小的比例，而一般的LDA/GGA甚至杂化泛函都很难正确描述这样的色散作用，因此需要修正，DMol3或CASTEP可选择的方案并不多，可以考虑DFT-D2（DFT-D的Grimme）或者M06L，前者修正的部分计算量是力场级别的，但后者计算量要大很多。但如果不进行修正可能本该物理吸附的分子就会被算成推斥状态，这样分子越飞越远自然也就很难计算成功。 4、如果打开自旋之后SCF难收敛可以去降低density mixing部分的charge。 最后，表面吸附而言，DMol3的计算速度要比CASTEP快得多，精度虽然略有损失但基本上在可接受的范围内，虽然各自有各自的缺点，但还是值得去尝试一下。 |
11楼2018-05-09 01:22:27
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3567.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
17楼2018-05-09 20:30:41
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
lidegui354: 金币+90 2018-05-08 17:55:17
感谢参与，应助指数 +1
lidegui354: 金币+90 2018-05-08 17:55:17
|
各位：大家好！ 如题，本人用CASTEP计算O2在Au(111)面的吸附能。 你確定嗎？ O2分子及Au(111)面已经会建了。谁能帮忙介绍一下，下面的操作流程吗？ 疑点三： 1. MS CASTEP help文件中关于CO吸附在Pd(110)面上的例子是在已经知道了C-Pd之间形成化学键（已知其键长）的基础上做的，而现在我不知道我的O2分子与Au的相互作用，这样怎么计算吸附能？ 沒關係，把O2放上去就行，自己就馳豫了 但放的位置，比較微妙 2.5-3A的距離吧 2. 怎么把O2放到Au(111)面上，直接copy进去吗？关键是这样放进去后怎么摆放O2分子，离Au(111)面多高，还是随意都行？ 如上 3. 把O2分子放到Au(111)面后，是直接进行几何优化吗？不用进行其他的运算吗？（比如动力学，因为用forcite模块进行吸附能计算的时候，弛豫过程是要进行动力学运算的，不懂CASTEP模块怎么样）。 先優化再說 |
3楼2018-05-08 12:32:02
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3567.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
6楼2018-05-08 19:10:16
卡开发发10楼2018-05-08 23:54:44
2楼2018-05-07 11:12:00
4楼2018-05-08 17:53:19
5楼2018-05-08 19:04:09
7楼2018-05-08 19:51:23
8楼2018-05-08 19:52:04
9楼2018-05-08 20:29:15
10