High Energy Physics - Lattice (hep-lat)

Off-lightcone Wilson-line operators are constructed using local operators connected by time-like or space-like Wilson lines, which ensure gauge invariance. Off-lightcone Wilson-line operators have broad applications in various contexts. For instance, space-like Wilson-line operators play a crucial role in determining quasi-distribution functions (quasi-PDFs), while time-like Wilson-line operators are essential for understanding quarkonium decay and production within the potential non-relativistic QCD (pNRQCD) framework. In this work, we establish a systematic approach for calculating the matching from the gradient-flow scheme to the $\overline{\rm MS}$ scheme in the limit of small flow time for off-lightcone Wilson-line operators. By employing the one-dimensional auxiliary-field formalism, we simplify the matching procedure, reducing it to the matching of local current operators. We provide one-loop level matching coefficients for these local current operators. For the case of hadronic matrix element related to the quark quasi-PDFs, we show at one-loop level that the finite flow time effect is very small as long as the flow radius is smaller than the physical distance $z$, which is usually satisfied in lattice gradient flow computations. Applications include lattice gradient flow computations of quark/gluon quasi-PDFs, gluonic correlators related to quarkonium decay and production in pNRQCD, and spin-dependent potentials in terms of chromoelectric and chromomagnetic field insertions into a Wilson loop.

We consider two-dimensional ${\cal N} = (2, 2)$ Yang--Mills theory with gauge group SU($N$) in Euclidean signature compactified on a torus with thermal fermion boundary conditions imposed on one cycle. We perform non-perturbative lattice analyses of this theory for large $12 \leq N \leq 20$. Although no holographic dual of this theory is yet known, we conduct numerical investigations to check for features similar to the two-dimensional ${\cal N} = (8, 8)$ Yang--Mills theory, which has a well-defined gravity dual. We perform lattice field theory calculations to determine the phase diagram, observing a spatial deconfinement transition similar to the maximally supersymmetric case. However, the transition does not continue to low temperature, implying the absence of a topology-changing transition between black hole geometries in any holographic dual for this four-supercharge theory.

We study nonequilibrium dynamics of relativistic $N$-component scalar field theories in Minkowski space-time in a classical-statistical regime, where typical occupation numbers of modes are much larger than unity. In this strongly correlated system far from equilibrium, the role of different phenomena such as nonlinear wave propagation and defect dynamics remains to be clarified. We employ persistent homology to infer topological features of the nonequilibrium many-body system for different numbers of field components $N$ via a hierarchy of cubical complexes. Specifically, we show that the persistent homology of local energy density fluctuations can give rise to signatures of self-similar scaling associated with topological defects, distinct from the scaling behaviour of nonlinear wave modes. This contributes to the systematic understanding of the role of topological defects for far-from-equilibrium time evolutions of nonlinear many-body systems.

To accelerate the HMC with field transformation, we consider a variant of the trivializing map, the decimation map, which can be regarded as a coarse-graining transformation. Using the 2D $U(1)$ pure gauge model, combined with the guided Monte Carlo algorithm, we show that the integrated autocorrelation time of the topological charge can be exponentially improved in the wall clock time. Our study indicates that incorporating renormalization group picture is a powerful and essential ingredient to accelerate the HMC at large $\beta$.

We derive an explicit expression for the Lellouch-Lüscher (LL) factor in the $K\to 3\pi$ decays at leading order (without derivative couplings). Several important technical details are addressed, like a proper decomposition into the isospin amplitudes, the choice of a minimal set of effective couplings and the renormalization, as well as the algorithm for the solution of the pertinent Faddeev equations in the infinite volume which is based on the contour deformation method. Most importantly, our numerical results demonstrate that the three-body force contributes very little to the LL factor. This result paves the way for the study of the $K\to 3\pi$ decays on the lattice.

The TMD soft function can be obtained by formulating the Wilson line in terms of auxiliary 1-dimensional fermion fields on the lattice. In this formulation, the directional vector of the auxiliary field in Euclidean space has the form $\tilde n = (in^0, \vec 0_\perp, n^3)$, where the time component is purely imaginary. The components of these complex directional vectors in the Euclidean space can be mapped directly to the rapidities of the Minkowski space soft function. We present the results of the one-loop calculation of the Euclidean space analog to the soft function using these complex directional vectors. As a result, we show that the calculation is valid only when the directional vectors obey the relation: $|r| = |n^3/n^0| > 1$, and that this result corresponds to a computation in Minkowski space with space-like directed Wilson lines. Finally, we show that a lattice calculable object can be constructed that has the desired properties of the soft function.

