Anisotropic fluctuations of momentum and angular momentum of heavy quarks in the pre-equilibrium stage of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>p</mml:mi> <mml:mi>A</mml:mi> </mml:mrow> </mml:ma

We simulate the real-time evolution of the SU(3)-glasma generated in the early stages of high-energy proton-nucleus collisions, employing classical lattice gauge theory techniques. Our setup incorporates a realistic modeling of the proton’s internal structure and includes longitudinal fluctuations in the initial state, enabling the study of genuinely nonboost-invariant collision dynamics. Focusing on the momentum and angular momentum anisotropies of heavy quarks in the infinite mass limit, we find that the system retains significant anisotropy well beyond the characteristic timescale <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"> <a:mi>τ</a:mi> <a:mo>=</a:mo> <a:mn>1</a:mn> <a:mo>/</a:mo> <a:msub> <a:mi>Q</a:mi> <a:mi>s</a:mi> </a:msub> </a:math> , with <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"> <c:msub> <c:mi>Q</c:mi> <c:mi>s</c:mi> </c:msub> </c:math> being the saturation energy of the glasma, the only scale in our model. This persistence of anisotropy is further confirmed in the more realistic, nonboost-invariant scenario, across a range of fluctuation amplitudes. These findings pave the way for future investigations involving dynamical heavy quarks and more quantitative initializations of the glasma.