Since the launch of LHC experiments it has been discovered that the high multiplicity trigger in pp, pA collisions finds events behaving differently from the typical (minimally biased) ones. In central pPb case it has been proven that those possess collective phenomena known as the radial, elliptic and triangular flows, similar to what is known in heavy ion (AA) collisions. In this paper we argue that at the ultra-high energies, E_lab ~ 10^{20} eV, of the observed cosmic rays this regime changes from a small-probability fluctuation to a dominant one. We estimate velocity of the transverse collective expansion for the light-light and heavy-light collisions, and find it comparable to what is observed at LHC for the central PbPb case. We argue that significant changes of spectra of various secondaries associated with this phenomenon should be important for the development of the cosmic ray cascades.

Inhomogeneities associated with the cosmological QCD and electroweak phase transitions produce hydrodynamical perturbations, longitudinal sounds and rotations. It has been demonstrated by Hindmarsh et al. that the sounds produce gravity waves (GW) well after the phase transition is over. We further argue that, under certain conditions, an inverse acoustic cascade may occur and move sound perturbations from the (UV) momentum scale at which the sound is originally produced to much smaller (IR) momenta. The weak turbulence regime of this cascade is studied via the Boltzmann equation, possessing stationary power and time-dependent self-similar solutions. We suggest certain indices for the strong turbulence regime as well, into which the cascade eventually proceeds. Finally, we point out that two on-shell sound waves can produce one on-shell gravity wave, and we evaluate the rate of the process using a standard sound loop diagram.

Light scalar fields very naturally appear in modern cosmological models, affecting such parameters of Standard Model as electromagnetic fine structure constant $\alpha$, dimensionless ratios of electron or quark mass to the QCD scale, $m_{e,q}/\Lambda_{QCD}$. Cosmological variations of these scalar fields should occur because of drastic changes of matter composition in Universe: the latest such event is rather recent (redshift $z\sim 0.5$), from matter to dark energy domination. In a two-brane model (we use as a pedagogical example) these modifications are due to changing distance to "the second brane", a massive companion of "our brane". Back from extra dimensions, massive bodies (stars or galaxies) can also affect physical constants. They have large scalar charge $Q_d$ proportional to number of particles which produces a Coulomb-like scalar field $\phi=Q_d/r$. This leads to a variation of the fundamental constants proportional to the gravitational potential, e.g. $\delta \alpha/ \alpha = k_\alpha \delta (GM/ r c^2)$. We compare different manifestations of this effect. The strongest limits $k_\alpha +0.17 k_e= (-3.5\pm 6) * 10^{-7}$ are obtained from the measurements of dependence of atomic frequencies on the distance from Sun (the distance varies due to the ellipticity of the Earth's orbit).

We address the emerging discrepancy between the Big Bang Nucleosynthesis data and standard cosmology, which asks for a bit longer evolution time. If this effect is real, one possible implication (in a framework of brane cosmology model) is that there is a ``dark radiation'' component which is negative and makes few percents of ordinary matter density. If so, all scales of this model can be fixed, provided brane-to-bulk leakage problem is solved.

We follow our previous paper on possible cosmological variation of weak scale (quark masses) and strong scale, inspired by data on cosmological variation of the electromagnetic fine structure constant from distant quasar (QSO) absorption spectra. In this work we identify the {\em strange quark mass} $m_s$ as the most important quantity, and the {\em sigma meson mass} as the ingredient of the nuclear forces most sensitive to it. As a result, we claim significantly stronger limits on ratio of weak/strong scale ($W=m_s/\Lambda_{QCD}$) variation following from our previous discussion of primordial Big-Bang Nucleosynthesis ($|\delta W/W|<0.006$) and Oklo natural nuclear reactor ($|\delta W/W|<1.2 \cdot 10^{-10}$; there is also a non-zero solution $\delta W/W=(-0.56 \pm 0.05) \cdot 10^{-9}$) .

Recent data on cosmological variation of the electromagnetic fine structure constant from distant quasar (QSO) absorption spectra have inspired a more general discussion of possible variation of other constants. We discuss variation of strong scale and quark masses. We derive the limits on their relative change from (i) primordial Big-Bang Nucleosynthesis (BBN); (ii) Oklo natural nuclear reactor, (iii) quasar absorption spectra, and (iv) laboratory measurements of hyperfine intervals.

We review recent work on the phase structure of QCD at very high baryon density. We introduce the phenomenon of color superconductivity and discuss how the quark masses and chemical potentials determine the structure of the superfluid quark phase. We comment on the possibility of kaon condensation at very high baryon density and study the competition between superfluid, density wave, and chiral crystal phases at intermediate density.