We present the discovery and high-cadence follow-up observations of SN 2018ivc, an unusual Type II supernova that exploded in NGC 1068 (D=10.1 Mpc). The light curve of SN 2018ivc declines piecewise-linearly, changing slope frequently, with four clear slope changes in the first 30 days of evolution. This rapidly changing light curve indicates that interaction between the circumstellar material and ejecta plays a significant role in the evolution. Circumstellar interaction is further supported by a strong X-ray detection. The spectra are rapidly evolving and dominated by hydrogen, helium, and calcium emission lines. We identify a rare high-velocity emission-line feature blueshifted at ~7800 km/s (in Ha, Hb, Pb, Pg, HeI, CaII), which is visible from day 18 until at least day 78 and could be evidence of an asymmetric progenitor or explosion. From the overall similarity between SN 2018ivc and SN 1996al, the \Ha{} equivalent width of its parent HII region, and constraints from pre-explosion archival Hubble Space Telescope images, we find that the progenitor of SN 2018ivc could be as massive as 52 Msun but is more likely <12 Msun. SN 2018ivc demonstrates the importance of the early discovery and rapid follow-up observations of nearby supernovae to study the physics and progenitors of these cosmic explosions.

We make use of the ASAGAO, deep 1.2 mm continuum observations of a 26 arcmin$^2$ region in the GOODS-South field obtained with ALMA, to probe dust-enshrouded star formation in $K$-band selected (i.e., stellar mass selected) galaxies, which are drawn from the ZFOURGE catalog. Based on the ASAGAO combined map, which was created by combining ASAGAO and ALMA archival data in the GOODS-South field, we find that 24 ZFOURGE sources have 1.2 mm counterparts with a signal-to-noise ratio $>$ 4.5 (1$\sigma\simeq$ 30 - 70 $\mu$Jy beam$^{-1}$ at 1.2 mm). Their median redshift is estimated to be $z_\mathrm{median}=$ 2.38 $\pm$ 0.14. They generally follow the tight relationship of the stellar mass versus star formation rate (i.e., the main sequence of star-forming galaxies). ALMA-detected ZFOURGE sources exhibit systematically larger infrared (IR) excess (IRX $\equiv L_\mathrm{IR}/L_\mathrm{UV}$) compared to ZFOURGE galaxies without ALMA detections even though they have similar redshifts, stellar masses, and star formation rates. This implies the consensus stellar-mass versus IRX relation, which is known to be tight among rest-frame-UV-selected galaxies, can not fully predict the ALMA detectability of stellar-mass-selected galaxies. We find that ALMA-detected ZFOURGE sources are the main contributors to the cosmic IR star formation rate density at $z$ = 2 - 3.

We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to estimate infrared luminosities and dust temperatures simultaneously. In conjunction with previous ALMA observations for six galaxies at z>6, we confirm that all the nine z=6-9 galaxies have high [OIII]/[CII] ratios of L[OIII]/L[CII]~3-20, ~10 times higher than z~0 galaxies. We also find a positive correlation between the [OIII]/[CII] ratio and the Lya equivalent width (EW) at the ~90% confidence level. We carefully investigate physical origins of the high [OIII]/[CII] ratios at z=6-9 using Cloudy, and find that high density of the interstellar medium, low C/O abundance ratio, and the cosmic microwave background attenuation are responsible to only a part of the z=6-9 galaxies. Instead, the observed high [OIII]/[CII] ratios are explained by 10-100 times higher ionization parameters or low photodissociation region (PDR) covering fractions of 0-10%, both of which are consistent with our [NII] observations. The latter scenario can be reproduced with a density bounded nebula with PDR deficit, which would enhance the Lya, Lyman continuum, and C+ ionizing photons escape from galaxies, consistent with the [OIII]/[CII]-Lya EW correlation we find.

We investigate the role of dense Mpc-scale environments in processing molecular gas of cluster galaxies as they fall into the cluster cores. We consider $\sim20$ luminous infrared galaxies (LIRGs) in intermediate-$z$ clusters, from the Hershel Lensing Survey and the Local Cluster Substructure Survey. They include MACS J0717.5+3745 at $z=0.546$ and Abell 697, 963, 1763, and 2219 at $z=0.2-0.3$. We have performed far infrared to ultraviolet spectral energy distribution modeling of the LIRGs, which span cluster-centric distances within $r/r_{200}\simeq0.2-1.6$. We have observed the LIRGs in CO(1$\rightarrow$0) or CO(2$\rightarrow$1) with the Plateau de Bure interferometer and its successor NOEMA, as part of five observational programs carried out between 2012 and 2017. We have compared the molecular gas to stellar mass ratio $M(H_2)/M_\star$, star formation rate (SFR), and depletion time ($\tau_{\rm dep}$) of the LIRGs with those of a compilation of cluster and field star forming galaxies. The targeted LIRGs have SFR, $M(H_2)/M_\star$, and $\tau_{\rm dep}$ that are consistent with those of both main sequence (MS) field galaxies and star forming galaxies from the comparison sample. However we find that the depletion time, normalized to the MS value, increases with increasing $r/r_{200}$, with a significance of $2.8\sigma$, which is ultimately due to a deficit of cluster core LIRGs with $\tau_{\rm dep}\gtrsim\tau_{\rm dep,MS}$. We suggest that a rapid exhaustion of the molecular gas reservoirs occurs in the cluster LIRGs and is effective in suppressing their star formation. This mechanism may explain the exponential decrease of the fraction of cluster LIRGs with cosmic time. The compression of the gas in LIRGs, possibly induced by intra-cluster medium shocks, may be responsible for the short depletion timescales, observed in a large fraction of cluster core LIRGs.

The NIRCam instrument on the upcoming James Webb Space Telescope (JWST) will offer an unprecedented view of the most distant galaxies. In preparation for future deep NIRCam extragalactic surveys, it is crucial to understand the color selection of high-redshift galaxies using the Lyman dropout technique. To that end, we have used the JAdes extraGalactic Ultradeep Artificial Realizations (JAGUAR) mock catalog to simulate a series of extragalactic surveys with realistic noise estimates. This enables us to explore different color selections and their impact on the number density of recovered high-redshift galaxies and lower-redshift interlopers. We explore how survey depth, detection signal-to-noise ratio, color selection method, detection filter choice, and the presence of the Ly$\alpha$ emission line affects the resulting dropout selected samples. We find that redder selection colors reduce the number of recovered high-redshift galaxies, but the overall accuracy of the final sample is higher. In addition, we find that methods that utilize two or three color cuts have higher accuracy because of their ability to select against low-redshift quiescent and faint dusty interloper galaxies. We also explore the near-IR colors of brown dwarfs and demonstrate that, while they are predicted to have low on-sky densities, they are most likely to be recovered in F090W dropout selection, but there are color cuts which help to mitigate this contamination. Overall, our results provide NIRCam selection methods to aid in the creation of large, pure samples of ultra high-redshift galaxies from photometry alone.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

We study six luminous Lyman-alpha emitters (LAEs) with very blue rest-frame UV continua at $5.7\le z \le 6.6$. These LAEs have previous HST and Spitzer IRAC observations. Combining our newly acquired HST images, we find that their UV-continuum slopes $\beta$ are in a range of $-3.4\le \beta \le -2.6$. Unlike previous, tentative detections of $\beta \simeq -3$ in photometrically selected, low-luminosity galaxies, our LAEs are spectroscopically confirmed and luminous ($M_{\rm UV}<-20$ mag). We model their broadband spectral energy distributions (SEDs), and find that two $\beta\simeq-2.6\pm0.2$ galaxies can be well fitted with young and dust-free stellar populations. However, it becomes increasingly difficult to fit bluer galaxies. We explore further interpretations by including non-zero LyC escape fraction $f_{\rm esc}$, very low metallicities, and/or AGN contributions. Assuming $f_{\rm esc}\simeq0.2$, we achieve the bluest slopes $\beta\simeq-2.7$ when nebular emission is considered. This can nearly explain the SEDs of two galaxies with $\beta\simeq-2.8$ and --2.9 ($\sigma_{\beta}=0.15$). Larger $f_{\rm esc}$ values and very low metallicities are not favored by the strong nebular line emission (evidenced by the IRAC flux) or the observed (IRAC 1 - IRAC 2) color. Finally, we find that the $\beta\simeq-2.9$ galaxy can potentially be well explained by the combination of a very young population with a high $f_{\rm esc}$ ($\ge0.5$) and an old, dusty population. We are not able to produce two $\beta \simeq -3.4 \pm0.4$ galaxies. Future deep spectroscopic observations are needed to fully understand these galaxies.

We explore the connection between dust and star formation, in the context of environmental effects on galaxy evolution. In particular, we exploit the susceptibility of dust to external processes to assess the influence of dense environment on star-forming galaxies. We have selected cluster Abell 1758 from the Local Cluster Substructure Survey (LoCuSS). Its complex dynamical state is an ideal test-bench to track dust removal and destruction in galaxies due to merger and accretion shocks. We present a systematic panchromatic study (from 0.15 $\rm \mu$m with GALEX to 500 $\rm \mu$m with Herschel) of spectroscopically confirmed star-forming cluster galaxies at intermediate redshift. We observe that the main subclusters (A1758N and A1758S) belong to two separate large-scale structures, with no overlapping galaxy members. Star-forming cluster members are distributed preferentially outside cluster central regions, and are not grouped in substructures. Rather, these galaxies are being funneled towards the main subclusters along separate accretion filaments. Additionally, we present the first study of dust-to-stellar (DTS) mass ratio used as indicator for local environmental influence on galaxy evolution. Star-forming cluster members show lower mean values (32% at 2.4$\rm \sigma$) of DTS mass ratio and lower levels of infrared emission from birth clouds with respect to coeval star-forming field galaxies. This picture is consistent with the majority of star-forming cluster members infalling in isolation. Upon accretion, star-formation is observed to decrease and warm dust is destroyed due to heating from the intracluster medium radiation, ram-pressure stripping and merger shocks.

A double source plane (DSP) system is a precious probe for the density profile of distant galaxies and cosmological parameters. However, these measurements could be affected by the surrounding environment of the lens galaxy. Thus, it is important to evaluate the cluster-scale mass for detailed mass modeling. We observed the {\it Eye of Horus}, a DSP system discovered by the Subaru HSC--SSP, with XMM--Newton. We detected two X-ray extended emissions, originating from two clusters, one centered at the {\it Eye of Horus}, and the other located $\sim100$ arcsec northeast to the {\it Eye of Horus}. We determined the dynamical mass assuming hydrostatic equilibrium, and evaluated their contributions to the lens mass interior of the Einstein radius. The contribution of the former cluster is $1.1^{+1.2}_{-0.5}\times10^{12}~M_{\odot}$, which is $21-76\%$ of the total mass within the Einstein radius. The discrepancy is likely due to the complex gravitational structure along the line of sight. On the other hand, the contribution of the latter cluster is only $\sim2\%$ on the {\it Eye of Horus}. Therefore, the influence associated with this cluster can be ignored.

The cold molecular gas in contemporary galaxies is structured in discrete cloud complexes. These giant molecular clouds (GMCs), with $10^4$-$10^7$ solar masses and radii of 5-100 parsecs, are the seeds of star formation. Highlighting the molecular gas structure at such small scales in distant galaxies is observationally challenging. Only a handful of molecular clouds were reported in two extreme submillimetre galaxies at high redshift. Here we search for GMCs in a typical Milky Way progenitor at z = 1.036. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we mapped the CO(4-3) emission of this gravitationally lensed galaxy at high resolution, reading down to 30 parsecs, which is comparable to the resolution of CO observations of nearby galaxies. We identify 17 molecular clouds, characterized by masses, surface densities and supersonic turbulence all of which are 10-100 times higher than present-day analogues. These properties question the universality of GMCs and suggest that GMCs inherit their properties from ambient interstellar medium. The measured cloud gas masses are similar to the masses of stellar clumps seen in the galaxy in comparable numbers. This corroborates the formation of molecular clouds by fragmentation of distant turbulent galactic gas disks, which then turn into stellar clumps ubiquitously observed in galaxies at cosmic noon.

We present the brightest cluster galaxies (BCGs) catalog for SPectroscoic IDentification of eROSITA Sources (SPIDERS) DR14 cluster program value-added catalog. We list the 416 BCGs identified as part of this process, along with their stellar mass, star formation rates, and morphological properties. We identified the BCGs based on the available spectroscopic data from SPIDERS and photometric data from SDSS. We computed stellar masses and SFRs of the BCGs on the basis of SDSS, WISE, and GALEX photometry using spectral energy distribution fitting. Morphological properties for all BCGs were derived by Sersic profile fitting using the software package SIGMA in different optical bands (g,r,i). We combined this catalog with the BCGs of galaxy groups and clusters extracted from the deeper AEGIS, CDFS, COSMOS, XMM-CFHTLS, and XMM-XXL surveys to study the stellar mass - halo mass relation using the largest sample of X-ray groups and clusters known to date. This result suggests that the mass growth of the central galaxy is controlled by the hierarchical mass growth of the host halo. We find a strong correlation between the stellar mass of BCGs and the mass of their host halos. This relation shows no evolution since z $\sim$ 0.65. We measure a mean scatter of 0.21 and 0.25 for the stellar mass of BCGs in a given halo mass at low ( $0.1<z < 0.3$ ) and high ( $0.3<z<0.65$ ) redshifts, respectively. We further demonstrate that the BCG mass is covariant with the richness of the host halos in the very X-ray luminous systems. We also find evidence that part of the scatter between X-ray luminosity and richness can be reduced by considering stellar mass as an additional variable.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

Recent instrumentation projects have allocated resources to develop codes for simulating astronomical images. Novel physics-based models are essential for understanding telescope, instrument, and environmental systematics in observations. A deep understanding of these systematics is especially important in the context of weak gravitational lensing, galaxy morphology, and other sensitive measurements. In this work, we present an adaptation of a physics-based ab initio image simulator: The Photon Simulator (PhoSim). We modify PhoSim for use with the Near-Infrared Camera (NIRCam) -- the primary imaging instrument aboard the James Webb Space Telescope (JWST). This photon Monte Carlo code replicates the observational catalog, telescope and camera optics, detector physics, and readout modes/electronics. Importantly, PhoSim-NIRCam simulates both geometric aberration and diffraction across the field of view. Full field- and wavelength-dependent point spread functions are presented. Simulated images of an extragalactic field are presented. Extensive validation is planned during in-orbit commissioning.

Using the Spitzer Space Telescope and Herschel Space Observatory, we have conducted a survey of infrared galaxies in the field of the galaxy cluster Abell S1063 (AS1063) at $z=0.347$, which is one of the most massive clusters known and a target of the HST CLASH and Frontier-Field surveys. The Spitzer/MIPS 24 $\mu$m and Herschel/PACS & SPIRE images revealed that the core of AS1063 is surprisingly devoid of infrared sources, showing only a few detectable sources within the central r$\sim1^{\prime}$. There is, however, one particularly bright source (2.3 mJy at 24 $\mu$m; 106 mJy at 160 $\mu$m), which corresponds to a background galaxy at $z=0.61$. The modest magnification factor (4.0$\times$) implies that this galaxy is intrinsically IR-luminous (L$_{\rm IR}=3.1\times10^{11}\ \rm L_{\odot}$). What is particularly interesting about this galaxy is that HST optical/near-infrared images show a remarkably bright and large (1 kpc) clump at one edge of the disk. Our follow-up optical/near-infrared spectroscopy shows Balmer (H$\alpha$-H8) and forbidden emission from this clump ([OII] $\lambda$3727, [OIII] $\lambda\lambda$4959,5007, [NII] $\lambda\lambda$6548,6583), indicating that it is a HII region. The HII region appears to have formed in-situ, as kinematically it is part of a rotating disk, and there is no evidence of nearby interacting galaxies. With an extinction correction of A$_{\rm V}=1.5$ mag, the star formation rate of this giant HII region is $\sim$10 M$_{\odot}$ yr$^{-1}$, which is exceptionally large, even for high redshift HII regions. Such a large and luminous HII region is often seen at $z\sim2$ but quite rare in the nearby Universe.

Galaxy formation depends on a complex interplay between gravitational collapse, gas accretion, merging, and feedback processes. Yet, after many decades of investigation, these concepts are poorly understood. This paper presents the argument that warm H$_2$ can be used as a tool to unlock some of these mysteries. Turbulence, shocks and outflows, driven by star formation, AGN activity or inflows, may prevent the rapid buildup of star formation in galaxies. Central to our understanding of how gas is converted into stars is the process by which gas can dissipate its mechanical energy through turbulence and shocks in order to cool. H$_2$ lines provide direct quantitative measurements of kinetic energy dissipation in molecular gas in galaxies throughout the Universe. Based on the detection of very powerful H$_2$ lines from z = 2 galaxies and proto-clusters at the detection limits of {\it Spitzer}, we are confident that future far-IR and UV H$_2$ observations will provide a wealth of new information and insight into galaxy evolution to high-z. Finally, at the very earliest epoch of star and galaxy formation, warm H$_2$ may also provide a unique glimpse of molecular gas collapse at 7 $<$ z $<$ 12 in massive dark matter (DM) halos on their way to forming the very first galaxies. Such measurements are beyond the reach of existing and planned observatories.

We report detections of two 1.2 mm continuum sources ($S_\mathrm{1.2mm}$ ~ 0.6 mJy) without any counterparts in the deep $H$- and/or $K$-band image (i.e., $K$-band magnitude $\gtrsim$ 26 mag). These near-infrared-dark faint millimeter sources are uncovered by ASAGAO, a deep and wide-field ($\simeq$ 26 arcmin$^2$) Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm survey. One has a red IRAC (3.6 and 4.5 $\mu$m) counterpart, and the other has been independently detected at 850 and 870 $\mu$m using SCUBA2 and ALMA Band 7, respectively. Their optical to radio spectral energy distributions indicate that they can lie at $z \gtrsim$ 3-5 and can be in the early phase of massive galaxy formation. Their contribution to the cosmic star formation rate density is estimated to be ~ 1 $\times$ 10$^{-3}$ $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ if they lie somewhere in the redshift range of $z$ ~ 3-5. This value can be consistent with, or greater than that of bright submillimeter galaxies ($S_\mathrm{870\mu m}>$ 4.2 mJy) at $z$ ~ 3-5. We also uncover 3 more candidates near-infrared-dark faint ALMA sources without any counterparts ($S_\mathrm{1.2mm}$ ~ 0.45-0.86 mJy). These results show that an unbiased ALMA survey can reveal the dust-obscured star formation activities, which were missed in previous deep optical/near-infrared surveys.

We quantify the star formation (SF) in the inner cores ($\mathcal{R}$/$R_{200}$$\leq$0.3) of 24 massive galaxy clusters at 0.2$\lesssim$$z$$\lesssim$0.9 observed by the $Herschel$ Lensing Survey and the Cluster Lensing and Supernova survey with $Hubble$. These programmes, covering the rest-frame ultraviolet to far-infrared regimes, allow us to accurately characterize stellar mass-limited ($\mathcal{M}_{*}$$>$$10^{10}$ $M_{\odot}$) samples of star-forming cluster members (not)-detected in the mid- and/or far-infrared. We release the catalogues with the photometry, photometric redshifts, and physical properties of these samples. We also quantify the SF displayed by comparable field samples from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. We find that in intermediate-$z$ cluster cores, the SF activity is suppressed with respect the field in terms of both the fraction ($\mathcal{F}$) of star-forming galaxies (SFG) and the rate at which they form stars ($\mathcal{SFR}$ and $s\mathcal{SFR} = \mathcal{SFR}/\mathcal{M}_{*}$). On average, the $\mathcal{F}$ of SFGs is a factor $\sim$$2$ smaller in cluster cores than in the field. Furthermore, SFGs present average $\mathcal{SFR}$ and $s\mathcal{SFR}$ typically $\sim$0.3 dex smaller in the clusters than in the field along the whole redshift range probed. Our results favour long time-scale quenching physical processes as the main driver of SF suppression in the inner cores of clusters since $z$$\sim$0.9, with shorter time-scale processes being very likely responsible for a fraction of the missing SFG population.

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI.

Deep MIR surveys have revealed numerous strongly star-forming galaxies at redshift z~<2. Their MIR fluxes are produced by a combination of continuum and PAH emission features. The PAH features can dominate the total MIR flux, but are difficult to measure without spectroscopy. We aim to study star-forming galaxies by using a blind spectroscopic survey at MIR wavelengths to understand evolution of their star formation rate (SFR) and specific SFR up to z~=0.5, by paying particular attention to their PAH properties. We conducted a low-resolution (R~=50) slitless spectroscopic survey at 5-13um of 9um flux-selected sources (>0.3 mJy) around the North Ecliptic Pole with the Infrared Camera (IRC) onboard AKARI. After removing 11 AGN candidates by using the IRC photometry, we identified 48 PAH galaxies with PAH 6.2, 7.7, and 8.6um features at z<0.5. The rest-frame optical-MIR spectral energy distributions (SEDs) based on CFHT and AKARI/IRC imaging covering 0.37-18um were produced, and analysed in conjunction with the PAH spectroscopy. We defined the PAH enhancement by using the luminosity ratio of the 7.7um PAH feature over the 3.5um stellar component of the SEDs. The rest-frame SEDs of all PAH galaxies have a universal shape with stellar and 7.7um bumps, except that the PAH enhancement significantly varies as a function of the PAH luminosities. We identified a PAH-enhanced population at z~>0.35, whose SEDs and luminosities are typical of luminous infrared galaxies. They show particularly larger PAH enhancement at high luminosity, implying that they are vigorous star-forming galaxies with elevated specific SFR. Our composite starburst model that combines a very young and optically very thick starburst with a very old population can successfully reproduce most of their SED characteristics, although we could not confirm this optically thick component from our spectral analysis.

We present the survey design, data reduction, construction of images, and source catalog of the Atacama Large Millimeter/submillimeter Array (ALMA) twenty-six arcmin^2 survey of GOODS-S at one-millimeter (ASAGAO). ASAGAO is a deep (1sigma ~ 61 uJy/beam for a 250 klambda-tapered map with a synthesized beam size of 0.51" x 0.45") and wide area (26 arcmin^2) survey on a contiguous field at 1.2 mm. By combining with ALMA archival data in the GOODS-South field, we obtained a deeper map in the same region (1sigma ~ 30 uJy/beam for a deep region with a 250 klambda-taper, and a synthesized beam size of 0.59" x 0.53"), providing the largest sample of sources (25 sources at >=5.0sigma, 45 sources at >=4.5sigma) among ALMA blank-field surveys to date. The number counts shows that 52(+11 -8)% of the extragalactic background light at 1.2 mm is resolved into discrete sources at S1.2m > 135 uJy. We create infrared (IR) luminosity functions (LFs) in the redshift range of z = 1-3 from the ASAGAO sources with KS-band counterparts, and constrain the faintest luminosity of the LF at 2.0 < z < 3.0. The LFs are consistent with previous results based on other ALMA and SCUBA-2 observations, which suggest a positive luminosity evolution and negative density evolution with increasing redshift. We find that obscured star-formation of sources with IR luminosities of log(L(IR)/Lsun)} ~> 11.8 account for ~~60%-90% of the z ~ 2 cosmic star-formation rate density.

We present an original phenomenological model to describe the evolution of galaxy number counts, morphologies, and spectral energy distributions across a wide range of redshifts (0.2<z<15) and stellar masses [Log10 M/Msun >6]. Our model follows observed mass and luminosity functions of both star-forming and quiescent galaxies, and reproduces the redshift evolution of colors, sizes, star-formation and chemical properties of the observed galaxy population. Unlike other existing approaches, our model includes a self-consistent treatment of stellar and photoionized gas emission and dust attenuation based on the BEAGLE tool. The mock galaxy catalogs generated with our new model can be used to simulate and optimize extragalactic surveys with future facilities such as the James Webb Space Telescope (JWST), and to enable critical assessments of analysis procedures, interpretation tools, and measurement systematics for both photometric and spectroscopic data. As a first application of this work, we make predictions for the upcoming JWST Advanced Deep Extragalactic Survey (JADES), a joint program of the JWST/NIRCam and NIRSpec Guaranteed Time Observations teams. We show that JADES will detect, with NIRCam imaging, thousands of galaxies at z>6, and tens at z>10 at m_AB<30 (5-sigma) within the 236 arcmin^2 of the survey. The JADES data will enable accurate constraints on the evolution of the UV luminosity function at z>8, and resolve the current debate about the rate of evolution of galaxies at z>8. Ready to use mock catalogs and software to generate new realizations are publicly available as the JAdes extraGalactic Ultradeep Artificial Realizations (JAGUAR) package.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We report fourteen and twenty-eight protocluster candidates at z=5.7 and 6.6 over 14 and 19 deg^2 areas, respectively, selected from 2,230 (259) Lya emitters (LAEs) photometrically (spectroscopically) identified with Subaru/Hyper Suprime-Cam (HSC) deep images (Keck, Subaru, and Magellan spectra and the literature data). Six out of the 42 protocluster candidates include 1-12 spectroscopically confirmed LAEs at redshifts up to z=6.574. By the comparisons with the cosmological Lya radiative transfer (RT) model reproducing LAEs with the reionization effects, we find that more than a half of these protocluster candidates are progenitors of the present-day clusters with a mass of > 10^14 M_sun. We then investigate the correlation between LAE overdensity delta and Lya rest-frame equivalent width EW_Lya^rest, because the cosmological Lya RT model suggests that a slope of EW_Lya^rest-delta relation is steepened towards the epoch of cosmic reionization (EoR), due to the existence of the ionized bubbles around galaxy overdensities easing the escape of Lya emission from the partly neutral intergalactic medium (IGM). The available HSC data suggest that the slope of the EW_Lya^rest-delta correlation does not evolve from the post-reionization epoch z=5.7 to the EoR z=6.6 beyond the moderately large statistical errors. There is a possibility that we would detect the evolution of the EW_Lya^rest - delta relation from z=5.7 to 7.3 by the upcoming HSC observations providing large samples of LAEs at z=6.6-7.3.

We report direct evidence of pre-processing of the galaxies residing in galaxy groups falling into galaxy clusters drawn from the Local Cluster Substructure Survey (LoCuSS). 34 groups have been identified via their X-ray emission in the infall regions of 23 massive ($\rm \langle M_{200}\rangle = 10^{15}\,M_{\odot}$) clusters at $0.15<z<0.3$. Highly complete spectroscopic coverage combined with 24 $\rm\mu$m imaging from Spitzer allows us to make a consistent and robust selection of cluster and group members including star forming galaxies down to a stellar mass limit of $\rm M_{\star} = 2\times10^{10}\,M_{\odot}$. The fraction $\rm f_{SF}$ of star forming galaxies in infalling groups is lower and with a flatter trend with respect to clustercentric radius when compared to the rest of the cluster galaxy population. At $\rm R\approx1.3\,r_{200}$ the fraction of star forming galaxies in infalling groups is half that in the cluster galaxy population. This is direct evidence that star formation quenching is effective in galaxies already prior to them settling in the cluster potential, and that groups are favourable locations for this process.

A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths.

Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the early phases of the dust-obscured Universe. Among the next generation of IR observatories, SPICA, observing in the 12-350 micron range, will be the only facility that can enable us to make the required leap forward in understanding the obscured star-formation rate and black-hole accretion rate densities (SFRD and BHARD, respectively) with respect to what Spitzer and Herschel achieved in the mid- and far-IR at z<3. In particular, SPICA will have the unique ability to trace the evolution of the obscured SFRD and BHARD over cosmic time, from the peak of their activity back to the reionisation epoch (i.e., 3<z<6-7), where its predecessors had severe limitations. Here we discuss the potential of both deep and shallow photometric surveys performed with the SPICA mid-IR instrument (SMI), enabled by the very low level of impact of dust obscuration in a band centred at 34 micron. These unique unbiased photometric surveys that SPICA will perform will be followed up by observations both with the SPICA spectrometers and with other facilities at shorter and longer wavelengths, with the aim to fully characterise the evolution of AGNs and star-forming galaxies after re-ionisation.

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed.

Galaxy clusters are expected to form hierarchically in a LCDM universe, growing primarily through mergers with lower mass clusters and the continual accretion of group-mass halos. Galaxy clusters assemble late, doubling their masses since z~0.5, and so the outer regions of clusters should be replete with infalling group-mass systems. We present an XMM-Newton survey to search for X-ray groups in the infall regions of 23 massive galaxy clusters at z~0.2, identifying 39 X-ray groups that have been spectroscopically confirmed to lie at the cluster redshift. These groups have mass estimates in the range 2x10^13-7x10^14Msun, and group-to-cluster mass ratios as low as 0.02. The comoving number density of X-ray groups in the infall regions is ~25x higher than that seen for isolated X-ray groups from the XXL survey. The average mass per cluster contained within these X-ray groups is 2.2x10^14Msun, or 19% of the mass within the primary cluster itself. We estimate that ~10^15Msun clusters increase their masses by 16% between z=0.223 and the present day due to the accretion of groups with M200>10^13.2Msun. This represents about half of the expected mass growth rate of clusters at these late epochs. The other half is likely to come from smooth accretion of matter not bound in halos. The mass function of the infalling X-ray groups appears significantly top-heavy with respect to that of field X-ray systems, consistent with expectations from numerical simulations, and the basic consequences of collapsed massive dark matter halos being biased tracers of the underlying large-scale density distribution.

We study a sample of 19 galaxy clusters in the redshift range $0.15<z<0.30$ with highly complete spectroscopic membership catalogues (to $K < K^{\ast}(\rm z)+1.5$) from the Arizona Cluster Redshift Survey (ACReS); individual weak-lensing masses and near-infrared data from the Local Cluster Substructure Survey (LoCuSS); and optical photometry from the Sloan Digital Sky Survey (SDSS). We fit the scaling relations between total cluster luminosity in each of six bandpasses (${\it grizJK}$) and cluster mass, finding cluster luminosity to be a promising mass proxy with low intrinsic scatter $\sigma_{\ln L|M}$ of only $\sim 10-20$ per cent for all relations. At fixed overdensity radius, the intercept increases with wavelength, consistent with an old stellar population. The scatter and slope are consistent across all wavelengths, suggesting that cluster colour is not a function of mass. Comparing colour with indicators of the level of disturbance in the cluster, we find a narrower variety in the cluster colours of 'disturbed' clusters than of 'undisturbed' clusters. This trend is more pronounced with indicators sensitive to the initial stages of a cluster merger, e.g. the Dressler Schectman statistic. We interpret this as possible evidence that the total cluster star formation rate is 'standardised' in mergers, perhaps through a process such as a system-wide shock in the intracluster medium.

