Strong Coupling Expansion Research Articles

Decomposing Imaginary-Time Feynman Diagrams Using Separable Basis Functions: Anderson Impurity Model Strong-Coupling Expansion

We present a deterministic algorithm for the efficient evaluation of imaginary-time diagrams based on the recently introduced discrete Lehmann representation (DLR) of imaginary-time Green’s functions. In addition to the efficient discretization of diagrammatic integrals afforded by its approximation properties, the DLR basis is separable in imaginary-time, allowing us to decompose diagrams into linear combinations of nested sequences of one-dimensional products and convolutions. Focusing on the strong-coupling bold-line expansion of generalized Anderson impurity models, we show that our strategy reduces the computational complexity of evaluating an Mth-order diagram at inverse temperature β and spectral width ωmax from O((βωmax)2M−1) for a direct quadrature to O(M(log(βωmax))M+1), with controllable high-order accuracy. We benchmark our algorithm using third-order expansions for multiband impurity problems with off-diagonal hybridization and spin-orbit coupling, presenting comparisons with exact diagonalization and quantum Monte Carlo approaches. In particular, we perform a self-consistent dynamical mean-field theory calculation for a three-band Hubbard model with strong spin-orbit coupling representing a minimal model of Ca2RuO4, demonstrating the promise of the method for modeling realistic strongly correlated multiband materials. For both strong and weak coupling expansions of low and intermediate order, in which diagrams can be enumerated, our method provides an efficient, straightforward, and robust blackbox evaluation procedure. In this sense, it fills a gap between diagrammatic approximations of the lowest order, which are simple and inexpensive but inaccurate, and those based on Monte Carlo sampling of high-order diagrams. Published by the American Physical Society 2024

Lattice Hamiltonian for adjoint QCD2

We introduce a Hamiltonian lattice model for the (1 + 1)-dimensional SU(Nc) gauge theory coupled to one adjoint Majorana fermion of mass m. The discretization of the continuum theory uses staggered Majorana fermions. We analyze the symmetries of the lattice model and find lattice analogs of the anomalies of the corresponding continuum theory. An important role is played by the lattice translation by one lattice site, which in the continuum limit involves a discrete axial transformation. On a lattice with periodic boundary conditions, the Hilbert space breaks up into sectors labeled by the Nc-ality p = 0, … Nc − 1. Our symmetry analysis implies various exact degeneracies in the spectrum of the lattice model. In particular, it shows that, for m = 0 and even Nc, the sectors p and p′ are degenerate if |p − p′| = Nc/2. In the Nc = 2 case, we explicitly construct the action of the Hamiltonian on a basis of gauge-invariant states, and we perform both a strong coupling expansion and exact diagonalization for lattices of up to 12 lattice sites. Upon extrapolation of these results, we find good agreement with the spectrum computed previously using discretized light-cone quantization. One of our new results is the first numerical calculation of the fermion bilinear condensate.

Tracy-Widom Distribution in Four-Dimensional Supersymmetric Yang-Mills Theories.

Various observables in different four-dimensional superconformal Yang-Mills theories can be computed exactly as Fredholm determinants of truncated Bessel operators. We exploit this relation to determine their dependence on the 't Hooft coupling constant. Unlike the weak coupling expansion, which has a finite radius of convergence, the strong coupling expansion is factorially divergent, necessitating the inclusion of nonperturbative, exponentially small corrections. We develop a method to systematically compute these corrections and discuss the resurgent properties of the resulting transseries.

The NNLO soft function for N-jettiness in hadronic collisions

We compute the N-jettiness soft function in hadronic collisions to next-to-next-to-leading order (NNLO) in the strong-coupling expansion. Our calculation is based on an extension of the SoftSERVE framework to soft functions that involve an arbitrary number of lightlike Wilson lines. We present numerical results for 1-jettiness and 2-jettiness, and illustrate that our formalism carries over to a generic number of jets by calculating a few benchmark points for 3-jettiness. We also perform a detailed analytic study of the asymptotic behaviour of the soft-function coefficients at the edges of phase space, where one of the jets becomes collinear to another jet or beam direction, and comment on previous calculations of the N-jettiness soft function.

Exact results for giant graviton four-point correlators

