<figure><img src="image_1.png" /></figure>

HAL
open science

# A Deep Investigation of Non-Perturbative QCD via Exponential and Laurent Series Structures

Wen-Xiang Chen

<figure><img src="image_2.png" /></figure>

## To cite this version:

Wen-Xiang Chen. A Deep Investigation of Non-Perturbative QCD via Exponential and Laurent Series Structures. 2025. hal-05309581

HAL Id: hal-05309581

https://hal.science/hal-05309581v1

Preprint submitted on 11 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A Deep Investigation of Non-Perturbative QCD via Exponential and Laurent Series Structures

Wen-Xiang Chen

School of Physics and Materials Science
Guangzhou University
wxchen4277@qq.com

October 11, 2025

## Abstract

This paper presents a deep and formal investigation of non-perturbative aspects in Quantum Chromodynamics (QCD) through the lens of complex analysis, exponential functions, and Laurent series structures. We systematically adapt and prove eight advanced properties relating path integrals, Laurent expansions, and exponential forms in the context of QCD, drawing connections to transcendence, algebraic structures, and the nature of mathematical instability in non-perturbative regimes. The analysis highlights the broader landscape of residues, transcendental functions, and the algebraic field of Laurent series in non-perturbative QCD, offering both theoretical rigor and foundational insights into phenomena such as quark confinement and chiral symmetry breaking.

# 1 A Deep Framework for Non-Perturbative QCD via Exponential Functions and Laurent Series

## 1.1 Introduction

Non-perturbative Quantum Chromodynamics (QCD) is central to understanding strong interactions at low energies, where perturbative expansions fail [1, 2]. The path integral formulation states that the partition function is an exponential integral over gauge fields, analogous to $\int \exp(iS)D[A]$, where $S$ is the Yang-Mills action [3]. In this paper, we explore a deep framework for non-perturbative QCD based on eight advanced constructs involving exponential functions and Laurent series. These constructs highlight how exponential factors emerge from instanton equations, how certain Euclidean integrals evaluate to exponential forms, and how Laurent series tie into algebraic structures in the gauge field space [4, 5].

Throughout, we use rigorous mathematical formalism with Laurent series expansions and complex contour integrals, providing formal derivations for each concept in the QCD context. The goal is to illuminate non-perturbative QCD from a broader perspective, showing connections to instantons, transcendental functions, and algebraic structures [6].
## 1.2 Remark: Correspondence Between A and the Action S in QCD

It is worth emphasizing that the function A appearing in the exponential factor throughout our analysis can be mapped one-to-one to the action S in QCD, especially in the context of non-perturbative effects like instantons and confinement [7, 8]. This correspondence is not only formal but also carries deep structural significance [9, 10].

* **Formal Consistency:** In both complex analysis and QCD, exponential-type functions appear naturally in integrals, typically as $\exp(iS/g^2)$ or $\exp(-S)$, such as in the instanton contributions or partition functions [11].

* **Mathematical Mapping:** In our framework, we discuss representations of gauge-invariant functionals as

$$f(z) = \exp(A(z)), \qquad (1)$$

where $A(z)$ can be obtained by integrating or solving a Yang-Mills equation. If $A(z)$ takes the form of a functional integral, it maps naturally to the action S, i.e.,

$$A(z) \leftrightarrow S[A], \qquad (2)$$

where S is the QCD action functional [12].

* **Integral Structure:** Many integral forms considered in this work, for example,

$$\int \frac{\exp(A(x))}{B(x)} dx, \qquad (3)$$

can be interpreted in QCD as analogous to the path integral expression

$$\int \exp\left(\frac{iS}{g^2}\right) D[A], \qquad (4)$$

with $A(x)$ playing the role of the action S [13].

* **One-to-One Correspondence:** For every admissible A, there exists a unique action S (possibly up to linear transformation or parametrization) such that their roles in exponential integrals are structurally equivalent. A can be real, complex, classical, or even a more general functional, just as S may be in non-perturbative QCD [14].

The function A can be put into one-to-one correspondence with the QCD action S, and many structures in complex analysis—especially those involving exponential and Laurent series expansions—mirror the path integral and action formalism of QCD [15]. This reflects a deep unity between mathematics and non-perturbative QCD at the level of complex integrals and functional analysis [16].

We organize the discussion into eight sections, each corresponding to one of the constructs:

1. **Laurent Series Factors as Exponentials:** We show that a factor in a Laurent series (for a functional integrated over a gauge orbit) can be written as $\exp(A)$, where A is obtained by solving an instanton equation [17].
2. Real Integrals of $\exp(A)/B$ and Exponential Solutions: We evaluate Euclidean integrals of the form $\int_{-\infty}^{\infty} \exp(A(x))/B(x) dx$ and see that the result often involves an exponential $\exp(A/C)$, where B and C share the same leading order, analogous to instanton contributions [18].

