应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Monte Carlo » 【求助】关于MC模拟中蛇行算法的能量计算
101/1返回列表
查看: 1058  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

nightwalker

金虫 (著名写手)

[交流] 【求助】关于MC模拟中蛇行算法的能量计算



各位用过Monte Carlo方法的虫友：
我最近在用该方法加入蛇行算法来计算相分离一类的工作，结果遇到一个程序中的错误。怎么都调试不过去。

我这个程序里面用到了蛇行算法，还用到了用于能量判断的Metropolics规则。处理的是亲水加疏水的三嵌段共聚物在溶液中的体系

结果发现：

不用能量判断的时候，蛇行算法没问题，程序可以运行且无错误。当然，这种情况下相当于无热体系。

然后，再加入能量判断但是取消蛇行的时候，系统可以分相。结果正确。

但是如果假如能量判断而且也有蛇行的时候后，就不再分相了。结果当然错误。检查发现这种处理方法下，能量也有波动，但是就是不分相。

百思不得其解，特来求助。

说明一下，我的能量判断方法也是通过计算相互接触对数的方法。

调试多日，无果，抓狂，跪求高人指点。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
奇文共欣赏，疑义相与析。
1楼 2010-08-24 10:42:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
★
lei0736(金币+1):谢谢 2010-08-26 12:42:46
nightwalker(金币+1): 2010-08-28 08:37:54
你这个跟我的那个 嵌段共聚物的程序差不多吧？
赞一下(1人) 回复此楼
好好学习，天天向上。
2楼2010-08-26 02:21:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
★ ★
lei0736(金币+2):谢谢 2010-08-26 12:42:51
nightwalker(金币+1): 2010-08-28 08:38:02
http://muchong.com/bbs/viewthread.php?tid=1921541&fpage=1
赞一下(1人) 回复此楼
好好学习，天天向上。
3楼2010-08-26 02:21:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)
★ ★
nightwalker(金币+1): 2010-08-28 08:38:06
lei0736(金币+2):谢谢 2010-08-29 14:54:53
引用回帖:
Originally posted by nightwalker at 2010-08-24 10:42:54:


各位用过Monte Carlo方法的虫友：
我最近在用该方法加入蛇行算法来计算相分离一类的工作，结果遇到一个程序中的错误。怎么都调试不过去。

我这个程序里面用到了蛇行算法，还用到了用于能量判断的Metrop ...

你所谓的蛇行是 REPTATION （剪断一端的单体或Bead，接到同一条链的另一端）？

如果是，那这种算法并不适用于嵌段共聚物，因为链的结构会被改变；如果不是，请先说明你的具体算法是怎样执行的。
赞一下(1人) 回复此楼
4楼2010-08-27 22:25:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nightwalker

金虫 (著名写手)
如楼上所说   确实是 reptation的方式

举例说明就是   加入原来的  5个链节  和一个空格 0 的顺序是

5  4  3  2  1   0

蛇行结束 之后  其顺序 就变成了

0  5  4  3  2  1


在此我请教楼上一句   所谓的改变链结构是什么意思

能否详细说明一下   

感谢帮助
赞一下 回复此楼
奇文共欣赏，疑义相与析。
5楼2010-08-28 08:41:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nightwalker

金虫 (著名写手)
sorry   刚才没仔细看   跟3楼 所说的 reptation  方式不一样

我采取的蛇行方式 是 头或尾 向前走一步  

后面的链节 依次跟进  没有改变链结构

所以我才认为  找不到错误  但是有确定结果是错误的
赞一下 回复此楼
奇文共欣赏，疑义相与析。
6楼2010-08-28 08:43:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nightwalker

金虫 (著名写手)
感谢 1楼 2楼 的交流程序  

我这就下来看看 学习一下   呵呵
赞一下 回复此楼
奇文共欣赏，疑义相与析。
7楼2010-08-28 08:45:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)
★ ★ ★
nightwalker(金币+1):谢谢参与
lei0736(金币+2):谢谢 2010-08-29 14:55:15
nightwalker(金币+1): 2010-08-31 11:04:33
引用回帖:
Originally posted by nightwalker at 2010-08-28 08:43:54:
sorry   刚才没仔细看   跟3楼 所说的 reptation  方式不一样

我采取的蛇行方式 是 头或尾 向前走一步  

后面的链节 依次跟进  没有改变链结构

所以我才认为  找不到错误  但是有确定结果是错误的

可以先设计个实验测试你的这个算法是否有错。譬如，如果该尝试MC运动被接受了，你可以计算新的能量$E_n = E_o + \delta E$，这里的$\delta E$就是用来判断MC运动是否被接受的能量变化。另外，再重新算体系的总能量$E_tot$，如果两者一样那这个算法才可能是正确的。
赞一下(2人) 回复此楼
8楼2010-08-28 17:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)
★ ★ ★
ghcacj(金币+3):谢谢 2010-08-31 17:45:21
nightwalker(金币+1): 2010-09-02 07:56:34
引用回帖:
Originally posted by nightwalker at 2010-08-24 10:42:54:


各位用过Monte Carlo方法的虫友：
我最近在用该方法加入蛇行算法来计算相分离一类的工作，结果遇到一个程序中的错误。怎么都调试不过去。

我这个程序里面用到了蛇行算法，还用到了用于能量判断的Metrop ...

一个问题：格子MC所考虑的体系不可能太大，而且MC运动的效率不会高。为何不考虑 用DPD呢？既可以处理连续的大体系，观察长时间的动力学过程，而且溶剂存在下的流体力学相互作用也可以考虑进来，更切合实际 。
赞一下(1人) 回复此楼
9楼2010-08-31 17:27:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nightwalker

金虫 (著名写手)
目前我只能先用MC方法了，DPD还没学会  呵呵
赞一下 回复此楼
奇文共欣赏，疑义相与析。
10楼2010-09-02 07:58:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nightwalker 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 13:22 by jzdxtiUmJ6PL
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 8/400 2026-08-04 13:15 by jzdxtiUmJ6PL
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 jRl3mE6ddZGq 2026-08-03 11/550 2026-08-04 12:58 by jzdxtiUmJ6PL
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 纯娱乐，不喜欢勿喷 +5 Tide man 2026-08-04 8/400 2026-08-04 11:43 by everxc
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 6/300 2026-08-04 11:34 by MWRYsWbJw4DE
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 11:31 by MWRYsWbJw4DE
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 4/200 2026-08-04 11:26 by MWRYsWbJw4DE
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[有机交流] 一个有机合成实验室都需要哪些设备？ 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 微信指数没变化，科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
关闭
请填处理意见
关闭