We present a first real-time study of hadronic scattering in a $(1+1)$-dimensional SU(2) lattice gauge theory with fundamental fermions using tensor-network techniques. Working in the gaugeless Hamiltonian formulation, we investigate scattering processes across sectors of fixed global baryon number $B = 0, 1, 2$, corresponding respectively to meson--meson, meson--baryon, and baryon--baryon collisions. At strong coupling, the $B = 0$ and $B = 2$ channels exhibit predominantly elastic dynamics closely resembling the U(1) Schwinger model. The mixed $B = 1$ sector displays qualitatively new behavior: meson and baryon wavepackets become entangled during the collision, with the slower state becoming spatially delocalized while the faster one propagates ballistically. We characterize these processes through local observables, entanglement entropy, and the information lattice.

We present the status of calculations of the form factors of the most relevant heavy-to-heavy and heavy-to-light decay channels. Using seven $N_f = 2+1+1$ HISQ ensembles, with lattice spacings ranging from 0.15 fm down to 0.06 fm, we calculate the form factors of the decays, including correlations among them. More than half of our ensembles feature physical pion masses, and the heavy quarks are simulated at their physical masses using the Wilson-clover action with the Fermilab interpretation. Even though we have recently seen huge qualitative and quantitative leaps in the characterization of heavy-to-heavy decays, these advances have failed to translate into improvements for the inclusive vs exclusive question, or the matter of the Lepton Flavor Universality ratios. In particular, in the $B\to D^{\ast}\ell\nu$ channel, the current situation of the lattice-QCD form factors is far from clear. Further, the latest lattice-QCD results on the heavy-to-light form factors display unexplained tensions that must urgently be resolved. The work presented here is an attempt to address these issues.

Lattice QCD (LQCD) calculations play a key role in the establishment of flavor anomalies. One of the most recent advancements in LQCD to this end has been the publication of several calculations of the $B\to D^\ast\ell\nu$ form factors, but despite all the anticipation, the LQCD results have been unable to give a final answer to the questions it was destined to answer. In this work I briefly review what is the current status of heavy-to-heavy and heavy-to-light semileptonic decays calculations in LQCD, and what we can expect for the near and not-so-near future.

The extraction of the static quark-antiquark potential from lattice QCD suffers from the poor signal-to-noise ratio of Wilson loops at large Euclidean times. To overcome this, smearing methods or the Coulomb gauge are used to improve the ground-state overlap with respect to the straight Wilson line trial state within the Wilson loop. To find excited states, complicated shapes are introduced to generate specific quantum numbers. Here, we introduce a neural-network parametrization of trial states, constructed with gauge-equivariant layers and optimized with a loss function that favors ground and excited states. In the quenched theory, we automatically obtain the interpolators for the ground and excited states.

We report our new lattice results for the sphaleron rate calculated within a thermal effective field theory of soft SU(N) gluons, where $N=2,3$, for a wide range of temperatures spanning from $0.6$-$10^{15}$ GeV at sufficiently large volumes. Comparing these results with sphaleron rates in a non-thermal SU(N) plasma where the infrared gluons are over-occupied, we estimate the typical thermalization time for these ultra-soft gluons during the early stages of reheating after inflation. We have also calculated the non-perturbative thermal axion production rate using lattice techniques which shows significant deviation from its perturbative estimate even at the electroweak scale.

We consider the problem of identifying the onset of the conformal window for QCD with $N_f$ massless flavors in the fundamental representation, and propose a new effective method to determine it from lattice simulations. This method is based on the investigation of the so-called Roberge-Weiss transition temperature $T_{RW}$, which is encountered at specific values of the imaginary baryon chemical potential, and can also be interpreted as the inverse of the critical spatial size at which charge conjugation is spontaneously broken in a finite box. Since $T_{RW}$ corresponds to a genuine phase transition for any value of the quark masses, it is a well-defined quantity; we argue that the critical $N_f$ at which $T_{RW}$ vanishes in the chiral limit coincides with the onset of the conformal window. We implement our proposal by investigating QCD with $N_f = 8$ flavors, discretized via stout improved staggered fermions and the tree-level improved Symanzik pure gauge action, at Euclidean temporal extents $N_t = 8, 10, 12, 16, 24$. In this case, we find evidence that $T_{RW}$ already vanishes in the chiral limit, indicating that $N_f = 8$ is already in the conformal window.

We present a solvable Hamiltonian that realizes an exact lattice chiral $U(1)_V \times U(1)_A$ symmetry. Nielsen-Ninomiya-type no-go theorems are evaded by using lattice bosons rather than fermions. The continuum limit is a compact boson field theory with an axion-like coupling. The $U(1)_V$ symmetry shifts the scalar, while $U(1)_A$ acts on local operators associated with short axion strings and is transmuted into a higher-form symmetry in the continuum limit. We demonstrate the chiral anomaly by showing that the lattice theta angle is shifted by an axial rotation when $U(1)_V$ is gauged. Gauging either $U(1)_V$ or $U(1)_A$ leads to lattice non-invertible and 2-group symmetries, respectively, matching the continuum picture.

