计算化学公社

标题: 扩散系数沿某个坐标的规律变化 [打印本页]

作者
Author: diaok 时间: 2016-7-1 18:09
标题: 扩散系数沿某个坐标的规律变化
RT
计算扩散系数沿某个坐标的分布是否可行？

我研究了某个体系中离子的扩散，体系中央带的电荷显然会对他产生影响。

最初的想法是算出所有时刻所有粒子的扩散系数，然后投影在距离上

然而msd的计算需要长时间以及大量粒子平均才能拟合，
很短的时间内msd和时间是线性相关，而没有平方的关系，
如果把体系分成很多块，粒子的穿梭太快，统计时间还是不够。

尝试了下速度分布，在热浴条件下好像看不出变化
用什么方法能够体现这种影响呢？

对应粒子的速度、实际坐标以及投影距离都拿到了，但是却没法分析。

作者
Author: sobereva 时间: 2016-7-2 01:59
这个应该不可行，相比之下讨论不同坐标区域内的粒子沿着某个一方向的平均速度之类倒是更可行

作者
Author: tobeant 时间: 2016-7-4 14:19
之前看过一篇文献，里面计算了固液界面法向方向的扩散系数分布，是把体系分块来计算的。
粒子穿梭太快的话，计算每个块内粒子穿过所需的时间，然后做个时间的分布，不知行不行？

作者
Author: sobereva 时间: 2016-7-4 16:08

tobeant 发表于 2016-7-4 14:19
之前看过一篇文献，里面计算了固液界面法向方向的扩散系数分布，是把体系分块来计算的。
粒子穿梭太快的话 ...

对于扩散系数，若令统计可行，就不得不把块取大点
如果是粒子明显是往单一方向运动，按你说的统计穿过时间好办得多

作者
Author: fhh2626 时间: 2016-7-4 23:33
不确定你是不是想要这样：利用harmonic restrain类似的方法限制体系的反应坐标，然后跑平衡，算出该坐标下的扩散系数

作者
Author: diaok 时间: 2016-7-6 11:49
非常感谢各位讨论提供思路
我先尝试下这种穿梭时间的方法

作者
Author: 萌新杨同学 时间: 2023-2-11 17:56
请问楼主解决了计算不同位置扩散系数的问题了吗？方便分享一下嘛

欢迎光临计算化学公社 (http://bbs.keinsci.com/)