⚛️ quantum physics

A nonrecursive method for computing the off-diagonal small-time heat kernel expansion

This paper presents a nonrecursive method for computing closed-form expressions of the off-diagonal small-time heat kernel expansion coefficients up to second order for the Klein-Gordon operator in flat space-time with electromagnetic fields, extending previous diagonal-only approaches and verifying accuracy against plane wave and constant field cases.

Original authors: Vyacheslav O. Guba, Alexey V. Reznichenko

Published 2026-09-11

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Vyacheslav O. Guba, Alexey V. Reznichenko

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the vast, invisible landscape of quantum physics, particles do not simply sit still; they constantly fluctuate, appearing and disappearing in a blur of probability. To understand how these particles behave, especially when they are subjected to forces like electromagnetism, physicists rely on a powerful mathematical tool known as the heat kernel. Despite its name, which evokes the warmth of a stove, this concept has nothing to do with temperature. Instead, it acts as a master map, describing how a particle's influence spreads out from one point in space to another over a tiny, almost instantaneous moment of time. This map is essential for calculating the energy and behavior of the universe at its most fundamental level, helping scientists predict how fields interact and how the vacuum itself responds to the presence of matter. However, drawing this map is notoriously difficult. While physicists have long known how to sketch the map for a particle starting and ending at the exact same spot, the task becomes exponentially harder when the particle begins at one location and ends at another. This "off-diagonal" scenario, where the start and finish points are separated, has remained a stubborn obstacle, forcing researchers to rely on complicated, step-by-step methods that often break down or become too unwieldy to use for general situations.

A team of researchers from the Budker Institute of Nuclear Physics and Novosibirsk State University has now cleared a path through this difficulty by developing a new, direct method to calculate these complex maps. Instead of building the solution piece by piece through a recursive process that requires solving one step to get to the next, they devised a streamlined approach that computes the entire picture in a single, non-recursive sweep. Their work focuses on a specific type of particle described by the Klein-Gordon equation, which governs how scalar particles move through flat space-time while interacting with electromagnetic fields. By treating the problem as a journey through a mathematical landscape where the start and end points are distinct, the team derived precise, closed-form expressions for the first few layers of the heat kernel's expansion. These expressions act as a detailed blueprint, allowing scientists to see exactly how the particle's behavior changes in the very first moments of its travel, without needing to approximate or simplify the distance between its origin and destination.

The power of this new method lies in its ability to handle the full complexity of the space-time environment. In previous approaches, calculating the interaction between two different points often required expanding the solution into a long, infinite series of terms, which was both tedious and prone to error. The new technique cuts through this by using a clever mathematical representation that keeps the distance between the two points fixed and constant throughout the calculation. This allows the researchers to isolate the specific effects of the electromagnetic field and the potential energy of the space, producing clear formulas for the coefficients that define the particle's behavior. The team successfully applied their method to find the solution up to the second order of time, a level of detail that captures subtle interactions previously difficult to pin down. To ensure their new map was accurate, they tested it against two well-known scenarios: a field that behaves like a plane wave, where the force ripples through space like a light wave, and a field that remains constant and uniform everywhere. In both cases, their new formulas perfectly matched the established, known results, confirming that the method works correctly for both simple and complex field configurations.

This achievement represents a significant refinement in the toolkit available to theoretical physicists. By providing a way to compute these expansions without getting bogged down in recursive loops, the method opens the door to more efficient calculations in quantum field theory. It allows researchers to explore the behavior of particles in arbitrary electromagnetic environments with a clarity that was previously out of reach. The work does not claim to solve every problem in the field, but it offers a robust, verified extension of existing techniques that makes the off-diagonal heat kernel accessible for a wider range of applications. For those studying the fundamental structure of matter and energy, having a reliable, direct way to chart the path of a particle between two distinct points is a crucial step forward, turning a previously tangled mathematical knot into a clear, navigable line.

Technical Summary: A Nonrecursive Method for Computing the Off-Diagonal Small-Time Heat Kernel Expansion

Problem Statement
The heat kernel, defined as the matrix element $h(\tau, x, y) = \langle x|e^{-iH\tau}|y\rangle$ of the operator $e^{-iH\tau}$ , is fundamental to quantum field theory (QFT) and spectral theory. For the Klein-Gordon operator in flat space-time coupled to electromagnetic fields, $H = D_\mu D^\mu + V(x)$ , exact analytical forms of the heat kernel are generally unavailable for arbitrary potentials $V(x)$ and vector potentials $A_\mu(x)$ . Consequently, researchers rely on the small-time asymptotic expansion ( $\tau \to 0$ ):
$h(\tau, x, y) = \frac{ie^{-iz^2/(4\tau)}}{(4\pi i\tau)^{d/2}} \sum_{n=0}^{\infty} a_n(x, y)\tau^n,$
where $z^\mu = x^\mu - y^\mu$ . While numerous methods exist to compute the diagonal coefficients $a_n(x, x)$ (crucial for one-loop renormalization), calculating the off-diagonal coefficients $a_n(x, y)$ is significantly more complex. Existing non-diagonal methods often rely on recursive relations or series expansions in $z^\mu$ , which can be cumbersome or limited in scope.

