Source author record

Mattia Dalla Brida

Mattia Dalla Brida appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Non-perturbative thermal QCD at all temperatures: the case of mesonic screening masses

We present a strategy based on the step-scaling technique to study non-perturbatively thermal QCD up to very high temperatures. As a first concrete application, we compute the flavour non-singlet meson screening masses at 12 temperatures covering the range from $T \sim 1$ GeV up to $\sim 160$ GeV in the theory with three massless quarks. The calculation is carried out by Monte Carlo simulations on the lattice by considering large spatial extensions in order to have negligible finite volume effects. For each temperature we have simulated 3 or 4 values of the lattice spacing, so as to perform the continuum limit extrapolation with confidence at a few permille accuracy. Chiral symmetry restoration manifests itself in our results through the degeneracy of the vector and the axial vector channels and of the scalar and the pseudoscalar ones. In the entire range of temperatures explored, the meson screening masses deviate from the free theory result, $2 πT$, by at most a few percent. These deviations, however, cannot be explained by the known leading term in the QCD coupling constant $g$ up to the highest temperature, where other contributions are still very relevant. In particular the vector-pseudoscalar mass splitting turns out to be of $O(g^4)$ in the entire range explored, and it remains clearly visible up to the highest temperature, where the two screening masses are still significantly different within our numerical precision. The pattern of different contributions that we have found explains why it has been difficult in the past to match non-perturbative lattice results at $T \sim 1$ GeV with the analytic behaviour at asymptotically high temperatures.

preprint2021arXiv

A novel computational paradigm for a precise determination of the hadronic contribution to $(g_μ-2)$ from lattice QCD

The hadronic contribution to the muon anomalous magnetic moment $a_μ=(g_μ-2)/2$ has to be determined at the per-mille level for the Standard Model prediction to match the expected final uncertainty of the ongoing E989 experiment. That is 3 times better than the current precision from the dispersive approach, and 5-15 times smaller than the uncertainty based on the purely theoretical determinations from lattice QCD. So far the stumbling-block is the large statistical error in the Monte Carlo evaluation of the required correlation functions which can hardly be tamed by brute force. In this talk we present our proposal to solve this problem by multi-level Monte Carlo integration, a technique which reduces the variance of correlators exponentially in the distance of the fields. We report the results of our feasibility tests for the computation of the Hadronic Vacuum Polarization on a lattice with a linear extension of 3~fm, a spacing of 0.065 fm, and a pion mass of 270 MeV. Indeed the two-level integration makes the contribution to the statistical error from long-distances de-facto negligible by accelerating its inverse scaling with the cost of the simulation. These findings establish multi-level Monte Carlo as a solid and efficient method for a precise lattice determination of the hadronic contribution to $a_μ$.

preprint2021arXiv

Past, present, and future of precision determinations of the QCD coupling from lattice QCD

Non-perturbative scale-dependent renormalization problems are ubiquitous in lattice QCD as they enter many relevant phenomenological applications. They require solving non-perturbatively the renormalization group equations for the QCD parameters and matrix elements of interest in order to relate their non-perturbative determinations at low energy to their high-energy counterparts needed for phenomenology. Bridging the large energy separation between the hadronic and perturbative regimes of QCD, however, is a notoriously difficult task. In this contribution we focus on the case of the QCD coupling. We critically address the common challenges that state-of-the-art lattice determinations have to face in order to be significantly improved. In addition, we review a novel strategy that has been recently put forward in order to solve this non-perturbative renormalization problem and discuss its implications for future precision determinations. The new ideas exploit the decoupling of heavy quarks to match $N_{\rm f}$-flavor QCD and the pure Yang-Mills theory. Through this matching the computation of the non-perturbative running of the coupling in QCD can be shifted to the computationally much easier to solve pure-gauge theory. We shall present results for the determination of the $Λ$-parameter of $N_{\rm f}=3$-flavor QCD where this strategy has been applied and proven successful. The results demonstrate that these techniques have the potential to unlock unprecedented precision determinations of the QCD coupling from the lattice. The ideas are moreover quite general and can be considered to solve other non-perturbative renormalization problems.

preprint2020arXiv

Non-perturbative definition of the QCD energy-momentum tensor on the lattice

