hiyangyue

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.5
帖子: 18
在線: 12小時(shí)
蟲(chóng)號(hào): 996292

[交流] 【求助】在哪里看Hellmann-Feynman stress?

vasp中，在哪里可以看到Hellmann-Feynman stress呢？INCAR里面改變什么參數(shù)，可以改變改力呢？初學(xué)者請(qǐng)多關(guān)照，非常感謝

回復(fù)此樓

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-02-16 17:11:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
帖子: 3782
在線: 873.8小時(shí)
蟲(chóng)號(hào): 1007127

★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
sunyang1988(金幣+1): 謝謝指點(diǎn) 2011-02-19 09:47:47

默認(rèn)的設(shè)置，會(huì)在輸出文件OUTCAR就輸出你提到的物理量的值，通過(guò)在outcar里面搜"Stress"或"stress"之類(lèi)的字眼。

贊一下(3人)

回復(fù)此樓

2樓2011-02-19 09:36:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hiyangyue

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.5
帖子: 18
在線: 12小時(shí)
蟲(chóng)號(hào): 996292

引用回帖:

Originally posted by valenhou001 at 2011-02-19 01:36:39:
默認(rèn)的設(shè)置，會(huì)在輸出文件OUTCAR就輸出你提到的物理量的值，通過(guò)在outcar里面搜"Stress"或"stress"之類(lèi)的字眼。

非常感謝，還是不是很清楚到底是哪個(gè)值。有的文獻(xiàn)會(huì)提到優(yōu)化到Hellmann-Feynman stress 直到0.1GPa, mei\沒(méi)發(fā)現(xiàn)哪里有Hellmann-Feynman stress參數(shù)啊

贊一下

回復(fù)此樓

3樓2011-02-20 11:14:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
帖子: 3782
在線: 873.8小時(shí)
蟲(chóng)號(hào): 1007127

★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
zzy870720z(金幣+2): 謝謝指教 2011-02-20 20:51:15

Hellmann-Feynman stress就是指的是通過(guò)Hellmann-Feynman定理來(lái)計(jì)算出來(lái)施加在原胞上的應(yīng)力，這個(gè)也就是vasp通常計(jì)算出來(lái)的應(yīng)力。
FORCE on cell =-STRESS in cart. coord.  units (eV/reduce length):
Direction X       Y       Z       XY    YZ    ZX
--------------------------------------------------------------------------------------
Alpha Z 2432.44 2432.44 2432.44
Ewald -7379.50  -7379.50  -7379.50    0.00    0.00    0.00
Hartree 851.73 851.73 851.73    0.00    0.00    0.00
  E(xc) -1102.88  -1102.88  -1102.88    0.00    0.00    0.00
Local    1686.25 1686.25 1686.25    0.00    0.00    0.00
n-local -602.33 -600.89 -599.96    -0.53    -0.16    -0.24
augment 927.05 927.05 927.05    0.00    0.00    0.00
Kinetic 3188.25 3187.42 3182.20    1.01    -0.65    0.09
-------------------------------------------------------------------------------------
Total    -0.02    -0.02    -0.02    0.00    0.00    0.00
in kB    -0.08    -0.08    -0.08    0.00    0.00    0.00
external pressure =    -0.08 kB  Pullay stress =       0.00 kB

在OUTCAR中找到類(lèi)似如上輸出。

是否要計(jì)算應(yīng)力，可以通過(guò)ISIF來(lái)控制的，默認(rèn)的是計(jì)算出應(yīng)力。

贊一下(2人)

回復(fù)此樓

4樓2011-02-20 16:40:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hiyangyue

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.5
帖子: 18
在線: 12小時(shí)
蟲(chóng)號(hào): 996292

引用回帖:

Originally posted by valenhou001 at 2011-02-20 08:40:26:
Hellmann-Feynman stress就是指的是通過(guò)Hellmann-Feynman定理來(lái)計(jì)算出來(lái)施加在原胞上的應(yīng)力，這個(gè)也就是vasp通常計(jì)算出來(lái)的應(yīng)力。
FORCE on cell =-STRESS in cart. coord. units (eV/reduce length):
Direc ...

非常非常感謝啦！
如果想畫(huà)應(yīng)力－應(yīng)變曲線，應(yīng)力就用in kB這行對(duì)應(yīng)的值就可以嗎？

贊一下

回復(fù)此樓

5樓2011-02-21 11:00:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
帖子: 3782
在線: 873.8小時(shí)
蟲(chóng)號(hào): 1007127