We study the four-point correlator O2O2DD\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\left\\langle {\\mathcal{O}}_2{\\mathcal{O}}_2\\mathcal{DD}\\right\\rangle $$\\end{document} in N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 super Yang-Mills theory (SYM) with SU(N) gauge group, where O2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\\mathcal{O}}_2 $$\\end{document} represents the superconformal primary operator with dimension two, while D\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{D} $$\\end{document} denotes a determinant operator of dimension N, which is holographically dual to a giant graviton D3-brane extending along S5. We analyse the integrated correlator associated with this observable, obtained after integrating out the spacetime dependence over a supersymmetric invariant measure. Similarly to other classes of integrated correlators in N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 SYM, this integrated correlator can be computed through supersymmetric localisation on the four-sphere. Employing matrix-model recursive techniques, we demonstrate that the integrated correlator can be reformulated as an infinite sum of protected three-point functions with known coefficients. This insight allows us to circumvent the complexity associated with the dimension-N determinant operator, significantly streamlining the large-N expansion of the integrated correlator. In the planar limit and beyond, we derive exact results for the integrated correlator valid for all values of the ’t Hooft coupling, and investigate the resurgent properties of their strong coupling expansion. Additionally, in the large-N expansion with fixed (complexified) Yang-Mills coupling, we deduce the SL(2, ℤ) completion of these results in terms of the non-holomorphic Eisenstein series. The proposed modular functions are confirmed by explicit instanton calculations from the matrix model, and agree with expectations from the holographic dual picture of known results in type IIB string theory.

The NNLO gluon beam function for jet-veto resummation

We compute the gluon beam function for jet-veto resummation to next-to-next-to-leading order (NNLO) in the strong-coupling expansion. Our calculation is based on an automated framework that was previously used for the computation of the respective quark beam function, and which we significantly extended for the present calculation. In particular, the perturbative matching kernels are directly calculated in momentum space, without the need to perform an additional Mellin transform. We present results for both gluon and quark-initiated processes, which we cross-checked with an independent semi-analytical method that exploits the similarity of the beam functions to the more familiar case of transverse-momentum resummation. Our computation is relevant for jet-veto resummations at NNLL′ accuracy.

Line operators, vortex statistics, and Higgs versus confinement dynamics

We study a 2+1D lattice gauge theory with fundamental representation scalar fields which has both Higgs and confining regimes with a spontaneously-broken U(1) 0-form symmetry. We show that the Higgs and confining regimes may be distinguished by a natural gauge invariant observable: the phase Ω of a correlation function of a vortex line operator linking with an electric Wilson line. We employ dualities and strong coupling expansions to analytically explore parameter regimes which were inaccessible in previous continuum calculations, and discuss possible implications for the phase diagram.

Integrated correlators at strong coupling in an orbifold of N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 SYM

We consider the 4dN\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 superconformal quiver gauge theory obtained by a ℤ2 orbifold of N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 super Yang-Mills (SYM). By exploiting supersymmetric localization, we study the integrated correlator of two Coulomb branch and two moment map operators and the integrated correlator of four moment map operators, determining exact expressions valid for any value of the ’t Hooft coupling in the planar limit. Additionally, for the second correlator, we obtain an exact expression also for the next-to-planar contribution. Then, we derive the leading terms of their strong-coupling expansions and outline the differences with respect to the N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 SYM theory.

Matrix product study of spin fractionalization in the one-dimensional Kondo insulator

The Kondo lattice is one of the classic examples of strongly correlated electronic systems. We conduct a controlled study of the Kondo lattice in one dimension, highlighting the role of excitations created by the composite fermion operator. Using time-dependent matrix product state methods, we compute various correlation functions and contrast them with both large-N mean-field theory and the strong-coupling expansion. We show that the composite fermion operator creates long-lived, charge-e and spin-1/2 excitations, which cover the low-lying single-particle excitation spectrum of the system. Furthermore, spin excitations can be thought to be composed of such fractionalized quasiparticles with a residual interaction which tend to disappear at weak Kondo coupling. Published by the American Physical Society 2024

Fermions, quantum gravity, and holography in two dimensions

We study a model comprising N flavors of Kähler Dirac fermion propagating on a triangulated two-dimensional disk which is constrained to have a negative average bulk curvature. Dirichlet boundary conditions are chosen for the fermions. Quantum fluctuations of the geometry are included by summing over all possible triangulations consistent with these constraints. We show in the limit N→∞ that the partition function is dominated by a regular triangulation of two-dimensional hyperbolic space. We use strong coupling expansions and Monte Carlo simulation to show that in this limit boundary correlators of the fermions have a power law dependence on boundary separation as one expects from holography. However, we argue that this behavior breaks down for any finite number of massive fields in the thermodynamic limit and quantum fluctuations of the bulk geometry drive the theory into a nonholographic phase. In contrast, for massless fermions, we find evidence that the boundary is conformal even for finite N. This is consistent with theoretical results in quantum Liouville theory. Published by the American Physical Society 2024