We study the most general Two Higgs Doublet Model with $SU(2)$ gauge fields on the lattice. The phase space is probed through the computation of gauge-invariant global observables serving as proxies for order parameters. In each phase, the spectrum of the theory is analysed for different combinations of bare couplings and different symmetry breaking patterns. The scale setting and determination of the running gauge coupling are performed through the Wilson flow computation of the action density.

Until recently the spin-flip processes in the deep inelastic scatterings are thought to be suppressed in the high energy. We found a positive intercept for the spin-flip generalized transverse momentum-dependent parton distribution (GTMDs) ${\rm Re}(F_{1,2})$ as, \begineqnarray \rm Re(F_1,2) ∼\left(\frac1x\right)^\bar \alpha_s\left(4\ln2-8/3\right) \left(\cos 3\phi_k∆ +\cos \phi_k∆\right). \nonumber \endeqnarray This is done by analytically solving the integro-differential evolution equation for ${\rm Re}(F_{1,2})$, recently proposed by Hatta and Zhou, in the dilute regime. Interestingly, the surviving solution corresponds to conformal spin $n=2$ and carries an explicit $\cos 3\phi_{k\Delta} + \cos \phi_{k\Delta}$ azimuthal dependence. As the imaginary part of $F_{1,2}$, is related to the spin-dependent odderon or gluon Siver function and scales as ${\rm Im}(F_{1,2}) \sim x^{0}$, the positive intercept for ${\rm Re}(F_{1,2})$, implies that it is expected to dominate over the gluon Siver function in the small-$x$ limit - and may directly impact the modeling of unpolarised GTMDs and associated spin-flip processes.

The phenomenon of unpaired Weyl fermions appearing on the sole 2n-dimensional boundary of a (2n+1)-dimensional manifold with massive Dirac fermions was recently analyzed in a companion paper by one of the authors. In this Letter we show that similar unpaired Weyl edge states can be seen on a finite lattice. In particular, we consider the discretized Hamiltonian for a Wilson fermion in (2+1) dimensions with a 1+1 dimensional boundary and continuous time. We demonstrate that the low lying boundary spectrum is indeed Weyl-like: it has a linear dispersion relation and definite chirality and circulates in only one direction around the boundary. We comment on how our results are consistent with Nielsen-Ninomiya theorem. This work removes one potential obstacle facing the program for regulating chiral gauge theories recently proposed by one of the authors.

Effective string theory has shown its universal power in the prediction of the spectrum of low-lying excited states of confining strings. Here we study confining flux tubes in $\mathbb{Z}_N$ gauge theories. For the $N=2$ theory, which corresponds to the 3d Ising gauge model, we compute the spectrum of low-lying excitations of confining strings and show that it agrees with the universal Nambu--Goto predictions except for an additional massive scalar resonance. This resonance, however, turns out to be a bulk glueball mixing with the flux tube excitations rather than a genuine string worldsheet state. In general $\mathbb{Z}_N$ gauge theories (dual to clock spin models), we observe a continuous phase transition for $N \geq 4$, while for $N > 5$ it is governed by the $O(2)$ universality class. The critical behavior of the string tension and mass gap is verified to be described by a dangerously irrelevant operator. At large $N$ the glueball spectrum is expected to approach the spectrum of U(1) gauge theory, which is confirmed by our lattice data.

Using an effective theory for an ensemble of center vortices, we observe the area law fall-off in the continuum limit for the SU($2$) gauge group in three-dimensional Euclidean space-time. The string tension is obtained in terms of the intrinsic properties of the vortices and the parameter describing their interactions. A good qualitative agreement between our results and the lattice ones is observed. In addition, we show that the repulsive force between the vortices increases with temperature. This behavior is expected due to the reduction of vortex configuration at higher temperatures, required for the deconfinement regime.

Reconstructed-correlator methods have been used to investigate thermal effects in mesonic correlation functions in a fit-independent manner. This technique has recently been extended to the baryonic sector. In this work different ways of implementing this approach for baryon correlators are examined. Using both real and synthetic data it is found that for heavy baryons, such as the $\Xi_{cc}(ccu)$ baryon, different choices are equivalent and that for the lighter nucleon the effect of different implementations is minimal. Further comparison to the so-called "double ratio" using the FASTSUM Generation 2L thermal ensembles shows that reconstructed-correlator ratios and double ratios contain nearly identical quantitative information.

