“…The method chooses the new configuration in the neighborhood of a current solution with unit probability if it has a lower energy and if it has higher energy than the current configuration, it will be accepted based on Boltzmann probability. For application to sophisticated systems with multimodal functionality and high-energy barriers between the potential energy minima in a system, several different modifications of the basic Metropolis MC simulations and simulated annealing have been devised (Beichl and Sullivan, 1999;Berg, 2003;Andricioaei and Straub, 1997;Iwamatsu and Okabe, 2000;Frantz, 2001;Ortiz et al, 2003). Most of these modifications incorporate appropriate changes in the basic acceptance probability function.…”