It is for the first time that quantum simulation for high-energy physics (HEP) is studied in the U.S. decadal particle-physics community planning, and in fact until recently, this was not considered a mainstream topic in the community. This fact speaks of a remarkable rate of growth of this subfield over the past few years, stimulated by the impressive advancements in quantum information sciences (QIS) and associated technologies over the past decade, and the significant investment in this area by the government and private sectors in the U.S. and other countries. High-energy physicists have quickly identified problems of importance to our understanding of nature at the most fundamental level, from tiniest distances to cosmological extents, that are intractable with classical computers but may benefit from quantum advantage. They have initiated, and continue to carry out, a vigorous program in theory, algorithm, and hardware co-design for simulations of relevance to the HEP mission. This Roadmap is an attempt to bring this exciting and yet challenging area of research to the spotlight, and to elaborate on what the promises, requirements, challenges, and potential solutions are over the next decade and beyond.

Chaos and complexity entail an entropic and computational obstruction to describing a system, and thus are intrinsically difficult to characterize. In this paper, we consider time evolution by Gaussian Unitary Ensemble (GUE) Hamiltonians and analytically compute out-of-time-ordered correlation functions (OTOCs) and frame potentials to quantify scrambling, Haar-randomness, and circuit complexity. While our random matrix analysis gives a qualitatively correct prediction of the late-time behavior of chaotic systems, we find unphysical behavior at early times including an $$ \mathcal{O}(1) $$ O 1 scrambling time and the apparent breakdown of spatial and temporal locality. The salient feature of GUE Hamiltonians which gives us computational traction is the Haar-invariance of the ensemble, meaning that the ensemble-averaged dynamics look the same in any basis. Motivated by this property of the GUE, we introduce k-invariance as a precise definition of what it means for the dynamics of a quantum system to be described by random matrix theory. We envision that the dynamical onset of approximate k-invariance will be a useful tool for capturing the transition from early-time chaos, as seen by OTOCs, to late-time chaos, as seen by random matrix theory.

We develop new tools for isolating CFTs using the numerical bootstrap. A “cutting surface” algorithm for scanning OPE coefficients makes it possible to find islands in high-dimensional spaces. Together with recent progress in large-scale semidefinite programming, this enables bootstrap studies of much larger systems of correlation functions than was previously practical. We apply these methods to correlation functions of charge-0, 1, and 2 scalars in the 3d O(2) model, computing new precise values for scaling dimensions and OPE coefficients in this theory. Our new determinations of scaling dimensions are consistent with and improve upon existing Monte Carlo simulations, sharpening the existing decades-old 8σ discrepancy between theory and experiment.

We prove that higher-dimension operators contribute positively to the entropy of a thermodynamically stable black hole at fixed mass and charge. Our results apply whenever the dominant corrections originate at tree level from quantum field theoretic dynamics. More generally, positivity of the entropy shift is equivalent to a certain inequality relating the free energies of black holes. These entropy inequalities mandate new positivity bounds on the coefficients of higher-dimension operators. One of these conditions implies that the charge-to-mass ratio of an extremal black hole asymptotes to unity from above for increasing mass. Consequently, large extremal black holes are unstable to decay to smaller extremal black holes and the weak gravity conjecture is automatically satisfied. Our findings generalize to arbitrary spacetime dimension and to the case of multiple gauge fields. The assumptions of this proof are valid across a range of scenarios, including string theory constructions with a dilaton stabilized below the string scale.

A recent proposal relates two dimensional holographic conformal field theories deformed by the integrable $$ T\overline{T} $$ T T ¯ flow to AdS₃ with a finite radial cutoff. We investigate this proposal by studying perturbative correlation functions on the two sides. For low point correlators of the stress tensor, we successfully match the deformed CFT results at large central charge to bulk results obtained in classical pure gravity. The deformed CFT also provides definite predictions for loop corrections in the bulk. We then include matter fields in the bulk. To reproduce the classical bulk two-point function of a scalar operator we show that the deformed CFT needs to be augmented with double trace scalar operators, with the $$ T\overline{T} $$ T T ¯ operator yielding corrections corresponding to loops in the bulk.

We study the 6j symbol for the conformal group, and its appearance in three seemingly unrelated contexts: the SYK model, conformal representation theory, and perturbative amplitudes in AdS. The contribution of the planar Feynman diagrams to the three-point function of the bilinear singlets in SYK is shown to be a 6j symbol. We generalize the computation of these and other Feynman diagrams to d dimensions. The 6j symbol can be viewed as the crossing kernel for conformal partial waves, which may be computed using the Lorentzian inversion formula. We provide closed-form expressions for 6j symbols in d = 1, 2, 4. In AdS, we show that the 6j symbol is the Lorentzian inversion of a crossing-symmetric tree-level exchange amplitude, thus efficiently packaging the doubletrace OPE data. Finally, we consider one-loop diagrams in AdS with internal scalars and external spinning operators, and show that the triangle diagram is a 6j symbol, while one-loop n-gon diagrams are built out of 6j symbols.