Reconciling Kubo and Keldysh Approaches to Fermi-Sea-Dependent Nonequilibrium Observables: Application to Spin Hall Current and Spin-Orbit Torque in Spintronics

Simao M. Joao,Marko D. Petrovic,J. M. Viana Parente Lopes,Aires Ferreira,Branislav K. Nikolic

2024-08-29

Abstract:Quantum transport studies of spin-dependent phenomena in solids commonly employ the Kubo or Keldysh formulas for the steady-state density matrix in the linear-response regime. Its trace with operators of interest -- such as, spin density, spin current density or spin torque -- gives expectation values of experimentally accessible observables. For such local quantities, these formulas require summing over the manifolds of {\em both} Fermi-surface and Fermi-sea quantum states. However, debates have been raging in the literature about vastly different physics the two formulations can apparently produce, even when applied to the same system. Here, we revisit this problem using a testbed of infinite-size graphene with proximity-induced spin-orbit and magnetic exchange effects. By splitting this system into semi-infinite leads and central active region, in the spirit of Landauer two-terminal setup for quantum transport, we prove the {\em numerically exact equivalence} of the Kubo and Keldysh approaches via the computation of spin Hall current density and spin-orbit torque in both clean and disordered limits. The key to reconciling the two approaches are the numerical frameworks we develop for: ({\em i}) evaluation of Kubo(-Bastin) formula for a system attached to semi-infinite leads, which ensure continuous energy spectrum and evade the need for phenomenological broadening in prior calculations; and ({\em ii}) proper evaluation of Fermi-sea term in the Keldysh approach, which {\em must} include the voltage drop across the central active region even if it is disorder free.

Mesoscale and Nanoscale Physics,Applied Physics,Computational Physics

What problem does this paper attempt to address?

This paper aims to solve the inconsistency problem of the Kubo and Keldysh methods when dealing with non - equilibrium observables that depend on the Fermi sea. Specifically, when these two methods are applied to the same system, they may produce completely different physical results, especially in the research of Spin Hall Current and Spin - Orbit Torque. By using an infinitely large graphene as a test platform and dividing it into semi - infinite leads and a central active region, the author proves the numerical exact equivalence of the Kubo and Keldysh methods in calculating the spin - Hall current density and spin - orbit torque.

### Main problems and solutions
1. **Problem background**:
- **Kubo method**: It is often used in linear response theory to obtain the density matrix in non - equilibrium states by calculating the contributions of the Fermi surface and the Fermi sea.
- **Keldysh method**: It is suitable for the non - linear response region and describes the electron distribution in non - equilibrium states by introducing additional Green functions.

2. **Existing controversies**:
- **Spin Hall Current**: The results obtained by the Kubo method and the Keldysh method in calculating the spin - Hall current are different.
- **Spin - Orbit Torque**: The Kubo method predicts the field - like (odd in magnetization) and damping - like (even in magnetization) spin - orbit torques, while the Keldysh method only predicts the field - like torque in the clean limit.

3. **Solutions**:
- **Improvement of the Kubo method**: A numerical framework has been developed so that the Kubo method can be applied to two - terminal settings, thus avoiding the ambiguity in the finite - size periodic lattice calculation.
- **Improvement of the Keldysh method**: The Fermi - sea term in the Keldysh method has been correctly constructed, especially considering the voltage drop on the central active region.

### Key technical details
- **Kubo method**:
- Use the Kubo - Bastin formula to calculate the contributions of the Fermi surface and the Fermi sea.
- The formulas are as follows:
\[
\hat{\rho}_{\text{neq}}^{\text{Kubo}}=\hat{\rho}_{\text{surf}}^{\text{Kubo}}+\hat{\rho}_{\text{sea}}^{\text{Kubo}}
\]
\[
\hat{\rho}_{\text{surf}}^{\text{Kubo}}=\frac{4eE_x}{\hbar}\text{Re}[\text{Im}\hat{G}\hat{v}_x\text{Im}\hat{G}]
\]
\[
\hat{\rho}_{\text{sea}}^{\text{Kubo}}=\frac{2eE_x}{\hbar}\int dE f(E)\text{Re}[(\hat{G}^r - \hat{G}^a)\hat{v}_x\partial_E\text{Re}\hat{G}]
\]

- **Keldysh method**:
- Use the Keldysh formula to calculate the contributions of the Fermi surface and the Fermi sea.
- The formulas are as follows:
\[
\hat{\rho}_{\text{neq}}^{\text{Keldysh}}=\hat{\rho}_{\text{surf}}^{\text{Keldysh}}+\hat{\rho}_{\text{sea}}^{\text{Keldysh}}
\]
\[
\hat{\rho}_{\text{surf}}