We present a numerical calculation of the Lee-Yang and Fisher zeros of the 2D Ising model using multi-point Padé approximants. We perform simulations for the 2D Ising model with ferromagnetic couplings both in the absence and in the presence of a magnetic field using a cluster spin-flip algorithm. We show that it is possible to extract genuine signature of Lee Yang and Fisher zeros of the theory through the poles of magnetization and specific heat, using multi-point Padé method. We extract the poles of magnetization using Padé approximants and compare their scaling with known results. We verify the circle theorem associated to the well known behaviour of Lee Yang zeros. We present our finite volume scaling analysis of the zeros done at $T=T_c$ for a few lattice sizes, extracting to a good precision the (combination of) critical exponents $\beta \delta$. The computation at the critical temperature is performed after the latter has been determined via the study of Fisher zeros, thus extracting both $\beta_c$ and the critical exponent $\nu$. Results already exist for extracting the critical exponents for the Ising model in 2 and 3 dimensions making use of Fisher and Lee Yang zeros. In this work, multi-point Padé is shown to be competitive with this respect and thus a powerful tool to study phase transitions.

Real-time quantities like spectral functions and transport coefficients are crucial for a proper understanding of the quark-gluon plasma created in relativistic heavy-ion collisions. Their numerical computation is plagued by a severe sign problem inherent in the real-time formulation of lattice field theories. In this letter, we present the first direct ab-initio computation of unequal-time correlation functions in non-Abelian lattice gauge theory, which are necessary to extract real-time quantities. We demonstrate non-trivial consistency relations among correlators, time-translation invariance, and agreement with Monte-Carlo results for thermal equilibrium in 3+1 dimensions by employing our stabilized complex Langevin method. Our work sets the stage to extract real-time observables in lattice gauge theory in a first-principles real-time framework.

A modified deep convolutional generative adversarial network (M-DCGAN) frame is proposed to study the N-dimensional (ND) topological quantities in lattice QCD based on the Monte Carlo (MC) simulations. We construct a new scaling structure including fully connected layers to support the generation of high-quality high-dimensional images for the M-DCGAN. Our results show that the M-DCGAN scheme of the Machine learning should be helpful for us to calculate efficiently the 1D distribution of topological charge and the 4D topological charge density compared with the case by the MC simulation alone.

We report recent progress on $\varepsilon_K$ evaluated directly from the standard model (SM) with lattice QCD inputs such as $\hat{B}_K$, $|V_{cb}|$, $|V_{us}|$, $|V_{ud}|$, $\xi_0$, $\xi_2$, $\xi_\text{LD}$, $f_K$, and $m_c$. We find that the standard model with exclusive $|V_{cb}|$ and lattice QCD inputs describes only 66\% of the experimental value of $|\varepsilon_K|$ and does not explain its remaining 34\%, which corresponds to a strong tension in $|\varepsilon_K|$ at the $4.9\sigma \sim 3.9\sigma$ level between the SM theory and experiment. We also find that this tension disappears when we use the inclusive value of $|V_{cb}|$ obtained using the heavy quark expansion based on the QCD sum rule approach.

In this paper we study the Chiral Separation Effect by means of first-principles lattice QCD simulations. For the first time in the literature, we determine the continuum limit of the associated conductivity using 2+1 flavors of dynamical staggered quarks at physical masses. The results reveal a suppression of the conductivity in the confined phase and a gradual enhancement toward the perturbative value for high temperatures. In addition to our dynamical setup, we also investigate the impact of the quenched approximation on the conductivity, using both staggered and Wilson quarks. Finally, we highlight the relevance of employing conserved vector and anomalous axial currents in the lattice simulations.

We present the first determination of the energy dependence of the $B$-$\bar{D}$ and $B^*$-$\bar{D}$ isospin-0, $S$-wave scattering amplitudes both below and above the thresholds using lattice QCD, which allows us to investigate rigorously whether mixed bottom-charm $\bar{b}\bar{c}ud$ tetraquarks exist as bound states or resonances. The scattering phase shifts are obtained using Lüscher's method from the energy spectra in two different volumes. To ensure that no relevant energy level is missed, we use large, symmetric $7 \times 7$ and $8 \times 8$ correlation matrices that include, at both source and sink, $B^{(*)}$-$\bar{D}$ scattering operators with the lowest three or four possible back-to-back momenta in addition to local $\bar{b}\bar{c}ud$ operators. We fit the energy dependence of the extracted scattering phase shifts using effective-range expansions. We observe sharp peaks in the $B^{(*)}$-$\bar{D}$ scattering rates close to the thresholds, which are associated with shallow bound states, either genuine or virtual, a few MeV or less below the $B^{(*)}$-$\bar{D}$ thresholds. In addition, we find hints for resonances with masses of order $100$ MeV above the thresholds and decay widths of order $200$ MeV.

