24小时热门版块排行榜    

查看: 4901  |  回复: 35
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

刘仕晨

木虫 (正式写手)

[求助] 新手求助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文件中没有定义我需要的元素,怎么修改阿,直接添加上就可以了么,还是怎么弄,有没有人修改过阿。
问题比较多,十分感谢!
回复此楼

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

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的回帖
查看全部 36 个回答

刘仕晨

木虫 (正式写手)

求各位大神帮忙啊,谢谢啦……
2楼2013-03-21 12:04:05
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 上海工程技术大学 激光智能制造课题组 2027级博士研究生招生 +5 oxidpl 2026-09-20 7/350 2026-09-27 10:43 by Equinoxhua
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:42 by 4IW0sJvtEqX8
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:26 by 4IW0sJvtEqX8
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:21 by 4IW0sJvtEqX8
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:15 by 4IW0sJvtEqX8
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:03 by 4IW0sJvtEqX8
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 oeyPlfmMUpOC 2026-09-26 3/150 2026-09-27 09:13 by 4IW0sJvtEqX8
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:42 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:39 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:44 by XgCC8uTcwILl
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 7/350 2026-09-25 11:18 by yhq9807
[基金申请] 师弟论文见刊大半年才想起来申请专利,还能抢救一下吗? +4 13108017953 2026-09-21 4/200 2026-09-24 10:06 by kudofaye
信息提示
请填处理意见