| 查看: 1580 | 回复: 6 | ||
[求助]
求教各位大神，关于今天看的一篇文献中，若干个关于计算参数设定的问题……
|
|
今天看了一篇石墨烯吸附K原子的文献，对计算方法一节里面的若干问题有所疑问，在这里开贴向各位大侠请教 问题1： 原文：The Brillouin zone was described using a Monkhorst-Pack scheme of special k-points to obtain a general convergence criteria of 5 × 10-3 eV for energies and 0.01 Å for structural variations between two different k-point grids. 这里0.01 Å是指POTIM的数值吗？ 问题2： 原文：Band diagrams and DOS analysis were obtained by fixing the Wigner-Seitz radius (rWIGS) for the support during integration over the number of electrons and then by setting rWIGS for the adsorbates within the radii of tangential spheres. This method allowed the accurate assignment of relevant atomic orbital attributions to a particular projected DOS peak. 这段提到的是RWIGS的设置，对于石墨烯吸附K，我的设置是：RWIGS=1.630 2.275(上述两个数值分别是C与K的RWIGS)，不知道这样设置对不对？文中的话没读懂是什么意思…… 问题3：The number of k-points used to describe the irreducible portion of the Brillouin zone was 12 (M&P mesh of 9 × 9 × 1). 这段话就实在是困扰我很久了，为什么9X9X1产生了12个点？？？怎么算的？？？ 请各位大侠多多指教 |
回复此楼» 猜你喜欢
面上提前没消息，有中的吗
已经有17人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有16人回复
-
一个有机合成实验室都需要哪些设备？
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教各位大神，别人文章的参考文献没有写标题，怎么找
已经有23人回复
- 关于投稿时作者简介的问题，请教各位大神。
已经有5人回复
- 【求助】今天做了次SBA-15，失败程度令人发指啊，求教各位
已经有27人回复
 |
2楼2013-10-25 19:05:20
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
奔跑的爷们: 金币+15, ★★★很有帮助 2013-10-25 23:52:35
liliangfang: 金币+3, 谢谢交流 2013-10-27 15:02:48
感谢参与，应助指数 +1
奔跑的爷们: 金币+15, ★★★很有帮助 2013-10-25 23:52:35
liliangfang: 金币+3, 谢谢交流 2013-10-27 15:02:48
|
1. 不是potim 就是两组不同的k点 分别做结构优化后 每个原子的位置变化都小于0.01 Angs 2. 他的意思应该是把C的rwigs固定不动，变K的rwigs. 作用可能就是说K对pdos的贡献可以算的更精确. 应该是这个意思，我没做过不好瞎说。。。离子的半径和原子差很多，所以rwigs是有一定任意性的。楼主可以试试bader analysis。这种方法可以比较“有道理”地决定离子的体积。就不用像文章里那样一遍遍的试了。代码在graeme henkleman的网页上有。用起来很方便。 3. 在outcar里找nkpts看等于几就知道了 |
3楼2013-10-25 21:50:45
4楼2013-10-25 23:52:26
5楼2013-10-26 00:24:41
6楼2013-10-26 09:33:47
7楼2013-10-26 11:07:02
50