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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 结构优化是否要DFT+U
121/2返回列表12下一页
查看: 7938  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者学渣进行中将赠送您 5 个金币

学渣进行中

新虫 (小有名气)

[求助] 结构优化是否要DFT+U 已有2人参与

在关于一个含Fe的三维晶体用VASP计算中，考虑DFT+U，那是在结构优化中就加吗？还是算能带和太密度时再加？计算弹性常数是否+U?U对弹性常数有影响吗？U来自于文献
回复此楼

» 猜你喜欢
1楼 2022-11-28 23:18:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

wangyikeco

禁虫 (正式写手)
感谢参与，应助指数 +1
本帖内容被屏蔽
2楼2022-11-30 01:33:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
5楼: Originally posted by 学渣进行中 at 2022-12-10 11:47:45
好的，谢谢
...

弹性常数可能有影响，不要说弹性常数，是否进行DFT+U对晶胞常数都可能会造成影响。
你可以参考LDAUTYPE=2方法作者Dudarev的文章10.1103/PhysRevB.57.1505。
赞一下(2人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
6楼2023-01-05 13:33:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
10楼: Originally posted by wangyikeco at 2023-01-11 05:18:04
10.1103/PhysRevB.57.1505文章发表于1998年，有时代局限性，他用的是LDA (Local Density Approximation)，不是更好GGA (Generalized Gradient Approximation) functionals (like PBE)的。

对于金属或合金，用GG ...

1、LDA确实有其时代局限性，因为不考虑梯度修正从而导致overbinding的问题，GGA核心来说是对这个问题的修正。如果仅考虑势能面形状相关的问题，包括导数相关的性质，比如声子或者模量，确实会有改善。但是，overbinding的改善之外，交换-关联泛函的考量程度还包括自作用偏差，不管是LDA或是GGA甚至是后续的metaGGA，这一点迄今仍然可能是缺陷。我印象当中应该也有相关的benchmark有表现出GGA对一些吸附问题仍然有很大的偏差，对于金属（d、f族）及其氧化物中的局域磁矩可能会低估，以及窄带绝缘体的gap问题。

2、DFT+U可以视作为自作用修正（SIC）的方式之一，这一点在M Cococcioni（比如05年的PRL）以及A, Zunger一些文章都有提到。从广义Koopmans定理可衍生出线性响应方法和cRPA来获取U值的方案，这个过程不见得不同化学环境中的同种原子就被一视同仁的处理，类似方法比较成熟，还有没有更好的方案正在发展也拭目以待。当然，这两种方案也只是给U做合理使用（不凑数）的一种方式。因此金属或合金到底是不是GGA就足够了，要看讨论的问题而定。

3、DFT+GW或DFT+DMFT确实很好，但也很贵。另外一般的DMFT可能当中的U还是要通过类似上述方式确定的（本质上还是Hubbard模型），以及GW可能会用G0W0或者GW0，对初始波函数或密度的选取一定程度也有依赖。

4、U近些年在一些文献的使用确实有凑数的可能性，我们不做过多评价。如同我在9楼所说，DFT+U有诸多的局限性和任意性，除此之外，上述LDA、GGA或是mGGA等半局域泛函本身构造也没考虑DFT+U下去拟合参数，是否引入DFT+U而在double counting中消除不了，这样的问题也很难说。之所以这么做还是杂化泛函或者GW的代价还是非常昂贵，只能做一些妥协。

5、说回LZ的问题，LZ的问题应该是两个层次：
（1）到底弹性常数要不要加，我的评价是，加了有改善的可能性。那么会不会更差，也有可能，9楼也对这个事情简单阐述了一下。
（2）加了对弹性常数有没有影响，这是我上面一直在强调的，因为Hamiltonian都改变了，这样一来势能面的极小点位置已经能量及其导数值都可能被改变，我也拿出文献做例子，逻辑上说与LDA或者GGA甚至mGGA也没关系。主要是对“计算弹性常数不要U, U对弹性常数没有影响”这句话无法认同。
赞一下(2人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
11楼2023-01-11 11:06:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyikeco

禁虫 (正式写手)
本帖内容被屏蔽
10楼2023-01-11 05:18:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

wangyikeco

禁虫 (正式写手)
本帖内容被屏蔽
3楼2022-12-01 02:13:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

学渣进行中

新虫 (小有名气)
引用回帖:
2楼: Originally posted by wangyikeco at 2022-11-30 01:33:57
结构优化不要DFT+U, the structure optimization is based on the equation of motion or Newtonian/modular dynamics. Classical mechanics is more than enough.

好的，谢谢

发自小木虫Android客户端
回复此楼
4楼2022-12-10 11:47:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

学渣进行中

新虫 (小有名气)
引用回帖:
3楼: Originally posted by wangyikeco at 2022-12-01 02:13:54
算Band Structure和DOS时再加 LDA+U or GGA+U.
计算弹性常数不要U, U对弹性常数没有影响

好的，谢谢

发自小木虫Android客户端
回复此楼
5楼2022-12-10 11:47:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rlafite

至尊木虫 (正式写手)

【答案】应助回帖

引用回帖:
6楼: Originally posted by 卡开发发 at 2023-01-05 13:33:51
弹性常数可能有影响，不要说弹性常数，是否进行DFT+U对晶胞常数都可能会造成影响。
你可以参考LDAUTYPE=2方法作者Dudarev的文章10.1103/PhysRevB.57.1505。...

U的数量级通常为,零点几个eV，而对称性等约束达几十个eV(甚至好几百eV), 加U通常也只导致百分之零点几的误差，晶胞常数优化本来就是很粗糙的过程, 所以DFT+U对晶胞常数进行优化毫无意义。
如果机时用不掉可以加着玩
赞一下 回复此楼
7楼2023-01-06 02:01:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
7楼: Originally posted by rlafite at 2023-01-06 02:01:58
U的数量级通常为,零点几个eV，而对称性等约束达几十个eV(甚至好几百eV), 加U通常也只导致百分之零点几的误差，晶胞常数优化本来就是很粗糙的过程, 所以DFT+U对晶胞常数进行优化毫无意义。
如果机时用不掉可以加着 ...

LZ问的是弹性模量的问题，从10.1103/PhysRevB.57.1505来说，DFT+U对弹性模量不止有影响，影响还不小，可不是百分之0.几的误差。正常来说，晶胞常数的优化可不是一个很粗糙的过程，尤其是模量、声子谱等可能还会更敏感。也许对于高通量的计算可以适当放宽条件，但是不表示一般的计算能够随意进行。

另外，如果按照0.几eV能量不重要，那么不同的交换关联泛函对相对能量影响都可能不超过这个范围，甚至基函数尺寸、k点测试都不重要，但实际情况显然不是这样，不过不慎重对待很可能得到一个不可靠的结果。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
8楼2023-01-06 06:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
这里额外做一些说法上的补充，并不是说DFT+U一定能保证正向改善结果，而是说，因为DFT+U改变了Hamilton，会影响到能量和核坐标的关系。U是个环境依赖的参数，如何合理设置大有讲究，尤其对于化学反应的过程，如何合理设置，迄今仍然是有学术争议的。除了U值本身有任意性，局域投影态、重复计入修正的方案也都不是唯一的。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
9楼2023-01-09 11:17:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学渣进行中 的主题更新
121/2返回列表12下一页
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +11 kulium 2026-08-21 14/700 2026-08-24 02:12 by 流流伤
[基金申请] 范进中举一文的中心思想 +4 炎黄贵胄 2026-08-22 5/250 2026-08-24 02:09 by 流流伤
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 让我中一个面上吧！ +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 什么时候开奖？ +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