| 6 | 1/1 | 返回列表 |
| 查看: 1427 | 回復(fù): 5 | ||||
[交流]
分子動(dòng)力學(xué)模擬中有可能出現(xiàn)混沌現(xiàn)象嗎?
|
|
我不是學(xué)物理的,甚至大學(xué)物理都沒有學(xué)過。關(guān)于混沌(CHAOS)我只是從一些科普書籍或是紀(jì)錄片中有所了解,我對(duì)它印象很深,因?yàn)檫@是個(gè)很“酷”的概念——完全決定性的系統(tǒng)也可以出現(xiàn)無法預(yù)測的現(xiàn)象。那么像分子動(dòng)力學(xué)這樣一個(gè)基于經(jīng)典力學(xué)的模擬過程,是否也可能出現(xiàn)結(jié)果迥異的混沌現(xiàn)象呢? 之所以想起這個(gè)問題是因?yàn)槲以谀M過程中發(fā)現(xiàn)同樣一個(gè)體系,同樣的模擬參數(shù),最后得出了完全不同的結(jié)果。我認(rèn)為這應(yīng)該是初始賦予的原子速度不一樣所致,但是會(huì)不會(huì)和混沌有什么關(guān)系呢?解答這個(gè)問題需要有很好的物理基礎(chǔ),同時(shí)又理解分子動(dòng)力學(xué)的過程,希望板上有這樣的大牛! [ 來自科研家族 化學(xué)生物學(xué) ] |
» 搶金幣啦!回帖就可以得到:
+1/268
+1/191
+1/183
+5/100
+1/91
+1/91
+2/62
+1/33
+1/26
+1/21
+1/11
+1/11
+1/11
+1/10
+1/8
+1/6
+1/5
+1/4
+1/4
+1/3
|
推薦一篇文章,http://dx.doi.org/10.1021/jp953409x 可能是第一篇講MD中的混沌現(xiàn)象。 MD中的混沌,比較接近于分叉現(xiàn)象 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +503 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
經(jīng)典力學(xué)中,只要是三體或者三體以上的問題,都是混沌(chaotic)的,也就是說沒有解析解,而數(shù)值解由于誤差的存在,求解的粒子的軌跡與粒子實(shí)際軌跡的誤差隨著時(shí)間的增加會(huì)指數(shù)級(jí)的增長(Lyapunov Instability,參見Understanding Molecular Simulation, Chapter 4.3.4)。所以對(duì)于MD,除非你只求解兩個(gè)粒子的運(yùn)動(dòng),否則體系必然是混沌的。 但MD并不關(guān)心粒子在相空間中的單條軌跡的精確求解,MD求的是系綜平均,也就是能夠把系綜所包含的相空間的點(diǎn)都走遍,即ergodicity。對(duì)于平衡態(tài)的性質(zhì),只要能實(shí)現(xiàn)ergodicity,所求的平均值就是準(zhǔn)確的,與你粒子軌跡算得準(zhǔn)確不準(zhǔn)確沒關(guān)系,所以對(duì)于平衡態(tài)的性質(zhì),理論上講MD和MC的結(jié)果是沒區(qū)別的,區(qū)別在于采樣的過程。對(duì)于涉及動(dòng)力學(xué)的性質(zhì)(比如擴(kuò)散)的計(jì)算,由于MD計(jì)算得到的軌跡和真實(shí)情況下的軌跡是有誤差的,所以計(jì)算的動(dòng)力學(xué)性質(zhì)也是有誤差的,但這個(gè)誤差并不是無限大,而是有界的,因而MD計(jì)算的動(dòng)力學(xué)性質(zhì)還是有實(shí)際意義的。 |
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 349求調(diào)劑 +5 | 李木子啊哈哈 2026-03-25 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +6 | WW.' 2026-03-23 | 8/400 |
|
|
[考研] 291求調(diào)劑 +4 | Y-cap 2026-03-29 | 5/250 |
|
|
[考研] 289求調(diào)劑 +13 | 新時(shí)代材料 2026-03-27 | 13/650 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|