Event-chain Monte Carlo with factor fields.

We study the dynamics of one-dimensional (1D) interacting particles simulated with the event-chain Monte Carlo algorithm (ECMC). We argue that previous versions of the algorithm suffer from a mismatch in the factor potential between different particle pairs (factors) and show that in 1D models, this mismatch is overcome by factor fields. ECMC with factor fields is motivated, in 1D, for the harmonic model, and validated for the Lennard-Jones model as well as for hard spheres. In 1D particle systems with short-range interactions, autocorrelation times generally scale with the second power of the system size for reversible Monte Carlo dynamics, and with its first power for regular ECMC and for molecular dynamics. We show, using simulations, that the autocorrelation time grows only with the square root of the system size for ECMC with factor fields. Mixing times, which bound the time to reach equilibrium from an arbitrary initial configuration, grow with the first power of the system size.

[1]  B. Alder,et al.  Decay of the Velocity Autocorrelation Function , 1970 .

[2]  A C Maggs,et al.  Multiscale Monte Carlo algorithm for simple fluids. , 2006, Physical review letters.

[3]  Koji Hukushima,et al.  Dynamics of One-Dimensional Ising Model without Detailed Balance Condition , 2013 .

[4]  K. Kimura,et al.  Anomalous diffusion analysis of the lifting events in the event-chain Monte Carlo for the classical XY models , 2017, 1709.01665.

[5]  Stefan Schaefer,et al.  Testing the event-chain algorithm in asymptotically free models , 2018, Physical Review D.

[6]  Werner Krauth,et al.  Event-chain Monte Carlo for classical continuous spin models , 2015, 1508.06541.

[7]  Werner Krauth,et al.  Cell-veto Monte Carlo algorithm for long-range systems. , 2016, Physical review. E.

[8]  E A J F Peters,et al.  Rejection-free Monte Carlo sampling for general potentials. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  Werner Krauth,et al.  All-atom computations with irreversible Markov chains. , 2018, The Journal of chemical physics.

[10]  Sharon C Glotzer,et al.  Hard-disk equation of state: first-order liquid-hexatic transition in two dimensions with three simulation methods. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[11]  P. Charbonneau,et al.  Clustering and assembly dynamics of a one-dimensional microphase former. , 2018, Soft matter.

[12]  Werner Krauth,et al.  Two-dimensional melting: from liquid-hexatic coexistence to continuous transitions. , 2014, Physical review letters.

[13]  W. Krauth.  Statistical Mechanics: Algorithms and Computations , 2006 .

[14]  Cyril Banderier,et al.  Local time for lattice paths and the associated limit laws , 2018, GASCom.

[15]  Werner Krauth,et al.  Irreversible Local Markov Chains with Rapid Convergence towards Equilibrium. , 2017, Physical review letters.

[16]  Werner Krauth,et al.  Two-step melting in two dimensions: first-order liquid-hexatic transition. , 2011, Physical review letters.

[17]  Feller William,et al.  An Introduction To Probability Theory And Its Applications , 1950 .

[18]  Michael Chertkov,et al.  Irreversible Monte Carlo Algorithms for Efficient Sampling , 2008, ArXiv.

[19]  D. Wallace,et al.  Thermoelastic Theory of Stressed Crystals and Higher-Order Elastic Constants , 1970 .

[20]  L. Tonks.  The Complete Equation of State of One, Two and Three-Dimensional Gases of Hard Elastic Spheres , 1936 .

[21]  Werner Krauth,et al.  Event-chain Monte Carlo algorithms for hard-sphere systems. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[22]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[23]  Ze Lei,et al.  Irreversible Markov chains in spin models: Topological excitations , 2017 .

[24]  Jan Kierfeld,et al.  Event-chain Monte Carlo algorithms for three- and many-particle interactions , 2016, 1611.09098.

[25]  P. Fearnhead,et al.  Piecewise deterministic Markov processes for scalable Monte Carlo on restricted domains , 2017, 1701.04244.

[26]  H. G. Petersen,et al.  Error estimates on averages of correlated data , 1989 .

[27]  Peter Winkler,et al.  Mixing Points on a Circle , 2005, APPROX-RANDOM.

[28]  Werner Krauth,et al.  Event-chain algorithm for the Heisenberg model: Evidence for z≃1 dynamic scaling. , 2015, Physical review. E, Statistical, nonlinear, and soft matter physics.

[29]  Ze Lei,et al.  Mixing and perfect sampling in one-dimensional particle systems , 2018, EPL (Europhysics Letters).

[30]  J. Wittmer,et al.  Scale-free center-of-mass displacement correlations in polymer melts without topological constraints and momentum conservation: a bond-fluctuation model study. , 2011, The Journal of chemical physics.

[31]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[32]  Peter Winkler,et al.  Mixing Points on an Interval , 2005, ALENEX/ANALCO.

[33]  Elizabeth L. Wilmer,et al.  Markov Chains and Mixing Times , 2008 .

[34]  B. Alder,et al.  Phase Transition for a Hard Sphere System , 1957 .

[35]  Koji Hukushima,et al.  Eigenvalue analysis of an irreversible random walk with skew detailed balance conditions. , 2015, Physical review. E.

[36]  T. R. Scavo,et al.  On the Geometry of Halley's Method , 1995 .

[37]  Werner Krauth,et al.  Generalized event-chain Monte Carlo: constructing rejection-free global-balance algorithms from infinitesimal steps. , 2013, The Journal of chemical physics.