Abstract: We unify two widely different approaches to understanding the infrared behavior of quantum chromodynamics (QCD), one essentially phenomenological, based on data, and the other computational, realized via quantum field equations in the continuum theory. Using the latter, we explain and calculate a process-independent running coupling for QCD, a new type of effective charge that is an analogue of the Gell-Mann-Low effective coupling in quantum electrodynamics. The result is almost identical to the process-dependent effective charge defined via the Bjorken sum rule, which provides one of the most basic constraints on our knowledge of nucleon spin structure. This reveals the Bjorken sum to be a near direct means by which to gain empirical insight into QCD's Gell-Mann-Low effective charge.

Abstract: The full quark-gluon vertex is a crucial ingredient for the dynamical generation of a constituent quark mass from the standard quark gap equation, and its nontransverse part may be determined exactly from the nonlinear Slav nov-Taylor identity that it satisfies. The resulting expression involves not only the quark propagator, but also the ghost dressing function and the quark-ghost kernel, and constitutes the non-abelian extension of the so-called “Ball-Chiu vertex,” known from QED. In the present work we carry out a detailed study of the impact of this vertex on the gap equation and the quark masses generated from it, putting particular emphasis on the contributions directly related with the ghost sector of the theory, and especially the quark-ghost kernel. In particular, we set up and solve the coupled system of six equations that determine the four form factors of the latter kernel and the two typical Dirac structures composing the quark propagator. Due to the incomplete implementation of the multiplicative renormalizability at the level of the gap equation, the correct anomalous dimension of the quark mass is recovered through the inclusion of a certain function, whose ultraviolet behavior is fixed, but its infrared completion is unknown; three particular Ansatze for this function are considered, and their effect on the quark mass and the pion decay constant is explored. The main results of this study indicate that the numerical impact of the quark-ghost kernel is considerable; the transition from a tree-level kernel to the one computed hem leads to a 20% increase in the value of the quark mass at the origin. Particularly interesting is the contribution of the fourth Ball-Chiu form factor, which, contrary to the Abelian case, is nonvanishing, and accounts for 10% of the total constituent quark mass.

Abstract: We present the detailed derivation of the longitudinal part of the three-gluon vertex from the Slavnov-Taylor identities that it satisfies, by means of a nonperturbative implementation of the Ball-Chiu construction; the latter, in its original form, involves the inverse gluon propagator, the ghost dressing function, and certain form factors of the ghost-gluon kernel. The main conceptual subtlety that renders this endeavor nontrivial is the infrared finiteness of the gluon propagator, and the resulting need to separate the vertex into two pieces, one that is intimately connected with the emergence of a gluonic mass scale, and one that satisfies the original set of Slavnov-Taylor identities, but with the inverse gluon propagator replaced by its “kinetic” term. The longitudinal form factors obtained by this construction are presented for arbitrary Euclidean momenta, as well as special kinematic configurations, parametrized by a single momentum. A particularly preeminent feature of the components comprising the tree-level vertex is their considerable suppression for momenta below 1 GeV, and the appearance of the characteristic “zero-crossing” in the vicinity of 100-200 MeV. Special combinations of the form factors derived with this method are compared with the results of recent large-volume lattice simulations, and are found to capture faithfully the rather complicated curves formed by the data. A similar comparison with results obtained from Schwinger-Dyson equations reveals a fair overall agreement, but with appreciable differences at intermediate energies. A variety of issues related to the distribution of the pole terms responsible for the gluon mass generation are discussed in detail, and their impact on the structure of the transverse parts is elucidated. In addition, a brief account of several theoretical and phenomenological possibilities involving these newly acquired results is presented.

Abstract: We present a novel analysis of the gluon gap equation, where its full nonlinear structure is duly taken into account. In particular, while in previous treatments the linearization of this homogeneous integral equation introduced an indeterminacy in the scale of the corresponding mass, the current approach determines it uniquely, once the value of the gauge coupling at a given renormalization point is used as input. A crucial ingredient for this construction is the “kinetic term” of the gluon propagator, whose form is not obtained from the complicated equation governing its evolution, but is rather approximated by suitable initial Ansatze, which are subsequently improved by means of a systematic iterative procedure. The multiplicative renormalization of the central equation is carried out following an approximate method, which is extensively employed in the studies of the standard quark gap equation. This approach amounts to the effective substitution of the vertex renormalization constants by kinematically simplified form factors of the three- and four-gluon vertices. The resulting numerical interplay, exemplified by the infrared suppression of the three-gluon vertex and the mild enhancement of the four-gluon vertex, is instrumental for obtaining positive-definite and monotonically decreasing running gluon masses. The resulting gluon propagators, put together from the gluon masses and kinetic terms obtained with this method, match rather accurately the data obtained from large-volume lattice simulations.

Abstract: We present a comprehensive study of the quark sector of 2 + 1 flavor QCD, based on a self-consistent treatment of the coupled system of Schwinger-Dyson equations for the quark propagator and the full quark-gluon vertex in the one-loop dressed approximation. The individual form factors of the quark-gluon vertex are expressed in a special tensor basis obtained from a set of gauge-invariant operators. The sole external ingredient used as input to our equations is the Landau gauge gluon propagator with 2 + 1 dynamical quark flavors, obtained from studies with Schwinger-Dyson equations, the functional renormalization group approach, and large volume lattice simulations. The appropriate renormalization procedure required in order to self-consistently accommodate external inputs stemming from other functional approaches or the lattice is discussed in detail, and the value of the gauge coupling is accurately determined at two vastly separated renormalization group scales. Our analysis establishes a clear hierarchy among the vertex form factors. We identify only three dominant ones, in agreement with previous results. The components of the quark propagator obtained from our approach are in excellent agreement with the results from Schwinger-Dyson equations, the functional renormalization group, and lattice QCD simulation, a simple benchmark observable being the chiral condensate in the chiral limit, which is computed as (245 MeV)(3). The present approach has a wide range of applications, including the self-consistent computation of bound-state properties and finite temperature and density physics, which are briefly discussed.