 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 786 | 回復(fù): 9 | |||
nightwalker金蟲 (著名寫手)
|
[交流]
【求助】關(guān)于MC模擬中蛇行算法的能量計(jì)算
|
 各位用過Monte Carlo方法的蟲友： 我最近在用該方法加入蛇行算法來計(jì)算相分離一類的工作，結(jié)果遇到一個(gè)程序中的錯(cuò)誤。怎么都調(diào)試不過去。 我這個(gè)程序里面用到了蛇行算法，還用到了用于能量判斷的Metropolics規(guī)則。處理的是親水加疏水的三嵌段共聚物在溶液中的體系 結(jié)果發(fā)現(xiàn)： 不用能量判斷的時(shí)候，蛇行算法沒問題，程序可以運(yùn)行且無錯(cuò)誤。當(dāng)然，這種情況下相當(dāng)于無熱體系。 然后，再加入能量判斷但是取消蛇行的時(shí)候，系統(tǒng)可以分相。結(jié)果正確。 但是如果假如能量判斷而且也有蛇行的時(shí)候后，就不再分相了。結(jié)果當(dāng)然錯(cuò)誤。檢查發(fā)現(xiàn)這種處理方法下，能量也有波動(dòng)，但是就是不分相。 百思不得其解，特來求助。 說明一下，我的能量判斷方法也是通過計(jì)算相互接觸對(duì)數(shù)的方法。 調(diào)試多日，無果，抓狂，跪求高人指點(diǎn)。 |
回復(fù)此樓» 猜你喜歡
一志愿華中科技大學(xué)071000，求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑院校信息
已經(jīng)有4人回復(fù)
-
085600材料與化工306
已經(jīng)有4人回復(fù)
-
286求調(diào)劑
已經(jīng)有10人回復(fù)
-
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有4人回復(fù)
-
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有8人回復(fù)
-
資源與環(huán)境 調(diào)劑申請(qǐng)(333分)
已經(jīng)有5人回復(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)物能量高于過渡態(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$就是用來判斷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ā)表 | |
|---|---|---|---|---|
|
[考研] 生物學(xué)調(diào)劑 +4 | Surekei 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +6 | 嘉年新程 2026-03-15 | 6/300 |
|
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|