This work discusses the main analogies and differences between the deterministic approach underlying most cosmological N-body simulations and the probabilistic interpretation of the problem that is often considered in mathematics and statistical mechanics. In practice, we advocate for averaging over an ensemble of $S$ independent simulations with $N$ particles each in order to study the evolution of the one-point probability density $\Psi$ of finding a particle at a given location of phase space $(\mathbf{x},\mathbf{v})$ at time $t$. The proposed approach is extremely efficient from a computational point of view, with modest CPU and memory requirements, and it provides an alternative to traditional N-body simulations when the goal is to study the average properties of N-body systems, at the cost of abandoning the notion of well-defined trajectories for each individual particle. In one spatial dimension, our results, fully consistent with those previously reported in the literature for the standard deterministic formulation of the problem, highlight the differences between the evolution of the one-point probability density $\Psi(x,v,t)$ and the predictions of the collisionless Boltzmann (Vlasov-Poisson) equation, as well as the relatively subtle dependence on the actual finite number $N$ of particles in the system. We argue that understanding this dependence with $N$ may actually shed more light on the dynamics of real astrophysical systems than the limit $N\to\infty$.

We present the results of prompt optical follow-up of the electromagnetic counterpart of the gravitational-wave event GW170817 by the Transient Optical Robotic Observatory of the South Collaboration (TOROS). We detected highly significant dimming in the light curves of the counterpart (Delta g=0.17+-0.03 mag, Delta r=0.14+-0.02 mag, Delta i=0.10 +- 0.03 mag) over the course of only 80 minutes of observations obtained ~35 hr after the trigger with the T80-South telescope. A second epoch of observations, obtained ~59 hr after the event with the EABA 1.5m telescope, confirms the fast fading nature of the transient. The observed colors of the counterpart suggest that this event was a "blue kilonova" relatively free of lanthanides.

The Latin American Giant Observatory (LAGO) is an extended cosmic ray observatory composed by a network of water-Cherenkov detectors spanning over different sites located at significantly different altitudes (from sea level up to more than $5000$\,m a.s.l.) and latitudes across Latin America, covering a huge range of geomagnetic rigidity cut-offs and atmospheric absorption/reaction levels. This detection network is designed to measure the temporal evolution of the radiation flux at ground level with extreme detail. The LAGO project is mainly oriented to perform basic research in three branches: high energy phenomena, space weather and atmospheric radiation at ground level. LAGO is built and operated by the LAGO Collaboration, a non-centralized collaborative union of more than 30 institutions from ten countries. These are the contributions of the LAGO Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands

This is the third of a series of papers of low X-ray luminosity galaxy clusters. In this work we present the weak lensing analysis of eight clusters, based on observations obtained with the Gemini Multi-Object Spectrograph in the $g'$, $r'$ and $i'$ passbands. For this purpose, we have developed a pipeline for the lensing analysis of ground-based images and we have performed tests applied to simulated data. We have determined the masses of seven galaxy clusters, six of them measured for the first time. For the four clusters with availably spectroscopic data, we find a general agreement between the velocity dispersions obtained via weak lensing assuming a Singular Isothermal Sphere profile, and those obtained from the redshift distribution of member galaxies. The correlation between our weak lensing mass determinations and the X-ray luminosities are suitably fitted by other observations of the $M-L_{X}$ relation and models.

We present the deepest near-infrared (ZJKs) photometry yet obtained of the Sagittarius dwarf spheroidal (Sgr dSph), using VISTA to survey 11 square degrees centred on its core. We list locations and ZJKs-band magnitudes for over 2.9 million sources in the field. We discuss the isolation of the Sgr dSph from the foreground and Galactic Bulge populations, identify the Sgr dSph's horizontal branch in the near-infrared for the first time, and map the density of the galaxy's stars. We present isochrones for the Sgr dSph and Bulge populations. These are consistent with the previously-reported properties of the Sgr dSph core: namely that it is dominated by a population between [Fe/H] ~ -1 dex and solar, with a significant [alpha/Fe] versus [Fe/H] gradient. While strong contamination from the Galactic Bulge prevents accurate measurement of the (Galactic) north side of the Sgr dSph, the dwarf galaxy can be well-approximated by a roughly ovaloid projection of characteristic size 4 x 2 degrees, beyond which the projected stellar density is less than half that of the region surrounding the core. The galaxy's major axis is perpendicular to the Galactic Plane, as in previous studies. We find slight evidence to confirm a metallicity gradient in the Sgr dSph and use isochrones to fit a distance of 24.3 +/- 2.3 kpc. We were unable to fully constrain the metallicity distribution of the Sgr dSph due to the Bulge contamination and strong correlation of [alpha/Fe] with metallicity, however we find that metal-poor stars ([Fe/H] <~ --1) make up <~29 per cent of the Sgr dSph's upper-RGB population. The Bulge population is best fit by a younger population with [Fe/H] ~ 0 and [alpha/Fe] ~ 0 or slightly higher. We find no evidence for a split, peanut- or X-shaped Bulge population in this line of sight (l = 5.6 +/- ~1 deg, b = -14.1 +/- ~3 deg).

We introduce a new method for the identification of galaxy systems in redshift surveys based on the halo model. This method is a modified version of the K-means identification algorithm developed by Yang et al (2005). We have calibrated and tested our algorithms using mock catalogs generated using the Millennium simulations (Springel et al. 2005) and applied them to the NYU-DR7 galaxy catalog (based on the SDSS datasets). Using this local sample of groups and clusters of galaxies we have measured the effect of gravitational redshift produced by their host dark matter haloes. Our results shows radial velocity decrements consistent with general relativity predictions and previous measurements by Wojtak et al (2011) in clusters of galaxies.

The nature of the dominant component of galaxies and clusters remains unknown. While the astrophysics community supports the cold dark matter (CDM) paradigm as a clue factor in the current cosmological model, no direct CDM detections have been performed. Faber and Visser 2006 have suggested a simple method for measuring the dark matter equation of state that combines kinematic and gravitational lensing data to test the widely adopted assumption of pressureless dark matter. Following this formalism, we have measured the dark matter equation of state for first time using improved techniques. We have found that the value of the equation of state parameter is consistent with pressureless dark matter within the errors. Nevertheless, the measured value is lower than expected because typically the masses determined with lensing are larger than those obtained through kinematic methods. We have tested our techniques using simulations and we have also analyzed possible sources of error that could invalidate or mimic our results. In the light of this result, we can now suggest that the understanding of the nature of dark matter requires a complete general relativistic analysis.