Observations have revealed massive (logM*/Msun>11) galaxies that were already dead when the universe was only ~2 Gyr. Given the short time before these galaxies were quenched, their past histories and quenching mechanism(s) are of particular interest. In this paper, we study star formation histories (SFHs) of 24 massive galaxies at 1.6<z<2.5. A deep slitless spectroscopy + imaging data set collected from multiple Hubble Space Telescope surveys allows robust determination of their spectral energy distributions and SFHs with no functional assumption on their forms. We find that most of our massive galaxies had formed > 50% of their extant masses by ~1.5 Gyr before the time of observed redshifts, with a trend where more massive galaxies form earlier. Their stellar-phase metallicities are already compatible with those of local early-type galaxies, with a median value of logZ*/Zsun=0.25 and scatter of ~0.15dex. In combination with the reconstructed SFHs, we reveal their rapid metallicity evolution from z~5.5 to ~2.2 at a rate of ~0.2dex/Gyr in log Z*/Zsun. Interestingly, the inferred stellar-phase metallicities are, when compared at half-mass time, ~0.25dex higher than observed gas-phase metallicities of star forming galaxies. While systematic uncertainties remain, this may imply that these quenched galaxies have continued low-level star formation, rather than abruptly terminating their star formation activity, and kept enhancing their metallicity until recently.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

We present a full data analysis of the pure-parallel Hubble Space Telescope (HST) imaging observations in the Brightest of Reionizing Galaxies Survey (BoRG[z9]) in Cycle 22. The medium-deep exposures with five HST/WFC3IR+UVIS filter bands from 79 independent sightlines (~370 arcmin^2) provide the least biased determination of number density for z>9 bright galaxies against cosmic variance. After a strict two-step selection for candidate galaxies, including dropout color and photometric redshift analyses, and revision of previous BoRG candidates, we identify one source at z~10 and two sources at z~9. The z~10 candidate shows evidence of line-of-sight lens magnification (mu~1.5), yet it appears surprisingly luminous (MUV ~ -22.6\pm0.3 mag), making it one of the brightest candidates at z > 8 known (~ 0.3 mag brighter than the z = 8.68 galaxy EGSY8p7, spectroscopically confirmed by Zitrin and collaborators). For z ~ 9 candidates, we include previous data points at fainter magnitudes and find that the data are well fitted by a Schechter luminosity function with alpha ~ -2.1, MUV ~ -21.5 mag, and log phi ~ -4.5 Mpc^-3mag^-1, for the first time without fixing any parameters. The inferred cosmic star formation rate density is consistent with unaccelerated evolution from lower redshift.

Local massive early-type galaxies are believed to have completed most of their star formation $\sim10$Gyr ago and evolved without having substantial star formation since. If so, their progenitors should have roughly solar stellar metallicities ($Z_*$), comparable to their values today. We report the discovery of two lensed massive ($\log M_*/M_\odot\sim11$), $z\sim2.2$ dead galaxies, that appear markedly metal deficient given this scenario. Using 17-band $HST$+$K_{s}$+$Spitzer$ photometry and deep $HST$ grism spectra from the GLASS and SN Refsdal follow-up campaigns covering features near $\lambda_{\rm rest}\sim4000$\AA, we find these systems to be dominated by A-type stars with $\log Z_*/Z_\odot=-0.40\pm0.02$ and $-0.49\pm0.03$ ($30$-$40\%$ solar) under standard assumptions. The second system's lower metallicity is robust to isochrone changes, though this choice can drive the first system's from $\log Z_*/Z_\odot=-0.6$ to 0.1. If these two galaxies are representative of larger samples, this finding suggests that evolutionary paths other than dry minor-merging are required for these massive galaxies. Future analyses with direct metallicity measurements-e.g., by the $James\ Webb\ Space\ Telescope$-will provide critical insight into the nature of such phenomena.

The selection of high redshift sources from broad-band photometry using the Lyman-break galaxy (LBG) technique is a well established methodology, but the characterization of its contamination for the faintest sources is still incomplete. We use the optical and near-IR data from four (ultra)deep Hubble Space Telescope legacy fields to investigate the contamination fraction of LBG samples at z~5-8 selected using a colour-colour method. Our approach is based on characterizing the number count distribution of interloper sources, that is galaxies with colors similar to those of LBGs, but showing detection at wavelengths shorter than the spectral break. Without sufficient sensitivity at bluer wavelengths, a subset of interlopers may not be properly classified, and contaminate the LBG selection. The surface density of interlopers in the sky gets steeper with increasing redshift of LBG selections. Since the intrinsic number of dropouts decreases significantly with increasing redshift, this implies increasing contamination from misclassified interlopers with increasing redshift, primarily by intermediate redshift sources with unremarkable properties (intermediate ages, lack of ongoing star formation and low/moderate dust content). Using Monte Carlo simulations, we estimate that the CANDELS deep data have contamination induced by photometric scatter increasing from ~2% at z~5 to ~6% at z~8 for a typical dropout color >1 mag, with contamination naturally decreasing for a more stringent dropout selection. Contaminants are expected to be located preferentially near the detection limit of surveys, ranging from 0.1 to 0.4 contaminants per arcmin2 at J=30, depending on the field considered. This analysis suggests that the impact of contamination in future studies of z>10 galaxies needs to be carefully considered.

Supernovae are important probes of the properties of stars at high redshifts because they can be detected at early epochs and their masses can be inferred from their light curves. Finding the first cosmic explosions in the universe will only be possible with the James Webb Space Telescope, the Wide-Field Infrared Survey Telescope and the next generation of extremely large telescopes. But strong gravitational lensing by massive clusters, like those in the Cluster Lensing and Supernova Survey with Hubble (CLASH), could reveal such events now by magnifying their flux by factors of 10 or more. We find that CLASH will likely discover at least 2 - 3 core-collapse supernovae at 5 < z < 12 and perhaps as many as ten. Future surveys of cluster lenses similar in scope to CLASH by the James Webb Space Telescope might find hundreds of these events out to z ~ 15 - 17. Besides revealing the masses of early stars, these ancient supernovae will also constrain cosmic star formation rates in the era of first galaxy formation.

The Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) enabled the search for the first galaxies observed at z ~ 8 - 11 (500 - 700 Myr after the Big Bang). To continue quantifying the number density of the most luminous galaxies (M_AB ~ -22.0) at the earliest epoch observable with HST, we search for z ~ 10 galaxies (F125W-dropouts) in archival data from the Brightest of Reionizing Galaxies (BoRG[z8]) survey, originally designed for detection of z ~ 8 galaxies (F098M-dropouts). By focusing on the deepest 293 arcmin^2 of the data along 62 independent lines of sight, we identify six z ~ 10 candidates satisfying the color selection criteria, detected at S/N > 8 in F160W with M_AB = -22.8 to -21.1 if at z = 10. Three of the six sources, including the two brightest, are in a single WFC3 pointing (~ 4 arcmin^2), suggestive of significant clustering, which is expected from bright galaxies at z ~ 10. However, the two brightest galaxies are too extended to be likely at z ~ 10, and one additional source is unresolved and possibly a brown dwarf. The remaining three candidates have m_AB ~ 26, and given the area and completeness of our search, our best estimate is a number density of sources that is marginally higher but consistent at 2{\sigma} with searches in legacy fields. Our study highlights that z ~ 10 searches can yield a small number of candidates, making tailored follow-ups of HST pure-parallel observations viable and effective.

We present a tally of Milky Way late-type dwarf stars in 68 WFC3 pure-parallel fields (227 arcmin^2) from the Brightest of Reionizing Galaxies (BoRG) survey for high-redshift galaxies. Using spectroscopically identified M-dwarfs in two public surveys, the CANDELS and the ERS mosaics, we identify a morphological selection criterion using the half-light radius (r50), a near-infrared J-H, G-J color region where M-dwarfs are found, and a V-J relation with M-dwarf subtype. We apply this morphological selection of stellar objects, color-color selection of M-dwarfs and optical-near-infrared color subtyping to compile a catalog of 274 M-dwarfs belonging to the disk of the Milky Way with a limiting magnitude of m_F125W < 24. Based on the M-dwarfs statistics, we conclude that (a) the previously identified North/South discrepancy in M-dwarf numbers persists in our sample; there are more M-dwarfs in the Northern fields on average than in Southern ones, (b) the Milky Way's single disk scale-height for M-dwarfs is 0.3-4 kpc, depending on sub-type, (c) {\bf ERRATUM:} we present corrected coordinates (AstroPy) and distances and find a constant $z_0$=600 pc for all types. (d) a second component is visible in the vertical distribution, with a different, much higher scale-height. We report the M-dwarf component of the Sagittarius stream in one of our fields with 11 confirmed M-dwarfs, 7 of which are at the stream's distance. The dwarf scale-height and the relative low incidence in our fields of L- and T-dwarfs in these fields makes it unlikely that these stars will be interlopers in great numbers in color-selected samples of high-redshift galaxies. The relative ubiquity of M-dwarfs however will make them ideal tracers of Galactic Halo substructure with EUCLID and reference stars for JWST observations.

We present the first results and design from the redshift z~9-10 Brightest of the Reionizing Galaxies {\it Hubble Space Telescope} survey BoRG[z9-10], aimed at searching for intrinsically luminous unlensed galaxies during the first 700 Myr after the Big Bang. BoRG[z9-10] is the continuation of a multi-year pure-parallel near-IR and optical imaging campaign with the Wide Field Camera 3. The ongoing survey uses five filters, optimized for detecting the most distant objects and offering continuous wavelength coverage from {\lambda}=0.35{\mu}m to {\lambda}=1.7{\mu}m. We analyze the initial ~130 arcmin$^2$ of area over 28 independent lines of sight (~25% of the total planned) to search for z>7 galaxies using a combination of Lyman break and photometric redshift selections. From an effective comoving volume of (5-25) $times 10^5$ Mpc$^3$ for magnitudes brighter than $m_{AB}=26.5-24.0$ in the $H_{160}$-band respectively, we find five galaxy candidates at z~8.3-10 detected at high confidence (S/N>8), including a source at z~8.4 with mAB=24.5 (S/N~22), which, if confirmed, would be the brightest galaxy identified at such early times (z>8). In addition, BoRG[z9-10] data yield four galaxies with $7.3 \lesssim z \lesssim 8$. These new Lyman break galaxies with m$\lesssim26.5$ are ideal targets for follow-up observations from ground and space based observatories to help investigate the complex interplay between dark matter growth, galaxy assembly, and reionization.

We give an overview of the Grism Lens Amplified Survey from Space (GLASS), a large Hubble Space Telescope program aimed at obtaining grism spectroscopy of the fields of ten massive clusters of galaxies at redshift z=0.308-0.686, including the Hubble Frontier Fields (HFF). The Wide Field Camera 3 yields near infrared spectra of the cluster cores, covering the wavelength range 0.81-1.69mum through grisms G102 and G141, while the Advanced Camera for Surveys in parallel mode provides G800L spectra of the infall regions of the clusters. The WFC3 spectra are taken at two almost orthogonal position angles in order to minimize the effects of confusion. After summarizing the scientific drivers of GLASS, we describe the sample selection as well as the observing strategy and data processing pipeline. We then utilize MACSJ0717.5+3745, a HFF cluster and the first one observed by GLASS, to illustrate the data quality and the high-level data products. Each spectrum brighter than H_AB=23 is visually inspected by at least two co-authors and a redshift is measured when sufficient information is present in the spectra. Furthermore, we conducted a thorough search for emission lines through all the GLASS WFC3 spectra with the aim of measuring redshifts for sources with continuum fainter than H_AB=23. We provide a catalog of 139 emission-line based spectroscopic redshifts for extragalactic sources, including three new redshifts of multiple image systems (one probable, two tentative). In addition to the data itself we also release software tools that are helpful to navigate the data.

The IRAC ultradeep field (IUDF) and IRAC Legacy over GOODS (IGOODS) programs are two ultradeep imaging surveys at 3.6{\mu}m and 4.5{\mu}m with the Spitzer Infrared Array Camera (IRAC). The primary aim is to directly detect the infrared light of reionization epoch galaxies at z > 7 and to constrain their stellar populations. The observations cover the Hubble Ultra Deep Field (HUDF), including the two HUDF parallel fields, and the CANDELS/GOODS-South, and are combined with archival data from all previous deep programs into one ultradeep dataset. The resulting imaging reaches unprecedented coverage in IRAC 3.6{\mu}m and 4.5{\mu}m ranging from > 50 hour over 150 arcmin^2, > 100 hour over 60 sq arcmin2, to 200 hour over 5 - 10 arcmin$^2$. This paper presents the survey description, data reduction, and public release of reduced mosaics on the same astrometric system as the CANDELS/GOODS-South WFC3 data. To facilitate prior-based WFC3+IRAC photometry, we introduce a new method to create high signal-to-noise PSFs from the IRAC data and reconstruct the complex spatial variation due to survey geometry. The PSF maps are included in the release, as are registered maps of subsets of the data to enable reliability and variability studies. Simulations show that the noise in the ultradeep IRAC images decreases approximately as the square root of integration time over the range 20 - 200 hours, well below the classical confusion limit, reaching 1{\sigma} point source sensitivities as faint as of 15 nJy (28.5 AB) at 3.6{\mu}m and 18 nJy (28.3 AB) at 4.5{\mu}m. The value of such ultradeep IRAC data is illustrated by direct detections of z = 7 - 8 galaxies as faint as HAB = 28.

We present a Bayesian framework to account for the magnification bias from both strong and weak gravitational lensing in estimates of high-redshift galaxy luminosity functions. We illustrate our method by estimating the $z\sim8$ UV luminosity function using a sample of 97 Y-band dropouts (Lyman break galaxies) found in the Brightest of Reionizing Galaxies (BoRG) survey and from the literature. We find the luminosity function is well described by a Schechter function with characteristic magnitude of $M^\star = -19.85^{+0.30}_{-0.35}$, faint-end slope of $\alpha = -1.72^{+0.30}_{-0.29}$, and number density of $\log_{10} \Psi^\star [\textrm{Mpc}^{-3}] = -3.00^{+0.23}_{-0.31}$. These parameters are consistent within the uncertainties with those inferred from the same sample without accounting for the magnification bias, demonstrating that the effect is small for current surveys at $z\sim8$, and cannot account for the apparent overdensity of bright galaxies compared to a Schechter function found recently by Bowler et al. (2014a,b) and Finkelstein et al. (2014). We estimate that the probability of finding a strongly lensed $z\sim8$ source in our sample is in the range $\sim 3-15 \%$ depending on limiting magnitude. We identify one strongly-lensed candidate and three cases of intermediate lensing in BoRG (estimated magnification $\mu>1.4$) in addition to the previously known candidate group-scale strong lens. Using a range of theoretical luminosity functions we conclude that magnification bias will dominate wide field surveys -- such as those planned for the Euclid and WFIRST missions -- especially at $z>10$. Magnification bias will need to be accounted for in order to derive accurate estimates of high-redshift luminosity functions in these surveys and to distinguish between galaxy formation models.

The discovery of 150 - 300 M$_{\odot}$ stars in the Local Group and pair-instability supernova candidates at low redshifts has excited interest in this exotic explosion mechanism. Realistic light curves for pair-instability supernovae at near-solar metallicities are key to identifying and properly interpreting these events as more are found. We have modeled pair-instability supernovae of 150 - 500 M$_{\odot}$ Z $\sim$ 0.1 - 0.4 Z$_{\odot}$ stars. These stars lose up to 80% of their mass to strong line-driven winds and explode as bare He cores. We find that their light curves and spectra are quite different from those of Population III pair-instability explosions, which therefore cannot be used as templates for low-redshift events. Although non-zero metallicity pair-instability supernovae are generally dimmer than their Population III counterparts, in some cases they will be bright enough to be detected at the earliest epochs at which they can occur, the formation of the first galaxies at $z \sim$ 10 - 15. Others can masquerade as dim, short duration supernovae that are only visible in the local universe and that under the right conditions could be hidden in a wide variety of supernova classes. We also report for the first time that some pair-instability explosions can create black holes with masses of $\sim$ 100 M$_{\odot}$.

Cosmological surface brightness dimming of the form $(1+z)^{-4}$ affects all sources. The strong dependence of surface brightness dimming on redshift z suggests the presence of a selection bias when searching for high-redshift galaxies, i.e. we tend to detect only those galaxies with a high surface brightness (SB). However, unresolved knots of emission are not affected by SB dimming, thus providing a way to test the clumpiness of high-z galaxies. Our strategy relies on the comparison of the total flux detected for the same source in surveys characterized by different depth. For all galaxies, deeper images permit the better investigation of low-SB features. Cosmological SB dimming makes these low-SB features hard to detect when going to higher and higher redshifts. We used the GOODS and HUDF Hubble Space Telescope legacy datasets to study the effect of SB dimming on low-SB features of high-redshift galaxies and compare it to the prediction for smooth sources. We selected a sample of Lyman-break galaxies at z~4 (i.e. B-band dropouts) detected in all of the datasets and found no significant trend when comparing the total magnitudes measured from images with different depth. Through Monte Carlo simulations we derived the expected trend for galaxies with different SB profiles. The comparison to the datahints at a compact distribution for most of the rest-frame ultraviolet light emitted from high-z galaxies.

Lyman-Break Galaxy (LBG) samples observed during reionization ($z\gtrsim6$) with Hubble Space Telescope's Wide Field Camera 3 are reaching sizes sufficient to characterize their clustering properties. Using a combined catalog from the Hubble eXtreme Deep Field and CANDELS surveys, containing $N=743$ LBG candidates at z>6.5 at a mean redshift of $z=7.2$, we detect a clear clustering signal in the angular correlation function (ACF) at $\sim4\sigma$, corresponding to a real-space correlation length $r_{0}=6.7^{+0.9}_{-1.0}h^{-1}$cMpc. The derived galaxy bias $b=8.6^{+0.9}_{-1.0}$ is that of dark-matter halos of $M=10^{11.1^{+0.2}_{-0.3}}$M$_{\odot}$ at $z=7.2$, and highlights that galaxies below the current detection limit ($M_{AB}\sim-17.7$) are expected in lower-mass halos ($M\sim10^{8}-10^{10.5}$M$_{\odot}$). We compute the ACF of LBGs at $z\sim3.8-z\sim5.9$ in the same surveys. A trend of increasing bias is found from $z=3.8$ ($b\sim3.0$) to $z=7.2$ ($b\sim8.6$), broadly consistent with galaxies at fixed luminosity being hosted in dark-matter halos of similar mass at $4<z<6$, followed by a slight rise in halo masses at $z\sim7$ ($\sim2\sigma$ confidence). Separating the data at the median luminosity of the $z=7.2$ sample ($M_{UV}=-19.4$) shows higher clustering at $z=5.9$ for bright galaxies ($r_{0}=5.5^{+1.4}_{-1.5}h^{-1}$cMpc, $b=6.2^{+1.2}_{-1.5}$) compared to faint galaxies ($r_{0}=1.9^{+1.1}_{-1.0}h^{-1}$cMpc, $b=2.7\pm1.2$) implying a constant mass-to-light ratio $\frac{dlogM}{dlogL}\sim1.2^{+1.8}_{-0.8}$. A similar trend is present in the $z=7.2$ sample with larger uncertainty. Finally, our bias measurements allow us to investigate the fraction of dark-matter halos hosting UV-bright galaxies (the duty-cycle, $\epsilon_{DC}$). At $z=7.2$ values near unity are preferred, which may be explained by the shortened halo assembly time at high-redshift.

We present Spitzer/IRAC observations of nine $z'$-band dropouts highly magnified (2<mu<12) by the Bullet Cluster. We combine archival imaging with our Exploratory program (SURFS UP), which results in a total integration time of ~30 hr per IRAC band. We detect (>3sigma) in both IRAC bands the brightest of these high-redshift galaxies, with [3.6]=23.80+-0.28 mag, [4.5]=23.78+-0.25 mag, and (H-[3.6])=1.17+-0.32 mag. The remaining eight galaxies are undetected to [3.6]~26.4 mag and [4.5]~26.0 mag with stellar masses of ~5x10^7 M_sol. The detected galaxy has an estimated magnification of mu=12+-4, which implies this galaxy has an ultraviolet luminosity of L_1500~0.3 L*_{z=7} --- the lowest luminosity individual source detected in IRAC at z>7. By modeling the broadband photometry, we estimate the galaxy has an intrinsic star-formation rate of SFR~1.3 M_sol/yr and stellar mass of M~2x10^9 M_sol, which gives a specific star-formation rate of sSFR~0.7 Gyr^-1. If this galaxy had sustained this star-formation rate since z~20, it could have formed the observed stellar mass (to within a factor of ~2), we also discuss alternate star-formation histories and argue the exponentially-increasing model is unlikely. Finally, based on the intrinsic star-formation rate, we estimate this galaxy has a likely [C II] flux of <f_[C II]> = 10^{-17} erg/s/cm2.

[Abbreviated] We present the largest search to date for $z\sim8$ Lyman break galaxies (LBGs) based on 350 arcmin$^2$ of HST observations in the V-, Y-, J- and H-bands from the Brightest of Reionizing Galaxies (BoRG) survey. The BoRG dataset includes $\sim$50 arcmin$^2$ of new data and deeper observations of two previous BoRG pointings, from which we present 9 new $z\sim8$ LBG candidates, bringing the total number of BoRG LBGs to 38 with $25.5\leqslant m_{J} \leqslant 27.6$ (AB system). We introduce a new Bayesian formalism for estimating the galaxy luminosity function (LF), which does not require binning (and thus smearing) of the data and includes a likelihood based on the formally correct binomial distribution as opposed to the often used approximate Poisson distribution. We demonstrate the utility of the new method on a sample of $97$ LBGs that combines the bright BoRG galaxies with the fainter sources published in Bouwens et al. (2012) from the HUDF and ERS programs. We show that the $z\sim8$ LF is well described by a Schechter function with a characteristic magnitude $M^\star = -20.15^{+0.29}_{-0.38}$, a faint-end slope of $\alpha = -1.87^{+0.26}_{-0.26}$, and a number density of $\log_{10} \phi^\star [\textrm{Mpc}^{-3}] = -3.24^{+0.25}_{-0.24}$. Integrated down to $M=-17.7$ this LF yields a luminosity density, $\log_{10} \epsilon [\textrm{erg}/\textrm{s/Hz/Mpc}^{3}] = 25.52^{+0.05}_{-0.05}$. Our LF analysis is consistent with previously published determinations within 1$\sigma$. We discuss the implication of our study for the physics of reionization. By assuming theoretically motivated priors on the clumping factor and the photon escape fraction we show that the UV LF from galaxy samples down to $M=-17.7$ can ionize only 10-50% of the neutral hydrogen at $z\sim8$. Full reionization would require extending the LF down to $M=-15$.

Early-type galaxies (ETGs) are observed to be more compact, on average, at $z \gtrsim 2$ than at $z\simeq 0$, at fixed stellar mass. Recent observational works suggest that such size evolution could reflect the similar evolution of the host dark matter halo density as a function of the time of galaxy quenching. We explore this hypothesis by studying the distribution of halo central velocity dispersion ($\sigma_0$) and half-mass radius ($r_{\rm h}$) as functions of halo mass $M$ and redshift $z$, in a cosmological $\Lambda$-CDM $N$-body simulation. In the range $0\lesssim z\lesssim 2.5$, we find $\sigma_0\propto M^{0.31-0.37}$ and $r_{\rm h}\propto M^{0.28-0.32}$, close to the values expected for homologous virialized systems. At fixed $M$ in the range $10^{11} M_\odot \lesssim M\lesssim 5.5 \times 10^{14} M_\odot$ we find $\sigma_0\propto(1+z)^{0.35}$ and $r_{\rm h}\propto(1+z)^{-0.7}$. We show that such evolution of the halo scaling laws is driven by individual haloes growing in mass following the evolutionary tracks $\sigma_0\propto M^{0.2}$ and $r_{\rm h}\propto M^{0.6}$, consistent with simple dissipationless merging models in which the encounter orbital energy is accounted for. We compare the $N$-body data with ETGs observed at $0\lesssim z\lesssim3$ by populating the haloes with a stellar component under simple but justified assumptions: the resulting galaxies evolve consistently with the observed ETGs up to $z \simeq 2$, but the model has difficulty reproducing the fast evolution observed at $z\gtrsim 2$. We conclude that a substantial fraction of the size evolution of ETGs can be ascribed to a systematic dependence on redshift of the dark matter haloes structural properties.

The Grism Lens-Amplified Survey from Space (GLASS) is a Hubble Space Telescope (HST) Large Program, which will obtain 140 orbits of grism spectroscopy of the core and infall regions of 10 galaxy clusters, selected to be among the very best cosmic telescopes. Extensive HST imaging is available from many sources including the CLASH and Frontier Fields programs. We introduce the survey by analyzing spectra of faint multiply-imaged galaxies and $z\gtrsim6$ galaxy candidates obtained from the first seven orbits out of fourteen targeting the core of the Frontier Fields cluster MACS0717.5+3745. Using the G102 and G141 grisms to cover the wavelength range 0.8-1.7$\mu$m, we confirm 4 strongly lensed systems by detecting emission lines in each of the images. For the 9 $z\gtrsim6$ galaxy candidates clear from contamination, we do not detect any emission lines down to a seven-orbit 1$\sigma$ noise level of $\sim$5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. Taking lensing magnification into account, our flux sensitivity reaches $\sim$0.2-5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. These limits over an uninterrupted wavelength range rule out the possibility that the high-$z$ galaxy candidates are instead strong line emitters at lower redshift. These results show that by means of careful modeling of the background - and with the assistance of lensing magnification - interesting flux limits can be reached for large numbers of objects, avoiding pre-selection and the wavelength restrictions inherent to ground-based multi-slit spectroscopy. These observations confirm the power of slitless HST spectroscopy even in fields as crowded as a cluster core.

Population III supernovae have been the focus of growing attention because of their potential to directly probe the properties of the first stars, particularly the most energetic events that can be seen at the edge of the observable universe. But until now pair-pulsation supernovae, in which explosive thermonuclear burning in massive stars fails to unbind them but can eject their outer layers into space, have been overlooked as cosmic beacons at the earliest redshifts. These shells can later collide and, like Type IIn supernovae, produce superluminous events in the UV at high redshifts that could be detected in the near infrared today. We present numerical simulations of a 110 M$_{\odot}$ pair-pulsation explosion done with the Los Alamos radiation hydrodynamics code RAGE. We find that collisions between consecutive pair pulsations are visible in the near infrared out to z $\sim$ 15 - 20 and can probe the earliest stellar populations at cosmic dawn.

We present an analysis of $\sim$1500 H160-selected photometric galaxies detected to a limiting magnitude of 27.8 in the HUDF, using imaging from the HST WFC3/IR camera in combination with archival UV, optical, and NIR imaging. We fit photometric redshifts and stellar population estimates for all galaxies with well-determined Spitzer IRAC fluxes, allowing for the determination of the cumulative mass function within the range $1<z<6$. By selecting samples of galaxies at a constant cumulative number density, we explore the co-evolution of stellar masses and star formation rates from z$\sim$6. We find a steady increase in the SFRs of galaxies at constant number density from z$\sim$6 to z$\sim$3. The peak epoch of star formation is found to shift to later times for galaxies with increasing number densities, in agreement with the expectations from cosmic downsizing. The observed SFRs can fully account for the mass growth to z$\sim$2 amongst galaxies with cumulative number densities greater than 10$^{-3.5}$ Mpc$^{-3}$. For galaxies with a lower constant number density we find the observed stellar masses are $\sim$3 times greater than that which may be accounted for by the observed star formation alone at late times, implying that growth from mergers plays an important role at $z<2$. We additionally observe a decreasing sSFR, equivalent to approximately one order of magnitude, from z$\sim$6 to z$\sim$2 amongst galaxies with number densities less than 10$^{-3.5}$ Mpc$^{-3}$ along with significant evidence that at any redshift the sSFR is higher for galaxies at higher number density. The combination of these findings can qualitatively explain the previous findings of a sSFR plateau at high redshift. Tracing the evolution of the fraction of quiescent galaxies for samples matched in cumulative number density over this redshift range, we find no unambiguous examples of quiescent galaxies at $z>4$.

Using new ultradeep Spitzer/IRAC photometry from the IRAC Ultradeep Field program (IUDF), we investigate the stellar populations of a sample of 63 Y-dropout galaxy candidates at z~8, only 650Myr after the Big Bang. The sources are selected from HST/ACS+WFC3/IR data over the Hubble Ultra Deep Field (HUDF), two HUDF parallel fields, and wide area data over the CANDELS/GOODS-South. The new Spitzer/IRAC data increase the coverage at 3.6 micron and 4.5 micron to ~120h over the HUDF reaching depths of ~28 (AB,1 sigma). The improved depth and inclusion of brighter candidates result in direct >3 sigma IRAC detections of 20/63 sources, of which 11/63 are detected at > 5 sigma. The average [3.6]-[4.5] colors of IRAC detected galaxies at z~8 are markedly redder than those at z~7, observed only 130Myr later. The simplest explanation is that we witness strong rest-frame optical emission lines (in particular [OIII]4959,5007+Hbeta) moving through the IRAC bandpasses with redshift. Assuming that the average rest-frame spectrum is the same at both z~7 and z~8 we estimate a rest-frame equivalent width of W([OIII]4959,5007+Hbeta) = 670 (+260,-170) Angstrom contributing 0.56 (+0.16,-0.11) mag to the 4.5 micron filter at z~8. The corresponding W(Halpha) = 430 (+160,-110) Angstrom implies an average specific star formation rate of sSFR = 11 (+11,-5) Gyr^-1 and a stellar population age of 100 (+100,-50) Myr. Correcting the spectral energy distribution for the contribution of emission lines lowers the average best-fit stellar masses and mass-to-light ratios by x3, decreasing the integrated stellar mass density to rho*(z=8,MUV<-18)=0.6 (+0.4,-0.3) x 10^6 Msun Mpc^-3.

The existence of supermassive black holes as early as z ~ 7 is one of the great unsolved problems in cosmological structure formation. One leading theory argues that they are born during catastrophic baryon collapse in z ~ 15 protogalaxies in strong Lyman-Werner UV backgrounds. Atomic line cooling in such galaxies fragments baryons into massive clumps that are thought to directly collapse to 10^4 - 10^5 solar-mass black holes. We have now discovered that some of these fragments can instead become supermassive stars that eventually explode as pair-instability supernovae with energies of ~ 10^55 erg, the most energetic explosions in the universe. We have calculated light curves and spectra for supermassive Pop III PI SNe with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible in NIR all-sky surveys by Euclid out to z ~ 10 - 15 and by WFIRST and WISH out to z ~ 15 - 20, perhaps revealing the birthplaces of the first quasars.