We report on the galaxy MACSJ0032-arc at z=3.6314 discovered during the Herschel Lensing snapshot Survey of massive galaxy clusters, and strongly lensed by the cluster MACSJ0032.1+1808. The successful detections of its rest-frame UV, optical, FIR, millimeter, and radio continua, and of its CO emission enable us to characterize, for the first time at such a high redshift, the stellar, dust, and molecular gas properties of a compact star-forming galaxy with a size smaller than 2.5 kpc, a fairly low stellar mass of 4.8x10^9 Msun, and a moderate IR luminosity of 4.8x10^11 Lsun. We find that the bulk of the molecular gas mass and star formation seems to be spatially decoupled from the rest-frame UV emission. About 90% of the total star formation rate is undetected at rest-frame UV wavelengths because of severe obscuration by dust, but is seen through the thermal FIR dust emission and the radio synchrotron radiation. The observed CO(4-3) and CO(6-5) lines demonstrate that high-J transitions, at least up to J=6, remain excited in this galaxy, whose CO spectral line energy distribution resembles that of high-redshift submm galaxies, even though the IR luminosity of MACSJ0032-arc is ten times lower. This high CO excitation is possibly due to the compactness of the galaxy. We find evidence that this high CO excitation has to be considered in the balance when estimating the CO-to-H2 conversion factor. The inferred depletion time of the molecular gas in MACSJ0032-arc supports the decrease in the gas depletion timescale of galaxies with redshift, although to a lesser degree than predicted by galaxy evolution models. Instead, the measured molecular gas fraction as high as 60-79% in MACSJ0032-arc favors the continued increase in the gas fraction of galaxies with redshift as expected, despite the plateau observed between z~1.5 and z~2.5.

Measurement of the gas velocity distribution in galaxy clusters provides insight into the physics of mergers, through which large scale structures form in the Universe. Velocity estimates within the intracluster medium (ICM) can be obtained via the Sunyaev-Zel'dovich (SZ) effect, but its observation is challenging both in term of sensitivity requirement and control of systematic effects, including the removal of contaminants. In this paper we report resolved observations, at 150 and 260 GHz, of the SZ effect toward the triple merger MACS J0717.5+3745 (z=0.55), using data obtained with the NIKA camera at the IRAM 30m telescope. Assuming that the SZ signal is the sum of a thermal (tSZ) and a kinetic (kSZ) component and by combining the two NIKA bands, we extract for the first time a resolved map of the kSZ signal in a cluster. The kSZ signal is dominated by a dipolar structure that peaks at -5.1 and +3.4 sigma, corresponding to two subclusters moving respectively away and toward us and coincident with the cold dense X-ray core and a hot region undergoing a major merging event. We model the gas electron density and line-of-sight velocity of MACS J0717.5+3745 as four subclusters. Combining NIKA data with X-ray observations from XMM-Newton and Chandra, we fit this model to constrain the gas line-of-sight velocity of each component, and we also derive, for the first time, a velocity map from kSZ data (i.e. that is model-dependent). Our results are consistent with previous constraints on the merger velocities, and thanks to the high angular resolution of our data, we are able to resolve the structure of the gas velocity. Finally, we investigate possible contamination and systematic effects with a special care given to radio and submillimeter galaxies. Among the sources that we detect with NIKA, we find one which is likely to be a high redshift lensed submillimeter galaxy.

The James Webb Space Telescope near-infrared camera (JWST NIRCam) has two 2.'2 $\times$ 2.'2 fields of view that are capable of either imaging or spectroscopic observations. Either of two $R \sim 1500$ grisms with orthogonal dispersion directions can be used for slitless spectroscopy over $\lambda = 2.4 - 5.0$ $\mu$m in each module, and shorter wavelength observations of the same fields can be obtained simultaneously. We present the latest predicted grism sensitivities, saturation limits, resolving power, and wavelength coverage values based on component measurements, instrument tests, and end-to-end modeling. Short wavelength (0.6 -- 2.3 $\mu$m) imaging observations of the 2.4 -- 5.0 $\mu$m spectroscopic field can be performed in one of several different filter bands, either in-focus or defocused via weak lenses internal to NIRCam. Alternatively, the possibility of 1.0 -- 2.0 $\mu$m spectroscopy (simultaneously with 2.4 -- 5.0 $\mu$m) using dispersed Hartmann sensors (DHSs) is being explored. The grisms, weak lenses, and DHS elements were included in NIRCam primarily for wavefront sensing purposes, but all have significant science applications. Operational considerations including subarray sizes, and data volume limits are also discussed. Finally, we describe spectral simulation tools and illustrate potential scientific uses of the grisms by presenting simulated observations of deep extragalactic fields, galactic dark clouds, and transiting exoplanets.

We report the serendipitous discovery of HSC J142449-005322, a double source plane lens system in the Hyper Suprime-Cam Subaru Strategic Program. We dub the system Eye of Horus. The lens galaxy is a very massive early-type galaxy with stellar mass of ~7x10^11 Msun located at z_L=0.795. The system exhibits two arcs/rings with clearly different colors, including several knots. We have performed spectroscopic follow-up observations of the system with FIRE on Magellan. The outer ring is confirmed at z_S2=1.988 with multiple emission lines, while the inner arc and counterimage is confirmed at z_S1=1.302. This makes it the first double source plane system with spectroscopic redshifts of both sources. Interestingly, redshifts of two of the knots embedded in the outer ring are found to be offset by delta_z=0.002 from the other knots, suggesting that the outer ring consists of at least two distinct components in the source plane. We perform lens modeling with two independent codes and successfully reproduce the main features of the system. However, two of the lensed sources separated by ~0.7 arcsec cannot be reproduced by a smooth potential, and the addition of substructure to the lens potential is required to reproduce them. Higher-resolution imaging of the system will help decipher the origin of this lensing feature and potentially detect the substructure.

The far-infrared fine-structure line [CII] at 1900.5\,GHz is known to be one of the brightest cooling lines in local galaxies, and therefore it has been suggested to be an efficient tracer for star-formation in very high-redshift galaxies. However, recent results for galaxies at $z>6$ have yielded numerous non-detections in star-forming galaxies, except for quasars and submillimeter galaxies. We report the results of ALMA observations of two lensed, star-forming galaxies at $z = 6.029$ and $z=6.703$. The galaxy A383-5.1 (star formation rate [SFR] of 3.2 M$_\odot$ yr$^{-1}$ and magnification of $\mu = 11.4\pm1.9$) shows a line detection with $L_{\rm [CII]} = 8.9\times10^{6}$ L$_\odot$, making it the lowest $L_{\rm [CII]}$ detection at $z>6$. For MS0451-H (SFR = 0.4 M$_\odot$ yr$^{-1}$ and $\mu = 100\pm20$) we provide an upper limit of $L_{\rm [CII]} < 3\times10^{5}$ L$_\odot$, which is 1\,dex below the local SFR-$L_{\rm [CII]}$ relations. The results are consistent with predictions for low-metallicity galaxies at $z>6$, however, other effects could also play a role in terms of decreasing $L_{\rm [CII]}$. The detection of A383-5.1 is encouraging and suggests that detections are possible, but much fainter than initially predicted.

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $\Lambda$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

We present a complete census of all Herschel-detected sources within the six massive lensing clusters of the HST Frontier Fields (HFF). We provide a robust legacy catalogue of 263 sources with Herschel fluxes, primarily based on imaging from the Herschel Lensing Survey (HLS) and PEP/HerMES Key Programmes. We optimally combine Herschel, Spitzer and WISE infrared (IR) photometry with data from HST, VLA and ground-based observatories, identifying counterparts to gain source redshifts. For each Herschel-detected source we also present magnification factor (mu), intrinsic IR luminosity and characteristic dust temperature, providing a comprehensive view of dust-obscured star formation within the HFF. We demonstrate the utility of our catalogues through an exploratory overview of the magnified population, including more than 20 background sub-LIRGs unreachable by Herschel without the assistance gravitational lensing.

We present a high-precision mass model of the galaxy cluster MACSJ1149.6+2223, based on a strong-gravitational-lensing analysis of Hubble Space Telescope Frontier Fields (HFF) imaging data and spectroscopic follow-up with Gemini/GMOS and VLT/MUSE. Our model includes 12 new multiply imaged galaxies, bringing the total to 22, comprised of 65 individual lensed images. Unlike the first two HFF clusters, Abell 2744 and MACSJ0416.1-2403, MACSJ1149 does not reveal as many multiple images in the HFF data. Using the Lenstool software package and the new sets of multiple images, we model the cluster with several cluster-scale dark-matter halos and additional galaxy-scale halos for the cluster members. Consistent with previous analyses, we find the system to be complex, composed of five cluster-scale halos. Their spatial distribution and lower mass, however, makes MACSJ1149 a less powerful lens. Our best-fit model predicts image positions with an RMS of 0.91". We measure the total projected mass inside a 200~kpc aperture as ($1.840\pm 0.006$)$\times 10^{14}$M$_{\odot}$, thus reaching again 1% precision, following our previous HFF analyses of MACSJ0416.1-2403 and Abell 2744. In light of the discovery of the first resolved quadruply lensed supernova, SN Refsdal, in one of the multiply imaged galaxies identified in MACSJ1149, we use our revised mass model to investigate the time delays and predict the rise of the next image between November 2015 and January 2016.

