| 5 | 1/1 | 返回列表 |
| 查看: 14137 | 回復(fù): 15 | |||||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||||
bjwang鐵桿木蟲 (正式寫手)
|
[求助]
VASP計算二維體系應(yīng)力中的單位問題:KB or N/M(牛/米)? 已有1人參與
|
||||||
|
利用改變晶格常數(shù)的方法計算應(yīng)力應(yīng)變的關(guān)系(stress-strain curve):epsilon = (a-a_{0})/a_{0}; 具體計算中:設(shè)置ISIF = 2,真空層取1nm。 驗證2D的石墨烯在5%的形變下的應(yīng)力(文獻中的計算值近似為150N/M或40GPa) 但我在OUTCAR中找到 “in kB -150.63161”,單位是KB。 1,換算成GPa 按照1KB=0.1GPa,我的結(jié)果150.63161KB=15.6GPa,與40GPa不符; 2,換算成N/M 150.63161KB=15.6GPa=15.6*10^9牛/米平方, 那么按照有關(guān)文獻的說法,算二維體系的應(yīng)力,考慮單位長度上的力,而不是單位面積上的,即單位要表示為牛/米, 牛/米平方要乘上真空層的厚度(這個不確定)? 但結(jié)果也不對:15.6*10^9牛/米平方*1nm=15.6N/M,和文獻中的150N/M剛好錯了一個數(shù)量級。 但有個帖子(http://www.gaoyang168.com/html/201406/7592923.html)說,在vasp計算石墨烯的楊氏模量中,單位就是N/M(為什么???) 那上面的應(yīng)力的單位也應(yīng)是N/M? 如果是這樣,結(jié)果就和文獻符合的很好,但是我上面的換算過程哪里有問題呢? 還是計算有問題? 謝謝各位! 附計算結(jié)果 FORCE on cell =-STRESS in cart. coord. units (eV): Direction XX YY ZZ XY YZ ZX -------------------------------------------------------------------------------------- Alpha Z 11.84580 11.84580 11.84580 Ewald 1223.83420 1231.72516 -1753.82933 0.00000 0.00000 0.00000 Hartree 1199.50298 1201.90716 -1177.06490 0.00000 0.00000 0.00000 E(xc) -57.55665 -57.59624 -59.87916 0.00000 0.00000 0.00000 Local -2603.83155 -2612.99666 2720.52797 0.00000 0.00000 0.00000 n-local -12.23578 -13.54492 -12.15478 0.08818 0.00000 0.00000 augment -1.40512 -1.40962 -1.31894 0.00000 0.00000 0.00000 Kinetic 229.49916 238.74199 271.32154 -4.37353 0.00000 0.00000 Fock 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ------------------------------------------------------------------------------------- Total -10.34696 -1.32733 -0.55181 0.00000 0.00000 0.00000 in kB -150.63161 -19.32335 -8.03326 0.00000 0.00000 0.00000 external pressure = -59.33 kB Pullay stress = 0.00 kB 附文獻stress-strain curve圖: N每M.png GPa.png |
第一性原理計算-tg | 表界面計算 | VASP計算資料收集byZF | vasp計算方法收集 |
第一性 |

鐵桿木蟲 (正式寫手)

|
嚴(yán)格來說,VASP 給出的3D 的應(yīng)力,你需要轉(zhuǎn)化成2D 的。 乘以thickness,但是這個做法是有爭議的,如果你看了一些文章的話。 另外,能把你上面提到的文獻題目貼出來嗎? 讓大家都學(xué)習(xí)下。 |

鐵桿木蟲 (正式寫手)


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩306英一數(shù)二 +8 | z1z2z3879 2026-03-16 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +8 | 曦熙兮 2026-03-15 | 8/400 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] (081700)化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|