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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP中怎么进行加U计算？
191/2返回列表12下一页
查看: 5719  |  回复: 18
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

niuxiaoqi721

铜虫 (小有名气)

[求助] VASP中怎么进行加U计算？

VASP中怎么进行加U计算？INCAR中我需要添加什么命令，加U影响对晶胞的结构优化吗？
回复此楼

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

vasp光学 第一性原理计算 DFT+U，HSE等

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
加油
1楼2012-06-28 15:41:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢回帖！ 2012-06-29 19:46:39
niuxiaoqi721: 金币+10, ★★★★★最佳答案, 哈哈 2012-07-05 09:17:17
The L(S)DA+U in VASP is switched on by means of the following tags


LDAU = .TRUE. Switches on the L(S)DA+U.
LDAUTYPE = 1|2|4 Type of L(S)DA+U (Default: LDAUTYPE = 2)
1
Rotationally invariant LSDA+U according to Liechtenstein et al.
4
Idem 1., but LDA+U instead of LSDA+U (i.e. no LSDA exchange splitting)
2
Dudarev's approach to LSDA+U (Default)
LDAUL = L .. -quantum number for which the on site interaction is added
(-1: no on site terms added, 1: p, 2: d, 3: f, Default: LDAUL = 2)
LDAUU = U .. Effective on site Coulomb interaction parameter
LDAUJ = J .. Effective on site Exchange interaction parameter
LDAUPRINT = 0|1|2 Controls verbosity of the L(S)DA+U module
(0: silent, 1: Write occupancy matrix to OUTCAR, 2: idem 1., plus potential matrix dumped to stdout, Default: LDAUPRINT = 0)
NB: LDAUL, LDAUU, and LDAUJ must be specified for all atomic species!
加U对晶格影响不会太大。
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

niuxiaoqi721
2楼2012-06-29 19:18:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★
xueht987: 金币+1, 应助指数+1, 多谢指教！ 2012-07-04 09:38:40
引用回帖:
6楼: Originally posted by niuxiaoqi721 at 2012-07-03 09:20:14
我需要添加的命令就是：
LDAU = .TRUE.
LDAUTYPE =2
LDAUL = 2
LDAUU =
LDAUJ =
LDAUPRINT = 0
按照说明书上的说法，就是我应该将上述命令加入到INCAR中，对吗？ 我计算的氧化铁，铁是d电子的，我现在只知 ...

从计算的最开始就加上吧，多一事不如少一事。
假设原子是Fe和O，就写成
LDAUL = 2 -1
LDAUU = 5 0
LDAUJ = 0 0
5eV的数值最好测试一下。通常是dos和实验值或HSE比较。

顺带说一句，如果想对某些Fe加U，某些不加（慎用！！！！一般是吸附原子和衬底），就把它们当成两种原子分别写出，POTCAR和POSCAR都这么做就行了。虽然说明书没说过，我测试了是可以的。
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

niuxiaoqi721 niuxiaoqi721
7楼2012-07-03 16:11:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★
xueht987: 金币+1, 多谢指教！ 2012-07-05 08:59:20
xueht987: 应助指数+1 2012-07-05 08:59:28
引用回帖:
16楼: Originally posted by wenyan0517 at 2012-07-04 21:23:24
请教这两个参数怎么确定呢？
LDAUU和LDAUJ  
凭经验，还是可以通过什么算出来？ 谢谢...

对LDAUTYPE = 2（常用）只有U-J才是有意义的。
经验值。三种方法：
1.实验上的gap多大，自己调到正好
2.测试电子结构、峰位移动，使+U结果靠近HSE
3.文献用多少，直接照用不解释
赞一下(1人) 回复此楼
17楼2012-07-04 22:49:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

umklapp

银虫 (小有名气)
有人说U对结果影响不大， 弛豫不用加U；不过我从来都是加U弛豫的；
赞一下 回复此楼
3楼2012-06-29 23:59:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2012-07-02 18:46:26
引用回帖:
3楼: Originally posted by umklapp at 2012-06-29 23:59:13
有人说U对结果影响不大， 弛豫不用加U；不过我从来都是加U弛豫的；

嗯，如果最终是要加U，弛豫最好就加上（免得静态计算出来比较大的受力烦心）。+U的计算量很少的。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

niuxiaoqi721
4楼2012-07-02 18:28:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niuxiaoqi721

铜虫 (小有名气)
送鲜花一朵
引用回帖:
4楼: Originally posted by mousekingadv at 2012-07-02 18:28:25
嗯，如果最终是要加U，弛豫最好就加上（免得静态计算出来比较大的受力烦心）。+U的计算量很少的。...

意思是我优化晶胞的时候可以不+U，当我进行切面，然后有固定层和驰域层的时候，此时优化结构就得进行+U，是吗
赞一下 回复此楼
加油
5楼2012-07-03 09:04:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niuxiaoqi721

铜虫 (小有名气)
送鲜花一朵
引用回帖:
2楼: Originally posted by mousekingadv at 2012-06-29 19:18:47
The L(S)DA+U in VASP is switched on by means of the following tags


LDAU = .TRUE. Switches on the L(S)DA+U.
LDAUTYPE = 1|2|4 Type of L(S)DA+U (Default: LDAUTYPE = 2)
1
Rotationally invaria ...

