The evolution of galaxy clusters can be affected by the repulsion described by the cosmological constant. This conclusion is reached within the modified weak-field General Relativity approach where the cosmological constant \Lambda enables to describe the common nature of the dark matter and the dark energy. Geometrical methods of theory of dynamical systems and the Ricci curvature criterion are used to reveal the difference in the instability properties of galaxy clusters which determine their evolutionary paths. Namely, it is shown that the clusters determined by the gravity with \Lambda-repulsion tend to become even more unstable than those powered only by Newtonian gravity, the effect to be felt at cosmological time scales.

A framework is developed enabling the global analysis of the stability of cosmological models using the local geometric characteristics of the infinite-dimensional superspace, i.e. using the generalised Jacobi equation reformulated for pseudo-Riemannian manifolds. We give a direct formalism for dynamical analysis in the superspace, the requisite equation pertinent for stability analysis of the universe by means of generalized covariant and Fermi derivative is derived. Then, the relevant definitions and formulae are retrieved for cosmological models with a scalar field.

An analytical method for investigation of the evolution of dynamical systems {\it with independent on time accuracy} is developed for perturbed Hamiltonian systems. The error-free estimation using of computer algebra enables the application of the method to complex multi-dimensional Hamiltonian and dissipative systems. It also opens principal opportunities for the qualitative study of chaotic trajectories. The performance of the method is demonstrated on perturbed two-oscillator systems. It can be applied to various non-linear physical and astrophysical systems, e.g. to the long-term planetary dynamics.

The recently developed method (Paper 1) enabling one to investigate the evolution of dynamical systems with an accuracy not dependent on time is developed further. The classes of dynamical systems which can be studied by that method are much extended, now including systems that are; (1) non-Hamiltonian, conservative; (2) Hamiltonian with time-dependent perturbation; (3) non-conservative (with dissipation). These systems cover various types of N-body gravitating systems of astrophysical and cosmological interest, such as the orbital evolution of planets, minor planets, artificial satellites due to tidal, non-tidal perturbations and thermal thrust, evolving close binary stellar systems, and the dynamics of accretion disks.

Extensive N-body simulations are among the key means for the study of numerous astrophysical and cosmological phenomena, so various schemes are developed for possibly higher accuracy computations. We demonstrate the principal possibility for revealing the evolution of a perturbed Hamiltonian system with an accuracy independent on time. The method is based on the Laplace transform and the derivation and analytical solution of an evolution equation in the phase space for the resolvent and using computer algebra.

In arxiv:1108.5354 the Kolmogorov-Smirnov (K-S) test and Kolmogorov stochasticity parameter (KSP) is applied to CMB data. Their interpretation of the KSP method, however, lacks essential elements. In addition, their main result on the Gaussianity of CMB was not a matter of debate in previous KSP-CMB studies which also included predictions on cold spots, point sources.

This short note is concerned with a recent paper by N{\ae}ss [arXiv:1105.5051]. We explain why the statements in the paper are absolutely irrelevant.

In arXiv:1105.1081, Eriksen and Wehus try to criticise the use of Kolmogorov's stochasticity parameter (KSP) for CMB data analysis. Their discussion is based on a serious misunderstanding of the randomness and Arnold's work. Their note includes numerous further inaccuracies and groundless statements. However, this short note is concerned with KSP method only.

The cosmic microwave background (CMB) radiation is characterized by well-established scales, the 2.7 K temperature of the Planckian spectrum and the $10^{-5}$ amplitude of the temperature anisotropy. These features were instrumental in indicating the hot and equilibrium phases of the early history of the Universe and its large scale isotropy, respectively. We now reveal one more intrinsic scale in CMB properties. We introduce a method developed originally by Kolmogorov, that quantifies a degree of randomness (chaos) in a set of numbers, such as measurements of the CMB temperature in some region. Considering CMB as a composition of random and regular signals, we solve the inverse problem of recovering of their mutual fractions from the temperature sky maps. Deriving the empirical Kolmogorov's function in the Wilkinson Microwave Anisotropy Probe's maps, we obtain the fraction of the random signal to be about 20 per cent, i.e. the cosmological sky is a weakly random one. The paper is dedicated to the memory of Vladimir Arnold (1937-2010).

The Kolmogorov stochasticity parameter is shown to act as a tool to detect point sources in the cosmic microwave background (CMB) radiation temperature maps. Kolmogorov CMB map constructed for the WMAP's 7-year datasets reveals tiny structures which in part coincide with point radio and Fermi/LAT gamma-ray sources. In the first application of this method, we identified several sources not present in the then available 0FGL Fermi catalog. Subsequently they were confirmed in the more recent and more complete 1FGL catalog, thus strengthening the evidence for the power of this methodology.

The power spectrum is obtained for the Kolmogorov stochasticity parameter map for WMAP's cosmic microwave background (CMB) radiation temperature datasets. The interest for CMB Kolmogorov map is that it can carry direct information about voids in the matter distribution, so that the correlations in the distribution of voids have to be reflected in the power spectrum. Although limited by the angular resolution of the WMAP, this analysis shows the possibility of acquiring this crucial information via CMB maps. Even the already obtained behavior, some of which is absent in the simulated maps, can influence the development of views on the void correlations at the large-scale web formation.