We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$\alpha$ emitters (LAEs) at $5.7<z<7$. These galaxies have deep optical and infrared images from Subaru, $HST$, and $Spitzer$/IRAC. We focus on a subset of 27 galaxies with IRAC detections, and characterize their stellar populations utilizing galaxy synthesis models based on the multi-band data and secure redshifts. By incorporating nebular emission estimated from the observed Ly$\alpha$ flux, we are able to break the strong degeneracy of model spectra between young galaxies with prominent nebular emission and older galaxies with strong Balmer breaks. The results show that our galaxies cover a wide range of ages from several to a few hundred million years (Myr), and a wide range of stellar masses from $\sim10^8$ to $\sim10^{11}$ $M_{\odot}$. These galaxies can be roughly divided into an `old' subsample and a `young' subsample. The `old' subsample consists of galaxies older than 100 Myr, with stellar masses higher than $10^9$ $M_{\odot}$. The galaxies in the `young' subsample are younger than $\sim$30 Myr, with masses ranging between $\sim10^8$ and $\sim3\times10^9$ $M_{\odot}$. Both subsamples display a correlation between stellar mass and star-formation rate (SFR), but with very different normalizations. The average specific SFR (sSFR) of the `old' subsample is 3--4 Gyr$^{-1}$, consistent with previous studies of `normal' star-forming galaxies at $z\ge6$. The average sSFR of the `young' subsample is an order of magnitude higher, likely due to starburst activity. Our results also indicate little or no dust extinction in the majority of the galaxies, as already suggested by their steep rest-frame UV slopes. Finally, LAEs and LBGs with strong Ly$\alpha$ emission are indistinguishable in terms of age, stellar mass, and SFR.

We use deep Hubble Space Telescope imaging of the Frontier Fields to accurately measure the galaxy rest-frame ultraviolet luminosity function (UV LF) in the redshift range $z \sim 6-8$. We combine observations in three lensing clusters A2744, MACS0416, MACS0717 and their associated parallels fields to select high-redshift dropout candidates. We use the latest lensing models to estimate the flux magnification and the effective survey volume in combination with completeness simulations performed in the source plane. We report the detection of 227 galaxy candidates at $z=6-7$ and 25 candidates at $z \sim 8$. While the total survey area is about 4 arcmin$^{2}$ in each parallel field, it drops to about 0.6 to 1 arcmin$^{2}$ in the cluster core fields because of the strong lensing. We compute the UV luminosity function at $z \sim 7$ using the combined galaxy sample and perform Monte Carlo simulations to determine the best fit Schechter parameters. We are able to reliably constrain the LF down to an absolute magnitude of $M_{UV}=-15.25$, which corresponds to 0.005$L^{\star}$. More importantly, we find that the faint-end slope remains steep down to this magnitude limit with $\alpha=-2.04_{-0.17}^{+0.13}$. Our results confirm the most recent results in deep blank fields but extend the LF measurements more than two magnitudes deeper. The UV LF at $z \sim 8$ is not very well constrained below $M_{UV}=-18$ due to the small number statistics and incompleteness uncertainties. To assess the contribution of galaxies to cosmic reionization we derive the UV luminosity density at $z\sim7$ by integrating the UV LF down to an observationally constrained limit of $M_{UV} = -15$. We show that our determination of Log($\rho_{UV}$)=$26.2\pm0.13$ (erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$) can be sufficient to maintain the IGM ionized.

Two millimeter observations of the MACS J1149.6+2223 cluster have detected a source that was consistent with the location of the lensed MACS1149-JD galaxy at z=9.6. A positive identification would have rendered this galaxy as the youngest dust forming galaxy in the universe. Follow up observation with the AzTEC 1.1 mm camera and the IRAM NOrthern Extended Millimeter Array (NOEMA) at 1.3 mm have not confirmed this association. In this paper we show that the NOEMA observations associate the 2 mm source with [PCB2012] 2882 ([PCB2012] 2882 is the NED-searchable name for this source.), source number 2882 in the Hubble Space Telescope ( HST) Cluster Lensing and Supernova (CLASH) catalog of MACS J1149.6+2223. This source, hereafter referred to as CLASH 2882, is a gravitationally lensed spiral galaxy at z=0.99. We combine the GISMO 2 mm and NOEMA 1.3 mm fluxes with other (rest frame) UV to far-IR observations to construct the full spectral energy distribution (SED) of this galaxy, and derive its star formation history, and stellar and interstellar dust content. The current star formation rate of the galaxy is 54/mu Msun yr-1, and its dust mass is about 5 10^7/mu Msun, where mu is the lensing magnification factor for this source, which has a mean value of 2.7. The inferred dust mass is higher than the maximum dust mass that can be produced by core collapse supernovae (CCSN) and evolved AGB stars. As with many other star forming galaxies, most of the dust mass in CLASH 2882 must have been accreted in the dense phases of the ISM.

We present a sample of 290 24-micron-selected active galactic nuclei (AGNs) mostly at z ~ 0.3 -- 2.5, within 5.2 square degrees distributed as 25' X 25' fields around each of 30 galaxy clusters in the Local Cluster Substructure Survey (LoCuSS). The sample is nearly complete to 1 mJy at 24 microns, and has a rich multi-wavelength set of ancillary data; 162 are detected by Herschel. We use spectral templates for AGNs, stellar populations, and infrared emission by star forming galaxies to decompose the spectral energy distributions (SEDs) of these AGNs and their host galaxies, and estimate their star formation rates (SFRs), AGN luminosities, and host galaxy stellar masses. The set of templates is relatively simple: a standard Type-1 quasar template; another for the photospheric output of the stellar population; and a far infrared star-forming template. For the Type-2 AGN SEDs, we substitute templates including internal obscuration, and some Type-1 objects require a warm component (T > 50 K). The individually Herschel- detected Type-1 AGNs and a subset of 17 Type-2 ones typically have luminosities > 10^{45} ergs/s, and supermassive black holes of ~ 3 X 10^8 Msun emitting at ~ 10% of the Eddington rate. We find them in about twice the numbers of AGN identified in SDSS data in the same fields, i.e., they represent typical high luminosity AGN, not an infrared-selected minority. These AGNs and their host galaxies are studied further in an accompanying paper.

We study the relation of active galactic nuclei (AGNs) to star formation in their host galaxies. Our sample includes 205 Type-1 and 85 Type-2 AGNs, 162 detected with Herschel, from fields surrounding 30 galaxy clusters in the Local Cluster Substructure Survey (LoCuSS). The sample is identified by optical line widths and ratios after selection to be brighter than 1 mJy at 24 microns. We show that Type-2 AGN [OIII]5007 line fluxes at high z can be contaminated by their host galaxies with typical spectrograph entrance apertures (but our sample is not compromised in this way). We use spectral energy distribution (SED) templates to decompose the galaxy SEDs and estimate star formation rates, AGN luminosities, and host galaxy stellar masses (described in an accompanying paper). The AGNs arise from massive black holes (~ 3 X 10^8 Msun) accreting at ~ 10% of the Eddington rate and residing in galaxies with stellar mass > 3 X 10^{10} Msun; those detected with Herschel have IR luminosity from star formation in the range of 10^{10} -- 10^{12} Lsun. We find that: 1.) the specific star formation rates in the host galaxies are generally consistent with those of normal star-forming (main sequence) galaxies; 2.) there is a strong correlation between the luminosities from star formation and the AGN; and 3.) however, the correlation may not result from a causal connection, but could arise because the black hole mass (and hence AGN Eddington luminosity) and star formation are both correlated with the galaxy mass.

Cosmological constraints from galaxy clusters rely on accurate measurements of the mass and internal structure of clusters. An important source of systematic uncertainty in cluster mass and structure measurements is the secure selection of background galaxies that are gravitationally lensed by clusters. This issue has been shown to be particular severe for faint blue galaxies. We therefore explore the selection of faint blue background galaxies, by reference to photometric redshift catalogs derived from the COSMOS survey and our own observations of massive galaxy clusters at z~0.2. We show that methods relying on photometric redshifts of galaxies in/behind clusters based on observations through five filters, and on deep 30-band COSMOS photometric redshifts are both inadequate to identify safely faint blue background galaxies. This is due to the small number of filters used by the former, and absence of massive galaxy clusters at redshifts of interest in the latter. We therefore develop a pragmatic method to combine both sets of photometric redshifts to select a population of blue galaxies based purely on photometric analysis. This sample yields stacked weak-lensing results consistent with our previously published results based on red galaxies. We also show that the stacked clustercentric number density profile of these faint blue galaxies is consistent with expectations from consideration of the lens magnification signal of the clusters. Indeed, the observed number density of blue background galaxies changes by ~10-30 per cent across the radial range over which other surveys assume it to be flat.

We present Early Science observations with the Large Millimeter Telescope, AzTEC 1.1 mm continuum images and wide bandwidth spectra (73-111 GHz) acquired with the Redshift Search Receiver, towards four bright lensed submillimetre galaxies identified through the Herschel Lensing Survey-snapshot and the SCUBA-2 Cluster Snapshot Survey. This pilot project studies the star formation history and the physical properties of the molecular gas and dust content of the highest redshift galaxies identified through the benefits of gravitational magnification. We robustly detect dust continuum emission for the full sample and CO emission lines for three of the targets. We find that one source shows spectroscopic multiplicity and is a blend of three galaxies at different redshifts (z=2.040, 3.252 and 4.680), reminiscent of previous high-resolution imaging follow-up of unlensed submillimetre galaxies, but with a completely different search method, that confirm recent theoretical predictions of physically unassociated blended galaxies. Identifying the detected lines as 12CO (J_up=2-5) we derive spectroscopic redshifts, molecular gas masses, and dust masses from the continuum emission. The mean H_2 gas mass of the full sample is (2.0 +- 0.2) x 10^11 M_sun/\mu, and the mean dust mass is (2.0+-0.2) x 10^9 M_sun/\mu, where \mu=2-5 is the expected lens amplification. Using these independent estimations we infer a gas-to-dust ratio of \delta_GDR=55-75, in agreement with other measurements of submillimetre galaxies. Our magnified high-luminosity galaxies fall on the same locus as other high-redshift submillimetre galaxies, extending the L'_CO - L_FIR correlation observed for local luminous and ultraluminous infrared galaxies to higher FIR and CO luminosities.

We present a high-precision mass model of galaxy cluster Abell 2744, based on a strong-gravitational-lensing analysis of the \emph{Hubble Space Telescope Frontier Fields} (HFF) imaging data, which now include both \emph{Advanced Camera for Surveys} and \emph{Wide-Field Camera 3} observations to the final depth. Taking advantage of the unprecedented depth of the visible and near-infrared data, we identify 34 new multiply imaged galaxies, bringing the total to 61, comprising 181 individual lensed images. In the process, we correct previous erroneous identifications and positions of multiple systems in the northern part of the cluster core. With the \textsc{Lenstool} software and the new sets of multiple images, we model the cluster using two cluster-scale dark matter halos plus galaxy-scale halos for the cluster members. Our best-fit model predicts image positions with an \emph{RMS} error of 0.69$\arcsec$, which constitutes an improvement by almost a factor of two over previous parametric models of this cluster. We measure the total projected mass inside a 200~kpc aperture as ($2.162\pm 0.005$)$\times 10^{14}M_{\odot}$, thus reaching 1\% level precision for the second time, following the recent HFF measurement of MACSJ0416.1-2403. Importantly, the higher quality of the mass model translates into an overall improvement by a factor of 4 of the derived magnification factor. % for the high-redshift lensed background galaxies. Together with our previous HFF gravitational lensing analysis, this work demonstrates that the HFF data enables high-precision mass measurements for massive galaxy clusters and the derivation of robust magnification maps to probe the early Universe.

The Spitzer-Cosmic Assembly Deep Near-Infrared Extragalactic Legacy Survey (S-CANDELS; PI G. Fazio) is a Cycle 8 Exploration Program designed to detect galaxies at very high redshifts (z > 5). To mitigate the effects of cosmic variance and also to take advantage of deep coextensive coverage in multiple bands by the Hubble Space Telescope Multi-Cycle Treasury Program CANDELS, S-CANDELS was carried out within five widely separated extragalactic fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the HST Deep Field North, and the Extended Groth Strip. S-CANDELS builds upon the existing coverage of these fields from the Spitzer Extended Deep Survey (SEDS) by increasing the integration time from 12 hours to a total of 50 hours but within a smaller area, 0.16 square degrees. The additional depth significantly increases the survey completeness at faint magnitudes. This paper describes the S-CANDELS survey design, processing, and publicly-available data products. We present IRAC dual-band 3.6+4.5 micron catalogs reaching to a depth of 26.5 AB mag. Deep IRAC counts for the roughly 135,000 galaxies detected by S-CANDELS are consistent with models based on known galaxy populations. The increase in depth beyond earlier Spitzer/IRAC surveys does not reveal a significant additional contribution from discrete sources to the diffuse Cosmic Infrared Background (CIB). Thus it remains true that only roughly half of the estimated CIB flux from COBE/DIRBE is resolved.

We present a study of the spatial distribution and kinematics of star-forming galaxies in 30 massive clusters at 0.15<z<0.30, combining wide-field Spitzer 24um and GALEX NUV imaging with highly-complete spectroscopy of cluster members. The fraction (f_SF) of star-forming cluster galaxies rises steadily with cluster-centric radius, increasing fivefold by 2r200, but remains well below field values even at 3r200. This suppression of star formation at large radii cannot be reproduced by models in which star formation is quenched in infalling field galaxies only once they pass within r200 of the cluster, but is consistent with some of them being first pre-processed within galaxy groups. Despite the increasing f_SF-radius trend, the surface density of star-forming galaxies actually declines steadily with radius, falling ~15x from the core to 2r200. This requires star-formation to survive within recently accreted spirals for 2--3Gyr to build up the apparent over-density of star-forming galaxies within clusters. The velocity dispersion profile of the star-forming galaxy population shows a sharp peak of 1.44-sigma_v at 0.3r500, and is 10--35% higher than that of the inactive cluster members at all cluster-centric radii, while their velocity distribution shows a flat, top-hat profile within r500. All of these results are consistent with star-forming cluster galaxies being an infalling population, but one that must also survive ~0.5--2Gyr beyond passing within r200. By comparing the observed distribution of star-forming galaxies in the stacked caustic diagram with predictions from the Millennium simulation, we obtain a best-fit model in which SFRs decline exponentially on quenching time-scales of 1.73\pm0.25 Gyr upon accretion into the cluster.

We present weak lensing and X-ray analysis of 12 low mass clusters from the CFHTLenS and XMM-CFHTLS surveys. We combine these systems with high-mass systems from CCCP and low-mass systems from COSMOS to obtain a sample of 70 systems, spanning over two orders of magnitude in mass. We measure core-excised Lx-Tx, M-Lx and M-Tx scaling relations and include corrections for observational biases. By providing fully bias corrected relations, we give the current limitations for Lx and Tx as cluster mass proxies. We demonstrate that Tx benefits from a significantly lower intrinsic scatter at fixed mass than Lx. By studying the residuals of the bias corrected relations, we show for the first time using weak lensing masses that galaxy groups seem more luminous and warmer for their mass than clusters. This implies a steepening of the M-Lx and M-Tx relations at low masses. We verify the inferred steepening using a different high mass sample from the literature and show that variance between samples is the dominant effect leading to discrepant scaling relations. We divide our sample into subsamples of merging and relaxed systems, and find that mergers may have enhanced scatter in lensing measurements, most likely due to stronger triaxiality and more substructure. For the Lx-Tx relation, which is unaffected by lensing measurements, we find the opposite trend in scatter. We also explore the effects of X-ray cross-calibration and find that Chandra calibration leads to flatter Lx-Tx and M-Tx relations than XMM-Newton.

Our objectives are to determine the properties of the interstellar medium (ISM) and of star-formation in typical star-forming galaxies at high redshift. Following up on our previous multi-wavelength observations with HST, Spitzer, Herschel, and the Plateau de Bure Interferometer (PdBI), we have studied a strongly lensed z=2.013 galaxy, the arc behind the galaxy cluster MACS J0451+0006, with ALMA to measure the [CII] 158 micron emission line, one of the main coolants of the ISM. [CII] emission from the southern part of this galaxy is detected at 10 $\sigma$. Taking into account strong gravitational lensing, which provides a magnification of $\mu=49$, the intrinsic lensing-corrected [CII]158 micron luminosity is $L(CII)=1.2 \times 10^8 L_\odot$. The observed ratio of [CII]-to-IR emission, $L(CII)/L(FIR) \approx (1.2-2.4) \times 10^{-3}$, is found to be similar to that in nearby galaxies. The same also holds for the observed ratio $L(CII)/L(CO)=2.3 \times 10^3$, which is comparable to that of star-forming galaxies and active galaxy nuclei (AGN) at low redshift. We utilize strong gravitational lensing to extend diagnostic studies of the cold ISM to an order of magnitude lower luminosity ($L(IR) \sim (1.1-1.3) \times 10^{11} L_\odot$) and SFR than previous work at high redshift. While larger samples are needed, our results provide evidence that the cold ISM of typical high redshift galaxies has physical characteristics similar to normal star forming galaxies in the local Universe.

To extend the molecular gas measurements to typical star-forming galaxies (SFGs) with SFR < 40 Msun yr^{-1} and M* < 2.5x10^{10} Msun at z~1.5-3, we have observed CO emission with the IRAM Plateau de Bure Interferometer and 30m telescope for five strongly-lensed galaxies selected from the Herschel Lensing Survey. These observations are combined with a compilation of CO measurements from the literature. We infer the luminosity correction factors r2,1 = 0.81+/-0.20 and r3,1 = 0.57+/-0.15 for the J=2 and J=3 CO transitions, respectively, valid for SFGs at z>1. The combined sample of CO-detected SFGs at z>1 shows a large spread in star formation efficiency (SFE), such that SFE extend beyond the low values of local spirals and overlap the distribution of z>1 sub-mm galaxies. We find that the spread in SFE (or equivalently in molecular gas depletion timescale) is due to primarily the specific star formation rate, but also stellar mass and redshift. Correlations of SFE with the offset from the main-sequence and the compactness of the starburst are less clear. The increase of the molecular gas depletion timescale with M* now revealed by low M* SFGs at z>1 and observed at z=0 is in contrast to the admitted constant molecular gas depletion timescale and the linear Kennicutt-Schmidt relation. We confirm an increase of the molecular gas fraction (fgas) from z~0.2 to z~1.2, followed by a very mild increase toward higher redshifts. At each redshift fgas shows a large dispersion due to the dependence of fgas on M*, producing a gradient of increasing fgas with decreasing M*. We provide the first measure of fgas of z>1 SFGs at the low-M* end (10^{9.4} < M*/Msun < 10^{9.9}), reaching a mean fgas = 0.69+/-0.18, which shows a clear fgas upturn. Finally, we find evidence for a non-universal dust-to-gas ratio among high-redshift SFGs and sub-mm galaxies, local spirals and ULIRGs with near-solar metallicities.

The CFHTLS presents a unique data set for weak lensing studies, having high quality imaging and deep multi-band photometry. We have initiated an XMM-CFHTLS project to provide X-ray observations of the brightest X-ray selected clusters within the wide CFHTLS area. Performance of these observations and the high quality of CFHTLS data, allows us to revisit the identification of X-ray sources, introducing automated reproducible algorithms, based on the multi-color red sequence finder. We have also introduced a new optical mass proxy. We provide the calibration of the red sequence observed in the CFHT filters and compare the results with the traditional single color red sequence and photoz. We test the identification algorithm on the subset of highly significant XMM clusters and identify 100% of the sample. We find that the integrated z-band luminosity of the red sequence galaxies correlates well with the X-ray luminosity with a surprisingly small scatter of 0.20 dex. We further use the multi-color red sequence to reduce spurious detections in the full XMM and RASS data sets, resulting in catalogs of 196 and 32 clusters, respectively. We made spectroscopic follow-up observations of some of these systems with HECTOSPEC and in combination with BOSS DR9 data. We also describe the modifications needed to the source detection algorithm in order to keep high purity of extended sources in the shallow X-ray data. We also present the scaling relation between X-ray luminosity and velocity dispersion.

We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dynamics, we complement this lensing analysis with a study of the intra-cluster gas using shallow Chandra data, and a three-dimensional model of the distribution and motions of cluster galaxies derived from over 100 spectroscopic redshifts. The multi-scale grid model obtained from our combined lensing analysis extends the high-precision strong-lensing mass reconstruction recently performed to cluster-centric distances of almost 1 Mpc. Our analysis detects the two well known mass concentrations in the cluster core. A pronounced offset between collisional and collisionless matter is only observed for the SW cluster component, while excellent alignment is found for the NE cluster. Both the lensing analysis and the distribution of cluster light strongly suggest the presence of a third massive structure, almost 2 arcmin SW of the cluster centre. Since no X-ray emission is detected in this region, we conclude that this structure is non-virialised and speculate that it might be part of a large-scale filament almost aligned with our line of sight. Combining all evidence from the distribution of dark and luminous matter, we propose two alternative scenarios for the trajectories of the components of MACSJ0416.1-2403. Upcoming deep X-ray observations that allow the detection of shock fronts, cold cores, and sloshing gas (all key diagnostics for studies of cluster collisions) will allow us to test, and distinguish between these two scenarios.

Exploiting the power of gravitational lensing, the Hubble Frontier Fields (HFF) program aims at observing six massive galaxy clusters to explore the distant Universe far beyond the depth limits of blank field surveys. Using the complete Hubble Space Telescope observations of the first HFF cluster Abell 2744, we report the detection of 50 galaxy candidates at $z \sim 7$ and eight candidates at $z \sim 8$ in a total survey area of 0.96 arcmin$^{2}$ in the source plane. Three of these galaxies are multiply-imaged by the lensing cluster. Using an updated model of the mass distribution in the cluster we were able to calculate the magnification factor and the effective survey volume for each galaxy in order to compute the ultraviolet galaxy luminosity function at both redshifts 7 and 8. Our new measurements extend the $z \sim 7$ UV LF down to an absolute magnitude of $M_{UV} \sim -15.5$. We find a characteristic magnitude of $M^{\star}_{UV}=-20.63^{+0.69}_{-0.56}$ mag and a faint-end slope $\alpha = -1.88^{+0.17}_{-0.20}$ close to previous determinations in blank fields. We show here for the first time that this slope remains steep down to very faint luminosities of 0.01$L^{\star}$. Although prone to large uncertainties, our results at $z \sim 8$ also seem to confirm a steep faint-end slope below 0.1$L^{\star}$. The HFF program is therefore providing an extremely efficient way to study the faintest galaxy populations at $z > 7$ that would otherwise be inaccessible with current instrumentation. The full sample of six galaxy clusters will provide yet better constraints on the build-up of galaxies at early epochs and their contribution to cosmic reionization.

We present ALMA observations of the [CII] line and far-infrared (FIR) continuum of a normally star-forming galaxy in the reionization epoch, the z=6.96 Ly-alpha emitter (LAE) IOK-1. Probing to sensitivities of sigma_line = 240 micro-Jy/beam (40 km/s channel) and sigma_cont = 21 micro-Jy/beam, we found the galaxy undetected in both [CII] and continuum. Comparison of UV - FIR spectral energy distribution (SED) of IOK-1, including our ALMA limit, with those of several types of local galaxies (including the effects of the cosmic microwave background, CMB, on the FIR continuum) suggests that IOK-1 is similar to local dwarf/irregular galaxies in SED shape rather than highly dusty/obscured galaxies. Moreover, our 3 sigma FIR continuum limit, corrected for CMB effects, implies intrinsic dust mass M_dust < 6.4 x 10^7 M_sun, FIR luminosity L_FIR < 3.7 x 10^{10} L_sun (42.5 - 122.5 micron), total IR luminosity L_IR < 5.7 x 10^{10} L_sun (8 - 1000 micron) and dust-obscured star formation rate (SFR) < 10 M_sun/yr, if we assume that IOK-1 has a dust temperature and emissivity index typical of local dwarf galaxies. This SFR is 2.4 times lower than one estimated from the UV continuum, suggesting that < 29% of the star formation is obscured by dust. Meanwhile, our 3 sigma [CII] flux limit translates into [CII] luminosity, L_[CII] < 3.4 x 10^7 L_sun. Locations of IOK-1 and previously observed LAEs on the L_[CII] vs. SFR and L_[CII]/L_FIR vs. L_FIR diagrams imply that LAEs in the reionization epoch have significantly lower gas and dust enrichment than AGN-powered systems and starbursts at similar/lower redshifts, as well as local star-forming galaxies.

Extending over three Hubble Space Telescope (HST) cycles, the Hubble Frontier Fields (HFF) initiative constitutes the largest commitment ever of HST time to the exploration of the distant Universe via gravitational lensing by massive galaxy clusters. We here present models of the mass distribution in the six HFF cluster lenses, derived from a joint strong- and weak-lensing analysis anchored by a total of 88 multiple-image systems identified in existing HST data. The resulting maps of the projected mass distribution and of the gravitational magnification effectively calibrate the HFF clusters as gravitational telescopes. Allowing the computation of search areas in the source plane, these maps are provided to the community to facilitate the exploitation of forthcoming HFF data for quantitative studies of the gravitationally lensed population of background galaxies. Our models of the gravitational magnification afforded by the HFF clusters allow us to quantify the lensing-induced boost in sensitivity over blank-field observations and predict that galaxies at $z>10$ and as faint as m(AB)=32 will be detectable, up to 2 magnitudes fainter than the limit of the Hubble Ultra Deep Field.

We present the first scaling relation between weak-lensing galaxy cluster mass, $M_{WL}$, and near-infrared luminosity, $L_K$. Our results are based on 17 clusters observed with wide-field instruments on Subaru, the United Kingdom Infrared Telescope, the Mayall Telescope, and the MMT. We concentrate on the relation between projected 2D weak-lensing mass and spectroscopically confirmed luminosity within 1Mpc, modelled as $M_{WL} \propto L_{K}^b$, obtaining a power law slope of $b=0.83^{+0.27}_{-0.24}$ and an intrinsic scatter of $\sigma_{lnM_{WL}|L_{K}}=10^{+8}_{-5}\%$. Intrinsic scatter of ~10% is a consistent feature of our results regardless of how we modify our approach to measuring the relationship between mass and light. For example, deprojecting the mass and measuring both quantities within $r_{500}$, that is itself obtained from the lensing analysis, yields $\sigma_{lnM_{WL}|L_{K}}=10^{+7}_{-5}\%$ and $b=0.97^{+0.17}_{-0.17}$. We also find that selecting members based on their (J-K) colours instead of spectroscopic redshifts neither increases the scatter nor modifies the slope. Overall our results indicate that near-infrared luminosity measured on scales comparable with $r_{500}$ (typically 1Mpc for our sample) is a low scatter and relatively inexpensive proxy for weak-lensing mass. Near-infrared luminosity may therefore be a useful mass proxy for cluster cosmology experiments.

We present a high-precision mass model of the galaxy cluster MACSJ0416.1-2403, based on a strong-gravitational-lensing analysis of the recently acquired Hubble Space Telescope Frontier Fields (HFF) imaging data. Taking advantage of the unprecedented depth provided by HST/ACS observations in three passbands, we identify 51 new multiply imaged galaxies, quadrupling the previous census and bringing the grand total to 68, comprising 194 individual lensed images. Having selected a subset of the 57 most securely identified multiply imaged galaxies, we use the Lenstool software package to constrain a lens model comprised of two cluster-scale dark-matter halos and 98 galaxy-scale halos. Our best-fit model predicts image positions with an $RMS$ error of 0.68'', which constitutes an improvement of almost a factor of two over previous, pre-HFF models of this cluster. We find the total projected mass inside a 200~kpc aperture to be $(1.60\pm0.01)\times 10^{14}\ M_\odot$, a measurement that offers a three-fold improvement in precision, reaching the percent level for the first time in any cluster. Finally, we quantify the increase in precision of the derived gravitational magnification of high-redshift galaxies and find an improvement by a factor of $\sim$2.5 in the statistical uncertainty. Our findings impressively confirm that HFF imaging has indeed opened the domain of high-precision mass measurements for massive clusters of galaxies.

We have discovered that Europa, Ganymede and Callisto are bright around 1.5 {\mu}m even when not directly lit by sunlight, based on observations from the Hubble Space Telescope and the Subaru Telescope. The observations were conducted with non-sidereal tracking on Jupiter outside of the field of view to reduce the stray light subtraction uncertainty due to the close proximity of Jupiter. Their eclipsed luminosity was $10^{-6}$-$10^{-7}$ of their uneclipsed brightness, which is low enough that this phenomenon has been undiscovered until now. In addition, Europa in eclipse was <1/10 of the others at 1.5 {\mu}m, a potential clue to the origin of the source of luminosity. Likewise, Ganymede observations were attempted at 3.6 {\mu}m by the Spitzer Space Telescope but it was not detected, suggesting a significant wavelength dependence. The reason why they are luminous even when in the Jovian shadow is still unknown, but forward-scattered sunlight by haze in the Jovian upper atmosphere is proposed as the most plausible candidate. If this is the case, observations of these Galilean satellites while eclipsed by the Jovian shadow provide us a new technique to investigate Jovian atmospheric composition, and investigating the transmission spectrum of Jupiter by this method is important for investigating the atmosphere of extrasolar giant planets by transit spectroscopy.

We present a catalog of 129 X-ray galaxy groups, covering a redshift range 0.04<z<1.23, selected in the ~3 square degree part of the CFHTLS W1 field overlapping XMM observations performed under the XMM-LSS project. We carry out a statistical study of the redshift evolution out to redshift one of the magnitude gap between the first and the second brightest cluster galaxies of a well defined mass-selected group sample. We find that the slope of the relation between the fraction of groups and the magnitude gap steepens with redshift, indicating a larger fraction of fossil groups at lower redshifts. We find that 22.2$\pm$6% of our groups at z$\leq$0.6 are fossil groups. We compare our results with the predictions of three semi-analytic models based on the Millennium simulation. The intercept of the relation between the magnitude of the brightest galaxy and the value of magnitude gap becomes brighter with increasing redshift. This trend is steeper than the model predictions which we attribute to the younger stellar age of the observed brightest cluster galaxies. This trend argues in favor of stronger evolution of the feedback from active galactic nuclei at z<1 compared to the models. The slope of the relation between the magnitude of the brightest cluster galaxy and the value of the gap does not evolve with redshift and is well reproduced by the models, indicating that the tidal galaxy stripping, put forward as an explanation of the occurrence of the magnitude gap, is both a dominant mechanism and is sufficiently well modeled.

We present a comprehensive study of star-forming (SF) galaxies in the HST Frontier Field recent cluster merger A2744 (z=0.308). Wide-field, ultraviolet-infrared (UV-IR) imaging enables a direct constraint of the total star formation rate (SFR) for 53 cluster galaxies, with SFR{UV+IR}=343+/-10 Msun/yr. Within the central 4 arcmin (1.1 Mpc) radius, the integrated SFR is complete, yielding a total SFR{UV+IR}=201+/-9 Msun/yr. Focussing on obscured star formation, this core region exhibits a total SFR{IR}=138+/-8 Msun/yr, a mass-normalised SFR{IR} of Sigma{SFR}=11.2+/-0.7 Msun/yr per 10^14 Msun and a fraction of IR-detected SF galaxies f{SF}=0.080(+0.010,-0.037). Overall, the cluster population at z~0.3 exhibits significant intrinsic scatter in IR properties (total SFR{IR}, Tdust distribution) apparently unrelated to the dynamical state: A2744 is noticeably different to the merging Bullet cluster, but similar to several relaxed clusters. However, in A2744 we identify a trail of SF sources including jellyfish galaxies with substantial unobscured SF due to extreme stripping (SFR{UV}/SFR{IR} up to 3.3). The orientation of the trail, and of material stripped from constituent galaxies, indicates that the passing shock front of the cluster merger was the trigger. Constraints on star formation from both IR and UV are crucial for understanding galaxy evolution within the densest environments.

We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect $5\times 10^{10}~\rm M_\odot$ of molecular gas within 10 kpc of the BCG. Its ensemble velocity profile width of $\sim 130 ~\rm km~s^{-1}$ FWHM is too narrow for the molecular cloud sto be supported in the galaxy by dynamic pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. Roughly $10^{10}~\rm M_\odot$ of molecular gas is projected $3-10 ~\rm kpc$ to the north-west and to the east of the nucleus with line of sight velocities lying between $-250 ~\rm km~s^{-1}$ to $+480 ~\rm km~s^{-1}$ with respect to the systemic velocity. The high velocity gas may be either inflowing or outflowing. However, the absence of high velocity gas toward the nucleus that would be expected in a steady inflow, and its bipolar distribution on either side of the nucleus, are more naturally explained as outflow. Star formation and radiation from the AGN are both incapable of driving an outflow of this magnitude. If so, the molecular outflow may be associated a hot outflow on larger scales reported by Kirkpatrick and colleagues. The molecular gas flow rate of approximately $200~\rm M_\odot ~yr^{-1}$ is comparable to the star formation rate of $100-180~\rm M_\odot ~yr^{-1}$ in the central disk. How radio bubbles would lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.

The Hubble Frontier Fields (HFF) program combines the capabilities of the Hubble Space Telescope (HST) with the gravitational lensing of massive galaxy clusters to probe the distant Universe to an unprecedented depth. Here, we present the results of the first combined HST and Spitzer observations of the cluster Abell 2744. We combine the full near-infrared data with ancillary optical images to search for gravitationally lensed high-redshift (z > 6) galaxies. We report the detection of 15 I814-dropout candidates at z ~ 6-7 and one Y105-dropout at z ~ 8 in a total survey area of 1.43 arcmin^2 in the source plane. The predictions of our lens model allow us to also identify five multiply-imaged systems lying at redshifts between z ~ 6 and z ~ 8. Thanks to constraints from the mass distribution in the cluster, we were able to estimate the effective survey volume corrected for completeness and magnification effects. This was in turn used to estimate the rest-frame ultraviolet luminosity function (LF) at z ~ 6-8. Our LF results are generally in agreement with the most recent blank field estimates, confirming the feasibility of surveys through lensing clusters. Although based on a shallower observations than what will be achieved in the final dataset including the full ACS observations, the LF presented here extends down to Muv ~ -18.5 at z ~ 7 with one identified object at Muv ~ -15 thanks to the highly-magnified survey areas. This early study forecasts the power of using massive galaxy clusters as cosmic telescopes and its complementarity to blank fields.

We present in this letter the first analysis of a z~8 galaxy candidate found in the Hubble and Spitzer imaging data of Abell 2744, as part of the Hubble Frontier Fields legacy program. We applied the most commonly-used methods to select exceptionally high-z galaxies by combining non-detection and color-criteria using seven HST bands. We used GALFIT on IRAC images for fitting and subtracting contamination of bright nearby sources. The physical properties have been inferred from SED-fitting using templates with and without nebular emission. This letter is focussed on the brightest candidate we found (m$_{F160W}$=26.2) over the 4.9 arcmin$^2$ field of view covered by the WFC3. It shows a non-detection in the ACS bands and at 3.6{\mu}m whereas it is clearly detected at 4.5{\mu}m with rather similar depths. This break in the IRAC data could be explained by strong [OIII]+H{\beta} lines at z~8 which contribute to the 4.5{\mu}m photometry. The best photo-z is found at z~8.0$^{+0.2}_{-0.5}$, although solutions at low-redshift (z~1.9) cannot be completely excluded, but they are strongly disfavoured by the SED-fitting work. The amplification factor is relatively small at {\mu}=1.49$\pm$0.02. The Star Formation Rate in this object is ranging from 8 to 60 Mo/yr, the stellar mass is in the order of M$_{\star}$=(2.5-10) x 10$^{9}$Mo and the size is r~0.35$\pm$0.15 kpc. This object is one of the first z~8 LBG candidates showing a clear break between 3.6{\mu}m and 4.5{\mu}m which is consistent with the IRAC properties of the first spectroscopically confirmed galaxy at a similar redshift. Due to its brightness, the redshift of this object could potentially be confirmed by near infrared spectroscopy with current 8-10m telescopes. The nature of this candidate will be revealed in the coming months with the arrival of new ACS and Spitzer data, increasing the depth at optical and near-IR wavelengths.

We present Submillimeter Array (SMA) [CII] 158um and Jansky Very Large Array (JVLA) $^{12}$CO(1-0) line emission maps for the bright, lensed, submillimeter source at $z=5.2430$ behind Abell 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple $^{12}$CO rotational transitions, [CI], water and [NII], providing some of the best constraints on the properties of the interstellar medium (ISM) in a galaxy at $z>5$. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8-1000um) = (1.6$\pm$0.1)x10^14 L_sun/mu, where the total magnification mu_total = 8.9$\pm$1.9, via a new lens model from the [CII] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components which correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source plane separation <1 kpc. The reddest dominates the continuum map (de-magnified L_FIR = (1.1$\pm$0.2)x10^13 L_sun), and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (dV ~ 500 km/s) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region, and photo-dissociation region (PDR) analysis suggests that it is lower density, cooler and forming stars less vigorously than the other components. Strikingly, it has very strong [NII] emission which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity.

We report our analysis of MACS J0717.5+3745 using 140 and 268 GHz Bolocam data collected at the Caltech Submillimeter Observatory. We detect extended Sunyaev-Zel'dovich (SZ) effect signal at high significance in both Bolocam bands, and we employ Herschel-SPIRE observations to subtract the signal from dusty background galaxies in the 268 GHz data. We constrain the two-band SZ surface brightness toward two of the sub-clusters of MACS J0717.5+3745: the main sub-cluster (named C), and a sub-cluster identified in spectroscopic optical data to have a line-of-sight velocity of +3200 km/s (named B). We determine the surface brightness in two separate ways: via fits of parametric models and via direct integration of the images. For both sub-clusters, we find consistent surface brightnesses from both analysis methods. We constrain spectral templates consisting of relativistically corrected thermal and kinetic SZ signals, using a jointly-derived electron temperature from Chandra and XMM-Newton under the assumption that each sub-cluster is isothermal. The data show no evidence for a kinetic SZ signal toward sub-cluster C, but they do indicate a significant kinetic SZ signal toward sub-cluster B. The model-derived surface brightnesses for sub-cluster B yield a best-fit line-of-sight velocity of v_z = +3450 +- 900 km/s, with (1 - Prob[v_z > 0]) = 1.3 x 10^-5 (4.2 sigma away from 0 for a Gaussian distribution). The directly integrated sub-cluster B SZ surface brightnesses provide a best-fit v_z = +2550 +- 1050 km/s, with (1 - Prob[v_z > 0]) = 2.2 x 10^-3 (2.9 sigma).

Multi-wavelength, optical to IR/sub-mm observations of 5 strongly lensed galaxies identified by the Herschel Lensing Survey, plus two well-studied lensed galaxies, MS1512-cB58 and the Cosmic Eye, for which we also provide updated Herschel measurements, are used to determine the physical properties of z~1.5-3 star-forming galaxies close to or below the detection limits of blank fields. We constrain their stellar and dust content, determine star formation rates and histories, dust attenuation and extinction laws, and other related properties. We perform SED-fits of the full photometry of each object as well for the optical and infrared parts separately, exploring various parameters, including nebular emission. The IR observations and emission line measurements, where available, are used a posteriori constraints on the models. Besides the various stellar population models we explore, we use the observed IR/UV ratio to estimate the extinction and create "energy conserving models", that constrain most accurately the physical properties of our sources. Our sample has a median lensing-corrected IR luminosity ~ 3e11 Lsun, stellar masses between 2e9 and 2e11 Msun, and IR/UV luminosity ratios spanning a wide range. The dust masses of our galaxies are in the range 2 to 17e7 Msun, extending previous studies at the same redshift down to lower masses. We do not find any particular trend of the dust temperature Tdust with IR luminosity, suggesting an overall warmer dust regime at our redshift regardless of luminosity. Lensing enables us to study the detailed physical properties of individual IR-detected z~1.5-3 galaxies up to a factor ~10 fainter than achieved with deep blank field observations. We demonstrate that multi-wavelength observations combining stellar and dust emission can constrain star formation histories and extinction laws of star-forming galaxies.

In the course of our 870um APEX/LABOCA follow up of the Herschel Lensing Survey we have detected a source in AS1063 (RXC J2248.7-4431), that has no counterparts in any of the Herschel PACS/SPIRE bands, it is a Herschel 'drop-out' with S_870/S_500>0.5. The 870um emission is extended and centered on the brightest cluster galaxy suggesting either a multiply imaged background source or substructure in the Sunyaev-Zel'dovich (SZ) increment due to inhomogeneities in the hot cluster gas of this merging cluster. We discuss both interpretations with emphasis on the putative lensed source. Based on the observed properties and on our lens model we find that this source could be the first SMG with a moderate far infrared luminosity (L_FIR<10^12 L_sol) detected so far at z>4. In deep HST observations we identified a multiply imaged z~6 source and we measured its spectroscopic redshift z=6.107 with VLT/FORS. This source could be associated with the putative SMG but it is most likely offset spatially by 10-30kpc and they could be interacting galaxies. With a FIR luminosity in the range [5-15]x10^{11} L_sol corresponding to a star formation rate in the range [80-260]M_sol/yr, this SMG would be more representative than the extreme starbursts usually detected at z>4. With a total magnification of ~25 it would open a unique window to the 'normal' dusty galaxies at the end of the epoch of reionization.

We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect 5E10 solar masses of molecular gas within 10 kpc of the BCG. Its velocity width of ~130 km/s FWHM is too narrow to be supported by dynamical pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. The disk is forming stars at a rate of 100-180 solar masses per year. Roughly 1E10 solar masses of molecular gas is projected 3-10 kpc to the north-west and to the east of the nucleus with line of sight velocities lying between -250 km/s to +480 km/s with respect to the systemic velocity. Although inflow cannot be ruled out, the rising velocity gradient with radius is consistent with a broad, bipolar outflow driven by radio jets or buoyantly rising X-ray cavities. The molecular outflow may be associated with an outflow of hot gas in Abell 1835 seen on larger scales. Molecular gas is flowing out of the BCG at a rate of approximately 200 solar masses per year, which is comparable to its star formation rate. How radio bubbles lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.

We present a detailed structural and morphological study of a large sample of spectroscopically-confirmed galaxies at z >= 6, using deep HST near-IR broad-band images and Subaru optical narrow-band images. The galaxy sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV continuum morphology in the HST images, from compact features to multiple component systems. The fraction of merging/interacting galaxies reaches 40% ~ 50% at the brightest end of M_1500 <= -20.5 mag. The intrinsic half-light radii r_{hl,in}, after correction for PSF broadening, are roughly between r_{hl,in} ~ 0.05" (0.3 kpc) and 0.3" (1.7 kpc) at M_1500 <= -19.5 mag. The median r_{hl,in} value is 0.16" (~0.9 kpc). This is consistent with the sizes of bright LAEs and LBGs at z >= 6 in previous studies. In addition, more luminous galaxies tend to have larger sizes, exhibiting a weak size-luminosity relation r_{hl,in} \propto L^{0.14} at M_1500 <= -19.5 mag. The slope of 0.14 is significantly flatter than those in fainter LBG samples. We discuss the morphology of z >= 6 galaxies with nonparametric methods, including the CAS system and the Gini and M_20 parameters, and demonstrate their validity through simulations. We search for extended Lyman-alpha emission halos around LAEs at z ~ 5.7 and 6.5, by stacking a number of narrow-band images. We do not find evidence of extended halos predicted by cosmological simulations. Such Lyman-alpha halos, if they exist, could be weaker than predicted. Finally, we investigate any positional misalignment between UV continuum and Lyman-alpha emission in LAEs. While the two positions are generally consistent, several merging galaxies show significant positional differences. This is likely caused by a disturbed ISM distribution due to merging activity.

We present an analysis of the levels and evolution of star formation activity in a representative sample of 30 massive galaxy clusters at 0.15<z<0.30 from the Local Cluster Substructure Survey (LoCuSS), combining wide-field Spitzer 24um data with extensive spectroscopy of cluster members. The specific-SFRs of massive (M>10^10 M_sun) star-forming cluster galaxies within r200 are found to be systematically 28% lower than their counterparts in the field at fixed stellar mass and redshift, a difference significant at the 8.7-sigma level. This is the unambiguous signature of star formation in most (and possibly all) massive star-forming galaxies being slowly quenched upon accretion into massive clusters, their SFRs declining exponentially on quenching time-scales in the range 0.7-2.0 Gyr. We measure the mid-infrared Butcher-Oemler effect over the redshift range 0.0-0.4, finding rapid evolution in the fraction (f_SF) of massive (M_K<-23.1) cluster galaxies within r200 with SFRs>3M_sun/yr, of the form f_SF (1+z)^7.6. We dissect the origins of the Butcher-Oemler effect, revealing it to be due to the combination of a ~3x decline in the mean specific-SFRs of star-forming cluster galaxies since z~0.3 with a ~1.5x decrease in number density. Two-thirds of this reduction in the specific-SFRs of star-forming cluster galaxies is due to the steady cosmic decline in the specific-SFRs among those field galaxies accreted into the clusters. The remaining one-third reflects an accelerated decline in the star formation activity of galaxies within clusters. The slow quenching of star-formation in cluster galaxies is consistent with a gradual shut down of star formation in infalling spiral galaxies as they interact with the intra-cluster medium via ram-pressure stripping or starvation mechanisms. We find no evidence for the build-up of cluster S0 bulges via major nuclear star-burst episodes.

We present deep HST near-IR and Spitzer mid-IR observations of a large sample of spectroscopically-confirmed galaxies at z >= 6. The sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in the 3.6um and 4.5um bands. Our galaxies also have deep optical imaging data from Subaru Suprime-Cam. We utilize the multi-band data and secure redshifts to derive their rest-frame UV properties. These galaxies have steep UV continuum slopes roughly between beta ~ -1.5 and -3.5, with an average value of beta ~ -2.3, slightly steeper than the slopes of LBGs in previous studies. The slope shows little dependence on UV continuum luminosity except for a few of the brightest galaxies. We find a statistically significant excess of galaxies with slopes around beta ~ -3, suggesting the existence of very young stellar populations with extremely low metallicity and dust content. Our galaxies have moderately strong rest-frame Lyman-alpha equivalent width (EW) in a range of ~10 to ~200 \AA. The star-formation rates are also moderate, from a few to a few tens solar masses per year. The LAEs and LBGs in this sample share many common properties, implying that LAEs represent a subset of LBGs with strong Lyman-alpha emission. Finally, the comparison of the UV luminosity functions between LAEs and LBGs suggests that there exists a substantial population of faint galaxies with weak Lyman-alpha emission (EW < 20 \AA) that could be the dominant contribution to the total ionizing flux at z >= 6.

We present deep spectroscopic observations of a Lyman-alpha emitter (LAE) candidate at z ~ 7.7 using the infrared spectrograph LUCI on the 2 x 8.4m Large Binocular Telescope (LBT). The candidate is the brightest among the four z ~ 7.7 LAE candidates found in a narrow-band imaging survey by Krug et al. 2012. Our spectroscopic data include a total of 7.5 hours of integration with LBT/LUCI and are deep enough to significantly (3.2-4.9 sigma) detect the Lyman-alpha emission line of this candidate, based on its Lyman-alpha flux 1.2 x 10^{-17} erg s^{-1} cm^{-2} estimated from the narrow-band photometry. However, we do not find any convincing signal at the expected position of its Lyman-alpha emission line, suggesting that this source is not an LAE at z ~ 7.7. The non-detection in this work, together with the previous studies of z ~ 7.7 LAEs, puts a strong constraint on the bright-end Lyman-alpha luminosity function (LF) at z ~ 7.7. We find a rapid evolution of the Lyman-alpha LF from z ~ 6.5 to 7.7: the upper limit of the z ~ 7.7 LF is more than 5 times lower than the z ~ 6.5 LF at the bright end (f > 1.0 x 10^{-17} erg s^{-1} cm^{-2}, or L > 6.9 x 10^{42} erg s^{-1}). This is likely caused by an increasing neutral fraction in the IGM that substantially attenuates Lyman-alpha emission at z ~ 7.7.

We present results from a Spitzer, mid-infrared spectroscopy study of a sample of 74 galaxies located in 23 Hickson Compact Groups, chosen to be at a dynamically-active stage of HI depletion. We find evidence for enhanced warm H2 emission (i.e. above that associated with UV excitation in star-forming regions) in 14 galaxies (~20%), with 8 galaxies having extreme values of L(H2 S(0)-S(3))/L(7.7micron PAH), in excess of 0.07. Such emission has been seen previously in the compact group HCG 92 (Stephan's Quintet), and was shown to be associated with the dissipation of mechanical energy associated with a large-scale shock caused when one group member collided, at high velocity, with tidal debris in the intragroup medium. Similarly, shock excitation or turbulent heating is likely responsible for the enhanced H2 emission in the compact group galaxies, since other sources of heating (UV or X-ray excitation from star formation or AGN) are insufficient to account for the observed emission. The group galaxies fall predominantly in a region of mid-infrared color-color space identified by previous studies as being connected to rapid transformations in HCG galaxy evolution. Furthermore, the majority of H2-enhanced galaxies lie in the optical "green valley" between the blue cloud and red-sequence, and are primarily early-type disk systems. We suggest that H2-enhanced systems may represent a specific phase in the evolution of galaxies in dense environments and provide new insight into mechanisms which transform galaxies onto the optical red sequence.

We have identified a very interesting Ly-alpha emitter, whose Ly-alpha emission line has an extremely large observed equivalent width of EW_0=436^{+422}_{-149}A, which corresponds to an extraordinarily large intrinsic rest-frame equivalent width of EW_0^{int}=872^{+844}_{-298}A after the average intergalactic absorption correction. The object was spectroscopically confirmed to be a real Ly-alpha emitter by its apparent asymmetric Ly-alpha line profile detected at z=6.538. The continuum emission of the object was definitely detected in our deep z'-band image; thus, its EW_0 was reliably determined. Follow-up deep near-infrared spectroscopy revealed emission lines of neither He II lambda1640 as an apparent signature of Population III, nor C IV lambda1549 as a proof of active nucleus. No detection of short-lived He II lambda1640 line is not necessarily inconsistent with the interpretation that the underlying stellar population of the object is dominated by Population III. We found that the observed extremely large EW_0 of the Ly-alpha emission and the upper limit on the EW_0 of the He II lambda1640 emission can be explained by population synthesis models favoring a very young age less than 2-4Myr and massive metal-poor (Z<10^{-5}) or even metal-free stars. The observed large EW_0 of Ly-alpha is hardly explained by Population I/II synthesis models with Z>10^{-3}. However, we cannot conclusively rule out the possibility that this object is composed of a normal stellar population with a clumpy dust distribution, which could enhance the Ly-alpha EW_0, though its significance is still unclear.

In the cores of some galaxy clusters the hot intracluster plasma is dense enough that it should cool radiatively in the cluster's lifetime, leading to continuous "cooling flows" of gas sinking towards the cluster center, yet no such cooling flow has been observed. The low observed star formation rates and cool gas masses for these "cool core" clusters suggest that much of the cooling must be offset by astrophysical feedback to prevent the formation of a runaway cooling flow. Here we report X-ray, optical, and infrared observations of the galaxy cluster SPT-CLJ2344-4243 at z = 0.596. These observations reveal an exceptionally luminous (L_2-10 keV = 8.2 x 10^45 erg/s) galaxy cluster which hosts an extremely strong cooling flow (dM/dt = 3820 +/- 530 Msun/yr). Further, the central galaxy in this cluster appears to be experiencing a massive starburst (740 +/- 160 Msun/yr), which suggests that the feedback source responsible for preventing runaway cooling in nearby cool core clusters may not yet be fully established in SPT-CLJ2344-4243. This large star formation rate implies that a significant fraction of the stars in the central galaxy of this cluster may form via accretion of the intracluster medium, rather than the current picture of central galaxies assembling entirely via mergers.

The AzTEC 1.1 mm survey of the SCUBA HAlf Degree Extragalactic Survey (SHADES) fields is the largest (0.7 deg2) blank-field millimetre-wavelength survey undertaken to date at a resolution of ~18" and a depth of ~1 mJy. We have used the deep optical-to-radio multi-wavelength data in the SHADES Lockman Hole East and SXDF/UDS fields to obtain galaxy identifications for ~64% (~80% with tentative identifications) of the 148 AzTEC-SHADES 1.1 mm sources, exploiting deep radio and 24 um data complemented by methods based on 8 um flux-density and red optical-infrared (i-K) colour. This unusually high identification rate can be attributed to the relatively bright millimetre-wavelength flux-density threshold, combined with the relatively deep supporting multi-frequency data. We have further exploited the optical-mid-infrared-radio data to derive a ~60% (~75% with tentative identifications) complete redshift distribution for the AzTEC-SHADES sources, yielding a median redshift of z~2.2, with a high-redshift tail extending to at least z~4. Despite the larger area probed by the AzTEC survey relative to the original SCUBA SHADES imaging, the redshift distribution of the AzTEC sources is consistent with that displayed by the SCUBA sources, and reinforces tentative evidence that the redshift distribution of mm/sub-mm sources in the Lockman Hole field is significantly different from that found in the SXDF/UDS field. Comparison with simulated surveys of similar scale extracted from semi-analytic models based on the Millennium simulation indicates that this is as expected if the mm/sub-mm sources are massive (M>10^11 Mo) star-forming galaxies tracing large-scale structures over scales of 10-20 Mpc. This confirms the importance of surveys covering several square degrees (as now underway with SCUBA2) to obtain representative samples of bright (sub)mm-selected galaxies.

Using far-infrared imaging from the "Herschel Lensing Survey", we derive dust properties of spectroscopically-confirmed cluster member galaxies within two massive systems at z~0.3: the merging Bullet Cluster and the more relaxed MS2137.3-2353. Most star-forming cluster sources (~90%) have characteristic dust temperatures similar to local field galaxies of comparable infrared (IR) luminosity (T_dust ~ 30K). Several sub-LIRG (L_IR < 10^11 L_sun) Bullet Cluster members are much warmer (T_dust > 37K) with far-infrared spectral energy distribution (SED) shapes resembling LIRG-type local templates. X-ray and mid-infrared data suggest that obscured active galactic nuclei do not contribute significantly to the infrared flux of these "warm dust" galaxies. Sources of comparable IR-luminosity and dust temperature are not observed in the relaxed cluster MS2137, although the significance is too low to speculate on an origin involving recent cluster merging. "Warm dust" galaxies are, however, statistically rarer in field samples (> 3sigma), indicating that the responsible mechanism may relate to the dense environment. The spatial distribution of these sources is similar to the whole far-infrared bright population, i.e. preferentially located in the cluster periphery, although the galaxy hosts tend towards lower stellar masses (M_* < 10^10 M_sun). We propose dust stripping and heating processes which could be responsible for the unusually warm characteristic dust temperatures. A normal star-forming galaxy would need 30-50% of its dust removed (preferentially stripped from the outer reaches, where dust is typically cooler) to recover a SED similar to a "warm dust" galaxy. These progenitors would not require a higher IR-luminosity or dust mass than the currently observed normal star-forming population.

We present a study of the distribution of X-ray AGN in a representative sample of 26 massive clusters at 0.15<z<0.30, combining Chandra observations with highly complete spectroscopy of cluster members down to M_K*+2. In total we identify 48 X-ray AGN among the cluster members, with luminosities 2x10^41-1x10^44erg/s. In the stacked caustic diagram, the X-ray AGN appear to preferentially lie along the caustics, suggestive of an infalling population. They also appear to avoid the region with lowest cluster-centric radii and relative velocities (r_proj<0.4 r_500; |v-<v>|/sigma_v<0.8), which is dominated by the virialized population of galaxies accreted earliest into the clusters. Moreover the velocity dispersion of the 48 X-ray AGN is 1.51x that of the overall cluster population, which is consistent with the sqrt(2) ratio expected by simple energetic arguments when comparing infalling versus virialized populations. This kinematic segregation is significant at the 4.66-sigma level. When splitting the X-ray AGN sample into two according to X-ray or infrared (IR) luminosity, both X-ray bright and IR-bright sub-samples show higher velocity dispersions than their X-ray dim and IR-dim counterparts at >2sigma significance. This is consistent with the nuclear activity responsible for the X-ray and IR emission being slowly shut down as the host galaxies are accreted into the cluster. Overall our results provide the strongest observational evidence to date that X-ray AGN found in massive clusters are an infalling population, and that the cluster environment very effectively suppresses radiatively-efficient nuclear activity in its member galaxies. These results are consistent with the view that for galaxies to host an X-ray AGN they should be the central galaxy within their dark matter halo and have a ready supply of cold gas.

We searched for z=7.3 Lya emitters (LAEs) behind two lensing clusters, Abell 2390 and CL 0024, with the Subaru Telescope Suprime-Cam and a narrowband NB1006 (FWHM ~ 21 nm centered at 1005 nm). We investigated if there exist objects consistent with the color of z=7.3 LAEs behind the clusters but could not detect any LAEs to the unlensed line limit F(Lya) ~ 6.9 x 10^{-18} erg/s/cm^2. Using several z=7 Lya luminosity functions (LFs) from the literature, we estimated and compared the expected detection numbers of z ~ 7 LAEs in lensing and blank field surveys in the case of using an 8m class ground based telescope. Given the steep bright-end slope of the LFs, when the detector field-of view (FOV) is comparable to the angular extent of a massive lensing cluster, imaging cluster(s) is more efficient in detecting z ~ 7 LAEs than imaging a blank field. However, the gain is expected to be modest, a factor of two at most and likely much less depending on the adopted LFs. The main advantage of lensing-cluster survey, therefore, remains to be the gain in depth and not necessarily in detection efficiency. For much larger detectors, the lensing effect becomes negligible and the efficiency of LAE detection is proportional to the instrumental FOV. We also inspected NB1006 images of three z ~ 7 z-dropouts previously detected in Abell 2390 and found that none of them are detected in NB1006. Two of them are consistent with predictions from the previous studies that they would be at lower redshifts. The other one has a photometric redshift of z ~ 7.3, and if it is at z=7.3, its unlensed Lya line flux would be very faint: F(Lya) < 4.4 x 10^{-18} erg/s/cm^2 (1 sigma upper limit) or rest frame equivalent width of W(Lya) < 26A. Its Lya emission might be attenuated by neutral hydrogen, as recent studies show that the fraction of Lyman break galaxies displaying strong Lya emission is lower at z ~ 7 than at z <~ 6.

We present a search for the [CII] 158micron fine structure line (a main cooling line of the interstellar medium) and the underlying far-infrared (FIR) continuum in three high-redshift (6.6<z<8.2) star-forming galaxies using the IRAM Plateau de Bure interferometer. We targeted two Lyman-Alpha-selected galaxies (Lyman-Alpha-Emitters, LAEs) with moderate UV-based star formation rates (SFR~20 M_sun/yr; Himiko at z=6.6 and IOK-1 at z=7.0) and a Gamma Ray Burst (GRB) host galaxy (GRB 090423 at z~8.2). Based on our 3 sigma rest-frame FIR continuum limits, previous (rest-frame) UV continuum measurements and spectral energy distribution (SED) fitting, we rule out SED shapes similar to highly obscured galaxies (e.g. Arp220, M82) and less extreme dust-rich nearby spiral galaxies (e.g. M51) for the LAEs. Conservatively assuming a SED shape typical of local spiral galaxies we derive upper limits for the FIR-based star formation rates (SFRs) of ~70 M_sun/yr, ~50 M_sun/yr and ~40 M_sun/yr for Himiko, IOK-1 and GRB 090423, respectively. For the LAEs these limits are only a factor ~3 higher than the published UV-based SFRs (uncorrected for extinction). This indicates that the dust obscuration in the z>6 LAEs studied here is lower by a factor of a few than what has recently been found in some LAEs at lower redshift (2<z<3.5) with similar UV-based SFRs. A low obscuration in our z>6 LAE sample is consistent with recent rest-frame UV studies of z~7 Lyman-Break-Galaxies (LBGs).

We report the discovery of a protocluster at z~6 containing at least eight cluster member galaxies with spectroscopic confirmations in the wide-field image of the Subaru Deep Field (SDF). The overdensity of the protocluster is significant at the 6 sigma level, based on the surface number density of i'-dropout galaxies. The overdense region covers ~36 sq. arcmin, and includes 30 i'-dropout galaxies. Follow-up spectroscopy revealed that 15 of these are real z~6 galaxies (5.7 < z < 6.3). Eight of the 15 are clustering in a narrow redshift range centered at z=6.01, corresponding to a seven-fold increase in number density over the average in redshift space. We found no significant difference in the observed properties, such as Ly-alpha luminosities and UV continuum magnitudes, between the eight protocluster members and the seven non-members. The velocity dispersion of the eight protocluster members is 647 km/s, which is about three times higher than that predicted by the standard cold dark matter model. This discrepancy could be attributed to the distinguishing three-dimensional distribution of the eight protocluster members. We discuss two possible explanations for this discrepancy: either the protocluster is already mature, with old galaxies at the center, or it is still immature and composed of three subgroups merging to become a larger cluster. In either case, this concentration of z=6.01 galaxies in the SDF may be one of the first sites of formation of a galaxy cluster in the universe.

We have discovered a 2.5 Mpc (projected) long filament of infrared-bright galaxies connecting two of the three ~5x10^14 Msun clusters making up the RCS 2319+00 supercluster at z=0.9. The filament is revealed in a deep Herschel Spectral and Photometric Imaging REceiver (SPIRE) map that shows 250-500um emission associated with a spectroscopically identified filament of galaxies spanning two X-ray bright cluster cores. We estimate that the total (8-1000um) infrared luminosity of the filament is Lir~5x10^12 Lsun, which, if due to star formation alone, corresponds to a total SFR 900 Msun/yr. We are witnessing the scene of the build-up of a >10^15 Msun cluster of galaxies, seen prior to the merging of three massive components, each of which already contains a population of red, passive galaxies that formed at z>2. The infrared filament demonstrates that significant stellar mass assembly is taking place in the moderate density, dynamically active circumcluster environments of the most massive clusters at high-redshift, and this activity is concomitant with the hierarchical build-up of large scale structure.

During our Herschel Lensing Survey (HLS) of massive galaxy clusters, we have discovered an exceptionally bright source behind the z=0.22 cluster Abell 773, which appears to be a strongly lensed submillimeter galaxy (SMG) at z=5.2429. This source is unusual compared to most other lensed sources discovered by Herschel so far, because of its higher submm flux (\sim 200mJy at 500\micron) and its high redshift. The dominant lens is a foreground z=0.63 galaxy, not the cluster itself. The source has a far-infrared (FIR) luminosity of L_FIR= 1.1 10^{14}/\mu Lo, where \mu is the magnification factor, likely \sim 11. We report here the redshift identification through CO lines with the IRAM-30m, and the analysis of the gas excitation, based on CO(7-6), CO(6-5), CO(5-4) detected at IRAM and the CO(2-1) at the EVLA. All lines decompose into a wide and strong red component, and a narrower and weaker blue component, 540\kms apart. Assuming the ultraluminous galaxy (ULIRG) CO-to-H2 conversion ratio, the H2 mass is 5.8 10^{11}/\mu Mo, of which one third is in a cool component. From the CI line we derive a CI/H2 number abundance of 6 10^{-5} similar to that in other ULIRGs. The H2O line is strong only in the red velocity component, with an intensity ratio I(H_2O)/I(CO) \sim 0.5, suggesting a strong local FIR radiation field, possibly from an active nucleus (AGN) component. We detect the [NII]205\mics line for the first time at high-z. It shows comparable blue and red components, with a strikingly broad blue one, suggesting strong ionized gas flows.

We present far-infrared (FIR) analysis of 68 Brightest Cluster Galaxies (BCGs) at 0.08 < z < 1.0. Deriving total infrared luminosities directly from Spitzer and Herschel photometry spanning the peak of the dust component (24-500um), we calculate the obscured star formation rate (SFR). 22(+6.2,-5.3)% of the BCGs are detected in the far-infrared, with SFR= 1-150 M_sun/yr. The infrared luminosity is highly correlated with cluster X-ray gas cooling times for cool-core clusters (gas cooling time <1 Gyr), strongly suggesting that the star formation in these BCGs is influenced by the cluster-scale cooling process. The occurrence of the molecular gas tracing Ha emission is also correlated with obscured star formation. For all but the most luminous BCGs (L_TIR > 2x10^11 L_sun), only a small (<0.4 mag) reddening correction is required for SFR(Ha) to agree with SFR_FIR. The relatively low Ha extinction (dust obscuration), compared to values reported for the general star-forming population, lends further weight to an alternate (external) origin for the cold gas. Finally, we use a stacking analysis of non-cool-core clusters to show that the majority of the fuel for star formation in the FIR-bright BCGs is unlikely to originate form normal stellar mass loss.

We present the analysis of Herschel SPIRE far-infrared (FIR) observations of the z = 2.515 lensed galaxy SMM J163554.2+661225. Combining new 250, 350, and 500 micron observations with existing data, we make an improved fit to the FIR spectral energy distribution (SED) of this galaxy. We find a total infrared (IR) luminosity of L(8--1000 micron) = 6.9 +/- 0.6x10^11 Lsol; a factor of 3 more precise over previous L_IR estimates for this galaxy, and one of the most accurate measurements for any galaxy at these redshifts. This FIR luminosity implies an unlensed star formation rate (SFR) for this galaxy of 119 +/- 10 Msol per yr, which is a factor of 1.9 +/- 0.35 lower than the SFR derived from the nebular Pa-alpha emission line (a 2.5-sigma discrepancy). Both SFR indicators assume identical Salpeter initial mass functions (IMF) with slope Gamma=2.35 over a mass range of 0.1 - 100 Msol, thus this discrepancy suggests that more ionizing photons may be necessary to account for the higher Pa-alpha-derived SFR. We examine a number of scenarios and find that the observations can be explained with a varying star formation history (SFH) due to an increasing star formation rate (SFR), paired with a slight flattening of the IMF. If the SFR is constant in time, then larger changes need to be made to the IMF by either increasing the upper-mass cutoff to ~ 200 Msol, or a flattening of the IMF slope to 1.9 +/- 0.15, or a combination of the two. These scenarios result in up to double the number of stars with masses above 20 Msol, which produce the requisite increase in ionizing photons over a Salpeter IMF with a constant SFH.

We present broad band photometry and photometric redshifts for 187611 sources located in ~0.5deg^2 in the Lockman Hole area. The catalog includes 389 X-ray detected sources identified with the very deep XMM-Newton observations available for an area of 0.2 deg^2. The source detection was performed on the Rc, z' and B band images and the available photometry is spanning from the far ultraviolet to the mid infrared, reaching in the best case scenario 21 bands. Astrometry corrections and photometric cross-calibrations over the entire dataset allowed the computation of accurate photometric redshifts. Special treatment is undertaken for the X-ray sources, the majority of which is active galactic nuclei. Comparing the photometric redshifts to the available spectroscopic redshifts we achieve for normal galaxies an accuracy of \sigma_{\Delta z/(1+z)}=0.036, with 12.7% outliers, while for the X-ray detected sources the accuracy is \sigma_{\Delta z/(1+z)}=0.069, with 18.3% outliers, where the outliers are defined as sources with |z_{phot}-z_{spec}|>0.15 (1+z_{spec})}. These results are a significant improvement over the previously available photometric redshifts for normal galaxies in the Lockman Hole, while it is the first time that photometric redshifts are computed and made public for AGN for this field.

We present Keck spectroscopic observations of z>6 Lyman-break galaxy (LBG) candidates in the Subaru Deep Field (SDF). The candidates were selected as i'-dropout objects down to z'=27 AB magnitudes from an ultra-deep SDF z'-band image. With the Keck spectroscopy we identified 19 LBGs with prominent Ly_alpha emission lines at 6< z < 6.4. The median value of the Ly_alpha rest-frame equivalent widths (EWs) is ~50 A, with four EWs >100 A. This well-defined spectroscopic sample spans a UV-continuum luminosity range of -21.8< M_{UV}<-19.5 (0.6~5 L*_{UV}) and a Ly_alpha luminosity range of (0.3~3) x 10^{43} erg s^{-1} (0.3~3 L*_ {Ly_alpha}). We derive the UV and Ly_alpha luminosity functions (LFs) from our sample at <z>~6.2 after we correct for sample incompleteness. We find that our measurement of the UV LF is consistent with the results of previous studies based on photometric LBG samples at 5<z<7. Our Ly_alpha LF is also generally in agreement with the results of Ly_alpha-emitter surveys at z~5.7 and 6.6. This study shows that deep spectroscopic observations of LBGs can provide unique constraints on both the UV and Ly_alpha LFs at z>6.

The spectral energy distributions (SED) of dusty galaxies at intermediate redshift may look similar to very high redshift galaxies in the optical/near infrared (NIR) domain. This can lead to the contamination of high redshift galaxy searches based on broad band optical/NIR photometry by lower redshift dusty galaxies as both kind of galaxies cannot be distinguished. The contamination rate could be as high as 50%. {This work shows how the far infrared (FIR) domain can help to recognize likely low-z interlopers in an optical/NIR search for high-z galaxies.} We analyse the FIR SEDs of two galaxies proposed as very high redshift ($z>7$) dropout candidates based on deep Hawk-I/VLT observations. The FIR SEDs are sampled with PACS/Herschel at 100 and 160\,$\mu$m, with SPIRE/Herschel at 250, 350 and 500\,$\mu$m and with LABOCA/APEX at 870\,$\mu$m. We find that redshifts $>7$ would imply extreme FIR SEDs (with dust temperatures $>100$\,K and FIR luminosities $>10^{13}$\,$L_{\odot}$). At z$\sim$2, instead, the SEDs of both sources would be compatible with that of typical ULIRGs/SMGs. Considering all the data available for these sources from visible to FIR we re-estimate the redshifts and we find $z\sim$1.6--2.5. Due to the strong spectral breaks observed in these galaxies, standard templates from the literature fail to reproduce the visible-near IR part of the SEDs even when additional extinction is included. These sources resemble strongly dust obscured galaxies selected in Spitzer observations with extreme visible-to-FIR colors, and the galaxy GN10 at $z=4$. Galaxies with similar SEDs could contaminate other high redshift surveys.

The He II \lambda 1640 emission line has been suggested as a direct probe of Population III (Pop III) stars at high-redshift, since it can arise from highly energetic ionizing photons associated with hot, metal free stars. We use the HST WFC3/F130N IR narrowband filter to probe He II \lambda 1640 emission in galaxy IOK-1 at z=6.96. The sensitivity of this measurement is >5x deeper than for previous measurements. From this deep narrowband imaging, combined with broadband observations in the F125W and F160W filters, we find the He II flux to be 1.2+/- 1.0x 10^-18 ergs/s/cm^2, corresponding to a 1\sigma upper limit on the Pop III star formation rate (SFR) of ~ 0.5 M_sun/yr for the case of a Salpeter IMF with 50-500M_sun and mass loss. Given that the broadband measurements can be fit with a UV continuum spectral flux density of ~ 4.85x 10^-10x \lambda ^-2.46 ergs/s/cm^2/A, which corresponds to an overall SFR of ~16+/-2.6 M_sun/yr, massive Pop III stars represent < 6% of the total star formation. This measurement places the strongest limit yet on metal-free star formation at high redshift, although the exact conversion from He II luminosity to Pop III SFR is highly uncertain due to the unknown IMF, stellar evolution, and photoionization effects. Although we have not detected He II \lambda 1640 at more than the 1.2\sigma level, our work suggests that a > 3\sigma level detection is possible with JWST.

We carried out extended spectroscopic confirmations of Ly-alpha emitters (LAEs) at z=6.5 and 5.7 in the Subaru Deep Field. Now, the total number of spectroscopically confirmed LAEs is 45 and 54 at z=6.5 and 5.7, respectively, and at least 81% (70%) of our photometric candidates at z=6.5 (5.7) have been spectroscopically identified as real LAEs. We made careful measurements of the Ly-alpha luminosity, both photometrically and spectroscopically, to accurately determine the Ly-alpha and rest-UV luminosity functions (LFs). The substantially improved evaluation of the Ly-alpha LF at z=6.5 shows an apparent deficit from z=5.7 at least at the bright end, and a possible decline even at the faint end, though small uncertainties remain. The rest-UV LFs at z=6.5 and 5.7 are in good agreement, at least at the bright end, in clear contrast to the differences seen in the Ly-alpha LF. These results imply an increase in the neutral fraction of the intergalactic medium from z=5.7 to 6.5. The rest-frame equivalent width (EW_0) distribution at z=6.5 seems to be systematically smaller than z=5.7, and it shows an extended tail toward larger EW_0. The bright end of the rest-UV LF can be reproduced from the observed Ly-alpha LF and a reasonable EW_0-UV luminosity relation. Integrating this rest-UV LF provides the first measurement of the contribution of LAEs to the photon budget required for reionization. The derived UV LF suggests that the fractional contribution of LAEs to the photon budget among Lyman break galaxies significantly increases towards faint magnitudes. Low-luminosity LAEs could dominate the ionizing photon budget, though this inference depends strongly on the uncertain faint-end slope of the Ly-alpha LF.

We investigate the nature and the physical properties of z, Y and J-dropout galaxies selected behind the lensing cluster A2667. This field is part of our project aimed at identifying z~7-10 candidates accessible to spectroscopic studies, based on deep photometry with ESO/VLT HAWK-I and FORS2 (zYJH and Ks-band images, AB(3 sigma)~26-27) on a sample of lensing clusters extracted from our multi-wavelength combined surveys with SPITZER, HST, and Herschel. In this paper we focus on the complete Y and J-dropout sample, as well as the bright z-dropouts fulfilling the selection criteria by Capak et al. (2011). 10 candidates are selected within the common field of ~33 arcmin2 (effective area once corrected for contamination and lensing dilution). All of them are detected in H and Ks bands in addition to J and/or IRAC 3.6/4.5, with H(AB)~23.4 to 25.2, and have modest magnification factors. Although best-fit photometric redshifts place all these candidates at high-z, the contamination by low-z interlopers is estimated at 50-75% level based on previous studies, and the comparison with the blank-field WIRCAM Ultra-Deep Survey (WUDS). The same result is obtained when photometric redshifts include a luminosity prior, allowing us to remove half of the original sample as likely z~1.7-3 interlopers with young stellar pulations and strong extinction. Two additional sources among the remaining sample could be identified at low-z based on a detection at 24 microns and on the HST z_850 band. These low-z interlopers are not well described by current templates given the large break, and cannot be easily identified based solely on optical and near-IR photometry. Given the estimated dust extinction and high SFRs, some of them could be also detected in the IR or sub-mm bands. After correction for likely contaminants, the observed counts at z>7.5 seem to be in agreement with an evolving LF. (abridged)

We report the discovery of a unique $z=6.027$ galaxy, multiply imaged by the cluster Abell 383 and detected in new Hubble Space Telescope ACS and WFC3 imaging, as well as in Warm Spitzer observations. This galaxy was selected as a pair of i-dropouts; its suspected high redshift was confirmed by the measurement of a strong Lyman-alpha line in both images using Keck/DEIMOS. Combining Hubble and Spitzer photometry after correcting for contamination by line emission (estimated to be a small effect), we identify a strong Balmer break of 1.5 magnitudes. Taking into account the magnification factor of 11.4+/-1.9 (2.65+/-0.17 mag) for the brightest image, the unlensed AB magnitude for the source is 27.2+/-0.05 in the H band, corresponding to a 0.4 L* galaxy, and 25.7+/-0.08 at 3.6 um. The UV slope is consistent with beta~2.0, and from the rest-frame UV continuum we measure a current star formation rate of 2.4+/-1.1 Msol/yr. The unlensed half-light radius is measured to be 300 pc, from which we deduce a star-forming surface density of ~10 Msol/yr/kpc2. The Lyman-alpha emission is found to be extended over ~3" along the slit, corresponding to ~5 kpc in the source plane. This can be explained by the presence of a much larger envelope of neutral hydrogen around the star-forming region. Finally, fitting the spectral energy distribution using 7 photometric data points with simple SED models, we derive the following properties: very little reddening, an inferred stellar mass of M*=6e9 Msol, and an inferred age of ~800 Myrs (corresponding to a redshift of formation of ~18). The star-formation rate of this object was likely much stronger in the past than at the time of observation, suggesting that we may be missing a fraction of galaxies at z~6 which have already faded in rest-frame UV wavelengths.

We present the observed-frame optical, near- and mid-infrared properties of X-ray selected AGN in the Lockman Hole. Using a likelihood ratio method on optical, near-infrared or mid-infrared catalogues, we assigned counterparts to 401 out of the 409 X-ray sources of the XMM-Newton catalogue. Accurate photometry was collected for all the sources from U to 24um. We used X-ray and optical criteria to remove any normal galaxies, galactic stars, or X-ray clusters among them and studied the multi-wavelength properties of the remaining 377 AGN. We used a mid-IR colour-colour selection to understand the AGN contribution to the optical and infrared emission. Using this selection, we identified different behaviours of AGN-dominated and host-dominated sources in X-ray-optical-infrared colour-colour diagrams. More specifically, the AGN dominated sources show a clear trend in the f_x/f_R vs. R-K and f_24um/f_R vs. R-K diagrams, while the hosts follow the behaviour of non X-ray detected galaxies. In the optical-near-infrared colour-magnitude diagram we see that the known trend of redder objects being more obscured in X-rays is stronger for AGN-dominated than for host-dominated systems. This is an indication that the trend is more related to the AGN contaminating the overall colours than any evolutionary effects. Finally, we find that a significant fraction (~30%) of the reddest AGN are not obscured in X-rays.

We have discovered an X-ray selected galaxy cluster with a spectroscopic redshift of 1.753. The redshift is of the brightest cluster galaxy (BCG), which is coincident with the peak of the X-ray surface brightness. We also have concordant photometric redshifts for seven additional candidate cluster members. The X-ray luminosity of the cluster is 3.68 +/- 0.70 x 10^43 erg s^-1 in the 0.1 - 2.4 keV band. The optical/IR properties of the BCG imply its formation redshift was ~5 if its stars formed in a short burst. This result continues the trend from lower redshift in which the observed properties of BCGs are most simply explained by a monolithic collapse at very high redshift instead of the theoretically preferred gradual hierarchical assembly at later times. However the models corresponding to different formation redshifts are more clearly separated as our observation epoch approaches the galaxy formation epoch. Although our infrared photometry is not deep enough to define a red sequence, we do identify a few galaxies at the cluster redshift that have the expected red sequence photometric properties.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We study the luminosity gap, dm12, between the first and second ranked galaxies in a sample of 59 massive galaxy clusters, using data from the Hale Telescope, HST, Chandra, and Spitzer. We find that the dm12 distribution, p(dm12), is a declining function of dm12, to which we fitted a straight line: p(dm12) propto -(0.13+/-0.02)dm12. The fraction of clusters with "large" luminosity gaps is p(dm12>=1)=0.37+/-0.08, which represents a 3sigma excess over that obtained from Monte Carlo simulations of a Schechter function that matches the mean cluster galaxy luminosity function. We also identify four clusters with "extreme" luminosity gaps, dm12>=2, giving a fraction of p(dm12>=2)=0.07+0.05-0.03. More generally, large luminosity gap clusters are relatively homogeneous, with elliptical/disky brightest cluster galaxies (BCGs), cuspy gas density profiles (i.e. strong cool cores), high concentrations, and low substructure fractions. In contrast, small luminosity gap clusters are heterogeneous, spanning the full range of boxy/elliptical/disky BCG morphologies, the full range of cool core strengths and dark matter concentrations, and have large substructure fractions. Taken together, these results imply that the amplitude of the luminosity gap is a function of both the formation epoch, and the recent infall history of the cluster. "BCG dominance" is therefore a phase that a cluster may evolve through, and is not an evolutionary "cul-de-sac". We also compare our results with semi-analytic model predictions based on the Millennium Simulation. None of the models are able to reproduce all of the observational results, underlining the inability of current models to match the empirical properties of BCGs. We identify the strength of AGN feedback and the efficiency with which cluster galaxies are replenished after they merge with the BCG in each model as possible causes of these discrepancies. [Abridged]

We present the AGN, star-forming, and morphological properties of a sample of 13 MIR-luminous (f(24) > 700uJy) IR-bright/optically-faint galaxies (IRBGs, f(24)/f(R) > 1000). While these z~2 sources were drawn from deep Chandra fields with >200 ks X-ray coverage, only 7 are formally detected in the X-ray and four lack X-ray emission at even the 2 sigma level. Spitzer IRS spectra, however, confirm that all of the sources are AGN-dominated in the mid-IR, although half have detectable PAH emission responsible for ~25% of their mid-infrared flux density. When combined with other samples, this indicates that at least 30-40% of luminous IRBGs have star-formation rates in the ULIRG range (~100-2000 Msun/yr). X-ray hardness ratios and MIR to X-ray luminosity ratios indicate that all members of the sample contain heavily X-ray obscured AGN, 80% of which are candidates to be Compton-thick. Furthermore, the mean X-ray luminosity of the sample, log L(2-10 keV)(ergs/s)=44.6, indicates that these IRBGs are Type 2 QSOs, at least from the X-ray perspective. While those sources most heavily obscured in the X-ray are also those most likely to display strong silicate absorption in the mid-IR, silicate absorption does not always accompany X-ray obscuration. Finally, ~70% of the IRBGs are merger candidates, a rate consistent with that of sub-mm galaxies (SMGs), although SMGs appear to be physically larger than IRBGs. These characteristics are consistent with the proposal that these objects represent a later, AGN-dominated, and more relaxed evolutionary stage following soon after the star-formation-dominated one represented by the SMGs.

We present extragalactic number counts and a lower limit estimate for the cosmic infrared background at 15 um from AKARI ultra deep mapping of the gravitational lensing cluster Abell 2218. This data is the deepest taken by any facility at this wavelength, and uniquely samples the normal galaxy population. We have de-blended our sources, to resolve photometric confusion, and de-lensed our photometry to probe beyond AKARI's blank-field sensitivity. We estimate a de-blended 5 sigma sensitivity of 28.7 uJy. The resulting 15 um galaxy number counts are a factor of three fainter than previous results, extending to a depth of ~ 0.01 mJy and providing a stronger lower limit constraint on the cosmic infrared background at 15 um of 1.9 +/- 0.5 nW m^-2 sr^-1.

We present wide-field Herschel/PACS observations of Abell 1689, a massive galaxy cluster at z=0.1832, from our Open Time Key Programme. We detect 39 spectroscopically confirmed 100micron-selected cluster members down to 1.5x10^10 Lsun. These galaxies are forming stars at rates in the range 1-10 Msun/yr, and appear to comprise two distinct populations: two-thirds are unremarkable blue, late-type spirals found throughout the cluster; the remainder are dusty red sequence galaxies whose star formation is heavily obscured with A(Halpha)~2 mag, and are found only in the cluster outskirts. The specific-SFRs of these dusty red galaxies are lower than the blue late-types, suggesting that the former are in the process of being quenched, perhaps via pre-processing, the unobscured star formation being terminated first. We also detect an excess of 100micron-selected galaxies extending ~6 Mpc in length along an axis that runs NE-SW through the cluster centre at >95% confidence. Qualitatively this structure is consistent with previous reports of substructure in X-ray, lensing, and near-infrared maps of this cluster, further supporting the view that this cluster is a dynamically active, merging system.

We present an analysis of an 8 arcminute diameter map of the area around the galaxy cluster Abell 1835 from jiggle map observations at a wavelength of 1.1 mm using the Bolometric Camera (Bolocam) mounted on the Caltech Submillimeter Observatory (CSO). The data is well described by a model including an extended Sunyaev-Zel'dovich (SZ) signal from the cluster gas plus emission from two bright background submm galaxies magnified by the gravitational lensing of the cluster. The best-fit values for the central Compton value for the cluster and the fluxes of the two main point sources in the field: SMM J140104+0252, and SMM J14009+0252 are found to be $y_{0}=(4.34\pm0.52\pm0.69)\times10^{-4}$, 6.5$\pm{2.0}\pm0.7$ mJy and 11.3$\pm{1.9}\pm1.1$ mJy, where the first error represents the statistical measurement error and the second error represents the estimated systematic error in the result. This measurement assumes the presence of dust emission from the cluster's central cD galaxy of $1.8\pm0.5$ mJy, based on higher frequency observations of Abell 1835. The cluster image represents one of the highest-significance SZ detections of a cluster in the positive region of the thermal SZ spectrum to date. The inferred central intensity is compared to other SZ measurements of Abell 1835 and this collection of results is used to obtain values for $y_{0} = (3.60\pm0.24)\times10^{-4}$ and the cluster peculiar velocity $v_{z} = -226\pm275$ km/s.

The Sunyaev-Zel'dovich (SZ) effect is a spectral distortion of the cosmic microwave background as observed through the hot plasma in galaxy clusters. This distortion is a decrement in the CMB intensity for lambda > 1.3 mm, an increment at shorter wavelengths, and small again by lambda ~250 um. As part of the Herschel Lensing Survey (HLS) we have mapped 1E0657-56 (the Bullet cluster) with SPIRE with bands centered at 250, 350 and 500 um and have detected the SZ effect at the two longest wavelengths. The measured SZ effect increment central intensities are Delta I_{0} = 0.097 +- 0.019 MJy sr^{-1} at 350 um and Delta I_{0} = 0.268 +- 0.031 MJy sr^{-1} at 500 um, consistent with the SZ effect spectrum derived from previous measurements at 2 mm. No other diffuse emission is detected. The presence of the finite temperature SZ effect correction is preferred by the SPIRE data at a significance of 2.1 sigma, opening the possibility that the relativistic SZ effect correction can be constrained by SPIRE in a sample of clusters. The results presented here have important ramifications for both sub-mm measurements of galaxy clusters and blank field surveys with SPIRE.

We present Herschel/PACS, MMT/Hectospec and XMM-Newton observations of Abell 1835, one of the brightest X-ray clusters on the sky, and the host of a strong cool core. Even though Abell 1835 has a prototypically "relaxed" X-ray morphology and no signs of ongoing merger activity in strong- and weak-lensing mass maps, it has a complex velocity distribution, suggesting that it is still accreting significant amounts of mass in the form of smaller satellite systems. Indeed, we find strong dynamical segregation of star-forming dusty galaxies from the optically selected cluster population. Most Herschel sources are found close to the virial radius of the cluster, and almost a third appear to be embedded within a filament feeding the cluster from the SW. We find that the most luminous infrared galaxies are likely involved in galaxy-galaxy interactions that may have triggered the current phase of star formation.

We present an early broad-brush analysis of Herschel/PACS observations of star-forming galaxies in 8 galaxy clusters drawn from our survey of 30 clusters at z~0.2. We define a complete sample of 192 spectroscopically confirmed cluster members down to L_TIR>3x10^10L_sun and L_K>0.25L_sun. The average K-band and bolometric infrared luminosities of these galaxies both fade by a factor of ~2 from clustercentric radii of 2r_200 to ~0.5r_200, indicating that as galaxies enter the clusters ongoing star-formation stops first in the most massive galaxies, and that the specific star-formation rate (SSFR) is conserved. On smaller scales the average SSFR jumps by 25%, suggesting that in cluster cores processes including ram pressure stripping may trigger a final episode of star-formation that presumably exhausts the remaining gas. This picture is consistent with our comparison of the Herschel-detected cluster members with the cluster mass distributions, as measured in our previous weak-lensing study of these clusters. For example, the spatial distribution of the Herschel sources is positively correlated with the structures in the weak-lensing mass maps at 5sigma significance, with the strongest signal seen at intermediate group-like densities. The strong dependence of the total cluster IR luminosity on cluster mass L_TIR propto M_virial^2 is also consistent with accretion of galaxies and groups of galaxies (i.e. the substructure mass function) driving the cluster IR luminosity. The most surprising result is that roughly half of the Herschel-detected cluster members have redder S_100/S_24 flux ratios than expected, based on the Rieke et al. models. On average cluster members are redder than non-members, and the fraction of red galaxies increases towards the cluster centers, both of which indicate that these colors are not attributable to systematic photometric errors. [Abridged]

The Herschel Lensing Survey (HLS) will conduct deep PACS and SPIRE imaging of ~40 massive clusters of galaxies. The strong gravitational lensing power of these clusters will enable us to penetrate through the confusion noise, which sets the ultimate limit on our ability to probe the Universe with Herschel. Here, we present an overview of our survey and a summary of the major results from our Science Demonstration Phase (SDP) observations of the Bullet Cluster (z=0.297). The SDP data are rich, allowing us to study not only the background high-redshift galaxies (e.g., strongly lensed and distorted galaxies at z=2.8 and 3.2) but also the properties of cluster-member galaxies. Our preliminary analysis shows a great diversity of far-infrared/submillimeter spectral energy distributions (SEDs), indicating that we have much to learn with Herschel about the properties of galaxy SEDs. We have also detected the Sunyaev-Zel'dovich (SZ) effect increment with the SPIRE data. The success of this SDP program demonstrates the great potential of the Herschel Lensing Survey to produce exciting results in a variety of science areas.

The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 micron bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 micron and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 micron selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities.This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.

We use deep, five band (100-500um) data from the Herschel Lensing Survey (HLS) to fully constrain the obscured star formation rate, SFR_FIR, of galaxies in the Bullet cluster (z=0.296), and a smaller background system (z=0.35) in the same field. Herschel detects 23 Bullet cluster members with a total SFR_FIR = 144 +/- 14 M_sun yr^-1. On average, the background system contains brighter far-infrared (FIR) galaxies, with ~50% higher SFR_FIR (21 galaxies; 207 +/- 9 M_sun yr^-1). SFRs extrapolated from 24um flux via recent templates (SFR_24) agree well with SFR_FIR for ~60% of the cluster galaxies. In the remaining ~40%, SFR_24 underestimates SFR_FIR due to a significant excess in observed S_100/S_24 (rest frame S_75/S_18) compared to templates of the same FIR luminosity.

We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint Herschel sources near the confusion noise. The method is based on the use of the positions of Spitzer/MIPS 24 um sources as priors. Our catalogs are able to reliably (5 sigma) recover galaxies with fluxes above 6 and 10 mJy in the PACS 100 and 160 um channels, respectively, and 12 to 18 mJy in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments.

The dust destruction timescales in the cores of clusters of galaxies are relatively short given their high central gas densities. However, substantial mid-infrared and sub-mm emission has been detected in many brightest cluster galaxies. In this letter we present Herschel PACS and SPIRE photometry of the brightest cluster galaxy in three strong cooling flow clusters, A1068, A2597 and Zw3146. This photometry indicates that a substantial mass of cold dust is present (>3 x 10^7 Mo) at temperatures significantly lower (20-28K) than previously thought based on limited MIR and/or sub-mm results. The mass and temperature of the dust appear to match those of the cold gas traced by CO with a gas-to-dust ratio of 80-120.

The question of how much gas cools in the cores of clusters of galaxies has been the focus of many, multiwavelength studies in the past 30 years. In this letter we present the first detections of the strongest atomic cooling lines, [C II], [O I] and [N I] in two strong cooling flow clusters, A1068 and A2597, using Herschel PACS. These spectra indicate that the substantial mass of cold molecular gas (>10^9 Mo) known to be present in these systems is being irradiated by intense UV radiation, most probably from young stars. The line widths of these FIR lines indicate that they share dynamics similar but not identical to other ionised and molecular gas traced by optical, near-infrared and CO lines. The relative brightness of the FIR lines compared to CO and FIR luminosity is consistent with other star-forming galaxies indicating that the properties of the molecular gas clouds in cluster cores and the stars they form are not unusual. These results provide additional evidence for a reservoir of cold gas that is fed by the cooling of gas in the cores of the most compact clusters and provide important diagnostics of the temperature and density of the dense clouds this gas resides in.

The most distant quasars known, at redshifts z=6, generally have properties indistinguishable from those of lower-redshift quasars in the rest-frame ultraviolet/optical and X-ray bands. This puzzling result suggests that these distant quasars are evolved objects even though the Universe was only seven per cent of its current age at these redshifts. Recently one z=6 quasar was shown not to have any detectable emission from hot dust, but it was unclear whether that indicated different hot-dust properties at high redshift or if it is simply an outlier. Here we report the discovery of a second quasar without hot-dust emission in a sample of 21 z=6 quasars. Such apparently hot-dust-free quasars have no counterparts at low redshift. Moreover, we demonstrate that the hot-dust abundance in the 21 quasars builds up in tandem with the growth of the central black hole, whereas at low redshift it is almost independent of the black hole mass. Thus z=6 quasars are indeed at an early evolutionary stage, with rapid mass accretion and dust formation. The two hot-dust-free quasars are likely to be first-generation quasars born in dust-free environments and are too young to have formed a detectable amount of hot dust around them.

Spitzer spectroscopy has revealed that ~80% of submm galaxies (SMGs) are starburst (SB) dominated in the mid-infrared. Here we focus on the remaining ~20% that show signs of harboring powerful active galactic nuclei (AGN). We have obtained Spitzer-IRS spectroscopy of a sample of eight SMGs which are candidates for harboring powerful AGN on the basis of IRAC color-selection (S8/S4.5>2; i.e. likely power-law mid-infrared SEDs). SMGs with an AGN dominating (>50%) their mid-infrared emission could represent `missing link' sources in an evolutionary sequence involving a major merger. First of all, we detect PAH features in all of the SMGs, indicating redshifts from 2.5-3.4, demonstrating the power of the mid-infrared to determine redshifts for these optically faint dusty galaxies. Secondly, we see signs of both star-formation (from the PAH features) and AGN activity (from continuum emission) in our sample: 62% of the sample are AGN-dominated in the mid-infrared with a median AGN content of 56%, compared with <30% on average for typical SMGs, revealing that our IRAC color selection has successfully singled out sources with proportionately more AGN emission than typical SB-dominated SMGs. However, we find that only about 10% of these AGN dominate the bolometric emission of the SMG when the results are extrapolated to longer infrared wavelengths, implying that AGN are not a significant power source to the SMG population overall, even when there is evidence in the mid-infrared for substantial AGN activity. When existing samples of mid-infrared AGN-dominated SMGs are considered, we find that S8/S4.5>1.65 works well at selecting mid-infrared energetically dominant AGN in SMGs, implying a duty cycle of ~15% if all SMGs go through a subsequent mid-infrared AGN-dominated phase in the proposed evolutionary sequence.

