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

回复此楼

» 猜你喜欢

博士读完未来一定会好吗已经有6人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
心脉受损已经有3人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复

坚持就是胜利！

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

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

souledge

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

★ ★
fzx2008: 金币+2, 专家考核, 谢谢指教 2012-04-28 13:14:36

引用回帖:

10楼: Originally posted by goldenfisher at 2012-04-28 12:55:02:
所谓SC，BCC和FCC，只是一个名称，何必拘泥于形式。关键是要搞清楚其中的物理。如果即用了BCC构型，又给出了不等价原子的坐标，大家想想会出现什么情况？

这个是显然，用SC、FCC等是为了方便建模和理解，但是没有也没关系，就算选择了ibrav=14，只要celldm(i)设置合理，照样可以组成其他各种结构。这个对于建模结果是没有影响的。
我的意思是，LZ可能查资料发现他要做的结构是BCC结构了，但是如果用ibrav=1建模，得到的是SC结构，跟已知的参考结构不符～
当然了，LZ也够那啥的，既然来问是否坐标设置有问题，那就该给出参考资料中的原始结构，而不是让我们在这里猜是否建模有问题。用LZ的参数建模，得到的结构的确比较奇怪：