Non-perturbative Renormalization Research Articles

Erratum: Nonperturbative renormalization in lattice QCD with three flavors of clover fermions: Using periodic and open boundary conditions [Phys. Rev. D 103 , 094511 (2021)

Erratum: Nonperturbative renormalization in lattice QCD with three flavors of clover fermions: Using periodic and open boundary conditions [Phys. Rev. D <b>103</b> , 094511 (2021)

Nonperturbative renormalization of the supercurrent in N=1 supersymmetric Yang-Mills theory

In this work, we study the nonperturbative renormalization of the supercurrent operator in $\mathcal{N}=1$ supersymmetric Yang-Mills theory, using a gauge-invariant renormalization scheme (GIRS). The proposed prescription addresses successfully the unwanted mixing of the supercurrent with other operators of equal or lower dimension, which respect the same global symmetries. This mixing is introduced by the unavoidable breaking of supersymmetry on the lattice. In GIRS all gauge-noninvariant operators, which mix with the supercurrent, are excluded from the renormalization procedure. The one remaining mixing operator is accessible by numerical simulations. We present results for the renormalization of the supercurrent using GIRS. We also compute at one-loop order the conversion matrix which relates the nonperturbative renormalization factors in GIRS to the reference scheme $\overline{\mathrm{MS}}$.

Low-temperature behavior of the O(N) models below two dimensions.

We investigate the critical behavior and the nature of the low-temperature phase of the O(N) models treating the number of field components N and the dimension d as continuous variables with a focus on the d≤2 and N≤2 quadrant of the (d,N) plane. We precisely chart a region of the (d,N) plane where the low-temperature phase is characterized by an algebraic correlation function decay similar to that of the Kosterlitz-Thouless phase but with a temperature-independent anomalous dimension η. We revisit the Cardy-Hamber analysis leading to a prediction concerning the nonanalytic behavior of the O(N) models' critical exponents and emphasize the previously not broadly appreciated consequences of this approach in d<2. In particular, we discuss how this framework leads to destabilization of the long-range order in favor of the quasi-long-range order in systems with d<2 and N<2. Subsequently, within a scheme of the nonperturbative renormalization group we identify the low-temperature fixed points controlling the quasi-long-range ordered phase and demonstrate a collision between the critical and the low-temperature fixed points upon approaching the lower critical dimension. We evaluate the critical exponents η(d,N) and ν^{-1}(d,N) and demonstrate a very good agreement between the predictions of the Cardy-Hamber type analysis and the nonperturbative renormalization group in d<2.

Improved renormalization scheme for nonlocal operators

In this paper we present an improved regularization-independent (RI)-type prescription appropriate for the nonperturbative renormalization of gauge-invariant nonlocal operators. In this prescription, the nonperturbative vertex function is improved by subtracting unwanted finite lattice spacing ($a$) effects, calculated in lattice perturbation theory. The method is versatile and can be applied to a wide range of fermion and gluon actions, as well as types of nonlocal operators. The presence of operator mixing can also be accommodated. Compared to the standard ${\mathrm{RI}}^{\ensuremath{'}}$ prescription, this variant can be recast as a supplementary finite renormalization, whose coefficients bring about corrections of higher order in $a$; consequently, it coincides with standard ${\mathrm{RI}}^{\ensuremath{'}}$ as $a\ensuremath{\rightarrow}0$---however, it can afford us a smoother and more controlled extrapolation to the continuum limit. In this proof-of-concept calculation we focus on nonlocal fermion bilinear operators containing a straight Wilson line. In the numerical implementation we use Wilson clover fermions and Iwasaki improved gluons. The finite-$a$ terms were calculated to one-loop level in lattice perturbation theory, and to all orders in $a$, using the same action as the nonperturbative vertex functions. We find that the method leads to significant improvement in the perturbative region indicated by small and intermediate values of the length of the Wilson line. This results in a robust extraction of the renormalization functions in that region. We have also applied the above method to operators with stout-smeared links. We show how to perform the perturbative correction for any number of smearing iterations and evaluate its effect on the power divergent renormalization coefficients.

Z_{4}-symmetric perturbations to the XY model from functional renormalization.

We employ the second order of the derivative expansion of the nonperturbative renormalization group to study cubic (Z_{4}-symmetric) perturbations to the classical XY model in dimensionality d∈[2,4]. In d=3 we provide accurate estimates of the eigenvalue y_{4} corresponding to the leading irrelevant perturbation and follow the evolution of the physical picture upon reducing spatial dimensionality from d=3 towards d=2, where we approximately recover the onset of the Kosterlitz-Thouless physics. We analyze the interplay between the leading irrelevant eigenvalues related to O(2)-symmetric and Z_{4}-symmetric perturbations and their approximate collapse for d→2. We compare and discuss different implementations of the derivative expansion in cases involving one and two invariants of the corresponding symmetry group.

Configuration space for quantum gravity in a locally regularized path integral

We discuss some aspects of the metric configuration space in quantum gravity in the background field formalism. We give a necessary and sufficient condition for the parametrization of Euclidean metric fluctuations such that (i) the signature of the metric is preserved in all configurations that enter the gravitational path integral, and (ii) the parametrization provides a bijective map between full Euclidean metrics and metric fluctuations about a fixed background. For the case of foliatable manifolds, we show how to parametrize fluctuations in order to preserve foliatability of all configurations. Moreover, we show explicitly that preserving the signature on the configuration space for the Lorentzian quantum gravitational path integral is most conveniently achieved by inequality constraints. We discuss the implementation of these inequality constraints in a nonperturbative renormalization group setup.

Non-perturbative renormalization of the QCD flavour-singlet local vector current

We compute non-perturbatively the renormalization constant of the flavour-singlet local vector current ZV in lattice QCD with 3 massless flavours. Gluons are discretized by the Wilson plaquette action while quarks by the O(a)-improved Wilson–Dirac operator. The constant ZV is fixed by comparing the expectation values (1-point functions) of the conserved and local vector currents at finite temperature in the presence of shifted boundary conditions and at non-zero imaginary chemical potential. Monte Carlo simulations with a moderate computational cost allow us to obtain ZV with an accuracy of about 8‰ for values of the inverse bare coupling constant β=6/g02 in the range 5.3≤β≤11.5.

Incompleteness of the large-N analysis of the O(N) models: Nonperturbative cuspy fixed points and their nontrivial homotopy at finite N.

We summarize the usual implementations of the large-N limit of O(N) models and show in detail why and how they can miss some physically important fixed points when they become singular in the limit N→∞. Using Wilson's renormalization group in its functional nonperturbative versions, we show how the singularities build up as N increases. In the Wilson-Polchinski version of the nonperturbative renormalization group, we show that the singularities are cusps, which become boundary layers for finite but large values of N. The corresponding fixed points being never close to the Gaussian, are out of reach of the usual perturbative approaches. We find four new fixed points and study them in all dimensions and for all N>0 and show that they play an important role for the tricritical physics of O(N) models. Finally, we show that some of these fixed points are bivalued when they are considered as functions of d and N thus revealing important and nontrivial homotopy structures. The Bardeen-Moshe-Bander phenomenon that occurs at N=∞ and d=3 is shown to play a crucial role for the internal consistency of all our results.

Asymptotically free and safe quantum gravity scenarios consistent with Hubble, laboratory, and inflation scale physics

To find out the possible scenarios for quantum gravity consistent with the observed universe, we numerically investigate the non-perturbative renormalization group equations of a general quadratic gravity theory recently derived by Sen, Wetterich and Yamada (\textit{JHEP} 03 (2022) 130). As boundary conditions, we impose consistency with the Hubble scale and the laboratory scale experiments, and the Starobinsky model of inflation. We find two kinds of trajectories which go to different regime at the trans-Planckian scales: i) a trajectory which flows to the asymptotically free regime, and ii) a trajectory which flows to the asymptotically safe regime. To determine the early-time cosmological scenario, an additional observational data from beyond the homogeneous and isotropic space-time is necessary.

