参考讲义中相应部分内容教程下载链接

认为制造一个相当大的缺陷，并进行优化之后，得到缺陷结构：

设置参数如下：

其中

• Frequency of fbMC steps，表示每计算多少步分子动力学，就开始进行fbMC计算；
• Number of fbMC steps，表示每次fbMC做多少步；
• Max atom displacement，表示fbMC过程中，原子每一步的位置变化虽然是随机的，但不得超过0.17埃半径的圆球的范围，这个值越大，fbMC的随机性越好，效果较好，但是如果过大，则可能导致整个体系完全崩毁，从而导致模拟无效，因此用户可以自行测试，最小值一般为0.1，可以适当增大。

40000步后，得到大部分自愈修复的结构：

通过观察Movie我们发现，很快就修复成这样了，但是后面很长一段时间内，就无法进一步完成修复了。因此我们在最后一帧结构基础上，升高温度到1098K，并增大ax atom displacement到0.2，其它条件不变，发现非常迅速地，最后一个缺陷也自动修复了：