24小时热门版块排行榜    

查看: 7985  |  回复: 11
【悬赏金币】回答本帖问题,作者学渣进行中将赠送您 5 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

学渣进行中

新虫 (小有名气)

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

在关于一个含Fe的三维晶体用VASP计算中,考虑DFT+U,那是在结构优化中就加吗?还是算能带和太密度时再加?计算弹性常数是否+U?U对弹性常数有影响吗?U来自于文献
回复此楼
已阅   回复此楼   关注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对弹性常数没有影响”这句话无法认同。
不一定挂在论坛,计算问题问题欢迎留言。
11楼2023-01-11 11:06:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

wangyikeco

禁虫 (正式写手)

感谢参与,应助指数 +1
本帖内容被屏蔽

2楼2022-11-30 01:33:57
已阅   回复此楼   关注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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jCd0dEvKHShX 2026-08-29 3/150 2026-08-30 19:26 by hZiFeudZyoGR
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 19:07 by xiangy672
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 18:51 by hZiFeudZyoGR
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 17:34 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:03 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 14/700 2026-08-30 17:01 by hZiFeudZyoGR
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
请填处理意见