We present a statistical analysis of a sample of 20 strong lensing clusters drawn from the Local Cluster Substructure Survey (LoCuSS), based on high resolution Hubble Space Telescope imaging of the cluster cores and follow-up spectroscopic observations using the Keck-I telescope. We use detailed parameterized models of the mass distribution in the cluster cores, to measure the total cluster mass and fraction of that mass associated with substructures within R<250kpc.These measurements are compared with the distribution of baryons in the cores, as traced by the old stellar populations and the X-ray emitting intracluster medium. Our main results include: (i) the distribution of Einstein radii is log-normal, with a peak and 1sigma width of <log(RE(z=2))>=1.16+/-0.28; (ii) we detect an X-ray/lensing mass discrepancy of <M_SL/M_X>=1.3 at 3 sigma significance -- clusters with larger substructure fractions displaying greater mass discrepancies, and thus greater departures from hydrostatic equilibrium; (iii) cluster substructure fraction is also correlated with the slope of the gas density profile on small scales, implying a connection between cluster-cluster mergers and gas cooling. Overall our results are consistent with the view that cluster-cluster mergers play a prominent role in shaping the properties of cluster cores, in particular causing departures from hydrostatic equilibrium, and possibly disturbing cool cores. Our results do not support recent claims that large Einstein radius clusters present a challenge to the CDM paradigm.

