| 查看: 2788 | 回复: 0 | ||||
394163597铁杆木虫 (著名写手)
|
[交流]
VASP练习：H在Cu(100)面的吸附能计算
|
|
按照《Density Functional Theory： A Practical Introduction》Chapter4 表面能及吸附能这部分，自己练习了H原子在Cu(100)面TOP、BRIDGE、HOLLOW三个位点的吸附能计算，计算结果和教材出入蛮大，想和算过的前辈比较比较结果。 以下是计算简介及结果： 1.H2分子能量计算：BOX大小8×8×8 angstrom，K点选取Gamma点，截断能350eV，优化得到H2分子键长0.75angstrom，总能-6.74eV 2.Cu(100)自由表面能量计算：p(1×1)，5层(100)面，10倍(100)层间距的真空层，固定底部两层、弛豫其余三层，K点Gamma 14×14×1，截断能350eV，开启偶极子矫正，优化得到Cu(100)自由表面总能-17.47eV 3.放置H原子在Cu(100)面上的TOP、BRIDGE、HOLLOW三个位点上，INCAR参数不变，优化得到相应的总能为-20.13eV、-20.77eV、-20.88eV 4.使用公式：吸附能=总能(H在100面)-1/2*总能(H2分子)-总能(100自由表面) 得到三个位点的相应吸附能为+0.71eV、+0.07eV、-0.04eV，与书中p104的计算结果+0.38eV、-0.11eV、-0.19eV的结果相差较大 我相信一部分原因在于参数选取不一致，如K mesh、Ecutoff、EDIFF(G)等等，但是结果相差也太大了吧  恳请前辈予以指点，我觉得我思路没错，也是对着教材和教程认真研究了的  如果也有练习的结果恳请交流交流 |
回复此楼» 收录本帖的淘帖专辑推荐
 Vasp |
» 猜你喜欢
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有188人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有120人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
