🔬 mesoscale physics

Unified contact-free formulation of linear-response transport theory

This paper introduces a unified, contact-free moment formulation of linear-response transport theory that utilizes scalar sources to calculate transport coefficients directly from retarded moment-moment correlators, thereby eliminating the need for separate magnetization corrections and providing a gauge-invariant framework applicable to particle, thermal, and spin transport.

Original authors: Ziqian Wang, Ji Feng

Published 2026-09-15

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Ziqian Wang, Ji Feng

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 hidden world of quantum materials, scientists are constantly trying to understand how heat and energy move through solids. While we can easily measure how electricity flows through a wire, tracking the flow of heat is far more elusive, especially when that heat moves sideways, perpendicular to the temperature difference. This sideways movement, known as the thermal Hall effect, has become a crucial tool for diagnosing the invisible geometric shapes that govern electrons and magnetic waves inside a material. For decades, the standard way to calculate this effect has been a two-step process. First, researchers calculate the direct flow of heat. Then, they must perform a separate, tricky correction to subtract out the "magnetization" of heat—a swirling, internal circulation of energy that stays trapped inside the material and never actually travels from one side to the other. This internal swirl is not a transport current, but because it is mathematically entangled with the real flow, it has to be carefully peeled away to get the right answer. This extra step has made the theory complex and dependent on arbitrary choices about how to define the local energy within the material.

A team of researchers at Peking University and the Hefei National Laboratory has now proposed a cleaner, unified way to handle these calculations. They developed a new mathematical framework that eliminates the need to subtract the internal heat swirl separately. Instead of treating the circulating heat as a nuisance to be removed later, their method builds the correction directly into the core calculation from the very beginning. By focusing on the "moments" of energy and spin—essentially measuring how these quantities are distributed across the material rather than just their total amount—they found a way to formulate the transport laws using only simple, scalar forces. In this new approach, the internal circulation and the real transport current are naturally balanced against each other before the final answer is taken, leaving a result that is free of the messy corrections that previously cluttered the theory.

The researchers demonstrated that this new method works by applying it to magnetic insulators, where heat is carried by collective waves of spin called magnons. In these materials, a temperature gradient can push these waves sideways, creating a thermal Hall effect, or even generate a sideways flow of pure spin without any net heat flow, known as the spin Nernst effect. Using their unified framework, the team showed that both of these phenomena emerge naturally from the same underlying geometric structure of the material's energy bands. The calculation automatically isolates the finite, measurable signal from the infinite, non-transporting background. This means that whether one is studying the flow of particles, heat, or spin, the same set of rules applies, and the result is independent of how the researcher chooses to assign energy to specific points in the crystal lattice.

The power of this new formulation lies in its ability to simplify the organization of complex transport problems. In the past, calculating the thermal Hall effect required introducing a fictitious gravitational field to mimic the temperature gradient, and then introducing a second, auxiliary field to account for the heat magnetization. The new method requires only the first, simple field. It treats the "contact" terms—the instantaneous interactions that usually require separate subtraction—as part of the natural response of the system. This removes the ambiguity of choosing a specific "gauge" or reference frame for the local current, ensuring that the final transport coefficient is a robust, physical observable that does not change based on mathematical conventions.

To test their theory, the authors applied it to a specific model of magnetic waves, known as Landau-Lifshitz spin waves. In this scenario, the mathematics revealed that the long-time behavior of the system's correlations naturally separates into two parts: a part that grows steadily over time, representing the true transport, and a part that remains bounded or grows too slowly to matter for the final measurement. The new formula acts as a filter, automatically discarding the non-transporting parts and keeping only the finite residue that corresponds to the actual Hall coefficient. When they calculated the thermal Hall conductivity and the spin Nernst coefficient for this model, their results matched the established, accepted formulas perfectly. This agreement serves as a strong validation, proving that their unified, contact-free approach captures the correct physics without the need for the traditional, cumbersome subtraction steps.

This work suggests a path forward for studying more complex, interacting systems where the simple models of non-interacting waves break down. Because the new method relies on the long-time behavior of correlations rather than on constructing specific local magnetization currents, it is well-suited for use in advanced numerical simulations of real, messy materials. It offers a way to calculate transport properties in systems with strong interactions or spin-orbit coupling, where the separation of current and magnetization has historically been a major hurdle. By providing a single, consistent operator framework for particle, thermal, and spin transport, the researchers have streamlined the theoretical toolkit available to physicists. The result is a clearer, more direct view of how quantum geometry drives the flow of energy and information through the solid state, stripping away the mathematical clutter to reveal the fundamental transport mechanisms underneath.

Technical Summary: Unified Contact-Free Formulation of Linear-Response Transport Theory

