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

24小时热门版块排行榜    

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

刘仕晨

木虫 (正式写手)
引用回帖:
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文件里的，上面那个的意思应该就是优化好了，然后就直接自动进行了一次自洽，我认为你应该就用优化好的结构参数把，为什么原子位置却有点出入呢？谢谢!
赞一下 回复此楼
31楼2013-05-07 19:37:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
31楼: Originally posted by 刘仕晨 at 2013-05-07 19:37:45
不好意思，又要麻烦您了，前几天在学习对体系进行full-relax,用vc-relax进行计算，在优化后它会自动进行一次自洽计算，但是奇怪的是，两次的原子位置却有点不大一样，我怎么看也不明白为什么，所以要麻烦您了，我给 ...

记住，alat units的意思是以你晶格常数为单位，而在最后一步优化中，晶格常数发生了变化，你应该用最后一步的原子坐标乘以最后一步的晶格常数，看看得到的是不是自洽计算的位置。
赞一下 回复此楼
32楼2013-05-08 08:24:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)
引用回帖:
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，这是我糊涂的地方，麻烦您了，谢谢！
赞一下 回复此楼
33楼2013-05-08 13:11:16
已阅   回复此楼   关注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的回帖

刘仕晨

木虫 (正式写手)
引用回帖:
34楼: Originally posted by gemucai at 2013-05-08 14:56:14
呵呵，自洽计算给出的不是分数坐标，它的总长度不是1，而是里面的a(1)，即1.016078，所以你看看0.5080391是不是它的一半。那么化成crystal坐标是不是就变成了0.5。...

明白了，看来还是以前就没有搞清楚，谢谢，谢谢，真不知道没有您，这个我该怎么理解这些东西，谢谢了！
赞一下 回复此楼
35楼2013-05-08 20:46:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

迷迷物理

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

我想请教:我算了一个162个原子的LDA+U体系，可是显示atomic wfc used for LDA+U projector are not orthogonalized是什么原因，该怎么修改？

发自小木虫Android客户端
赞一下 回复此楼
36楼2018-04-01 09:05:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 刘仕晨 的主题更新
364/4返回列表上一页1234
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +13 一样YWY 2026-04-04 13/650 2026-04-10 11:07 by mattzhming
[考研] 材料工程085601，270求调剂 +27 @ASDF1234 2026-04-08 29/1450 2026-04-10 10:48 by 紫翼精灵
[考研] 机械专硕270求调剂，接受跨专业 +10 老师看看我吧aba 2026-04-09 11/550 2026-04-10 10:41 by cal0306
[考研] 求调剂 +3 电气300求调剂不 2026-04-08 3/150 2026-04-10 10:19 by LHGeng
[考研] 08工学 309分求调剂 +6 Yin DY 2026-04-08 6/300 2026-04-10 09:18 by Delta2012
[考研] 297求调剂 +8 Kwgyz 2026-04-09 8/400 2026-04-09 23:22 by may_新宇
[考研] 087100初试311求调剂 +3 任雅琴 2026-04-09 3/150 2026-04-09 22:42 by lbsjt
[考研] 305求调剂 +4 77Qi 2026-04-07 4/200 2026-04-09 17:27 by wp06
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 085404 293求调剂 +7 勇远库爱314 2026-04-08 7/350 2026-04-09 16:02 by 猪会飞
[考研] 0860004 求调剂 309分 +7 Yin DY 2026-04-08 7/350 2026-04-09 14:06 by ditto77778
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +11 Mr. Z 2026-04-05 12/600 2026-04-08 16:15 by luoyongfeng
[考研] 277求调剂 +4 考研调剂lxh 2026-04-06 6/300 2026-04-08 10:40 by 逆水乘风
[考研] 调剂 +3 电气300求调剂不 2026-04-08 6/300 2026-04-08 09:39 by 电气300求调剂不
[考研] 319分085702安全工程求调剂 +6 rious 2026-04-05 6/300 2026-04-07 09:42 by jp9609
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 考研调剂 +7 15615482637 2026-04-04 7/350 2026-04-06 22:56 by chenzhimin
[考研] 一志愿北交大材料工程总分358求调剂 +10 cs0106 2026-04-05 12/600 2026-04-06 19:41 by 无际的草原
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 本科211，专业085404，293分请求调剂 +5 莲菜就是藕吧 2026-04-04 5/250 2026-04-04 14:08 by 这是一个无聊的�
信息提示
关闭
请填处理意见
关闭