 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 787 | 回復(fù): 9 | |||
nightwalker金蟲 (著名寫手)
|
[交流]
【求助】關(guān)于MC模擬中蛇行算法的能量計(jì)算
|
 各位用過(guò)Monte Carlo方法的蟲友： 我最近在用該方法加入蛇行算法來(lái)計(jì)算相分離一類的工作，結(jié)果遇到一個(gè)程序中的錯(cuò)誤。怎么都調(diào)試不過(guò)去。 我這個(gè)程序里面用到了蛇行算法，還用到了用于能量判斷的Metropolics規(guī)則。處理的是親水加疏水的三嵌段共聚物在溶液中的體系 結(jié)果發(fā)現(xiàn)： 不用能量判斷的時(shí)候，蛇行算法沒問題，程序可以運(yùn)行且無(wú)錯(cuò)誤。當(dāng)然，這種情況下相當(dāng)于無(wú)熱體系。 然后，再加入能量判斷但是取消蛇行的時(shí)候，系統(tǒng)可以分相。結(jié)果正確。 但是如果假如能量判斷而且也有蛇行的時(shí)候后，就不再分相了。結(jié)果當(dāng)然錯(cuò)誤。檢查發(fā)現(xiàn)這種處理方法下，能量也有波動(dòng)，但是就是不分相。 百思不得其解，特來(lái)求助。 說(shuō)明一下，我的能量判斷方法也是通過(guò)計(jì)算相互接觸對(duì)數(shù)的方法。 調(diào)試多日，無(wú)果，抓狂，跪求高人指點(diǎn)。 |
回復(fù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有6人回復(fù)
-
考研調(diào)劑
已經(jīng)有3人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有16人回復(fù)
-
318求調(diào)劑
已經(jīng)有5人回復(fù)
-
285求調(diào)劑
已經(jīng)有5人回復(fù)
-
306求調(diào)劑
已經(jīng)有3人回復(fù)
-
354求調(diào)劑
已經(jīng)有10人回復(fù)
-
303求調(diào)劑
已經(jīng)有5人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 相互作用的昨天和今天：從力、能量到信息 [當(dāng)然是轉(zhuǎn)載]
已經(jīng)有6人回復(fù)
- 求回歸方程參數(shù)估計(jì)的最優(yōu)化算法matlab代碼
已經(jīng)有9人回復(fù)
- 關(guān)于遺傳算法中交叉的程序
已經(jīng)有5人回復(fù)
- 【求助】關(guān)于Gaussian計(jì)算中的能量問題
已經(jīng)有14人回復(fù)
- 【轉(zhuǎn)帖】蒙特卡羅算法計(jì)算圓周率
已經(jīng)有12人回復(fù)
- 【求助】產(chǎn)物能量高于過(guò)渡態(tài)能量如何解釋啊
已經(jīng)有8人回復(fù)
- 【討論】MS中sorption模塊的MC計(jì)算
已經(jīng)有24人回復(fù)
- 【求助】MC方法計(jì)算吸附量
已經(jīng)有14人回復(fù)
- 【求助】關(guān)于吉布斯自由能的算法
已經(jīng)有17人回復(fù)
2樓2010-08-26 02:21:07
★ ★
lei0736(金幣+2):謝謝 2010-08-26 12:42:51
nightwalker(金幣+1): 2010-08-28 08:38:02
lei0736(金幣+2):謝謝 2010-08-26 12:42:51
nightwalker(金幣+1): 2010-08-28 08:38:02
3樓2010-08-26 02:21:50
yahoohoo
鐵桿木蟲 (著名寫手)
- 模擬EPI: 10
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 1.55
- 金幣: 7632.4
- 散金: 251
- 紅花: 4
- 帖子: 1176
- 在線: 167小時(shí)
- 蟲號(hào): 74894
- 注冊(cè): 2005-06-15
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
4樓2010-08-27 22:25:02
nightwalker
金蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3062.2
- 散金: 1095
- 紅花: 3
- 帖子: 1703
- 在線: 207.1小時(shí)
- 蟲號(hào): 1014530
- 注冊(cè): 2010-05-09
- 性別: GG
- 專業(yè): 高分子物理與高分子物理化
5樓2010-08-28 08:41:10
nightwalker
金蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3062.2
- 散金: 1095
- 紅花: 3
- 帖子: 1703
- 在線: 207.1小時(shí)
- 蟲號(hào): 1014530
- 注冊(cè): 2010-05-09
- 性別: GG
- 專業(yè): 高分子物理與高分子物理化
6樓2010-08-28 08:43:54
nightwalker
金蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3062.2
- 散金: 1095
- 紅花: 3
- 帖子: 1703
- 在線: 207.1小時(shí)
- 蟲號(hào): 1014530
- 注冊(cè): 2010-05-09
- 性別: GG
- 專業(yè): 高分子物理與高分子物理化
7樓2010-08-28 08:45:55
yahoohoo
鐵桿木蟲 (著名寫手)
- 模擬EPI: 10
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 1.55
- 金幣: 7632.4
- 散金: 251
- 紅花: 4
- 帖子: 1176
- 在線: 167小時(shí)
- 蟲號(hào): 74894
- 注冊(cè): 2005-06-15
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
★ ★ ★
nightwalker(金幣+1):謝謝參與
lei0736(金幣+2):謝謝 2010-08-29 14:55:15
nightwalker(金幣+1): 2010-08-31 11:04:33
nightwalker(金幣+1):謝謝參與
lei0736(金幣+2):謝謝 2010-08-29 14:55:15
nightwalker(金幣+1): 2010-08-31 11:04:33
|
可以先設(shè)計(jì)個(gè)實(shí)驗(yàn)測(cè)試你的這個(gè)算法是否有錯(cuò)。譬如，如果該嘗試MC運(yùn)動(dòng)被接受了，你可以計(jì)算新的能量$E_n = E_o + \delta E$，這里的$\delta E$就是用來(lái)判斷MC運(yùn)動(dòng)是否被接受的能量變化。另外，再重新算體系的總能量$E_tot$，如果兩者一樣那這個(gè)算法才可能是正確的。 |
8樓2010-08-28 17:25:25
yahoohoo
鐵桿木蟲 (著名寫手)
- 模擬EPI: 10
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 1.55
- 金幣: 7632.4
- 散金: 251
- 紅花: 4
- 帖子: 1176
- 在線: 167小時(shí)
- 蟲號(hào): 74894
- 注冊(cè): 2005-06-15
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
9樓2010-08-31 17:27:43
nightwalker
金蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3062.2
- 散金: 1095
- 紅花: 3
- 帖子: 1703
- 在線: 207.1小時(shí)
- 蟲號(hào): 1014530
- 注冊(cè): 2010-05-09
- 性別: GG
- 專業(yè): 高分子物理與高分子物理化
10樓2010-09-02 07:58:57
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 nightwalker 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑院校信息 +4 | CX 330 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +13 | 冬十三 2026-03-15 | 14/700 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +5 | 張vvvv 2026-03-15 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|