We present a strategy to define non-perturbatively the energy-momentum tensor in Quantum Chromodynamics (QCD) which satisfies the appropriate Ward identities and has the right trace anomaly. The tensor is defined by regularizing the theory on a lattice, and by fixing its renormalization constants non-perturbatively by suitable Ward identities associated to the Poincare' invariance of the continuum theory. The latter are derived in thermal QCD with a non-zero imaginary chemical potential formulated in a moving reference frame. A renormalization group analysis leads to simple renormalization-group-invariant definitions of the gluonic and fermionic contributions to either the singlet or the non-singlet components of the tensor, and therefore of their form factors among physical states. The lattice discussion focuses on the Wilson discretization of quark fields but the strategy is general. Specific to that case, we also carry out the analysis for the on-shell O(a)-improvement of the energy-momentum tensor. The renormalization and improvement programs profit from the fact that, as shown here, the thermal theory enjoys de-facto automatic O(a)-improvement at finite temperature. The validity of the proposal is scrutinized analytically by a study to 1-loop order in lattice perturbation theory with shifted and twisted (for quarks only) boundary conditions. The latter provides also additional useful insight for a precise non-perturbative calculation of the renormalization constants. The strategy proposed here is accessible to Monte Carlo computations, and in this sense it provides a practical way to define non-perturbatively the energy-momentum tensor in QCD.

preprint2019arXiv

Non-perturbative renormalization by decoupling

We propose a new strategy for the determination of the QCD coupling. It relies on a coupling computed in QCD with $N_{\rm f} \geq 3$ degenerate heavy quarks at a low energy scale $μ_{\rm dec}$, together with a non-perturbative determination of the ratio $Λ/μ_{\rm dec}$ in the pure gauge theory. We explore this idea using a finite volume renormalization scheme for the case of $N_{\rm f} = 3$ QCD, demonstrating that a precise value of the strong coupling $α_s$ can be obtained. The idea is quite general and can be applied to solve other renormalization problems, using finite or infinite volume intermediate renormalization schemes.

preprint2016arXiv

Determination of the QCD $Λ$-parameter and the accuracy of perturbation theory at high energies

We discuss the determination of the strong coupling $α_\mathrm{\overline{MS}}^{}(m_\mathrm{Z})$ or equivalently the QCD $Λ$-parameter. Its determination requires the use of perturbation theory in $α_s(μ)$ in some scheme, $s$, and at some energy scale $μ$. The higher the scale $μ$ the more accurate perturbation theory becomes, owing to asymptotic freedom. As one step in our computation of the $Λ$-parameter in three-flavor QCD, we perform lattice computations in a scheme which allows us to non-perturbatively reach very high energies, corresponding to $α_s = 0.1$ and below. We find that (continuum) perturbation theory is very accurate there, yielding a three percent error in the $Λ$-parameter, while data around $α_s \approx 0.2$ is clearly insufficient to quote such a precision. It is important to realize that these findings are expected to be generic, as our scheme has advantageous properties regarding the applicability of perturbation theory.

preprint2016arXiv

Importance of closed quark loops for lattice QCD studies of tetraquarks

To investigate the light scalar tetraquark candidate $a_0(980)$ (quantum numbers $J^P =0^+$), a correlation matrix including a variety of two- and four-quark interpolating operators has to be computed. We discuss efficient techniques to compute the elements of this correlation matrix, in particular diagrams with closed quark loops. Furthermore, we present evidence that such diagrams are not negligible given our precision, and their contribution is essential to obtain physically meaningful results. In particular, we find indications of the existence of an "additional" state around the two-particle thresholds of $K+ \bar K$ and $η+π$, which could correspond to the $a_0(980)$ meson.

preprint2016arXiv

Perturbative renormalization of $ΔF = 2$ four-fermion operators with the chirally rotated Schrödinger functional

The chirally rotated Schrödinger functional ($χ$SF) renders the mechanism of automatic $O(a)$ improvement compatible with Schrödinger functional (SF) renormalization schemes. Here we define a family of renormalization schemes based on the $χ$SF for a complete basis of $ΔF = 2$ parity-odd four-fermion operators. We compute the corresponding scale-dependent renormalization constants to one-loop order in perturbation theory and obtain their NLO anomalous dimensions by matching to the $\overline{\textrm{MS}}$ scheme. Due to automatic $O(a)$ improvement, once the $χ$SF is renormalized and improved at the boundaries, the step scaling functions (SSF) of these operators approach their continuum limit with $O(a^{2})$ corrections without the need of operator improvement.

preprint2016arXiv

The chirally rotated Schrödinger functional: theoretical expectations and perturbative tests