We study the effects of taste degeneracy on the continuum scaling of the localization properties of the staggered Dirac operator in high-temperature QCD using numerical simulations on the lattice, focusing in particular on the position of the mobility edge separating localized and delocalized modes at the low end of the spectrum. We find that, if the continuum limit is approached at fixed spatial volume, the restoration of taste symmetry leads to sizeable systematic effects on estimates for the mobility edge obtained from spectral statistics, which become larger and larger as the lattice spacing is decreased. Such systematics, however, are found to decrease if the volume is increased at fixed lattice spacing. We argue that spectral statistics estimate correctly the position of the mobility edge in the thermodynamic limit at fixed spacing, and support this with an independent numerical analysis based directly on the properties of the Dirac eigenvectors, that are unaffected by taste degeneracy. We then provide a theoretical argument justifying the exchange of the thermodynamic and continuum limits when studying localization. This allows us to use spectral statistics to determine the position of the mobility edge, and to obtain a controlled continuum extrapolation of the mobility edge for the first time.

We derive an expression for the energy-momentum tensor in the discrete lattice formulation of pure glue QCD. The resulting expression satisfies the continuity equation for energy conservation up to numerical errors with a symmetric procedure for the time discretization. In the case of the momentum conservation equation, we obtain an expression that is of higher accuracy in lattice spacing ($\mathcal{O}(a^2)$) than the naive discretization where fields in the continuum expressions are replaced by discretized counterparts. The improvements are verified by performing numerical tests on the derived expressions using classical real-time lattice gauge theory simulations. We demonstrate substantial reductions in relative error of one to several orders of magnitude compared to a naive discretization for both energy and momentum conservation equations. We expect our formulation to have applications in the area of pre-equilibrium dynamics in ultrarelativistic heavy ion collisions, in particular for the extraction of transport coefficients such as shear viscosity.

Representing the parton distribution functions (PDFs) of the proton and other hadrons through flexible, high-fidelity parametrizations has been a long-standing goal of particle physics phenomenology. This is particularly true since the chosen parametrization methodology can play an influential role in the ultimate PDF uncertainties as extracted in QCD global analyses; these, in turn, are often determinative of the reach of experiments at the LHC and other facilities to non-standard physics, including at large $x$, where parametrization effects can be significant. In this study, we explore a series of encoder-decoder machine-learning (ML) models with various neural-network topologies as efficient means of reconstructing PDFs from meaningful information stored in an interpretable latent space. Given recent effort to pioneer synergies between QCD analyses and lattice-gauge calculations, we formulate a latent representation based on the behavior of PDFs in Mellin space, i.e., their integrated moments, and test the ability of various models to decode PDFs from this information faithfully. We introduce a numerical package, $\texttt{PDFdecoder}$, which implements several encoder-decoder models to reconstruct PDFs with high fidelity and use this tool to explore strengths and pitfalls of neural-network approaches to PDF parametrization. We additionally dissect patterns of learned correlations between encoded Mellin moments and reconstructed PDFs which suggest opportunities for further improvements to ML-based approaches to PDF parametrizations and uncertainty quantification.

We propose a method to prepare Gaussian wave packets with momentum on top of the interacting ground state of a fermionic Hamiltonian. Using Givens rotation, we show how to efficiently obtain expectation values of observables throughout the evolution of the wave packets on digital quantum computers. We demonstrate our technique by applying it to the staggered lattice formulation of the Thirring model and studying the scattering of two wave packets. Monitoring the the particle density and the entropy produced during the scattering process, we characterize the phenomenon and provide a first step towards studying more complicated collision processes on digital quantum computers. In addition, we perform a small-scale demonstration on IBM's quantum hardware, showing that our method is suitable for current and near-term quantum devices.

Using a simplified lattice version of the electroweak sector of the standard model, with dynamical fermions excluded, we determine at fixed Weinberg angle the transition line between the confined phase and the Higgs phase, the latter defined as the region where the global center subgroup of the gauge group is spontaneously broken, and "separation of charge" confinement disappears. We then search, via lattice Monte Carlo simulations, for possible neutral vector bosons in the Higgs region, apart from the photon and Z. There are numerical indications of a "light Z" in the lattice data (along with the photon and the Z), but a lack of the expected scaling of the light mass particle excludes any firm conclusions about the physical spectrum.

