Schwinger-Dyson Research Articles

Tensor approximation of functional differential equations.

Functional differential equations(FDEs) play a fundamental role in many areas of mathematical physics, including fluid dynamics (Hopf characteristic functional equation), quantum field theory (Schwinger-Dyson equations), and statistical physics. Despite their significance, computing solutions to FDEs remains a longstanding challenge in mathematical physics. In this paper we address this challenge by introducing approximation theory and high-performance computational algorithms designed for solving FDEs on tensor manifolds. Our approach involves approximating FDEs using high-dimensional partial differential equations(PDEs), and then solving such high-dimensional PDEs on a low-rank tensor manifold leveraging high-performance (parallel) tensor algorithms. The effectiveness of the proposed approach is demonstrated through its application to the Burgers-Hopf FDE, which governs the characteristic functional of the stochastic solution to the Burgers equationevolving from a random initial state.

QCD equation of state and thermodynamic observables from computationally minimal Dyson-Schwinger equations

We study the QCD equation of state and other thermodynamic observables including the isentropic trajectories and the speed of sound. These observables are of eminent importance for the understanding of experimental results in heavy ion collisions and also provide a QCD input for studies of the timeline of heavy-ion-collisions with hydrodynamical simulations. They can be derived from the quark propagator whose gap equation is solved within a minimal approximation to the Dyson-Schwinger equations (DSEs) of QCD at finite temperature and density. This minimal approximation aims at a combination of computational efficiency and simplification of the truncation scheme while maintaining quantitative precision. This minimal DSE scheme is confronted and benchmarked with results for correlation functions and observables from lattice QCD at vanishing density and quantitative functional approaches at finite density. Published by the American Physical Society 2024

Correlation Functions Involving Dirac Fields from Homotopy Algebras II: The Interacting Theory

Abstract We extend the formula for correlation functions of free scalar field theories and Dirac field theories in terms of quantum $A_{\infty }$ algebras presented in arXiv:2305.11634 to general scalar-Dirac systems. We obtain the result that the same formula as in the previous paper holds in this case. We show that correlation functions from our formula satisfy the Schwinger–Dyson equations. We therefore confirm that correlation functions from our formula express correlation functions from the ordinary approach of quantum field theory.

Four-gluon vertex in collinear kinematics

To date, the four-gluon vertex is the least explored component of the QCD Lagrangian, mainly due to the vast proliferation of Lorentz and color structures required for its description. In this work we present a nonperturbative study of this vertex, based on the one-loop dressed Schwinger–Dyson equation obtained from the 4PI effective action. A vast simplification is brought about by resorting to “collinear” kinematics, where all momenta are parallel to each other, and by appealing to the charge conjugation symmetry in order to eliminate certain color structures. Out of the fifteen form factors that comprise the transversely-projected version of this vertex, two are singled out and studied in detail; the one associated with the classical tensorial structure is moderately suppressed in the infrared regime, while the other diverges logarithmically at the origin. Quite interestingly, both form factors display the property known as “planar degeneracy” at a rather high level of accuracy. With these results we construct an effective charge that quantifies the strength of the four-gluon interaction, and compare it with other vertex-derived charges from the gauge sector of QCD.

Universal Phenomenology at Critical Exceptional Points of Nonequilibrium O(N) Models

In thermal equilibrium the dynamics of phase transitions is largely controlled by fluctuation-dissipation relations: On the one hand, friction suppresses fluctuations, while on the other hand, the thermal noise is proportional to friction constants. Out of equilibrium, this balance dissolves and one can have situations where friction vanishes due to antidamping in the presence of a finite noise level. We study a wide class of O(N) field theories where this situation is realized at a phase transition, which we identify as a critical exceptional point. In the ordered phase, antidamping induces a continuous limit cycle rotation of the order parameter with an enhanced number of 2N−3 Goldstone modes. Close to the critical exceptional point, however, fluctuations diverge so strongly due to the suppression of friction that in dimensions d<4 they universally either destroy a preexisting static order or give rise to a fluctuation-induced first-order transition. This is demonstrated within a full resummation of loop corrections via Dyson-Schwinger equations for N=2, and a generalization for arbitrary N, which can be solved in the long wavelength limit. We show that in order to realize this physics it is not necessary to drive a system far out of equilibrium: Using the peculiar protection of Goldstone modes, the transition from an xy magnet to a ferrimagnet is governed by an exceptional critical point once weakly perturbed away from thermal equilibrium. Published by the American Physical Society 2024

