We propose a cross-correlation method for the searches of ultra-light fields, in particular, with a space network of atomic sensors. The main motivation of the approach is cancellation of uncorrelated noises in the observation data and unique pattern the fields leave on the cross-spectrum, depending on their nature (i.e., scalar, vector or tensor). In particular, we analytically derive a dependence of the cross-spectrum on the angle between two pairs of detectors. We then confirm obtained angular curves with a numerical simulation. We imply application of the method to the detection of dark matter and gravitational waves.

We analyze possible effects of the dark matter environment on the atomic clock stability measurements. The dark matter is assumed to exist in a form of waves of ultralight scalar fields or in a form of topological defects (monopoles and strings). We identify dark matter signal signatures in clock Allan deviation plots that can be used to constrain the dark matter coupling to the Standard Model fields. The existing data on the Al+/Hg+ clock comparison are used to put new limits on the dilaton dark matter in the region of masses m > 10^{-15} eV. We also estimate the sensitivities of future atomic clock experiments in space, including the cesium microwave and strontium optical clocks aboard the International Space Station, as well as a potential nuclear clock. These experiments are expected to put new limits on the topological dark matter in the range of masses 10^{-10} eV < m < 10^{-6} eV.

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.

The (absence of detecting) lensing dispersion of Supernovae type Ia (SNIa) can be used as a novel and extremely efficient probe of cosmology. In this preliminary example we analyze its consequences for the primordial power spectrum. The main setback is the knowledge of the power spectrum in the non-linear regime, 1 Mpc^{-1} < k < 10^2-10^3 Mpc^{-1} up to redshift of about unity. By using the lensing dispersion and conservative estimates in this regime of wavenumbers, we show how the current upper bound sigma_mu(z=1) < 0.12on existing data gives strong indirect constraints on the primordial power spectrum. The probe extends our handle on the spectrum to a total of 12-15 inflation e-folds. These constraints are so strong that they are already ruling out a large portion of the parameter space allowed by PLANCK for running alpha = d n_s / d ln k and running of running beta = d^2 n_s / d ln k^2. The bounds follow a linear relation to a very good accuracy. A conservative bound disfavors any enhancement above the line beta(k_0)=0.036-0.42 alpha(k_0) and a realistic estimate disfavors any enhancement above the line beta(k_0)=0.022-0.44 alpha(k_0).