Nonabelian gauge theories are chaotic in the classical limit. We discuss new evidence from SU(2) lattice gauge theory that they are also chaotic at the quantum level. We also describe possible future studies aimed at discovering the consequences of this insight. Based on a lecture presented by the first author at the Particles and Plasmas Symposium 2024.

A framework for computing the superfluid gap in an effective field theory (EFT) of fermions interacting via momentum independent contact forces is developed. The leading universal corrections in the EFT are one-loop in-medium effects at the Fermi surface, and reproduce the well-known Gor’kov-Melik-Barkhudarov result. The complete subleading universal corrections are presented here, and include one-loop effects away from the Fermi surface, two-loop in-medium effects, as well as modifications to the fermion propagator. Together, these effects are found to reduce the gap at low densities. Applications to neutron superfluidity in neutron stars are discussed.

We explore the utility of qutrits and qubits for simulating the flavor dynamics of dense neutrino systems. The evolution of such systems impacts some important astrophysical processes, such as core-collapse supernovae and the nucleosynthesis of heavy nuclei. Many-body simulations require classical resources beyond current computing capabilities for physically relevant system sizes. Quantum computers are therefore a promising candidate to efficiently simulate the many-body dynamics of collective neutrino oscillations. Previous quantum simulation efforts have primarily focused on properties of the two-flavor approximation due to their direct mapping to qubits. Here, we present new quantum circuits for simulating three-flavor neutrino systems on qutrit- and qubit-based platforms, and demonstrate their feasibility by simulating systems of two, four and eight neutrinos on IBM and Quantinuum quantum computers.

Rapidly growing capabilities of quantum simulators to probe quantum many-body phenomena require new methods to characterize increasingly complex states. We present a protocol that constrains quantum states by experimentally measured correlation functions which only scales polynomially with system size. This method enables measurement of a quantum state’s entanglement structure, opening a new route to study entanglement-related phenomena. Our approach extends Gaussian state parameterizations by systematically incorporating higher-order correlations. We show the protocol’s usefulness in conjunction with current and forthcoming experimental capabilities, focusing on weakly interacting fermions as a proof of concept. Here, the lowest non-trivial expansion quantitatively predicts early time thermalization dynamics, including signaling the on-set of quantum chaos indicated by the entanglement Hamiltonian.

This work is supported by the DOE, Office of Science, Office of Nuclear Physics, IQuS (\url{https://iqus.uw.edu}), via the program on Quantum Horizons: QIS Research and Innovation for Nuclear Science under Award DE-SC0020970 and by the European Union’s Horizon Europe research and innovation program under Grant Agreement No. 101113690 (PASQuanS2.1), the ERC Starting grant QARA (Grant No.~101041435), the EU-QUANTERA project TNiSQ (N-6001), and by the Austrian Science Fund (FWF): COE 1 and quantA. This work was enabled, in part, by the use of advanced computational, storage and networking infrastructure provided by the Hyak supercomputer system at the University of Washington

Simulations of scattering processes are essential in understanding the physics of our universe. Computing relevant scattering quantities from ab initio methods is extremely difficult on classical devices because of the substantial computational resources needed. This work reports the development of an algorithm that makes it possible to obtain phase shifts for generic non-relativistic elastic scattering processes on a quantum computer. Such algorithm is based on extracting phase shifts from the direct implementation of the real-time evolution. The algorithm is improved by a variational procedure, making it more accurate and resistant to the noise of quantum . The reliability of the algorithm is first demonstrated by means of classical numerical simulations for different potentials, and later tested on existing quantum hardware, specifically on IBM quantum processors.

We study the spin-dependent in-medium dynamics of quarkonia by using the potential nonrelativistic QCD (pNRQCD) and the open quantum system framework. We consider the pNRQCD Lagrangian valid up to the order. r/M^0=r and r^0/M=1/M in the double power counting. By considering the Markovian condition and applying the Wigner transformation upon the diagonal spin components of the quarkonium density matrix with the semiclassical expansion, we systematically derive the Boltzmann transport equation for quarkonia with polarization dependence in the quantum optical limit. Unlike the spin-independent collision terms governed by certain chromoelectric field correlators, new gauge invariant correlators of chromomagnetic fields determine the recombination and dissociation terms with polarization dependence at the order we are working. We also derive a Lindblad equation describing the in-medium transitions between spin-singlet and spin-triplet heavy quark-antiquark pairs in the quantum Brownian motion limit. The Lindblad equation is governed by new transport coefficients defined in terms of the chromomagnetic field correlators. Our formalism is generic and valid for both weakly-coupled and strongly-coupled quark gluon plasmas. It can be further applied to study spin alignment of vector quarkonia in heavy ion collisions.