Center for Theoretical Physics Publications

Hamiltonian simulation is a fundamental problem at the heart of quantum computation, and the associated simulation algorithms are useful building blocks for designing larger quantum algorithms. In order to be successfully concatenated into a larger quantum algorithm, a Hamiltonian simulation algorithm must succeed with arbitrarily high success probability \(1-\delta\) while only requiring a single copy of the initial state, a property which we call fully-coherent. Although optimal Hamiltonian simulation has been achieved by quantum signal processing (QSP), with query complexity linear in time \(t\) and logarithmic in inverse error \(\ln(1/\epsilon)\), the corresponding algorithm is not fully-coherent as it only succeeds with probability close to \(1/4\). While this simulation algorithm can be made fully-coherent by employing amplitude amplification at the expense of appending a \(\ln(1/\delta)\) multiplicative factor to the query complexity, here we develop a new fully-coherent Hamiltonian simulation algorithm that achieves a query complexity additive in \(\ln(1/\delta)\): \(\Theta\big( \|\mathcal{H}\| |t| + \ln(1/\epsilon) + \ln(1/\delta)\big)\). We accomplish this by compressing the spectrum of the Hamiltonian with an affine transformation, and applying to it a QSP polynomial that approximates the complex exponential only over the range of the compressed spectrum. We further numerically analyze the complexity of this algorithm and demonstrate its application to the simulation of the Heisenberg model in constant and time-dependent external magnetic fields. We believe that this efficient fully-coherent Hamiltonian simulation algorithm can serve as a useful subroutine in quantum algorithms where maintaining coherence is paramount.

Recent FPGA accelerator cards promise large acceleration factors for some specific computational tasks. In the context of Lattice QCD calculations, we investigate the possible gain of moving the \(SU(3)\) gauge field smearing routine to such accelerators. We study Xilinx Alveo U280 cards and use the associated Vitis high-level synthesis framework. We discuss the possible pros and cons of such a solution based on the gathered benchmarks.

We study 't Hooft anomalies and the related anomaly inflow for subsystem global symmetries. These symmetries and anomalies arise in a number of exotic systems, including models with fracton order such as the X-cube model. As is the case for ordinary global symmetries, anomalies for subsystem symmetries can be canceled by anomaly inflow from a bulk theory in one higher dimension; the corresponding bulk is therefore a non-trivial subsystem symmetry protected topological (SSPT) phase. We demonstrate these phenomena in several examples with continuous and discrete subsystem global symmetries, as well as time-reversal symmetry. For each example we describe the boundary anomaly, and present classical continuum actions for the corresponding bulk SSPT phases, which describe the response of background gauge fields associated with the subsystem symmetries. Interestingly, we show that the anomaly does not uniquely specify the bulk SSPT phase. In general, the latter may also depend on how the symmetry and the associated foliation structure on the boundary are extended into the bulk.

We study the phase structure of effective models of finite-density QCD using analytic and lattice simulation techniques developed for the study of non-Hermitian and \(\mathcal{PT}\)-symmetric QFTs. Finite-density QCD is symmetric under the combined operation of the charge and complex conjugation operators \(\mathcal{CK}\), which falls into the class of so-called generalized \(\mathcal{PT}\) symmetries. We show that \(\mathcal{PT}\)-symmetric quantum field theories can support patterned ground-state field configurations in the vicinity of a critical endpoint. We apply our methods to a lattice heavy quark model at nonzero chemical potential that displays patterning behavior for a range of parameters. We derive a simple approximate criterion for the formation of these patterns, which can be used with lattice results.