We present photometry, photometric redshifts and extra galactic number counts for ultra deep 15 micron mapping of the gravitational lensing cluster Abell 2218 (A2218), which is the deepest image taken by any facility at this wavelength. This data resolves the cosmic infrared background (CIRB) beyond the 80% that blank field AKARI surveys aim to achieve. To gain an understanding of galaxy formation and evolution over the age of the Universe a necessary step is to fully resolve the CIRB, which represents the dust-shrouded cosmic star formation history. Observing through A2218 gives magnifications of up to a factor of 10, thus allowing the sampling of a more representative spread of high redshift galaxies, which comprise the bulk of the CIRB. 19 pointed observations were taken by AKARIs IRC MIR-L channel, and a final combined image with an area of 122.3 square arcminutes and effective integration time of 8460 seconds was achieved. The 5 sigma sensitivity limit is estimated at 41.7 uJy. An initial 5 sigma catalogue of 565 sources was extracted giving 39 beams per source, which shows the image is confusion limited. Our 15 micron number counts show strong evolution consistent with galaxy evolution models that incorporate downsizing in star formation.

We present the results of a search for galaxy clusters in Subaru-XMM Deep Field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 2x10^{-15} ergs cm^{-2} s^{-1} over one of the widest XMM-Newton contiguous raster surveys, covering an area of 1.3 square degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. The red sequence technique allows us to identify 57 cluster candidates. We report on the progress with the cluster spectroscopic follow-up and derive their properties based on the X-ray luminosity and cluster scaling relations. In addition, 3 sources are identified as X-ray counterparts of radio lobes, and in 3 further sources, X-ray counterpart of radio lobes provides a significant fraction of the total flux of the source. In the area covered by NIR data, our identification success rate achieves 86%. We detect a number of radio galaxies within our groups and for a luminosity-limited sample of radio galaxies we compute halo occupation statistics using a marked cluster mass function. We compare the cluster detection statistics in the SXDF with the predictions of concordance cosmology and current knowledge of the X-ray cluster properties, concluding that a reduction of concordance sigma_8 value by 5% is required in order to match the prediction of the model and the data. This conclusion still needs verification through the completion of cluster follow-up.

We present the identification of a bright submm source, SMMJ163555.5+661300, detected in the lensing cluster Abell2218, for which we have accurately determined the position using observations from the Submillimeter Array (SMA). The identified optical counterpart has a spectroscopic redshift of z=4.044+-0.001 if we attribute the single emission line detected at lambda=6140AA to Lyman-alpha. This redshift identification is in good agreement with the optical/near-infrared photometric redshift as well as the submm flux ratio S_450/S_850~1.6, the radio-submm flux ratio S_1.4/S_850 < 0.004, and the 24um to 850um flux ratio S_24/S_850 < 0.005. Correcting for the gravitational lensing amplification of ~5.5, we find that the source has a far-infrared luminosity of 1.3x10^12 Lsun, which implies a star formation rate of 230 Msun/yr. This makes it the lowest-luminosity SMG known at z>4 to date. Previous CO(4-3) emission line obserations yielded a non-detection, for which we derived an upper limit of the CO line luminosity of L'_CO = 0.3x10^10 K km/s/pc^2, which is not inconsistent with the L'_CO - L_FIR relation for starburst galaxies. The best fit model to the optical and near-infrared photometry give a stellar population with an age of 1.4 Gyr and a stellar mass of 1.6x10^10 Msun. The optical morphology is compact and in the source plane the galaxy has an extent of ~6kpc x 3kpc with individual star forming knots of <500 pc in size. J163556 is not resolved in the SMA data and we place a strict upper limit on the size of the starburst region of 8kpc x 3kpc, which implies a lower limit on the star formation rate surface density of 12 Msun/yr/kpc^2. The redshift of J163556 extends the redshift distribution of faint, lensed SMGs, and we find no evidence that these have a different redshift distribution than bright SMGs.

We report the detection of the Paschen-alpha emission line in the z=2.515 galaxy SMM J163554.2+661225 using Spitzer spectroscopy. SMM J163554.2+661225 is a sub-millimeter-selected infrared (IR)-luminous galaxy maintaining a high star-formation rate (SFR), with no evidence of an AGN from optical or infrared spectroscopy, nor X-ray emission. This galaxy is lensed gravitationally by the cluster Abell 2218, making it accessible to Spitzer spectroscopy. Correcting for nebular extinction derived from the H-alpha and Pa-alpha lines, the dust-corrected luminosity is L(Pa-alpha) = (2.57+/-0.43) x 10^43 erg s^-1, which corresponds to an ionization rate, Q = (1.6+/-0.3) x 10^55 photons s^-1. The instantaneous SFR is 171+/-28 solar masses per year, assuming a Salpeter-like initial mass function. The total IR luminosity derived using 70, 450, and 850 micron data is L(IR) = (5-10) x 10^11 solar luminosities, corrected for gravitational lensing. This corresponds to a SFR=90-180 solar masses per year, where the upper range is consistent with that derived from the Paschen-alpha luminosity. While the L(8 micron) / L(Pa-alpha) ratio is consistent with the extrapolated relation observed in local galaxies and star-forming regions, the rest-frame 24 micron luminosity is significantly lower with respect to local galaxies of comparable Paschen-alpha luminosity. Thus, SMM J163554.2+661225 arguably lacks a warmer dust component (T ~ 70 K), which is associated with deeply embedded star formation, and which contrasts with local galaxies with comparable SFRs. Rather, the starburst is consistent with star-forming local galaxies with intrinsic luminosities, L(IR) ~ 10^10 solar luminosities, but "scaled-up" by a factor of 10-100.

We present the first results from the largest deep extragalactic millimetre-wavelength survey undertaken to date. These results are derived from maps covering over 0.7 deg^2, made at 1.1mm, using the AzTEC continuum camera mounted on the James Clerk Maxwell Telescope. The maps were made in the two fields originally targeted at 0.85mm with SCUBA in the SHADES project, namely the Lockman Hole East (mapped to a depth of 0.9-1.3 mJy rms) and the Subaru XMM Deep Field (1.0-1.7 mJy rms). The wealth of existing and forthcoming deep multi-frequency data in these two fields will allow the bright mm source population revealed by these images to be explored in detail in subsequent papers. Here we present the maps themselves, a catalogue of 114 high-significance sub-millimetre galaxy detections, and a thorough statistical analysis leading to the most robust determination to date of the 1.1mm source number counts. Through careful comparison, we find that both the COSMOS and GOODS North fields, also imaged with AzTEC, contain an excess of mm sources over the new 1.1mm source-count baseline established here. In particular, our new AzTEC/SHADES results indicate that very luminous high-redshift dust enshrouded starbursts (S_{1.1} > 3 mJy) are 25-50% less common than would have been inferred from these smaller surveys, thus highlighting the potential roles of cosmic variance and clustering in such measurements. We compare number count predictions from recent models of the evolving mm/sub-mm source population to these SMG surveys, which provide important constraints for the ongoing refinement of semi-analytic and hydrodynamical models of galaxy formation, and find that all recent models over-predict the number of bright sub-millimetre galaxies found in this survey.

We present 20 band photometry from the far-IR to X-ray in the Spitzer IRAC dark field. The bias for the near-IR camera on Spitzer is calibrated by observing a ~20 arcminute diameter "dark" field near the north ecliptic pole roughly every two-to-three weeks throughout the mission duration of Spitzer. The field is unique for its extreme depth, low background, high quality imaging, time-series information, and accompanying photometry including data taken with Akari, Palomar, MMT, KPNO, Hubble, and Chandra. This serendipitous survey contains the deepest mid-IR data taken to date. This dataset is well suited for studies of intermediate redshift galaxy clusters, high redshift galaxies, the first generation of stars, and the lowest mass brown dwarfs, among others. This paper provides a summary of the data characteristics and catalog generation from all bands collected to date as well as a discussion of photometric redshifts and initial and expected science results and goals. To illustrate the scientific potential of this unique dataset, we also present here IRAC color color diagrams.

We used the large binocular camera (LBC) mounted on the large binocular telescope (LBT) to observe the Lockman Hole in the U, B, and V bands. Our observations cover an area of 925 sq.arcmin. We reached depths of 26.7, 26.3, and 26.3 mag(AB) in the three bands, respectively, in terms of 50% source detection efficiency, making this survey the deepest U-band survey and one of the deepest B and V band surveys with respect to its covered area. We extracted a large number of sources (~89000), detected in all three bands and examined their surface density, comparing it with models of galaxy evolution. We find good agreement with previous claims of a steep faint-end slope of the luminosity functions, caused by late-type and irregular galaxies at z>1.5. A population of dwarf star-forming galaxies at 1.5<z<2.5 is needed to explain the U-band number counts. We also find evidence of strong supernova feedback at high redshift. This survey is complementary to the r, i, and z Lockman Hole survey conducted with the Subaru telescope and provides the essential wavelength coverage to derive photometric redshifts and select different types of sources from the Lockman Hole for further study.

We study the mid-infrared (MIR) properties of galaxies in 30 massive galaxy clusters at 0.02<z<0.40, using panoramic Spitzer/MIPS 24micron and NIR data. This is the largest sample of clusters to date with MIR data covering not only the cluster cores, but extending into the infall regions. We revisit the Butcher-Oemler effect, measuring the fraction of massive infrared-luminous galaxies (K<K*+1.5, L_IR>5x10^10L_sun) within r_200, finding a steady increase in the fraction with redshift from ~3% at z=0.02 to ~10% by z=0.30, and an rms cluster-to-cluster scatter about this trend of 0.03. The best-fit redshift evolution model is of the form f_SF ~ (1+z)^5.7, which is stronger redshift evolution than that of L*_IR in both clusters and the field. We find that, statistically, this excess is associated with galaxies found at large cluster-centric radii, implying that the MIR Butcher-Oemler effect can be explained by a combination of both the global decline in star-formation in the universe since z~1 and enhanced star formation in the infall regions of clusters at intermediate redshifts. This picture is supported by a simple infall model based on the Millennium Simulation semi-analytic galaxy catalogs, whereby star-formation in infalling galaxies is instantaneously quenched upon their first passage through the cluster, in that the observed radial trends of f_SF trace those inferred from the simulations. We also find that f_SF does not depend on simple indicators of the dynamical state of clusters, including the offset between the brightest cluster galaxy and the peak of the X-ray emission. This is consistent with the picture described above in that most new star-formation in clusters occurs in the infall regions, and is thus not sensitive to the details of cluster-cluster mergers in the core regions.

We present the results from an optical and near-infrared spectroscopic study of the ultraviolet-luminous z = 2.73 galaxy, the 8 o'clock arc. Due to gravitational lensing, this galaxy is magnified by a factor of > 10, allowing in-depth measurements which are usually unfeasible at such redshifts. In the optical spectra, we measured the systemic redshift of the galaxy, z = 2.7322 +/- 0.0012, using stellar photospheric lines. This differs from the redshift of absorption lines in the interstellar medium, z = 2.7302 +/- 0.0006, implying gas outflows on the order of 160 km/s. With H and K-band near-infrared spectra, we have measured nebular emission lines of Halpha, Hbeta, Hgamma, [N II] and [O III], which have a redshift z = 2.7333 +/- 0.0001, consistent with the derived systemic redshift. From the Balmer decrement, we measured the dust extinction to be A_5500 = 1.17 +/- 0.36 mag. Correcting Halpha for dust extinction and the assumed lensing factor, we measure a star-formation rate of ~ 270 Msol/yr, which is higher than ~ 85% of star-forming galaxies at z ~ 2-3. Using combinations of all detected emission lines, we find that the 8 o'clock arc has a gas-phase metallicity of ~ 0.8 Zsol, showing that enrichment at high-redshift is not rare, even in blue, star-forming galaxies. Studying spectra from two of the arc components separately, we find that one component dominates the dust extinction and star-formation rate, although the metallicities between the two components are similar. We derive the mass via stellar population modeling, and find that the arc has a total stellar mass of ~ 4 x 10^11 Msol, which falls on the mass-metallicity relation at z ~ 2. Finally, we estimate the total gas mass, and find it to be only ~ 12% of the stellar mass, implying that the 8 o'clock arc is likely nearing the end of a starburst.

We analyze a sample of galaxies chosen to have F24>0.5 mJy and satisfy a certain IRAC color criterion. IRS spectra yield redshifts, spectral types, and PAH luminosities, to which we add multi-wavelength broadband photometry. Stellar population modeling and IRS spectra together demonstrate that the double criteria used to select this sample have efficiently isolated massive star-forming galaxies at z~1.9. This is the first starburst-dominated ULIRG sample at high redshift with total infrared luminosity measured directly from FIR and millimeter photometry, and as such gives us the first accurate view of broadband SEDs for starburst galaxies at extremely high luminosity and at all wavelengths. Similar broadband data are assembled for three other galaxy samples -- local starburst galaxies, local AGN/ULIRGS, and a second 24mu-luminous z~1.9 sample dominated by AGN. L(PAH)/L(IR) for the new z~1.9 starburst sample is the highest ever seen, some three times higher than in local starbursts, whereas in AGNs this ratio is depressed below the starburst trend, often severely. Several pieces of evidence imply that AGNs exist in this starburst dominated sample, except two of which even host very strong AGN, while they still have very strong PAH emission. The ACS images show most objects have very extended morphologies in the rest-frame UV band, thus extended distribution of PAH molecules. Such an extended distribution prevents further destruction PAH molecules by central AGNs. We conclude that objects in this sample are ULIRGs powered mainly by starburst; and the total infrared luminosity density contributed by this type of objects is 0.9-2.6x 10^7 Lsol/Mpc^3.

We present the first galaxy evolution results from the Local Cluster Substructure Survey (LoCuSS), a multi-wavelength survey of 100 X-ray selected galaxy clusters at 0.15<z<0.30. LoCuSS combines far-UV through far-IR observations of cluster galaxies with gravitational lensing analysis and X-ray data to investigate the interplay between the hierarchical assembly of clusters and the evolution of cluster galaxies. Here we present new panoramic Spitzer/MIPS 24micron observations of the merging cluster Abell 1758 at z=0.279 spanning 6.5x6.5Mpc and reaching a 90% completeness limit of 400uJy. We estimate a global cluster SFR of 910\pm320 M_sun/yr within 3 Mpc of the cluster centre, originating from 42 galaxies with L_IR > 5x10^10 L_sun. The obscured activity in A1758 is therefore comparable with that in Cl 0024+1654, the most active cluster previously studied at 24um. The obscured galaxies faithfully trace the cluster potential as revealed by the weak-lensing mass map of the cluster, including numerous mass peaks at R~2-3Mpc that are likely associated with infalling galaxy groups and filamentary structures. However the core (R<500kpc) of A1758N is 2x more active in the IR than that of A1758S, likely reflecting differences in the recent dynamical history of the two clusters. The 24micron results from A1758 therefore suggest that dust-obscured cluster galaxies are common in merging clusters and suggests that obscured activity in clusters is triggered by both the details of cluster-cluster mergers and processes that operate at larger radii including those within in-falling groups. Our ongoing far-UV through far-IR observations of a large sample of clusters should allow us to disentangle the different physical processes responsible for triggering obscured star formation in clusters.

