| 24小時(shí)熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 14138 | 回復(fù): 15 | |||||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||||||
bjwang鐵桿木蟲(chóng) (正式寫(xiě)手)
|
[求助]
VASP計(jì)算二維體系應(yīng)力中的單位問(wèn)題:KB or N/M(牛/米)? 已有1人參與
|
||||||
|
利用改變晶格常數(shù)的方法計(jì)算應(yīng)力應(yīng)變的關(guān)系(stress-strain curve):epsilon = (a-a_{0})/a_{0}; 具體計(jì)算中:設(shè)置ISIF = 2,真空層取1nm。 驗(yàn)證2D的石墨烯在5%的形變下的應(yīng)力(文獻(xiàn)中的計(jì)算值近似為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)文獻(xiàn)的說(shuō)法,算二維體系的應(yīng)力,考慮單位長(zhǎng)度上的力,而不是單位面積上的,即單位要表示為牛/米, 牛/米平方要乘上真空層的厚度(這個(gè)不確定)? 但結(jié)果也不對(duì):15.6*10^9牛/米平方*1nm=15.6N/M,和文獻(xiàn)中的150N/M剛好錯(cuò)了一個(gè)數(shù)量級(jí)。 但有個(gè)帖子(http://www.gaoyang168.com/html/201406/7592923.html)說(shuō),在vasp計(jì)算石墨烯的楊氏模量中,單位就是N/M(為什么???) 那上面的應(yīng)力的單位也應(yīng)是N/M? 如果是這樣,結(jié)果就和文獻(xiàn)符合的很好,但是我上面的換算過(guò)程哪里有問(wèn)題呢? 還是計(jì)算有問(wèn)題? 謝謝各位! 附計(jì)算結(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 附文獻(xiàn)stress-strain curve圖: N每M.png GPa.png |
第一性原理計(jì)算-tg | 表界面計(jì)算 | VASP計(jì)算資料收集byZF | vasp計(jì)算方法收集 |
第一性 |

鐵桿木蟲(chóng) (正式寫(xiě)手)

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

鐵桿木蟲(chóng) (正式寫(xiě)手)


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +9 | 冬十三 2026-03-15 | 10/500 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] 求材料調(diào)劑 +5 | 隔壁陳先生 2026-03-12 | 5/250 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
[考研] 【考研調(diào)劑求收留】 +3 | Ceciilia 2026-03-11 | 3/150 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|