| 查看: 2327 | 回复: 15 | ||||
王金荣1987铁虫 (初入文坛)
|
[求助]
关于PWscf中的位置坐标设置问题
|
|
我在用PWscf计算一个氮化物的时的脚本是这样的 cd /home/ghz/wjr/cscl for i in `seq -w 4.11 0.1 9.31` do cd /home/ghz/wjr/cscl mkdir $i cd $i cat > RhN.in < calculation = 'scf' restart_mode='from_scratch', prefix='RhN', tstress = .true. tprnfor = .true. forc_conv_thr=1.0D-4 etot_conv_thr=1.0D-6 pseudo_dir = '/home/ghz/wjr/espresso-4.3-examples/pseudo/', outdir='/home/ghz/wjr/tmp/' / &system ibrav= 3, celldm(1) =$i, nat= 2, ntyp= 2,ecutwfc=100, occupations='smearing', smearing='marzari-vanderbilt', degauss=0.05 / &electrons diagonalization='david' mixing_mode = 'plain' mixing_beta = 0.7 conv_thr = 1.0d-8 / ATOMIC_SPECIES Rh 102.906 Rh.pbe-rrkjus_lb.UPF N 14.007 N.pbe-paw_kj.UPF ATOMIC_POSITIONS crystal Rh 0.00 0.00 0.00 N 0.50 0.50 0.50 K_POINTS automatic 24 24 24 0 0 0 EOF /home/ghz/wjr/espresso-4.3/bin/pw.x sed -n '/! total energy/p' RhNk=$i.out | awk '{print $5}' >> ../total_energy.txt sed -n '/unit-cell volume =/p' RhNk=$i.out | awk '{print $4}' >> ../cell_volume.txt cd .. done 但是算出来的结果用EOS拟合之后得到的体模量，晶格常数和文献上别人用其他软件算的差了很远。我怀疑是我的位置坐标设的不对，算的这个晶格类型是体心立方，ibrav=3。这个问题已经纠结我两个星期了，我是新手，请高手指点啊，急死了。 |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有12人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- PWSCF中如何填写基本的参数
已经有9人回复
- PWscf计算homo lumo或者valence conduction band，问题已解决
已经有14人回复
- PWscf中CP的温度单位
已经有5人回复
- 【求助】pwscf计算中的warning要紧吗？如何消除呢？
已经有9人回复
- 【原创】PWSCF中的nbnd和degauss
已经有11人回复
- 【求助】pwscf中的问题
已经有4人回复
- 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算
已经有23人回复
- 【求助】PWSCF输出的projwfc文件中，哪些文件是dxy、dxz、dyz轨道？
已经有10人回复
王金荣1987
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 143.6
- 帖子: 16
- 在线: 10.6小时
- 虫号: 1434566
- 注册: 2011-10-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2012-04-27 13:07:45
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
3楼2012-04-27 15:24:26
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
你不觉得LZ的赝势用得蛮乱么？第一个应该是rrkj形式的超软，第二个则是PAW～～～～ 4楼2012-04-27 19:46:26
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
5楼2012-04-28 09:40:38
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
6楼2012-04-28 12:03:07
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
7楼2012-04-28 12:09:14
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
★ ★
fzx2008: 金币+2, 专家考核, 谢谢回帖！ 2012-04-28 13:08:12
fzx2008: 金币+2, 专家考核, 谢谢回帖！ 2012-04-28 13:08:12
|
话是这么说，手册中也很清楚得写着不同类型的赝势可以混用，比如NCPP和USPP混用之类的。不过试了几次，总觉得怪怪的，大概可能因为用的NCPP是GTH的，属于解析赝势，而USPP是Vanderbilt的…… 不过，话说，如果用了不同XC泛函的，然后用input_dft参数强制XC泛函的话可以不可以呢？看着CASTEP似乎就那么个玩法，整个GGA的赝势就一套赝势XX_pbe00.ups，然后PBE，PW91，WC等泛函都是外加的…… 顺便，感觉这个模型有点纠结，如果是ibrav=1，那就不是BCC结构了，而是SC……因为原点和体心的原子类型不同，不符合BCC的对称性……搜了ICSD，没看到有BCC结构的RhN，只见到了四方结构的RhN2…… |
8楼2012-04-28 12:34:32
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
9楼2012-04-28 12:50:13
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
10楼2012-04-28 12:55:02
