rho_epsilon^1(t) = sum q e for integral R^2 j m_epsilon'(t-tau mu) (del p(epsilon,epsilon')) (j,tau)