Additive Renormalization Research Articles

Low-energy theorem revisited and OPE in massless QCD

We revisit a low-energy theorem (LET) of NSVZ type in SU(N) QCD with Nf massless quarks derived in [1] by implementing it in dimensional regularization. The LET relates n-point correlators in the l.h.s. to n + 1-point correlators with the extra insertion of TrF2 at zero momentum in the r.h.s. We demonstrate that, for 2-point correlators of an operator O in the l.h.s., the LET implies that, in general, the integrated 3-point correlator in the r.h.s. needs in perturbation theory an infinite additive renormalization in addition to the multiplicative one. We relate the above counterterm to a corresponding divergent contact term in a certain coefficient of the OPE of TrF2 with O in the momentum representation, thus extending to any operator O an independent argument that first appeared for O = TrF2 in [2]. Finally, we demonstrate that in the asymptotically free phase of QCD the aforementioned counterterm in the LET is actually finite nonperturbatively after resummation to all perturbative orders. We also briefly recall the implications of the LET in the gauge-invariant framework of dimensional regularization for the perturbative and nonperturbative renormalization in large-N QCD. The implications of the LET inside and above the conformal window of SU(N) QCD with Nf massless quarks will appear in a forthcoming paper.

Classical and quantum computing of shear viscosity for (2+1)D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for (2+1)-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density ηs is consistent with a well-known holographic result 14π at several temperatures on a 4×4 honeycomb lattice with the local electric representation truncated at jmax=12. We also find the ratio of the spectral function and frequency ρxy(ω)ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond jmax=12 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including jmax truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

Coexistence of Interacting Charge Density Waves in a Layered Semiconductor.

Coexisting orders are key features of strongly correlated materials and underlie many intriguing phenomena from unconventional superconductivity to topological orders. Here, we report the coexistence of two interacting charge-density-wave (CDW) orders in EuTe_{4}, a layered crystal that has drawn considerable attention owing to its anomalous thermal hysteresis and a semiconducting CDW state despite the absence of perfect Fermi surface nesting. By accessing unoccupied conduction bands with time- and angle-resolved photoemission measurements, we find that monolayers and bilayers of Te in the unit cell host different CDWs that are associated with distinct energy gaps. The two gaps display dichotomous evolutions following photoexcitation, where the larger bilayer CDW gap exhibits less renormalization and faster recovery. Surprisingly, the CDW in the Te monolayer displays an additional momentum-dependent gap renormalization that cannot be captured by density-functional theory calculations. This phenomenon is attributed to interlayer interactions between the two CDW orders, which account for the semiconducting nature of the equilibrium state. Our findings not only offer microscopic insights into the correlated ground state of EuTe_{4} but also provide a general nonequilibrium approach to understand coexisting, layer-dependent orders in a complex system.

An observation on Feynman diagrams with axial anomalous subgraphs in dimensional regularization with an anticommuting γ5

Through the calculation of the matrix element of the singlet axial-current operator between the vacuum and a pair of gluons in dimensional regularization with an anti-commuting γ5 defined in a Kreimer-scheme variant, we find that additional renormalization counter-terms proportional to the Chern-Simons current operator are needed starting from O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O} $$\\end{document}(αs2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\\alpha}_s^2 $$\\end{document}) in QCD. This is in contrast to the well-known purely multiplicative renormalization of the singlet axial-current operator defined with a non-anticommuting γ5. Consequently, without introducing compensation terms in the form of additional renormalization, the Adler-Bell-Jackiw anomaly equation does not hold automatically in the bare form in this kind of schemes. We determine the corresponding (gauge-dependent) coefficient to O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O} $$\\end{document}(αs3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\\alpha}_s^3 $$\\end{document}) in QCD, using a variant of the original Kreimer prescription which is implemented in our computation in terms of the standard cyclic trace together with a constructively-defined γ5. Owing to the factorized form of these divergences, intimately related to the axial anomaly, we further performed a check, using concrete examples, that with γ5 treated in this way, the axial-current operator needs no more additional renormalization in dimensional regularization but only for non-anomalous amplitudes in a perturbatively renormalizable theory. To be complete, we provide a few additional ingredients needed for a proposed extension of the algorithmic procedure formulated in the above analysis to potential applications to a renormalizable anomaly-free chiral gauge theory, i.e. the electroweak theory.