Schwinger–Dyson Equation on the Complex Plane: A Four-Fermion Interaction Model at Finite Temperature

Abstract We extend the Schwinger–Dyson equation (SDE) on the complex plane, which was treated in our previous research, to finite temperature. As a simple example, we solve the SDE for a model with four-fermion interactions in (1+1) space-time dimensions in the strong coupling region. We investigate the properties of the effective mass and energy for the fermions, especially near the phase transition temperature.

Strebel Differentials and String Field Theory

Abstract A closed string worldsheet of genus g with n punctures can be presented as a contact interaction in which n semi-infinite cylinders are glued together in a specific way via the Strebel differential on it, if $n\ge 1,\ 2g-2+n\gt 0$. We construct a string field theory of closed strings such that all the Feynman diagrams are represented by such contact interactions. In order to do so, we define off-shell amplitudes in the underlying string theory using the combinatorial Fenchel–Nielsen coordinates to describe the moduli space and derive a recursion relation satisfied by them. Utilizing the Fokker–Planck formalism, we construct a string field theory from which the recursion relation can be deduced through the Schwinger–Dyson equation. The Fokker–Planck Hamiltonian consists of kinetic terms and three-string interaction terms.

Effective potential of composite operators in the first order region of the QCD phase transition

We propose a method to determine the effective potential of QCD from the gap equation by introducing the homotopy method between the solutions of the equation of motion. Via this method, the effective potential beyond the bare vertex approximation is obtained, which then generalizes the Cornwall-Jackiw-Tomboulis effective potential for the bilocal composite operators. Moreover, the extended effective potential is set to be a function of self-energy instead of the propagator, which is the key point for the potential to be bounded from below. We then investigate the extended effective potential in the cases of phase transition of the QCD vacuum with a small current quark mass, the first order phase transition of QCD at finite temperature, and high baryon chemical potential. In the former case, the effective potential shows as an inflection point at the critical mass where the multiple solutions of the Dyson-Schwinger equation (DSE) vanishes, which is consistent with that obtained by solving the DSE directly. For the latter case, the in-medium properties, such as the latent heat and the difference of trace anomaly, of QCD is obtained. Published by the American Physical Society 2024

BPZ equations for higher degenerate fields and nonperturbative Dyson-Schwinger equations

In the two-dimensional Liouville conformal field theory, correlation functions involving a degenerate field satisfy partial differential equations due to the decoupling of the null descendant field. On the other hand, the instanton partition function of a four-dimensional N=2 supersymmetric theory in the Ω-background at a special point of the parameter space also satisfies a partial differential equation resulting from the constraints of the gauge field configurations. This partial differential equation can be proved using the nonperturbative Dyson-Schwinger equations. We show for the next-to-simplest case that the partial differential equations obtained from two different perspectives can be identified, thereby confirming an assertion of the Bogomol’nyi-Prasad-Sommerfield/conformal field theory correspondence. Published by the American Physical Society 2024

First radial excitations of mesons and diquarks in a contact interaction

We present a calculation for the masses of the first radially excited states of 40 mesons and diquarks made up of u, d, s, c, and b quarks, including states that contain one or both heavy quarks. To this end, we employ a combined analysis of the Bethe-Salpeter and Schwinger-Dyson equations within a self-consistent and symmetry-preserving vector-vector contact interaction. The same set of parameters describes ground and excited states of mesons and their diquark partners. The wave function of the first radial excitation contains a zero whose location is correlated with an additional parameter dF which is a function of dressed quark masses. Our results satisfy the equal spacing rules given by the Gell-Mann Okubo mass relations. Wherever possible, we make comparisons of our findings with known experimental observations as well as theoretical predictions of several other models and approaches including lattice quantum chromodynamics, finding a very good agreement. We report predictions for a multitude of radial excitations not yet observed in experiments. Published by the American Physical Society 2024

Ghost-gluon vertex in the presence of the Gribov horizon: General kinematics

