版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

bjwang

鐵桿木蟲 (正式寫手)

應助: 4 (幼兒園)
金幣: 5860.3
散金: 94
紅花: 6
帖子: 792
在線: 556.3小時
蟲號: 250077
注冊: 2006-05-08
專業(yè): 凝聚態(tài)物性 II ：電子結構

米）? 已有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計算方法收集
第一性

» 猜你喜歡

287求調(diào)劑已經(jīng)有7人回復
325求調(diào)劑已經(jīng)有5人回復
343求調(diào)劑已經(jīng)有4人回復
求調(diào)劑推薦材料 304 已經(jīng)有4人回復
316求調(diào)劑已經(jīng)有4人回復
351求調(diào)劑已經(jīng)有4人回復
349求調(diào)劑已經(jīng)有5人回復
材料學碩333求調(diào)劑已經(jīng)有7人回復
342求調(diào)劑已經(jīng)有3人回復
調(diào)劑求收留已經(jīng)有7人回復

» 本主題相關商家推薦: (我也要在這里推廣)

» 本主題相關價值貼推薦，對您同樣有幫助:

vasp加電場弛豫原子位置很怪異已經(jīng)有5人回復
vasp計算能帶用origin處理遇到的問題請教已經(jīng)有0人回復
計算輸運的免費軟件已經(jīng)有6人回復
計算含有過渡金屬的氧化物二維體系如何設置自旋情況已經(jīng)有6人回復
求助關于用vasp計算二維石墨烯納米帶所得到的能帶結果與報道對不上已經(jīng)有7人回復

It is a shame for the soul to be first to give way in this life,when thy body does not give way.

1樓 2016-03-27 12:27:45

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

bjwang

鐵桿木蟲 (正式寫手)

應助: 4 (幼兒園)
金幣: 5860.3
散金: 94
紅花: 6
帖子: 792
在線: 556.3小時
蟲號: 250077
注冊: 2006-05-08
專業(yè): 凝聚態(tài)物性 II ：電子結構

利用vasp計算低維體系的應力單位問題（主題有誤，以石墨烯拉5%為例，現(xiàn)更新匯總如下）：
a，對二維體系:
以N/M為單位，結果要乘以真空層厚度
轉換如下： 150.63161KB*1nm=15.06GPa*1nm=15.06*10^9牛/米平方*1nm=15.06N/M
和文獻PRL 102, 235502 (2009)還是可以相比的。
以Gpa為單位，結果要先乘以真空層厚度算出單位長度上的力，再除以石墨烯的實際厚度（這里取石墨的層間距3.4埃），在上面基礎上： 15.06N/M除以3.4埃=44.29Gpa，基本和文獻一致。
b，對一維體系（參考上面的思路，驗證了一些一維體系，至少數(shù)量級是對的）：
算出kb的結果要轉化為GPa：計算值（kb）*真空層面積/(帶寬*層間距)
以上處理方法是否合適（先不考慮精度）？歡迎批評指正

。
另外，也看了一些結合vasp來算應力的，比如PHYSICAL REVIEW B 85, 125428 (2012)，但只有理論和結果，不知具體如何操作的

，哪位有需要或興趣共同探討一下，先謝過了。

贊一下(8人)

回復此樓

It is a shame for the soul to be first to give way in this life,when thy body does not give way.

6樓2016-03-29 13:22:05

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

q965240010

銅蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 216.8
紅花: 2
帖子: 209
在線: 191.3小時
蟲號: 3292577
注冊: 2014-06-25
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

6樓: Originally posted by bjwang at 2016-03-29 13:22:05
利用vasp計算低維體系的應力單位問題（主題有誤，以石墨烯拉5%為例，現(xiàn)更新匯總如下）：
a，對二維體系:
以N/M為單位，結果要乘以真空層厚度
轉換如下： 150.63161KB*1nm=15.06GPa*1nm=15.06*10^9牛/米平方*1n ...

你這樓太好了，自己的帖子解釋的還這么詳細，感謝�。。�

贊一下(2人)

回復此樓

7樓2016-03-29 19:08:10

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

gaozhibin

銅蟲 (正式寫手)

應助: 9 (幼兒園)
金幣: 100.9
散金: 202
紅花: 7
帖子: 391
在線: 149.5小時
蟲號: 2617861
注冊: 2013-08-26
性別: MM
專業(yè): 凝聚態(tài)物性I:結構、力學和

引用回帖:

3樓: Originally posted by bjwang at 2016-03-27 20:59:15
上面文獻的題目為：Graphene Mechanics: Current Status and Perspectives...

