| 查看: 2094 | 回复: 12 | ||||
| 当前主题已经存档。 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
fo3mt2金虫 (小有名气)
|
[交流]
【求助】关于DFT+U功能的使用问题！
|
|||
|
大家好： 我有一个问题想向大家请教。我正在使用Materials Studio 4.4中的CASTEP工具作半导体材料的计算。我在计算二氧化钛（是Rutile晶型）能带。由于钛原子是4d周期的过度金属元素，我使用了DFT+U的功能来计算。我把我的操作过程简诉如下：首先选中名为Rutile.xsd文件中的一个钛原子，然后点击Modify菜单下的Electronic Configuration，选择U tab，在d电子一栏里我输入U参数值为5.00，在Spin tab下我选了Formal Spin中的Up为+4,Down为-4，然后我关闭了Electronic Configuration。在Modules菜单下，我选了CASTEP工具，在Task一栏里我选了Energy,Functional一栏里我选了GGA-WC，然后在Spin polarized,Use fromal as initial和Use LDA+U三个框里打了钩，将Charge数值设定为0。我将Energy cutoff,SCF tolerance和k-point set都设定成Ultra-fine，Pseudopotentials选择的是Ultrasoft。在Properties一栏里我选了Band structure,Density of states和Electron density difference三个任务，最后将任务保存并上传到服务器上。 可是问题来了，当我打开Rutile.cell文件时，我看到HUBBARD_U BLOCK里不仅钛原子的d电子的U参数值被设定成5.00，而且氧原子的d电子的U参数值也被设定成了5.00。可是氧原子的价电子里是没有d电子的。为什么U校正参数值里会有氧原子的d电子数值（5.00）呢？是不是我的操作过程有问题呢？直接这样执行能量的计算，会不会出错呢？我把Rutile.cell，Rutile.param，Rutile.xsd和Rutile.castep四个文件用wordPad记录了下来（因为小木虫论坛的商船附件功能不允许直接上传.cell，.param，.xsd和.castep格式的文件）并分别传上来了。能否请大家帮忙看看，是不是我的操作有问题呢？多谢指教！ |
回复此楼» 收录本帖的淘帖专辑推荐
 模拟 |
» 猜你喜欢
极限拉扯，如何抉择
已经有15人回复
-
申请科研项目纯靠实力能成功吗？
已经有8人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
freshgirl(金币+1,VIP+0):谢谢积极参与~ 7-24 23:12
wuli8(金币+1,VIP+0):http://emuch.net/bbs/viewthread.php?tid=1354235 7-25 16:29
wuli8(金币+1,VIP+0):3ks,(*^__^*) …… 7-25 16:29
小木虫(金币+0.5):给个红包，谢谢回帖交流
freshgirl(金币+1,VIP+0):谢谢积极参与~ 7-24 23:12
wuli8(金币+1,VIP+0):http://emuch.net/bbs/viewthread.php?tid=1354235 7-25 16:29
wuli8(金币+1,VIP+0):3ks,(*^__^*) …… 7-25 16:29
|
楼主，我按照你的方法试了下，没有出现你说的问题啊。。。你最好重新再建立模型，重新算一下。 另外，我觉得你的自选设置的可能有点问题，你的晶胞自旋是+4，我觉得应该是0。当然，如果你在SCF中选了optimisize spin的话这样应该没影响。 又看了下你的castep文件，感觉结果应该能用 LDA+U occupation numbers ------------------------ Species Ion l spin occupancies ============================== Ti 1 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 2 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 3 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 4 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 5 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 6 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 7 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 8 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 9 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 10 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 11 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 12 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 13 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 14 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 15 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 Ti 16 2 up: 0.23 0.26 0.28 0.51 0.52 dn: 0.23 0.26 0.28 0.51 0.52 ================================== 这里没有说O的事，只有Ti的LDA+U [ Last edited by stou on 2009-7-24 at 10:33 ] |
10楼2009-07-24 10:29:38
fo3mt2
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1332.6
- 散金: 54
- 红花: 4
- 帖子: 129
- 在线: 158.9小时
- 虫号: 809980
- 注册: 2009-07-16
- 专业: 凝聚态物性 II ：电子结构
2楼2009-07-22 10:05:05
fo3mt2
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1332.6
- 散金: 54
- 红花: 4
- 帖子: 129
- 在线: 158.9小时
- 虫号: 809980
- 注册: 2009-07-16
- 专业: 凝聚态物性 II ：电子结构
3楼2009-07-22 10:05:45
fo3mt2
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1332.6
- 散金: 54
- 红花: 4
- 帖子: 129
- 在线: 158.9小时
- 虫号: 809980
- 注册: 2009-07-16
- 专业: 凝聚态物性 II ：电子结构
4楼2009-07-22 10:06:16
200