QCD transition at the physical point, and its scaling window from twisted mass Wilson fermions

We study the scaling properties of the finite temperature QCD phase transition, for light quark masses ranging from the heavy quark regime to their physical values. The lattice results are obtained in the fixed scale approach from simulations of Nf=2+1+1 flavours of Wilson fermions at maximal twist. We identify an order parameter free from the linear contributions in mass due to additive renormalization and regular terms in the Equation of State, which proves useful for the assessment of the hypothesized universal behaviour. We find compatibility with the 3D O(4) universality class for the physical pion mass and temperatures 120 MeV ≲T≲300 MeV. We discuss violation of scaling at larger masses and a possible cross-over to mean field behaviour. The chiral extrapolation T0=134−4+6 MeV of the pseudocritical temperature is robust against predictions of different universality classes and consistent with its estimate from the O(4) Equation of State for the physical pion mass.

Renormalization of the topological charge density in QCD with dimensional regularization

To all orders of perturbation theory, the renormalization of the topological charge density in dimensionally regularized QCD is shown to require no more than an additive renormalization proportional to the divergence of the flavour-singlet axial current. The proof is based on the standard BRS analysis of the QCD vertex functional in the background gauge and exploits the special algebraic properties of the charge density through the Stora–Zumino chain of descent equations.

Varieties and properties of central-branch Wilson fermions

We focus on the four-dimensional central-branch Wilson fermion, which makes good use of six species at the central branch of the Wilson Dirac spectrum and possesses the extra $U(1)_{\overline V}$ symmetry. With introducing new insights we discuss the prohibition of additive mass renormalization for all the six species, SSB of $U(1)_{\overline V}$ in strong-coupling QCD, the absence of the sign problem, and the usefulness for many-flavor QCD simulation. We then construct several varieties of the central-branch fermions and study their properties. In particular, we investigate the two-flavor version, where the Dirac spectrum has seven branches and two species live at the central branch. Although the hypercubic symmetry is broken, the other symmetries are the same as those of the original one. We study this setup in terms of lattice perturbation theory, strong-coupling QCD, the absence of sign problem, and the parameter tuning for Lorentz symmetry restoration. By comparing the properties of the original and two-flavor version, we find that the existence of hypercubic symmetry as well as $U(1)_{\overline V}$ is essential for the absence of additive mass renormalization of all the central-branch species. As the other two-flavor version, we investigate the central-branch staggered-Wilson fermion, which is obtained from the eight-flavor central-branch Wilson fermion via spin diagonalization. We argue that it is free from any additive mass renormalization and is regarded as a minimally doubled fermion with less symmetry breaking.

Local interstellar spectra and solar modulation of cosmic ray electrons and positrons

Low energy cosmic rays are modulated by the solar activity when they propagation in the heliosphere, leading to ambiguities in understanding their acceleration at sources and propagation in the Milky Way. By means of the precise measurements of the e−,e+,e−+e+, and e+/(e−+e+) spectra by AMS-02 near the Earth, as well as the very low energy measurements of the e−+e+ fluxes by Voyager-1 far away from the Sun, we derive the local interstellar spectra (LIS) of e− and e+ components individually. Our method is based on a non-parametric description of the LIS of e− and e+ and a force-field solar modulation model. We then obtain the evolution of the solar modulation parameters based on the derived LIS and the monthly fluxes of cosmic ray e− and e+ measured by AMS-02. To better fit the monthly data, additional renormalization factors for e− and e+ have been multiplied to the modulated fluxes. We find that the inferred solar modulation parameters of positrons are in good agreement with that of cosmic ray nuclei, and the time evolutions of the solar modulation parameters of electrons and positrons differ after the reversal of the heliosphere magnetic field polarity, which shows clearly the charge-sign dependent modulation effect.

Monte Carlo simulations of overlap Majorana fermions

Supersymmetric Yang-Mills (SYM) theories in four dimensions exhibit many interesting non-perturbative phenomena that can be studied by means of Monte Carlo lattice simulations. However, the lattice regularization breaks supersymmetry explicitly, and in general a fine tuning of a large number of parameters is required to correctly extrapolate the theory to the continuum limit. From this perspective, it is important to preserve on the lattice as many symmetries of the original continuum action as possible. Chiral symmetry for instance prevents an additive renormalization of the fermion mass. A (modified) version of chiral symmetry can be preserved exactly if the Dirac operator fulfills the Ginsparg-Wilson relation. In this contribution, we present an exploratory non-perturbative study of N=1 supersymmetric Yang-Mills theory using the overlap formalism to preserve chiral symmetry at non-zero lattice spacings. N=1 SYM is an ideal benchmark toward the extension of our studies to more complex supersymmetric theories, as the only parameter to be tuned is the gluino mass. Overlap fermions allow therefore to simulate the theory without fine-tuning. We compare our approach to previous investigations of the same theory, and we present clear evidences for gluino condensation.

