| 查看: 1263 | 回復(fù): 16 | |||
fannieyang金蟲 (小有名氣)
|
[交流]
分子動學(xué)模擬中的一個問題
|
||
|
做了一個孤立系統(tǒng)分子動力學(xué)模擬測試(即沒有熱浴、沒有耗散力,只有帶電粒子之間的庫倫作用),發(fā)現(xiàn)系統(tǒng)的動能、勢能隨時間都在變大,直到模擬崩掉。感到有這么幾個困惑: (1)程序應(yīng)該沒有問題,因為加上熱浴后的計算結(jié)果就是正確的。 (2)假定程序沒有問題,就從物理基本原理方面理解:系統(tǒng)是孤立的,那么根據(jù)熵增原理,系統(tǒng)狀態(tài)應(yīng)該趨向無序演化,直到平衡態(tài)。因此我能否這樣理解:也許我的程序在崩掉的時候還沒有達到平衡態(tài)。。。 可是系統(tǒng)的勢能、動能隨時間都增大顯然不符合能量守恒,這怎么解釋呢。。。 (3) 無奈,俺只能懷疑求解運動方程中數(shù)值誤差造成的能量不守恒,于是減小了時間步長測試,果然可以計算的時間延長了,但是動能、勢能依然較為緩慢的增加。是不是一定要使用精度更高的數(shù)值解法,才能滿足能量守恒呢? 各位大俠怎么看呢?求指教。 |
| 祝福 |
|
|
|
|
至尊木蟲 (文壇精英)

|
|
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 329求調(diào)劑 +7 | 想上學(xué)吖吖 2026-03-19 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 265求調(diào)劑 +8 | 梁梁校校 2026-03-17 | 8/400 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研]
|
.6lL 2026-03-18 | 7/350 |
|
|
[考研] 304求調(diào)劑 +5 | 曼殊2266 2026-03-18 | 5/250 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|