We study the effect of the cubic torus topology of the Universe on scalar cosmological perturbations which define the gravitational potential. We obtain three alternative forms of the solution for both the gravitational potential produced by point-like masses, and the corresponding force. The first solution includes the expansion of delta-functions into Fourier series, exploiting periodic boundary conditions. The second one is composed of summed solutions of the Helmholtz equation for the original mass and its images. Each of these summed solutions is the Yukawa potential. In the third formula, we express the Yukawa potentials via Ewald sums. We show that for the present Universe, both the bare summation of Yukawa potentials and the Yukawa-Ewald sums require smaller numbers of terms to yield the numerical values of the potential and the force up to desired accuracy. Nevertheless, the Yukawa formula is yet preferable owing to its much simpler structure.

We study a three-component universe filled with dust-like matter in the form of discrete inhomogeneities (e.g., galaxies) and perfect fluids characterized by linear and nonlinear equations of state. Within the cosmic screening approach, we develop the theory of scalar and vector perturbations. None of the energy density contrasts associated with the distinct components is treated as small. Consequently, the derived equations are valid at both sub- and super-horizon scales and enable simulations for a variety of cosmological models.

We investigate the influence of the chimney topology $T\times T\times R$ of the Universe on the gravitational potential and force that are generated by point-like massive bodies. We obtain three distinct expressions for the solutions. One follows from Fourier expansion of delta functions into series using periodicity in two toroidal dimensions. The second one is the summation of solutions of the Helmholtz equation, for a source mass and its infinitely many images, which are in the form of Yukawa potentials. The third alternative solution for the potential is formulated via the Ewald sums method applied to Yukawa-type potentials. We show that, for the present Universe, the formulas involving plain summation of Yukawa potentials are preferable for computational purposes, as they require a smaller number of terms in the series to reach adequate precision.

Within the extension of the $$\varLambda $$Λ CDM model, allowing for the presence of neutrinos or warm dark matter, we develop the analytical cosmological perturbation theory. It covers all spatial scales where the weak gravitational field regime represents a valid approximation. Discrete particles – the sources of the inhomogeneous gravitational field – may be relativistic. Similarly to the previously investigated case of nonrelativistic matter, the Yukawa interaction range is naturally incorporated into the first-order scalar metric corrections.