Search

Agent Based Models (ABMs) have emerged as a powerful tool for investigating complex social interactions, particularly in the context of public health and infectious disease investigation. In an effort to enhance the conventional ABM, enabling automated model calibration and reducing the computational resources needed for scaling up the model, we have developed a tensorized and differentiable agent-based model by coupling Graph Neural Network (GNN) and Long Short-Term Memory (LSTM) network. The model was employed to investigate the 2019 measles outbreak occurred in New Zealand, demonstrating a promising ability to accurately simulate the outbreak dynamics, particularly during the peak period of repeated cases. This paper shows that by leveraging the latest Artificial Intelligence (AI) technology and the capabilities of traditional ABMs, we gain deeper insights into the dynamics of infectious disease outbreaks. This, in turn, helps us make more informed decision when developing effective strategies that strike a balance between managing outbreaks and minimizing disruptions to everyday life., Comment: 11 pages, 9 figures

We present zELDA(redshift Estimator for Line profiles of Distant Lyman-Alpha emitters), an open source code to fit Lyman-Alpha (Lya) line profiles. The main motivation is to provide the community with an easy to use and fast tool to analyze Lya line profiles uniformly to improve the understating of Lya emitting galaxies. zELDA is based on line profiles of the commonly used 'shell-model' pre-computed with the full Monte Carlo radiative transfer code LyaRT. Via interpolation between these spectra and the addition of noise, we assemble a suite of realistic Lya spectra which we use to train a deep neural network. We show that the neural network can predict the model parameters to high accuracy (e.g.,.0.34 dex HI column density for R=12000) and thus allows for a significant speedup over existing fitting methods. As a proof of concept, we demonstrate the potential of zELDA by fitting 97 observed Lya line profiles from the LASD data base. Comparing the fitted value with the measured systemic redshift of these sources, we find that Lya determines their rest frame Lya wavelength with a remarkable good accuracy of 0.3A (75 km/s). Comparing the predicted outflow properties and the observed Lya luminosity and equivalent width, we find several possible trends. For example, we find an anticorrelation between the Lya luminosity and the outflow neutral hydrogen column density, which might be explained by the radiative transfer process within galaxies, Comment: 34 pages, 23 figures. zELDA's tutorial and installation guide: https://lya-zelda.readthedocs.io/en/latest/

The next generation of spectroscopic surveys will target emission-line galaxies (ELGs) to produce constraints on cosmological parameters. We study the large scale structure traced by ELGs using a combination of a semi-analytical model of galaxy formation, a code that computes the nebular emission from HII regions using the properties of the interstellar medium, and a large-volume, high-resolution N-body simulation. We consider fixed number density samples where galaxies are selected by either their H$\alpha$, [OIII]$\lambda 5007$ or [OII]$\lambda \lambda 3727-3729$ emission line luminosities. We investigate the assembly bias signatures of these samples, and compare them to those of stellar mass and SFR selected samples. Interestingly, we find that the [OIII]- and [OII]-selected samples display scale-dependent bias on large scales and that their assembly bias signatures are also scale-dependent. Both these effects are more pronounced for lower number density samples. The [OIII] and [OII] emitters that contribute most to the scale dependence tend to have a low gas-phase metallicity and are preferentially found in low-density regions. We also measure the baryon acoustic oscillation (BAO) feature and the $\beta$ parameter related to the growth rate of overdensities. We find a slight tendency for the BAO peak to shift toward smaller scales for [OII] emitters and that $\beta$ is scale-dependent at large scales. Our results suggest that ELG samples include environmental effects that should be modelled in order to remove potential systematic errors that could affect the estimation of cosmological parameters., Comment: 15 pages, 12 figures

We present the photometric determination of the bright-end (L_Lya>10^43.5 erg/s) of the Lya luminosity function (LF) within four redshifts windows in the interval 2.210^43.3 erg/s) Lya-emitting candidates. We test our selection with two spectroscopic follow-up programs at the GTC telescope, confirming ~89% of the targets as line-emitting sources, with ~64% being genuine z~2.2 QSOs. We extend the 2.22 AGN/QSOs (likely dominant in our samples) are described by a structurally different LF than z>2 star-forming LAEs, namely with L^*_QSOs ~ 100 L^*_LAEs and Phi^*_QSOs ~ 10^-3 Phi^*_LAEs. Finally, our method identifies very efficiently as high-z line-emitters sources without previous spectroscopic confirmation, currently classified as stars (~2000 objects in each redshift bin, on average). Assuming a large predominance of Lya-emitting AGN/QSOs in our samples, this supports the scenario by which these are the most abundant class of z>2 Lya emitters at L_Lya>10^43.3 erg/s., Comment: 27 pages (6 in appendix), 21 figures (5 in appendix), submitted to A&A

