We review and compare two different CMB dipole estimators discussed in the literature, and assess their performances through Monte Carlo simulations. The first method amounts to simple template regression with partial sky data, while the second method is an optimal Wiener filter (or Gibbs sampling) implementation. The main difference between the two methods is that the latter approach takes into account correlations with higher-order CMB temperature fluctuations that arise from non-orthogonal spherical harmonics on an incomplete sky, which for recent CMB data sets (such as Planck) is the dominant source of uncertainty. For an accepted sky fraction of 81% and an angular CMB power spectrum corresponding to the best-fit Planck 2018 $\Lambda$CDM model, we find that the uncertainty on the recovered dipole amplitude is about six times smaller for the Wiener filter approach than for the template approach, corresponding to 0.5 and 3$~\mu$K, respectively. Similar relative differences are found for the corresponding directional parameters and other sky fractions. We note that the Wiener filter algorithm is generally applicable to any dipole estimation problem on an incomplete sky, as long as a statistical and computationally tractable model is available for the unmasked higher-order fluctuations. The methodology described in this paper forms the numerical basis for the most recent determination of the CMB solar dipole from Planck, as summarized by arXiv:2007.04997.

We constrain the spectral index of polarized synchrotron emission, $\beta_s$, by correlating the recently released 2.3 GHz S-Band Polarization All Sky Survey (S-PASS) data with the 23 GHz 9-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We sub-divide the S-PASS field, which covers the Southern Ecliptic hemisphere, into $15^{\circ}\times 15^{\circ}$ regions, and estimate the spectral index of polarized synchrotron emission within each region using a simple but robust T-T plot technique. Three different versions of the S-PASS data are considered, corresponding to either no correction for Faraday rotation; Faraday correction based on the rotation measure model presented by the S-PASS team; or Faraday correction based on a rotation measure model presented by Hutschenreuter and En{\ss}lin. We find that the correlation between S-PASS and WMAP is strongest when applying the S-PASS model. Adopting this correction model, we find that the mean spectral index of polarized synchrotron emission gradually steepens from $\beta_s\approx-2.8$ at low Galactic latitudes to $\beta_s\approx-3.2$ at high Galactic latitudes, in good agreement with previously published results. Finally, we consider two special cases defined by the BICEP2 and SPIDER fields, and obtain mean estimates of $\beta_{BICEP2}=-3.22\pm0.06$ and $\beta_{SPIDER}=-3.21\pm0.03$, respectively. Adopting the WMAP 23 GHz sky map bandpass filtered to including angular scales only between $2^{\circ}$ and $10^{\circ}$ as a spatial template, we constrain the root-mean-square synchrotron polarization amplitude to be less than $0.03\mu K$ ($0.009\mu K$) at 90 GHz (150 GHz) for the BICEP2 field, corresponding roughly to a tensor-to-scalar ratio of $r\lesssim0.02$ ($r\lesssim0.005$), respectively. Very similar constraints are obtained for the SPIDER field.

We study the late-time evolution of the Universe where dark energy (DE) is parametrized by a modified generalized Chaplygin gas (mGCG) on top of cold dark matter (CDM). We also take into account the radiation content of the Universe. In this context, the late stage of the evolution of the universe refers to the epoch where CDM is already clustered into inhomogeneously distributed discrete structures (galaxies, groups and clusters of galaxies). Under these conditions, the mechanical approach is an adequate tool to study the Universe deep inside the cell of uniformity. To be more accurate, we study scalar perturbations of the Friedmann-Lema\^itre-Robertson-Walker metric due to inhomogeneities of CDM as well as fluctuations of radiation and mGCG, the later driving the late-time acceleration of the universe. Our analysis applies as well to the case where mGCG plays the role of DM and DE. We select the sets of parameters of the mGCG that are compatible with the mechanical approach. These sets define prospective mGCG models. By comparing the selected sets of models with some of the latest observational data results, we conclude that the mGCG is in tight agreement with those observations particularly for a mGCG playing the role of DE and DM.

