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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 压力下的电子结构及弹性计算
51/1返回列表
查看: 925  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

x夏日的约定

铜虫 (小有名气)

[求助] 压力下的电子结构及弹性计算

小弟有个问题想请教各位虫友，请问计算压力下的电子结构和弹性时，在晶体结构优化完后，进行态密度、能带和弹性性质计算时，要不要加压力？希望获得各位虫友的回答，谢谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-08-24 19:00:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

x夏日的约定

铜虫 (小有名气)
引用回帖:
7楼: Originally posted by xueht987 at 2012-08-26 09:37:47
这个可能是你应变为零时的结构没有弛豫好，重新弛豫，之后在此基础上逐步加压或者应变。...

可是我仔细对比我的弹性计算结果，在没有应变的情况下计算出的体系自由能与在计算弹性之前的结构优化后的体系自由能的数据是相同的，这又是怎么回事？下面是我的计算文件设置，麻烦大家帮我看一下：
（1）在压力15GPa下晶格优化,得到：
free  energy   TOTEN  =       -27.116931 eV
energy  without entropy= -27.116931  energy(sigma->0) =  -27.116931
（2）弹性常数C11、C12计算
script.sh文件参数：
ENCUT=520
ISTART=0
ICHARG=2
ISMEAR=0;SIGMA=0.2
NSW=80;IBRION=2
POTIM=0.2
IALGO=48
EDIFF=5E-06
ISIF=2
PREC=Normal
LWAVE=.FALSE.
LCHARG=.FALSE.
!

echo "delta=$i"; mpirun -np 2 vasp5.2

cp CONTCAR pos.$i
cp CONTCAR POSCAR

cat > INCAR < SYSTEM=
ENCUT=520
ISTART=0
ICHARG=2
ISMEAR=0
IALGO=48
EDIFF=5E-6
PREC=Normal
LWAVE=.FALSE.
LCHARG=.FALSE.

计算结果：
-0.018 1 F= -.26572346E+02 E0= -.26572346E+02
-0.015 1 F= -.26674069E+02 E0= -.26674069E+02
-0.012 1 F= -.26770758E+02 E0= -.26770758E+02
-0.009 1 F= -.26863577E+02 E0= -.26863577E+02
-0.006 1 F= -.26952164E+02 E0= -.26952164E+02
-0.003 1 F= -.27036589E+02 E0= -.27036589E+02
0.00   1 F= -.27116931E+02 E0= -.27116931E+02
0.003  1 F= -.27193287E+02 E0= -.27193287E+02
0.006  1 F= -.27265715E+02 E0= -.27265715E+02
当应变超过0时，体系总能仍在下降，这是为什么？明明在应变为0的时候，体系的总能和优化的总能相一致，可是随着应变的变化，体系的总能却仍然降低。
赞一下 回复此楼
9楼2012-08-26 14:24:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

xueht987

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
uuv2010: 金币+2, 多谢详细指导 2012-09-29 10:26:54
结构优化完之后再算性质的话，可以不用在输入文件加压了，因为性质计算的时候结构没有变化，也即此时算的就是处于压力下结构的性质。不放心的话，可以找个简单的例子试一试，加与不加应该不会有很大差别。
赞一下 回复此楼
2楼2012-08-25 21:41:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hwceng0816

木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
这是需要的，几何优化后的计算和加压后的计算式有区别的。
赞一下 回复此楼
3楼2012-08-26 06:51:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

x夏日的约定

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by xueht987 at 2012-08-25 21:41:44
结构优化完之后再算性质的话，可以不用在输入文件加压了，因为性质计算的时候结构没有变化，也即此时算的就是处于压力下结构的性质。不放心的话，可以找个简单的例子试一试，加与不加应该不会有很大差别。

哦，但是我计算C11和C12的时候出现一个情况就是在我给定的应变下，体系的总能不断的减小，而并不是在应变为零时达到最低点，然后再增加。我尝试将K点增大，但结果还是不好，这是怎么回事呢？
赞一下 回复此楼
4楼2012-08-26 08:52:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