We explore how to mitigate the clustering distortions in Lyman-$\alpha$ emitters (LAEs) samples caused by the miss-identification of the Lyman-$\alpha$ (Ly$\alpha$) wavelength in their Ly$\alpha$ line profiles. We use the Ly$\alpha$ line profiles from our previous LAE theoretical model that includes radiative transfer in the interstellar and intergalactic mediums. We introduce a novel approach to measure the systemic redshift of LAEs from their Ly$\alpha$ line using neural networks. In detail, we assume that, for a fraction of the whole LAE population their systemic redshift is determined precisely through other spectral features. We then use this subset to train a neural network that predicts the Ly$\alpha$ wavelength given a Ly$\alpha$ line profile. We test two different training sets: i) the LAEs are selected homogeneously and ii) only the brightest LAEs are selected. In comparison with previous approaches in the literature, our methodology improves significantly both accuracy and precision in determining the Ly$\alpha$ wavelength. In fact, after applying our algorithm in ideal Ly$\alpha$ line profiles, we recover the clustering unperturbed down to 1cMpc/h. Then, we test the performance of our methodology in realistic Ly$\alpha$ line profiles by downgrading their quality. The machine learning techniques work well even if the Ly$\alpha$ line profile quality is decreased considerably. We conclude that LAE surveys such as HETDEX would benefit from determining with high accuracy the systemic redshift of a subpopulation and applying our methodology to estimate the systemic redshift of the rest of the galaxy sample., Comment: 24 Pages, 16 figures, a lot of fun

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1

We present a synthetic galaxy lightcone specially designed for narrow-band optical photometric surveys. To reduce time-discreteness effects, unlike previous works, we directly include the lightcone construction in the \texttt{L-Galaxies} semi-analytic model applied to the subhalo merger trees of the {\tt Millennium} simulation. Additionally, we add a model for the nebular emission in star-forming regions, which is crucial for correctly predicting the narrow/medium-band photometry of galaxies. Explicitly, we consider, individually for each galaxy, the contribution of 9 different lines: $\rm Ly{\alpha}$ (1216\AA), \Hb (4861\AA), \Ha (6563\AA), {\oii} (3727\AA, 3729\AA), {\oiii} (4959\AA, 5007\AA), $\rm [\ion{Ne}{III}]$ (3870\AA), {\oi} (6300\AA), $\rm [\ion{N}{II}]$ (6548\AA, 6583\AA), and $\rm [\ion{S}{II}]$ (6717\AA, 6731\AA). We validate our lightcone by comparing galaxy number counts, angular clustering, and \Ha, \Hb, {\oii} and {\oiiiFd} luminosity functions to a compilation of observations. As an application of our mock lightcones, we generate catalogues tailored for J-PLUS, a large optical galaxy survey featuring 5 broad and 7 medium band filters. We study the ability of the survey to correctly identify, with a simple \textit{three filter method}, a population of emission-line galaxies at various redshifts. We show that the $4000\AA$ break in the spectral energy distribution of galaxies can be misidentified as line emission. However, all significant excess (larger than 0.4 magnitudes) can be correctly and unambiguously attributed to emission line galaxies. Our catalogues are publicly released to facilitate their use in interpreting narrow-band surveys and for quantifying the impact of line emission in broad band photometry., Comment: Published in A&A; Mock catalogue available at https://www.j-plus.es/ancillarydata/mock\_galaxy\_lightcone

