qyw 发表于 2023-10-5 09:49
请问这个MSD为啥会降低呢

没想明白，可能体系太复杂了

请问这个MSD为啥会降低呢

sobereva 发表于 2022-10-27 04:53
165原子的生物质分子扩散运动当然难了。而且体系里又只有这么一个这种分子，采样误差也容易比较大。
采样 ...

明白了，谢谢sob老师

165原子的生物质分子扩散运动当然难了。而且体系里又只有这么一个这种分子，采样误差也容易比较大。
采样不是说你动力学模拟的步长，而是计算MSD时候考虑的“有效”的数据量，诸如跑很多次取平均，或者一个很大的体系里有多种这样的分子做统计平均，都属于。轨迹保存频率（保存间隔）设得过高也没什么意义，相邻帧过于接近的话增加保存频率不会显著增加有效数据。
对于本来内在可重复性就差的数据（比如一个难扩散分子的MSD，体系里还就这一个分子，没法做同种分子的平均），不可能要求结果的可重复性好，你只能通过增加采样（诸如反复跑多次取平均，或者体系内纳入多种这种分子）来让你要研究的信息的统计误差较小、能说明问题。

sobereva 发表于 2022-10-26 15:26
跟体系很有关，没有具体描述没法说
诸如像小分子液体的情况，小分子扩散运动很容易，就很容易得到较平直的 ...

谢谢sob老师解答，我是算生物质分子在50%的乙醇溶液的扩散系数，该分子有165个原子，总共12648个原子，扩散难度不算难吧.
采样的话，有算过步长为0.5fs、1fs、2fs的，总时长3ns，每1ps输出；也算过步长为1fs，总时长3ns，每10fs、1ps、10ps输出的，最后跑出来的MSD重复性10fs相对较好一点。
sob老师，以我这样的体系，以步长为1fs，取样在三万步，即每100fs输出的话，这样的采样程度算是充分吗？

跟体系很有关，没有具体描述没法说
诸如像小分子液体的情况，小分子扩散运动很容易，就很容易得到较平直的MSD曲线，而且可重复性较好
如果是扩散困难的情况，诸如磷脂膜中的磷脂，MSD的可重现性和平直程度就没那么理想。但这也和被统计分子的数目、计算MSD的方式（是否考虑时间平均）有很大关系，采样不充分的情况此问题会很显著。在计算方式确实得当的前提下，遇到可重复性差的情况，也可以跑多次取统计平均。