The nature of the Fermi gamma-ray Galactic Center Excess (GCE) has remained a persistent mystery for over a decade. Although the excess is broadly compatible with emission expected due to dark matter annihilation, an explanation in terms of a population of unresolved astrophysical point sources e.g., millisecond pulsars, remains viable. The effort to uncover the origin of the GCE is hampered in particular by an incomplete understanding of diffuse emission of Galactic origin. This can lead to spurious features that make it difficult to robustly differentiate smooth emission, as expected for a dark matter origin, from more "clumpy" emission expected for a population of relatively bright, unresolved point sources. We use recent advancements in the field of simulation-based inference, in particular density estimation techniques using normalizing flows, in order to characterize the contribution of modeled components, including unresolved point source populations, to the GCE. Compared to traditional techniques based on the statistical distribution of photon counts, our machine learning-based method is able to utilize more of the information contained in a given model of the Galactic Center emission, and in particular can perform posterior parameter estimation while accounting for pixel-to-pixel spatial correlations in the gamma-ray map. This makes the method demonstrably more resilient to certain forms of model misspecification. On application to Fermi data, the method generically attributes a smaller fraction of the GCE flux to unresolved point sources when compared to traditional approaches. We nevertheless infer such a contribution to make up a non-negligible fraction of the GCE across all analysis variations considered, with at least \(38^{+9}_{-19}\%\) of the excess attributed to unresolved points sources in our baseline analysis.

In holographic duality an eternal AdS black hole is described by two copies of the boundary CFT in the thermal field double state. This identification has many puzzles, including the boundary descriptions of the event horizons, the interiors of the black hole, and the singularities. Compounding these mysteries is the fact that, while there is no interaction between the CFTs, observers from them can fall into the black hole and interact. We address these issues in this paper. In particular, (i) we reformulate the meeting-behind-the-horizon puzzle by introducing a concept called causal connectability for any two quantum systems (which can in principle interact with each other), and show that this simple bulk experiment cannot be described by the boundary system at finite \(N\); (ii) we present a resolution of the puzzle by arguing there are emergent type III\(_1\) von Neumann subalgebras in the operator algebra of the boundary system in the large \(N\) limit; (iii) the type III\(_1\) structure leads to an explicit construction in the boundary theory of the evolution operator for a bulk in-falling observer, making manifest the emergence of the black hole horizons, the interiors, and the associated causal structure.

We review two algorithmic advances that bring us closer to reliable quantum simulations of model systems in high energy physics and beyond on noisy intermediate-scale quantum (NISQ) devices. The first method is the dimensional expressivity analysis of quantum circuits, which allows for constructing minimal but maximally expressive quantum circuits. The second method is an efficient mitigation of readout errors on quantum devices. Both methods can lead to significant improvements in quantum simulations, e.g., when variational quantum eigensolvers are used.

We study the mixed-state entanglement structure of chaotic quantum many-body systems at late times using the recently developed \(\textit{equilibrium approximation}\). A rich entanglement phase diagram emerges when we generalize this technique to evaluate the logarithmic negativity for various universality classes of macroscopically thermalized states. Unlike in the infinite temperature case, when we impose energy constraints at finite temperature, the phase diagrams for the logarithmic negativity and the mutual information become distinct. In particular, we identify a regime where the negativity is extensive but the mutual information is sub-extensive, indicating a large amount of \(\textit{bound entanglement}\). When applied to evaporating black holes, these results imply that there is quantum entanglement within the Hawking radiation long before the Page time, although this entanglement may not be distillable into EPR pairs.

Astrometry -- the precise measurement of positions and motions of celestial objects -- has emerged as a promising avenue for characterizing the dark matter population in our Galaxy. By leveraging recent advances in simulation-based inference and neural network architectures, we introduce a novel method to search for global dark matter-induced gravitational lensing signatures in astrometric datasets. Our method based on neural likelihood-ratio estimation shows significantly enhanced sensitivity to a cold dark matter population and more favorable scaling with measurement noise compared to existing approaches based on two-point correlation statistics, establishing machine learning as a powerful tool for characterizing dark matter using astrometric data.

The pion light-cone distribution amplitude (LCDA) is a central non-perturbative object of interest for high-energy exclusive processes in quantum chromodynamics. In this article, the second Mellin moment of the pion LCDA is determined as a proof-of-concept calculation for the first numerical implementation of the heavy-quark operator product expansion (HOPE) method. The resulting value for the second Mellin moment, determined in quenched QCD at a pion mass of \(m_\pi=550\) MeV at a factorization scale of 2 GeV is \( \langle \xi^2 \rangle = 0.210 \pm 0.013\text{ (stat.)} \pm 0.034\text{ (sys.)}\). This result is compatible with those from previous determinations of this quantity.

