Parton model description of multiparticle azimuthal correlations in pA collisions

In K. Dusling et al. [arXiv: 1705.00745.], an initial-state parton model of quarks scattering off a dense nuclear target was shown to qualitatively reproduce the systematics of multiparticle azimuthal anisotropy cumulants measured in proton/deuteron-nucleus (pA) collisions at the Relativistic Heavy Ion Collider and the LHC. The systematics included (i) the behavior of the four-particle cumulant c(2){4}, which generates a real four-particle second Fourier harmonic v(2){4}; (ii) the ordering v(2){2} > v(2){4} approximate to v(2){6} approximate to v(2){8} for two-, four-, six-, and eight-particle Fourier harmonics; and (iii) the behavior of so-called symmetric cumulants SC(2,3) and SC(2,4). These features of azimuthal multiparticle cumulants were previously interpreted as a signature of hydrodynamic flow; our results challenge this interpretation. We expand here upon our previous study and present further details and novel results on the saturation scale and transverse momentum (p(perpendicular to)) dependence of multiparticle azimuthal correlations. We find that the dependence of v(2){2} and v(2){4} on the number of color domains in the target varies with the p(perpendicular to) window explored. We extend our prior discussion of symmetric cumulants and compute as yet unmeasured symmetric cumulants. We investigate the N-c dependence of v(2){2} and v(2){4}. We contrast our results, which include multiple scatterings of each quark off the target, to the Glasma graph approximation, where each quark suffers at most two-gluon exchanges with the target. We find that coherent multiple scattering is essential to obtain a positive definite v(2){4}. We provide an algorithm to compute expectation values of arbitrary products of the dipole lightlike Wilson line correlators.