Search

We investigate the effects of boost invariance breaking on the isotropization of pressure in the glasma, using the $3+1$D glasma simulation. The breaking is attributed to spatial fluctuations of the classical color charge density along the collision axis. We present numerical results for pressure and energy density at mid-rapidity and across a wider rapidity region. It is found that, despite varying longitudinal correlation lengths, the behaviors of the pressure isotropizations are qualitatively similar. The numerical results suggest that, in the initial stage, longitudinal color electromagnetic fields develop, similar to those in the boost invariant glasma. Subsequently, these fields evolve into a dilute glasma, expanding longitudinally in a manner akin to a dilute gas. We also show that the energy density at mid-rapidity exhibits a $1/\tau$ decay in the dilute glasma stage., Comment: 16 page, 6 figures

We propose a new numerical method for $3+1$D glasma simulation using Milne coordinates. We formulate the classical Yang-Mills field and $3$D classical color current on a lattice at the initial proper time, specified as a moment just before the collision of the two nuclei. By solving the evolution equations, we extract observables of the $3$D glasma at later times. We demonstrate the efficiency of our method in terms of numerical cost and apply it to the central collisions of Au-Au. We also discuss possible further improvements of our method., Comment: 14 pages, 7 figures

We investigate efficiency of a gauge-covariant neural network and an approximation of the Jacobian in optimizing the complexified integration path toward evading the sign problem in lattice field theories. For the construction of the complexified integration path, we employ the path optimization method. The $2$-dimensional $\text{U}(1)$ gauge theory with the complex gauge coupling constant is used as a laboratory to evaluate the efficiency. It is found that the gauge-covariant neural network, which is composed of the Stout-like smearing, can enhance the average phase factor, as the gauge-invariant input does. For the approximation of the Jacobian, we test the most drastic case in which we perfectly drop the Jacobian during the learning process. It reduces the numerical cost of the Jacobian calculation from ${\cal O}(N^3)$ to ${\cal O}(1)$, where $N$ means the number of degrees of freedom of the theory. The path optimization using this Jacobian approximation still enhances the average phase factor at expense of a slight increase of the statistical error., Comment: 8 pages, 5 figures; accepted version

We investigate the possible thermalization process of the highly occupied and weakly coupled Yang-Mills fields expanding along the beam axis through an evaluation of the entropy, particle number, and pressure anisotropy. The time evolution of the system is calculated by solving the equation of motion for the Wigner function in the semiclassical approximation with initial conditions mimicking the glasma. For the evaluation of the entropy, we adopt the Husimi-Wehrl (HW) entropy, which is obtained by using the Husimi function, a positive semidefinite quantum distribution function given by smearing the Wigner function. By numerical calculations at $g=0.1$ and $0.2$, the entropy production is found to occur together with the particle creation in two distinct stages: In the first stage, the particle number and the entropy at low longitudinal momenta grow rapidly. In the second stage, the particle number and the entropy of higher longitudinal momentum modes show slower increase. The pressure anisotropy remains in our simulation and implies that the system is still out-of-equilibrium., Comment: 28 pages, 15 figures

Real-time evolution of replicas of classical field is proposed as an approximate simulator of real-time quantum field dynamics at finite temperatures. We consider $N$ classical field configurations dubbed as replicas which interact with each other via the $\tau$-derivative terms and evolve with the classical equation of motion. The partition function of replicas is found to be proportional to that of quantum field in the imaginary time formalism. As the replica index $\tau$ can be regarded as the imaginary time index, the replica evolution is technically the same as the molecular dynamics part of the hybrid Monte-Carlo sampling and the replica configurations should reproduce the correct quantum equilibrium distribution after the long-time evolution. At the same time, evolution of the replica-index average of field variables is described by the classical equation of motion when the fluctuations are small. In order to examine the real-time propagation properties of replicas, we first discuss replica evolution in quantum mechanics. Statistical averages of observables are precisely obtained by the initial condition average of replica evolution, and the time evolution of the unequal-time correlation function, $\langle x(t) x(t')\rangle$, in a harmonic oscillator is also described well by the replica evolution in the range $T/\omega > 0.5$. Next, we examine the statistical and dynamical properties of the $\phi^4$ theory in the 4+1 dimensional spacetime, which contains three spatial, one replica index or the imaginary time, and one real-time. We note that the Rayleigh-Jeans divergence can be removed in replica evolution with $N \geq 2$ when the mass counterterm is taken into account. We also find that the thermal mass obtained from the unequal-time correlation function at zero momentum grows as a function of the coupling as in the perturbative estimate in the small coupling region., Comment: 27 pages, 11 figures (text and Figs. 6 and 7 are updated in v2)

We investigate the shear viscosity $\eta$ of the classical Yang-Mills (CYM) field on a lattice by using the Green-Kubo formula, where the shear viscosity is calculated from the time-correlation function of the energy-momentum tensor in equilibrium. Dependence of the shear viscosity $\eta(g,T)$ on the coupling $g$ and temperature $T$ is represented by a scaling function $f_\eta(g^2T)$ as $\eta(g,T)=Tf_\eta(g^2T)$ due to the scaling-invariant property of the CYM. The explicit functional form of $f_\eta(g^2T)$ is successfully determined from the calculated shear viscosity: It turns out that $\eta(g,T)$ of the CYM field is proportional to $1/g^{1.10-1.88}$ at weak coupling, which is a weaker dependence on $g$ than that in the leading-order perturbation theory but consistent with that of the "anomalous viscosity" $\eta\propto 1/g^{1.5}$ under the strong disordered field. The obtained shear viscosity is also found to be roughly consistent with that estimated through the analysis of the anisotropy of the pressure of the CYM dynamics in the expanding geometry with recourse to a hydrodynamic equation., Comment: 10pages, 11figures

We investigate the shear viscosity of massless classical scalar fields in the $\phi^4$ theory on a lattice by using the Green-Kubo formula. Based on the scaling property of the classical field, the shear viscosity is represented using a scaling function. Equilibrium expectation value of the time-correlation function of the energy-momentum tensor is evaluated as the ensemble average of the classical field configurations, whose time evolution is obtained by solving the classical equation of motion starting from the initial condition in thermal equilibrium. It is found that there are two distinct damping time scales in the time-correlation function, which is found to show damped oscillation behavior in the early stage around a slow monotonous decay with an exponential form, and the slow decay part is found to dominate the shear viscosity in the massless classical field theory. This kind of slow decay is also known to exist in the molecular dynamics simulation, then it may be a generic feature of dense matter., Comment: 21 pages, 11 figures