The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

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.2<z<3.3. Our work is based on the Javalambre Photometric Local Universe Survey (J-PLUS) first data-release, which provides multiple narrow-band measurements over ~1000 deg^2, with limiting magnitude r~22. The analysis of high-z Lya-emitting sources over such a wide area is unprecedented, and allows to select a total of ~14,500 hyper-bright (L_Lya>10^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.2<z<3.3 Lya LF for the first time above L_Lya~10^44 erg/s and down to densities of ~10^-8 Mpc^-3. Our results unveil with high detail the Schechter exponential-decay of the brightest-end of the Lya LF, complementing the power-law component of previous LF determinations at 43.3<Log_10(L_Lya / [erg/s])<44. We measure Phi^*=(3.33+-0.19)x10^-6, Log(L^*)=44.65+-0.65 and alpha=-1.35+-0.84 as an average over the redshifts we probe. These values are significantly different than the typical Schechter parameters measured for the Lya LF of high-z star-forming LAEs. This suggests that z>2 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.

We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6<z<1.2$. We compare model predictions with different observational data sets, including DEEP2--Firefly galaxies with absolute magnitudes. We estimate the [OII] luminosity, L[OII], using simple relations derived both from the models and observations and also using a public code. This code ideally uses as input instantaneous star formation rates (SFRs), which are only provided by one of the SAMs under consideration. We use this SAM to study the feasibility of inferring galaxies' L[OII] for models that only provide average SFRs. We find that the post-processing computation of L[OII] from average SFRs is accurate for model galaxies with dust attenuated L[OII]$\lesssim10^{42.2}$erg s$^{-1}$ ($<5\%$ discrepancy). We also explore how to derive the [OII] luminosity from simple relations using global properties usually output by SAMs. Besides the SFR, the model L[OII] is best correlated with the observed-frame $u$ and $g$ broad-band magnitudes. These correlations have coefficients (r-values) above 0.64 and a dispersion that varies with L[OII]. We use these correlations and an observational one based on SFR and metallicity to derive L[OII]. These relations result in [OII] luminosity functions and halo occupation distributions with shapes that vary depending on both the model and the method used. Nevertheless, for all the considered models, the amplitude of the clustering at scales above 1$h^{-1}$Mpc remains unchanged independently of the method used to derive L[OII].

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.

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.

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.

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<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope.

We estimate the constraining power of J-PAS for parameters of an interacting dark energy cosmology. The survey is expected to map several millions of luminous red galaxies, emission line galaxies and quasars in an area of thousands of square degrees in the northern sky with precise photometric redshift measurements. Forecasts for the DESI and Euclid surveys are also evaluated and compared to J-PAS. With the Fisher matrix approach, we find that J-PAS can place constraints on the interaction parameter comparable to those from DESI, with an absolute uncertainty of about $0.02$, when the interaction term is proportional to the dark matter energy density, and almost as good, of about $0.01$, when the interaction is proportional to the dark energy density. For the equation of state of dark energy, the constraints from J-PAS are slightly better in the two cases (uncertainties $0.04$ - $0.05$ against $0.05$ - $0.07$ around the fiducial value $-1$). Both surveys stay behind Euclid but follow it closely, imposing comparable constraints in all specific cases considered.

Ultracool dwarfs (UCDs) are objects with spectral types equal or later than M7. Most of them have been discovered using wide-field imaging surveys. The Virtual Observatory (VO) has proven to be of great utility to efficiently exploit these astronomical resources. We aim to validate a VO methodology designed to discover and characterize UCDs in the J-PLUS photometric survey. J-PLUS is a multiband survey carried out with the wide angle T80Cam optical camera mounted on the 0.83-m telescope JAST/T80 in the Observatorio Astrof\'isico de Javalambre. In this work we make use of the Internal Data Release (IDR) covering 528 deg$^2$. We complement J-PLUS photometry with other catalogues in the optical and IR using VOSA, a VO tool that estimates physical parameters from the spectral energy distribution fitting to collections of theoretical models. Objects identified as UCDs are distinguished from background M giants and highly reddened stars using parallaxes and proper motions from Gaia DR2. We identify 559 UCDs, ranging from i=16.2 to 22.4 mag, of which 187 are candidate UCDs not previously reported in the literature. This represents an increase in the number of known UCDs of about 50% in the studied region of the sky, particularly at the faint end of our sensitivity, which is interesting as reference for future wide and deep surveys such as Euclid. Three candidates constitute interesting targets for exoplanet surveys because of their proximity (<40 pc). We also analyze the kinematics of UCDs in our catalogue and find evidence that it is consistent with a Galactic thin-disk population, except for 6 objects that might be members of the thick disk. The results obtained validate the proposed methodology, which will be used in future J-PLUS and J-PAS releases. Considering the region of the sky covered by the IDR used, we foresee to discover 3,000-3,500 new UCDs at the end of the J-PLUS project.

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.

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<z<4 using the same tracer, carried out using massively parallel wide field multi-object slit spectroscopy from space, will provide definitive dark energy measurements with minimal observational systematics by design. The HD GRS Wide will illuminate the nature of dark energy, and lead to revolutionary advances in particle physics and cosmology. It will also trace the cosmic web of dark matter and provide key insight into large-scale structure in the Universe. The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

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.

We present a new methodology for the estimation of stellar atmospheric parameters from narrow- and intermediate-band photometry of the Javalambre Photometric Local Universe Survey (J-PLUS), and propose a method for target pre-selection of low-metallicity stars for follow-up spectroscopic studies. Photometric metallicity estimates for stars in the globular cluster M15 are determined using this method. By development of a neural-network-based photometry pipeline, we aim to produce estimates of effective temperature, $T_{\rm eff}$, and metallicity, [Fe/H], for a large subset of stars in the J-PLUS footprint. The Stellar Photometric Index Network Explorer, SPHINX, is developed to produce estimates of $T_{\rm eff}$ and [Fe/H], after training on a combination of J-PLUS photometric inputs and synthetic magnitudes computed for medium-resolution (R ~ 2000) spectra of the Sloan Digital Sky Survey. This methodology is applied to J-PLUS photometry of the globular cluster M15. Effective temperature estimates made with J-PLUS Early Data Release photometry exhibit low scatter, \sigma($T_{\rm eff}$) = 91 K, over the temperature range 4500 < $T_{\rm eff}$ (K) < 8500. For stars from the J-PLUS First Data Release with 4500 < $T_{\rm eff}$ (K) < 6200, 85 $\pm$ 3% of stars known to have [Fe/H] <-2.0 are recovered by SPHINX. A mean metallicity of [Fe/H]=-2.32 $\pm$ 0.01, with a residual spread of 0.3 dex, is determined for M15 using J-PLUS photometry of 664 likely cluster members. We confirm the performance of SPHINX within the ranges specified, and verify its utility as a stand-alone tool for photometric estimation of effective temperature and metallicity, and for pre-selection of metal-poor spectroscopic targets.

In the present paper we aim to validate a methodology designed to extract the Halpha emission line flux from J-PLUS photometric data. J-PLUS is a multi narrow-band filter survey carried out with the 2 deg2 field of view T80Cam camera, mounted on the JAST/T80 telescope in the OAJ, Teruel, Spain. The information of the twelve J-PLUS bands, including the J0660 narrow-band filter located at rest-frame Halpha, is used over 42 deg2 to extract de-reddened and [NII] decontaminated Halpha emission line fluxes of 46 star-forming regions with previous SDSS and/or CALIFA spectroscopic information. The agreement of the inferred J-PLUS photometric Halpha fluxes and those obtained with spectroscopic data is remarkable, with a median comparison ratio R = 1.05 +- 0.25. This demonstrates that it is possible to retrieve reliable Halpha emission line fluxes from J-PLUS photometric data. With an expected area of thousands of square degrees upon completion, the J-PLUS dataset will allow the study of several star formation science cases in the nearby universe, as the spatially resolved star formation rate of nearby galaxies at z < 0.015, and how it is influenced by the environment, morphology or nuclear activity. As an illustrative example, the close pair of interacting galaxies NGC3994 and NGC3995 is analyzed, finding an enhancement of the star formation rate not only in the center, but also in outer parts of the disk of NGC3994.

The spatial variations of stellar population properties within a galaxy are intimately related to their formation process. Therefore, spatially resolved studies of galaxies are essential to uncover their formation and assembly. The Javalambre Photometric Local Universe Survey (J-PLUS) is a dedicated multi-filter designed to observed ~8500 deg2 using twelve narrow-, intermediate- and broad-band filters in the optical range. In this study, we test the potential of the multi-filter observation carried out with J-PLUS to investigate the properties of spatially-resolved nearby galaxies. We present detailed 2D maps of stellar population properties (age, metallicity, extinction, and stellar mass surface density) for two early-type galaxies observed in both, J-PLUS and CALIFA surveys: NGC 5473 and NGC 5485. Radial structures are also compared and luminosity- and mass-weighted profiles are derived. We use MUFFIT to process the J-PLUS observations, and two different techniques (STARLIGHT and STECKMAP) to analyze IFU CALIFA data. We demonstrate that this novel technique delivers radial stellar population gradients in good agreement with the IFU technique CALIFA/STECKMAP although comparison of the absolute values reveals the existence of intrinsic systematic differences. Radial stellar population gradients differ when CALIFA/STARLIGHT methodology is used. Age and metallicity radial profiles derived from J-PLUS/MUFFIT are very similar when luminosity- or mass-weighted properties are used, suggesting that the contribution of a younger component is small. Comparison between the three methodologies reveals some discrepancies suggesting that the specific characteristics of each method causes important differences. We conclude that the ages, metallicities and extinction derived for individual galaxies not only depend on the chosen models but also depend on the method used.

J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrof\'isico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a combination of broad, medium and narrow-band filters, optimally designed to extract the rest-frame spectral features (the 3700-4000\AA\ Balmer break region, H$\delta$, Ca H+K, the G-band, the Mgb and Ca triplets) that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each pixel of the sky observed. With a typical depth of AB $\sim 21.25$ mag per band, this filter set thus allows for an indiscriminate and accurate characterization of the stellar population in our Galaxy, it provides an unprecedented 2D photo-spectral information for all resolved galaxies in the local universe, as well as accurate photo-z estimates ($\Delta\,z\sim 0.01-0.03$) for moderately bright (up to $r\sim 20$ mag) extragalactic sources. While some narrow band filters are designed for the study of particular emission features ([OII]/$\lambda$3727, H$\alpha$/$\lambda$6563) up to $z < 0.015$, they also provide well-defined windows for the analysis of other emission lines at higher redshifts. As a result, J-PLUS has the potential to contribute to a wide range of fields in Astrophysics, both in the nearby universe (Milky Way, 2D IFU-like studies, stellar populations of nearby and moderate redshift galaxies, clusters of galaxies) and at high redshifts (ELGs at $z\approx 0.77, 2.2$ and $4.4$, QSOs, etc). With this paper, we release $\sim 36$ sq.deg of J-PLUS data, containing about $1.5\times 10^5$ stars and $10^5$ galaxies at $r<21$ mag.

Our goal is to morphologically classify the sources identified in the images of the J-PLUS early data release (EDR) into compact (stars) or extended (galaxies) using a suited Bayesian classifier. J-PLUS sources exhibit two distinct populations in the r-band magnitude vs. concentration plane, corresponding to compact and extended sources. We modelled the two-population distribution with a skewed Gaussian for compact objects and a log-normal function for the extended ones. The derived model and the number density prior based on J-PLUS EDR data were used to estimate the Bayesian probability of a source to be star or galaxy. This procedure was applied pointing-by-pointing to account for varying observing conditions and sky position. Finally, we combined the morphological information from g, r, and i broad bands in order to improve the classification of low signal-to-noise sources. The derived probabilities are used to compute the pointing-by-pointing number counts of stars and galaxies. The former increases as we approach to the Milky Way disk, and the latter are similar across the probed area. The comparison with SDSS in the common regions is satisfactory up to r ~ 21, with consistent numbers of stars and galaxies, and consistent distributions in concentration and (g - i) colour spaces. We implement a morphological star/galaxy classifier based on PDF analysis, providing meaningful probabilities for J-PLUS sources to one magnitude deeper (r ~ 21) than a classical boolean classification. These probabilities are suited for the statistical study of 150k stars and 101k galaxies with 15 < r < 21 present in the 31.7 deg2 of the J-PLUS EDR. In a future version of the classifier, we will include J-PLUS colour information from twelve photometric bands.

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.

Emission line galaxies (ELGs) are used in several ongoing and upcoming surveys (SDSS-IV/eBOSS, DESI) as tracers of the dark matter distribution. Using a new galaxy formation model, we explore the characteristics of [OII] emitters, which dominate optical ELG selections at $z\simeq 1$. Model [OII] emitters at $0.5<z<1.5$ are selected to mimic the DEEP2, VVDS, eBOSS and DESI surveys. The luminosity functions of model [OII] emitters are in reasonable agreement with observations. The selected [OII] emitters are hosted by haloes with $M_{\rm halo}\geq 10^{10.3}h^{-1}{\rm M}_{\odot}$, with ~90% of them being central star-forming galaxies. The predicted mean halo occupation distributions of [OII] emitters has a shape typical of that inferred for star-forming galaxies, with the contribution from central galaxies, $\langle N \rangle_{\left[OII\right]\, cen}$, being far from the canonical step function. The $\langle N \rangle_{\left[OII\right]\, cen}$ can be described as the sum of an asymmetric Gaussian for disks and a step function for spheroids, which plateaus below unity. The model [OII] emitters have a clustering bias close to unity, which is below the expectations for eBOSS and DESI ELGs. At $z\sim 1$, a comparison with observed g-band selected galaxy, which are expected to be dominated by [OII] emitters, indicates that our model produces too few [OII] emitters that are satellite galaxies. This suggests the need to revise our modelling of hot gas stripping in satellite galaxies.

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.

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.

We present a new color-selection technique, based on the Bruzual & Charlot models convolved with the bands of the ALHAMBRA survey, and the redshifted position of the Balmer jump to select star-forming galaxies in the redshift range 0.5 < z < 1.5. These galaxies are dubbed Balmer jump Galaxies BJGs. We apply the iSEDfit Bayesian approach to fit each detailed SED and determine star-formation rate (SFR), stellar mass, age and absolute magnitudes. The mass of the haloes where these samples reside are found via a clustering analysis. Five volume-limited BJG sub-samples with different mean redshifts are found to reside in haloes of median masses $\sim 10^{12.5 \pm 0.2} M_\odot$ slightly increasing toward z=0.5. This increment is similar to numerical simulations results which suggests that we are tracing the evolution of an evolving population of haloes as they grow to reach a mass of $\sim 10^{12.7 \pm 0.1} M_\odot$ at z=0.5. The likely progenitors of our samples at z$\sim$3 are Lyman Break Galaxies, which at z$\sim$2 would evolve into star-forming BzK galaxies, and their descendants in the local Universe are elliptical galaxies.Hence, this allows us to follow the putative evolution of the SFR, stellar mass and age of these galaxies. From z$\sim$1.0 to z$\sim$0.5, the stellar mass of the volume limited BJG samples nearly does not change with redshift, suggesting that major mergers play a minor role on the evolution of these galaxies. The SFR evolution accounts for the small variations of stellar mass, suggesting that star formation and possible minor mergers are the main channels of mass assembly.

We present a survey of $z\sim3$ Ly$\alpha$ emitters (LAEs) within the fields of the VLT LBG Redshift Survey. The data encompasses 5 independent survey fields co-spatial with spectroscopic LBG data and covering a larger total area than previously analysed for LAE number counts and clustering. This affords an improved analysis over previous work by minimising the effects of cosmic variance and allowing the cross-clustering analysis of LAEs and LBGs. Our photometric sample consists of $\approx600$ LAE candidates, over an area of 1.07~deg$^2$, with equivalent widths of $\gtrsim65$~\AA\ and a flux limit of $\approx2\times10^{-17}$~erg~cm$^{-2}$~s$^{-1}$. From spectroscopic follow-up, we measured a success rate of $78\pm18\%$. We find the $R$-band continuum luminosity function to be $\sim10\times$ lower than the luminosity function of LBGs at this redshift, consistent with previous studies. Exploiting the large area of the survey, we estimate the LAE auto-correlation function and find a clustering length of $r_0=2.86\pm0.33~h^{-1}$~Mpc, low compared to the $z\sim3$ LBG population, but somewhat higher than previous LAE measurements. This corresponds to a median halo mass of $M_{\rm DM}=10^{11.0\pm0.3}~h^{-1}~$M$_{\odot}$. We present an analysis of clustering length versus continuum magnitude and find that the measurements for LAEs and LBGs are consistent at faint magnitudes. Our combined dataset of LAEs and LBGs allows us to measure, for the first time, the LBG-LAE cross-correlation, finding a clustering length of $r_0=3.29\pm0.57~h^{-1}$~Mpc and a LAE halo mass of $10^{11.1\pm0.4}~h^{-1}$~M$_{\odot}$. Overall, we conclude that LAEs inhabit primarily low mass halos, but form a relatively small proportion of the galaxy population found in such halos.

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.

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

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.

We make use of a semi-analytical model of galaxy formation to investigate the origin of the observed correlation between [a/Fe] abundance ratios and stellar mass in elliptical galaxies. We implement a new galaxy-wide stellar initial mass function (Top Heavy Integrated Galaxy Initial Mass Function, TH-IGIMF) in the semi-analytic model SAG and evaluate its impact on the chemical evolution of galaxies. The SFR-dependence of the slope of the TH-IGIMF is found to be key to reproducing the correct [a/Fe]-stellar mass relation. Massive galaxies reach higher [a/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 low as 5 solar masses and 2, respectively. A mild 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 [a/Fe] values. Hence, the time-scale of galaxy formation alone cannot explain the slope of the [a/Fe]-galaxy mass relation, but is responsible for the big dispersion of [a/Fe] abundance ratios at fixed stellar mass.We further test the hyphothesis 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.

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.

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.

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<z<7 are predicted using two outflow geometries: a Shell of neutral gas and a Wind ejecting material, both expanding at constant velocity. We characterise the differences in the Ly-alpha line profiles predicted by the two outflow geometries in terms of their width, asymmetry and shift from the line centre for a set of outflows with different hydrogen column densities, expansion velocities and metallicities. In general, the Ly-alpha line profile of the Shell geometry is broader and more asymmetric, and the Ly-alpha escape fraction is lower than with the Wind geometry for the same set of parameters. In order to implement the outflow geometries in the semi-analytical model GALFORM, a number of free parameters in the outflow model are set by matching the luminosity function of Ly-alpha emitters over the whole observed redshift range. The models are consistent with the observationally inferred Ly-alpha escape fractions, equivalent width distributions and with the shape of the Ly-alpha line from composite spectra. Interestingly, our predicted UV luminosity function of Ly-alpha emitters and the fraction of Ly-alpha emitters in Lyman-break galaxy samples at high redshift are in partial agreement with observations. Attenuation of the Ly-alpha line by the presence of a neutral intergalactic medium at high redshift could be responsible for this disagreement. We predict that Ly-alpha emitters constitute a subset of the galaxy population with lower metallicities, lower instantaneous star formation rates and larger sizes than the overall population at the same UV luminosity.

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.

We use a model for the evolution of galaxies in the far-IR based on the LambdaCDM cosmology to make detailed predictions for upcoming cosmological surveys with the Herschel Space Observatory. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al. (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF appears necessary to explain the number counts and redshifts of faint sub-mm galaxies. In this paper, we present predictions for galaxy luminosity functions, number counts and redshift distributions in the Herschel imaging bands. We find that source confusion will be a serious problem in the deepest planned surveys. We also show predictions for physical properties such as star formation rates and stellar, gas and halo masses, together with fluxes at other wavelengths (from the far-UV to the radio) relevant for multi-wavelength follow-up observations. We investigate what fraction of the total IR emission from dust and of the high-mass star formation over the history of the Universe should be resolved by planned surveys with Herschel, and find a fraction ~30-50%, depending on confusion. Finally, we show that galaxies in Herschel surveys should be significantly clustered.

Future galaxy redshift surveys aim to measure cosmological quantities from the galaxy power spectrum. A prime example is the detection of baryonic acoustic oscillations (BAOs), providing a standard ruler to measure the dark energy equation of state, w(z), to high precision. The strongest practical limitation for these experiments is how quickly accurate redshifts can be measured for sufficient galaxies to map the large-scale structure. A promising strategy is to target emission-line (i.e. star-forming) galaxies at high-redshift (z~0.5-2); not only is the space density of this population increasing out to z~2, but also emission-lines provide an efficient method of redshift determination. Motivated by the prospect of future dark energy surveys targeting H-alpha emitters at near-infrared wavelengths (i.e. z>0.5), we use the latest empirical data to model the evolution of the H-alpha luminosity function out to z~2, and thus provide predictions for the abundance of H-alpha emitters for practical limiting fluxes. We caution that the estimates presented in this work must be tempered by an efficiency factor, epsilon, giving the redshift success rate from these potential targets. For a range of practical efficiencies and limiting fluxes, we provide an estimate of nP_{0.2}, where n is the 3D galaxy number density and P_{0.2} is the galaxy power spectrum evaluated at k=0.2h/Mpc. Ideal surveys must provide nP_{0.2}>1 in order to balance shot-noise and cosmic variance errors. We show that a realistic emission-line survey (epsilon=0.5) could achieve nP_{0.2}=1 out to z~1.5 with a limiting flux of 10^{-16} erg/s/cm^{-2}. If the limiting flux is a factor 5 brighter, then this goal can only be achieved out to z~0.5, highlighting the importance of survey depth and efficiency in cosmological redshift surveys.

Future galaxy surveys will map the galaxy distribution in the redshift interval $0.5<z<2$ using near-infrared cameras and spectrographs. The primary science goal of such surveys is to constrain the nature of the dark energy by measuring the large-scale structure of the Universe. This requires a tracer of the underlying dark matter which maximizes the useful volume of the survey. We investigate two potential survey selection methods: an emission line sample based on the \ha line and a sample selected in the H-band. We present predictions for the abundance and clustering of such galaxies, using two published versions of the \galform galaxy formation model. Our models predict that \ha selected galaxies tend to avoid massive dark matter haloes and instead trace the surrounding filamentary structure; H-band selected galaxies, on the other hand, are found in the highest mass haloes. This has implications for the measurement of the rate at which fluctuations grow due to gravitational instability. We use mock catalogues to compare the effective volumes sampled by a range of survey configurations. To give just two examples: a redshift survey down to $H_{\rm AB}=22$ samples an effective volume that is $\sim 5-10$ times larger than that probed by an \ha survey with $\logfha > -15.4$; a flux limit of at least $\logfha = -16$ is required for an \ha sample to become competitive in effective volume.

We combine a semi-analytical model of galaxy formation with a very large simulation which follows the growth of large scale structure in a LambdaCDM universe to predict the clustering of Ly-alpha emitters. We find that the clustering strength of Ly-alpha emitters has only a weak dependence on Ly-alpha luminosity but a strong dependence on redshift. With increasing redshift, Ly-alpha emitters trace progressively rarer, higher density regions of the universe. Due to the large volume of the simulation, over 100 times bigger than any previously used for this application, we can construct mock catalogues of Ly-alpha emitters and study the sample variance of current and forthcoming surveys. We find that the number and clustering of Ly-alpha emitters in our mock catalogues are in agreement with measurements from current surveys, but that there is a considerable scatter in these quantities. We argue that a proposed survey of emitters at z=8.8 should be extended significantly in solid angle to allow a robust measurement of Ly-alpha emitter clustering.

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.

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.

We describe the results of an extremely deep, 0.28 deg^2 survey for z = 3.1 Ly-alpha emission-line galaxies in the Extended Chandra Deep Field South. By using a narrow-band 5000 Anstrom filter and complementary broadband photometry from the MUSYC survey, we identify a statistically complete sample of 162 galaxies with monochromatic fluxes brighter than 1.5 x 10^-17 ergs cm^-2 s^-1 and observers frame equivalent widths greater than 80 Angstroms. We show that the equivalent width distribution of these objects follows an exponential with a rest-frame scale length of w_0 = 76 +/- 10 Angstroms. In addition, we show that in the emission line, the luminosity function of Ly-alpha galaxies has a faint-end power-law slope of alpha = -1.49 +/- 0.4, a bright-end cutoff of log L^* = 42.64 +/- 0.2, and a space density above our detection thresholds of 1.46 +/- 0.12 x 10^-3 h70^3 galaxies Mpc^-3. Finally, by comparing the emission-line and continuum properties of the LAEs, we show that the star-formation rates derived from Ly-alpha are ~3 times lower than those inferred from the rest-frame UV continuum. We use this offset to deduce the existence of a small amount of internal extinction within the host galaxies. This extinction, coupled with the lack of extremely-high equivalent width emitters, argues that these galaxies are not primordial Pop III objects, though they are young and relatively chemically unevolved.

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.