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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助：vasp 计算固体相变，能量。
41/1返回列表
查看: 1850  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fangyongxinxi

新虫 (正式写手)

[求助] 求助：vasp 计算固体相变，能量。

自己初次用vasp，在计算一个相变，遇到三个问题：
1. 该物质有两个相，A，B，具有相同原子数，但不同的晶格常数（体积当然也不一样），现在需要比较两个体系的能量。自己的做法是：对A进行结构优化，在计算完后，采用 grep entropy OUTCAR , 提取 energy  without entropy 后面的数值，当做A体系的能量（Ea）；同样的方法对 B 结构进行操作，得到能量(Eb)。
问题：Ea - Eb 就是两个相的能量差吗 ？（这里不考虑过渡态问题，只考虑两个相绝对能量差），是不是还需要有一个有关体积的能量项，比如，大气压*（AB 体积差）？
2. 计算单个体系的能量，直接采用优化结构后OUTCAR中的 “采用 grep entropy OUTCAR , 提取 energy  without entropy 后面的数值”，这样对吗？还是需要将优化好的结构CONTCAR拿出来，在进行能量计算？（自己属于初级阶段，具体操作过程还不太清楚）。
3. 补充问一个K点的问题。以前问过如果一个体系是8个原子，K点设置为8*8*8；如果体系扩大为2*2*2倍，变成64个原子体系，则K点相应变为4*4*4.
现在的问题是：8个原子体系为7个原子加一个杂质原子，K点设置为8*8*8；当体系变成64个原子时候（63个原子+1个杂质原子），这时候K点还是4*4*4 吗？

期待高手解答，谢谢啦~\(≧▽≦)/~啦啦啦

[ Last edited by fangyongxinxi on 2011-11-15 at 10:28 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-11-15 10:19:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

fangyongxinxi(金币+14): 2011-11-15 14:20:11
引用回帖:
1楼: Originally posted by fangyongxinxi at 2011-11-15 10:19:38:
自己初次用vasp，在计算一个相变，遇到三个问题：
1. 该物质有两个相，A，B，具有相同原子数，但不同的晶格常数（体积当然也不一样），现在需要比较两个体系的能量。自己的做法是：对A进行结构优化，在计算完后， ...

1， 在两者晶胞大小一致的时候可以，但是原子数如果不相等就要乘以一个原子数的比
2，在T=0时，这两个能量相差很小，
如：
  free  energy   TOTEN  =       -15.732155 eV
  energy  without entropy=      -15.732155  energy(sigma->0) =      -15.732155

取前面的值吧

3，掺杂后k点还是一样的。

另外， 你优化后，加大k点提高精度，使用ismear=-5重新计算一次能量，不要直接取优化的用
赞一下 回复此楼
嘟啦啦嘟啦啦嘟
2楼2011-11-15 13:43:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangyongxinxi

新虫 (正式写手)
引用回帖:
2楼: Originally posted by 贺仪 at 2011-11-15 13:43:24:
1， 在两者晶胞大小一致的时候可以，但是原子数如果不相等就要乘以一个原子数的比
2，在T=0时，这两个能量相差很小，
如：
  free  energy   TOTEN  =       -15.732155 eV
  energy  without entropy=      ...

第一个问题的答案我还没完全明白。
版主的意思是，我算的Ea-Eb 就是我需要的两个体系的绝对能量差，不需要考虑体积变化引起的能量部分？
谢谢。
赞一下 回复此楼
3楼2011-11-15 14:23:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)
free energy 中已包含PV
赞一下 回复此楼
嘟啦啦嘟啦啦嘟
4楼2011-11-15 14:34:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fangyongxinxi 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +7 大米饭！ 2026-03-15 7/350 2026-03-16 10:25 by 了了了了。。
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 复试调剂 +3 呼呼？~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
信息提示
关闭
请填处理意见
关闭