We present methods to constrain fermionic condensates on the level of the path integral, which grant access to the quantum effective potential in the infinite volume limit. In the case of a spontaneously broken symmetry, this potential possesses a manifestly flat region, which is inaccessible to the standard approach on the lattice. However, by constraining the appropriate order parameters such as the chiral condensate, one is then able to probe the flat region. We demonstrate our method of constraining fermionic condensates in the 2-dimensional Gross-Neveu model, which exhibits a spontaneously broken chiral symmetry. We show how the potential flattens for increasing volume and that the flat region is dominated by inhomogeneous field configurations.

We discuss a novel effective-field-theory-based approach for extracting two-body scattering information from finite volume energies, serving as an alternative to Lüscher's method. By explicitly incorporating one-pion exchange, we overcome the challenging left-hand cut problem in Lüscher's method and can handle finite volume energy levels both below and above the left-hand cut. Applied to the lattice data for $DD^*$ scattering at a pion mass of 280 MeV, as an illustrative example, our results reveal the significant impact of the one-pion exchange on P-wave and S-wave phase shifts. The pole position of the $T_{cc}(3875)^+$ state, extracted from the finite-volume energy levels at this pion mass while taking into account left-hand cut effects, range corrections and partial-wave mixing, is consistent with a near-threshold resonance. This study demonstrates, for the first time, that two-body scattering information can be reliably extracted from lattice spectra including the left-hand cut.

We propose a lattice-QCD-suitable framework for computing the two-photon long-distance contribution to the complex $K_{\rm L}\rightarrow\mu^+\mu^-$ decay amplitude, where QED is treated perturbatively in the continuum and infinite-volume. We provide preliminary numerical results on the quark-connected diagrams on one ensemble at physical pion mass from this method, with well-controlled systematic errors. The successful application of this method will allow the determination of the dispersive part of the aforementioned contribution from first-principles and enable a meaningful comparison between the Standard-Model prediction and experiment.

Preliminary lattice QCD results for $D\pi$ scattering in isospin $I=\frac{1}{2}$ channel are presented. Utilizing the $N_f=2+1$ Wilson-Clover configuration at two volumes ($L^3 \times T=32^3 \times 96$ and $48^3 \times 96$) with the same lattice spacing ($a=0.07746(18)$ fm) and pion mass ($m_\pi \approx 303$ MeV), various two-particle operators in both the COM and the moving frames are constructed and the corresponding finite-volume spectra are determined from their correlation functions. The $S$ and $P$-wave scattering phase shifts are then extracted using the Lüscher approach, assuming negligible contributions from higher partial waves. A virtual state associated with the $D_0^*$ is also identified.

As the study of three-hadron physics from lattice QCD matures, it is necessary to develop proper analysis tools in order to reliably study a variety of phenomena, including resonance spectroscopy and nuclear structure. Reconstructing the three-particle scattering amplitude requires solving integral equations, which can be written in terms of data-constrained dynamical functions and physical on-shell quantities. The driving term in these equations is the so-called one-particle exchange, which leads to a kinematic divergence for particles on-mass-shell. A vital component in defining three-particle amplitudes with definite parity and total angular momentum, which are used in spectroscopic studies, is to project the one-particle exchange into definite partial waves. We present a general procedure to construct exact analytic partial wave projections of the one-particle exchange contribution for any system composed of three spinless hadrons. Our result allows one full control over the analytic structure of the projection, which we explore for some low-lying partial waves with applications to three pions.

Dilepton decays of pseudoscalar mesons have been drawing particular interest, thanks to their sensitivity to both the QCD dynamics at low energy and also signals beyond the Standard Model. In this context, we shortly review our recent study on an improved Standard-Model prediction for the rare decay $\pi^0\to e^+e^-$, and compare it with the first determination on the lattice that predicted also the $\pi^0\to \gamma\gamma $ decay width as a byproduct. In addition, we discuss our recent work on $K_L\to\ell^+\ell^-$ decays and its connection to lattice QCD. We comment on the current uncertainty estimates and discuss how they could be improved profiting from future experiments and progress in lattice QCD.

