 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 14144 | 回復(fù): 15 | |||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||||
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。 驗(yàn)證2D的石墨烯在5%的形變下的應(yīng)力（文獻(xiàn)中的計算值近似為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)的說法，算二維體系的應(yīng)力，考慮單位長度上的力，而不是單位面積上的，即單位要表示為牛/米， 牛/米平方要乘上真空層的厚度（這個不確定）？ 但結(jié)果也不對：15.6*10^9牛/米平方*1nm=15.6N/M，和文獻(xiàn)中的150N/M剛好錯了一個數(shù)量級。 但有個帖子（http://www.gaoyang168.com/html/201406/7592923.html）說，在vasp計算石墨烯的楊氏模量中，單位就是N/M（為什么？？？） 那上面的應(yīng)力的單位也應(yīng)是N/M? 如果是這樣，結(jié)果就和文獻(xiàn)符合的很好，但是我上面的換算過程哪里有問題呢？ 還是計算有問題？ 謝謝各位！ 附計算結(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 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一性原理計算-tg | 表界面計算 | VASP計算資料收集byZF | vasp計算方法收集 |
| 第一性 |
» 猜你喜歡
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿武理材料305分求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有3人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有7人回復(fù)
-
311求調(diào)劑
已經(jīng)有3人回復(fù)
-
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線
已經(jīng)有7人回復(fù)
-
295求調(diào)劑
已經(jīng)有5人回復(fù)
-
354求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- vasp加電場弛豫原子位置很怪異
已經(jīng)有5人回復(fù)
- vasp計算能帶用origin處理遇到的問題請教
已經(jīng)有0人回復(fù)
- 計算輸運(yùn)的免費(fèi)軟件
已經(jīng)有6人回復(fù)
- 計算含有過渡金屬的氧化物二維體系如何設(shè)置自旋情況
已經(jīng)有6人回復(fù)
- 求助關(guān)于用vasp計算二維石墨烯納米帶所得到的能帶結(jié)果與報道對不上
已經(jīng)有7人回復(fù)
bjwang
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 5860.3
- 散金: 94
- 紅花: 6
- 帖子: 792
- 在線: 556.3小時
- 蟲號: 250077
- 注冊: 2006-05-08
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2016-03-27 20:59:15
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
bjwang: 金幣+10, ★★★很有幫助, 謝謝！乘以厚度的爭議我知道，但上面乘以真空層厚度還是數(shù)量級不對呀？另外，不知通過3D如何轉(zhuǎn)化為2D的應(yīng)力？ 2016-03-27 20:57:42
感謝參與，應(yīng)助指數(shù) +1
bjwang: 金幣+10, ★★★很有幫助, 謝謝！乘以厚度的爭議我知道，但上面乘以真空層厚度還是數(shù)量級不對呀？另外，不知通過3D如何轉(zhuǎn)化為2D的應(yīng)力？ 2016-03-27 20:57:42
|
嚴(yán)格來說，VASP 給出的3D 的應(yīng)力，你需要轉(zhuǎn)化成2D 的。 乘以thickness，但是這個做法是有爭議的，如果你看了一些文章的話。 另外，能把你上面提到的文獻(xiàn)題目貼出來嗎？ 讓大家都學(xué)習(xí)下。 |
2樓2016-03-27 20:22:02
bjwang4樓2016-03-28 14:49:47
bjwang
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 5860.3
- 散金: 94
- 紅花: 6
- 帖子: 792
- 在線: 556.3小時
- 蟲號: 250077
- 注冊: 2006-05-08
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2016-03-29 12:36:15
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +7 | SUICHILD 2026-03-12 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 344求調(diào)劑 +5 | knight344 2026-03-16 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +9 | yangfz 2026-03-17 | 9/450 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實(shí)用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|