Alexander Polyakov, a theoretical physicist now at Princeton University, caught a glimpse of the future of quantum theory in 1981. A range of mysteries, from the wiggling of strings to the binding of quarks into protons, demanded a new mathematical tool whose silhouette he could just make out.

“There are methods and formulae in science which serve as master keys to many apparently different problems,” he wrote in the introduction to a now famous four-page letter in *Physics Letters B*. “At the present time we have to develop an art of handling sums over random surfaces.”

Polyakov’s proposal proved powerful. In his paper he sketched out a formula that roughly described how to calculate averages of a wildly chaotic type of surface, the “Liouville field.” His work brought physicists into a new mathematical arena, one essential for unlocking the behavior of theoretical objects called strings and building a simplified model of quantum gravity.

Years of toil would lead Polyakov to breakthrough solutions for other theories in physics, but he never fully understood the mathematics behind the Liouville field.

Over the last seven years, however, a group of mathematicians has done what many researchers thought impossible. In a trilogy of landmark publications, they have recast Polyakov’s formula using fully rigorous mathematical language and proved that the Liouville field flawlessly models the phenomena Polyakov thought it would.

“It took us 40 years in math to make sense of four pages,” said Vincent Vargas, a mathematician at the French National Center for Scientific Research and co-author of the research with Rémi Rhodes of Aix-Marseille University, Antti Kupiainen of the University of Helsinki, François David of the French National Center for Scientific Research, and Colin Guillarmou of Paris-Saclay University.

The three papers forge a bridge between the pristine world of mathematics and the messy reality of physics — and they do so by breaking new ground in the mathematical field of probability theory. The work also touches on philosophical questions regarding the objects that take center stage in the leading theories of fundamental physics: quantum fields.

“This is a masterpiece in mathematical physics,” said Xin Sun, a mathematician at the University of Pennsylvania.

## Infinite Fields

In physics today, the main actors in the most successful theories are fields — objects that fill space, taking on different values from place to place.

In classical physics, for example, a single field tells you everything about how a force pushes objects around. Take Earth’s magnetic field: The twitches of a compass needle reveal the field’s influence (its strength and direction) at every point on the planet.

Fields are central to quantum physics, too. However, the situation here is more complicated due to the deep randomness of quantum theory. From the quantum perspective, Earth doesn’t generate one magnetic field, but rather an infinite number of different ones. Some look almost like the field we observe in classical physics, but others are wildly different.

But physicists still want to make predictions — predictions that ideally match, in this case, what a mountaineer reads on a compass. Assimilating the infinite forms of a quantum field into a single prediction is the formidable task of a “quantum field theory,” or QFT. This is a model of how one or more quantum fields, each with their infinite variations, act and interact.

Driven by immense experimental support, QFTs have become the basic language of particle physics. The Standard Model is one such QFT, depicting fundamental particles like electrons as fuzzy bumps that emerge from an infinitude of electron fields. It has passed every experimental test to date (although various groups may be on the verge of finding the first holes).

Physicists play with many different QFTs. Some, like the Standard Model, aspire to model real particles moving through the four dimensions of our universe (three spatial dimensions plus one dimension of time). Others describe exotic particles in strange universes, from two-dimensional flatlands to six-dimensional uber-worlds. Their connection to reality is remote, but physicists study them in the hopes of gaining insights they can carry back into our own world.

Polyakov’s Liouville field theory is one such example.

## Gravity’s Field

The Liouville field, which is based on an equation from complex analysis developed in the 1800s by the French mathematician Joseph Liouville, describes a completely random two-dimensional surface — that is, a surface, like Earth’s crust, but one in which the height of every point is chosen randomly. Such a planet would erupt with mountain ranges of infinitely tall peaks, each assigned by rolling a die with infinite faces.

Such an object might not seem like an informative model for physics, but randomness is not devoid of patterns. The bell curve, for example, tells you how likely you are to randomly pass a seven-foot basketball player on the street. Similarly, bulbous clouds and crinkly coastlines follow random patterns, but it’s nevertheless possible to discern consistent relationships between their large-scale and small-scale features.

Liouville theory can be used to identify patterns in the endless landscape of all possible random, jagged surfaces. Polyakov realized this chaotic topography was essential for modeling strings, which trace out surfaces as they move. The theory has also been applied to describe quantum gravity in a two-dimensional world. Einstein defined gravity as space-time’s curvature, but translating his description into the language of quantum field theory creates an infinite number of space-times — much as the Earth produces an infinite collection of magnetic fields. Liouville theory packages all those surfaces together into one object. It gives physicists the tools to measure the curvature —and hence, gravitation — at every location on a random 2D surface.

“Quantum gravity basically means random geometry, because quantum means random and gravity means geometry,” said Sun.

Polyakov’s first step in exploring the world of random surfaces was to write down an expression defining the odds of finding a particular spiky planet, much as the bell curve defines the odds of meeting someone of a particular height. But his formula did not lead to useful numerical predictions.

To solve a quantum field theory is to be able to use the field to predict observations. In practice, this means calculating a field’s “correlation functions,” which capture the field’s behavior by describing the extent to which a measurement of the field at one point relates, or correlates, to a measurement at another point. Calculating correlation functions in the photon field, for instance, can give you the textbook laws of quantum electromagnetism.

Polyakov was after something more abstract: the essence of random surfaces, similar to the statistical relationships that make a cloud a cloud or a coastline a coastline. He needed the correlations between the haphazard heights of the Liouville field. Over the decades he tried two different ways of calculating them. He started with a technique called the Feynman path integral and ended up developing a workaround known as the bootstrap. Both methods came up short in different ways, until the mathematicians behind the new work united them in a more precise formulation.

## Add ’Em Up

You might imagine that accounting for the infinitely many forms a quantum field can take is next to impossible. And you would be right. In the 1940s Richard Feynman, a quantum physics pioneer, developed one prescription for dealing with this bewildering situation, but the method proved severely limited.

Take, again, Earth’s magnetic field. Your goal is to use quantum field theory to predict what you’ll observe when you take a compass reading at a particular location. To do this, Feynman proposed summing all the field’s forms together. He argued that your reading will represent some average of all the field’s possible forms. The procedure for adding up these infinite field configurations with the proper weighting is known as the Feynman path integral.

It’s an elegant idea that yields concrete answers only for select quantum fields. No known mathematical procedure can meaningfully average an infinite number of objects covering an infinite expanse of space in general. The path integral is more of a physics philosophy than an exact mathematical recipe. Mathematicians question its very existence as a valid operation and are bothered by the way physicists rely on it.

“I’m disturbed as a mathematician by something which is not defined,” said Eveliina Peltola, a mathematician at the University of Bonn in Germany.

Physicists can harness Feynman’s path integral to calculate exact correlation functions for only the most boring of fields — free fields, which do not interact with other fields or even with themselves. Otherwise, they have to fudge it, pretending the fields are free and adding in mild interactions, or “perturbations.” This procedure, known as perturbation theory, gets them correlation functions for most of the fields in the Standard Model, because nature’s forces happen to be quite feeble.

But it didn’t work for Polyakov. Although he initially speculated that the Liouville field might be amenable to the standard hack of adding mild perturbations, he found that it interacted with itself too strongly. Compared to a free field, the Liouville field seemed mathematically inscrutable, and its correlation functions appeared unattainable.

## Up by the Bootstraps

Polyakov soon began looking for a workaround. In 1984, he teamed up with Alexander Belavin and Alexander Zamolodchikov to develop a technique called the bootstrap — a mathematical ladder that gradually leads to a field’s correlation functions.

To start climbing the ladder, you need a function which expresses the correlations between measurements at a mere three points in the field. This “three-point correlation function,” plus some additional information about the energies a particle of the field can take, forms the bottom rung of the bootstrap ladder.

From there you climb one point at a time: Use the three-point function to construct the four-point function, use the four-point function to construct the five-point function, and so on. But the procedure generates conflicting results if you start with the wrong three-point correlation function in the first rung.

Polyakov, Belavin and Zamolodchikov used the bootstrap to successfully solve a variety of simple QFT theories, but just as with the Feynman path integral, they couldn’t make it work for the Liouville field.

Then in the 1990s two pairs of physicists — Harald Dorn and Hans-Jörg Otto, and Zamolodchikov and his brother Alexei — managed to hit on the three-point correlation function that made it possible to scale the ladder, completely solving the Liouville field (and its simple description of quantum gravity). Their result, known by their initials as the DOZZ formula, let physicists make any prediction involving the Liouville field. But even the authors knew they had arrived at it partially by chance, not through sound mathematics.

“They were these kind of geniuses who guessed formulas,” said Vargas.

Educated guesses are useful in physics, but they don’t satisfy mathematicians, who afterward wanted to know where the DOZZ formula came from. The equation that solved the Liouville field should have come from some description of the field itself, even if no one had the faintest idea how to get it.

“It looked to me like science fiction,” said Kupiainen. “This is never going to be proven by anybody.”

## Taming Wild Surfaces

In the early 2010s, Vargas and Kupiainen joined forces with the probability theorist Rémi Rhodes and the physicist François David. Their goal was to tie up the mathematical loose ends of the Liouville field — to formalize the Feynman path integral that Polyakov had abandoned and, just maybe, demystify the DOZZ formula.

As they began, they realized that a French mathematician named Jean-Pierre Kahane had discovered, decades earlier, what would turn out to be the key to Polyakov’s master theory.

“In some sense it’s completely crazy that Liouville was not defined before us,” Vargas said. “All the ingredients were there.”

The insight led to three milestone papers in mathematical physics completed between 2014 and 2020.

They first polished off the path integral, which had failed Polyakov because the Liouville field interacts strongly with itself, making it incompatible with Feynman’s perturbative tools. So instead, the mathematicians used Kahane’s ideas to recast the wild Liouville field as a somewhat milder random object known as the Gaussian free field. The peaks in the Gaussian free field don’t fluctuate to the same random extremes as the peaks in the Liouville field, making it possible for the mathematicians to calculate averages and other statistical measures in sensible ways.

“Somehow it’s all just using the Gaussian free field,” Peltola said. “From that they can construct everything in the theory.”

In 2014, they unveiled their result: a new and improved version of the path integral Polyakov had written down in 1981, but fully defined in terms of the trusted Gaussian free field. It’s a rare instance in which Feynman’s path integral philosophy has found a solid mathematical execution.

“Path integrals can exist, do exist,” said Jörg Teschner, a physicist at the German Electron Synchrotron.

With a rigorously defined path integral in hand, the researchers then tried to see if they could use it to get answers from the Liouville field and to derive its correlation functions. The target was the mythical DOZZ formula — but the gulf between it and the path integral seemed vast.

“We’d write in our papers, just for propaganda reasons, that we want to understand the DOZZ formula,” said Kupiainen.

The team spent years prodding their probabilistic path integral, confirming that it truly had all the features needed to make the bootstrap work. As they did so, they built on earlier work by Teschner. Eventually, Vargas, Kupiainen and Rhodes succeeded with a paper posted in 2017 and another in October 2020, with Colin Guillarmou. They derived DOZZ and other correlation functions from the path integral and showed that these formulas perfectly matched the equations physicists had reached using the bootstrap.

“Now we’re done,” Vargas said. “Both objects are the same.”

The work explains the origins of the DOZZ formula and connects the bootstrap procedure —which mathematicians had considered sketchy — with verified mathematical objects. Altogether, it resolves the final mysteries of the Liouville field.

“It’s somehow the end of an era,” said Peltola. “But I hope it’s also the beginning of some new, interesting things.”

## New Hope for QFTs

Vargas and his collaborators now have a unicorn on their hands, a strongly interacting QFT perfectly described in a nonperturbative way by a brief mathematical formula that also makes numerical predictions.

Now the literal million-dollar question is: How far can these probabilistic methods go? Can they generate tidy formulas for all QFTs? Vargas is quick to dash such hopes, insisting that their tools are specific to the two-dimensional environment of Liouville theory. In higher dimensions, even free fields are too irregular, so he doubts the group’s methods will ever be able to handle the quantum behavior of gravitational fields in our universe.

But the fresh minting of Polyakov’s “master key” will open other doors. Its effects are already being felt in probability theory, where mathematicians can now wield previously dodgy physics formulas with impunity. Emboldened by the Liouville work, Sun and his collaborators have already imported equations from physics to solve two problems regarding random curves.

Physicists await tangible benefits too, further down the road. The rigorous construction of the Liouville field could inspire mathematicians to try their hand at proving features of other seemingly intractable QFTs — not just toy theories of gravity but descriptions of real particles and forces that bear directly on the deepest physical secrets of reality.

“[Mathematicians] will do things that we can’t even imagine,” said Davide Gaiotto, a theoretical physicist at the Perimeter Institute.