24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

刘仕晨

木虫 (正式写手)

应助: 30 (小学生)
金币: 3276.8
散金: 6
红花: 11
帖子: 319
在线: 179.2小时
虫号: 2201333
注册: 2012-12-23
性别: GG
专业: 凝聚态物性 II ：电子结构

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

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有41人回复
回收溶剂求助已经有7人回复
硝基苯如何除去已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
三无产品还有机会吗已经有6人回复

1楼 2013-03-20 19:12:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

应助: 30 (小学生)
金币: 3276.8
散金: 6
红花: 11
帖子: 319
在线: 179.2小时
虫号: 2201333
注册: 2012-12-23
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

32楼: Originally posted by gemucai at 2013-05-08 08:24:51
记住，alat units的意思是以你晶格常数为单位，而在最后一步优化中，晶格常数发生了变化，你应该用最后一步的原子坐标乘以最后一步的晶格常数，看看得到的是不是自洽计算的位置。...

我明白您的意思，因为我曾经问过您这个问题，但是自洽计算的晶格常数和full-relax后的晶格常数是一样的，那么Cr的位置应该是（0.5 0.5 0.3353）对吧，而不是自洽中的( 0.5080391 0.5080391 0.3406890 )我不明白的是，为什么0.5的位置变成了0.5080391，这是我糊涂的地方，麻烦您了，谢谢！