你第一幅圖形是來自一篇PRL（他們不是用的VASP，而用的連續(xù)彈性理論和緊束縛模擬），所以這個數(shù)據(jù)是不能和VASP 比較的。
你第二幅圖形是來自一篇 journal of applied physics，他里面用的分子動力學模擬或者MC；

方法不一樣，肯定結果沒有辦法對比。

希望對你有幫助。

贊一下(1人)

回復此樓

» 本帖已獲得的紅花（最新10朵）

bjwang

加油努力，堅持爬行的小烏龜

4樓2016-03-28 14:49:47

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

往事往亦

木蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 1868.1
紅花: 1
帖子: 680
在線: 213.6小時
蟲號: 3159289
注冊: 2014-04-24
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

除以3.4這個石墨層間距，我感覺應該是不對的，另外我通過修改了constr_cell_relax.F源代碼（只優(yōu)化x方向-石墨烯），就是對Y方向施加拉伸應變（拉伸x%后固定，并用的ISIF=3），這樣做后得到的應力-應變數(shù)據(jù)為：
in kB -3.23705 -70.14611 -2.60080 0.00000 0.00000 0.00000 1.02
in kB -0.15934 -27.23946 -4.29416 0.00000 0.00000 0.00000 1.04
in kB 0.13776 -57.52228 -4.20366 0.00000 0.00000 0.00000 1.06
in kB -0.00191 -84.64178 -4.00855 0.00000 0.00000 0.00000 1.08
in kB 0.03199 -108.67303 -3.94575 0.00000 0.00000 0.00000 1.10
in kB -0.13334 -128.35889 -3.79733 0.00000 0.00000 0.00000 1.12
in kB 0.00870 -145.29313 -3.87864 0.00000 0.00000 0.00000 1.14
in kB 0.16160 -160.36230 -3.73108 0.00000 0.00000 0.00000 1.16
in kB 0.23975 -172.93999 -3.73327 0.00000 0.00000 0.00000 1.18
in kB -0.13232 -181.16550 -3.70646 0.00000 0.00000 0.00000 1.20
也發(fā)現(xiàn)了一些問題：
（1）我想要的效果是，石墨烯的晶格常數(shù)b，z值不變，但允許a值變化，并且原子坐標都可以變化，（沒達到,,改過代碼后，晶格常數(shù)b,z不變達到了，a值變化也達到了，但是原子坐標貌似也只優(yōu)化了x方向的分量，，y,z方向的分量沒變化，，）
（2）計算的數(shù)據(jù)中1.02怎么出現(xiàn)了不遵循從小到大的變化了，，（好奇怪）
（3）我的真空層取得是上下各10A的真空厚度，都乘以20A后，，（我感覺像石墨烯這種高對稱，我做的是Y方向拉伸，但仍是應該和你的數(shù)據(jù)符合哎，但怎么會有這么的差別呢，奇怪，

）

贊一下(1人)

回復此樓

青春不留遺憾。

14樓2017-05-05 17:59:25

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

gaozhibin

銅蟲 (正式寫手)

應助: 9 (幼兒園)
金幣: 100.9
散金: 202
紅花: 7
帖子: 391
在線: 149.5小時
蟲號: 2617861
注冊: 2013-08-26
性別: MM
專業(yè): 凝聚態(tài)物性I:結構、力學和

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
bjwang: 金幣+10, ★★★很有幫助, 謝謝！乘以厚度的爭議我知道，但上面乘以真空層厚度還是數(shù)量級不對呀？另外，不知通過3D如何轉化為2D的應力？ 2016-03-27 20:57:42

嚴格來說，VASP 給出的3D 的應力，你需要轉化成2D 的。乘以thickness，但是這個做法是有爭議的，如果你看了一些文章的話。

另外，能把你上面提到的文獻題目貼出來嗎？讓大家都學習下。

贊一下

回復此樓

加油努力，堅持爬行的小烏龜

2樓2016-03-27 20:22:02

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

bjwang

鐵桿木蟲 (正式寫手)

應助: 4 (幼兒園)
金幣: 5860.3
散金: 94
紅花: 6
帖子: 792
在線: 556.3小時
蟲號: 250077
注冊: 2006-05-08
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

2樓: Originally posted by gaozhibin at 2016-03-27 20:22:02
嚴格來說，VASP 給出的3D 的應力，你需要轉化成2D 的。乘以thickness，但是這個做法是有爭議的，如果你看了一些文章的話。

另外，能把你上面提到的文獻題目貼出來嗎？讓大家都學習下。

上面文獻的題目為：Graphene Mechanics: Current Status and Perspectives

贊一下

回復此樓

It is a shame for the soul to be first to give way in this life,when thy body does not give way.