Methodology
The authors present a nonrecursive method to compute the off-diagonal expansion coefficients $a_n(x, y)$ up to the second order in $\tau$ . The approach extends the nonrecursive technique described by Nepomechie [8], which was previously restricted to diagonal values.

The core of the method involves deriving a specific integral representation for the heat kernel that separates the Gaussian dependence on the proper time $\tau$ from the operator dynamics. The derivation proceeds as follows:

Integral Representation: Starting from the delta function representation, the authors utilize the operator identity $e^{-ikx}D_\mu e^{ikx} = D_\mu + ik_\mu$ to rewrite the heat kernel in terms of an integral over an auxiliary momentum variable $q$ .
Operator Ordering: To handle the non-commutativity of operators within the exponential, the authors introduce an auxiliary parameter $\gamma$ and define an operator $L(\gamma)$ satisfying a differential equation. This allows the use of a path-ordered exponential (P-exp) to disentangle the operators.
Gauge and Shift Handling: The method explicitly accounts for the gauge potential $A_\mu$ and the coordinate shift $x \to x-z$ through the definition of a specific operator $B(\gamma)$ , which absorbs the gauge phase.
Final Representation: The heat kernel is expressed as:
$h(\tau, x, y) = \frac{ie^{-iz^2/(4\tau)}}{(4\pi i\tau)^{d/2}} e^{-iez_\mu \int_0^1 d\theta A^\mu(x-\theta z)} e^{-z\cdot\partial} \mathcal{P} \exp\left\{ -i\tau \left( \int_0^1 d\theta H_x(\theta) - \int_0^1 d\theta_1 D_\mu(\theta_1) \int_0^1 d\theta_2 D^\mu(\theta_2) \right) \right\} 1.$
Here, the operators act on the constant function $\psi(x)=1$ . The expansion in powers of $\tau$ is obtained by expanding the ordered exponential and performing the resulting differential operations.

Key Contributions and Results
The primary contribution of this work is the derivation of closed-form expressions for the expansion coefficients $b_n(x, y)$ (where $a_n$ corresponds to the terms in the expansion of the bracketed term in Eq. 28) up to second order ( $n=0, 1, 2$ ).

Zeroth Order: $b_0(x, y) = 1$ .
First Order ( $b_1$ ): The authors provide a closed-form expression involving integrals of the potential $V$ , the field strength tensor $F_{\mu\nu}$ , and their derivatives along the path connecting $x$ and $y$ . The expression is gauge invariant and explicitly depends on the separation vector $z^\mu$ through path integrals of the field strength (e.g., $B_\mu(s) \equiv z^\alpha \int_0^s d\theta F_{\alpha\mu}$ ).
Second Order ( $b_2$ ): A complex but explicit closed-form expression is derived, containing terms quadratic in the potentials and field strengths, as well as higher-order derivatives. These expressions retain the explicit dependence on $z^\mu$ via the path-ordered integrals of the field strength tensor.

The paper explicitly states that the method is not limited to these orders; subsequent coefficients can be computed via straightforward manipulations of differential operators, though the algebraic complexity increases.

Verification
To validate the derived formulas, the authors applied their general expressions to two specific cases where exact solutions are known:

Plane Wave Fields: For a plane wave background $A_\mu(x) = \epsilon_\mu f(n \cdot x)$ , the authors demonstrated that their general formulas for $b_1$ and $b_2$ reduce exactly to the known small-time asymptotics of the plane wave heat kernel.
Constant Electromagnetic Fields: For constant fields $F_{\mu\nu}$ and constant potential $V=m^2$ , the derived coefficients were shown to reproduce the expansion of the known Schwinger proper-time result for constant fields.

Significance
The paper claims that this work provides a nonrecursive alternative for computing off-diagonal heat kernel coefficients, extending the diagonal-only method of Nepomechie [8] to arbitrary space-time points $x$ and $y$ . The significance lies in providing a systematic, closed-form approach that avoids the recursive relations or series expansions in $z^\mu$ often required by other methods. The authors verify that their results correctly reproduce established asymptotics for specific field configurations, confirming the validity of the general formalism. The method is presented as a tool for obtaining expansion coefficients for general potentials without the need for recursive step-by-step derivation.

More like this