Uncovering a chirally suppressed mechanism of 0νββ decay with LHC searches

∆L = 2 lepton number violation (LNV) at the TeV scale could provide an alternative interpretation of positive signal(s) in future neutrinoless double beta (0νββ) decay experiments. An interesting class of models from this point of view are those that at low energies give rise to dimension-9 vector operators and a dimension-7 operator, both of whose 0νββ-decay rates are “chirally suppressed”. We study and compare the sensitivities of 0νββ-decay experiments and LHC searches to a simplified model in this class of TeV-scale LNV that is also SU(2)L × U(1)Y gauge invariant. The searches for 0νββ decay, which are here diluted by a chiral suppression of the vector operators, are found to be less constraining than LHC searches whose reach is increased by the assumed kinematic accessibility of the mediator particles. For the chirally suppressed dimension-7 operator generated by TeV-scale mediators, in contrast, 0νββ-decay searches place strong constraints on the size of the new Yukawa coupling. Signals of this model at the LHC and 0νββ-decay experiments are entirely uncorrelated with the observed neutrinos masses, as these new sources of LNV give negligible contributions to the latter. We find the prospects for the high-luminosity LHC and ton-scale 0νββ-decay experiments to uncover the chirally-suppressed mechanism with TeV-scale LNV to be promising. We also comment on the sensitivity of the 0νββ-decay lifetime to certain unknown low-energy constants that in the case of dimension-9 scalar operators are expected to be large due to non-perturbative renormalization.

