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

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 .

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.

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.

With the ever increasing size and complexity of fully self-consistent simulations of galaxy formation within the framework of the cosmic web, the demands upon object finders for these simulations has simultaneously grown. To this extent we initiated the Halo Finder Comparison Project that gathered together all the experts in the field and has so far led to two comparison papers, one for dark matter field haloes (Knebe et al. 2011), and one for dark matter subhaloes (Onions et al. 2012). However, as state-of-the-art simulation codes are perfectly capable of not only following the formation and evolution of dark matter but also account for baryonic physics (e.g. hydrodynamics, star formation, feedback) object finders should also be capable of taking these additional processes into consideration. Here we report on a comparison of codes as applied to the Constrained Local UniversE Simulation (CLUES) of the formation of the Local Group which incorporates much of the physics relevant for galaxy formation. We compare both the properties of the three main galaxies in the simulation (representing the MW, M31, and M33) as well as their satellite populations for a variety of halo finders ranging from phase-space to velocity-space to spherical overdensity based codes, including also a mere baryonic object finder. We obtain agreement amongst codes comparable to (if not better than) our previous comparisons, at least for the total, dark, and stellar components of the objects. However, the diffuse gas content of the haloes shows great disparity, especially for low-mass satellite galaxies. This is primarily due to differences in the treatment of the thermal energy during the unbinding procedure. We acknowledge that the handling of gas in halo finders is something that needs to be dealt with carefully, and the precise treatment may depend sensitively upon the scientific problem being studied.

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.

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.