Masanori Hanada

Understanding the phases of quantum chromodynamics (QCD) at finite temperature and baryon density is crucial for describing matter in heavy-ion collisions and the interior of neutron stars. While the transition from the confined phase to the deconfined phase has been extensively studied at vanishing chemical potential, recent theoretical work suggests the existence of intermediate phases with partial deconfinement, where only a subset of the color degrees of freedom is deconfined. In this paper, we study the generalization of partial deconfinement in QCD in the Veneziano large-Nc limit (Nf/Nc fixed) to finite baryon chemical potential. We find that the partially deconfined phase has finer structure, and hence that there are four phases in QCD at finite temperature and moderately large chemical potential: complete confinement, complete deconfinement, and two kinds of partial deconfinement. A key ingredient is the refined understanding of the meaning of the ‘Gross-Witten-Wadia point’, i.e., the opening of a gap in the distribution of Polyakov line phases: either string condensation or baryon condensation causes the GWW transition. As a by-product, we observe the emergence of the QCD critical point from the interplay between baryon condensation and partial deconfinement. While our approach is necessarily qualitative due to the fermion sign problem, it provides a unified theoretical framework for understanding the rich phase structure of dense QCD matter and offers new perspectives on the location and nature of the QCD critical point.

Quantum computers promise to simulate complex quantum field theories that are intractable for classical computers, potentially revealing new physics in strongly interacting systems. Current approaches for simulating Yang-Mills gauge theories face significant technical barriers due to complex group structures and complicated truncation schemes that vary drastically between different theories. Here we show that the orbifold lattice formulation provides a universal framework for quantum simulation of Yang-Mills theories with arbitrary gauge groups and dimensions. Our approach reduces all theories to the same simple Hamiltonian form, enabling implementation with standard quantum gates regardless of system complexity. We demonstrate explicit quantum circuits using only controlled-NOT and single-qubit operations, with concrete resource estimates for time evolution algorithms. This universal framework simplifies quantum simulation of gauge theories and enables systematic scaling to larger systems on fault-tolerant quantum computers.

The Sachdev-Ye-Kitaev (SYK) model may provide us with a good starting point for the experimental study of quantum chaos and holography in the laboratory. Still, the four-local interaction of fermions makes quantum simulation challenging, and it would be good to search for simpler models that keep the essence. In this paper, we argue that the four-local interaction may not be important by introducing a few models that have two-local interactions. The first model is a generalization of the spin-SYK model, which is obtained by replacing the spin variables with SU(d) generators. Simulations of this class of models might be straightforward on qudit-based quantum devices. We study the case of d=3,4,5,6 numerically and observe quantum chaos already for two-local interactions in a wide energy range. We also introduce modifications of spin-SYK and SYK models that have similar structures as the SU(d) model (e.g., H=∑_{p,q}J_{pq}χ_{p}χ_{p+1}χ_{q}χ_{q+1} instead of the original SYK Hamiltonian H=∑_{p,q,r,s}J_{pqrs}χ_{p}χ_{q}χ_{r}χ_{s}), which shows strongly chaotic features although the interaction is essentially two-local. These models may be a good starting point for the quantum simulation of the original SYK model.

We apply the variational Monte Carlo method based on neural network quantum states, using a neural autoregressive flow architecture as our ansatz, to determine the ground state wave function of the bosonic SU ( N ) Yang-Mills-type two-matrix model at strong coupling. Previous literature hinted at the inaccuracy of such an approach at strong coupling. In this work, the accuracy of the results is tested using lattice Monte Carlo simulations: we benchmark the expectation value of the energy of the ground state for system sizes N that are beyond brute-force exact diagonalization methods. We observe that the variational method with neural network states reproduces the right ground state energy when the width of the network employed in this work is sufficiently large. We confirm that the correct result is obtained for N = 2 and 3, while obtaining a precise value for N = 4 requires more resources than the amount available for this work.

A bstract: For matrix model and QFT, we discuss how dual gravitational geometry emerges from matrix degrees of freedom (specifically, adjoint scalars in super Yang-Mills theory) and how operator algebra that describes an arbitrary region of the bulk geometry can be constructed. We pay attention to the subtle difference between the notions of wave packets that describe low-energy excitations: QFT wave packet associated with the spatial dimensions of QFT, matrix wave packet associated with the emergent dimensions from matrix degrees of freedom, and bulk wave packet which is a combination of QFT and matrix wave packets. In QFT, there is an intriguing interplay between QFT wave packet and matrix wave packet that connects quantum entanglement and emergent geometry. We propose that the bulk wave packet is the physical object in QFT that describes the emergent geometry from entanglement. This proposal sets a unified view on two seemingly different mechanisms of holographic emergent geometry: one based on matrix eigenvalues and the other based on quantum entanglement. Further intuition comes from the similarity to a traversable wormhole discussed as the dual description of the coupled SYK model by Maldacena and Qi: the bulk can be seen as an eternal traversable wormhole connecting boundary regions.

Linear confinement with Casimir scaling of the string tension in confining gauge theories is a consequence of a certain property of the Polyakov loop related to random matrices. This mechanism does not depend on the details of the theories (neither the gauge group nor dimensions) and explains approximate Casimir scaling below string-breaking length. In this paper, we study 3d SU(2) pure Yang-Mills theory numerically and find the same random-matrix behavior for rectangular Wilson loops. We conjecture that this is a universal feature of strongly coupled confining gauge theories.

We construct a model of Pauli spin operators with all-to-all 4-local interactions by replacing Majorana fermions in the SYK model with spin operators. Equivalently, we replace fermions with hard-core bosons. We study this model numerically and compare the properties with those of the SYK model. We observe a striking quantitative coincidence between the spin model and the SYK model, which suggests that this spin model is strongly chaotic and, perhaps, can play some role in holography. We also discuss the path-integral approach with multi-local fields and the possibility of quantum simulations. This model may be an interesting target for quantum simulations because Pauli spins are easier to implement than fermions on qubit-based quantum devices.

We propose an orbifold lattice formulation of QCD suitable for quantum simulations. We show explicitly how to encode gauge degrees of freedom into qubits using noncompact variables, and how to write down a simple truncated Hamiltonian in the coordinate basis. We show that SU(3) gauge group variables and quarks in the fundamental representation can be implemented straightforwardly on qubits, for arbitrary truncation of the gauge manifold.

We describe how the general mechanism of partial deconfinement applies to large-N QCD and the partially-deconfined phase inevitably appears between completely-confined and completely-deconfined phases. Furthermore, we propose how the partial deconfinement can be observed in the real-world QCD with SU(3) gauge group. We propose the relationship between the behaviors of the Polyakov loop and other quantities. We test our proposal against lattice simulation data and find a nontrivial matching.

To simulate bosons on a qubit- or qudit-based quantum computer, one has to regularize the theory by truncating infinite-dimensional local Hilbert spaces to finite dimensions. In the search for practical quantum applications, it is important to know how big the truncation errors can be. In general, it is not easy to estimate errors unless we have a good quantum computer. In this paper, we show that traditional sampling methods on classical devices, specifically Markov Chain Monte Carlo, can address this issue for a rather generic class of bosonic systems with a reasonable amount of computational resources available today. As a demonstration, we apply this idea to the scalar field theory on a two-dimensional lattice, with a size that goes beyond what is achievable using exact diagonalization methods. This method can be used to estimate the resources needed for realistic quantum simulations of bosonic theories, and also, to check the validity of the results of the corresponding quantum simulations.