Absence versus Presence of Dissipative Quantum Phase Transition in Josephson Junctions.

Dissipative quantum phase transition has been widely believed to occur in a Josephson junction coupled to a resistor despite a lack of concrete experimental evidence. Here, on the basis of both numerical and analytical nonperturbative renormalization group analyses, we reveal breakdown of previous perturbative arguments and defy the common wisdom that the transition always occurs at the quantum resistance R_{Q}=h/(4e^{2}). We find that renormalization group flows in nonperturbative regimes induce nonmonotonic renormalization of the charging energy and lead to a qualitatively different phase diagram, where the insulator phase is strongly suppressed to the deep charge regime (Cooper pair box), while the system is always superconducting in the transmon regime. We identify a previously overlooked dangerously irrelevant term as an origin of the failure of conventional understandings. Our predictions can be tested in recent experiments realizing high-impedance long superconducting waveguides and would provide a solution to the long-standing controversy about the fate of dissipative quantum phase transition in the resistively shunted Josephson junction.

Renormalization of Transverse-Momentum-Dependent Parton Distribution on the Lattice.

To calculate the transverse-momentum-dependent parton distribution functions (TMDPDFs) from lattice QCD, an important goal yet to be realized, it is crucial to establish a viable nonperturbative renormalization approach for linear divergences in the corresponding Euclidean quasi-TMDPDF correlators in large-momentum effective theory. We perform a first systematic study of the renormalization property of the quasi-TMDPDFs by calculating the relevant matrix elements in a pion state at five lattice spacings ranging from 0.03fm to 0.12fm. We demonstrate that the square root of the Wilson loop combined with the short distance hadron matrix element provides a successful method to remove all ultraviolet divergences of the quasi-TMD operator, and thus provides the necessary justification to perform a continuum limit calculation of TMDPDFs. In contrast, the popular regularization independent momentum subtraction renormalization (RI/MOM) scheme fails to eliminate all linear divergences.

Lattice study of a magnetic contribution to heavy quark momentum diffusion