In the near future galaxy surveys will target Lyman alpha emitting galaxies (LAEs) to unveil the nature of the dark energy. It has been suggested that the observability of LAEs is coupled to the large scale properties of the intergalactic medium. Such coupling could introduce distortions into the observed clustering of LAEs, adding a new potential difficulty to the interpretation of upcoming surveys. We present a model of LAEs that incorporates Lyman-alpha radiative transfer processes in the interstellar and intergalactic medium. The model is implemented in the GALFORM semi-analytic model of galaxy of formation and evolution. We find that the radiative transfer inside galaxies produces selection effects over galaxy properties. In particular, observed LAEs tend to have low metallicities and intermediate star formation rates. At low redshift we find no evidence of a correlation between the spatial distribution of LAEs and the intergalactic medium properties. However, at high redshift the LAEs are linked to the line of sight velocity and density gradient of the intergalactic medium. The strength of the coupling depends on the outflow properties of the galaxies and redshift. This effect modifies the clustering of LAEs on large scales, adding non linear features. In particular, our model predicts modifications to the shape and position of the baryon acoustic oscillation peak. This work highlights the importance of including radiative transfer physics in the cosmological analysis of LAEs., Comment: 25 pages, 18 figures , a lot of love

Galaxies form and evolve in the context of their local and large-scale environments. Their baryonic content that we observe with imaging and spectroscopy is intimately connected to the properties of their dark matter halos, and to their location in the "cosmic web" of large-scale structure. Very large spectroscopic surveys of the local universe (e.g., SDSS and GAMA) measure galaxy positions (location within large-scale structure), statistical clustering (a direct constraint on dark matter halo masses), and spectral features (measuring physical conditions of the gas and stars within the galaxies, as well as internal velocities). Deep surveys with the James Webb Space Telescope (JWST) will revolutionize spectroscopic measurements of redshifts and spectral properties for galaxies out to the epoch of reionization, but with numerical statistics and over cosmic volumes that are too small to map large-scale structure and to constrain halo properties via clustering. Here, we consider advances in understanding galaxy evolution that would be enabled by very large spectroscopic surveys at high redshifts: very large numbers of galaxies (outstanding statistics) over large co-moving volumes (large-scale structure on all scales) over broad redshift ranges (evolution over most of cosmic history). The required observational facility can be established as part of the probe portfolio by NASA within the next decade., Comment: 8 pages (including cover page and references), 3 figures. Science white paper submitted to Astro2020. arXiv admin note: substantial text overlap with arXiv:1802.01539

The nature of dark energy remains a profound mystery 20 years after the discovery of cosmic acceleration. A very high number density galaxy redshift survey over a wide area (HD GRS Wide) spanning the redshift range of 0.5

We investigate the spatial distribution of Lyman alpha emitting galaxies (LAEs) at $z \approx 2.67$, selected from the NOAO Deep Wide-Field Survey (NDWFS), using two-point statistics and topological diagnostics adopted from network science. We measure the clustering length, $r_0 \approx 4 h^{-1}$ Mpc, and the bias, $b_{LAE} = 2.2^{+0.2}_{-0.1}$. Fitting the clustering with halo occupation distribution (HOD) models results in two disparate possibilities: (1) where the fraction of central galaxies is $<$1% in halos of mass $>10^{12}$$M_\odot$; and (2) where the fraction is $\approx$20%. We refer to these two scenarios as the `Dusty Core Scenario' for Model#1 since most of central galaxies in massive halos are dead in Ly$\alpha$ emission, and the `Pristine Core Scenario' for Model#2 since the central galaxies are bright in Ly$\alpha$ emission. Traditional two-point statistics cannot distinguish between these disparate models given the current data sets. To overcome this degeneracy, we generate mock catalogs for each HOD model using a high resolution $N$-body simulation and adopt a network statistics approach, which provides excellent topological diagnostics for galaxy point distributions. We find three topological anomalies from the spatial distribution of observed LAEs, which are not reproduced by the HOD mocks. We find that Model#2 matches better all network statistics than Model#1, suggesting that the central galaxies in $> 10^{12} h^{-1} M_\odot$ halos at $z \approx 2.67$ need to be less dusty to be bright as LAEs, potentially implying some replenishing channels of pristine gas such as the cold mode accretion., Comment: 23 pages, 18 figures, accepted by MNRAS

