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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » MoF3用 LDA+U ，U和J值该设置多少
202/2返回列表上一页12
查看: 1937  |  回复: 19
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

umklapp

银虫 (小有名气)
★ ★ ★ ★
sunyang1988: 金币+2, 谢谢交流 2012-07-22 14:26:23
franch: 金币+2, 谢谢回帖交流。 2012-07-28 18:11:45
引用回帖:
7楼: Originally posted by aylixmut at 2012-07-21 14:44:43
什么是Open Core calculation，具体怎么计算？，还请多多指点一下，谢谢。...

就是把的d或者f 电子放到core 里去按照U的定义去计算（两个电离能之差）， 这样就相当于关闭了电子的 hybridization， U值比较可靠。
但是
1） 要注意你的离子上 有多少电子占据 d/f态;
2） VASP应该不能实现 Open Core, 需要其它code;
赞一下 回复此楼
11楼2012-07-22 13:05:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylixmut

铜虫 (小有名气)
引用回帖:
11楼: Originally posted by umklapp at 2012-07-22 13:05:52
就是把的d或者f 电子放到core 里去按照U的定义去计算（两个电离能之差）， 这样就相当于关闭了电子的 hybridization， U值比较可靠。
但是
1） 要注意你的离子上 有多少电子占据 d/f态;
2） VASP应该不能实现  ...

我们是用vasp算，原来vasp这么弱，各位大牛介绍的方法，都难实现。悲催啊！
赞一下 回复此楼
12楼2012-07-22 14:41:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

itpfeng

木虫 (著名写手)
★
sunyang1988: 金币+1, 谢谢交流 2012-07-24 00:06:24
引用回帖:
5楼: Originally posted by souledge at 2012-07-21 14:39:34
我用EXX算个简单的AlAs，结果算了两天还没出结果，果断中止了……以前也试过用HSE计算一个简单体系，跑了一夜一个SCF步都没跑完……
话说，+U和Hybrid都是被看作为trick的~还是直接GW吧...

虽然简单但+U还是有一定的物理思想在内的，hybrid则是因为拟合不好而纯正在凑。如果+U被认为是trick的话，那GW就也是trick了
赞一下 回复此楼
13楼2012-07-23 22:42:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★ ★
franch: 金币+2, 谢谢回帖交流。 2012-07-28 18:12:08
引用回帖:
13楼: Originally posted by itpfeng at 2012-07-23 22:42:34
虽然简单但+U还是有一定的物理思想在内的，hybrid则是因为拟合不好而纯正在凑。如果+U被认为是trick的话，那GW就也是trick了...

也是~至少对于+U问题，还有非整数占据和屏蔽库仑力之类的解释~Hybrid基本纯粹是因为HF（高估）和DFT（低估）的中和一下可以得到好的带隙~~~
赞一下 回复此楼
思想重于技巧，内涵重于表象
14楼2012-07-24 00:53:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmr0125

木虫 (正式写手)

【答案】应助回帖

★ ★
franch: 金币+2, 鼓励交流 2012-07-27 13:02:05
查文献去看，具体不好说。我算的MoO3的U-J=6.3.U的值不一定，J在0.9左右就可以的。可以去查一些文献，具体的设置可以参照手册。自己看懂还是很容易的。
赞一下 回复此楼
人生如此！
15楼2012-07-27 09:38:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)
Mo 一般不用加U吧
赞一下 回复此楼
16楼2012-07-28 08:42:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylixmut

铜虫 (小有名气)
★ ★
franch: 金币+2, 鼓励交流 2012-07-28 18:12:21
引用回帖:
16楼: Originally posted by jghe at 2012-07-28 08:42:27
Mo 一般不用加U吧

不加U最好，可是为什么不加U啊，
Mo有d电子还不满，
可我怎么优化晶格参数都是偏实验值小很多。不知道该怎么办。
可能得搞明白这个结构优化是如何优化，我们设得是ＩＳＩＦ＝３
赞一下 回复此楼
17楼2012-07-28 16:30:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunny、小芳

新虫 (小有名气)
引用回帖:
15楼: Originally posted by xmr0125 at 2012-07-27 09:38:11
查文献去看，具体不好说。我算的MoO3的U-J=6.3.U的值不一定，J在0.9左右就可以的。可以去查一些文献，具体的设置可以参照手册。自己看懂还是很容易的。

您好，您算MoO3的时候，U=6.4，J=0.9吗？计算得到的带隙是多少？我也查文献看到大家都用6.3，但是我计算得到的带隙还是比较小……跟实验值差很多。
赞一下 回复此楼
18楼2019-03-31 10:07:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
18楼: Originally posted by sunny、小芳 at 2019-03-31 10:07:48
您好，您算MoO3的时候，U=6.4，J=0.9吗？计算得到的带隙是多少？我也查文献看到大家都用6.3，但是我计算得到的带隙还是比较小……跟实验值差很多。...

你计算参数与文献相同吗？你使用的程序是？如果是VASP的话，我遇到的一种可能是对d、f电子系统需要分别指定LMAXMIX=4、6，在DFT+U的情况下这东西不止会影响数值稳定性。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
19楼2019-03-31 12:23:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
前面的帖子我说下自己的观点：
1、11楼的做法是有的，我没记错的话可能可以是用在LMTO或LAPW上面。
2、线性响应方法QE能做，实际上VASP也能做，可以用LDAUTYPE=3来实现，也有人写过教程http://grandcentral.apam.columbi ... sponse_u/index.html，SIESTA在4.1-b2过后的版本也可以使用这种方法，如果不是很熟悉的话有机会我可以单独开帖子说。其他可能还要程序能做，只是我不太清楚。除了上面两种其实还有通过Constrainted RPA的方法，但我不太熟悉。
3、关于DFT+U的意义，+U或者说Hubbard模型本身确实是个模型而不是严格的理论，但不能说+U没有物理意义。从自作用修正的角度，Cococcioni的关于线性响应计算U的文章已经解释U应当看做DFT在电子增减过程中不服从广义Koopmans定理的补偿，大体上说杂化泛函在理论上也是相同的思路；从GW的角度来看+U可以看成是零频近似，这在Anisimov的提DFT+U的文章也提到过。当然，DFT+U有他的缺陷，比如：1、U值选取实际上操作是个复杂的过程；2、其实某些部分还是有任意性，将能带在投影到局域轨道的时候，局域轨道的选取不是唯一的，还有double count的部分也不是唯一的；3、+U主要修正的还是单粒子能级，对能量、能垒方面的修正到底是好是坏是个未知数，Truhlar做过这方面的测试，结果显示+U对能垒不一定是改善的。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
20楼2019-03-31 12:46:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aylixmut 的主题更新
202/2返回列表上一页12
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