Consider a quantum field theory with action I[ϕ], and suppose we 'turn on' a set of classical currents Jr(x) coupled to the fields ϕr(x) of the theory. The complete vacuum vacuum amplitude in the presence of these currents is then
Z[J]≡⟨VAC, out|VAC, in⟩J=∫[∏s,ydϕs(y)]exp(iI[ϕ]+i∫d4xϕr(x)Jr(x)+ϵ terms).
The Feynman rules for calculating Z[J] are just the same as for calculating the vacuum-vacuum amplitude Z[0] in the absence of the external current, except that the Feynman diagrams now contain vertices of a new kind, to which a single ϕr-line is attached. Such a vertex labelled with a coordinate x contributes a 'coupling' factor iJr(x) to the integrand of the position-space Feynman amplitude. Equivalently, we could say that in the expansion of Z[J] in powers of J, the coefficient of the term proportional to iJr(x)iJs(y)⋯ is just the sum of diagrams with external lines (including propagators) corresponding to the fields ϕr(x),ϕs(y,) etc. In particular, the first derivative gives the vacuum matrix element of the quantum mechanical operator Φr(x) corresponding to ϕr(x):
[δδJr(y)Z[J]]J=0=∫[∏r,xdϕr(x)]ϕr(y)exp{iI[ϕ]+ϵ terms}=i⟨VAC, out|Φr(x)|VAC, in⟩J=0.
What's the difference between quantum mechanical operator Φr(x) and its corresponding quantum fields ϕr(x) ?