We point out a weak side of the commonly used determination of scalar cosmological perturbations lying in the fact that their average values can be nonzero for some matter distributions. It is shown that introduction of the finite-range gravitational potential instead of the infinite-range one resolves this problem. The concrete illustrative density profile is investigated in detail in this connection.

We study the late-time evolution of the Universe where dark energy (DE) is presented by a barotropic fluid on top of cold dark matter (CDM). We also take into account the radiation content of the Universe. Here by the late stage of the evolution we refer to the epoch where CDM is already clustered into inhomogeneously distributed discrete structures (galaxies, groups and clusters of galaxies). Under this condition the mechanical approach is an adequate tool to study the Universe deep inside the cell of uniformity. More precisely, we study scalar perturbations of the FLRW metric due to inhomogeneities of CDM as well as fluctuations of radiation and DE. For an arbitrary equation of state for DE we obtain a system of equations for the scalar perturbations within the mechanical approach. First, in the case of a constant DE equation of state parameter $w$, we demonstrate that our method singles out the cosmological constant as the only viable dark energy candidate. Then, we apply our approach to variable equation of state parameters in the form of three different linear parametrizations of $w$, e.g., the Chevallier-Polarski-Linder perfect fluid model. We conclude that all these models are incompatible with the theory of scalar perturbations in the late Universe.

We consider lattice Universes with spatial topologies $T\times T\times T$, $\; T\times T\times R\; $ and $\; T\times R\times R$. In the Newtonian limit of General Relativity, we solve the Poisson equation for the gravitational potential in the enumerated models. In the case of point-like massive sources in the $T\times T\times T$ model, we demonstrate that the gravitational potential has no definite values on the straight lines joining identical masses in neighboring cells, i.e. at points where masses are absent. Clearly, this is a nonphysical result since the dynamics of cosmic bodies is not determined in such a case. The only way to avoid this problem and get a regular solution at any point of the cell is the smearing of these masses over some region. Therefore, the smearing of gravitating bodies in $N$-body simulations is not only a technical method but also a physically substantiated procedure. In the cases of $\; T\times T\times R\; $ and $\; T\times R\times R$ topologies, there is no way to get any physically reasonable and nontrivial solution. The only solutions we can get here are the ones which reduce these topologies to the $T\times T\times T$ one.

In this paper we consider the Universe at the late stage of its evolution and deep inside the cell of uniformity. At these scales, the Universe is filled with inhomogeneously distributed discrete structures (galaxies, groups and clusters of galaxies). Supposing that a small fraction of colored objects escaped hadronization and survived up to now in the form of quark-gluon nuggets (QNs), and also taking into account radiation, we investigate scalar perturbations of the FRW metrics due to inhomogeneities of dustlike matter as well as fluctuations of QNs and radiation. In particular, we demonstrate that the nonrelativistic gravitational potential is defined by the distribution of inhomogeneities/fluctuations of both dustlike matter and QNs. Consequently, QNs can be distributed around the baryonic inhomogeneities (e.g., galaxies) in such a way that it can solve the problem of the flatness of the rotation curves. We also show that the fluctuations of radiation are caused by both the inhomogeneities in the form of galaxies and the fluctuations of quark-gluon nuggets. Therefore, if QNs exist, the CMB anisotropy should contain also the contributions from QNs. Additionally, the spatial distribution of the radiation fluctuations is defined by the gravitational potential. All these results look physically reasonable.

We calculate the present expansion of our Universe endowed with relict colored objects - quarks and gluons - that survived hadronization either as isolated islands of quark-gluon "nuggets", or spread uniformly in the Universe. In the first scenario, the QNs can play the role of dark matter. In the second scenario, we demonstrate that uniform colored objects can play the role of dark energy providing the late-time accelerating expansion of the Universe.