Correlation functions are important probes for the behavior of quantum field theories. Already at tree-level, the refined Gribov-Zwanziger (RGZ) effective action for Yang-Mills theories provides a good approximation for the gluon propagator, as compared to that calculated by nonperturbative methods such as lattice field theory and Dyson-Schwinger equations. However, the study of higher correlation functions of the RGZ theory is still at its beginning. In this work we evaluate the ghost-gluon vertex function in Landau gauge at one-loop level, in d=4 space-time dimensions for the gauge groups SU(2) and SU(3). More precisely, we extend the analysis conducted in [1] for the soft-gluon limit to an arbitrary kinematic configuration. We introduce renormalization group effects by means of a toy model for the running coupling and investigate the impact of such a model in the ultraviolet tails of our results. We find that RGZ results match fairly closely those from lattice simulations, Schwinger-Dyson equations and the Curci-Ferrari model for three different kinematic configurations. This is compatible with RGZ being a feasible theory for the strong interaction in the infrared regime. Published by the American Physical Society 2024

Bound states from the spectral Bethe-Salpeter equation

We compute the bound state properties of three-dimensional scalar ϕ4 theory in the broken phase. To this end, we extend the recently developed technique of spectral Dyson-Schwinger equations to solve the Bethe-Salpeter equation and determine the bound state spectrum. We employ consistent truncations for the two-, three- and four-point functions of the theory that recover the scaling properties in the infinite coupling limit. Our result for the mass of the lowest-lying bound state in this limit agrees very well with lattice determinations. Published by the American Physical Society 2024

Coloured combinatorial maps and quartic bi-tracial 2-matrix ensembles from noncommutative geometry

We compute the first twenty moments of three convergent quartic bi-tracial 2-matrix ensembles in the large N limit. These ensembles are toy models for Euclidean quantum gravity originally proposed by John Barrett and collaborators. A perturbative solution is found for the first twenty moments using the Schwinger-Dyson equations and properties of certain bi-colored unstable maps associated to the model. We then apply a result of Guionnet et al. to show that the perturbative and convergent solution coincide for a small neighbourhood of the coupling constants. For each model we compute an explicit expression for the free energy, critical points, and critical exponents in the large N limit. In particular, the string susceptibility is found to be γ = 1/2, hinting that the associated universality class of the model is the continuous random tree.

Correlation functions of scalar field theories from homotopy algebras

We present expressions for correlation functions of scalar field theories in perturbation theory using quantum A∞ algebras. Our expressions are highly explicit and can be used for theories both in Euclidean space and in Minkowski space including quantum mechanics. Correlation functions at a given order of perturbation theory can be calculated algebraically without using canonical quantization or the path integral, and we demonstrate it explicitly for φ3 theory. We show that the Schwinger-Dyson equations are satisfied as an immediate consequence of the form of the expressions based on quantum A∞ algebras.

Gravitational bound waveforms from amplitudes

With the aim of computing bound waveforms from scattering amplitudes, we explore gravitational two-body dynamics using the Schwinger-Dyson equations and Bethe-Salpeter recursion. We show that the tree-level scattering waveform admits a natural analytic continuation, in rapidity, to the bound waveform, which we confirm from an independent calculation, in the Post-Newtonian expansion, of the time-domain multipoles at large eccentricity. We demonstrate consistency of this scattering-to-bound map with the Damour-Deruelle prescription for orbital elements in the quasi-Keplerian parametrization (which enters into the evaluation of the multipoles) and with the analytic continuation, in the binding energy, of radiated energy and angular momentum at 3PM.

Dependence of the Landau gauge ghost-gluon-vertex on the number of flavors

The gauge-boson, ghost and fermion propagators as well as the gauge-boson-ghost vertex function are studied for SU(N), Sp(2N), and SO(N) gauge groups. We solve a set of coupled Dyson-Schwinger equations in Landau gauge for a variable, fractional number Nf of massless fermions in the fundamental representation. For large Nf we find a phase transition from a chirally broken into a chirally symmetric phase that is consistent with the behavior expected inside the conformal window. Even in the presence of fermions the gauge-boson-ghost-vertex dressing remains small. In the conformal window this vertex shows the expected power law behavior. It does not assume its tree-level value in the far infrared, but the respective dressing function is a constant greater than 1. Published by the American Physical Society 2024

Renormalisable Non-Local Quark–Gluon Interaction: Mass Gap, Chiral Symmetry Breaking and Scale Invariance

