We study the effects associated with nonlinearity of $f(R)$ gravity and of the background perfect fluid manifested in the Kaluza-Klein model with spherical compactification. The background space-time is perturbed by a massive gravitating source which is pressureless in the external space but has an arbitrary equation of state (EoS) parameter in the internal space. As characteristics of a nonlinear perfect fluid, the squared speeds of sound are not equal to the background EoS parameters in the external and internal spaces. In this setting, we find exact solutions to the linearized Einstein equations for the perturbed metric coefficients. For nonlinear models with $f^{\prime\prime}(R_0)\neq0$, we show that these coefficients acquire correction terms in the form of two summed Yukawa potentials and that in the degenerated case, the solutions are reduced to a single Yukawa potential with some "corrupted" prefactor (in front of the exponential function), which, in addition to the standard $1/r$ term, contains a contribution independent of the three-dimensional distance $r$. In the linear $f''(R)=0$ model, we generalize the previous studies to the case of an arbitrary nonlinear perfect fluid. We also investigate the particular case of the nonlinear background perfect fluid with zero speed of sound in the external space and demonstrate that a non-trivial solution exists only in the case of $f''(R_0)=0$.

Two distinct perturbative approaches have been recently formulated within General Relativity, arguing for the screening of gravity in the $\Lambda$CDM Universe. In this paper we compare them and show that the offered screening concepts, each characterized by its own interaction range, can peacefully coexist. Accordingly, we advance a united scheme, determining the gravitational potential at all scales, including regions of nonlinear density contrasts, by means of a simple Helmholtz equation with the effective cosmological screening length. In addition, we claim that cosmic structures may not grow at distances above this Yukawa range and confront its current value with dimensions of the largest known objects in the Universe.

It is shown that the holographic principle of quantum gravity (in the hologram of the Universe, and therefore in our Galaxy, and of course on the Sun!), in which the conflict between the theory of gravitation and quantum mechanics disappears, gives rise to the Babcock-Leighton holographic mechanism. Unlike the solar dynamo models, it generates a strong toroidal magnetic field by means of the thermomagnetic Ettingshausen-Nernst (EN) effect in the tachocline. Hence, it can be shown that with the help of the thermomagnetic EN effect, a simple estimate of the magnetic pressure of an ideal gas in the tachocline of e.g. the Sun can indirectly prove that by using the holographic principle of quantum gravity, the repulsive toroidal magnetic field of the tachocline ($B_{tacho}^{Sun} = 4.1 \cdot 10^7 ~G = - B_{core}^{Sun}$) precisely "neutralizes" the magnetic field in the Sun core, since the projections of the magnetic fields in the tachocline and the core have equal values but opposite directions. The basic problem is a generalized problem of the antidynamo model of magnetic flux tubes (MFTs), where the nature of both holographic effects (the thermomagnetic EN~effect and Babcock-Leighton holographic mechanism), including magnetic cycles, manifests itself in the modulation of asymmetric dark matter (ADM) and, consequently, the solar axion in the Sun interior.

Within the cosmic screening approach, we obtain the exact formulas for the velocity-independent gravitational potentials produced by matter in the form of discrete sources distributed in the open and closed Universes. These formulas demonstrate that spatial curvature of the Universe considerably affect the form of the potentials and forces. While in the open Universe the gravitational force undergoes exponential suppression at cosmological distances, in the closed Universe the force induced by an individual mass is equal to zero at the antipodal point with respect to this mass. The derived formulas are applicable for investigations of the motion of astrophysical objects (e.g., galaxies) in the open and closed Universes, and for simulations of the large scale structure formation.

We study the second-order scalar perturbations in the conventional $\Lambda$CDM Universe within the cosmic screening approach. The analytic expressions for both the velocity-independent and velocity-dependent second-order scalar perturbations are derived. We demonstrate how the Yukawa screening effect, which is inherent in the first-order metric corrections, manifests itself in the second-order ones. It is shown that the obtained formulas for the second-order perturbations are reduced to the known post-Newtonian expressions at distances much smaller than the Yukawa screening length. In the era of precision cosmology, these analytic formulas play an important role since the second-order metric corrections may affect the interpretation of observational data (e.g., the luminosity-redshift relation, gravitational lensing, baryon acoustic oscillations).

We investigate cosmological perturbations for nonlinear $f(R)$ models within the cosmic screening approach. Matter is considered both in the form of a set of discrete point-like massive bodies and in the form of a continuous pressureless perfect fluid. We perform full relativistic analysis of the first-order theory of scalar perturbations for arbitrary nonlinear $f(R)$ models and demonstrate that scalar potentials $\Phi(t,\mathbf{r})$ and $\Psi(t,\mathbf{r})$ are determined by a system of only two master equations. Our equations are applicable at all spatial scales as long as the approximation $\delta R/\bar R \ll 1$ (which is usually assumed in studies devoted to cosmological perturbations in $f(R)$ models) works.

