Archive for Glauber dynamics

Natural quantum MCMC

Posted in Books, Statistics, University life with tags , , , , , , , , , , , on November 10, 2025 by xi'an

Quite unexpectedly, I opened Nature of 16 October 2025 only to realise it featured a paper by Chi-Fang Chen et al. on a valid MCMC algorithm for quantum computing. The motivation is for simulating quantum multi-body physics, with details that escape me. The paper claims this is the first valid quantum MCMC quantum thermal simulation, with further application similar to the MCMC revolution in Bayesian inference. (An earlier proposal by Davies in 1974 cannot be implemented for many-bodies problems.)

“Although early proposals for simulating thermalization directly considered simulating the global system–bath Hamiltonian evolution, more tractable approaches use a master equation, or Lindbladian, to capture the effect of a large bath on the small system using a continuous-time Markovian process. However, previous approaches inherited issues of the prototypical Davies generator, which worked well in quantum optics but has unphysical features in noncommuting many-body systems because of the exponentially small level spacing. Our nature-inspired algorithm takes precisely the form of a Lindbladian, inherits the locality from the physical model, exactly satisfies detailed balance and resembles the interactions we expect from weak coupling to a Markovian thermal bath.”

The setup is one of a continuous-time quantum Markov chain (or Lindbladian) over a finite set of configurations. The arguments for validating the algorithm (into a completely spelled out theorem!) involve quantum detailed balance (illustrated by the above cartoon) and locality (which I understand as a form of Markovianity in the updating rule). The transition attached to the algorithm however remains incomprehensible to me, involving a Fourier transform of a Hermitian jump operator (both notions escaping me in this context).  And the illustration on an Ising model does not seem particularly quantum related, albeit the “Pauli operators” may be unrelated with the original spin binaries… There is also no MCMC reject step.

“Detailed balance (…) gives a coherent Gibbs sampler, in which we prepare the purified Gibbs state by a natural adiabatic path parametrized by the inverse temperature β. Of course, the purified Gibbs state reduces to the (mixed) Gibbs state once we trace out the purifying replica system, but the purification opens doors to advanced quantum algorithmic tools, such as verification (for example, through the swap test or measuring the energy of the parent Hamiltonian) or faster mean estimation.”

If I understand correctly the above “Gibbs sampler” relates to the Gibbs distribution as in the original 1984 paper of Geman and Geman! It may thus be that the current paper is similarly annunciating a new era in scientific computing, although I remain in the dark as to how to link it with a practical problem.

Martin Hairer gets Breakthrough Prize (and $3M)

Posted in Books, University life with tags , , , , , , , , , on September 14, 2020 by xi'an

Just heard the news that Fields Medallist Martin Hairer (formerly U of Warwick) got the 2021 Breakthrough Prize in Mathematics for his unification theory of stochastic partial differential equations, which he likens to a form of Taylor expansion in the massive Inventiones paper describing this breakthrough. (Looking at the previous winners of the prize, who also made its selection committee, this represents a break from focussing primarily on algebraic geometry! If not from sticking to male recipients…)

We introduce a new notion of “regularity structure” that provides an algebraic framework allowing to describe functions and/or distributions via a kind of “jet” or local Taylor expansion around each point. The main novel idea is to replace the classical polynomial model which is suitable for describing smooth functions by arbitrary models that are purpose-built for the problem at hand. In particular, this allows to describe the local behaviour not only of functions but also of large classes of distributions. We then build a calculus allowing to perform the various operations (multiplication, composition with smooth functions, integration against singular kernels) necessary to formulate fixed point equations for a very large class of semi-linear PDEs driven by some very singular (typically random) input. This allows, for the first time, to give a mathematically rigorous meaning to many interesting stochastic PDEs arising in physics. The theory comes with convergence results that allow to interpret the solutions obtained in this way as limits of classical solutions to regularised problems, possibly modified by the addition of diverging counterterms. These counterterms arise naturally through the action of a “renormalisation group” which is defined canonically in terms of the regularity structure associated to the given class of PDEs. Our theory also allows to easily recover many existing results on singular stochastic PDEs (KPZ equation, stochastic quantisation equations, Burgers-type equations) and to understand them as particular instances of a unified framework. One surprising insight is that in all of these instances local solutions are actually “smooth” in the sense that they can be approximated locally to arbitrarily high degree as linear combinations of a fixed family of random functions/distributions that play the role of “polynomials” in the theory. As an example of a novel application, we solve the long-standing problem of building a natural Markov process that is symmetric with respect to the (finite volume) measure describing the \Phi^4_ 3 Euclidean quantum field theory. It is natural to conjecture that the Markov process built in this way describes the Glauber dynamic of 3-dimensional ferromagnets near their critical temperature.