We derive a Nambu–Jona-Lasinio (NJL) model from a non-local gauge theory and show that it has confining properties at low energies. In particular, we present an extended approach to non-local QCD and a complete revision of the technique of Bender, Milton and Savage applied to non-local theories, providing a set of Dyson–Schwinger equations in differential form. In the local case, we obtain closed-form solutions in the simplest case of the scalar field and extend it to the Yang–Mills field. In general, for non-local theories, we use a perturbative technique and a Fourier series and show how higher-order harmonics are heavily damped due to the presence of the non-local factor. The spectrum of the theory is analysed for the non-local Yang–Mills sector and found to be in agreement with the local results on the lattice in the limit of the non-locality mass parameter running to infinity. In the non-local case, we confine ourselves to a non-locality mass that is sufficiently large compared to the mass scale arising from the integration of the Dyson–Schwinger equations. Such a choice results in good agreement, in the proper limit, with the spectrum of the local theory. We derive a gap equation for the fermions in the theory that gives some indication of quark confinement in the non-local NJL case as well. Confinement seems to be a rather ubiquitous effect that removes some degrees of freedom in the original action, favouring the appearance of new observable states, as seen, e.g., for quantum chromodynamics at lower energies.

Random finite noncommutative geometries and topological recursion

In this paper, we investigate a model for quantum gravity on finite noncommutative spaces using the theory of blobbed topological recursion. The model is based on a particular class of random finite real spectral triples {(\mathcal{A}, \mathcal{H}, D, \gamma, J)} , called random matrix geometries of type {(1,0)} , with a fixed fermion space {(\mathcal{A}, \mathcal{H}, \gamma, J)} and a distribution of the form {e^{- \mathcal{S} (D)}\, {\mathrm{d}} D} over the moduli space of Dirac operators. The action functional {\mathcal{S} (D)} is considered to be a sum of terms of the form {\prod_{i=1}^s \mathrm{Tr} ({D^{n_i}})} for arbitrary {s \geqslant 1} . The Schwinger–Dyson equations satisfied by the connected correlators {W_n} of the corresponding multi-trace formal 1-Hermitian matrix model are derived by a differential geometric approach. It is shown that the coefficients {W_{g,n}} of the large N expansion of {W_n} ’s enumerate discrete surfaces, called stuffed maps, whose building blocks are of particular topologies. The spectral curve {( {\Sigma, \omega_{0,1}, \omega_{0,2}} )} of the model is investigated in detail. In particular, we derive an explicit expression for the fundamental symmetric bidifferential {\omega_{0,2}} in terms of the formal parameters of the model.

Systematic Incorporation of Ionic Hard-Core Size into the Debye-Hückel Theory via the Cumulant Expansion of the Schwinger-Dyson Equations.

The Debye-Hückel (DH) formalism of bulk electrolytes equivalent to the Gaussian-level closure of the electrostatic Schwinger-Dyson identities without the interionic hard-core (HC) coupling is extended via the cumulant treatment of these equations augmented by HC interactions. By comparing the monovalent ion activity and pressure predictions of our cumulant-corrected DH (CCDH) theory with hypernetted-chain results and Monte Carlo simulations from the literature, we show that this rectification extends the accuracy of the DH formalism from submolar to molar salt concentrations. In the case of internal energies or the general case of divalent electrolytes mainly governed by charge correlations, the improved accuracy of the CCDH theory is limited to submolar ion concentrations. Comparison with experimental data from the literature shows that, via the adjustment of the hydrated ion radii, CCDH formalism can equally reproduce the nonuniform effect of salt increment on the ionic activity coefficients up to molar concentrations. The inequality satisfied by these HC sizes coincides with the cationic branch of the Hofmeister series.

QCD anomalies in electromagnetic processes: A solution to the γ→3π puzzle

In this work, the γ→3π form factor is calculated within the Dyson-Schwinger equations framework using a contact interaction model within the so-called modified rainbow ladder truncation. The present calculation takes into account the pseudovector component in the pion Bethe-Salpeter amplitude (BSA) and π−π scattering effects, producing a γ→3π anomaly that is 1+6Rπ2 larger than the low energy prediction. Here Rπ is the relative ratio of the pseudovector and pseudoscalar components in the pion BSA; with our parameters’ input, this correction raises the γ→3π anomaly by around 10%. The main outcome of this work is the unveiling of the origin of such a correction, which could be a possible explanation of the discrepancy between the existing experimental data and the low energy prediction. Moreover, it is highlighted how the magnitude of the anomaly is affected in effective theories that require an irremovable ultraviolet cutoff. We find that for both the anomalous processes π→2γ and γ→3π, the missing contribution to the anomaly can be compensated by the additional structures related with the quark anomalous magnetic moment. Published by the American Physical Society 2024