| 5 | 1/1 | 返回列表 |
| 查看: 14139 | 回復: 15 | |||||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||||
bjwang鐵桿木蟲 (正式寫手)
|
[求助]
VASP計算二維體系應力中的單位問題:KB or N/M(牛/米)? 已有1人參與
|
||||||
|
利用改變晶格常數(shù)的方法計算應力應變的關系(stress-strain curve):epsilon = (a-a_{0})/a_{0}; 具體計算中:設置ISIF = 2,真空層取1nm。 驗證2D的石墨烯在5%的形變下的應力(文獻中的計算值近似為150N/M或40GPa) 但我在OUTCAR中找到 “in kB -150.63161”,單位是KB。 1,換算成GPa 按照1KB=0.1GPa,我的結果150.63161KB=15.6GPa,與40GPa不符; 2,換算成N/M 150.63161KB=15.6GPa=15.6*10^9牛/米平方, 那么按照有關文獻的說法,算二維體系的應力,考慮單位長度上的力,而不是單位面積上的,即單位要表示為牛/米, 牛/米平方要乘上真空層的厚度(這個不確定)? 但結果也不對:15.6*10^9牛/米平方*1nm=15.6N/M,和文獻中的150N/M剛好錯了一個數(shù)量級。 但有個帖子(http://www.gaoyang168.com/html/201406/7592923.html)說,在vasp計算石墨烯的楊氏模量中,單位就是N/M(為什么???) 那上面的應力的單位也應是N/M? 如果是這樣,結果就和文獻符合的很好,但是我上面的換算過程哪里有問題呢? 還是計算有問題? 謝謝各位! 附計算結果 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計算方法收集 |
第一性 |

鐵桿木蟲 (正式寫手)


鐵桿木蟲 (正式寫手)


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +9 | 呼呼呼。。。。 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
[考研] 【0856】化學工程(085602)313 分,本科學科評估A類院校化學工程與工藝,誠求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
[考研] 085600材料與化工 +4 | 安全上岸! 2026-03-16 | 4/200 |
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學上- 2026-03-12 | 5/250 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|