Second-order dc conductivity in the velocity-gauge Keldysh formalism: gauge-invariant decomposition into nonlinear Drude, Berry-curvature-dipole, and quantum-metric responses

Abstract: We derive a gauge-invariant clean-limit decomposition of the second-order dc nonlinear conductivity in multiband tight-binding systems within the velocity-gauge Keldysh Green's function formalism. In the constant-relaxation-time approximation, the dc response separates into four contributions with distinct lifetime $τ$ scalings and physical origins: the nonlinear Drude term $σ^{\mathrm{ND}}_{ijk}\proptoτ^{2}$, the Berry-curvature-dipole term $σ^{\mathrm{BCD}}_{ijk}\proptoτ$, the intraband quantum-metric-dipole term $σ^{\mathrm{intra\text{-}QMD}}_{ijk}\proptoτ^{0}$, and the interband quantum-metric-dipole term $σ^{\mathrm{inter\text{-}QMD}}_{ijk}\proptoτ^{0}$. The intraband term is a Fermi-surface dipole of the ordinary band quantum metric, while the interband term is written, in the present representation, as a Fermi-sea-type response involving a band-normalized quantum metric. Working entirely within the velocity-gauge Keldysh--Kubo framework, we show that all connection-dependent commutator terms generated in the band-basis expansion cancel exactly between the covariant-quantum-connection sector $σ^{\mathcal{C}}_{ijk}$ and the three-Berry-connection sector $σ^{\mathcal{T}}_{ijk}$, making the role of the Peierls contact velocity vertices $V_{ij}$ and $V_{ijk}$ explicit; a complementary projector-based derivation appears in Ulrich et al., Phys. Rev. B 113, L201107 (2026), and our Fermi-surface dc-limit expression agrees with that reference after accounting for index and convention differences. As a diagnostic illustration, we introduce a real two-band model in which the Berry curvature and hence the BCD response vanish identically while the intraband quantum-metric dipole remains finite, establishing a practical route to quantum-metric dc responses not reducible to the Berry-curvature-dipole mechanism.

• The paper rigorously derives a gauge-invariant formulation for second-order DC conductivity in multiband systems employing the velocity-gauge Keldysh formalism.
• It decomposes the nonlinear current into distinct contributions—nonlinear Drude (ND), Berry curvature dipole (BCD), and quantum metric dipole (QMD) terms—each with unique lifetime scaling.
• The methodology offers clear diagnostics for separating geometric and topological effects in nonlinear transport, paving the way for both theoretical and experimental advances.

Gauge-Invariant Decomposition of Second-Order DC Conductivity in the Velocity-Gauge Keldysh Formalism

Introduction and Context

This work provides a comprehensive derivation of the second-order dc nonlinear conductivity in multiband tight-binding systems, employing the velocity-gauge Keldysh Green’s function formalism. Beyond offering gauge invariance, the approach yields a clean decomposition of the response into physically meaningful sectors: nonlinear Drude (ND), Berry curvature dipole (BCD), and two distinct quantum-metric-dipole (QMD) contributions—one intraband (intra-QMD) and one interband (inter-QMD). These contributions are organized according to their scaling with the electronic quasiparticle lifetime $\tau$: respectively, $\mathcal{O}(\tau^2)$, $\mathcal{O}(\tau)$, and $\mathcal{O}(\tau^0)$ for the ND, BCD, and QMD terms.

The paper situates its contribution within the ongoing development of geometric response theory, where both the Berry curvature and the quantum metric of Bloch bands act as geometric sources of transport phenomena. While the BCD mechanism has been established as a leading nonlinear source term, quantum metrics have more recently been recognized as distinct contributors to nonlinear responses, especially in systems with vanishing Berry curvature.

Formalism: Peierls Substitution and the Velocity-Gauge Keldysh Framework

The model considers a multiband tight-binding Hamiltonian, coupled to uniform electromagnetic fields via the Peierls substitution, $$H_{\bm{k}} \rightarrow H_{\bm{k} + (e/\hbar)\bm{A}(t)}$$, in the velocity gauge. Expanding the Hamiltonian in powers of the vector potential $\bm{A}$, this yields first, second, and third order velocity-vertex operators ($V_i$, $V_{ij}$, $V_{ijk}$) that encode the system’s band structure via momentum derivatives.

The current operator derives from the functional derivative of the Hamiltonian with respect to $\bm{A}$. Crucially, higher-order contact (diamagnetic) current vertices arising from the Peierls expansion (such as $\mathcal{O}(\tau^2)$0 and $\mathcal{O}(\tau^2)$1) are necessary to enforce gauge invariance and remove spurious divergences in the dc limit.

Figure 1: Diagrammatic representation of the four second-order current contributions in the Peierls velocity gauge: current vertices and Peierls contact vertices underpin distinct contributions.

All nonlinear response calculations are formulated within the Keldysh nonequilibrium Green’s function formalism, providing a systematic perturbative expansion in $\mathcal{O}(\tau^2)$2.

Second-Order DC Conductivity: Lifetime Hierarchy and Geometric Decomposition

Within the constant-relaxation-time approximation, the second-order nonlinear current separates into diagrammatically distinct terms corresponding to different sequences of vertex insertions. After reducing to the DC limit and eliminating apparent $\mathcal{O}(\tau^2)$3 singularities using gauge symmetry and the Ward–Takahashi identity, the response admits a universal decomposition: $\mathcal{O}(\tau^2)$4 Each term is associated with distinct physical origins and $\mathcal{O}(\tau^2)$5 dependences.

Nonlinear Drude (ND) Term ($\mathcal{O}(\tau^2)$6):

Originates entirely from single-band intraband processes and depends only on band dispersion. It corresponds to the conventional semiclassical Drude picture, generalized to nonlinear order.

Berry Curvature Dipole (BCD) Term ($\mathcal{O}(\tau^2)$7):

Arises from the imaginary part of the quantum geometric tensor (the Berry curvature), combined with a pinch singularity in the Keldysh kernel. This is the clean metal counterpart to the nonlinear Hall effect associated with Berry curvature dipoles [Sodemann & Fu PRL 2015, (2606.22359)].

Quantum Metric Dipole (QMD) Terms ($\mathcal{O}(\tau^2)$8):

Intrinsic "geometric" contributions controlled by the quantum metric. The intra-QMD term is explicitly a Fermi-surface dipole of the quantum metric, while the inter-QMD is a Fermi-sea–weighted contribution involving derivatives of a band-normalized metric. They persist even if the BCD vanishes due to symmetries (e.g., in PT-symmetric systems or real two-band models).

Figure 2: Two-model summary of band structures and nonlinear geometric diagnostics, contrasting BCD-dominated (Dirac) and QMD-dominated (real) two-band models.

The entire analysis is performed directly in the velocity gauge, making explicit how Peierls contact vertices ensure gauge invariance at each diagrammatic level.

Gauge Invariance and Cancellation of Connection-Dependent Terms

A central technical achievement is the explicit demonstration of cancellation of gauge-dependent "connection" terms (covariant derivatives of Berry connections and three-Berry-connection loops) in the final physical observable. Although such terms arise naturally in the band-basis expansion of the Keldysh-Kubo response (and have sometimes been individually attributed physical meaning), the physical current is shown to depend only on gauge-invariant band geometric quantities (Berry curvature for BCD, quantum metric for QMD), consistent with the projector calculus developed in [Mitscherling et al. 2025; Ulrich et al. 2026].

This cancellation reflects the broader geometric principle asserted in the literature: only projectors and their covariant derivatives enter true physical quantities, not arbitrary choices of Bloch-state phases or Berry connections.

Model Diagnostics: Distinguishing Geometric and Berry Contributions

Minimal two-band models are used to clarify the distinct origins of BCD and QMD responses:

• Tilted Massive Dirac Model (BCD-dominated): Exhibits finite BCD and QMD responses, with the former dominating at high $\mathcal{O}(\tau^2)$9. The geometric origin can be mapped directly onto Berry curvature gradients.
• Real Two-Band Model (QMD-dominated): By constructing a model with real Bloch functions (and thus identically zero Berry curvature), BCD is forced to vanish, and nonlinear response arises purely from the quantum metric.

  Figure 3: Momentum-resolved intraband-QMD integrand $\mathcal{O}(\tau)$0 in the Dirac and real models, contrasting BCD and QMD contributions at the Fermi surface.

These examples provide practical diagnostic tools for distinguishing between nonlinear responses of geometric (QMD) and topological (BCD) origin, and are relevant for experimental interpretation—such as in topological antiferromagnets where QMD signatures are isolated [Wang et al. Nature 2023].

Theoretical and Practical Implications

The explicit decomposition into ND, BCD, intra-QMD, and inter-QMD terms has several important implications:

• It provides a rigorous basis for interpreting nonlinear dc transport experiments in terms of band geometry, enabling extrication of quantum-metric contributions in the absence of Berry curvature.
• The lifetime scaling clarifies the roles of disorder and dissipation, distinguishing between extrinsic (disorder-dominated) and intrinsic (geometry-dominated) regimes.
• The gauge-invariant approach sidesteps ambiguities of connection-dependent quantities that have complicated prior analyses and poses a template for formulating observables in higher-order nonlinear response.
• The methods and decompositions here generalize naturally to the formulation of third-order and higher nonlinear responses in crystals, with expected analogs in the emergence of higher quantum-geometric tensors and multipoles [Mandal et al. PRB 2024, Zhang et al. PRB 2023].

Outlook and Prospects

This unified framework provides a foundation for future theoretical investigations and materials discovery targeting nonlinear geometric transport phenomena. Systematic inclusion of disorder, explicit vertex corrections, and generalization beyond nondegenerate bands are natural directions. The velocity-gauge approach also suggests numerical strategies for first-principles calculations of geometric response tensors in complex materials. Experimental efforts can leverage these diagnostics to design or identify materials—in particular PT-symmetric or antiferromagnetic compounds—where nonlinear quantum-metric responses can be isolated and maximized.

Conclusion

This work delivers a rigorous, gauge-invariant decomposition of the second-order dc conductivity for multiband tight-binding systems, articulated within the velocity-gauge Keldysh formalism. The transport response is shown to consist of nonlinear Drude, Berry-curvature-dipole, and both intraband and interband quantum-metric-dipole contributions, each with well-defined physical origin, gauge structure, and disorder scaling. The methodology provides both a conceptual and computational blueprint for dissecting and predicting nonlinear transport phenomena in topological and quantum-geometric materials.

Figure 4: Diagrammatic representation of the zeroth- and first-order current contributions, emphasizing the necessity of Peierls contact vertices for gauge invariance.

Figure 5: Diagrammatic representation of the four second-order current contributions, further underscoring the interplay of paramagnetic and contact terms in DC response.

---

References:

The technical results and further discussion can be found in "Second-order dc conductivity in the velocity-gauge Keldysh formalism: gauge-invariant decomposition into nonlinear Drude, Berry-curvature-dipole, and quantum-metric responses" (2606.22359).