Strongly interacting matter exhibits new phases under extreme conditions. Matter was exposed to such extremes not only in the Early Universe, but also today in the cores of neutron stars, as well as in laboratory experiments at a much smaller scale. We study the underlying theory, Quantum Chromodynamics (QCD) with the methods of statistical physics and explore the various phases we may encounter in experiment, such as the Quark Gluon Plasma. We briefly summarize the experimental evidence for the new forms of matter and review the theoretical efforts to embed these findings in the broader context of quantum field theory, with special attention to exact and broken symmetries and critical behaviour.

We combine the recently introduced Non-Abelian Quasi-Particle Model (NAQPM) for gluons with an ideal Fermi gas of three quark species with the aim to describe the equation of state (energy density vs. temperature) of $2+1$ - flavour Lattice-QCD at zero chemical potential. Allowing temperature dependent masses for the fermions, we show that above a critical temperature $T_c$ the quark mass has to drop rapidly in order to obtain energy density values compatible with the Lattice-QCD results. Within this framework, thus, the restoration of chiral symmetry in the system is observed. Furthermore, we demonstrate that the gluon variance -- which is a fundamental quantity of the NAQPM -- is strongly correlated to the fermion mass and decreases by orders of magnitude through the transition. The high temperature phenomenological characteristics of the gluon appear consistent to properties of the perturbative QCD gluon. The model indicates that color deconfinement and chiral symmetry restoration are interrelated and classical configurations of the QCD dynamics play an important role to the criticality of the system.

We consider the case of a free real massive bulk scalar in D=4 dimensions, and embed two parallel plates as interfaces on which we impose non-linear boundary conditions, either Dirichlet- or Neumann-like, parameterized by a new coupling constant g. This mimics a non-Abelian gauge theory supplemented with boundary conditions on surfaces embedded in the bulk. We present the first evidence for a non-perturbative 1/g^2 boundary mass generation and its ensuing correction to the standard Casimir energy. This becomes possible by incorporating dynamical corrections to the effective boundary fields, which are used to build in the boundary conditions directly at the action level.

In this work, we present a first-principles lattice-QCD calculation of the unpolarized quark PDF for the pion and the kaon. The lattice data rely on matrix elements calculated for boosted mesons coupled to non-local operators containing a Wilson line. The calculations on this lattice ensemble correspond to two degenerate light, a strange, and a charm quark ($N_f=2+1+1$), using maximally twisted mass fermions with a clover term. The lattice volume is $32^3\times 64$, with a lattice spacing of 0.0934 fm, and a pion mass of 260 MeV. Matrix elements are calculated for hadron boosts of $|P_3| = 0,~0.41,~0.83,~1.25,~1.66,$ and 2.07 GeV. To match lattice QCD results to their light-cone counterparts, we employ two complementary frameworks: the large-momentum effective theory (LaMET) and the short-distance factorization (SDF). Using these approaches in parallel, we also test the lattice data to identify methodology-driven systematics. Results are presented for the standard quark PDFs, as well as the valence sector. Beyond obtaining the PDFs, we also explore the possibility of extracting information on SU(3) flavor-symmetry-breaking effects. For LaMET, we also parametrize the momentum dependence to obtain the infinite-momentum PDFs.

Finite density lattice QCD usually relies on extrapolations in baryon chemical potential ($\mu_B$), be it Taylor expansion, T' expansion (\cite{Borsanyi:2021sxv}) or analytical continuation. However, their range of validity is difficult to control. In the canonical formulation, the baryon density is the parameter of the system, not $\mu_B$. Here we demonstrate that we can access finite density QCD in the canonical formulation with physical quark masses. We present first results with both the strangeness ($n_S$) and baryon ($n_B$) densities as parameters. Specifically, we compute the QCD pressure and chemical potentials as functions of $n_B$ and $n_S$. Our computations rely on high-statistics simulations with 2+1 4HEX-staggered fermions.

Simulations of energy loss and hadronization are essential for understanding a range of phenomena in non-equilibrium strongly-interacting matter. We establish a framework for performing such simulations on a quantum computer and apply it to a heavy quark moving across a modest-sized 1+1D SU(2) lattice of light quarks. Conceptual advances with regard to simulations of non-Abelian versus Abelian theories are developed, allowing for the evolution of the energy in light quarks, of their local non-Abelian charge densities, and of their multi-partite entanglement to be computed. The non-trivial action of non-Abelian charge operators on arbitrary states suggests mapping the heavy quarks to qubits alongside the light quarks, and limits the heavy-quark motion to discrete steps among spatial lattice sites. Further, the color entanglement among the heavy quarks and light quarks is implemented using hadronic operators, and Domain Decomposition is shown to be effective in quantum state preparation. Scalable quantum circuits that account for the heterogeneity of non-Abelian charge sectors across the lattice are used to prepare the interacting ground-state wavefunction in the presence of heavy quarks. The discrete motion of heavy quarks between adjacent spatial sites is implemented using fermionic SWAP operations. Quantum simulations of the dynamics of a system on $L=3$ spatial sites are performed using IBM's ${\tt ibm\_pittsburgh}$ quantum computer using 18 qubits, for which the circuits for state preparation, motion, and one second-order Trotter step of time evolution have a two-qubit depth of 398. A suite of error mitigation techniques are used to extract the observables from the simulations, providing results that are in good agreement with classical simulations. The framework presented here generalizes straightforwardly to other non-Abelian groups, including SU(3) for quantum chromodynamics.