Lattice gauge theory for the Haldane conjecture and central-branch Wilson fermion

Abstract We develop a $(1+1)$D lattice $U(1)$ gauge theory in order to define the two-flavor massless Schwinger model, and discuss its connection with the Haldane conjecture. We propose to use the central-branch Wilson fermion, which is defined by relating the mass, $m$, and the Wilson parameter, $r$, by $m+2r=0$. This setup gives two massless Dirac fermions in the continuum limit, and it turns out that no fine-tuning of $m$ is required because the extra $U(1)$ symmetry at the central branch, $U(1)_{\overline{V}}$, prohibits additive mass renormalization. Moreover, we show that the Dirac determinant is positive semi-definite and this formulation is free from the sign problem, so a Monte Carlo simulation of the path integral is possible. By identifying the symmetry at low energy, we show that this lattice model has a mixed ’t Hooft anomaly between $U(1)_{\overline{V}}$, lattice translation, and lattice rotation. We discuss its relation to the anomaly of half-integer anti-ferromagnetic spin chains, so our lattice gauge theory is suitable for numerical simulation of the Haldane conjecture. Furthermore, it gives a new and strict understanding on the parity-broken phase (Aoki phase) of the $2$D Wilson fermion.

Effects of K–K1 Transitions on Kaon-Field Renormalization in the NJL Model Framework

An additional kaon-field renormalization dictated by transitions between pseudoscalar and axial-vector fields is derived in the Nambu–Jona-Lasinio model framework. Possible mixings between the $${{K}_{{1A}}}$$ and $${{K}_{{1B}}}$$ axial-vector fields with quantum numbers of $${{J}^{{PC}}} = {{1}^{{ + + }}}$$ and $${{J}^{{PC}}} = {{1}^{{ + - }}}$$, respectively, are taken into account. We demonstrate that the $${{a}_{1}}(1260)$$ and $${{b}_{1}}(1235)$$ do not mix within the SU(2) $$ \times $$ SU(2) chiral scheme. However, a nonzero mixing between the $${{K}_{{1A}}}$$ and $${{K}_{{1B}}}$$ strange axial-vector mesons arises as soon as the strange-quark mass is not equal to those of the u and d quarks, whereby the chiral symmetry is essentially violated. As a result, the $${{K}_{1}}(1270)$$ and $${{K}_{1}}(1400)$$ mesons emerge as physical states. These are taken into account when formulating the additional renormalization of the kaon field.

Full O(a) improvement in electrostatic QCD

EQCD is a 3D bosonic theory containing SU(3) and an adjoint scalar, which efficiently describes the infrared, nonperturbative sector of hot QCD and which is highly amenable to lattice study. We improve the matching between lattice and continuum EQCD by determining the final unknown coefficient in the O(a) matching, an additive scalar mass renormalization. We do this numerically by using the symmetry-breaking phase transition point of the theory as a line of constant physics. This prepares the ground for a precision study of the transverse momentum diffusion coefficient C(qperp) within this theory. As a byproduct, we provide an updated version of the EQCD phase diagram.

Mass spectrum of 2-dimensional mathcal{N}=left(2,2right) super Yang-Mills theory on the lattice

