The creation of a 3D map of the bulge using RRLyrae (RRL) is one of the main goals of the VVV(X) surveys. The overwhelming number of sources under analysis request the use of automatic procedures. In this context, previous works introduced the use of Machine Learning (ML) methods for the variable star classification. Our goal is the development and analysis of an automatic procedure, based on ML, for the identification of RRLs in the VVV Survey. This procedure will be use to generate reliable catalogs integrated over several tiles in the survey. After the reconstruction of light-curves, we extract a set of period and intensity-based features. We use for the first time a new subset of pseudo color features. We discuss all the appropriate steps needed to define our automatic pipeline: selection of quality measures; sampling procedures; classifier setup and model selection. As final result, we construct an ensemble classifier with an average Recall of 0.48 and average Precision of 0.86 over 15 tiles. We also make available our processed datasets and a catalog of candidate RRLs. Perhaps most interestingly, from a classification perspective based on photometric broad-band data, is that our results indicate that Color is an informative feature type of the RRL that should be considered for automatic classification methods via ML. We also argue that Recall and Precision in both tables and curves are high quality metrics for this highly imbalanced problem. Furthermore, we show for our VVV data-set that to have good estimates it is important to use the original distribution more than reduced samples with an artificial balance. Finally, we show that the use of ensemble classifiers helps resolve the crucial model selection step, and that most errors in the identification of RRLs are related to low quality observations of some sources or to the difficulty to resolve the RRL-C type given the date.

Machine learning algorithms are highly useful for the classification of time series data in astronomy in this era of peta-scale public survey data releases. These methods can facilitate the discovery of new unknown events in most astrophysical areas, as well as improving the analysis of samples of known phenomena. Machine learning algorithms use features extracted from collected data as input predictive variables. A public tool called Feature Analysis for Time Series (FATS) has proved an excellent workhorse for feature extraction, particularly light curve classification for variable objects. In this study, we present a major improvement to FATS, which corrects inconvenient design choices, minor details, and documentation for the re-engineering process. This improvement comprises a new Python package called "feets", which is important for future code-refactoring for astronomical software tools.

We present the study of nineteen low X-ray luminosity galaxy clusters (L$_X \sim$ 0.5--45 $\times$ $10^{43}$ erg s$^{-1}$), selected from the ROSAT Position Sensitive Proportional Counters (PSPC) Pointed Observations (Vikhlinin et al. 1998) and the revised version of Mullis et al. (2003) in the redshift range of 0.16 to 0.7. This is the introductory paper of a series presenting the sample selection, photometric and spectroscopic observations and data reduction. Photometric data in different passbands were taken for eight galaxy clusters at Las Campanas Observatory; three clusters at Cerro Tololo Interamerican Observatory; and eight clusters at the Gemini Observatory. Spectroscopic data were collected for only four galaxy clusters using Gemini telescopes. With the photometry, the galaxies were defined based on the star-galaxy separation taking into account photometric parameters. For each galaxy cluster, the catalogues contain the PSF and aperture magnitudes of galaxies within the 90\% completeness limit. They are used together with structural parameters to study the galaxy morphology and to estimate photometric redshifts. With the spectroscopy, the derived galaxy velocity dispersion of our clusters ranged from 507 km~s$^{-1}$ for [VMF98]022 to 775 km~s$^{-1}$ for [VMF98]097 with signs of substructure. Cluster membership has been extensively discussed taking into account spectroscopic and photometric redshift estimates. In this sense, members are the galaxies within a projected radius of 0.75 Mpc from the X-ray mission peak and with cluster centric velocities smaller than the cluster velocity dispersion or 6000 km~s$^{-1}$, respectively. These results will be used in forthcoming papers to study, among the main topics, the red cluster sequence, blue cloud and green populations; the galaxy luminosity function and cluster dynamics.

