| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 14207 | 回復(fù): 15 | |||||||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||||||
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ě)手)
|
利用vasp計(jì)算低維體系的應(yīng)力單位問(wèn)題(主題有誤,以石墨烯拉5%為例,現(xiàn)更新匯總?cè)缦拢?br />
a,對(duì)二維體系: 以N/M為單位,結(jié)果要乘以真空層厚度 轉(zhuǎn)換如下: 150.63161KB*1nm=15.06GPa*1nm=15.06*10^9牛/米平方*1nm=15.06N/M 和文獻(xiàn)PRL 102, 235502 (2009)還是可以相比的。 以Gpa為單位,結(jié)果要先乘以真空層厚度算出單位長(zhǎng)度上的力,再除以石墨烯的實(shí)際厚度(這里取石墨的層間距3.4埃),在上面基礎(chǔ)上: 15.06N/M除以3.4埃=44.29Gpa,基本和文獻(xiàn)一致。 b,對(duì)一維體系(參考上面的思路,驗(yàn)證了一些一維體系,至少數(shù)量級(jí)是對(duì)的): 算出kb的結(jié)果要轉(zhuǎn)化為GPa:計(jì)算值(kb)*真空層面積/(帶寬*層間距) 以上處理方法是否合適(先不考慮精度)?歡迎批評(píng)指正 。另外,也看了一些結(jié)合vasp來(lái)算應(yīng)力的,比如PHYSICAL REVIEW B 85, 125428 (2012),但只有理論和結(jié)果,不知具體如何操作的 ,哪位有需要或興趣共同探討一下,先謝過(guò)了。 |

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

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


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +4 | 篤行致遠(yuǎn). 2026-03-28 | 4/200 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專(zhuān)業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 調(diào)劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 08開(kāi)頭275求調(diào)劑 +4 | 拉誰(shuí)不重要 2026-03-26 | 4/200 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生(材料、化學(xué)、電化學(xué),環(huán)境) +4 | 我愛(ài)學(xué)電池 2026-03-23 | 4/200 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問(wèn)津 2026-03-23 | 5/250 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專(zhuān)碩) +4 | allen-yin 2026-03-23 | 6/300 |
|