24小时热门版块排行榜    

查看: 7247  |  回复: 101
本帖产生 3 个 模拟EPI ,点击这里进行查看

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+2):谢谢 2010-11-16 16:29:49
非常感谢,碰到高手了,其实我知道你以前是LXH老师的学生。你说到的NH3吸附在Cu-BTC上我也想尝试,请问你可否把文献列出来,我们一起做做看,碰到高手了,自然要多学习,这样才能成长的快。
71楼2010-11-14 19:47:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by ghcacj at 2010-11-14 19:47:31:
非常感谢,碰到高手了,其实我知道你以前是LXH老师的学生。你说到的NH3吸附在Cu-BTC上我也想尝试,请问你可否把文献列出来,我们一起做做看,碰到高手了,自然要多学习,这样才能成长的快。

此贴非常值得关注。。。
好好学习,天天向上。
72楼2010-11-16 10:52:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

引用回帖:
Originally posted by zyj8119 at 2010-11-16 10:52:48:

此贴非常值得关注。。。

没有下文了?
好好学习,天天向上。
73楼2010-11-22 13:58:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主


小木虫(金币+0.5):给个红包,谢谢回帖交流
是啊,没有下文了,可能LZ最近比较忙
74楼2010-11-22 15:39:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

引用回帖:
Originally posted by ghcacj at 2010-11-22 15:39:10:
是啊,没有下文了,可能LZ最近比较忙

三个老师在组里访问. 过去的一周, 没有做成任何事情....

Next week will be better? :-)
Life, Love, Laugh.
75楼2010-11-23 00:43:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
ghcacj(金币+10):精彩 2010-11-23 09:19:17
引用回帖:
Originally posted by zyj8119 at 2010-11-22 13:58:47:

没有下文了?

小update一下.

常常看到文献上有这样的说法:

"At least 20 million trials were used in the simulations, among which the first half was used for equilibrium, and the last half was used to calculate the ensemble averages."

那么这里的"20 million trials" 和 "50%" 这两个关键词是怎样来的?

同样用我们这里的例子, 同样用类似前面某次update提到的体系能量变化图, 我们再仔细看看:

上图!



大家可以看到, 我这里也是跑了 20 million, 也是将前50%扔掉, 后50%系综平均.

但是很明显, 从大图(黑线)可以看出, 体系的总能量在(10E6~14E6)这个阶段还是有很明显的降低. 从小图1(红线, 18E6~20E6)看出, 就算是将总迭代(2E7)的80%扔掉, 体系的能量在最后的20%那里还是在降低.

从小图2(蓝线, 19E6~20E6) 看起来, 似乎体系还是在波动, 但是考虑到波动的绝对值小于 20KJ (注意, 这里是 20KJ for the whole system, NOT 20KJ/mol), 再结合大图(黑线)最后 19E6~20E6处体现出来的'平坦', 我们有理由说服自己, 在19E6步以后, 体系已经达到平衡态了. 差不多可以安全地做系综平均了.

但是, 为了百分百放心, 最好在这个基础上, 再往后算10E6步, 然后看看这个区间 19E6~30E6 的能量波动, 如果能量的波动像小图2(蓝线)所示范围震荡, 那么OK, 放心大胆系综平均, 得出你的结论吧. 反之, 如果体系能量还是像小图1(红线)所示, 有一路下降的趋势(哪怕斜率平坦), 那么对不起, 体系还没有平衡, 继续跑吧......


以上小小update, 可以看作如何选择迭代总步数, 如果选择系综平均(到底扔掉前面的百分之多少)的一些依据.

[ Last edited by qphll on 2010-11-23 at 07:02 ]
Life, Love, Laugh.
76楼2010-11-23 06:56:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by qphll at 2010-11-23 06:56:08:



小update一下.

常常看到文献上有这样的说法:

"At least 20 million trials were used in the simulations, among which the first half was used for equilibrium, and the last half was used  ...

谢谢楼主继续跟进。。。
好好学习,天天向上。
77楼2010-11-23 08:30:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+3):谢谢 2010-11-23 12:34:19
引用回帖:
Originally posted by qphll at 2010-11-23 06:56:08:



小update一下.

常常看到文献上有这样的说法:

"At least 20 million trials were used in the simulations, among which the first half was used for equilibrium, and the last half was used  ...

不过目前我看到过的文献,计算步数最大一般是4E7了,似乎很难有更大的,这个感觉不光和计算体系有关,也和硬件配置相关,几年前几百万步一样可以被接受,现在似乎一定要1E7以上了。这个就和计算能接受的误差有关,对于贴出这个图,如果能接受的误差比较大,那么1KJ的能量波动也可以接受,当然有人更喜欢精确的数据。

[ Last edited by ghcacj on 2010-11-23 at 10:17 ]
78楼2010-11-23 09:20:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

★ ★ ★
zh1987hs(金币+3): 谢谢 2010-11-23 12:46:03
引用回帖:
Originally posted by ghcacj at 2010-11-23 09:20:27:


不过目前我看到过的文献,计算步数最大一般是4E7了,似乎很难有更大的,这个感觉不光和计算体系有关,也和硬件配置相关,几年前几百万步一样可以被接受,现在似乎一定要1E7以上了。这个就和计算能接受的误差有 ...

几年前, 实验对于表面,界面和受限孔内的现象几乎束手无策, 理论计算怎么说怎么有理. 现在仪器和实验技术的发展, 要求理论计算能提供更加精细的信息, 不单单满足于往左还是往右这样的趋势, 而且要求能说出往左多少度这样的信息. 所以计算精度要求的提高,也是势在必行.
Life, Love, Laugh.
79楼2010-11-23 12:42:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主


小木虫(金币+0.5):给个红包,谢谢回帖交流
很好,LZ上次说模拟Cu-BTC吸附NH3的例子,有文献吗?
80楼2010-11-23 13:56:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qphll 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +16 基诺咪客 2026-08-17 20/1000 2026-08-20 12:44 by 201120518022
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 言之凿凿,8月21号(明天)放榜 +12 医学老男孩 2026-08-20 14/700 2026-08-20 10:53 by ran2233
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 21/1050 2026-08-20 10:37 by Vivilian
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见