A high-energy photon polarimeter for astrophysics studies in the energy range from 10 MeV to 800 MeV is considered. The proposed concept uses a stack of silicon micro-strip detectors where they play the roles of both a converter and a tracker. The purpose of this paper is to outline the parameters of such a polarimeter and to estimate the productivity of measurements. Our study supported by a Monte Carlo simulation shows that with a one-year observation period the polarimeter will provide 6% accuracy of the polarization degree for photon energies above 100 MeV, which would be a significant advance relative to the currently explored energy range of a few MeV. The proposed polarimeter design could easily be adjusted to the specific photon energy range to maximize efficiency if needed.

We consider 5D brane world models with broken global 4D Poincar\a'{e} invariance (4D part of the spacetime metric is not conformal to the Minkowski spacetime). The bulk is filled with the negative cosmological constant and may contain a perfect fluid. In the case of empty bulk (the perfect fluid is absent), it is shown that one brane solution always has either a physical or a coordinate singularity in the bulk. We cut off these singularities in the case of compact two brane model and obtain regular exact solutions for both 4D Poincar\a'{e} broken and restored invariance. When the perfect fluid is present in the bulk, we get the master equation for the metric coefficients in the case of arbitrary bulk perfect fluid equation of state (EoS) parameters. In two particular cases of EoS, we obtain the analytic solutions for thin and thick branes. First one generalizes the well known Randall-Sundrum model with one brane to the case of the bulk anisotropic perfect fluid. In the second solution, the 4D Poincar\a'{e} invariance is restored. Here, the spacetime goes asymptotically to the anti-de Sitter one far from the thick brane.

We investigate dark and visible matter distribution in the Coma cluster in the case of the Navarro-Frenk-White (NFW) profile. A toy model where all galaxies in the cluster are concentrated inside a sphere of an effective radius $R_{\mathrm{eff}}$ is considered. It enables us to obtain the mean velocity dispersion as a function of $R_{\mathrm{eff}}$. We show that, within the observation accuracy of the NFW parameters, the calculated value of $R_{\mathrm{eff}}$ can be rather close to the observable cutoff of the galaxy distribution. Moreover, the comparison of our toy model with the observable data and simulations leads to the following preferable NFW parameters for the Coma cluster: $R_{200} \approx 1.77\,h^{-1} \, \mathrm{Mpc} = 2.61\, \mathrm{Mpc}$, $c=3\div 4$ and $M_{200}= 1.29h^{-1}\times10^{15}M_{\odot}$. In the Coma cluster the most of galaxies are concentrated inside a sphere of the effective radius $R_{\mathrm{eff}}\sim 3.7$ Mpc and the line-of-sight velocity dispersion is $1004\, \mathrm{km}\, \mathrm{s}^{-1}$.

We study a universe filled with cold dark matter in the form of discrete inhomogeneities (e.g., galaxies) and dark energy in the form of a continuous perfect fluid. We develop a first-order scalar perturbation theory in the weak gravity limit around a spatially flat Friedmann universe. Our approach works at all cosmic scales and incorporates linear and nonlinear effects with respect to energy density fluctuations. The gravitational potential can be split into individual contributions from each matter source. Each potential is characterized by a Yukawa interaction with the same range, which is of the order of 3700 Mpc at the present time. The derived equations can form the theoretical basis for numerical simulations for a wide class of modern cosmological models.

We make a generalization of a self-consistent first-order perturbation scheme, being suitable for all (sub-horizon and super-horizon) scales, which has been recently constructed for the concordance cosmological model and discrete presentation of matter sources, to the case of extended models with extra perfect fluids and continuous presentation. Namely, we derive a single equation determining the scalar perturbation and covering the whole space as well as define the corresponding universal Yukawa interaction range. We also demonstrate explicitly that the structure growth is suppressed at distances exceeding this fundamental range.

In the $\Lambda$CDM framework, presenting nonrelativistic matter inhomogeneities as discrete massive particles, we develop the second-order cosmological perturbation theory. Our approach relies on the weak gravitational field limit. The derived equations for the second-order scalar, vector and tensor metric corrections are suitable at arbitrary distances including regions with nonlinear contrasts of the matter density. We thoroughly verify fulfilment of all Einstein equations as well as self-consistency of order assignments. In addition, we achieve logical positive results in the Minkowski background limit. Feasible investigations of the cosmological backreaction manifestations by means of relativistic simulations are also outlined.

Without exceeding the limits of the conventional $\Lambda$CDM paradigm, we argue for Yukawa law of interparticle interaction as the law of gravitation in the real expanding inhomogeneous Universe. It covers the whole space and comes up to take place of Newtonian gravity, which is restricted exclusively to sub-horizon distances. The large-scale screening of gravitational interaction between every two nonrelativistic massive particles is ensured by the homogeneous cosmological background (specifically, by the nonzero average rest mass density of nonrelativistic matter). We take advantage of the uniform matter distribution case (i.e. the homogeneous Universe limit) to demonstrate superiority of Yukawa gravity. Attention is also devoted to the concrete particular case of inhomogeneity.

