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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 用侯博士计算弹性常数的方法，计算NITI B2相的C11和C12出现了问题。求助高手！
71/1返回列表
查看: 2238  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 环境领域全国重点实验室招收博士1-2名 +3 QGZDSYS 2026-03-13 3/150 2026-03-18 00:13 by PY-Hzb
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +7 yangfz 2026-03-17 7/350 2026-03-17 23:57 by 星空星月
[考研] 085601专硕，总分342求调剂，地区不限 +4 share_joy 2026-03-16 4/200 2026-03-17 21:32 by hmn_wj
[考研] 材料与化工求调剂 +6 为学666 2026-03-16 6/300 2026-03-17 20:15 by peike
[考研] 08工科 320总分 求调剂 +4 梨花珞晚风 2026-03-17 4/200 2026-03-17 13:38 by houyaoxu
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] （081700）化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
信息提示
关闭
请填处理意见
关闭