[求助] 关于vasp+phonony计算thermal properties的问题

请教两个关于vasp+phonony计算thermal properties的问题，计算结果示例： #      T [K]      F [kJ/mol]    S [J/K/mol]  C_v [J/K/mol]     E [kJ/mol]         0.000      14.1618956      0.0000000      0.0000000     14.1618956       10.000      14.1618243      0.0352697      0.1253221     14.1621770       20.000      14.1601534      0.3905657      1.3436077     14.1679647 ...... 问题1：单位里[/mol]是指1mol原子，还是每元胞原子数量*1mol，还是每晶胞原子数量*1mol，还是每超胞原子数量*1mol？ 问题2：vasp计算的初始元胞的基态能量小于零（E≈-14ev），扩包2*2*2计算声子谱，T=0K时，F为什么大于零？我知道这个F是Helmholtz free energy，manual里面也给出的数值计算公式。但是我想：F=U-TS，0温0压下，内能U和vasp计算的能量E相差很小，但是phonony计算的结果却是31.32kJ/mol（≈0.15eV/atom），这个F还包含了什么能量呢？ 问题3：我想计算相同材料FCC→BCT的相变温度，一般是用Gibbs  free energy作为两相平衡点判据，在零压下，是不是可以直接用这里的F（Helmholtz free energy）作为判据呢？ 谢谢各位大能

回复此楼

» 猜你喜欢

• FileCode能看出啥？ 已经有23人回复
• 应该是93bebmhtak前后十一个字符比较关键 已经有36人回复
• 关于Filecode分析方法 已经有9人回复
• 2019年青年基金涵评意见，大家看看几个A，几个B? 已经有10人回复
• 不应该看fileCode 已经有6人回复
• 好奇怪的filecode 已经有6人回复
• 长年满屏的广告，版主太不责任了。基金也等的急 已经有4人回复
• filecode 已经有15人回复
• 职称评审，求友友推荐见刊最快的期刊 已经有4人回复
• 奇怪，两个人的filecode固定段从头到尾一模一样 已经有10人回复

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