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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 用侯博士计算弹性常数的方法，计算NITI B2相的C11和C12出现了问题。求助高手！
71/1返回列表
查看: 2239  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

losenq

新虫 (小有名气)

[求助] 用侯博士计算弹性常数的方法，计算NITI B2相的C11和C12出现了问题。求助高手！

我计算的Niti B2结构，用了两种方法来计算弹性常数，一种是通过VASP直接计算，IBRION=6， ISIF=3， NFREE=4

计算出来C11=194GPA，C12=178GPA，C44=53GPA。与其他参考文献对照了一下，发现C11和C12相差稍远，C44倒是有很接近的。

于是又用侯博士的那种方法，具体方法小木虫上面应该有挺多帖子。DELTA变化范围从-0.02到0.02，间隔0.0025。

第一步算C44，E=（0，0，0，DELTA，DELTA，DELTA）
得到的结果是49GPA，这与参考文献相差不多，参考值46GPA。

第二步算C11，C12。
分别用E=（DELTA，DELTA，0，0，0，0）
和
E=（DELTA，DELTA，DELTA，0，0，0）
计算（E-E0）/V0 ~DELTA，然后进行拟和，联合求解计算C11，C12，发现这两个值居然相当接近！
C11=159GPA
C12=160GPA

这与文献相当不符合阿。文献里面一般是C11-C12〉0

我检查了好几遍觉得程序没什么问题，就用的侯博士的程序~~

现在不知道该怎么办了。。。。求助各位高手阿！！！


请版主代扣金币吧。。。我还不太熟悉怎么弄啊。。。谢谢了。。。
回复此楼

» 猜你喜欢

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

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+3, 专家考核, 谢谢指教！ 2012-09-07 16:37:28
losenq: 金币+5, ★★★很有帮助 2012-10-14 20:07:17
对于立方体系，为什么不试试M. J. Mehl的方法？先做EOS拟合，得到体模量，且体模量 B0 = (C11 + 2*C12) / 3。
然后做C11-C12的拟合，拟合方程，DeltaE = (C11-C12)*V0*d^2
其中d表示应变，对应的应变矩阵为：
d  0  0
0  -d  0
0  0  d^2/(1-d^2)
然后联立体模量求解C11和C12。
其实LZ的方法中，
d  0  0
0  d  0
0  0  d
的应变也是求解体模量的应变了，但是拟合的是简单的能量响应，而不是EOS，所以可能结果上会有点出入。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

losenq losenq
思想重于技巧，内涵重于表象
2楼2012-09-07 15:07:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)
送鲜花一朵
引用回帖:
2楼: Originally posted by souledge at 2012-09-07 15:07:08
对于立方体系，为什么不试试M. J. Mehl的方法？先做EOS拟合，得到体模量，且体模量 B0 = (C11 + 2*C12) / 3。
然后做C11-C12的拟合，拟合方程，DeltaE = (C11-C12)*V0*d^2
其中d表示应变，对应的应变矩阵为：
d  ...

十分感谢！那我现在用这种方法试一试。
赞一下 回复此楼
3楼2012-09-10 09:24:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)
引用回帖:
2楼: Originally posted by souledge at 2012-09-07 15:07:08
对于立方体系，为什么不试试M. J. Mehl的方法？先做EOS拟合，得到体模量，且体模量 B0 = (C11 + 2*C12) / 3。
然后做C11-C12的拟合，拟合方程，DeltaE = (C11-C12)*V0*d^2
其中d表示应变，对应的应变矩阵为：
d  ...

您好，我在做C11-C12的时候，在文件defvect.f中定义strain(3)=delta^2/(1-delta^2)后，在gfortran -o defvect.x defvect.f时遇到问题：

错误：（1）语句无法归类。

请问在defvect.f中是不是不能进行计算运行，还是说语句、符号有错误？

谢谢了~
赞一下 回复此楼
4楼2012-09-12 15:28:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★
sunyang1988: 金币+1, 专家考核, 谢谢交流 2012-09-13 10:39:05
引用回帖:
4楼: Originally posted by losenq at 2012-09-12 15:28:37
您好，我在做C11-C12的时候，在文件defvect.f中定义strain(3)=delta^2/(1-delta^2)后，在gfortran -o defvect.x defvect.f时遇到问题：

错误：（1）语句无法归类。

请问在defvect.f中是不是不能进行计算运行 ...

Fortran中，N次方的计算用x**N表示，而不是常用的^表示。
顺便，如果不是很熟悉Fortran，为了防止在代码中出现自己看不到的问题，强烈建议还是先手动修改晶格进行计算。
不知道哪个defvect.f是个怎样的代码，其实用VASP求解弹性常数，自己写脚本（然后Origin之类的软件拟合），或者直接用VASP提供的直接法求，都是不错的选择。
赞一下 回复此楼
思想重于技巧，内涵重于表象
5楼2012-09-12 18:48:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)
引用回帖:
5楼: Originally posted by souledge at 2012-09-12 18:48:03
Fortran中，N次方的计算用x**N表示，而不是常用的^表示。
顺便，如果不是很熟悉Fortran，为了防止在代码中出现自己看不到的问题，强烈建议还是先手动修改晶格进行计算。
不知道哪个defvect.f是个怎样的代码，其实 ...

用VASP直接算得的弹性常数与参考文献给出值相差较大。。。。所以后来才改用这种方法算的。嗯，发现了，确实是FORTRAN里二次方是**不是^。谢谢啦！
赞一下 回复此楼
6楼2012-09-13 19:35:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)
送鲜花一朵
引用回帖:
2楼: Originally posted by souledge at 2012-09-07 15:07:08
对于立方体系，为什么不试试M. J. Mehl的方法？先做EOS拟合，得到体模量，且体模量 B0 = (C11 + 2*C12) / 3。
然后做C11-C12的拟合，拟合方程，DeltaE = (C11-C12)*V0*d^2
其中d表示应变，对应的应变矩阵为：
d  ...

您好，我使用了您说的这种方法，拟合\delta E=(c11-c12) \delta^2

发现得出的数据很奇怪啊。delta E~\delta^2曲线不是平滑曲线。。。。有震荡起伏。。。。好悲催啊。

已调整过很多参数了，还是无法解决这个问题啊。。。不知道哪里出错了。。侯博士说应变\delta的值不能太过接近，否则计算出的能量有误差，但是我的delta取值就是-0.02，-0.015，-0.01，-0.005，0， 0.005（两边对称）。

现在不知道该用什么招了，求助大侠啊！多谢了！
赞一下 回复此楼
7楼2012-10-14 20:15:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 losenq 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 281求调剂（0805） +3 烟汐忆海 2026-03-16 9/450 2026-03-18 07:23 by 烟汐忆海
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 08工科 320总分 求调剂 +4 梨花珞晚风 2026-03-17 4/200 2026-03-17 13:38 by houyaoxu
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 一志愿华中师范071000，325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 328化工专硕求调剂 +4 。，。，。，。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
关闭
请填处理意见
关闭