24小时热门版块排行榜

返回列表

» 猜你喜欢

博士读完未来一定会好吗已经有27人回复
投稿精细化工已经有4人回复
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！已经有4人回复
导师想让我从独立一作变成了共一第一已经有9人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:

求晶格常数优化小程序（a，b固定，只优化c轴）已经有9人回复
求教，为何我优化出石墨烯的晶格常数是2.4683，而文献里都是2.47？？？已经有25人回复
优化后晶格常数变化很大原因已经有10人回复
VASP做结构优化时用ISIF=3只想对晶格常数中的一个维度放开经行优化已经有7人回复
结构优化对于计算精度到底有多大影响已经有6人回复
计算锗烯的晶格常数和起皱高度电场E=0.01 已经有4人回复
无法优化出文献中的晶格常数已经有7人回复
声子谱虚频求助已经有10人回复
vasp计算三角晶系转换成primitive cell优化后，晶格常数和实验值怎么比较？急求已经有4人回复
求教各位大侠，优化石墨烯的晶格常数时，IBRION与ISIF是如何设置的呢？已经有4人回复
优化完成，发现晶格常数跟实验的数据相差太大已经有8人回复
Castep优化结果的晶格常数最终结果怎么看？已经有13人回复
VASP结构优化ISIF=3，为什么晶格常数不变已经有12人回复
怎么确定气体如（氩气）固体的晶格常数已经有5人回复
求石墨烯介电常数已经有6人回复
用LDA赝势优化出来的晶格常数一定比实验值小吗？谢谢！已经有10人回复
【求助】castep 怎么预测晶格常数已经有8人回复
【求助】怎么计算镁合金的晶格常数？已经有7人回复

1楼2014-06-26 19:28:43

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

liqizuiyang

木虫 (著名写手)

应助: 59 (初中生)
金币: 6859.5
散金: 261
红花: 94
帖子: 2093
在线: 708.9小时
虫号: 2366686
注册: 2013-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
chuanghua304: 金币+15, ★★★很有帮助, 谢谢，真速度。 2014-06-26 22:47:35
franch: 金币+2, 谢谢回帖交流。。。 2014-07-08 09:49:18

#! /bin/bash

# ncores
NCORE=24

# dft flavor
dft="pz-vbc"

# celldm related vars
alist="2.20 2.25 2.30 2.35 2.40 2.45 2.50 2.55 2.60 2.65 2.70 2.75"
a2b="1.889726124546135"
vacc="18.0"

for a in $alist ; do
celldm1=$(awk 'BEGIN{printf("%11.9f", '$a' * '$a2b')}')
celldm3=$(awk 'BEGIN{printf("%11.9f", '$vacc' / '$a'  )}')
cat > $a.in << EOF
&CONTROL
  prefix          = '$a'
  calculation    = 'relax'
  restart_mode    = 'from_scratch'
  pseudo_dir    = './'
  outdir          = './'
  wfcdir          = './'
  verbosity       = 'high'
  wf_collect    = .TRUE.
  nstep          = 250
  etot_conv_thr = 7.3d-5
  forc_conv_thr = 3.8d-4
  tprnfor       = .TRUE.
/
&SYSTEM
  ibrav          = 4
  celldm(1)       = $celldm1
  celldm(3)       = $celldm3
  ntyp          = 1
  nat             = 2
  nbnd          = 8
  occupations    = 'smearing'
  smearing       = 'mp'
  degauss       = 0.015
  ecutwfc       = 60.0
/
&ELECTRONS
  electron_maxstep = 250
  conv_thr       = 1.0d-8
  mixing_mode    = 'plain'
  mixing_beta    = 0.7
  mixing_ndim    = 8
  diagonalization  = 'david'
  diago_david_ndim = 4
/
&IONS
  ion_dynamics    = 'bfgs'
  bfgs_ndim       = 1
/
ATOMIC_SPECIES
C  12.011 C.$dft.UPF
ATOMIC_POSITIONS {crystal}
C 0.000000000 0.000000000 0.500000000
C 0.666666667 0.333333333 0.500000000
K_POINTS {automatic}
12 12 1 0 0 0
EOF
mpirun -np $NCORE pw.x -in $a.in > $a.out
wait
rm -r $a.save
done