• ### We are not the 99 percent: quantifying asphericity in the distribution of Local Group satellites(1805.03188)

May 8, 2018 astro-ph.GA
We use simulations to build an analytic probability distribution for the asphericity in the satellite distribution around Local Group (LG) type galaxies in the Lambda Cold Dark Matter (LCDM) paradigm. We use this distribution to estimate the atypicality of the satellite distributions in the LG even when the underlying simulations do not have enough systems fully resembling the LG in terms of its typical masses, separation and kinematics. We demonstrate the method using three different simulations (Illustris-1, Illustris-1-Dark and ELVIS) and a number of satellites ranging from 11 to 15. Detailed results differ greatly among the simulations suggesting a strong influence of the typical DM halo mass, the number of satellites and the simulated baryonic effects. However, there are three common trends. First, at most $2\%$ of the pairs are expected to have satellite distributions with the same asphericity as the LG; second, at most $80\%$ of the pairs have a halo with a satellite distribution as aspherical as in M31; and third, at most $4\%$ of the pairs have a halo with satellite distribution as planar as in the MW. These quantitative results place the LG at the level of a $3\sigma$ outlier in the LCDM paradigm. We suggest that understanding the reasons for this atypicality requires quantifying the asphericity probability distribution as a function of halo mass and large scale environment. The approach presented here can facilitate that kind of study and other comparisons between different numerical setups and choices to study satellites around LG pairs in simulations.
• ### Unbiased clustering estimates with the DESI fibre assignment(1805.00951)

May 2, 2018 astro-ph.CO
The Emission Line Galaxy survey made by the Dark Energy Spectroscopic Instrument (DESI) survey will be created from five passes of the instrument on the sky. On each pass, the constrained mobility of the ends of the fibres in the DESI focal plane means that the angular-distribution of targets that can be observed is limited. Thus, the clustering of samples constructed using a limited number of passes will be strongly affected by missing targets. In two recent papers, we showed how the effect of missing galaxies can be corrected when calculating the correlation function using a weighting scheme for pairs. Using mock galaxy catalogues we now show that this method provides an unbiased estimator of the true correlation function for the DESI survey after any number of passes. We use multiple mocks to determine the expected errors given one to four passes, compared to an idealised survey observing an equivalent number of randomly selected targets. On BAO scales, we find that the error is a factor 2 worse after one pass, but that after three or more passes, the errors are very similar. Thus we find that the fibre assignment strategy enforced by the design of DESI will not affect the cosmological measurements to be made by the survey, and can be removed as a potential risk for this experiment.
• ### Modelling the gas kinematics of an atypical Lyman-alpha emitting compact dwarf galaxy(1710.05534)

Oct. 16, 2017 astro-ph.GA
Star-forming Compact Dwarf Galaxies (CDGs) resemble the expected pristine conditions of the first galaxies in the Universe and are the best systems to test models on primordial galaxy formation and evolution. Here we report on one of such CDGs, Tololo 1214-277, which presents a broad, single peaked, highly symmetric Ly$\alpha$ emission line that had evaded theoretical interpretation so far. In this paper we reproduce for the first time these line features with two different physically motivated kinematic models: an interstellar medium composed by outflowing clumps with random motions and an homogeneous gaseous sphere undergoing solid body rotation. The multiphase model requires a clump velocity dispersion of $54.3\pm 0.6$ km s$^{-1}$ with outflows of $54.3\pm 5.1$ km s$^{-1}$, while the bulk rotation velocity is constrained to be $348^{+75}_{-48}$ km s$^{-1}$. We argue that the results from the multiphase model provide a correct interpretation of the data. In that case the clump velocity dispersion implies a dynamical mass of $2\times 10^{9}$ M$_{\odot}$, ten times its baryonic mass. If future kinematic maps of Tololo 1214-277 confirm the velocities suggested by the multiphase model, it would provide additional support to expect such kinematic state in primordial galaxies, opening the opportunity to use the models and methods presented in this paper to constrain the physics of star formation and feedback in the early generation of Ly-$\alpha$ emitting galaxies.
• ### Tracing the cosmic web(1705.03021)

May 8, 2017 astro-ph.CO
The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.
• DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.
• DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.
• ### Boosting Lya and HeII 1640A Line Fluxes from Pop III Galaxies: Stochastic IMF Sampling and Departures from Case-B(1609.02150)

Sept. 7, 2016 astro-ph.GA
We revisit calculations of nebular hydrogen Lya and HeII 1640A line strengths for population III galaxies, undergoing continuous and bursts of star formation. We focus on initial mass functions (IMFs) motivated by recent theoretical studies, which generally span a lower range of stellar masses than earlier works. We also account for case-B departures and the stochastic sampling of the IMF. In agreement with previous works, we find that departures from case-B can enhance the Lya flux by a factor of a few, but we argue that this enhancement is driven mainly by collisional excitation and ionization, and not due to photoionization from the n = 2 state of atomic hydrogen. The increased sensitivity of the Lya flux to the high-energy end of the galaxy spectrum makes it more subject to stochastic sampling of the IMF. The latter introduces a dispersion in the predicted nebular line fluxes around the deterministic value by as much as a factor of ~4. In contrast, the stochastic sampling of the IMF has less impact on the emerging Lyman Werner (LW) photon flux. When case-B departures and stochasticity effects are combined, nebular line emission from population III galaxies can be up to one order of magnitude brighter than predicted by 'standard' calculations that do not include these effects. This enhances the prospects for detection with future facilities such as JWST and large, groundbased telescopes.
• ### Impact of Cosmic Variance on the Galaxy-Halo Connection for Lyman-$\alpha$ Emitters(1606.07443)

June 23, 2016 astro-ph.GA
• ### Target allocation yields for massively multiplexed spectroscopic surveys with fibers(1406.1787)

June 30, 2014 astro-ph.IM
We present Simulated Annealing fiber-to-target allocation simulations for the proposed DESI and 4MOST massively multiplexed spectroscopic surveys, and for both Poisson and realistically clustered mock target samples. We simulate both Echidna and theta-phi actuator designs, including the restrictions caused by the physical actuator characteristics during repositioning. For DESI, with theta-phi actuators, used in 5 passes over the sky for a mock ELG/LRG/QSO sample, with matched fiber and target densities, a total target allocation yield of 89.3% was achieved, but only 83.7% for the high-priority Ly-alpha QSOs. If Echidna actuators are used with the same pitch and number of passes, the yield increases by 5.7% and 16% respectively. Echidna also allows a factor-of-two increase in the number of close Ly-alpha QSO pairs that can be observed. Echidna spine tilt causes a variable loss of throughput, with average loss being the same as the loss at the rms tilt. With a natural tilt minimization scheme, we find an rms tilt always close to 0.58 x maximum. There is an additional but much smaller defocus loss, equivalent to an average defocus of 30microns. These tilt losses offset the gains in yield for Echidna, but because the survey strategy is driven by the higher priority targets, a clear survey speed advantage remains. For 4MOST, high and low latitude sample mock catalogs were supplied by the 4MOST team, and allocations were carried out with the proposed Echidna-based positioner geometry. At high latitudes, the resulting target completeness was 85.3% for LR targets and 78.9% for HR targets. At low latitude, the target completeness was 93.9% for LR targets and 71.2% for HR targets.
• ### The kinematics of the Local Group in a cosmological context(1303.2690)

March 11, 2013 astro-ph.CO
Recent observations constrained the tangential velocity of M31 with respect to the Milky Way (MW) to be v_tan<34.4 km/s and the radial velocity to be in the range v_rad=-109+/- 4.4 km/s (van der Marel et al. 2012). In this study we use a large volume high resolution N-body cosmological simulation (Bolshoi) together with three constrained simulations to statistically study this kinematics in the context of the LCDM. The comparison of the ensembles of simulated pairs with the observed LG at the 1-sigma level in the uncertainties has been done with respect to the radial and tangential velocities, the reduced orbital energy (e_tot), angular momentum (l_orb) and the dimensionless spin parameter, lambda. Our main results are: (i) the preferred radial and tangential velocities for pairs in LCDM are v_rad=-80+/-20 km/s, v_tan=50+/-10 km/s, (ii) pairs around that region are 3 to 13 times more common than pairs within the observational values, (iii) 15%to 24% of LG-like pairs in LCDM have energy and angular momentum consistent with observations while (iv) 9% to 13% of pairs in the same sample show similar values in the inferred dimensionless spin parameter. It follows that within current observational uncertainties the quasi-conserved quantities that characterize the orbit of the LG, i.e. e_tot, r_orb and lambda, do not challenge the standard LCDM model, but the model is in tension with regard to the actual values of the radial and tangential velocities. This might hint to a problem of the LCDM model to reproduce the observed LG.
• ### Effects of Star Formation Stochasticity on the Ly-alpha & Lyman Continuum Emission from Dwarf Galaxies during Reionization(1206.0726)