We present FLaREON (Fast Lyman-Alpha Radiative Escape from Outflowing Neutral gas), a public python package that delivers fast and accurate Lyman alpha escape fractions and line profiles over a wide range of outflow geometries and properties. The code incorporates different algorithms, such as interpolation and machine learning to predict Lyman alpha line properties from a pre-computed grid of outflow configurations based on the outputs of a Monte Carlo radiative transfer code. Here we describe the algorithm, discuss its performance and illustrate some of its many applications. Most notably, FLaREON can be used to infer the physical properties of the outflowing medium from an observed Lyman alpha line profile, including the escape fraction, or it can be run over millions of objects in a galaxy formation model to simulate the escape of Lyman alpha photons in a cosmological volume., Comment: 7 pages, 5 figures

Lyman-{\alpha} emitters (LAEs) are a promising target to probe the large scale structure of the Universe at high redshifts, $z\gtrsim 2$. However, their detection is sensitive to radiative transfer effects that depend on local astrophysical conditions. Thus, modeling the bulk properties of this galaxy population remains challenging for theoretical models. Here we develop a physically-motivated scheme to predict LAEs in cosmological simulations. The escape of Ly{\alpha} photons is computed using a Monte Carlo radiative transfer code which outputs a Ly{\alpha} escape fraction. To speed-up the process of assigning escape fractions to individual galaxies, we employ fitting formulae that approximate the full Monte Carlo results within an accuracy of 10% for a broad range of column densities, gas metallicities and gas bulk velocities. We apply our methodology to the semi-analytical model GALFORM on a large N-body simulation. The Ly{\alpha} photons escape through an outflowing neutral gas medium, implemented assuming different geometries. This results in different predictions for the typical column density and outflow velocities of the LAE population. To understand the impact of radiative transfer on our predictions, we contrast our models against a simple abundance matching assignment. Our full models populate LAEs in less massive haloes than what is obtained with abundance matching. Overall, radiative transfer effects result in better agreement when confronting the properties of LAEs against observational measurements. This suggest that incorporating the effects of Ly{\alpha} radiative transfer in the analysis of this galaxy population, including their clustering, can be important for obtaining an unbiased interpretation of future datasets., Comment: 25 pages, 18 figures. Submitted to MNRAS. Comments are welcome

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS Probe science spans four broad categories: (1) Revolutionizing galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Opening a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probing the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Exploring the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS Probe is a 1.5m telescope with a field of view of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000 over 1-4 microns, and a spectroscopic multiplex factor >5,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology. ATLAS Probe will lead to transformative science over the entire range of astrophysics: from galaxy evolution to the dark Universe, from Solar System objects to the dusty regions of the Milky Way., Comment: 46 pages. Accepted by PASA

We use the semi-analytic model of galaxy formation SAG to study the relevance of mass and environmental quenching on satellite galaxies. We find that environmental processes dominate the star formation (SF) quenching of low-mass satellites ($M_{\star} \lesssim 10^{10.5}\, {\rm M}_{\odot}$), whereas high-mass galaxies typically quench as centrals. High-mass galaxies that remain actively forming stars while being accreted are found to be mainly affected by mass quenching after their first infall. For a given stellar mass, our model predicts SF quenching to be less efficient in low-mass haloes both before and after infall, in contradiction with common interpretations of observational data. Our model supports a two-stage scenario to explain the SF quenching. Initially, the SF of satellites resembles that of centrals until the gas cooling rate is reduced to approximately half its value at infall. Then, the SF fades through secular processes that exhaust the cold gas reservoir. This reservoir is not replenished efficiently due to the action of either ram-pressure stripping (RPS) of the hot gas in low-mass satellites, or feedback from the active galactic nucleus (AGN) in high-mass satellites. The delay times for the onset of SF quenching are found to range from $\approx 3\,{\rm Gyr}$ to $\approx 1\,{\rm Gyr}$ for low-mass ($M_{\star} \approx 10^{10}\, {\rm M}_{\odot}$) and high-mass ($M_{\star} \approx 10^{11}\, {\rm M}_{\odot}$) satellites, respectively. SF fades in $\approx 1.5\,{\rm Gyr}$, largely independent of stellar mass. We find that the SF quenching of low-mass satellites supports the so-called delay-then-rapid quenching scenario. However, the SF history of $z=0$ passive satellites of any stellar mass is better described by a delay-then-fade quenching scenario., Comment: Moderate changes introduced, added figures, revised discussion on quenching scenario. Submitted to 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 \gtrsim 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_{\star}\gtrsim 10^{10}\,{\rm M}_{\odot}$)., Comment: New figure added (Fig. 5); section 5 modified accordingly. Published in MNRAS

