24小时热门版块排行榜    

查看: 1234  |  回复: 8

byin

铜虫 (正式写手)


[交流] 用E-V曲线优化结构的时候,能量取哪个?

VASP提供了三种能量,F,F without Entropy,F sigma->0。
个人觉得既然第一原理计算的是0K时的状态,就应该用F=U+TS=U去做结构优化。VASP中的所谓自由能F,是由展宽引入的“温度”,只是为了帮助收敛,没有物理意义。

但是侯博士的教程上写是使用F,不知道大家用什么?为什么?
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★
byin(金币+1): 谢谢参与
liliangfang: 金币+2, 谢谢指教 2012-03-25 09:26:53
fzx2008: 金币+1, 谢谢指教 2012-03-26 14:16:45
F实际是总能U。  F without Entropy 是U-TS。  
你说的对,这里的展宽的引入是为了帮助收敛。   尽管如此,判断sigma的使用值的标准就是:  F和F without Entropy的差小于1meV/atom or  1meV/f.u.
所以其实用哪个都可以。
但如果追求物理上的正确, 还是应该用F。因此候博士的使用是非常正确的。

[ Last edited by gavinliu7390 on 2012-3-26 at 14:28 ]
2楼2012-03-25 08:48:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)


引用回帖:
2楼: Originally posted by gavinliu7390 at 2012-03-25 08:48:56:
F实际是总能U。  F without Entropy 是U+TS。  
你说的对,这里的展宽的引入是为了帮助收敛。   尽管如此,判断sigma的使用值的标准就是:  F和F without Entropy的差小于1meV/atom or  1meV/f.u.
所以其实用哪 ...

你好,谢谢回复。

有两点问题:
1,F without Entropy 从字面理解是从F里面去掉Entropy,那就是U,为什么你解释为U+TS?

2、这里的F不是真正意义上的自由能,为什么还说他是物理的?
3楼2012-03-26 13:14:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★
WDD880227: 金币+2, 非常感谢您的详细解答 2012-03-26 15:43:10
fzx2008: 金币+2, 谢谢指教 2012-03-26 21:43:42
1, 这里vasp的OUTCAR写的是 free energy  TOTEN
energy without entropy
energy(sigma->0)
实际在源代码里写的是
WRITE(TIU6,7261)TOTEN,TOTEN-E%EENTROPY,TOTEN-E%EENTROPY/(2+NORDER)
所以  free energy  TOTEN 就是总能, energy without entropy 是自由能。
2, 因为计算的是0K  所以总能和自由能是相等的, 所以这个Free energy TOTEN是总能,又是没有处理过的自由能。
4楼2012-03-26 14:33:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)


引用回帖:
4楼: Originally posted by gavinliu7390 at 2012-03-26 14:33:57:
1, 这里vasp的OUTCAR写的是 free energy  TOTEN
energy without entropy
energy(sigma->0)
实际在源代码里写的是
WRITE(TIU6,7261)TOTEN,TOTEN-E%EENTROPY,TOTEN-E%EENTROPY/(2+NORDER)
所以  free en ...

十分感谢!

但是如何从他的输出语句判断谁是自由能,谁是内能?可否详细解释一下?万分感激!
5楼2012-03-26 21:36:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
御剑江湖: 金币+3, 谢谢回帖交流 2012-03-28 19:41:29
Free energy TOTEN是内能。  note:尽管它是内能,但是我们做的是0K计算,它还是自由能。
energy without entropy 是自由能。  note: 这里的熵不一定是真是的电子熵,这里熵是用展宽sigma算出来的。
6楼2012-03-27 11:30:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)


引用回帖:
6楼: Originally posted by gavinliu7390 at 2012-03-27 11:30:16:
Free energy TOTEN是内能。  note:尽管它是内能,但是我们做的是0K计算,它还是自由能。
energy without entropy 是自由能。  note: 这里的熵不一定是真是的电子熵,这里熵是用展宽sigma算出来的。

谢谢,我完全明白你的意思。

我现在的疑问是你是如何从输出语句里面得出上述结论的?还请指教!
7楼2012-03-28 18:02:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008: 金币+1, 谢谢指教 2012-03-28 19:01:24
我列的代码不全,如果想完全理解,需要继续读代码,就能弄清楚了。

» 本帖已获得的红花(最新10朵)

8楼2012-03-28 18:55:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)


送鲜花一朵
引用回帖:
8楼: Originally posted by gavinliu7390 at 2012-03-28 18:55:11:
我列的代码不全,如果想完全理解,需要继续读代码,就能弄清楚了。

ok thanks again!
9楼2012-03-28 20:00:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 byin 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 0703化学调剂 ,六级已过,有科研经历 +14 曦熙兮 2026-03-15 14/700 2026-03-21 19:12 by ColorlessPI
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
请填处理意见