Oct. 18, 2012 astro-ph.CO, astro-ph.GA
Observations of distant galaxies play a key role in improving our understanding of the Epoch of Reionization (EoR). The observed Ly-alpha emission line strength - quantified by its restframe equivalent width (EW) - provides a valuable diagnostic of stellar populations and dust in galaxies during and after the EoR. In this paper we quantify the effects of star formation stochasticity on the predicted Ly-alpha EW in dwarf galaxies, using the publicly available code SLUG ('Stochastically Light Up Galaxies'). We compute the number of hydrogen ionizing photons, as well as flux in the Far UV for a set of models with star formation rates (SFR) in the range 10-3-1 Msol/yr. From these fluxes we compute the luminosity, L-alpha, and the EW of the Ly-alpha line. We find that stochasticity alone induces a broad distribution in L-alpha and EW at a fixed SFR, and that the widths of these distributions decrease with increasing SFR. We parameterize the EW probability density function (PDF) as an SFR-dependent double power law. We find that it is possible to have EW as low as ~EW0/4 and as high as ~3 times the EW0, where EW0 denotes the expected EW in the absence of stochasticity. We argue that stochasticity may therefore be important when linking drop-out and narrow-band selected galaxies, when identifying population III galaxies, and that it may help to explain the large EW (EW > 100 - 200 A) observed for a fraction of Ly-alpha- selected galaxies. Finally, we show that stochasticity can also affect the inferred escape fraction of ionizing photons from galaxies. In particular, we argue that stochasticity may simultaneously explain the observed anomalous ratios of the Lyman continuum flux density to the (non-ionizing) UV continuum density in so-called Lyman-Bump galaxies at z = 3.1, as well as the absence of such objects among a sample of z = 1.3 drop-out galaxies.
• ### Towards a panchromatic picture of galaxy evolution during the reionization epoch(1208.0819)

Aug. 3, 2012 astro-ph.CO
There are thousands of confirmed detections of star forming galaxies at high redshift (z > 4). These observations rely primarily on the detection of the spectral Lyman Break and the Lyman-alpha emission line. Theoretical modelling of these sources helps to interpret the observations in the framework of the standard cosmological paradigm. We present results from the High-z MareNostrum Project, aimed at constructing a panchromatic picture of the high redshift galaxy evolution that will improve our understanding of young star forming galaxies. Our simulation successfully reproduces the observational constraints from Lyman Break Galaxies and Lyman-alpha emitters at 5 < z < 7 . Based on this model we make predictions on the expected Far Infrared (FIR) emission that should be observed for LAEs. These predictions will help to settle down the question on the dust content of massive high-z galaxies, an issue that will be feasible to probe observationally with the Atacama Large Millimetre Array (ALMA).
• ### A kinematic classification of the cosmic web(1201.3367)

July 16, 2012 astro-ph.CO
A new approach for the classification of the cosmic web is presented. In extension of the previous work of Hahn et al. (2007) and Forero-Romero et al. (2009) the new algorithm is based on the analysis of the velocity shear tensor rather than the gravitational tidal tensor. The procedure consists of the construction of the the shear tensor at each (grid) point in space and the evaluation of its three eigenvectors. A given point is classified to be either a void, sheet, filament or a knot according to the number of eigenvalues above a certain threshold, 0, 1, 2, or 3 respectively. The threshold is treated as a free parameter that defines the web. The algorithm has been applied to a dark matter only, high resolution simulation of a box of side-length 64$h^{-1}$Mpc and N = $1024^3$ particles with the framework of the WMAP5/LCDM model. The resulting velocity based cosmic web resolves structures down to <0.1$h^{-1}$Mpc scales, as opposed to the ~1$h^{-1}$Mpc scale of the tidal based web. The under-dense regions are made of extended voids bisected by planar sheets, whose density is also below the mean. The over-dense regions are vastly dominated by the linear filaments and knots. The resolution achieved by the velocity based cosmic web provides a platform for studying the formation of halos and galaxies within the framework of the cosmic web.
• ### Visualising Matter and Cosmologies: A Transhistorical Example(1110.6772)

