Recent classical-statistical numerical simulations have established the "bottom-up" thermalization scenario of Baier et al. as the correct weak coupling effective theory for thermalization in ultrarelativistic heavy-ion collisions. We perform a parametric study of photon production in the various stages of this bottom-up framework to ascertain the relative contribution of the off-equilibrium "Glasma" relative to that of a thermalized Quark-Gluon Plasma. Taking into account the constraints imposed by the measured charged hadron multiplicities at RHIC and the LHC, we find that Glasma contributions are important especially for large values of the saturation scale at both energies. These non-equilibrium effects should therefore be taken into account in studies where weak coupling methods are employed to compute photon yields.

We investigate axial charge production in two-color QCD out of equilibrium. We compute the real-time evolution starting with spatially homogeneous strong gauge fields, while the fermions are in vacuum. The idealized class of initial conditions is motivated by glasma flux tubes in the context of heavy-ion collisions. We focus on axial charge production at early times, where important aspects of the anomalous dynamics can be derived analytically. This is compared to real-time lattice simulations. Quark production at early times leading to anomalous charge generation is investigated using Wilson fermions. Our results indicate that coherent gauge fields can transiently produce significant amounts of axial charge density, while part of the induced charges persist to be present even well beyond characteristic decoherence times. The comparisons to analytic results provide stringent tests of real-time representations of the axial anomaly on the lattice.

We revisit the problem of quark production in high energy heavy ion collisions, at leading order in $\alpha_s$ in the color glass condensate framework. In this first paper, we setup the formalism and express the quark spectrum in terms of a set of mode functions of the Dirac equation. We determine analytically their initial value in the Fock-Schwinger gauge on a proper time surface $Q_s\tau_0\ll 1$, in a basis which makes manifest the boost invariance properties of this problem. We also describe a statistical algorithm to perform the sampling of the mode functions.

In this paper, we show how classical statistical field theory techniques can be used to efficiently perform the numerical evaluation of the non-perturbative Schwinger mechanism of particle production by quantum tunneling. In some approximation, we also consider the back-reaction of the produced particles on the external field, as well as the self-interactions of the produced particles.