The eXtreme Deep Field (XDF) combines data from ten years of observations with the HST Advanced Camera for Surveys (ACS) and the Wide-Field Camera 3 Infra-Red (WFC3/IR) into the deepest image of the sky ever in the optical/near-IR. Since the initial observations on the Hubble Ultra-Deep Field (HUDF) in 2003, numerous surveys and programs, including supernova followup, HUDF09, CANDELS, and HUDF12 have contributed additional imaging data across the HUDF region. Yet these have never been combined and made available as one complete ultra-deep optical and near-infrared image dataset. We do so now for the eXtreme Deep Field (XDF) program. Our new and improved processing techniques provide higher quality reductions of the total dataset. All WFC3 near-IR and optical ACS data sets have been fully combined and accurately matched, resulting in the deepest imaging ever taken at these wavelengths ranging from 29.1 to 30.3 AB mag (5sigma in a 0.35" diameter aperture) in 9 filters. The gains in the optical for the four filters done in the original ACS HUDF correspond to a typical improvement of 0.15 mag, with gains of 0.25 mag in the deepest areas. Such gains are equivalent to adding ~130 to ~240 orbits of ACS data to the HUDF. Improved processing alone results in a typical gain of ~0.1 mag. Our 5sigma (optical+near-IR) SExtractor catalogs reveal about 14140 sources in the full field and about 7121 galaxies in the deepest part of the XDF (the HUDF09 region). The XDF is the deepest image of the universe ever taken, reaching, in the combined image for a flat f_nu source, to 31.2 AB mag 5sigma (32.9 at 1sigma) in a 0.35" diameter aperture.

The first stars are the key to the formation of primitive galaxies, early cosmological reionization and chemical enrichment, and the origin of supermassive black holes. Unfortunately, in spite of their extreme luminosities, individual Population III stars will likely remain beyond the reach of direct observation for decades to come. However, their properties could be revealed by their supernova explosions, which may soon be detected by a new generation of NIR observatories such as JWST and WFIRST. We present light curves and spectra for Pop III pair-instability supernovae calculated with the Los Alamos radiation hydrodynamics code RAGE. Our numerical simulations account for the interaction of the blast with realistic circumstellar envelopes, the opacity of the envelope, and Lyman absorption by the neutral IGM at high redshift, all of which are crucial to computing the NIR signatures of the first cosmic explosions. We find that JWST will detect pair-instability supernovae out to z > 30, WFIRST will detect them in all-sky surveys out to z ~ 15 - 20 and LSST and Pan-STARRS will find them at z ~ 7 - 8. The discovery of these ancient explosions will probe the first stellar populations and reveal the existence of primitive galaxies that might not otherwise have been detected.

We present MOSFIRE spectroscopy of 13 candidate z~8 galaxies selected as Y-dropouts as part of the BoRG pure parallel survey. We detect no significant lya emission (our median 1-sigma rest frame equivalent width sensitivity is in the range 2-16 AA). Using the Bayesian framework derived in a previous paper, we perform a rigorous analysis of a statistical subsample of non-detections for ten Y-dropouts, including data from the literature, to study the cosmic evolution of the lya emission of Lyman Break Galaxies. We find that lya emission is suppressed at z~8 by at least a factor of three with respect to z~6 continuing the downward trend found by previous studies of z-dropouts at z~7. This finding suggests a dramatic evolution in the conditions of the intergalactic or circumgalactic media in just 300 Myrs, consistent with the onset of reionization or changes in the physical conditions of the first generations of starforming regions.

We present a new technique to quantify the light contribution coming from the faint high redshift ($z\sim6$) galaxies below the detection threshold of imaging data, set conventionally at S/N=4.5. We illustrate the technique with an application to Hubble Space Telescope Advanced Camera for Surveys images in the F775W and F850LP filters of the Ultra Deep Field parallel field NICP12. The aim of this analysis is to extend by a few magnitudes the faint end of the luminosity function at $z\sim6$. After masking all the detected sources in the field we apply a Fast Fourier Transform to obtain the spatial power spectrum of the background signal. The power spectrum permits us to separate the background noise signal, the residuals due to the data reduction of the wide field, and the overall signal produced by faint galaxies. The ratio of the signal in the i_775 and z_850 bands is used to estimate the contribution of the faint i-dropout objects. We rely on extensive Monte Carlo simulations to characterize various sources of uncertainty and quantify the number of faint i-dropout galaxies in the field. The analysis allows us to put constraints on the luminosity function at $z\sim6$ down to z_850= 30 mag, 2.5 mag fainter than with standard techniques on the same data. The data are consistent with a faint end slope of the luminosity function of $\alpha = -1.9$. Assuming a specific set of values for the clumping factor, escape fraction, and spectral energy distribution, we find that the $z\sim6$ undetected galaxies down to z_850=30 mag could have driven cosmic reionization.

The detection of Pop III supernovae could directly probe the primordial IMF for the first time, unveiling the properties of the first galaxies, early chemical enrichment and reionization, and the seeds of supermassive black holes. Growing evidence that some Pop III stars were less massive than 100 solar masses may complicate prospects for their detection, because even though they would have been more plentiful they would have died as core-collapse supernovae, with far less luminosity than pair-instability explosions. This picture greatly improves if the SN shock collides with a dense circumstellar shell ejected during a prior violent LBV type eruption. Such collisions can turn even dim SNe into extremely bright ones whose luminosities can rival those of pair-instability SNe. We present simulations of Pop III Type IIn SN light curves and spectra performed with the Los Alamos RAGE and SPECTRUM codes. Taking into account Lyman-alpha absorption in the early universe and cosmological redshifting, we find that 40 solar mass Pop III Type IIn SNe will be visible out to z ~ 20 with JWST and out to z ~ 7 with WFIRST. Thus, even low mass Pop III SNe can be used to probe the primeval universe.

Understanding the properties of Pop III stars is prerequisite to elucidating the nature of primeval galaxies, the chemical enrichment and reionization of the early IGM, and the origin of supermassive black holes. While the primordial IMF remains unknown, recent evidence from numerical simulations and stellar archaeology suggests that some Pop III stars may have had lower masses than previously thought, 15 - 50 \Ms in addition to 50 - 500 \Ms. The detection of Pop III supernovae by JWST, WFIRST or the TMT could directly probe the primordial IMF for the first time. We present numerical simulations of 15 - 40 \Ms Pop III core-collapse SNe done with the Los Alamos radiation hydrodynamics code RAGE. We find that they will be visible in the earliest galaxies out to z ~ 10 - 15, tracing their star formation rates and in some cases revealing their positions on the sky. Since the central engines of Pop III and solar-metallicity core-collapse SNe are quite similar, future detection of any Type II supernovae by next-generation NIR instruments will in general be limited to this epoch.

We present a comprehensive analysis of z>8 galaxies based on ultra-deep WFC3/IR data. We constrain the evolution of the UV luminosity function (LF) and luminosity densities from z~11 to z~8 by exploiting all the WFC3/IR data over the Hubble Ultra-Deep Field from the HUDF09 and the new HUDF12 program, in addition to the HUDF09 parallel field data, as well as wider area WFC3/IR imaging over GOODS-South. Galaxies are selected based on the Lyman Break Technique in three samples centered around z~9, z~10 and z~11, with seven z~9 galaxy candidates, and one each at z~10 and z~11. We confirm a new z~10 candidate (with z=9.8+-0.6) that was not convincingly identified in our first z~10 sample. The deeper data over the HUDF confirms all our previous z>~7.5 candidates as genuine high-redshift candidates, and extends our samples to higher redshift and fainter limits (H_160~29.8 mag). We perform one of the first estimates of the z~9 UV LF and improve our previous constraints at z~10. Extrapolating the lower redshift UV LF evolution should have revealed 17 z~9 and 9 z~10 sources, i.e., a factor ~3x and 9x larger than observed. The inferred star-formation rate density (SFRD) in galaxies above 0.7 M_sun/yr decreases by 0.6+-0.2 dex from z~8 to z~9, in good agreement with previous estimates. The low number of sources found at z>8 is consistent with a very rapid build-up of galaxies across z~10 to z~8. From a combination of all current measurements, we find a best estimate of a factor 10x decrease in the SFRD from z~8 to z~10, following (1+z)^(-11.4+-3.1). Our measurements thus confirm our previous finding of an accelerated evolution beyond z~8, and signify a rapid build-up of galaxies with M_UV<-17.7 within only ~200 Myr from z~10 to z~8, in the heart of cosmic reionization.

The first stars ended the cosmic Dark Ages and created the first heavy elements necessary for the formation of planets and life. The properties of these stars remain uncertain, and it may be decades before individual Pop III stars are directly observed. Their masses, however, can be inferred from their supernova explosions, which may soon be found in both deep-field surveys by JWST and in all-sky surveys by WFIRST. We have performed radiation hydrodynamical simulations of the near infrared signals of Pop III pair-instability supernovae in realistic circumstellar environments with Lyman absorption by the neutral intergalactic medium. We find that JWST and WFIRST will detect these explosions out to z ~ 30 and 20, respectively, unveiling the first generation of stars in the universe.

We report the discovery of 33 Lyman break galaxy (LBG) candidates at z~8 detected in Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging as part of the Brightest of Reionizing Galaxies (BoRG) pure-parallel survey. The ongoing BoRG survey currently has the largest area (274 arcmin^2) with Y_098 (or Y_105), J_125, and H_160 band coverage needed to search for z~8 galaxies, about three times the current CANDELS area, and slightly larger than what will be the final CANDELS wide component with Y_105 data. Our sample of 33 relatively bright Y_098-dropout galaxies have J_125 band magnitudes between 25.5 and 27.4 mag. This is the largest sample of bright (J_125 <~ 27.4) z~8 galaxy candidates presented to date. Combining our dataset with the Hubble Ultra-Deep Field (HUDF09) dataset, we constrain the rest-frame ultraviolet galaxy luminosity function at z~8 over the widest dynamic range currently available. The combined datasets are well fitted by a Schechter function, i.e. \phi(L) = \phi_* (L/L_*)^{\alpha}\ e^{-(L/L_*)}, without evidence for an excess of sources at the bright end. At 68% confidence, we derive \phi_* = (4.3^{+3.5}_{-2.1}) \times 10^{-4} Mpc^{-3}, M_* = -20.26^{+0.29}_{-0.34}, and a very steep faint-end slope \alpha = -1.98^{+0.23}_{-0.22}. While the best-fit parameters still have a strong degeneracy, especially between \phi_* and M_*, our improved coverage at the bright end has reduced the uncertainty of the faint-end power-law slope at z~8 compared to the best previous determination at +/-0.4. With a future expansion of the BoRG survey, combined with planned ultradeep WFC3/IR observations, it will be possible to further reduce this uncertainty and clearly demonstrate the steepening of the faint-end slope compared to measurements at lower redshift, thereby confirming the key role played by small galaxies in the reionization of the universe.

The timing and duration of the reionization epoch is crucial to the emergence and evolution of structure in the universe. The relative roles that star-forming galaxies, active galactic nuclei and quasars play in contributing to the metagalactic ionizing background across cosmic time remains uncertain. Deep quasar counts provide insights into their role, but the potentially crucial contribution from star-formation is highly uncertain due to our poor understanding of the processes that allow ionizing radiation to escape into the intergalactic medium (IGM). The fraction of ionizing photons that escape from star-forming galaxies is a fundamental free parameter used in models to "fine-tune" the timing and duration of the reionization epoch that occurred somewhere between 13.4 and 12.7 Gyrs ago (redshifts between 12 > z > 6). However, direct observation of Lyman continuum (LyC) photons emitted below the rest frame \ion{H}{1} ionization edge at 912 \AA\ is increasingly improbable at redshifts z > 3, due to the steady increase of intervening Lyman limit systems towards high z. Thus UV and U-band optical bandpasses provide the only hope for direct, up close and in depth, observations of the types of environment that favor LyC escape. By quantifying the evolution over the past 11 billion years (z < 3) of the relationships between LyC escape and local and global parameters ..., we can provide definitive information on the LyC escape fraction that is so crucial to answering the question of, how did the universe come to be ionized? Here we provide estimates of the ionizing continuum flux emitted by "characteristic" (L_{uv}^*) star-forming galaxies as a function of look back time and escape fraction, finding that at z = 1 (7.6 Gyrs ago) L_{uv}^* galaxies with an escape fraction of 1% have a flux of 10^{-19} ergs cm^{-2} s^{-1} \AA^{-1}.

We present the first results of our spectroscopic follow-up of 6.5 < z < 10 candidate galaxies behind clusters of galaxies. We report the spectroscopic confirmation of an intrinsically faint Lyman break galaxy (LBG) identified as a z 850LP-band dropout behind the Bullet Cluster. We detect an emission line at {\lambda} = 9412 {\AA} at >5{\sigma} significance using a 16 hr long exposure with FORS2 VLT. Based on the absence of flux in bluer broadband filters, the blue color of the source, and the absence of additional lines, we identify the line as Ly{\alpha} at z = 6.740 \pm 0.003. The integrated line flux is f = (0.7 \pm 0.1 \pm 0.3) \times 10^{-17} erg^{-1} s^{-1} cm^{-2} (the uncertainties are due to random and flux calibration errors, respectively) making it the faintest Ly{\alpha} flux detected at these redshifts. Given the magnification of {\mu} = 3.0 \pm 0.2 the intrinsic (corrected for lensing) flux is f^int = (0.23 \pm 0.03 \pm 0.10 \pm 0.02) \times 10^{-17} erg^{-1} s^{-1} cm^{-2} (additional uncertainty due to magnification), which is ~2-3 times fainter than other such measurements in z ~ 7 galaxies. The intrinsic H 160W-band magnitude of the object is m^int(H_160W)=27.57 \pm 0.17, corresponding to 0.5 L* for LBGs at these redshifts. The galaxy is one of the two sub-L* LBG galaxies spectroscopically confirmed at these high redshifts (the other is also a lensed z = 7.045 galaxy), making it a valuable probe for the neutral hydrogen fraction in the early universe.

The HUDF09 data are the deepest near-IR observations ever, reaching to 29.5 mag. Luminosity functions (LF) from these new HUDF09 data for 132 z\sim7 and z\sim8 galaxies are combined with new LFs for z\sim5-6 galaxies and the earlier z\sim4 LF to reach to very faint limits (<0.05 L*(z=3)). The faint-end slopes alpha are steep: -1.79+/-0.12 (z\sim5), -1.73+/-0.20 (z\sim6), -2.01+/-0.21 (z\sim7), and -1.91+/-0.32 (z\sim8). Slopes alpha\lesssim-2 lead to formally divergent UV fluxes, though galaxies are not expected to form below \sim-10 AB mag. These results have important implications for reionization. The weighted mean slope at z\sim6-8 is -1.87+/-0.13. For such steep slopes, and a faint-end limit of -10 AB mag, galaxies provide a very large UV ionizing photon flux. While current results show that galaxies can reionize the universe by z\sim6, matching the Thomson optical depths is more challenging. Extrapolating the current LF evolution to z>8, taking alpha to be -1.87+/-0.13 (the mean value at z\sim6-8), and adopting typical parameters, we derive Thomson optical depths of 0.061_{-0.006}^{+0.009}. However, this result will change if the faint-end slope alpha is not constant with redshift. We test this hypothesis and find a weak, though uncertain, trend to steeper slopes at earlier times (dalpha/dz\sim-0.05+/-0.04), that would increase the Thomson optical depths to 0.079_{-0.017}^{+0.063}, consistent with recent WMAP estimates (tau=0.088+/-0.015). It may thus not be necessary to resort to extreme assumptions about the escape fraction or clumping factor. Nevertheless, the uncertainties remain large. Deeper WFC3/IR+ACS observations can further constrain the ionizing flux from galaxies.

We use imaging obtained with the Hubble Space Telescope Wide Field Camera 3 to search for z_850 dropouts at z~7 and J_110 dropouts at z~9 lensed by the Bullet Cluster. In total we find 10 z_850 dropouts in our 8.27 arcmin^2 field. Using magnification maps from a combined weak and strong lensing mass reconstruction of the Bullet Cluster and correcting for estimated completeness levels, we calculate the surface density and luminosity function of our z_850 dropouts as a function of intrinsic (accounting for magnification) magnitude. We find results consistent with published blank field surveys, despite using much shallower data, and demonstrate the effectiveness of cluster surveys in the search for z~7 galaxies.

Spectroscopic confirmation of galaxies at z~7 and above has been extremely difficult, owing to a drop in intensity of Ly-alpha emission in comparison with samples at z~6. This crucial finding could potentially signal the ending of cosmic reionization. However it is based on small datasets, often incomplete and heterogeneous in nature. We introduce a flexible Bayesian framework, useful to interpret such evidence. Within this framework, we implement two simple phenomenological models: a smooth one, where the distribution of Ly-alpha is attenuated by a factor \es with respect to z~6; a patchy one where a fraction \ep is absorbed/non-emitted while the rest is unabsorbed. From a compilation of 39 observed z~7 galaxies we find \es=0.69+-0.12 and \ep=0.66+-0.16. The models can be used to compute fractions of emitters above any equivalent width W. For W>25\AA, we find X^{25}_{z=7}=0.37+-0.11 (0.14+-0.06) for galaxies fainter (brighter) than M_{UV}=-20.25 for the patchy model, consistent with previous work, but with smaller uncertainties by virtue of our full use of the data. At z~8 we combine new deep (5-\sigma flux limit 10^{-17}ergs^{-1}cm^{-2}) Keck-NIRSPEC observations of a bright Y-dropout identified by our BoRG Survey, with those of three objects from the literature and find that the inference is inconclusive. We compute predictions for future near-infrared spectroscopic surveys and show that it is challenging but feasible to constrain the distribution of Ly-alpha emitters at z~8 and distinguish between models.

Theoretical and numerical modeling of dark-matter halo assembly predicts that the most luminous galaxies at high redshift are surrounded by overdensities of fainter companions. We test this prediction with HST observations acquired by our Brightest of Reionizing Galaxies (BoRG) survey, which identified four very bright z~8 candidates as Y-dropout sources in four of the 23 non-contiguous WFC3 fields observed. We extend here the search for Y-dropouts to fainter luminosities (M_* galaxies with M_AB\sim-20), with detections at >5sigma confidence (compared to >8sigma confidence adopted earlier) identifying 17 new candidates. We demonstrate that there is a correlation between number counts of faint and bright Y-dropouts at >99.84% confidence. Field BoRG58, which contains the best bright z\sim8 candidate (M_AB=-21.3), has the most significant overdensity of faint Y-dropouts. Four new sources are located within 70arcsec (corresponding to 3.1 comoving Mpc at z=8) from the previously known brighter z\sim8 candidate. The overdensity of Y-dropouts in this field has a physical origin to high confidence (p>99.975%), independent of completeness and contamination rate of the Y-dropout selection. We modeled the overdensity by means of cosmological simulations and estimate that the principal dark matter halo has mass M_h\sim(4-7)x10^11Msun (\sim5sigma density peak) and is surrounded by several M_h\sim10^11Msun halos which could host the fainter dropouts. In this scenario, we predict that all halos will eventually merge into a M_h>2x10^14Msun galaxy cluster by z=0. Follow-up observations with ground and space based telescopes are required to secure the z\sim8 nature of the overdensity, discover new members, and measure their precise redshift.

We identify 73 z~7 and 59 z~8 candidate galaxies in the reionization epoch, and use this large 26-29.4 AB mag sample of galaxies to derive very deep luminosity functions to <-18 AB mag and the star formation rate density at z~7 and z~8. The galaxy sample is derived using a sophisticated Lyman-Break technique on the full two-year WFC3/IR and ACS data available over the HUDF09 (~29.4 AB mag, 5 sigma), two nearby HUDF09 fields (~29 AB mag, 14 arcmin) and the wider area ERS (~27.5 AB mag) ~40 arcmin**2). The application of strict optical non-detection criteria ensures the contamination fraction is kept low (just ~7% in the HUDF). This very low value includes a full assessment of the contamination from lower redshift sources, photometric scatter, AGN, spurious sources, low mass stars, and transients (e.g., SNe). From careful modelling of the selection volumes for each of our search fields we derive luminosity functions for galaxies at z~7 and z~8 to <-18 AB mag. The faint-end slopes alpha at z~7 and z~8 are uncertain but very steep at alpha = -2.01+/-0.21 and alpha=-1.91+/-0.32, respectively. Such steep slopes contrast to the local alpha<~-1.4 and may even be steeper than that at z~4 where alpha=-1.73+/-0.05. With such steep slopes (alpha<~-1.7) lower luminosity galaxies dominate the galaxy luminosity density during the epoch of reionization. The star formation rate densities derived from these new z~7 and z~8 luminosity functions are consistent with the trends found at later times (lower redshifts). We find reasonable consistency, with the SFR densities implied from reported stellar mass densities, being only ~40% higher at z<7. This suggests that (1) the stellar mass densities inferred from the Spitzer IRAC photometry are reasonably accurate and (2) that the IMF at very high redshift may not be very different from that at later times.

In this paper, we present a derivation of the rest-frame 1400A luminosity function (LF) at redshift six from a new application of the maximum likelihood method by exploring the five deepest HST/ACS fields, i.e., the HUDF, two UDF05 fields, and two GOODS fields. We work on the latest improved data products, which makes our results more robust than those of previous studies. We use un-binned data and thereby make optimal use of the information contained in the dataset. We focus on the analysis to a magnitude limit where the completeness is larger than 50% to avoid possibly large errors in the faint end slope that are difficult to quantify. We also take into account scattering in and out of the dropout sample due to photometric errors by defining for each object a probability that it belongs to the dropout sample. We find the best fit Schechter parameters to the z~6 LF are: alpha = 1.87 +/- 0.14, M* = -20.25 +/- 0.23, and phi*=1.77^{+0.62}_{-0.49} * 10^{-3} Mpc^{-3}. Such a steep slope suggests that galaxies, especially the faint ones, are possibly the main sources of ionizing photons in the universe at redshift six. We also combine results from all studies at z~6 to reach an agreement in 95% confidence level that -20.45<M*<-20.05 and -1.90<alpha<-1.55. The luminosity density has been found not to evolve significantly between z~6 and z~5, but considerable evolution is detected from z~6 to z~3.

