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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 新手求助pwscf的进行磁性运算和LDA+U运算
364/4返回列表上一页1234
查看: 4159  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[考研] 0856调剂，是学校就去 +6 sllhht 2026-03-19 7/350 2026-03-19 19:50 by 制度的
[考研] 288求调剂，一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 085600材料与化工 +5 安全上岸！ 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 复试调剂 +3 呼呼？~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
信息提示
关闭
请填处理意见
关闭