# Structure of longitudinal chromomagnetic fields in high energy collisions

###### Abstract

We compute expectation values of spatial Wilson loops in the forward light cone of high-energy collisions. We consider ensembles of gauge field configurations generated from a classical Gaussian effective action as well as solutions of high-energy renormalization group evolution with fixed and running coupling. The initial fields correspond to a color field condensate exhibiting domain-like structure over distance scales of order the saturation scale. At later times universal scaling emerges at large distances for all ensembles, with a nontrivial critical exponent. Finally, we compare the results for the Wilson loop to the two-point correlator of magnetic fields.

###### pacs:

24.85.+p,25.75.-q,12.38.Mh, 12.38.Lg^{†}

^{†}preprint:

## I Introduction

Heavy ion collisions at high energies involve non-linear dynamics of strong QCD color fields (1). These soft fields correspond to gluons with light-cone momentum fractions x\ll 1, which can be described in the “Color Glass Condensate” (CGC) framework. Because of the high gluon occupation number the gluon field can be determined from the classical Yang-Mills equations with a static current on the light cone (2); *McLerran:1994ka; *McLerran:1994vd. It consists of gluons with a transverse momentum on the order of the density of valence charges per unit transverse area, Q_{\mathrm{s}}^{2} (5). Parametrically, the saturation momentum scale Q_{\mathrm{s}} separates the regime of non-linear color field interactions from the perturbative (linear) regime. It is commonly defined using a two-point function of electric Wilson lines, the “dipole scattering amplitude” evaluated in the field of a single hadron or nucleus (6) as described below.

Before the collision the individual fields of projectile and target are two dimensional pure gauges; in light cone gauge,

\alpha^{i}_{m}=\frac{i}{g}\,V_{m}\,\partial^{i}V_{m}^{\dagger} | (1) |

where m=1,\,2 labels the projectile and target, respectively. Here V_{m} are light-like SU({N_{\mathrm{c}}}) Wilson lines, which correspond to the eikonal phase of a high energy projectile passing through the classical field shockwave (7); (8).

The field in the forward light cone after the collision up to the formation of a thermalized plasma is commonly called the “glasma” (9). Immediately after the collision longitudinal chromo-electric and magnetic fields E_{z},~{}B_{z}\sim 1/g dominate (10); *Fries:2006pv; (9). They fluctuate according to the random local color charge densities of the valence sources. The magnitude of the color charge fluctuations is related to the saturation scale Q_{\mathrm{s}}^{2}. The transverse gauge potential at proper time \tau\equiv\sqrt{t^{2}-z^{2}}\to 0, is given by (12)

A^{i}=\alpha_{1}^{i}+\alpha_{2}^{i}~{}. | (2) |

Note that while the fields of the individual projectiles \alpha_{m}^{i} are pure gauges, for a non-Abelian gauge theory A^{i} is not. Hence, spatial Wilson loops evaluated in the field A^{i} are not equal to 1. The field at later times is then obtained from the classical Yang-Mills equations of motion, which can be solved either analytically in an expansion in the field strength (12); (13) or numerically on a lattice (14); (15); *Krasnitz:2003jw; *Lappi:2003bi; (18). The Wilson loop, and the magnetic field correlator, provide an explicitly gauge-invariant method to study the nonperturbative dynamics of these fields, complementary to studies of the gluon spectrum (19).

Spatial Wilson loops at very early times \tau have recently been studied numerically in Ref. (20), using the MV model (2); *McLerran:1994ka; *McLerran:1994vd for the colliding color charge sheets. It was observed that the loops effectively satisfy area law scaling for radii \,\,\vbox{\hbox{$\buildrel\displaystyle>\over{\sim}$}}\,\,1/Q_{\mathrm{s}}, up to a few times this scale. Furthermore, Ref. (21) found that two-point correlators of B_{z} over distances \lesssim 1/Q_{\mathrm{s}} correspond to two dimensional screened propagators with a magnetic screening mass a few times Q_{\mathrm{s}}. This indicates that the initial fields exhibit structure such that magnetic flux does not spread uniformly over the transverse plane (like in a Coulomb phase) but instead is concentrated in small domains.

