С заданным фермионным заряженным полем.
С ЗАДАННЫМ, т.е. с коммутирующим (в гайзенберговском представлении!) с ЭМ полем -- это и есть с классическим током. Если фермионное поле не зависит от электромагнитного, с ним коммутирует, то все просто, заменяем ток, выраженный через фермионное поле, просто буквой

и приходим к предыдущему.
Собственно, можно вообще ничего не делать. Просто понимать, что фотонный пропагатор -- это функция Грина (с смысле диффуравнений), выражающая поле через ток. Ну а перейти от интегрального уравнения к дифференциальному... Подействовали даламберианом и бац,

в знаменателе как не бывало

С пропагатором в фейнмановской калибровке так вообще все сразу.
В общем тут есть разные варианты. Можно просто убедиться, что оператор 4-потенциала тождественно удовлетворяет уравнениям Максвелла (но это без тока). Но квантовать удобнее свободное поле, а взаимодействие с током включить потом, в гамильтониан. Тогда понадобится коммутатор. И вообще так понятнее физически (классический ток порождает поле в глауберовском состоянии, тогда не ноль среднее поле и проч. бла-бла-бла). Можно сразу написать операторные выражения для

и

и их усреднять, так, пожалуй, даже проще.