We study a universe filled with dust-like matter in the form of discrete inhomogeneities (e.g., galaxies and their groups and clusters) and two sets of perfect fluids with linear and nonlinear equations of state, respectively. The background spacetime geometry is defined by the FLRW metric. In the weak gravitational field limit, we develop the first-order scalar and vector cosmological perturbation theory. Our approach works at all cosmological scales (i.e. sub-horizon and super-horizon ones) and incorporates linear and nonlinear effects with respect to energy density fluctuations. We demonstrate that the scalar perturbation (i.e. the gravitational potential) as well as the vector perturbation can be split into individual contributions from each matter source. Each of these contributions satisfies its own equation. The velocity-independent parts of the individual gravitational potentials are characterized by a finite time-dependent Yukawa interaction range being the same for each individual contribution. We also obtain the exact form of the gravitational potential and vector perturbation related to the discrete matter sources. The self-consistency of our approach is thoroughly checked. The derived equations can form the theoretical basis for numerical simulations for a wide class of cosmological models.

In the framework of the concordance cosmological model the first-order scalar and vector perturbations of the homogeneous background are derived in the weak gravitational field limit without any supplementary approximations. The sources of these perturbations (inhomogeneities) are presented in the discrete form of a system of separate point-like gravitating masses. The found expressions for the metric corrections are valid at all (sub-horizon and super-horizon) scales and converge at all points except at locations of the sources. The average values of these metric corrections are zero (thus, first-order backreaction effects are absent). Both the Minkowski background limit and the Newtonian cosmological approximation are reached under certain well-defined conditions. An important feature of the velocity-independent part of the scalar perturbation is revealed: up to an additive constant this part represents a sum of Yukawa potentials produced by inhomogeneities with the same finite time-dependent Yukawa interaction range. The suggested connection between this range and the homogeneity scale is briefly discussed along with other possible physical implications.

It is shown that the hypothesis of the axion mechanism for solar luminosity variations suggesting that the solar axion particles are born in the core of the Sun and may be efficiently converted back into $\gamma$-quanta in the magnetic field of the solar overshoot tachocline is physically relevant. As a result, it is also shown that the intensity variations of the $\gamma$-quanta of axion origin, induced by the magnetic field variations in the tachocline via the thermomagnetic Ettingshausen-Nernst effect, directly cause the Sun luminosity variations and eventually characterize the active and quiet states of the Sun. Within the framework of this mechanism estimations of the strength of the axion coupling to a photon ($g_{a \gamma} = 4.4 \cdot 10^{-11} GeV^{-1}$) and the hadronic axion particle mass ($m_a \sim 3.2 \cdot 10^{-2} eV$) have been obtained. It is also shown that the claimed axion parameters do not contradict any known experimental and theoretical model-independent limitations.

The conventional $\Lambda$CDM cosmological model supplemented by the inflation concept describes the Universe very well. However, there are still a few concerns: new Planck data impose constraints on the shape of the inflaton potential, which exclude a lot of inflationary models; dark matter is not detected directly, and dark energy is not understood theoretically on a satisfactory level. In this brief sketch we investigate an alternative cosmological model with spherical spatial geometry and an additional perfect fluid with the constant parameter $\omega=-1/3$ in the linear equation of state. It is demonstrated explicitly that in the framework of such a model it is possible to satisfy the supernovae data at the same level of accuracy as within the $\Lambda$CDM model and at the same time suppose that the observed cosmic microwave background (CMB) radiation originates from a very limited space region. This is ensured by introducing an additional condition of light propagation between the antipodal points during the age of the Universe. Consequently, the CMB uniformity can be explained without the inflation scenario. The corresponding drawbacks of the model with respect to its comparison with the CMB data are also discussed.

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.

Scalar cosmological perturbations are investigated in the framework of a model with interacting dark energy and dark matter. In addition to these constituents, the inhomogeneous Universe is supposed to be filled with the standard noninteracting constituents corresponding to the conventional $\Lambda$CDM model. The interaction term is chosen in the form of a linear combination of dark sector energy densities with evolving coefficients. The methods of discrete cosmology are applied, and strong theoretical constraints on the parameters of the model are derived. A brief comparison with observational data is performed.

In the framework of a cosmological model of the Universe filled with a nonrelativistic particle soup, we easily reproduce inflation due to the quantum potential. The lightest particles in the soup serve as a driving force of this simple, natural and promising mechanism. It is explicitly demonstrated that the appropriate choice of their mass and fraction leads to reasonable numbers of e-folds. Thus, the direct introduction of the quantum potential into cosmology of the earliest Universe gives ample opportunities of successful reconsideration of the modern inflationary theory.

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 the Universe contains also the cosmological constant and a perfect fluid with a negative constant equation of state (EoS) parameter $\omega$ (e.g., quintessence, phantom or frustrated network of topological defects), we investigate scalar perturbations of the FRW metrics due to inhomogeneities. Our analysis shows that, to be compatible with the theory of scalar perturbations, this perfect fluid, first, should be clustered and, second, should have the equation of state parameter $\omega=-1/3$. In particular, this value corresponds to the frustrated network of cosmic strings. Therefore, the frustrated network of domain walls with $\omega =-2/3$ is ruled out. A perfect fluid with $\omega =-1/3$ neither accelerates nor decelerates the Universe. We also obtain the equation for the nonrelativistic gravitational potential created by a system of inhomogeneities. Due to the perfect fluid with $\omega = -1/3$, the physically reasonable solutions take place for flat, open and closed Universes. This perfect fluid is concentrated around the inhomogeneities and results in screening of the gravitational potential.