The planar limit of the N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 E-theory: numerical calculations and the large λ expansion

We study correlation functions of local operators and Wilson loop expectation values in the planar limit of a 4d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 superconformal SU(N) YM theory with hypermultiplets in the symmetric and antisymmetric representations of the gauge group. This so called E theory is closely related to N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 SYM and has a holographic description in terms of a ℤ2 orientifold of AdS5 × S5. Using recent matrix model results based on supersymmetric localization we develop efficient numerical methods to calculate two- and three-point functions of certain single trace operators as well as 1/2-BPS Wilson loop expectation values as a function of the ’t Hooft coupling λ. We use our numerical results to arrive at simple analytic expressions for these correlators valid up to sixth order in the λ−1/2 strong coupling expansion. These results provide explicit field theory predictions for the α′ corrections to the supergravity approximation of type IIB string theory on the AdS5 × S5/ℤ2 orientifold.

The ABJM Hagedorn Temperature from Integrability

We use the quantum spectral curve to compute the Hagedorn temperature for ABJM theory in terms of the interpolating function h(λ). At weak coupling we compute this temperature up to eight-loop order, showing that it matches the known tree-level and two-loop results. At strong coupling we compute the dependence numerically, showing that it is consistent with expectations from supergravity and the plane-wave limit for the four leading terms in the strong coupling expansion, up to an overall shift of the zero-point energy for type IIA string theory on AdS4 × ℂP3. We conjecture an analytic form for this shift to leading order that is consistent with our numerical results.

Higher order corrections to the Hagedorn temperature at strong coupling

We propose a general formula for higher order corrections to the value of the Hagedorn temperature of a class of holographic confining gauge theories in the strong coupling expansion. Inspired by recent proposals in the literature, the formula combines the sigma-model string expansion with an effective approach. In particular, it includes the sigma-model contributions to the Hagedorn temperature at next-to-next-to leading order, which are computed in full generality. For mathcal{N} = 4 SYM on S3 our result agrees with numerical estimates with excellent precision. We use the general formula to predict the value of the Hagedorn temperature for ABJM on S2 and for the dual of purely RR global AdS3.

Defect correlators in a mathcal{N} = 2 SCFT at strong coupling

We study the correlation function between one single-trace scalar operator and a circular Wilson loop in the 4d mathcal{N} = 2 superconformal field theory with gauge group SU(N) and matter transforming in the symmetric and anti-symmetric representations. By exploiting supersymmetric localization, we resum the perturbative expansion of this correlator in the large-N ’t Hooft limit. Furthermore, using both analytical and numerical techniques, we provide a prediction for the leading term of its strong coupling expansion and we compare this prediction to numerical Padé resummations of the perturbative series.

High-order parametrization of the hypergeometric-Meijer approximants

In this work, we introduce an extension to the hypergeometric algorithm we developed before for the resummation of divergent series. The extension overcome the time-consuming problem we face in the parametrization process of the hypergeometric approximants. In the previous version, the parametrization process results in a set of non-linear equations which is hard to solve (might be impossible) for relatively high orders (greater than six) using normal PCs with Mathematica software for instance. To solve this problem, we formulate an equivalent (order by order) linear set of equations which is easy to solve in an appropriate time using normal PCs. We also show that such extension of the hypergeometric resummation algorithm is able to employ non-perturbative information like strong-coupling and large-order asymptotic data which are always used to accelerate the convergence. We applied the algorithm for different orders (up to O(29)) of the ground state energy of the x4 anharmonic oscillator with and without the non-perturbative information. We also considered the available 20 orders for the ground state energy of the PT−symmetric ix3 anharmonic oscillator as well as the given 20 orders of its strong-coupling expansion or equivalently the Yang–Lee model. For high order weak-coupling parametrization, accurate results have been obtained for the ground state energy and the non-perturbative parameters describing strong-coupling and large-order asymptotic behaviors. The employment of the non-perturbative data accelerated the convergence very clearly. The High temperature expansion for the susceptibility within the SQ lattice has been also considered and led to accurate prediction for the critical exponent and critical temperature.

