| 查看: 4803 | 回复: 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文件中没有定义我需要的元素，怎么修改阿，直接添加上就可以了么，还是怎么弄，有没有人修改过阿。 问题比较多，十分感谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理和电化学 | 第一性原理 |
» 猜你喜欢
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
能否退出参与的面上项目解除限项
已经有17人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助pwscf中tools/pwi2xsf.sh 的用法 谢谢
已经有6人回复
- [求助]PWscf安装
已经有11人回复
- 【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？
已经有14人回复
- 【求助】pwscf错误提示diagonalization(ZHEGV*) falied
已经有8人回复
- 【求助】求助pwscf一些参数的物理意义和设定经验
已经有13人回复
- 【求助】pwscf结构优化结果
已经有4人回复
- 【求助】pwscf中的问题
已经有4人回复
- 【求助】pwscf结构优化中结果出的怪异啊
已经有16人回复
- 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算
已经有23人回复
- 【求助】PWSCF输出的projwfc文件中，哪些文件是dxy、dxz、dyz轨道？
已经有10人回复
- 【求助】pwscf怎样计算介电常数和波恩有效电荷
已经有13人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
- 【求助】酬谢20金币紧急求助pwscf中用ld1.x生成K赝势
已经有11人回复
刘仕晨
木虫 (正式写手)
- 应助: 30 (小学生)
- 金币: 3276.8
- 散金: 6
- 红花: 11
- 帖子: 319
- 在线: 179.2小时
- 虫号: 2201333
- 注册: 2012-12-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2013-03-25 21:50:14
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
18楼2013-03-28 10:23:11
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
34楼2013-05-08 14:56:14
刘仕晨
木虫 (正式写手)
- 应助: 30 (小学生)
- 金币: 3276.8
- 散金: 6
- 红花: 11
- 帖子: 319
- 在线: 179.2小时
- 虫号: 2201333
- 注册: 2012-12-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2013-03-21 12:04:05
weixp2008
木虫 (著名写手)
- 应助: 15 (小学生)
- 金币: 9742.1
- 散金: 75
- 红花: 2
- 帖子: 1190
- 在线: 459.7小时
- 虫号: 753583
- 注册: 2009-04-21
- 专业: 凝聚态物性 II ：电子结构
3楼2013-03-21 12:34:23
刘仕晨
木虫 (正式写手)
- 应助: 30 (小学生)
- 金币: 3276.8
- 散金: 6
- 红花: 11
- 帖子: 319
- 在线: 179.2小时
- 虫号: 2201333
- 注册: 2012-12-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2013-03-21 13:31:24
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
5楼2013-03-21 21:03:17
刘仕晨
木虫 (正式写手)
- 应助: 30 (小学生)
- 金币: 3276.8
- 散金: 6
- 红花: 11
- 帖子: 319
- 在线: 179.2小时
- 虫号: 2201333
- 注册: 2012-12-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2013-03-22 14:27:45
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4953.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
刘仕晨: 金币+10, ★★★★★最佳答案 2013-03-25 21:40:56
刘仕晨: 金币+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
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4953.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
迷迷物理9楼2013-03-25 22:13:02
刘仕晨
木虫 (正式写手)
- 应助: 30 (小学生)
- 金币: 3276.8
- 散金: 6
- 红花: 11
- 帖子: 319
- 在线: 179.2小时
- 虫号: 2201333
- 注册: 2012-12-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2013-03-25 22:21:40
20