24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1937  |  回复: 19

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电子还不满,
可我怎么优化晶格参数都是偏实验值小很多。不知道该怎么办。
可能得搞明白这个结构优化是如何优化,我们设得是ISIF=3
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 的主题更新
信息提示
请填处理意见