Simulation of a <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>3</mml:mn><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mi mathvariant="norma

We apply the 3D glasma simulation method using Milne coordinates, proposed in our previous work [H. Matsuda and X.-G. Huang, [] ], to the early stage of the Au-Au collisions at <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:msqrt><a:msub><a:mi>s</a:mi><a:mrow><a:mi>NN</a:mi></a:mrow></a:msub></a:msqrt><a:mo>=</a:mo><a:mn>200</a:mn><a:mtext> </a:mtext><a:mtext> </a:mtext><a:mi>GeV</a:mi></a:math>. The nucleus model prior to the collisions, which offers the initial condition for the 3D glasma simulation, is constructed to account for the longitudinal structure of the nucleus, the finite thickness of nucleons, and their random positions along the collision axis. We investigate rapidity profiles for a wide range of physical quantities of the glasma, including energy, pressure, fluctuations of topological charge, eccentricity, and angular momentum. In particular, we elucidate the behavior of eccentricity and angular momentum, which are physical quantities dependent on the geometric shape of the glasma, across a wide range of impact parameter regions. Published by the American Physical Society 2024