★
zzy870720z(金幣+1): 謝謝回復(fù) 2011-02-21 15:58:25
hiyangyue(金幣+5): 非常感謝 2011-02-23 14:45:14

可以的。

回復(fù)此樓

6樓2011-02-21 15:40:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hiyangyue

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.5
帖子: 18
在線: 12小時(shí)
蟲(chóng)號(hào): 996292

引用回帖:

Originally posted by valenhou001 at 2011-02-21 15:40:42:
可以的。

我看有的文獻(xiàn)上算應(yīng)力是用一個(gè)公式算出來(lái)的，不是直接用in kB這行的數(shù)值。我不知道用這兩種方法有什么區(qū)別？還是一樣，哪種方法都行呢？

贊一下

回復(fù)此樓

7樓2011-02-25 10:34:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
帖子: 3782
在線: 873.8小時(shí)
蟲(chóng)號(hào): 1007127

★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
gzqdyouxia(金幣+2): 謝謝指教 2011-02-25 20:52:45

貼出你的文獻(xiàn)，指出你所提到的公式，大家才知道你談的是什么，以解釋疑惑。
一種：晶體的應(yīng)變能（晶格應(yīng)變之后的總能減去未應(yīng)變時(shí)的總能）－－－應(yīng)變量，通過(guò)應(yīng)變能與彈性常數(shù)和應(yīng)變量的關(guān)系擬合。
另外一種：應(yīng)變之后體系所收的應(yīng)力－－應(yīng)變量，也可以得到彈性常數(shù)。

贊一下(2人)

回復(fù)此樓

8樓2011-02-25 10:46:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hiyangyue

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.5
帖子: 18
在線: 12小時(shí)
蟲(chóng)號(hào): 996292

引用回帖:

Originally posted by valenhou001 at 2011-02-25 10:46:15:
貼出你的文獻(xiàn)，指出你所提到的公式，大家才知道你談的是什么，以解釋疑惑。
一種：晶體的應(yīng)變能（晶格應(yīng)變之后的總能減去未應(yīng)變時(shí)的總能）－－－應(yīng)變量，通過(guò)應(yīng)變能與彈性常數(shù)和應(yīng)變量的關(guān)系擬合。
另外一種：應(yīng) ...

http://g.zhubajie.com/urllink.php?id=106799792oclivomvqlf44pn
這篇論文中的第48頁(yè)（或論文本身頁(yè)碼的第36頁(yè)）中，應(yīng)力是用一個(gè)公式算出來(lái)的，我不知道是應(yīng)該算呢？還是直接從輸出文件中看？

贊一下

回復(fù)此樓

9樓2011-02-25 18:05:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
帖子: 3782
在線: 873.8小時(shí)
蟲(chóng)號(hào): 1007127

★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
sunyang1988(金幣+2): 謝謝侯老師 2011-02-26 11:07:40

可以的話，能否將找個(gè)文獻(xiàn)的pdf版本發(fā)到valenhou@gmail.com。謝謝了。

贊一下(2人)

回復(fù)此樓

10樓2011-02-25 23:26:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
帖子: 3782
在線: 873.8小時(shí)
蟲(chóng)號(hào): 1007127

★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
zzy870720z(金幣+1): 謝謝分享 2011-02-26 08:28:14

有關(guān)第一原理計(jì)算應(yīng)力的文獻(xiàn)：
O. H. Nielsen and Richard M. Martin,, First-Principles Calculation of Stress Phys. Rev. Lett. 50, 697–700 (1983)

贊一下(3人)

回復(fù)此樓

11樓2011-02-25 23:36:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +21	rare12345 2026-03-18	21/1050	2026-03-20 14:31 by 無(wú)懈可擊111
[基金申請(qǐng)] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+3	babangida 2026-03-19	7/350	2026-03-20 14:13 by hejicker
[考研] 本人考085602 化學(xué)工程專(zhuān)碩 +18	不知道叫什么！ 2026-03-15	20/1000	2026-03-20 13:52 by danney002
[考研] 279分求調(diào)劑一志愿211 +9	chaojifeixia 2026-03-19	10/500	2026-03-20 12:29 by lature00
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3	葵梓衛(wèi)隊(duì) 2026-03-18	5/250	2026-03-19 19:35 by 給你你注意休息
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 287求調(diào)劑 +3	晨昏線與星海 2026-03-19	4/200	2026-03-19 12:32 by peike
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 材料專(zhuān)業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 304求調(diào)劑 +5	素年祭語(yǔ) 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike