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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 压力下的电子结构及弹性计算
51/1返回列表
查看: 998  |  回复: 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289求调剂 +15 硕星赴 2026-03-23 15/750 2026-03-26 02:16 by BruceLiu320
[考研] 一志愿上海交大生物与医药专硕324分，求调剂 +6 jiajunX 2026-03-22 6/300 2026-03-25 23:05 by licg0208
[考研] 308求调剂 +5 墨墨漠 2026-03-25 5/250 2026-03-25 22:19 by 544594351
[考研] 网络空间安全0839招调剂 +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 303求调剂 +6 元夕元 2026-03-20 7/350 2026-03-25 12:00 by edmund7
[考研] B区考研调剂 +4 yqdszhdap－ 2026-03-22 5/250 2026-03-25 08:51 by baoball
[考研] 289材料与化工（085600）B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 求调剂 +5 林之夕 2026-03-24 5/250 2026-03-24 17:16 by dick_runner
[考研] 279分求调剂 一志愿211 +18 chaojifeixia 2026-03-19 20/1000 2026-03-24 10:34 by dolphin_ycj
[考研] 291求调剂 +8 hhhhxn.. 2026-03-23 8/400 2026-03-23 23:15 by peike
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 一志愿070300浙大化学358分，求调剂！ +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
信息提示
关闭
请填处理意见
关闭