3樓2016-03-27 20:59:15

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

bjwang

鐵桿木蟲 (正式寫手)

應助: 4 (幼兒園)
金幣: 5860.3
散金: 94
紅花: 6
帖子: 792
在線: 556.3小時
蟲號: 250077
注冊: 2006-05-08
專業(yè): 凝聚態(tài)物性 II ：電子結構

送紅花一朵

引用回帖:

4樓: Originally posted by gaozhibin at 2016-03-28 14:49:47
你第一幅圖形是來自一篇PRL（他們不是用的VASP，而用的連續(xù)彈性理論和緊束縛模擬），所以這個數(shù)據(jù)是不能和VASP 比較的。
你第二幅圖形是來自一篇 journal of applied physics，他里面用的分子動力學模擬或者MC； ...

謝謝！

回復此樓

It is a shame for the soul to be first to give way in this life,when thy body does not give way.

5樓2016-03-29 12:36:15

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

bjwang

鐵桿木蟲 (正式寫手)

應助: 4 (幼兒園)
金幣: 5860.3
散金: 94
紅花: 6
帖子: 792
在線: 556.3小時
蟲號: 250077
注冊: 2006-05-08
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

7樓: Originally posted by q965240010 at 2016-03-29 19:08:10
你這樓太好了，自己的帖子解釋的還這么詳細，感謝�。�！...

主要是想和大家交流一下

贊一下

回復此樓

It is a shame for the soul to be first to give way in this life,when thy body does not give way.

8樓2016-03-30 19:51:43

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

bjwang

鐵桿木蟲 (正式寫手)

應助: 4 (幼兒園)
金幣: 5860.3
散金: 94
紅花: 6
帖子: 792
在線: 556.3小時
蟲號: 250077
注冊: 2006-05-08
專業(yè): 凝聚態(tài)物性 II ：電子結構

正如2樓所說，這種處理方法可能不精確，有爭議，哪位算過應力應變關系的同仁，能否出面科普一下如何結合vasp更精確的得到應力應變關系？看了些文獻大部分都有理論推導和結果，不知道具體細節(jié)，是有對應的代碼或程序處理嗎？有了解的指點一二，謝謝！

贊一下

回復此樓

It is a shame for the soul to be first to give way in this life,when thy body does not give way.

9樓2016-03-30 19:58:05

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

q965240010

銅蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 216.8
紅花: 2
帖子: 209
在線: 191.3小時
蟲號: 3292577
注冊: 2014-06-25
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

想問下您這個真空層面積指的是什么~

贊一下

回復此樓

10樓2016-04-12 10:31:33

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 bjwang 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 求調(diào)劑推薦材料 304 +4	荷包蛋hyj 2026-03-26	4/200	2026-03-27 00:39 by wxiongid
[考研] 一志愿鄭州大學，080500學碩，總分317分求調(diào)劑 +4	舉個栗子oi 2026-03-24	5/250	2026-03-26 23:15 by 不吃魚的貓
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 292求調(diào)劑 +9	鵝鵝鵝額額額額?/a> 2026-03-25	10/500	2026-03-26 16:27 by 不吃魚的貓
[考研] 一志愿西北大學總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	7/350	2026-03-26 16:07 by 不吃魚的貓
[考研] 279 分求調(diào)劑 +3	睡個好覺_16 2026-03-24	3/150	2026-03-26 15:43 by zzll406
[考研] 總分293求調(diào)劑 +6	加一一九 2026-03-25	8/400	2026-03-26 13:30 by yujianx
[考研] 315分求調(diào)劑 +5	26考研上岸版26 2026-03-26	5/250	2026-03-26 12:11 by laoshidan
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4	gdlf9999 2026-03-24	4/200	2026-03-25 23:01 by boxking200
[考研] 299求調(diào)劑 +7	某某某某位 2026-03-21	8/400	2026-03-25 20:34 by 熱情沙漠
[考研] 321求調(diào)劑 +3	璞玉~~ 2026-03-25	3/150	2026-03-25 19:07 by Zhanglab-TJU
[考研] 0854人工智能方向招收調(diào)劑 +4	章小魚567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 0703化學求調(diào)劑 +6	奶油草莓. 2026-03-22	7/350	2026-03-25 10:00 by shangxh
[有機交流] 有機合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 材料考研調(diào)劑生 +3	黃粱一夢千年 2026-03-24	3/150	2026-03-24 17:00 by barlinike
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt
[考研] 298求調(diào)劑 +8	上岸6666@ 2026-03-20	8/400	2026-03-23 11:02 by laoshidan
[考研] 276求調(diào)劑 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在風里
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike