We present a description of the new complex installation for the study of extensive air showers which was created at the Tien Shan mountain cosmic ray station, as well as the results of the test measurements made there in 2014-2016. At present, the system for registration of electromagnetic shower component consists of $\sim$100 detector points built on the basis of plastic scintillator plates with the sensitive area of 0.25m$^2$ and 1m$^2$, spread equidistantly over $\sim$10$^4$m$^2$ space. The dynamic range of scintillation amplitude measurements is currently about $(3-7)\cdot 10^4$, and there is a prospect of it being extended up to $\sim$10$^6$. The direction of shower arrival is defined by signal delays from a number of the scintillators placed cross-wise at the periphery of the detector system. For the investigation of nuclear active shower components there was created a multi-tier 55m$^2$ ionization-neutron calorimeter with a sum absorber thickness of $\sim$1000g/cm$^2$, typical spatial resolution of the order of 10cm, and dynamic range of ionization measurement channel about $\sim$10$^5$. Also, the use of saturation-free neutron detectors is anticipated for registration of the high- and low-energy hadron components in the region of shower core. A complex of underground detectors is designed for the study of muonic and penetrative nuclear-active components of the shower. The full stack of data acquisition, detector calibration, and shower parameters restoration procedures are now completed, and the newly obtained shower size spectrum and lateral distribution of shower particles occur in agreement with conventional data. Future studies in the field of $10^{14}-10^{17}$eV cosmic ray physics to be held at the new shower installation are discussed.

Experimental complex (EC) NEVOD includes a number of unique experimental facilities for studies of main components of cosmic rays on the Earth's surface. The complex is used for the basic research of CR flux characteristics and their interactions in the energy range 10^15 - 10^19 eV, and for applied investigations directed to the development of methods of the muon diagnostics of the atmosphere and the Earth's magnetosphere and near-terrestrial space. To extend the experimental capabilities and raising the status of the installation to the Mega Science level, nowadays new large-scale detectors: array for the EAS registration - NEVOD-EAS, detector of atmospheric neutrons - URAN, and large-area coordinate-tracking detector - TREK, are being deployed around EC NEVOD. The description of new detectors and a common trigger system to ensure the joint operation together with other detectors of EC NEVOD are presented.

We applied the method of continuous wavelet-transform to high-quality time-frequency analysis to the sets of observations of relative sunspot numbers. Wavelet analysis of these data reveals the following pattern: at the same time there are several activity cycles whose periods vary widely from the quasi biennial up to the centennial period. These relatively low-frequency periodic variations of the solar activity gradually change the values of periods of different cycles in time. This phenomenon can be observed in every cycle of activity.

This paper is concerned with the figures of equilibrium of a self-gravitating ideal fluid with density stratification and a steady-state velocity field. As in the classical setting, it is assumed that the figures or their layers uniformly rotate about an axis fixed in space.

The statistical analysis and the spherical wavelet analysis of the SDSS DR7 quasars distribution and of the WMAP CMB anisotropy are performed. They revealed the qualitative agreement between the angular power spectrum of CMB and the angular power spectrum of the quasar distribution on the celestial sphere. The angular correlation function and the angular power spectrum of the quasar distribution may be described by the power laws. The large quasar groups are discovered and they form the fractal set: the relation between their angular size and a number of quasar groups with this size is characterized by a power-law with fractal dimension 2.08.

We use 1-dimensional numerical simulations to study spherical collapse in the f(R) gravity models. We include the nonlinear coupling of the gravitational potential to the scalar field in the theory and use a relaxation scheme to follow the collapse. We find an unusual enhancement in density near the virial radius which may provide observable tests of gravity. We also use the estimated collapse time to calculate the critical overdensity $\delta_c$ used in calculating the mass function and bias of halos. We find that analytical approximations previously used in the literature do not capture the complexity of nonlinear spherical collapse.

Modifications of general relativity provide an alternative explanation to dark energy for the observed acceleration of the universe. We calculate quasilinear effects in the growth of structure in f(R) models of gravity using perturbation theory. We find significant deviations in the bispectrum that depend on cosmic time, length scale and triangle shape. However the deviations in the reduced bispectrum Q for f(R) models are at the percent level, much smaller than the deviations in the bispectrum itself. This implies that three-point correlations can be predicted to a good approximation simply by using the modified linear growth factor in the standard gravity formalism. Our results suggest that gravitational clustering in the weakly nonlinear regime is not fundamentally altered, at least for a class of gravity theories that are well described in the Newtonian regime by the parameters G_eff and Phi/Psi. This approximate universality was also seen in the N-body simulation measurements of the power spectrum by Stabenau and Jain (2006), and in other recent studies based on simulations. Thus predictions for such modified gravity models in the regime relevant to large-scale structure observations may be less daunting than expected on first principles. We discuss the caveats that apply to such predictions.