I look at propagation of jets that act as a probe of the medium in the phenomenologically relevant case of a short lived dilute Quark Gluon Plasma(QGP) created during heavy ion collisions. Working in the regime where the lifetime of the medium is of the order or smaller than the formation time of the energetic jet, I derive a factorization formula with a manifest separation of scales for an illustrative jet substructure observable using the Open Quantum system Effective Field Theory(EFT) approach developed in \cite{Vaidya:2020lih}, in terms of a medium structure function and a medium induced jet function . The medium structure function or the "PDF of the medium" that describes the observable independent physics of the QGP remains identical to that in a dilute long lived medium considered in \cite{Vaidya:2021vxu}, while the medium induced jet function is modified and incorporates quantum interference between the hard interaction that creates the jet and subsequent medium interactions. The medium jet function continues to enjoy a BFKL rapidity evolution just as for a long lived medium, albeit with a modified renormalization scale that now depends on the medium lifetime.

It is well known that in Anderson localized systems, starting from a random product state the entanglement entropy remains bounded at all times. However, we show that adding a single boundary term to an otherwise Anderson localized Hamiltonian leads to unbounded growth of entanglement. Our results imply that Anderson localization is not a local property. One cannot conclude that a subsystem has Anderson localized behavior without looking at the whole system, as a term that is arbitrarily far from the subsystem can affect the dynamics of the subsystem in such a way that the features of Anderson localization are lost.

We present exact results on a novel kind of emergent random matrix universality that quantum many-body systems at infinite temperature can exhibit. Specifically, we consider an ensemble of pure states supported on a small subsystem, generated from projective measurements of the remainder of the system in a local basis. We rigorously show that the ensemble, derived for a class of quantum chaotic systems undergoing quench dynamics, approaches a universal form completely independent of system details: it becomes uniformly distributed in Hilbert space. This goes beyond the standard paradigm of quantum thermalization, which dictates that the subsystem relaxes to an ensemble of quantum states that reproduces the expectation values of local observables in a thermal mixed state. Our results imply more generally that the distribution of quantum states themselves becomes indistinguishable from those of uniformly random ones, i.e. the ensemble forms a quantum state-design in the parlance of quantum information theory. Our work establishes bridges between quantum many-body physics, quantum information and random matrix theory, by showing that pseudo-random states can arise from isolated quantum dynamics, opening up new ways to design applications for quantum state tomography and benchmarking.

One of the most important properties of neural networks is the clustering of local minima of the loss function near the global minimum, enabling efficient training. Though generative models implemented on quantum computers are known to be more expressive than their traditional counterparts, it has empirically been observed that these models experience a transition in the quality of their local minima. Namely, below some critical number of parameters, all local minima are far from the global minimum in function value; above this critical parameter count, all local minima are good approximators of the global minimum. Furthermore, for a certain class of quantum generative models, this transition has empirically been observed to occur at parameter counts exponentially large in the problem size, meaning practical training of these models is out of reach. Here, we give the first proof of this transition in trainability, specializing to this latter class of quantum generative model. We use techniques inspired by those used to study the loss landscapes of classical neural networks. We also verify that our analytic results hold experimentally even at modest model sizes.

We prove a positive volume theorem for asymptotically AdS spacetimes: the maximal volume slice has nonnegative vacuum-subtracted volume, and the vacuum-subtracted volume vanishes if and only if the spacetime is identically pure AdS. Under the Complexity=Volume proposal, this constitutes a positive holographic complexity theorem. The result features a number of parallels with the positive energy theorem, including the assumption of an energy condition that excludes false vacuum decay (the AdS weak energy condition). Our proof is rigorously established in broad generality in four bulk dimensions, and we provide strong evidence in favor of a generalization to arbitrary dimensions. Our techniques also yield a holographic proof of Lloyd's bound for a class of bulk spacetimes. We further establish a partial rigidity result for wormholes: wormholes with a given throat size are more complex than AdS-Schwarzschild with the same throat size.