The chirally rotated Schrödinger functional ($χ$SF) with massless Wilson-type fermions provides an alternative lattice regularization of the Schrödinger functional (SF), with different lattice symmetries and a common continuum limit expected from universality. The explicit breaking of flavour and parity symmetries needs to be repaired by tuning the bare fermion mass and the coefficient of a dimension 3 boundary counterterm. Once this is achieved one expects the mechanism of automatic O($a$) improvement to be operational in the $χ$SF, in contrast to the standard formulation of the SF. This is expected to significantly improve the attainable precision for step-scaling functions of some composite operators. Furthermore, the $χ$SF offers new strategies to determine finite renormalization constants which are traditionally obtained from chiral Ward identities. In this paper we consider a complete set of fermion bilinear operators, define corresponding correlation functions and explain the relation to their standard SF counterparts. We discuss renormalization and O($a$) improvement and then use this set-up to formulate the theoretical expectations which follow from universality. Expanding the correlation functions to one-loop order of perturbation theory we then perform a number of non-trivial checks. In the process we obtain the action counterterm coefficients to one-loop order and reproduce some known perturbative results for renormalization constants of fermion bilinears. By confirming the theoretical expectations, this perturbative study lends further support to the soundness of the $χ$SF framework and prepares the ground for non-perturbative applications.

preprint2016arXiv

The determination of $α_s$ by the ALPHA collaboration

We review the ALPHA collaboration strategy for obtaining the QCD coupling at high scale. In the three-flavor effective theory it avoids the use of perturbation theory at $α> 0.2$ and at the same time has the physical scales small compared to the cutoff $1/a$ in all stages of the computation. The result $Λ_\overline{MS}^{(3)}=332(14)$~MeV is translated to $α_\overline{MS}(m_Z)=0.1179(10)(2)$ by use of (high order) perturbative relations between the effective theory couplings at the charm and beauty quark "thresholds". The error of this perturbative step is discussed and estimated as $0.0002$.

preprint2016arXiv

The gradient flow coupling from numerical stochastic perturbation theory

Perturbative calculations of gradient flow observables are technically challenging. Current results are limited to a few quantities and, in general, to low perturbative orders. Numerical stochastic perturbation theory is a potentially powerful tool that may be applied in this context. Precise results using these techniques, however, require control over both statistical and systematic uncertainties. In this contribution, we discuss some recent algorithmic developments that lead to a substantial reduction of the cost of the computations. The matching of the ${\overline{\rm MS}}$ coupling with the gradient flow coupling in a finite box with Schrödinger functional boundary conditions is considered for illustration.

preprint2015arXiv

A status update on the determination of $Λ_{\overline{\rm MS}}^{N_{\rm f}=3}$ by the ALPHA collaboration

The ALPHA collaboration aims to determine $α_s(m_Z)$ with a total error below the percent level. A further step towards this goal can be taken by combining results from the recent simulations of 2+1-flavour QCD by the CLS initiative with a number of tools developed over the years: renormalized couplings in finite volume schemes, recursive finite size techniques, two-loop renormalized perturbation theory and the (improved) gradient flow on the lattice. We sketch the strategy, which involves both the standard SF coupling in the high energy regime and a gradient flow coupling at low energies. This implies the need for matching both schemes at an intermediate switching scale, $L_{\rm swi}$, which we choose roughly in the range 2-4 GeV. In this contribution we present a preliminary result for this matching procedure, and we then focus on our almost final results for the scale evolution of the SF coupling from $L_{\rm swi}$ towards the perturbative regime, where we extract the $N_{\rm f} = 3$ $Λ$-parameter, $Λ_{\overline{\rm MS}}^{N_{\rm f}=3}$, in units of $L_{\rm swi}$ . Connecting $L_{\rm swi}$ and thus the $Λ$-parameter to a hadronic scale such as $F_K$ requires 2 further ingredients: first, the connection of $L_{\rm swi}$ to $L_{\rm max}$ using a few steps with the step-scaling function of the gradient flow coupling, and, second, the continuum extrapolation of $L_{\rm max} F_K$.

preprint2015arXiv

Computation of correlation matrices for tetraquark candidates with $J^P = 0^+$ and flavor structure $q_1 \bar{q_2} q_3 \bar{q}_3$

