Search

Galaxies in cosmic voids have been reported with properties related to a delayed evolution with respect to the Universe in general. These characteristics reflect the interaction of galaxies with the environment. However, it is not clear the degree of influence of the large-scale structure on the properties of void galaxies or, if these are only influenced by the low local density around them typical of these regions. In this article we identified cosmic voids in the SDSS-DR16 and studied various properties of galaxies, such as g-r colour, star formation rate, and concentration. To characterise the local environment, we have identified groups of galaxies and studied their properties as a function of their dark matter and stellar masses, analysing separately those found in voids and in the general sample. Our results show that galaxies that inhabit haloes of a given mass (below \sim 10^13.5 M_\dot ), are bluer, have a higher star formation rate and are less concentrated when the host halo is inside voids compared to other regions. For larger halo masses, the trend disappears. We also analyse whether the properties of galaxies are sensitive to the type of voids that inhabit. This is done by separating voids embedded in overdense regions (S-type) from those that asymptotically converge to the average density of the universe (R-type). We found that galaxies in R-type voids are bluer, with higher SFR and less concentration than in S-type voids. Our results indicate some degree of correlation of galaxy properties with the large-scale environment provided by voids, suggesting possible second-order mechanisms in galaxy evolution., 10 pages, 9 figures. Acccepted for publication in MNRAS

All estimators of the two-point correlation function are based on a random catalogue, a set of points with no intrinsic clustering following the selection function of a survey. High-accuracy estimates require the use of large random catalogues, which imply a high computational cost. We propose to replace the standard random catalogues by glass-like point distributions or glass catalogues, which are characterized by a power spectrum $P(k)\propto k^4$ and exhibit significantly less power than a Poisson distribution with the same number of points on scales larger than the mean inter-particle separation. We show that these distributions can be obtained by iteratively applying the technique of Zeldovich reconstruction commonly used in studies of baryon acoustic oscillations (BAO). We provide a modified version of the widely used Landy-Szalay estimator of the correlation function adapted to the use of glass catalogues and compare its performance with the results obtained using random samples. Our results show that glass-like samples do not add any bias with respect to the results obtained using Poisson distributions. On scales larger than the mean inter-particle separation of the glass catalogues, the modified estimator leads to a significant reduction of the variance of the Legendre multipoles $\xi_\ell(s)$ with respect to the standard Landy-Szalay results with the same number of points. The size of the glass catalogue required to achieve a given accuracy in the correlation function is significantly smaller than when using random samples. Even considering the small additional cost of constructing the glass catalogues, their use could help to drastically reduce the computational cost of configuration-space clustering analysis of future surveys while maintaining high-accuracy requirements., Comment: 9 pages, 7 figures, minor changes to match version accepted by MNRAS

Associations of dwarf galaxies are loose systems composed exclusively of dwarf galaxies. These systems were identified in the Local Volume for the first time more than 30 yr ago. We study these systems in the cosmological framework of the λ cold dark matter (λCDM) model.We consider the Small MultiDark Planck simulation and populate its dark matter haloes by applying the semi-analytic model of galaxy formation SAG. We identify galaxy systems using a friends-of-friends algorithm with a linking length equal to b = 0.4Mpc h-1to reproduce the size of dwarf galaxy associations detected in the Local Volume. Our samples of dwarf systems are built up removing those systems that have one or more galaxies with stellar mass larger than a maximum thresholdMmax.We analyse three different samples defined by log10(Mmax[M⊙ h-1]) = 8.5, 9.0, and 9.5. On average, our systems have typical sizes of ∼ 0.2Mpc h-1, velocity dispersion of ∼ 30km s-1, and estimated total mass of ∼ 1011M⊙ h-1. Such large typical sizes suggest that individual members of a given dwarf association reside in different dark matter haloes and are generally not substructures of any other halo. Indeed, in more than 90 per cent of our dwarf systems their individual members inhabit different dark matter haloes, while only in the remaining 10 per cent members do reside in the same halo. Our results indicate that the λCDM model can naturally reproduce the existence and properties of dwarf galaxies' associations without much difficulty. Fil: Yaryura, Claudia Yamila. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; Argentina Fil: Abadi, Mario Gabriel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; Argentina Fil: Gottlöber, Stefan. Leibniz Universitat Hannover; Alemania Fil: Libeskind, Noam I.. Leibniz Universitat Hannover; Alemania Fil: Cora, Sofia Alejandra. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Astrofísica La Plata. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas. Instituto de Astrofísica La Plata; Argentina Fil: Ruiz, Andrés Nicolás. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; Argentina. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba; Argentina Fil: Vega Martínez, Cristian Antonio. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Astrofísica La Plata. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas. Instituto de Astrofísica La Plata; Argentina. Universidad de La Serena; Chile Fil: Yepes, Gustavo. Universidad Autónoma de Madrid; España Fil: Behroozi, Peter. George Mason University. School Of Physics. Astronomy And Computational Sciences; Estados Unidos

