[交流] 【求助】分子動力學中體系溫度不斷增加? 已有11人參與

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有8人回復
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名已經(jīng)有3人回復
材料專碩306英一數(shù)二已經(jīng)有10人回復
301求調(diào)劑已經(jīng)有6人回復
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復
302求調(diào)劑已經(jīng)有6人回復
26博士申請已經(jīng)有3人回復
268求調(diào)劑已經(jīng)有3人回復
311求調(diào)劑已經(jīng)有10人回復
被我言中：新模板不強調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

誰能告訴我一下用分子動力學怎么計算能量和溫度，謝謝了啊已經(jīng)有7人回復
國內(nèi)有開發(fā)分子動力學勢函數(shù)的課題組嗎？已經(jīng)有12人回復
【討論】誰能幫我講講分子動力學中的parrinello-Rahman方法？已經(jīng)有4人回復
【討論】大家好，加水后跑完分子動力學，在log中發(fā)現(xiàn)了問題已經(jīng)有15人回復
【求助】用分子動力學觀點分析晶體熔化為何除保持一定溫度外還需提供一定熱量？已經(jīng)有6人回復
【求助】分子動力學中退火求助已經(jīng)有14人回復

1樓 2010-09-12 13:37:05

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yahoohoo

鐵桿木蟲 (著名寫手)

模擬EPI: 10
應(yīng)助: 0 (幼兒園)
貴賓: 1.55
金幣: 7632.4
散金: 251
紅花: 4
帖子: 1176
在線: 167小時
蟲號: 74894
注冊: 2005-06-15
性別: GG
專業(yè): 理論和計算化學

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zh1987hs(金幣+1):謝謝 2010-09-12 15:11:45

引用回帖:

Originally posted by weizhiyong at 2010-09-12 13:37:05:
在MD中,體系的動能和勢能之和不守恒,而且體系的溫度不斷增加,是什么原因啊?

先交待下你研究的體系，以及系綜，然后是怎么做MD的。這樣才能分析原因。

2樓2010-09-12 14:44:24

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

weizhiyong

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1365.8
散金: 11
紅花: 1
帖子: 356
在線: 178.7小時
蟲號: 585996
注冊: 2008-08-08
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

Originally posted by yahoohoo at 2010-09-12 14:44:24:

先交待下你研究的體系，以及系綜，然后是怎么做MD的。這樣才能分析原因。

我是先按照波爾滋蔓分布給定原子的初始速度之后開始給體系調(diào)溫,經(jīng)過一段時間的調(diào)溫后,體系在NVE系綜下運行.

不管調(diào)溫的步數(shù)有多少,只要停止調(diào)溫,系統(tǒng)的溫度就會慢慢升高,而且整個體系的能量頁步守恒.不知道什么原因回導致這樣的后果?

贊一下

3樓2010-09-12 14:59:49

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yahoohoo

鐵桿木蟲 (著名寫手)

模擬EPI: 10
應(yīng)助: 0 (幼兒園)
貴賓: 1.55
金幣: 7632.4
散金: 251
紅花: 4
帖子: 1176
在線: 167小時
蟲號: 74894
注冊: 2005-06-15
性別: GG
專業(yè): 理論和計算化學

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
lei0736(金幣+2):謝謝 2010-09-12 18:16:44

引用回帖:

Originally posted by weizhiyong at 2010-09-12 14:59:49:

我是先按照波爾滋蔓分布給定原子的初始速度之后開始給體系調(diào)溫,經(jīng)過一段時間的調(diào)溫后,體系在NVE系綜下運行.

不管調(diào)溫的步數(shù)有多少,只要停止調(diào)溫,系統(tǒng)的溫度就會慢慢升高,而且整個體系的能量頁步守恒 ...

如果時間步長較大，可以試試小步長，如果還是不行那就是程序的問題

4樓2010-09-12 16:51:42

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wyding

木蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 3579.2
散金: 633
紅花: 7
帖子: 370
在線: 386.3小時
蟲號: 1088971
注冊: 2010-09-04
性別: GG
專業(yè): 物理學II

