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

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.