We report the discovery of a Compton-thick Active Galactic Nucleus (AGN) at z=4.76 in the 4 Ms Chandra Deep Field South. This object was selected as a V-band dropout in HST/ACS images and previously recognized as an AGN from optical spectroscopy. The 4 Ms Chandra observations show a significant (~4.2sigma) X-ray detection at the V-band dropout position. The X-ray source displays a hardness ratio of HR=0.23+-0.24, which, for a source at z~5, is highly suggestive of Compton-thick absorption. The source X-ray spectrum is seen above the background level in the energy range of ~0.9-4 keV, i.e., in the rest-frame energy range of ~5-23 keV. When fixing the photon index to Gamma=1.8, the measured column density is N_H=1.4^{+0.9}_{-0.5} x 10^{24} cm^{-2}, which is Compton-thick. To our knowledge, this is the most distant heavily obscured AGN, confirmed by X-ray spectral analysis, discovered so far. The intrinsic (de-absorbed), rest-frame luminosity in the 2-10 keV band is ~2.5 x 10^{44} erg s^{-1}, which places this object among type-2 quasars. The Spectral Energy Distribution shows that massive star formation is associated with obscured black hole accretion. This system may have then been caught during a major co-eval episode of black hole and stellar mass assembly at early times. The measure of the number density of heavily obscured AGN at high redshifts will be crucial to reconstruct the BH/galaxy evolution history from the beginning.

Searches for very-high-redshift galaxies over the past decade have yielded a large sample of more than 6,000 galaxies existing just 900-2,000 million years (Myr) after the Big Bang (redshifts 6 > z > 3; ref. 1). The Hubble Ultra Deep Field (HUDF09) data have yielded the first reliable detections of z ~ 8 galaxies that, together with reports of a gamma-ray burst at z ~ 8.2 (refs 10, 11), constitute the earliest objects reliably reported to date. Observations of z ~ 7-8 galaxies suggest substantial star formation at z > 9-10. Here we use the full two-year HUDF09 data to conduct an ultra-deep search for z ~ 10 galaxies in the heart of the reionization epoch, only 500 Myr after the Big Bang. Not only do we find one possible z ~ 10 galaxy candidate, but we show that, regardless of source detections, the star formation rate density is much smaller (~10%) at this time than it is just ~200 Myr later at z ~ 8. This demonstrates how rapid galaxy build-up was at z ~ 10, as galaxies increased in both luminosity density and volume density from z ~ 8 to z ~ 10. The 100-200 Myr before z ~ 10 is clearly a crucial phase in the assembly of the earliest galaxies.

We present the first results on the search for very bright (M_AB -21) galaxies at redshift z~8 from the Brightest of Reionizing Galaxies (BoRG) survey. BoRG is a Hubble Space Telescope Wide Field Camera 3 pure-parallel survey that is obtaining images on random lines of sight at high Galactic latitudes in four filters (F606W, F098M, F125W, F160W), with integration times optimized to identify galaxies at z>7.5 as F098M-dropouts. We discuss here results from a search area of approximately 130 arcmin^2 over 23 BoRG fields, complemented by six other pure-parallel WFC3 fields with similar filters. This new search area is more than two times wider than previous WFC3 observations at z~8. We identify four F098M-dropout candidates with high statistical confidence (detected at greater than 8sigma confidence in F125W). These sources are among the brightest candidates currently known at z~8 and approximately ten times brighter than the z=8.56 galaxy UDFy-38135539. They thus represent ideal targets for spectroscopic followup observations and could potentially lead to a redshift record, as our color selection includes objects up to z~9. However, the expected contamination rate of our sample is about 30% higher than typical searches for dropout galaxies in legacy fields, such as the GOODS and HUDF, where deeper data and additional optical filters are available to reject contaminants.

We present UV luminosity functions (LFs) at 1500 Angstrom derived from the HST Early Release Science WFC3/UVIS data acquired over ~50 arcmin^2 of the GOODS-South field. The LFs are determined over the entire redshift range z=0.75-2.5 using two methods, similar to those used at higher redshifts for Lyman Break Galaxies (LBGs): (1) 13-band UV+optical+NIR photometric redshifts to study galaxies in the range z=0.5-2 in three bins of dz=0.5, and (2) dropout samples in three redshift windows centered at z~1.5, z~1.9, and z~2.5. The characteristic luminosity dims by ~1.5 mag from z=2.5 to z=0.75, consistent with earlier work. However, the other Schechter function parameters, the faint-end slope and the number density, are found to be remarkably constant over the range z=0.75-2.5. Using these LF determinations we find the UV luminosity density to increase by ~1.4 dex according to (1+z)^{2.58+-0.15} from z~0 to its peak at z~2.5. Strikingly, the inferred faint-end slopes for our LFs are all steeper than alpha=-1.5, in agreement with higher-redshift LBG studies. Since the faint-end slope in the local universe is found to be much flatter with alpha~=-1.2, this poses the question as to when and how the expected flattening occurs. Despite relatively large uncertainties, our data suggest alpha~=-1.7 at least down to z~1. These new results from such a shallow early dataset demonstrate very clearly the remarkable potential of WFC3/UVIS for the thorough characterization of galaxy evolution over the full redshift range z~0.5 to z~3.

We study the clustering properties of the first galaxies formed in the Universe. We find that, due to chemical enrichment of the inter-stellar medium by isolated Population III stars formed in mini-halos at redshift z>30, the (chronologically) first galaxies are composed of metal-poor Population II stars and are highly clustered on small scales. In contrast, chemically pristine galaxies in halos with mass M~10^8 M_sun may form at z<20 in relatively underdense regions of the Universe. This occurs once self-enrichment by Population III in mini-halos is quenched by the build-up of an $H_2$ photo-dissociating radiative background in the Lyman-Werner bands. We find that these chemically pristine galaxies are spatially uncorrelated. Thus, we expect that deep fields with the James Webb Space Telescope may detect clusters of chemically enriched galaxies but individual chemically pristine objects. We predict that metal-free galaxies at 10 <= z <= 15$ have surface densities of about 80 per square arcmin and per unit redshift but most of them will be too faint even for JWST. However, the predicted density makes these objects interesting targets for searches behind lensing clusters.

We investigate the Spitzer/IRAC properties of 36 z~7 z-dropout galaxies and 3 z~8 Y-dropout galaxies derived from deep/wide-area WFC3/IR data of the Early Release Science, the ultradeep HUDF09, and wide-area NICMOS data. We fit stellar population synthesis models to the SEDs to derive mean redshifts, stellar masses, and ages. The z~7 galaxies are best characterized by substantial ages (>100 Myr) and M/L_V~0.2. The main trend with decreasing luminosity is that of bluing of the far-UV slope from beta~-2.0 to beta~-3.0. This can be explained by decreasing metallicity, except for the lowest luminosity galaxies (0.1 L*{z=3}), where low metallicity and smooth SFHs fail to match the blue far-UV and moderately red H-[3.6] color. Such colors may require episodic SFHs with short periods of activity and quiescence ("on-off" cycles) and/or a contribution from emission lines. The stellar mass of our sample of z~7 star forming galaxies correlates with SFR according to log M* = 8.70 (+-0.09) + 1.06(+-0.10) log SFR, implying star formation may have commenced at z>10. No galaxies are found with SFRs much higher or lower than the past averaged SFR suggesting that the typical star formation timescales are probably a substantial fraction of the Hubble time. We report the first IRAC detection of Y-dropout galaxies at z~8. The average rest-frame U-V=0.3 (AB) of the 3 galaxies are similar to faint z~7 galaxies, implying similar M/L. The stellar mass density to M_{UV,AB}<-18 is rho*(z=8) = 1.8^{+0.7}_{-1.0} x 10^6 M_sun Mpc^-3, following log rho*(z)= 10.6(+-0.6) - 4.4(+-0.7) log(1+z) [M_sun Mpc^-3] over 3<z< 8.

The new Wide Field Camera 3/IR observations on the Hubble Ultra-Deep Field started investigating the properties of galaxies during the reionization epoch. To interpret these observations, we present a novel approach inspired by the conditional luminosity function method. We calibrate our model to observations at z=6 and assume a non-evolving galaxy luminosity versus halo mass relation. We first compare model predictions against the luminosity function measured at z=5 and z=4. We then predict the luminosity function at z>=7 under the sole assumption of evolution in the underlying dark-matter halo mass function. Our model is consistent with the observed z>6.5 galaxy number counts in the HUDF survey and suggests a possible steepening of the faint-end slope of the luminosity function: alpha(z>8)< -1.9 compared to alpha=-1.74 at z=6. Although we currently see only the brightest galaxies, a hidden population of lower luminosity objects (L/L_{*}> 10^{-4}) might provide >75% of the total reionizing flux. Assuming escape fraction f_{esc}~0.2, clumping factor C~5, top heavy-IMF and low metallicity, galaxies below the detection limit produce complete reionization at z>8. For solar metallicity and normal stellar IMF, reionization finishes at z>6, but a smaller C/f_{esc} is required for an optical depth consistent with the WMAP measurement. Our model highlights that the star formation rate in sub-L_* galaxies has a quasi-linear relation to dark-matter halo mass, suggesting that radiative and mechanical feedback were less effective at z>6 than today.

We use the ultra-deep WFC3/IR data over the HUDF and the Early Release Science WFC3/IR data over the CDF-South GOODS field to quantify the broadband spectral properties of candidate star-forming galaxies at z~7. We determine the UV-continuum slope beta in these galaxies, and compare the slopes with galaxies at later times to measure the evolution in beta. For luminous L*(z=3) galaxies, we measure a mean UV-continuum slope beta of -2.0+/-0.2, which is comparable to the beta~-2 derived at similar luminosities at z~5-6. However, for the lower luminosity 0.1L*(z=3) galaxies, we measure a mean beta of -3.0+/-0.2. This is substantially bluer than is found for similar luminosity galaxies at z~4, just 800 Myr later, and even at z~5-6. In principle, the observed beta of -3.0 can be matched by a very young, dust-free stellar population, but when nebular emission is included the expected beta becomes >~-2.7. To produce these very blue beta's (i.e., beta~-3), extremely low metallicities and mechanisms to reduce the red nebular emission are likely required. For example, a large escape fraction (i.e., f_{esc}>~0.3) could minimize the contribution from this red nebular emission. If this is correct and the escape fraction in faint z~7 galaxies is >~0.3, it may help to explain how galaxies reionize the universe.

We utilize the newly-acquired, ultra-deep WFC3/IR observations over the HUDF to search for star-forming galaxies at z~8-8.5, only 600 million years from recombination, using a Y_{105}-dropout selection. The new 4.7 arcmin**2 WFC3/IR observations reach to ~28.8 AB mag (5 sigma) in the Y_{105}J_{125}H_{160} bands. These remarkable data reach ~1.5 AB mag deeper than the previous data over the HUDF, and now are an excellent match to the HUDF optical ACS data. For our search criteria, we use a two-color Lyman-Break selection technique to identify z~8-8.5 Y_{105}-dropouts. We find 5 likely z~8-8.5 candidates. The sources have H_{160}-band magnitudes of ~28.3 AB mag and very blue UV-continuum slopes, with a median estimated beta of <~-2.5 (where f_{\lambda}\propto \lambda^{\beta}). This suggests that z~8 galaxies are not only essentially dust free but also may have very young ages or low metallicities. The observed number of Y_{105}-dropout candidates is smaller than the 20+/-6 sources expected assuming no evolution from z~6, but is consistent with the 5 expected extrapolating the Bouwens et al. 2008 LF results to z~8. These results provide evidence that the evolution in the LF seen from z~7 to z~3 continues to z~8. The remarkable improvement in the sensitivity of WFC3/IR has enabled HST to cross a threshold, revealing star-forming galaxies at z~8-9.

We study the Spitzer Infrared Array Camera (IRAC) mid-infrared (rest-frame optical) fluxes of 14 newly WFC3/IR-detected z=7 z_{850}-dropout galaxies and 5 z=8 Y_{105}-dropout galaxies. The WFC3/IR depth and spatial resolution allow accurate removal of contaminating foreground light, enabling reliable flux measurements at 3.6 micron and 4.5 micron. None of the galaxies are detected to [3.6]=26.9 (AB, 2 sigma), but a stacking analysis reveals a robust detection for the z_{850}-dropouts and an upper limit for the Y_{105}-dropouts. We construct average broadband SEDs using the stacked ACS, WFC3, and IRAC fluxes and fit stellar population synthesis models to derive mean redshifts, stellar masses, and ages. For the z_{850}-dropouts, we find z=6.9^{+0.1}_{-0.1}, (U-V)_{rest}=0.4, reddening A_V=0, stellar mass M*=1.2^{+0.3}_{-0.6} x 10^9 M_sun (Salpeter IMF). The best-fit ages ~300Myr, M/L_V=0.2, and SSFR=1.7Gyr^{-1} are similar to values reported for luminous z=7 galaxies, indicating the galaxies are smaller but not younger. The sub-L* galaxies observed here contribute significantly to the stellar mass density and under favorable conditions may have provided enough photons for sustained reionization at 7<z<11. In contrast, the z=8.3^{+0.1}_{-0.2} Y_{105}-dropouts have stellar masses that are uncertain by 1.5 dex due to the near-complete reliance on far-UV data. Adopting the 2 sigma upper limit on the M/L(z=8), the stellar mass density to M_{UV,AB} < -18 declines from rho*(z=7)=3.7^{+1.0}_{-1.8} x 10^6 M_sun Mpc^{-3} to rho*(z=8) < 8 x 10^5 M_sun Mpc^{-3}, following (1+z)^{-6} over 3<z<8. Lower masses at z=8 would signify more dramatic evolution, which can be established with deeper IRAC observations, long before the arrival of the James Webb Space Telescope.

We present a sample of 16 robust z~7 z-drop galaxies detected by the newly installed WFC3/IR camera on the Hubble Space Telescope. Our analysis is based on the first epoch data of the HUDF09 program covering the Hubble Ultra Deep Field with 60 orbits of Y, J, and H observations. These remarkable data cover 4.7 arcmin^2 and are the deepest NIR images ever taken, reaching to ~29 mag AB (5 sigma). The 16 z~6.5-7.5 galaxies have been identified based on the Lyman Break technique utilizing (z-Y) vs. (Y-J) colors. They have magnitudes J = 26.0-29.0 (AB), an average apparent half-light radius of ~0.16 arcsec (<~1 kpc), and show very blue colors (some even beta<~-2.5), in particular at low luminosities. The WFC3/IR data confirms previous NICMOS detections indicating that the dropout selection at z~7 is very reliable. Our data allow a first determination of the faint end slope of the z~7 luminosity function, reaching down to M_UV ~ -18, a full magnitude fainter than previous measurements. When fixing phi*=1.4e-3/Mpc^3/mag to the value previously measured at z~6, we find a best-fit value of alpha=-1.77+-0.20, with a characteristic luminosity of M_*=-19.91+-0.09. This steep slope is similar to what is seen at z~2-6 and indicates that low luminosity galaxies could potentially provide adequate flux to reionize the universe. The remarkable depth and resolution of these new images provide insights into the coming power of JWST.

We present a first morphological study of z~7-8 Lyman Break galaxies (LBGs) from Oesch et al. 2009 and Bouwens et al. 2009 detected in ultra-deep near-infrared imaging of the Hubble Ultra Deep field (HUDF) by the HUDF09 program. With an average intrinsic size of 0.7+-0.3 kpc these galaxies are found to be extremely compact having an average observed surface brightness of mu_J ~= 26 mag arcsec^(-2), and only two out of the full sample of 16 z~7 galaxies show extended features with resolved double cores. By comparison to lower redshift LBGs it is found that only little size evolution takes place from z~7 to z~6, while galaxies between z~4-5 show more extended wings in their apparent profiles. The average size scales as (1+z)^(-m) with m=1.12+-0.17 for galaxies with luminosities in the range (0.3-1)L*_{z=3} and with m=1.32+-0.52 for (0.12-0.3)L*_{z=3}, consistent with galaxies having constant comoving sizes. The peak of the size distribution changes only slowly from z~7 to z~4. However, a tail of larger galaxies (>~ 1.2 kpc) is gradually built up towards later cosmic times, possibly via hierarchical build-up or via enhanced accretion of cold gas. Additionally, the average star-formation surface density of LBGs with luminosities (0.3-1)L*_{z=3} is nearly constant at Sigma_{SFR}=1.9 Msun/yr/kpc^2 over the entire redshift range z~4-7 suggesting similar star-formation efficiencies at these early epochs. The above evolutionary trends seem to hold out to z~8 though the sample is still small and possibly incomplete.

We present a study of the properties of the star-cluster systems around pseudo-bulges of late-type spiral galaxies using a sample of 11 galaxies with distances from 17 to 37 Mpc. Star clusters are identified from multiband HST ACS and WFPC2 imaging data by combining detections in 3 bands (F435W and F814W with ACS and F606W with WFPC2). The photometric data are then compared to population synthesis models to infer the masses and ages of the star clusters. Photometric errors and completeness are estimated by means of artificial source Monte Carlo simulations. Dust extinction is estimated by considering F160W NICMOS observations of the central regions of the galaxies, augmenting our wavelength coverage. In all galaxies we identify star clusters with a wide range of ages, from young (age < 8 Myr) blue clusters, with typical mass of 10^3 Msun to older (age > 100-250 Myr), more massive, red clusters. Some of the latter might likely evolve into objects similar to the Milky Way's globular clusters. We compute the specific frequencies for the older clusters with respect to the galaxy and bulge luminosities. Specific frequencies relative to the galaxy light appear consistent with the globular cluster specific frequencies of early-type spirals. We compare the specific frequencies relative to the bulge light with the globular cluster specific frequencies of dwarf galaxies, which have a surface-brightness profile that is similar to that of the pseudo-bulges in our sample. The specific frequencies we derive for our sample galaxies are higher than those of the dwarf galaxies, supporting an evolutionary scenario in which some of the dwarf galaxies might be the remnants of harassed late-type spiral galaxies which hosted a pseudo-bulge.