We propose a connection between the visualisation of cosmic matter and structure formation in the Cartesian tradition and that used by contemporary astrophysics. More precisely, we identify cosmological simulations of large scale structure in the Universe with the system of vortices in Descartes physics. This connection operates at different levels of the images: their representational purpose; the theoretical systems behind their use; and, finally, their function and materiality as visual productions. A skilled use of image analysis is necessary to stress the continuities and peculiarities between different epochs and disciplines.
• ### Predictability in Semi-Analytic Models of Galaxy Formation(0801.1953)

Jan. 13, 2008 astro-ph
We propose a general framework to scrutinize the performance of semi-analytic codes of galaxy formation. The approach is based on the analysis of the outputs from the model after a series of perturbations in the input parameters controlling the baryonic physics. The perturbations are chosen in a way that they do not change the results in the luminosity function or mass function of the galaxy population. We apply this approach on a particular semi-analytic model called GalICS. We chose to perturb the parameters controlling the efficiency of star formation and the efficiency of supernova feedback. We keep track of the baryonic and observable properties of the central galaxies in a sample of dark matter halos with masses ranging from 10^{10} M_sol to 10^{13} M_sol. We find very different responses depending on the halo mass. For small dark matter halos its central galaxy responds in a highly predictable way to small perturbation in the star formation and feedback efficiency. For massive dark matter halos, minor perturbations in the input parameters can induce large fluctuations on the properties of its central galaxy, at least $\sim 0.1$ in (B-V) color or $\sim 0.5$ mag in U or r filter, in a seemingly random fashion. We quantify this behavior through an objective scalar function we call predictability. We argue that finding the origin of this behavior needs additional information from other approximations and different semi-analytic codes. Furthermore, the implementation of an scalar objective function, such as the predictability, opens the door to quantitative benchmarking of semi-analytic codes based on its numerical performance.
• ### The Coarse Geometry of Merger Trees in \Lambda CDM(0801.1955)

Jan. 13, 2008 astro-ph
We introduce the contour process to describe the geometrical properties of merger trees. The contour process produces a one-dimensional object, the contour walk, which is a translation of the merger tree. We portray the contour walk through its length and action. The length is proportional to to the number of progenitors in the tree, and the action can be interpreted as a proxy of the mean length of a branch in a merger tree. We obtain the contour walk for merger trees extracted from the public database of the Millennium Run and also for merger trees constructed with a public Monte-Carlo code which implements a Markovian algorithm. The trees correspond to halos of final masses between 10^{11} h^{-1} M_sol and 10^{14} h^{-1} M_sol. We study how the length and action of the walks evolve with the mass of the final halo. In all the cases, except for the action measured from Markovian trees, we find a transitional scale around 3 \times 10^{12} h^{-1} M_sol. As a general trend the length and action measured from the Markovian trees show a large scatter in comparison with the case of the Millennium Run trees.
• ### LeMoMaF: Lensed Mock Map Facility(astro-ph/0611932)

Nov. 30, 2006 astro-ph
We present the Lensed Mock Map Facility (LeMoMaF), a tool designed to perform mock weak lensing measurements on numerically simulated chunks of the universe. Coupling N-body simulations to a semi-analytical model of galaxy formation, LeMoMaF can create realistic lensed images and mock catalogues of galaxies, at wavelengths ranging from the UV to the submm. To demonstrate the power of such a tool we compute predictions of the source-lens clustering effect on the convergence statistics, and quantify the impact of weak lensing on galaxy counts in two different filters. We find that the source-lens clustering effect skews the probability density function of the convergence towards low values, with an intensity which strongly depends on the redshift distribution of galaxies. On the other hand, the degree of enhancement or depletion in galaxy counts due to weak lensing is independent of the source-lens clustering effect. We discuss the impact on the two-points shear statistics to be measured by future missions like SNAP and LSST. The source-lens clustering effect would bias the estimation of sigma_8 from two point statistics by 2% -5%. We conclude that accurate photometric redshifts for individual galaxies are necessary in order to quantify and isolate the source-lens clustering effect.