The present paper extends this previous work as follows. We perform lattice measurements of spatial Wilson loops over a much broader range of radii to analyze their behavior at short (R\ll 1/Q_{\mathrm{s}}) and long (R\gg 1/Q_{\mathrm{s}}) distances. We also implement the so-called JIMWLK (5); (22); *Jalilian-Marian:1997gr; *Jalilian-Marian:1997dw; *JalilianMarian:1998cb; *Iancu:2000hn; *Iancu:2001md; *Ferreiro:2001qy; *Iancu:2001ad; (30); *Mueller:2001uk high-energy functional renormalization group evolution which resums observables to all orders in \alpha_{\mathrm{s}}\log(1/x). High-energy evolution modifies the classical ensemble of gauge field configurations (4), (5) to account for nearly boost invariant quantum fluctuations at rapidities far from the sources. Finally, we also solve the Yang-Mills equations in the forward light cone to study the time evolution of magnetic flux loops.

The calculation of the initial conditions and the numerical solution
of the classical boost-invariant^{1}^{1}The YM equations are solved
in terms of the coordinates \tau=\sqrt{t^{2}-z^{2}},
\eta=\frac{1}{2}\ln\frac{t+z}{t-z} and {\mathbf{x}_{T}}; hence \,\mathrm{d}s^{2}=\,\mathrm{d}\tau^{2}-\tau^{2}\,\mathrm{d}\eta^{2}-\,\mathrm{%
d}{\mathbf{x}_{T}}^{2}. Yang-Mills fields in the
initial stages of a heavy ion collision have been documented in
the references given below, so here we will only describe
them very briefly in Sec. II before moving on
to show our results in Secs. III and IV.

## II Lattice implementation

We work on a two dimensional square lattice of N_{\perp}^{2} points with periodic boundary conditions and consider color sources that fill the whole transverse plane. The lattice spacing is denoted as a, thus the area of the lattice in physical units is L^{2}=N_{\perp}^{2}a^{2}. The calculations are performed for {N_{\mathrm{c}}}=3 colors. In this work we only consider symmetric collisions, where the color charges of both colliding nuclei are taken from the same probability distribution.

In this work we compare three different initial conditions for the classical Yang-Mills equations: the classical MV model (5) as well as fixed and running coupling JIMWLK evolution. We define the saturation scale Q_{\mathrm{s}}(Y) at rapidity Y through the expectation value of the dipole operator as

\frac{1}{{N_{\mathrm{c}}}}\left\langle\,\mathrm{Tr}\,V^{\dagger}({\mathbf{x}_{% T}})V({\mathbf{y}_{T}})\right\rangle_{Y,|{\mathbf{x}_{T}}-{\mathbf{y}_{T}}|=% \sqrt{2}/Q_{\mathrm{s}}}=e^{-1/2}. | (3) |

Throughout this paper we shall use Q_{\mathrm{s}} defined in this way from the light-like Wilson lines V({\mathbf{x}_{T}}) in the fundamental representation. The saturation scale is the only scale in the problem and we attempt to construct the various initial conditions in such a way that the value of Q_{\mathrm{s}}a is similar, to ensure a similar dependence on discretization effects.

In the MV model the Wilson lines are obtained from a classical color charge density \rho as

V({\mathbf{x}_{T}})=\mathbb{P}\exp\left\{i\int\,\mathrm{d}x^{-}g^{2}\frac{1}{% \boldsymbol{\nabla}_{T}^{2}}\rho^{a}({\mathbf{x}_{T}},x^{-})\right\}, | (4) |

where \mathbb{P} denotes path-ordering in x^{-}. The color charge density is a random variable with a local Gaussian probability distribution

P[\rho^{a}]\sim\exp\left\{-\int\,\mathrm{d}^{2}{\mathbf{x}_{T}}\,\mathrm{d}x^{% -}\frac{\rho^{a}({\mathbf{x}_{T}},x^{-})\rho^{a}({\mathbf{x}_{T}},x^{-})}{2\mu% ^{2}(x^{-})}\right\}, | (5) |

The total color charge \int\,\mathrm{d}x^{-}\mu^{2}(x^{-})\sim Q_{\mathrm{s}}^{2} is proportional to the thickness of a given nucleus.

