 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 2605 | 回復(fù): 3 | |||
[求助]
lammps中關(guān)于pressure的問題 已有1人參與
|
一直不理解這個問題：在lammps中，首先做結(jié)構(gòu)優(yōu)化，使用fix box/relax 和 minimize 命令，其中fix box/relax 命令設(shè)定壓強為0 bar，最后一般得到的模型的壓強接近0 bar，但是接著加入溫度后，在nvt系綜下做弛豫，壓強就會變很大，幾千bar， 遠遠大于標(biāo)準(zhǔn)氣壓，我很是不能理解，這種情況下構(gòu)型應(yīng)該不存在了吧？但查閱了很多計算，比如體積熱容和熱脹系數(shù)的計算，都有這個情況，不知道這里面的物理意義是什么？新手求解答， |
回復(fù)此樓» 猜你喜歡
求調(diào)劑院校信息
已經(jīng)有6人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有17人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑
已經(jīng)有3人回復(fù)
-
260求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過
已經(jīng)有7人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有7人回復(fù)
-
275求調(diào)劑
已經(jīng)有7人回復(fù)
-
070300化學(xué)求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 關(guān)于Green-Kubo方法模擬水的粘度
已經(jīng)有5人回復(fù)
- 求助lammps輸出與輸入的問題
已經(jīng)有0人回復(fù)
- 關(guān)于lammps中expected floating point parameter in input
已經(jīng)有13人回復(fù)
- 關(guān)于LAMMPS 如果定義charge
已經(jīng)有0人回復(fù)
- 【求助】關(guān)于LAMMPS建模
已經(jīng)有9人回復(fù)
- 關(guān)于LAmmps的windows下源程序編譯調(diào)試問題
已經(jīng)有7人回復(fù)
- 關(guān)于 run 0 command 重要性
已經(jīng)有0人回復(fù)
- 求助各位前輩關(guān)于在LAMMPS中angle_style的設(shè)置問題。
已經(jīng)有2人回復(fù)
- 關(guān)于towhee 回答zhang_jaj的問題
已經(jīng)有44人回復(fù)
- 一個關(guān)于非周期性邊界條件的問題
已經(jīng)有5人回復(fù)
- 【分享】mdbbs發(fā)布了三年優(yōu)質(zhì)資源匯總電子書
已經(jīng)有48人回復(fù)
- 【分享】mdbbs發(fā)布了三年優(yōu)質(zhì)資源匯總電子書
已經(jīng)有40人回復(fù)
alundilong
木蟲 (正式寫手)
- 應(yīng)助: 129 (高中生)
- 金幣: 2019.1
- 紅花: 18
- 帖子: 309
- 在線: 56小時
- 蟲號: 595468
- 注冊: 2008-09-08
- 性別: GG
- 專業(yè): 傳熱傳質(zhì)學(xué)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
從宏觀角度看，nvt，密度溫度確定。假設(shè)體系熱平衡，那么兩個狀態(tài)量定下來，所有其他狀態(tài)都定了，包括壓力。然而現(xiàn)在壓力變化很大，這個應(yīng)該是讓樓主困惑的地方吧？ 這個問題不好回答吧。 可能找本統(tǒng)計力學(xué)的書看看，會有收獲。 在分子動力學(xué)中，我覺得這里的壓力用應(yīng)力代替更加貼切。 這里是計算應(yīng)力的公式 http://lammps.sandia.gov/doc/compute_stress_atom.html 尤其這一段話： compute peratom all stress/atom NULL compute p all reduce sum c_peratom[1] c_peratom[2] c_peratom[3] variable press equal -(c_p[1]+c_p[2]+c_p[3])/(3*vol) thermo_style custom step temp etotal press v_press 說明了壓力就是應(yīng)力張量對角元素的均值（假設(shè)不可壓，這是成立的） 這個compute計算出的壓力就是lammps顯示的壓力。 仔細看看看公式，很大的壓力，其實來自于virial項（動能變化不大），導(dǎo)致了內(nèi)部應(yīng)力變化劇烈，表現(xiàn)出來的就是很大的壓力變化。 另外一個分析角度就是尺度了： 分子動力學(xué)計算的空間，時間尺度極小。 時間上：分子松弛時間的倍數(shù)（可以很多很多倍，然而還是不夠大） 空間上：分子自由程的倍數(shù)（可以很多很多倍，然而還是不夠） 平時所講的壓力，其實是個統(tǒng)計量（統(tǒng)計力學(xué)）。統(tǒng)計需要足夠的統(tǒng)計樣本，換言之，時間要足夠長，空間也要足夠大�，F(xiàn)在分子動力學(xué)很難達到這兩點。那么，在樣本不夠的情況下，壓力變化很大也是可以理解的。 以上為本人拙見，希望有所幫助。 |
2樓2015-10-21 01:56:14
|
你好，謝謝你的回答。您的回答很有啟發(fā)，關(guān)于壓強變動劇烈，主要來自于vriial項的說法，我很贊同，我在mail-list上看到steve也有類似的回復(fù)。其實，我的主要問題還不是壓強變動劇烈，我想換個說法可能更容易表述：比如說，我想測1bar（標(biāo)準(zhǔn)大氣壓），300k室溫下的材料熱導(dǎo)，應(yīng)該怎么做呢？我的做法是，先做minimize和fix_relax 讓pressure為1bar，然后依次做nvt ( 300 k，進入這個過程時，之前0 bar 的系統(tǒng) pressure 陡增到很大），npt (1 bar，體積變化，pressure變動劇烈)，nvt (300 k)，然后再測熱導(dǎo)，不知道這樣做對不對呢？我看到lammps的例子中并沒有考慮得這么復(fù)雜，直接nvt平衡后就開始測熱導(dǎo)了。希望能幫助解答，謝謝。 |
3樓2015-10-23 21:05:53
liuxun200376
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 425.2
- 散金: 100
- 紅花: 1
- 帖子: 156
- 在線: 157.4小時
- 蟲號: 1296409
- 注冊: 2011-05-15
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
4樓2015-10-26 10:23:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 未曾彷徨 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+3
|
酥酥魚.. 2026-03-21 | 3/150 |
|
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[論文投稿]
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|