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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 压力下的电子结构及弹性计算
111/2返回列表12下一页
查看: 915  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

x夏日的约定

铜虫 (小有名气)

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

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

» 猜你喜欢

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

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的回帖

x夏日的约定

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by hwceng0816 at 2012-08-26 06:51:32
这是需要的，几何优化后的计算和加压后的计算式有区别的。

嗯，谢谢，我也计算过在弹性计算时，在script.sh中输入压力参数后，通过对比分析无应变时的能量，发现体系的能量与在计算压力下晶格参数优化后的自由能有很大差距。你能提供一些关于加压后的计算公式吗？谢谢！
赞一下 回复此楼
5楼2012-08-26 08:56:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

x夏日的约定

铜虫 (小有名气)
我还想请教一个问题：就是计算弹性的时候，自由能和应变之间的关系可不可以用公式转换成自由能与应变对应的压力的公式？如果有，可以提供一点资料给我吗？谢谢了
赞一下 回复此楼
6楼2012-08-26 09:21:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xueht987

木虫 (正式写手)

【答案】应助回帖

★
uuv2010: 金币+1, 多谢补充 2012-09-29 10:27:04
引用回帖:
4楼: Originally posted by x夏日的约定 at 2012-08-26 08:52:56
哦，但是我计算C11和C12的时候出现一个情况就是在我给定的应变下，体系的总能不断的减小，而并不是在应变为零时达到最低点，然后再增加。我尝试将K点增大，但结果还是不好，这是怎么回事呢？...

这个可能是你应变为零时的结构没有弛豫好，重新弛豫，之后在此基础上逐步加压或者应变。
赞一下 回复此楼
7楼2012-08-26 09:37:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★
xueht987: 金币+1, 是的 2012-08-27 09:55:55
引用回帖:
7楼: Originally posted by xueht987 at 2012-08-26 09:37:47
这个可能是你应变为零时的结构没有弛豫好，重新弛豫，之后在此基础上逐步加压或者应变。...

压力下的性质计算，需要先在加压下做结果优化吗？
赞一下 回复此楼
8楼2012-08-26 10:18:15
已阅   回复此楼   关注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的回帖

x夏日的约定

铜虫 (小有名气)
引用回帖:
8楼: Originally posted by cenwanglai at 2012-08-26 10:18:15
压力下的性质计算，需要先在加压下做结果优化吗？...

因为我想了解随着压力的变化，体系的弹性性质的变化
赞一下 回复此楼
10楼2012-08-26 14:25:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 x夏日的约定 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