In the numerical calculation the MV model initial conditions have been constructed as described in Ref. (18), discretizing the longitudinal coordinate Y in N_{y}=100 steps. For the calculations using the MV model directly for the initial conditions (1), (2) we have performed simulations on lattices of two different sizes: N_{\perp}=1024, with the MV model color charge parameter g^{2}\mu L=156 which translates into Q_{\mathrm{s}}a=0.119; and with N_{\perp}=2048, using g^{2}\mu L=550, which results in Q_{\mathrm{s}}a=0.172.

The MV model also provides the configurations used as the initial condition for quantum evolution in rapidity via the JIMWLK renormalization group equation, starting at Y=\log x_{0}/x=0. Performing a step \Delta Y in rapidity opens phase space for radiation of additional gluons which modify the classical action (4), (5). This process can be expressed as a “random walk” in the space of light-like Wilson lines V({\mathbf{x}_{T}}) (30); *Mueller:2001uk; (32); (33):

\displaystyle\partial_{Y}V({\mathbf{x}_{T}})=V({\mathbf{x}_{T}})\frac{i}{\pi}% \int\,\mathrm{d}^{2}{\mathbf{u}_{T}}\frac{({\mathbf{x}_{T}}-{\mathbf{u}_{T}})^% {i}\eta^{i}({\mathbf{u}_{T}})}{({\mathbf{x}_{T}}-{\mathbf{u}_{T}})^{2}} | (6) | ||

\displaystyle-\frac{i}{\pi}\int\,\mathrm{d}^{2}{\mathbf{v}_{T}}V({\mathbf{v}_{% T}})\frac{({\mathbf{x}_{T}}-{\mathbf{v}_{T}})^{i}\eta^{i}({\mathbf{v}_{T}})}{(% {\mathbf{x}_{T}}-{\mathbf{v}_{T}})^{2}}V^{\dagger}({\mathbf{v}_{T}})V({\mathbf% {x}_{T}}), |

where the Gaussian white noise \eta^{i}=\eta^{i}_{a}t^{a} satisfies \langle\eta^{a}_{i}({\mathbf{x}_{T}})\rangle=0 and, for fixed coupling,

\langle\eta^{a}_{i}({\mathbf{x}_{T}})\;\eta^{b}_{j}({\mathbf{y}_{T}})\rangle=% \alpha_{\mathrm{s}}\delta^{ab}\delta_{ij}\delta^{(2)}({\mathbf{x}_{T}}-{% \mathbf{y}_{T}}). | (7) |

Here the equation is written in the left-right symmetric form introduced in (34); (33).

The fixed coupling JIMWLK equation is solved using the numerical
method developed in (32); (35); (33).
For the smaller lattice size N_{\perp}=1024 we start with the MV model
with g^{2}\mu L=31 and without a mass regulator, which corresponds
to an initial Q_{\mathrm{s}}a=0.0218. After \Delta y=1.68/\alpha_{\mathrm{s}} units of
evolution in rapidity^{2}^{2}For fixed coupling the evolution
variable is \alpha_{\mathrm{s}}y, so we do not need to specify a particular value
of \alpha_{\mathrm{s}} separately. this results in Q_{\mathrm{s}}a=0.145. For a
N_{\perp}=2048-lattice we again start with g^{2}\mu L=31,
corresponding to Q_{\mathrm{s}}a=0.0107, and after \Delta y=1.8/\alpha_{\mathrm{s}}
units of evolution end up with Q_{\mathrm{s}}a=0.141.

For running coupling the evolution is significantly slower. We use the running coupling prescription introduced in (33), where the scale of the coupling is taken as the momentum conjugate to the distance in the noise correlator in Eq. (7). For the smaller N_{\perp}=1024 lattice we again start with g^{2}\mu L=31, i.e. Q_{\mathrm{s}}a=0.0218 and evolve for \Delta Y=10 units in rapidity, arriving at Q_{\mathrm{s}}a=0.118. For the larger N_{\perp}=2048 lattice we test a configuration that is farther from the IR cutoff, starting the JIMWLK evolution with g^{2}\mu L=102.4, i.e. Q_{\mathrm{s}}a=0.0423 and evolve for \Delta Y=10 units in rapidity, arriving at Q_{\mathrm{s}}a=0.172. In the rc-JIMWLK simulations the QCD scale is taken as \Lambda_{\mathrm{QCD}}a=0.00293 and the coupling is frozen to a value \alpha_{0}=0.76 in the infrared below 2.5\Lambda_{\mathrm{QCD}}.