We use the semi-analytical model of galaxy formation GALFORM to characterise an indirect signature of AGN feedback in the environment of radio galaxies at high redshifts. The predicted environment of radio galaxies is denser than that of radio-quiet galaxies with the same stellar mass. This is consistent with observational results from the CARLA survey. Our model shows that the differences in environment are due to radio galaxies being hosted by dark matter haloes that are ~1.5 dex more massive than those hosting radio-quiet galaxies with the same stellar mass. By running a control-simulation in which AGN feedback is switched-off, we identify AGN feedback as the primary mechanism affecting the build-up of the stellar component of radio galaxies, thus explaining the different environment in radio galaxies and their radio-quiet counterparts. The difference in host halo mass between radio loud and radio quiet galaxies translates into different galaxies populating each environment. We predict a higher fraction of passive galaxies around radio loud galaxies compared to their radio-quiet counterparts. Furthermore, such a high fraction of passive galaxies shapes the predicted infrared luminosity function in the environment of radio galaxies in a way that is consistent with observational findings. Our results suggest that the impact of AGN feedback at high redshifts and environmental mechanisms affecting galaxies in high halo masses can be revealed by studying the environment of radio galaxies, thus providing new constraints on galaxy formation physics at high redshifts., Comment: 12 pages, 7 Figures, submitted to MNRAS

Galaxy surveys aim to map the large-scale structure of the Universe and use redshift space distortions to constrain deviations from general relativity and probe the existence of massive neutrinos. However, the amount of information that can be extracted is limited by the accuracy of theoretical models used to analyze the data. Here, by using the L-Galaxies semi-analytical model run over the MXXL N-body simulation, we assess the impact of galaxy formation on satellite kinematics and the theoretical modelling of redshift-space distortions. We show that different galaxy selection criteria lead to noticeable differences in the radial distributions and velocity structure of satellite galaxies. Specifically, whereas samples of stellar mass selected galaxies feature satellites that roughly follow the dark matter, emission line satellite galaxies are located preferentially in the outskirts of halos and display net infall velocities. We demonstrate that capturing these differences is crucial for modelling the multipoles of the correlation function in redshift space, even on large scales. In particular, we show how modelling small scale velocities with a single Gaussian distribution leads to a poor description of the measure clustering. In contrast, we propose a parametrization that is flexible enough to model the satellite kinematics, and that leads to and accurate description of the correlation function down to sub-Mpc scales. We anticipate that our model will be a necessary ingredient in improved theoretical descriptions of redshift space distortions, which together could result in significantly tighter cosmological constraints and a more optimal exploitation of future large datasets., Comment: 16 pages, 15 figures, typos corrected and moderate revision after referee report

The Large Synoptic Survey Telescope (LSST) will survey the southern sky from 2022--2032 with unprecedented detail. Since the observing strategy can lead to artifacts in the data, we investigate the effects of telescope-pointing offsets (called dithers) on the $r$-band coadded 5$\sigma$ depth yielded after the 10-year survey. We analyze this survey depth for several geometric patterns of dithers (e.g., random, hexagonal lattice, spiral) with amplitude as large as the radius of the LSST field-of-view, implemented on different timescales (per season, per night, per visit). Our results illustrate that per night and per visit dither assignments are more effective than per season. Also, we find that some dither geometries (e.g., hexagonal lattice) are particularly sensitive to the timescale on which the dithers are implemented, while others like random dithers perform well on all timescales. We then model the propagation of depth variations to artificial fluctuations in galaxy counts, which are a systematic for large-scale structure studies. We calculate the bias in galaxy counts caused by the observing strategy, accounting for photometric calibration uncertainties, dust extinction, and magnitude cuts; uncertainties in this bias limit our ability to account for structure induced by the observing strategy. We find that after 10 years of the LSST survey, the best dither strategies lead to uncertainties in this bias smaller than the minimum statistical floor for a galaxy catalog as deep as $r$$<$27.5. A few of these strategies bring the uncertainties close to the statistical floor for $r$$<$25.7 after only one year of survey., Comment: 16 pages, 11 figures, 1 table. Updated to match the version accepted by ApJ

