24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

fo3mt2

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1329.1
散金: 54
红花: 4
帖子: 129
在线: 158.7小时
虫号: 809980
注册: 2009-07-16
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】关于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格式的文件）并分别传上来了。能否请大家帮忙看看，是不是我的操作有问题呢？多谢指教！

回复此楼

» 收录本帖的淘帖专辑推荐

模拟

» 猜你喜欢

最失望的一年已经有17人回复
拟解决的关键科学问题还要不要写已经有8人回复
为什么nbs上溴没有产物点出现呢已经有6人回复
求推荐博导已经有4人回复
存款400万可以在学校里躺平吗已经有34人回复
求助一下有机合成大神已经有4人回复
求推荐英文EI期刊已经有5人回复
26申博已经有3人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
疑惑？已经有5人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2009-07-22 10:03:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

stou

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 471.6
红花: 2
帖子: 164
在线: 27小时
虫号: 738274
注册: 2009-04-02
性别: GG
专业: 其他无机非金属材料

★ ★ ★ ★
小木虫(金币+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 ]