While the average metallicity of the intergalactic medium rises above Z~10^{-3} Zsun by the end of the reionization, pockets of metal-free gas can still exist at later times. We quantify the presence of a long tail in the formation rate of metal-free halos during late stages of reionization (redshift z~6), which might offer the best window to detect Population III stars. Using cosmological simulations for the growth of dark matter halos, coupled with analytical recipes for the metal enrichment of their interstellar medium, we show that pockets of metal-free gas exist at z~6 even under the assumption of high efficiency in metal pollution via winds. A comoving metal-free halo formation rate d^2n/dtdV > 10^{-9} Mpc^{-3}yr^{-1} is expected at z=6 for halos with virial temperature T_{vir}~10^4 K (mass ~10^8 Msun), sufficient to initiate cooling even with strong negative radiative feedback. Under the assumption of a single Population III supernova formed per metal-free halo, we expect an observed supernova rate of 2.6x10^{-3} deg^{-2}yr^{-1} in the same redshift range. These metal-free stars and their supernovae will be isolated and outside galaxies (at distances >150 h^{-1} kpc) and thus significantly less biased than the general population of ~10^8 Msun halos at z~6. Supernova searches for metal-free explosions must thus rely on large area surveys. If metal-free stars produce very luminous supernovae, like SN2006gy, then a multi-epoch survey reaching m_AB =27 at 1 micron is sufficient for detecting them at z=6. While the Large Synoptic Survey Telescope will not reach this depth in the z band, it will be able to detect several tens of Population III supernovae in the i and r bands at z <5.5, when their observed rate is down to 3-8x10^{-4} deg^{-2} yr^{-1}.

[abridged] Theoretical investigations have suggested the presence of Intermediate Mass Black Holes (IMBHs, with masses in the 100-10000 Msun range) in the cores of some Globular Clusters (GCs). In this paper we present the first application of a new technique to determine the presence or absence of a central IMBH in globular clusters that have reached energy equipartition via two-body relaxation. The method is based on the measurement of the radial profile for the average mass of stars in the system, using the fact that a quenching of mass segregation is expected when an IMBH is present. Here we measure the radial profile of mass segregation using main-sequence stars for the globular cluster NGC 2298 from resolved source photometry based on HST-ACS data. The observations are compared to expectations from direct N-body simulations of the dynamics of star clusters with and without an IMBH. The mass segregation profile for NGC 2298 is quantitatively matched to that inferred from simulations without a central massive object over all the radial range probed by the observations, that is from the center to about two half-mass radii. Profiles from simulations containing an IMBH more massive than ~ 300-500 Msun (depending on the assumed total mass of NGC 2298) are instead inconsistent with the data at about 3 sigma confidence, irrespective of the IMF and binary fraction chosen for these runs. While providing a null result in the quest of detecting a central black hole in globular clusters, the data-model comparison carried out here demonstrates the feasibility of the method which can also be applied to other globular clusters with resolved photometry in their cores.

One of the main themes in extragalactic astronomy for the next decade will be the evolution of galaxies over cosmic time. Many future observatories, including JWST, ALMA, GMT, TMT and E-ELT will intensively observe starlight over a broad redshift range, out to the dawn of the modern Universe when the first galaxies formed. It has, however, become clear that the properties and evolution of galaxies are intimately linked to the growth of their central black holes. Understanding the formation of galaxies, and their subsequent evolution, will therefore be incomplete without similarly intensive observations of the accretion light from supermassive black holes (SMBH) in galactic nuclei. To make further progress, we need to chart the formation of typical SMBH at z>6, and their subsequent growth over cosmic time, which is most effectively achieved with X-ray observations. Recent technological developments in X-ray optics and instrumentation now bring this within our grasp, enabling capabilities fully matched to those expected from flagship observatories at longer wavelengths.

[abridged] The First Stars in the Universe form out of pristine primordial gas clouds that have been radiatively cooled to a few hundreds of degrees Kelvin either via molecular or atomic (Lyman-Alpha) hydrogen lines. This primordial mode of star formation is eventually quenched once radiative and/or chemical (metal enrichment) feedbacks mark the transition to Population II stars. In this paper we present a model for the formation rate of Population III stars based on Press-Schechter modeling coupled with analytical recipes for gas cooling and radiative feedback. Our model also includes a novel treatment for metal pollution based on self-enrichment due to a previous episode of Population III star formation in progenitor halos. With this model we derive the star formation history of Population III stars, their contribution to the re-ionization of the Universe and the time of the transition from Population III star formation in minihalos to that in more massive halos where atomic hydrogen cooling is also possible. We consider a grid of models highlighting the impact of varying the values for the free parameters used, such as star formation and feedback efficiency. The most critical factor is the assumption that only one Population III star is formed in a halo. In this scenario, metal free stars contribute only to a minor fraction of the total number of photons required to re-ionize the universe. In addition, metal free star formation is primarily located in minihalos and chemically enriched halos become the dominant locus of star formation very early in the life of the Universe, at redshift z~25. If instead multiple metal free stars are allowed to form out of a single halo, then there is an overall boost of Population III star formation, with a consequent significant contribution to the re-ionizing radiation budget.

We present a sample of $i_{775}$-dropout candidates identified in five Hubble Advanced Camera for Surveys fields centered on Sloan Digital Sky Survey QSOs at redshift $z\sim 6$. Our fields are as deep as the Great Observatory Origins Deep Survey (GOODS) ACS images which are used as a reference field sample. We find them to be overdense in two fields, underdense in two fields, and as dense as the average density of GOODS in one field. The two excess fields show significantly different color distributions from that of GOODS at the 99% confidence level, strengthening the idea that the excess objects are indeed associated with the QSO. The distribution of $i_{775}$-dropout counts in the five fields is broader than that derived from GOODS at the 80% to 96% confidence level, depending on which selection criteria were adopted to identify $i_{775}$-dropouts; its width cannot be explained by cosmic variance alone. Thus, QSOs seem to affect their environments in complex ways. We suggest the picture where the highest redshift QSOs are located in very massive overdensities and are therefore surrounded by an overdensity of lower mass halos. Radiative feedback by the QSO can in some cases prevent halos from becoming galaxies, thereby generating in extreme cases an underdensity of galaxies. The presence of both enhancement and suppression is compatible with the expected differences between lines of sight at the end of reionization as the presence of residual diffuse neutral hydrogen would provide young galaxies with shielding from the radiative effects of the QSO.

[Abridged] We detect three (plus one less certain) z-dropout sources in two separate fields of our UDF05 HST NICMOS images. These z~7 Lyman-Break Galaxy (LBG) candidates allow us to constrain the Luminosity Function (LF) of the star forming galaxy population at those epochs. By assuming a change in only M* and adopting a linear evolution in redshift, anchored to the measured values at z~6, the best fit evolution coefficient is found to be 0.43+-0.19 mag per unit redshift (0.36+-0.18, if including all four candidates), which provides a value of M*(z=7.2)=-19.7+-0.3. This implies a steady evolution for the LBG LF out to z~7, at the same rate that is observed throughout the z~3 to 6 period. This puts a strong constraint on the star-formation histories of z~6 galaxies, whose ensemble star-formation rate density must be lower by a factor 2 at ~170 Myr before the epoch at which they are observed. In particular, a large fraction of stars in the z~6 LBG population must form at redshifts well above z~7. Extrapolating this steady evolution of the LF out to higher redshifts, we estimate that galaxies would be able to reionize the universe by z~6, provided that the faint-end slope of the z>7 LF steepens to alpha~-1.9, and that faint galaxies, with luminosities below the current detection limits, contribute a substantial fraction of the required ionizing photons. This scenario gives however an integrated optical depth to electron scattering that is ~2sigma below the WMAP-5 measurement. Therefore, altogether, our results indicate that, should galaxies be the primary contributors to reionization, either the currently detected evolution of the galaxy population slows down at z>7, or the LF evolution must be compensated by a decrease in metallicity and a corresponding increase in ionization efficiency at these early epochs.

Near-infrared (hereafter NIR) data may provide complementary information to the traditional optical population synthesis analysis of unresolved stellar populations because the spectral energy distribution of the galaxies in the 1-2.5\mum range is dominated by different types of stars than at optical wavelengths. Furthermore, NIR data are subjected to less absorption and hence could constrain the stellar populations in dust-obscured galaxies. We want to develop observational constraints on the stellar populations of unresolved stellar systems in the NIR. To achieve this goal we need a benchmark sample of NIR spectra of ``simple'' early-type galaxies, to be used for testing and calibrating the outputs of population synthesis models. We obtained low-resolution (R~1000) long-slit spectra between 1.5 and 2.4\mum for 14 nearby early-type galaxies using SofI at NTT and higher resolution (R~3000) long-slit spectra, centered at the MgI at ~1.51\mum for a heterogeneous sample of 5 nearby galaxies observed with ISAAC at VLT. We defined spectral indices for CO, NaI, CaI and MgI features and measured the strengths of these features in the sample galaxies. We defined a new global NIR metallicity index, suitable for abundance measurements in low-resolution spectra. Finally, we present an average NIR spectrum of an early-type galaxy, built from a homogenized subset of our sample. The NIR spectra of the sample galaxies show great similarity and the strength of some features does correlate with the iron abundance [Fe/H] and optical metal features of the galaxies. The data suggest that the NIR metal features, in combination with a hydrogen absorption feature may be able to break the age-metallicity degeneracy just like the Mg and Fe features in the optical wavelength range.

A number of recent numerical investigations concluded that the remnants of rare structures formed at very high redshift, such as the very first stars and bright redshift z~6 QSOs, are preferentially located at the center of the most massive galaxy clusters at redshift z=0. In this paper we readdress this question using a combination of cosmological simulations of structure formation and extended Press-Schechter formalism and we show that the typical remnants of Population III stars are instead more likely to be found in a group environment, that is in dark matter halos of mass ~2x10^{13} h^{-1}M_sun. Similarly, the descendants of the brightest z~6 QSOs are expected to be in medium-sized clusters (mass of a few 10^{14} h^{-1}M_sun), rather than in the most massive superclusters (M>10^{15} h^{-1}M_sun) found within the typical 1 Gpc^3 cosmic volume where a bright z~6 QSO lives. The origin of past claims that the most massive clusters preferentially host these remnants is rooted in the numerical method used to initialize their numerical simulations: Only a small region of the cosmological volume of interest was simulated with sufficient resolution to identify low-mass halos at early times, and this region was chosen to host the most massive halo in the cosmological volume at late times. The conclusion that the earliest structures formed in the entire cosmological volume evolve into the most massive halo at late times was thus arrived at by construction. We demonstrate that, to the contrary, the first structures to form in a cosmological region evolve into relatively typical objects at later times. We propose alternative numerical methods for simulating the earliest structures in cosmological volumes.

We explore the design of a space mission, Project Lyman, which has the goal of quantifying the ionization history of the universe from the present epoch to a redshift of z ~ 3. Observations from WMAP and SDSS show that before a redshift of z >~ 6 the first collapsed objects, possibly dwarf galaxies, emitted Lyman continuum (LyC) radiation shortward of 912 A, reionizing most of the universe. How LyC escapes from galactic environments, whether it induces positive or negative feedback on the local and global collapse of structures, and the role played by clumping, molecules, metallicity and dust are major unanswered theoretical questions, requiring observational constraint. Numerous intervening Lyman limit systems, which frustrate the detection of LyC from high z objects, thin below z ~ 3 where there are a few objects with apparently very high fesc. At low z there are only controversial detections and a handful of upper limits. A wide-field multi-object spectroscopic survey with moderate spectral and spatial resolution can quantify fesc within diverse spatially resolved galactic environments over redshifts with significant evolution in galaxy assemblage and quasar activity. It can also calibrate LyC escape against Ly-alpha escape, providing an essential tool to JWST for probing the beginnings of reionization. We present calculations showing the evolution of the characteristic apparent magnitude of star-forming galaxy luminosity functions at 900 A, as a function of redshift and assumed escape fraction to determine the required aperture for detecting LyC. We review our efforts to build a pathfinding dual order multi-object spectro/telescope with a (0.5deg)^2 field-of-view, using a GSFC microshutter array, and crossed delay-line micro-channel plate detector.

In many theoretical scenarios it is expected that intermediate-mass black holes (IMBHs, with masses M ~ 100-10000 solar masses) reside at the centers of some globular clusters. However, observational evidence for their existence is limited. Several previous numerical investigations have focused on the impact of an IMBH on the cluster dynamics or brightness profile. Here we instead present results from a large set of direct N-body simulations including single and binary stars. These show that there is a potentially more detectable IMBH signature, namely on the variation of the average stellar mass between the center and the half-light radius. We find that the existence of an IMBH quenches mass segregation and causes the average mass to exhibit only modest radial variation in collisionally relaxed star clusters. This differs from when there is no IMBH. To measure this observationally requires high resolution imaging at the level of that already available from the Hubble Space Telescope (HST) for the cores of a large sample of galactic globular clusters. With a modest additional investment of HST time to acquire fields around the half-light radius, it will be possible to identify the best candidate clusters to harbor an IMBH. This test can be applied only to globulars with a half-light relaxation time less than or equal to 1 Gyr, which is required to guarantee efficient energy equipartition due to two-body relaxation.

We describe the scientific motivations, the mission concept and the instrumentation of SPACE, a class-M mission proposed for concept study at the first call of the ESA Cosmic-Vision 2015-2025 planning cycle. SPACE aims to produce the largest three-dimensional evolutionary map of the Universe over the past 10 billion years by taking near-IR spectra and measuring redshifts for more than half a billion galaxies at 0<z<2 down to AB~23 over 3\pi sr of the sky. In addition, SPACE will also target a smaller sky field, performing a deep spectroscopic survey of millions of galaxies to AB~26 and at 2<z<10+. These goals are unreachable with ground-based observations due to the ~500 times higher sky background. To achieve the main science objectives, SPACE will use a 1.5m diameter Ritchey-Chretien telescope equipped with a set of arrays of Digital Micro-mirror Devices (DMDs) covering a total field of view of 0.4 deg2, and will perform large-multiplexing multi-object spectroscopy (e.g. ~6000 targets per pointing) at a spectral resolution of R~400 as well as diffraction-limited imaging with continuous coverage from 0.8mum to 1.8mum.

I will focus on the study of the Cosmic Dark ages and in particular on the properties of the galaxies responsible for the reionization of Hydrogen and on the possibility of detecting the first generation of galaxies and the first stars. I will make the case that there is a large and widening gap between our present and future capabilities in imaging and what we can follow-up spectroscopically even with the NIRSpec spectrograph on the James Webb Space Telescope because of its moderate field of view. A high throughput, high multiplexing near-IR spectrometer on the Extremely Large Telescope would be important to contain this gap. Studying the Lyman alpha forest of QSOs at z>=7 and measuring the kinematics of small mass galaxies at high-z are also problems unsuited to JWST and where a high-throughput, medium-high resolution near-IR spectrograph on the ELT would be essential.

We study cosmic variance in deep high redshift surveys and its influence on the determination of the luminosity function for high redshift galaxies. For several survey geometries relevant for HST and JWST instruments, we characterize the distribution of the galaxy number counts. This is obtained by means of analytic estimates via the two point correlation function in extended Press-Schechter theory as well as by using synthetic catalogs extracted from N-body cosmological simulations of structure formation. We adopt a simple luminosity - dark halo mass relation to investigate the environment effects on the fitting of the luminosity function. We show that in addition to variations of the normalization of the luminosity function, a steepening of its slope is also expected in underdense fields, similarly to what is observed within voids in the local universe. Therefore, to avoid introducing artificial biases, caution must be taken when attempting to correct for field underdensity, such as in the case of HST UDF i-dropout sample, which exhibits a deficit of bright counts with respect to the average counts in GOODS. A public version of the cosmic variance calculator based on the two point correlation function integration is made available on the web.

We discuss the link between dark matter halos hosting the first PopIII stars formed at redshift z > 40 and the rare, massive, halos that are generally considered to host bright z~6 quasars. We show that within the typical volume occupied by one bright high-z QSO the remnants of the first several thousands PopIII stars formed do not end up in the most massive halos at z~6, but rather live in a large variety of environments. The black hole seeds planted by these very first PopIII stars can easily grow to M > 10^{9.5} Msun by z=6 assuming Eddington accretion with radiative efficiency epsilon~0.1. Therefore quenching of the accretion is crucial to avoid an overabundance of supermassive black holes. We implement a simple feedback model for the growth of the seeds planted by PopIII stars and obtain a z~6 BH mass function consistent with the observed QSO luminosity function.

We present the UDF05 project, a HST Large Program of deep ACS (F606W, F775W, F850LP) and NICMOS (F110W, F160W) imaging of three fields, two of which coincide with the NICP1-4 NICMOS parallel observations of the Hubble Ultra Deep Field (HUDF). In this first paper we use the ACS data for the NICP12 field, as well as the original HUDF ACS data, to measure the UV Luminosity Function (LF) of z~5 Lyman Break Galaxies (LBGs) down to very faint levels. Specifically, based on a V-i, i-z selection criterion, we identify a sample of 101 and 133 candidate z~5 galaxies down to z_{850}=28.5 and 29.25 magnitudes in the NICP12 and in the HUDF fields, respectively. Using an extensive set of Monte Carlo simulations we derive corrections for observational biases and selection effects, and construct the rest-frame 1400 A LBG LF over the range M_{1400}=[-21.4, -17.1], i.e. down to ~0.04 L* at z~5, and complement it with data from the Subaru Deep Field (SDF) from Yoshida et al. (2006) to extend it to the brighter end (M_{1400}>-22.2). We show that: (i) Different assumptions regarding the SED distribution of the LBG population, dust properties and intergalactic absorption result in a 25% variation in the number density of LBGs at z~5; (ii) Under consistent assumptions for dust properties and intergalactic absorption, the HUDF is about 30% under-dense in z~5 LBGs relative to the NICP12 field, a variation which is well explained by cosmic variance; (iii) The faint-end slope of the LF is independent of the specific assumptions for the input physical parameters, and has a value of alpha ~ -1.6, similar to the faint-end slope of the LF that has been measured for LBGs at z~3 and z~6.

