| 查看: 1030 | 回复: 10 | ||
[求助]
压力下的电子结构及弹性计算
|
| 小弟有个问题想请教各位虫友，请问计算压力下的电子结构和弹性时，在晶体结构优化完后，进行态密度、能带和弹性性质计算时，要不要加压力？希望获得各位虫友的回答，谢谢！ |
回复此楼» 猜你喜欢
0831生医工第一轮调剂失败求助
已经有17人回复
-
291求调剂
已经有11人回复
-
297，工科调剂?河南农业大学本科
已经有14人回复
-
0854求调剂
已经有20人回复
-
26药学专硕105500求调剂
已经有7人回复
-
一志愿华中农业071010，320求调剂
已经有15人回复
-
2026 WR青拔
已经有7人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有9人回复
-
271求调剂
已经有37人回复
-
申博自荐
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 厦大研究生招生1、化学软件2、电子结构理论与计算方法研究方向的区别
已经有22人回复
- 【求助】VASP中，“结构优化”和“电子自洽计算”的关系？
已经有8人回复
- 【讨论】电子结构能量E与自由能G，如何判断化学反应走向
已经有16人回复
- 【求助】非零压下弹性常数计算的初始结构
已经有9人回复
- 【求助】Ecut, k点数等设置需随压力变化吗?
已经有10人回复
- 管理理论每日谈更新
已经有65人回复
- 纱线毛羽cv(%)值是什么意思？
已经有2人回复
xueht987
木虫 (正式写手)
- 应助: 112 (高中生)
- 贵宾: 0.081
- 金币: 3514.2
- 散金: 106
- 红花: 6
- 帖子: 593
- 在线: 273.4小时
- 虫号: 1059317
- 注册: 2010-07-17
- 性别: GG
- 专业: 金属材料表面科学与工程
2楼2012-08-25 21:41:44
hwceng0816
木虫 (著名写手)
- 应助: 130 (高中生)
- 金币: 4121
- 散金: 244
- 红花: 46
- 帖子: 1914
- 在线: 1123.9小时
- 虫号: 1529751
- 注册: 2011-12-08
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2012-08-26 06:51:32
4楼2012-08-26 08:52:56
5楼2012-08-26 08:56:03
6楼2012-08-26 09:21:18
xueht987
木虫 (正式写手)
- 应助: 112 (高中生)
- 贵宾: 0.081
- 金币: 3514.2
- 散金: 106
- 红花: 6
- 帖子: 593
- 在线: 273.4小时
- 虫号: 1059317
- 注册: 2010-07-17
- 性别: GG
- 专业: 金属材料表面科学与工程
7楼2012-08-26 09:37:47
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
8楼2012-08-26 10:18:15
|
可是我仔细对比我的弹性计算结果，在没有应变的情况下计算出的体系自由能与在计算弹性之前的结构优化后的体系自由能的数据是相同的，这又是怎么回事？下面是我的计算文件设置，麻烦大家帮我看一下： （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
10楼2012-08-26 14:25:46
