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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 新手求助pwscf的进行磁性运算和LDA+U运算
361/4返回列表1234下一页
查看: 3999  |  回复: 35
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

刘仕晨

木虫 (正式写手)

[求助] 新手求助pwscf的进行磁性运算和LDA+U运算

组里一直用的vasp，现在导师让学习pwscf，自己在看example和小木虫，自己捣鼓，但是最近在计算费米面的时候遇到一些问题，望各位大神能指导指导，谢谢！
1，考虑磁性的时候我是画CrO2的费米面，对Cr加了铁磁，在starting_magnetization搞了好久，因为一直没有从vasp转过来，这里是自旋极化率，值在-1到1之间，刚开始纠结好了好久，怎么会在-1到1之间呢，那么问题来了，这个值我到底怎么确定呢，如果是按计算公式，在费米能级上P=/,那绝缘体怎么办，这个值是多少呢？
2，计算的out中就给出了absolute magnetization到底是怎么回事，total magnetization这个我已经能理解，关键的是pwscf不能像vasp那样给出每个离子的磁矩，不知道每个离子的磁矩要怎么算，看到一个帖子里是这样写的，但是不是很懂，“tot_magnetization=N_total majority spin charge - N_minority spin charge
                     absolute_magnetization=|N_total majority spin charge -N_minority spin charge|
                    例如体系为铁磁，tot_magnetization=absolute_magnetization
                    体系为反铁磁性，tot_magnetization=0, absolute_magnetization= 2 * N_total majority spin charge”
对absolute magnetization还是没有感觉，也不知道单个离子的磁矩要怎么得到？
3，关于LDA+U，看了侯老子的博客，感觉是不是只要设置lda_plus_u=.true. ，再指定Hubbard_U(I)的值就可以了，因为vasp里还需要指定了是在在d或者f上加U，vasp里加U的同时要设置J的值，pwscf是不是不需要设置阿？
4，侯老师同时还给出了需要可选的参数，那些参数到底在什么时候需要考虑阿，比如：Hubbard_alpha(I)，starting_ns_eigenvalue(m,ispin,I) ，
5，如果PW/set_hubbard_l.f90 和PW/tabd.f90文件中没有定义我需要的元素，怎么修改阿，直接添加上就可以了么，还是怎么弄，有没有人修改过阿。
问题比较多，十分感谢！
回复此楼

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

第一性原理和电化学 第一性原理

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-03-20 19:12:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

刘仕晨

木虫 (正式写手)
引用回帖:
7楼: Originally posted by huazhorg at 2013-03-25 20:42:33
1. starting_magnetization只是个初始猜测，【-1，1】以内随便设一个，当然不能为0， 一般最终会收敛到正确的结果上去的。
2. total magnetization是矢量和，absolute magnetization可以认为是取摸求和，一般大于或 ...

1.第4个问题，我不是很懂，我想问下您的意思是，考虑强关联的时候只需要设置U ，不需要考虑其他那些参数，还是在强关联的时候要考虑U，我没有很懂你的意思，就是其他那些可选的参数都选默认呢，还是在什么时候是需要设置的？
2.现在还有一个新的问题，我在计算CrO2的时候，我是用的LDA，模守恒的势，ecutwfc=30.不知道是不是小了点，因为我看到example里的Ni就放了24，就没放多大，但是当我变成100的时候，发现fermi energy 和总能量都有很大的变化，包括磁矩都在变化，是不是ecutwfc太小的原因……
麻烦了谢谢！
赞一下(1人) 回复此楼
8楼2013-03-25 21:50:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
14楼: Originally posted by 刘仕晨 at 2013-03-27 15:39:27
您好，我在测试ecut和K点的时候遇到了新的问题，我是这样做的，分别设ecut=30 40 50 60 70 80 90 100 110 120,然后在每个ecut值下面慢慢增加K点，但是奇怪的是，在ecut=50 60和70的时候自洽运算会出现对称性操作丢 ...

你用8.3503乘以0.3294852，看看跟5.5026乘以0.5是不是一样的，就是给你换了坐标而已。
赞一下(1人) 回复此楼
18楼2013-03-28 10:23:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
33楼: Originally posted by 刘仕晨 at 2013-05-08 13:11:16
我明白您的意思，因为我曾经问过您这个问题，但是自洽计算的晶格常数和full-relax后的晶格常数是一样的，那么Cr的位置应该是（0.5 0.5 0.3353）对吧，而不是自洽中的(   0.5080391   0.5080391   0.3406890  )我不 ...

