Using strong-disorder renormalization group, numerical exact diagonalization, and quantum Monte Carlo methods, we revisit the random antiferromagnetic XXZ spin-1/2 chain focusing on the long-length and ground-state behavior of the average time-indepe
ndent spin-spin correlation function C(l)=upsilon l^{-eta}. In addition to the well-known universal (disorder-independent) power-law exponent eta=2, we find interesting universal features displayed by the prefactor upsilon=upsilon_o/3, if l is odd, and upsilon=upsilon_e/3, otherwise. Although upsilon_o and upsilon_e are nonuniversal (disorder dependent) and distinct in magnitude, the combination upsilon_o + upsilon_e = -1/4 is universal if C is computed along the symmetric (longitudinal) axis. The origin of the nonuniversalities of the prefactors is discussed in the renormalization-group framework where a solvable toy model is considered. Moreover, we relate the average correlation function with the average entanglement entropy, whose amplitude has been recently shown to be universal. The nonuniversalities of the prefactors are shown to contribute only to surface terms of the entropy. Finally, we discuss the experimental relevance of our results by computing the structure factor whose scaling properties, interestingly, depend on the correlation prefactors.
