24小时热门版块排行榜

返回列表

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

甲斐之虎

银虫 (正式写手)

应助: 60 (初中生)
金币: 302.3
散金: 220
红花: 7
帖子: 970
在线: 465小时
虫号: 2273045
注册: 2013-02-03
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 关于晶格常数弛豫，一直有个疑问已有3人参与

比如，对于石墨烯，假设PBE泛函能量最小那点为晶格常数=2.45A
按照，侯博士那个测试脚本可以进行测试，比如设置i = 2.44 2.45 2.46 2.47 2.48……
但是，如果不那么做，而是一开始建立模型取实验值a=2.46A，然后对石墨烯进行全弛豫（ISIF=3），是否一样可以得到a=2.45A这个结果呢？
另外，如果这个方法可行，那么这个方法对于其他二维材料是否也同样可行呢？
不知道有没有表达清楚，请各位前辈不吝赐教，先行谢过了

附侯博士测试脚本：
#!/bin/sh
rm WAVECAR
for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30
do
cat > POSCAR <<!
Al-fcc
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1
Direct
0.0 0.0 0.0
!
echo / a = $i angstrom /; time vasp
E=‘grep /TOTEN0 OUTCAR | tail -1 | awk -{printf /%12.6f \n0, $5 }.‘
V=‘grep /volume0 OUTCAR | tail -1 | awk -{printf /%12.4f \n0 , $5}.‘
echo $V $E >>EtVo.dat
done

回复此楼