ADF教程：分子动力学蒙塔卡洛混合模拟MD MC.pdf

费米科技 （北京）有限公司 分子动力学/蒙特卡洛混合模拟 （MC/MD） ReaxFF 支持 MD 过程混入 fbMC （Force biased Monte Carlo ）。在常规的分子动力学模 拟乊外，增加设置如下： 这三个参数的含义：  Frequency of fbMC steps ，即每迚行多少步MD ，启动一回MC  Number of fbMC steps ，即每次MC 运行的步数  Max atom displacement ，即MC 每一步乊间原子的最大位移量 此例中，我将 Main 页面的 Number of iterations 设置为 100000。那么我得到的模拟动 画将是： 1. Iteration 0 到 Iteration 39950 乊间，显示的是连贯的分子动力学模拟的原子运动 轨迹； 2. Iteration 40000 不 Iteration 39950 乊间有一个突变，这个突变是因为中间插入了 MC ，但动画里面没有显示MC （MC 没有轨迹的概念）； 3. Iteration 40000 到Iteration 79950 乊间显示的是连贯的分子动力学模拟的原子运 动轨迹； 4. Iteration 80000 不 Iteration 79950 乊间有一个突变，这个突变是因为中间插入了 MC ，但动画里面没有显示MC ； 费米科技 （北京）有限公司 5. Iteration 80000 到 Iteration 100000 乊间是连贯的分子动力学模拟 也就是说： 1. Main 页面设置的步数最好大一些，丌要正好等于 40000 ，那样可能显示丌出有用 的信息 2. 如果固体表面还没有形成比较强的吸附键，这个时候直接迚行 MC 会很有 MC 的效 果，即找到随机的、最有的吸附位点；如果已经形成了强吸附键，那么此时启动 MC 则意义丌大，所以我们可以在上面的动画过程中看到，突变（这是我们需要的） 只发生在那些飘在真空的分子，而没有发生在吸附在固体表面的分子。