A striking coincidence of revolution periods of S-stars orbiting a supermassive black hole at the Galactic Center of the Milky Way and oscillation periods of such solar and terrestrial observables as the sunspot number, the geomagnetic field Y-component and the global temperature is established on basis of the corresponding experimental data. Rejecting randomness of this discovered coincidence, we put forward a hypothesis that modulation of dark matter flows in the Milky Way by the S-stars is responsible for such a frequency transfer from the Galactic Center to the Solar System.

This paper is devoted to different modifications of two standard softenings of the gravitational attraction (namely the Plummer and Hernquist softenings), which are commonly used in cosmological simulations based on the particle-particle (PP) method, and their comparison. It is demonstrated that some of the proposed alternatives lead to almost the same accuracy as in the case of the pure Newtonian interaction, even despite the fact that the force resolution is allowed to equal half the minimum interparticle distance. The revealed way of precision improvement gives an opportunity to succeed in solving Gurzadyan's Problem 5 and bring modern computer codes up to a higher standard.

It is shown that the hypothesis of the axion mechanism of Sun luminosity suggesting that the solar axion particles are born in the core of the Sun and may be efficiently converted back into $\gamma$-quanta in the magnetic field of the solar overshoot tachocline is physically relevant. As a result, it is also shown that the intensity variations of the $\gamma$-quanta of axion origin, induced by the magnetic field variations in the tachocline, directly cause the Sun luminosity %and total solar irradiance (TSI) variations and eventually characterize the active and quiet states of the Sun. Within the framework of this mechanism estimations of the strength of the axion coupling to a photon ($g_{a \gamma} = 3.6 \cdot 10^{-11} GeV^{-1}$) and the hadronic axion particle mass ($m_a \sim 2.3 \cdot 10^{-2} eV$) have been obtained. It is also shown that the claimed axion parameters do not contradict any known experimental and theoretical model-independent limitations.

In this paper we study scalar perturbations of the metric for nonlinear $f(R)$ models. We consider the Universe at the late stage of its evolution and deep inside the cell of uniformity. We investigate the astrophysical approach in the case of Minkowski spacetime background and two cases in the cosmological approach, the large scalaron mass approximation and the quasi-static approximation, getting explicit expressions for scalar perturbations for both these cases. In the most interesting quasi-static approximation, the scalar perturbation functions depend on both the nonlinearity function $f(R)$ and the scale factor $a$. Hence, we can study the dynamical behavior of the inhomogeneities (e.g., galaxies and dwarf galaxies) including into consideration their gravitational attraction and the cosmological expansion, and also taking into account the effects of nonlinearity. Our investigation is valid for functions $f(R)$ which have stable de Sitter points in future with respect to the present time, that is typical for the most popular $f(R)$ models.

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 consider Kaluza-Klein (KK) models where internal spaces are compact Einstein spaces. These spaces are stabilized by background matter (e.g., monopole form-fields). We perturb this background by a compact matter source (e.g., the system of gravitating masses) with the zero pressure in the external/our space and an arbitrary pressure in the internal space. We show that the Einstein equations are compatible only if the matter source is smeared over the internal space and perturbed metric components do not depend on coordinates of extra dimensions. The latter means the absence of KK modes corresponding to the metric fluctuations. Maybe, the absence of KK particles in LHC experiments is explained by such mechanism.

We consider the Universe deep inside the cell of uniformity. At these scales, the Universe is filled with inhomogeneously distributed discrete structures (galaxies, groups and clusters of galaxies), which perturb the background Friedmann model. Here, the mechanical approach (Eingorn & Zhuk, 2012) is the most appropriate to describe the dynamics of the inhomogeneities which is defined, on the one hand, by gravitational potentials of inhomogeneities and, on the other hand, by the cosmological expansion of the Universe. In this paper, we present additional arguments in favor of this approach. First, we estimate the size of the cell of uniformity. With the help of the standard methods of statistical physics and for the galaxies of the type of the Milky Way and Andromeda, we get that it is of the order of 190 Mpc which is rather close to observations. Then, we show that the nonrelativistic approximation (with respect to the peculiar velocities) is valid for $z \lesssim 10$, i.e. approximately for 13 billion years from the present moment. We consider scalar perturbations and, within the $\Lambda$CDM model, justify the main equations. Moreover, we demonstrate that radiation can be naturally incorporated into our scheme. This emphasizes the viability of our approach. This approach gives a possibility to analyze different cosmological models and compare them with the observable Universe. For example, we indicate some problematic aspects of the spatially flat models. Such models require a rather specific distribution of the inhomogeneities to get a finite potential at any points outside gravitating masses. We also criticize the application of the Schwarzschild-de Sitter solution to the description of the motion of test bodies on the cosmological background.