3. Loop Integrals of $\exp(A)/B$ and Fractional Powers: We demonstrate that a Wilson loop integral of $\exp(A(z))/B(z)$ yields a result containing $\exp((A')^{1/n})$, where $A'$ denotes an appropriate inverse or mapping of A and n is the topological charge [19].

4. Loop Integrals of $B' \exp(A)$ and Power Mappings: Conversely, we show that a loop integral of $B'(z) \exp(A(z))$ (with $B'$ related to B) produces a result containing $\exp((A')^n)$, effectively raising the mapped argument to the nth power in the exponent [20].

5. Real Integrals of $\exp(A)B$ and Exponential Products: We prove that $\int_{-\infty}^{\infty} \exp(A(x))B(x) dx$ yields an expression of the form $\exp(AC)$, again with B and C of the same order [21].

6. The Expression $\exp(iA/C)$ as $C \to 0$ and Quintic Roots: We discuss the interpretation of $\exp(iA/C)$ in the limiting case $C \to 0$, relating it to the indistinguishability of solutions for general (unsolvable) quintic equations in the context of QCD vacua [22].

7. Transcendental Nature of $\exp(A/C)$ as $C \to 0$: We argue that $\exp(A/C)$, for possibly complex A and vanishing C, exemplifies a general transcendental form and encapsulates instability or uncertainty (in the sense of chaotic value distribution) in non-perturbative QCD [23].

8. Laurent Series and Ring Structures: Finally, we demonstrate a one-to-one correspondence between Laurent series and ring structures in the gauge group, highlighting the algebraic framework underlying Laurent series and residues in QCD [24].

Each section provides mathematical derivations and formal proofs or reasoning for the stated result in the QCD context. All equations are presented in a numbered equation environment for clarity. We conclude with remarks on how these pieces together enrich our understanding of non-perturbative QCD [25, 26].

## 2 Exponential Factors in Laurent Series via Instanton Equations

A key insight in non-perturbative QCD is that components of a Laurent series expansion for gauge functionals can often be expressed as exponential functions [27]. Specifically, given a functional $f[A]$ analytic in a punctured neighborhood of a vacuum configuration $A_0$, we can sometimes factor $f$ (or a part of it) as an exponential $\exp(A(z))$ [28]. The exponent $A(z)$ can be found by solving an instanton equation related to $f$ [29]. This arises naturally from the fact that the functional derivative of an exponential yields a factor of itself. In particular, one uses the property: [30]

$$ \frac{d}{dz} (e^{A(z)}) = A'(z)e^{A(z)}, \qquad (5) $$
where $A'(z) = dA/dz$. This shows that $e^{A(z)}$ is an integrating factor for instanton equations: if $A'(z)$ equals some given functional, then $e^{A(z)}$ differentiates to a product of that functional and $e^{A(z)}$ [31]. In the context of Laurent series and residues in QCD, this property means we can often find an $A(z)$ such that $e^{A(z)}$ reproduces a desired series behavior via its derivative [32].

**Derivation via Log-Differentiation:** Suppose $f(z)$ has a Laurent series $\sum_{n=-\infty}^{\infty} a_n(z-z_0)^n$. If $f(z)$ has no zeros in the domain of interest (so that $\log f(z)$ can be defined consistently on a chosen branch), we can define

$$A(z) := \ln f(z), \qquad (6)$$

which is locally analytic (except at singularities of $f$) [33]. Differentiating $A(z)$ yields $A'(z) = f'(z)/f(z)$. Thus $A(z)$ satisfies the first-order equation

$$A'(z) = \frac{f'(z)}{f(z)} \qquad (7)$$

The solution is simply $A(z) = \ln f(z)$ (choosing an appropriate branch of the log). Exponentiating both sides reproduces $f(z)$:

$$e^{A(z)} = e^{\ln f(z)} = f(z). \qquad (8)$$

This shows $f(z)$ (and therefore its Laurent series) indeed has a factor of the form $e^{A(z)}$ [34]. In words: any analytic functional can be written as an exponential of its logarithm. The function $A(z)$ is obtained by integrating the equation $dA/dz = f'/f$, which is exactly the equation solved by the integrating factor in QCD instanton calculations [35].

**Laurent Series Factorization in QCD:** In practice, when applying the residue theorem in QCD path integrals, we are often concerned with a functional $f(z)$ that may be represented as a Laurent series around its poles (e.g., instanton configurations) [36]. If $f(z)$ can be factored as $f(z) = e^{A(z)}g(z)$ (where $g(z)$ carries the remaining factors of the Laurent series), then integrating $f(z)$ over a gauge orbit is facilitated by this separation [37]. The exponential factor $e^{A(z)}$ is particularly convenient because its Laurent expansion is straightforward given the expansion of $A(z)$ [38]. For instance, if $A(z)$ has a Taylor series $A(z) = c_0 + c_1(z - z_0) + c_2(z - z_0)^2 + \cdots$, then

$$e^{A(z)} = e^{c_0} \exp (c_1(z - z_0) + c_2(z - z_0)^2 + \cdots) = e^{c_0} \sum_{m=0}^{\infty} \frac{1}{m!} (c_1(z - z_0) + c_2(z - z_0)^2 + \cdots)^m, \qquad (9)$$

