In this highly-cited paper (or its pdf) on photoemission, Eq.(3) gives the current density in terms of $G^<$ $$ {\bf j}({\bf r},t) = 2\hbar\left( \frac{e\hbar}{2m} (\partial_{\bf r'} - \partial_{\bf r}) + \frac{ie^2}{m}{\bf A}({\bf r},t)\right) G^<({\bf r},t;{\bf r'},t)\big|_{{\bf r'}={\bf r}}. $$
It simply follows from Eq.(2) the electromagnetic coupling Hamiltonian $$H_1 = \int d{\bf r} \psi^\dagger({\bf r},t) \left( \frac{ie\hbar}{m}{\bf A}({\bf r},t)\cdot\partial_{\bf r} + \frac{e^2}{2m}{\bf A}^2({\bf r},t)\right)\psi({\bf r},t).$$ I think one can calculate ${\bf j}=\langle{\delta S/\delta {\bf A}\rangle}$ where the action $S$ contains the coupling part $H_1$ and the expectation value leads to some Green's function.
But why is it the lesser one $G^<$? This is not obvious to me.