We use the GALFORM semi-analytical model to study high density regions traced by radio galaxies and quasars at high redshifts. We explore the impact that baryonic physics has upon the properties of galaxies in these environments. Star-forming emission-line galaxies (Ly{\alpha} and H{\alpha} emitters) are used to probe the environments at high redshifts. Radio galaxies are predicted to be hosted by more massive haloes than quasars, and this is imprinted on the amplitude of galaxy overdensities and cross-correlation functions. We find that Ly{\alpha} radiative transfer and AGN feedback indirectly affect the clustering on small scales and also the stellar masses, star- formation rates and gas metallicities of galaxies in dense environments. We also investigate the relation between protoclusters associated with radio galaxies and quasars, and their present- day cluster descendants. The progenitors of massive clusters associated with radio galaxies and quasars allow us to determine an average protocluster size in a simple way. Overdensities within the protoclusters are found to correlate with the halo descendant masses. We present scaling relations that can be applied to observational data. By computing projection effects due to the wavelength resolution of modern spectrographs and narrow-band filters we show that the former have enough spectral resolution to map the structure of protoclusters, whereas the latter can be used to measure the clustering around radio galaxies and quasars over larger scales to determine the mass of dark matter haloes hosting them., Comment: 13 pages, 9 figures,accepted for publication after minor revision

We exploit wide-field Ly$alpha$ imaging with Subaru to probe the environment around TN J1338-1942, a powerful radio galaxy with a >100 kpc Ly$\alpha$ halo at z=4.11. We used a sample of Ly$\alpha$ emitters (LAEs) down to $\log(L_{\rm Ly\alpha} [erg\, s^{-1}])\sim 42.8$ to measure the galaxy density around TNJ1338, compared to a control sample from a blank field taken with the same instrument. We found that TNJ1338 resides in a region with a peak overdensity of $\delta_{\rm LAE}=2.8\pm 0.5$ on scales of $8\, h^{-1}\rm Mpc$ (on the sky) and $112\, h^{-1}\rm Mpc$ (line of sight) in comoving coordinates. Adjacent to this overdensity, we found a strong underdensity where virtually no LAEs are detected. We used a semi-analytical model of LAEs derived from the Millennium Simulation to compare our results with theoretical predictions. While the theoretical density distribution is consistent with the blank field, overdense regions such as that around TNJ1338 are very rare, with a number density of $6.4\times 10^{-8}\rm Mpc^{-3}$ (comoving), corresponding to the densest < 0.4 percentile at $z\simeq 4.1$. We also found that the Ly$\alpha$ luminosity function in the TNJ1338 field differs from that in the blank field: the number of bright LAEs ($\log(L_{\rm Ly\alpha}[erg\,s^{-1}]) \gtrsim 43.3$) is enhanced, while the number of fainter LAEs is relatively suppressed. These results suggest that some powerful radio galaxies associated with Ly$\alpha$ nebulae reside in extreme overdensities on $\sim 3$--$6\, \rm Mpc$ scales, where star-formation and AGN activity may be enhanced via frequent galaxy mergers or high rates of gas accretion from the surroundings., Comment: 20 papges, 9 figures (incl. appendix), accepted for publication in MNRAS

Galaxy surveys targeting emission lines are characterising 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 photoionisation code \mapp. We characterise the interstellar medium (ISM) of galaxies by relating the ionisation parameter of gas in galaxies to their cold gas metallicity, obtaining a reasonable agreement with the observed \ha, \oii, \oiii\ luminosity functions, and the the BPT diagram for local star-forming galaxies. The average ionisation 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 $M_{\rm halo} \lesssim 10^{12} {[\rm h^{-1} M_{\odot}]}$. Finally, we exploit our modelling technique to predict galaxy number counts up to $z\sim 10$ by targeting far-infrared (FIR) emission lines detectable with submillimetre facilities, Comment: 16 pages, 9 figures. Published in MNRAS after minor revision

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 welcome