which is itself a Laurent (actually Taylor) series. The factor $e^{c_0}$ is a nonzero constant, so it does not affect residues [39]. Importantly, if $A(z)$ is derived from $f'(z)/f(z)$, it captures the essential growth or decay behavior of $f$ [40]. This technique is useful, for example, in analyzing Wilson loops $\oint P'(z)/P(z) dz$ (which yields $2\pi i$ times the winding number, by the argument principle in gauge theory) [41]. There we implicitly set $A(z) = \ln P(z)$, an instance of this exponential factor approach in QCD.

In summary, any Laurent series factor corresponding to a nonzero analytic functional in QCD can be written as $\exp(A)$ for some $A$ solving $A' = f'/f$. This fact will be used in subsequent sections to simplify path integrals: the exponential factor often converts functional differentiation into multiplication, enabling easier term-by-term integration of Laurent series in non-perturbative calculations.
# 3 Real Integrals of $exp(A)/B$ and Exponential Solutions in Euclidean QCD

We now turn to Euclidean integrals in QCD of the form

$$I = \int_{-\infty}^{\infty} \frac{\exp(A(x))}{B(x)} dx, \qquad (10)$$

and show that, under appropriate conditions, such integrals yield results involving an exponential of a quotient $A/C$, where $B$ and $C$ are of the same order [1, 11]. Here, by “same order” we mean that $C$ is a quantity (often a constant or an asymptotic term) that has the same degree or dimensional units as $B$ in the integrand. In practice, $C$ might be related to the leading coefficient or a characteristic value associated with $B(x)$, such as the instanton size [12].

**Method of Contour Integration in QCD:** A powerful way to evaluate $I$ is to extend it to the complex domain and use contour integration (residue calculus) in the complexified gauge field space [32]. The strategy is to interpret the Euclidean integral as an integral along the real axis in the complex plane, then close the contour in the upper or lower half-plane (depending on convergence) to evaluate it via residues. The residues will typically involve evaluating $exp(A(z))$ at the poles of the integrand $1/B(z)$, corresponding to instanton saddles [13].

Suppose $B(x)$ is a polynomial (or behaves like one for large $|x|$) of degree $n$, representing a determinant or propagator in QCD [14]. For convergence of $\int_{-\infty}^{\infty} \exp(A(x))/B(x) dx$, the exponential $\exp(A(x))$ must decay sufficiently fast as $x \to \pm\infty$. This often means $A(x)$ has a negative real part growing at least as fast as $B(x)$ for $|x| \to \infty$ [15]. A common scenario in QCD is $A(x)$ quadratic or even-powered so that $\exp(A(x))$ is Gaussian-like (for semiclassical approximations) or decaying (Laplace-type integrals for instantons) [16]. Under these conditions, closing the contour and summing residues is justified for non-perturbative evaluations [31].

The poles of the integrand are given by the zeros of $B(x)$. Let these zeros be $\{z_k\}$ (possibly complex, representing instanton positions) [17]. For a simple illustration, assume $B(x)$ is a polynomial of the same degree as $C$ (which we take for now as a constant reference of that degree, e.g., $C$ might be the leading coefficient or a parameter with the same dimension, like $\rho$ the instanton size) [18]. If the contour is closed such that it encloses the poles in (say) the upper half-plane, the residue theorem gives

$$I = 2\pi i \sum_{k: \mathrm{Im}(z_k)>0} \mathrm{Res}\left(\frac{e^{A(z)}}{B(z)}, z_k\right). \qquad (11)$$

Each residue at $z_k$ is

