A crucial step in extracting physical predictions from lattice QCD simulations is the scale setting, i.e. the determination of the lattice spacing ($a$) in physical units. Herein, the relative scale setting for different $\beta$'s is discussed, using the Landau gauge gluon propagator computed with large statistical ensembles. After setting the relative scales, finite size effects are observed in the ultraviolet regime and handled in an effective description, inspired in perturbation theory. The new devised procedure is efficient in handling the finite size effects, linking the lattice simulations with continuum perturbation theory for the high momenta regime. Furthermore, the procedure can be extended to handle other Green functions computed within lattice QCD simulations.

Although ensemble generation remains a central challenge in lattice field theory simulations, recent advances in generative modeling may offer a path to accelerated sampling in these contexts. In this work, we implement a framework for efficiently training diffusion models acting on ${\rm SU}(N)$ degrees of freedom, adapting the traditional score matching technique to the group manifold. We demonstrate that our models can effectively reproduce several target densities, resulting in precise unbiased expectation values. These results mark a step for diffusion models towards modeling full ${\rm SU}(N)$ lattice field theories, including lattice Quantum Chromodynamics.

We introduce a learning method for recovering action parameters in lattice field theories. Our method is based on the minimization of a convex loss function constructed using the Schwinger-Dyson relations. We show that score matching, a popular learning method, is a special case of our construction of an infinite family of valid loss functions. Importantly, our general Schwinger-Dyson-based construction applies to gauge theories and models with Grassmann-valued fields used to represent dynamical fermions. In particular, we extend our method to realistic lattice field theories including quantum chromodynamics.

We propose a practical formulation of the overlap Dirac operator in lattice QCD that employs the diagonal Kenney-Laub rational iterates - expressed via their partial fraction decomposition - to approximate the matrix sign function. We investigate this approximation using the Brillouin operator as kernel, in addition to the standard Wilson Dirac operator. Numerical results show improved chiral symmetry preservation and computational efficiency compared to the Chebyshev polynomial approach.

We present the first principle calculation of charmed baryon decay constants employing 2+1 flavor gauge ensembles with lattice spacings ranging from 0.05 to 0.1 fm and pion masses between 136 and 310 MeV. Under $SU(3)$ flavor symmetry, we construct the charmed baryon interpolating operators and compute the corresponding hadronic matrix elements to extract the bare decay constants for each ensemble. The non-perturbative renormalization is performed via the symmetric momentum-subtraction scheme. After performing systematic chiral and continuum extrapolations, we obtain the decay constants with a precision of $8\sim 16\%$ from first principles.

We present a key algorithmic improvement to the generalized combinatorial Feynman--Vdovichenko method for calculating the critical temperature of the Ising model on Sierpiński carpets $SC_k(a,b)$, originally introduced in {\tt arXiv:1505.02699}. By reformulating the method in terms of purely real-valued transfer matrices, we substantially reduce their dimension. This optimization, together with modern computational resources, enables us to reach generation $k=10$ for the canonical $SC_k(3,1)$ carpet. Extrapolation from these data yields the most accurate estimate to date of the critical temperature $T_c^{(3,1)} = 1.4782927(26)$. We further extend the analysis to additional members of the $SC_k(a,b)$ family and report their corresponding critical temperatures.

We determine for the first time the two-loop renormalization-group (RG) equation for the nucleon light-cone distribution amplitude, which constitutes the last missing ingredient for the complete next-to-leading-logarithmic corrections to the nucleon form factors in the hard-collinear factorization framework. Applying the conformal expansion for this fundamental nucleon distribution amplitude then enables us to construct an analytic solution that captures the desired scale dependence of phenomenologically interesting series coefficients. Importantly, the two-loop RG evolutions of these central hadronic quantities can bring about noticeable impacts on the corresponding leading-logarithmic results for three sample models of the nucleon distribution amplitude.

In this chapter we introduce the $\theta$-dependence and the topological properties of QCD, features of the strongly interacting sector which give rise to the strong CP problem in the more general context of the Standard Model of particle physics. We discuss the analytical approaches that can be used to obtain qualitative, or in some cases quantitative, information on the $\theta$-dependence of QCD and QCD-like models, discussing their range of validity and comparing their predictions with the numerical results obtained by means of lattice simulations.