Probing the growth of structure from the epoch of hydrogen recombination to the formation of the first stars and galaxies is one of the most important uncharted areas of observational cosmology. Far-IR spectroscopy covering $\lambda$ 100-500 microns from space, and narrow partial transmission atmospheric bands available from the ground, opens up the possibility of probing the molecular hydrogen and metal fine-structure lines from primordial clouds from which the first stars and galaxies formed at 6 < z $<$ 15. Building on Spitzer observations of unexpectedly powerful H2 emission from shocks, we argue that next-generation far-IR space telescopes may open a new window into the main cloud cooling processes and feedback effects which characterized this vital, but unexplored epoch. Without this window, we are essential blind to the dominant cloud cooling which inevitably led to star formation and cosmic reionization.

We report the discovery of a giant Lya emitter (LAE) with a Spitzer/IRAC counterpart near the reionization epoch at z=6.595. The giant LAE is found from the extensive 1 deg^2 Subaru narrow-band survey for z=6.6 LAEs in the Subaru/XMM-Newton Deep Survey (SXDS) field, and subsequently identified by deep spectroscopy of Keck/DEIMOS and Magellan/IMACS. Among our 207 LAE candidates, this LAE is not only the brightest narrow-band object with L(Lya) = 3.9+/-0.2 x 10^43 erg/s in our survey volume of 10^6 Mpc^3, but also a spatially extended Lya nebula with the largest isophotal area whose major axis is at least ~3''. This object is more likely to be a large Lya nebula with a size of >~17 kpc than to be a strongly-lensed galaxy by a foreground object. Our Keck spectrum with medium-high spectral and spatial resolutions suggests that the velocity width is v(FWHM)=251+/-21 km/s, and that the line-center velocity changes by ~60 km/s in a 10-kpc range. The stellar mass and star-formation rate are estimated to be 0.9-5.0 x 10^10 Mo and >34 Mo/yr, respectively, from the combination of deep optical to infrared images of Subaru, UKIDSS-Ultra Deep Survey, and Spitzer/IRAC. Although the nature of this object is not yet clearly understood, this could be an important object for studying cooling clouds accreting onto a massive halo, or forming-massive galaxies with significant outflows contributing to cosmic reionization and metal enrichment of inter-galactic medium.

We present the results of a Spitzer IRAC and MIPS 24 micron study of extended Lyman-alpha clouds (or Lyman-alpha Blobs, LABs) within the SSA22 filamentary structure at z = 3.09. We detect 6/26 LABs in all IRAC filters, four of which are also detected at 24 micron, and find good correspondence with the 850 micron measurements of Geach et al. 2005. An analysis of the rest-frame ultraviolet, optical, near- and mid-infrared colors reveals that these six systems exhibit signs of nuclear activity (AGN)and/or extreme star formation. Notably, they have properties that bridge galaxies dominated by star formation (Lyman-break galaxies; LBGs) and those with AGNs (LBGs classified as QSOs). The LAB systems not detected in all four IRAC bands, on the other hand, are, as a group, consistent with pure star forming systems, similar to the majority of the LBGs within the filament. These results indicate that the galaxies within LABs do not comprise a homogeneous population, though they are also consistent with scenarios in which the gas halos are ionized through a common mechanism such as galaxy-scale winds driven by the galaxies within them, or gravitational heating of the collapsing cloud itself.

We present a high signal-to-noise spectrum of a bright galaxy at z = 4.9 in 14 h of integration on VLT FORS2. This galaxy is extremely bright, i_850 = 23.10 +/- 0.01, and is strongly-lensed by the foreground massive galaxy cluster Abell 1689 (z=0.18). Stellar continuum is seen longward of the Ly-alpha emission line at ~7100 \AA, while intergalactic H I produces strong absorption shortward of Ly-alpha. Two transmission spikes at ~6800 Angstroms (A) and ~7040 A are also visible, along with other structures at shorter wavelengths. Although fainter than a QSO, the absence of a strong central ultraviolet flux source in this star forming galaxy enables a measurement of the H I flux transmission in the intergalactic medium (IGM) in the vicinity of a high redshift object. We find that the effective H I optical depth of the IGM is remarkably high within a large 14 Mpc (physical) region surrounding the galaxy compared to that seen towards QSOs at similar redshifts. Evidently, this high-redshift galaxy is located in a region of space where the amount of H I is much larger than that seen at similar epochs in the diffuse IGM. We argue that observations of high-redshift galaxies like this one provide unique insights on the nascent stages of baryonic large-scale structures that evolve into the filamentary cosmic web of galaxies and clusters of galaxies observed in the present universe.