We perform a detailed analysis of the fermionic sign problem in a series of one dimensional integrals, that are achieved as extreme (one-site) limits of genuine physics models. Altogether we studied a Hubbard-like, a Gross-Neveu-like, a Thirring-like and a Chern-Simons-like integral. We compare the Lefschetz-thimble structure for these integrals with contours obtained with the holomorphic flow equations at different flow-times and with numerically optimized continuous integration contours, defined by a maximal value of the expectation values of the phases. With the holomorphic flow equation, we perform the large flow-time limit, so that the average phase corresponds to its value on the thimbles. In some of these integrals (the Hubbard-, Gross-Neveu-, and Chern-Simons-like integrals), we observe that the convergence to this value is not monotonic, meaning that there is an optimal flow-time where the sign problem is weaker than it is on the thimbles. Furthermore, we find that for all of these toy models, numerical optimization can find continuous contours on which the sign problem is considerably weaker than it is both on the thimbles and at flowed integration contours at the optimal flow-time.

We revisit quantum false vacuum decay for the one-dimensional Ising model, focusing on the real-time nucleation and growth of true vacuum bubbles. Via matrix product state simulations, we demonstrate that for a wide range of parameters, the full time-dependent quantum state is well described by a Gaussian ansatz in terms of domain wall operators, with the associated vacuum bubble wave function evolving according to the linearized time-dependent variational principle. The emerging picture shows three different stages of evolution: an initial nucleation of small bubbles, followed by semi-classical bubble growth, which in turn is halted by the lattice phenomenon of Bloch oscillations. Furthermore, we find that the resonant bubble only plays a significant role in a certain region of parameter-space. However, when significant, it does lead to an approximately constant decay rate during the intermediate stage. Moreover, this rate is in quantitative agreement with the analytical result of Rutkevich (Phys. Rev. B 60, 14525) for which we provide an independent derivation based on the Gaussian ansatz.

False vacuum decay in scalar quantum field theory (QFT) is a cornerstone of early Universe cosmology and high energy physics, yet its real-time dynamics is essentially inaccessible to classical computation due to its non-perturbative, highly entangled dynamics. We introduce a general Hamiltonian framework for simulating full interacting QFTs, using a spatial lattice of continuos-variable ``qumodes'' -- bosonic local oscillators whose high-dimensional local Hilbert space faithfully captures interacting field dynamics. This construction is rooted in continuous-variable quantum computing (CVQC), and provides a unified platform spanning efficient classical tensor-network methods and emerging photonic quantum hardware. The first key advance of this work is a robust and scaleable method for preparing the QFT in its correct initial vacuum state. We develop an imaginary-time preparation algorithm tailored to qumode lattices, that efficiently projects onto the vacuum even in strongly coupled regimes. This provides a controllable starting point for studying nonperturbative dynamics such as tunnelling and real-time decay. Building on this, we use a time-evolving block decimation algorithm to capture the real-time dynamics of the scalar field. Our second key advance is the identification and excitation of the negative fluctuation mode of the bounce configuration on the qumode lattice. A small displacement along this mode produces the expected tachyonic growth, driving fully coherent bubble nucleation without requiring classically supercritical seeds. This demonstrates that the qumode lattice captures non-perturbative quantum dynamics that lie beyond the classical treatments. Our results establish the qumode network as a scalable framework for non-equilibrium scalar QFT phenomena and pave the way for higher-dimensional studies and continuous-variable quantum computing implementations.

The Sachdev-Ye-Kitaev (SYK) model, fundamental to quantum chaos, many-body physics, and holographic duality, exhibits both maximal chaos and exact solvability, yet its entanglement entropy has never been measured experimentally. Here, we present the first measurement of entanglement entropy growth under the SYK Hamiltonian using IBM's superconducting quantum computer as a programmable quantum platform. We implement an optimized swap-based many-body interference protocol tailored for limited qubit connectivity to measure entanglement entropy. Additionally, we enhance the randomized measurement protocol for entropy measurement by expanding the unitary ensemble, which increases the circuit volume and pushes the capability of the near-term devices. However, we tackle this challenge by employing a quantum multi-programming strategy that parallelizes circuit execution. This work paves the way for exploring chaotic quantum dynamics and establishes scalable methods for investigating complex entanglement growth on current quantum platforms.

