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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 求助审稿人提出的有关晶格优化的问题
182/2返回列表上一页12
查看: 1242  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

iamkz

至尊木虫 (小有名气)

寻识小生
引用回帖:
9楼: Originally posted by babaleo at 2013-03-03 10:09:41
不妨就算几个体积优化, 然后B-M状态方程拟和一下得出零压结果, 回复时把两者结果比较一下就应该搞定了吧....

刚才我查了一下以前的帖子，我看到你在别人帖子里回复怎样进行 体积优化 了，谢谢你。
我计算的是六方结构，按照你在别的帖子里的方法，“建模时用想要的体积V, 估个a值, 用六方晶系体积公式算出c,然后优化.”  在CASTEP优化的时候“选optimize cell，到xxxxxx.cell 文件里加上:"fix_vol : true"然后,这样可以同时优化原子位置和c/a但是保持体积不变”
是不是这不算出来得到E-V曲线，然后就用你说的B-M状态方程进行拟和，得到E最小时的V值，从而得到晶格常数a和c。

这里有个问题，这么做的话是不是六方结构的角度没有改变呢？谢谢
赞一下 回复此楼
落花焉知流水之意
11楼2013-03-04 17:12:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)
引用回帖:
11楼: Originally posted by iamkz at 2013-03-04 17:12:14
刚才我查了一下以前的帖子，我看到你在别人帖子里回复怎样进行 体积优化 了，谢谢你。
我计算的是六方结构，按照你在别的帖子里的方法，“建模时用想要的体积V, 估个a值, 用六方晶系体积公式算出c,然后优化.”  在 ...

模型还是六方(不要降P1), 所以只有a,c 和原子位置是变量, 其他都被对称性固定了.
赞一下 回复此楼
12楼2013-03-04 21:44:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamkz

至尊木虫 (小有名气)

寻识小生
引用回帖:
12楼: Originally posted by babaleo at 2013-03-04 21:44:00
模型还是六方(不要降P1), 所以只有a,c 和原子位置是变量, 其他都被对称性固定了....

按你这么计算 “用你说的B-M状态方程进行拟和，得到E最小时的V值，从而得到晶格常数a和c。” 可是计算过程中c/a同样也是优化的，这样最终得到V能推断出a和c吗？谢谢
还有我昨晚查找了一些文献，B-M状态方程的表达式好像有很多种，对于六方结构用哪种比较好呢？
麻烦你了
赞一下 回复此楼
落花焉知流水之意
13楼2013-03-05 09:43:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)

【答案】应助回帖

★ ★
iamkz: 金币+2, ★★★很有帮助, 谢谢 2013-03-05 12:32:32
引用回帖:
13楼: Originally posted by iamkz at 2013-03-05 09:43:49
按你这么计算 “用你说的B-M状态方程进行拟和，得到E最小时的V值，从而得到晶格常数a和c。” 可是计算过程中c/a同样也是优化的，这样最终得到V能推断出a和c吗？谢谢
还有我昨晚查找了一些文献，B-M状态方程的表达 ...

CASTEP里你如果只是固定体积, 优化过程中c, a和原子位置一起都优化了, 所以只要找到E最小时的V值,在那个V优化时, a 和c 都应该有被优化. 你不妨用两组不同a,c, 但是V相同做起始模型, 自己验证一下. 优化后的结果应该一样.
至于B-M EOS, 不管什么结构都适用, 都是针对块体均匀材料的.
赞一下 回复此楼
14楼2013-03-05 11:09:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamkz

至尊木虫 (小有名气)

寻识小生
引用回帖:
14楼: Originally posted by babaleo at 2013-03-05 11:09:47
CASTEP里你如果只是固定体积, 优化过程中c, a和原子位置一起都优化了, 所以只要找到E最小时的V值,在那个V优化时, a 和c 都应该有被优化. 你不妨用两组不同a,c, 但是V相同做起始模型, 自己验证一下. 优化后的结果应 ...

哦，这下我明白了，谢谢你，我知道该怎么计算了。
可是对于B-M EOS还有有点疑问，我在文献上看到两种不同的表达式，具体用哪种呢？谢谢

  

不知道图片上传得对不对
赞一下 回复此楼
落花焉知流水之意
15楼2013-03-05 12:37:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamkz

至尊木虫 (小有名气)

寻识小生
引用回帖:
14楼: Originally posted by babaleo at 2013-03-05 11:09:47
CASTEP里你如果只是固定体积, 优化过程中c, a和原子位置一起都优化了, 所以只要找到E最小时的V值,在那个V优化时, a 和c 都应该有被优化. 你不妨用两组不同a,c, 但是V相同做起始模型, 自己验证一下. 优化后的结果应 ...

查看附件吧，不好意思，不太会传图片

1_副本.jpg



图片1.jpg
赞一下 回复此楼
落花焉知流水之意
16楼2013-03-05 12:38:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)
Equation of state: Birch-Murnaghan
E = E0 + (9/16)*B*V0*[(eta**2-1)**3*BP + (eta**2-1)**2*(6-4*eta**2)]
     --> eta = (V0/V)**(1/3)
1 eV/Angstrom^3 = 160.2176487 GPa
赞一下 回复此楼
17楼2013-03-05 21:55:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamkz

至尊木虫 (小有名气)

寻识小生
引用回帖:
17楼: Originally posted by babaleo at 2013-03-05 21:55:27
Equation of state: Birch-Murnaghan
E = E0 + (9/16)*B*V0*
     --> eta = (V0/V)**(1/3)
1 eV/Angstrom^3 = 160.2176487 GPa

谢谢你，我试着拟合看看
那我在文章里引用的话，应该是引用我贴出来的哪个公式比较好呢？
赞一下 回复此楼
落花焉知流水之意
18楼2013-03-06 08:01:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 iamkz 的主题更新
182/2返回列表上一页12
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