Neutrinoless double-beta decay ($0\nu\beta\beta$) is a rare hypothesised process that, if discovered, would establish that the neutrino is Majorana, that is, it is its own antiparticle. Interpretation of experimental results relies on knowledge of nuclear matrix elements, whose large model uncertainty is the limiting factor in comparing measured (bounds on) half-lives to the neutrino mass. Nuclear effective field theory and lattice QCD have the potential to compute these matrix elements with better control over uncertainties, enhancing the discovery potential of next-generation $0\nu\beta\beta$ experiments. This work will survey various lattice QCD double-beta decay calculations and discuss their implications.

Three-hadron spectroscopy is a key frontier in our understanding of the hadron spectrum. In recent years, significant formal and numerical advances have paved the way for studying three-hadron processes directly from lattice QCD, with outstanding applications including the Roper resonance and the doubly charmed tetraquark. This requires theoretical frameworks that relate finite-volume energies to infinite-volume three-particle scattering amplitudes. In this contribution, I discuss recent progress in formulating such frameworks for generic three-hadron systems, and present numerical results for three-meson systems at maximal isospin with physical quark masses, as well as our recent investigation of the three-body dynamics of the doubly charmed tetraquark, $T_{\rm cc}$.

We will report recent progress on the QCD phase diagram at finite temperature and density. In particular, we discuss the universal scaling of the chiral transition in the limit of two massless quarks and one strange quark. We also discuss influence of other control parameter as chemical potentials, external magnetic field strength and number of quark flavors on the chiral transition. From calculations of Taylor expansion coefficients of the pressure w.r.t the baryon chemical potential and at imaginary chemical potential, we discuss estimates of the QCD critical point. Those estimates make use of the universal scaling ansatz of the Lee-Yang edge singularity.

Numerical studies of phase transitions in statistical and quantum lattice models provide crucial insights into the corresponding Conformal Field Theories (CFTs). In higher dimensions, comparing finite-volume numerical results to infinite-volume CFT data is facilitated by choosing the sphere $S^{d-1}$ as the spatial manifold. Recently, the fuzzy sphere regulator in Ref. [Zhu et al, Phys. Rev. X 13 021009 (2023)] has enabled such studies with exact rotational invariance, yielding impressive agreement with known 3D Ising CFT predictions, as well as new results. However, systematic improvements and a deeper understanding of finite-size corrections remain essential. In this work, we revisit the fuzzy sphere regulator, focusing on the original Ising model, with two main goals. First, we assess the robustness of this approach using Conformal Perturbation Theory (CPT), to which we provide a detailed guidebook. We demonstrate how CPT provides a unified framework for determining the critical point, the speed of light, and residual deviations from CFT predictions. Applying this framework, we study finite-size corrections and clarify the role of tuning the model in minimizing these effects. Second, we develop a novel method for extracting Operator Product Expansion (OPE) coefficients from fuzzy sphere data. This method leverages the sensitivity of energy levels to detuning from criticality, providing new insights into level mixing and avoided crossings in finite systems. Our work also includes validation of CPT in a 1+1D Ising model away from the integrable limit.