24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1249  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289求调剂 +11 新时代材料 2026-03-27 11/550 2026-03-28 18:09 by YK。。
[考研] 一志愿北京工业大学,324分求调剂 +4 零八# 2026-03-28 4/200 2026-03-28 15:01 by 17865157980
[考研] 0856求调剂 +11 zhn03 2026-03-25 12/600 2026-03-28 13:32 by 唐沐儿
[考研] 求调剂 +6 芦lty 2026-03-25 7/350 2026-03-28 13:13 by 唐沐儿
[考研] 291求调剂 +15 hhhhxn.. 2026-03-23 21/1050 2026-03-28 11:26 by self2008
[考研] 308求调剂 +7 墨墨漠 2026-03-27 7/350 2026-03-28 07:43 by 热情沙漠
[考研] 295求调剂 +5 1428151015 2026-03-27 6/300 2026-03-28 04:04 by fmesaito
[考研] 265求调剂 +8 小木虫085600 2026-03-27 8/400 2026-03-27 22:16 by 无际的草原
[考研] 285求调剂 +4 AZMK 2026-03-27 7/350 2026-03-27 20:59 by AZMK
[考研] 一志愿南师大0703化学 275求调剂 +4 Ripcord上岸 2026-03-27 4/200 2026-03-27 17:00 by zhyzzh
[考研] 材料求调剂 +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by 不吃魚的貓
[考研] 359求调剂 +4 王了个楠 2026-03-25 4/200 2026-03-27 08:43 by 不吃魚的貓
[考研] 材料调剂 +8 匹克i 2026-03-23 8/400 2026-03-27 08:11 by hypershenger
[考研] 总分293求调剂 +6 加一一九 2026-03-25 8/400 2026-03-26 13:30 by yujianx
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 0854电子信息求调剂 +7 α____ 2026-03-22 9/450 2026-03-25 13:37 by α____
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 361求调剂 +3 Glack 2026-03-22 3/150 2026-03-23 22:03 by fuyu_
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
信息提示
请填处理意见