Heavy quarks placed within a hot QCD medium undergo Brownian motion, characterized by specific transport coefficients. Their determination can be simplified by expanding them in T/M, where T is the temperature and M is a heavy quark mass. The leading term in the expansion originates from the colour-electric part of a Lorentz force, whereas the next-to-leading order involves the colour-magnetic part. We measure a colour-magnetic 2-point correlator in quenched QCD at T ∼ (1.2 − 2.0)Tc. Employing multilevel techniques and non-perturbative renormalization, a good signal is obtained, and its continuum extrapolation can be estimated. Modelling the shape of the corresponding spectral function, we subsequently extract the momentum diffusion coefficient, κ. For charm (bottom) quarks, the magnetic contribution adds ∼ 30% (10%) to the electric one. The same increases apply also to the drag coefficient, η. As an aside, the colour-magnetic spectral function is computed at NLO.

Proton momentum and angular momentum decompositions with overlap fermions

We present a calculation of the proton momentum and angular momentum decompositions using overlap fermions on a $2+1$-flavor RBC/UKQCD domain-wall lattice at 0.143 fm with a pion mass of 171 MeV which is close to the physical one. A complete determination of the momentum and angular momentum fractions carried by up, down, strange, and glue inside the proton has been done with valence pion masses varying from 171 to 391 MeV. We have utilized fast Fourier transform on the stochastic-sandwich method for connected-insertion parts and the cluster-decomposition error reduction technique technique for disconnected-insertion parts has been used to reduce statistical errors. The full nonperturbative renormalization and mixing between the quark and glue operators are carried out. The final results are normalized with the momentum and angular momentum sum rules and reported at the physical valence pion mass at $\overline{\mathrm{MS}}(\ensuremath{\mu}=2\text{ }\text{ }\mathrm{GeV})$. The renormalized momentum fractions for the quarks and glue are $⟨x{⟩}^{q}=0.491(20)(23)$ and $⟨x{⟩}^{g}=0.509(20)(23)$, respectively, and the renormalized total angular momentum fractions for quarks and glue are $2{J}^{q}=0.539(22)(44)$ and $2{J}^{g}=0.461(22)(44)$, respectively. The quark spin fraction is $\mathrm{\ensuremath{\Sigma}}=0.405(25)(37)$ from our previous work and the quark orbital angular momentum fraction is deduced from $2{L}^{q}=2{J}^{q}\ensuremath{-}\mathrm{\ensuremath{\Sigma}}$ to be 0.134(22)(44).

Nonperturbative renormalization of the lattice Sommerfield vector model

The lattice Sommerfield model, describing a massive vector gauge field coupled to a light fermion in 2d, is an ideal candidate to verify perturbative conclusions. In contrast with continuum exact solutions, we prove that there is no infinite field renormalization, implying the reduction of the degree of the ultraviolet divergence, and that anomalies are non renormalized. Such features are the counterpart of analogue properties at the basis of the Standard model perturbative renormalizability. The results are non-perturbative, in the sense that the averages of gauge invariant observables are expressed in terms of convergent expansions uniformly in the lattice and volume.

Nonperturbative renormalization of quantum thermodynamics from weak to strong couplings

By solving the exact master equation of open quantum systems, we formulate the quantum thermodynamics from weak to strong couplings. The open quantum systems exchange matters, energies and information with their reservoirs through quantum particles tunnelings that are described by the generalized Fano-Anderson Hamiltonians. We find that the exact solution of the reduced density matrix of these systems approaches a Gibbs-type state in the steady-state limit for both the weak and strong system-reservoir coupling strengths. When the couplings become strong, thermodynamic quantities of the system must be renormalized. The renormalization effects are obtained nonperturbatively after exactly traced over all reservoir states through the coherent state path integrals. The renormalized system Hamiltonian is characterized by the renormalized system energy levels and interactions, corresponding to the quantum work done by the system. The renormalized temperature is introduced to characterize the entropy production counting the heat transfer between the system and the reservoir. We further find that only with the renormalized system Hamiltonian and other renormalized thermodynamic quantities, can the exact steady state of the system be expressed as the standard Gibbs state. Consequently, the corresponding exact steady-state particle occupations in the renormalized system energy levels obey the Bose-Einstein and the Fermi-Dirac distributions for bosonic and fermionic systems, respectively. Thus, the conventional statistical mechanics and thermodynamics are thereby rigorously deduced from quantum dynamical evolution.