呵呵，自洽计算给出的不是分数坐标，它的总长度不是1，而是里面的a(1)，即1.016078，所以你看看0.5080391是不是它的一半。那么化成crystal坐标是不是就变成了0.5。
赞一下(1人) 回复此楼
34楼2013-05-08 14:56:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

刘仕晨

木虫 (正式写手)
求各位大神帮忙啊，谢谢啦……
赞一下 回复此楼
2楼2013-03-21 12:04:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weixp2008

木虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
fzx2008: 金币+1, 鼓励交流！ 2013-03-21 19:42:22
磁性计算，个人以为不太靠谱
赞一下 回复此楼
3楼2013-03-21 12:34:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)
引用回帖:
3楼: Originally posted by weixp2008 at 2013-03-21 12:34:23
磁性计算，个人以为不太靠谱

这个我也是初学者不是很懂，但是我们用vasp算的时候都加了磁性的，现在用pwscf画费米面的话，应该也是要加磁性加U吧，这个磁性计算有多少不靠谱呢……
赞一下 回复此楼
4楼2013-03-21 13:31:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
我这两天帮你问问那些算强关联的人，我一直在算很轻的元素。
赞一下 回复此楼
5楼2013-03-21 21:03:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)
引用回帖:
5楼: Originally posted by gemucai at 2013-03-21 21:03:17
我这两天帮你问问那些算强关联的人，我一直在算很轻的元素。

谢谢，谢谢，麻烦您了！
赞一下 回复此楼
6楼2013-03-22 14:27:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
刘仕晨: 金币+10, ★★★★★最佳答案 2013-03-25 21:40:56
1. starting_magnetization只是个初始猜测，【-1，1】以内随便设一个，当然不能为0， 一般最终会收敛到正确的结果上去的。
2. total magnetization是矢量和，absolute magnetization可以认为是取摸求和，一般大于或等于前者。
3.只需要设U，不需要设J。U的值，可以根据文献的值，也可以自己测试，我认为，能够使结果得到正确的磁矩的最小U值就够了。
4. 涉及到过渡金属活强关联体系，在你想得到正确的磁性性质或能带结构，就需要设置U。
5. 最近版的espresso里比较全，一般的过渡金属都考虑进去了，Cr自然也考虑进去了。如果万一没考虑，结果会有出错提示，让你修改什么文件，然后重新编译。
赞一下 回复此楼
7楼2013-03-25 20:42:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)
引用回帖:
8楼: Originally posted by 刘仕晨 at 2013-03-25 21:50:14
1.第4个问题，我不是很懂，我想问下您的意思是，考虑强关联的时候只需要设置U ，不需要考虑其他那些参数，还是在强关联的时候要考虑U，我没有很懂你的意思，就是其他那些可选的参数都选默认呢，还是在什么时候是需 ...

强关联体系，最好加U，另外要开启spin极化，无论这个体系有没有磁性，开启spin都不会影响其结果；相反，如果体系有磁性但没有开启spin的话，结果肯定就是有问题的了。所以一般spin +U就够了。如果要用模守恒赝势，30肯定小了，一般在90以上。你设100发现结果相差很大，就是这个原因：之前的参数设小了，结果根本是错的。
赞一下 回复此楼

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

迷迷物理
9楼2013-03-25 22:13:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)
引用回帖:
9楼: Originally posted by huazhorg at 2013-03-25 22:13:02
强关联体系，最好加U，另外要开启spin极化，无论这个体系有没有磁性，开启spin都不会影响其结果；相反，如果体系有磁性但没有开启spin的话，结果肯定就是有问题的了。所以一般spin +U就够了。如果要用模守恒赝势， ...

谢谢，谢谢，受益匪浅，对我这个新手来说收获很多，我现在怀疑我的K点也设的不够，所以我现在在测试K点，画出K点和能量的关系，然后拿那个K点再来加大ecut，谢谢
那在加U同时开启spin，那这些参数呢Hubbard_alpha(I)，starting_ns_eigenvalue(m,ispin,I) 等等的，这些参数怎么设置阿，还是用默认值阿，还是不用管？什么时候需要设阿，谢谢
赞一下 回复此楼
10楼2013-03-25 22:21:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 刘仕晨 的主题更新
361/4返回列表1234下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