Lammps模拟,velocity中随机数变化,结果为什么也发生变化
velocity all create 300.0 123456789 rot yes dist gaussian
随机种子用12345678和123456789运行后结果不一样 返回小木虫查看更多
今日热帖
velocity all create 300.0 123456789 rot yes dist gaussian
随机种子用12345678和123456789运行后结果不一样 返回小木虫查看更多
沙发
不同随机数肯定会有点差别啊,相空间取样的问题
既然有差别,那随机数设置还有什么意义
最终得到的系统统计平均值比如温度,压力这些一样就可以了啊。分子动力学计算你看单个原子的结果是没意义的
但是我的模拟中随机数设定不同,结果完全不同
你贴的那条语句的含义是,给系统原子赋予随机初速度,使得统计意义下系统的初始温度为300K。不同的随机数生成的原子初速度可能不会相同,但是它们的体系温度肯定都会是300K,不信你让Lammps输出体系温度比对下。当然有的问题需要明确指定原子的初速度。比如做冲击相变模拟,需要给左端五层原子一定初速度让它去压缩右边原子,从而产生冲击波的传播。这时你要选定的原子作为一个region,然后直接施加速度就完了。所以,要么你的问题依赖于具体原子的初速度,那你就直接给原子指定初速度。如果只依赖于体系的温度,那么那样根据初始温度来给原子生成随机初速度没问题
,