我需要添加的命令就是：
LDAU = .TRUE.
LDAUTYPE =2
LDAUL = 2
LDAUU =
LDAUJ =
LDAUPRINT = 0
按照说明书上的说法，就是我应该将上述命令加入到INCAR中，对吗？ 我计算的氧化铁，铁是d电子的，我现在只知道U-J为5ev就可以，可中间的命令怎么填，还有就是不是所有原子都要定义吗？氧原子怎么办？我确实新手，一点不懂
赞一下 回复此楼
加油
6楼2012-07-03 09:20:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niuxiaoqi721

铜虫 (小有名气)
送鲜花一朵
引用回帖:
7楼: Originally posted by mousekingadv at 2012-07-03 16:11:03
从计算的最开始就加上吧，多一事不如少一事。
假设原子是Fe和O，就写成
LDAUL = 2 -1
LDAUU = 5 0
LDAUJ = 0 0
5eV的数值最好测试一下。通常是dos和实验值或HSE比较。

顺带说一句，如果想对某些Fe加U，某些 ...

+U时这3个参数是对所有种类原子设置吗？那LDAU = .TRUE.  LDAUTYPE =2   LDAUPRINT = 0  就不用对各个原子分别设置了哦；还有一点就是+U是针对磁性原子才使用吗？如果我吸附一个汞原子，那么就是LDAUL = 2 -1 -1   LDAUU = 5 0 0  LDAUJ = 0 0 0，这样对吗？
赞一下 回复此楼
加油
8楼2012-07-03 17:07:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fyh987

金虫 (小有名气)
引用回帖:
4楼: Originally posted by mousekingadv at 2012-07-02 18:28:25
嗯，如果最终是要加U，弛豫最好就加上（免得静态计算出来比较大的受力烦心）。+U的计算量很少的。...

请问一下什么是弛豫啊？为什么我加U计算时时间特别长啊，老是不收敛，LZ计算的是什么体系
赞一下 回复此楼
努力挑战自己
9楼2012-07-03 17:15:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niuxiaoqi721

铜虫 (小有名气)
引用回帖:
9楼: Originally posted by fyh987 at 2012-07-03 17:15:48
请问一下什么是弛豫啊？为什么我加U计算时时间特别长啊，老是不收敛，LZ计算的是什么体系...

驰域就是不固定这些原子啊，在poscar中可以设置，至于+U方面的计算我完全不懂，正在咨询
赞一下 回复此楼
加油
10楼2012-07-03 17:24:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 niuxiaoqi721 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 江南大学到瑞士招聘，称取消非升即走，改预聘+长聘 +22 babu2015 2024-06-18 24/1200 2024-06-21 10:37 by hktk1987
[基金申请] 青年和面上，哪个上会难度更大 +6 今晚推荐22 2024-06-21 6/300 2024-06-21 10:18 by migrant09
[基金申请] 江苏省333人才工程 出校后被刷的概率大不大？ +4 maxbirdzhang 2024-06-19 5/250 2024-06-21 08:57 by maxbirdzhang
[考博] 一名额难求啊 +9 pinbo拼搏 2024-06-19 9/450 2024-06-21 08:54 by 朔月要疯
[精细化工] 试剂价格 +7 旭必上岸 2024-06-17 10/500 2024-06-20 23:35 by 相与处于陆的鱼
[找工作] 初始合伙人来啦！（生物试剂耗材标准品） +16 欢快的小科研人 2024-06-15 27/1350 2024-06-20 22:19 by zzw1992127
[基金申请] 我标书代码变了 +63 学员NHuqdk 2024-06-16 86/4300 2024-06-20 22:03 by 841790061
[有机交流] 想要用氢化钠拔掉吲哚N上的氢取代酰氯 50+3 光敏剂 2024-06-19 4/200 2024-06-20 18:41 by HF111001
[论文投稿] 投稿求助 +4 平凡的日子 2024-06-19 5/250 2024-06-20 16:24 by yueyueyue@
[教师之家] 每次骚扰女学生的都是院系领导，而不是普通教师，小编们要注意措辞正确 +11 zju2000 2024-06-15 13/650 2024-06-20 15:19 by ase123456
[精细化工] 切削液交流群 +5 SZ吴彦祖 2024-06-17 5/250 2024-06-20 12:10 by 毛爱波BBB
[基金申请] F03青年基金函评结果 +3 暨阳一只柴 2024-06-19 4/200 2024-06-20 11:15 by 暨阳一只柴
[催化] 镍负载氧化铝的保存问题 8+3 lwn0130 2024-06-15 6/300 2024-06-20 09:00 by lwn0130
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中！ +10 恐龙爸爸 2024-06-14 10/500 2024-06-20 08:42 by jsguo4127
[论文投稿] 审稿 +5 香瓜木香 2024-06-19 6/300 2024-06-19 17:44 by xli1984
[基金申请] 博后基金刷到的BUG,图片来的更直观 +15 carolloo 2024-06-17 16/800 2024-06-19 09:42 by msjy
[高分子] 烧瓶内合成聚酯 +3 大帝国乐 2024-06-17 7/350 2024-06-18 11:10 by 大帝国乐
[有机交流] 求教最后一步问题出在哪 20+3 苏州浩恒 2024-06-15 4/200 2024-06-18 08:33 by like_2005
[硕博家园] 博士毕业高校和就业的相关问题 +7 SCITOPPP 2024-06-14 11/550 2024-06-18 07:51 by yinxing1995
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
信息提示
关闭
请填处理意见
关闭