We carry out hydrodynamical simulation of the evolution of fluid in relativistic heavy-ion collisions with random initial fluctuations. The time evolution of power spectrum of momentum anisotropies shows very strong correspondence with the physics of cosmic microwave anisotropies as was earlier predicted by some of us. In particular our results demonstrate suppression of superhorizon fluctuations and the correspondence between the location of the first peak in the power spectrum of momentum anisotropies and the length scale of fluctuations and expected freezeout time scale (more precisely, the sound horizon size at freezeout).

We propose a simple technique to detect any anisotropic expansion stage in the history of the universe starting from the inflationary stage to the surface of last scattering from the CMBR data. We use the property that any anisotropic expansion in the universe would deform the shapes of the primordial density perturbations and this deformation can be detected in a shape analysis of superhorizon fluctuations in CMBR. Using this analysis we obtain the constraint on any previous anisotropic expansion of the universe to be less than about 35%.

It is known that the presence of background magnetic field in cosmic plasma distorts the acoustic peaks in CMBR. This primarily results from different types of waves in the plasma with velocities depending on the angle between the magnetic field and the wave vector. We consider the consequences of these effects in relativistic heavy-ion collisions where very strong magnetic fields arise during early stages of the plasma evolution. We show that flow coefficients can be significantly affected by these effects when the magnetic field remains strong during early stages due to strong induced fields in the conducting plasma. In particular, the presence of magnetic field can lead to enhancement in the elliptic flow coefficient $v_2$.

Recently we have shown that there are crucial similarities in the physics of cosmic microwave background radiation (CMBR) anisotropies and the flow anisotropies in relativistic heavy-ion collision experiments (RHICE). We also argued that, following CMBR anisotropy analysis, a plot of root-mean square values of the flow coefficients, calculated in a lab fixed frame for RHICE, can yield important information about the nature of initial state anisotropies and their evolution. Here we demonstrate the strength of this technique by showing that elliptic flow for non-central collisions can be directly determined from such a plot without any need for the determination of event-plane.

We focus on the initial state spatial anisotropies, originating at the thermalization stage, for central collisions in relativistic heavy-ion collisions. We propose that a plot of the root mean square values of the flow coefficients $\sqrt{\bar {v_n^2}} \equiv v_n^{rms}$, calculated in a lab fixed coordinate system, for a large range of $n$, from 1 to about 30, can give non-trivial information about the initial stages of the system and its evolution. We also argue that for all wavelengths $\lambda$ of the anisotropy (at the surface of the plasma region) much larger than the acoustic horizon size $H_s^{fr}$ at the freezeout stage, the resulting values of $v_n^{rms}$ should be suppressed by a factor of order $2H_s^{fr}/\lambda$. With initial flow being zero, we discuss the possibility that the resulting flow could show imprints of coherent oscillations in the plot of $v_n^{rms}$ for sub-horizon modes. For gold-gold collision at 200 GeV, these features are expected to occur for n $\ge$ 5, with $n < 4$ modes showing supression. This has strong similarities with the physics of the anisotropies of the cosmic microwave background radiation (CMBR) resulting from inflationary density fluctuations in the universe. It seems possible that the statistical fluctuations due to finite multiplicity may not be able to mask such features in the flow data, or, at least a non-trivial overall shape of the plot of $v_n^{rms}$ may be inferred. In that case, the successes of analysis of CMBR anisotropy power spectrum to get cosmological parameters can be applied for relativistic heavy-ion collisions to learn about various relevant parameters at the early stages of the evolving system.