We present an evaluation of the glueball spectrum for configurations produced with $N_f=1$ dynamical fermions as a function of the $m_{\rm PCAC}$ mass. We obtained masses of states that fall into the irreducible representations of the octahedral group of rotations in combination with the quantum numbers of charge conjugation $C$ and parity $P$. Due to the low signal to noise ratio, practically, we can only extract masses for the irreducible representations $R^{PC}=$ $A_1^{++}$, $E^{++}$, $T_2^{++}$ as well as $A_1^{-+}$. We make use of the Generalized Eigenvalue Problem (GEVP) with an operator basis consisting only of gluonic operators. Throughout this work we are aiming towards the identification of the effects of light dynamical quarks on the glueball spectrum and how this compares to the statistically more precise spectrum of SU(3) pure gauge theory. We used large gauge ensembles which consist of ${\sim {~\cal O}}(10 {\rm K})$ configurations. Our findings demonstrate that the low-lying spectrum of the scalar, tensor as well as pseudo-scalar glueballs receive negligible contributions from the inclusion of $N_f=1$ dynamical fermions.

Neutrino oscillation experiments require accurate reconstructions of neutrino energies, which depend in part on a theoretical understanding of the axial $N \rightarrow \Delta$ transition form factors. Lattice QCD studies of this transition require construction of all hadronic states with energies up to the mass of the $\Delta$ resonance, which includes $N\pi$ and $N\pi\pi$ for physical quark mass values. Building interpolating operators from sparse grids of source and sink points is a versatile method of approximating all-to-all quark propagators that has been successfully used in other multi-hadron calculations. This work will discuss applications of this method to $N\pi$ and $N\pi\pi$ systems and present preliminary results.

We set up a numerical S-matrix bootstrap problem to rigorously constrain bound state couplings given by the residues of poles in elastic amplitudes. We extract upper bounds on these couplings that follow purely from unitarity, crossing symmetry, and the Roy equations within their proven domain of validity. First we consider amplitudes with a single spin 0 or spin 2 bound state, both with or without a self-coupling. Subsequently we investigate amplitudes with the spectrum of bound states corresponding to the estimated glueball masses of pure SU(3) Yang-Mills. In the latter case the 'glue-hedron', the space of allowed couplings, provides a first-principles constraint for future lattice estimates.

The two-particle finite-volume scattering formalism derived by Lüscher and generalized in many subsequent works does not hold for energies far enough below the two-particle threshold to reach the nearest left-hand cut. The breakdown of the formalism is signaled by the fact that a real scattering amplitude is predicted in a regime where it should be complex. In this work, we address this limitation by deriving an extended formalism that includes the nearest branch cut, arising from single particle exchange. We focus on two-nucleon ($NN \to NN$) scattering, for which the cut arises from pion exchange, but give expressions for any system with a single channel of identical particles. The new result takes the form of a modified quantization condition that can be used to constrain an intermediate K-matrix in which the cut is removed. In a second step, integral equations, also derived in this work, must be used to convert the K-matrix to the physical scattering amplitude. We also show how the new formalism reduces to the standard approach when the $N \to N \pi$ coupling is set to zero.

We describe a procedure to determine moments of parton distribution functions of any order in lattice quantum chromodynamics (QCD). The procedure is based on the gradient flow for fermion and gauge fields. The flowed matrix elements of twist-2 operators renormalize multiplicatively, and the matching with the physical matrix elements can be obtained using continuum symmetries and the irreducible representations of Euclidean 4-dimensional rotations. We calculate the matching coefficients at one-loop in perturbation theory for moments of any order in the flavor nonsinglet case. We also give specific examples of operators that could be used in lattice QCD computations. It turns out that it is possible to choose operators with identical Lorentz indices and still have a multiplicative matching. One can thus use twist-2 operators exclusively with temporal indices, thus substantially improving the signal-to-noise ratio in the computation of the hadronic matrix elements.

TEK reduction is a well-established technique that allows single-site simulations of Yang-Mills theory in the large-$N_c$ limit by exploiting volume reduction induced by twisted boundary conditions. We performed simulations for $SU(841)$ for several gauge couplings and applied standard Wilson flow techniques combined with a tree-level improvement methodology to set the lattice scale. The wide range of gauge couplings covered by our simulations allows us to explore a region in the coupling space where our data exhibits asymptotic scaling and perturbation theory could be used to analyze the behaviour of the $\beta$-function. In this talk, I will review the methodology used and go through the main results we obtained, including a determination of the $\Lambda$-parameter of Yang-Mills theory at large-$N_c$ in $\overline{\text{MS}}$-scheme.

In this work we consider strongly interacting dark matter candidates as composite states of $N_f=2$ fermions charged under a dark $Sp(4)$ gauge group in the fundamental representation. We give expressions that allow the calculation of correlation functions of two pseudo-Nambu-Goldstone-bosons with lattice field theory and present first results on the scattering phase shift in the isospin-2 channel in the theory from first principles. We give a lower limit on the dark matter particle mass by comparing our results with astrophysical constraints on the cross-section.

