| 5 | 1/1 | 返回列表 |
| 查看: 1796 | 回復: 11 | ||||
| 本帖產(chǎn)生 2 個 1ST強帖 ,點擊這里進行查看 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
輕風龍騎士新蟲 (初入文壇)
|
[交流]
【求助】siesta 做NOSE動力學出錯,請高手指教 已有4人參與
|
|||
| 在做nose的溫控動力學的時候,我設置的初始和目標溫度分別是2700K和2600K,做降溫的動力學模擬,但是我的systemlable.MDE文件中怎么顯示的每一步的溫度有問題。涸诘谌降臅r候就已經(jīng)到了2620K左右,之后溫度更低最低只有1400K左右,然后又開始升高,最高高達3000K以上,這是為什么?是不是我那個參數(shù)錯了,請高手指教 |
siesta問題 |

新蟲 (初入文壇)

木蟲 (著名寫手)
葉落鷹飛

銅蟲 (初入文壇)
|
我沒有用過nose模塊,不過我在用verletMD計算時出現(xiàn)過類似的溫度起伏,甚至是劇烈起伏的問題。 我想可能有如下的原因給你參考。 首先宏觀體系的溫度就如二樓所說是個統(tǒng)計平均值,實際上是對分子熱運動的一個平均描述。 到了微觀層次,按照能量均分定理,每個原子,每個自由度能量是1/2kBT,簡化考慮就可以采用1/2kBT=1/2mv**2。意思是在一維無轉(zhuǎn)動的情況下,這一個自由度的能量對應于粒子平動的動能。在程序中就是這么計算得到的溫度。輸出文件中有一個systemlable.XV里面有各個原子不同方向的速度。通過公式就可以把各個速度對應的能量轉(zhuǎn)化為溫度。 你在輸入文件中給定了初始的溫度,程序就給各個原子不同的速度讓其震動起來。但是原子一旦發(fā)生震動,勢能必然變化,那么各個原子的動能(速度)也隨之變化,因為總能量要保持守恒。這樣就導致了體系溫度的變化。這種漲落是不可避免的,不算是錯誤,但是可能與我們想要的所謂恒定溫度情況有較大的偏差。 |
新蟲 (初入文壇)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +7 | 李李不服輸 2026-03-26 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +6 | 小蠟新筆 2026-03-29 | 6/300 |
|
|
[考研] 299求調(diào)劑 +9 | 15188958825 2026-03-25 | 9/450 |
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 一志愿中南大學化學0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
[考研] 一志愿南京航空航天大學材料學碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[有機交流]
高溫高壓反應求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研] 生物學 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
[考研] 材料科學與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|