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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 用侯博士计算弹性常数的方法，计算NITI B2相的C11和C12出现了问题。求助高手！
51/1返回列表
查看: 2273  |  回复: 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的回帖

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的回帖
查看全部 7 个回答

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  ...

您好，我在做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的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂 +24 不逢春 2026-04-07 25/1250 2026-04-10 09:39 by 690616278
[考研] 22408 366分，本科211，一志愿西工大 +3 Rubt 2026-04-09 3/150 2026-04-10 09:38 by liuhuiying09
[考研] 材料与化工371求调剂 +17 陪琳看海 2026-04-04 18/900 2026-04-10 07:51 by 314126402
[考研] 求调剂 +5 archer.. 2026-04-09 7/350 2026-04-09 22:18 by lbsjt
[考研] 0703化学 +31 妮妮ninicgb 2026-04-04 35/1750 2026-04-09 21:06 by zhouxiaoyu
[硕博家园] 有没有学校材料专业收跨调(一志愿085410) +5 momo(上岸版) 2026-04-06 8/400 2026-04-09 15:07 by only周
[考研] 22408 270分 +9 sanjin020722 2026-04-08 10/500 2026-04-09 14:11 by 诗与自由
[考研] 一志愿985初试354分生物调剂 +3 031001 2026-04-06 3/150 2026-04-09 00:30 by Evan_Liu
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 一志愿哈工大，初试329，求环境科学与工程调剂！ +11 余未辛 2026-04-06 11/550 2026-04-08 15:21 by screening
[考研] 274求调剂求调剂 +10 Jachenbingoo 2026-04-06 13/650 2026-04-08 14:25 by zhq0425
[考研] 284求调剂 +17 梵@@ 2026-04-06 17/850 2026-04-08 11:35 by 1shin_ichi
[考研] 295求调剂 +18 xndjjj 2026-04-04 19/950 2026-04-07 11:02 by wangjy2002
[考研] 机械专硕274求调剂，不挑专业学校 +6 泛泛2333 2026-04-05 8/400 2026-04-06 18:06 by 泛泛2333
[考研] 341求调剂 +3 学无止境，冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 一志愿南农090401，268，求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 调剂 +5 asdasdassda 2026-04-03 6/300 2026-04-03 20:27 by 岸上的一条鱼
信息提示
关闭
请填处理意见
关闭