24小时热门版块排行榜    

CyRhmU.jpeg
查看: 949  |  回复: 7

identation

金虫 (著名写手)

[交流] 【求助】PWSCF如何进行thermal dynamics的计算已有4人参与

请问PWSCF如何计算材料热容随温度的变化以及自由能随温度的变化?用什么模块呢?可否有输入文件的例子和计算步骤提供呢?求高人解答,感激不尽!
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

★ ★
identation(金币+5): 2010-05-20 19:20:18
zzy870720z(金币+2):感谢交流 2010-05-20 22:02:45
引用回帖:
Originally posted by identation at 2010-05-20 17:12:11:
请问PWSCF如何计算材料热容随温度的变化以及自由能随温度的变化?用什么模块呢?可否有输入文件的例子和计算步骤提供呢?求高人解答,感激不尽!

计算声子谱, 可以求得热容和自由能随温度变化。只是这部分没有考虑声子的非谐部分。
真理是一点点接近的!
2楼2010-05-20 19:08:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

引用回帖:
Originally posted by gavinliu7390 at 2010-05-20 19:08:12:


计算声子谱, 可以求得热容和自由能随温度变化。只是这部分没有考虑声子的非谐部分。

那是用ph.x模块吗?那输入的in文件例子你有么?
3楼2010-05-20 19:22:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★
gavinliu7390(金币+3):谢谢提示! 2010-05-20 20:09:57
identation(金币+20):谢谢老师 2010-05-20 23:15:39
参考我这里的blog,并阅读相关的手册部分,去掉计算中去掉la2F=.true.相关的东西就可以了。

http://old.blog.edu.cn/user1/11542/archives/2007/1648431.shtml
4楼2010-05-20 19:53:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)


identation(金币+5): 2010-05-22 09:27:49
identation(金币+100):先奉上金币100,QHA计算还有诸多不懂,望指教,感激不尽!! 2010-05-22 09:30:44
gavinliu7390(金币+1):谢谢交流! 2010-05-22 10:41:24
用最新版本里面的QHA模块,可以处理这些基本性质。
5楼2010-05-20 23:44:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

引用回帖:
Originally posted by goldenfisher at 2010-05-20 23:44:23:
用最新版本里面的QHA模块,可以处理这些基本性质。

QHA里面的examples有si、al和alas,那真正计算的时候是按readme把edit me和环境变量设定好了,然后QHA.x运行run me么?但是赝势文件存放在哪呢?需不需要预先进行自洽计算之类的呢?
6楼2010-05-22 09:29:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★
gavinliu7390(金币+3):感谢提示! 2010-05-22 10:42:10
identation(金币+60):谢谢老师 2010-05-22 11:29:39
QHA.x只是最后一步。

大致的计算过程:

i)给一个晶格常数,进行pw.x模块的计算,
ii)进行ph.x模块的计算,得到力学常数和动力学矩阵元,(ph.x, q2r.x,等)
iii)再就是qha.x,在这个里面设置温度。
然后再给一个晶格常数,重复i), ii), iii)。
7楼2010-05-22 09:46:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★
identation(金币+1):谢谢参与
gavinliu7390(金币+1):谢谢交流! 2010-05-22 10:43:00
正如楼上所说,步骤就是这些。但是必须注意这些性质都是在简谐近似下得到的,温度太高的时候不太准确。
8楼2010-05-22 10:00:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见