--- title: Graphical Functions in Quantum Field Theory url: https://www.emergentmind.com/topics/graphical-functions type: topic --- # Graphical Functions in Quantum Field Theory Searching arXiv for the core literature on graphical functions in quantum field theory to ground the article in the cited papers. Searching for foundational papers on scalar graphical functions, parametric representations, even-dimensional theory, and spin generalizations. Graphical functions are special massless three-point position-space Feynman integrals associated with a connected graph carrying three distinguished external vertices, usually labeled \(0\), \(1\), and \(z\). After fixing the external configuration by a Euclidean motion so that \(x_0=0\), \(x_1=1\), and \(x_z=z\in\mathbb C\setminus\{0,1\}\), the resulting integral becomes a function on the punctured complex plane. In even dimensions \(d=2\lambda+2>2\), and under the usual infrared and ultraviolet power-counting conditions, scalar graphical functions define single-valued real-analytic functions on \(\mathbb C\setminus\{0,1\}\); in the scalar case they are moreover positive functions \(\mathbb C\setminus\{0,1\}\to\mathbb R_+\) [1509.07296][2105.05015]. Their theory sits at the intersection of position-space quantum field theory, Symanzik-polynomial techniques, Möbius symmetry, single-valued multiple polylogarithms, and high-loop period calculations, and it has recently been generalized from scalar theories to theories with spin by incorporating numerator structures and spin-dependent differential operators [1302.6445][2504.05850]. ## 1. Definition and normalization conventions Let \(G\) be a connected graph with three distinguished external vertices and internal vertex set \(V_{\mathrm{int}}\). In the scalar setting, each edge \(e=\{u,v\}\) carries a weight \(\nu_e>0\), and one sets \[ Q_e=|x_u-x_v|^2. \] In Euclidean dimension \(d=2\lambda+2>2\), the position-space Feynman integral is \[ f_G^{(\lambda)}(z) = \prod_{v\in V_{\mathrm{int}}}\int_{\mathbb R^d} d^d x_v\,\pi^{-d/2} \prod_{e\in E_G} Q_e^{-\lambda \nu_e}, \] with the external points fixed at \(0\), \(1\), and \(z\in\mathbb C\setminus\{0,1\}\) [1509.07296]. In the four-dimensional normalization emphasized in later work, one equivalently writes \[ A_G(z_0,z_1,z_2) = \int_{(\mathbb R^d)^{|V_{\mathrm{int}}|}} \prod_{i\in V_{\mathrm{int}}}\frac{d^d x_i}{\pi^{d/2}} \prod_{e=\{u,v\}\in E(G)}\frac{1}{|x_u-x_v|^{2\nu_e}}, \] and then pulls back to a function of \(z\) by fixing \(z_0=0\), \(z_1=1\), \(z_2=z\) [2105.05015]. The literature uses more than one divergence convention. In the parametric treatment of Panzer and Schnetz, the superficial degree of divergence is \[ M_G=\lambda\sum_e \nu_e-\frac d2 |V_{\mathrm{int}}|, \] whereas in the spin-generalized formulation one writes \[ N_G=\sum_e \nu_e-\frac{\lambda+1}{\lambda}V_{\mathrm{int}} \] and factors out an overall scale so that the graphical function depends only on the conformal ratio \(z\) [1509.07296][2504.05850]. These are notationally different but serve the same structural purpose: they isolate the dimensionally correct three-point object. A basic point is that graphical functions are “graphical” in the sense that, after the external frame is fixed, the integral depends only on the combinatorics of the graph and the assigned edge weights. In scalar theories with convergent integrals, they are special three-point massless Feynman integrals and form a distinguished subclass of position-space amplitudes [1509.07296]. ## 2. Parametric representation and convergence A central structural result is the projective Schwinger-parametric representation. Introducing one Schwinger parameter \(\alpha_e>0\) for each edge, one defines the first and second Symanzik polynomials \[ \Psi_G(\alpha)=\sum_{T\text{ spanning tree of }G}\prod_{e\notin T}\alpha_e, \] and \[ \Phi_G(\alpha;z,\bar z) = \sum_{i\frac d2 |V_{\mathrm{int}}(g)|, \] and ultraviolet convergence requires that for every induced \(g\) with at most one external vertex, \[ \lambda\sum_{e\in g}\nu_e<\frac d2 (|V(g)|-1). \] Under these hypotheses the scalar integral converges absolutely [1509.07296]. The parametric form also yields analyticity. The key observation is that \(\Phi_G(\alpha;z,\bar z)\) is affine-linear in the squared distances \(s_{ij}=|x_i-x_j|^2\), with positive real parts in the appropriate domain. Combined with combinatorial degree bounds such as \[ \deg_g \Psi_G\ge |V_{\mathrm{int}}(g)|, \qquad \max_g \deg \Psi_G\le |V(g)|-1, \] and analogous estimates for forest polynomials, this permits differentiation under the integral and proves real analyticity on \(\mathbb C\setminus\{0,1\}\) via Morera’s theorem or dominated convergence for holomorphic parameters [1509.07296]. ## 3. Analytic structure, singularities, and symmetries In even dimensions \(d=2X+2\ge 4\), the scalar graphical function \(f_G(z)\) is real-valued on \(\mathbb C\setminus\{0,1\}\), invariant under complex conjugation, and single-valued real-analytic there [2105.05015]. In the scalar theory with general even dimension \(D=2\lambda+2\), Panzer–Schnetz theory further states that \(f_G(z)\) is a single-valued real-analytic function with possible logarithmic singularities at \(0\), \(1\), and \(\infty\) [2407.17133]. A refined local statement is available in \(d=4\) with all \(\nu_e=1\). At each singular point \(s\in\{0,1,\infty\}\), \(f_G(z)\) admits a single-valued log-Laurent expansion, with an explicit bound \[ M_s=\max\Bigl\{-\deg\bigl(G[V_{\rm int}\cup\{s\}]\bigr),\,1-\nu_s\Bigr\}, \] where \(\nu_s\) is the sum of edge weights attached to \(s\) [2105.05015]. This places graphical functions naturally on the punctured Riemann sphere. Their symmetry group is larger than complex conjugation. After “completion,” which adjoins the point \(\infty\) and balances external weighted valences, a permutation \(\sigma\) of \(\{0,1,z,\infty\}\) induces a Möbius transformation \(\phi(\sigma)\) of the cross-ratio variable, and one has \[ f_G(z)=f_{\sigma(G)}(\phi(\sigma)(z)). \] This \(S_4\)-cross-ratio symmetry is one of the main organizing principles of the theory [1302.6445][2105.05015]. Graphical functions also satisfy exact factorization identities. If removing the external vertices disconnects the graph into components, then the graphical function factorizes as the product of the graphical functions of the components [2105.05015]. In the spin-generalized setting the same factorization persists in tensorial form, \[ f_G^\alpha(z)=f_{G_1}^{\beta_1}(z)\,f_{G_2}^{\beta_2}(z), \] when the graph disconnects upon removal of \(\{0,1,z\}\) [2407.17133]. ## 4. Single-valued polylogarithms and recursive computation A major advance in the subject is the identification of graphical functions with explicit spaces of single-valued polylogarithmic functions. Brown’s single-valued multiple polylogarithms \(P_w(z)\) are built from ordinary iterated integrals \(L_w(z)\) by imposing trivial monodromy around \(0\), \(1\), and \(\infty\), with \(P_\emptyset=1\) and \(P_{0^n}(z)=(\ln z)^n/n!\) [1302.6445]. Many graphical functions can then be represented in the form \[ f_G(z)=\frac{g(z)}{z-\bar z}, \qquad g(z)\in\mathcal A, \] where \(\mathcal A\) is the shuffle algebra generated by the \(P_w(z)\) up to suitable weight [1302.6445]. This representation underlies an effective recursive calculus. In even dimensions, if \(G_1\) is obtained from \(G\) by attaching an edge of weight \(1\) at the external point \(z\), then the two graphical functions satisfy a Laplace-type equation \[ A_z^{(X-1)} f_{G_1}(z)=\frac{1}{(X-1)!}\,f_G(z), \] with \[ A_z^{(X-1)} = (z-\bar z)\,\partial_z\partial_{\bar z} +(X-1)\bigl(\partial_z(z-\bar z)-\partial_{\bar z}(z-\bar z)\bigr). \] Solving this PDE inside the single-valued function space yields the append-edge algorithm [2105.05015]. The algorithmic version for constructible graphs proceeds by completion, decomposition after removing \(\{0,1,z,\infty\}\), sequential-function building by repeated edge appending, evaluation of the period integral by the residue theorem inside the polylogarithmic algebra, and reduction to a standard MZV basis; the Maple package `polylog_procedures` implements this workflow [1302.6445]. Later expositions place the same strategy in the broader algebra of generalized single-valued hyperlogarithms and emphasize automated implementations such as `HyperlogProcedures` and `HyperFORM` [2604.25739]. The one-loop three-star already displays the characteristic analytic structure. In \(d=4\), \[ f_G^{(4)}(z) = \int\frac{d^4x}{\pi^2}\, \frac{1}{|x|^2\,|x-1|^2\,|x-z|^2} = \frac{2\,\mathrm{Im}\bigl(\mathrm{Li}_2(z)+\ln(1-z)\ln|z|\bigr)} {|z-\bar z|\,|z-1|^2}, \] equivalently \[ f_G^{(4)}(z)=\frac{D(z)}{|z-\bar z|\,(z-1)(\bar z-1)}, \] where \(D(z)\) is the Bloch–Wigner dilogarithm [2105.05015]. A 2026 review further formulates period extraction in terms of “anti-residues” at \(0\), \(1\), and \(\infty\), making Stokes-theoretic reduction explicit [2604.25739]. ## 5. Planar duality and momentum-space correspondences For externally planar graphs there is an exact planar-duality theorem. If \(G\) embeds in the plane with \(0\), \(1\), and \(z\) on the outer face, one constructs its dual \(G^\star\) by adjoining an extra vertex connected to \(0\), \(1\), and \(z\), dualizing, and then removing the dual triangle. When \[ M_G=\frac d2, \qquad \lambda\nu_{e^\star}=\frac d2-\lambda\nu_e, \] the dual graph also satisfies \(M_{G^\star}=d/2\), and the parametric polynomials obey \[ \Psi_G(\alpha)=\Phi_{G^\star}(\alpha), \qquad \Phi_G(\alpha)=\Psi_{G^\star}(\alpha). \] Consequently, \[ f_{G^\star}^{(\lambda)}(z) = f_G^{(\lambda)}(z)\, \prod_{e\in E_G}\frac{\Gamma(\lambda\nu_e)}{\Gamma(\lambda\nu_{e^\star})}. \] This theorem generalizes the familiar momentum-position duality for off-shell three-point integrals and is especially effective in four dimensions, where many dual graphs reduce to ladder or completed graphs with known closed forms [1509.07296]. A related formulation in the general even-dimensional theory is Fourier or position/momentum duality for planar graphs. If each propagator \(1/|x_u-x_v|^{2\nu_e}\) is replaced by its Fourier-transformed mate with weight \(\tilde\nu_e\) satisfying \[ \nu_e+\tilde\nu_e=\frac d2+1, \] then the graphical function becomes proportional to that of the planar dual graph \(G^\ast\); in parametric language this is associated with the Cremona transformation on Schwinger parameters [2105.05015]. Recent expository work pushes this correspondence further. The one-loop triangle is described as self-dual, and the general identity \[ I(\nu_1,\nu_2,\nu_3;z) = I(\lambda+1-\nu_1,\lambda+1-\nu_2,\lambda+1-\nu_3;z) \] is presented as the three-point manifestation of momentum-space/position-space self-duality [2604.25739]. ## 6. Periods, exact evaluations, and high-loop applications Graphical functions are a principal tool for computing Feynman periods. For a completed primitive \(\phi^4\) graph with four external edges \(0\), \(1\), \(z\), and \(\infty\), the period is \[ P(\Gamma)=\int f_\Gamma(z)\,\frac{d^2 z}{\pi}, \] and if a graph has only two external points, or one isolated external point, the corresponding graphical function is constant and coincides with the usual Feynman period [1302.6445][2105.05015]. The period theory is especially developed for sequential or zig-zag families. For the \(n\)-loop zig-zag graph \(Z_n\), the period was reduced to special values of sequential functions, then to combinations of \(\zeta(2^a,3,2^b)\), and finally evaluated in closed form by applying Zagier’s theorem: \[ P(Z_n) = 4\,\frac{(2n-2)!}{n!(n-1)!} \Bigl[1-\frac{1-( -1)^n}{2^{2n-3}}\Bigr]\zeta(2n-3). \] This establishes the Broadhurst–Kreimer zig-zag conjecture [1302.6445]. The same framework yields exact formulas modulo products for two further \(\phi^4\) families \(A_{m,n}\) and \(B_{m,n}\), and it extends algorithmically to a larger class of constructible graphs [1302.6445]. In even dimensions the method also covers fishnet ladder integrals \(G_{1,n}\), whose graphical functions in \(d=4\) are expressible in terms of single-valued classical polylogarithms of weight \(2n\) [2105.05015]. These structural results translate directly into perturbative calculations. Graphical-function methods have been used to compute renormalization constants and anomalous dimensions in four-dimensional \(\phi^4\) theory to seven and eight loops, and in six-dimensional \(\phi^3\) theory to five loops [2105.05015]. A later review presents graphical functions as a framework that also handles singular cases, dimensional regularization in non-integer dimension, and automated high-loop examples [2604.25739]. ## 7. Generalization to theories with spin The scalar theory has been generalized to theories with spin by introducing numerator propagators carrying multi-indices. In dimension \(D=2\lambda+2>2\), the spin-\(k\) numerator propagator is \[ Q^\alpha_\nu(x,y) = \frac{(y^{\alpha_1}-x^{\alpha_1})\cdots(y^{\alpha_k}-x^{\alpha_k})} {\|x-y\|^{2\lambda\nu+k}}, \qquad |\alpha|=k, \] with equivalent differential form \[ Q_{\nu;\alpha}(x,y) = \partial_y^{\alpha_1}\cdots\partial_y^{\alpha_k} \frac{1}{\|x-y\|^{2\lambda\nu-k}}. \] A general spinful Feynman integral is obtained by assigning such propagators to edges and contracting repeated indices with the Euclidean metric \(g^{\mu\nu}=\delta^{\mu\nu}\) [2407.17133][2504.05850]. In the three-point case one defines \[ A_G^\alpha(0,z_1,z_2)=\|z_1\|^{-2\lambda N_G}f_G^\alpha(z), \] and expands the tensor structure as \[ f_G^\alpha(z) = \sum_{\pi\in\Pi_2^\alpha} f_{G,\pi}(z)\,g^{\pi^0}\,\hat z_1^{\pi^1}\,\hat z_2^{\pi^2}, \] where \(\Pi_2^\alpha\) partitions the spin indices into metric-pairs and indices carried by the external unit vectors \(\hat z_1\) and \(\hat z_2\). The components \(f_{G,\pi}(z)\) are single-valued functions of \(z\) with at worst logarithmic singularities at \(0\), \(1\), and \(\infty\) [2407.17133]. The formalism has two notable algebraic consequences. First, nondegenerate bilinear pairings in the Brauer algebra allow spin-\(\alpha\) periods and two-point amplitudes to be reduced to linear combinations of scalar periods of augmented graphs [2407.17133]. Second, the action of differentiation with respect to the external point \(z_2\) is encoded by a block-triangular effective Laplace operator \(-2\,\mathrm{ex}_\lambda^\alpha\). Appending a scalar edge of weight \(1\) at \(z_2\) produces the differential equation \[ -2\,\mathrm{ex}_\lambda^\alpha\, f_{G_1}^\alpha(z) = -\frac{1}{\Gamma(\lambda)}\,f_G^\alpha(z), \] which can be inverted recursively sector by sector in the number of \(\hat z_2\)-factors, imposing single-valuedness and excluding poles of order \(\ge 2\lambda\) at \(z=1\) [2407.17133]. The function space remains remarkably stable under this generalization. Three-point graphical functions with numerators still lie in the same space of single-valued combinations of generalized hyperlogarithms, which preserves the feasibility of high-loop computations without Laporta IBP [2407.17133]. The 2025 full spin treatment applies these ideas to primitive periods in Yukawa–\(\phi^4\) theory up to eight loops; one completed ten-loop graph, corresponding to eight loops before completion, was evaluated numerically as \[ P_{\overline G}=2.7071150367306\ldots \] in \(15\) minutes on a single desktop core using the Maple package `HyperlogProcedures` [2504.05850]. This places spinful graphical functions within the same computational paradigm that had previously proved effective for scalar \(\phi^4\) and \(\phi^3\) theories. Source: https://www.emergentmind.com/topics/graphical-functions