As already mentioned above, RG evolution in rapidity resums quantum corrections to the fields \alpha_{m}^{\mu} of the individual charge sheets to all orders in \alpha_{s}\log 1/x, with leading logarithmic accuracy. In other words, the effective action at Y is modified from that at Y=0, written in Eq. (5).

Once an ensemble of Wilson lines V({\mathbf{x}_{T}}) at a rapidity Y is constructed, separately for both projectile and target, these configurations define \alpha_{1}^{i} and \alpha_{2}^{i} in light-cone gauge as written in Eq. (1); the initial field A^{i} of produced soft gluons at proper time \tau=+0 corresponds to their sum, Eq. (2). The evolution to \tau>0 follows from the real-time Hamiltonian evolution described in Ref. (14). This has been used in many classical field calculations, e.g. in Refs. (15); *Krasnitz:2003jw; *Lappi:2003bi, or more recently for the first study of the effects of JIMWLK-evolution on the gluon spectrum (19), and in the IP-glasma model for the initial conditions for hydrodynamics (36). On the N_{\perp}=2048 lattices we evolve the fields up to Q_{\mathrm{s}}\tau=5 and on the smaller N_{\perp}=1024 ones to Q_{\mathrm{s}}\tau=10. In this study, the nuclei are taken to fill the whole transverse lattice, with periodic boundary conditions.

## III Wilson loop

In the continuum the spatial (magnetic) Wilson loop is defined as the trace of a path ordered exponential of the gauge field around a closed path of area A in the transverse plane:

W(A)=\frac{1}{{N_{\mathrm{c}}}}\left<\,\mathrm{Tr}\,\mathbb{P}\exp\left\{ig% \oint_{\partial A}\,\mathrm{d}{\mathbf{x}_{T}}\cdot\mathbf{A}_{T}\right\}% \right>. | (8) |

On the lattice this is easily discretized as the product of link matrices around a square of area A. For N_{c}\geq 3 colors any particular Wilson loop is complex but the ensemble average is real.

We have measured the expectation value of the Wilson loop in the glasma field, with different initial conditions and at different times Q_{\mathrm{s}}\tau. The results of the calculation are shown in Fig. 1. As expected, the magnetic flux through a loop generically increases with its area. Focusing first on the curves corresponding to the initial time \tau=0 we observe that the resummation of quantum fluctuations (JIMWLK evolution) increases the flux through small loops of area AQ_{\mathrm{s}}^{2}<1. This can be understood intuitively as due to emission of additional virtual soft gluons in the pure gauge fields of the colliding charge sheets. On the other hand, the flux through large loops, AQ_{\mathrm{s}}^{2}\,\,\vbox{\hbox{$\buildrel\displaystyle>\over{\sim}$}}\,\,2, decreases. This indicates uncorrelated fluctuations of magnetic flux over such areas and is consistent with the suggestion that the flux is “bundled” in domains with a typical area \sim 1/Q_{\mathrm{s}}^{2} (20). Accordingly, loops of area \sim 1.5Q_{\mathrm{s}}^{2} are invariant under high-energy evolution.

Moving on to finite times we see that the flux through loops of fixed area decreases with \tau. This is, of course, a consequence of the decreasing field strength in an expanding metric. For small loops the ordering corresponding to the different initial conditions (MV, rc-JIMWLK, fc-JIMWLK) is preserved even at later times. However, for large loops one observes a striking “universality” emerging at Q_{\mathrm{s}}\tau\sim 5 as the curves for all initial conditions fall on top of each other.

The data from Fig. 1 shows an approximately linear dependence of \ln(-\ln W) on \ln(AQ_{\mathrm{s}}^{2}), with different slopes in the regime of small AQ_{\mathrm{s}}^{2}\ll 1 vs. large AQ_{\mathrm{s}}^{2}\gg 1. Based on this observation we fit the data to

W(A)=\exp\left\{-(\sigma A)^{\gamma}\right\}, | (9) |

with separate parametrizations for the IR and UV regimes:

\displaystyle\mathrm{IR:} | \displaystyle e^{0.5}<AQ_{\mathrm{s}}^{2}<e^{5} | (10) | |||

\displaystyle\mathrm{UV:} | \displaystyle e^{-3.5}<AQ_{\mathrm{s}}^{2}<e^{-0.5}~{}. | (11) |

In addition to limiting the fits to the quoted ranges we also restrict
them to the region where W>0.01 and the statistical error
on W is less than 0.2W; beyond these limits the data exhibits too large
fluctuations for a meaningful fit. Figures 2 and
3 show the time dependence of the exponents
\gamma in the IR and UV regions. The “string tension” \sigma
naturally decreases as \sim 1/\tau because of the longitudinal
expansion of the glasma, which leads to B_{z}\sim 1/\sqrt{\tau}. We
therefore show, in Figs. 4
and 5, the time dependence of the combination
\tau\sigma/Q_{\mathrm{s}}, where this leading effect is scaled out. The values
of \sigma/Q_{\mathrm{s}}^{2} for Q_{\mathrm{s}}\tau=0 are given in the
captions^{3}^{3}For the MV model, at \tau=0 we find
\sigma/Q_{\mathrm{s}}^{2}=0.44 in the IR region which is about four times
larger than the value reported in ref. (20). Our
present results refer to N_{c}=3 colors while
ref. (20) considered N_{c}=2; also, our current
definition of Q_{\mathrm{s}} via Eq. (3) leads to
smaller values for this quantity than the definition used
in (20). Finally, \sigma is extracted from fits
over a somewhat different range of areas..

The results in the “UV”-regime probed by Wilson loops of small area are shown in Figs. 2 and 4. They are easily understood from the differences in the initial condition. The gluon spectrum in the MV model falls steeply as a function of k_{T}, leading to a steep dependence of short-distance correlators on the distance. Our result for the UV exponent in Fig. 2 is close to the A^{2}-scaling obtained analytically in a weak field expansion (37). The difference is probably due to a combination of logarithmic corrections and lattice UV cutoff effects. For the JIMWLK ensembles the gluon spectrum is much harder (19), especially for fixed coupling. This manifests itself in smaller values of both \gamma and \sigma. In addition, the UV exponents are remarkably time independent at Q_{\mathrm{s}}\tau>1: this is consistent with the expectation that at such time the UV modes can be viewed as noninteracting gluons whose spectrum is close to the expectation from a perturbative k_{T}-factorized calculation (38).

The behavior in the IR regime (Figs. 3, 5) probed by large Wilson loops points to a very different picture. At \tau=0 the exponents \gamma and, to a lesser extent, the values of \sigma depend very much on the initial conditions. As already alluded to above, the scaling exponents \gamma_{\rm IR}<1 obtained for the JIMWLK fields indicate that quantum emissions increase magnetic flux fluctuations at the scale \sim 1/Q_{\mathrm{s}}, much smaller than the area of the loop. It is interesting to note that for the rather strong fixed-coupling evolution the initial scaling exponent is not too far above \gamma_{\rm IR}=1/2 corresponding to perimeter scaling.

At times Q_{\mathrm{s}}\tau\gtrsim 3, however, one observes a remarkable universality in the IR as the curves corresponding to different initial conditions collapse onto a single curve in Fig. 1. The string tensions in Fig. 5 are within 10% of each other at late Q_{\mathrm{s}}\tau, and the exponents \gamma in Fig. 3 are very close to each other, with values around \gamma_{\rm IR}\approx 1.2\dots 1.3. The exponent gradually decreases with \tau, potentially approaching the area law \gamma=1 at late times. The initial evolution points at a rapid rearrangement of “magnetic hot spots” to some universal field configurations at later time, Q_{\mathrm{s}}\tau\gtrsim 3.

We stress that the *universal* behavior of
large magnetic loops, characterized by a nontrivial power-law dependence
on the loop area, sets in at rather early time scales of a few times
1/Q_{\mathrm{s}}, independent of initial conditions.
Actual area law scaling \gamma=1 is approached only later.
This behavior
mirrors a similar universality between MV and JIMWLK results seen in
the IR part of the gluon spectrum (determined from correlators of
gauge fixed fields) in Ref. (19).
Since the structure of the fields does not seem to depend on the initial
conditions, we infer that this universality in due to stong interactions
in the glasma phase.
This universal behavior of
the Wilson loop for different initial conditions at Q_{\mathrm{s}}\tau\gtrsim 3 and AQ_{\mathrm{s}}^{2}\gg 1 is the main result of this paper.