In the present work we analyse mathcal{N}=left(2,2right) supersymmetric Yang-Mills (SYM) theory with gauge group SU(2) in two dimensions by means of lattice simulations. The theory arises as dimensional reduction of mathcal{N}=1 SYM theory in four dimensions. As in other gauge theories with extended supersymmetry, the classical scalar potential has flat directions which may destabilize numerical simulations. In addition, the fermion determinant need not be positive and this sign-problem may cause further problems in a stochastic treatment. We demonstrate that mathcal{N}=left(2,2right) super Yang-Mills theory has actually no sign problem and that the flat directions are lifted and thus stabilized by quantum corrections. Only the bare masses of the scalars experience a finite additive renormalization in this finite theory. On various lattices with different lattice constants we determine the scalar masses and hopping parameters for which the supersymmetry violating terms are minimal. By studying four Ward identities and by monitoring the π-mass we show that supersymmetry is indeed restored in the continuum limit. In the second part we calculate the masses of the low-lying bound states. We find that in the infinite-volume and supersymmetric continuum limit the Veneziano-Yankielowicz super-multiplet becomes massless and the Farrar-Gabadadze-Schwetz super-multiplet decouples from the theory. In addition, we estimate the masses of the excited mesons in the Veneziano-Yankielowicz multiplet. We observe that the gluino-glueballs have comparable masses to the excited mesons.

Meson masses in electromagnetic fields with Wilson fermions

We determine the light meson spectrum in QCD in the presence of background magnetic fields using quenched Wilson fermions. Our continuum extrapolated results indicate a monotonous reduction of the connected neutral pion mass as the magnetic field grows. The vector meson mass is found to remain nonzero, a finding relevant for the conjectured $\rho$-meson condensation at strong magnetic fields. The continuum extrapolation was facilitated by adding a novel magnetic field-dependent improvement term to the additive quark mass renormalization. Without this term, sizable lattice artifacts that would deceptively indicate an unphysical rise of the connected neutral pion mass for strong magnetic fields are present. We also investigate the impact of these lattice artifacts on further observables like magnetic polarizabilities and discuss the magnetic field-induced mixing between $\rho$-mesons and pions. We also derive Ward-Takashi identities for QCD+QED both in the continuum formulation and for (order $a$-improved) Wilson fermions.

Renormalization of Supersymmetric QCD on the Lattice

We perform a pilot study of the perturbative renormalization of a Supersymmetric gauge theory with matter fields on the lattice. As a specific example, we consider Supersymmetric N=1 QCD (SQCD). We study the self-energies of all particles which appear in this theory, as well as the renormalization of the coupling constant. To this end we compute, perturbatively to one-loop, the relevant two-point and three-point Green’s functions using both dimensional and lattice regularizations. Our lattice formulation involves theWilson discretization for the gluino and quark fields; for gluons we employ the Wilson gauge action; for scalar fields (squarks) we use naive discretization. The gauge group that we consider is SU(Nc), while the number of colors, Nc, the number of flavors, Nf, and the gauge parameter, α, are left unspecified. We obtain analytic expressions for the renormalization factors of the coupling constant (Zg) and of the quark (ZΨ), gluon (Zu), gluino (Zλ), squark (ZA±), and ghost (Zc) fields on the lattice. We also compute the critical values of the gluino, quark and squark masses. Finally, we address the mixing which occurs among squark degrees of freedom beyond tree level: we calculate the corresponding mixing matrix which is necessary in order to disentangle the components of the squark field via an additional finite renormalization.

Large thermopower from dressed quasiparticles in the layered cobaltates and rhodates

The origin of the large thermopower in ${\mathrm{Na}}_{x}{\mathrm{CoO}}_{2}$ is complicated by correlation phenomena. To disentangle the effects from multiple interactions, we use angle-resolved photoemission to study ${\mathrm{K}}_{x}{\mathrm{RhO}}_{2}$, an isostructural analogy of ${\mathrm{Na}}_{x}{\mathrm{CoO}}_{2}$ with large thermopower and weak electron correlation. Using the experimentally measured electronic structure, we demonstrate that the thermopower in ${\mathrm{K}}_{x}{\mathrm{RhO}}_{2}$ can be quantitatively explained within the quasiparticle framework after including an electron-phonon mass enhancement effect. Extending the analysis to the cobaltate, we find the doubling in thermopower is well accounted for by additional band renormalization from electron correlation. As such, the large thermopower emerges from the itinerant quasiparticles dressed by hierarchical electron-phonon and electron-electron interactions.

Supersymmetric QCD on the lattice: An exploratory study

