Multi-messenger astronomy received a great boost following the discovery of kilonova AT2017gfo, the optical counterpart of the gravitational wave source GW170817 associated with the short gamma-ray burst GRB 170817A. AT2017gfo was the first kilonova that could be extensively monitored in time both photometrically and spectroscopically. Previously, only few candidates have been observed against the glare of short GRB afterglows. In this work, we aim to search the fingerprints of AT2017gfo-like kilonova emissions in the optical/NIR light curves of 39 short GRBs with known redshift. For the first time, our results allow us to study separately the range of luminosity of the blue and red components of AT2017gfo-like kilonovae in short GRBs. In particular, the red component is similar in luminosity to AT2017gfo, while the blue kilonova can be more than 10 times brighter. Finally, we find further evidence to support all the claimed kilonova detections and we exclude an AT2017gfo-like kilonova in GRBs 050509B and 061201.

Proper motion studies of stars in the centre of the Milky Way have been typically limited to the Arches and Quintuplet clusters and to the central parsec. Here, we present the first results of a large-scale proper motion study of stars within several tens of parsecs of Sagittarius A* based on our $0.2''$ angular resolution GALACTICNUCLEUS survey (epoch 2015) combined with NICMOS/HST data from the Paschen-$\alpha$ survey (epoch 2008). This study will be the first extensive proper motion study of the central $\sim 36' \times 16'$ of the Galaxy, which is not covered adequately by any of the existing astronomical surveys such as Gaia because of its extreme interstellar extinction ($A_{V} \gtrsim 30$ mag). Proper motions can help us to disentangle the different stellar populations along the line-of-sight and interpret their properties in combination with multi-wavelength photometry from GALACTICNUCLEUS and other sources. It also allows us to infer the dynamics and interrelationship between the different stellar components of the Galactic Centre (GC). In particular, we use proper motions to detect co-moving groups of stars which can trace low mass or partially dissolved young clusters in the GC that can hardly be discovered by any other means. Our pilot study in this work is on a field in the nuclear bulge associated by HII regions that show the presence of young stars. We detect the first group of co-moving stars coincident with an HII region. Using colour-magnitude diagrams, we infer that the co-moving stars are consistent with being the post-main sequence stars with ages of few Myrs. Simulations show that this group of stars is a real group that can indicate the existence of a dissolving or low to intermediate mass young cluster. A census of these undiscovered clusters will ultimately help us to constrain star formation at the GC in the past few ten Myrs.

The Radiative Transfer Model (RTM) and the retrieval algorithm, incorporated in the SCIATRAN 2.2 software package developed at the Institute of Remote Sensing/Institute of Enviromental Physics of Bremen University (Germany), allows to simulate, among other things, radiance/irradiance spectra in the 2400-24 000 {\AA} range. In this work we present applications of RTM to two case studies. In the first case the RTM was used to simulate direct solar irradiance spectra, with different water vapor amounts, for the study of the water vapor content in the atmosphere above Sierra Nevada Observatory. Simulated spectra were compared with those measured with a spectrometer operating in the 8000-10 000 {\AA} range. In the second case the RTM was used to generate telluric model spectra to subtract the atmospheric contribution and correct high-resolution stellar spectra from atmospheric water vapor and oxygen lines. The results of both studies are discussed.

The "Extreme starbursts in the local universe" workshop was held at the Insituto de Astrofisica de Andalucia in Granada, Spain on 21-25 June 2010. Bearing in mind the advent of a new generation of facilities such as JWST, Herschel, ALMA, eVLA and eMerlin, the aim of the workshop was to bring together observers and theorists to review the latest results. The purpose of the workshop was to address the following issues: what are the main modes of triggering extreme starbursts in the local Universe? How efficiently are stars formed in extreme starbursts? What are the star formation histories of local starburst galaxies? How well do the theoretical simulations model the observations? What can we learn about starbursts in the distant Universe through studies of their local counterparts? How important is the role of extreme starbursts in the hierarchical assembly of galaxies? How are extreme starbursts related to the triggering of AGN in the nuclei of galaxies? Overall, 41 talks and 4 posters with their corresponding 10 minutes short talks were presented during the workshop. In addition, the workshop was designed with emphasis on discussions, and therefore, there were 6 discussion sessions of up to one hour during the workshop. Here is presented a summary of the purposes of the workshop as well as a compilation of the abstracts corresponding to each of the presentations. The summary and conclusions of the workshop along with a description of the future prospects by Sylvain Veilleux can be found in the last section of this document. A photo of the assistants is included.

Over the past several years, numerous examples of X-ray cavities coincident with radio sources have been observed in so-called "cool core" clusters of galaxies. Motivated by these observations, we explore the evolution and the effect of cavities on a cooling intracluster medium (ICM) numerically, adding relevant physics step by step. In this paper we present a first set of hydrodynamical, high resolution (1024^3 effective grid elements), three-dimensional simulations, together with two-dimensional test cases. The simulations follow the evolution of radio cavities, modeled as bubbles filled by relativistic plasma, in the cluster atmosphere while the ICM is subject to cooling. We find that the bubble rise retards the development of a cooling flow by inducing motions in the ICM which repeatedly displace the material in the core. Even bubbles initially set significantly far from the cluster center affect the cooling flow, although much later than the beginning of the simulation. The effect is, however, modest: the cooling time is increased by at most only 25%. As expected, the overall evolution of pure hydrodynamic bubbles is at odds with observations, showing that some additional physics has to be considered in order to match the data.

We review the general properties of the intracluster medium (ICM) in clusters that host a cooling flow, and in particular the effects on the ICM of the injection of hot plasma by a powerful active galactic nucleus (AGN). It is observed that, in some cases, the hot plasma produces cavities in the ICM that finally detach and rise, perhaps buoyantly. The gas dynamics induced by the rising bubbles can help explain the absence of a cooled gas component in clusters with a cooling flow. This scenario is explored using numerical simulations.

Numerical hydro-N-body simulations are very important tools for making theoretical predictions for the formation of galaxy clusters. They show that the atmospheres of clusters of galaxies have quite complex angular and thermal structures. The full understanding of the physical processes behind these features can be only achieved by direct comparison of observations to hydro-N-body simulations. Although simple in principle, these comparisons are not always trivial. In fact, real data are convolved with the instrument response which may substantially influence the apparent properties of the studied features. To overcome this problem we build the software package X-MAS devoted to simulate X-ray observations of galaxy clusters obtained from hydro-N-body simulations. In this paper we present how this software package works and discuss its application to the simulation of Chandra ACIS-S3 observations. We compare some of the main physical properties of the input data to the ones derived from simulated observations after performing a standard imaging and spectral analysis. We show that, if the thermal structure of the cluster along a particular line of sight is quite complex, the projected spectroscopic temperature obtained from the observation is significantly lower than the emission-weighed value inferred directly from hydrodynamical simulation. This implies that much attention must be paid in the theoretical interpretation of observational temperatures.

A software package able to simulate imaging observations of galaxy clusters by the Chandra X-ray telescope is here presented. We start from high resolution N-body hydrodynamical simulations of galaxy clusters and assign to each gas particle a spectrum of emissivity, after assuming the MeKaL model. We then construct spatial images of the source differential flux which are used to create lists of incoming X-ray photons, preserving information on photon direction and energy. The photon lists are passed on to the Chandra simulator (MARX) to produce the final observation events. Background events are added to complete the simulation. Data analysis is currently in progress and simulated observations by other telescopes will become available in the future.

We examine how the statistics of the quadrupoles of (projected) cluster masses can discriminate between flat cold dark matter (CDM) universes with or without a cosmological constant term. Even in the era of high precision cosmology that cosmic microwave background experiments should open soon, it is important to devise self consistency tests of cosmogonic theories tuned at the matter radiation decoupling epoch using data from the non--linear evolved universe. We build cluster catalogs from two large volume simulations of a ``tilted'' CDM model and a $\Lambda$CDM model with cosmic density parameter $\Omega_m=0.35$ and cosmological constant contribution $\Omega_\Lambda=0.65$. From the projected mass distribution of the clusters we work out the quadrupoles Q and examine their dependence on cluster mass and the cosmological model. We find that TCDM clusters have systematically larger quadrupoles than their \lcdm counterpart. The effect is mass dependent: massive clusters ($M\gsim 10^{15}\Msunh$) have quadrupoles differing by more than 30% in the two models, while for $M\lsim 4\times10^{14}\Msunh$ the difference rapidly drops to $\sim 1$%. Performing a K-S test of the Q distributions, we estimate that using just the 15 most massive clusters in the simulation volume ($360\Mpch$ a side) we can discriminate between TCDM and \lcdm at a confidence level better than 99.9%. In the volume probed by exhisting observations, there are potentially several hundred clusters with masses above the threshold for which the differences in the quadrupoles become relevant. Should weak lensing data become available for this whole set, a quadrupole analysis may be expected to discriminate among different values of $\Lambda$.

We evaluate the dependence of the cluster correlation length r_c on the mean intercluster separation D_c, for three models with critical matter density, vanishing vacuum energy (Lambda = 0) and COBE normalized: a tilted CDM (tCDM) model (n=0.8) and two blue mixed models with two light massive neutrinos yielding Omega_h = 0.26 and 0.14 (MDM1 and MDM2, respectively). All models approach the observational value of sigma_8 (and, henceforth, the observed cluster abundance) and are consistent with the observed abundance of Damped Lyman_alpha systems. Mixed models have a motivation in recent results of neutrino physics; they also agree with the observed value of the ratio sigma_8/sigma_25, yielding the spectral slope parameter Gamma, and nicely fit LCRS reconstructed spectra. We use parallel AP3M simulations, performed in a wide box (side 360/h Mpc) and with high mass and distance resolution, enabling us to build artificial samples of clusters, whose total number and mass range allow to cover the same D_c interval inspected through APM and Abell cluster clustering data. We find that the tCDM model performs substantially better than n=1 critical density CDM models. Our main finding, however, is that mixed models provide a surprisingly good fit of cluster clustering data.

The passage from angular to spatial correlations, in the case of spatial clustering length depending on the average distance between nearby objects is studied. We show that, in a number of cases, the scaling law of angular correlation amplitudes, which holds for constant spatial clustering length, is still true also for a luminosity dependent spatial correlation. If the Limber equation is then naively used to obtain `the' spatial clustering length from the angular function amplitude, a quantity close to the average object separation is obtained. The case of cluster clustering is explicitly considered.

We study the cluster mass function in mixed dark matter (MDM) models, using two COBE normalized simulations with Omega_h = 0.26 and n=1.2, and Omega_h = 0.14 and n = 1.05, both with 2 massive nu's (MDM1 and MDM2, respectively). For the sake of comparison, we also simulate a CDM model with spectral index n=0.8 (TCDM), also COBE normalized. We argue that, in our non--hydro simulations, where CDM particles describe both actual CDM and baryons, the galaxy distribution traces CDM particles. Therefore, we use them to define clusters and their velocities to work out cluster masses. As CDM particles are more clustered than HDM and therefore have, in average, greater velocities, this leads to significant differences from PS predictions. Such predictions agree with simulations if both HDM and CDM are used to define clusters. Clusters defined through CDM in MDM models, instead, are less numerous than PS estimates, by a factor ~0.3, at the low mass end; the factor becomes \~0.6-0.8, depending on the mix, on intermediate mass scales (~4-5, h^-1 10^14 Msun) and almost vanishes on the high mass end. Therefore: (i) MDM models expected to overproduce clusters over intermediate scales are viable; (ii) the greater reduction factor at small scales agrees with the observational data dependence on the cluster mass M (which, however, may be partially due to sample incompleteness); (iii) the higher spectral normalization is felt at large scales, where MDM models produce more objects (hence, large clusters) than CDM. MDM1 even exceeds Donahue et al. (1998) findings, while MDM2 is consistent with them. (abridged)

We compare the distribution of internal velocity dispersions of galaxy clusters for an observational sample to those obtained from a set of N-body simulations of seven COBE-normalised cosmological scenarios: the standard CDM (SCDM) and a tilted (n=0.85) CDM (TCDM) model, a CHDM model with 25% of massive neutrinos, two low-density LCDM models with Omega_0=0.3 and 0.5, two open OCDM models with Omega_0=0.4 and 0.6. Simulated clusters are observed in projection so as to reproduce the main observational biases and are analysed by applying the same algorithm for interlopers removal and velocity dispersion estimate as for the reference observational sample. Velocity dispersions for individual clusters can be largely affected by observational biases in a model-dependent way: models in which clusters had less time to virialize show larger discrepancies between 3D and projected velocity dispersions. From the comparison with real clusters we find that both SCDM and TCDM largely overproduce clusters. The CHDM model marginally overproduces clusters and requires a somewhat larger sigma_8 than a purely CDM model in order to produce the same cluster abundance. The LCDM model with Omega_0=0.3 agrees with data, while the open model with Omega_0=0.4 and 0.6 underproduces and marginally overproduces clusters, respectively.