# Monte Carlo simulations

MathOverflow Asked by Alekk on November 7, 2020

I was wondering what were the models of statistical physics that are still considered difficult/slow to simulate (exactly, or approximately) with the current technology of Monte Carlo approaches. I can think of growing long self-avoiding walks, simulate very dense hard spheres configurations, etc… though I may be wrong since I am not aware of the new approaches to tackle these problems. It could be nice to have a brief description of the way people attack these difficult simulations.

Lattice gauge theory is indeed still very demanding, although the last decade has seen tremendous progress, both in terms of algorithms and of software implementations that make good use of modern computer hardware (GPUs & Co.). To some extent, lattice gauge theory simulations now serve as benchmarks for some HPC systems.

A hard problem is still posed by lattice gauge theories with fermions in the presence of a chemical potential (i.e. finite-density QCD), where there exists a "sign problem", i.e. where the determinant of the Dirac operator appearing in the path integral to be evaluated is no longer positive, whence the ordinary Markov-chain MC methods no longer apply. There are some attempts to bypass this using complex Langevin methods (but the theoretical foundation isn't 100% settled) or imaginary chemical potentials (but the continuation to real chemical potential only gets you so far).

Answered by gmvh on November 7, 2020

Here is a bit of a braindump, ranging from things physicists actually write down to theory that might never be seen by anyone in a lab. Please comment if you want some more specific responses, and I can try to hunt something down.

The large particle physics projects use intensive MC simulations, and these are not computationally quick or easy, taking days or weeks to run (my only direct experience is with SNO, and indir.ectly CERN because of that, but this must be true for most projects). These often seem very specific (i.e. they model much of the world, including where dirt might come from and so on, not something like 'hard discs' or other abstract systems)

As Steve Huntsman mentions, lattice guage theory is demanding (though there have been many huge successes!), both in terms of the theory and in terms of the crazy calculations people have succeeded in doing.

Outside of my direct experience, but within recent reading, random surfaces and the 'surface integrals' (which are in some sense analogous to 'path integrals', though I don't know enough physics to know why) are so hard that we are very far from making any sort of reasonable success. There are many good survey articles by physicists on the subject (google searching will get some), but as a mathematician it is worthwhile to remember that these models are still fairly mathematically intractable.

Statistical physics (of the graph-dynamics/spin-system type) is, as far as I know, a mixed bag right now. Until recently, 'all' of it was I think quite hard. Some recent advances (especially coupling from the past) have allowed sampling from the Ising model on a lattice and a number of other formerly intractable problems, while leaving other nearby problems sort of untouched.

I don't have any great references for the above, and am not involved in that research, so take it with a grain of salt. One source you should really look at is David Wilson's spectacular archive on MCMC in general. He seems to have a strong interest in physical models: http://dimacs.rutgers.edu/~dbwilson/exact.html/

Another general reference (it is not primarily about physics, but has physics references at the end, and was an intro to the area for me) is http://www.ams.org/bull/2009-46-02/S0273-0979-08-01238-X/S0273-0979-08-01238-X.pdf. In it, he suggests reading "Statistical Mechanics: Algorithms and Computations" by Krauth (2006). I have only looked at it briefly, but it certainly touches upon many modern problems, including those that you mention, and is probably fairly up-to-date.

EDIT: I'll add something that is probably obvious, just in case you are new, since it wasn't obvious to me when I first looked at this subject. The 'hard' examples you mention (like solid-balls-in-boxes) are the ones that are just on the borderline between being tractable math and total messes. The Ising model has now crossed over into 'tractable' land, as have some types of graph-colouring problems. On the other hand, the first thing I mention (use of MC for large particle physics experiments) is often completely intractable, and nobody can turn these sorts of things into math problems - physicists just do something that seems reasonable and run the simulation for a while. There's nothing wrong with that, but you should be aware that the types of 'mathy' hard sampling problems you mention are ones where there is some hope of rigorously showing that your sample is close to correct, and there are many more problems where rigorous and/or sharp analysis is essentially impossible.

Answered by user2282 on November 7, 2020

Lattice gauge theory is extremely demanding (four dimensions, small space and time scales, extremely complicated dynamics, etc.), and application-specific supercomputing infrastructure has been used for it. Monte Carlo simulations are used to demonstrate QCD phenomenology.

Answered by Steve Huntsman on November 7, 2020

## Related Questions

### triviality of homology with local coefficients

1  Asked on February 17, 2021

### Generators of sandpile groups of wheel graphs

1  Asked on February 16, 2021 by castor

### Conditions for pointwise convergence of indicators precomposed with uniformly continuous sequence

1  Asked on February 16, 2021 by bernard_karkanidis

### (Weakly) connected sets with large (out-)boundary

0  Asked on February 15, 2021

### Global minimum of sum of a non-convex and convex function, where minima of the non-convex function can be found

0  Asked on February 15, 2021 by proof-by-wine

### Hilbert class field of Quadratic fields

1  Asked on February 14, 2021 by sina

### The variety induced by an extension of a field

2  Asked on February 12, 2021 by federico-fallucca

### Why is this algebra called the q-Weyl algebra?

1  Asked on February 12, 2021 by gmra

### Cohomology of a simplicial abelian group $X_bullet$, where $S_n$ acts on $X_n$

1  Asked on February 12, 2021 by patrick-nicodemus

### Hyperplane arrangements whose regions all have the same shape

1  Asked on February 12, 2021 by christian-gaetz

### Twisted winding number

1  Asked on February 9, 2021 by jack-l

### Generalized moment problem for discrete distributions

1  Asked on February 8, 2021 by harharkh

### Convergence of Riemann’s Product representation of Xi

2  Asked on February 7, 2021 by mustafa-said

### Positive subharmonic functions with constant integral blowing up at boundary

1  Asked on February 7, 2021 by fozz

### Symmetry in Hardy-Littlewood k-tuple conjecture

1  Asked on February 6, 2021 by sylvain-julien

### Looking up the Mordell-Weil rank and generator(s) of a Weierstrass Equation

4  Asked on February 5, 2021 by john-r-ramsden

### A question on moduli space of Hitchin’s equations

1  Asked on February 5, 2021

### Writing papers in pre-LaTeX era?

28  Asked on February 3, 2021 by psihodelia

### Checking if Hochschild cohomology $mathit{HH}^2(A)=0$

2  Asked on February 2, 2021 by serge

### On a degenerate SDE in the unit ball

0  Asked on February 2, 2021