We perform a pilot study of the perturbative renormalization of a Supersymmetric gauge theory with matter fields on the lattice. As a specific example, we consider Supersymmetric ${\cal N}{=}1$ QCD (SQCD). We study the self-energies of all particles which appear in this theory, as well as the renormalization of the coupling constant. To this end we compute, perturbatively to one-loop, the relevant two-point and three-point Green's functions using both dimensional and lattice regularizations. Our lattice formulation involves the Wilson discretization for the gluino and quark fields; for gluons we employ the Wilson gauge action; for scalar fields (squarks) we use na\"ive discretization. The gauge group that we consider is $SU(N_c)$, while the number of colors, $N_c$, the number of flavors, $N_f$, and the gauge parameter, $\alpha$, are left unspecified. We obtain analytic expressions for the renormalization factors of the coupling constant ($Z_g$) and of the quark ($Z_\psi$), gluon ($Z_u$), gluino ($Z_\lambda$), squark ($Z_{A_\pm}$), and ghost ($Z_c$) fields on the lattice. We also compute the critical values of the gluino, quark and squark masses. Finally, we address the mixing which occurs among squark degrees of freedom beyond tree level: we calculate the corresponding mixing matrix which is necessary in order to disentangle the components of the squark field via an additional finite renormalization.

Energy Level Alignment at Hybridized Organic–Metal Interfaces: The Role of Many-Electron Effects

Hybridized molecule/metal interfaces are ubiquitous in molecular and organic devices. The energy level alignment (ELA) of frontier molecular levels relative to the metal Fermi level (EF) is critical to the conductance and functionality of these devices. However, a clear understanding of the ELA that includes many-electron self-energy effects is lacking. Here, we investigate the many-electron effects on the ELA using state-of-the-art, benchmark GW calculations on prototypical chemisorbed molecules on Au(111), in eleven different geometries. The GW ELA is in good agreement with photoemission for monolayers of benzene-diamine on Au(111). We find that in addition to static image charge screening, the frontier levels in most of these geometries are renormalized by additional screening from substrate-mediated intermolecular Coulomb interactions. For weakly chemisorbed systems, such as amines and pyridines on Au, this additional level renormalization (~1.5 eV) comes solely from static screened exchange energy, allowing us to suggest computationally more tractable schemes to predict the ELA at such interfaces. However, for more strongly chemisorbed thiolate layers, dynamical effects are present. Our ab initio results constitute an important step towards the understanding and manipulation of functional molecular/organic systems for both fundamental studies and applications.

Lattice simulations withNf=2+1improved Wilson fermions at a fixed strange quark mass

The explicit breaking of chiral symmetry of the Wilson fermion action results in additive quark mass renormalization. Moreover, flavour singlet and non-singlet scalar currents acquire different renormalization constants with respect to continuum regularization schemes. This complicates keeping the renormalized strange quark mass fixed when varying the light quark mass in simulations with $N_f=2+1$ sea quark flavours. Here we present and validate our strategy within the CLS (Coordinated Lattice Simulations) effort to achieve this in simulations with non-perturbatively order-$a$ improved Wilson fermions. We also determine various combinations of renormalization constants and improvement coefficients.

One-shot calculation of temperature-dependent optical spectra and phonon-induced band-gap renormalization

Recently, Zacharias et al [Phys. Rev. Lett. 115, 177401 (2015)] developed a new ab initio theory of temperature-dependent optical absorption spectra and band gaps in semiconductors and insulators. In that work the zero-point renormalization and the temperature dependence were obtained by sampling the nuclear wavefunctions using a stochastic approach. In the present work, we show that the stochastic sampling can be replaced by fully deterministic supercell calculations based on a single optimal configuration of the atomic positions. We demonstrate that a single calculation is able to capture the temperature-dependent band gap renormalization including quantum nuclear effects in direct and indirect-gap semiconductors, as well as phonon-assisted optical absorption in indirect-gap semiconductors. In order to demonstrate this methodology we calculate from first principles the temperature-dependent optical absorption spectra and the renormalization of direct and indirect band gaps in Si, C, and GaAs, and we obtain good agreement with experiment and with previous calculations. In this work we also establish the formal connection between the Williams-Lax theory of optical transitions and the related theories of indirect absorption by Hall, Bardeen, and Blatt, and of temperature-dependent band structures by Allen and Heine. Furthermore, we identify an additional band gap renormalization that arises in the case of degenerate band extrema, and which has not been considered in previous work. The present methodology enables systematic ab initio calculations of optical absorption spectra at finite temperature, including both direct and indirect transitions. This feature will be useful for high-throughput calculations of optical properties at finite temperature, and for calculating temperature-dependent optical properties using high-level theories such as GW and Bethe-Salpeter approaches.