24小时热门版块排行榜    

查看: 2547  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xiaoxiao45

至尊木虫 (正式写手)

[求助] berkeleygw 输入文件问题 已有1人参与

有哪位大神指导一下berkeleygw中的WFN.in文件lattice vector怎么计算的吗?为什么kgrid.x运行后bravais lattice和crystal 的space group总是对不上?谢谢了

另外知道晶体的point group: R3c
a=b=c= 5.667
alfa=beta=gama=59.17
原子坐标 1 Bi,2 Fe,3 O
1  0.00000000  0.00000000  0.00000000
1  0.50000000  0.50000000  0.50000000
2  0.22120000  0.22120000  0.22120000
2  0.72120000  0.72120000  0.72120000
3  0.53500000  0.94000000  0.39400000
3  0.94000000  0.39400000  0.53500000
3  0.03500000  0.89400000  0.44000000
3  0.39400000  0.53500000  0.94000000
3  0.89400000  0.44000000  0.03500000
3  0.44000000  0.03500000  0.89400000
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

xiaoxiao45

至尊木虫 (正式写手)

我的数据是用的是cta Materialia 59 (2011) 4229–4234
原子分数坐标是按照Wyckoff position设置的,对称性为R3c,然后我把所有的位置数据拷贝到WFN.in中,只是不知道lattice vector (WFN.in 文件需要输入的)怎么计算,
然后用你所说的方法先算scf(pw.x),进而把计算结果带到WFN.in中计算,得到的WFN.log文件,文件中有布拉菲对称和晶体对称性,结果得不到R3c的结果,因为我感觉如果得不到所需要的结果,那么后面得到约化的K-point是不是也不对。 那么接来的wfn_q和wfn_co就无法计算了。

所以,这几天一直纠结着这问题,头发又掉不少,不知道liqizuiyang遇到过此类问题没?
求帮忙指条明路
9楼2014-10-06 18:07:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

liqizuiyang

木虫 (著名写手)

★ ★ ★ ★
franch: 金币+4, 谢谢回帖交流。。。 2014-10-06 23:36:30
1. kgrid.x中的基矢坐标和原子坐标均为在同一坐标系下的直角坐标,而不是分数坐标;
2. 如果是简单晶格,布拉伐格子对称性和晶体对称性才能对上号;复式晶格一般对不上。
2楼2014-10-05 16:34:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoxiao45

至尊木虫 (正式写手)

请问能不能给个详细的计算步骤,万分感谢,可以加送金币

[ 发自手机版 http://muchong.com/3g ]
3楼2014-10-05 23:17:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

★ ★ ★ ★ ★
franch: 金币+5, 谢谢回帖交流。。。 2014-10-06 23:36:41
如果用QE算DFT部分,有个简便方法。

先进行scf计算,输出文件中这一部分就是基矢坐标:
     crystal axes: (cart. coord. in units of alat)
               a(1) = (   4.582581   4.582581   0.000000 )  
               a(2) = (   0.000000   4.582581   4.582581 )  
               a(3) = (   4.582581   0.000000   4.582581 )  

这一部分就是同一坐标系下的原子坐标:
   Cartesian axes

     site n.     atom                  positions (alat units)
         1           O   tau(   1) = (   2.3756976   2.2068836   2.3756976  )
         2           O   tau(   2) = (   4.6669876   4.4981749   2.3756989  )
         3           O   tau(   3) = (   2.3756989   4.4981749   4.6669876  )
         4           O   tau(   4) = (   4.6669876   6.7894636   4.6669876  )
         5           O   tau(   5) = (   4.4981749   2.3756989   4.4981749  )
         6           O   tau(   6) = (   6.7894636   4.6669876   4.4981749  )
         7           O   tau(   7) = (   4.4981749   4.6669876   6.7894636  )
         8           O   tau(   8) = (   6.7894648   6.9582788   6.7894648  )
         9           Mg  tau(   9) = (   2.2912906   0.0000000   2.2912906  )
        10           Mg  tau(  10) = (   4.5825812   2.2912906   2.2912906  )
        11           Mg  tau(  11) = (   2.2912906   2.2912906   4.5825812  )
        12           Mg  tau(  12) = (   4.5825812   4.5825812   4.5825812  )
        13           Ca  tau(  13) = (   0.0000000   0.0000000   0.0000000  )
        14           Ca  tau(  14) = (   2.2912906   2.2912906   0.0000000  )
        15           Ca  tau(  15) = (   0.0000000   2.2912906   2.2912906  )
        16           Ca  tau(  16) = (   2.2912906   4.5825812   2.2912906  )
4楼2014-10-06 10:00:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:25 by mibUvS8DDCwf
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:57 by vZfe6xYu34yj
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:33 by vZfe6xYu34yj
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:00 by vZfe6xYu34yj
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:24 by vZfe6xYu34yj
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 QUjhNVAcOSff 2026-09-13 7/350 2026-09-14 21:21 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
请填处理意见