After a brief overview of electroweak scale baryogenesis, I present a method for a controlled computation of the CP-violating sources relevant for baryogenesis at a first order phase transition. The method is based on a weak coupling and gradient expansions of the Dyson-Schwinger equations for chiral fermions coupled to bosons and gauge fields. The sources appear both in the flow and collision terms of the kinetic equation. A comparison of our results with source computations in literature will also be carried out. The sources we compute can then be used for baryogenesis computations in various extensions of the Standard Model.