We discuss the link between dark matter halos hosting the first PopIII stars and the rare, massive, halos that are generally considered to host bright quasars at high redshift z~6. The main question that we intend to answer is whether the super-massive black holes powering these QSOs grew out from the seeds planted by the first intermediate massive black holes created in the universe. This question involves a dynamical range of 10^13 in mass and we address it by combining N-body simulations of structure formation to identify the most massive halos at z~6 with a Monte Carlo method based on linear theory to obtain the location and formation times of the first light halos within the whole simulation box. We show that the descendants of the first ~10^6 Msun virialized halos do not, on average, end up in the most massive halos at z~6, but rather live in a large variety of environments. The oldest PopIII progenitors of the most massive halos at z~6, form instead from density peaks that are on average one and a half standard deviations more common than the first PopIII star formed in the volume occupied by one bright high-z QSO. The intermediate mass black hole seeds planted by the very first PopIII stars at z>40 can easily grow to masses m_BH>10^9.5 Msun by z=6 assuming Eddington accretion with radiative efficiency \epsilon~0.1. Quenching of the black hole accretion is therefore crucial to avoid an overabundance of supermassive black holes at lower redshift. This can be obtained if the mass accretion is limited to a fraction \eta~6*10^{-3} of the total baryon mass of the halo hosting the black hole. The resulting high end slope of the black hole mass function at z=6 is \alpha ~ -3.7, a value within the 1\sigma error bar for the bright end slope of the observed quasar luminosity function at z=6.

ABRIDGED: We use HSTACS and NICMOS imaging to study the structure and colors of a sample of nine late-type spirals. We find: (1) A correlation between bulge and disks scale-lengths, and a correlation between the colors of the bulges and those of the inner disks. Our data show a trend for bulges to be more metal-enriched than their surrounding disks, but otherwise no simple age-metallicity connection between these systems; (2) A large range in bulge stellar population properties, and, in particular, in stellar ages. Specifically, in about a half of the late-type bulges in our sample the bulk of the stellar mass was produced recently. Thus, in a substantial fraction of the z=0 disk-dominated bulged galaxies, bulge formation occurs after the formation/accretion of the disk; (3) In about a half of the late-type bulges in our sample, however, the bulk of the stellar mass was produced at early epochs; (4) Even these "old" late-type bulges host a significant fraction of stellar mass in a young(er) c component; (5) A correlation for bulges between stellar age and stellar mass, in the sense that more massive late-type bulges are older than less massive late-type bulges. Since the overall galaxy luminosity (mass) also correlates with the bulge luminosity (mass), it appears that the galaxy mass regulates not only what fraction of itself ends up in the bulge component, but also "when" bulge formation takes place. We show that dynamical friction of massive clumps in gas-rich disks is a plausible disk-driven mode for the formation of "old" late-type bulges. If disk evolutionary processes are responsible for the formation of the entire family of late-type bulges, CDM simulations need to produce a similar number of initially bulgeless disks in addition to the disk galaxies that are observed to be bulgeless at z=0.

This paper presents the Hubble Ultra Deep Field (HUDF), a one million second exposure of an 11 square minute-of-arc region in the southern sky with the Hubble Space Telescope. The exposure time was divided among four filters, F435W (B435), F606W (V606), F775W (i775), and F850LP (z850), to give approximately uniform limiting magnitudes mAB~29 for point sources. The image contains at least 10,000 objects presented here as a catalog. Few if any galaxies at redshifts greater than ~4 resemble present day spiral or elliptical galaxies. Using the Lyman break dropout method, we find 504 B-dropouts, 204 V-dropouts, and 54 i-dropouts. Using these samples that are at different redshifts but derived from the same data, we find no evidence for a change in the characteristic luminosity of galaxies but some evidence for a decrease in their number densities between redshifts of 4 and 7. The ultraviolet luminosity density of these samples is dominated by galaxies fainter than the characteristic luminosity, and the HUDF reveals considerably more luminosity than shallower surveys. The apparent ultraviolet luminosity density of galaxies appears to decrease from redshifts of a few to redshifts greater than 6. The highest redshift samples show that star formation was already vigorous at the earliest epochs that galaxies have been observed, less than one billion years after the Big Bang.

We estimate the optical depth distribution of dust present in absorption systems along the line of sight of high redshift galaxies and the resulting reddening. We characterize the probability distribution of the transmission to a given redshift and the shape of the effective mean extinction law by means of analytical estimates and Monte Carlo simulations. We present our results in a format useful for applications to present samples of high redshift galaxies and discuss the implications for observations with the James Webb Space Telescope. Our most realistic model takes into account the metallicity evolution of Damped Lyman Alpha absorbers and predicts that the effects of dust absorption are modest: at redshift z \gtrsim 5 the transmission is above 0.8 at an emitted wavelength \lambda_e = 0.14 \mu m with probability 90%. Therefore dust obscuration along the line of sight will affect only marginally observations at very high redshift.

We estimate the intrinsic neutral gas density in Damped Lyman Alpha systems ($\Omega_{HI}^{(DLA)}$) in the redshift range $ 2.2 \lesssim z \lesssim 5$ from the DLA SDSS DR_3 sample of optically selected quasars. We take into account self-consistently the obscuration on background quasars due to the dust present in Damped Lyman Alpha systems. We model the column density and redshift distribution of these systems by using both a non-parametric and a parametric approach. Under conservative assumptions on the dust content of Damped Lyman $\alpha$ systems, we show that selection effects lead to underestimating the intrinsic neutral gas density by at least $15\%$ with respect to the observed neutral gas density. Over the redshift range $[2.2;5.5]$ we find $\Omega_{HI}^{(DLA)}=0.97^{+0.08+0.28}_{-0.06-0.15} \cdot 10^{-3}$, where the first set of error bars gives the $1\sigma$ random errors and the second set gives the modeling uncertainty dependent on the fraction of metals in dust - from 0\% to 50\%. This value compares with $\Omega_{HI}^{(DLA)}=0.82^{+0.05}_{-0.05}$ ($1\sigma$ error bars), which is obtained when no correction for dust is introduced. In the model with half of the metals mass in dust we cannot constraint $\Omega_{HI}^{(DLA)}$ at a confidence level higher than $90\%$. In this case there is indeed a probability of about $10\%$ that the intrinsic column density distribution of DLA systems is a power law $f(N_{HI}) \propto 1/N_{HI}^{~1.95}$. In contrast, with $25 \%$ of the metals in dust - the most realistic estimate - a power law is ruled out at $99.5\%$ of confidence level.

The James Webb Space Telescope (JWST) is a large (6.6m), cold (50K), infrared-optimized space observatory that will be launched early in the next decade. The observatory will have four instruments: a near-infrared camera, a near-infrared multi-object spectrograph, and a tunable filter imager will cover the wavelength range, 0.6 to 5.0 microns, while the mid-infrared instrument will do both imaging and spectroscopy from 5.0 to 29 microns. The JWST science goals are divided into four themes. The End of the Dark Ages: First Light and Reionization theme seeks to identify the first luminous sources to form and to determine the ionization history of the early universe. The Assembly of Galaxies theme seeks to determine how galaxies and the dark matter, gas, stars, metals, morphological structures, and active nuclei within them evolved from the epoch of reionization to the present day. The Birth of Stars and Protoplanetary Systems theme seeks to unravel the birth and early evolution of stars, from infall on to dust-enshrouded protostars to the genesis of planetary systems. The Planetary Systems and the Origins of Life theme seeks to determine the physical and chemical properties of planetary systems including our own, and investigate the potential for the origins of life in those systems. To enable these observations, JWST consists of a telescope, an instrument package, a spacecraft and a sunshield. The telescope consists of 18 beryllium segments, some of which are deployed. The segments will be brought into optical alignment on-orbit through a process of periodic wavefront sensing and control. The JWST operations plan is based on that used for previous space observatories, and the majority of JWST observing time will be allocated to the international astronomical community through annual peer-reviewed proposal opportunities.

We present results from a search for high-redshift J--band ``dropout'' galaxies in the portion of the GOODS southern field that is covered by extremely deep imaging from the Hubble Ultradeep Field (HUDF).Using observations at optical, near-infrared and mid-infrared wavelengths from the Hubble and Spitzer Space Telescopes and the ESO-VLT, we search for very massive galaxies at high redshifts and find one particularly remarkable candidate. Its spectral energy distribution is consistent with a galaxy at z ~ 6.5 and a stellar mass of 6x10e11 M(sun) (for a Salpeter IMF). We interpret a prominent photometric break between the near-infrared and Spitzer bandpasses as the 3646A Balmer discontinuity. The best-fitting models have low reddening and ages of several hundred Myr, placing the formation of the bulk of the stars at z > 9. Alternative models of dusty galaxies at z ~ 2.5 are possible but provide significantly poorer fits. The object is detected with Spitzer at 24 micron. This emission originats from an obscured active nucleus or star formation. We present optical and near-infrared spectroscopy which has, thus far, failed to detect any spectral features. This helps limit the solution in which the galaxy is a starburst or active galaxy at z ~ 2.5. If the high-redshift interpretation is correct, this object would be an example of a galaxy that formed by a process strongly resembling traditional models of monolithic collapse, in a way which a very large mass of stars formed within a remarkably short period of time, at very high redshift.

We investigate the central regions of 23 spiral galaxies using archival NICMOS imaging and STIS spectroscopy. The sample is taken from our program to determine the masses of central massive black holes (MBH) in 54 nearby spiral galaxies. Stars are likely to contribute significantly to any dynamical central mass concentration that we find in our MBH program and this paper is part of a series to investigate the nuclear properties of these galaxies. We use the Nuker law to fit surface brightness profiles, derived from the NICMOS images, to look for nuclear star clusters and find possible extended sources in 3 of the 23 galaxies studied (13 per cent). The fact that this fraction is lower than that inferred from optical Hubble Space Telescope studies is probably due to the greater spatial resolution of those studies. Using R-H and J-H colors and equivalent widths of H-alpha emission (from the STIS spectra) we investigate the nature of the stellar population with evolutionary models. Under the assumption of hot stars ionizing the gas, as opposed to a weak AGN, we find that there are young stellar populations (~10-20 Myr) however these data do not allow us to determine what percentage of the total nuclear stellar population they form. Also, in an attempt to find any unknown AGN we use [N II] and [S II] line flux ratios (relative to H-alpha) and find tentative evidence for weak AGN in NGC 1300 and NGC 4536.

We report on a complete sample of 7 luminous early-type galaxies in the Hubble Ultra Deep Field (UDF) with spectroscopic redshifts between 1.39 and 2.47 and to K<23 AB. Using the BzK selection criterion we have pre-selected a set of objects over the UDF which fulfill the photometric conditions for being passively evolving galaxies at z>1.4. Low-resolution spectra of these objects have been extracted from the HST+ACS grism data taken over the UDF by the GRAPES project. Redshift for the 7 galaxies have been identified based on the UV feature at rest frame 2640<lambda<2850 AA. This feature is mainly due to a combination of FeII, MgI and MgII absorptions which are characteristic of stellar populations dominated by stars older than about 0.5 Gyr. The redshift identification and the passively evolving nature of these galaxies is further supported by the photometric redshifts and by the overall spectral energy distribution (SED), with the ultradeep HST+ACS/NICMOS imaging revealing compact morphologies typical of elliptical/early-type galaxies. From the SED we derive stellar masses of 10^{11}Msun or larger and ages of about 1 Gyr. Their space density at < z >=1.7 appears to be roughly a factor of 2--3 smaller than that of their local counterparts, further supporting the notion that such massive and old galaxies are already ubiquitous at early cosmic times. Much smaller effective radii are derived for some of the objects compared to local massive ellipticals, which may be due to morphological K corrections, evolution, or the presence of a central point-like source. Nuclear activity is indeed present in a subset of the galaxies, as revealed by them being hard X-ray sources, hinting to AGN activity having played a role in discontinuing star formation.

We present grism spectra taken with the Advanced Camera for Surveys to identify 29 red sources with (i'-z') > 0.9 in the Hubble Ultra Deep Field (HUDF). Of these 23 are found to be galaxies at redshifts between z=5.4 and 6.7, identified by the break at 1216 angstrom due to IGM absorption; two are late type dwarf stars with red colors; and four are galaxies with colors and spectral shape similar to dust reddened or old galaxies at redshifts z~1-2. This constitutes the largest uniform, flux-limited sample of spectroscopically confirmed galaxies at such faint fluxes z'<27.5. Many are also among the most distant spectroscopically confirmed galaxies (at redshifts up to z=6.7). We find a significant overdensity of galaxies at redshifts z=5.9 +/- 0.2. Nearly two thirds of the galaxies in our sample (15/23) belong to this peak. Taking into account the selection function and the redshift sensitivity of the survey, we get a conservative overdensity of at least a factor of two along the line-of-sight. The galaxies found in this redshift peak are also localized in the plane of the sky in a non-random manner, occupying about half of the ACS chip. Thus the volume overdensity is a factor of four. The star-formation rate derived from detected sources in this overdense region is sufficient to reionize the local IGM.

We present Space Telescope Imaging Spectrograph emission-line spectra of the central regions of the spiral galaxies NGC 1300 and NGC 2748. From the derived kinematics of the nuclear gas we have found evidence for central supermassive black holes in both galaxies. The estimated mass of the black hole in NGC 1300 is 6.6 (+6.3, -3.2) x 10^7 solar masses and in NGC 2748 is 4.4 (+3.5, -3.6) x 10^7 solar masses (both at the 95% confidence level). These two black hole mass estimates contribute to the poorly sampled low-mass end of the nuclear black hole mass spectrum.

We present tentative evidence for primordial clustering, manifested as an excess of color-selected objects in the field of the QSO SDSS J1030+0524 at redshift z=6.28. We have selected objects red in i_{775}-z_{850} on the basis of Hubble Space Telescope Advanced Camera for Surveys imaging of a field centered on the QSO. Compared to data at comparable depth obtained by the GOODS survey, we find an excess of objects with (i_{775}-z_{850}) \geq 1.5 in the QSO field. The significance of the detection is estimated to be ~97% on the basis of the counts alone and increases to 99.4% if one takes into account the color distribution. If confirmed this would represent the highest redshift example of galaxy clustering and would have implications on models for the growth of structure. Bias-driven clustering of first luminous objects forming in the highest peaks of the primordial density field is expected in most models of early structure formation. The redshift of one of the candidates has been found to be z=5.970 by our spectroscopy with Keck I/LRIS, confirming the validity of our color selection.

We have discovered an extended Lyman alpha plume associated with a compact source at redshift 5.4 in slitless spectroscopic data from the Grism ACS Program for Extragalactic Science (GRAPES) project. The spatial extent of the emission is about 6 x 1.5 kpc (1 x 0.25 arcsec). Combining our grism data and the broadband images from the Hubble UltraDeep Field (UDF) images, we find a Lyman alpha line flux of 2e-17 erg/cm2/s and surface brightness 7e-17 erg/cm2/s/arcsec2. The UDF images show diffuse continuum emission associated with the Lyman alpha plume (hereafter UDF 5225), with three embedded knots. The morphology of UDF 5225 is highly suggestive of a galaxy in assembly. It is moreover possible that the prominent Lyman alpha emission from this object is due to an active nucleus, and that we are seeing the simultaneous growth through accretion of a galaxy and its central black hole. Followup observations at higher spectral resolution could test this hypothesis.

Using N-body simulations we show that the process of formation of the brightest cluster galaxy through dissipationless galactic cannibalism can affect the inner cluster dark matter density profile. In particular, we use as realistic test case the dynamical evolution of the galaxy cluster C0337-2522 at redshift z=0.59, hosting in its centre a group of five elliptical galaxies which are likely to be the progenitor of a central giant elliptical. After the formation of the brightest cluster galaxy, the inner cluster dark matter density profile is significantly flatter (logarithmic slope 0.49<beta<0.90) than the original cusp (beta=1), as a consequence of dynamical friction heating of the massive galaxies against the diffuse cluster dark matter. In our simulations we have assumed that the cluster galaxies are made of stars only. We also show that the presence of galactic dark matter haloes can steepen the cluster central density profile. We conclude that galactic cannibalism could be a viable physical mechanism to reconcile - at least at the cluster scale - the flat dark matter haloes inferred observationally in some galaxy clusters with the steep haloes predicted by cosmological simulations.

We present surface photometry for the central regions of a sample of 48 spiral galaxies (mostly unbarred and barred of types Sbc or Sc) observed with the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. Surface brightness profiles were derived and modeled with a Nuker law. We also analyzed archival Wide Field Planetary Camera~2 images with a larger field of view, available for 18 galaxies in our sample. We modeled the extracted bulge surface brightness profiles with an exponential, a de Vaucouleurs or a Sersic profile. In agreement with previous studies, we find that bulges of Sbc galaxies fall into two categories: bulges well described by an exponential profile and those well described by an de Vaucouleurs profile. Only one galaxy requires the use of a more general Sersic profile to properly describe the bulge. Nuclear photometrically distinct components are found in ~55% of the galaxies. For those that we classify as star clusters based on their resolved extent we find absolute magnitudes that are brighter on average than those previously identified in spiral galaxies. This might be due to a bias in our sample toward star forming galaxies, combined with a trend for star forming galaxies to host brighter central clusters.

