Given a polynomially bounded
functional F over the field configurations, then, for any
state vector (which is a solution of the QFT), |\psi\rangle, we have :\left\langle\psi\left|\mathcal{T}\left\{\frac{\delta}{\delta\varphi}F[\varphi]\right\}\right|\psi\right\rangle = -i\left\langle\psi\left|\mathcal{T}\left\{F[\varphi]\frac{\delta}{\delta\varphi}S[\varphi]\right\}\right|\psi\right\rangle where \frac{\delta}{\delta \varphi} is the
functional derivative with respect to \varphi, S is the
action functional and \mathcal{T} is the
time ordering operation. Equivalently, in the
density state formulation, for any (valid) density state \rho, we have :\rho\left(\mathcal{T}\left\{\frac{\delta}{\delta\varphi}F[\varphi]\right\}\right) = -i\rho\left(\mathcal{T}\left\{ F[\varphi] \frac{\delta}{\delta\varphi}S[\varphi]\right\}\right). This
infinite set of equations can be used to solve for the correlation functions
nonperturbatively. To make the connection to diagrammatic techniques (like
Feynman diagrams) clearer, it is often convenient to split the action S as :S[\varphi] = \tfrac{1}{2}\varphi^{i}D^{-1}_{ij}\varphi^{j} + S_{\text{int}}[\varphi], where the first term is the quadratic part and D^{-1} is an invertible symmetric (antisymmetric for fermions) covariant tensor of rank two in the
deWitt notation whose inverse, D is called the bare propagator and S_{\text{int}}[\varphi] is the "interaction action". Then, we can rewrite the SD equations as :\left\langle\psi\left\vert\mathcal{T}\left\{F \varphi^j\right\}\right\vert\psi\right\rangle=\left\langle\psi\left\vert\mathcal{T}\left\{iF_{,i}D^{ij}-FS_{\text{int},i}D^{ij}\right\}\right\vert\psi\right\rangle. If F is a functional of \varphi, then for an
operator K, F[K] is defined to be the operator which substitutes K for \varphi. For example, if :F[\varphi]=\frac{\partial^{k_1}}{\partial x_1^{k_1}}\varphi(x_1)\cdots \frac{\partial^{k_n}}{\partial x_n^{k_n}}\varphi(x_n) and G is a functional of J, then :F\left[-i\frac{\delta}{\delta J}\right]G[J]=(-i)^n \frac{\partial^{k_1}}{\partial x_1^{k_1}}\frac{\delta}{\delta J(x_1)} \cdots \frac{\partial^{k_n}}{\partial x_n^{k_n}}\frac{\delta}{\delta J(x_n)} G[J]. If we have an "
analytic" (a function that is locally given by a convergent power series)
functional Z (called the
generating functional) of J (called the
source field) satisfying :\frac{\delta^n Z}{\delta J(x_1) \cdots \delta J(x_n)}[0]=i^n Z[0] \langle\varphi(x_1)\cdots \varphi(x_n)\rangle, then, from the properties of the functional integrals :{\left \langle \frac{\delta \mathcal{S}}{\delta \varphi(x)}\left[\varphi \right] + J(x)\right\rangle}_J=0, the Schwinger–Dyson equation for the generating functional is :\frac{\delta S}{\delta \varphi(x)}\left[-i \frac{\delta}{\delta J} \right] Z[J] + J(x)Z[J]=0. If we expand this equation as a
Taylor series about J = 0, we get the entire set of Schwinger–Dyson equations. ==An example:
φ4==