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

回复此楼

» 猜你喜欢

药学硕士，第一、第二作者已发表6 篇 SCI，药理方向及相关方向2026年/2027年博士申请已经有4人回复
26年博士申请自荐-电催化已经有3人回复
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业已经有6人回复
收到国自然专家邀请后几年才会有本子送过来评已经有4人回复
考博已经有5人回复
26年申博自荐-计算机视觉已经有4人回复
药化及相关博士的申请已经有3人回复
一篇MDPI论文改变了学习工作和生活已经有4人回复

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 ：电子结构

引用回帖:

29楼: Originally posted by gemucai at 2013-04-27 09:04:46
当然两个都可以，就看你喜欢用什么样的坐标系了...

不好意思，又要麻烦您了，前几天在学习对体系进行full-relax,用vc-relax进行计算，在优化后它会自动进行一次自洽计算，但是奇怪的是，两次的原子位置却有点不大一样，我怎么看也不明白为什么，所以要麻烦您了，我给您把数据贴出来：
End of BFGS Geometry Optimization

   Final enthalpy = -167.8961699933 Ry
Begin final coordinates
   new unit-cell volume = 409.58929 a.u.^3 ( 60.69487 Ang^3 )

CELL_PARAMETERS (alat=  8.35030000)
1.016078133 0.000000000 0.000000000
0.000000000 1.016078133 0.000000000
0.000000000 0.000000000 0.681378067

ATOMIC_POSITIONS (crystal)
Cr    0.000000000 0.000000000  -0.000000000
Cr    0.500000000 0.500000000 0.500000000
O       0.304047718 0.304047718  -0.000000000
O       0.695952282 0.695952282 0.000000000
O       0.804047718 0.195952282 0.500000000
O       0.195952282 0.804047718 0.500000000
End final coordinates

   A final scf calculation at the relaxed structure.
   The G-vectors are recalculated for the final unit cell
   Results may differ from those at the preceding step.

   Parallelization info
   --------------------
   sticks: dense  smooth    PW    G-vecs: dense smooth    PW
   Min       858    858 234             35945 35945 5075
   Max       859    859 236             35946 35946 5076
   Sum       6869 6869 1877             287563 287563 40605

   bravais-lattice index    =          0
   lattice parameter (alat)  =    8.3503  a.u.
   unit-cell volume       =    409.5893 (a.u.)^3
   number of atoms/cell    =          6
   number of atomic types =          2
   number of electrons    =       36.00
   number of Kohn-Sham states=          22
   kinetic-energy cutoff    =    300.0000  Ry
   charge density cutoff    = 1200.0000  Ry
   convergence threshold    =    3.8E-12
   mixing beta             =    0.7000
   number of iterations used =          8  plain    mixing
   Exchange-correlation    = SLA-PZ-NOGX-NOGC ( 1 1 0 0 0)
   EXX-fraction             =       0.00

   celldm(1)= 8.350300  celldm(2)= 0.000000  celldm(3)= 0.000000
   celldm(4)= 0.000000  celldm(5)= 0.000000  celldm(6)= 0.000000

   crystal axes: (cart. coord. in units of alat)
            a(1) = ( 1.016078 0.000000 0.000000 )
            a(2) = ( 0.000000 1.016078 0.000000 )
            a(3) = ( 0.000000 0.000000 0.681378 )

   reciprocal axes: (cart. coord. in units 2 pi/alat)
            b(1) = (  0.984176  0.000000  0.000000 )
            b(2) = (  0.000000  0.984176  0.000000 )
            b(3) = (  0.000000  0.000000  1.467614 )

   PseudoPot. # 1 for Cr read from file:
   /home/liushichen/tools/codes/espresso-4.3.2/pseudo/Cr.pz-hgh.UPF
   MD5 check sum: 1ecfd84ed12b05c881064968d61b62c6
   Pseudo is Norm-conserving, Zval =  6.0
   Generated in analytical, separable form
   Using radial grid of 1183 points,  6 beta functions with:
            l(1) = 0
            l(2) = 0
            l(3) = 0
            l(4) = 1
            l(5) = 1
            l(6) = 2

   PseudoPot. # 2 for  O read from file:
   /home/liushichen/tools/codes/espresso-4.3.2/pseudo/O.pz-hgh.UPF
   MD5 check sum: 9c3213c38e11527cd347dae851a6ebc4
   Pseudo is Norm-conserving, Zval =  6.0
   Generated in analytical, separable form
   Using radial grid of 1095 points,  1 beta functions with:
            l(1) = 0

   atomic species valence mass    pseudopotential
      Cr          6.00 52.00000    Cr( 1.00)
      O             6.00 16.00000    O( 1.00)

   Starting magnetic structure
   atomic species magnetization
      Cr          0.500
      O          0.000

   LDA+U calculation, Hubbard_lmax = 2
   atomic species  L Hubbard U  Hubbard alpha
      Cr          2 0.220496 0.000000

   16 Sym.Ops. (with inversion)

Cartesian axes

   site n.    atom                positions (alat units)
      1          Cr  tau( 1) = ( 0.0000000 0.0000000  -0.0000000  )
      2          Cr  tau( 2) = ( 0.5080391 0.5080391 0.3406890  )
      3          O tau( 3) = ( 0.3089362 0.3089362  -0.0000000  )
      4          O tau( 4) = ( 0.7071419 0.7071419 0.0000000  )
      5          O tau( 5) = ( 0.8169753 0.1991028 0.3406890  )
      6          O tau( 6) = ( 0.1991028 0.8169753 0.3406890  )
这是在vc-relax的out文件里的，上面那个的意思应该就是优化好了，然后就直接自动进行了一次自洽，我认为你应该就用优化好的结构参数把，为什么原子位置却有点出入呢？谢谢!