Nearly accurate solutions for Ising-like models using Maximal Entropy Random Walk

12/31/2019
by   Jarek Duda, et al.
0

While one-dimensional Markov processes are well understood, going to higher dimensions there are only a few analytically solved Ising-like models, in practice requiring to use costly, uncontrollable and inaccurate Monte-Carlo methods. There is discussed analytical approach for approximated problem exploiting Hammersley-Clifford theorem, which allows to realize generation of random Markov field through scanning line-by-line using local statistical model as in lossless image compression. While conditional distributions for such statistical model could be found with Monte-Carlo methods, there is discussed use of Maximal Entropy Random Walk (MERW) to obtain it from approximation of lattice as infinite in one direction and finite in the remaining. Specifically, in the finite directions we build alphabet of all patterns, find transition matrix containing energy for all pairs of such patterns, from its dominant eigenvector getting probability distribution of pairs in their infinite 1D Markov process, which can be translated to statistical model for line-by-line scan. Such inexpensive models, requiring seconds on a laptop for attached implementation and directly providing probability distributions of patterns, were tested for mean entropy and energy, getting maximal ≈ 0.02 error from analytical solution near critical point, which quickly improves to extremely accurate e.g. ≈ 10^-10 error for coupling constant J≈ 0.2.

READ FULL TEXT
research
07/21/2022

Metropolis Monte Carlo sampling: convergence, localization transition and optimality

Among random sampling methods, Markov Chain Monte Carlo algorithms are f...
research
08/20/2019

A Fast Method to Calculate Hitting Time Distribution for a Random Walk on Connected Undirected Graph

With the advent of increasingly larger graphs, we need to find a quick a...
research
12/19/2016

Random Walk Models of Network Formation and Sequential Monte Carlo Methods for Graphs

We introduce a class of network models that insert edges by connecting t...
research
02/26/2022

A Log-Gaussian Cox Process with Sequential Monte Carlo for Line Narrowing in Spectroscopy

We propose a statistical model for narrowing line shapes in spectroscopy...
research
08/23/2021

Conditional sequential Monte Carlo in high dimensions

The iterated conditional sequential Monte Carlo (i-CSMC) algorithm from ...
research
03/22/2019

An n-dimensional Rosenbrock Distribution for MCMC Testing

The Rosenbrock function is an ubiquitous benchmark problem for numerical...

Please sign up or login with your details

Forgot password? Click here to reset