We search for evidence of diffuse Ly-alpha emission from extended neutral hydrogen surrounding Ly-alpha emitting galaxies (LAEs) using deep narrow-band images of the Extended Chandra Deep Field South. By stacking the profiles of 187 LAEs at z = 2.06, 241 LAEs at z = 3.10, and 179 LAEs at z = 3.12, and carefully performing low-surface brightness photometry, we obtain mean surface brightness maps that reach 9.9, 8.7, and 6.2 * 10^{-19} ergs cm^{-2} s^{-1} arcsec^{-2} in the emission line. We undertake a thorough investigation of systematic uncertainties in our surface brightness measurements, and find that our limits are 5--10 times larger than would be expected from Poisson background fluctuations; these uncertainties are often underestimated in the literature. At z ~ 3.1, we find evidence for extended halos with small scale lengths of 5--8 kpc in some, but not all of our sub-samples. We demonstrate that sub-samples of LAEs with low equivalent widths and brighter continuum magnitudes are more likely to possess such halos. At z ~ 2.1, we find no evidence of extended Ly-alpha emission down to our detection limits. Through Monte-Carlo simulations, we also show that we would have detected large diffuse LAE halos if they were present in our data sets. We compare these findings to other measurements in the literature, and discuss possible instrumental and astrophysical reasons for the discrepancies., Comment: 18 pages, 9 figures. Replaced with accepted version to be published in the Astrophysical Journal. Changes include expanded introduction and discussion of deep surface photometry techniques as well as additional numerical tests on the profiles. Primary conclusions remain unchanged

We study the properties of Ly-alpha emitters in a cosmological framework by computing the escape of Ly-alpha photons through galactic outflows. We combine the GALFORM semi-analytical model of galaxy formation with a Monte Carlo Ly-alpha radiative transfer code. The properties of Ly-alpha emitters at 0

