| 查看: 2335 | 回复: 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。这个问题已经纠结我两个星期了，我是新手，请高手指点啊，急死了。 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有69人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有11人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
3楼2012-04-27 15:24:26
王金荣1987
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 143.6
- 帖子: 16
- 在线: 10.6小时
- 虫号: 1434566
- 注册: 2011-10-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2012-04-27 13:07:45
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