★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zh1987hs(金幣+4):謝謝指點 2010-09-16 11:06:38

正如樓上，可能是算法的問題。只要時間步長足夠小（不能太小，否則模擬時間長度過�。琋VE系統(tǒng)能量是守恒的（上下有些波動，畢竟是數(shù)值模擬嘛）！NVT的步數(shù)也要足夠大，否則不能達到平衡！這樣NVE后溫度有10％之內(nèi)的波動都是容許的！還有算溫度的時候，各粒子速度要減去質(zhì)心的速度，MD要保持質(zhì)心不動，否則溫度沒有意義！就像運動火車中杯子里的開水，采用統(tǒng)計學算溫度時要減去火車的速度，就是這個道理！

[ Last edited by wyding on 2010-9-16 at 11:01 ]

向強者看齊，即使高度不夠也是昂首挺胸！

5樓2010-09-16 10:58:24

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

qphll

金蟲 (正式寫手)

模擬EPI: 10
應(yīng)助: 18 (小學生)
金幣: 2708.7
散金: 2294
紅花: 8
帖子: 554
在線: 227.1小時
蟲號: 89654
注冊: 2005-08-29
性別: GG
專業(yè): 碳素材料與超硬材料

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流

引用回帖:

Originally posted by wyding at 2010-09-16 10:58:24:
正如樓上，可能是算法的問題。只要時間步長足夠小（不能太小，否則模擬時間長度過小），NVE系統(tǒng)能量是守恒的（上下有些波動，畢竟是數(shù)值模擬嘛）！NVT的步數(shù)也要足夠大，否則不能達到平衡！這樣NVE后溫度有10％之 ...

這個回復很到位. 贊一下.

Life, Love, Laugh.

6樓2010-12-03 10:00:36

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zsuzwr

木蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
貴賓: 0.15
金幣: 2811.7
散金: 10
帖子: 554
在線: 275小時
蟲號: 209571
注冊: 2006-03-06
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zh1987hs(金幣+2): 謝謝 2011-02-20 21:31:49

我也碰到類似問題，出自摩擦項

Lala

7樓2011-02-18 15:28:48

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

weizhiyong

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1365.8
散金: 11
紅花: 1
帖子: 356
在線: 178.7小時
蟲號: 585996
注冊: 2008-08-08
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

★
御劍江湖(金幣+1): 謝謝回復 2011-04-09 09:55:35

引用回帖:

Originally posted by zsuzwr at 2011-02-18 07:28:48:
我也碰到類似問題，出自摩擦項

在NVE中，是不可能出現(xiàn)這個現(xiàn)象的，因為系統(tǒng)是孤立的，和外界絕熱。最后我發(fā)現(xiàn)是程序的問題。

8樓2011-02-20 03:31:15

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

relonfbq

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1294.3
散金: 1140
紅花: 2
帖子: 301
在線: 99.7小時
蟲號: 744333
注冊: 2009-04-09
專業(yè): 金屬材料的微觀結(jié)構(gòu)

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流

我的系綜是設(shè)想npt系統(tǒng)，通過p-r恒壓模型，再通過標定速度的方法標定溫度，可是標定了幾千步下來，停止標定溫度，溫度漲落還是很厲害？這是怎么一回事�。�

9樓2011-04-07 22:10:07

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

weizhiyong

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1365.8
散金: 11
紅花: 1
帖子: 356
在線: 178.7小時
蟲號: 585996
注冊: 2008-08-08
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

★
御劍江湖(金幣+1): 謝謝回復 2011-04-09 09:55:15

引用回帖:

Originally posted by relonfbq at 2011-04-07 22:10:07:
我的系綜是設(shè)想npt系統(tǒng)，通過p-r恒壓模型，再通過標定速度的方法標定溫度，可是標定了幾千步下來，停止標定溫度，溫度漲落還是很厲害？這是怎么一回事��？

據(jù)我現(xiàn)在所知，一般是勢函數(shù)錯了的問題，原子由晶體變成非晶體轉(zhuǎn)變所致