In this chapter, we introduce and explain one of the fundamental tools in the study of dynamic critical phenomena, namely dynamic perturbation theory. Inevitably, large parts of Chapter 4 need to be rather technical. Complementary to the straightforward iterative method for the solution of non-linear Langevin equations presented in Section 3.2.3, we now describe the more elegant and efficient field-theoretic techniques. Yet both for the derivation of general properties of the perturbation series and for the sake of practical calculations, the response functional and the following elaborations on Feynman diagrams, cumulants, and vertex functions prove indispensable. As is the case with every efficient formalism, once one has become acquainted with the Janssen–De Dominicis response functional for the construction of dynamic field theory, and therefrom the perturbation expansion in terms of vertex functions, it serves to save a considerable amount of rather tedious work. For the purely relaxational O(n)-symmetric models A and B, we derive the fluctuation-dissipation theorem within this formalism, and then systematically construct the perturbation series and its diagrammatic representation in terms of Feynman graphs. In order to reduce the required efforts to a minimum, the generating functionals first for the cumulants, and then, motivated by Dyson's equation for the propagator, for the one-particle irreducible vertex functions are examined. As an example, the relevant vertex functions for models A and B are explicitly evaluated to two-loop order. Furthermore, alternative formulations of the perturbation expansion are discussed, and the Feynman rules are given in both the frequency and the time domain. The results of this chapter provide the foundation for the renormalization group treatment in the subsequent Chapter 5.