Tensor renormalization group approach to ( 1+1 )-dimensional SU(2) principal chiral model at finite density

We apply the tensor renormalization group method to the (1+1)-dimensional SU(2) principal chiral model at finite chemical potential with the use of the Gauss-Legendre quadrature to discretize the SU(2) Lie group. The internal energy at vanishing chemical potential $\mu=0$ shows good consistency with the prediction of the strong and weak coupling expansions. This indicates an effectiveness of the Gauss-Legendre quadrature for the partitioning of the SU(2) Lie group. In the finite density region with $\mu\ne 0$ at the strong coupling we observe the Silver-Blaze phenomenon for the number density.

Perturbative boundaries of quantum advantage: Real-time evolution for digitized λϕ4 lattice models

The real time evolution of quantum field theory models can be calculated order by order in perturbation theory. For $\ensuremath{\lambda}{\ensuremath{\phi}}^{4}$ models, the perturbative series have a zero radius of convergence which in part motivated the design of digitized versions suitable for quantum computing. In agreement with general arguments suggesting that a large field cutoff modifies Dyson's reasoning and improves convergence properties, we show that the harmonic digitizations of $\ensuremath{\lambda}{\ensuremath{\phi}}^{4}$ lattice field theories lead to weak coupling expansions with a finite radius of convergence. Similar convergence properties are found for strong coupling expansions. We compare the resources needed to calculate the real-time evolution of the digitized models with perturbative expansions to those needed to do so with universal quantum computers. Unless new approximate methods can be designed to calculate long perturbative series for large systems efficiently, it appears that the use of universal quantum computers with digitizations involving a few qubits per site has the potential for more efficient calculations of the real-time evolution for large systems at intermediate coupling.

Ferromagnetism in the SU( N ) Kondo lattice model: SU( N ) double exchange and supersymmetry

We study the ground-state properties of the SU(N)-generalization of the Kondo-lattice model in one dimension when the Kondo coupling J_K (both ferromagnetic and antiferromagnetic) is sufficiently strong. Both cases can be realized using alkaline-earth-like cold gases in optical lattices. Specifically, we first carry out the strong-coupling expansion and identify two insulating phases (one of which is the SU(N)-analogue of the well-known gapped Kondo singlet phase). We then rigorously establish that the ground state in the low-density (for J_K<0) or the high-density (for J_K>0) region is ferromagnetic. The results are accounted for by generalizing the double-exchange mechanism to SU(N) "spins". Possible realizations of Bose-Fermi supersymmetry SU(N|1) in the (generalized) SU(N) Kondo-lattice model are discussed as well.

Strong coupling expansions in mathcal{N} = 2 quiver gauge theories

We study the 3-point functions of gauge-invariant scalar operators in four dimensional mathcal{N} = 2 superconformal quiver theories using supersymmetric localization in the planar limit of a large number of colors. By exploiting a web of nontrivial relations, we show that the 3-point functions can be expressed in terms of the 2-point functions through exact Ward-like identities that are valid for all values of the coupling constant. In this way, using recent results about the 2-point functions, we are able to obtain the asymptotic strong-coupling expansion of the 3-point functions and of the corresponding structure constants in the planar limit. Our results extend to sub-leading orders what has been recently found at leading order, where a precise match with calculations within the AdS/CFT correspondence at the supergravity level is possible. Therefore, our findings can be interpreted also as a prediction for the sub-leading string corrections to these holographic calculations.

Vertex-Based Diagrammatic Treatment of Light-Matter-Coupled Systems.

We propose a diagrammatic MonteCarlo approach for quantum impurity models, which can be regarded as a generalization of the strong-coupling expansion for fermionic impurity models. The algorithm is based on a self-consistently computed three-point vertex and a stochastically sampled four-point vertex, and it allows one to obtain numerically exact results in a wide parameter regime. The performance of the algorithm is demonstrated with applications to a spin-boson model representing an emitter in a waveguide. As a function of the coupling strength, the spin exhibits a delocalization-localization crossover at low temperatures, signaling a qualitative change in the real-time relaxation. In certain parameter regimes, the response functions of the emitter coupled to the electromagnetic continuum can be described by an effective Rabi model with appropriately defined parameters. We also discuss the spatial distribution of the photon density around the emitter.