## IV Magnetic field correlator

In this section we analyze gauge-invariant
two-point magnetic field correlators of the
form^{4}^{4}We include the factor g^{2} for convenience, because the
quantity that appears naturally in the classical lattice formulation
is actually gB.

C_{B}(r)\equiv 2g^{2}\,\mathrm{Tr}\,\left<B_{z}({\mathbf{x}_{T}})U_{{\mathbf{x% }_{T}}\to{\mathbf{y}_{T}}}B_{z}({\mathbf{y}_{T}})U^{\dagger}_{{\mathbf{x}_{T}}% \to{\mathbf{y}_{T}}}\right>~{}. | (12) |

The points {\mathbf{x}_{T}} and {\mathbf{y}_{T}} are separated in the x or y direction by a distance r=|{\mathbf{x}_{T}}-{\mathbf{y}_{T}}|, and the Wilson line U_{{\mathbf{x}_{T}}\to{\mathbf{y}_{T}}} is the ordered product of links along the straight line separating these points.

The magnetic field B_{z}=t^{a}B_{z}^{a} on the lattice is defined as the traceless antihermitian part of the plaquette as

gB_{z}^{a}({\mathbf{x}_{T}})=2\,\mathrm{Re}\,\mathrm{Tr}\,t^{a}U_{x,y}({% \mathbf{x}_{T}})~{}, | (13) |

where the transverse plaquette is

U_{i,j}({\mathbf{x}_{T}})=U_{i}({\mathbf{x}_{T}})U_{j}({\mathbf{x}_{T}}+% \mathbf{i}_{T})U^{\dagger}_{i}({\mathbf{x}_{T}}+\mathbf{j}_{T})U^{\dagger}_{j}% ({\mathbf{x}_{T}})~{}. | (14) |

Here U_{i}({\mathbf{x}_{T}}) denotes the link matrix in the i-direction based at {\mathbf{x}_{T}} and \mathbf{i}_{T},\mathbf{j}_{T} are unit vectors.

The resulting magnetic field correlator r\;C_{B}(r) is plotted in
Fig. 6. We have multiplied by r to better expose the
behavior around r\sim 1/Q_{\mathrm{s}}. At the initial time there is a
significant anticorrelation at intermediate distances. It shows the
domain structure of the field such that B_{z} is likely to flip
sign^{5}^{5}Recall that B transforms homogeneously. Hence, unlike
the links U_{i}({\mathbf{x}_{T}}), the magnetic field B_{z}({\mathbf{x}_{T}}) can be
diagonalized everywhere by a suitable gauge transformation. (or
direction) over distances of order 1/Q_{\mathrm{s}}. This structure then
changes very rapidly: already at time Q_{\mathrm{s}}\tau\sim 2 the fields
have rearranged such that the anticorrelation has disappeared. Also,
the subsequent time evolution results in damping of the fluctuations
at Q_{\mathrm{s}}r\,\,\vbox{\hbox{$\buildrel\displaystyle>\over{\sim}$}}\,\,1 which are present in the initial field
configurations. On the other hand, the strong short-distance
correlations around the peak are not affected much by the time
evolution beyond Q_{\mathrm{s}}\tau\sim 2, aside from a decrease in
magnitude. In particular, no “infrared diffusion” of the peak
towards larger distances is observed.

Given the clear scaling behavior of the Wilson loop one might expect to see a similar phenomenon for the magnetic field correlator. A very naive scaling argument would assume that if C(r)\sim r^{-\alpha} then the Wilson loop should scale as

\displaystyle-\ln W\sim\int\limits_{A}\,\mathrm{d}^{2}{\mathbf{x}_{T}}\,% \mathrm{d}^{2}{\mathbf{y}_{T}}\,C(|{\mathbf{x}_{T}}-{\mathbf{y}_{T}}|)\\ \displaystyle\sim R^{4-\alpha}\sim A^{\frac{4-\alpha}{2}}=A^{\gamma}~{}. | (15) |

