| 5 | 1/1 | 返回列表 |
| 查看: 1786 | 回復(fù): 11 | ||||
| 本帖產(chǎn)生 2 個 1ST強帖 ,點擊這里進行查看 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
輕風(fēng)龍騎士新蟲 (初入文壇)
|
[交流]
【求助】siesta 做NOSE動力學(xué)出錯,請高手指教 已有4人參與
|
|||
| 在做nose的溫控動力學(xué)的時候,我設(shè)置的初始和目標溫度分別是2700K和2600K,做降溫的動力學(xué)模擬,但是我的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ā)生震動,勢能必然變化,那么各個原子的動能(速度)也隨之變化,因為總能量要保持守恒。這樣就導(dǎo)致了體系溫度的變化。這種漲落是不可避免的,不算是錯誤,但是可能與我們想要的所謂恒定溫度情況有較大的偏差。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料292調(diào)劑 +4 | 橘頌思美人 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料專碩找調(diào)劑 +3 | 哈哈哈吼吼吼哈 2026-03-23 | 3/150 |
|
|
[考研] 327求調(diào)劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +3 | Promise-jyl 2026-03-23 | 3/150 |
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強芒? 2026-03-21 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 268求調(diào)劑 +9 | 簡單點0 2026-03-17 | 9/450 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|