We present recent results from the FASTSUM collaboration, using anisotropic lattice QCD to study spectral properties of heavy quarkonia and open heavy flavour systems at high temperature. For heavy quarkonium, our results using a number of different methods suggest a small but significant and robust negative mass shift as well as an increasing thermal width. We present the first lattice results for masses and spectral functions of B mesons at high temperature, and preliminary results for a high-precision calculation of the static quark potential.

The Lorentzian type IIB matrix model is a promising candidate for a nonperturbative formulation of superstring theory. In this model, the eigenvalue distribution of the $N\times N$ bosonic matrices $A_\mu$ $(\mu = 0 , \ldots , 9)$ represents an emergent spacetime, which is determined by the dynamics of the model in the large-$N$ limit. Here we perform numerical simulations of the model overcoming the sign problem by the complex Langevin method with the matrix size $N$ up to $128$. In order to avoid the singular drift problem due to the Pfaffian, which appears after integrating out the fermionic matrices, we deform the model in a manner inspired by the supersymmetric deformation, which is used to define the ``polarized type IIB matrix model'' in the Euclidean case. We find that the deformed model exhibits a phase in which (3+1)-dimensional expanding spacetime emerges with both space and time being smooth and real.

Discovery of instantons in colliders will provide experimental evidence for the topological properties of the QCD vacuum. In this work, we propose jet correlation observables that can unambiguously discriminate between instanton-induced processes and perturbative hard scattering events in pp collisions at LHC energies. By calculating the instanton sizes and their separations in 2+1 flavor QCD with physical quark masses, we provide constraints on the center-of-mass energies of the produced hadrons in an instanton-induced process. Our proposal is directly applicable for future ep measurements at the Electron-Ion Collider, offering a cleaner environment to probe instanton-induced processes.

Local gauge symmetry underlies fundamental interactions and strongly correlated quantum matter, yet existing machine-learning approaches lack a general, principled framework for learning under site-dependent symmetries, particularly for intrinsically nonlocal observables. Here we introduce a gauge-equivariant graph neural network that embeds non-Abelian symmetry directly into message passing via matrix-valued, gauge-covariant features and symmetry-compatible updates, extending equivariant learning from global to fully local symmetries. In this formulation, message passing implements gauge-covariant transport across the lattice, allowing nonlocal correlations and loop-like structures to emerge naturally from local operations. We validate the approach across pure gauge, gauge-matter, and dynamical regimes, establishing gauge-equivariant message passing as a general paradigm for learning in systems governed by local symmetry.

The nature of the finite temperature phase transition of QCD depends on the particle density and the mass of the dynamical quarks. We discuss the properties of the phase transition at high density, considering an effective theory describing the high-density heavy-quark limit of QCD. This effective theory is a simple model in which the Polyakov loop is a dynamical variable, and the quark Boltzmann factor is controlled by only one parameter, $C(\mu,m_q)$, which is a function of the quark mass $m_q$ and the chemical potential $\mu$. The Polyakov loop is an order parameter of $Z_3$ symmetry, and the fundamental properties of the phase transition are thought to be determined by the $Z_3$ symmetry broken by the phase transition. By replacing the Polyakov loop with $Z_3$ spin, we find that the effective model becomes a three-dimensional three-state Potts model ($Z_3$ spin model) with a complex external field term. We investigate the phase structure of the Potts model and discuss QCD in the heavy-quark region. As the density varies from $\mu=0$ to $\mu=\infty$, we find that the phase transition is first order in the low-density region, changes to a crossover at the critical point, and then becomes first-order again. This strongly suggests the existence of a first-order phase transition in the high density heavy-quark region of QCD.

We present a determination of the scalar and tensor $\Lambda\to p$ transition form factors using lattice QCD. These form factors are relevant for semileptonic hyperon decays in the presence of extensions of the Standard Model that include scalar and tensor interactions. The calculation is carried out using a gauge ensemble of twisted mass fermions at the physical pion mass, following the same strategy as our recent study on vector and axial form factors for the same transition. We provide the complete set of form factors as functions of $q^2$ employing a model-independent parametrization. We examine their impact on searches for non-standard charged-current interactions via the muon-to-electron decay-rate ratio $R^{\mu e}=\Gamma(\Lambda\to p\mu\bar\nu_\mu)/\Gamma(\Lambda\to pe\bar\nu_e)$, where scalar and tensor contributions enter linearly and are helicity-enhanced relative to the electron channel. We compare this first-principles prediction for the decay-rate ratio with recent experimental measurements, thereby enabling improved constraints on non-standard charged-current interactions.