In this paper, we consider a system of gravitating bodies in Kaluza-Klein models with toroidal compactification of extra dimensions. To simulate the astrophysical objects (e.g., our Sun and pulsars) with energy density much greater than pressure, we suppose that these bodies are pressureless in the external/our space. At the same time, they may have nonzero parameters \omega_{(\bar\alpha -3)} \, (\bar\alpha =4,\ldots,D) of the equations of state in the extra dimensions. We construct the Lagrange function of this many-body system for any value of \Sigma =\sum_{\bar\alpha} \omega_{(\bar\alpha -3)}. Moreover, the gravitational tests (PPN parameters, perihelion and periastron advances) require negligible deviation from the latent soliton value \Sigma =-(D-3)/2. However, the presence of pressure/tension in the internal space results necessarily in the smearing of the gravitating masses over the internal space and in the absence of the KK modes. This looks very unnatural from the point of view of quantum physics.

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.

We consider a system of gravitating bodies in Kaluza-Klein models with toroidal compactification of extra dimensions. To simulate the astrophysical objects (e.g., our Sun and pulsars) with the energy density much greater than the pressure, we suppose that these bodies are pressureless in the external/our space. At the same time, they may have nonzero parameters \omega_{(\bar\alpha -3)} \, (\bar\alpha =4,...,D) of the equations of state in the extra dimensions. We construct the Lagrange function of this many-body system for any value of \Sigma =\sum_{\bar\alpha} \omega_{(\bar\alpha -3)}. Moreover, the gravitational tests (PPN parameters, perihelion/periastron advance) require negligible deviation from the latent soliton value \Sigma =-(D-3)/2. However, the presence of pressure/tension in the internal space results necessarily in the smearing of the gravitating masses over the internal space and in the absence of the KK modes. This looks very unnatural from the point of quantum physics.

We show the existence of the strong inverse correlation between the temporal variations of the toroidal component of the magnetic field in the solar tachocline (the bottom of the convective zone) and the Earth magnetic field (the Y-component). The possibility that the hypothetical solar axions, which can transform into photons in external electric or magnetic fields (the inverse Primakoff effect), can be the instrument by which the magnetic field of the Sun convective zone modulates the magnetic field of the Earth is considered. We propose the axion mechanism of Sun luminosity and "solar dynamo -- geodynamo" connection, where the energy of one of the solar axion flux components emitted in M1 transition in 57Fe nuclei is modulated at first by the magnetic field of the solar tachocline zone (due to the inverse coherent Primakoff effect) and after that is resonantly absorbed in the core of the Earth, thereby playing the role of the energy modulator of the Earth magnetic field. Within the framework of this mechanism estimations of the strength of the axion coupling to a photon (g_{a \gamma} = 7.07e-11 1/GeV), the axion-nucleon coupling (g_{an} = 3.20e-7), the axion-electron coupling (g_{ae} = 5.28e-11) and the axion mass (m_a = 17eV) have been obtained. It is also shown that the claimed axion parameters do not contradict to any known experimental and theoretical model-independent limitations. We consider the effect of dark matter in the form of 17eV axions on the extragalactic back-ground light (EBL) and find that the theoretical spectral intensity of the EBL from decaying axions as a function of the observed wavelength is in good agreement with the known experimental data for the near ultraviolet, optical and near infrared bands (1500-20000 angstrom).

We consider Kaluza-Klein models with background matter in the form of a multicomponent perfect fluid. This matter provides spherical compactification of the internal space with an arbitrary number of dimensions. The gravitating source has the dust-like equation of state in the external/our space and an arbitrary equation of state (with the parameter $\Omega$) in the internal space. In the single-component case, tension ($\Omega=-1/2$) is the necessary condition to satisfy both the gravitational tests in the solar system and the thermodynamical observations. In the multicomponent case, we propose two models satisfying both of these observations. One of them also requires tension $\Omega=-1/2$, but the second one is of special interest because is free of tension, i.e. $\Omega=0$. To get this result, we need to impose certain conditions.