We present the results of a systematic search for gravitationally-lensed continuum Lyman break `drop-outs' beyond a redshift 7 conducted via very deep imaging through six foreground clusters undertaken with the Hubble and Spitzer Space Telescopes. The survey has yielded 10 z-band and 2 J-band drop-out candidates to photometric limits of J_110~=26.2 AB (5sigma). Taking into account the magnifications afforded by our clusters (1-4 magnitudes), we probe the presence of z>7 sources to unlensed limits of J_{110}~=30 AB, fainter than those charted in the Hubble Ultradeep Field. To verify the fidelity of our candidates we conduct a number of tests for instrumental effects which would lead to spurious detections, and carefully evaluate the likelihood of foreground contamination by considering photometric uncertainties in the drop-out signature, the upper limits from stacked IRAC data and the statistics of multiply-imaged sources. Overall, we conclude that we can expect about half of our sample of z-band drop-outs are likely to be at high redshift. An ambitious infrared spectroscopic campaign undertaken with the NIRSPEC spectrograph at the WM Keck Observatory for seven of the most promising candidates failed to detect any Lyman-alpha emission highlighting the challenge of making further progress in this field. While the volume density of high redshift sources will likely remain uncertain until more powerful facilities are available, our data provides the first potentially interesting constraints on the UV luminosity function at z~=7.5 at intrinsically faint limits. We discuss the implications of our results in the context of the hypothesis that the bulk of the reionizing photons in the era 7<z<12 arise in low luminosity galaxies undetected by conventional surveys.

In this article, we apply strong lensing techniques in Abell 1703, a massive X-ray luminous galaxy cluster at z=0.28. Our analysis is based on imaging data both from space and ground in 8 bands, complemented with a spectroscopic survey. Abell 1703 looks rather circular from the general shape of its multiply imaged systems and present a dominant giant elliptical cD galaxy in its centre. This cluster exhibits a remarkable bright 'central ring' formed by 4 bright images at z_{spec}=0.888 located very close to the cD galaxy, providing observational constraints that are potentially very interesting to probe the central mass distribution. The stellar contribution from the cD galaxy (~1.25 10^{12} M_{sun} within 7") is accounted for in our parametric mass modelling, and the underlying smooth dark matter component distribution is described using a generalized NFW profile parametrized with a central logarithmic slope \alpha. We find that within the range where observational constraints are present (from ~5" to ~50"), the slope of the dark matter distribution in Abell1703 is equal to 1.09^{+0.05}_{-0.11} (3\sigma confidence level). The concentration parameter is equal to c_{200} ~ 3.5, and the scale radius is constrained to be larger than the region where observational constraints are available. Within this radius, the 2D mass is equal to M(50")=2.4 10^{14} M_{\sun}. We cannot draw any conclusions on cosmological models at this point since we lack results from realistic numerical simulations containing baryons to make a proper comparison. We advocate the need for a sample of observed and simulated unimodal relaxed galaxy clusters in order to make reliable comparisons, and potentially provide a test of cosmological models.

Using a sample of ~28,000 sources selected at 3.6-4.5 microns with Spitzer observations of the HDF-N, the CDF-S, and the Lockman Hole (surveyed area: ~664 arcmin^2), we study the evolution of the stellar mass content of the Universe at 0<z<4. We calculate stellar masses and photometric redshifts, based on ~2,000 templates built with stellar and dust emission models fitting the UV-to-MIR SEDs of galaxies with spectroscopic redshifts. We estimate stellar mass functions for different redshift intervals. We find that 50% of the local stellar mass density was assembled at 0<z<1 (average SFR:0.048 M_sun/yr/Mpc^3), and at least another 40% at 1<z<4 (average SFR: 0.074 M_sun/yr/Mpc^3). Our results confirm and quantify the ``downsizing'' scenario of galaxy formation. The most massive galaxies (M>10^12.0 M_sun) assembled the bulk of their stellar content rapidly (in 1-2 Gyr) beyond z~3 in very intense star formation events (producing high specific SFRs). Galaxies with 10^11.5<M/M_sun<10^12.0 assembled half of their stellar mass before z~1.5, and more than 90% of their mass was already in place at z~0.6. Galaxies with M<10^11.5 M_sun evolved more slowly (presenting smaller specific SFRs), assembling half of their stellar mass below z~1. About 40% of the local stellar mass density of 10^9.0<M/M_sun<10^11.0 galaxies was assembled below z~0.4, most probably through accretion of small satellites producing little star formation. The cosmic stellar mass density at z>2.5 is dominated by optically faint (R>25) red galaxies (Distant Red Galaxies or BzK sources) which account for ~30% of the global population of galaxies, but contribute at least 60% to the cosmic stellar mass density. Bluer galaxies (e.g., Lyman Break Galaxies) are more numerous but less massive, contributing less than 50% to the global stellar mass density at high redshift.

Quillen et al.(2007) presented an imaging survey with the {\it Spitzer Space Telescope} of 62 brightest cluster galaxies with optical line emission located in the cores of X-ray luminous clusters. They found that at least half of these sources have signs of excess infrared emission. Here we discuss the nature of the IR emission and its implications for cool core clusters. The strength of the mid-IR excess emission correlates with the luminosity of the optical emission lines. Excluding the four systems dominated by an AGN, the excess mid-infrared emission in the remaining brightest cluster galaxies is likely related to star formation. The mass of molecular gas (estimated from CO observations) is correlated with the IR luminosity as found for normal star forming galaxies. The gas depletion time scale is about 1 Gyr. The physical extent of the infrared excess is consistent with that of the optical emission line nebulae. This supports the hypothesis that the star formation occurs in molecular gas associated with the emission line nebulae and with evidence that the emission line nebulae are mainly powered by ongoing star formation. We find a correlation between mass deposition rates (${\dot M}_X$) estimated from the X-ray emission and the star formation rate estimated from the infrared luminosity. The star formation rates are 1/10 to 1/100 of the mass deposition rates suggesting that the re-heating of the ICM is generally very effective in reducing the amount of mass cooling from the hot phase but not eliminating it completely.

We present a comparison between the SCUBA Half Degree Extragalactic Survey (SHADES) at 450 and 850 microns in the Lockman Hole East with a deep Spitzer Space Telescope survey at 3.6-24 microns conducted in Guaranteed Time. Using stacking analyses we demonstrate a striking correspondence between the galaxies contributing the submm extragalactic background light, with those likely to dominate the backgrounds at Spitzer wavelengths. Using a combination BRIzK plus Spitzer photometric redshifts, we show that at least a third of the Spitzer-identified submm galaxies at 1<z<1.5 appear to reside in overdensities when the density field is smoothed at 0.5-2 Mpc comoving diameters, supporting the high-redshift reversal of the local star formation - galaxy density relation. We derive the dust-shrouded cosmic star formation history of galaxies as a function of assembled stellar masses. For model stellar masses <10^11 Msun, this peaks at lower redshifts than the ostensible z~2.2 maximum for submm point sources, adding to the growing consensus for ``downsizing'' in star formation. Our surveys are also consistent with ``downsizing'' in mass assembly. Both the mean star formation rates <dM/dt> and specific star formation rates <(1/M)dM/dt> are in striking disagreement with some semi-analytic predictions from the Millenium simulation. The discrepancy could either be resolved with a top-heavy initial mass function, or a significant component of the submm flux heated by the interstellar radiation field.

We present estimates of the photometric redshifts, stellar masses and star formation histories of sources in the SCUBA HAlf Degree Extragalactic Survey (SHADES). This paper describes the 60 SCUBA sources detected in the Lockman Hole covering an area of ~320 square arcmin. Using photometry spanning the B band to 8um, we find that the average SCUBA source forms a significant fraction of its stars in an early period of star formation and that most of the remainder forms in a shorter more intense burst around the redshift it is observed. This trend does not vary significantly with source redshift but the exact ratio of early to late mass is quite sensitive to the way extinction is treated in the modelling. However, the sources show a clear increase in stellar mass with redshift, consistent with downsizing. In terms of SED types, only two out of the 51 sources we have obtained photometric redshifts for are best fit by a quasar-like spectral energy distribution, with approximately 80 per cent of the sources being best fit with late-type spectra (Sc, Im and starburst). By including photometry at 850um, we conclude that the average SCUBA source is forming stars at a rate somewhere between 6 and 30 times the rate implied from the rest-frame optical in a dust obscured burst and that this burst creates 15-65 per cent of the total stellar mass. Using a simplistic calculation, we estimate from the average star formation history that between one in five and one in 15 bright (L_* +2 mag < L_optical < L_* -1 mag) galaxies in the field over the interval 0 < z < 3 will at some point in their lifetime experience a similar energetic dusty burst of star formation. Finally, we compute the evolution of the star formation rate density and find it peaks around z=2.

We present a reconstruction of the mass distribution of galaxy cluster Abell 1689 at z = 0.18 using detected strong lensing features from deep HST/ACS observations and extensive ground based spectroscopy. Earlier analyses have reported up to 32 multiply imaged systems in this cluster, of which only 3 were spectroscopically confirmed. In this work, we present a parametric strong lensing mass reconstruction using 24 multiply imaged systems with newly determined spectroscopic redshifts, which is a major step forward in building a robust mass model. In turn, the new spectroscopic data allows a more secure identification of multiply imaged systems. The resultant mass model enables us to reliably predict the redshifts of additional multiply imaged systems for which no spectra are currently available, and to use the location of these systems to further constrain the mass model. In particular, we have detected 5 strong galaxy-galaxy lensing systems just outside the Einstein ring region, further constraining the mass profile. Our strong lensing mass model is consistent with that inferred from our large scale weak lensing analysis derived using CFH12k wide field images. Thanks to a new method for reliably selecting a well defined background lensed galaxy population, we resolve the discrepancy found between the strong and weak lensing mass models reported in earlier work. [ABRIDGED]

We present Spitzer/IRS mid-infrared spectra for 15 gravitationally lensed, 24 micron--selected galaxies, and combine the results with 4 additional very faint galaxies with IRS spectra in the literature. The median intrinsic 24 micron flux density of the sample is 130 microJy, enabling a systematic survey of the spectral properties of the very faint 24 micron sources that dominate the number counts of Spitzer cosmological surveys. Six of the 19 galaxy spectra (32%) show the strong mid-IR continuua expected of AGN; X-ray detections confirm the presence of AGN in three of these cases, and reveal AGNs in two other galaxies. These results suggest that nuclear accretion may contribute more flux to faint 24 micron--selected samples than previously assumed. Almost all the spectra show some aromatic (PAH) emission features; the measured aromatic flux ratios do not show evolution from z=0. In particular, the high S/N mid-IR spectrum of SMM J163554.2+661225 agrees remarkably well with low--redshift, lower--luminosity templates. We compare the rest-frame 8 micron and total infrared luminosities of star--forming galaxies, and find that the behavior of this ratio with total IR luminosity has evolved modestly from z=2 to z=0. Since the high aromatic--to--continuum flux ratios in these galaxies rule out a dominant contribution by AGN, this finding implies systematic evolution in the structure and/or metallicity of infrared sources with redshift. It also has implications for the estimates of star forming rates inferred from 24 micron measurements, in the sense that at z ~2, a given observed frame 24 micron luminosity corresponds to a lower bolometric luminosity than would be inferred from low-redshift templates of similar luminosity at the corresponding rest wavelength.

The gravitationally lensed Lyman-alpha emitting galaxy, HCM6A, detected by Hu et al. (2002) at z=6.56 behind the Abell 370 cluster was observed with the MAMBO-2 array of bolometers at 1.2mm wavelength. The galaxy was not detected down to 1.08 mJy (3 sigma), but the depth of the observations and the lens amplification allow us to improve by approximately one order of magnitude previously published upper limits on far infrared emission of Lyman-alpha emitting galaxies at this redshift. The following upper limits are derived from our observations assuming typical dust parameters: dust mass <5.3x10^7 Msun, IR luminosity <2.1x10^{11} Lsun, and star formation rate, SFR<35 Msun/yr. The observed restframe UV--optical--IR spectral energy distribution (SED) of this galaxy is compatible with that of normal spiral galaxies or blue compact dwarf galaxies. SEDs of prototypical ULIRGs, such as Arp 220, are clearly excluded. Finally, we obtain an upper limit of < 2.1x10^{-2} Msun/yr/Mpc^{-3} for the dust-obscured SFR density of Lyman-alpha selected galaxies at z~6.6.

We study the properties and nature of extremely red galaxies (ERO, R-K>5.6) found behind two lensing clusters and compare them with other known galaxy populations. New HST/ACS observations, Spitzer IRAC and MIPS, and Chandra/ACIS observations of the two lensing clusters Abell 1835 and AC114 contemplate our earlier optical and near-IR observations and have been used to study extremely red objects (EROs) in these deep fields. We have found 6 and 9 EROs in Abell 1835 and AC114. Several (7) of these objects are undetected up to the I and/or z band, and are hence ``optical'' drop-out sources. The photometric redshifts of most of our sources (80%) are z~0.7-1.5. According to simple colour-colour diagrams the majority of our objects would be classified as hosting old stellar populations. However, there are clear signs of dusty starbursts for several among them. These objects correspond to the most extreme ones in R-K colour. We estimate a surface density of (0.97+-0.31) arcmin-2 for EROs with (R-K>5.6) at K<20.5. Among our 15 EROs 6 (40 %) also classify as distant red galaxies (DRGs). 11 of 13 EROs with available IRAC photometry also fulfil the selection criteria for IRAC selected EROs (IEROs) of Yan et al. (2004). SED modelling shows that ~ 36 % of the IEROs in our sample are luminous or ultra-luminous infrared galaxies ((U)LIRG). Some very red DRGs are found to be very dusty starbursts, even (U)LIRGs, as also supported by their mid-IR photometry. No indication for AGNs is found, although faint activity cannot be excluded for all objects. From mid-IR and X-ray data 5 objects are clearly classified as starbursts. The derived properties are quite similar to those of DRGs and IEROs, except for 5 extreme objects in terms of colours, for which a very high extinction (Av>3) is found.

Determining an accurate position for a submm galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2-D clustering - to an assessment of their detailed, multi-wavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or IR counterparts, and hence accurate positions, for over two thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-um and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parameterisation appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6", or 15-50/(sin i) kpc at z ~ 2, consistent with early bursts seen in merger simulations.

Historically, variable and transient sources have both surprised astronomers and provided new views of the heavens. Here we report the discovery of an optical transient in the outskirts of the lenticular galaxy Messier 85 in the Virgo Cluster. With a peak absolute R magnitude of -12 this event is distinctly brighter than novae, but fainter than type Ia supernovae (expected from a population of old stars in lenticular galaxies). Archival images of the field do not show a luminous star at that position with an upper limit of ~-4.1, so it is unlikely to be a giant eruption from a luminous blue variable star. Over a two month period the transient emitted radiation energy of almost 10E47 erg and subsequently faded in the opical sky. It is similar to, but more luminous at peak by a factor of 6 than, an enigmatic transient in the galaxy M31. A possible origin of M85 OT2006-1 is a stellar merger. If so, searches for similar events in nearby galaxies will not only allow to study the physics of hyper-Eddington sources, but also probe an important phase in the evolution of stellar binary systems.

We have compiled the largest sample of multiwavelength spectral energy distributions (SEDs) of Broad Absorption Line (BAL) quasars to date, from the radio to the X-ray. We present new Spitzer MIPS (24, 70, and 160 micron) observations of 38 BAL quasars in addition to data from the literature and public archives. In general, the mid-infrared properties of BAL quasars are consistent with those of non-BAL quasars of comparable luminosity. In particular, the optical-to-mid-infrared luminosity ratios of the two populations are indistinguishable. We also measure or place upper limits on the contribution of star formation to the far-infrared power. Of 22 (57%) upper limits, seven quasars have sufficiently sensitive constraints to conclude that star formation likely contributes little (<20%) to their far-infrared power. The 17 BAL quasars (45%) with detected excess far-infrared emission likely host hyperluminous starbursts with L_fir,SF=10^{13-14} L_sun. Mid-infrared through X-ray composite BAL quasar SEDs are presented, incorporating all of the available photometry. Overall, we find no compelling evidence for inherent differences between the SEDs of BAL vs. non-BAL quasars of comparable luminosity. Therefore a ``cocoon'' picture of a typical BAL quasar outflow whereby the wind covers a large fraction of the sky is not supported by the mid-infrared SED comparison with normal quasars, and the disk-wind paradigm with a typical radio-quiet quasar hosting a BAL region remains viable.

The absolute calibration and characterization of the Multiband Imaging Photometer for Spitzer (MIPS) 70 micron coarse- and fine-scale imaging modes are presented based on over 2.5 years of observations. Accurate photometry (especially for faint sources) requires two simple processing steps beyond the standard data reduction to remove long-term detector transients. Point spread function (PSF) fitting photometry is found to give more accurate flux densities than aperture photometry. Based on the PSF fitting photometry, the calibration factor shows no strong trend with flux density, background, spectral type, exposure time, or time since anneals. The coarse-scale calibration sample includes observations of stars with flux densities from 22 mJy to 17 Jy, on backgrounds from 4 to 26 MJy sr^-1, and with spectral types from B to M. The coarse-scale calibration is 702 +/- 35 MJy sr^-1 MIPS70^-1 (5% uncertainty) and is based on measurements of 66 stars. The instrumental units of the MIPS 70 micron coarse- and fine-scale imaging modes are called MIPS70 and MIPS70F, respectively. The photometric repeatability is calculated to be 4.5% from two stars measured during every MIPS campaign and includes variations on all time scales probed. The preliminary fine-scale calibration factor is 2894 +/- 294 MJy sr^-1 MIPS70F^-1 (10% uncertainty) based on 10 stars. The uncertainty in the coarse- and fine-scale calibration factors are dominated by the 4.5% photometric repeatability and the small sample size, respectively. The 5-sigma, 500 s sensitivity of the coarse-scale observations is 6-8 mJy. This work shows that the MIPS 70 micron array produces accurate, well calibrated photometry and validates the MIPS 70 micron operating strategy, especially the use of frequent stimulator flashes to track the changing responsivities of the Ge:Ga detectors.

We determine the nature, redshift, stellar populations and dust properties of optically faint or non-detected extremely red objects (ERO) found from our survey of the lensing clusters A1835 and AC114. We perform SED fitting and use deep optical, HST, VLT, Spitzer data, and for some objects and sub-mm data. For most of the lensed EROs we find photometric redshifts showing a strong degeneracy between "low-z" (z~1-3) and high-z (z~6-7). Although formally best fits are often found at high-z, their resulting bright absolute magnitudes, the number density of these objects, and in some cases Spitzer photometry or longer wavelength observations, suggest strongly that all of these objects are at "low-z". The majority of these objects are best fitted with young (<~ 0.5-0.7 Gyr) and dusty starbursts. Indications for strong extinction, with A_V~2.4-4, are found in some objects. For dusty objects star formation rates (SFR) have been estimated from the bolometric luminosity determined after fitting of semi-empirical starburst, ERO, and ULIRG templates. Typically we find SFR~(1-18) Msun/yr. Again, SMMJ14009+0252 stands out as a LIRG with SFR 1000 Msun/yr. Concerning the comparison objects, we argue that the massive post-starburst z~6.5 galaxy candidate HUDF-J2 showing observed properties very similar to our EROs, is more likely a dusty starburst at z~2.3-2.6. This interpretation also naturally explains the observed 24 micron emission. Both empirically and from our SED fits we find that the IRAC selectec EROs from Yan et al (2004) show very similar properties to our lensed EROs. Reasonable fits are found for most of them with relatively young and dusty stellar populations.

We describe two peculiar galaxies falling into the massive galaxy clusters Abell 1689 (z~0.18) and 2667 (z~0.23) respectively. Hubble Space Telescope images show extraordinary trails composed of bright blue knots (-16.5<M<-11.5 mag) and stellar streams associated with each of these systems. Combining optical, near and mid-infrared and radio observations we prove that while both galaxies show similar extended trails of star-forming knots, their recent star formation histories are different. One (~L*) is experiencing a strong burst of star formation, appearing as a rare example of a luminous infrared cluster galaxy. In comparison, the other (~ 0.1 L*) has recently ceased its star formation activity. Our model suggests that the morphologies and star formation in these galaxies have been influenced by the combined action of tidal interaction (likely with the cluster potential) and of ram pressure with the intracluster medium. These results can be used to gain more insights to the origin of S0s, dwarf and ultra-compact dwarf (UCD) cluster galaxies.

We summarise the current status of our project to identify and study z~6-10 galaxies thanks to strong gravitational lensing. Building on the detailed work from Richard et al. (2006), we present results from new follow-up observations (imaging) undertaken with ACS/HST and the Spitzer Space Telescope and compare our results with findings from the Hubble Ultra-Deep Field (UDF). These new observations are in agreement with the high-z nature for the vast majority of the candidates presented in Richard et al. (2006). We also discuss the properties of other optical dropout sources found in our searches and related objects (EROs, sub-mm galaxies,...) from other surveys.

We have measured the contribution of submillimeter and mid-infrared sources to the extragalactic background radiation at 70 and 160um. Specifically, we have stacked flux in 70 and 160um Spitzer Space Telescope (Spitzer) observations of the Canada-UK Deep Sub-millimeter Survey 14h field at the positions of 850um sources detected by SCUBA and also 8 and 24um sources detected by Spitzer. We find that per source, the SCUBA galaxies are the strongest and the 8um sources the weakest contributors to the background flux at both 70 and 160um. Our estimate of the contribution of the SCUBA sources is higher than previous estimates. However, expressed as a total contribution, the full 8um source catalogue accounts for twice the total 24um source contribution and ~10 times the total SCUBA source contribution. The 8um sources account for the majority of the background radiation at 160um with a flux of 0.87+/-0.16 MJy/sr and at least a third at 70um with a flux of 0.103+/-0.019 MJy/sr. These measurements are consistent with current lower limits on the background at 70 and 160um. Finally, we have investigated the 70 and 160um emission from the 8 and 24um sources as a function of redshift. We find that the average 70um flux per 24um source and the average 160um flux per 8 and 24um source is constant over all redshifts, up to z~4. In contrast, the low-redshift half (z<1) of the of 8um sample contributes approximately four times the total 70um flux of the high-redshift half. These trends can be explained by a single non-evolving SED.

We present 0.5 -160 micron Spectral Energy Distributions (SEDs) of galaxies, detected at 70microns with the Multiband Imaging Photometer for Spitzer (MIPS), using broadband imaging data from Spitzer and ground-based telescopes. Spectroscopic redshifts, in the range 0.2<z<1.5, have been measured as part of the Deep Extragalactic Evolutionary Probe2 (DEEP2) project. Based on the SEDs we explore the nature and physical properties of the sources. Using the optical spectra we derive Hbeta and [OII]-based Star Formation Rates (SFR) which are 10-100 times lower than SFR estimates based on IR and radio. The median offset in SFR between optical and IR is reduced by a factor of ~3 when we apply a typical extinction corrections. We investigate mid-to-far infrared correlations for low redshift (>0.5) and high redshift (0.5<z<1.2) bins. Using this unique ``far-infrared'' selected sample we derive an empirical mid to far-infrared relationship that can be used to estimate the infrared energy budget of galaxies in the high-redshift universe. Our sample can be used as a template to translate far-infrared luminosities into bolometric luminosities for high redshift objects.

We discuss images of M87 from 3.6 to 160 um obtained with Spitzer. As found previously, there is an excess in the far infrared over a simple power law interpolation from the radio to the resolved nonthermal features in the mid-infrared and optical. We show that this excess is most likely warm dust in the galaxy itself, and that the properties of this emission component are similar to the far infrared emission of normal giant elliptical galaxies. The new observations let us determine the spectrum of the jet and surrounding lobes of nonthermal emission. We find that even in the lobes the synchrotron break frequency is in the optical, probably requiring in situ particle acceleration not only in the jet but in the lobes as well.

We present the Spitzer/IRS mid-infrared spectrum of the infrared-luminous (L_{IR}=4e11 L_sun) brightest cluster galaxy (BCG) in the X-ray-luminous cluster Z3146 (z=0.29). The spectrum shows strong aromatic emission features, indicating that the dominant source of the infrared luminosity is star formation. The most striking feature of the spectrum, however, is the exceptionally strong molecular hydrogen (H2) emission lines, which seem to be shock-excited. The line luminosities and inferred warm H2 gas mass (~1e10 M_sun) are 6 times larger than those of NGC 6240, the most H2-luminous galaxy at z <~ 0.1. Together with the large amount of cold H2 detected previously (~1e11 M_sun), this indicates that the Z3146 BCG contains disproportionately large amounts of both warm and cold H2 gas for its infrared luminosity, which may be related to the intracluster gas cooling process in the cluster core.

Observing the first galaxies formed during the reionisation epoch, i.e. approximately within the first billion years after the Big Bang, remains one of the challenges of contemporary astrophysics. Several efforts are being undertaken to search for such remote objects. Combining the near-IR imaging power of the VLT and the natural effect of strong gravitational lensing our pilot program has allowed us to identify several galaxy candidates at redshift 6 <~ z <~ 10. The properties of these objects and the resulting constraints on the star formation rate density at high redshift are discussed. Finally we present the status of follow-up observations (ISAAC spectroscopy, HST and Spitzer imaging) and discuss future developments.

As part of the "All Wavelength Extended Groth Strip International Survey" (AEGIS), we describe the panchromatic characterization of an X-ray luminous active galactic nucleus (AGN) in a merging galaxy at z=1.15. This object is detected at infrared (8mic, 24mic, 70mic, 160mic), submillimeter (850mic) and radio wavelengths, from which we derive a bolometric luminosity L_bol ~ 9x10^12 Lsol. We find that the AGN clearly dominates the hot dust emission below 40mic but its total energetic power inferred from the hard X-rays is substantially less than the bolometric output of the system. About 50% of the infrared luminosity is indeed produced by a cold dust component that probably originates from enshrouded star formation in the host galaxy. In the context of a coeval growth of stellar bulges and massive black holes, this source might represent a ``transition'' object sharing properties with both quasars and luminous starbursts. Study of such composite galaxies will help address how the star formation and disk-accretion phenomena may have regulated each other at high redshift and how this coordination may have participated to the build-up of the relationship observed locally between the masses of black holes and stellar spheroids.

We report the detection of rest--frame 6.2 and 7.7 \micron emission features arising from Polycyclic Aromatic Hydrocarbons (PAH) in the Spitzer/IRS spectrum of an infrared-luminous Lyman break galaxy at z=3.01. This is currently the highest redshift galaxy where these PAH emission features have been detected. The total infrared luminosity inferred from the MIPS 24 \micron and radio flux density is 2$\times10^{13}$ L$_{\odot}$, which qualifies this object as a so--called hyperluminous infrared galaxy (HyLIRG). However, unlike local HyLIRGs which are generally associated with QSO/AGNs and have weak or absent PAH emission features, this HyLIRG has very strong 6.2 and 7.7 \micron PAH emission. We argue that intense star formation dominates the infrared emission of this source, although we cannot rule out the presence of a deeply obscured AGN. This LBG appears to be a distorted system in the HST ACS F606W and F814W images, possibly indicating that a significant merger or interaction is driving the large IR luminosity.

A new 6-cm survey of almost 0.6 square degrees to a limit of 0.55-mJy/beam (10-sigma) finds 37 isolated radio sources and 7 radio source pairs (not necessarily physical companions). IRAC counterparts are identified for at least 92% of the radio sources within the area of deep IRAC coverage, which includes 31 isolated sources and 6 pairs. This contrasts with an identification rate of <74% to R<23.95 in visible light. Eight of the IRAC galaxies have power law spectral energy distributions, implying that the mid-infrared emission comes from a powerful AGN. The remaining 26 IRAC galaxies show stellar emission in the mid-infrared, probably in most of these galaxies because the stellar emission is bright enough to outshine an underlying AGN. The infrared colors suggest that the majority of these galaxies are bulge-dominated and have redshifts between approximately 0.5 and 1.0. Visible spectra from the DEEP2 redshift survey, available for 11 galaxies, are consistent with this suggestion. The IRAC galaxies fall into two distinct groups in a color-magnitude diagram, one group (the "stripe") includes all the AGN. The other group (the "blue clump") has blue 3.6 to 8 micron colors and a small range of 8 micron magnitudes. This separation should be useful in classifying galaxies found in other radio surveys.

We describe the spectral energy distributions (SEDs) of Lyman Break Galaxies (LBGs) at z~3 using deep mid-infrared and optical observations of the Extended Groth Strip, obtained with IRAC and MIPS on board Spitzer and from the ground, respectively. We focus on LBGs with detections at all four IRAC bands, in particular the 26 galaxies with IRAC 8 micron band (rest--frame K-band) detections. We use stellar population synthesis models and probe the stellar content of these galaxies. Based on best--fit continuous star-formation models we derive estimates of the stellar mass for these LBGs. As in previous studies, we find that a fraction of LBGs have very red colors and large estimated stellar masses (M_stellar > 5 x 10 ^{10} M_solar): the present Spitzer data allow us, for the first time, to study these massive LBGs in detail. We discuss the link between these LBGs and submm-luminous galaxies.We find that the number density of these massive LBGs at high redshift is higher than predicted by current semi-analytic models of galaxy evolution.

We identify 17 possible 8.0 micron-selected counterparts to the submillimeter galaxies in the CUDSS 14-hour field, derived from deep imaging carried out with the IRAC and MIPS instruments aboard the Spitzer Space Telescope. Ten of the 17 counterparts are not the same as those previously identified at shorter wavelengths. We argue that 8.0 micron selection offers a better means for identifying counterparts to submillimeter galaxies than near-infrared or optical selection. Based on the panchromatic SEDs, most counterparts appear to be powered by ongoing star formation. Power-law fits to the SEDs suggest that five objects in the 8.0 micron-selected counterpart sample harbor dominant AGNs; a sixth object is identified as a possible AGN. The 3.6 to 8.0 micron colors of the infrared-selected counterparts are significantly redder than the general IRAC galaxy population in the CUDSS 14-hour field.

We present the identification of a bright submillimeter galaxy (SMG) in the field of Abell 2218. The galaxy has a spectroscopic redshift of ~4, and is currently the highest redshift SMG known. It is detected at all wavelengths from optical to submm, including the Spitzer IRAC bands. We discuss the properties of this galaxy, which is undergoing intense star formation at a rate~600 Msun/yr. We also compare the properties to those of radio-preselected submm-bright galaxies. The z~4 result extends the redshift distribution of SMGs.

We have studied the infrared properties of the brightest cluster galaxies (BCGs) located in the cores of X-ray-luminous clusters at 0.15 < z < 0.35. The majority of the BCGs are not particularly infrared-luminous compared with other massive early-type galaxies, suggesting that the cluster environment has little influence on the infrared luminosities of the BCGs. The exceptions, however, are the BCGs in the three X-ray-brightest clusters in the sample, A1835, Z3146, and A2390. These BCGs have a prominent far-infrared peak in their spectral energy distributions (SEDs), and two of them (those in A1835 and Z3146) can be classified as luminous infrared galaxies (LIRGs: L_{IR} > 10^{11} L_{sun}). Although radio AGNs are found to be prevalent among the BCGs, the infrared luminosities of these three BCGs, judged from the infrared SED signatures, are likely to be powered by star formation. Considering the overall trend that clusters with shorter radiative gas cooling times harbor more infrared-luminous BCGs, the enhanced star formation may be caused by the cooling cluster gas accreting onto the BCGs.

We are currently undertaking a Spitzer GTO program to image ~30 massive lensing clusters at moderate redshift with both IRAC and MIPS. By taking advantage of the gravitatinoal lensing power of these clusters, we will study the population of faint galaxies that are below the nominal Spitzer detection limits. Here, we present a few examples of our science programs.

We present 350 micron photometry of all 17 galaxy candidates in the Lockman Hole detected in a 1.1 mm Bolocam survey. Several of the galaxies were previously detected at 850 microns, at 1.2 mm, in the infrared by Spitzer, and in the radio. Nine of the Bolocam galaxy candidates were detected at 350 microns and two new candidates were serendipitously detected at 350 microns (bringing the total in the literature detected in this way to three). Five of the galaxies have published spectroscopic redshifts, enabling investigation of the implied temperature ranges and a comparison of photometric redshift techniques. Lambda = 350 microns lies near the spectral energy distribution peak for z = 2.5 thermally emitting galaxies. Thus, luminosities can be measured without extrapolating to the peak from detection wavelengths of lambda > 850 microns. Characteristically, the galaxy luminosities lie in the range 1.0 - 1.2 x 10^13 L_solar, with dust temperatures in the range of 40 K to 70 K, depending on the choice of spectral index and wavelength of unit optical depth. The implied dust masses are 3 - 5 x 10^8 M_solar. We find that the far-infrared to radio relation for star-forming ULIRGs systematically overpredicts the radio luminosities and overestimates redshifts on the order of Delta z ~ 1, whereas redshifts based on either on submillimeter data alone or the 1.6 micron stellar bump and PAH features are more accurate.

Gravitational lensing by massive galaxy clusters is a powerful tool for the discovery and study of high redshift galaxies, including those at z>=6 likely responsible for cosmic re-ionization. Pello et al. recently used this technique to discover a candidate gravitationally magnified galaxy at z=10 behind the massive cluster lens Abell 1835 (z=0.25). We present new Keck (LRIS) and Spitzer Space Telescope (IRAC) observations of the z=10 candidate (hereafter #1916) together with a re-analysis of archival optical and near-infrared imaging from the Hubble Space Telescope and VLT respectively. Our analysis therefore extends from the atmospheric cut-off at lambda_obs=0.35um out to lambda_obs=5um with Spitzer/IRAC. The z=10 galaxy is not detected in any of these data, including an independent reduction of Pello et al.'s discovery H- and K-band imaging. We conclude that there is no statistically reliable evidence for the existence of #1916. We also assess the implications of our results for ground-based near-infrared searches for gravitationally magnified galaxies at z>~7. The broad conclusion is that such experiments remain feasible, assuming that space-based optical and mid-infrared imaging are available to break the degeneracy with low redshift interlopers (e.g. z~2-3) when fitting spectral templates to the photometric data.

Mid-infrared observations of Active Galactic Nuclei (AGN) are important for understanding of the physical conditions around the central accretion engines. Chandra and XMM-Newton X-ray observations of a 300 arcmin^2 region in the Extended Groth Strip are used to select a sample of~150 AGN. The Spitzer instruments IRAC and MIPS detect 68-80% of these sources, which show a wide range of mid-infrared properties. About 40% of the sources have red power-law spectral energy distributions (f_nu ~ nu^alpha, alpha<0) in the 3.6-8 um IRAC bands. In these sources the central engine dominates the emission at both X-ray and IR wavelengths. Another 40% of the sources have blue mid-IR spectral energy distributions (alpha>0) with their infrared emission dominated by the host galaxy; the remaining 20% are not well-fit by a power law. Published IRAC color criteria for AGN select most of the red sources, but only some of the blue sources. As with all other known methods, selecting AGN with mid-IR colors will not produce a sample that is simultaneously complete and reliable. The IRAC SED type does not directly correspond to X-ray spectral type (hard/soft). The mid-IR properties of X-ray-detected Lyman-break, radio, submillimeter, and optically-faint sources vary widely and, for the most part, are not distinct from those of the general X-ray/infrared source population. X-ray sources emit 6-11% of the integrated mid--IR light, making them significant contributors to the cosmic infrared background.

We investigate the redshift dependence of the contribution to the extragalactic far-infrared/sub-millimeter background from galaxies detected by the Spitzer Space Telescope at 8um and 24um. Using seven-band optical to mid-infrared photometry, we estimate photometric redshifts for the Spitzer sources which appear to be mostly L* galaxies at a median redshift of z=1.0. These sources, extracted from deep 8um and 24um mosaics of the CUDSS 14-hour field with 5sigma limits of 5.8uJy and 70uJy respectively, exhibit significant 850um and 450um emission as observed by SCUBA. At 850um, after removing >=4sigma sources and those securely identified in our previous cross-matching paper, we measure stacked flux at the significance level of 4.4sigma and 2.9sigma from the full 8um and 24um galaxy catalogue respectively. At 450um, flux is detected from all 8um galaxies at the level of 3.5sigma, while there is no significant emission from the 24um galaxies. We find that the 850um flux is emitted almost exclusively at z>~1.3 from the Spitzer sources with 0.44mJy (4.7sigma) per 8um source and 0.51mJy (2.8sigma) per 24um source. This corresponds to a contribution of (16+/-3)% toward the 850um extra-galactic background from the 8um sources and (5.0+/-1.8)% from the 24um sources. At 450um, only the 8um sources within the redshift interval 1<z<2 exhibit significant emission with an average flux per source of 3.35mJy (3.0sigma). This is a contribution of (37+/-12)% to the 450um background.

We present Spitzer 24micron imaging of 1.5 < z < 2.5 Distant Red Galaxies (DRGs) in the 10arcmin by 10arcmin Extended Hubble Deep Field South of the Multiwavelength Survey by Yale-Chile. We detect 65% of the DRGs with K_AB < 23.2 mag at S_24micron > 40uJy, and conclude that the bulk of the DRG population are dusty active galaxies. A mid-infrared (MIR) color analysis with IRAC data suggests that the MIR fluxes are not dominated by buried AGN, and we interpret the high detection rate as evidence for a high average star formation rate of <SFR> = 130+/-30 M/yr. From this, we infer that DRGs are important contributors to the cosmic star formation rate density at z ~ 2, at a level of \~0.02 M/yr/Mpc^3 to our completeness limit of K_AB = 22.9 mag

We investigate the nature of a sample of 92 Spitzer/MIPS 24 micron selected galaxies in the CDFS, showing power law-like emission in the Spitzer/IRAC 3.6-8 micron bands. The main goal is to determine whether the galaxies not detected in X-rays (47% of the sample) are part of the hypothetical population of obscured AGN not detected even in deep X-ray surveys. The majority of the IR power-law galaxies are ULIRGs at z>1, and those with LIRG-like IR luminosities are usually detected in X-rays. The optical to IR spectral energy distributions (SEDs) of the X-ray detected galaxies are almost equally divided between a BLAGN SED class (similar to an optically selected QSO) and a NLAGN SED (similar to the BLAGN SED but with an obscured UV/optical continuum). A small fraction of SEDs resemble warm ULIRG galaxies (e.g., Mrk231). Most galaxies not detected in X-rays have SEDs in the NLAGN+ULIRG class as they tend to be optically fainter, and possibly more obscured. Moreover, the IR power-law galaxies have SEDs significantly different from those of high-z (z_sp>1) IR (24 micron) selected and optically bright (VVDS I_AB<=24) star-forming galaxies whose SEDs show a very prominent stellar bump at 1.6 micron. The galaxies detected in X-rays have 2-8 keV rest-frame luminosities typical of AGN. The galaxies not detected in X-rays have global X-ray to mid-IR SED properties that make them good candidates to contain IR bright X-ray absorbed AGN. If all these sources are actually obscured AGN, we would observe a ratio of obscured to unobscured 24 micron detected AGN of 2:1, whereas models predict a ratio of up to 3:1. Additional studies using Spitzer to detect X-ray-quiet AGN are likely to find more such obscured sources.

We investigate the properties of massive galaxies at z=1-3.5 using HST observations, ground-based near-IR imaging, and Spitzer Space Telescope observations at 3-24 micron. We identify 153 distant red galaxies (DRGs) with J-K > 2.3 mag (Vega) in the southern GOODS field. This sample is approximately complete in stellar mass for passively evolving galaxies above 10^11 solar masses and z < 3. The galaxies identified by this selection are roughly split between objects whose optical and near-IR rest-frame light is dominated by evolved stars combined with ongoing star formation, and galaxies whose light is dominated by heavily reddened starbursts. Very few of the galaxies (< 10%) have no indication of current star formation. Using SFR estimates that include the reradiated IR emission, the DRGs at z=1.5-3 with stellar masses > 10^11 solar masses have specific SFRs (SFRs per unit stellar mass) ranging from 0.2 to 10 Gyr^-1, with a mean value of ~2.4 Gyr^-1. The DRGs with stellar masses > 10^11 solar masses and 1.5 < z < 3 have integrated specific SFRs greater the global value over all galaxies. In contrast, we find that galaxies at z = 0.3-0.75 with these stellar masses have integrated specific SFRs less than the global value, and more than an order of magnitude lower than that for massive DRGs at z = 1.5-3. At z < 1, lower-mass galaxies dominate the overall cosmic mass assembly. This suggests that the bulk of star formation in massive galaxies occurs at early cosmic epochs and is largely complete by z~1.5. [Abridged]

In order to examine the relative importance of powerful starbursts and Compton-thick AGNs in NGC 6240, we have obtained mid-infrared images and low-resolution spectra of the galaxy with sub-arcsecond spatial resolution using the Keck Telescopes. Despite the high spatial resolution (~200 pc) of our data, no signature of the hidden AGNs has been detected in the mid-infrared. The southern nucleus, which we show provides 80-90% of the total 8-25 um luminosity of the system, has a mid-infrared spectrum and a mid-/far-infrared spectral energy distribution consistent with starbursts. At the same time, however, it is also possible to attribute up to 60% of the bolometric luminosity to an AGN, consistent with X-ray observations, if the AGN is heavily obscured and emits mostly in the far-infrared. This ambiguity arises because the intrinsic variation of properties among a given galaxy population (e.g., starbursts) introduces at least a factor of a few uncertainty even into the most robust AGN-starburst diagnostics. We conclude that with present observations it is not possible to determine the dominant power source in galaxies when AGN and starburst luminosities are within a factor of a few of each other.

A deep mid- and far-infrared survey in the Extended Groth Strip (EGS) area gives 3.6 to 8\micron flux densities or upper limits for 253 Lyman Break Galaxies (LBGs). The LBGs are a diverse population but with properties correlated with luminosity. The LBGs show a factor of 30 range in indicated stellar mass and a factor of 10 in apparent dust content relative to stellar mass. About 5% of LBGs are luminous at all wavelengths with powerful emission at rest 6\micron. In the rest 0.9 to 2\micron spectral range these galaxies have stellar spectral slopes with no sign of an AGN power law component, suggesting that their emission is mainly powered by intensive star formation. Galaxies in this luminous population share the infrared properties of cold SCUBA sources: both are massive and dusty starburst galaxies at $2<z<3$; their stellar mass is larger than $10^{11} M_{\odot}$. We suggest that these galaxies are the progenitors of present-day giant elliptical galaxies, with a substantial fraction of their stars already formed at $z \approx 3$.

We analyze a sample of ~2600 MIPS/Spitzer 24mic sources brighter than ~80muJy and located in the Chandra Deep Field South to characterize the evolution of the comoving infrared (IR) energy density of the Universe up to z~1. Using published ancillary optical data we first obtain a nearly complete redshift determination for the 24mic objects associated with R<24 counterparts at z<1. We find that the 24mic population at 0.5<z<1 is dominated by ``Luminous Infrared Galaxies'' (i.e., 10^11 L_sol < L_IR < 10^12 L_sol), the counterparts of which appear to be also luminous at optical wavelengths and tend to be more massive than the majority of optically-selected galaxies. We finally derive 15mic and total IR luminosity functions (LFs) up to z~1. In agreement with the previous results from ISO and SCUBA and as expected from the MIPS source number counts, we find very strong evolution of the contribution of the IR-selected population with lookback time. Pure evolution in density is firmly excluded by the data, but we find considerable degeneracy between strict evolution in luminosity and a combination of increases in both density and luminosity (L*_IR prop. to (1+z)^{3.2_{-0.2}^{+0.7}}, Phi*_IR prop. to (1+z)^{0.7_{-0.6}^{+0.2}}). Our results imply that the comoving IR energy density of the Universe evolves as (1+z)^(3.9+/-0.4) up to z~1 and that galaxies luminous in the infrared (i.e., L_IR > 10^11 L_IR) are responsible for 70+/-15% of this energy density at z~1. Taking into account the contribution of the UV luminosity evolving as (1+z)^~2.5, we infer that these IR-luminous sources dominate the star-forming activity beyond z~0.7. The uncertainties affecting these conclusions are largely dominated by the errors in the k-corrections used to convert 24mic fluxes into luminosities.

Two images of Cassiopeia A obtained at 24 micrometer with the Spitzer Space Telescope over a one year time interval show moving structures outside the shell of the supernova remnant to a distance of more than 20 arcmin. Individual features exhibit apparent motions of 10 to 20 arcsec per year, independently confirmed by near-infrared observations. The observed tangential velocities are at roughly the speed of light. It is likely that the moving structures are infrared echoes, in which interstellar dust is heated by the explosion and by flares from the compact object near the center of the remnant.

We use a 24 micron selected sample containing more than 8,000 sources to study the evolution of star-forming galaxies in the redshift range from z=0 to z~3. We obtain photometric redshifts for most of the sources in our survey using a method based on empirically-built templates spanning from ultraviolet to mid-infrared wavelengths. The accuracy of these redshifts is better than 10% for 80% of the sample. The derived redshift distribution of the sources detected by our survey peaks at around z=0.6-1.0 (the location of the peak being affected by cosmic variance), and decays monotonically from z~1 to z~3. We have fitted infrared luminosity functions in several redshift bins in the range 0<z<~3. Our results constrain the density and/or luminosity evolution of infrared-bright star-forming galaxies. The typical infrared luminosity (L*) decreases by an order of magnitude from z~2 to the present. The cosmic star formation rate (SFR) density goes as (1+z)^{4.0\pm0.2} from z=0 to z=0.8. From z=0.8 to z~1.2, the SFR density continues rising with a smaller slope. At 1.2<z<3, the cosmic SFR density remains roughly constant. The SFR density is dominated at low redshift (z<0.5) by galaxies which are not very luminous in the infrared (L_TIR<1.e11 L_sun, where L_TIR is the total infrared luminosity, integrated from 8 to 1000 micron). The contribution from luminous and ultraluminous infrared galaxies (L_TIR>1.e11 L_sun) to the total SFR density increases steadily from z~0 up to z~2.5, forming at least half of the newly-born stars by z~1.5. Ultraluminous infrared galaxies (L_TIR>1.e12 L_sun) play a rapidly increasing role for z>~1.3.

We report MIPS photometry of 20 radio-loud quasars and galaxies at 24 and 70 um (and of five at 160 um). We combine this sample with additional sources detected in the far infrared by IRAS and ISO for a total of 47 objects, including 23 steep spectrum Type I AGNs: radio-loud quasars and broad line radio galaxies; and 24 Type II AGNs: narrow line and weak line FR II radio galaxies. Of this sample, the far infrared emission of all but 3C 380 appears to be dominated by emission by dust heated by the AGN and by star formation. The AGN appears to contribute more than 50% of the far infrared luminosity in most of sources. It is also expected that the material around the nucleus is optically thin in the far infrared. Thus, the measurements at these wavelengths can be used to test the orientation-dependent unification model. As predicted by the model, the behavior of the sources is consistent with the presence of an obscuring circumnuclear torus; in fact, we find it may still have significant optical depth at 24 um.In addition, as expected for the radio-loud quasars, there is a significant correlation between the low frequency radio (178 MHz) and the 70 um emission, two presumably isotropic indicators of nuclear activity. This result is consistent with the simple unified scheme. However, there is a population of radio galaxies that are underluminous at 70 um compared with the radio-loud quasars and hence are a challenge to the simple unified model.

We present a first analysis of deep 24 micron observations with the Spitzer Space Telescope of a sample of nearly 1500 galaxies in a thin redshift slice, 0.65<z<0.75. We combine the infrared data with redshifts, rest-frame luminosities, and colors from COMBO-17, and with morphologies from Hubble Space Telescope images collected by the GEMS and GOODS projects. To characterize the decline in star-formation rate (SFR) since z~0.7, we estimate the total thermal infrared (IR) luminosities, SFRs, and stellar masses for the galaxies in this sample. At z~0.7, nearly 40% of intermediate and high-mass galaxies (with stellar masses >2x10^10 solar masses) are undergoing a period of intense star formation above their past-averaged SFR. In contrast, less than 1% of equally-massive galaxies in the local universe have similarly intense star formation activity. Morphologically-undisturbed galaxies dominate the total infrared luminosity density and SFR density: at z~0.7, more than half of the intensely star-forming galaxies have spiral morphologies, whereas less than \~30% are strongly interacting. Thus, a decline in major-merger rate is not the underlying cause of the rapid decline in cosmic SFR since z~0.7. Physical properties that do not strongly affect galaxy morphology - for example, gas consumption and weak interactions with small satellite galaxies - appear to be responsible.

We describe the data reduction algorithms for the Multiband Imaging Photometer for Spitzer (MIPS) instrument. These algorithms were based on extensive preflight testing and modeling of the Si:As (24 micron) and Ge:Ga (70 and 160 micron) arrays in MIPS and have been refined based on initial flight data. The behaviors we describe are typical of state-of-the-art infrared focal planes operated in the low backgrounds of space. The Ge arrays are bulk photoconductors and therefore show a variety of artifacts that must be removed to calibrate the data. The Si array, while better behaved than the Ge arrays, does show a handful of artifacts that also must be removed to calibrate the data. The data reduction to remove these effects is divided into three parts. The first part converts the non-destructively read data ramps into slopes while removing artifacts with time constants of the order of the exposure time. The second part calibrates the slope measurements while removing artifacts with time constants longer than the exposure time. The third part uses the redundancy inherit in the MIPS observing modes to improve the artifact removal iteratively. For each of these steps, we illustrate the relevant laboratory experiments or theoretical arguments along with the mathematical approaches taken to calibrate the data. Finally, we describe how these preflight algorithms have performed on actual flight data.

The first six months of flight data from the Multiband Imaging Photometer for Spitzer (MIPS) were used to test MIPS reduction algorithms based on extensive preflight laboratory data and modeling. The underlying approach for the preflight algorithms has been found to be sound, but some modifications have improved the performance. The main changes are scan mirror dependent flat fields at 24 micron, hand processing to remove the time dependent stim flash latents and fast/slow response variations at 70 micron, and the use of asteroids and other sources instead of stars for flux calibration at 160 micron due to a blue ``leak.'' The photometric accuracy of flux measurpements is currently 5%, 10%, and 20% at 24, 70, and 160 micron, respectively. These numbers are expected to improve as more flight data are analyzed and data reduction algorithms refined.

Using deep observations of the Chandra Deep Field South obtained with MIPS at 24mic, we present our preliminary estimates on the evolution of the infrared (IR) luminosity density of the Universe from z=0 to z~1. We find that a pure density evolution of the IR luminosity function is clearly excluded by the data. The characteristic luminosity L_IR* evolves at least by (1+z)^3.5 with lookback time, but our monochromatic approach does not allow us to break the degeneracy between a pure evolution in luminosity or an evolution in both density and luminosity. Our results imply that IR luminous systems (L_IR > 10^11 L_sol) become the dominant population contributing to the comoving IR energy density beyond z~0.5-0.6. The uncertainties affecting our measurements are largely dominated by the poor constraints on the spectral energy distributions that are used to translate the observed 24mic flux into luminosities.

We report the detection of a z~7 galaxy strongly lensed by the massive galaxy cluster Abell 2218 (z=0.175) at 3.6 and 4.5 um using the Spitzer Observatory and at 1.1 um using the Hubble Space Telescope. The new data indicate a refined photometric redshift in the range of 6.6-6.8 depending on the presence of Ly-alpha emission. The spectral energy distribution is consistent with having a significant Balmer break, suggesting that the galaxy is in the poststarburst stage with an age of at least ~50 Myr and quite possibly a few hundred Myr. This suggests the possibility that a mature stellar population is already in place at such a high redshift. Compared with typical Lyman break galaxies at z~3-4, the stellar mass is an order of magnitude smaller (~10^{9} Msun), but the specific star formation rate (star formation rate/M_{star}) is similarly large (> 10^{-9} yr^{-1}), indicating equally vigorous star-forming activity.

New surveys with the Spitzer space telescope identify distant star-forming and active galaxies by their strong emission at far-infrared wavelengths, which provides strong constraints on these galaxies' bolometric energy. Using early results from Spitzer surveys at 24 micron, we argue that the faint sources correspond to the existence of a population of infrared-luminous galaxies at z > 1 that are not expected from predictions based on previous observations from ISO and IRAS. Combining Spitzer images with deep ground-based optical and Hubble Space Telescope imaging, we discuss the properties of galaxies selected at 24 micron in the region of the Chandra Deep Field South, including redshift and morphological distributions. Galaxies with z < 1 constitute roughly half of the faint 24 micron sources. Infrared-luminous galaxies at these redshifts span a wide variety of normal to strongly interacting/merging morphologies, which suggests that a range of mechanisms produce infrared activity. Large-area, joint surveys between Spitzer and HST are needed to understand the complex relation between galaxy morphology, structure, environment and activity level, and how this evolves with cosmic time. We briefly discuss strategies for constructing surveys to maximize the legacy of these missions.

We derive galaxy source counts at 70 and 160 microns using the Multiband Imaging Photometer for Spitzer (MIPS) to map the Chandra Deep Field South (CDFS) and other fields. At 70 microns, our observations extend upwards about 2 orders of magnitude in flux density from a threshold of 15 mJy, and at 160 microns they extend about an order of magnitude upward from 50 mJy. The counts are consistent with previous observations on the bright end. Significant evolution is detected at the faint end of the counts in both bands, by factors of 2-3 over no-evolution models. This evolution agrees well with models that indicate most ofthe faint galaxies lie at redshifts between 0.7 and 0.9. The new Spitzer data already resolve about 23% of the Cosmic Far Infrared Background at 70 microns and about 7% at 160 microns.

We present Spitzer IRAC (3.6, 4.5, 5.8, 8.0 um) and MIPS (24 um) observations of the SCUBA submillimeter sources and uJy VLA radio sources in a 5'x5' area in the Lockman Hole East region. Out of the ~40 SCUBA/VLA sources in the field, Spitzer counterparts were detected for nearly all except for the few low-weight SCUBA detections. We show that the majority (80-90%) of the detected sources are cold (i.e., starburst-like) infrared-luminous galaxies (L_{IR} > 10^{11} Lsun) at redshift 0.5<z<3.5, whose star-formation rate density (SFRD) is comparable to that of the optically-selected star-forming galaxies.

We present images of the galaxy NGC 55 at 24, 70, and 160 micron obtained with the Multiband Imaging Photometer for Spitzer (MIPS) instrument aboard the Spitzer Space Telescope. The new images display the far infrared emission in unprecedented detail and demonstrate that the infrared morphology differs dramatically from that at shorter wavelengths. The most luminous emission region in the galaxy is marginally resolved at 24 micron and has a projected separation of nearly 520 pc from the peak emission in the optical and near infrared. This region is responsible for ~9% of the total emission at 24 micron and is likely a young star formation region. We show that this and other compact sources account for more than 1/3 of the total 24 micron emission. We compute a total infrared luminosity for NGC 55 of 1.2*10^9 L_sun. The star formation rate implied by our measurements is 0.22 M_sun/yr. We demonstrate that the cold dust is more extended than the warm dust in NGC 55--the minor-axis scale heights are 0.32, 0.43, and 0.49 kpc at 24, 70 and 160 micron, respectively. The dust temperature map shows a range of temperatures that are well-correlated with the 24 micron surface brightness, from 20 K in low-surface-brightness regions to 26 K in high-surface-brightness regions.

We present Spitzer observations in five wavebands between 3.6 and 24um of an unbiased sample of 9 luminous, dusty galaxies selected at 1200um by the MAMBO camera on the IRAM 30-m telescope, a population akin to the well-known submm or `SCUBA' galaxies (hereafter SMGs). Owing to the coarse resolution of submm/mm instrumentation, SMGs have traditionally been difficult to identify at other wavelengths. We compare our multi-wavelength catalogs to show that the overlap between 24 and 1200um must be close to complete at these flux levels. We find that all (4/4) of the most secure >=4sigma SMGs have robust >=4sigma counterparts at 1.4GHz, while the fraction drops to 7/9 using all >=3sigma SMGs. We show that combining mid-IR and marginal (>=3sigma) radio detections provides plausible identifications in the remaining cases, enabling us to identify the complete sample. Accretion onto an obscured central engine is betrayed by the shape of the mid-IR continuum emission for several sources, confirming Spitzer's potential to weed out active galaxies. We demonstrate the power of a S(24um)/S(8um) vs S(8um)/S(4.5um) color-color plot as a diagnostic for this purpose. However, we conclude that the majority (~75%) of SMGs have rest-frame mid-/far-IR SEDs commensurate with obscured starbursts. Sensitive 24-um observations are clearly a useful route to identify and characterize reliable counterparts to high-redshift far-IR-bright galaxies, complementing what is possible via deep radio imaging.

