force-bias Monte Carlo (fbMC)

The Uniform-acceptance force-bias Monte Carlo (fbMC) method (Phys. Rev. B 85, 134301 (2012)) has been implemented in ADF ReaxFF by the group of Erik Neyts.

Usually, ReaxFF optimizations are alternated with fbMC steps to accelerate events such as reactions. The acceleration (MC step) is controlled by the temperature and maximum displacement, which should be balanced to keep close to the real dynamics. fbMC does not scale time uniformly (J. Chem. Phys. 141, 204104 (2014)). Diffusion coefficients and reaction rates should probably use quartic mass scaling.

In the developers version of ADF ReaxFF, fbMC mass scaling can be set. It defaults to quadratic scaling as in the original fbMC implementation. Setting imcroo to 4 would mean quartic mass scaling.

Input keywords

These keywords should be set at the end of the control file. fbMC control parameters can also be set in the GUI.

  • drmax – maximum displacement
  • imcfrq – frequency of fbMC steps
  • imcstp – number of fbMC steps

A worked out example for fbMC versus unaccelerated ReaxFF MD for healing graphene has been made for the 2-day ReaxFF workshop.