 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 2604 | 回復(fù): 3 | ||
[求助]
lammps中關(guān)于pressure的問題 已有1人參與
|
一直不理解這個問題：在lammps中，首先做結(jié)構(gòu)優(yōu)化，使用fix box/relax 和 minimize 命令，其中fix box/relax 命令設(shè)定壓強(qiáng)為0 bar，最后一般得到的模型的壓強(qiáng)接近0 bar，但是接著加入溫度后，在nvt系綜下做弛豫，壓強(qiáng)就會變很大，幾千bar， 遠(yuǎn)遠(yuǎn)大于標(biāo)準(zhǔn)氣壓，我很是不能理解，這種情況下構(gòu)型應(yīng)該不存在了吧？但查閱了很多計算，比如體積熱容和熱脹系數(shù)的計算，都有這個情況，不知道這里面的物理意義是什么？新手求解答， |
回復(fù)此樓» 猜你喜歡
招收博士1-2人
已經(jīng)有4人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有7人回復(fù)
-
生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑
已經(jīng)有7人回復(fù)
-
0856材料專碩353求調(diào)劑
已經(jīng)有4人回復(fù)
-
山東省面上項目限額評審
已經(jīng)有4人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有5人回復(fù)
-
生物學(xué)071000 329分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿華中科技大學(xué)071000，求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑院校信息
已經(jīng)有4人回復(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顯示的壓力。 仔細(xì)看看看公式，很大的壓力，其實來自于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é)很難達(dá)到這兩點。那么，在樣本不夠的情況下，壓力變化很大也是可以理解的。 以上為本人拙見，希望有所幫助。 |
2樓2015-10-21 01:56:14
|
你好，謝謝你的回答。您的回答很有啟發(fā)，關(guān)于壓強(qiáng)變動劇烈，主要來自于vriial項的說法，我很贊同，我在mail-list上看到steve也有類似的回復(fù)。其實，我的主要問題還不是壓強(qiáng)變動劇烈，我想換個說法可能更容易表述：比如說，我想測1bar（標(biāo)準(zhǔn)大氣壓），300k室溫下的材料熱導(dǎo)，應(yīng)該怎么做呢？我的做法是，先做minimize和fix_relax 讓pressure為1bar，然后依次做nvt ( 300 k，進(jìn)入這個過程時，之前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ā)表 | |
|---|---|---|---|---|
|
[考研] 生物學(xué)調(diào)劑 +4 | Surekei 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +5 | 生物工程調(diào)劑 2026-03-17 | 9/450 |
|
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
|
[考研] 298求調(diào)劑 +4 | 上岸6666@ 2026-03-20 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
|
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗） +3 | 1911236844 2026-03-17 | 3/150 |
|
|
|
[考研]
 271材料工程求調(diào)劑
+8
|
.6lL 2026-03-18 | 8/400 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|