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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 新手求助pwscf的进行磁性运算和LDA+U运算
364/4返回列表上一页1234
查看: 4803  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 没有任何消息-是不是就凉了 +6 图啦图啦 2026-08-24 6/300 2026-08-24 17:44 by 小小红帆
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 我面上完蛋了 +8 且听虎啸 2026-08-20 9/450 2026-08-24 13:52 by zzuzxg
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 明天应该可查了！？ +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 什么时候开奖？ +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