The conjecture that several recently observed mesons have a structure, which is not dominated by an ordinary quark-antiquark pair, but by a four-quark structure, is being actively investigated both theoretical and experimentally. Such a state may be characterized as a mesonic molecule or as a diquark-antidiquark pair. Lattice QCD provides a theoretically sound framework to study such states. To quantitatively investigate the internal structure of such mesons, one needs to precisely compute correlation matrices containing several interpolating operators including two and four quarks. Here we discuss certain technical aspects of such correlation matrices suited to study tetraquark candidates with $J^P = 0^+$ and flavor structure $q_1 \bar{q_2} q_3 \bar{q}_3$, e.g.\ the $a_0(980)$ meson, the $D_{s0}^\ast$ meson and some of the charged $c \bar{c}$ $X$ states. Some numerical results for the $a_0(980)$ meson are presented.

preprint2014arXiv

A dynamical study of the chirally rotated Schrödinger functional in QCD

The chirally rotated Schrödinger functional for Wilson-fermions allows for finite-volume, mass-independent renormalization schemes compatible with automatic O($a$) improvement. So far, in QCD, the set-up has only been studied in the quenched approximation. Here we present first results for $N_{\rm f} = 2$ dynamical quark-flavours for several renormalization factors of quark-bilinears. We discuss how these renormalization factors can be easily obtained from simple ratios of two-point functions, and show how automatic O($a$) improvement is at work. As a by-product of this investigation the renormalization of the non-singlet axial current, $Z_A$, is determined very precisely.

preprint2014arXiv

Investigation of the tetraquark candidate $a_0(980)$: technical aspects and preliminary results

We discuss technical aspects and first results of a lattice QCD study of the $a_0(980)$ state. We employ various interpolating operators of quark-antiquark, mesonic molecule, diquark-antidiquark and two-meson type. Both connected and disconnected contributions including diagrams with closed fermion loops are computed. To keep statistical errors small, it is essential to optimize the computation of these diagrams by choosing that combination of techniques most appropriate for each type of diagram from the correlation matrix of interpolating operators. We illustrate, how this can be done, by discussing certain diagrams in detail. We also present preliminary results corresponding to a $4\times 4$ submatrix computed with 2+1 flavors of clover fermions.

preprint2013arXiv

Numerical Stochastic Perturbation Theory in the Schrödinger Functional

The Schrödinger functional (SF) is a powerful and widely used tool for the treatment of a variety of problems in renormalization and related areas. Albeit offering many conceptual advantages, one major downside of the SF scheme is the fact that perturbative calculations quickly become cumbersome with the inclusion of higher orders in the gauge coupling and hence the use of an automated perturbation theory framework is desirable. We present the implementation of the SF in numerical stochastic perturbation theory (NSPT) and compare first results for the running coupling at two loops in pure SU(3) Yang-Mills theory with the literature.

preprint2013arXiv

Scalar mesons and tetraquarks from twisted mass lattice QCD

We study light scalar mesons with particular focus on the a_0(980) using lattice QCD with 2+1+1 dynamical quark flavors. To investigate the structure of these scalar mesons and to identify, whether a sizeable tetraquark component is present, we use a large set of operators, including diquark-antidiquark, mesonic molecule and two-meson operators. We find that the low-lying states overlap essentially exclusively with two-meson states. This indicates that in the channels investigated no tightly bound four-quark states of either molecular or diquark-antidiquark type exist.

preprint2012arXiv

Lattice investigation of the tetraquark candidates a0(980) and kappa

It is a long discussed issue whether light scalar mesons have sizeable four-quark components. We present an exploratory study of this question using N_f=2+1+1 twisted mass lattice QCD. A mixed action approach ignoring disconnected contributions is used to calculate correlator matrices consisting of mesonic molecule, diquark-antidiquark and two-meson interpolating operators with quantum numbers of the scalar mesons a0(980) (1(0^++)) and kappa (1/2(0+)). The correlation matrices are analyzed by solving the generalized eigenvalue problem. The theoretically expected free two-particle scattering states are identified, while no additional low lying states are observed. We do not observe indications for bound four-quark states in the channels investigated.

preprint2012arXiv

Scalar mesons and tetraquarks by means of lattice QCD

We study the light scalar mesons a_0(980) and kappa using N_f = 2+1+1 flavor lattice QCD. In order to probe the internal structure of these scalar mesons, and in particular to identify, whether a sizeable tetraquark component is present, we use a large set of operators, including diquark-antidiquark, mesonic molecule and two-meson operators. The inclusion of disconnected diagrams, which are technically rather challenging, but which would allow us to extend our work to e.g. the f_0(980) meson, is introduced and discussed.