The chaos in stellar systems is studied using the theory of dynamical systems and the Van Kampen stochastic differential equation approach. The exponential instability (chaos) of spherical N-body gravitating systems, already known previously, is confirmed. The characteristic timescale of that instability is estimated confirming the collective relaxation time obtained by means of the Maupertuis principle.

We study nonrelativistic gravity using the Hamiltonian formalism. For the dynamics of general relativity (relativistic gravity) the formalism is well known and called the Arnowitt-Deser-Misner (ADM) formalism. We show that if the lapse function is constrained correctly, then nonrelativistic gravity is described by a consistent Hamiltonian system. Surprisingly, nonrelativistic gravity can have solutions identical to relativistic gravity ones. In particular, (anti-)de Sitter black holes of Einstein gravity and IR limit of Horava gravity are locally identical.

The properties of geodesics flow are studied in a Friedmann-Robertson-Walker metric perturbed due to the inhomogeneities of matter. The basic, averaged Jacobi equation is derived, which reveals that the low density regions (voids) are able to induce hyperbolicity, even if the global curvature of the Universe is zero or slightly positive. It is shown that the energy independence is a characteristic property of these geometric effects. The importance of these conclusions is determined by the temperature independent ellipticity of excursion sets and regions of different randomness found in Kolmogorov CMB maps.

The plane-mirror symmetry previously noticed in the Cosmic Microwave Background (CMB) temperature anisotropy maps of Wilkinson Microwave Anisotropy Probe is shown to possess certain anomalous properties. The degree of the randomness determined by the Kolmogorov stochasticity parameter in the both symmetry regions appears to have identical values which, however, essentially differ from the corresponding values for other sky regions. If the mirroring were of cosmological origin, this would imply either additional randomizing properties in those directions of the Universe or their different line-of-sight depth. This analysis also provides a way to test the hypothesis of a link between the nature of dark energy and inhomogeneities.

A new map of the sky representing the degree of randomness in the cosmic microwave background (CMB) temperature has been obtained. The map based on estimation of the Kolmogorov stochasticity parameter clearly distinguishes the contribution of the Galactic disk from the CMB and reveals regions of various degrees of randomness that can reflect the properties of inhomogeneities in the Universe. For example, among the high randomness regions is the southern non-Gaussian anomaly, the Cold Spot, with a stratification expected for the voids. Existence of its counterpart, a Northern Cold Spot with almost identical randomness properties among other low-temperature regions is revealed. By its informative power, Kolmogorov's map can be complementary to the CMB temperature and polarization sky maps.

The degree of randomness in the cosmic microwave background (CMB) maps has been shown to be measurable and indicative of the hyperbolic effect of voids on the propagation of the photon beams. In terms of an introduced porosity parameter, we now obtain the criterion for a hyperbolicity due to domination of voids in the line-of-sight direction. The criterion seems to be supported by the parameters of 30 Mpc scale voids revealed by the galactic redshift surveys. The importance of the geometrical effect of the voids is also due to their possible role in the temperature independent ellipticity of the excursion sets, as well as in the nature of dark energy. The role of possible larger scale inhomogeneities, from superclusters to semi-Hubble scale, in this effect remains yet unclear.

The Kolmogorov stochasticity parameter (KSP) is applied to quantify the degree of randomness (stochasticity) in the temperature maps of the Cosmic Microwave Background radiation maps. It is shown that, the KSP for the WMAP5 maps is about twice higher than that of the simulated maps for the concordance \LambdaCDM cosmological model, implying the existence of a randomizing effect not taken into account in the model. As was revealed recently, underdense regions in the large scale matter distributions, i.e the voids, possess hyperbolic and hence randomizing properties. The degree of randomness for the Cold Spot appears to be about twice higher than the average over the mean temperature level spots in the sky, which supports the void nature of the Cold Spot. Kolmogorov's parameter then acts as a quantitative tracer of the voids via CMB.

We propose a geometrical approach to the investigation of Hamiltonian systems on (Pseudo) Riemannian manifolds. A new geometrical criterion of instability and chaos is proposed. This approach is more generic than well known reduction to the geodesic flow. It is applicable to various astrophysical and cosmological problems.

As was shown before (Gurzadyan and Kocharyan, 1992, 1993ab) the statistical properties (exponential mixing) of motion of CMB photon beams in Friedmann Universe with negative curvature can have definite observable consequences including the decrease of the CMB anisotropy amplitude after the last scattering, distortion of images on CMB sky maps. Here we consider one more consequence concerning the behaviour of temperature autocorrelation function C(\theta,\beta), where \theta is the modulation (separation) angle, and \beta the observing beam effective angle. We prove that in the Universe with negative curvature any spectrum of perturbations at last scattering epoch becomes a scale-invariant by present time. It is also predicted that the temperature autocorrelation function should increase while obtained from measurements with small beam sizes within some mixing beam angle.