Multiloop calculations with parametric integration in critical dynamics: the four-loop analytic study of model A of $φ^4$ theory

Abstract: We perform an analytical four loop calculation of exponent $z$ in model A of critical dynamics in $d=4-2\varepsilon$ dimensions. This is the first time such a large order of perturbation theory has been calculated analytically for models of critical dynamics. To do this, we apply the modern method of parametrical integration with hyperlogaritms. We discuss in detail peculiarities of application of this method to critical dynamics, e.g. the problem of linear-irreducible diagrams already present in four loop (contrary to statics where the first linear-irreducible diagram appears in six loop).

• The paper presents the first analytic four-loop evaluation of the dynamic critical exponent z in model A of φ⁴ theory.
• It leverages advanced parametric integration with hyperlogarithms and a novel stream decomposition method for precise diagram reduction.
• The results validate high-precision numerical studies and offer new insights into mitigating linear-irreducibility in dynamic RG computations.

Four-loop Analytic Computation of the Dynamic Exponent $z$ in Model A of $\phi^4$ Theory via Parametric Integration

---

Introduction

The paper "Multiloop calculations with parametric integration in critical dynamics: the four-loop analytic study of model A of $\phi^4$ theory" (2512.10591) presents the first analytic four-loop evaluation of the dynamic critical exponent $z$ in the context of $O(n)$-symmetric model A, formulated in $d=4-2\varepsilon$ dimensions. This model encompasses a dissipative generalization of the classic static $\phi^4$ theory, crucial for describing dynamic critical phenomena such as the slowdown of relaxation time near criticality in a wide variety of systems—ranging from magnetic ordering, alloys, multiferroics, and superfluid transitions to collective behavior in biological and social systems.

The work addresses the longstanding computational barrier in analytic multiloop treatments of dynamic models, which notoriously lag behind static counterparts due to the complexity of time-correlated propagator structures and diagram proliferation. The study leverages modern Feynman parametric integration methods—specifically, hyperlogarithms/Goncharov polylogarithms—and introduces new algorithmic techniques for diagram reduction and handling linear-irreducibility, achieving analytic control over previously numerically intractable terms.

---

Renormalization Framework and Diagrammatic Structure

The action for model A is renormalized in the MS scheme with explicit account of the time derivative and spatial relaxation terms. Owing to the underlying symmetry and dissipative nature of the model, the renormalization constants for the static fields and couplings coincide with those of the static $\phi^4$ theory, while $Z_\lambda$ fully encodes the intrinsic dynamic contributions, extractable from the one-particle irreducible $\langle \psi' \psi' \rangle$ function at zero frequency.

Perturbative expansions up to four loops necessitate precise cancellation and organization of $\varepsilon$-pole structures, requiring detailed analysis of diagrammatic divergences, especially those associated with dynamic propagators $\sim (-i\omega+\lambda k^2)^{-1}$ (enforcing time-causality and leading to "time cuts" in diagrammatic notation). The exponential time dependence dramatically increases the number of time-ordered diagrams and enforces frequency-momentum dependence at every loop order.

---

Modified Diagram Reduction and Renormalization

Building on previous numeric sector-decomposition approaches for dynamic RG calculations [AIKV_4lSD17, AEHIKKZ_5l_2_22], this study implements an enhanced analytic diagram reduction scheme. Unlike in static models—where integration by parts (IBP) and counting arguments yield manageable sets of master integrals—dynamic models retain a vast array of time-topologies, many with complicated subdivergences.

The authors systematize the reduction of diagrams, exploiting equivalence between renormalization constants in the static and dynamic sectors to target the elimination of dynamic subgraphs wherever possible. This procedure is not simply a combinatorial optimization, but is tailored to enable analytic subtraction of dynamic subdivergencies at the integrand level. Notably, for four-loop diagrams, the reduction algorithm reduces the number of irreducible dynamic subgraphs from 17 to just one, facilitating parametric methods.

For the remaining problematic subgraphs, they employ a BPHZ-like scheme with IR rearrangement to render divergent integrals convergent, ensuring compatibility with parametric integration and regularization requirements.

---

Parametric Integration and Linear Reducibility

Feynman parametric representation is then employed, recasting loop diagrams into multidimensional integrals over simplex-constrained Schwinger (Feynman) parameters. The Maple-based HyperInt package [Panzer15] automates the symbolic integration as long as the integrand exhibits linear reducibility: each step must allow representation in terms of hyperlogarithms/GPLs, with non-trivial dependence determined by the diagram's spanning tree structure and the algebraic complexity of dynamic time cuts.

While most diagrams at four loops are linearly reducible, the $C_9$ diagram reveals the onset of genuine linear-irreducibility, manifested in quadratic dependencies on dynamic Feynman parameters. Standard rationalizing change-of-variables techniques are inefficient; instead, the authors introduce a "stream decomposition" approach—splitting the integrand into partitions that are individually linearly reducible in different integration orders and compensating with explicit analytic counterterms for subdivergent parts. This enables the analytic extraction of $\varepsilon$-pole coefficients for $Z_1$ at four loops.

All analytic expressions for diagram contributions are cross-validated against sector-decomposition numerical results for high-precision agreement.

---

Dynamic Exponent $z$: Analytic Results and Comparison

Employing the computed coefficients, the dynamic exponent $z$ is constructed via the $\varepsilon$ expansion, to $\mathcal{O}(\varepsilon^4)$,

$$z(\varepsilon,n) = 2 + \varepsilon^2 \frac{2(n+2)}{(n+8)^2}[6\ln(4/3)-1] + \varepsilon^3 \text{(higher polylog terms)} + \varepsilon^4 \text{(numerical coefficients)} + \mathcal{O}(\varepsilon^5)$$

with explicit expressions for three-loop contributions in terms of dilogarithms, and the four-loop coefficient given numerically due to the complexity of the $C_9$ diagram. The numerical value for $z$ at the Ising ($n=1$) universality class and $d=3$ is in full consistency with established five-loop numerical RG studies [AEHIKKZ_5l_2_22], e.g., $z = 2.0236(8)$.

Partial analytic cross-checks are performed via 1/$n$ expansions, confirming correspondence to previous results [HHM_2l2], with further simplification of GPL expressions shown possible (Appendices).

---

Practical and Theoretical Implications

This paper demonstrates the feasibility and limitations of analytic parametric-integration approaches for dynamic critical models at high perturbative orders. The onset of linear-irreducibility due to time cuts in dynamic diagrams is shown to be an intrinsic obstruction, likely to proliferate at higher loops and in more complex models. The stream decomposition method provides a template for circumventing this barrier for certain cases, but the complexity of resulting GPL alphabets and the challenge of further algebraic reduction remain open problems.

From a practical perspective, the work enables direct computation of physical observables (critical slowdown exponents) relevant to phase transition kinetics, magnetic relaxation, and emergent behavior in complex systems. On a theoretical level, the methodology advances analytic multiloop RG approaches in non-equilibrium statistical field theory, with scope for further development in both symbolic integration tools and diagram reduction algorithms.

---

Conclusion

This work establishes analytic four-loop control over the dynamic critical exponent $z$ in model A of $\phi^4$ theory, employing parametric integration with hyperlogarithms and advanced diagram reduction techniques. The analytic procedure, augmented with stream decomposition for linearly irreducible integrals, robustly reproduces and supports prior high-precision numerical results, clarifies the nature of diagrammatic complexity in dynamic critical theories, and sets the stage for further progress at higher loops and in broader classes of models. The strategies developed herein will inform developments in both symbolic analytic tools and RG analysis of critical dynamics for years to come.