Dust plays an important role in the evolution of a galaxy, since it is one of the main ingredients for efficient star formation. Dust grains are also a sink/source of metals when they are created/destroyed, and, therefore, a self-consistent treatment is key in order to correctly model chemical evolution. In this work, we discuss the implementation of dust physics into our current multiphase model, which also follows the evolution of atomic, ionised and molecular gas. Our goal is to model the conversion rates among the different phases of the interstellar medium, including the creation, growth and destruction of dust, based on physical principles rather than phenomenological recipes inasmuch as possible. We first present the updated set of differential equations and then discuss the results. We calibrate our model against observations of the Milky Way Galaxy and compare its predictions with extant data. Our results are broadly consistent with the observed data for intermediate and high metallicities, but the models tend to produce more dust than observed in the low metallicity regime.

We present a 2-dimensional chemical evolution code applied to a Milky Way type galaxy, incorporating the role of spiral arms in shaping azimuthal abundance variations, and confront the predicted behaviour with recent observations taken with integral field units. To the usual radial distribution of mass, we add the surface density of the spiral wave and study its effect on star formation and elemental abundances. We compute five different models: one with azimuthal symmetry which depends only on radius, while the other four are subjected to the effect of a spiral density wave. At early times, the imprint of the spiral density wave is carried by both the stellar and star formation surface densities; conversely, the elemental abundance pattern is less affected. At later epochs, however, differences among the models are diluted, becoming almost indistinguishable given current observational uncertainties. At the present time, the largest differences appear in the star formation rate and/or in the outer disc (R$\ge$ 18\,kpc). The predicted azimuthal oxygen abundance patterns for $t \le 2$\,Gyr are in reasonable agreement with recent observations obtained with VLT/MUSE for NGC 6754

We study how the use of several stellar subtraction methods and line fitting approaches can affect the derivation of the main kinematic parameters (velocity and velocity dispersion fields) of the ionized gas component. The target of this work is the nearby galaxy NGC 2906, observed with the MUSE instrument at Very Large Telescope. A sample of twelve spectra is selected from the inner (nucleus) and outer (spiral arms) regions, characterized by different ionization mechanisms. We compare three different methods to subtract the stellar continuum (FIT3D, STARLIGHT and pPXF), combined with one of the following stellar libraries: MILES, STELIB and GRANADA+MILES. The choice of the stellar subtraction method is the most important ingredient affecting the derivation of the gas kinematics, followed by the choice of the stellar library and by the line fitting approach. In our data, typical uncertainties in the observed wavelength and width of the H\alpha and [NII] lines are of the order of <\delta\lambda>_rms \sim 0.1\AA\ and <\delta\sigma>_rms \sim 0.2\AA\ (\sim 5 and 10km/s, respectively). The results obtained from the [NII] line seem to be slightly more robust, as it is less affected by stellar absorption than H\alpha. All methods considered yield statistically consistent measurements once a mean systemic contribution \Delta\bar\lambda=\Delta\bar\sigma=0.2xDelta_{MUSE} is added in quadrature to the line fitting errors, where \Delta_{MUSE} = 1.1\AA\ \sim 50 km/s denotes the instrumental resolution of the MUSE spectra. Although the subtraction of the stellar continuum is critical in order to recover line fluxes, any method (including none) can be used in order to measure the gas kinematics, as long as an additional component of 0.2 x Delta_MUSE is added to the error budget.

We study the evolution of oxygen abundance radial gradients as a function of time for the Milky Way Galaxy obtained with our {\sc Mulchem} chemical evolution model. We review the recent data of abundances for different objects observed in our Galactic disc. We analyse with our models the role of the growth of the stellar disc, as well as the effect of infall rate and star formation prescriptions, or the pre-enrichment of the infall gas, on the time evolution of the oxygen abundance radial distribution. We compute the radial gradient of abundances within the {\sl disk}, and its corresponding evolution, taking into account the disk growth along time. We compare our predictions with the data compilation, showing a good agreement. Our models predict a very smooth evolution when the radial gradient is measured within the optical disc with a slight flattening of the gradient from $\sim -0.057$\,dex\,kpc$^{-1}$ at $z=4$ until values around $\sim -0.015$\,dex\,kpc$^{-1}$ at $z=1$ and basically the same gradient until the present, with small differences between models. Moreover, some models show a steepening at the last times, from $z=1$ until $z=0$ in agreement with data which give a variation of the gradient in a range from $-0.02$ to $-0.04$\,de\,kpc$^{-1}$ from $t=10$\,Gyr until now. The gradient measured as a function of the normalized radius $R/R_{\rm eff}$ is in good agreement with findings by CALIFA and MUSE, and its evolution with redshift falls within the error bars of cosmological simulations.

The Advanced Telescope for High Energy Astrophysics (Athena) is the X-ray observatory large mission selected by the European Space Agency (ESA), within its Cosmic Vision 2015-2025 programme, to address the "Hot and Energetic Universe" scientific theme (Nandra et al. 2013), and it is provisionally due for launch in the early 2030s. The Square Kilometer Array (SKA) is the next generation radio observatory and consists of two telescopes, one comprised of dishes operating at mid frequencies (SKA1-MID) and located in South Africa, and the other comprised of Log-Periodic antennas operating at low radio frequencies (SKA1-LOW), which will be located in Australia (Braun et al. 2017). The scientific commissioning of the radio telescope is planned to begin in 2021-2022. The SKA-Athena Synergy Team (SAST) has been tasked to single out the potential scientific synergies between Athena and SKA. The astrophysical community was involved in this exercise primarily through a dedicated SKA-Athena Synergy Workshop, which took place on April 24-25, 2017 at SKAO, Jodrell Bank, Manchester. The final result of the synergy exercise, this White Paper, describes in detail a number of scientific opportunities that will be opened up by the combination of Athena and SKA, these include: 1. the Cosmic Dawn; 2. the Evolution of black holes and galaxies; 3. Active galaxy feedback in galaxy clusters; 4. Non-thermal phenomena in galaxy clusters; 5. Detecting the cosmic web; 6. Black-hole accretion physics and astrophysical transients; 7. Galactic astronomy: stars, planets, pulsars and supernovae.

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$.

In this work we present a set of synthetic observations that mimic the properties of the Integral Field Spectroscopy (IFS) survey CALIFA, generated using radiative transfer techniques applied to hydrodynamical simulations of galaxies in a cosmological context. The simulated spatially-resolved spectra include stellar and nebular emission, kinematic broadening of the lines, and dust extinction and scattering. The results of the radiative transfer simulations have been post-processed to reproduce the main properties of the CALIFA V500 and V1200 observational setups. The data has been further formatted to mimic the CALIFA survey in terms of field of view size, spectral range and sampling. We have included the effect of the spatial and spectral Point Spread Functions affecting CALIFA observations, and added detector noise after characterizing it on a sample of 20 galaxies. The simulated datacubes are suited to be analyzed by the same algorithms used on real IFS data. In order to provide a benchmark to compare the results obtained applying IFS observational techniques to our synthetic datacubes, and test the calibration and accuracy of the analysis tools, we have computed the spatially-resolved properties of the simulations. Hence, we provide maps derived directly from the hydrodynamical snapshots or the noiseless spectra, in a way that is consistent with the values recovered by the observational analysis algorithms. Both the synthetic observations and the product datacubes are public and can be found in the collaboration website http://astro.ft.uam.es/selgifs/data_challenge/.

Spiral arms are the most singular features in disc galaxies. These structures can exhibit different patterns, namely grand design and flocculent arms, with easily distinguishable characteristics. However, their origin and the mechanisms shaping them are unclear. The overall role of spirals in the chemical evolution of disc galaxies is another unsolved question. In particular, it has not been fully explored if the \hii\,regions of spiral arms present different properties from those located in the interarm regions. Here we analyse the radial oxygen abundance gradient of the arm and interarm star forming regions of 63 face-on spiral galaxies using CALIFA Integral Field Spectroscopy data. We focus the analysis on three characteristic parameters of the profile: slope, zero-point, and scatter. The sample is morphologically separated into flocculent versus grand design spirals and barred versus unbarred galaxies. We find subtle but statistically significant differences between the arm and interarm distributions for flocculent galaxies, suggesting that the mechanisms generating the spiral structure in these galaxies may be different to those producing grand design systems, for which no significant differences are found. We also find small differences in barred galaxies, not observed in unbarred systems, hinting that bars may affect the chemical distribution of these galaxies but not strongly enough as to be reflected in the overall abundance distribution. In light of these results, we propose bars and flocculent structure as two distinct mechanisms inducing differences in the abundance distribution between arm and interarm star forming regions.

In our grid of multiphase chemical evolution models (Moll\'a & D\'iaz, 2005), star formation in the disk occurs in two steps: first, molecular gas forms, and then stars are created by cloud-cloud collisions or interactions of massive stars with the surrounding molecular clouds. The formation of both molecular clouds and stars are treated through the use of free parameters we refer to as efficiencies. In this work we modify the formation of molecular clouds based on several new prescriptions existing in the literature, and we compare the results obtained for a chemical evolution model of the Milky Way Galaxy regarding the evolution of the Solar region, the radial structure of the Galactic disk, and the ratio between the diffuse and molecular components, HI/H$_2$. Our results show that the six prescriptions we have tested reproduce fairly consistent most of the observed trends, differing mostly in their predictions for the (poorly-constrained) outskirts of the Milky Way and the evolution in time of its radial structure. Among them, the model proposed by Ascasibar et al. (2017), where the conversion of diffuse gas into molecular clouds depends on the local stellar and gas densities as well as on the gas metallicity, seems to provide the best overall match to the observed data.

Galaxies in the local Universe are known to follow bimodal distributions in the global stellar populations properties. We analyze the distribution of the local average stellar-population ages of 654,053 sub-galactic regions resolved on ~1-kpc scales in a volume-corrected sample of 394 galaxies, drawn from the CALIFA-DR3 integral-field-spectroscopy survey and complemented by SDSS imaging. We find a bimodal local-age distribution, with an old and a young peak primarily due to regions in early-type galaxies and star-forming regions of spirals, respectively. Within spiral galaxies, the older ages of bulges and inter-arm regions relative to spiral arms support an internal age bimodality. Although regions of higher stellar-mass surface-density, mu*, are typically older, mu* alone does not determine the stellar population age and a bimodal distribution is found at any fixed mu*. We identify an "old ridge" of regions of age ~9 Gyr, independent of mu*, and a "young sequence" of regions with age increasing with mu* from 1-1.5 Gyr to 4-5 Gyr. We interpret the former as regions containing only old stars, and the latter as regions where the relative contamination of old stellar populations by young stars decreases as mu* increases. The reason why this bimodal age distribution is not inconsistent with the unimodal shape of the cosmic-averaged star-formation history is that i) the dominating contribution by young stars biases the age low with respect to the average epoch of star formation, and ii) the use of a single average age per region is unable to represent the full time-extent of the star-formation history of "young-sequence" regions.

This paper describes the Bayesian Technique for Multi-image Analysis (BaTMAn), a novel image-segmentation technique based on Bayesian statistics that characterizes any astronomical dataset containing spatial information and performs a tessellation based on the measurements and errors provided as input. The algorithm iteratively merges spatial elements as long as they are statistically consistent with carrying the same information (i.e. identical signal within the errors). We illustrate its operation and performance with a set of test cases including both synthetic and real Integral-Field Spectroscopic data. The output segmentations adapt to the underlying spatial structure, regardless of its morphology and/or the statistical properties of the noise. The quality of the recovered signal represents an improvement with respect to the input, especially in regions with low signal-to-noise ratio. However, the algorithm may be sensitive to small-scale random fluctuations, and its performance in presence of spatial gradients is limited. Due to these effects, errors may be underestimated by as much as a factor of two. Our analysis reveals that the algorithm prioritizes conservation of all the statistically-significant information over noise reduction, and that the precise choice of the input data has a crucial impact on the results. Hence, the philosophy of BaTMAn is not to be used as a `black box' to improve the signal-to-noise ratio, but as a new approach to characterize spatially-resolved data prior to its analysis. The source code is publicly available at http://astro.ft.uam.es/SELGIFS/BaTMAn .

The aim of this work is to present our new series of chemical evolution models computed for spiral and low mass galaxies of different total masses and star formation efficiencies. We analyze the results of models, in particular the evolution of the radial gradient of oxygen abundance. Furthermore, we study the role of the infall rate and of the star formation history on the variations of this radial gradient. The relations between the O/H radial gradient and other spiral galaxies characteristics as the size or the stellar mass are also shown. We find that the radial gradient is mainly a scale effect which basically does not change with the redshift (or time) if it is measured within the optical radius. Moreover, when it is measured as a function of a normalized radius, show a similar value for all galaxies masses, showing a correlation with a dispersion around an average value which is due to the differences star formation efficiencies, in agreement with the idea of an universal O/H radial gradient

The velocity function is a fundamental observable statistic of the galaxy population, similarly impor- tant as the luminosity function, but much more difficult to measure. In this work we present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km/s for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy population (60 < v_circ < 170 km/s), we use the HIPASS velocity function. For the massive galaxy population (170 < v_circ < 320 km/s), we use stellar circular velocities from the Calar Alto Legacy Integral Field Area Survey (CALIFA). In earlier work we obtained the measurements of circular velocity at the 80% light radius for 226 galaxies and demonstrated that the CALIFA sample can produce volume- corrected galaxy distribution functions. The CALIFA velocity function includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies and has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both velocity functions can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed velocity function is compared to velocity functions derived from cosmological simulations of the z = 0 galaxy population. We find that dark matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic simulations fare better, but still do not fully reproduce observations.

This paper describes the Third Public Data Release (DR3) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. Science-grade quality data for 667 galaxies are made public, including the 200 galaxies of the Second Public Data Release (DR2). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto Observatory. Three different spectral setups are available, i) a low-resolution V500 setup covering the wavelength range 3749-7500 AA (4240-7140 AA unvignetted) with a spectral resolution of 6.0 AA (FWHM), for 646 galaxies, ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 AA (3650-4620 AA unvignetted) with a spectral resolution of 2.3 AA (FWHM), for 484 galaxies, and iii) the combination of the cubes from both setups (called COMBO), with a spectral resolution of 6.0 AA and a wavelength range between 3700-7500 AA (3700-7140 AA unvignetted), for 446 galaxies. The Main Sample, selected and observed according to the CALIFA survey strategy covers a redshift range between 0.005 and 0.03, spans the color-magnitude diagram and probes a wide range of stellar mass, ionization conditions, and morphological types. The Extension Sample covers several types of galaxies that are rare in the overall galaxy population and therefore not numerous or absent in the CALIFA Main Sample. All the cubes in the data release were processed using the latest pipeline, which includes improved versions of the calibration frames and an even further improved im- age reconstruction quality. In total, the third data release contains 1576 datacubes, including ~1.5 million independent spectra. It is available at http://califa.caha.es/DR3.

The aim of this paper is to characterize the radial structure of the star formation rate (SFR) in galaxies in the nearby Universe as represented by the CALIFA survey. The sample under study contains 416 galaxies observed with IFS, covering a wide range of Hubble types and stellar masses. Spectral synthesis techniques are applied to obtain radial profiles of the intensity of the star formation rate in the recent past, and the local sSFR. To emphasize the behavior of these properties for galaxies that are on and off the main sequence of star formation (MSSF) we stack the individual radial profiles in bins of galaxy morphology and stellar masses. Our main results are: a) The intensity of SFR shows declining profiles that exhibit very little differences between spirals. The dispersion between the profiles is significantly smaller in late type spirals. This confirms that the MSSF is a sequence of galaxies with nearly constant intensity of SFR b) sSFR values scale with Hubble type and increase radially outwards, with a steeper slope in the inner 1 HLR. This behavior suggests that galaxies are quenched inside-out, and that this process is faster in the central, bulge-dominated part than in the disks. c) As a whole, and at all radii, E and S0 are off the MSSF. d) Applying the volume-corrections for the CALIFA sample, we obtain a density of star formation in the local Universe of 0.0105 Msun/yr/Mpc^{-3}. Most of the star formation is occurring in the disks of spirals. e) The volume averaged birthrate parameter, b'=0.39, suggests that the present day Universe is forming stars at 1/3 of its past average rate. E, S0, and the bulge of early type spirals contribute little to the recent SFR of the Universe, which is dominated by the disks of later spirals. f) There is a tight relation between the intensity of the SFR and stellar mass, defining a local MSSF relation with a logarithmic slope of 0.8.

In this work, we investigate the contribution of dust scattering to the diffuse H-alpha emission observed in nearby galaxies. As initial conditions for the spatial distribution of HII regions, gas, and dust, we take three Milky Way-like galaxies from state-of-the-art cosmological hydrodynamical simulations that implement different prescriptions for star formation, feedback, and chemical enrichment. Radiative transfer has been solved a posteriori, using the publicly-available Monte Carlo code Sunrise to take into account dust absorption and scattering of the H-alpha photons, originating exclusively from the HII regions. No contribution from recombinations in the diffuse ionized gas (DIG) component is explicitly or implicitly included in our model. Our main result is that the flux arising from scattered light is of the order of 1-2 per cent of the H-alpha flux coming directly from the HII regions. Building upon previous studies, we conclude that the DIG contributes lass than 50 per cent of the total H-alpha emission.

The "main sequence of galaxies" $-$ defined in terms of the total star formation rate $\psi$ vs. the total stellar mass $M_*$ $-$ is a well-studied tight relation that has been observed at several wavelengths and at different redshifts. All earlier studies have derived this relation from integrated properties of galaxies. We recover the same relation from an analysis of spatially-resolved properties, with integral field spectroscopic (IFS) observations of 306 galaxies from the CALIFA survey. We consider the SFR surface density in units of log(M$_{\odot}$ yr$^{-1}$ Kpc$^{-2}$) and the stellar mass surface density in units of log(M$_{\odot}$ Kpc$^{-2}$) in individual spaxels which probe spatial scales of 0.5-1.5 Kpc. This local relation exhibits a high degree of correlation with small scatter ($\sigma = 0.23$ dex), irrespective of the dominant ionisation source of the host galaxy or its integrated stellar mass. We highlight: $(i)$ the integrated star formation main sequence formed by galaxies whose dominant ionisation process is related to star formation, for which we find a slope of 0.81 $\pm 0.02$; (ii) the spatially-resolved relation obtained with the spaxel analysis, we find a slope of 0.72 $\pm 0.04$; (iii) for the integrated main sequence we identified also a sequence formed by galaxies that are dominated by an old stellar population, which we have called the retired galaxies sequence.

We present Pipe3D, an analysis pipeline based on the FIT3D fitting tool, devel- oped to explore the properties of the stellar populations and ionized gas of Integral Field Spectroscopy data. Pipe3D was created to provide with coherent, simple to distribute, and comparable dataproducts, independently of the origin of the data, focused on the data of the most recent IFU surveys (e.g., CALIFA, MaNGA, and SAMI), and the last generation IFS instruments (e.g., MUSE). Along this article we describe the different steps involved in the analysis of the data, illustrating them by showing the dataproducts derived for NGC 2916, observed by CALIFA and P-MaNGA. As a practical use of the pipeline we present the complete set of dataproducts derived for the 200 datacubes that comprises the V500 setup of the CALIFA Data Release 2 (DR2), making them freely available through the network (ftp://ftp.caha.es/CALIFA/dataproducts/DR2/Pipe3D). Finally, we explore the hypothesis that the properties of the stellar populations and ionized gas of galaxies at the effective radius are representative of the overall average ones, finding that this is indeed the case.

We measured the gas abundance profiles in a sample of 122 face-on spiral galaxies observed by the CALIFA survey and included all spaxels whose line emission was consistent with star formation. This type of analysis allowed us to improve the statistics with respect to previous studies, and to properly estimate the oxygen distribution across the entire disc to a distance of up to 3-4 disc effective radii (r$_e$). We confirm the results obtained from classical HII region analysis. In addition to the general negative gradient, an outer flattening can be observed in the oxygen abundance radial profile. An inner drop is also found in some cases. There is a common abundance gradient between 0.5 and 2.0 r$_e$ of $\alpha_{O/H} = -\,0.075\,\rm{dex}/r_e$ with a scatter of $\sigma = 0.016\,\rm{dex}/r_e$ when normalising the distances to the disc effective radius. By performing a set of Kolmogorov-Smirnov tests, we determined that this slope is independent of other galaxy properties, such as morphology, absolute magnitude, and the presence or absence of bars. In particular, barred galaxies do not seem to display shallower gradients, as predicted by numerical simulations. Interestingly, we find that most of the galaxies in the sample with reliable oxygen abundance values beyond $\sim 2$ effective radii (57 galaxies) present a flattening of the abundance gradient in these outer regions. This flattening is not associated with any morphological feature, which suggests that it is a common property of disc galaxies. Finally, we detect a drop or truncation of the abundance in the inner regions of 27 galaxies in the sample; this is only visible for the most massive galaxies.

We present an improved version of FIT3D, a fitting tool for the analysis of the spectroscopic properties of the stellar populations and the ionized gas derived from moderate resolution spectra of galaxies. FIT3D is a tool developed to analyze Integral Field Spectroscopy data and it is the basis of Pipe3D, a pipeline already used in the analysis of datasets like CALIFA, MaNGA, and SAMI. We describe the philosophy behind the fitting procedure, and in detail each of the different steps in the analysis. We present an extensive set of simulations in order to estimate the precision and accuracy of the derived parameters for the stellar populations. In summary, we find that using different stellar population templates we reproduce the mean properties of the stellar population (age, metallicity, and dust attenuation) within ~0.1 dex. A similar approach is adopted for the ionized gas, where a set of simulated emission- line systems was created. Finally, we compare the results of the analysis using FIT3D with those provided by other widely used packages for the analysis of the stellar population (Starlight, Steckmap, and analysis based on stellar indices) using real high S/N data. In general we find that the parameters for the stellar populations derived by FIT3D are fully compatible with those derived using these other tools.

Models of galaxy formation predict that gas accretion from the cosmic web is a primary driver of star formation over cosmic history. Except in very dense environments where galaxy mergers are also important, model galaxies feed from cold streams of gas from the web that penetrate their dark matter haloes. Although these predictions are unambiguous, the observational support has been indirect so far. Here we report spectroscopic evidence for this process in extremely metal-poor galaxies (XMPs) of the local Universe, taking the form of localized starbursts associated with gas having low metallicity. Detailed abundance analyses based on Gran Telescopio Canarias (GTC) optical spectra of ten XMPs show that the galaxy hosts have metallicities around 60 % solar on average, while the large star-forming regions that dominate their integrated light have low metallicities of some 6 % solar. Because gas mixes azimuthally in a rotation timescale (a few hundred Myr), the observed metallicity inhomogeneities are only possible if the metal-poor gas fell onto the disk recently. We analyze several possibilities for the origin of the metal-poor gas, favoring the metal-poor gas infall predicted by numerical models. If this interpretation is correct, XMPs trace the cosmic web gas in their surroundings, making them probes to examine its properties.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

This work investigates the main mechanism(s) that regulate the specific star formation rate (SSFR) in nearby galaxies, cross-correlating two proxies of this quantity -- the equivalent width of the \Ha\ line and the $(u-r)$ colour -- with other physical properties (mass, metallicity, environment, morphology, and the presence of close companions) in a sample of $\sim82500$ galaxies extracted from the Sloan Digital Sky Survey (SDSS). The existence of a relatively tight `ageing sequence' in the colour-equivalent width plane favours a scenario where the secular conversion of gas into stars (i.e. `nature') is the main physical driver of the instantaneous SSFR and the gradual transition from a `chemically primitive' (metal-poor and intensely star-forming) state to a `chemically evolved' (metal-rich and passively evolving) system. Nevertheless, environmental factors (i.e. `nurture') are also important. In the field, galaxies may be temporarily affected by discrete `quenching' and `rejuvenation' episodes, but such events show little statistical significance in a probabilistic sense, and we find no evidence that galaxy interactions are, on average, a dominant driver of star formation. Although visually classified mergers tend to display systematically higher EW(H$\alpha$) and bluer $(u-r)$ colours for a given luminosity, most galaxies with high SSFR have uncertain morphologies, which could be due to either internal or external processes. Field galaxies of early and late morphological types are consistent with the gradual `ageing' scenario, with no obvious signatures of a sudden decrease in their SSFR. In contrast, star formation is significantly reduced and sometimes completely quenched on a short time scale in dense environments, where many objects are found on a `quenched sequence' in the colour-equivalent width plane.

This work studies the relation between gas-phase oxygen abundance and stellar-to-gas fraction in nearby galaxies. We first derive the theoretical prediction, and argue that this relation is fundamental, in the sense that it must be verified regardless of the details of the gas accretion and star formation histories. Moreover, it should hold on "local" scales, i.e. in regions of the order of 1 kpc. These predictions are then compared with a set of spectroscopic observations, including both integrated and resolved data. Although the results depend somewhat on the adopted metallicity calibration, observed galaxies are consistent with the predicted relation, imposing tight constraints on the mass-loading factor of (enriched) galactic winds. The proposed parametrization of the star fraction-metallicity relation is able to describe the observed dependence of the oxygen abundance on gas mass at fixed stellar mass. However, the "local" mass-metallicity relation also depends on the relation between stellar and gas surface densities.

This paper describes the Second Public Data Release (DR2) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. The data for 200 objects are made public, including the 100 galaxies of the First Public Data Release (DR1). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5 m telescope at the Calar Alto observatory. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the wavelength range 3745-7500 \AA\ with a spectral resolution of 6.0 \AA\ (FWHM), and (ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 \AA\ with a spectral resolution of 2.3 \AA\ (FWHM). The sample covers a redshift range between 0.005 and 0.03, with a wide range of properties in the Color-Magnitude diagram, stellar mass, ionization conditions, and morphological types. All released cubes were reduced with the latest pipeline, including improved spectrophotometric calibration, spatial registration and spatial resolution. The spectrophotometric calibration is better than 6% and the median spatial resolution is 2.5". Altogether the second data release contains over 1.5 million spectra. It is available at http://califa.caha.es/DR2.

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

We investigate the production of electrons and positrons in the Milky Way within the context of dark matter annihilation. Upper limits on the relevant cross-section are obtained by combining observational data at different wavelengths (from Haslam, WMAP, and Fermi all-sky intensity maps) with recent measurements of the electron and positron spectra in the solar neighbourhood by PAMELA, Fermi, and HESS. We consider synchrotron emission in the radio and microwave bands, as well as inverse Compton scattering and final-state radiation at gamma-ray energies. According to our results, the dark matter annihilation cross-section into electron-positron pairs should not be higher than the canonical value for a thermal relic if the mass of the dark matter candidate is smaller than a few GeV. In addition, we also derive a stringent upper limit on the inner logarithmic slope, alpha, of the density profile of the Milky Way dark matter halo (alpha < 1 if m_dm < 5 GeV, alpha < 1.3 if m_dm < 100 GeV and alpha < 1.5 if m_dm < 2 TeV) assuming that cross-section = 3 x 10^(-26) cm^3 s(-1). A logarithmic slope steeper than alpha about 1.5 is hardly compatible with a thermal relic lighter than about 1 TeV, regardless of the dominant annihilation channel.

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

This work studies a simplified model of the gravitational instability of an initially homogeneous infinite medium, represented by $\TT^d$, based on the approximation that the mean fluid velocity is always proportional to the local acceleration. It is shown that, mathematically, this assumption leads to the restricted Patlak-Keller-Segel model considered by J\"ager and Luckhaus or, equivalently, the Smoluchowski equation describing the motion of self-gravitating Brownian particles, coupled to the modified Newtonian potential that is appropriate for an infinite mass distribution. We discuss some of the fundamental properties of a non-local generalization of this model where the effective pressure force is given by a fractional Laplacian with $0<\alpha<2$, and illustrate them by means of numerical simulations. Local well-posedness in Sobolev spaces is proven, and we show the smoothing effect of our equation, as well as a \emph{Beale-Kato-Majda}-type criterion in terms of $\rhomax$. It is also shown that the problem is ill-posed in Sobolev spaces when it is considered backward in time. Finally, we prove that, in the critical case (one conservative and one dissipative derivative), $\rhomax(t)$ is uniformly bounded in terms of the initial data for sufficiently large pressure forces.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

Extremely metal-poor (XMP) galaxies are chemically, and possibly dynamically, primordial objects in the local Universe. Our objective is to characterize the HI content of the XMP galaxies as a class, using as a reference the list of 140 known local XMPs compiled by Morales-Luis et al. (2011). We have observed 29 XMPs, which had not been observed before at 21 cm, using the Effelsberg radio telescope. This information was complemented with HI data published in literature for a further 53 XMPs. In addition, optical data from the literature provided morphologies, stellar masses, star-formation rates and metallicities. Effelsberg HI integrated flux densities are between 1 and 15 Jy km/s, while line widths are between 20 and 120 km/s. HI integrated flux densities and line widths from literature are in the range 0.1 - 200 Jy km/s and 15 - 150 km/s, respectively. Of the 10 new Effelsberg detections, two sources show an asymmetric double-horn profile, while the remaining sources show either asymmetric (7 sources) or symmetric (1 source) single-peak 21 cm line profiles. An asymmetry in the HI line profile is systematically accompanied by an asymmetry in the optical morphology. Typically, the g-band stellar mass-to-light ratios are ~0.1, whereas the HI gas mass-to-light ratios may be up to 2 orders of magnitude larger. Moreover, HI gas-to-stellar mass ratios fall typically between 10 and 20, denoting that XMPs are extremely gas-rich. We find an anti-correlation between the HI gas mass-to-light ratio and the luminosity, whereby fainter XMPs are more gas-rich than brighter XMPs, suggesting that brighter sources have converted a larger fraction of their HI gas into stars. The dynamical masses inferred from the HI line widths imply that the stellar mass does not exceed 5% of the dynamical mass, while the \ion{H}{i} mass constitutes between 20 and 60% of the dynamical mass. (abridged)

While various codes exist to systematically and robustly find haloes and subhaloes in cosmological simulations (Knebe et al., 2011, Onions et al., 2012), this is the first work to introduce and rigorously test codes that find tidal debris (streams and other unbound substructure) in fully cosmological simulations of structure formation. We use one tracking and three non-tracking codes to identify substructure (bound and unbound) in a Milky Way type simulation from the Aquarius suite (Springel et al., 2008) and post-process their output with a common pipeline to determine the properties of these substructures in a uniform way. By using output from a fully cosmological simulation, we also take a step beyond previous studies of tidal debris that have used simple toy models. We find that both tracking and non-tracking codes agree well on the identification of subhaloes and more importantly, the {\em unbound tidal features} associated with them. The distributions of basic properties of the total substructure distribution (mass, velocity dispersion, position) are recovered with a scatter of $\sim20%$. Using the tracking code as our reference, we show that the non-tracking codes identify complex tidal debris with purities of $\sim40%$. Analysing the results of the substructure finders, we find that the general distribution of {\em substructures} differ significantly from the distribution of bound {\em subhaloes}. Most importantly, both bound and unbound {\em substructures} together constitute $\sim18%$ of the host halo mass, which is a factor of $\sim2$ higher than the fraction in self-bound {\em subhaloes}. However, this result is restricted by the remaining challenge to cleanly define when an unbound structure has become part of the host halo. Nevertheless, the more general substructure distribution provides a more complete picture of a halo's accretion history.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles.We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important.

Recent measurements of the temperature of the sky in the radio band, combined with literature data, have convincingly shown the existence of a cosmic radio background with an amplitude of $\sim 1$ K at 1 GHz and a spectral energy distribution that is well described by a power law with index $\alpha \simeq -0.6$. The origin of this signal remains elusive, and it has been speculated that it could be dominated by the contribution of star-forming galaxies at high redshift \change{if the far infrared-radio correlation $q(z)$ evolved} in time. \change{We fit observational data from several different experiments by the relation $q(z) \simeq q_0 - \beta \log(1+z)$ with $q_0 = 2.783 \pm 0.024$ and $\beta = 0.705 \pm 0.081$ and estimate the total radio emission of the whole galaxy population at any given redshift from the cosmic star formation rate density at that redshift. It is found that} star-forming galaxies can only account for $\sim$13 percent of the observed intensity of the cosmic radio background.

Using analytical models and cosmological N-body simulations, we study the free-free radio emission from ionized gas in clusters and groups of galaxies. The results obtained with the simulations are compared with analytical predictions based on the mass function and scaling relations. Earlier works based on analytical models have shown that the average free-free signal from small haloes (galaxies) during and after the reionization time could be detected with future experiments as a distortion of the CMB spectrum at low frequencies ($\nu <$ 5 GHz). We focus on the period after the reionization time (from redshift $z=0$ up to $z=7$) and on haloes that are more massive than in previous works (groups and clusters). We show how the average signal from haloes with $M > 10^{13} h^{-1} M_{\odot}$ is less than 10% the signal from the more abundant and colder smaller mass haloes. However, the individual signal from the massive haloes could be detected with future experiments opening the door for a new window to study the intracluster medium.

Due to its proximity to the Earth and its nearly edge-on geometry, the Orion Bar provides an excellent testbed for detailed models of the structure of HII regions and the surrounding photon-dominated regions. In the present study, a self-consistent model of the structure of the Orion Nebula in the vicinity of the Bar is built under the assumption of approximate ionization, thermal, and hydrostatic equilibrium. It is found that a fairly simple geometry is able to describe the surface brightness profiles of the emission lines tracing the ionized HII region with a remarkable accuracy, independent of the prescription adopted to set the magnetic field or the population of cosmic rays. Although we consider different scenarios for these non-thermal components, none of the models is able to provide a fully satisfactory match to the observational data for the atomic layer, and the predicted column densities of several molecular species are always well above the measured abundances. Contrary to previous studies, we conclude that a more elaborate model is required in order to match all the available data.

We identify a spectroscopic sequence of galaxies, analogous to the Hubble sequence of morphological types, based on the Automatic Spectroscopic K-means (ASK) classification. Considering galaxy spectra as multidimensional vectors, the majority of the spectral classes are distributed along a well defined curve going from the earliest to the latest types, suggesting that the optical spectra of normal galaxies can be described in terms of a single affine parameter. Optically-bright active galaxies, however, appear as an independent, roughly orthogonal branch that intersects the main sequence exactly at the transition between early and late types.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

The rotation curves of spiral galaxies are one of the basic predictions of the cold dark matter paradigm, and their shape in the innermost regions has been hotly debated over the last decades. The present work shows that dark matter annihilation into electron-positron pairs may affect the observed rotation curve by a significant amount. We adopt a model-independent approach, where all the electrons and positrons are injected with the same initial energy E_0 ~ m_dm*c^2 in the range from 1 MeV to 1 TeV and the injection rate is constrained by INTEGRAL, Fermi, and HESS data. The pressure of the relativistic electron-positron gas is determined by solving the diffusion-loss equation, considering inverse Compton scattering, synchrotron radiation, Coulomb collisions, bremsstrahlung, and ionization. For values of the gas density and magnetic field that are representative of the Milky Way, it is estimated that pressure gradients are strong enough to balance gravity in the central parts if E_0 < 1 GeV. The exact value depends somewhat on the astrophysical parameters, and it changes dramatically with the slope of the dark matter density profile. For very steep slopes, as those expected from adiabatic contraction, the rotation curves of spiral galaxies would be affected on ~kpc scales for most values of E_0. By comparing the predicted rotation curves with observations of dwarf and low surface brightness galaxies, we show that the pressure from dark matter annihilation may improve the agreement between theory and observations in some cases, but it also imposes severe constraints on the model parameters (most notably, the inner slope of halo density profile, as well as the mass and the annihilation cross-section of dark matter particles into electron-positron pairs).

The Field Estimator for Arbitrary Spaces (FiEstAS) computes the continuous probability density field underlying a given discrete data sample in multiple, non-commensurate dimensions. The algorithm works by constructing a metric-independent tessellation of the data space based on a recursive binary splitting. Individual, data-driven bandwidths are assigned to each point, scaled so that a constant "mass" M0 is enclosed. Kernel density estimation may then be performed for different kernel shapes, and a combination of balloon and sample point estimators is proposed as a compromise between resolution and variance. A bias correction is evaluated for the particular (yet common) case where the density is computed exactly at the locations of the data points rather than at an uncorrelated set of locations. By default, the algorithm combines a top-hat kernel with M0=2.0 with the balloon estimator and applies the corresponding bias correction. These settings are shown to yield reasonable results for a simple test case, a two-dimensional ring, that illustrates the performance for oblique distributions, as well as for a six-dimensional Hernquist sphere, a fairly realistic model of the dynamical structure of stellar bulges in galaxies and dark matter haloes in cosmological N-body simulations. Results for different parameter settings are discussed in order to provide a guideline to select an optimal configuration in other cases. Source code is available upon request.

We report sub-arcsecond scale radio continuum observations of a field of 35 by 22 arcmin centred in NGC3351 obtained with the Multi-Element Radio Linked Interferometer Network (MERLIN). We found 23 radio sources in this field, 6 of which are projected within the D$_{25}$ isophote of the galaxy, and 3 are located inside the central 100 arcsec in radius. Two of these three are significantly extended, while the third one is relatively compact. This one is the only source with a previously detected counterpart at other wavelengths and could constitute the radio counterpart of a young supernova remnant. The other two are probably related to jets from a background AGN. We are not able to detect individual supernovae or SNRs in the central region ($r<600$ pc) of the galaxy. This could imply that the ionising populations of the circumnuclear star-forming regions (CNSFRs) are too young (less than a few Myr) to host supernovae. Also the diffusion length of the relativistic electrons in the ISM associated with the SN from previous events of star formation seems to be larger than our maximum resolution of 50 pc in radius. Detecting the thermal bremsstrahlung emission from the circumnuclear HII regions probably requires deeper observations.

We present a method to include the effects of gravity in the plasma physics code Cloudy. More precisely, a term is added to the desired gas pressure in order to enforce hydrostatic equilibrium, accounting for both the self-gravity of the gas and the presence of an optional external potential. As a test case, a plane-parallel model of the vertical structure of the Milky Way disk near the solar neighbourhood is considered. It is shown that the gravitational force determines the scale height of the disk, and it plays a critical role in setting its overall chemical composition. However, other variables, such as the shape of incident continuum and the intensity of the Galactic magnetic field, strongly affect the predicted structure.

We report the detection of a high Galactic latitude, large scale, 7-sigma signal in WMAP 5yr and spatially correlated with the ecliptic plane. Two possible candidates are studied, namely unresolved sources and Zodiacal light emission. We determine the strength of the Zodiacal light emission at WMAP frequencies and estimate the contribution from unresolved extragalactic sources. Neither the standard Zodiacal light emission nor the unresolved sources alone seem to be able to explain the observed signal. Other possible interpretations like Galactic foregrounds and diffuse Sunyaev-Zel'dovich effect also seem unlikely. We check if our findings could affect the low-l anomalies which have been reported in the WMAP data. Neither Zodiacal light emission nor unresolved point source residuals seem to affect significantly the quadrupole and octupole measurements. However, a signal with a quasi-blackbody spectrum and with a spatial distribution similar to the Zodiacal light emission, could explain both the anomalous signal and the low-ell anomalies. Future data (Planck) will be needed in order to explain the origin of this signal.

This paper describes a new algorithm for Monte Carlo integration, based on the Field Estimator for Arbitrary Spaces (FiEstAS). The algorithm is discussed in detail, and its performance is evaluated in the context of Bayesian analysis, with emphasis on multimodal distributions with strong parameter degeneracies. Source code is available upon request.

WMAP observations at mm wavelengths are sensitive to the Sunyaev-Zel'dovich effect in galaxy clusters. Among all the objects in the sky, the Virgo cluster is expected to provide the largest integrated signal. Based on models compatible with the X-ray emission observed in the ROSAT All Sky Survey, we predict a two-sigma detection of the SZ effect from Virgo in the WMAP 3-year data. Our analysis reveals a 3-sigma signal on scales of 5 degrees, although the frequency dependence deviates from the theoretical expectation for the SZ effect. The main sources of uncertainty are instrumental noise, and most importantly, possible contamination from point sources and diffuse back/foregrounds. In particular, a population of unresolved extragalactic sources in Virgo would explain the observed intensity and frequency dependence. In order to resolve this question one needs to wait for experiments like Planck to achieve the required accuracy.

Thanks to the ever increasing computational power and the development of more sophisticated algorithms, numerical N-body simulations are now uncovering several phenomenological relations between the physical properties of dark matter haloes in position and velocity space. It is the aim of the present work to investigate in detail the dynamical structure of dark matter haloes, as well as its possible dependence on mass and its evolution with redshift up to z=5. We use high-resolution cosmological simulations of individual objects to compute the radially-averaged profiles of several quantities, scaled by the radius Rmax at which the circular velocity attains its maximum value, Vmax. No systematic dependence on mass or cosmic epoch are found within Rmax, and all the different radial profiles are well fit by simple analytical models. However, our results suggest that several properties are not `universal' outside this radius. [Abridged]

We present a simple model to describe the dark matter density, the gas density, and the gas temperature profiles of galaxy clusters. Analytical expressions for these quantities are given in terms of only five free parameters with a clear physical meaning: the mass M of the dark matter halo (or the characteristic temperature T_0), the characteristic scale radius a, the cooling radius in units of a (0<alpha<1), the central temperature in units of T_0 (0<t<1), and the asymptotic baryon fraction in units of the cosmic value (f~1). It is shown that our model is able to reproduce the three-dimensional density and temperature profiles inferred from X-ray observations of real clusters within a 20 per cent accuracy over most of the radial range. Some possible applications are briefly discussed.

In the current paradigm of cosmic structure formation, dark matter plays a key role on the formation and evolution of galaxies through its gravitational influence. On microscopic scales, dark matter particles are expected to annihilate amongst themselves into different products, with some fraction of the energy being transferred to the baryonic component. It is the aim of the present work to show that, in the innermost regions of dark matter halos, heating by dark matter annihilation may be comparable to the cooling rate of the gas. We use analytical models of the dark matter and gas distributions in order to estimate the heating and cooling rates, as well as the energy available from supernova explosions. Depending on the model parameters and the precise nature of dark matter particles, the injected energy may be enough to balance radiative cooling in the cores of galaxy clusters. On galactic scales, it would inhibit star formation more efficiently than supernova feedback. Our results suggest that dark matter annihilation prevents gas cooling and star formation within at least $0.01-1$ per cent of the virial radius.

We test the Secondary Infall Model (SIM) by direct comparison with the results of N-body simulations. Eight cluster-size and six galactic-size dark matter haloes have been selected at $z=0$ and re-simulated with high resolution. Based on their density profiles at the initial redshift, we compute their evolution by the SIM, assuming a simple prescription for the angular momentum. A comparison of the density profiles obtained by the SIM and the numerical experiments at $z=5$, 1 and 0 shows that, for most of the haloes at most epochs, the SIM reproduces the simulated mater distribution with a typical fractional deviation of less than 40 per cent over more than six order of magnitudes in the density. It is also found that, within the SIM framework, most of the diversity in the shape of the density profiles at $z=0$ arises from the scatter in the primordial initial conditions rather than the scatter in the angular momentum distribution. A crude optimization shows that a similar degree of agreement is obtained for galactic and cluster haloes, but the former seem to require slightly higher amounts of angular momentum. Our main conclusion is that the SIM provides a viable dynamical model for predicting the structure and evolution of the density profile of dark matter haloes.

Chandra X-ray observations revealed the presence of cold fronts (sharp contact discontinuities between gas regions with different temperatures and densities) in the centers of many, if not most, relaxed clusters with cool cores. We use high-resolution simulations of idealized cluster mergers to show that they are due to sloshing of the cool gas in the central gravitational potential, which is easily set off by minor mergers and can persist for gigayears. The only necessary condition is a steep entropy profile, as observed in cooling flow clusters. Even if the infalling subcluster has no gas during core passage, the gravitational disturbance sets the main mass peak (gas and dark matter together) in motion relative to the surrounding gas. A rapid change in the direction of motion causes a change in ram pressure, which pushes the cool gas away from the dark matter peak and triggers sloshing. For nonzero impact parameters, the cool gas acquires angular momentum, resulting in a characteristic spiral pattern of cold fronts. There is little visible disturbance outside the cool core in such a merger. If the subcluster retains its gas during core passage, the cool central gas of the main cluster is more easily decoupled from the dark matter peak. Subsequently, some of that gas, and often the cool gas from the subcluster, falls back to the center and starts sloshing. However, in such a merger, global disturbances are readily visible in X-rays for a long time. We conclude that cold fronts at the centers of relaxed clusters, often spiral or concentric-arc in shape, are probably caused by encounters with small subhalos stripped of all their gas at the early infall stages.

About one year ago, it was speculated that decaying or annihilating Light Dark Matter (LDM) particles could explain the flux and extension of the 511 keV line emission in the galactic centre. Here we present a thorough comparison between theoretical expectations of the galactic positron distribution within the LDM scenario and observational data from INTEGRAL/SPI. Unlike previous analyses, there is now enough statistical evidence to put tight constraints on the shape of the dark matter halo of our galaxy, if the galactic positrons originate from dark matter. For annihilating candidates, the best fit to the observed 511 keV emission is provided by a radial density profile with inner logarithmic slope gamma=1.03+-0.04. In contrast, decaying dark matter requires a much steeper density profile, gamma>1.5, rather disfavoured by both observations and numerical simulations. Within the annihilating LDM scenario, a velocity-independent cross-section would be consistent with the observational data while a cross-section purely proportional to v^2 can be rejected at a high confidence level. Assuming the most simplistic model where the galactic positrons are produced as primaries, we show that the LDM candidate should be a scalar rather than a spin-1/2 particle and obtain a very stringent constraint on the value of the positron production cross-section to explain the 511 keV emission. One consequence is that the value of the fine structure constant should differ from that recommended in the CODATA. This is a very strong test for the LDM scenario and an additional motivation in favour of experiments measuring alpha directly. Our results finally indicate that an accurate measurement of the shape of the dark halo profile could have a tremendous impact on the determination of the origin of the 511 keV line and vice versa.

The aim of the present work is to show that, contrary to popular belief, galaxy clusters are **not** expected to be self-similar, even when the only energy sources available are gravity and shock-wave heating. In particular, we investigate the scaling relations between mass, luminosity and temperature of galaxy groups and clusters in the absence of radiative processes. Theoretical expectations are derived from a polytropic model of the intracluster medium and compared with the results of high-resolution adiabatic gasdynamical simulations. It is shown that, in addition to the well-known relation between the mass and concentration of the dark matter halo, the effective polytropic index of the gas also varies systematically with cluster mass, and therefore neither the dark matter nor the gas profiles are exactly self-similar. It is remarkable, though, that the effects of concentration and polytropic index tend to cancel each other, leading to scaling relations whose logarithmic slopes roughly match the predictions of the most basic self-similar models. We provide a phenomenological fit to the relation between polytropic index and concentration, as well as a self-consistent scheme to derive the non-linear scaling relations expected for any cosmology and the best-fit normalizations of the M-T, L-T and F-T relations appropriate for a Lambda-CDM universe. The predicted scaling relations reproduce observational data reasonably well for massive clusters, where the effects of cooling and star formation are expected to play a minor role.

[Abridged] We present a novel technique, dubbed FiEstAS, to estimate the underlying density field from a discrete set of sample points in an arbitrary multidimensional space. FiEstAS assigns a volume to each point by means of a binary tree. Density is then computed by integrating over an adaptive kernel. As a first test, we construct several Monte Carlo realizations of a Hernquist profile and recover the particle density in both real and phase space. At a given point, Poisson noise causes the unsmoothed estimates to fluctuate by a factor ~2 regardless of the number of particles. This spread can be reduced to about 1 dex (~26 per cent) by our smoothing procedure. [...] We conclude that our algorithm accurately measure the phase-space density up to the limit where discreteness effects render the simulation itself unreliable. Computationally, FiEstAS is orders of magnitude faster than the method based on Delaunay tessellation that Arad et al. employed, making it practicable to recover smoothed density estimates for sets of 10^9 points in 6 dimensions.

We present results from a set of high (512^3 effective resolution), and ultra-high (1024^3) SPH adiabatic cosmological simulations of cluster formation aimed at studying the internal structure of the intracluster medium (ICM). We derive a self-consistent analytical model of the structure of the intracluster medium (ICM). We discuss the radial structure and scaling relations expected from purely gravitational collapse, and show that the choice of a particular halo model can have important consequences on the interpretation of observational data. The validity of the approximations of hydrostatic equilibrium and a polytropic equation of state are checked against results of our simulations. The properties of the ICM are fully specified when a 'universal' profile is assumed for either the dark or the baryonic component. We also show the first results from an unprecedented large-scale simulation of 500 Mpc/h and 2 times 512^3 gas and dark matter particles. This experiment will make possible a detailed study of the large-scale distribution of clusters as a function of their X-ray properties.

Simple self-consistent models of galaxy groups and clusters are tested against the results of high-resolution adiabatic gasdynamical simulations. We investigate two models based on the existence of a 'universal' dark matter density profile and two versions of the beta-model. The mass distribution of relaxed clusters can be fitted by phenomenological formulae proposed in the literature. Haloes that have experienced a recent merging event are systematically less concentrated and show steeper profiles than relaxed objects near the centre. The hot X-ray emitting gas is found to be in approximate hydrostatic equilibrium with the dark matter potential, and it is well described by a polytropic equation of state. Analytic formulae for the gas density and temperature can be derived from these premises. Though able to reproduce the X-ray surface brightness, the beta-model is shown to provide a poor description of our numerical clusters. We find strong evidence of a 'universal' temperature profile that decreases by a factor of 2-3 from the centre to the virial radius, whereas baryon fraction and entropy are monotonically increasing functions. Numerical resolution and entropy conservation play a key role in the shapes of the profiles at small radii.

The radial mass distribution of dark matter haloes is investigated within the framework of the spherical infall model. We present a new formulation of spherical collapse including non-radial motions, and compare the analytical profiles with a set of high-resolution N-body simulations ranging from galactic to cluster scales. We argue that the dark matter density profile is entirely determined by the initial conditions, which are described by only two parameters: the height of the primordial peak and the smoothing scale. These are physically meaningful quantities in our model, related to the mass and formation time of the halo. Angular momentum is dominated by velocity dispersion, and it is responsible for the shape of the density profile near the centre. The phase-space density of our simulated haloes is well described by a power-law profile, rho/sigma^3 = 10^{1.46\pm0.04} (rho_c/Vvir^3) (r/Rvir)^{-1.90\pm0.05}. Setting the eccentricity of particle orbits according to the numerical results, our model is able to reproduce the mass distribution of individual haloes.

We present some of the results of an ongoing collaboration to sudy the dynamical properties of galaxy clusters by means of high resolution adiabatic SPH cosmological simulations. Results from our numerical clusters have been tested against analytical models often used in X-ray observations: $\beta$ model (isothermal and polytropic) and those based on universal dark matter profiles. We find a universal temperature profile, in agreement with AMR gasdynamical simulations of galaxy clusters. Temperature decreases by a factor 2-3 from the center to virial radius. Therefore, isothermal models (e.g. $\beta$ model) give a very poor fit to simulated data. Moreover, gas entropy profiles deviate from a power law near the center, which is also in very good agreement with independent AMR simulations. Thus, if future X-ray observations confirm that gas in clusters has an extended isothermal core, then non-adiabatic physics would be required in order to explain it.

Clusters of galaxies are studied from a theoretical point of view, comparing with observational results whenever possible. The problem is approached both analytically as well as by means of high-resoultion numerical simulations. The dark matter halo, the hot intracluster gas, and the stellar component are investigated separately. Numerical clusters are consistent with a relatively simple scenario, in which these objects form around local maxima of the primordial density field, smoothed on Mpc sacales. Hot diffuse gas is in approximate hydrostatic equilibrium with the dark matter potential, and it is well described by a polytropic equation of state. Global X-ray properties are closely interrelated, and deviations from self-similarity are expected even when radiative processes are not considered. Star formation is severely reduced in the inner regions of clusters, due to the fact that infalling galaxies loose their hot gas reservoirs during the first orbit.

The cosmic star formation history in Cold Dark Matter dominated cosmological scenarios is studied by means of hydrodynamical numerical simulations. In particular, we explore a low density model with a Lambda-term and two high density models with different power spectra, all of them being spatially flat. Our simulations employ a fully nonlinear N-body and Eulerian hydrodynamics algorithm with a model for star formation and supernovae feedback that depends on two phenomenological parameters determined in agreement with previous papers. We find a nearly constant star formation rate beyond z=1, and we discuss which facts may determine the decrease in the SFR from z=1 to the present epoch. The Lambda-term cosmology with realistic parameters for star formation and feedback best reproduces the observed star formation history.

From a set of 3D cosmological simulations which incorporate a self-consistent model of the star-gas interactions, we have been able to obtain a statistically significant sample of galaxy-type halos with observational properties, like colors and luminosities. Using this data-base, we have studied general relations of galaxies, such as the Tully-Fisher relation, luminosity functions, environmental dependences, or cosmic star formation density.