In this paper we use Euclidean gravity methods to show that charged black holes which are sufficiently close to extremality must be able to decay. The argument proceeds by showing that Euclidean gravity would otherwise imply a violation of charge quantization. As this is the assumption which leads to the weak gravity conjecture, our argument gives a derivation of that conjecture. We use a small negative cosmological constant as an infrared regulator, but our argument applies to near-extremal black holes which are arbitrarily small compared to the \(AdS\) curvature scale. We also give a universal formula for the density of black hole microstates which transform in each irreducible representation of any finite gauge group. Since each representation appears with nonzero fraction, this gives a new proof of the completeness hypothesis for finite gauge fields. Based on these observations we make two conjectures about many-body quantum physics: we propose a lower bound on the critical temperature for the instability of a semi-local quantum liquid, and we propose that our formula for the density of black hole microstates in each representation of a finite gauge group also applies at high energy to any quantum field theory with a finite group global symmetry.

We work out the effective gravitational action for 2D massive Euclidean fermions in a small mass expansion. Besides the leading Liouville action, the order \(m^2\) gravitational action contains a piece characteristic of the Mabuchi action, much as for 2D massive scalars, but also several non-local terms involving the Green's functions and Green's functions at coinciding points on the manifold.

These lectures, presented at the 2021 Les Houches Summer School on Dark Matter, provide an introduction to key methods and tools of indirect dark matter searches, as well as a status report on the field circa summer 2021. Topics covered include the possible effects of energy injection from dark matter on the early universe, methods to calculate both the expected energy distribution and spatial distribution of particles produced by dark matter interactions, an outline of theoretical models that predict diverse signals in indirect detection, and a discussion of current constraints and some claimed anomalies. These notes are intended as an introduction to indirect dark matter searches for graduate students, focusing primarily on intuition-building estimates and useful concepts and tools.

The low-energy spectrum and scattering of two-nucleon systems are studied with lattice quantum chromodynamics using a variational approach. A wide range of interpolating operators are used: dibaryon operators built from products of plane-wave nucleons, hexaquark operators built from six localized quarks, and quasi-local operators inspired by two-nucleon bound-state wavefunctions in low-energy effective theories. Sparsening techniques are used to compute the timeslice-to-all quark propagators required to form correlation-function matrices using products of these operators. Projection of these matrices onto irreducible representations of the cubic group, including spin-orbit coupling, is detailed. Variational methods are applied to constrain the low-energy spectra of two-nucleon systems in a single finite volume with quark masses corresponding to a pion mass of 806 MeV. Results for S- and D-wave phase shifts in the isospin singlet and triplet channels are obtained under the assumption that partial-wave mixing is negligible. Tests of interpolating-operator dependence are used to investigate the reliability of the energy spectra obtained and highlight both the strengths and weaknesses of variational methods. These studies and comparisons to previous studies using the same gauge-field ensemble demonstrate that interpolating-operator dependence can lead to significant effects on the two-nucleon energy spectra obtained using both variational and non-variational methods, including missing energy levels and other discrepancies. While this study is inconclusive regarding the presence of two-nucleon bound states at this quark mass, it provides robust upper bounds on two-nucleon energy levels that can be improved in future calculations using additional interpolating operators and is therefore a step toward reliable nuclear spectroscopy from the underlying Standard Model of particle physics.

Motivated by questions related to the landscape of flux compactifications, we combine new and existing techniques into a systematic, streamlined approach for computing vertical fluxes and chiral matter multiplicities in 4D F-theory models. A central feature of our approach is the conjecturally resolution-independent intersection pairing of the vertical part of the integer middle cohomology of smooth elliptic CY fourfolds, relevant for computing chiral indices and related aspects of 4D F-theory flux vacua. We illustrate our approach by analyzing vertical flux backgrounds for F-theory models with simple, simply-laced gauge groups and generic matter content, as well as models with U(1) gauge factors. We explicitly analyze resolutions of these F-theory models in which the elliptic fiber is realized as a cubic in \(\mathbb P^2\) over an arbitrary (e.g., not necessarily toric) smooth base, and confirm the resolution-independence of the intersection pairing of the vertical part of the middle cohomology. In each model we study, we find that vertical flux backgrounds can produce nonzero multiplicities for all anomaly-free chiral matter field combinations, suggesting that F-theory geometry imposes no additional linear constraints beyond those implied by anomaly cancellation.