Using a cosmological N-body numerical simulation of the formation of a galaxy cluster- sized halo, we analyze the temporal evolution of its globular cluster population. We follow the dynamical evolution of 38 galactic dark matter halos orbiting in a galaxy cluster that at redshift z=0 has a virial mass of 1.71 * 10 ^14 Msol h^-1. In order to mimic both "blue" and "red" populations of globular clusters, for each galactic halo we select two different sets of particles at high redshift (z ~ 1), constrained by the condition that, at redshift z=0, their average radial density profiles are similar to the observed profiles. As expected, the general galaxy cluster tidal field removes a significant fraction of the globular cluster populations to feed the intracluster population. On average, halos lost approximately 16% and 29% of their initial red and blue globular cluster populations, respectively. Our results suggest that these fractions strongly depend on the orbital trajectory of the galactic halo, specifically on the number of orbits and on the minimum pericentric distance to the galaxy cluster center that the halo has had. At a given time, these fractions also depend on the current clustercentric distance, just as observations show that the specific frequencyof globular clusters S_N depends on their clustercentric distance.

Recently alternative approaches in cosmology seeks to explain the nature of dark matter as a direct result of the non-linear spacetime curvature due to different types of deformation potentials. In this context, a key test for this hypothesis is to examine the effects of deformation on the evolution of large scales structures. An important requirement for the fine analysis of this pure gravitational signature (without dark matter elements) is to characterize the position of a galaxy during its trajectory to the gravitational collapse of super clusters at low redshifts. In this context, each element in an gravitational N-body simulation behaves as a tracer of collapse governed by the process known as chaotic advection (or lagrangian turbulence). In order to develop a detailed study of this new approach we develop the COsmic LAgrangian TUrbulence Simulator (COLATUS) to perform gravitational N-body simulations based on Compute Unified Device Architecture (CUDA) for graphics processing units (GPUs). In this paper we report the first robust results obtained from COLATUS.

Although the theoretical understanding of the nonlinear gravitational clustering has greatly advanced in the last decades, in particular by the outstanding improvement on numerical N-body simulations, the physics behind this process is not fully elucidated. The main goal of this work is the study of the possibility of a turbulent-like physical process in the formation of structures, galaxies and clusters of galaxies, by the action of gravity alone. We use simulation data from the Virgo Consortium, in ten redshift snapshots (from 0 to 10). From this we identify galaxy-sized and cluster-sized dark matter haloes, by using a FoF algorithm and applying a boundedness criteria, and study the gravitational potential energy spectra. We find that the galaxy-sized haloes energy spectrum follows closely a Kolmogorov's power law, similar to the behaviour of dynamically turbulent processes in fluids. This means that the gravitational clustering of dark matter may admit a turbulent-like representation.

In this paper we present a study of the cosmic flow from the point of view of how clusterings at different dynamical regimes in an expanding universe evolve according to a `coarse-grained' partitioning of their ranked energy distribution. By analysing a Lambda-CDM cosmological simulation from the Virgo Project, we find that cosmic flows evolve in an orderly sense, when tracked from their coarse-grained energy cells, even when nonlinearities are already developed. We show that it is possible to characterize scaling laws for the Pairwise Velocity Distribution in terms of the energy cells, generally valid at the linear and nonlinear clustering regimes.

We present an analysis of the behaviour of the `coarse-grained' (`mesoscopic') rank partitioning of the mean energy of collections of particles composing virialized dark matter halos in a Lambda-CDM cosmological simulation. We find evidence that rank preservation depends on halo mass, in the sense that more massive halos show more rank preservation than less massive ones. We find that the most massive halos obey Arnold's theorem (on the ordering of the characteristic frequencies of the system) more frequently than less massive halos. This method may be useful to evaluate the coarse-graining level (minimum number of particles per energy cell) necessary to reasonably measure signatures of `mesoscopic' rank orderings in a gravitational system.

We investigate two important questions about the use of the nonextensive thermostatistics (NETS) formalism in the context of nonlinear galaxy clustering in the Universe. Firstly, we define a quantitative criterion for justifying nonextensivity at different physical scales. Then, we discuss the physics behind the ansatz of the entropic parameter $q(r)$. Our results suggest the approximate range where nonextensivity can be justified and, hence, give some support to the applicability of NETS to the study of large scale structures.

There has been a trend in the past decade to describe the large-scale structures in the Universe as a (multi)fractal set. However, one of the main objections raised by the opponents of this approach deals with the transition to homogeneity. Moreover, they claim there is not enough sampling space to determine a scaling index which characterizes a (multi)fractal set. In this work we propose an alternative solution to this problem, using the generalized thermostatistics formalism. We show that applying the idea of nonextensivity, intrinsic to this approach, it is possible to derive an expression for the correlation function, describing the scaling properties of large-scale structures in the Universe and the transition to homogeneity, which is in good agreement with observational data.