We compare the available catalogs of $z\approx6$ galaxies in the Hubble Ultra-Deep Field (UDF) and in the Great Observatories Origins Deep Survey (GOODS) with the expected properties of the sources of cosmological reionization from our previous theoretical study. Our approach is based on the mean surface brightness of the sources required for reionization and depends on relatively few undetermined parameters. We find that the observed mean surface brightness of galaxies at $z \approx 6$ is sufficient for reionization, provided that the sources are composed of hot metal-free or metal-poor stars, regardless of whether reionization occurs over a short or long interval of redshift. The broad agreement between the new observations and our predictions suggests that we may have detected the sources responsible for some or even all of the reionization of hydrogen.

We present deep unbiased spectroscopy of the Hubble Ultra Deep Field (UDF) carried out using the slitless grism spectroscopy mode of the Advance Camera for Surveys on the Hubble Space Telescope (HST). The GRIsm ACS Program for Extragalactic Science (GRAPES) achieves continuum detection as faint as $z_{AB}=27.2$ using 40 orbits ($9.2 \times 10^4$ seconds) on HST. The data were taken at four orientation angles to correct for the overlap of spectra. GRAPES data provide a unique, uninterrupted, low resolution (R=100) spectral coverage for $5500\AA < \lambda < 10500\AA$, and allow us to detect high redshift galaxies at $4 < z < 7$ whether they have $\lya$ lines or just show the Lyman Break, as well as find low luminosity AGNs in an unbiased fashion. This paper describes in detail the observations and the data reduction, and examines the quality of the extracted spectra. Subsequent papers will deal with the analysis of the data. The extracted and calibrated GRAPES spectra will be available from MAST at STScI.

The X-ray cluster C0337-2522 at redshift $z=0.59$ hosts in its core a group of five elliptical galaxies. Using N-body simulations we show that a multiple merging event among the five galaxies is expected to take place in the next few Gyrs, forming a central brightest cluster galaxy. We also find indications that dynamical friction heating associated with this event is likely to modify the central slope of the cluster dark matter density profile.

Motivated by recent evidence that the epoch of reionization of hydrogen may have ended at a redshift as low as z~6 we consider the detectability of the sources responsible for this reionization. The main idea is that reionization places limits on the mean surface brightness of the population of reionization sources. Reducing the number of model-dependent assumptions to a minimum, we discuss the observability of these sources with existing and planned telescopes. We define a family of models characterized by two parameters: the Lyman continuum escape fraction f_c from the sources, and the clumpiness parameter C of the intergalactic medium. The minimum surface brightness model corresponds to a value of unity for both parameters. We find that the detection of the non-ionizing UV continuum of the reionization sources will be difficult to accomplish before the launch of JWST if these sources have a mean surface brightness close to the minimum value. However, if the values of f_c and C are more realistic, the reionization sources may well be detected by HST before JWST, perhaps by exploiting gravitational lensing amplification by foreground clusters of galaxies. Instead of a detection in the continuum, one may attempt to detect the Lyalpha emission line by narrow-band imaging. Present Lyalpha searches at z>6 suggest that either the typical sources are fainter than dwarf galaxies, or that the escape fraction of ionizing photons is much higher than 50%, so that there are no bright compact HII regions formed around the ionizing sources.

Limiting the number of model-dependent assumptions to a minimum, we discuss the detectability of the sources responsible for reionization with existing and planned telescopes. We conclude that if reionization sources are UV-efficient, minimum luminosity sources, then it may be difficult to detect them before the advent of the James Webb Space Telescope (JWST). The best approach before the launch of JWST may be either to exploit gravitational lensing by clusters of galaxies, or to search for strong Ly-alpha sources by means of narrow-band excess techniques or slitless grism spectroscopy.

According to the galactic cannibalism model, cD galaxies are formed in the center of galaxy clusters by merging of massive galaxies and accretion of smaller stellar systems: however, observational examples of the initial phases of this process are lacking. We have identified a strong candidate for this early stage of cD galaxy formation: a group of five elliptical galaxies in the core of the X-ray cluster C0337-2522 at redshift z=0.59. With the aid of numerical simulations, in which the galaxies are represented by N-body systems, we study their dynamical evolution up to z=0; the cluster dark matter distribution is also described as a N-body system. We find that a multiple merging event in the considered group of galaxies will take place before z=0 and that the merger remnant preserves the Fundamental Plane and the Faber-Jackson relations, while its behavior with respect to the Mbh-sigma relation is quite sensitive to the details of black hole merging [abridged].

As part of the Hubble Deep Field South program, a set of shorter 2-orbit observations were obtained of the area adjacent to the deep fields. The WFPC2 flanking fields cover a contiguous solid angle of 48 square arcminutes. Parallel observations with the STIS and NICMOS instruments produce a patchwork of additional fields with optical and near-infrared (1.6 micron) response. Deeper parallel exposures with WFPC2 and NICMOS were obtained when STIS observed the NICMOS deep field. These deeper fields are offset from the rest, and an extended low surface brightness object is visible in the deeper WFPC2 flanking field. In this data paper, which serves as an archival record of the project, we discuss the observations and data reduction, and present SExtractor source catalogs and number counts derived from the data. Number counts are broadly consistent with previous surveys from both ground and space. Among other things, these flanking field observations are useful for defining slit masks for spectroscopic follow-up over a wider area around the deep fields, for studying large-scale structure that extends beyond the deep fields, for future supernova searches, and for number counts and morphological studies, but their ultimate utility will be defined by the astronomical community.

We present HST/STIS spectra of the Sbc spiral galaxy NGC 4041 which were used to map the velocity field of the gas in its nuclear region. We detect the presence of a compact (r~0.4" ~40 pc), high surface brightness, rotating nuclear disk co-spatial with a nuclear star cluster. The disk is characterized by a rotation curve with a peak to peak amplitude of ~40 km/s and is systematically blueshifted by 10 - 20 km/s with respect to the galaxy systemic velocity. With the standard assumption of constant mass-to-light ratio and with the nuclear disk inclination taken from the outer disk, we find that a dark point mass of 1(-0.7;+0.6) 10^7 Msun is needed to reproduce the observed rotation curve. However the observed blueshift suggests the possibility that the nuclear disk could be dynamically decoupled. Following this line of reasoning we relax the standard assumptions and find that the kinematical data can be accounted for by the stellar mass provided that either the central mass-to-light ratio is increased by a factor of ~2 or that the inclination is allowed to vary. This model results in a 3 sigma upper limit of 6 10^6 Msun on the mass of any nuclear black hole. Overall, our analysis only allows us to set an upper limit of 2 10^7 Msun on the mass of the nuclear black hole. If this upper limit is taken in conjunction with an estimated bulge B magnitude of -17.7 and with a central stellar velocity dispersion of 95 km/s, then these results are not inconsistent with both the MBH-Lsph and the MBH-sigma correlations. Constraints on black hole masses in spiral galaxies of types as late as Sbc are still very scarce and therefore the present result adds an important new datapoint to our understanding of black hole demography.

We present the first results of a study of the expected properties of the first stellar generations in the Universe. In particular, we consider and discuss a series of properties that, on the basis of the emission from associated HII regions, permit one to discern {\it bona fide} primeval stellar generations from the ones formed after pollution from supernova explosions. The expected performance of NGST for the study and the characterization of primordial sources is also discussed.

This paper reviews our current understanding of the process of re-ionization of the Universe, focusing especially on those models where re-ionization is caused by UV radiation from massive stars. After reviewing the expected properties of stars at zero metallicity, I discuss the properties of primordial HII regions and their observability.

We discuss dynamical mass measurements of high-z galaxies with the Next Generation Space Telescope (NGST). In particular, we review some of the observational limits with the current instrument/telescope generation, we discuss the redshift limits and caveats for absorption and emission lines studies with NGST, and the existence of suitable targets at high redshift. We also briefly summarize strengths and weaknesses of proposed NGST instruments for dynamical studies.

The evolution of substructure embedded in non-dissipative dark halos is studied through N-body simulations of isolated systems, both in and out of initial equilibrium, complementing cosmological simulations of the growth of structure. We determine by both analytic calculations and direct analysis of the N-body simulations the relative importance of various dynamical processes acting on the clumps, such as the removal of material by global tides, clump-clump heating, clump-clump merging and dynamical friction. Our comparison between merging and disruption processes implies that spiral galaxies cannot be formed in a proto-system that contains a few large clumps, but can be formed through the accretion of many small clumps; elliptical galaxies form in a more clumpy environment than do spiral galaxies. Our results support the idea that the central cusp in the density profiles of dark halos is the consequence of self-limiting merging of small, dense halos. This implies that the collapse of a system of clumps/substructure is not sufficient to form a cD galaxy, with an extended envelope; plausibly subsequent accretion of large galaxies is required. Persistent streams of material from disrupted clumps can be found in the outer regions of the final system, and at an overdensity of around 0.75, can cover 10% to 30% of the sky.

We report evidence for nuclear stellar disks in 3 early-type spirals, namely NGC 1425, NGC 3898 and NGC 4698, revealed by WFPC2/F606W images out of a sample of 38 spiral galaxies, selected from the Hubble Space Telescope Data Archive. Adopting the photometric method introduced by Scorza & Bender we derived their central surface brightness and scalelength by assuming them to be infinitesimally thin exponential disks. No nuclear disk was found in barred galaxies or galaxies of Hubble type later than Sb. The external origin of the disk in NGC 4698 is strongly suggested by its orthogonal geometrical decoupling with respect to the host galaxy.

We have measured the Fundamental Plane (FP) parameters for a sample of 30 field early-type galaxies (E/S0) in the redshift range 0.1<z<0.66. We find that: i) the FP is defined and tight out to the highest redshift bin; ii) the intercept \gamma evolves as d\gamma/dz=0.58+0.09-0.13 (for \Omega=0.3, \Omega_{\Lambda}=0.7), or, in terms of average effective mass to light ratio, as d\log(M/L_B)/dz=-0.72+0.11-0.16, i.e. faster than is observed for cluster E/S0 -0.49+-0.05. In addition, we detect [OII] emission >5\AA in 22% of an enlarged sample of 42 massive E/S0 in the range 0.1<z<0.73, in contrast with the quiescent population observed in clusters at similar z. We interpret these findings as evidence that a significant fraction of massive field E/S0 experiences secondary episodes of star-formation at z<1.

We present surface brightness profiles for 56 of the 78 spiral galaxies observed in the HST/{\tt NICMOS/Camera-2} F160W snapshot survey introduced in Carollo et al. (2001; paper I), as well as surface brightness profiles for 23 objects out of the 41 that were also observed in the F110W filter. We fit these surface brightness profiles with the Nuker law of Lauer et al. (1995), and use the smooth analytical descriptions of the data to compute the average nuclear stellar cusp slopes $<\gamma>$ in the 0.1$''-0.5''$ radial range. Our main result is the startling similarity between the nuclear stellar cusp slopes $<\gamma>$ in the near-infared compared to those derived in the visual passband. This similarity has several implications: (1.) There are typically little or no optical-NIR {\it global} color gradients, and thus no global stellar population variations, inside $\sim 50$-100 pc from the nucleus in nearby spirals. (2.) The large observed range of the strength of the nuclear stellar cusps seen in the HST optical study of spiral galaxies reflects a physical difference between galaxies, and is not an artifact caused by nuclear dust and/or recent star formation. (3.) The dichotomy between $R^{1/4}$ bulges, with steep nuclear stellar cusps $<\gamma > \sim 1$, and exponential bulges, with shallow nuclear stellar cusps $<\gamma > < 0.3$, is also not an artifact of the effects of dust or recent star formation. (4) A surrounding massive disk appears to have no effect on the rise of the stellar density distribution within the innermost hundred pc of the

