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

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

新蟲 (初入文壇)

新蟲 (初入文壇)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京理工大學本科211材料工程294求調(diào)劑 +7 | mikasa的圍巾 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 291求調(diào)劑 +4 | Y-cap 2026-03-29 | 5/250 |
|
|
[考研] 329求調(diào)劑 +7 | 星野? 2026-03-26 | 7/350 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 295求調(diào)劑 +4 | wei-5 2026-03-26 | 4/200 |
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 一志愿中南大學化學0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 0703化學338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 打過很多競賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|