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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Abinit » 重金请教采用准谐近似求热膨胀系数的方法
51/1返回列表
查看: 5157  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 2 个 1ST强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

mazuju028

木虫 (正式写手)

[求助] 重金请教采用准谐近似求热膨胀系数的方法

在abinit教程中只看到用简谐近似计算热动力学的例子(trf2_7.in)，可是在文献中看到很多都是采用准谐近似，两者应该是有所区别的吧。通过文献阅读，我罗列了一下基本计算步骤，以立方晶系为例，
1，基于平衡晶格参数a0,罗列一组晶格常数：a1,a2,a3,a4……
2，对上面给出的一组结构进行结构优化，可以得到V1，V2，V3……及其对应的能量E1，E2，E3……
3，分别计算V1，V2，V3……所对应结构的声子谱；
4，分别计算V1，v2，V3……所对应的结构的热动力学量（不同温度下的声子自由能，声子内能，声子比热）；
5，问题就是，怎么基于上面给出的数据得到热膨胀系数：1/V*(dV/dT)?
还有个疑问就是上面步骤是否有问题？
谢谢指教。

[ Last edited by mazuju028 on 2011-5-12 at 09:12 ]
回复此楼

» 收录本帖的淘帖专辑推荐

仿真建模与计算 MS资源 培训、总结 心路历程

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
锐气藏于胸，和气浮于脸，才气见于事，义气施与人
1楼 2011-05-12 08:51:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
引用回帖:
Originally posted by valenhou001 at 2011-05-14 22:51:56:
大致可采用这样的步骤：
1）在实验或理论计算出来的晶格常数或平衡体积(V0)的左右各取6～7体积点，使得所要计算的体积点的个数在10个左右，而且这些体积的范围大致在[(1- x) V0,  (1+x) V0]，x为小于0.1~0.3的正 ...

呵呵， 说的挺好！
回复此楼
真理是一点点接近的！
10楼2011-05-15 09:36:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

gavinliu7390

木虫 (著名写手)

叶落鹰飞

【答案】应助回帖

★ ★ ★
zzy870720z(金币+3, 1ST强帖+1): 谢谢详细指教 2011-05-12 15:15:38
工具，vasp+phonopy/fropho/phonon 或者也可以用线性响应的办法(abinit or quantum espresso)，计算材料是 Zr。
公式是膨胀系数=1/v*(Δv/ΔT).
正常在 ab initio  的计算中，是无法计算热膨胀的。需要考虑声子
振动才能进行计算。计算这个还是稍微有一点点麻烦。例如，我想计
算 0 压的热膨胀系数。需要计算对应 0 压下各个温度的体积(这里的 0
压主要指的是声子是 0 压的)。进而带入公式求得热膨胀系数。
有 F（v，T）=Min（Uab+Uph+PV）是 Uab 是 ab inito 的能量，Uph
是声子振动能量。Uab和 Uph都是 V 和 T 的函数。这里由于我想计算
的是 0 压强，所以 PV=0.
所以我需要改变体积，计算不同体积下的声子谱，进而求得不同
体积下的声子能量。再加上该体积下 vasp 计算的能量，就是总能，
总能最小即使是平衡体积。这个平衡体积即是在这个温度下的真实体
积。这主要是由于考虑了离子振动的能量。之所以麻烦，是因为每一
个温度，都需要做物态方程拟合，求平衡体积。
赞一下(2人) 回复此楼
真理是一点点接近的！
2楼2011-05-12 15:04:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)
引用回帖:
Originally posted by gavinliu7390 at 2011-05-12 15:04:48:
工具，vasp+phonopy/fropho/phonon 或者也可以用线性响应的办法(abinit or quantum espresso)，计算材料是 Zr。
公式是膨胀系数=1/v*(Δv/ΔT).
正常在 ab initio  的计算中，是无法计算热膨胀的。需要考虑声子 ...

你好，非常感谢你的细致回答，我用的是abinit线性响应的方法。你的说法中我有一点疑问，你刚开始说的就是0压，可是在后面需要求不同体积下的声子能量时，这里不同的体积已经对应不同的压强了啊。和前面说的0压是否冲突了？
还有一个疑问就是针对不同体积需要结构优化吗？我用abinit尝试过发现如果是全优化，所有的体积都会优化到同一个值，就是0压下的平衡体积，不知道你是怎么处理的？

呵呵，还有你再回帖子时填上应助回帖，不然给不了你金币，呵呵。

[ Last edited by mazuju028 on 2011-5-13 at 22:42 ]
赞一下 回复此楼
锐气藏于胸，和气浮于脸，才气见于事，义气施与人
3楼2011-05-13 08:22:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)
我对二楼的问题，有谁知道吗？顶一顶了，哈哈。
赞一下 回复此楼
锐气藏于胸，和气浮于脸，才气见于事，义气施与人
4楼2011-05-13 22:31:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