A space-based galaxy redshift survey would have enormous power in constraining dark energy and testing general relativity, provided that its parameters are suitably optimized. We study viable space-based galaxy redshift surveys, exploring the dependence of the Dark Energy Task Force (DETF) figure-of-merit (FoM) on redshift accuracy, redshift range, survey area, target selection, and forecast method. Fitting formulae are provided for convenience. We also consider the dependence on the information used: the full galaxy power spectrum P(k, P(k) marginalized over its shape, or just the Baryon Acoustic Oscillations (BAO). We find that the inclusion of growth rate information (extracted using redshift space distortion and galaxy clustering amplitude measurements) leads to a factor of ~ 3 improvement in the FoM, assuming general relativity is not modified. This inclusion partially compensates for the loss of information when only the BAO are used to give geometrical constraints, rather than using the full P(k) as a standard ruler. We find that a space-based galaxy redshift survey covering ~20,000 deg^2 over 0.5< z < 2 with \sigma_z/(1+z)<= 0.001 exploits a redshift range that is only easily accessible from space, extends to sufficiently low redshifts to allow both a vast 3-D map of the universe using a single tracer population, and overlaps with ground-based surveys to enable robust modeling of systematic effects. We argue that these parameters are close to their optimal values given current instrumental and practical constraints., Comment: 13 pages, 5 figures. Version with slightly modified text. Accepted by MNRAS

The clustering of galaxies observed in future redshift surveys will provide a wealth of cosmological information. Matching the signal at different redshifts constrains the dark energy driving the acceleration of the expansion of the Universe. In tandem with these geometrical constraints, redshift-space distortions (RSD) depend on the build up of large-scale structure. As pointed out by many authors measurements of these effects are intrinsically coupled. We investigate this link, and argue that it strongly depends on the cosmological assumptions adopted when analysing data. Using representative assumptions for the parameters of the Euclid survey in order to provide a baseline future experiment, we show how the derived constraints change due to different model assumptions. We argue that even the assumption of a Friedman-Robertson-Walker (FRW) space-time is sufficient to reduce the importance of the coupling to a significant degree. Taking this idea further, we consider how the data would actually be analysed and argue that we should not expect to be able to simultaneously constrain multiple deviations from the standard $\Lambda$CDM model. We therefore consider different possible ways in which the Universe could deviate from the $\Lambda$CDM model, and show how the coupling between geometrical constraints and structure growth affects the measurement of such deviations., Comment: 11 pages, 5 figures. Matches published version. For the fisher matrix code email Lado Samushia

Future galaxy surveys will map the galaxy distribution in the redshift interval $0.5 -15.4$; a flux limit of at least $\logfha = -16$ is required for an \ha sample to become competitive in effective volume., Comment: submitted to MNRAS

We studied the clustering properties and multiwavelength spectral energy distributions of a complete sample of 162 Ly Alpha-Emitting (LAE) galaxies at z=3.1 discovered in deep narrow-band MUSYC imaging of the Extended Chandra Deep Field South. LAEs were selected to have observed frame equivalent widths >80A and emission line fluxes >1.5E-17 erg/cm^2/s. Only 1% of our LAE sample appears to host AGN. The LAEs exhibit a moderate spatial correlation length of r_0=3.6+0.8-1.0 Mpc, corresponding to a bias factor b=1.7+0.3-0.4, which implies median dark matter halo masses of log10(M_med) = 10.9+0.5-0.9 M_sun. Comparing the number density of LAEs, (1.5+-0.3)E-3/Mpc^3, with the number density of these halos finds a mean halo occupation ~1-10%. The evolution of galaxy bias with redshift implies that most z=3.1 LAEs evolve into present-day galaxies with L<2.5L*, whereas other z>3 galaxy populations typically evolve into more massive galaxies. Halo merger trees show that z=0 descendants occupy halos with a wide range of masses, with a median descendant mass close to that of L*. Only 30% of LAEs have sufficient stellar mass (>~3E9 M_sun) to yield detections in deep Spitzer-IRAC imaging. A two-population SED fit to the stacked UBVRIzJK+[3.6,4.5,5.6,8.0]micron fluxes of the IRAC-undetected objects finds that the typical LAE has low stellar mass (1.0+0.6-0.4 E9 M_sun), moderate star formation rate (2+-1 M_sun/yr), a young component age of 20+30-10 Myr, and little dust (A_V<0.2). The best fit model has 20% of the mass in the young stellar component, but models without evolved stars are also allowed., Comment: ApJ, in press, 7 pages including 4 color figures

Context: Many current and future surveys aim to detect the highest redshift (z >~ 7) sources through their Lyman-alpha (Ly-alpha) emission, using the narrow-band imaging method. However, to date the surveys have only yielded non-detections and upper limits as no survey has reached the necessary combination of depth and area to detect these very young star forming galaxies. Aims: We aim to calculate model luminosity functions and mock surveys of Ly-alpha emitters at z >~ 7 based on a variety of approaches. Methods: We calculate model luminosity functions at different redshifts based on three different approaches: a semi-analytical model based on CDM, a simple phenomenological model, and an extrapolation of observed Schechter functions at lower redshifts. The results of the first two models are compared with observations made at redshifts z ~ 5.7 and z ~ 6.5, and they are then extrapolated to higher redshift. Results: We present model luminosity functions for redshifts between z = 7 - 12.5 and give specific number predictions for future planned or possible narrow-band surveys for Ly-alpha emitters. We also investigate what constraints future observations will be able to place on the Ly-alpha luminosity function at very high redshift. Conclusion: It should be possible to observe z = 7 - 10 Ly-alpha emitters with present or near-future instruments if enough observing time is allocated. In particular, large area surveys such as ELVIS (Emission Line galaxies with VISTA Survey) will be useful in collecting a large sample. However, to get a large enough sample to constrain well the z >= 10 Ly-alpha luminosity function, instruments further in the future, such as an ELT, will be necessary., Comment: 8 pages, 4 figures, 3 tables, Accepted in A&A, best printed in colour

In the last decade, the technique of finding Ly-alpha emitters through narrow-band imaging has become a promising method of detecting high redshift galaxies. Ly-alpha emitters have been found from redshifts z ~ 2, up to the highest redshift source known to date at z = 6.96. Several surveys are also underway to find z = 7 - 9 sources. But these very high redshift sources are too faint to be studied in great detail, and more information can be found from studying the same class of objects at lower redshifts. Here we present our survey strategy to determine the nature of Ly-alpha emitters at lower redshifts, through multi-wavelength surveys, and our plans to extend the survey to redshift z = 8.8., Comment: 6 pages, 3 figures, to appear in the proceedings of 'At the Edge of the Universe', October 2006, Sintra, Portugal