Charmonia are often supposed to provide simple hydrogen-like ``atomic'' systems that can be used to obtain insights into heavier-quark QCD. We use continuum Schwinger function methods to analyse this hypothesis in connection with ground-state scalar and pseudoscalar charmonia and find that a more complex picture of these states may be necessary. For instance, considering orbital angular momentum, the $\chi_{c0}$ is not a simple $P$-wave system; similarly, the $\eta_c$ wave function contains more than merely $S$-wave contributions. The distribution amplitudes (DAs) and distribution functions (DFs) of these mesons are also nontrivial. For instance, the $\chi_{c0}$ DA is not positive definite: owing to QCD symmetries, it possesses domains of balanced negative and positive support. This feature is also expressed in the $\chi_{c0}$ DF, but differences between $\chi_{c0}$ and $\eta_c$ DFs diminish under scale evolution. Notably, the light-front momentum fraction carried by glue is the same in both states: it is 10\% less than the in-pion glue momentum fraction. Whilst experimental confirmation of the predictions herein is unlikely, our results should serve as benchmarks for complementary theory attempts to understand local and global structural features of heavier-quark hadrons.

We compute the spectrum of light mesons, composed by up, down, and strange quarks, using a symmetry-preserving approximation that permits the inclusion of fully-dressed quark-gluon vertices in the key dynamical equations. This method is characterized by the use of the standard symmetric kinematic configuration as a seed in the corresponding Schwinger-Dyson equation, yielding finally the full kinematic dependence of all eight form factors composing the transversely-projected quark-gluon vertex. The extension of this approach to the case of distinct nonvanishing current quark masses is discussed, and the compatibility with the fundamental Ward-Takahashi identities demonstrated. The corresponding Bethe-Salpeter kernel is composed by three different diagrammatic structures, which may be deduced from the attendant quark gap equation by applying the standard "cutting" rules. The masses of the light mesons are computed by first determining the eigenvalue of the Bethe-Salpeter equation as a function of Euclidean momenta, and then using the Schlessinger extrapolation method to determine the Minkowski momentum for which this eigenvalue becomes unity. The resulting meson masses are in good agreement with experimental values, and substantially improve upon predictions from the rainbow-ladder approximation.

In this work, we study exotic theta terms in the 2+1d $\phi$-theory, which provides a continuum description of the XY-plaquette model. The $\phi$-theory can be viewed as a fractonic analogue of the 1+1d compact boson and exhibits momentum and winding subsystem symmetries. In this theory, discontinuous field configurations play a crucial role. Although such configurations spoil the naive topology of the field, they induce nontrivial backreactions that give rise to new topological terms. We study two types of theta terms, which we call the bulk theta term and the foliated theta term. The foliated theta term is constructed by coupling winding currents on neighboring leaves of a foliation. Remarkably, the corresponding theta angle can vary spatially without affecting the classical equations of motion. Both theta terms lead to generalized Witten effects, in which vortex operators carrying winding subsystem charge acquire momentum subsystem charge. In the case of the foliated theta angle, the Witten effect exhibits a more intricate structure: vortex operators acquire a quadrupolar momentum charge. We demonstrate these features using lattice realizations based on the modified Villain formulation.

We study the quantum Ising model on (2+1)-dimensional anti-de Sitter space using Matrix Product States (MPS) and Matrix Product Operators (MPOs). We explore the bulk phase diagram of the theory on regular tessellations of hyperbolic space with coordination number seven and find disordered and ordered phases separated by a phase transition. We find that the boundary-boundary spin correlation function exhibits power law scaling deep in the disordered phase of the Ising model consistent with holography. At the critical point, we find the boundary entanglement entropy scales logarithmically with subsystem size but away from this, we see a linear scaling. In comparison, the full system exhibits a volume law scaling, which is expected in chaotic and/or highly connected systems. We also measure Out of time Ordered Correlators (OTOCs) to explore the scrambling behavior of the theory.

The pole structure of the $\Lambda(1405)$ has been a topic of debate for a long time. Chiral perturbation theory predicts that its experimental spectrum may be explained by a two pole structure originating in the $SU(3)$ chiral dynamics of the baryon-meson interaction. The $SU(3)$-symmetric flavor point is readily accessible in lattice QCD, in this work we study the baryon-meson states directly at this point. We construct interpolation operators that belong to the irreducible representations of $SU(3)$ that are attractive in the channel with the quantum numbers of the (singlet and two octets). The extracted energy levels can be used as input for chiral perturbation theory to find the poles associated with each representation. The relevant correlation functions are computed on $SU(3)$-symmetric ensembles with $M_{\pi}\approx 714$ MeV using the distillation technique.