24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

王金荣1987

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 143.6
帖子: 16
在线: 10.6小时
虫号: 1434566
注册: 2011-10-10
性别: GG
专业: 凝聚态物性I:结构、力学和

[求助] 关于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。这个问题已经纠结我两个星期了，我是新手，请高手指点啊，急死了。

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有5人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有11人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有16人回复
初秋的晨风已经有3人回复
面上没中，邀请各位路过的虫友分析一下分数已经有11人回复
Ei源刊怎么投已经有4人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有6人回复
两块石头已经有5人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复

坚持就是胜利！

1楼 2012-04-25 21:43:06

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

gemucai

木虫 (正式写手)

应助: 16 (小学生)
金币: 1774.6
散金: 157
红花: 4
帖子: 693
在线: 201.7小时
虫号: 816478
注册: 2009-07-27
专业: 凝聚态物性 II ：电子结构

引用回帖:

15楼: Originally posted by souledge at 2012-04-28 14:58:51:
Si是金刚石结构，应该是
ibrav = 2,
Si1 0 0 0
Si2 1/4 1/4 1/4
的才是～
一个顶点一个体心的是NaCl结构，也叫做B1结构，所以NaCl的建模是：
ibrav = 2,
Na 0 0 0
Cl 1/2 1/2 1/2

时间长了不记得了