Abstract: The Colour Glass Condensate effective field theory is a useful framework for studying heavy ion collisions at ultrarelativistic energies, such as those reached at the Large Hadron Collider at CERN. In this framework, we study the correlations of Wilson lines that appear explicitly in cross section expressions. Of particular importance is the evolution of these correlations as a function of rapidity. We use two nonlinear evolution equations to study this behaviour. In the first part of the talk, I will discuss JIMWLK (Jalilian-Marian–Iancu–McLerran–Weigert–Leonidov–Kovner) evolution formulated as a Langevin equation. The Langevin picture is particularly useful for numerical studies and is well suited for studying correlators between particles at unequal rapidities. It also offers an alternative interpretation of evolution as a stochastic process. I will consider inclusive quark-gluon production as a concrete example to i) give a diagrammatic interpretation of the evolution equations that result and ii) go to the dilute limit to study the emergence of BFKL (Balitsky-Fadin-Kuraev-Lipatov) dynamics from this formulation. In the second part of the talk, I will discuss the next-to-leading order BK (Balitsky-Kovchegov) equation, which describes the evolution of 2-point correlators. The 6-point correlators that appear in this equation are typically calculated in the large-Nc limit. I will present a fully analytic calculation of these correlators in the finite Nc case, using the Gaussian Truncation. The results can be used for numerical studies to find the relative importance of finite Nc corrections to the next-to-leading order evolution equation.