What progress has been made towards sampling from the 2D lattice Ising model with the following Hamiltonian:
$H=-J\sum_{\langle i,j \rangle}S_iS_j - \sum_i b_iS_i$
Where the first sum runs over all nearest-neighbour positions on a lattice. Especially, my emphasis is on very large lattice sizes where Markov-chain Monte Carlo methods become unwieldy.