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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 关于PWscf中的位置坐标设置问题
51/1返回列表
查看: 1888  |  回复: 15
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

王金荣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 < &control
    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 RhNk=$i.out
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。这个问题已经纠结我两个星期了，我是新手，请高手指点啊，急死了。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
坚持就是胜利！
1楼2012-04-25 21:43:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
4楼: Originally posted by souledge at 2012-04-27 19:46:26:
你不觉得LZ的赝势用得蛮乱么？第一个应该是rrkj形式的超软，第二个则是PAW～～～～

对啊，能不能混用也是个问题。
还有一点，超软赝势要设置ecutrho=1000，如果默认的话就太小了。然而另一个PAW赝势，不太清楚要设置为多少。。。
赞一下 回复此楼
5楼2012-04-28 09:40:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

王金荣1987

铁虫 (初入文坛)
跪求各位大侠相救啊。
赞一下 回复此楼
坚持就是胜利！
2楼2012-04-27 13:07:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
★ ★ ★
fzx2008: 金币+1, 鼓励交流 2012-04-27 16:19:23
WDD880227: 金币+2, 感谢交流给出建议！ 2012-04-28 08:58:11
我不是高手，只能给一点小建议：
1.鉴于你的能量截断和K点已经足够高，检查赝势。
2.你用GGA与别人用LDA比的话软很多是正常的。
3.暂时没想到
赞一下 回复此楼
3楼2012-04-27 15:24:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★
fzx2008: 金币+1, 专家考核, 谢谢回帖！ 2012-04-28 13:07:24
引用回帖:
3楼: Originally posted by gemucai at 2012-04-27 15:24:26:
我不是高手，只能给一点小建议：
1.鉴于你的能量截断和K点已经足够高，检查赝势。
2.你用GGA与别人用LDA比的话软很多是正常的。
3.暂时没想到

你不觉得LZ的赝势用得蛮乱么？第一个应该是rrkj形式的超软，第二个则是PAW～～～～
赞一下 回复此楼
思想重于技巧，内涵重于表象
4楼2012-04-27 19:46:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