Schwinger Keldysh field theory is a widely used paradigm to study nonequilibrium dynamics of quantum many-body systems starting from a thermal state. We extend this formalism to describe nonequilibrium dynamics of quantum systems starting from arbitrary initial many-body density matrices. We show how this can be done for both bosons and fermions, and for both closed and open quantum systems, using additional sources coupled to bilinears of the fields at the initial time, calculating Green's functions in a theory with these sources, and then taking appropriate set of derivatives of these Green's functions with respect to initial sources to obtain physical observables. The set of derivatives depends on the initial density matrix. The physical correlators in a dynamics with arbitrary initial conditions do not satisfy Wick's theorem, even for noninteracting systems. However, our formalism constructs intermediate ``$n\ensuremath{-}\mathrm{particle}$ Green's functions'' that obey Wick's theorem and provide a prescription to obtain physical correlation functions from them. This allows us to obtain analytic answers for all physical many-body correlation functions of a noninteracting system even when it is initialized to an arbitrary density matrix. We use these exact expressions to obtain an estimate of the violation of Wick's theorem and relate it to presence of connected multiparticle initial correlations in the system. We illustrate this new formalism by calculating density and current profiles in many-body fermionic and bosonic open quantum systems initialized to nontrivial density matrices. We have also shown how this formalism can be extended to interacting many-body systems.
We extend the self-consistent Green's functions formalism to take into account three-body interactions. We analyze the perturbative expansion in terms of Feynman diagrams and define effective one- and two-body interactions, which allows for a substantial reduction of the number of diagrams. The procedure can be taken as a generalization of the normal ordering of the Hamiltonian to fully correlated density matrices. We give examples up to third order in perturbation theory. To define nonperturbative approximations, we extend the equation-of-motion method in the presence of three-body interactions. We propose schemes that can provide nonperturbative resummation of three-body interactions. We also discuss two different extensions of the Koltun sum rule to compute the ground state of a many-body system.