This is the first of two papers where we present the analysis of an HST NICMOS-Cam2 near-infared (NIR) snapshot survey in the F160W (H) filter for a sample of 78 spiral galaxies selected from the UGC and ESOLV catalogs. For 69 of these objects we provide nuclear color information derived by combining the H data either with additional NICMOS F110W (J) images or with V WFPC2/HST data. Here we present the NIR images and the optical-NIR color maps. We focus our attention on the properties of the photometrically-distinct `nuclei' which are found embedded in most of the galaxies, and provide measurements of their half-light radii and magnitudes in the H (and when available, in the J) band. We find that: (i) In the NIR, the nuclei embedded in the bright early- to intermediate-type galaxies span a much larger range in brightness than the nuclei which are typically found embedded in bulgeless late-type disks; (ii) Nuclei are found in both non-barred and barred hosts, in large-scale ($\gta$ 1 kpc) as well as in nuclear (up to a few 100pc) bars; (iii) There is a significant increase in half-light radius with increasing luminosity of the nucleus in the early/intermediate types; (iv) The nuclei of early/intermediate-type spirals cover a large range of optical-NIR colors, from $V-H\approx$ -0.5 to 3. Some nuclei are bluer and others redder than the surrounding galaxy, indicating the presence of activity or reddening by dust in many of these systems; (v) Some early/intermediate nuclei are elongated and/or slightly offset from the isophotal center of the host galaxy. On average, however, these nuclei appear as centered, star-cluster-like structures similar to those which are found in the late-type disks. [abridged]

We have carried out a search for Lyman alpha emission from galaxies at z~=2.4 over a field of 1200 sq. arcmin using the CFH12K camera at the CFHT and a custom medium band filter. The search has uncovered 58 candidates, corresponding to a completeness-corrected source density of ~0.07 sq. arcmin^-2. Our sources have red colors (B-I~=1.8) which imply either that a large fraction of the light is highly reddened and we are detecting Lyman alpha through special lines of sight, or that these objects contain an underlying older stellar population. While for each individual object we cannot discriminate between these alternatives, we conclude that most of the objects actually contain an older component because the star formation rates inferred from the picture based on reddening, applied to all candidates, would imply an exceedingly high star formation rate, i.e. more than two orders of magnitude above the peak cosmic star formation rate (e.g. Lilly et al. 1996).

We report on the study of a sample of 25 field early-type galaxies, in the redshift range z~0.1-0.5, selected on the basis of colours and morphology from the HST-MDS. Field early-type galaxies define a tight Fundamental Plane (FP) out to z~0.4, with scatter unchanged with respect to local samples, within the observational errors. The intermediate redshift FP is offset with respect to the FP of the Coma Cluster. The offset of the FP is found to increase with redshift. The evolution of the FP is studied quantitatively with a Bayesian-Montecarlo technique. By applying this technique, we find that the offset of the intercept of the FP (\Delta \gamma) with respect to the local FP increases as \Delta \gamma = \tau z with the following 68 per cent limits: 0.33<\tau<0.44 (for \Omega=1, \Omega_{\Lambda}=0) or 0.44<\tau<0.56 (for \Omega=0.3,\Omega_{\Lambda}=0.7). In addition, we interpret the results in terms of the evolution of the stellar populations, under the assumption of passive evolution. In a single-burst scenario, the observed properties are consistent with those of a stellar population formed at z>2 (for \Omega=1, \Omega_{\Lambda}=0, H_0=50 \kms Mpc^{-1}) or 0.8<z<1.6 (for \Omega=0.3, \Omega_{\Lambda}=0.7, H_0=65 \kms Mpc^{-1}). If a small fraction of the stellar mass is formed in a secondary burst, the primary burst may have occurred at higher z. Finally, the intercept and scatter of the FP found for field early-type galaxies and for cluster data at z~0.4 are mutually consistent, within the observational errors. If higher redshift (up to z=0.83) cluster data are considered, the ages of the stellar populations of field early-type galaxies inferred from a single-burst scenario are found to be marginally smaller than the ages derived for the cluster galaxies. [shortened]

A sample of field early-type galaxies (E/S0) at intermediate redshift ($z\sim0.1-0.6$) is selected, based on morphology and colours from HST-WFPC2 parallel images. Photometric structural parameters (effective radius $R_{\tx{e}}$ and effective surface brightness $SB_{\tx{e}}$) are derived through the F606W and F814W filters, using luminosity profile fitting and two-dimensional fitting techniques. The combined parameter that enters the Fundamental Plane ($\log R_{\tx{e}}-\beta SB_{\tx{e}}$, with $\beta\approx0.32$) is shown to suffer from significantly smaller uncertainties (r.m.s. 0.03) than the individual structural parameters (e.g. $\sim 15$ per cent r.m.s. on the effective radius). High signal-to-noise intermediate resolution spectra, taken at the ESO-3.6m, yield redshifts for 35 galaxies and central velocity dispersions for 22 galaxies. Central velocity dispersions are derived using a library of stellar templates covering a wide range of spectral types, in order to study the effects of templates mismatches. The average random error on the central velocity dispersion is found to be 8 per cent and the average systematic error due to template mismatch is found to be 5 per cent. The errors on velocity dispersion measurement and the effects of template mismatches are studied by means of extensive Montecarlo simulations. In addition, we investigate whether the determination of the velocity dispersion is sensitive to the spectral range used, finding that the value of velocity dispersion is unchanged when the spectral regions that include the absorption features Ca HK and NaD are masked out during the fit.

The dynamical friction timescale for globular clusters to sink to the center of a dwarf elliptical galaxy (dE) is significantly less than a Hubble time if the halos have King-model or isothermal profiles and the globular clusters formed with the same radial density profile as the underlying stellar population. We examine the summed radial distribution of the entire globular cluster systems and the bright globular cluster candidates in 51 Virgo and Fornax Cluster dEs for evidence of dynamical friction processes. We find that the summed distribution of the entire globular cluster population closely follows the exponential profile of the underlying stellar population. However, there is a deficit of bright clusters within the central regions of dEs (excluding the nuclei), perhaps due to the orbital decay of these massive clusters into the dE cores. We also predict the magnitude of each dE's nucleus assuming the nuclei form via dynamical friction. The observed trend of decreasing nuclear luminosity with decreasing dE luminosity is much stronger than predicted if the nuclei formed via simple dynamical friction processes. We find that the bright dE nuclei could have been formed from the merger of orbitally decayed massive clusters, but the faint nuclei are several magnitudes fainter than expected. These faint nuclei are found primarily in M_V > -14 dEs which have high globular cluster specific frequencies and extended globular cluster systems. In these galaxies, supernovae-driven winds, high central dark matter densities, extended dark matter halos, the formation of new star clusters, or tidal interactions may act to prevent dynamical friction from collapsing the entire globular cluster population into a single bright nucleus.

We describe an on-going HST program aimed at determining the relationship between the nuclear black hole mass and bulge mass in spiral galaxies. We have selected a volume limited sample of 54 nearby spiral galaxies for which we already have ground based emission line rotation curves, CCD surface photometry and radio maps. We are now obtaining HST/STIS longslit observations of each of the galaxies in the sample in order to determine the nuclear Halpha rotation curve at high (~0.1") spatial resolution. We will use these data to measure the unresolved dark mass concentration at the nucleus of each object. Here we show the first results from observations of objects in the sample.

We derive the light profiles for a sample of 25 dwarf elliptical galaxies observed by us with HST/WFPC2 in F555W and F814W. These profiles are fitted with Nuker, R^1/4, exponential, and Sersic laws, and are also used to derive the nuclear cusp slopes gamma. We discuss the correlation of nuclear cusp slope with galactic luminosity, presence of a nucleus, and type of light profile. The results are compared to those found in the literature for elliptical galaxies and the bulges of spiral galaxies. We find that as a class the nuclear regions of dwarf ellipticals are very similar to those of the exponential bulges of spiral galaxies, and have nuclear cusp slopes shallower than those of bulges that were well fitted by a de Vaucouleurs R^1/4 profile with the same luminosity. For the 14 nucleated galaxies in our sample this conclusion is less certain than for the 11 non-nucleated objects since it relies on an extrapolation of galaxy light under the nucleus. In terms of their light profiles and nuclear properties, most spheroidal stellar systems can be broadly divided into two subclasses: the exponential shallow cusp objects and the R^1/4, steep cusp objects. Membership of a class does not appear to correlate with the presence of a massive stellar disk.

The Hubble Deep Field-South observations targeted a high-galactic-latitude field near QSO J2233-606. We present WFPC2 observations of the field in four wide bandpasses centered at roughly 300, 450, 606, and 814 nm. Observations, data reduction procedures, and noise properties of the final images are discussed in detail. A catalog of sources is presented, and the number counts and color distributions of the galaxies are compared to a new catalog of the HDF-N that has been constructed in an identical manner. The two fields are qualitatively similar, with the galaxy number counts for the two fields agreeing to within 20%. The HDF-S has more candidate Lyman-break galaxies at z > 2 than the HDF-N. The star-formation rate per unit volume computed from the HDF-S, based on the UV luminosity of high-redshift candidates, is a factor of 1.9 higher than from the HDF-N at z ~ 2.7, and a factor of 1.3 higher at z ~ 4.

We briefly discuss our survey of Extremely Red Objects (EROs) with HST/NICMOS and deep ground-based images and summarize its results. We use the high-resolution NICMOS images to sub-divide EROs into morphological classes, tentatively arguing that about 60 % of EROs (18 out of 30) are early-type galaxies. About 20 % (6 out of 30 objects) appear to be flattened disk systems. Another 10 % (3 out of 30) have an irregular morphology. Thus, combining disk and irregular systems, our data are compatible with 30 % of EROs being dusty starbursts. A small fraction (3 out of 30) appear point-like even on the HST images and could therefore be either compact high redshift galaxies or stars in the Milky Way halo.

We have analysed $V$, $H$ and $J$ HST images for a sample of early- to late-type spiral galaxies, and reported elsewhere the statistical frequency of $R^{1/4}$-law and exponential bulges in our sample as a function of Hubble type, and the frequency of occurrence and structural properties of the resolved central nuclei hosted by intermediate- to late-type bulges and disks (see references in the text). Here we use these data to show that: {\it (i)} The $V-H$ color distribution of the $R^{1/4}$ bulges peaks around $<V-H> \sim 1.3$, with a sigma $\Delta (V-H) \sim 0.1$ magnitudes. Assuming a solar metallicity, these values correspond to stellar ages of $\approx 6 \pm 3$ Gyrs. In contrast, the $V-H$ color distribution of the exponential bulges peaks at $<V-H> \sim 0.9$ and has a sigma $\Delta (V-H) \sim 0.4$ mags. This likely implies significantly smaller ages and/or lower metallicities for (a significant fraction of the stars in) the exponential bulges compared to the $R^{1/4}$-law spheroids. {\it (ii)} Most of the central nuclei hosted by the exponential bulges have $V-H$ and $J-H$ colors which are compatible with relatively unobscured stellar populations. Assuming no or little dust effects, ages $\gta 1$ Gyrs are suggested for these nuclei, which in turn imply masses of about a few 10$^6$ to a few 10$^7$ M$_\odot$, i.e., sufficient to dissolve progenitor bars with masses consistent with those inferred for the exponential bulges by their luminosities. {\it (iii)} While different bulge-nucleus pairs cover a large range of $V-H$ colors, each bulge-nucleus pair has quite similar $V-H$ colors, and thus possibly similar stellar populations. (abridged)

We present preliminary results on the evolution of the stellar populations of field early-type galaxies (E/S0) from z=0.4 to z=0. The diagnostic tool used in this study is the Fundamental Plane (FP), a tight empirical correlation between their central velocity dispersion (sigma), effective radius (Re), and effective surface brightness (SBe), which is observed to hold in the local Universe. Using HST-WFPC2 archive images and spectra obtained at the ESO-3.6m telescope we measured the FP parameters for a sample of ~30 field E/S0s at z=0.2-0.4. Remarkably, field E/S0s at intermediate redshift also define a tight FP, with scatter unchanged with respect to that of local samples. The intermediate redshift FP is offset from the local one, in the sense that, for given RE and sigma, galaxies are brighter at z=0.4 than at z=0. The implication of the offset of the FP in terms of passive evolution of the stellar population depends on its star formation history. In a single burst scenario, the stellar populations of field E/S0s were formed at z=0.8-1.6 (Omega=0.3; Omega_{Lambda}=0.7; H_0=50 km/s/Mpc). Alternatively, the bulk of stars (90% in mass) can be formed at high redshift (z~3), and the rest in a secondary burst occurred more recently (z~0.5-0.8).

We have collected data for 69 Damped Lyman-alpha (DLA) systems, to investigate the chemical evolution of galaxies in the redshift interval 0.0 < z < 4.4. In doing that, we have adopted the most general approach used so far to correct for dust depletion. The best solution, obtained through chi^2 minimization, gives as output parameters the global DLA metallicity and the dust-to-metals ratio. Clear evolution of the metallicity vs. redshift is found (99.99% significance level), with average values going from ~1/30 solar at z~4.1 to ~3/5 solar at z~0.5. We also find that the majority of DLAs (~60%) shows dust depletion patterns which most closely resemble that of the warm halo clouds in the Milky Way, and have dust-to-metals ratios very close to warm halo clouds.

We have collected optical follow-up observations from the ground and the Hubble Space Telescope (HST) for 13.74 square ' of archival images taken with the Near Infrared Camera and Multi Object Spectrograph (NICMOS). We use two criteria to select E/S0 galaxy candidates at z>~1 and z>~1.5 (HizECs) based on colors and infrared morphology. The observed surface density of HizECs is significantly smaller than what is predicted assuming the local luminosity function (LF) of E/S0, constant comoving density, and pure passive evolution with a high redshift of formation of the stellar populations (z>~ 5). On the other hand, assuming a low redshift (z<~2) of formation, not enough HizECs are predicted. The data are very well matched assuming that a substantial fraction of present-day E/S0s (10-66%, depending on cosmology and on the local LF of E/S0s) formed at z = 3 or higher and then evolved passively. The rest of the current E/S0 population may have formed at lower redshift, or may have been rejected by the selection criteria because of interactions or more recent episodes of star-formation (disturbed morphology and/or bluer colors).

We present measurements of the Fundamental Plane (FP) parameters (the effective radius, the mean effective surface brightness, and the central velocity dispersion) of six field elliptical galaxies at intermediate redshift. The imaging is taken from the Medium Deep Survey of the Hubble Space Telescope, while the kinematical data are obtained from long-slit spectroscopy using the 3.6-m ESO telescope. The Fundamental Plane appears well defined in the field even at redshift $\approx$ 0.3. The data show a shift in the FP zero point with respect to the local relation, possibly indicating modest evolution, consistent with the result found for intermediate redshift cluster samples. The FP slopes derived for our field data, plus other cluster ellipticals at intermediate redshift taken from the literature, differ from the local ones, but are still consistent with the interpretation of the FP as a result of homology, of the virial theorem and of the existence of a relation between luminosity and mass, $L \propto M^{\eta}$. We also derive the surface brightness vs. effective radius relation for nine galaxies with redshift up to $z \approx0.6$, and data from the literature; the evolution that can be inferred is consistent with what is found using the FP.

The quasar in the Hubble Deep Field South (HDFS), J2233-606 (z=2.23) has been exhaustively observed by ground based telescopes and by the STIS spectrograph on board the Hubble Space Telescope at low, medium and high resolution in the spectral interval from 1120 A to 10000 A. This very large base-line represents a unique opportunity to study in detail the distribution of clouds associated with emitting structures in the field of the quasar and in nearby fields already observed as part of the HDFS campaign. Here we report the main properties of the Lyman-alpha clouds in the intermediate redshift range 1.20-2.20, where our present knowledge has been complicated by the difficulty in producing good data. The number density is shown to be higher than what is expected by extrapolating the results from both lower and higher redshifts: 63\pm8 lines with log N_{HI}\geq14.0 are found (including metal systems) at <z>=1.7, to be compared with ~40 lines predicted by extrapolating from previous studies. The redshift distribution of the Lyman-alpha clouds shows a region spanning z=1.383-1.460 (comoving size of 94 h^{-1}_{65} Mpc, Omega_o=1) with a low density of absorption lines; we detect 5 lines in this region, compared with the 16 expected from an average density along the line of sight. The two point correlation function shows a positive signal up to scales of about 3 h^{-1}_{65} Mpc and an amplitude that is larger for larger HI column densities. The average Doppler parameter is about 27 km/s, comparable to the mean value found at z > 3, thus casting doubts on the temperature evolution of the Lyman-alpha clouds.

The combined use of the ESO Very Large Telescope (VLT) UT1 Science Verification (SV) images and of the Hubble Space Telescope (HST) Hubble Deep Field South observations allows us to strengthen the identification as a candidate elliptical galaxy of the Extremely Red Object HDFS 223251-603910 previously identified by us on the basis of NICMOS and Cerro Tololo Interamerican Observatory imaging. The photometry presented here includes VLT data in U, B, V, R, I, a STIS unfiltered image, NICMOS J, H, and K band data, thus combining the 16.5 hours of VLT SV exposures with 101 hours of HST observing. The object is detected in all images except the VLT U band and is one of the reddest known with B-K=9.7+-0.5. We consider a wide range of models with different ages, metallicities, star formation histories and dust content, and conclude that the observed spectral energy distribution agrees best with that of an old elliptical galaxy at redshift just below 2. Alternative possibilities are discussed in light of their likelihood and of the perspective of spectroscopic confirmation.

The specific globular cluster frequencies (S_N) for 24 dwarf elliptical (dE) galaxies in the Virgo and Fornax Clusters and the Leo Group imaged with the Hubble Space Telescope are presented. Combining all available data, we find that for nucleated dEs --- which are spatially distributed like giant ellipticals in galaxy clusters --- S_N(dE,N)=6.5 +- 1.2 and S_N increases with M_V, while for non-nucleated dEs --- which are distributed like late-type galaxies --- S_N(dE,noN)=3.1 +- 0.5 and there is little or no trend with M_V. The S_N values for dE galaxies are thus on average significantly higher than those for late-type galaxies, which have S_N < 1. This suggests that dE galaxies are more akin to giant Es than to late-type galaxies. If there are dormant or stripped irregulars hiding among the dE population, they are likely to be among the non-nucleated dEs. Furthermore, the similarities in the properties of the globular clusters and in the spatial distributions of dE,Ns and giant Es suggest that neither galaxy mass or galaxy metallicity is responsible for high values of S_N. Instead, most metal-poor GCs may have formed in dwarf-sized fragments that merged into larger galaxies.

We discuss how laboratory experiments can be used to place constraints on possible variations of the fine structure constant alpha in the observationally relevant redshift interval z ~= 0 - 5, within a rather general theoretical framework. We find a worst case upper limit for Delta alpha / alpha of 8 x 10^-6 for z <= 5 and Delta alpha / alpha of 0.9 x 10^-6 for z <= 1.6. The derived limits are at variance with the recent findings by Webb et al. (1998), who claim an observed variation of Delta alpha/alpha = -2.6 +- 0.4 x 10^-5 at 1<z<1.6.

We report the serendipitous discovery of an extremely red object in the Hubble Deep Field South (HDFS) Test NICMOS (Near Infrared Camera and Multi Object Spectrograph) field of view. The object is resolved in the NICMOS image and has a light profile very well described by an r^{1/4} law with effective radius r_e=0."20+-0."05 and H_{AB}=21.7+-0.1 magnitudes. In contrast, the galaxy is undetected in the R and I band ground based images taken at the CTIO 4 m Blanco Telescope, giving a lower limit to the color of (R-H)_{AB}>3.9 and (I-H)_{AB}>3.5 at the 95% confidence level. The colors of a range of synthetic galactic spectra are computed, showing that the object is likely to be an ``old'' elliptical galaxy at redshift z>~1.7. Alternatively the colors can be reproduced by an ``old'' elliptical galaxy at somewhat lower redshift (z>~1) with significant amount of dust, or by a younger galaxy at higher redshift. This object represents a very interesting target for future VLT observations.

Globular clusters (GCs) provide one way of determining when the most important periods of star formation in dwarf elliptical (dE) galaxies occurred. Thus, they are also useful for comparing the star formation histories of dwarf and giant galaxies. We present results on the GC systems of 24 nearby dE galaxies imaged with the Hubble Space Telescope. We find that the specific GC frequency, S_N, increases with M_V. Also, S_N in dEs is more like the values found in giant ellipticals rather than in spirals, and the mean S_N for nucleated dEs is roughly a factor of 2 higher than for non-nucleated dEs. The colors suggest that most of the GCs formed before the bulk of the stars in the galaxies.

In this paper, the third of a series dedicated to the investigation of the nuclear properties of spiral galaxies, we have (i) modelled the WFPC2 F606W nuclear surface brightness profiles of 41 spiral galaxies presented in Carollo et al. 1997c, 1998 with the analytical law introduced by Lauer et al. 1995, and (ii) deprojected these surface brightness profiles and their analytical fits, so as to estimate the nuclear stellar densities of bulges of spiral galaxies. We find that the nuclear stellar cusps (quantified by the average logarithmic slope of the surface brightness profiles within 0.1''-0.5'') are significantly different for R^1/4-law and exponential bulges. The former have nuclear properties similar to those of early-type galaxies, i.e. similar values of nuclear cusps for comparable luminosities, and increasingly steeper stellar cusps with decreasing luminosity. By contrast, exponential bulges have (underlying the light contribution from photometrically distinct, central compact sources) comparative shallower stellar cusps, and likely lower nuclear densities, than R^1/4-law bulges.

We report the analysis of HST WFPC2 F606W images of 40 spirals belonging to the sample introduced in Carollo et al. (paper I), where 35 other targets were discussed. We describe the optical morphology of the new 40 galaxies, derive the surface brightness profiles for 25 of them, and present the results of photometric decompositions of these profiles into a ``bulge'' (R^1/4 or exponential) and a disk component. The analysis of the sample of 75 galaxies puts on a statistically more solid ground the main results presented in paper I. Furthermore, the analysis of the sample of 75 objects shows that: A. Several of the non-classical inner structures are well fitted by an exponential profile. These ``exponential bulges'' are fainter than R^1/4 bulges for given total galaxy luminosity and (catalog) Hubble type later than Sab. B. Irregular/exponential bulges typically host central compact sources. C. The central sources are present in all types of disk galaxies, starting with systems as early as S0a. About 60 % of Sb to Sc galaxies host a central compact source. Many of the galaxies which host compact sources contain a barred structure. D. Galaxies with apparent nuclear star formation and the brightest compact sources, are preferentially the early- and intermediate- type (S0a-Sb) systems. E. None of the features depends on environment: isolated and not isolated galaxies show indistinguishable properties. Independently from the physical nature of the non-classical inner structures, our main conclusion is that a significant fraction of galaxies classified from the ground as relatively early-type shows a rich variety of central properties, and little or no morphological/photometric evidence for a smooth, R^1/4-law bulge.

I present galaxy models resulting from violent relaxation in the presence of a pre-existing black hole. The models are computed by maximizing the entropy of the stellar dynamical system. I show that their properties are very similar to those of adiabatic growth models for a suitable choice of parameters. This suggests that observations of nuclear light profiles and kinematics alone may not be sufficient to discriminate between scenarios where a black hole grows adiabatically in the core of a galaxy, and scenarios where the black hole formation preceeds galaxy formation.

We present WFPC2 images in the F450W, F606W and F814W filters of the interacting pair of galaxies NGC 454. Our data indicate that the system is in the early stages of interaction. A population of young star-clusters has formed around the late component, and substantial amounts of gas have sunk into the center of the earlier component, where it has not yet produced significant visible star formation or nuclear activity. We have photometric evidence that the star-clusters have strong line emission, which indicate the presence of a substantial component of hot, massive stars which formed less than 5-10 Myrs ago.

We address the problem of deriving, from the observed projected metallicity gradients, the intrinsic metallicity distribution of elliptical galaxies as a function of their integrals of motion. The method is illustrated by an application to anisotropic spherical Hernquist models. We also compare the derived metallicity distribution with those expected from two very simple models of galaxy formation and find that the more dissipative scheme agrees better with the typical metallicity distribution of ellipticals.