Recent works have shown that the properties of galaxy populations in dark matter halos vary with large-scale environments. These results suggest a variation in the halo occupation distribution (HOD) in extreme density environments. To analyse these effects, we identify cosmic voids and future virialised structures (FVS) in the SDSS-DR12 and estimate the HOD within these superstructures using group catalogues as dark matter halo proxies. Our goal is to use observational galaxy data to characterise the HOD within voids and FVS, explore the different properties of these galaxies populations and compare them with the general results outside these superstructures. Using a galaxy group catalogue we compute the HOD within both types of superstructures. We also study the dependence on the results on the main void and FVS properties. We also analysed the mean stellar age of the galaxies inside these regions. In all cases, we compare the results with those derived from the Field sample. Inside voids, we find a strong decrease in HOD concerning the Field results. The mean number of satellites fall to 50%. Inside FVS, the HOD shows a significant increase to the Field, with a 40% excess in the mean number of satellites. In both regions, the differences with respect to the Field increases for the extreme values of the density environments. We obtain no signs of variations related to intrinsic characteristics of voids and FVS. We find that the cumulative distribution of the mean age of stars of the central galaxy also varies in these regions. Finally, we explore the HOD for the 25% youngest (oldest) galaxies. We find that for the low-mass groups the youngest galaxies are only present inside voids. On the other hand, for the high-mass groups the FVS environments show the same increase in the HOD concerning the Field. We find that cosmic voids lack of oldest galaxies., Comment: 9 pages, 9 figures. First version. Sent to refereed in A&A

We connect galaxy properties with their orbital classification by analysing a sample of galaxies with stellar mass $M_{\star} \geq 10^{8.5}h^{-1}M_\odot$ residing in and around massive and isolated galaxy clusters with mass $M_{200} > 10^{15}h^{-1}M_\odot$ at redshift $z=0$. The galaxy population is generated by applying the semi-analytic model of galaxy formation SAG on the cosmological simulation MultiDark Planck 2. We classify galaxies considering their real orbits (3D) and their projected phase-space position using the ROGER code (2D). We define five categories: cluster galaxies, galaxies that have recently fallen into a cluster, backsplash galaxies, infalling galaxies, and interloper galaxies. For each class, we analyse the $g-r$ colour, the specific star formation rate (sSFR), and the stellar age, as a function of the stellar mass. For the 3D classes, we find that cluster galaxies have the lowest sSFR, and are the reddest and the oldest, as expected from environmental effects. Backsplash galaxies have properties intermediate between the cluster and recent infaller galaxies. For each 2D class, we find an important contamination by other classes. We find it necessary to separate the galaxy populations in red and blue to perform a more realistic analysis of the 2D data. For the red population, the 2D results are in good agreement with the 3D predictions. Nevertheless, when the blue population is considered, the 2D analysis only provides reliable results for recent infallers, infalling galaxies and interloper galaxies., 5 figures, 12 pages, accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

The environment within dark matter haloes can quench the star formation of galaxies. However, environmental effects beyond the virial radius of haloes ($\gtrsim$ 1 Mpc) are less evident. An example is the debated correlation between colour or star formation in central galaxies and neighbour galaxies in adjacent haloes at large separations of several Mpc, referred to as two-halo galactic conformity. We use two galaxy catalogues generated from different versions of the semi-analytic model SAG applied to the MDPL2 cosmological simulation and the IllustrisTNG300 cosmological hydrodynamical simulation to study the two-halo conformity by measuring the quenched fraction of neighbouring galaxies as a function of the real-space distance from central galaxies. We find that low-mass central galaxies in the vicinity of massive systems ($M_{\rm 200c}$ $\geq$ 10$^{13}$ $h^{-1}~\rm M_{\odot}$) out to 5 $h^{-1}$ Mpc are preferentially quenched compared to other central galaxies at fixed stellar mass $M_{\star}$ or fixed host halo mass $M_{\rm 200c}$ at $z$ ~ 0. In all the galaxies catalogues is consistent that the low-mass ($M_{\star} < 10^{10}$ $h^{-1}~\rm M_{\odot}$ or $M_{\rm 200c} < 10^{11.8}$ $h^{-1}~\rm M_{\odot}$) central galaxies in the vicinity of clusters and, especially, groups of galaxies mostly produce the two-halo galactic conformity. On average, the quenched low-mass central galaxies are much closer to massive haloes than star-forming central galaxies of the same mass (by a factor of ~5). Our results agree with other works regarding the environmental influence of massive haloes that can extend beyond the virial radius and affect nearby low-mass central galaxies., 15 pages and 10 figures without appendix. Accepted for publication in MNRAS

We present a new approach to describe statistics of the non-linear matter density field that exploits a degeneracy in the impact of different cosmological parameters on the linear dimensionless matter power spectrum, $\Delta^2_{\rm L}(k)$. We classify all cosmological parameters into two groups, shape parameters, which determine the shape of $\Delta^2_{\rm L}(k)$, and evolution parameters, which only affect its amplitude at any given redshift. With this definition, the time evolution of $\Delta^2_{\rm L}(k)$ in models with identical shape parameters but different evolution parameters can be mapped from one to the other by relabelling the redshifts that correspond to the same clustering amplitude, which we characterize by the linear mass fluctuation in spheres of radius $12\,{\rm Mpc}$, $\sigma_{12}(z)$. We use N-body simulations to show that the same evolution mapping relation gives a good description of the non-linear power spectrum, the halo mass function, or the full density field. The deviations from the exact degeneracy are the result of the different structure formation histories experienced by each model to reach the same clustering amplitude and can be accurately described in terms of differences in the suppression factor $g(a) = D(a)/a$. These relations can be used to drastically reduce the number of parameters required to describe the cosmology dependence of the power spectrum. We show how this can help to speed up the inference of parameter constraints from cosmological observations. We also present a new design of an emulator of the non-linear power spectrum whose predictions can be adapted to an arbitrary choice of evolution parameters and redshift., Comment: 13 pages, 11 figures, replaced to match version accepted by MNRAS

There are hints suggesting that properties of galaxy populations in dark matter haloes may depend on their large-scale environment. Recent works point out that very low-density environments influence halo occupation distribution (HOD), however there is not a similar analysis focused on high-density environments. Here we use a simulated set of future virialized superstructures (FVS) to analyse the occupation of galaxies in haloes within these high globally dense regions. We use a publicly available simulated galaxy set constructed with a semi-analytical model to identify FVS in the simulation. Then, we computed the HOD within these superstructures for different absolute magnitude thresholds and make several analysis including the comparison to the global HOD results. We study the dependence on the results on properties of the FVS such as density and volume as well as consider the morphology of galaxies. We also analysed the properties of the stellar content of galaxies and the formation time of the haloes inside FVS. We find a significant increase in the HOD inside FVS. This result is present for all absolute magnitude thresholds explored. The effect is larger in the densest regions of FVS, but does not depend on the volume of the superstructure. We also find that the stellar-mass content of galaxies considerably differs inside the superstructures. Low mass haloes have their central and satellite galaxies with a higher stellar mass content (50%), and exhibit mean star ages (20%) older than average. For massive haloes in FVS we find that only the stellar mass of satellite galaxies varies considerably corresponding to a decrease of 50%. We find a significant statistical difference between the formation times of haloes in FVS and the average population. Haloes residing in superstructures formed earlier, a fact that leads to several changes in the HOD and their member galaxy properties., 12 pages, 15 figure. Published version (by the A&A)

Context. The extragalactic vision we have through the Milky Way is very unclear. There is significant extinction of the optical emission from objects located in the region called the zone of avoidance (ZOA). NIR wavelengths are less affected by extinction, and therefore the infrared surveys in this zone are a potential source of astronomical discoveries. Nevertheless, these observations need to be compared with cosmological simulations in order to carry out high−accuracy studies. Aims. Our aim is to identify extragalactic sources in the ZOA, using infrared images of the VVV survey. We consider mock galaxy catalogues in order to interpret observational results. Methods. We studied a region of 1.636 square degrees corresponding to the VVV tile b204. Using SExtractor, we analysed photometric data generating a catalogue of extended sources in this area. In order to confirm these sources as galaxy candidates we visually inspected RGB images looking for typical galaxy features. Using 2MASX and GCMW catalogued sources we tested completeness and contamination of our catalogue and define suitable colour cuts to select galaxies. We also compared the observational results with those obtained from two semi-analytical models on Dark Matter simulations. One galaxy catalogue was constructed with the SAG semi-analytic model of galaxy formation, and the other one was constructed with the L-Galaxies semi-analytic model. Results. By adopting CLASS_STAR < 0.5, r1/2 > 0.7 arcsec and specific colour cuts (J − Ks > 0.97, J − H > 0 and H − Ks > 0) we generated an automatic catalogue of extended sources. After visual inspection we identified 624 sources with 10, Instituto de Astrofísica de La Plata