The transition from a false vacuum to the true vacuum is a real-time phenomenon of interest in many contexts. It represents a special challenge for strongly interacting non-Abelian gauge theories because standard spacetime lattices incorporate imaginary time from the outset. To attain real-time phenomena, Hamiltonian lattice methods are being developed for quantum computers. The present work considers SU(2) gauge theory on a minimal lattice in three spatial dimensions, where round-the-world strings called torelons can travel and interact. This minimal 3D lattice has a triamond structure whose properties are elucidated by calculating the spectrum of torelon states. Then, by introducing a twisted boundary condition, real-time evolution is used to demonstrate the decay of a false vacuum. Calculations in the present work are done on classical computers except for one benchmark study of imaginary time evolution that ran on the ibm_brisbane quantum computer.

We propose a method to construct the initial tensor representation of partition functions and observables for the tensor renormalization group (TRG). The TRG is a numerical calculation technique that utilizes a tensor network representations of physical quantities to investigate physical properties without encountering the sign problem. To apply the TRG, it is essential to construct a locally connected tensor network suitable for recursive coarse-graining. We present a systematic approach for translating a general tensor representation of the partition function to this form. Furthermore, we show the dependence of TRG algorithms on the choice of the initial tensor network representation and propose an improvement of TRG algorithms in this respect

A deep understanding of the intricate interactions between particles within a system is a key approach to revealing the essential characteristics of the system, whether it is an in-depth analysis of molecular properties in the field of chemistry or the design of new materials for specific performance requirements in materials science. To this end, we propose Graph Attention Hamiltonian Neural Network (GAHN), a neural network method that can understand the underlying structure of lattice Hamiltonian systems solely through the dynamic trajectories of particles. We can determine which particles in the system interact with each other, the proportion of interactions between different particles, and whether the potential energy of interactions between particles exhibits even symmetry or not. The obtained structure helps the neural network model to continue predicting the trajectory of the system and further understand the dynamic properties of the system. In addition to understanding the underlying structure of the system, it can be used for detecting lattice structural abnormalities, such as link defects, abnormal interactions, etc. These insights benefit system optimization, design, and detection of aging or damage. Moreover, this approach can integrate other components to deduce the link structure needed for specific parts, showcasing its scalability and potential. We tested it on a challenging molecular dynamics dataset, and the results proved its ability to accurately infer molecular bond connectivity, highlighting its scientific research potential.

We provide a universal framework for the quantum simulation of SU(N) Yang--Mills theories on fault-tolerant digital quantum computers adopting the orbifold lattice formulation. As warm-up examples, we also consider simple models, including scalar field theory and the Yang--Mills matrix model, to illustrate the universality of our formulation, which shows up in the fact that the truncated Hamiltonian can be expressed in the same simple form for any N, any dimension, and any lattice size, in stark contrast to the popular approach based on the Kogut--Susskind formulation. In all these cases, the truncated Hamiltonian can be programmed on a quantum computer using only standard tools well-established in the field of quantum computation. As a concrete application of this universal framework, we consider Hamiltonian time evolution by Suzuki--Trotter decomposition. This turns out to be a straightforward task due to the simplicity of the truncated Hamiltonian. We also provide a simple circuit structure that contains only CNOT and one-qubit gates, independent of the details of the theory investigated.

We compute the particle creation of a harmonic oscillator using the classical stochastic method. This recently constructed method reproduces all the vacuum expectation values in quantum theory. We prepare the vacuum state at the initial time and evolve it over time using Langevin equations of motion. By averaging over the ensemble, we compute the energy of the state at the final time and determine the amount of particles created. We verify that the particle creation agrees with predictions from quantum theory and thus give an evidence that this method can really compute the quantum dynamics.

We propose a method to construct a tensor network representation of partition functions without singular value decompositions nor series expansions. The approach is demonstrated for one- and two-dimensional Ising models and we study the dependence of the tensor renormalization group (TRG) on the form of the initial tensors and their symmetries. We further introduce variants of several tensor renormalization algorithms. Our benchmarks reveal a significant dependence of various TRG algorithms on the choice of initial tensors and their symmetries. However, we show that the boundary TRG technique can eliminate the initial tensor dependence for all TRG methods. The numerical results of TRG calculations can thus be made significantly more robust with only a few changes in the code. Furthermore, we study a three-dimensional $\mathbb{Z}_2$ gauge theory without gauge-fixing and confirm the applicability of the initial tensor construction. Our method can straightforwardly be applied to systems with longer range and multi-site interactions, such as the next-nearest neighbor Ising model.