| 查看: 1959 | 回复: 7 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
王金荣1987铁虫 (初入文坛)
|
[求助]
如何用PWSCF计算弹性常数啊？
|
|||
我用的是将基矢的微形变，放到脚本里面算出不同形变下的能量，再用origin和形变参数伽马进行拟合得到能量关于伽马的函数的泰勒展开式的二阶系数B2,在用下面图中的公式算出弹性常数。但是怎么算都不对啊，恶心死了，这个东西整了一两个月了。就是整不对啊！！！！ 要是那个高手能帮我解决了这个问题我把所有的积分都给他，我都愿意的。  弹性常数计算1.jpg [ Last edited by fzx2008 on 2012-12-4 at 19:35 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 good resources |
» 猜你喜欢
一志愿A区211，22408 321求调剂
已经有8人回复
-
295分求调剂
已经有6人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有13人回复
-
085600材料与化工调剂
已经有6人回复
-
期刊推荐
已经有5人回复
-
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业
已经有7人回复
-
又一批高校组建人工智能学院 师资行吗 不是骗人吗
已经有4人回复
-
337求调剂
已经有3人回复
-
申博/考博
已经有5人回复
-
申博
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求专家解惑：pwscf并行计算cpu分配问题
已经有3人回复
- 新手利用pwscf 优化 ZnO的晶格常数老是出问题， 求助！
已经有9人回复
- 求PWSCF使用说明书
已经有3人回复
- 【求助】pwscf计算得到的能带图很乱
已经有16人回复
- 【求助】利用abinit或pwscf软件在单机上可以计算半导体材料的哪些性质?
已经有8人回复
- 【求助】求Fortan90或Fortran95编译器，安装PWscf用的。
已经有6人回复
- 【求助】关于PWSCF介电常数和声子谱的计算
已经有7人回复
- 【求助】PWSCF计算介电极化时为何出现如下提示？
已经有19人回复
- 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算
已经有23人回复
- 【求助】pwscf怎样计算介电常数和波恩有效电荷
已经有13人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
- 【求助】pwscf能计算材料的一些宏观性质吗？
已经有9人回复
王金荣1987
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 143.6
- 帖子: 16
- 在线: 10.6小时
- 虫号: 1434566
- 注册: 2011-10-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
送红花一朵 
8楼2013-04-23 09:50:59
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
franch: 金币+5, 专家考核, 谢谢回帖交流， 2012-12-03 21:46:08
王金荣1987: 金币+30, ★有帮助, 这个我大致有点明白，但是具体的还是有点模糊，那个M J Mehl的文章我确实没有看，我是学着VASP的做法做的，你能把那个文献给我一下吗？ 2012-12-04 12:31:36
感谢参与，应助指数 +1
franch: 金币+5, 专家考核, 谢谢回帖交流， 2012-12-03 21:46:08
王金荣1987: 金币+30, ★有帮助, 这个我大致有点明白，但是具体的还是有点模糊，那个M J Mehl的文章我确实没有看，我是学着VASP的做法做的，你能把那个文献给我一下吗？ 2012-12-04 12:31:36
|
咋看LZ的计算方法这么纠结呢？有没有看M J Mehl的文章呢？对于立方晶（看LZ只需要C11，C12和C44），只需要计算三步： 1，使用EOS拟合计算体模量：B0 = （ C11 +2C12 ）／ 3； 2，假设形变为x，则有形变矩阵：e1=x, e2=-x, e3=x^2/(1-x^2)，其他应变量为零。再把x-E（x）关系拟合方程dE=V0 （C11-C12）x^2，其中V0是平衡位置体积，dE是形变后总能与平衡位置总能的差。然后联立该结果和体模量得到C11和C12。 3，对于C44，进行形变为：e3=x^2/(4-x^2), e6=x，其他向量为零。拟合方程 dE = 0.5 V0 C44 x^2。 |
» 本帖已获得的红花（最新10朵）
2楼2012-12-03 20:19:04
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
3楼2012-12-04 05:15:59
zhoufeng1982
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 2431.1
- 散金: 1
- 红花: 3
- 帖子: 2691
- 在线: 82.6小时
- 虫号: 1193591
- 注册: 2011-01-21
- 专业: 物理化学
4楼2012-12-04 08:21:53
5