
Recent classicalstatistical numerical simulations have established the
"bottomup" thermalization scenario of Baier et al. as the correct weak
coupling effective theory for thermalization in ultrarelativistic heavyion
collisions. We perform a parametric study of photon production in the various
stages of this bottomup framework to ascertain the relative contribution of
the offequilibrium "Glasma" relative to that of a thermalized QuarkGluon
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 nonequilibrium 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 twocolor QCD out of equilibrium.
We compute the realtime 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 heavyion
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 realtime 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 realtime 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 FockSchwinger 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 nonperturbative
Schwinger mechanism of particle production by quantum tunneling. In some
approximation, we also consider the backreaction of the produced particles on
the external field, as well as the selfinteractions of the produced particles.