| 8 | 1/1 | 返回列表 |
| 查看: 1033 | 回復(fù): 7 | ||||
chuntao118新蟲 (著名寫手)
|
[求助]
Noneqilibrium dynamics和 equilibrium dynamics 有什么差別? 已有1人參與
|
| Noneqilibrium dynamics和 equilibrium dynamics 有什么差別?在進行MD模擬時又有哪些具體的不同?有誰能夠具體詳細(xì)的解釋一下,或者提供個比較詳細(xì)可靠的資料!謝謝!當(dāng)然,大家也可以發(fā)表自己的見解,和研究心得!謝謝!謝謝!金幣是不是問題。只要回答的可信,有道理。都會奉上金幣! |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +49 |
新蟲 (著名寫手)
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +49 |
|
equilibrium dynamics MD 依賴于系綜那? 是的。只要在平衡態(tài)MD里計算的性質(zhì)都是系綜平均。 是不是說Noneqilibrium dynamics MD 與系綜就沒有關(guān)系? 這個不好說。只能說NEMD不是在一個系綜里進行采樣計算,但其計算的性質(zhì)還是和之前沒有被擾動的系綜是有關(guān)系的。 舉個例子,為了用NEMD計算體系在某個溫度下的剪切粘度,對系統(tǒng)添加一個擾動來產(chǎn)生剪切力,進而通過形變計算粘度。整個模擬的過程是非平衡態(tài)的(系統(tǒng)一直在演化),但模擬得到的粘度,對應(yīng)于擾動之前的系統(tǒng)的粘度。很顯然,擾動之前的處于平衡態(tài)的系統(tǒng),對應(yīng)一個系綜。 系綜這個概念本來就是和平衡態(tài)捆綁在一起的,因為根據(jù)統(tǒng)計物理的遍歷假設(shè),物理量的系綜平均=時間平均=平衡態(tài)的熱力學(xué)值。平衡態(tài)意味著體系各個熱力學(xué)性質(zhì)不隨時間變化,而對于隨時間演化的非平衡態(tài)過程,沒有一個與之對應(yīng)的系綜。 有時對于NEMD模擬,會存在所謂的穩(wěn)態(tài)(steady state),比如體系中存在一個穩(wěn)定的熱流,或者粒子流,但體系的其他性質(zhì)像體積什么的都不變。這種處于穩(wěn)態(tài)的非平衡態(tài),和與之相對應(yīng)的平衡態(tài)的系綜有什么關(guān)聯(lián),是當(dāng)前非平衡態(tài)統(tǒng)計物理研究的一項課題之一, 一個比較典型的例子是處于穩(wěn)態(tài)的self-propelled particle。 |
新蟲 (著名寫手)
送紅花一朵 |
"舉個例子,為了用NEMD計算體系在某個溫度下的剪切粘度,對系統(tǒng)添加一個擾動來產(chǎn)生剪切力,進而通過形變計算粘度。整個模擬的過程是非平衡態(tài)的(系統(tǒng)一直在演化),但模擬得到的粘度,對應(yīng)于擾動之前的系統(tǒng)的粘度。很顯然,擾動之前的處于平衡態(tài)的系統(tǒng),對應(yīng)一個系綜“ 回答的很專業(yè),但是還有點問題想請教一下:就一上面您舉的例子吧?要是計算體系的在某個溫度下的剪切粘度吧?以Poiseuille flow來說吧!在采用力驅(qū)動時,應(yīng)該是NEMD吧(原因:因為體系受外力驅(qū)動)?但是要保證體系的溫度是某個確定的溫度是不是就要采用恒溫器來調(diào)整溫度那?謝謝! |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +49 |
|
我沒有用NEMD計算粘度的實際操作經(jīng)驗,只是對原理略知一二,以下回答可能不準(zhǔn)確。 的確,在lammps用NEMD計算粘度的說明中,當(dāng)使用fix deform命令施加擾動時,需要同時把體系耦合一個恒溫器來防止體系被加熱。 類似的,在某些施加外加電場計算電導(dǎo)率的NEMD模擬中,系統(tǒng)也會耦合一個恒溫器,來防止外加電場對離子的不斷加速產(chǎn)生的升溫。 所以,采用恒溫器在NEMD模擬中控制溫度是說得通的。 需要說明一點,像溫度還有壓力這樣的熱力學(xué)量,它們的定義是只和平衡態(tài)相關(guān)的: 對于溫度,定義就是平衡態(tài)能量對熵的偏導(dǎo)數(shù),在分子模擬中瞬時溫度的計算具體體現(xiàn)為對粒子平均動能的計算,粒子動能的系綜平均就是溫度乘以玻爾茲曼常數(shù)。 對于壓力,定義就是平衡態(tài)能量對體積的偏導(dǎo)數(shù),在分子模擬中瞬時壓力的計算具體體現(xiàn)為對應(yīng)力張量(stress tensor)或者Virial的計算,應(yīng)力張量或者Virial的系綜平均就對應(yīng)于壓力。 對于非平衡態(tài),你沒法把不斷演化的系統(tǒng)和一個系綜聯(lián)系起來(因為做不到遍歷),所以在這種情況下平衡態(tài)下定義的溫度和壓力對于非平衡態(tài)嚴(yán)格來講是失效的。雖然在NEMD模擬中,仍然可以通過計算平均動能和應(yīng)力張量來計算所謂的瞬時溫度和瞬時壓力,但這樣得到的所謂的壓力和溫度的物理意義,目前從理論上來講還不是很清晰,有一些人在做這方面的工作。 |
新蟲 (著名寫手)
木蟲 (正式寫手)
材料大數(shù)據(jù)

| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 332求調(diào)劑 +5 | Zz版 2026-03-13 | 5/250 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +3 | 梨花珞晚風(fēng) 2026-03-17 | 3/150 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +3 | zbzihdhd 2026-03-15 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +3 | Du.11 2026-03-16 | 3/150 |
|
|
[考研] 化學(xué)調(diào)劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[基金申請]
有必要更換申報口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 341求調(diào)劑 +3 | 番茄頭--- 2026-03-10 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | zjy101327 2026-03-11 | 6/300 |
|
|
[考研] 求材料調(diào)劑 +5 | 隔壁陳先生 2026-03-12 | 5/250 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|