I was asked to find conductivity $ \sigma (\omega, T) $ using methods of qft, or more exactly using Matsubara Green's function, for the following system: qed with massless dirac fermions in case $ T \neq 0 $. As I understand, coductivity is a linear response coefficient to applying electric field. $$ j_i (t) = \int_0^{\infty} \sigma_{ij} (\tau, T) E_j (t - \tau) d \tau $$ So I need to use Kubo formula: $$ \sigma (\omega, T) = \frac{i}{\hbar} \int_{0}^{\infty} < [\hat{X}(t), \hat{Y}(0)] > e^{i \omega t} dt $$
I think $ \hat{X} $ should be current operator: $$ \hat{X_{\mu}} = e \bar{\psi} \gamma_{\mu} \psi $$ I can write excitation to electromagnetic field as: $$ \hat{H}_{ext} = i e \bar{\psi} A_{\mu} \gamma_{\mu} \psi $$ But if I use this as $ \hat{Y} $, I will find linear response to the vector potential, not electric field. Maybe, it is still correct and I just need to recalculate conductivity from this linear response coefficient in someway? So that is the first question: what operator I should use as $ \hat{Y} $?
After that, I want to use Wick theorem to express $ <[\hat{X}(t), \hat{Y}(t)]> $ through Matsubara Green's function. Can I use free Green's function $ S(i\omega_n, p) = \frac{1}{i\omega_n \gamma_0 + \bar{p} \bar{\gamma} + i \epsilon} $, or I need to calculate first loop correction?