青青青 发表于 2019-6-30 16:57 取样距离 |
k64_cc 发表于 2019-3-29 15:05 请问您说的dt是什么呢？ |
|
本帖最后由 青青青 于 2019-6-26 14:01 编辑 终于看到有写MSD的单位的了 |
| 可以用速度自相关函数来算啊 |
青青青 发表于 2019-3-28 16:41 比如说你有100个frame，dt取80 frame，是不是就只有20个data让你average了…… |
本帖最后由 青青青 于 2019-4-1 11:49 编辑 k64_cc 发表于 2018-8-16 16:59 您好，请问为什么MSD后段sample少呢？后段统计误差为啥会大呢？谢谢解答 |
|
本帖最后由 k64_cc 于 2018-8-16 17:01 编辑 等下你dt都100ns了……总数据量不得几倍于此，可怕。 MSD后段sample本来就少，不要也罢。前半段不是挺棒的嘛…… |
winterzen 发表于 2018-8-15 00:29 你说的很正确。其实之前组里的人用基于TST的方法做过类似体系，我这里先试了最简单的MD是因为我的材料里笼之间的窗口并不小（和之前的工作相比），所以我认为可能能垒不那么高，可以先用MD试试。然而看来我还是低估了烷烃的吸附能力，正己烷这样大的分子（虽然非极性）还是会比较强地吸附在特定位点上。 |
sobereva 发表于 2018-8-14 18:14 我也觉得和统计误差有关，因为我发现吸附量高时线性很好，而吸附量极低（整个系统内只有几个烷烃分子）时几乎毫无线性可言，这种介于两者之间的情况不太好确定。谢谢sob |
xylz6188 发表于 2018-8-14 17:43 300ps？我这已经100ns了，你是说300ns吧... |
|
从你的MSD曲线估测你的体系的扩散系数在1E-11 m^2/s的数量级上，是比较慢的扩散体系；另外你的多孔材料的孔结构中是不是有些大笼，笼之间通过一些狭小窗口连接？分子在其中扩散时可能大部分时间都处于稳定的大笼中的吸附位上，扩散的发生需要分子翻过自由能垒高的窗口跳跃到其他吸附位来进行，这种扩散方式也可能导致MSD不够平滑，不容易通过Einstein方程计算扩散系数。有一些基于TST处理rare event的方法可以处理这类慢扩散问题，可以参考 Phys. Rev. Lett., 93, 248301 Phys. Rev. Lett., 96, 044501 JCP, 122, 224712 JPCC, 120, pp 1110–1120 |
| 参与人数Participants 1 | eV +2 | 收起 理由Reason |
|---|---|---|
|
| + 2 |
| MSD曲线末端的统计误差比较大，往往因此趋势也和前面不同，统计前70%~80%左右就行了。当前这样也还可以了，如果觉得图像不理想，可以重新跑，直到跑出比较理想的曲线 |
| 在此基础上模拟一个300ps试一试 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-18 09:59 , Processed in 0.177139 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites