Subhalo accretion through filaments

We track subhalo orbits of galaxy and group sized halos in cosmological simulations. We identify filamentary structures around halos and we use these to define a sample of subhalos accreted from filaments as well as a control sample of subhalos accreted from other directions. We use these samples to study differences in satellite orbits produced by filamentary accretion. Our results depend on host halo mass. We find that for low masses, subhalos accreted from filaments show $\sim10\%$ shorter lifetimes compared to the control sample, they show a tendency towards more radial orbits, reach halo central regions earlier, and are more likely to merge with the host. For higher mass halos this lifetime difference dissipates and even reverses for cluster sized halos. This behavior appears to be connected to the fact that more massive hosts are connected to stronger filaments with higher velocity coherence and density, with slightly more radial subhalo orbits. Because subhalos tend to follow the coherent flow of the filament, it is possible that such thick filaments are enough to shield the subhalo from the effect of dynamical friction at least during their first infall. We also identify subhalo pairs/clumps which merge with one another after accretion. They survive as a clump for only a very short time, which is even shorter for higher subhalo masses, suggesting that the Magellanic Clouds and other Local group satellite associations, may have entered the MW virial radius very recently and probably are in their first infall. Filaments boost the accretion of satellite associations.Comment: Accepted for publication in Ap

Dark Matter annihilation energy output and its effects on the high-z IGM

We study the case of DM self annihilation, in order to assess its importance as an energy injection mechanism, to the IGM in general, and to the medium within particular DM haloes. We consider thermal relic WIMP particles with masses of 10GeV and 1TeV and we analyse in detail the clustering properties of DM in a $\Lambda$CDM cosmology, on all hierarchy levels, from haloes and their mass function, to subhaloes and the DM density profiles within them, considering adiabatic contraction by the presence of a SMBH. We then compute the corresponding energy output, concluding that DM annihilation does not constitute an important feedback mechanism. We also calculate the effects that DM annihilation has on the IGM temperature and ionization fraction, and we find that assuming maximal energy absorption, at z ~ 10, for the case of a 1TeV WIMP, the ionization fraction could be raised to $6 \times 10^{-4}$ and the temperature to 10K, and in the case of a 10GeV WIMP, the IGM temperature could be raised to 200K and the ionization fraction to $8 \times 10^{-3}$. We conclude that DM annihilations cannot be regarded as an alternative reionization scenario. Regarding the detectability of the WIMP through the modifications to the 21 cm differential brightness temperature signal ($\delta$Tb), we conclude that a thermal relic WIMP with mass of 1TeV is not likely to be detected from the global signal alone, except perhaps at the 1-3mK level in the frequency range 30MHz < $\nu$ < 35MHz corresponding to 40 < z < 46. However, a 10GeV mass WIMP may be detectable at the 1-3mK level in the frequency range 55MHz < $\nu$ < 119MHz corresponding to 11 < z < 25, and at the 1-10mK level in the frequency range 30MHz < $\nu$ < 40MHz corresponding to 35 < z < 46.Comment: 23 pages, 12 figures, accepted for publication in MNRA

The Intrinsic Shape of Galaxies in SDSS/Galaxy Zoo

By modelling the axis ratio distribution of SDSS DR8 galaxies we find the intrinsic 3D shapes of spirals and ellipticals. We use morphological information from the Galaxy Zoo project and assume a non-parametric distribution intrinsic of shapes, while taking into account dust extinction. We measure the dust extinction of the full sample of spiral galaxies and find a smaller value than previous estimations, with an edge-on extinction of $E_0 = 0.284^{+0.015}_{-0.026}$ in the SDSS r band. We also find that the distribution of minor to major axis ratio has a mean value of $0.267 \pm 0.009$, slightly larger than previous estimates mainly due to the lower extinction used; the same affects the circularity of galactic discs, which are found to be less round in shape than in previous studies, with a mean ellipticity of $0.215 \pm 0.013$. For elliptical galaxies, we find that the minor to major axis ratio, with a mean value of $0.584 \pm 0.006$, is larger than previous estimations due to the removal of spiral interlopers present in samples with morphological information from photometric profiles. These interlopers are removed when selecting ellipticals using Galaxy Zoo data. We find that the intrinsic shapes of galaxies and their dust extinction vary with absolute magnitude, colour and physical size. We find that bright elliptical galaxies are more spherical than faint ones, a trend that is also present with galaxy size, and that there is no dependence of elliptical galaxy shape with colour. For spiral galaxies we find that the reddest ones have higher dust extinction as expected, due to the fact that this reddening is mainly due to dust. We also find that the thickness of discs increases with luminosity and size, and that brighter, smaller and redder galaxies have less round discs.Comment: 15 pages, 16 figures, submitted to MNRA

Future virialized structures: An analysis of superstructures in SDSS-DR7

We construct catalogues of present superstructures that, according to a LCDM scenario, will evolve into isolated, virialized structures in the future. We use a smoothed luminosity density map derived from galaxies in SDSS-DR7 data and separate high luminosity density peaks. The luminosity density map is obtained from a volume-limited sample of galaxies in the spectroscopic galaxy catalogue, within the SDSS-DR7 footprint area and in the redshift range 0.04 < z < 0.12. Other two samples are constructed for calibration and testing purposes, up to z = 0.10 and z = 0.15. The luminosity of each galaxy is spread using an Epanechnikov kernel of 8Mpc/h radius, and the map is constructed on a 1 Mpc/h cubic cells grid. Future virialized structures (FVS) are identified as regions with overdensity above a given threshold, calibrated using a LCDM numerical simulation, and the criteria presented by D\"unner et al. (2006). We assume a constant mass-to-luminosity ratio and impose the further condition of a minimum luminosity of 10^{12}Lsol. According to our calibrations with a numerical simulation, these criteria lead to a negligible contamination by less overdense (non FVS) superstructures.We present a catalogue of superstructures in the SDSS-DR7 area within redshift 0.04 < z < 0.12 and test the reliability of our method by studying different subsamples as well as a mock catalogue.We compute the luminosity and volume distributions of the superstructures finding that about 10% of the luminosity (mass) will end up in future virialized structures. The fraction of groups and X-ray clusters in these superstructures is higher for groups/clusters of higher mass, suggesting that future cluster mergers will involve the most massive systems. We also analyse known structures in the present Universe and compare with our catalogue of FVS.Comment: 14 pages, 11 figures, modified to match accepted version in MNRAS. PDF with high resolution colour figures is available at http://www.oac.uncor.edu/apache2-default/adminweb/html/WEB/preprints/2011.01/FVS-DR7.pd

Galaxy groups in the 2dF galaxy redshift survey: Large Scale Structure with Groups

We use the 2dF Galaxy Group Catalogue constructed by Merch\'an & Zandivarez to study the large scale structure of the Universe traced by galaxy groups. We concentrate on the computation of the power spectrum and the two point correlation function. The resulting group power spectrum shows a similar shape to the galaxy power spectrum obtained from the 2dF Galaxy Redshift Survey by Percival et al., but with a higher amplitude quantified by a relative bias in redshift space of $b_s(k) \sim 1.5$ . The group two point correlation function for the total sample is well described by a power law with correlation length s_0=8.9 \pm 0.3 \mpc and slope $\gamma=-1.6 \pm 0.1$ on scales s < 20 \mpc. In order to study the dependence of the clustering properties on group mass we split the catalogue in four subsamples defined by different ranges of group virial masses. These computations allow a fair estimate of the relation described by the correlation length $s_0$ and the mean intergroup separation $d_c$ for galaxy systems of low mass. We also extend our study to the redshift space distortions of galaxy groups, where we find that the anisotropies in the clustering pattern of the 2dF group catalogue are consistent with gravitational instability, with a flattening of the redshift-space correlation function contours in the direction of the line of sight.Comment: 11 pages, 9 figures, resubmitted to MNRAS after revisio

Angular momentum-Large-scale structure alignments in LCDM models and the SDSS

We study the alignments between the angular momentum of individual objects and the large-scale structure in cosmological numerical simulations and real data from the Sloan Digital Sky Survey, Data Release 6. To this end we measure anisotropies in the two point cross-correlation function around simulated halos and observed galaxies, studying separately the 1- and 2-halo regimes. The alignment of the angular momentum of dark-matter haloes in LCDM simulations is found to be dependent on scale and halo mass. At large distances (2-halo regime), the spins of high mass haloes are preferentially oriented in the direction perpendicular to the distribution of matter; lower mass systems show a weaker trend that may even reverse to show an angular momentum in the plane of the matter distribution. In the 1-halo term regime, the angular momentum is aligned in the direction perpendicular to the matter distribution; the effect is stronger than for the 1-halo term and increases for higher mass systems. On the observational side, we focus our study on galaxies in the Sloan Digital Sky Survey, Data Release 6 (SDSS-DR6) with elongated apparent shapes, and study alignments with respect to the major semi-axis. We find an excess of structure in the direction of the major semi-axis for all samples; the red sample shows the highest alignment (2.7+-0.08%) and indicates that the angular momentum of flattened spheroidals tends to be perpendicular to the large-scale structure. (Abridged)Comment: 10 pages, 6 figures, accepted for publication in MNRAS; the definitive version is available at www.blackwell-synergy.co

Calibration of semi-analytic models of galaxy formation using Particle Swarm Optimization

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$CDM 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 less 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.Comment: 11 pages, 4 figures, 1 table. Accepted for publication in ApJ. Comments are welcom

Automated detection of filaments in the large scale structure of the universe

We present a new method to identify large scale filaments and apply it to a cosmological simulation. Using positions of haloes above a given mass as node tracers, we look for filaments between them using the positions and masses of all the remaining dark-matter haloes. In order to detect a filament, the first step consists in the construction of a backbone linking two nodes, which is given by a skeleton-like path connecting the highest local dark matter (DM) density traced by non-node haloes. The filament quality is defined by a density and gap parameters characterising its skeleton, and filament members are selected by their binding energy in the plane perpendicular to the filament. This membership condition is associated to characteristic orbital times; however if one assumes a fixed orbital timescale for all the filaments, the resulting filament properties show only marginal changes, indicating that the use of dynamical information is not critical for the method. We test the method in the simulation using massive haloes($M>10^{14}$h$^{-1}M_{\odot}$) as filament nodes. The main properties of the resulting high-quality filaments (which corresponds to $\simeq33$ of the detected filaments) are, i) their lengths cover a wide range of values of up to $150$h$^{-1}$Mpc, but are mostly concentrated below 50h$^{-1}$Mpc; ii) their distribution of thickness peaks at $d=3.0$h$^{-1}$Mpc and increases slightly with the filament length; iii) their nodes are connected on average to $1.87\pm0.18$ filaments for $\simeq 10^{14.1}M_{\odot}$ nodes; this number increases with the node mass to $\simeq 2.49\pm0.28$ filaments for $\simeq 10^{14.9}M_{\odot}$ nodes.Comment: 17 pages, 13 figures, MNRAS Accepte

How accurate is it to update the cosmology of your halo catalogues?

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 <10 percent when moving away from the base simulation by up to 2 times the uncertainty in the WMAP7 cosmological parameters. The merger tree properties related to the final galaxy population in haloes also show small biases; the time since the last major merger, the assembly time-scale, and a time-scale related to the stellar ages show correlated biases which indicate that the spectral shapes of galaxies would only be affected by global age changes of ~150 Myr. We show some of these biases for different separations in the cosmological parameters with respect to the desired cosmology so that these can be used to estimate the expected accuracy of the resulting halo population. We also present a way to construct grids of simulations to provide stable accuracy across the Omega_m vs sigma_8 parameter space.Comment: 14 pages, 2 tables, 10 figures. Accepted for publication in MNRA