| 查看: 1602 | 回复: 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格式的文件）并分别传上来了。能否请大家帮忙看看，是不是我的操作有问题呢？多谢指教！ |
回复此楼» 收录本帖的淘帖专辑推荐
 模拟 |
» 猜你喜欢
化工学硕294分，求导师收留
已经有37人回复
-
求调剂
已经有11人回复
-
260求调剂
已经有4人回复
-
一志愿华中农业071010，320求调剂
已经有19人回复
-
304求调剂
已经有7人回复
-
求博导｜生物质基多孔碳/超级电容方向，已有相关成果，寻能源材料/碳材料方向老师
已经有3人回复
-
二苯甲酮酸类衍生物
已经有6人回复
-
接受任何调剂
已经有4人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有12人回复
-
320求调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
fo3mt2
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1329.1
- 散金: 54
- 红花: 4
- 帖子: 129
- 在线: 158.7小时
- 虫号: 809980
- 注册: 2009-07-16
- 专业: 凝聚态物性 II ：电子结构
2楼2009-07-22 10:05:05
fo3mt2
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1329.1
- 散金: 54
- 红花: 4
- 帖子: 129
- 在线: 158.7小时
- 虫号: 809980
- 注册: 2009-07-16
- 专业: 凝聚态物性 II ：电子结构
3楼2009-07-22 10:05:45
fo3mt2
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1329.1
- 散金: 54
- 红花: 4
- 帖子: 129
- 在线: 158.7小时
- 虫号: 809980
- 注册: 2009-07-16
- 专业: 凝聚态物性 II ：电子结构
4楼2009-07-22 10:06:16
fuzp
铁杆木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 5745.3
- 散金: 139
- 红花: 12
- 帖子: 555
- 在线: 72.8小时
- 虫号: 422943
- 注册: 2007-07-19
- 性别: GG
- 专业: 无机非金属类光电信息与功
5楼2009-07-22 21:31:10
fo3mt2
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1329.1
- 散金: 54
- 红花: 4
- 帖子: 129
- 在线: 158.7小时
- 虫号: 809980
- 注册: 2009-07-16
- 专业: 凝聚态物性 II ：电子结构
6楼2009-07-23 04:38:12
crante
木虫 (正式写手)
登高
- 应助: 33 (小学生)
- 金币: 4032.2
- 散金: 129
- 红花: 1
- 帖子: 665
- 在线: 164.8小时
- 虫号: 218001
- 注册: 2006-03-12
- 性别: GG
- 专业: 半导体材料
7楼2009-07-23 12:06:56
crante
木虫 (正式写手)
登高
- 应助: 33 (小学生)
- 金币: 4032.2
- 散金: 129
- 红花: 1
- 帖子: 665
- 在线: 164.8小时
- 虫号: 218001
- 注册: 2006-03-12
- 性别: GG
- 专业: 半导体材料
8楼2009-07-23 12:09:32
fo3mt2
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1329.1
- 散金: 54
- 红花: 4
- 帖子: 129
- 在线: 158.7小时
- 虫号: 809980
- 注册: 2009-07-16
- 专业: 凝聚态物性 II ：电子结构
|
谢谢你的回复，我确实试过的你说的把氧原子的U参数值改为零。操作过程如下，选中钛原子，把它的U参数值改成5.00，然后选中氧原子，把它的U参数值改成0.00，然后其它设置都没变，保存计算任务之后，再重新选中被保存任务中的Rutile.xsd文件中的氧原子，在察看Electronic configuration,可是显示的却是5.00，我真的不知道该怎么做了。。。真的求各位朋友帮忙啦。。。 至于Rutile.cell文件我是上传到服务器上之后，在linux系统下打开我保存并上传上来的任务时看到的，你可以看我上传上来的附件里，有这个文件的内容。。。 还是在此求大家帮忙啦，万分感谢！ |
9楼2009-07-23 14:34:06
★ ★ ★ ★
小木虫(金币+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
50