Problem Statement
Conventional linear-response theory for thermal Hall transport relies on a formulation that combines a current–current Kubo response with an energy-magnetization correction. This correction is typically derived by introducing an auxiliary pseudogravitational magnetic field, a method established by Cooper et al. and refined by Qin et al. While physically correct, this approach introduces a second fictitious probe and separates the final transport observable into two distinct parts: a Kubo term and a magnetization term. A significant drawback of this conventional route is that both individual pieces depend on the specific choice of local energy and current representatives (the "local-current gauge"), even though only their sum constitutes a physical transport observable. Furthermore, this separation complicates the unified treatment of particle, thermal, and spin transport channels.

Methodology
The authors develop a unified, contact-free moment formulation for linear response to scalar sources. The core of the methodology involves:

Definition of Density Moments: Instead of working directly with local currents or densities, the authors define finite- $q$ moments of density operators ( $\hat{A}_{\alpha}$ ) and source operators ( $\hat{\Sigma}^a_{\alpha}$ ). These moments are related to integrated currents via continuity equations, specifically $\hat{J}^a_\alpha = \dot{\hat{A}}_\alpha - \hat{\Sigma}^a_\alpha$ . This formulation avoids ill-defined position operators under periodic boundary conditions.
Scalar Source Response: The response is calculated using only scalar sources (e.g., Luttinger's gravitational potential for thermal transport, Zeeman fields for spin). The perturbation to the grand Hamiltonian is expressed in terms of these scalar sources and their conjugate moments.
Contact-Free Derivation: The linear-response coefficient is derived by analyzing the retarded moment–moment correlator. A key step involves an integration by parts of the susceptibility, which generates an equal-time endpoint term (the commutator of moments). The authors demonstrate that this endpoint term exactly cancels the explicit "contact correction" arising from the field dependence of the current operator before the DC limit ( $\omega \to 0$ ) is taken.
Low-Frequency Extraction: The transport coefficient is identified as the low-frequency residue of the retarded moment–moment correlator. Specifically, for channels where source contributions vanish in the DC limit, the coefficient is proportional to $\lim_{\omega \to 0} (\omega^+) \chi_{A_\alpha, B_\beta}(\omega)$ . This isolates the finite DC residue from the ballistic ( $t^2$ ) sector, which is regularized by the finite- $q$ definition.

Key Contributions

Unified Framework: The formulation organizes particle, grand-energy, and specified spin-density channels within a single source–moment response matrix. It treats thermal and spin transport on equal footing without requiring separate magnetization subtractions.
Gauge and Redundancy Independence: The resulting transport coefficient is:
- Independent of the local-current gauge (transverse circulating currents).
- Invariant under admissible redistributions of bond or interaction energy (density improvements), provided these redistributions do not overlap with conserved polarization modes.
- Formulated without an auxiliary pseudogravitational magnetic field.
Contact-Free Operator Formalism: The method incorporates the equal-time endpoint contributions (which contain the magnetization physics) directly into the moment response, eliminating the need for explicit additive magnetization terms in the final formula.

Results
The authors apply this formalism to quadratic Landau–Lifshitz spin waves to demonstrate its efficacy:

Magnon Thermal Hall: By using the grand-energy moment in the response channel, the authors derive the magnon thermal Hall conductivity. The calculation isolates the finite DC residue from the ballistic sector, yielding the standard result involving the Berry curvature of the magnon bands and Bose thermal weights.
Magnon Spin Nernst: By replacing the grand-energy moment with a longitudinal spin-density moment (while retaining the thermal source), the authors derive the magnon spin Nernst coefficient.
Geometric Structure: Both results emerge from the same paraunitary band geometry, confirming that different transported densities are handled within a single unified framework. The secular growth of polarization-moment correlations is shown to isolate the finite Hall residue, while regular contributions are annihilated by the $\omega^2$ prefactor in the low-frequency limit.

Significance and Claims
The paper claims to provide a "unified operator framework" that simplifies the organization of particle, thermal, and spin transport. Its primary significance lies in streamlining calculations by eliminating the need for separate contact and magnetization corrections and by removing the dependence on specific local-current gauges.

The authors emphasize that their formulation does not discard heat magnetization; rather, it incorporates the corresponding endpoint contribution naturally within the scalar-source Kubo problem. They note that the operator formulation itself does not rely on quadratic or non-interacting approximations, suggesting that the method is applicable to interacting systems, provided one can evaluate the long-time behavior of polarization-moment correlators. The work serves as a theoretical unification that recovers established formulas for magnon thermal Hall and spin Nernst effects while offering a more robust and gauge-invariant derivation.

More like this