24小时热门版块排行榜    

查看: 1581  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +7 wangmingqi 2026-08-28 8/400 2026-08-28 18:42 by zhuzg0628
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +11 Kittylucky 2026-08-27 12/600 2026-08-28 16:26 by lyxiaomc
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见