Understanding the intricate three-dimensional internal structure of the nucleon has been a long-standing challenge. The main quantitative tool to map out this structure are the generalized parton distributions (GPDs). In these proceedings, we present our exploratory results from applying the pseudo-distribution approach on the lattice, showing the unpolarized isovector $E$ GPD as an example. We use one ensemble of $N_f=2+1+1$ twisted mass fermions with a clover term, at a non-physical pion mass of $260 \textrm{ MeV}$ and a lattice spacing of $0.093 \textrm{ fm}$.

The numerical sign problem poses a seemingly insurmountable barrier to the simulation of many fascinating systems. We apply neural networks to deform the region of integration, mitigating the sign problem of systems with strongly correlated electrons. In this talk we present our latest architectural developments as applied to contour deformation. We also demonstrate its applicability to real systems, namely perylene.

The extraction of spectral densities from Euclidean correlators evaluated on the lattice is an important problem, as these quantities encode physical information on scattering amplitudes, finite-volume spectra, inclusive decay rates, and transport coefficients. In this contribution, we show that the Bayesian approach to this "inverse" problem, based on Gaussian processes, can be reformulated in a way that yields a solution equivalent, up to statistical uncertainties, to the one obtained in a Backus-Gilbert approach. After discussing this equivalence, we point out its implications for a reliable determination of spectral densities from lattice simulations.

