We generalize the rejection-free event-chain Monte Carlo algorithm from many particle systems with pairwise interactions to systems with arbitrary three- or many-particle interactions. We introduce generalized lifting probabilities between particles and obtain a general set of equations for lifting probabilities, the solution of which guarantees maximal global balance. We validate the resulting three-particle event-chain Monte Carlo algorithms on three different systems by comparison with conventional local Monte Carlo simulations: (i) a test system of three particles with a three-particle interaction that depends on the enclosed triangle area; (ii) a hard-needle system in two dimensions, where needle interactions constitute three-particle interactions of the needle end points; (iii) a semiflexible polymer chain with a bending energy, which constitutes a three-particle interaction of neighboring chain beads. The examples demonstrate that the generalization to many-particle interactions broadens the applicability of event-chain algorithms considerably.
[1]
Berend Smit,et al.
Understanding molecular simulation: from algorithms to applications
,
1996
.
[2]
R. Lipowsky,et al.
Stretching of semiflexible polymers with elastic bonds
,
2004,
The European physical journal. E, Soft matter.
[3]
H. Kleinert.
Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets
,
2006
.