24小时热门版块排行榜

返回列表

2longriver

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 17.2
帖子: 35
在线: 57.6小时
虫号: 4730075
注册: 2016-05-27
专业: 金属功能材料

[求助] 关于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）作为判据呢？

谢谢各位大能