We investigate the luminous X-ray sources in the Lockman Hole (LH) and the Extended Groth Strip (EGS) detected at 24microns using MIPS and also with IRAC on board Spitzer. We assemble optical/infrared spectral energy distributions (SEDs) for 45 X-ray/24micron sources in the EGS and LH. Only about 1/4 of the hard X-ray/24micron sources show pure type 1 AGN SEDs. More than half of the X-ray/24micron sources have stellar-emission-dominated or obscured SEDs, similar to those of local type 2 AGN and spiral/starburst galaxies. One-third of the sources detected in hard X-rays do not have a 24micron counterpart. Two such sources in the LH have SEDs resembling those of S0/elliptical galaxies. The broad variety of SEDs in the optical-to-Spitzer bands of X-ray selected AGN means that AGN selected according to the behavior in the optical/infrared will have to be supplemented by other kinds of data (e.g., X-ray) to produce unbiased samples of AGN.

We present preliminary results on 24micron detections of luminous infrared galaxies at z>1 with the Multiband Imaging Photometer for Spitzer (MIPS). Observations were performed in the Lockman Hole and the Extended Groth Strip (EGS), and were supplemented by data obtained with the Infrared Array Camera (IRAC) between 3 and 9microns. The positional accuracy of ~2arcsec for most MIPS/IRAC detections provides unambiguous identifications of their optical counterparts. Using spectroscopic redshifts from the Deep Extragalactic Evolutionary Probe survey, we identify 24micron sources at z>1 in the EGS, while the combination of the MIPS/IRAC observations with $BVRIJHK$ ancillary data in the Lockman Hole also shows very clear cases of galaxies with photometric redshifts at 1<z<2.5. The observed 24micron fluxes indicate infrared luminosities greater than 10^11 L_sol, while the data at shorter wavelengths reveal rather red and probably massive (M>=M*) galaxy counterparts. It is the first time that this population of luminous objects is detected up to z~2.5 in the infrared. Our work demonstrates the ability of the MIPS instrument to probe the dusty Universe at very high redshift, and illustrates how the forthcoming Spitzer deep surveys will offer a unique opportunity to illuminate a dark side of cosmic history not explored by previous infrared experiments.

IRAC imaging of a 4'7x4'7 area in the Lockman Hole detected over 400 galaxies in the IRAC 3.6 micron and 4.5 micron bands, 120 in the 5.8 micron, and 80 in the 8 micron bandin 30 minutes of observing time. Color-color diagrams suggest that about half of these galaxies are at redshifts 0.6<z<1.3 with about a quarter at higher redshifts (z>1.3). We also detect IRAC counterparts for 6 of the 7 SCUBA sources and all 9 XMM sources in this area. The detection of the counterparts of the SCUBA sources and galaxies at z>1.3 demonstrates the ability of IRAC to probe the universe at very high redshifts.

We present an updated phenomenological galaxy evolution model to fit the Spitzer 24, 70 and 160 microns number counts as well as all the previous mid and far infrared observations. Only a minor change of the co-moving luminosity density distribution in the previous model (Lagache, Dole, Puget 2003), combined with a slight modification of the starburst template spectra mainly between 12 and 30 microns, are required to fit all the data available. We show that the peak in the MIPS 24 micron counts is dominated by galaxies with redshift between 1 and 2, with a non negligible contribution from the z>2 galaxies (~30% at S=0.2 mJy). The very close agreement between the model and number counts at 15 and 24 microns strikingly implies that (1) the PAHs (Policyclic Aromatic Hydrocarbons) features remain prominent in the redshift band 0.5 to 2.5 and (2) the IR energy output has to be dominated by ~3 10^11 Lo to ~3 10^12 Lo galaxies from redshift 0.5 to 2.5. Combining Spitzer with the Infrared Space Observatory (ISO) deep cosmological surveys gives for the first time an unbiased view of the infrared Universe from z=0 to z=2.5.

We use the source counts measured with the Multiband Imaging Photometer for Spitzer (MIPS) at 24, 70, and 160 microns to determine the 5-sigma confusion limits due to extragalactic sources: 56 micro-Jy, 3.2 and 40 mJy at 24, 70 and 160 microns, respectively. We also make predictions for confusion limits for a number of proposed far infrared missions of larger aperture (3.5 to 10m diameter).

Galaxy source counts in the infrared provide strong constraints on the evolution of the bolometric energy output from distant galaxy populations. We present the results from deep 24 micron imaging from Spitzer surveys, which include approximately 50,000 sources to an 80% completeness of 60 uJy. The 24 micron counts rapidly rise at near-Euclidean rates down to 5 mJy, increase with a super-Euclidean rate between 0.4 - 4 mJy, and converge below 0.3 mJy. The 24 micron counts exceed expectations from non-evolving models by a factor >10 at 0.1 mJy. The peak in the differential number counts corresponds to a population of faint sources that is not expected from predictions based on 15 micron counts from ISO. We argue that this implies the existence of a previously undetected population of infrared-luminous galaxies at z ~ 1-3. Integrating the counts to 60 uJy, we derive a lower limit on the 24 micron background intensity of 1.9 +/- 0.6 nW m-2 sr-1 of which the majority (~ 60%) stems from sources fainter than 0.4 mJy. Extrapolating to fainter flux densities, sources below 60 uJy contribute 0.8 {+0.9/-0.4} nW m-2 sr-1 to the background, which provides an estimate of the total 24 micron background of 2.7 {+1.1/-0.7} nW m-2 sr-1.

We present sub-millimeter statistical detections of galaxies discovered in the 5'x5' Spitzer Early Release Observations (to 4-15 microJy 5 sigma at 3.6-8 microns, 170 microJy at 24 microns) through a stacking analysis of our reanalysed SCUBA 8mJy survey maps, and a Spitzer identification of a new sub-millimeter point source in the 8mJy survey region. For sources detected at 5.8 or 8 microns (154 and 111 sources respectively), we detect positive skews in the sub-millimeter flux distributions at 99.2-99.8% confidence using Kolmogorov-Smirnov tests, at both 850 microns and 450 microns. We also marginally detect the Spitzer 24 micron galaxies at 850 microns at 97% confidence, and place limits on the mean sub-millimeter fluxes of the 3.6 and 4.5 micron sources. Integrating the sub-millimeter fluxes of the Spitzer populations, we find the 5.8 micron galaxies contribute 0.12 +/- 0.05 nW/m^2/sr to the 850 micron background, and 2.4 +/- 0.7 nW/m^2/sr to the 450 micron background; similar contributions are made by the 8 micron-selected sample. We infer that the populations dominating the 5.8 and 8 micron extragalactic background light also contribute around a quarter of the 850 micron background and the majority of the 450 micron background.

We investigate Extremely Red Objects (EROs) using near- and mid-infrared observations in five passbands (3.6 to 24 micron) obtained from the Spitzer Space Telescope, and deep ground-based R and K imaging. The great sensitivity of the IRAC camera allows us to detect 64 EROs in only 12 minutes of IRAC exposure time, by means of an R-[3.6] color cut (analogous to the traditional red R-K cut). A pure infrared K-[3.6] red cut detects a somewhat different population and may be more effective at selecting z > 1.3 EROs. We find 17% of all galaxies detected by IRAC at 3.6 or 4.5 micron to be EROs. These percentages rise to about 40% at 5.8 micron, and about 60% at 8.0 micron. We utilize the spectral bump at 1.6 micron to divide the EROs into broad redshift slices using only near-infrared colors (2.2/3.6/4.5 micron). We conclude that two-thirds of all EROs lie at redshift z > 1.3. Detections at 24 micron imply that at least 11% of 0.6 < z < 1.3 EROs and at least 22% of z > 1.3 EROs are dusty star-forming galaxies.

We report the discovery of two new gravitationally lensed quasars, WFI J2026-4536 and WFI J2033-4723, at respective source redshifts of z=2.23 and z=1.66. Both systems are quadruply imaged and have similar PG1115-like image configurations. WFI J2026-4536 has a maximum image separation of 1.4", a total brightness of g = 16.5, and a relatively simple lensing environment, while WFI J2033-4723 has a maximum image separation of 2.5", an estimated total brightness of g = 17.9, and a more complicated environment of at least six galaxies within 20". The primary lensing galaxies are detected for both systems after PSF subtraction. Several of the broadband flux ratios for the two lenses show a strong (0.1-0.4 mags) trend with wavelength, suggesting either microlensing or differential extinction through the lensing galaxy. For WFI J2026-4536, the total quasar flux has dimmed by 0.1 mag in the blue but only half as much in the red over three months, suggestive of microlensing-induced variations. For WFI J2033-4723, resolved spectra of some of the quasar components reveal emission line flux ratios that agree better with the macromodel predictions than either the broadband or continuum ratios, also indicative of microlensing. The predicted differential time delays for WFI J2026-4536 are short, ranging from 1-2 weeks for the long delay, but are longer for WFI J2033-4723, ranging from 1-2 months. Both systems hold promise for future monitoring campaigns aimed at microlensing or time delay studies.

The Palomar Integral Field Spectrograph was used to probe a variety of environments in nine nearby galaxies that span a range of morphological types, luminosities, metallicities, and infrared-to-blue ratios. For the first time, near-infrared spectroscopy was obtained for nuclear or bright HII regions in star-forming galaxies over two spatial dimensions (5.7"x10.0") in the [FeII] (1.257um), [FeII] (1.644um), Paschen beta (1.282um), H2 (2.122um), and Brackett gamma (2.166um) transition lines. These data yield constraints on various characteristics of the star-forming episodes in these regions, including their strength, maturity, spatial variability, and extinction. The HII regions stand out from the nuclei. Unlike observations of nuclear regions, HII region near-infrared observations do not show a spatial coincidence of the line and continuum emission; the continuum and line maps of HII regions usually show distinct and sometimes spatially-separated morphologies. Gauging from Paschen beta and Brackett gamma equivalent widths and luminosities, the HII regions have younger episodes of star formation than the nuclei and more intense radiation fields. Near-infrared line ratio diagnostics suggest that HII regions have "purer" starbursting properties. The correlation between ionizing photon density and mid-infrared color is consistent with the star formation activity level being higher for HII regions than for nuclei. And though the interpretation is complicated, on a purely empirical basis the HII regions show lower Fe1+ abundances than nuclei by an order of magnitude.

Images at 12.5 microns of nuclei of three nearby Seyfert galaxies -- NGC 1275, NGC 4151 and NGC 7469 -- have been obtained with the Keck 10-m Telescope. NGC 7469 is resolved and deconvolution delineates a structure <0.04"x0.08" or <13x26 pc at a position angle of 135deg. From a comparison with structure seen at millimeter wavelengths, this structure is interpreted as a disk aligned with the molecular gas in the central few hundred parsecs of the galaxy. NGC 1275 and NGC 4151 are not resolved; limits on the sizes of these nuclei are 0.08" and 0.16", corresponding to physical spatial scales of 28 and 10 pc. The lower limits to the brightness temperatures implied by these size limits and the measured flux densities are within ~50K of the 12 micron to 25 micron color temperatures of these systems as inferred from IRAS observations. The angular size limits are within a factor of 2.5 of the sizes required to spatially resolve thermal emission from dust heated by a central luminosity source. These sizes preclude significant contributions to the nuclear infrared emission from star forming regions.

We present near-infrared imaging and spectroscopic observations of two FR II high-redshift radio galaxies (HzRGs), 4C 40.36 (z=2.3) and 4C 39.37 (z=3.2), obtained with the Hubble, Keck, and Hale Telescopes. High resolution images were taken with filters both in and out of strong emission lines, and together with the spectroscopic data, the properties of the line and continuum emissions were carefully analyzed. Our analysis of 4C 40.36 and 4C 39.37 shows that strong emission lines (e.g., [O III] 5007 A and H alpha+[N II]) contribute to the broad-band fluxes much more significantly than previously estimated (80% vs. 20-40%), and that when the continuum sources are imaged through line-free filters, they show an extremely compact morphology with a high surface brightness. If we use the R^1/4-law parametrization, their effective radii (r(e)) are only 2-3 kpc while their restframe B-band surface brightnesses at r(e) are I(B) ~ 18 mag/arcsec^2. Compared with z ~ 1 3CR radio galaxies, the former is x3-5 smaller, while the latter is 1-1.5 mag brighter than what is predicted from the I(B)-r(e) correlation. Although exponential profiles produce equally good fits for 4C 40.36 and 4C 39.37, this clearly indicates that with respect to the z~1 3CR radio galaxies, the light distribution of these two HzRGs is much more centrally concentrated. Spectroscopically, 4C 40.36 shows a flat (fnu=const) continuum while 4C 39.37 shows a spectrum as red as that of a local giant elliptical galaxy. Although this difference may be explained in terms of a varying degree of star formation, the similarities of their surface brightness profiles and the submillimeter detection of 4C 39.37 might suggest that the intrinsic spectra is equally blue (young stars or an AGN), and that the difference is the amount of reddening.

We present revised infrared (2.2 micron) astrometry of the field containing the Galactic microquasar GRS 1758-258, using observations at the Keck I 10-m telescope. We find three candidates for the microquasar within a 3 sigma error circle, but none within 2 sigma. We show that if the 18.4 day X-ray period of GRS 1758-258 is due to a binary orbit, then only one of the three candidates, an early K-type giant, is large enough to power the microquasar via Roche lobe overflow. We therefore identify this star as the infrared counterpart of GRS 1758-258, which we classify as a low mass X-ray binary. Long term infrared monitoring of this source should provide further information about the microquasar system, including a confirmation of the X-ray period and an estimate of the compact object's mass.

Low spectral resolution (delta lambda/lambda ~ 50) mid-infrared observations with high angular resolution (0.3--0.5'') using the Long Wavelength Spectrometer on the Keck I Telescope are reported of the nuclei of five highly luminous infrared bright galaxies. Spectra of eight distinct nuclei, ranging in luminosity from 10^11 L_sun to >10^12 L_sun have been obtained. Four of the nuclei show the characteristic PAH emission features, i.e. 11.3 micron emission as well as the 8.6 micron shoulder of the 7.7 micron band. The other nuclei show either weak PAH emission bands or no evidence for these bands. The high spatial resolution of the observations reveals extended emission in the 11.3 micron PAH band associated with several of the compact nuclear sources. When proper account is taken of the diffuse PAH emission, most of the compact sources show little or no directly associated PAH emission. The diffuse PAH emission is extended over spatial scales of 100--500 pc; its presence shows that there is significant circumnuclear UV/optical emission exciting the aromatic bands, most likely associated with circumnuclear starbursts. After the spectra of the nuclear sources are corrected for the spectrum of the diffuse PAH emission, the peak apparent silicate optical depth at 9.7 microns can be as great as 15, corresponding to > 150 magnitudes of visible light extinction. Because of the large silicate optical depths, mid-infrared spectra are not probing the nature of the true nuclei in the most opaque compact sources.

Observations for seven infrared luminous starburst galaxies are reported in the mid-infrared from 8 - 18 microns using the Keck Telescopes with spatial resolution approaching the diffraction limit. All of the galaxies observed show evidence of strong interactions based on optical morphologies. For these galaxies, a substantial fraction, usually more than 50%, of the infrared luminosity is generated in regions ranging in sizes from 100pc -- 1 Kpc. Nuclear starbursts often dominate the infrared luminosity, but this is not always true. In some galaxies, most notably NGC 6090, substantial infrared luminosity greatly in excess of the nuclear luminosity is generated in regions associated with the physical interaction between two galaxies. The radio emission is a good tracer of the location of high luminosity young stars. The visual/ultraviolet radiation output of the nearby star forming galaxies is dominated by emission from regions that are generally not producing the copious infrared luminosity of the systems. The regions of high infrared luminosity in local starburst galaxies are significantly smaller than the galaxies as a whole. The integrated spectral energy distributions (SEDs) of these galaxies are very different from the SEDs of the regions of star formation. If the SEDs of star-forming regions in these galaxies reflect the SEDs found in forming galaxies at high redshift, the distant galaxies should be dominated by the mid and far-infrared luminosity output far more than the integrated luminous output of nearby starburst galaxies would suggest.

We present deep infrared ($2.2 \mu$m) imaging of the Galactic microquasars 1E1740-2942 and GRS 1758-258 using the Keck-I 10-meter telescope in June 1998. The observations were taken under excellent seeing conditions ($\sim 0.45 \arcsec$ full-width half-maximum), making them exceptionally deep for these crowded fields. We used the USNO-A2.0 catalog to astrometrically calibrate the infrared images (along with an optical CCD image in the case of GRS 1758-258), providing independent frame ties to the known radio positions of the objects. For 1E1740-2942, we confirm potential candidates for the microquasar previously identified by Marti et al., and show that none of the objects near the microquasar have varied significantly from 1998 to 1999. For GRS 1758-258, our astrometry indicates a position shifted from previous reports of candidates for the microquasar. We find no candidates inside our 90% confidence radius to a $2 \sigma$ limiting magnitude of $K_s = 20.3$ mag. We discuss the implications of these results for the nature of the microquasar binary systems.

Mid-infrared observations of the central source of NGC 1068 have been obtained with a spatial resolution in the deconvolved image of 0.1" (~ 7pc). The central source is extended by ~1" in the north-south direction but appears unresolved in the east-west direction over most of its length. About two-thirds of its flux can be ascribed to a core structure which is itself elongated north-south and does not show a distinct unresolved compact source. The source is strongly asymmetric, extending significantly further to the north than to the south. The morphology of the mid-infrared emission appears similar to that of the radio jet, and has features which correlate with the images in [O III]. Its 12.5-24.5 micron color temperature ranges from 215 to 260 K and does not decrease smoothly with distance from the core. Silicate absorption is strongest in the core and to the south and is small in the north. The core, apparently containing two-thirds of the bolometric luminosity of the inner 4" diameter area, may be explained by a thick, dusty torus near the central AGN viewed at an angle of ~65 deg to its plane. There are, however, detailed difficulties with existing models, especially the narrow east-west width of the thin extended mid-infrared "tongue" to the north of the core. We interpret the tongue as re-processed visual and ultraviolet radiation that is strongly beamed and that originates in the AGN.

We present Keck high-resolution near-IR (2.2 microns; FWHM~0.15") and mid-IR (12.5 microns; FWHM~0.4") images of APM08279+5255, a z=3.91 IR-luminous BALQSO with a prodigious apparent bolometric luminosity of 5x10^{15} Lsun, the largest known in the universe. The K-band image shows that this system consists of three components, all of which are likely to be the gravitationally lensed images of the same background object, and the 12.5 micron image shows a morphology consistent with such an image configuration. Our lens model suggests that the magnification factor is ~100 from the restframe UV to mid-IR, where most of the luminosity is released. The intrinsic bolometric luminosity and IR luminosity of APM08279+5255 are estimated to be 5x10^{13} Lsun and 1x10^{13} Lsun, respectively. This indicates that APM 08279+5255 is intriniscally luminous, but it is not the most luminous object known. As for its dust contents, little can be determined with the currently available data due to the uncertainties associated with the dust emissivity and the possible effects of differential magnification. We also suggest that the lensing galaxy is likely to be a massive galaxy at z~3.

Observations of ultraluminous infrared galaxies (ULIRGs) with an achieved resolution approaching the diffraction limit in the mid-infrared from 8 - 25 $\mu$m using the Keck Telescopes are reported. We find extremely compact structures, with spatial scales of $< 0.3''$ (diameter) in six of the seven ULIRGs observed. These compact sources emit between 30% and 100% of the mid-infrared energy from these galaxies. We have utilized the compact mid-infrared structures as a diagnostic of whether an AGN or a compact (100 -- 300 pc) starburst is the primary power source in these ULIRGs. In Markarian 231, the upper limit on the diameter of the 12.5 $\mu$m source, 0.13$''$, shows that the size of the infrared source must increase with increasing wavelength, consistent with AGN models. In IRAS 05189-2524 and IRAS 08572+3915 there is strong evidence that the source size increases with increasing wavelength. This suggests heating by a central source rather than an extended luminosity source, consistent with the optical classification as an AGN. The compact mid-infrared sources seen in the other galaxies cannot be used to distinguish the ultimate luminosity source. If these ULIRGs are powered by compact starbursts, the star formation rates seen in the central few hundred parsecs far exceed the global rates seen in nearby starburst galaxies, and approach the surface brightness of individual clusters in nearby starburst galaxies.

Using NICMOS on HST, we have imaged the emission-line nebulae and the line-free continuum in 4C 40.36, a ultra-steep spectrum FR II radio galaxy at z=2.269. The line-free continuum was found to be extremely compact and symmetric while the emission-line nebulae seen in H-alpha+[N II] show very clumpy structures spreading almost linearly over 16 kpc. However, this linear structure is clearly misaligned from the radio axis. The SED of the line-free continuum is very flat, suggesting that if the continuum emission is produced by a single source, it is likely to be a young bursting stellar population or scattered AGN light. However, because of the lack of a line-free optical image with a comparable spatial resolution, we cannot exclude the possibility that the observed SED is a composite of a young blue population and an old red population.

Images of Arp 220 from 3.45 to 24.5 $\mu$m with 0.5'' resolution are presented which clearly separate the nucleus into at least two components. The western component is about three times more luminous than the eastern component, but the silicate absorption in the fainter, eastern component is roughly 50% greater than the absorption in the western component. Each component is marginally resolved. The two components seen at 24.5 $\mu$m are identified with the two radio components. The western source most likely coincides with the high extinction disk previously suggested to exist in Arp 220, while the eastern nucleus is identified with a faint highly reddened source seen in HST 2.2 $\mu$m NICMOS images. The two nuclei together account for essentially all of the measured 24.5 $\mu$m flux density. Two models are presented, both of which fit the observations. In one, the majority of the total luminosity is produced in an extended star formation region and in the other most of the luminosity is produced in the compact but extincted regions associated with the two nuclei seen at 24.5 $\mu$m. In both pictures, substantial luminosity at 100 $\mu$m emerges from a component having a diameter of 2--3'' ($\sim$1kpc).

We present a deep BVrIK multicolor catalog of galaxies in the field of the high redshift (z=4.7) quasar BR 1202-0725. Reliable colors have been measured for galaxies selected down to R=25. The choice of the optical filters has been optimized to define a robust multicolor selection of galaxies at 3.8<z< 4.5. Within this interval the surface density of galaxy candidates with z~4 in this field is 1 arcmin^{-2}. Photometric redshifts have been derived for the galaxies in the field with the maximum likelihood analysis using the GISSEL library of ~10^6 synthetic spectra. The accuracy of the method used has been discussed and tested using galaxies in the Hubble Deep Field (HDF) with known spectroscopic redshifts and accurate photometry. A peak in the redshift distribution is present at z~0.6 with relatively few galaxies at z>1.5. At variance with brighter surveys (I<22.5) there is a tail in the distribution towards high redshifts up to z~4. The luminosity function at z~0.6 shows a steepening for M_B>-19. The observed cosmological ultraviolet luminosity density is computed in the overall redshift interval z=0.3-4.5 reaching a value 2x10^{19} W/Hz/Mpc^3 at z~0.8. We have derived in a homogeneous way, using the GISSEL libraries, the physical parameters connected with the fitted spectral energy distributions. The bulk of the blue intermediate redshift population with z=0.4-1 mostly consists of very young star-forming galaxies with a median starburst age of the order of a few 10^8 yr and typical mass in luminous stars ~2x10^8 Mo.

We present an ISO SWS and ISOPHOT-S, mid-infrared spectroscopic survey of 15 ultra-luminous IRAS galaxies. We combine the survey results with a detailed case study, based on near-IR and mm imaging spectroscopy, of one of the sample galaxies (UGC 5101). We compare the near- and mid-IR characteristics of these ultra-luminous galaxies to ISO and literature data of thirty starburst and active galactic nuclei (AGN), template galaxies. We find that 1) 70-80% of the ultra-luminous IRAS galaxies in our sample are predominantly powered by recently formed massive stars. 20-30% are powered by a central AGN. These conclusions are based on a new infrared 'diagnostic diagram' involving the ratio of high to low excitation mid-IR emission lines on the one hand, and on the strength of the 7.7um PAH feature on the other hand. 2) at least half of the sources probably have simultaneously an active nucleus and starburst activity in a 1-2 kpc diameter circum-nuclear disk/ring. 3) the mid-infrared emitting regions are highly obscured. After correction for these extinctions, we estimate that the star forming regions in ULIRGs have ages between 10^7 and 10^8 years, similar to but somewhat larger than those found in lower luminosity starburst galaxies. 4) in the sample we have studied there is no obvious trend for the AGN component to dominate in the most compact, and thus most advanced mergers. Instead, at any given time during the merger evolution, the time dependent compression of the circum-nuclear interstellar gas, the accretion rate onto the central black hole and the associated radiation efficiency may determine whether star formation or AGN activity dominates the luminosity of the system.

A large deep and nearly complete B<24.5 redshift sample is used to measure the change in distribution function of the stellar mass production rate in individual galaxies with redshift. The evolution of the star formation rate distribution with redshift is interpreted in terms of the history of spiral galaxy formation, with the disk component modelled as a single evolving entity, and the characteristic timescales, luminosities, and epochs varying according to galaxy type. The more massive forming galaxies seen at z=1 to 3 are identified as earlier type spirals, whose star formation rates are initially high and then decline rapidly at z<1, while for later type spirals and smaller mass irregulars the mass formation rates at z<1 are lower, and the formation process persists to redshifts much closer to the present epoch. We find that these models can be consistent with the data and fit well into a broad picture of other recent results if qnought=0.02 and many of the disks begin their growth at z<<3, but that they predict too many bright star formers at high z in flat universes.

We present the results of ground-based and HST imaging studies targeted on z>4.5 quasar fields. High-redshift galaxies identified in deep narrow-band Lyman alpha images of the fields surrounding the quasars BR1202-0725 (z=4.694) and BR2237-0607 (z=4.558) have been confirmed with follow-up spectroscopy using the LRIS spectrograph on the Keck 10m telescope. These high-redshift galaxies are typically only moderately luminous objects and their sub-L* magnitudes imply star formation rates of only a few solar masses per year. For BR1202-0725 tip-tilt imaging in Lyman alpha at the UH 2.2m telescope and deep narrow-band IR imaging in [OII] at CFHT reveals a complex and structured system which may be merging to produce the host galaxy of the quasar. For BR2237-0607 the Lyman-alpha-emitting galaxies show the strong Lyman alpha emission and weak continuum expected for objects before substantial chemical enrichment and consequent dust formation and extinction has occurred. These objects appear to be galaxies in their first outburst of star formation.

Using the OH-airglow suppressor spectrograph at the University of Hawaii 2.2m telescope and the CGS4 spectrometer at the United Kingdom Infrared Telescope, we have found exceptionally large Balmer decrements in two unusual high-z QSOs, Hawaii 167 (z=2.36, Ha/Hb = 13) and Q0059-2735 (z=1.59, Ha/Hb = 7.6), the latter being a so-called low-ionization broad absorption line QSO (BALQSO). We argue that these objects are young QSOs heavily enshrouded by dust. In fact, the internal reddening might be so large as to completely extinguish the QSO light in the restframe UV, allowing us to see the underlying stellar population. Our possible detection of the 4000 A break in Hawaii 167 supports this idea. Its small amplitude indicates a very young age for the population, ~15 Myrs. To explain the properties of these QSOs, we propose a model in which a young QSO is surrounded by a shell of young massive stars mixed with significant amounts of dust. We predict that as the QSO emerges from this dust cocoon, it will eventually take on the appearance of a normal BALQSO.

We report the detection of a nearby emission-line companion to the z=4.695 quasar BR1202-0725. Deep narrow-band exposures on this field from the UH 2.2 m show a Ly alpha flux of 1.5\times\ten{-16} ergs cm^{-2} s^{-1}. High-resolution HST WFC2 imaging in the F814W filter band shows continuum structure near the emission position, at 2.6" NW of the quasar, corresponding to a projected separation of ~7.5 h^{-1} kpc for \qnought=0.5, where h=\hnought/100\kms Mpc^{-1}. We discuss possible explanations for the combined line and color properties. The ionization is most likely produced by the quasar, but if due to underlying star formation would require a star formation rate of ~7 h^2 solar masses per yr.