The Workshop on results of the Project Kosmomikrofizyka-2 (Astroparticle Physics) of the National Academy of Sciences (NAS) of Ukraine "Astrophysical and cosmological problems of invisible mass and dark energy in the Universe" was held on November 21-22, 2012 in the Institute for Nuclear Research, Kyiv, Ukraine (http://lpd.kinr.kiev.ua/kmf12). This Project was carried out during three years (2010-2012) by scientists from various universities and institutes of the National Academy of Sciences of Ukraine; it was a logical continuation of the previous scientific program of the NAS of Ukraine "Researches of structure and composition of the Universe, hidden mass and dark energy (Kosmomikrofizyka)" in 2007-2009. These programs were devoted to theoretical and experimental investigations in astronomy, astrophysics, cosmology, physics of atomic nuclei and particle physics, which are related with the problems of dark matter and dark energy in the Universe.

In this paper we consider dynamical behavior of astrophysical objects (galaxies and dwarf galaxies) taking into account both gravitational attraction between them and cosmological expansion of the Universe. First, we obtain the general system of equations and apply them to some abstract systems of galaxies. Then we investigate the collision between the Milky Way and Andromeda in future. We distinguish two models. For the first one, we do not take into account the influence of Intra-Group Matter (IGrM). In this case we demonstrate that for currently known parameters of this system the collision is hardly plausible because of the angular momentum. These galaxies will approach the minimum distance of about 290 Kpc in 4.44 Gyr from present and then begin to run away irreversibly from each other. For the second model, we take into account dynamical friction due to IGrM. We find a characteristic value of the IGrM particle velocity dispersion $\tilde\sigma=2.306$. For $\tilde\sigma\leq2.306$ the merger will take place, but for the bigger values of $\tilde\sigma$ the merger can be problematic. If the temperature of the IGrM particles is $10^5$ K, this characteristic value of $\tilde\sigma$ corresponds to the IGrM particle mass 17 MeV. Therefore, for the IGrM particles with masses less than 17 MeV the merger becomes problematic. We also define the region in the vicinity of our Local Group where the formation of Hubble flows starts. For such processes the zero-acceleration surface (where gravitational attraction is balanced by cosmological accelerated expansion) plays the crucial role. We show that such surface is absent for the Local Group. Instead, we find two points and one circle with zero acceleration. Nevertheless, there is a nearly closed area around the MW and M31 where the absolute value of the acceleration is approximately equal to zero. Hubble flows are formed outside of this area.

In this letter, we consider the six-dimensional Kaluza-Klein models with spherical compactification of the internal space. Here, we investigate the case of bare gravitating compact objects with the dustlike equation of state $\hat p_0=0$ in the external (our) space and an arbitrary equation of state $\hat p_1=\Omega \hat\varepsilon$ in the internal space, where $\hat \varepsilon$ is the energy density of the source. This gravitating mass is spherically symmetric in the external space and uniformly smeared over the internal space. In the weak field approximation, the conformal variations of the internal space volume generate the admixture of the Yukawa potential to the usual Newton's gravitational potential. For sufficiently large Yukawa masses, such admixture is negligible and the metric coefficients of the external spacetime coincide with the corresponding expressions of General Relativity. Then, these models satisfy the classical gravitational tests. However, we show that gravitating masses acquire effective relativistic pressure in the external space. Such pressure contradicts the observations of compact astrophysical objects (e.g., the Sun). The equality $\Omega =-1/2$ (i.e. tension) is the only possibility to preserve the dustlike equation of state in the external space. Therefore, in spite of agreement with the gravitational experiments for an arbitrary value of $\Omega$, tension ($\Omega=-1/2$) plays a crucial role for the considered models.

In this paper, we consider the Universe deep inside of the cell of uniformity. At these scales, the Universe is filled with inhomogeneously distributed discrete structures (galaxies, groups and clusters of galaxies), which disturb the background Friedmann model. We propose mathematical models with conformally flat, hyperbolic and spherical spaces. For these models, we obtain the gravitational potential for an arbitrary number of randomly distributed inhomogeneities. In the cases of flat and hyperbolic spaces, the potential is finite at any point, including spatial infinity, and valid for an arbitrary number of gravitating sources. For both of these models, we investigate the motion of test masses (e.g., dwarf galaxies) in the vicinity of one of the inhomogeneities. We show that there is a distance from the inhomogeneity, at which the cosmological expansion prevails over the gravitational attraction and where test masses form the Hubble flow. For our group of galaxies, it happens at a few Mpc and the radius of the zero-acceleration sphere is of the order of 1 Mpc, which is very close to observations. Outside of this sphere, the dragging effect of the gravitational attraction goes very fast to zero.

We investigate Kaluza-Klein models in the case of spherical compactification of the internal space with an arbitrary number of dimensions. The gravitating source has the dust-like equation of state in the external/our space and an arbitrary equation of state (with the parameter $\Omega$) in the internal space. We get the perturbed (up to $O(1/c^2)$) metric coefficients. For the external space, these coefficients consist of two parts: the standard general relativity expressions plus the admixture of the Yukawa interaction. This admixture takes place only for some certain condition which is equivalent to the condition for the internal space stabilization. We demonstrate that the mass of the Yukawa interaction is defined by the mass of the gravexciton/radion. In the Solar system, the Yukawa mass is big enough for dropping the admixture of this interaction and getting good agreement with the gravitational tests for any value of $\Omega$. However, the gravitating body acquires the effective relativistic pressure in the external space which vanishes only in the case of tension $\Omega=-1/2$ in the internal space.

In these remarks, we clarify the problematic aspects of gravitational interaction in a weak-field limit of Kaluza-Klein models. We explain why some models meet the classical gravitational tests, while others do not. We show that variation of the total volume of the internal spaces generates the fifth force. This is the main reason of the problem. It happens for all considered models (linear with respect to the scalar curvature and nonlinear $f(R)$, with toroidal and spherical compactifications). We explicitly single out the contribution of the fifth force to nonrelativistic gravitational potentials. In the case of models with toroidal compactification, we demonstrate how tension (with and without effects of nonlinearity) of the gravitating source can fix the total volume of the internal space, resulting in the vanishing fifth force and consequently in agreement with the observations. It takes place for latent solitons, black strings and black branes. We also demonstrate a particular example where non-vanishing variations of the internal space volume do not contradict the gravitational experiments. In the case of spherical compactification, the fifth force is replaced by the Yukawa interaction for models with the stabilized internal space. For large Yukawa masses, the effect of this interaction is negligibly small, and considered models satisfy the gravitational tests at the same level of accuracy as general relativity.

In the six-dimensional Kaluza-Klein model with the multidimensional cosmological constant $\Lambda_6$, we obtain the black brane with spherical compactification of the internal space. The matter source for this exact solution consists of two parts. First, it is a fine-tuned homogeneous perfect fluid which provides spherical compactification of the internal space. Second, it is a gravitating massive body with the dustlike equation of state in the external space and tension $\hat p_1=-(1/2)\hat\varepsilon$ in the internal space. This solution exists both in the presence and absence of $\Lambda_6$. In the weak-field approximation, we also get solutions of the linearized Einstein equations for the model with spherical compactification. Here, the gravitating matter source has the dustlike equation of state in the external space and an arbitrary equation of state $\hat p_1=\Omega \hat\varepsilon$ in the internal space. In the case $\Lambda_6>0$ and $\Omega\neq -1/2$, these approximate solutions tend asymptotically to the weak-field limit of the exact black brane solution. Both the exact and asymptotic black branes satisfy the gravitational experiments at the same level of accuracy as general relativity.

We investigate the classical gravitational tests for the six-dimensional Kaluza-Klein model with spherical (of a radius $a$) compactification of the internal space. The model contains also a bare multidimensional cosmological constant $\Lambda_6$. The matter, which corresponds to this ansatz, can be simulated by a perfect fluid with the vacuum equation of state in the external space and an arbitrary equation of state with the parameter $\omega_1$ in the internal space. For example, $\omega_1=1$ and $\omega_1=2$ correspond to the monopole two-forms and the Casimir effect, respectively. In the particular case $\Lambda_6=0$, the parameter $\omega_1$ is also absent: $\omega_1=0$. In the weak-field approximation, we perturb the background ansatz by a point-like mass. We demonstrate that in the case $\omega_1>0$ the perturbed metric coefficients have the Yukawa type corrections with respect to the usual Newtonian gravitational potential. The inverse square law experiments restrict the parameters of the model: $a/\sqrt{\omega_1}\lesssim 6\times10^{-3}\ {{cm}}$. Therefore, in the Solar system the parameterized post-Newtonian parameter $\gamma$ is equal to 1 with very high accuracy. Thus, our model satisfies the gravitational experiments (the deflection of light and the time delay of radar echoes) at the same level of accuracy as General Relativity. We demonstrate also that our background matter provides the stable compactification of the internal space in the case $\omega_1>0$. However, if $\omega_1=0$, then the parameterized post-Newtonian parameter $\gamma=1/3$, which strongly contradicts the observations.

We investigate classical gravitational tests for the Kaluza-Klein model with spherical compactification of the internal two-dimensional space. In the case of the absence of a multidimensional bare cosmological constant, the only matter which corresponds to the proposed metric ansatz is a perfect fluid with the vacuum equation of state in the external space and the dust-like equation of state in the internal space. We perturb this background by a compact massive source with the dust-like equation of state in both external and internal spaces (e.g., a point-like mass), and obtain the metric coefficients in the weak-field approximation. It enables to calculate the parameterized post-Newtonian parameter $\gamma$. We demonstrate that $\gamma=1/3$ which strongly contradicts the observations.

We investigate nonlinear f(R) theories in the Kaluza-Klein models with toroidal compactification of extra dimensions. A point-like matter source has the dust-like equation of state in our three dimensions and nonzero equations of state in the extra dimensions. We obtain solutions of linearized Einstein equations with this matter source taking into account effects of nonlinearity of the model. There are two asymptotic regions where solutions satisfy the gravitational tests at the same level of accuracy as General Relativity. According to these asymptotic regions, there are two classes of solutions. We call these solutions asymptotic latent solitons. The asymptotic latent solitons from the first class generalize the known result of the linear theory. The asymptotic black strings and black branes are particular cases of these asymptotic solutions. The second class of asymptotic solitons exists only in multidimensional nonlinear models. The main feature for both of these classes of solutions is that the matter sources have tension in the extra dimensions.

In the Randall-Sundrum model with one brane, we found the approximate and exact solutions for gravitational potentials and accelerations of test bodies in these potentials for different geometrical configurations. We applied these formulas for calculation of the gravitational interaction between two spheres and found the approximate and exact expressions for the relative force corrections to the Newton's gravitational force. We demonstrated that the difference between relative force corrections for the approximate and exact cases increases with the parameter $l$ (for the fixed distance $r$ between centers of the spheres). On the other hand, this difference increases with decreasing of the distance between the centers of the spheres (for the fixed curvature scale parameter $l$). We got the upper limit for the curvature scale parameter $l\lesssim 10\, \mu$m. For these values of $l$, the difference between the approximate and exact solutions is negligible.

We investigate a point-like massive source in non-linear f(R) theories in the case of arbitrary number of spatial dimensions D\geq 3. If D>3 then extra dimensions undergo toroidal compactification. We consider a weak-field approximation with Minkowski and de Sitter background solutions. In both these cases point-like massive sources demonstrate good agreement with experimental data only in the case of ordinary three-dimensional (D=3) space. We generalize this result to the case of perfect fluid with dust-like equations of state in the external and internal spaces. This perfect fluid is uniformly smeared over all extra dimensions and enclosed in a three-dimensional sphere. In ordinary three dimensional (D=3) space, our formulas are useful for experimental constraints on parameters of f(R) models.

In Kaluza-Klein models with an arbitrary number of toroidal internal spaces, we investigate soliton solutions which describe the gravitational field of a massive compact object. We single out the physically interesting solution corresponding to a point-like mass. For the general solution we obtain equations of state in the external and internal spaces. These equations demonstrate that the point-like mass soliton has dust-like equations of state in all spaces. We also obtain the PPN parameters, which give the possibility to obtain the formulas for perihelion shift, deflection of light and time delay of radar echoes. Additionally, the gravitational experiments lead to a strong restriction on the parameter of the model: $\tau = -(2.1\pm 2.3)\times 10^{-5}$. The point-like mass solution contradicts this restriction. The condition $\tau=0$ satisfies the experimental limitation and defines a new class of solutions which are indistinguishable from general relativity. We call such solutions latent solitons. Black strings and black branes belong to this class. Moreover, the condition of stability of the internal spaces singles out black strings/branes from the latent solitons and leads uniquely to the black string/brane equations of state $p_i=-\epsilon/2$, in the internal spaces and to the number of the external dimensions $d_0=3$. The investigation of multidimensional static spherically symmetric perfect fluid with dust-like equation of state in the external space confirms the above results.

In Kaluza-Klein models, we investigate soliton solutions of Einstein equation. We obtain the formulas for perihelion shift, deflection of light, time delay of radar echoes and PPN parameters. We find that the solitonic parameter k should be very big: |k|\geq 2.3\times10^4. We define a soliton solution which corresponds to a point-like mass source. In this case the soliton parameter k=2, which is clearly contrary to this restriction. Similar problem with the observations takes place for static spherically symmetric perfect fluid with the dust-like equation of state in all dimensions. The common for both of these models is the same equations of state in our three dimensions and in the extra dimensions. All dimensions are treated at equal footing. To be in agreement with observations, it is necessary to break the symmetry between the external/our and internal spaces. It takes place for black strings which are particular examples of solitons with k\to \infty. For such k, black strings are in concordance with the observations. Moreover, we show that they are the only solitons which are at the same level of agreement with the observations as in general relativity. Black strings can be treated as perfect fluid with dust-like equation of state p_0=0 in the external/our space and very specific equation of state p_1=-(1/2)\epsilon in the internal space. The latter equation is due to negative tension in the extra dimension. We also demonstrate that dimension 3 for the external space is a special one. Only in this case we get the latter equation of state. We show that the black string equations of state satisfy the necessary condition of the internal space stabilization. Therefore, black strings are good candidates for a viable model of astrophysical objects (e.g., Sun) if we can provide a satisfactory explanation of negative tension for particles constituting these objects.

It is found the exact solution of the Poisson equation for the multidimensional space with topology $M_{3+d}=\mathbb{R}^3\times T^d$. This solution describes smooth transition from the newtonian behavior $1/r_3$ for distances bigger than periods of tori (the extra dimension sizes) to multidimensional behavior $1/r^{1+d}_{3+d}$ in opposite limit. In the case of one extra dimension $d=1$, the gravitational potential is expressed via compact and elegant formula. It is shown that the corrections to the gravitational constant in the Cavendish-type experiment can be within the measurement accuracy of Newton's gravitational constant $G_N$. It is proposed models where the test masses are smeared over some (or all) extra dimensions. In 10-dimensional spacetime with 3 smeared extra dimensions, it is shown that the size of 3 rest extra dimensions can be enlarged up to submillimeter for the case of 1TeV fundamental Planck scale $M_{Pl(10)}$. In the models where all extra dimensions are smeared, the gravitational potential exactly coincides with the newtonian one. Nevertheless, the hierarchy problem can be solved in these models.

In the case of one extra dimension, well known Newton's potential $\phi (r_3)=-G_N m/r_3$ is generalized to compact and elegant formula $\phi(r_3,\xi)=-(G_N m/r_3)\sinh(2\pi r_3/a)[\cosh(2\pi r_3/a)-\cos(2\pi\xi/a)]^{-1}$ if four-dimensional space has topology $\mathbb{R}^3\times T$. Here, $r_3$ is magnitude of three-dimensional radius vector, $\xi$ is extra dimension and $a$ is a period of a torus $T$. This formula is valid for full range of variables $r_3 \in [0,+\infty)$ and $\xi\in [0,a]$ and has known asymptotic behavior: $\phi \sim 1/r_3$ for $r_3>>a$ and $\phi \sim 1/r_4^2$ for $r_4=\sqrt{r_3^2+\xi^2}<<a$. Obtained formula is applied to an infinitesimally thin shell, a shell, a sphere and two spheres to show deviations from the newtonian expressions. Usually, these corrections are very small to observe at experiments. Nevertheless, in the case of spatial topology $\mathbb{R}^3\times T^{d}$, experimental data can provide us with a limitation on maximal number of extra dimensions.