The area integrals in the first line extend over |{\mathbf{x}_{T}}|,|{\mathbf{y}_{T}}|<R. Thus, \gamma=1.225 extracted from the Wilson loop at Q_{\mathrm{s}}\tau=10 would give C_{B}(r)\sim r^{-1.55}. On a logarithmic scale C_{B}(r) does indeed qualitatively resemble such behavior as shown in Fig. 7. However, this kind of scaling is less conclusive than for the Wilson loop (see also appendix); this could be an indication for the presence of higher cumulants in the expansion of the spatial Wilson loop (39); *DiGiacomo:2000va.

## V Summary

In this paper we have provided some insight into the fields produced initially in a high-energy collision of dense color charge sheets. We have focused, in particular, on the structure of the longitudinal magnetic field B_{z}.

We consider both purely classical as well as JIMWLK RG evolved gauge field ensembles on which we measure expectation values of spatial Wilson loops and two-point correlation functions of B_{z}. These show that the initial fields exhibit domain-like structure over distance scales of the order of the inverse saturation scale 1/Q_{\mathrm{s}}. Classical YM evolution to later times leads to universal scaling, for all ensembles, of the magnetic loop with area, with a nontrivial critical exponent. Also, the anti-correlation of B_{z}({\mathbf{x}_{T}}) over distances \sim 1/Q_{\mathrm{s}} disappears, which we interpret as rearrangement, possibly accompanied by transverse expansion, of the magnetic field domains.

The emergence of a color field condensate in high-energy collisions of dense hadrons or nuclei is a very interesting phenomenon, and its dynamics remains to be understood in more detail. In closing we only draw attention to recent arguments that the presence of such a condensate might have important implications for the process of (pre-) thermalization in high multiplicity collisions (41).

## Appendix A Relation between the Wilson loop and magnetic field correlator

In an abelian theory there is a simple relation between the Wilson loop and the magnetic field due to Stokes’ theorem:

\oint_{\partial A}\,\mathrm{d}{\mathbf{x}}\cdot\mathbf{A}=\int_{A}\,\mathrm{d}% ^{2}{\mathbf{x}}\,B_{z}({\mathbf{x}})~{}. | (16) |

If we assume that in the nonabelian case the magnetic field in each color channel a is independent, and that it consists of uncorrelated domains which are much smaller than the area A and distributed as Gaussian random variables, we obtain the following estimate for the Wilson loop:

\displaystyle\frac{1}{{N_{\mathrm{c}}}}\,\mathrm{Tr}\,\exp\left\{ig\oint_{% \partial A}\,\mathrm{d}{\mathbf{x}}\cdot\mathbf{A}\right\}\\ \displaystyle\approx\exp\left\{-\frac{g^{2}}{2{N_{\mathrm{c}}}}\left<\,\mathrm% {Tr}\,\left[\int_{A}\,\mathrm{d}^{2}{\mathbf{x}}\,B_{z}({\mathbf{x}})\right]^{% 2}\right>\right\}\\ \displaystyle=\exp\left\{-\frac{1}{4{N_{\mathrm{c}}}}\int_{A}\,\mathrm{d}^{2}{% \mathbf{x}}\,\mathrm{d}^{2}{\mathbf{y}}\,C_{B}({\mathbf{x}}-{\mathbf{y}})% \right\}. | (17) |

In Fig. 8 we compare the result of a numerical integration of the r.h.s. of Eq. (17) using the measured magnetic field correlator, to the direct measurement of the Wilson loop. It can be seen that the two are in a relatively good agreement. This consistency check supports the interpretation of B_{z} as independent field domains of area \sim 1/Q_{\mathrm{s}}^{2}.

## Acknowledgements

T. L. is supported by the Academy of Finland, projects 133005, 267321 and 273464. This work was done using computing resources from CSC – IT Center for Science in Espoo, Finland. A. D. acknowledges support by the DOE Office of Nuclear Physics through Grant No. DE-FG02-09ER41620 and from The City University of New York through the PSC-CUNY Research Award Program, grant 66514-0044. The authors thank the Yukawa Institute for Theoretical Physics, Kyoto University, where part of this work was done during the YITP-T-13-05 workshop on “New Frontiers in QCD”.

10