应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】关于DFT+U功能的使用问题！
131/2返回列表12下一页
查看: 1489  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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格式的文件）并分别传上来了。能否请大家帮忙看看，是不是我的操作有问题呢？多谢指教！
回复此楼

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

模拟

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2009-07-22 10:03:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fo3mt2

金虫 (小有名气)
这个是Rutile param文件。
赞一下 回复此楼
2楼2009-07-22 10:05:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fo3mt2

金虫 (小有名气)
这个是Rutile castep文件。
赞一下 回复此楼
3楼2009-07-22 10:05:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fo3mt2

金虫 (小有名气)
这个是Rutile xsd文件。
赞一下 回复此楼
4楼2009-07-22 10:06:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fuzp

铁杆木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuchenwf(金币+2,VIP+0):谢谢回帖交流 7-24 00:36
你在Rutile.xsd文件中选氧原子，在Electronic Configuration中看看它的U值，
.cell里应该没有氧的U值的，也许你有没提到的误操作。
赞一下(11人) 回复此楼
5楼2009-07-22 21:31:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fo3mt2

金虫 (小有名气)
非常感谢你的回复。按照你说的方式做过了，可是Electronic Configuration中氧原子的U参数值仍然是5.00。我真的把我的全部操作过程都写在上面了，可还是不知道为什么会这样。我想问一下，如果我直接用vi编辑器将Rutile.cell文件中氧原子的U参数值删除，然后计算可以吗?那样做的话，会得到正确结果吗?希望大家予以指教，在此先谢过了!
赞一下 回复此楼
6楼2009-07-23 04:38:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crante

木虫 (正式写手)

登高
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuchenwf(金币+1,VIP+0):谢谢回帖交流 7-24 00:37
你把Electronic Configuration中氧原子的U参数值改为0 啊.
赞一下(11人) 回复此楼
努力才能成功
7楼2009-07-23 12:06:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crante

木虫 (正式写手)

登高
问下: 你提到Rutile.cell，Rutile.param，Rutile.xsd和Rutile.castep四个文件
其中Rutile.cell 在那看? 我只看到后面3个文件
赞一下 回复此楼
努力才能成功
8楼2009-07-23 12:09:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fo3mt2

金虫 (小有名气)
谢谢你的回复，我确实试过的你说的把氧原子的U参数值改为零。操作过程如下，选中钛原子，把它的U参数值改成5.00，然后选中氧原子，把它的U参数值改成0.00，然后其它设置都没变，保存计算任务之后，再重新选中被保存任务中的Rutile.xsd文件中的氧原子，在察看Electronic configuration,可是显示的却是5.00，我真的不知道该怎么做了。。。真的求各位朋友帮忙啦。。。
至于Rutile.cell文件我是上传到服务器上之后，在linux系统下打开我保存并上传上来的任务时看到的，你可以看我上传上来的附件里，有这个文件的内容。。。
还是在此求大家帮忙啦，万分感谢！
赞一下 回复此楼
9楼2009-07-23 14:34:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stou

银虫 (小有名气)
★ ★ ★ ★
小木虫(金币+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 ]
赞一下(48人) 回复此楼
10楼2009-07-24 10:29:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fo3mt2 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