Au cours d’une simulation monte carlo, le système se déplace dans l’espace des phases en obéissant à la probabilité de transition W(x_i® x_j)=e^{-b D H} où D H=H_i-H_j et x_i représente une certaine configuration. Or, selon l’algorithme de hubb2_3.f, deux configurations consécutives ne peuvent différer que par l’orientation d’un spin. Ainsi, il est possible qu’un sous-espace F de l’espace de phase W soit tel que W(x_i® x_j)£ e " x_iÎ F et x_jÎ W \F où x_i et x_j ne diffèrent que par un spin, avec e < < 1. Dans ce cas, le système restera coincé dans le sous-espace F ce qui biaisera les statistiques.

Afin de remédier à ce problème, nous avons recours à l’option de mouvement globaux qui, après un certain nombre de mesures, flip le spin des électrons sur un site et cela pour toutes les tranches de temps, quittant ainsi probablement le sous-espace F .

Retour à la page de thermalisation