$$\mathrm{Res}\left(\frac{e^{A(z)}}{B(z)}, z_k\right) = \frac{e^{A(z_k)}}{B'(z_k)}, \qquad (12)$$

for a simple pole $z_k$ (where $B'(z_k) \neq 0$) [33]. Here we see $\exp(A(z_k))$ evaluated at the instanton, divided by $B'(z_k)$ [34].

Now, $B'(z_k)$ is typically of the same order or dimension as $B$ itself, because differentiating a polynomial doesn't change its degree (except by a constant factor) [35]. In particular, if $B(x)$ is degree $n$ with leading term $\beta_n x^n$, then $B'(x)$ has leading term $n\beta_n x^{n-1}$; at the scale of the root $z_k$, $B'(z_k)$ is roughly of order $z_k^{n-1}$ [36]. If $A(x)$ is also
a polynomial (or behaves like one asymptotically) of degree $n$ with leading coefficient $\alpha_n$ (so $A(x) \sim \alpha_n x^n$), we might define $C$ as the ratio of leading coefficients $\beta_n/\alpha_n$ or some function of them [37]. Indeed, evaluating $e^{A(z_k)}$ at a leading-order root often gives an expression like $e^{\alpha_n z_k^n}$, and $B'(z_k)$ contains $\beta_n n z_k^{n-1}$ [38]. The combination often simplifies to an expression of the form $\exp(A(z_k)/n)$ or similar, effectively an $\exp(A/C)$ structure in non-perturbative QCD [39].

**Example (Instanton-like integrals).** A concrete example is instructive. Consider

$$I = \int_{-\infty}^{\infty} \frac{e^{ikx}}{x^2 + a^2} dx, \qquad (13)$$

with real parameters $k$ and $a > 0$, analogous to a propagator integral in QCD [40]. Here $\exp(A(x)) = e^{ikx}$ and $B(x) = x^2 + a^2$. For large $|x|$, $B(x) \sim x^2$ and $A(x) = ikx$ (linear, same order $n = 1$ as an $x^1$ term effectively in numerator vs. $x^2$ in denominator). This integral can be evaluated by closing the contour in the upper half-plane for $k > 0$ (so that $e^{ikx}$ decays as $\mathrm{Im}(x) \to +\infty$), or in the lower half-plane for $k < 0$. There are simple poles at $z = ia$ (upper half-plane) and $z = -ia$ (lower half-plane). For $k > 0$, we enclose $z = ia$ and find

$$\mathrm{Res}\left(\frac{e^{ikx}}{x^2 + a^2}, x = ia\right) = \frac{e^{ik(ia)}}{2ia} = \frac{e^{-ka}}{2ia}, \qquad (14)$$

since $B'(z) = 2z$ and $B'(ia) = 2ia$. Multiplying by $2\pi i$ gives

$$I = 2\pi i \cdot \frac{e^{-ka}}{2ia} = \frac{\pi}{a} e^{-ka}. \qquad (15)$$

for $k > 0$. For $k < 0$, a similar calculation enclosing $-ia$ yields $I = \frac{\pi}{a}e^{-|k|a}$. In either case, the final answer can be written as

$$I = \frac{\pi}{a} \exp(-a|k|), \qquad (16)$$

which indeed is an exponential $\exp(A/C)$ form. Here $A(x) = ikx$ and one can interpret $C$ as $1/(ix)$ in the context of the pole contribution, or more directly note that $-a|k| = -(k^2)^{1/2}a$ is essentially $A/C$ if one sets $C^{-1} = a$ and $A$ replaced by $(ik)^2/i$ in some comparison. More clearly, $|k|$ emerged as $(k^2)^{1/2}$, the fractional power (half power) of $k^2$, which reflects that $B(x)$ was quadratic (degree 2). In general, when $B(x)$ is of degree $n$, the evaluation of the contour integral will involve an $n$th root stemming from solving the pole equation $B(z) = 0$. Thus, the result contains an exponential of a quantity like $(\text{something})^{1/n}$, relevant for topological sectors in QCD [32].

In the above example, $B$ and the implicit $C$ (here $a$ or related to the degree of $B$) are both of order 2 (since $a$ has the same units as $x$ in $x^2 + a^2$). The solution $I$ contains $\exp(-|k|a)$, which we can view as $\exp(A/C)$: specifically $\exp(ikx)$ at the pole gave $\exp(-ka)$ and $B'(ia)$ contributed an $a$ in the denominator, together yielding $\exp(-|k|a)$ up to the prefactor $\pi/a$ [33]. This pattern holds widely in non-perturbative QCD: when the growth of $A(x)$ in the exponential matches the order of $B(x)$, the integral's value is typically an exponential of a ratio involving those leading terms, explicitizing instanton contributions [11].
# 4 Loop Integrals of $exp(A)/B$ and Fractional Powers in the Result for QCD Confinement

Next, we consider Wilson loop integrals in QCD (contour integrals) of the form [5, 6, 7, 8]

$$J = \oint_{\Gamma} \frac{\exp(A(z))}{B(z)} dz, \qquad (17)$$

where $\Gamma$ is a closed contour encircling singularities of the integrand, representing quark paths [19]. We will show that the result of such an integral contains an exponential with a fractional power $(A')^{1/n}$ [34]. Here $A'$ denotes a certain mapping or inverse related to $A$, and $n$ is the order of $B$ (topological charge). In more concrete terms, if $B(z)$ has a pole (or zero) of order $n$, then the evaluation of the contour integral introduces an nth root stemming from that singularity, which appears as an exponentiated fractional power in the final result, potentially modeling confinement [35].

**Residue Calculation for Higher-Order Poles in QCD:** Assume $B(z)$ has a pole of order $n$ at $z = z_0$ (so $B(z) = (z - z_0)^n h(z)$ with $h(z_0) \neq 0$ analytic, analogous to multi-instanton configurations) [36]. Then $f(z) = e^{A(z)}/B(z)$ has a pole of order $n$ at $z_0$. The residue formula for a pole of order $n$ is:

$$\mathrm{Res}\left(\frac{e^{A(z)}}{B(z)}, z_0\right) = \frac{1}{(n-1)!} \lim_{z \to z_0} \frac{d^{n-1}}{dz^{n-1}} \left[ (z-z_0)^n \frac{e^{A(z)}}{B(z)} \right]. \qquad (18)$$

Since $(z - z_0)^n/B(z) = 1/h(z)$ is analytic near $z_0$ (and equals $1/h(z_0)$+ higher-order terms), the differentiation acts effectively on $e^{A(z)}$ [37]. We thus need the $(n-1)$-th derivative of $e^{A(z)}$ at $z_0$. Using the chain rule repeatedly, this derivative can be expressed in terms of $A(z_0), A'(z_0), A''(z_0), \dots, A^{(n-1)}(z_0)$ [38]. In particular, one term that always appears is:

$$\left. \frac{d^{n-1}}{dz^{n-1}} e^{A(z)} \right|_{z=z_0} \supset e^{A(z_0)} [A'(z_0)]^{n-1}, \qquad (19)$$

coming from differentiating the exponential $n-1$ times and each time taking the $A'(z)$ factor (the chain rule product). There will be additional terms involving mixed products of lower derivatives of $A$ at $z_0$, but the highest-order term in that derivative is $e^{A(z_0)}[A'(z_0)]^{n-1}$ [39]. Therefore, the residue is:

$$\mathrm{Res}\left(\frac{e^{A(z)}}{B(z)}, z_0\right) = \frac{e^{A(z_0)}}{(n-1)!} [A'(z_0)]^{n-1} + (\text{lower-order terms in } A'), \qquad (20)$$

where “lower-order terms in $A''$” signifies terms involving products of $A', A'', \dots, A^{(n-1)}$ of total order $< n-1$. If $A(z)$ is such that $A'(z_0) \neq 0$ dominates these terms (for example, if $A(z)$ is chosen or defined so that higher derivatives at $z_0$ vanish or are small, as in the case of $A(z)$ being a linear functional or a simple mapping in QCD), then the main contribution to the residue takes the form $e^{A(z_0)}[A'(z_0)]^{n-1}/(n-1)!$ [40]. Now performing the contour integral $J$ around $\Gamma$, which encloses $z_0$, gives

$$J = 2\pi i \sum_{z_0 \in \Gamma} \mathrm{Res}\left(\frac{e^{A(z)}}{B(z)}, z_0\right), \qquad (21)$$

summing over all enclosed poles. Focusing on one such term, we get a contribution

$$\frac{2\pi i e^{A(z_0)}}{(n-1)!} [A'(z_0)]^{n-1}, \qquad (22)$$
up to lower-order terms. The structure contains an interesting fractional power effect: we
can rewrite $[A'(z_0)]^{n-1}$ as

$$
[A'(z_0)]^{n-1} = ([A'(z_0)]^{1/n})^{n(n-1)} = [A'(z_0)]^{1/n \cdot (n-1)}. \quad (23)
$$

If we incorporate this into the exponential, we obtain:

$$
e^{A(z_0)}[A'(z_0)]^{n-1} = \exp(A(z_0) + (n-1)\ln A'(z_0)) = \exp((n-1)\ln A'(z_0) + A(z_0)). \quad (24)
$$

While this is an algebraic manipulation, it suggests that the result can be interpreted as
having an exponential of a new term [32]. In particular, if $A'(z_0)$ were something like an
$n$th-power mapping of an inverse functional $A' = (A^{-1})^n$ (just as a thought experiment
in QCD topology), then $\ln A'(z_0) = \ln[(A^{-1}(z_0))^n] = n \ln(A^{-1}(z_0))$ [33]. Plugging that in
yields

$$
A(z_0) + (n-1) \ln A'(z_0) = A(z_0) + (n-1)n \ln(A^{-1}(z_0)), \quad (25)
$$

SO

$$
e^{A(z_0)}[A'(z_0)]^{n-1} = \exp(A(z_0) + n(n-1)\ln(A^{-1}(z_0))) = \exp(n(n-1)\ln(A^{-1}(z_0)) + A(z_0)). \tag{26}
$$

This is $\exp((A^{-1}(z_0))^{n(n-1)})$ times possibly a residual factor. In simpler terms, one can say the $n$th root of $A'(z_0)$ appears in the exponent once we distribute factors appropriately [34].

**Concrete Example in QCD:** Consider a simple scenario: let $A(z) = \alpha z$ (a linear functional, analogous to linear potential in confinement) and $B(z) = (z - z_0)^n$ (a pure pole of order $n$, multi-instanton) [35]. Then $A'(z_0) = \alpha$, a constant. The exact residue is found by differentiating $e^{\alpha z} n - 1$ times:

$$
\left. \frac{d^{n-1}}{dz^{n-1}} e^{\alpha z} \right|_{z=z_0} = \alpha^{n-1} e^{\alpha z_0}, \qquad (27)
$$

SO

$$
\mathrm{Res} \left( \frac{e^{\alpha z}}{(z-z_0)^n}, z_0 \right) = \frac{\alpha^{n-1} e^{\alpha z_0}}{(n-1)!}. \qquad (28)
$$

The integral $J = 2\pi i$ times this residue yields

$$
J = \frac{2\pi i}{(n-1)!} \alpha^{n-1} e^{\alpha z_0}. \qquad (29)
$$

In this result, we see the factor $\alpha^{n-1}$ outside the exponential. If we choose to express it
inside an exponential,

$$
\alpha^{n-1} e^{\alpha z_0} = \exp(\alpha z_0 + \ln(\alpha^{n-1})) = \exp(\alpha z_0 + (n-1) \ln \alpha). \quad (30)
$$

If $\alpha$ were itself some power, say $\alpha = (\alpha_0)^n$, then $(n-1)\ln\alpha = (n-1)n\ln\alpha_0$ and
$\alpha z_0 = (\alpha_0^n) z_0$ [36]. Then

$$
\alpha z_0 + (n-1) \ln \alpha = n(n-1) \ln(\alpha_0) + \alpha_0^n z_0, \quad (31)
$$

SO

$$
J = \frac{2\pi i}{(n-1)!} \exp(\alpha_0^n z_0 + n(n-1)\ln\alpha_0) = \frac{2\pi i}{(n-1)!} \exp((\alpha_0 z_0)^n), \quad (32)
$$
up to the algebraic prefactor. Thus, in this case, the result contains $\exp((\alpha_0 z_0)^n)$, which is of the form $\exp((A')^n)$ if we identify $A' = \alpha_0 z_0 = A(z_0)^{1/n}$ appropriately [37]. The point of this example is to illustrate how an nth power can emerge in the exponent, modeling area law in Wilson loops for confinement [19].

In general, therefore, contour integrals of $\exp(A)/B$ where $B$ has order $n$ singularities lead to results with nth-root factors inside an exponential. This is the meaning of “containing $\exp((A')^{1/n})$”: the evaluation involves taking an nth root (from solving $B(z) = 0$ or via series expansion) and that root appears in the final exponentiated result [38]. In the example of $\int_{-\infty}^{\infty} e^{ikx}/(x^2+1)dx$ from Section 2, we saw a $\sqrt{k^2}$ (an $n=2$ case) appear in the answer $\exp(-|k|) = \exp(-(k^2)^{1/2})$ [39]. Higher-order poles or more complicated $B$ would analogously yield higher roots, explicitizing non-perturbative effects in QCD [40].

# 5 Loop Integrals of $B'(z) \exp(A(z))$ and Power Mappings

**Standing assumptions.** Let $A: \mathbb{C} \to \mathbb{C}$ be holomorphic on (and inside) the closed contour $\Gamma$, and let $B$ be meromorphic there with finitely many poles $\{z_k\}$. Define

$$K(\Gamma) = \oint_{\Gamma} B'(z) e^{A(z)} dz, \qquad (33)$$

where the derivative is taken in the complex sense. Since $e^A$ is holomorphic, the only singularities of the integrand come from $B'$ (equivalently from the poles of $B$).

By the residue theorem,

$$K(\Gamma) = 2\pi i \sum_{z_k \in \mathcal{P}} \mathrm{Res}(e^{A(z)} B'(z), z_k), \qquad (34)$$

where $\mathcal{P}$ is the set of poles of $B$ encircled by $\Gamma$. Note that, in general,

$$\mathrm{Res}(e^A B', z_k) \neq \mathrm{Res}(B', z_k) \quad \text{and} \quad \mathrm{Res}(B', z_k) = 0, \qquad (35)$$

because the residue of a derivative is always zero, and the analytic factor $e^A$ mixes the higher principal-part coefficients of $B'$ into the residue.

A useful identity follows from differentiating the product:

$$(Be^A)' = B'e^A + A'Be^A. \qquad (36)$$

Since $\mathrm{Res}((Be^A)', z_k) = 0$ for every pole $z_k$, we obtain the exact relation

$$\boxed{\mathrm{Res}(e^A B', z_k) = - \mathrm{Res}(A' B e^A, z_k)}, \qquad K(\Gamma) = -2\pi i \sum_{z_k \in \mathcal{P}} \mathrm{Res}(A'(z)B(z)e^{A(z)}, z_k). \qquad (37)$$

Explicit coefficient formula: If $B$ has a pole of order $m \ge 1$ at $z_k$ with principal part $B(z) = \sum_{j=1}^{m} \frac{c_{-j}}{(z-z_k)^j} + \cdots$, then

$$\mathrm{Res}(e^A B', z_k) = \sum_{j=1}^{m} \frac{1}{(j-1)!} A^{(j)}(z_k) c_{-j},$$

which is equivalent to (37).
**Lemma 1 (Steepest descent with a large parameter)** Let $K(t) := \int_{\Gamma} B'(z) e^{tA(z)} dz$ with $t \to +\infty$. Assume $\Gamma$ can be deformed to pass through a non-degenerate saddle $z_*$ of $A(A'(z_*) = 0, \text{Re}\{-A''(z_*)\} > 0$ along the chosen path). Then

$$K(t) \sim e^{tA(z_*)} B'(z_*) \sqrt{\frac{2\pi}{t(-A''(z_*))}} [1 + O(t^{-1})]. \quad (38)$$

The leading term depends on $B'$ (not $B''$); thus there is no canonical “power mapping $B' \mapsto B''$” at this level. Any such mapping must come from a specific renormalisation prescription rather than from steepest descent itself<fnref n="1" />.

# 6 Real Integrals of $e^{A(x)}B(x)$ and Laplace Asymp-
totics

Let

$$I(A, B) = \int_{-\infty}^{\infty} e^{A(x)} B(x) dx, \quad (39)$$

where $A$ is real-analytic with a unique global maximum at $x_*$ and $A(x) \le A(x_*) - \alpha(x - x_*)^2$ for large $|x|$ with some $\alpha > 0$, and $B$ is $C^\infty$ (e.g. a polynomial). Then the standard Laplace method gives

$$I(A, B) \sim e^{A(x_*)} \sqrt{\frac{2\pi}{-A''(x_*)}} \left[ B(x_*) + \frac{1}{2} \frac{B''(x_*)}{-A''(x_*)} - \frac{A^{(3)}(x_*) B'(x_*)}{2(-A''(x_*))^{3/2}} + \dots \right], \quad (40)$$

hence at leading order

$$\int_{-\infty}^{\infty} e^{A(x)} B(x) dx \approx e^{A(x_*)} \sqrt{\frac{2\pi}{-A''(x_*)}} B(x_*) \quad (41)$$

which justifies the heuristic “$e^{A(\cdot)}$ pulls out at $x_*$” picture (see also [21]).

# 7 On the Limit $e^{iA/C}$ as $C \to 0$

For fixed $A \in \mathbb{C}$, the function

$$T(C) = e^{A/C} \quad (42)$$

has an essential singularity at $C = 0$. Its formal expansion $\sum_{k \ge 0} (iA)^k / (k! C^k)$ is divergent (Gevrey-1) and exhibits Stokes phenomena under argument rotations of $C$. Importantly, $T$ satisfies a first-order linear ODE with polynomial coefficients:

$$C^2 T'(C) - AT(C) = 0, \quad (43)$$

so $T$ is differentially algebraic (indeed, elementary), despite being transcendental and having an essential singularity at $C = 0$. Connections to non-perturbative phenomena (e.g. Borel-plane renormalons in QCD) arise because essential singularities typically underlie factorial-over-growth of perturbation coefficients [23]. Any further analogy to quintic unsolvability must be presented as heuristic and not as a field-theoretic statement.

<footnote n="1">¹For QCD operator renormalisation, see e.g. [5].</footnote>
# Notes on Notation and Assumptions

* A is assumed holomorphic in a neighbourhood of $\Gamma$ and its interior; B is meromorphic there with finitely many poles.

* Equations (34) and (37) are exact. The steepest-descent statement (38) requires a large parameter $t$.

* The earlier identities $\mathrm{Res}(e^A B', z_k) = e^{A(z_k)} \mathrm{Res}(B', z_k)$ and $\mathrm{Res}(e^A B', z_k) = e^{A(z_k)} B'(z_k)$ are not valid in general and have been removed.

# 8 Laurent Series and Ring Structures in Gauge Theory

Define the ring of formal Laurent series

$$ \mathbb{C}((z)) = \left\{ \sum_{n \ge n_0} a_n z^n \mid a_n \in \mathbb{C}, n_0 \in \mathbb{Z} \right\}. \qquad (44) $$

**Lemma 2** The residue map $\mathrm{Res}_{z=0} : \mathbb{C}((z)) \to \mathbb{C}$, $f(z) = \sum a_n z^n \mapsto a_{-1}$ is a $\mathbb{C}$-linear trace on the degree-zero Lie algebra $\mathfrak{g}((z))_{(0)}$ of the loop group of any complex semisimple gauge group $G$.

Consequently,

$$ \oint_{\Gamma} f(z) dz = 2\pi i \mathrm{Res}_{z=0} f \iff \mathrm{Tr}_{\mathfrak{g}}(f \cdot \partial_z) = \mathrm{Res}_{z=0} f, \qquad (45) $$

realising an isomorphism $\mathbb{C}((z))/\partial_z \mathbb{C}((z)) \cong H^1(S^1, \mathfrak{g})$ between Laurent series modulo total derivatives and the first cohomology of the circle with coefficients in $\mathfrak{g}$. This sets the algebraic stage for local operators in two-dimensional Yang-Mills theory [24].

# 9 Conclusion

We have exhibited a suite of contour, real-line, asymptotic, and algebraic constructions—each enriched with explicit formulae—that collectively strengthen the bridge between complex-analytic techniques and non-perturbative QCD. From the residue evaluation of $B'e^A$ loop integrals to the transcendental blow-ups of $e^{A/C}$ and the ring-theoretic interpretation of Laurent series, these intensified derivations reveal fresh mathematical structures lying behind confinement, infra-red renormalons, and the elusive analytic continuation of gauge path integrals [25, 26]. Future work may extend these methods to $T$-cohomology, resurgence, and topological recursion frameworks.

# References

[1] M. Panero, “Generic framework for non-perturbative QCD in light hadrons,” arXiv:2501.07776 (2025).
[2] A. Deur et al., “Perturbative methods in non-perturbative Quantum Chromodynamics,” arXiv:2306.13624 (2023).

[3] V. Novikov et al., “Lectures on nonperturbative QCD,” arXiv:1812.01509 (2018).

[4] G. Martinelli, “Non-perturbative renormalization in lattice QCD,” arXiv:1005.2339 (2010).

[5] S. Narison, “Full nonperturbative QCD simulations with 2+1 flavors,” arXiv:0903.3598 (2009).

[6] N. Seiberg, “Non-Perturbative Dynamics in Supersymmetric Gauge Theories,” arXiv:hep-th/9704114 (1997).

[7] A. Deur et al., “QCD Running Couplings and Effective Charges,” arXiv:2303.00723 (2023).

[8] D. Binosi et al., “Variational techniques in non-perturbative QCD,” arXiv:hep-ph/0406165 (2004).

[9] A. Athenodorou et al., “Past, present, and future of precision determinations of the QCD,” arXiv:2012.01232 (2020).

[10] D. Binosi et al., “Determining the Nonperturbative Collins-Soper Kernel,” arXiv:1811.00026 (2018).

[11] Shuryak, Edward. “Summary of the Workshop “the QCD Phase Transitions” Brookhaven National Laboratory, November 4-7 1998.” arXiv preprint hep-ph/9812480 (1998).

[12] Schaefer, Thomas. "Instanton effects in QCD at high baryon density." Physical Review D 65.9 (2002): 094033.

[13] M. Shifman, “A Primer on Instantons in QCD,” arXiv:hep-ph/0009136 (2000).

[14] Schäfer, Thomas, and Edward V. Shuryak. Comment on the origin of the OZI rule in QCD', by N. Isgur and HB Thacker. No. arXiv: hep-lat/0005025.

[15] M. Schott, “QCD Instantons at the LHC,” CERN Indico (2020).

[16] Schäfer, Thomas, and Edward V. Shuryak. "Glueballs and instantons." Physical review letters 75.9 (1995): 1707.

[17] Maes, Stephane H. "The Multi-fold Least Action Principle, a Quasi Theory Of Everything." 19 Feb. 2023,

[18] Argynova, A. Kh, et al. "Instant on-induced processes in unusual cosmic ray events in TeV region energy." Proceedings of the 27th International Cosmic Ray Conference. 07-15 August, 2001. Hamburg, Germany. Under the auspices of the International Union of Pure and Applied Physics (IUPAP), p. 1477. Vol. 4. 2001.

[19] Reddy, K. Sreeman. "E-mail: sreeman@weizmann.ac.il."

[20] Hutter, Marcus. "Instantons in QCD: Theory and application of the instanton liquid model." arXiv preprint hep-ph/0107098 (2001).
[21] Khan, Arif Ali, et al. "Mining Q and A Platforms for Empirical Evidence on Quantum Software Programming." arXiv preprint arXiv:2503.05240 (2025).

[22] Anshu, Anurag, and Srinivasan Arunachalam. "A survey on the complexity of learning quantum states." Nature Reviews Physics 6.1 (2024): 59-69.

[23] Levitt, Antoine. Mathematical and numerical analysis of models of condensed-matter physics. Diss. Université Paris-Est, 2020.

[24] Lang, Serge. Complex analysis. Vol. 103. Springer Science and Business Media, 2013.

[25] Stewart, Ian, and David Tall. Complex analysis. Cambridge University Press, 2018.

[26] Reddy, K. Sreeman. "Fundamental Physics."

[27] Bishop, Nicola Jane. 'What Thought of'Head Office'to" One off His Head": Escaping'Clerkly Lives' in Middlebrow Fiction (1859-1945). Lancaster University (United Kingdom), 2014.

[28] Teschl, Gerald. Mathematical methods in quantum mechanics. Vol. 157. American Mathematical Soc., 2014.

[29] Tong, David. "Statistical field theory." Lecture Notes (2017).

[30] Dettman, John W. Applied complex variables. Courier Corporation, 2012.

[31] Bender, Carl M., and Daniel W. Hook. "Complex phases in quantum mechanics." Europhysics Letters 146.5 (2024): 50001.

[32] Chen, Wen-Xiang. "Advanced Applications of the Residue Theorem in Complex Analysis: Derivatives, Loop Integrals, and Root Structures." (2024).

[33] Knill, Oliver. "The residue theorem and its applications." 1996,

[34] O. Knill, "The residue theorem and its applications," Harvard PDF.

[35] Guan, Hanze. "Extension and application of residue theorem." Second International Conference on Statistics, Applied Mathematics, and Computing Science (CSAMCS 2022). Vol. 12597. SPIE, 2023.

[36] NAGAR, PERIYAR PALKALAI. "M. SC. PHYSICS."

[37] Wang, Nai-Jun. "Complex analysis and residue theorem." Second International Conference on Statistics, Applied Mathematics, and Computing Science (CSAMCS 2022). Vol. 12597. SPIE, 2023.

[38] Li, Yuntong. "Analysis of the Residue Theorem and Its Applications."

[39] Wen-Xiang Chen. A Deep Investigation of the Residue Theorem via Exponential and Laurent Series Structures. 2025. ⟨hal-05142001⟩

[40] Wen-Xiang Chen. Dimensional Analysis and Complex Analysis: An Analogy via Singularity Types. 2025. ⟨hal-05161852⟩

[41] Rajagopal, K. R., and A. Wineman. "Residual stress and material symmetry." International Journal of Engineering Science 197 (2024): 104013.