We present one-loop perturbative results of the renormalization functions for a complete set of nonlocal quark bilinear operators containing an asymmetric staple-shaped Wilson line, using a family of improved lattice actions. This study is relevant for the nonperturbative investigations regarding the renormalization of the unpolarized, helicity and transversity transverse-momentum dependent parton distribution functions (TMDPDFs) in lattice QCD. We employ a number of different versions of regularization-independent (RI$'$) renormalization prescriptions which address the power and logarithmic divergences of such nonlocal operators, the pinch-pole singularities at infinite Wilson-line lengths, as well as the mixing among operators of different Dirac structures, as dictated by discrete symmetries. All cancelations of divergences and admixtures are confirmed by our results at one-loop level. We compare all the different prescriptions and we provide the conversion matrices at one-loop order which relate the matrix elements of the staple operators in RI$'$ to the reference scheme $\overline{\rm MS}$.

Lattice gauge-equivariant convolutional neural networks (L-CNNs) can be used to form arbitrarily shaped Wilson loops and can approximate any gauge-covariant or gauge-invariant function on the lattice. Here we use L-CNNs to describe fixed point (FP) actions which are based on renormalization group transformations. FP actions are classically perfect, i.e., they have no lattice artifacts on classical gauge-field configurations satisfying the equations of motion, and therefore possess scale invariant instanton solutions. FP actions are tree-level Symanzik-improved to all orders in the lattice spacing and can produce physical predictions with very small lattice artifacts even on coarse lattices. We find that L-CNNs are much more accurate at parametrizing the FP action compared to older approaches. They may therefore provide a way to circumvent critical slowing down and topological freezing towards the continuum limit.

We compute the screening masses of fields with nucleon quantum numbers for a wide range of temperatures between $T \sim 1$ GeV and $T\sim 160$ GeV. The computation has been performed by means of Monte Carlo simulations of lattice QCD with $N_f=3$ flavors of $O(a)$-improved Wilson fermions: we exploit a novel strategy which has recently allowed to determine for the first time non-singlet mesonic screening masses up to extremely high temperatures. The baryonic screening masses are measured with a few per-mille precision in the continuum limit, and percent deviations from the free theory result $3\pi T$ are clearly visible even at the highest temperatures. The observed degeneracy of the positive and negative parity state's screening mass, expected from Ward identities associated to non-singlet axial transformations, provides further evidence for the restoration of chiral symmetry in the high temperature regime of QCD.

We demonstrate the efficiency of the bond weighting method for the Grassmann tensor renormalization group (TRG). Benchmarking with the two-dimensional Gross-Neveu model with the Wilson fermion at finite density, we show that the bond weighting method improves the accuracy of the original Grassmann TRG. We also provide a sample code of the bond-weighted TRG that can be applied to the two-dimensional models including fermions on a square lattice.

In this study, we investigate the renormalization of a complete set of gauge-invariant non-local gluon operators up to one-loop in lattice perturbation theory. Our computations have been performed in both dimensional and lattice regularizations, using the Wilson gluon action, leading to the renormalization functions in the modified Minimal Subtraction $(\overline{\text{MS}})$ scheme, as well as conversion factors from the modified regularization invariant $(RI')$ scheme to $\overline{\text{MS}}$.

We study interactions between the S-wave octet pseudo-scalar (PS) meson and octet baryon in the flavor SU(3) limit using the HAL QCD method at the PS meson mass $m_M\approx 670~\textrm{MeV}$. We focus on the singlet and two octet channels, where the poles corresponding to $\Lambda(1405)$ have been predicted in the chiral unitary model. For calculations with $\Lambda$-baryon source operators with zero momentum, we employ the conventional stochastic calculation combined with the covariant-approximation averaging to calculate the all-to-all propagators. Due to a zero of the R-correlator (a kind of wave function), the leading order (LO) potential obtained by the single channel analysis has a singular point in all channels, which makes it difficult to obtain reliable binding energies. To overcome this problem, we take a linear combination of two octet R-correlators with a relative weight such that it does not cross zero, as two octet channels are suggested to couple to the same low-energy states with different weights. The potential calculated from such the linear combination shows strong attraction without singularities, though its shape depends on the relative weight. Our estimation for the binding energy in the octet channel is $E^{8_{s(a)}}_{\textrm{bind}}=163(7)(^{+16} _{-64})~\textrm{MeV}$, which is consistent with 156(8) MeV estimated from the two-point correlation function within errors.

The Landau gauge four-gluon vertex is studied using high statistical lattice simulations for several momentum configurations. Furthermore, the outcome of the lattice QCD simulations is compared with calculations performed with continuum Schwinger-Dyson equations.

The gradient-flow formalism proves to be a useful tool in lattice calculations of quantum chromodynamics. For example, it can be used as a scheme to renormalize composite operators by inverting the short-flow-time expansion of the corresponding flowed operators. In this paper, we consider the short-flow-time expansion of five quark bilinear operators, the scalar, pseudoscalar, vector, axialvector, and tensor currents, and compute the matching coefficients through next-to-next-to-leading order QCD. Among other applications, our results constitute one ingredient for calculating bag parameters of mesons within the gradient-flow formalism on the lattice.

We present a preliminary result of the kaon semileptonic form factor calculated at the smallest lattice spacing in the PACS10 configurations, whose physical volumes are more than (10 fm)$^4$ at the physical point. The configurations were generated using the Iwasaki gauge action and $N_f=2+1$ stout-smeared nonperturbatively $O(a)$ improved Wilson quark action at the three lattice spacings, 0.085, 0.063, and 0.041 fm. The value of $|V_{us}|$ in the continuum limit is estimated from our results including the preliminary one. We compare our result of $|V_{us}|$ with the previous results and those through the kaon leptonic decay.

In this work we perform calculations in order to determine the renormalization factors and the mixing coefficients of the Yukawa and the quartic couplings in $\mathcal{N} = 1$ Supersymmetric QCD. The Yukawa couplings describe the interactions between gluino, quark and squark fields whereas the quartic couplings describe four-squark interactions. We discretize the action on a Euclidean lattice using the Wilson formulation for the gluino, quark and gluon fields; for squark fields (scalar fields) we employ naïve discretization. At the quantum level Yukawa and quartic interactions suffer from mixing with other operators which have the same transformation properties. Exploiting parity and charge conjugation symmetries of the Supersymmetric QCD action, we reduce the allowed mixing patterns. We compute, perturbatively to one-loop and to the lowest order in the lattice spacing, the relevant three-point Green's functions so as to fine tune the Yukawa couplings and the relevant four-point Green's functions to fine tune the quartic couplings. We use both dimensional and lattice regularizations as required for implementing the Modified Minimal Subtraction scheme ($\overline{\rm MS}$).

This paper introduces a transverse-momentum dependent (TMD) factorization scheme designed to unify both large and small Bjorken-x regimes. We compute the next-to-leading order (NLO) quantum chromodynamics (QCD) corrections to the gluon TMD operator for an unpolarized hadron within this proposed scheme. This leads to the emergence of a new TMD evolution, incorporating those in transverse momentum, rapidity, and Bjorken-x. When matched to the collinear factorization scheme, our factorization scheme faithfully reproduces the well-established Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) and Collins-Soper-Sterman (CSS) evolutions. Conversely, matching with high-energy factorization not only yields the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution but also reveals distinctive signatures of CSS logarithms. The development of this novel TMD factorization scheme, capable of seamlessly reconciling disparate Bjorken-x regimes and faithfully reproducing established QCD evolution equations, has the potential to significantly advance our comprehension of high-energy processes and three-dimensional parton structures of hadrons.