Proton decay matrix elements on the lattice at physical pion mass

Proton decay is a major prediction of grand-unified theories (GUT) and its observation would indicate baryon number violation that is required for baryogenesis. Many decades of searching for proton decay have constrained its rate and ruled out some of the simplest GUT models. Apart from the baryon number-violating interactions, this rate also depends on transition amplitudes between the protons and mesons or leptons produced in the decay, which are matrix elements of three-quark operators. We report the nonperturbative calculation of these matrix elements for the most studied two-body decay channels into a meson and antilepton done on a lattice with physical light and strange quark masses and lattice spacings $a\ensuremath{\approx}0.14\text{ }\text{ }\mathrm{fm}$ and 0.20 fm. We perform nonperturbative renormalization and excited state analysis to control associated systematic effects. Our results largely agree with previous lattice calculations done with heavier quark masses and thus remove ambiguity in ruling out some simple GUT theories due to quark mass dependence of hadron structure.

Operator product expansion coefficients from the nonperturbative functional renormalization group

Using the nonperturbative functional renormalization group (FRG) within the Blaizot-M\'endez-Galain-Wschebor approximation, we compute the operator product expansion (OPE) coefficient $c_{112}$ associated with the operators $\mathcal{O}_1\sim\varphi$ and $\mathcal{O}_2\sim\varphi^2$ in the three-dimensional $\mathrm{O}(N)$ universality class and in the Ising universality class ($N=1$) in dimensions $2 \leq d \leq 4$. When available, exact results and estimates from the conformal bootstrap and Monte-Carlo simulations compare extremely well to our results, while FRG is able to provide values across the whole range of $d$ and $N$ considered.

Importance truncation for the in-medium similarity renormalization group

Ab initio nuclear many-body frameworks require extensive computational resources, especially when targeting heavier nuclei. Importance-truncation (IT) techniques allow one to significantly reduce the dimensionality of the problem by neglecting unimportant contributions to the solution of the many-body problem. In this work, we apply IT methods to the nonperturbative in-medium similarity renormalization group (IMSRG) approach and investigate the induced errors for ground-state energies in different mass regimes based on different nuclear Hamiltonians. We study various importance measures, which define the IT selection, and identify two measures that perform best, resulting in only small errors to the full IMSRG(2) calculations even for sizable compression ratios. The neglected contributions are accounted for in a perturbative way and serve as an estimate of the IT-induced error. Overall we find that the IT-IMSRG(2) performs well across all systems considered, while the largest compression ratios for a given error can be achieved when using soft Hamiltonians and for large single-particle bases.

Gaussian information bottleneck and the non-perturbative renormalization group

The renormalization group (RG) is a class of theoretical techniques used to explain the collective physics of interacting, many-body systems. It has been suggested that the RG formalism may be useful in finding and interpreting emergent low-dimensional structure in complex systems outside of the traditional physics context, such as in biology or computer science. In such contexts, one common dimensionality-reduction framework already in use is information bottleneck (IB), in which the goal is to compress an ‘input’ signal X while maximizing its mutual information with some stochastic ‘relevance’ variable Y. IB has been applied in the vertebrate and invertebrate processing systems to characterize optimal encoding of the future motion of the external world. Other recent work has shown that the RG scheme for the dimer model could be ‘discovered’ by a neural network attempting to solve an IB-like problem. This manuscript explores whether IB and any existing formulation of RG are formally equivalent. A class of soft-cutoff non-perturbative RG techniques are defined by families of non-deterministic coarsening maps, and hence can be formally mapped onto IB, and vice versa. For concreteness, this discussion is limited entirely to Gaussian statistics (GIB), for which IB has exact, closed-form solutions. Under this constraint, GIB has a semigroup structure, in which successive transformations remain IB-optimal. Further, the RG cutoff scheme associated with GIB can be identified. Our results suggest that IB can be used to impose a notion of ‘large scale’ structure, such as biological function, on an RG procedure.