| 查看: 4755 | 回复: 25 | ||||
[交流]
【求助】一定压强下的弹性常数计算方法 已有10人参与
|
|
我计算的思路是这样的： 1。改变晶胞的体积，由BM方程拟合出该体积下的压强。 2。选出拟合出来的某个压强下的晶胞体积，此体积的晶胞对应某个压强，然后按侯博士的方法，施加应变，拟合二次曲线。 可是给一定压强下的晶胞，施加应变后，能量-应变关系基本就是一条单调递减的直线，指向平衡体积的方向。能量选用的是VASP输出的能量，没有考虑压强进去。 请问各位在计算高压下的弹性常数时，是怎么操作的呢？谢谢大家了。 |
回复此楼» 收录本帖的淘帖专辑推荐
 vasp压力 |
» 猜你喜欢
回收溶剂求助
已经有7人回复
-
职称评审没过，求安慰
已经有40人回复
-
硝基苯如何除去
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 弹性常数计算出现诡异。。。。来看看喽
已经有5人回复
- 加压后的体系怎么计算弹性常数？（或其它性质）
已经有6人回复
- 【求助】CASTEP计算弹性常数得到应力应变曲线
已经有17人回复
- 【求助】非零压下弹性常数计算的初始结构
已经有9人回复
- 【求助】弹性常数与压强
已经有8人回复
- 【求助】MS弹性常数计算分析问题
已经有5人回复
- 【求助】CASTEP服务器弹性常数的计算
已经有18人回复
- 【讨论】用vasp计算过弹性常数的请进！！更新中！！
已经有51人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★
xh512(金币+10):谢谢侯老师 2010-10-26 21:11:17
zzy870720z(金币+1):谢谢指教 2010-11-15 12:05:50
xh512(金币+10):谢谢侯老师 2010-10-26 21:11:17
zzy870720z(金币+1):谢谢指教 2010-11-15 12:05:50
|
参考这篇文献： J. Phys.: Condens. Matter 21 (2009) 115401 http://iopscience.iop.org/0953-8984/21/11/115401/ |
2楼2010-10-26 20:56:16
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指教 2010-11-13 14:47:28
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指教 2010-11-13 14:47:28
|
最近也一直学这个, 得到多位虫友帮助, 思路如下: (1)高精度ISIF=3优化结构,压强通过PSTRESS设定 (2)根据材料的对称性, 选定几个应变,比如 (100000), 等,每个选定应变方式下取6-8个不同幅度应变值优化原子位置然后算静态能(参考侯博士教程) (3)从静态能结果OUTCAR里找出各应变幅度对应的应力值,根据材料对称性, 决定每个应力对应的应变, 比如说(1000000) 中F1=C11*应变, 等等. $>grep "in kB" OUTCAR 会显示 F1, F2, F3, F4, F5 F6 (单位kbar ) (4)对6-8个不同幅度应变值下的得到的不同 F-Cij 进行拟和 就可以得到Cij |
3楼2010-11-13 11:59:00
sg18408926
至尊木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 10945.9
- 散金: 246
- 红花: 1
- 帖子: 1075
- 在线: 431.4小时
- 虫号: 727956
- 注册: 2009-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2010-11-15 11:57:34
5楼2010-11-15 14:47:31
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
6楼2010-11-16 09:10:59
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
7楼2010-11-16 09:45:59
8楼2010-11-17 23:15:56
fengshiquan
铁杆木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 5830.3
- 散金: 80
- 红花: 5
- 帖子: 599
- 在线: 1569.1小时
- 虫号: 884725
- 注册: 2009-10-26
- 专业: 原子和分子物理
9楼2011-03-31 10:14:22
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
10楼2011-04-01 07:29:37