Un algorithme de simulation de l'interaction électrostatique est introduit, fondé sur un champ électrique fluctuant contraint par la loi de Gauss. Sa localité lui confère un coût linéaire en Monte Carlo.<br /><br />L'utilisation de conditions aux bords périodiques ne requiert pas de terme correctif, en contraste avec les méthodes usuelles. La simulation de systèmes quasi-bidimensionnels en est facilitée. Des surfaces quelconques à potentiel fixé peuvent être introduites.<br /><br />La loi de Gauss couple la dynamique des charges et du champ. La présence de charges libres accélère la relaxation du champ. La mobilité des charges est freinée par une « traînée de champ » qu'elles laissent en se déplaçant ; la modélisation quantitative de cet effet conduit à une méthode systématique pour le supprimer.