This is the second part of a thorough investigation of the redshift-space effects that affect void properties and the impact they have on cosmological tests. Here, we focus on the void-galaxy cross-correlation function, specifically, on the projected versions that we developed in a previous work. The pillar of the analysis is the one-to-one relationship between real and redshift-space voids above the shot-noise level identified with a spherical void finder. Under this mapping, void properties are affected by three effects: (i) a systematic expansion as a consequence of the distortions induced by galaxy dynamics, (ii) the Alcock-Paczynski volume effect, which manifests as an overall expansion or contraction depending on the fiducial cosmology, and (iii) a systematic off-centring along the line of sight as a consequence of the distortions induced by void dynamics. We found that correlations are also affected by an additional source of distortions: the ellipticity of voids. This is the first time that distortions due to the off-centring and ellipticity effects are detected and quantified. With a simplified test, we verified that the Gaussian streaming model is still robust provided all these effects are taken into account, laying the foundations for improvements in current models in order to obtain unbiased cosmological constraints from spectroscopic surveys. Besides this practical importance, this analysis also encodes key information about the structure and dynamics of the Universe at the largest scales. Furthermore, some of the effects constitute cosmological probes by themselves, as is the case of the void ellipticity., Comment: 14 pages, 9 figures, published in MNRAS, accepted after minor comments

Studies of large-scale structures in the Universe, such as superstructures or cosmic voids, have been widely used to characterize the properties of the cosmic web through statistical analyses. On the other hand, the 2-point correlation function of large-scale tracers such as galaxies or halos provides a reliable statistical measure. However, this function applies to the spatial distribution of point-like objects, and therefore it is not appropriate for extended large structures which strongly depart from spherical symmetry. Here we present an analysis based on the standard correlation function formalism that can be applied to extended objects exhibiting arbitrary shapes. Following this approach, we compute the probability excess $\Xi$ of having spheres sharing parts of cosmic structures with respect to a realization corresponding to a distribution of the same structures in random positions. For this aim, we identify superstructures defined as Future Virialized Structures (FVSs) in semi-anaytic galaxies on the MPDL2 MultiDark simulation. We have also identified cosmic voids to provide a joint study of their relative distribution with respect to the superstructures. Our analysis suggests that $\Xi$ provides useful characterizations of the large scale distribution, as suggested from an analysis of sub-sets of the simulation. Even when superstructure properties may exhibit negligible variations across the sub-sets, $\Xi$ has the sensitivity to statistically distinguish sub-boxes that departs from the mean at larger scales. Thus, our methods can be applied in analysis of future surveys to provide characterizations of large-scale structure suitable to distinguish different theoretical scenarios., Comment: 8 pages, 6 figures. Accepted for publication in MNRAS

We present the ROGER (Reconstructing Orbits of Galaxies in Extreme Regions) code, which uses three different machine learning techniques to classify galaxies in, and around, clusters, according to their projected phase-space position. We use a sample of 34 massive, $M_{200}>10^{15} h^{-1} M_{\odot}$, galaxy clusters in the MultiDark Planck 2 (MDLP2) simulation at redshift zero. We select all galaxies with stellar mass $M_{\star} \ge 10^{8.5} h^{-1}M_{\odot}$, as computed by the semi-analytic model of galaxy formation SAG, that are located in, and in the vicinity of, the clusters and classify them according to their orbits. We train ROGER to retrieve the original classification of the galaxies out of their projected phase-space positions. For each galaxy, ROGER gives as output the probability of being a cluster galaxy, a galaxy that has recently fallen into a cluster, a backsplash galaxy, an infalling galaxy, or an interloper. We discuss the performance of the machine learning methods and potential uses of our code. Among the different methods explored, we find the K-Nearest Neighbours algorithm achieves the best performance., Acceptep for its publication in the MNRAS Journal. Code available at github repository

Evidence shows properties of dark matter haloes may vary with large-scale environment. By studying the halo occupation distribution in cosmic voids it is possible to obtain useful information that can shed light on the subject. The history of the formation of the haloes and galaxies residing in these regions is likely to differ from the global behaviour given their extreme environment. We use two public access simulated galaxy catalogues constructed with different methods: a semi-analytical model and a hydrodynamic simulation. In both, we identify cosmic voids and we measure the halo occupation distribution inside these regions for different absolute magnitude thresholds. We compare these determinations with the overall results and we study the dependence of different characteristics of the voids. Also, we analyze the stellar content and the formation time of the haloes inside voids and confront the general halo population results. Inside the voids, we find a significantly different halo occupation distribution with respect to the general results. This is present in all absolute magnitude ranges explored. We obtain no signs of variation related to void characteristics indicating that the effects depend only on the density of the large-scale environment. Additionally, we find that the stellar mass content also differs within voids, which host haloes with less massive central galaxies (10%) as well as satellites with significantly lower stellar mass content (30%). Finally, we find a slight difference between the formation times of the haloes which are younger in voids than the average population. These characteristics indicate that haloes populating voids have had a different formation history, inducing significant changes on the halo occupation distribution., 11 pages and 12 figures. Submitted to A&A. Comments are welcome

We present a new cosmological test using the distribution of galaxies around cosmic voids without assuming a fiducial cosmology. The test is based on a physical model for the void-galaxy cross-correlation function projected along and perpendicular to the line of sight. We treat correlations in terms of void-centric angular distances and redshift differences between void-galaxy pairs, hence it is not necessary to assume a fiducial cosmology. This model reproduces the coupled dynamical (Kaiser effect, RSD) and geometrical (Alcock-Paczynski effect, GD) distortions that affect the correlation measurements. It also takes into account the scale mixing due to the projection ranges in both directions. The model is general, so it can be applied to an arbitrary cylindrical binning scheme, not only in the case of the projected correlations. It primarily depends on two cosmological parameters: $\Omega_m$, the matter fraction of the Universe today (sensitive to GD), and $\beta$, the ratio between the growth rate factor of density perturbations and the tracer bias (sensitive to RSD). In the context of the new generation of galaxy spectroscopic surveys, we calibrated the test using the Millennium XXL simulation for different redshifts. The method successfully recovers the cosmological parameters. We studied the effect of measuring with different projection ranges, finding robust results up to wide ranges. The resulting data covariance matrices are relatively small, which reduces the noise in the Gaussian likelihood analysis and will allow the usage of a smaller number of mock catalogues. The performance evaluated in this work indicates that the developed method is a promising test to be applied on real data., Comment: 13 pages. 7 Figures. Accepted in MNRAS

We present results from the semi-analytic model of galaxy formation SAG applied on the MULTIDARK simulation MDPL2. SAG features an updated supernova (SN) feedback scheme and a robust modelling of the environmental effects on satellite galaxies. This incorporates a gradual starvation of the hot gas halo driven by the action of ram pressure stripping (RPS), that can affect the cold gas disc, and tidal stripping (TS), which can act on all baryonic components. Galaxy orbits of orphan satellites are integrated providing adequate positions and velocities for the estimation of RPS and TS. The star formation history and stellar mass assembly of galaxies are sensitive to the redshift dependence implemented in the SN feedback model. We discuss a variant of our model that allows to reconcile the predicted star formation rate density at z ≳ 3 with the observed one, at the expense of an excess in the faint end of the stellar mass function at z= 2. The fractions of passive galaxies as a function of stellar mass, halo mass, and the halo-centric distances are consistent with observational measurements. The model also reproduces the evolution of the main sequence of star-forming central and satellite galaxies. The similarity between them is a result of the gradual starvation of the hot gas halo suffered by satellites, in which RPS plays a dominant role. RPS of the cold gas does not affect the fraction of quenched satellites but it contributes to reach the right atomic hydrogen gas content for more massive satellites (M* ≳1010M⊙)., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MULTIDARK simulations (i.e. MDPL2, with a volume of (1 h-1 Gpc)3 and mass resolution of 1.5 × 109 h-1 M⊙) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite-mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by ongoing and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available., Facultad de Ciencias Astronómicas y Geofísicas, Instituto de Astrofísica de La Plata

We study the stellar populations of bulges of MilkyWay-like (MW-like) galaxies with the aim of identifying the physical processes involved in the formation of the bulge of our Galaxy. We use the semi-analytic model of galaxy formation and evolution SAG adapted to this aim; these kind of models can trace the properties of galaxies and their components like stellar discs, bulges and haloes, but resolution limits prevent them from reaching the scale of stellar populations (SPs). Properties of groups of stars formed during single star formation events are stored and tracked in the model and results are compared with observations of stars in the galactic bulge. MW-like galaxies are selected using two different criteria. One of them considers intrinsic photo-metric properties and the second is focused on the cosmological context of the local group (LG) of galaxies. We compare our model results with spectroscopic and photometric stellar metallicity distributions. We find that 87 per cent of stars in bulges of MW-type galaxies in our model are accreted and formed in starbursts during disc instability events. Mergers contribute to 13 per cent of the mass budget of the bulge and are responsible for the low metallicity tail of the distribution. Abundance ratios of α elements with respect to iron, [α/Fe], are measured in SPs of model galaxies. The patterns found in the model for SPs with different origins help to explain the lack of a gradient of [α/Fe] ratios in observed stars along the minor axis of the bulge., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

We compute the bulk motions of cosmic voids, using a $\Lambda$CDM numerical simulation considering the mean velocities of the dark matter inside the void itself and that of the haloes in the surrounding shell. We find coincident values of these two measures in the range $\sim$ 300-400 km/s, not far from the expected mean peculiar velocities of groups and galaxy clusters. When analysing the distribution of the pairwise relative velocities of voids, we find a remarkable bimodal behaviour consistent with an excess of both systematically approaching and receding voids. We determine that the origin of this bimodality resides in the void large scale environment, since once voids are classified into void-in-void (R-type) or void-in-cloud (S-type), R-types are found mutually receding away, while S-types approach each other. The magnitude of these systematic relative velocities account for more than 100 km/s, reaching large coherence lengths of up to 200 h$^{-1}$ Mpc . We have used samples of voids from the Sloan Digital Sky Survey Data Release 7 (SDSS-DR7) and the peculiar velocity field inferred from linear theory, finding fully consistent results with the simulation predictions. Thus, their relative motion suggests a scenario of a sparkling Universe, with approaching and receding voids according to their local environment., Comment: 6 pages, 4 Figures, accepted for publication in MNRAS Letters

We analyse the clustering of cosmic voids using a numerical simulation and the main galaxy sample from the Sloan Digital Sky Survey. We take into account the classification of voids into two types that resemble different evolutionary modes: those with a rising integrated density profile (void-in-void mode, or R-type) and voids with shells (void-in-cloud mode, or S-type). The results show that voids of the same type have stronger clustering than the full sample. We use the correlation analysis to define void clumps, associations with at least two voids separated by a distance of at most the mean void separation. In order to study the spatial configuration of void clumps, we compute the minimal spanning tree and analyse their multiplicity, maximum length and elongation parameter. We further study the dynamics of the smaller sphere that encloses all the voids in each clump. Although the global densities of void clumps are different according to their member-void types, the bulk motions of these spheres are remarkably lower than those of randomly placed spheres with the same radii distribution. In addition, the coherence of pairwise void motions does not strongly depend on whether voids belong to the same clump. Void clumps are useful to analyse the large-scale flows around voids, since voids embedded in large underdense regions are mostly in the void-in-void regime, were the expansion of the larger region produces the separation of voids. Similarly, voids around overdense regions form clumps that are in collapse, as reflected in the relative velocities of voids that are mostly approaching., Comment: Accepted for publication on MNRAS, 9 pages, 5 figures

The expanding complex pattern of filaments, walls and voids build the evolving cosmic web with material flowing from underdense onto high density regions. Here we explore the dynamical behaviour of voids and galaxies in void shells relative to neighboring overdense superstructures, using the Millenium Simulation and the main galaxy catalogue in Sloan Digital Sky Survey data. We define a correlation measure to estimate the tendency of voids to be located at a given distance from a superstructure. We find voids-in-clouds (S-types) preferentially located closer to superstructures than voids-in-voids (R-types) although we obtain that voids within $\sim40~\mathrm{Mpc}\,\mathrm{h}^{-1}$ of superstructures are infalling in a similar fashion independently of void type. Galaxies residing in void shells show infall towards the closest superstructure, along with the void global motion, with a differential velocity component depending on their relative position in the shell with respect to the direction to the superstructure. This effect is produced by void expansion and therefore is stronger for R-types. We also find that galaxies in void shells facing the superstrucure flow towards the overdensities faster than galaxies elsewere at the same relative distance to the superstructure. The results obtained for the simulation are also reproduced for the SDSS data with a linearized velocity field implementation., 10 pages, 9 figures. Accepted for publication in MNRAS

The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs., La lista completa de autores se encuentra en el documento., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

We test and present the application of the full rescaling method by Angulo & White (2010) to change the cosmology of halo catalogues in numerical simulations for cosmological parameter search using semi-analytic galaxy properties. We show that a reduced form of the method can be applied in small simulations with box side of ~50/h Mpc. We perform statistical tests on the accuracy of the properties of rescaled individual haloes, and also on the rescaled population as a whole. We find that individual positions and velocities are recovered with almost no detectable biases. The dispersion in the recovered halo mass does not seem to depend on the resolution of the simulation. Regardless of the halo mass, the individual accretion histories, spin parameter evolution and fraction of mass in substructures are well recovered. The mass of rescaled haloes can be underestimated (overestimated) for negative (positive) variations of either sigma_8 or Omega_m, in a way that does not depend on the halo mass. Statistics of abundances and correlation functions of haloes show also small biases of

We perform a statistical study of the global motion of cosmic voids using both a numerical simulation and observational data. We analyse their relation to large-scale mass flows and the physical effects that drive those motions. We analyse the bulk motions of voids, defined by the mean velocity of haloes in the surrounding shells in the numerical simulation, and by galaxies in the Sloan Digital Sky Survey Data Release 7. We find void mean bulk velocities close to 400 km/s, comparable to those of haloes (~ 500-600 km/s), depending on void size and the large-scale environment. Statistically, small voids move faster than large ones, and voids in relatively higher density environments have higher bulk velocities than those placed in large underdense regions. Also, we analyze the mean mass density around voids finding, as expected, large-scale overdensities (underdensities) along (opposite to) the void motion direction, suggesting that void motions respond to a pull-push mechanism. This contrasts with massive cluster motions who are mainly governed by the pull of the large-scale overdense regions. Our analysis of void pairwise velocities shows how their relative motions are generated by large-scale density fluctuations. In agreement with linear theory, voids embedded in low (high) density regions mutually recede (attract) each other, providing the general mechanism to understand the bimodal behavior of void motions. In order to compare the theoretical results and the observations we have inferred void motions in the SDSS using linear theory, finding that the estimated observational void motions are in statisticalagreement with the results of the simulation. Regarding large-scale flows, our results suggest a scenario of galaxies and galaxy systems flowing away from void centers with the additional, and morerelevant, contribution of the void bulk motion to the total velocity., 10 pages, 10 figures, Accepted for publication in MNRAS

We present a new technique for modelling submillimetre galaxies (SMGs): the 'Count Matching' approach. Using light cones drawn from a semi-analytic model of galaxy formation, we choose physical galaxy properties given by the model as proxies for their submillimetre luminosities, assuming a monotonic relationship. As recent interferometric observations of the Extended Chandra Deep Field-South show that the brightest sources detected by single-dish telescopes are comprised by emission from multiple fainter sources, we assign the submillimetre fluxes so that the combined Large APEX BOlometer CAmera (LABOCA) plus bright-end Atacama Large Millimeter/submillimeter Array observed number counts for this field are reproduced. After turning the model catalogues given by the proxies into submillimetre maps, we perform a source extraction to include the effects of the observational process on the recovered counts and galaxy properties. We find that for all proxies, there are lines of sight giving counts consistent with those derived from LABOCA observations, even for input sources with randomized positions in the simulated map. Comparing the recovered redshift, stellar mass and host halo mass distributions for model SMGs with observational data, we find that the best among the proposed proxies is that in which the submillimetre luminosity increases monotonically with the product between dust mass and star formation rate (SFR). This proxy naturally reproduces a positive trend between SFR and bolometric IR luminosity. The majority of components of blended sources are spatially unassociated., Instituto de Astrofísica de La Plata

Inspired on the well known dynamical dichotomy predicted in voids, where some underdense regions expand whereas others collapse due to overdense surrounding regions, we explored the interplay between the void inner dynamics and its large scale environment. The environment is classified depending on its density as in previous works. We analyse the dynamical properties of void-centered spherical shells at different void-centric distances depending on this classification. The above dynamical properties are given by the angular distribution of the radial velocity field, its smoothness, the field dependence on the tracer density and shape, and the field departures from linear theory. We found that the velocity field in expanding voids follows more closely the linear prediction, with a more smooth velocity field. However when using velocity tracers with large densities such deviations increase. Voids with sizes around $18\,h^{-1}\,Mpc$ are in a transition regime between regions with expansion overpredicted and underpredicted from linear theory. We also found that velocity smoothness increases as the void radius, indicating the laminar flow dominates the expansion of larger voids (more than $18\,h^{-1}\,Mpc$). The correlations observed suggest that nonlinear dynamics of the inner regions of voids could be dependent on the evolution of the surrounding structures. These also indicate possible scale couplings between the void inner expansion and the large scale regions where voids are embedded. These results shed some light to the origin of nonlinearities in voids, going beyond the fact that voids just quickly becomes nonlinear as they become emptier., Comment: 12 pages, 8 Figures, accepted for publication in MNRAS

We present a fast and accurate method to select an optimal set of parameters in semi-analytic models of galaxy formation and evolution (SAMs). Our approach compares the results of a model against a set of observables applying a stochastic technique called Particle Swarm Optimization (PSO), a self-learning algorithm for localizing regions of maximum likelihood in multidimensional spaces that outperforms traditional sampling methods in terms of computational cost. We apply the PSO technique to the SAG semi-analytic model combined with merger trees extracted from a standard Lambda Cold Dark Matter N-body simulation. The calibration is performed using a combination of observed galaxy properties as constraints, including the local stellar mass function and the black hole to bulge mass relation. We test the ability of the PSO algorithm to find the best set of free parameters of the model by comparing the results with those obtained using a MCMC exploration. Both methods find the same maximum likelihood region, however, the PSO method requires one order of magnitude fewer evaluations. This new approach allows a fast estimation of the best-fitting parameter set in multidimensional spaces, providing a practical tool to test the consequences of including other astrophysical processes in SAMs., Facultad de Ciencias Astronómicas y Geofísicas, Instituto de Astrofísica de La Plata

The angular momentum of galactic discs in semi-analytic models of galaxy formation is usually updated in time as material is accreted to the disc by adopting a constant dimensionless spin parameter and little attention is paid to the effects of accretion with misaligned angular momenta. These effects are the subject of this paper, where we adopt a Monte Carlo simulation for the changes in the direction of the angular momentum of a galaxy disc as it accretes matter based on accurate measurements from dark-matter haloes in the Millennium II simulation. In our semi-analytic model implementation, the flips seen in the dark-matter haloes are assumed to be the same for the cold baryons; however, we also assume that in the latter the flip also entails a difficulty for the disc to increase its angular momentum which causes the disc to become smaller relative to a no-flip case. This makes star formation to occur faster, especially in low-mass galaxies at all redshifts allowing galaxies to reach higher stellar masses faster. We adopt a new condition for the triggering of starbursts during mergers. As these produce the largest flips it is natural to adopt the disc instability criterion to evaluate the triggering of bursts in mergers instead of one based on mass ratios as in the original model. The new implementation reduces the average lifetimes of discs by a factor of ∼2, while still allowing old ages for the present-day discs of large spiral galaxies. It also provides a faster decline of star formation in massive galaxies and a better fit to the bright end of the luminosity function at z = 0., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

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., 15 pages, 11 figures, 2 tables. 1 figure removed to reduce paper extension, shorter discussion, accepted for publication in MNRAS

Galaxy surveys targeting emission lines are characterizing the evolution of star-forming galaxies, but there is still little theoretical progress in modelling their physical properties. We predict nebular emission from star-forming galaxies within a cosmological galaxy formation model. Emission lines are computed by combining the semi-analytical model SAG with the photoionization code MAPPINGS-III. We characterize the interstellar medium of galaxies by relating the ionization parameter of gas in galaxies to their cold gas metallicity, obtaining a reasonable agreement with the observed Hα, [O II] λ 3727, [O III] λ 5007 luminosity functions, and the BPT diagram for local star-forming galaxies. The average ionization parameter is found to increase towards low star formation rates and high redshifts, consistent with recent observational results. The predicted link between different emission lines and their associated star formation rates is studied by presenting scaling relations to relate them. Our model predicts that emission-line galaxies have modest clustering bias, and thus reside in dark matter haloes of masses below Mhalo ≲ 1012 [h−1 M⊙]. Finally, we exploit our modelling technique to predict galaxy number counts up to z ∼ 10 by targeting far-infrared emission lines detectable with submillimetre facilities., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

We make use of a semi-analytical model of galaxy formation to investigate the origin of the observed correlation between [α/Fe] abundance ratios and stellar mass in elliptical galaxies.We implement a new galaxy-wide stellar initial mass function (IMF; top-heavy integrated galaxy initial mass function, TH-IGIMF) in the semi-analytic model SAG (acronym for semi-analytic galaxies) and evaluate its impact on the chemical evolution of galaxies. The star formation rate (SFR) dependence of the slope of the TH-IGIMF is found to be key to reproducing the correct [α/Fe]-stellar mass relation. Massive galaxies reach higher [α/Fe] abundance ratios because they are characterized by more top-heavy IMFs as a result of their higher SFR. As a consequence of our analysis, the value of the minimum embedded star cluster mass and of the slope of the embedded cluster mass function, which are free parameters involved in the TH-IGIMF theory, are found to be as lowas 5 and 2M⊙, respectively. Amild downsizing trend is present for galaxies generated assuming either a universal IMF or a variable TH-IGIMF.We find that, regardless of galaxy mass, older galaxies (with formation redshifts ≳2) are formed in shorter time-scales (≲2 Gyr), thus achieving larger [α/Fe] values. Hence, the time-scale of galaxy formation alone cannot explain the slope of the [α/Fe]-galaxy mass relation, but is responsible for the big dispersion of [α/Fe] abundance ratios at fixed stellar mass. We further test the hypothesis of a TH-IGIMF in elliptical galaxies by looking into mass-to-light ratios, and luminosity functions. Models with a TH-IGIMF are also favoured by these constraints. In particular, mass-to-light ratios agree with observed values for massive galaxies while being overpredicted for less massive ones; this overprediction is present regardless of the IMF considered., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

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., Comment: 28 pages containing 13 figures & 4 tables + 9 pages appendix containing another 4 tables + 4 pages of references, accepted for publication in MNRAS

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., 10 pages, 14 figures and 3 tables. Accepted for publication in MNRAS