We present the discovery of the second binary with a Roche lobe-filling hot subdwarf transferring mass to a white dwarf (WD) companion. This 56 minute binary was discovered using data from the Zwicky Transient Facility. Spectroscopic observations reveal an He-sdOB star with an effective temperature of $T_{\rm eff}=33,700\pm1000$ K and a surface gravity of $log(g)=5.54\pm0.11$. The GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the He-sdOB star and shows an eclipse of the He-sdOB by an accretion disk as well as a weak eclipse of the WD. We infer a He-sdOB mass of $M_{\rm sdOB}=0.41\pm0.04$ M$_\odot$ and a WD mass of $M_{\rm WD}=0.68\pm0.05$ M$_\odot$. The weak eclipses imply a WD black-body temperature of $63,000\pm10,000$ K and a radius $R_{\rm WD}=0.0148\pm0.0020$ M$_\odot$ as expected for a WD of such high temperature. The He-sdOB star is likely undergoing hydrogen shell burning and will continue transferring mass for $\approx1$ Myrs at a rate of $10^{-9} M_\odot {\rm yr}^{-1}$ which is consistent with the high WD temperature. The hot subdwarf will then turn into a WD and the system will merge in $\approx30$ Myrs. We suggest that Galactic reddening could bias discoveries towards preferentially finding Roche lobe-filling systems during the short-lived shell burning phase. Studies using reddening corrected samples should reveal a large population of helium core-burning hot subdwarfs with $T_{\rm eff}\approx25,000$ K in binaries of 60-90 minutes with WDs. Though not yet in contact, these binaries would eventually come into contact through gravitational wave emission and explode as a sub-luminous thermonuclear supernova or evolve into a massive single WD.

Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

Helium is expected to be present in the ejecta of some hydrogen-poor superluminous supernovae (SLSN-I). However, so far only one event, PTF10hgi has been identified with He features in its photospheric spectra (Quimby et al. 2018). At $z=0.0866$, ZTF19aawfbtg (SN2019hge) has more than 10 optical spectra at phases from $-41$ to $+103$ days. Most of these spectra match well with that of PTF10hgi, suggesting SN2019hge as the second case of a He-rich SLSN-I. Confirmation comes from its near-IR spectrum taken at $+34$ days, revealing He I features with P-Cygni profiles at 1.083 and 2.058$\mu$m. Using the optical spectra of PTF10hgi and SN2019hge as templates, we examine 70 SLSN-I discovered by ZTF in the first two years of operation and found additional five SLSN-I with distinct He-features, yielding a lower limit to the He-rich fraction of $\sim9$\%. It is suggested that the excitation of He I atoms requires non-thermal radiation. These He-rich events can not be explained by traditional $^{56}$Ni mixing model because of their blue spectra, high peak luminosity and slow rising time scales. Magnetar models offer a possible solution since pulsar winds naturally generate high energy particles as sources of non-thermal excitation. However, how the magnetar energy is transported to the outer layers is unclear. An alternative model is ejecta interaction with H-poor CSM which may be supported by the observed light curve undulations. These six SLSN-Ib appear to have relatively low-peak luminosities, around $-20.06\pm0.16$ mag ($g$-band, rest-frame), in comparison to the full range of $-20$ to $-23$ mag covered by SLSN-I.

We report the discovery and panchromatic followup observations of the young Type Ic supernova, SN2020oi, in M100, a grand design spiral galaxy at a mere distance of $14$ Mpc. We followed up with observations at radio, X-ray and optical wavelengths from only a few days to several months after explosion. The optical behaviour of the supernova is similar to those of other normal Type Ic supernovae. The event was not detected in the X-ray band but our radio observation revealed a bright mJy source ($L_{\nu} \approx 1.2 \times 10^{27} {\rm erg\,s}^{-1} {\rm Hz}^{-1}$). Given, the relatively small number of stripped envelope SNe for which radio emission is detectable, we used this opportunity to perform a detailed analysis of the comprehensive radio dataset we obtained. The radio emitting electrons initially experience a phase of inverse Compton cooling which leads to steepening of the spectral index of the radio emission. Our analysis of the cooling frequency points to a large deviation from equipartition at the level of $\epsilon_e/\epsilon_B \gtrsim 200$, similar to a few other cases of stripped envelope SNe. Our modeling of the radio data suggests that the shockwave driven by the SN ejecta into the circumstellar matter (CSM) is moving at $\sim 3\times 10^{4}\,{\rm km\,s}^{-1}$. Assuming a constant mass-loss from the stellar progenitor, we find that the mass-loss rate is $\dot{M} \approx 1.4\times 10^{-4}\,{M}_{\odot}\,{\rm yr}^{-1}$, for an assumed wind velocity of $1000\,{\rm km\,s}^{-1}$. The temporal evolution of the radio emission suggests a radial CSM density structure steeper than the standard $r^{-2}$.

We present ZTF20aajnksq (AT2020blt), a fast-fading ($\Delta r=2.4$ mag in $\Delta t=1.3$ days) red ($g-r\approx0.6$ mag) and luminous ($M_{1626}=-25.9$) optical transient at $z=2.9$ discovered by the Zwicky Transient Facility (ZTF). AT2020blt shares several features in common with afterglows to long-duration gamma-ray bursts (GRBs): (1) an optical light curve well-described by a broken power-law with a break at $t_\mathrm{j}=1$ day (observer-frame); (2) a luminous $(L_X = 10^{46}$ $\mathrm{erg}$ $\mathrm{s}^{-1})$ X-ray counterpart; and (3) luminous ($L_\nu = 4 \times 10^{31}$ $\mathrm{erg}$ $\mathrm{sec}^{-1}$ $\mathrm{Hz}^{-1}$ at 10 GHz) radio emission. However, no GRB was detected in the 0.74d between the last ZTF non-detection ($r > 20.64$) and the first ZTF detection ($r = 19.57$), with an upper limit on the isotropic-equivalent gamma-ray energy release of $E_{\gamma,\mathrm{iso}} < 7 \times 10^{52}$ erg. AT2020blt is thus the third afterglow-like transient discovered without a detected GRB counterpart (after PTF11agg and ZTF19abvizsw) and the second (after ZTF19abvizsw) with a redshift measurement. We conclude that the properties of AT2020blt are consistent with a classical (initial Lorentz factor $\Gamma_0 \gtrsim 100$) on-axis GRB that was missed by high-energy satellites. Furthermore, by estimating the rate of transients with light curves similar to that of AT2020blt in ZTF high-cadence data, we agree with previous results that there is no evidence for an afterglow-like phenomenon that is significantly more common than classical GRBs. We conclude by discussing the status and future of fast-transient searches in wide-field high-cadence optical surveys.

We present observations of ZTF18abfcmjw (SN2019dge), a helium-rich supernova with a fast-evolving light curve indicating an extremely low ejecta mass ($\approx 0.3\,M_\odot$) and low kinetic energy ($\approx 1.2\times 10^{50}\,{\rm erg}$). Early-time (<4 d after explosion) photometry reveal evidence of shock cooling from an extended helium-rich envelope of $\sim0.1\,M_\odot$ located at $\sim 3\times 10^{12}\,{\rm cm}$ from the progenitor. Early-time He II line emission and subsequent spectra show signatures of interaction with helium-rich circumstellar material, which extends from $\gtrsim 5\times 10^{13}\,{\rm cm}$ to $\gtrsim 2\times 10^{16}\,{\rm cm}$. We interpret SN2019dge as a helium-rich supernova from an ultra-stripped progenitor, which originates from a close binary system consisting of a mass-losing helium star and a low-mass main sequence star or a compact object (i.e., a white dwarf, a neutron star, or a black hole). We infer that the local volumetric birth rate of 19dge-like ultra-stripped SNe is in the range of 1400--8200$\,{\rm Gpc^{-3}\, yr^{-1}}$ (i.e., 2--12% of core-collapse supernova rate). This can be compared to the observed coalescence rate of compact neutron star binaries that are not formed by dynamical capture.

Cosmic neutrinos provide a unique window into the otherwise-hidden mechanism of particle acceleration in astrophysical objects. A flux of high-energy neutrinos was discovered in 2013, and the IceCube Collaboration recently associated one high-energy neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. However a combined analysis of many similar active galaxies revealed no excess from the broader population, leaving the vast majority of the cosmic neutrino flux unexplained. Here we present the association of a radio-emitting tidal disruption event (AT2019dsg) with another high-energy neutrino, identified as part of our systematic search for optical counterparts to high-energy neutrinos with the Zwicky Transient Facility (ZTF). The probability of finding any radio-emitting tidal disruption event by chance is 0.5%, while the probability of finding one as bright in bolometric energy flux as AT2019dsg is 0.2%. Our electromagnetic observations can be explained through a multi-zone model, with radio analysis revealing a central engine, embedded in a UV photosphere, that powers an extended synchrotron-emitting outflow. This provides an ideal site for PeV neutrino production. The association suggests that tidal disruption events contribute to the cosmic neutrino flux. Unlike previous work which considered the rare subset of tidal disruption events with relativistic jets, our observations of AT2019dsg suggest an empirical model with a mildly-relativistic outflow.

Colors of Type Ia supernovae in the first few days after explosion provide a potential discriminant between different models. In this paper, we present $g-r$ colors of 65 Type Ia supernovae discovered within 5 days from first light by the Zwicky Transient Facility in 2018, a sample that is about three times larger than that in the literature. We find that $g-r$ colors are intrinsically rather homogeneous at early phases, with about half of the dispersion attributable to photometric uncertainties ($\sigma_\mathrm{noise}\sim\sigma_\mathrm{int}\sim$ 0.18 mag). Colors are nearly constant starting from 6 days after first light ($g-r\sim-0.15$ mag), while the time evolution at earlier epochs is characterized by a continuous range of slopes, from events rapidly transitioning from redder to bluer colors (slope of $\sim-0.25$ mag day$^{-1}$) to events with a flatter evolution. The continuum in the slope distribution is in good agreement both with models requiring some amount of $^{56}$Ni mixed in the outermost regions of the ejecta and with "double-detonation" models having thin helium layers ($M_\mathrm{He}=0.01\,M_\odot$) and varying carbon-oxygen core masses. At the same time, six events show evidence for a distinctive "red bump" signature predicted by "double-detonation" models with larger helium masses. We finally identify a significant correlation between the early-time $g-r$ slopes and supernova brightness, with brighter events associated to flatter color evolution (p-value=0.006). The distribution of slopes, however, is consistent with being drawn from a single population, with no evidence for two components as claimed in the literature based on $B-V$ colors.

We present optical, radio, and X-ray observations of SN2020bvc (=ASASSN20bs; ZTF20aalxlis), a nearby ($z=0.0252$; $d$=114 Mpc) broad-lined (BL) Type Ic supernova (SN). Our observations show that SN2020bvc shares several properties in common with the Ic-BL SN2006aj, which was associated with the low-luminosity gamma-ray burst (LLGRB) 060218. First, the 10 GHz radio light curve is on the faint end of LLGRB-SNe ($L_\mathrm{radio} \approx 10^{37}$erg/s): we model our VLA observations (spanning 13-43 d) as synchrotron emission from a mildly relativistic ($v \gtrsim 0.3c$) forward shock. Second, with Swift and Chandra we detect X-ray emission ($L_X \approx 10^{41}$erg/s) that is not naturally explained as inverse Compton emission or as part of the same synchrotron spectrum as the radio emission. Third, high-cadence ($6\times$/night) data from the Zwicky Transient Facility (ZTF) shows a double-peaked optical light curve, the first peak from shock-cooling emission from extended low-mass material (mass $M<10^{-2} M_\odot$ at radius $R>10^{12}$cm) and the second peak from the radioactive decay of Ni-56. SN2020bvc is the first confirmed double-peaked Ic-BL SN discovered without a GRB trigger, and shows X-ray and radio emission similar to LLGRB-SNe: this is consistent with models in which the same mechanism produces both the LLGRB and the shock-cooling emission. For four of the five other nearby ($z\lesssim0.05$) Ic-BL SNe with ZTF high-cadence data, we rule out a first peak like that seen in SN2006aj and SN2020bvc, i.e. that lasts $\approx 1$d and reaches a peak luminosity $M \approx -18$. X-ray and radio follow-up observations of future such events will establish whether double-peaked optical light curves are indeed predictive of LLGRB-like X-ray and radio emission.

(Abridged) Using the Zwicky Transient Facility alert stream, we are conducting a large campaign to spectroscopically classify all transients occurring in galaxies in the Census of the Local Universe (CLU) catalog. The aim of the experiment is to construct a spectroscopically complete, volume-limited sample of transients coincident within 100" of CLU galaxies out to 200 Mpc, and to a depth of 20 mag. We describe the survey design and spectroscopic completeness from the first 16 months of operations. We present results from a systematic search for Calcium rich gap transients in the sample of 22 low luminosity (peak absolute magnitude $M > -17$), hydrogen poor events found in the experiment (out of 754 spectroscopically classified SNe). We report the detection of eight Calcium rich gap transients, and constrain their volumetric rate to be at least $\approx 15\pm5$% of the SN Ia rate. Combining this sample with ten events from the literature, we find a likely continuum of spectroscopic properties ranging from events with SN Ia-like features (Ca-Ia objects) to SN Ib/c-like features (Ca-Ib/c objects) at peak light. Within the Ca-Ib/c events, we find two populations of events distinguished by their red ($g - r \approx 1.5$ mag) or green ($g - r \approx 0.5$ mag) spectral colors at $r$-band peak, wherein redder events show strong line blanketing signatures, slower light curves, weaker He lines and lower [Ca II]/[O I] in the nebular phase. Together, we find that the spectroscopic continuum, volumetric rates and striking old environments are consistent with the explosive burning of He shells on low mass white dwarfs. We posit that Ca-Ia and red Ca-Ib/c objects are consistent with the double detonation of He shells with high He burning efficiency, while green Ca-Ib/c objects could arise from less efficient He burning scenarios such as detonations in low density He shells or He shell deflagrations.

We present ZTF18abvkwla (the "Koala"), a fast blue optical transient discovered in the Zwicky Transient Facility (ZTF) One-Day Cadence (1DC) Survey. ZTF18abvkwla has a number of features in common with the groundbreaking transient AT2018cow: blue colors at peak ($g-r\approx-0.5$ mag), a short rise time from half-max of under two days, a decay time to half-max of only three days, a high optical luminosity ($M_{g,\mathrm{peak}}\approx-20.6$mag), a hot ($\gtrsim 40,000$K) featureless spectrum at peak light, and a luminous radio counterpart. At late times ($\Delta t>80$d) the radio luminosity of ZTF18abvkwla ($\nu L_\nu \gtrsim 10^{40}$erg/s at 10 GHz, observer-frame) is most similar to that of long-duration gamma-ray bursts (GRBs). The host galaxy is a dwarf starburst galaxy ($M\approx5\times10^{8}M_\odot$, $\mathrm{SFR}\approx7 M_\odot$/yr) that is moderately metal-enriched ($\log\mathrm{[O/H]} \approx 8.5$), similar to the hosts of GRBs and superluminous supernovae. As in AT2018cow, the radio and optical emission in ZTF18abvkwla likely arise from two separate components: the radio from fast-moving ejecta ($\Gamma \beta c >0.38c$) and the optical from shock-interaction with confined dense material ($<0.07M_\odot$ in $\sim 10^{15}$cm). Compiling transients in the literature with $t_\mathrm{rise} <5$d and $M_\mathrm{peak}<-20$mag, we find that a significant number are engine-powered, and suggest that the high peak optical luminosity is directly related to the presence of this engine. From 18 months of the 1DC survey, we find that transients in this rise-luminosity phase space are at least two to three orders of magnitude less common than CC SNe. Finally, we discuss strategies for identifying such events with future facilities like the Large Synoptic Survey Telescope, and prospects for detecting accompanying X-ray and radio emission.

We present ZTF18aaqjovh (SN 2018bvw), a high-velocity ("broad-lined") stripped-envelope (Type Ic) supernova (Ic-BL SN) discovered in the Zwicky Transient Facility one-day cadence survey. ZTF18aaqjovh shares a number of features in common with engine-driven explosions: the photospheric velocity and the shape of the optical light curve are very similar to that of the Type Ic-BL SN 1998bw, which was associated with a low-luminosity gamma-ray burst (LLGRB) and had relativistic ejecta. However, the radio luminosity of ZTF18aaqjovh is almost two orders of magnitude fainter than that of SN 1998bw at the same velocity phase, and the shock velocity is at most mildly relativistic (v=0.06-0.4c). A search of high-energy catalogs reveals no compelling GRB counterpart to ZTF18aaqjovh, and the limit on the prompt GRB luminosity of $L_{\gamma,\mathrm{iso}} \approx 1.6 \times 10^{48}$ erg/sec excludes a classical GRB but not an LLGRB. Altogether, ZTF18aaqjovh represents another transition event between engine-driven SNe associated with GRBs and "ordinary" Ic-BL SNe.

We report the discovery of the first short period binary in which a hot subdwarf star (sdOB) fills its Roche lobe and started mass transfer to its companion. The object was discovered as part of a dedicated high-cadence survey of the Galactic Plane named the Zwicky Transient Facility and exhibits a period of $P_{\rm orb}=39.3401(1)$ min, making it the most compact hot subdwarf binary currently known. Spectroscopic observations are consistent with an intermediate He-sdOB star with an effective temperature of $T_{\rm eff}=42,400\pm300$ K and a surface gravity of $\log(g)=5.77\pm0.05$. A high-signal-to noise GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the sdOB star and an eclipse of the sdOB by an accretion disk. We infer a low-mass hot subdwarf donor with a mass $M_{\rm sdOB}=0.337\pm0.015$ M$_\odot$ and a white dwarf accretor with a mass $M_{\rm WD}=0.545\pm0.020$ M$_\odot$. Theoretical binary modeling indicates the hot subdwarf formed during a common envelope phase when a $2.5-2.8$ M$_\odot$ star lost its envelope when crossing the Hertzsprung Gap. To match its current $P_{\rm orb}$, $T_{\rm eff}$, $\log(g)$, and masses, we estimate a post-common envelope period of $P_{\rm orb}\approx150$ min, and find the sdOB star is currently undergoing hydrogen shell burning. We estimate that the hot subdwarf will become a white dwarf with a thick helium layer of $\approx0.1$ M$_\odot$ and will merge with its carbon/oxygen white dwarf companion after $\approx17$ Myr and presumably explode as a thermonuclear supernova or form an R CrB star.

We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) in the Ultra-Violet (UV) and visible light. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical circumstellar material (CSM), we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% on the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modelling of these events.

While it is clear that Type Ia supernovae (SNe) are the result of thermonuclear explosions in C/O white dwarfs (WDs), a great deal remains uncertain about the binary companion that facilitates the explosive disruption of the WD. Here, we present a comprehensive analysis of a unique, and large, data set of 127 SNe Ia with exquisite coverage by the Zwicky Transient Facility (ZTF). High-cadence (6 observations per night) ZTF observations allow us to measure the SN rise time and examine its initial evolution. We develop a Bayesian framework to model the early rise as a power-law in time, which enables the inclusion of priors in our model. For a volume-limited subset of normal SNe Ia, we find the mean power-law index is consistent with 2 in the $r_\mathrm{ztf}$-band ($\alpha_r = 2.01\pm0.02$), as expected in the expanding fireball model. There are, however, individual SNe that are clearly inconsistent with $\alpha_r=2$. We estimate a mean rise time of 18.5$\,$d (with a range extending from $\sim$15$-$22$\,$d), though this is subject to the adopted prior. We identify an important, previously unknown, bias whereby the rise times for higher redshift SNe within a flux-limited survey are systematically underestimated. This effect can be partially alleviated if the power-law index is fixed to $\alpha=2$, in which case we estimate a mean rise time of 21.0$\,$d (with a range from $\sim$18$-$23$\,$d). The sample includes a handful or rare and peculiar SNe Ia. Finally, we conclude with a discussion of lessons learned from the ZTF sample that can eventually be applied to Large Synoptic Survey Telescope observations.

While tidal disruption events (TDEs) have long been heralded as laboratories for the study of quiescent black holes, the small number of known TDEs and uncertainties in their emission mechanism have hindered progress towards this promise. Here present 17 new TDEs that have been detected recently by the Zwicky Transient Facility along with Swift UV and X-ray follow-up observations. Our homogeneous analysis of the optical/UV light curves, including 22 previously known TDEs from the literature, reveals a clean separation of light curve properties with spectroscopic class. The TDEs with Bowen fluorescence features in their optical spectra have smaller blackbody radii, as well as longer rise times and higher disruption rates compared to the rest of the sample. The Bowen fluorescence mechanism requires a high density which can be reached at smaller radii, which in turn yields longer diffusion timescales. Thus, the difference in rise times suggests the pre-peak TDE light curves are governed not by the fallback timescale, but instead by the diffusion of photons through the tidal debris. The small subset of TDEs that show only helium emission lines in their spectra have the longest rise times, the highest luminosities and the lowest rates. We also report, for the first time, the detection of soft X-ray flares from a TDE on day timescales. Based on the fact the flares peak at a luminosity similar to the optical/UV blackbody luminosity, we attribute them to brief glimpses through a reprocessing layer that otherwise obscures the inner accretion flow.

We present detailed observations of ZTF18abukavn (SN2018gep), discovered in high-cadence data from the Zwicky Transient Facility as a rapidly rising ($1.4\pm0.1$ mag/hr) and luminous ($M_{g,\mathrm{peak}}=-20$ mag) transient. It is spectroscopically classified as a broad-lined stripped-envelope supernova (Ic-BL SN). The high peak luminosity ($L_{\mathrm{bol}} \gtrsim 3 \times 10^{44}$ erg $\mathrm{sec}^{-1}$), the short rise time ($t_{\mathrm{rise}}= 3$ days in $g$-band), and the blue colors at peak ($g-r\sim-0.4$) all resemble the high-redshift Ic-BL iPTF16asu, as well as several other unclassified fast transients. The early discovery of SN2018gep (within an hour of shock breakout) enabled an intensive spectroscopic campaign, including the highest-temperature ($T_{\mathrm{eff}}\gtrsim40,000$ K) spectra of a stripped-envelope SN. A retrospective search revealed luminous ($M_g \sim M_r \approx -14\,$mag) emission in the days to weeks before explosion, the first definitive detection of precursor emission for a Ic-BL. We find a limit on the isotropic gamma-ray energy release $E_\mathrm{\gamma,iso}<4.9 \times 10^{48}$ erg, a limit on X-ray emission $L_{\mathrm{X}} < 10^{40}\,$erg sec$^{-1}$, and a limit on radio emission $\nu L_\nu \lesssim 10^{37}\,$erg sec$^{-1}$. Taken together, we find that the early ($<10\,$days) data are best explained by shock breakout in a massive shell of dense circumstellar material ($0.02\,M_\odot$) at large radii ($3 \times 10^{14}\,$cm) that was ejected in eruptive pre-explosion mass-loss episodes. The late-time ($>10$ days) light curve requires an additional energy source, which could be the radioactive decay of Ni-56.

Comets are primitive objects that formed in the protoplanetary disk, and have been largely preserved over the history of the Solar System. However, they are not pristine, and surfaces of cometary nuclei do evolve. In order to understand the extent of their primitive nature, we must define the mechanisms that affect their surfaces and comae. We examine the lightcurve of comet 240P/NEAT over three consecutive orbits, and investigate three events of significant brightening ($\Delta m \sim -2$ mag). Unlike typical cometary outbursts, each of the three events are long-lived, with enhanced activity for at least 3 to 6 months. The third event, observed by the Zwicky Transient Facility, occurred in at least two stages. The anomalous behavior appears to have started after the comet was perturbed by Jupiter in 2007, reducing its perihelion distance from 2.53 to 2.12 au. We suggest that the brightening events are temporary transitions to a higher baseline activity level, brought on by the increased insolation, which has warmed previously insulated sub-surface layers. The new activity is isolated to one or two locations on the nucleus, indicating that the surface or immediate sub-surface is heterogeneous. Further study of this phenomenon may provide insight into cometary outbursts, the structure of the near-surface nucleus, and cometary nucleus mantling.

The Zwicky Transient Facility (ZTF) is performing a three-day cadence survey of the visible Northern sky (~3$\pi$). The transient candidates found in this survey are announced via public alerts. As a supplementary product ZTF is also conducting a large spectroscopic campaign: the ZTF Bright Transient Survey (BTS). The goal of the BTS is to spectroscopically classify all extragalactic transients brighter than 18.5 mag at peak brightness and immediately announce those classifications to the public. Extragalactic discoveries from ZTF are predominantly Supernovae (SNe). The BTS is the largest flux-limited SN survey to date. Here we present a catalog of the761 SNe that were classified during the first nine months of the survey (2018 Apr. 1 to 2018 Dec. 31). The BTS SN catalog contains redshifts based on SN template matching and spectroscopic host galaxy redshifts when available. Based on this data we perform an analysis of the redshift completeness of local galaxy catalogs, dubbed as the Redshift Completeness Fraction (RCF; the number of SN host galaxies with known spectroscopic redshift prior to SN discovery divided by the total number of SN hosts). In total, we identify the host galaxies of 512 Type Ia supernovae, 227 of which have known spectroscopic redshifts, yielding an RCF estimate of $44\% \pm1\%$. We find a steady decrease in the RCF with increasing distance in the local universe. For z<0.05, or ~200 Mpc, we find RCF=0.6, which has important ramifications when searching for multimessenger astronomical events. Prospects for dramatically increasing the RCF are limited to new multi-fiber spectroscopic instruments, or wide-field narrowband surveys. We find that existing galaxy redshift catalogs are only $50\%$ complete at $r\approx16.9$ mag. Pushing this limit several magnitudes deeper will pay huge dividends when searching for electromagnetic counterparts to gravitational wave events.

Early-time observations of Type Ia supernovae (SNe Ia) are essential to constrain their progenitor properties. In this paper, we present high-quality light curves of 127 SNe Ia discovered by the Zwicky Transient Facility (ZTF) in 2018. We describe our method to perform forced point spread function (PSF) photometry, which can be applied to other types of extragalactic transients. With a planned cadence of six observations per night ($3g+3r$), all of the 127 SNe Ia are detected in both $g$ and $r$ band more than 10\,d (in the rest frame) prior to the epoch of $g$-band maximum light. The redshifts of these objects range from $z=0.0181$ to 0.165; the median redshift is 0.074. Among the 127 SNe, 50 are detected at least 14\,d prior to maximum light (in the rest frame), with a subset of 9 objects being detected more than 17\,d before $g$-band peak. This is the largest sample of young SNe Ia collected to date; it can be used to study the shape and color evolution of the rising light curves in unprecedented detail. We discuss six peculiar events in this sample, including one 02cx-like event ZTF18abclfee (SN\,2018crl), one Ia-CSM SN ZTF18aaykjei (SN\,2018cxk), and four objects with possible super-Chandrasekhar mass progenitors: ZTF18abhpgje (SN\,2018eul), ZTF18abdpvnd (SN\,2018dvf), ZTF18aawpcel (SN\,2018cir) and ZTF18abddmrf (SN\,2018dsx).

We present photometry and spectroscopy of four hydrogen-poor luminous supernovae discovered during the two-month science commissioning and early operations of the Zwicky Transient Facility (ZTF) survey. Three of these objects, SN2018bym (ZTF18aapgrxo), SN2018avk (ZTF18aaisyyp) and SN2018bgv (ZTF18aavrmcg) resemble typical SLSN-I spectroscopically, while SN2018don (ZTF18aajqcue) may be an object similar to SN2007bi experiencing considerable host galaxy reddening, or an intrinsically long-lived, luminous and red SN Ic. We analyze the light curves, spectra, and host galaxy properties of these four objects and put them in context of the population of SLSN-I. SN2018bgv stands out as the fastest-rising SLSN-I observed to date, with a rest-frame g-band rise time of just 10 days from explosion to peak -- if it is powered by magnetar spin-down, the implied ejecta mass is only ~1 M$_{\odot}$. SN2018don also displays unusual properties -- in addition to its red colors and comparatively massive host galaxy, the light curve undergoes some of the strongest light curve undulations post-peak seen in a SLSN-I, which we speculate may be due to interaction with circumstellar material. We discuss the promises and challenges of finding SLSNe in large-scale surveys like ZTF given the observed diversity in the population.

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

We report the discovery of six active galactic nuclei (AGN) caught "turning on" during the first nine months of the Zwicky Transient Facility (ZTF) survey. The host galaxies were classified as LINERs by weak narrow forbidden line emission in their archival SDSS spectra, and detected by ZTF as nuclear transients. In five of the cases, we found via follow-up spectroscopy that they had transformed into broad-line AGN, reminiscent of the changing-look LINER iPTF 16bco. In one case, ZTF18aajupnt/AT2018dyk, follow-up HST UV and ground-based optical spectra revealed the transformation into a narrow-line Seyfert 1 (NLS1) with strong [Fe VII, X, XIV] and He II 4686 coronal lines. Swift monitoring observations of this source reveal bright UV emission that tracks the optical flare, accompanied by a luminous soft X-ray flare that peaks ~60 days later. Spitzer follow-up observations also detect a luminous mid-infrared flare implying a large covering fraction of dust. Archival light curves of the entire sample from CRTS, ATLAS, and ASAS-SN constrain the onset of the optical nuclear flaring from a prolonged quiescent state. Here we present the systematic selection and follow-up of this new class of changing-look LINERs, compare their properties to previously reported changing-look Seyfert galaxies, and conclude that they are a unique class of transients well-suited to test the uncertain physical processes associated with the LINER accretion state.

We report the discovery of a mid-infrared variable AGN which is hosted by an ultraluminous infrared galaxy (ULIRG) candidate in the Sloan Stripe 82 field. \textit{WISE} \textit{J030654.88+010833.6} is a red, extended galaxy, which we estimate to be at a photometric redshift of 0.28 $\leq$ z $\leq$ 0.31, based on its optical and near-infrared spectral energy distribution (SED). The factor of two variability over 8 years seen in the \textit{WISE} 3.4 and 4.6 $\mu$m wavelength channels is not clearly correlated with optical variability in archival data. Based on our estimation of the physical parameters of the host galaxy, \textit{J030654.88+010833.6} is possibly a composite AGN/starburst ULIRG in a phase where high star formation $\sim$ 70 M$_{\odot}$ year$^{-1}$ is occurring. Our estimate of the black hole mass to stellar mass ratio also appears to be consistent with that of broad-line AGN in the local universe. The long-term variability of \textit{J030654.88+010833.6} as seen in the \textit{WISE} \textit{W1} and \textit{W2} light curves is likely due to variations in the accretion rate, with the energy being reprocessed by dust in the vicinity of the AGN.

General relativity predicts that short orbital period binaries emit significant gravitational radiation, and the upcoming Laser Interferometer Space Antenna (LISA) is expected to detect tens of thousands of such systems; however, few have been identified, and only one is eclipsing--the double white dwarf binary SDSS J065133.338+284423.37, which has an orbital period of 12.75 minutes. Here, we report the discovery of an eclipsing double white dwarf binary system with an orbital period of only 6.91 minutes, ZTF J153932.16+502738.8. This system has an orbital period close to half that of SDSS J065133.338+284423.37 and an orbit so compact that the entire binary could fit within the diameter of the planet Saturn. The system exhibits a deep eclipse, and a double-lined spectroscopic nature. We observe rapid orbital decay, consistent with that expected from general relativity. ZTF J153932.16+502738.8 is a significant source of gravitational radiation close to the peak of LISA's sensitivity, and should be detected within the first week of LISA observations.

We explore methods for the identification of stellar flare events in irregularly sampled data of ground-based time domain surveys. In particular, we describe a new technique for identifying flaring stars, which we have implemented in a publicly available Python module called "PyVAN". The approach uses the Differential Evolution algorithm to optimize parameters of empirically derived light-curve templates for different types of stars to fit a candidate light-curve. The difference of the likelihoods that these best-fit templates produced the observed data is then used to delineate targets that are well explained by a flare template but simultaneously poorly explained by templates of common contaminants. By testing on light-curves of known identity and morphology, we show that our technique is capable of recovering flaring status in $69\%$ of all light-curves containing a flare event above thresholds drawn to include $\lt1\%$ of any contaminant population. By applying to Palomar Transient Factory data, we show consistency with prior samples of flaring stars, and identify a small selection of candidate flaring G-type stars for possible follow-up.

Using high-cadence observations from the Zwicky Transient Facility at low Galactic latitudes, we have discovered a new class of pulsating, hot, compact stars. We have found four candidates, exhibiting blue colors ($g-r\leq-0.1$ mag), pulsation amplitudes of $>5\%$, and pulsation periods of $200 - 475$ sec. Fourier transforms of the lightcurves show only one dominant frequency. Phase-resolved spectroscopy for three objects reveals significant radial velocity, T$_{\rm eff}$ and log(g) variations over the pulsation cycle, consistent with large amplitude radial oscillations. The mean T$_{\rm eff}$ and log(g) for these stars are consistent with hot subdwarf B (sdB) effective temperatures and surface gravities. We calculate evolutionary tracks using MESA and adiabatic pulsations using GYRE for low-mass helium-core pre-white dwarfs and low mass helium-burning stars. Comparison of low-order radial oscillation mode periods with the observed pulsation periods show better agreement with the pre-white dwarf models. Therefore, we suggest that these new pulsators and Blue Large-Amplitude Pulsators (BLAPs) could be members of the same class of pulsators, composed of young $\approx0.25-0.35$ M$_\odot$ helium-core pre-white dwarfs.

We present the Census of the Local Universe (CLU) narrow-band survey to search for emission-line (\ha) galaxies. CLU-\ha~has imaged $\approx$3$\pi$ of the sky (26,470~deg$^2$) with 4 narrow-band filters that probe a distance out to 200~Mpc. We have obtained spectroscopic follow-up for galaxy candidates in 14 preliminary fields (101.6~deg$^2$) to characterize the limits and completeness of the survey. In these preliminary fields, CLU can identify emission lines down to an \ha~flux limit of $10^{-14}$~$\rm{erg~s^{-1}~cm^{-2}}$ at 90\% completeness, and recovers 83\% (67\%) of the \ha~flux from catalogued galaxies in our search volume at the $\Sigma$=2.5 ($\Sigma$=5) color excess levels. The contamination from galaxies with no emission lines is 61\% (12\%) for $\Sigma$=2.5 ($\Sigma$=5). Also, in the regions of overlap between our preliminary fields and previous emission-line surveys, we recover the majority of the galaxies found in previous surveys and identify an additional $\approx$300 galaxies. In total, we find 90 galaxies with no previous distance information, several of which are interesting objects: 7 blue compact dwarfs, 1 green pea, and a Seyfert galaxy; we also identified a known planetary nebula. These objects show that the CLU-\ha~survey can be a discovery machine for objects in our own Galaxy and extreme galaxies out to intermediate redshifts. However, the majority of the CLU-\ha~galaxies identified in this work show properties consistent with normal star-forming galaxies. CLU-\ha~galaxies with new redshifts will be added to existing galaxy catalogs to focus the search for the electromagnetic counterpart to gravitational wave events.

We describe ZStreak, a semi-real-time pipeline specialized in detecting small, fast-moving near-Earth asteroids (NEAs) that is currently operating on the data from the newly-commissioned Zwicky Transient Facility (ZTF) survey. Based on a prototype originally developed by Waszczak et al. (2017) for the Palomar Transient Factory (PTF), the predecessor of ZTF, ZStreak features an improved machine-learning model that can cope with the $10\times$ data rate increment between PTF and ZTF. Since its first discovery on 2018 February 5 (2018 CL), ZTF/ZStreak has discovered $45$ confirmed new NEAs over a total of 232 observable nights until 2018 December 31. Most of the discoveries are small NEAs, with diameters less than $\sim100$ m. By analyzing the discovery circumstances, we find that objects having the first to last detection time interval under 2 hr are at risk of being lost. We will further improve real-time follow-up capabilities, and work on suppressing false positives using deep learning.

We report the discovery of non-stellar hydrogen Balmer and metastable helium absorption lines accompanying a transient, high-velocity (0.05$c$) broad absorption line (BAL) system in the optical spectra of the tidal disruption event (TDE) AT2018zr ($z=0.071$). In the HST UV spectra, absorption of high- and low-ionization lines are also present at this velocity, making AT2018zr resemble a low-ionization broad absorption line (LoBAL) QSO. We conclude that these transient absorption features are more likely to arise in fast outflows produced by the TDE than absorbed by the unbound debris. In accordance with the outflow picture, we are able to reproduce the flat-topped H$\alpha$ emission in a spherically expanding medium, without invoking the typical prescription of an elliptical disk. We also report the appearance of narrow ($\sim$1000~km~s$^{-1}$) NIII$\lambda$4640, HeII$\lambda4686$, H$\alpha$, and H$\beta$, emission in the late-time optical spectra of AT2018zr, which may be a result of UV continuum hardening at late time as observed by Swift. Including AT2018zr, we find a high association rate (3 out of 4) of BALs in the UV spectra of TDEs. This suggests that outflows may be ubiquitous among TDEs and may be less sensitive to viewing angle effects compared to QSO outflows.

The Zwicky Transient Facility is a large optical survey in multiple filters producing hundreds of thousands of transient alerts per night. We describe here various machine learning (ML) implementations and plans to make the maximal use of the large data set by taking advantage of the temporal nature of the data, and further combining it with other data sets. We start with the initial steps of separating bogus candidates from real ones, separating stars and galaxies, and go on to the classification of real objects into various classes. Besides the usual methods (e.g., based on features extracted from light curves) we also describe early plans for alternate methods including the use of domain adaptation, and deep learning. In a similar fashion we describe efforts to detect fast moving asteroids. We also describe the use of the Zooniverse platform for helping with classifications through the creation of training samples, and active learning. Finally we mention the synergistic aspects of ZTF and LSST from the ML perspective.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ~3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5sigma in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The realtime pipeline employs a novel image-differencing algorithm, optimized for the detection of point source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. The reconstructed astrometric accuracy per science image with respect to Gaia is typically 45 to 85 milliarcsec. This is the RMS per axis on the sky for sources extracted with photometric S/N >= 10. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae, rare flux transients, variable stars, eclipsing binaries, variability from active galactic nuclei, counterparts to gravitational wave sources, a more complete census of Type Ia supernovae, and Solar System objects.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

The detonation of a helium shell on a white dwarf has been proposed as a possible explosion triggering mechanism for Type Ia supernovae. Here, we report ZTF 18aaqeasu (SN 2018byg/ATLAS 18pqq), a peculiar Type I supernova, consistent with being a helium-shell double-detonation. With a rise time of $\approx 18$ days from explosion, the transient reached a peak absolute magnitude of $M_R \approx -18.2$ mag, exhibiting a light curve akin to sub-luminous SN 1991bg-like Type Ia supernovae, albeit with an unusually steep increase in brightness within a week from explosion. Spectra taken near peak light exhibit prominent Si absorption features together with an unusually red color ($g-r \approx 2$ mag) arising from nearly complete line blanketing of flux blue-wards of 5000 \AA. This behavior is unlike any previously observed thermonuclear transient. Nebular phase spectra taken at and after $\approx 30$ days from peak light reveal evidence of a thermonuclear detonation event dominated by Fe-group nucleosynthesis. We show that the peculiar properties of ZTF 18aaqeasu are consistent with the detonation of a massive ($\approx 0.15$ M$_\odot$) helium shell on a sub-Chandrasekhar mass ($\approx 0.75$ M$_\odot$) white dwarf after including mixing of $\approx 0.2$ M$_\odot$ of material in the outer ejecta. These observations provide evidence of a likely rare class of thermonuclear supernovae arising from detonations of massive helium shells.

Although dwarf carbon (dC) stars are thought universally to be binaries to explain the presence of $C_2$ in their spectra while still near main sequence luminosity, direct observational evidence for binarity is remarkably scarce. Here we report the detection of a 2.92 d periodicity in both photometry and radial velocity of SDSS J125017.90+252427.6, an $r=16.4$ dC star. This is the first photometric binary dC, and only the second dC spectroscopic binary. The relative phase of the photometric period to the spectroscopic observations suggests that the photometric variations are a reflection effect due to heating from an unseen companion. The observed radial velocity amplitude of the dC component ($K = 98.8\pm10.7$ km/s) is consistent with a white dwarf companion, presumably the evolved star that earlier donated the carbon to the dC, although substantial orbital evolution must have occurred. Large synoptic photometric surveys such as the Palomar Transient Factory, used for this work, may prove useful for identifying binaries among the shorter period dC stars.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

We report the discovery and analysis of 36 new eclipsing EL CVn-type binaries, consisting of a core helium-composition pre-white dwarf and an early-type main-sequence companion, more than doubling the known population of these systems. We have used supervised machine learning methods to search 0.8 million lightcurves from the Palomar Transient Factory, combined with SDSS, Pan-STARRS and 2MASS colours. The new systems range in orbital periods from 0.46-3.8 d and in apparent brightness from ~14-16 mag in the PTF $R$ or $g^{\prime}$ filters. For twelve of the systems, we obtained radial velocity curves with the Intermediate Dispersion Spectrograph at the Isaac Newton Telescope. We modelled the lightcurves, radial velocity curves and spectral energy distributions to determine the system parameters. The radii (0.3-0.7 $\mathrm{R_{\odot}}$) and effective temperatures (8000-17000 K) of the pre-He-WDs are consistent with stellar evolution models, but the masses (0.12-0.28 $\mathrm{M_{\odot}}$) show more variance than models predicted. This study shows that using machine learning techniques on large synoptic survey data is a powerful way to discover substantial samples of binary systems in short-lived evolutionary stages.

The Zwicky Transient Facility is a new robotic-observing program, in which a newly engineered 600-MP digital camera with a pioneeringly large field of view, 47~square degrees, will be installed into the 48-inch Samuel Oschin Telescope at the Palomar Observatory. The camera will generate $\sim 1$~petabyte of raw image data over three years of operations. In parallel related work, new hardware and software systems are being developed to process these data in real time and build a long-term archive for the processed products. The first public release of archived products is planned for early 2019, which will include processed images and astronomical-source catalogs of the northern sky in the $g$ and $r$ bands. Source catalogs based on two different methods will be generated for the archive: aperture photometry and point-spread-function fitting.

The Spectral and Photometric Imaging Receiver (SPIRE) was launched as one of the scientific instruments on board of the space observatory Herschel. The SPIRE photometer opened up an entirely new window in the Submillimeter domain for large scale mapping, that up to then was very difficult to observe. There are already several catalogs that were produced by individual Herschel science projects. Yet, we estimate that the objects of only a fraction of these maps will ever be systematically extracted and published by the science teams that originally proposed the observations. The SPIRE instrument performed its standard photometric observations in an optically very stable configuration, only moving the telescope across the sky, with variations in its configuration parameters limited to scan speed and sampling rate. This and the scarcity of features in the data that require special processing steps made this dataset very attractive for producing an expert reduced catalog of point sources that is being described in this document. The Catalog was extracted from a total of 6878 unmodified SPIRE scan map observations. The photometry was obtained by a systematic and homogeneous source extraction procedure, followed by a rigorous quality check that emphasized reliability over completeness. Having to exclude regions affected by strong Galactic emission, that pushed the limits of the four source extraction methods that were used, this catalog is aimed primarily at the extragalactic community. The result can serve as a pathfinder for ALMA and other Submillimeter and Far-Infrared facilities. 1,693,718 sources are included in the final catalog, splitting into 950688, 524734, 218296 objects for the 250\mu m, 350\mu m, and 500\mu m bands, respectively. The catalog comes with well characterized environments, reliability, completeness, and accuracies, that single programs typically cannot provide.

We report the discovery and analysis of PTF1 J085713+331843, a new eclipsing post common-envelope detached white-dwarf red-dwarf binary with a 2.5h orbital period discovered by the Palomar Transient Factory. ULTRACAM multicolour photometry over multiple orbital periods reveals a light curve with a deep flat-bottomed primary eclipse and a strong reflection effect. Phase-resolved spectroscopy shows broad Balmer absorption lines from the DA white dwarf and phase-dependent Balmer emission lines originating on the irradiated side of the red dwarf. The temperature of the DA white dwarf is $T_\mathrm{WD} = 25700 \pm 400\,$K and the spectral type of the red dwarf is M3-5. A combined modelling of the light curve and the radial velocity variations results in a white dwarf mass of $M_\mathrm{WD} = 0.61^{+0.18}_{-0.17}\, \mathrm{M_{\odot}}$ and radius of $R_\mathrm{WD} = 0.0175^{+0.0012}_{-0.0011}\, \mathrm{R_{\odot}}$, and a red dwarf mass and radius of $M_\mathrm{RD} = 0.19^{+0.10}_{-0.08}\, \mathrm{M_{\odot}}$ and $R_\mathrm{RD} = 0.24^{+0.04}_{-0.04}\, \mathrm{R_{\odot}}$. The system is either a detached cataclysmic variable or has emerged like from the common envelope phase at nearly its current orbital period. In $\sim70\,$Myr, this system will become a cataclysmic variable in the period gap.

We present the results of a Spitzer/Herschel infrared photometric analysis of the largest (716) and highest-redshift (z=1.8) sample of Brightest Cluster Galaxies (BCGs), those from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS). Given the tension that exists between model predictions and recent observations of BCGs at z<2, we aim to uncover the dominant physical mechanism(s) guiding the stellar-mass buildup of this special class of galaxies, the most massive in the Universe uniquely residing at the centres of galaxy clusters. Through a comparison of their stacked, broadband, infrared spectral energy distributions (SEDs) to a variety of SED model templates in the literature, we identify the major sources of their infrared energy output, in multiple redshift bins between 0 < z < 1.8. We derive estimates of various BCG physical parameters from the stacked {\nu}L{\nu} SEDs, from which we infer a star-forming, as opposed to a 'red and dead' population of galaxies, producing tens to hundreds of solar masses per year down to z=0.5. This discovery challenges the accepted belief that BCGs should only passively evolve through a series of gas-poor, minor mergers since z~4 (De Lucia & Blaizot 2007), but agrees with the improved semi-analytic model of hierarchical structure formation of Tonini et al. (2012), which predicts star-forming BCGs throughout the epoch considered. We attribute the star formation inferred from the stacked infrared SEDs to both major and minor 'wet' (gas-rich) mergers, based on a lack of key signatures (to date) of the cluster cooling flows to which BCG star formation is typically attributed, as well as a number of observational and simulation-based studies that support this scenario.

Far-infrared (FIR) images and photometry are presented for 201 Luminous and Ultraluminous Infrared Galaxies [LIRGs: log$(L_{\rm IR}/L_\odot) = 11.00 - 11.99$, ULIRGs: log$(L_{\rm IR}/L_\odot) = 12.00 - 12.99$], in the Great Observatories All-Sky LIRG Survey (GOALS) based on observations with the $Herschel$ $Space$ $Observatory$ Photodetector Array Camera and Spectrometer (PACS) and the Spectral and Photometric Imaging Receiver (SPIRE) instruments. The image atlas displays each GOALS target in the three PACS bands (70, 100, and 160 $\mu$m) and the three SPIRE bands (250, 350, and 500 $\mu$m), optimized to reveal structures at both high and low surface brightness levels, with images scaled to simplify comparison of structures in the same physical areas of $\sim$$100\times100$ kpc$^2$. Flux densities of companion galaxies in merging systems are provided where possible, depending on their angular separation and the spatial resolution in each passband, along with integrated system fluxes (sum of components). This dataset constitutes the imaging and photometric component of the GOALS Herschel OT1 observing program, and is complementary to atlases presented for the Hubble Space Telescope (Evans et al. 2017, in prep.), Spitzer Space Telescope (Mazzarella et al. 2017, in prep.), and Chandra X-ray Observatory (Iwasawa et al. 2011, 2017, in prep.). Collectively these data will enable a wide range of detailed studies of AGN and starburst activity within the most luminous infrared galaxies in the local Universe.

We present the discovery of the hot subdwarf B star (sdB) binary PTF1 J082340.04+081936.5. The system has an orbital period P$_{\rm orb}=87.49668(1)$ min (0.060761584(10) days), making it the second-most compact sdB binary known. The lightcurve shows ellipsoidal variations. Under the assumption that the sdB primary is synchronized with the orbit, we find a mass $M_{\rm sdB}=0.45^{+0.09}_{-0.07}$ M$_\odot$, a companion white dwarf mass $M_{\rm WD}=0.46^{+0.12}_{-0.09}$ M$_\odot$ and a mass ratio $q = \frac{M_{\rm WD}}{M_{\rm sdB}}=1.03^{+0.10}_{-0.08}$. The future evolution was calculated using the MESA stellar evolution code. Adopting a canonical sdB mass of $M_{\rm sdB}=0.47$ M$_\odot$, we find that the sdB still burns helium at the time it will fill its Roche lobe if the orbital period was less than 106 min at the exit from the last common envelope phase. For longer common envelope exit periods the sdB will have stopped burning helium and turned into a C/O white dwarf at the time of contact. Comparing the spectroscopically derived log(g) and $T_{\rm eff}$ with our MESA models, we find that an sdB model with a hydrogen envelope mass of $5\times10^{-4} M_\odot$ matches the measurements at a post-common envelope age of 94 Myr, corresponding to a post-common envelope orbital period of 109 min which is close to the limit to start accretion while the sdB is still burning helium.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We describe the near real-time transient-source discovery engine for the intermediate Palomar Transient Factory (iPTF), currently in operations at the Infrared Processing and Analysis Center (IPAC), Caltech. We coin this system the IPAC/iPTF Discovery Engine (or IDE). We review the algorithms used for PSF-matching, image subtraction, detection, photometry, and machine-learned (ML) vetting of extracted transient candidates. We also review the performance of our ML classifier. For a limiting signal-to-noise ratio of 4 in relatively unconfused regions, "bogus" candidates from processing artifacts and imperfect image subtractions outnumber real transients by ~ 10:1. This can be considerably higher for image data with inaccurate astrometric and/or PSF-matching solutions. Despite this occasionally high contamination rate, the ML classifier is able to identify real transients with an efficiency (or completeness) of ~ 97% for a maximum tolerable false-positive rate of 1% when classifying raw candidates. All subtraction-image metrics, source features, ML probability-based real-bogus scores, contextual metadata from other surveys, and possible associations with known Solar System objects are stored in a relational database for retrieval by the various science working groups. We review our efforts in mitigating false-positives and our experience in optimizing the overall system in response to the multitude of science projects underway with iPTF.

We used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,\mu$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,\mu$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $\Phi_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $\Phi_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,\mu$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.

We present a five-band Herschel study (100-500um) of three galaxy clusters at z~1.2 from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS). With a sample of 120 spectroscopically-confirmed cluster members, we investigate the role of environment on galaxy properties utilizing the projected cluster phase space (line-of-sight velocity versus clustercentric radius), which probes the time-averaged galaxy density to which a galaxy has been exposed. We divide cluster galaxies into phase-space bins of (r/r200) x (v/sigma_v), tracing a sequence of accretion histories in phase space. Stacking optically star-forming cluster members on the Herschel maps, we measure average infrared star formation rates, and, for the first time in high-redshift galaxy clusters, dust temperatures for dynamically distinct galaxy populations---namely, recent infalls and those that were accreted onto the cluster at an earlier epoch. Proceeding from the infalling to virialized (central) regions of phase space, we find a steady decrease in the specific star formation rate and increase in the stellar age of star-forming cluster galaxies. We perform a probability analysis to investigate all acceptable infrared spectral energy distributions within the full parameter space and measure a ~4 sigma drop in the average dust temperature of cluster galaxies in an intermediate phase-space bin, compared to an otherwise flat trend with phase space. We suggest one plausible quenching mechanism which may be consistent with these trends, invoking ram-pressure stripping of the warmer dust for galaxies within this intermediate accretion phase.

The Herschel Space Observatory was the fourth cornerstone mission in the European Space Agency (ESA) science programme with excellent broad band imaging capabilities in the sub-mm and far-infrared part of the spectrum. Although the spacecraft finished its observations in 2013, it left a large legacy dataset that is far from having been fully scrutinised and still has a large potential for new scientific discoveries. This is specifically true for the photometric observations of the PACS and SPIRE instruments. Some source catalogues have already been produced by individual observing programs, but there are many observations that risk to remain unexplored. To maximise the science return of the SPIRE and PACS data sets, we are in the process of building the Herschel Point Source Catalogue (HPSC) from all primary and parallel mode observations. Our homogeneous source extraction enables a systematic and unbiased comparison of sensitivity across the different Herschel fields that single programs will generally not be able to provide. The catalogue will be made available online through archives like the Herschel Science Archive (HSA), the Infrared Science Archive (IRSA), and the Strasbourg Astronomical Data Center (CDS).

We present the results of a MIPS-24um study of the Brightest Cluster Galaxies (BCGs) of 535 high-redshift galaxy clusters. The clusters are drawn from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS), which effectively provides a sample selected on total stellar mass, over 0.2 < z < 1.8 within the Spitzer Wide-Area Infrared Extragalactic (SWIRE) Survey fields. 20%, or 106 clusters have spectroscopically confirmed redshifts, and the rest have redshifts estimated from the color of their red sequence. A comparison with the public SWIRE images detects 125 individual BCGs at 24um > 100uJy, or 23%. The luminosity-limited detection rate of BCGs in similar richness clusters (Ngal> 12) increases rapidly with redshift. Above z ~ 1, an average of ~20\% of the sample have 24um-inferred infrared luminosities of LIR > 10^12 Lsun, while the fraction below z ~ 1 exhibiting such luminosities is < 1 \%. The Spitzer-IRAC colors indicate the bulk of the 24um-detected population is predominantly powered by star formation, with only 7/125 galaxies lying within the color region inhabited by Active Galactic Nuclei (AGN). Simple arguments limit the star-formation activity to several hundred million years and this may therefore be indicative of the timescale for AGN feedback to halt the star formation. Below redshift z ~ 1 there is not enough star formation to significantly contribute to the overall stellar mass of the BCG population, and therefore BCG growth is likely dominated by dry-mergers. Above z~ 1, however, the inferred star formation would double the stellar mass of the BCGs and is comparable to the mass assembly predicted by simulations through dry mergers. We cannot yet constrain the process driving the star formation for the overall sample, though a single object studied in detail is consistent with a gas-rich merger.

The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard is, however, showing its age. Developed in the late 1970s the FITS authors made a number of implementation choices for the format that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not appreciate the challenges which we would be facing today in astronomical computing. Difficulties we now face include, but are not limited to, having to address the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets and the need to capture significantly more complex metadata and data relationships. There are members of the community today who find some (or all) of these limitations unworkable, and have decided to move ahead with storing data in other formats. This reaction should be taken as a wakeup call to the FITS community to make changes in the FITS standard, or to see its usage fall. In this paper we detail some selected important problems which exist within the FITS standard today. It is not our intention to prescribe specific remedies to these issues; rather, we hope to call attention of the FITS and greater astronomical computing communities to these issues in the hopes that it will spur action to address them.

The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard, however, is showing its age. Developed in the late 1970s, the FITS authors made a number of implementation choices that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not anticipate the challenges which we are facing today in astronomical computing. Difficulties we now face include, but are not limited to, addressing the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets, and capturing significantly more complex metadata and data relationships. There are members of the community today who find some or all of these limitations unworkable, and have decided to move ahead with storing data in other formats. If this fragmentation continues, we risk abandoning the advantages of broad interoperability, and ready archivability, that the FITS format provides for astronomy. In this paper we detail some selected important problems which exist within the FITS standard today. These problems may provide insight into deeper underlying issues which reside in the format and we provide a discussion of some lessons learned. It is not our intention here to prescribe specific remedies to these issues; rather, it is to call attention of the FITS and greater astronomical computing communities to these problems in the hope that it will spur action to address them.

The photometer section of SPIRE is one of the key instruments on board of Herschel. Its legacy depends very much on how well the scanmap observations that it carried out during the Herschel mission can be converted to high quality maps. In order to have a comprehensive assessment on the current status of SPIRE map-making, as well as to provide guidance for future development of the SPIRE scan-map data reduction pipeline, we carried out a test campaign on SPIRE map-making. In this report, we present results of the tests in this campaign.

The different algorithms appropriate for point source photometry on data from the SPIRE instrument on-board the Herschel Space Observatory, within the Herschel Interactive Processing Environment (HIPE) are compared. Point source photometry of a large ensemble of standard calibration stars and dark sky observations is carried out using the 4 major methods within HIPE: SUSSEXtractor, DAOphot, the SPIRE Timeline Fitter and simple Aperture Photometry. Colour corrections and effective beam areas as a function of the assumed source spectral index are also included to produce a large number of photometric measurements per individual target, in each of the 3 SPIRE bands (250, 350, 500um), to examine both the accuracy and repeatability of each of the 4 algorithms. It is concluded that for flux densities down to the level of 30mJy that the SPIRE Timeline Fitter is the method of choice. However, at least in the 250 and 350um bands, all 4 methods provide photometric repeatability better than a few percent down to at approximately 100mJy. The DAOphot method appears in many cases to have a systematic offset of ~8% in all SPIRE bands which may be indicative of a sub-optimal aperture correction. In general, aperture photometry is the least reliable method, i.e. largest scatter between observations, especially in the longest wavelength band. At the faintest fluxes, <30mJy, SUSSEXtractor or DAOphot provide a better alternative to the Timeline Fitter.

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

We have observed four massive galaxy clusters with the SPIRE instrument on the Herschel Space Observatory and measure a deficit of surface brightness within their central region after subtracting sources. We simulate the effects of instrumental sensitivity and resolution, the source population, and the lensing effect of the clusters to estimate the shape and amplitude of the deficit. The amplitude of the central deficit is a strong function of the surface density and flux distribution of the background sources. We find that for the current best fitting faint end number counts, and excellent lensing models, the most likely amplitude of the central deficit is the full intensity of the cosmic infrared background (CIB). Our measurement leads to a lower limit to the integrated total intensity of the CIB of I(250 microns) > 0.69_(-0.03)^(+0.03) (stat.)_(-0.06)^(+0.11) (sys.) MJy/sr, with more CIB possible from both low-redshift sources and from sources within the target clusters. It should be possible to observe this effect in existing high angular resolution data at other wavelengths where the CIB is bright, which would allow tests of models of the faint source component of the CIB.

We describe the procedure used to flux calibrate the three-band submillimetre photometer in the Spectral and Photometric Imaging REceiver (SPIRE) instrument on the Herschel Space Observatory. This includes the equations describing the calibration scheme, a justification for using Neptune as the primary calibration source, a description of the observations and data processing procedures used to derive flux calibration parameters (for converting from voltage to flux density) for every bolometer in each array, an analysis of the error budget in the flux calibration for the individual bolometers, and tests of the flux calibration on observations of primary and secondary calibrators. The procedure for deriving the flux calibration parameters is divided into two parts. In the first part, we use observations of astronomical sources in conjunction with the operation of the photometer internal calibration source to derive the unscaled derivatives of the flux calibration curves. To scale the calibration curves in Jy/beam/V, we then use observations of Neptune in which the beam of each bolometer is mapped using Neptune observed in a very fine scan pattern. The total instrumental uncertainties in the flux calibration for the individual bolometers is ~0.5% for most bolometers, although a few bolometers have uncertainties of ~1-5% because of issues with the Neptune observations. Based on application of the flux calibration parameters to Neptune observations performed using typical scan map observing modes, we determined that measurements from each array as a whole have instrumental uncertainties of 1.5%. This is considerably less than the absolute calibration uncertainty associated with the model of Neptune, which is estimated at 4%.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

The infrared (IR) emission of M_* galaxies (10^{10.4} < M_{star} < 10^{11.0} M_\sun) in galaxy pairs, derived using data obtained in Herschel (PEP/HerMES) and Spitzer (S-COSMOS) surveys, is compared to that of single disk galaxies in well matched control samples to study the cosmic evolution of the star-formation enhancement induced by galaxy-galaxy interaction. Both the mean IR SED and mean IR luminosity of star-forming galaxies (SFGs) in SFG+SFG (S+S) pairs in the redshift bin of 0.6 < z < 1 are consistent with no star-formation enhancement. SFGs in S+S pairs in a lower redshift bin of 0.2 < z < 0.6 show marginal evidence for a weak star-formation enhancement. Together with the significant and strong sSFR enhancement shown by SFGs in a local sample of S+S pairs (obtained using previously published Spitzer observations), our results reveal a trend for the star-formation enhancement in S+S pairs to decrease with increasing redshift. Between z=0 and z=1, this decline of interaction-induced star-formation enhancement occurs in parallel with the dramatic increase (by a factor of ~10) of the sSFR of single SFGs, both can be explained by the higher gas fraction in higher z disks. SFGs in mixed pairs (S+E pairs) do not show any significant star-formation enhancement at any redshift. The difference between SFGs in S+S pairs and in S+E pairs suggests a modulation of the sSFR by the inter-galactic medium IGM in the dark matter halos (DMH) hosting these pairs.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

We use the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) to explore the specific star-formation activity of galaxies and their evolution near the peak of the cosmic far-infrared (FIR) background at 70 and 160um. We use a stacking analysis to determine the mean FIR properties of well defined subsets of galaxies at flux levels well below the FIR catalogue detection limits of SWIRE and other Spitzer surveys. We tabulate the contribution of different subsets of galaxies to the FIR background at 70um and 160um. These long wavelengths provide a good constraint on the bolometric, obscured emission. The large area provides good constraints at low z and in finer redshift bins than previous work. At all redshifts we find that the specific FIR Luminosity (sLFIR) decreases with increasing mass, following a trend L_FIR/M* propto M_* ^beta with beta =-0.38\pm0.14. This is a more continuous change than expected from the {Delucia2007} semi-analytic model suggesting modifications to the feedback prescriptions. We see an increase in the sLFIR by about a factor of ~100 from 0<z<2 and find that the sLFIR evolves as (1+z)^alpha with alpha=4.4\pm0.3 for galaxies with 10.5 < log M*/Msun < 12. This is considerably steeper than the {Delucia2007} semi-analytic model (alpha \sim 2.5). When separating galaxies into early and late types on the basis of the optical/IR spectral energy distributions we find that the decrease in sLFIR with stellar mass is stronger in early type galaxies (beta ~ -0.46), while late type galaxies exhibit a flatter trend (beta \sim -0.15). The evolution is strong for both classes but stronger for the early type galaxies. The early types show a trend of decreasing strength of evolution as we move from lower to higher masses while the evolution of the late type galaxies has little dependence on stellar mass. We suggest that in late-type galaxies we are seeing a consistently declining sSFR..

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

We describe the production and verification of sky maps of the five SPIRE fields observed as part of the Herschel Multi-tiered Extragalactic Survey (HerMES) during the Science Demonstration Phase (SDP) of the Herschel mission. We have implemented an iterative map-making algorithm (SHIM; The SPIRE-HerMES Iterative Mapper) to produce high fidelity maps that preserve extended diffuse emission on the sky while exploiting the repeated observations of the same region of the sky with many detectors in multiple scan directions to minimize residual instrument noise. We specify here the SHIM algorithm and outline the various tests that were performed to determine and characterize the quality of the maps and verify that the astrometry, point source flux and power on all relevant angular scales meets the needs of the HerMES science goals. These include multiple jackknife tests, determination of the map transfer function and detailed examination of the power spectra of both sky and jackknife maps. The map transfer function is approximately unity on scales from one arcminute to one degree. Final maps (v1.0), including multiple jackknives, as well as the SHIM pipeline, have been used by the HerMES team for the production of SDP papers.

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

We want to study the mid-infrared properties and the starburst and AGN contributions, of 24um sources at z~2, through analysis of mid-infrared spectra combined with millimeter, radio, and infrared photometry. Mid-infrared spectroscopy allows us to recover accurate redshifts. A complete sample of 16 Spitzer-selected sources (ULIRGs) believed to be starbursts at z~2 ("5.8um-peakers") was selected in the (0.5 sq.deg.) J1064+56 SWIRE Lockman Hole field. These sources have S(24um)>0.5mJy, a stellar emission peak redshifted to 5.8um, and r'(Vega)>23. The entire sample was observed with the low resolution units of the Spitzer/IRS infrared spectrograph. These sources have 1.2mm observations with IRAM 30m/MAMBO and very deep 20cm observations from the VLA. Nine of our sources also benefit from 350um observation and detection from CSO/SHARC-II. The entire sample shows good quality IRS spectra dominated by strong PAH features. The main PAH features at 6.2, 7.7, 8.6, and 11.3um have high S/N average luminosities of 2.90, 10.38, 3.62, and 2.29x10^{10}Lsun, respectively. We derived accurate redshifts spanning from 1.75 to 2.28. The average of these redshifts is 2.017. This result confirms that the selection criteria of "5.8um-peakers" associated with a strong detection at 24um are reliable to select sources at z~2. We have analyzed the different correlations between PAH emission and infrared, millimeter, and radio emission. Practically all our sources are strongly dominated by starburst emission. We have also defined two subsamples based on the equivalent width at 7.7um to investigate AGN contributions. Our sample contains strong starbursts and represents a particularly 24um-bright class of SMGs. The very good correlation between PAH and far-IR luminosities is now confirmed in high-z starburst ULIRGs. These sources show a small AGN contribution to the mid-IR, around ~20% in most cases.

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.

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

SPIRE, the Spectral and Photometric Imaging Receiver, is the Herschel Space Observatory's submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 {\mu}m, and an imaging Fourier transform spectrometer (FTS) covering 194-671 {\mu}m (447-1550 GHz). In this paper we describe the initial approach taken to the absolute calibration of the SPIRE instrument using a combination of the emission from the Herschel telescope itself and the modelled continuum emission from solar system objects and other astronomical targets. We present the photometric, spectroscopic and spatial accuracy that is obtainable in data processed through the "standard" pipelines. The overall photometric accuracy at this stage of the mission is estimated as 15% for the photometer and between 15 and 50% for the spectrometer. However, there remain issues with the photometric accuracy of the spectra of low flux sources in the longest wavelength part of the SPIRE spectrometer band. The spectrometer wavelength accuracy is determined to be better than 1/10th of the line FWHM. The astrometric accuracy in SPIRE maps is found to be 2 arcsec when the latest calibration data are used. The photometric calibration of the SPIRE instrument is currently determined by a combination of uncertainties in the model spectra of the astronomical standards and the data processing methods employed for map and spectrum calibration. Improvements in processing techniques and a better understanding of the instrument performance will lead to the final calibration accuracy of SPIRE being determined only by uncertainties in the models of astronomical standards.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

[Abridged]We exploit a large homogeneous dataset to derive a self-consistent picture of IR emission based on the time-dependent 24, 15, 12 and 8micron monochromatic and bolometric IR luminosity functions (LF) over the 0<z<2.5 redshift range. Our analysis is based on the combination of data from deep Spitzer surveys in the VVDS-SWIRE and GOODS areas. To our limiting flux of S(24)=400microJy our derived sample in VVDS-SWIRE includes 1494 sources, and 666 and 904 sources brighter than S(24)=80microJy are catalogued in GOODS-S and GOODS-N, respectively, for a total area of ~0.9 square degs. We obtain reliable optical identifications and redshifts, providing us a rich and robust dataset for our luminosity function determination. Based on the multi-wavelength information available, we constrain the LFs at 8, 12, 15 and 24micron. We also extrapolate total IR luminosities from our best-fit to the observed SEDs of each source, and use this to derive the bolometric LF and comoving volume emissivity up to z~2.5. In the 0<z<1 interval, the bolometric IR luminosity density evolves as (1+z)^3.8+/-0.4. Although more uncertain at higher-z, our results show a flattening of the IR luminosity density at z>1. The mean redshift of the peak in the source number density shifts with luminosity: the brighest IR galaxies appear to be forming stars earlier in cosmic time (z>1.5), while the less luminous ones keep doing it at more recent epochs (z~1 for L(IR)<10^11L_sun). Our results suggest a rapid increase of the galaxy IR comoving volume emissivity back to z~1 and a constant average emissivity at z>1. We also seem to find a difference in the evolution rate of the source number densities as a function of luminosity, a downsizing evolutionary pattern similar to that reported from other samples of cosmic sources.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

We present a study of a 20cm selected sample in the Deep SWIRE VLA Field, reaching a limiting flux density of ~13.5 uJy at the image center. In a 0.6x0.6 square degrees field, we are able to assign an optical/IR counterpart to 97% of the radio sources. Up to 11 passbands from the NUV to 4.5um are then used to sample the spectral energy distribution (SED) of these counterparts in order to investigate the nature of the host galaxies. By means of an SED template library and stellar population synthesis models we estimate photometric redshifts, stellar masses, and stellar population properties, dividing the sample in three sub-classes of quiescent, intermediate and star-forming galaxies. We focus on the radio sample in the redshift range 0.3<z<1.3 where we estimate to have a redshift completeness higher than 90%, and study the properties and redshift evolution of these sub-populations. We find that, as expected, the relative contributions of AGN and star-forming galaxies to the uJy population depend on the flux density limit of the sample. At all flux levels a significant population of "green-valley" galaxies is observed. While the actual nature of these sources is not definitely understood, the results of this work may suggest that a significant fraction of faint radio sources might be composite (and possibly transition) objects, thus a simple "AGN vs star-forming" classification might not be appropriate to fully understand what faint radio populations really are.

A complete sample of 33 sources believed to be starbursts ("5.8micron-peakers") was selected in the (0.5 sq. deg.) J1046+56 field with selection criteria F_(24micron)>400muJy, the presence of a redshifted stellar emission peak at 5.8um, and r'(Vega)>23. The field, part of the SWIRE Lockman Hole field, benefits from very deep VLA/GMRT 20cm, 50cm, and 90cm radio data (all 33 sources are detected at 50cm), and deep 160micron and 70micron Spitzer data. The 33 sources, with photometric redshifts ~1.5-2.5, were observed at 1.2mm with IRAM-30m/MAMBO to an rms ~0.7-0.8mJy in most cases. Their millimeter, radio, 7-band Spitzer, and near-IR properties were jointly analyzed. The entire sample of 33 sources has an average 1.2mm flux density of 1.56+/-0.22mJy and a median of 1.61mJy, so the majority of the sources can be considered SMGs. Four sources have confirmed 4sigma detections, and nine were tentatively detected at the 3sigma level. Because of its 24micron selection, our sample shows systematically lower F_(1.2mm)/F_(24micron) flux ratios than classical SMGs, probably because of enhanced PAH emission. A median FIR SED was built by stacking images at the positions of 21 sources in the region of deepest Spitzer coverage. Its parameters are T_(dust)=37+/-8K, L_(FIR)=2.5x10^{12}Lo, and SFR=450Mo/yr. The FIR-radio correlation provides another estimate of L_(FIR) for each source, with an average value of 4.1x10^{12}Lo; however, this value may be overestimated because of some AGN contribution. Most of our targets are also luminous star-forming BzK galaxies which constitute a significant fraction of weak SMGs at 1.7<z<2.3.

We report a moderate-depth (70 ksec), contiguous 0.7 sq.deg, Chandra survey, in the Lockman Hole Field of the Spitzer/SWIRE Legacy Survey coincident with a completed, ultra-deep VLA survey with deep optical and near-infrared imaging in-hand. The primary motivation is to distinguish starburst galaxies and AGN, including the significant, highly obscured (log N_H >23) subset. Chandra has detected 775 X-ray sources to a limiting broad band (0.3-8 keV) flux ~4E-16 erg/cm^2/s. We present the X-ray catalog, fluxes, hardness ratios and multi-wavelength fluxes. The log N vs. log S agrees with those of previous surveys covering similar flux ranges. The Chandra and Spitzer flux limits are well matched: 771 (99%) of the X-ray sources have IR or optical counterparts, and 333 have MIPS 24 um detections. There are 4 optical-only X-ray sources and 4 with no visible optical/IR counterpart. The very deep (~2.7 uJy rms) VLA data yields 251 (> 4 sigma) radio counterparts, 44% of the X-ray sources in the field. We confirm that the tendency for lower X-ray flux sources to be harder is primarily due to absorption. As expected, there is no correlation between observed IR and X-ray flux. Optically bright, Type 1 and red AGN lie in distinct regions of the IR vs X-ray flux plots, demonstrating the wide range of SEDs in this sample and providing the potential for classification/source selection. The classical definition of radio-loudness, R_L, is not effective at selecting strong radio sources for red and relatively optically faint AGN. Using the 24 um to radio flux ratio (q_24) results in the usual ~10% radio-loud AGN (13 of 147).

We present Spitzer 70 and 160 micron observations of the COSMOS Spitzer survey (S-COSMOS). The data processing techniques are discussed for the publicly released products consisting of images and source catalogs. We present accurate 70 and 160 micron source counts of the COSMOS field and find reasonable agreement with measurements in other fields and with model predictions. The previously reported counts for GOODS-North and the extragalactic First Look Survey are updated with the latest calibration, and counts are measured based on the large area SWIRE survey to constrain the bright source counts. We measure an extragalactic confusion noise level of sigma_c = 9.4+/-3.3 mJy (q=5) for the MIPS 160-micron band based on the deep S-COSMOS data and report an updated confusion noise level of sigma_c = 0.35+/-0.15 mJy (q=5) for the MIPS 70-micron band.

We present results on low-resolution mid-infrared (MIR) spectra of 70 infrared-luminous galaxies obtained with the Infrared Spectrograph (IRS) onboard Spitzer. We selected sources from the European Large Area Infrared Survey (ELAIS) with S15 > 0.8 mJy and photometric or spectroscopic z > 1. About half of the sample are QSOs in the optical, while the remaining sources are galaxies, comprising both obscured AGN and starbursts. We classify the spectra using well-known infrared diagnostics, as well as a new one that we propose, into three types of source: those dominated by an unobscured AGN (QSOs), obscured AGN, and starburst-dominated sources. Starbursts concentrate at z ~ 0.6-1.0 favored by the shift of the 7.7-micron PAH band into the selection 15 micron band, while AGN spread over the 0.5 < z < 3.1 range. Star formation rates (SFR) are estimated for individual sources from the luminosity of the PAH features. An estimate of the average PAH luminosity in QSOs and obscured AGN is obtained from the composite spectrum of all sources with reliable redshifts. The estimated mean SFR in the QSOs is 50-100 Mo yr^-1, but the implied FIR luminosity is 3-10 times lower than that obtained from stacking analysis of the FIR photometry, suggesting destruction of the PAH carriers by energetic photons from the AGN. The SFR estimated in obscured AGN is 2-3 times higher than in QSOs of similar MIR luminosity. This discrepancy might not be due to luminosity effects or selection bias alone, but could instead indicate a connection between obscuration and star formation. However, the observed correlation between silicate absorption and the slope of the near- to mid-infrared spectrum is compatible with the obscuration of the AGN emission in these sources being produced in a dust torus.

We report on--off pointed MAMBO observations at 1.2 mm of 61 Spitzer-selected star-forming galaxies from the SWIRE survey. The sources are selected on the basis of bright 24um fluxes (f_24um>0.4mJy) and of stellar dominated near-infrared spectral energy distributions in order to favor z~2 starburst galaxies. The average 1.2mm flux for the whole sample is 1.5+/-0.2 mJy. Our analysis focuses on 29 sources in the Lockman Hole field where the average 1.2mm flux (1.9+/-0.3 mJy) is higher than in other fields (1.1+/-0.2 mJy). The analysis of the sources multi-wavelength spectral energy distributions indicates that they are starburst galaxies with far-infrared luminosities ~10^12-10^13.3 Lsun, and stellar masses of ~0.2-6 x10^11 M_sun. Compared to sub-millimeter selected galaxies (SMGs), the SWIRE-MAMBO sources are among those with the largest 24um/millimeter flux ratios. The origin of such large ratios is investigated by comparing the average mid-infrared spectra and the stacked far-infrared spectral energy distributions of the SWIRE-MAMBO sources and of SMGs. The mid-infrared spectra exhibit strong PAH features, and a warm dust continuum. The warm dust continuum contributes to ~34% of the mid-infrared emission, and is likely associated with an AGN component. This constribution is consistent with what is found in SMGs. The large 24um/1.2mm flux ratios are thus not due to AGN emission, but rather to enhanced PAH emission compared to SMGs. The analysis of the stacked far-infrared fluxes yields warmer dust temperatures than typically observed in SMGs. Our selection favors warm ultra-luminous infrared sources at high-z, a class of objects that is rarely found in SMG samples. Our sample is the largest Spitzer-selected sample detected at millimeter wavelengths currently available.

We derive the properties of dusty tori in Active Galactic Nuclei (AGN) from the comparison of observed Spectral Energy Distributions (SEDs) of SDSS quasars and a precomputed grid of torus models. The observed SEDs comprise SDSS photometry, 2MASS J, H, and K data, whenever available and mid-Infrared (MIR) data from the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. The adopted model is that of Fritz et al., 2006. The fit is performed by standard chi^2 minimisation, the model however can be multi-component comprising a stellar and a starburst components, whenever necessary. Models with low equatorial optical depth, tau_9.7, were allowed as well as ``traditional'' models with tau_9.7 > 1.0, corresponding to A_V > 22 and the results were compared. Fits using high optical depth tori models only produced dust more compactly distributed than in the configuration where all tau_9.7 models were permitted. Tori with decreasing dust density with the distance from the centre were favoured while there was no clear preference for models with or without angular variation of the dust density. The computed outer radii of the tori are of some tens of parsecs large but can reach, in a few cases, a few hundreds of parsecs. The mass of dust, M_Dust, and infrared luminosity, L_IR, integrated in the wavelength range between 1 and 1000 micron, do not show significant variations with redshift, once the observational biases are taken into account. Objects with 70 micron detections, representing 25% of the sample, are studied separately and the starburst contribution (whenever present) to the IR luminosity can reach, in the most extreme but very few cases, 80%.

We present the broad-band Spectral Energy Distributions (SEDs) of the largest available highly (72%) complete spectroscopic sample of mid-infrared (MIR) selected galaxies and AGN at intermediate redshift. The sample contains 203 extragalactic sources from the 15-micron survey in the ELAIS-SWIRE field S1, all with measured spectroscopic redshift. Most of these sources have full multi-wavelength coverage from the far-UV to the far-infrared and lie in the redshift range 0.1<z<1.3. Due to its size, this sample allows us for the first time to characterise the spectral properties of the sources responsible for the strong evolution observed in the MIR. Based on SED-fitting technique we have classified the MIR sources, identifying AGN signatures in about 50% of them. This fraction is significantly higher than that derived from optical spectroscopy (~29%) and is due in particular to the identification of AGN activity in objects spectroscopically classified as galaxies. It is likely that in most of our objects, the AGN is either obscured or of low-luminosity, and thus it does not dominate the energetic output at any wavelength, except in the MIR, showing up just in the range where the host galaxy SED has a minimum. The fraction of AGN strongly depends on the flux density, with that derived through the SED-fitting being about 20% at S(15)~0.5-1 mJy and gradually increasing up to 100% at S(15)>10 mJy, while that obtained from optical spectroscopy never being >30%, even at the higher flux densities. The results of this work will be very useful for updating all the models aimed at interpreting the deep infrared survey data and, in particular, for constraining the nature and the role of dust-obscured systems in the intermediate/high-redshift Universe.

We present the SWIRE Photometric Redshift Catalogue, 1025119 redshifts of unprecedented reliability and accuracy. Our method is based on fixed galaxy and QSO templates applied to data at 0.36-4.5 mu, and on a set of 4 infrared emission templates fitted to infrared excess data at 3.6-170 mu. The code involves two passes through the data, to try to optimize recognition of AGN dust tori. A few carefully justified priors are used and are the key to supression of outliers. Extinction, A_V, is allowed as a free parameter. We use a set of 5982 spectroscopic redshifts, taken from the literature and from our own spectroscopic surveys, to analyze the performance of our method as a function of the number of photometric bands used in the solution and the reduced chi^2. For 7 photometric bands the rms value of (z_{phot}-z_{spec})/(1+z_{spec}) is 3.5%, and the percentage of catastrophic outliers is ~1%. We discuss the redshift distributions at 3.6 and 24 mu. In individual fields, structure in the redshift distribution corresponds to clusters which can be seen in the spectroscopic redshift distribution. 10% of sources in the SWIRE photometric redshift catalogue have z >2, and 4% have z>3, so this catalogue is a huge resource for high redshift galaxies. A key parameter for understanding the evolutionary status of infrared galaxies is L_{ir}/L_{opt}, which can be interpreted as the specific star-formation rate for starbursts. For dust tori around Type 1 AGN, L_{tor}/L_{opt} is a measure of the torus covering factor and we deduce a mean covering factor of 40%.

We present the optical-to-submm spectral energy distributions for 33 radio & mid-IR identified submillimetre galaxies discovered via the SHADES 850 micron SCUBA imaging in the Subaru-XMM Deep Field (SXDF). Optical data for the sources comes from the Subaru-XMM Deep Field (SXDF) and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6 and 4.5 micron fluxes. We then fit spectral energy distribution (SED) templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The infrared template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submm sources and on the SHADES SXDF field. Our template fitting finds that AGN, present in about 10% of sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fit by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submm galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.

This paper presents galaxy source counts at 24 microns in the six Spitzer Wide-field InfraRed Extragalactic (SWIRE) fields. The source counts are compared to counts in other fields, and to model predictions that have been updated since the launch of Spitzer. This analysis confirms a very steep rise in the Euclidean-normalized differential number counts between 2 mJy and 0.3 mJy. Variations in the counts between fields show the effects of sample variance in the flux range 0.5-10 mJy, up to 100% larger than Poisson errors. Nonetheless, a "shoulder" in the normalized counts persists at around 3 mJy. The peak of the normalized counts at 0.3 mJy is higher and narrower than most models predict. In the ELAIS N1 field, the 24 micron data are combined with Spitzer-IRAC data and five-band optical imaging, and these bandmerged data are fit with photometric redshift templates. Above 1 mJy the counts are dominated by galaxies at z less than 0.3. By 300 microJy, about 25% are between z ~ 0.3-0.8, and a significant fraction are at z ~ 1.3-2. At low redshifts the counts are dominated by spirals, and starbursts rise in number density to outnumber the spirals' contribution to the counts below 1 mJy.

We investigate the clustering of galaxies selected in the 3.6 micron band of the Spitzer Wide-area Infrared Extragalactic (SWIRE) legacy survey. The angular two-point correlation function is calculated for eleven samples with flux limits of S_3.6 > 4-400 mujy, over an 8 square degree field. The angular clustering strength is measured at >5-sigma significance at all flux limits, with amplitudes of A=(0.49-29)\times10^{-3} at one degree, for a power-law model, A\theta^{-0.8}. We estimate the redshift distributions of the samples using phenomological models, simulations and photometric redshifts, and so derive the spatial correlation lengths. We compare our results with the GalICS (Galaxies In Cosmological Simulations) models of galaxy evolution and with parameterized models of clustering evolution. The GalICS simulations are consistent with our angular correlation functions, but fail to match the spatial clustering inferred from the phenomological models or the photometric redshifts. We find that the uncertainties in the redshift distributions of our samples dominate the statistical errors in our estimates of the spatial clustering. At low redshifts (median z<0.5) the comoving correlation length is approximately constant, r_0=6.1\pm0.5h^{-1} Mpc, and then decreases with increasing redshift to a value of 2.9\pm0.3h^{-1} Mpc for the faintest sample, for which the median redshift is z=1. We suggest that this trend can be attributed to a decrease in the average galaxy and halo mass in the fainter flux-limited samples, corresponding to changes in the relative numbers of early- and late-type galaxies. However, we cannot rule out strong evolution of the correlation length over 0.5<z<1.

We use a simple optical/infrared (IR) photometric selection of high-redshift QSOs that identifies a Lyman Break in the optical photometry and requires a red IR color to distinguish QSOs from common interlopers. The search yields 100 z~3 (U-dropout) QSO candidates with 19<r'<22 over 11.7 deg^2 in the ELAIS-N1 (EN1) and ELAIS-N2 (EN2) fields of the Spitzer Wide-area Infrared Extragalactic (SWIRE) Legacy Survey. The z~3 selection is reliable, with spectroscopic follow-up of 10 candidates confirming they are all QSOs at 2.83<z<3.44. We find that our z~4$ (g'-dropout) sample suffers from both unreliability and incompleteness but present 7 previously unidentified QSOs at 3.50<z<3.89. Detailed simulations show our z~3 completeness to be ~80-90% from 3.0<z<3.5, significantly better than the ~30-80% completeness of the SDSS at these redshifts. The resulting luminosity function extends two magnitudes fainter than SDSS and has a faint end slope of beta=-1.42 +- 0.15, consistent with values measured at lower redshift. Therefore, we see no evidence for evolution of the faint end slope of the QSO luminosity function. Including the SDSS QSO sample, we have now directly measured the space density of QSOs responsible for ~70% of the QSO UV luminosity density at z~3. We derive a maximum rate of HI photoionization from QSOs at z~3.2, Gamma = 4.8x10^-13 s^-1, about half of the total rate inferred through studies of the Ly-alpha forest. Therefore, star-forming galaxies and QSOs must contribute comparably to the photoionization of HI in the intergalactic medium at z~3.

(Abridged) We selected high-z massive galaxies at 5.8 microns, in the SWIRE ELAIS-S1 field (1 sq. deg.). Galaxies with the 1.6 microns stellar peak redshifted into the IRAC bands (z~1-3, called ``IR-peakers'') were identified. Stellar masses were derived by means of spectro-photometric fitting and used to compute the stellar mass function (MF) at z=1-2 and 2-3. A parametric fit to the MF was performed, based on a Bayesian formalism, and the stellar mass density of massive galaxies above z=2 determined. We present the first systematic study of the very-massive tail of the galaxy stellar mass function at high redshift. A total of 326 sources were selected. The majority of these galaxies have stellar masses in excess of 1e11 Msun and lie at z>1.5. The availability of mid-IR data turned out to be a valuable tool to constrain the contribution of young stars to galaxy SEDs, and thus their M(stars)/L ratio. The influence of near-IR data and of the chosen stellar library on the SED fitting are also discussed. A significant evolution is found not only for galaxies with M~1e11 Msun, but also in the highest mass bins considered. The comoving number density of these galaxies was lower by more than a factor of 10 at z=2-3, with respect to the local estimate. SWIRE 5.8 micron peakers more massive than 1.6x1e11 Msun provide 30-50% of the total stellar mass density in galaxies at z=2-3.

By combining the VVDS with the SWIRE data, we have built the currently largest spectroscopic sample of galaxies selected in the rest-frame near-infrared. These allow us to investigate, for the first time using spectroscopic redshifts, the clustering evolution of galaxies selected from their rest-frame near-infrared luminosity in the redshift range 0.2<z<2.1. Therefore we use the projected two-point correlation function wp(rp) to study the three dimensional clustering properties of galaxies detected at 3.6 and 4.5 microns with IRAC and for which we have spectroscopic redshifts from first epoch VVDS. We find that in the flux limited samples at 3.6 and 4.5 microns, the clustering length does not change from redshift ~2 to the present. In addition, we find that the clustering lengths are systematically higher when galaxy samples are selected from increasingly redder wavelengths. These values are larger than the typical values found for I-band selected galaxies in the same redshift range, but we find that the difference in clustering length between I-band and 3.6/4.5 microns selected samples is decreasing with increasing redshift to become comparable at z~1.5. We interpret this as evidence that galaxies with older stellar populations and galaxies actively forming stars reside in comparably overdense environments at epochs earlier than z~1.5. The increasing difference in clustering length observed between rest-frame UV-optical and infrared selected samples could then be an indication that star formation is gradually shifting to lower density regions as cosmic time increases.

We study the ultraviolet to far-infrared (hereafter UV-to-IR) SEDs of a sample of intermediate redshift (0.2 < z < 0.7) UV-selected galaxies from the ELAIS-N1 and ELAIS-N2 fields by fitting a multi-wavelength dataset to a library of GRASIL templates. Star formation related properties of the galaxies are derived from the library of models by using the Bayesian statistics. We find a decreasing presence of galaxies with low attenuation and low total luminosity as redshift decreases, which does not hold for high total luminosity galaxies. In addition the dust attenuation of low mass galaxies increases as redshift decreases, and this trend seems to disappear for galaxies with M* > 10^11 M_sun. This result is consistent with a mass dependent evolution of the dust to gas ratio, which could be driven by a mass dependent efficiency of star formation in star forming galaxies. The specific star formation rates (SSFR) decrease with increasing stellar mass at all redshifts, and for a given stellar mass the SSFR decreases with decreasing redshift. The differences in the slope of the M*--SSFR relation found between this work and others at similar redshift could be explained by the adopted selection criteria of the samples which, for a UV selected sample, favours blue, star forming galaxies.

We present measurements of the spatial clustering of ultraluminous infrared galaxies in two redshift intervals, 1.5<z<2.0 and 2<z<3. Both samples cluster strongly, with r_0=14.40+/-1.99 h^-1 Mpc for the 2<z<3 sample, and r_0=9.40+/-2.24 h^-1 Mpc for the 1.5<z<2.0 sample, making them among the most biased galaxies at these epochs. These clustering amplitudes are consistent with both populations residing in dark matter haloes with masses of ~6x10^13 solar masses. We infer that a minimum dark matter halo mass is an important factor for all forms of luminous, obscured activity in galaxies at z>1. Adopting plausible models for the growth of DM haloes with redshift, then the haloes hosting the 2<z<3 sample will likely host the richest clusters of galaxies at z=0, whereas the haloes hosting the 1.5<z<2.0 sample will likely host poor to rich clusters at z=0.

We present the stellar calibrator sample and the conversion from instrumental to physical units for the 24 micron channel of the Multiband Imaging Photometer for Spitzer (MIPS). The primary calibrators are A stars, and the calibration factor based on those stars is 4.54*10^{-2} MJy sr^{-1} (DN/s)^{-1}, with a nominal uncertainty of 2%. We discuss the data-reduction procedures required to attain this accuracy; without these procdures, the calibration factor obtained using the automated pipeline at the Spitzer Science Center is 1.6% +/- 0.6% lower. We extend this work to predict 24 micron flux densities for a sample of 238 stars which covers a larger range of flux densities and spectral types. We present a total of 348 measurements of 141 stars at 24 micron. This sample covers a factor of ~460 in 24 micron flux density, from 8.6 mJy up to 4.0 Jy. We show that the calibration is linear over that range with respect to target flux and background level. The calibration is based on observations made using 3-second exposures; a preliminary analysis shows that the calibration factor may be 1% and 2% lower for 10- and 30-second exposures, respectively. We also demonstrate that the calibration is very stable: over the course of the mission, repeated measurements of our routine calibrator, HD 159330, show a root-mean-square scatter of only 0.4%. Finally, we show that the point spread function (PSF) is well measured and allows us to calibrate extended sources accurately; Infrared Astronomy Satellite (IRAS) and MIPS measurements of a sample of nearby galaxies are identical within the uncertainties.

We present the spectral energy distributions (SEDs) of a hard X-ray selected sample. The sample contains 136 sources with F(2-10 keV)>10^-14 erg/cm^2/s and 132 are AGNs. The sources are detected in a 1 square degree area of the XMM-Newton-Medium Deep Survey where optical data from the VVDS, CFHTLS surveys, and infrared data from the SWIRE survey are available. Based on a SED fitting technique we derive photometric redshifts with sigma(1+z)=0.11 and 6% of outliers and identify AGN signatures in 83% of the objects. This fraction is higher than derived when a spectroscopic classification is available. The remaining 17+9-6% of AGNs shows star-forming galaxy SEDs (SF class). The sources with AGN signatures are divided in two classes, AGN1 (33+6-1%) and AGN2 (50+6-11). The AGN1 and AGN2 classes include sources whose SEDs are fitted by type 1 and type 2 AGN templates, respectively. On average, AGN1s show soft X-ray spectra, consistent with being unabsorbed, while AGN2s and SFs show hard X-ray spectra, consistent with being absorbed. The analysis of the average SEDs as a function of X-ray luminosity shows a reddening of the IR SEDs, consistent with a decreasing contribution from the host galaxy at higher luminosities. The AGNs in the SF classes are likely obscured in the mid-infrared, as suggested by their low L(3-20micron)/Lcorr(0.5-10 keV) ratios. We confirm the previously found correlation for AGNs between the radio luminosity and the X-ray and the mid-infrared luminosities. The X-ray-radio correlation can be used to identify heavily absorbed AGNs. However, the estimated radio fluxes for the missing AGN population responsible for the bulk of the background at E>10 keV are too faint to be detected even in the deepest current radio surveys.

Our goal is to probe the populations of obscured and unobscured AGN investigating their optical-IR and X-ray properties as a function of X-ray flux, luminosity and redshift within a hard X-ray selected sample of 136 X-ray sources in the XMM Medium Deep Survey (XMDS) with wide multiwavelength coverage. The XMDS area is covered with optical photometry from the VVDS and CFHTLS surveys and infrared Spitzer data. Based on the X-ray luminosity and X-ray to optical ratio, 132 sources are likely AGN, of which 122 have unambiguous optical - IR identification. The observed optical and IR spectral energy distributions of sources are fitted with AGN/galaxy templates in order to classify them and compute photometric redshifts. 70% of the AGN are fitted by a type 2 AGN or a star forming galaxy template and are grouped together in a single class of ``optically obscured'' AGN. They have ``red'' optical colors and generally show significant X-ray absorption from X-ray spectra or hardness ratios (N$_H > 10^{22}$ cm$^{-2}$). Sources with SEDs typical of type 1 AGN have ``blue'' optical colors and exhibit X-ray absorption in about 30% of cases. We performed a stacking analysis for obscured and type 1 AGN. The stacked X-ray spectrum of obscured AGN is flatter than that of type 1 AGN and has an average spectral slope of Gamma = 1.6. The subsample of objects fitted by a galaxy template has an even harder stacked spectrum, with Gamma = 1.2 - 1.3. The obscured fraction is larger at lower fluxes, lower redshifts and lower luminosities. X-ray absorption is less common than ``optical'' obscuration and its incidence is nearly constant with redshift and luminosity. This implies that X-ray absorption is not necessarily related to optical obscuration.

(Abridged) High-redshift ultra luminous infrared galaxies contribute the bulk of the cosmic IR background and are the best candidates for very massive galaxies in formation at z>1.5. We present Keck/LRIS optical spectroscopy of 35 z>1.4 luminous IR galaxies in the Spitzer Wide-area Infra-Red Extragalactic survey (SWIRE) northern fields (Lockman Hole, ELAIS-N1, ELAIS-N2). The primary targets belong to the ``IR-peak'' class of galaxies, having the 1.6 micron (restframe) stellar feature detected in the IRAC Spitzer channels.The spectral energy distributions of the main targets are thoroughly analyzed, by means of spectro-photometric synthesis and multi-component fits (stars + starburst dust + AGN torus). The IR-peak selection technique is confirmed to successfully select objects above z=1.4, though some of the observed sources lie at lower redshift than expected. Among the 16 galaxies with spectroscopic redshift, 62% host an AGN component, two thirds being type-1 and one third type-2 objects. The selection, limited to r'<24.5, is likely biased to optically-bright AGNs. The SEDs of non-AGN IR-peakers resemble those of starbursts (SFR=20-500 Msun/yr) hosted in massive (M>1e11 Msun) galaxies. The presence of an AGN component provides a plausible explanation for the spectroscopic/photometric redshift discrepancies, as the torus produces an apparent shift of the peak to longer wavelengths. These sources are analyzed in IRAC and optical-IR color spaces. In addition to the IR-peak galaxies, we present redshifts and spectral properties for 150 objects, out of a total of 301 sources on slits.

We study dust attenuation and stellar mass of $\rm z\sim 0.6$ star-forming galaxies using new SWIRE observations in IR and GALEX observations in UV. Two samples are selected from the SWIRE and GALEX source catalogs in the SWIRE/GALEX field ELAIS-N1-00 ($\Omega = 0.8$ deg$^2$). The UV selected sample has 600 galaxies with photometric redshift (hereafter photo-z) $0.5 \leq z \leq 0.7$ and NUV$\leq 23.5$ (corresponding to $\rm L_{FUV} \geq 10^{9.6} L_\sun$). The IR selected sample contains 430 galaxies with $f_{24\mu m} \geq 0.2$ mJy ($\rm L_{dust} \geq 10^{10.8} L_\sun$) in the same photo-z range. It is found that the mean $\rm L_{dust}/L_{FUV}$ ratios of the z=0.6 UV galaxies are consistent with that of their z=0 counterparts of the same $\rm L_{FUV}$. For IR galaxies, the mean $\rm L_{dust}/L_{FUV}$ ratios of the z=0.6 LIRGs ($\rm L_{dust} \sim 10^{11} L_\sun$) are about a factor of 2 lower than local LIRGs, whereas z=0.6 ULIRGs ($\rm L_{dust} \sim 10^{12} L_\sun$) have the same mean $\rm L_{dust}/L_{FUV}$ ratios as their local counterparts. This is consistent with the hypothesis that the dominant component of LIRG population has changed from large, gas rich spirals at z$>0.5$ to major-mergers at z=0. The stellar mass of z=0.6 UV galaxies of $\rm L_{FUV} \leq 10^{10.2} L_\sun$ is about a factor 2 less than their local counterparts of the same luminosity, indicating growth of these galaxies. The mass of z=0.6 UV lunmous galaxies (UVLGs: $\rm L_{FUV} > 10^{10.2} L_\sun$) and IR selected galaxies, which are nearly exclusively LIRGs and ULIRGs, is the same as their local counterparts.

Spectra have been obtained with the Infrared Spectrograph (IRS) on the Spitzer Space Telescope for 20 sources in the Lockman Hole field of the SWIRE survey. The sample is divided between sources with indicators of an obscured AGN, based primarily on X-ray detections of optically-faint sources, and sources with indicators of a starburst, based on optical and near-infrared spectral energy distributions (SEDs) which show a luminosity peak from stellar photospheric emission. Ten of the 11 AGN sources have IRS spectra which show silicate absorption or are power laws; only one AGN source shows PAH emission features. All 9 of the sources showing starburst SEDs in the near-infrared show PAH emission features in the IRS spectra. Redshifts are determined from the IRS spectra for all 9 starbursts (1.0 < z < 1.9) and 8 of the 11 AGN (0.6 < z < 2.5). Classification as AGN because of an X-ray detection, the classification as AGN or starburst derived from the photometric SED, and the IRS spectroscopic classification as AGN (silicate absorption) or starburst (PAH emission) are all consistent in 18 of 20 sources. The surface density for starbursts which are most luminous in the mid-infrared is less than that for the most luminous AGN within the redshift interval 1.7 < z < 1.9. This result implies that mid-infrared source counts at high redshift are dominated by AGN for f(24micron) > 1.0 mJy.

We present the rest-frame optical and infrared colours of a complete sample of 1114 z<0.3 galaxies from the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) and the Sloan Digital Sky Survey (SDSS). We discuss the optical and infrared colours of our sample and analyse in detail the contribution of dusty star-forming galaxies and AGN to optically selected red sequence galaxies. We propose that the optical (g-r) colour and infrared log(L_{24}/L_{3.6}) colour of galaxies in our sample are determined primarily by a bulge-to-disk ratio. The (g-r) colour is found to be sensitive to the bulge-to-disk ratio for disk-dominated galaxies, whereas the log(L_{24}/L_{3.6}) colour is more sensitive for bulge-dominated systems. We identify ~18% (195 sources) of our sample as having red optical colours and infrared excess. Typically, the infrared luminosities of these galaxies are found to be at the high end of star-forming galaxies with blue optical colours. Using emission line diagnostic diagrams, 78 are found to have an AGN contribution, and 117 are identified as star-forming systems. The red (g-r) colour of the star-forming galaxies could be explained by extinction. However, their high optical luminosities cannot. We conclude that they have a significant bulge component. The number densities of optically red star-forming galaxies are found to correspond to ~13% of the total number density of our sample. In addition, these systems contribute ~13% of the total optical luminosity density, and 28% of the total infrared luminosity density of our SWIRE/SDSS sample. These objects may reduce the need for "dry-mergers".

We construct rest-frame luminosity functions at 3.6, 4.5, 5.8, 8 and 24 microns over the redshift range 0<z<2 for galaxies and 0<z<4 for optical QSOs, using optical and infrared data from the Spitzer Wide-area InfraRed Extragalactic survey. The 3.6 and 4.5 micron galaxy LFs show evidence for moderate positive luminosity evolution up to z~1.5, consistent with the passive ageing of evolved stellar populations. Their comoving luminosity density was found to evolve passively, gradually increasing out to z~0.5-1 but flattening, or even declining, at higher redshift. Conversely, the 24 micron galaxy LF, which is more sensitive to obscured star formation and/or AGN activity, undergoes strong positive evolution, with the derived IR energy density and SFR density proportional to (1+z)^gamma with gamma=4.5^{+0.7}_{-0.6} and the majority of this evolution occurring since z~1. Optical QSOs, however, show positive luminosity evolution in all bands, out to the highest redshifts (3<z<4). Modelling as L* proportional to (1+z)^gamma gave gamma=1.3^{+0.1}_{-0.1} at 3.6 micron, gamma=1.0^{+0.1}_{-0.1} at 4.5 micron and stronger evolution at the longer wavelengths (5.8, 8 and 24 micron), of gamma~3. Comparison of the galaxy LFs to predictions from a semi-analytic model based on CDM indicate that an IMF skewed towards higher mass star formation in bursts compared to locally is preferred. As a result the currently inferred massive star formation rates in distant sub-mm sources may require substantial downwards revision.

We fit a parametric model comprising a mixture of multi-dimensional Gaussian functions to the 3.6 to 8um colour and optical photo-z distribution of galaxy populations in the ELAIS-N1 and Lockman Fields of SWIRE. For 16,698 sources in ELAIS-N1 we find our data are best modelled (in the sense of the Bayesian Information Criterion) by the sum of four Gaussian distributions or modes (C_a, C_b, C_c and C_d). We compare the fit of our empirical model with predictions from existing semi-analytic and phenomological models. We infer that our empirical model provides a better description of the mid-infrared colour distribution of the SWIRE survey than these existing models. This colour distribution test is thus a powerful model discriminator and complementary to comparisons of number counts. We use our model to provide a galaxy classification scheme and explore the nature of the galaxies in the different modes of the model. C_a consists of dusty star-forming systems such as ULIRG's. Low redshift late-type spirals are found in C_b, where PAH emission dominates at 8um. C_c consists of dusty starburst systems at intermediate redshifts. Low redshift early-type spirals and ellipticals dominate C_d. We thus find a greater variety of galaxy types than one can with optical photometry alone. Finally we develop a new technique to identify unusual objects, and find a selection of outliers with very red IRAC colours. These objects are not detected in the optical, but have very strong detections in the mid-infrared. These sources are modelled as dust-enshrouded, strongly obscured AGN, where the high mid-infrared emission may either be attributed to dust heated by the AGN or substantial star-formation. These sources have z_ph ~ 2-4, making them incredibly infrared luminous, with a L_IR ~ 10^(12.6-14.1) L_sun.

We present measurements of the spatial clustering of galaxies with stellar masses >10^11Msun, infrared luminosities >10^12 Lsun, and star formation rates >200Msun per year in two redshift intervals; 1.5<z<2.0 and 2<z<3. Both samples cluster very strongly, with spatial correlation lengths of r_0=14.40+/-1.99 h^-1Mpc for the 2<z<3 sample, and r_0=9.40+/-2.24 h^-1Mpc for the 1.5<z<2.0 sample. These clustering amplitudes are consistent with both populations residing in dark matter haloes with masses of ~6x10^13Msun, making them among the most biased galaxies at these epochs. We infer, from this and previous results, that a minimum dark matter halo mass is an important factor for all forms of luminous, obscured activity in galaxies at z>1, both starbursts and AGN. Adopting plausible models for the growth of DM haloes with redshift, then the haloes hosting the 2<z<3 sample will likely host the richest clusters of galaxies at z=0, whereas the haloes hosting the 1.5<z<2.0 sample will likely host poor to rich clusters at z=0. We conclude that ULIRGs at z>1 signpost stellar buildup in galaxies that will reside in clusters at z=0, with ULIRGs at increasing redshifts signposting the buildup of stars in galaxies that will reside in increasingly rich clusters.

We use a simple optical/infrared photometric selection of high redshift QSOs which identifies a Lyman Break in the optical and requires a red IRAC color to distinguish QSOs from common interlopers. We find 100 U-dropout (z=3) QSO candidates with r'<22 within 11.2 deg^2 in the ELAIS-N1 & ELAIS-N2 fields in the Spitzer Wide-Area Infrared Extragalactic (SWIRE) Legacy Survey. Spectroscopy of 10 candidates shows that they are all QSOs with 2.83<z<3.44. We use detailed simulations which incorporate variations in QSO SEDs, IGM transmission and imaging depth to derive a completeness of 85-90% between 3.0<z<3.4. The resulting luminosity function is two magnitudes fainter than SDSS and, when combined with those data, gives a faint end slop \beta = 1.62 \pm 0.18, consistent with measurements at z<2 and steeper than initial measurements at the same redshift.

We discuss spectral energy distributions, photometric redshifts, redshift distributions, luminosity functions, source-counts and the far infrared to optical luminosity ratio for sources in the SWIRE Legacy Survey. The spectral energy distributions of selected SWIRE sources are modelled in terms of a simple set of galaxy and quasar templates in the optical and near infrared, and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst, and AGN dust torus) in the mid infrared. The optical data, together with the IRAC 3.6 and 4.5 mu data, have been used to determine photometric redshifts. For galaxies with known spectroscopic redshifts there is a notable improvement in the photometric redshift when the IRAC data are used, with a reduction in the rms scatter from 10% in (1+z) to 5%. While further spectroscopic data are needed to confirm this result, the prospect of determining good photometric redshifts for the 2 million extragalactic objects in SWIRE is excellent. The distribution of the different infrared sed types in the L{ir}/L{opt} versus L{ir} plane, where L{ir} and L{opt} are the infrared and optical bolometric luminosities, is discussed. Source-counts at 24, 70 and 160 mu are discussed, and luminosity functions at 3.6 and 24 mu are presented.

We present the reduction of the 24 micron data obtained during the first cosmological survey performed by the Spitzer Space Telescope (First Look Survey, FLS). The survey consists of a shallow observation of 2.5x2 sq deg centered at 17h18m +59d30m (main survey) and a deeper observation of 1x0.5 sq deg centered at 17h17m +59d45m(verification survey). Issues with the reduction of the 24 micron MIPS data are discussed and solutions to attenuate instrumental effects are proposed and applied to the data. Approximately 17000 sources are extracted with a SNR greater than five. The photometry of the point sources is evaluated through PSF fitting using an empirical PSF derived from the data. Aperture corrections and the absolute calibration have been checked using stars in the field. Astrometric and photometric errors depend on the SNR of the source varying between 0.35-1 arcsec and 5-15%, respectively, for sources detected at 20-5 sigma. The flux of the 123 extended sources have been estimated through aperture photometry. The extended sources cover less than 0.3% of the total area of the survey. Based on simulations, the main and verification surveys are 50% complete at 0.3 and 0.15 mJy, respectively. Counterparts have been searched for in optical and radio catalogs. More than 80% of the 24 micron sources have a reliable optical counterpart down to R=25.5. 16% of the sources have a 20 cm counterpart down to 0.1 mJy and ~ 80% of the radio-infrared associations have a reliable optical counterpart. A residual map is obtained by subtracting point sources detected at the 3-sigma level and interpolating the regions occupied by extended sources. Several galactic clouds with low and intermediate velocities are identified by comparison with neutral Hydrogen data from this field.

This work focuses on the properties of dusty tori in active galactic nuclei (AGN) derived from the comparison of SDSS type 1 quasars with mid-Infrared (MIR) counterparts and a new, detailed torus model. The infrared data were taken by the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. Basic model parameters are constraint, such as the density law of the graphite and silicate grains, the torus size and its opening angle. A whole variety of optical depths is supported. The favoured models are those with decreasing density with distance from the centre, while there is no clear tendency as to the covering factor, ie small, medium and large covering factors are almost equally distributed. Based on the models that better describe the observed SEDs, properties such as the accretion luminosity, the mass of dust, the inner to outer radius ratio and the hydrogen column density are computed.

Using the large multi-wavelength data set in the chandra/SWIRE Survey (0.6 square degrees in the Lockman Hole), we show evidence for the existence of highly obscured (Compton-thick) AGN, estimate a lower limit to their surface density and characterize their multi-wavelength properties. Two independent selection methods based on the X-ray and infrared spectral properties are presented. The two selected samples contain 1) 5 X-ray sources with hard X-ray spectra and column densities > 10^24 cm-2, and 2) 120 infrared sources with red and AGN-dominated infrared spectral energy distributions (SEDs). We estimate a surface density of at least 25 Compton-thick AGN per square degree detected in the infrared in the chandra/SWIRE field of which ~40% show distinct AGN signatures in their optical/near-infrared SEDs, the remainings being dominated by the host-galaxy emission. Only ~33% of all Compton-thick AGN are detected in the X-rays at our depth (F(0.3-8 keV)>10^-15 erg/cm2/s. We report the discovery of two sources in our sample of Compton-thick AGN, SWIRE_J104409.95+585224.8 (z=2.54) and SWIRE_J104406.30+583954.1 (z=2.43), which are the most luminous Compton-thick AGN at high-z currently known. The properties of these two sources are discussed in detail with an analysis of their spectra, SEDs, luminosities and black-hole masses.

The ESO-Spitzer extragalactic Imaging Survey (ESIS) is the optical follow up of the Spitzer Wide-Area InfraRed Extragalactic (SWIRE) survey in the ELAIS-S1 area. This paper presents B, V, R Wide Field Imager observations of the first 1.5 square degree of the ESIS survey. Data reduction is described including astrometric calibration, illumination and color corrections, completeness and photometric accuracy estimates. Number counts and color distributions are compared to literature observational and theoretical data, including non-evolutionary, PLE, evolutionary and semi-analytic Lambda-CDM galaxy models, as well as Milky Way stellar predictions. ESIS data are in good agreement with previous works and are best reproduced by evolutionary and hierarchical Lambda-CDM scenarios. The ELAIS-S1 area benefits from extensive follow-up from X-ray to radio frequencies: some potential uses of the multi-wavelength observations are illustrated. Optical-Spitzer color-color plots promise to be very powerful tools to disentangle different classes of sources (e.g. AGNs, starbursts, quiescent galaxies). Ultraviolet GALEX data are matched to optical and Spitzer samples, leading to a discussion of galaxy properties in the UV-to-24 microns color space. The spectral energy distribution of a few objects, from the X-rays to the far-IR are presented as examples of the multi-wavelength study of galaxy emission components in different spectral domains.

We address the question of how to deal with confusion limited surveys in the mid-infrared domain by using informations from higher frequency observations over the same sky regions. Such informations, once applied to apparently extended mid-infrared sources, which are indeed ``blends'' of two or more different sources, allow us to disentangle the single counterparts and to split the measured flux density into different components. We present the application of this method to the 24 micron Spitzer archival data in the GOODS-EN1 test field, where apparently extended, ``blended'' sources constitute about 20% of a reliable sample of 983 sources detected above the 5-sigma threshold down to 23 microJy. As higher frequency data-set we have considered the public IRAC images and catalogues on the same field. We show that the 24 micron sample is almost unbiased down to ~40 microJy and the careful application of the deblending procedure does not require any statistical completeness correction (at least at the flux level considered). This is probed by direct comparison of our results with those of Chary et al. (2004), who analysed the same data-set through extensive Monte Carlo simulations. The deblending procedure reduces of about 30% the confusion limit of the MIPS 24 micron survey, allowing one to obtain reliable source counts down to ~40 microJy. The extrapolation of the source counts down to fainter fluxes suggests that our 24 micron sample is able to resolve ~62% of the cosmic background down to a flux level of 38 microJy.

We present ground-based optical and Spitzer infrared imaging observations of the interacting galaxy UGC 10214, the "Tadpole Galaxy" (z = 0.0310), focusing on the star formation activity in the nuclear, disk, spiral arms and tidal tail regions. The major findings of this study are that the Tadpole is actively forming stars in the main disk outside of the nucleus and in the tidal plume, with an estimated mean star formation rate of ~2 to 4 M_sun/yr. The most prominent sites of mid-infrared emission define a "ring" morphology that, combined with the overall morphology of the system, suggest the interaction may belong to the rare class of off-center collisional ring systems that form both shock-induced rings of star formation and tidal plumes. The nuclear emission is solely powered by older stars, with little evidence for ongoing star formation at the center of the Tadpole. Extra-nuclear star formation accounts for >50% of the total star formation in the disk and spiral arms, featuring infrared-bright 'hot spots' that exhibit strong PAH emission, whose band strength is comparable to that of late-type star-forming disk galaxies. The tidal tail, which extends 2 arcmin (~75 kpc) into the intergalactic medium, is populated by super massive star clusters likely triggered by the galaxy-galaxy interaction that has distorted UGC 10214 into its current "tadpole" shape.

We present Spitzer 70um and 160um observations of the Spitzer extragalactic First Look Survey (xFLS). The data reduction techniques and the methods for producing co-added mosaics and source catalogs are discussed. Currently, 26% of the 70um sample and 49% of the 160um-selected sources have redshifts. The majority of sources with redshifts are star-forming galaxies at z<0.5, while about 5% have infrared colors consistent with AGN. The observed infrared colors agree with the spectral energy distribution (SEDs) of local galaxies previously determined from IRAS and ISO data. The average 160um/70um color temperature for the dust is Td~= 30+/-5 K, and the average 70um/24um spectral index is alpha~= 2.4+/-0.4. The observed infrared to radio correlation varies with redshift as expected out to z~1 based on the SEDs of local galaxies. The xFLS number counts at 70um and 160um are consistent within uncertainties with the models of galaxy evolution, but there are indications that the current models may require slight modifications. Deeper 70um observations are needed to constrain the models, and redshifts for the faint sources are required to measure the evolution of the infrared luminosity function.

We present results of a morphological analysis of a small subset of the Spitzer Wide-area InfraRed Extragalactic survey (SWIRE) galaxy population. The analysis is based on public ACS data taken inside the SWIRE N1 field, which are the deepest optical high-resolution imaging available within the SWIRE fields as of today. Our reference sample includes 156 galaxies detected by both ACS and SWIRE. Among the various galaxy morphologies, we disentangle two main classes, spheroids (or bulge-dominated galaxies) and disk-dominated ones, for which we compute the number counts as a function of flux. We then limit our sample to objects with IRAC fluxes brighter than 10 microJy, estimated ~90% completeness limit of the SWIRE catalogues, and compare the observed counts to model predictions. We find that the observed counts of the spheroidal population agree with the expectations of a hierarchical model while a monolithic scenario predicts steeper counts. Both scenaria, however, under-predict the number of late-type galaxies. These observations show that the large majority (close to 80 per cent) of the 3.6 and 4.5 micron galaxy population, even at these moderately faint fluxes, is dominated by spiral and irregular galaxies or mergers.

We present Infrared Array Camera (IRAC) data and source catalogs from the Spitzer Space Telescope Extragalactic First Look Survey. The data were taken in four broad bands centered at nominal wavelengths of 3.6, 4.5, 5.8 and 8.0 microns. A set of mosaics and catalogs have been produced which are ~80% complete and ~99% reliable to their chosen flux density limits. The main field survey covers 3.8 deg^2, and has flux density limits of 20muJy, 25muJy, 100muJy and 100muJy at wavelengths of 3.6, 4.5, 5.8 and 8.0 microns,respectively. The deeper ``verification'' survey covers 0.25 deg^2 with limits of 10muJy, 10muJy, 30muJy and 30muJy, respectively. We also include deep data in the ELAIS-N1 field which covers 0.041deg^2 with limits of 4muJy, 3muJy, 10muJy and 10muJy, respectively, but with only two wavelength coverage at a given sky position. The final bandmerged catalogs contain 103193 objects in the main field, 12224 in the verification field and 5239 in ELAIS-N1. Flux densities of high signal-to-noise objects are accurate to about 10%, and the residual systematic error in the absolute flux density scale is ~2-3%. We have successfully extracted sources at source densities as high as 100000 deg^-2 in our deepest 3.6 and 4.5 micron data. The mosaics and source catalogs will be made available through the Spitzer Science Center archive and the Infrared Science Archive.

We exploit deep combined observations with Spitzer and Chandra of the SWIRE survey in the ELAIS-N1 region, to investigate the nature of the faint X-ray and IR sources in common, to identify AGN/starburst diagnostics, and to study the sources of the X-ray and IR cosmic backgrounds. In the 17'x17' area of the Chandra ACIS-I image there are 3400 SWIRE near-IR sources with 4 sigma detections in at least 2 IRAC bands and 988 sources detected at 24micron with MIPS brighter than 0.1 mJy. Of these, 102 IRAC and 59 MIPS sources have Chandra counterparts, out of a total of 122 X-ray sources present in the area with S(0.5-8 kev)>10^(-15) erg/cm^2/s. We have constructed SEDs for each source using data from the 4 IRAC wavebands, Chandra fluxes, and optical follow-up data in the wavebands U, g', r', i', Z, and H. We fit a number of spectral templates to the SEDs at optical and infrared wavelengths to determine photometric redshifts and spectral categories, and also make use of diagnostics based on the X-ray luminosities, hardness ratios, X-ray to infrared spectral slopes and optical morphologies. Although we have spectroscopic redshifts for only a minority of the Chandra sources, the available SEDs constrain the redshifts for most of the sample sources, which turn out to be typically at 0.5<z<2. We find that 39% of the Chandra sources are dominated by type-1 AGN emission, 23% display optical/IR spectra typical of type-2 AGNs, while the remaining 38% fraction show starburst-like or even normal galaxy spectra. Since we prove that all these galaxies are dominated by AGN emission in X-rays this brings the fraction of type-1 AGNs to be 80% of the type-2: even assuming that all the Chandra sources undetected by Spitzer are type-2 AGNs, the type-1 fraction would exceed 1/3 of the total population (abridged).

We discuss optical associations, spectral energy distributions and photometric redshifts for SWIRE sources in the ELAIS-N1 area and the Lockman Validation Field. The band-merged IRAC (3.6, 4.5, 5.8 and 8.0 mu) and MIPS (24, 70, 160 mu) data have been associated with optical UgriZ data from the INT Wide Field Survey in ELAIS-N1, and with our own optical Ugri data in Lockman-VF. The spectral energy distributions of selected ELAIS sources in N1 detected by SWIRE, most with spectroscopic redshifts, are modelled in terms of a simple set of galaxy and quasar templates in the optical and near infrared, and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst, and AGN dust torus) in the mid infrared. The optical data, together with the IRAC 3.6 and 4.5 mu data, have been used to determine photometric redshifts. For galaxies with known spectroscopic redshifts there is a notable improvement in the photometric redshift when the IRAC data are used, with a reduction in the rms scatter from 10% in (1+z) to 7%. The photometric redshifts are used to derive the 3.6 and 24 mu redshift distribution and to compare this with the predictions of models. For those sources with a clear mid infrared excess, relative to the galaxy starlight model used for the optical and near infrared, the mid and far infrared data are modelled in terms of the same dust emission templates. The proportions found of each template type are: cirrus 31%, M82 29%, Arp 220 10%, AGN dust tori 29%. The distribution of the different infrared sed types in the L_{ir}/L_{opt} versus L_{ir} plane, where L_{ir} and L_{opt} are the infrared and optical bolometric luminosities, is discussed.

We present a mid-infrared analysis of 35 quasars with spectroscopic redshifts selected from the Spitzer Wide-area InfraRed Extragalactic Survey (SWIRE). We discuss their optical and mid-infrared (MIR) colors, and show that these quasars occupy well defined regions in MIR color-color space. We examine the issue of type-I AGN candidate selection in detail and propose new selection methods based on mid-IR colors. The available multi-band data allows us to construct two new, well-sampled quasar templates, covering wavelengths from the ultraviolet to the MIR.

We present the first results of the angular auto-correlation function of the galaxies detected by the Infrared Array Camera (IRAC) instrument in the First Look Survey (FLS) of the Spitzer Space Telescope. We detect significant signals of galaxy clustering within the survey area. The angular auto-correlation function of the galaxies detected in each of the four IRAC instrument channels is consistent with a power-law form $w(\theta)=A\theta^{1-\gamma}$ out to $\theta = 0.2\arcdeg$, with the slope ranging from $\gamma = 1.5$ to 1.8. We estimate the correlation amplitudes $A$ to be $2.95 \times 10^{-3}$, $2.03 \times 10^{-3}$, $4.53 \times 10^{-3}$, and $2.34 \times 10^{-3}$ at $\theta=1\arcdeg$ for galaxies detected in the IRAC 3.6$\mu$m, 4.5$\mu$m, 5.8$\mu$m, and 8.0$\mu$m instrument channels, respectively. We compare our measurements at 3.6$\mu$m with the previous K-band measurements, and discuss the implications of these results.

We report on submillimeter and infrared observations of 28 radio-selected galaxies in the Spitzer First Look Survey Verification field (FLSV). All of the radio-selected galaxies that show evidence for emission at 850um with SCUBA have Spitzer counterparts at 24um, while only half of the radio-selected galaxies without 850um emission have detectable counterparts at 24um. The data show a wide range of infrared colors (S70/S24 < 5--30, S8/S3.6 < 0.3--4), indicative of a mixture of infrared-warm AGN and cooler starburst dominated sources. The galaxies showing 850um emission have Spitzer flux densities and flux density ratios consistent with the range of values expected for high-redshift (z=1--4) ultraluminous infrared galaxies.

We characterize the SWIRE galaxy populations in the SWIRE validation field within the Lockman Hole, based on the 3.6-24$\mu$ Spitzer data and deep U,g',r',r' optical imaging within an area ~1/3 sq. deg for ~16,000 Spitzer-SWIRE sources. The entire SWIRE survey will discover over 2.3 million galaxies at 3.6$\mu$m and almost 350,000 at 24$\mu$m; ~70,000 of these will be 5-band 3.6-24$\mu$ detections. The colors cover a broad range, generally well represented by redshifted spectral energy distributions of known galaxy populations, however significant samples of unusually blue objects in the [3.6-4.5]$\mu$m color are found, as well as many objects very red in the 3.6-24$\mu$m mid-IR. Nine of these are investigated and are interpreted as star-forming systems, starbursts and AGN from z=0.37 to 2.8, with luminosities from L$_{IR}$=10$^{10.3}$ to 10$^{13.7}$ L$_{\odot}$

Using the SPITZER SPACE TELESCOPE and the VLA, we present the first DIRECT evidence that the well-known far-infrared/radio correlation is valid to cosmologically significant redshift. We also confirm, with improved statistics compared with previous surveys, a similar result for the Mid-IR/radio correlation. We explore the dependence of monochromatic q_{24} and q_{70} on z. The results were obtained by matching Spitzer sources at 24 and 70 microns with VLA 1.4 GHz micro-Jy radio sources obtained for the SPITZER FLS. Spectroscopic redshifts have been obtained for over 500 matched IR/radio sources using observations at WIYN, Keck and archival SDSS data extending out to z > 2. We find that q_{24} shows significantly more dispersion than q_{70}. By comparing the observed fluxes at 70, 24 and 4.5 microns with a library of SED templates, we find that the larger dispersion in q_{24} is predictable in terms of systematic variations in SED shape throughout the population. Although the models are not able to encompass the full range of observed behavior (both the presence of either extremely flat or extremely steep IR SEDs), the fitting parameters were used to `k-correct' the higher-z galaxies which resulted in a reduced scatter in q. For comparison, we also corrected these data using the SED for M82. The results for 24 and 70 microns provide strong consistent evidence for the universality of the mid-IR/radio and far-IR/radio correlations out to redshifts of at least z = 1.

We present the Spitzer MIPS 24 micron source counts in the Extragalactic First Look Survey main, verification and ELAIS-N1 fields. Spitzer's increased sensitivity and efficiency in large areal coverage over previous infrared telescopes, coupled with the enhanced sensitivity of the 24 micron band to sources at intermediate redshift, dramatically improve the quality and statistics of number counts in the mid-infrared. The First Look Survey observations cover areas of, respectively, 4.4, 0.26 and 0.015 sq.deg. and reach 3-sigma depths of 0.11, 0.08 and 0.03 mJy. The extragalactic counts derived for each survey agree remarkably well. The counts can be fitted by a super-Euclidean power law of index alpha=-2.9 from 0.2 to 0.9 mJy, with a flattening of the counts at fluxes fainter than 0.2 mJy. Comparison with infrared galaxy evolution models reveals a peak's displacement in the 24 micron counts. This is probably due to the detection of a new population of galaxies with redshift between 1 and 2, previously unseen in the 15 micron deep counts.

In this Letter, we present the initial characterization of extragalactic 24um sources in the Spitzer First Look Survey (FLS) by examining their counterparts at 8um and R-band. The color-color diagram of 24-to-8 vs. 24-to-0.7um is populated with 18,734 sources brighter than the 3sigma flux limit of 110uJy, over an area of 3.7sq.degrees. The 24-to-0.7um colors of these sources span almost 4 orders of magnitudes, while the 24-to-8um colors distribute at least over 2 orders of magnitudes. In addition to identifying ~30% of the total sample with infrared quiescent, mostly low redshift galaxies, we also found that: (1) 23% of the 24um sources (~1200/sq.degrees) have very red 24-to-8 and 24-to-0.7 colors and are probably infrared luminous starbursts with L(IR)>3x10^(11)Lsun at z>1. In particular, 13% of the sample (660/sq.degrees) are 24um detected only, with no detectable emission in either 8um or R-band. These sources are the candidates for being ULIRGs at z>2. (2) 2% of the sample (85/sq.degrees) have colors similar to dust reddened AGNs, like Mrk231 at z~0.6-3. (3) We anticipate that some of these sources with extremely red colors may be new types of sources, since they can not be modelled with any familiar type of spectral energy distribution. We find that 17% of the 24um sources have no detectable optical counterparts brighter than R limit of 25.5mag. Optical spectroscopy of these optical extremely faint 24um sources would be very difficult, and mid-infrared spectroscopy from the Spitzer would be critical for understanding their physical nature (Abridged).

We carried out the direct measurement of the fraction of dusty sources in a sample of extremely red galaxies with (R - Ks) > 5.3mag and Ks < 20.2mag, using 24um data from the Spitzer Space Telescope. Combining deep 24um, Ks- and R-band data over an area of 64sq.arcmin in ELAIS N1, we find that 50%+-6% of our ERO sample have measurable 24um flux above the 3sigma flux limit of 40uJy. This flux limit corresponds to a SFR of 12Msun/yr at z~1. The 24um-detected EROs have 24-to-2.2 and 24-to-0.7um flux ratios consistent with infrared luminous, dusty sources at z>1, and an order of magnitude too red to be explained by an infrared quiescent spiral or a pure old stellar population at any redshift. Keck optical spectroscopy of a sample of similarly selected EROs in the FLS field suggests that most of the EROs in ELAIS N1 are probably at z~1. The mean 24um flux (167uJy) of the 24um-detected ERO sample roughly corresponds to the total infrared luminosity of (3-10)x10^(11)Lsun at z=1-1.5, which is similar to that of local LIRGs and ULIGs. The corresponding SFR is ~(50-170)Msun/yr. If the time scale of this starbursting phase is on the order of 10^8yr as inferred for the local LIRGs and ULIGs, the lower limit on the masses of these 24um-detected EROs is (5-20)x10^9Msun. It is plausible that some of the starburst EROs are in the midst of violent transformation to become massive early type galaxies at the epoch of z ~ 1-2.

Selection of active galactic nuclei (AGN) in the infrared allows the discovery of AGN whose optical emission is extinguished by dust. In this paper, we use the Spitzer Space Telescope First Look Survey (FLS) to assess what fraction of AGN with mid-infrared luminosities comparable to quasars are missed in optical quasar surveys due to dust obscuration. We begin by using the Sloan Digital Sky Survey (SDSS) database to identify 54 quasars within the 4 deg^2 extragalactic FLS. These quasars occupy a distinct region in mid-infrared color space by virtue of their strong, red, continua. This has allowed us to define a mid-infrared color criterion for selecting AGN candidates. About 2000 FLS objects have colors consistent with them being AGN, but most are much fainter in the mid-infrared than the SDSS quasars, which typically have 8 micron flux densities, S(8.0), ~1 mJy. We have investigated the properties of the 43 objects with S(8.0) >= 1 mJy satisfying our AGN color selection. This sample should contain both unobscured quasars, and AGN which are absent from the SDSS survey due to extinction in the optical. After removing 16 known quasars, three probable normal quasars, and eight spurious or confused objects from the initial sample of 43, we are left with 16 objects which are likely to be obscured quasars or luminous Seyfert-2 galaxies. This suggests the numbers of obscured and unobscured AGN are similar in samples selected in the mid-infrared at S(8.0)~1 mJy.

We present the first analysis of large-scale clustering from the Spitzer Wide-area InfraRed Extragalactic legacy survey: SWIRE. We compute the angular correlation function of galaxies selected to have 3.6 micron fluxes brighter than 32 muJy in three fields totaling two square degrees in area. In each field we detect clustering with a high level of significance. The amplitude and slope of the correlation function is consistent between the three fields and is modeled as $w(theta)=A\theta^{1-\gamma}$ with $A=(0.6\pm0.3)\times10^{-3}, \gamma=2.03\pm0.10$. With a fixed slope of $\gamma=1.8$, we obtain an amplitude of $A=(1.7\pm0.1)\times10^{-3}$. Assuming an equivalent depth of K~18.7 mag we find our errors are smaller but our results are consistent with existing clustering measurements in K-band surveys and with stable clustering models. We estimate our median redshift $z\simeq 0.75$ and this allows us to obtain an estimate of the three-dimensional correlation function $\xi(r)$, for which we find $r_0=4.4\pm0.1 h^{-1}$ Mpc.

The SIRTF Wide-area InfraRed Extragalactic survey (SWIRE), the largest SIRTF Legacy program, is a wide-area, imaging survey to trace the evolution of dusty, star-forming galaxies, evolved stellar populations, and AGN as a function of environment, from redshifts z~3 to the current epoch. SWIRE will survey 7 high-latitude fields, totaling 60 - 65 sq. deg. in all 7 SIRTF bands: IRAC 3.6, 4.5, 5.6, 8 microns and MIPS 24, 70, 160 microns. The Legacy Extragalactic Catalog may contain in excess of 2 million IR-selected galaxies, dominated by (1) ~150,000 luminous infrared galaxies (LIRGs: L{FIR}>10^11 L_sun), ~7000 of these with z>2; (2) 1 million early-type galaxies, ~10,000 with z>2; and (3) \~20,000 classical AGN, plus significantly more dust-obscured QSO/AGN among the LIRGs. SWIRE will provide an unprecedented view of the evolution of galaxies, structure, and AGN. The key scientific goals of SWIRE are: (1) to determine the evolution of actively star-forming and passively evolving galaxies in order to understand the history of galaxy formation in the context of cosmic structure formation; (2) to determine the evolution of the spatial distribution and clustering of evolved galaxies, starbursts and AGN in the key redshift range, 0.5<z<3, over which much of cosmic evolution has occurred; (3) to determine the evolutionary relationship between ``normal galaxies'' and AGN, and the contribution of AGN accretion energy vs stellar nucleosynthesis to the cosmic backgrounds. SWIRE's large area is important to create statistically significant population samples over enough volume cells that we can resolve the star formation history as a function of epoch and environment The large volume is also optimised for finding rare objects.

Phenomenological models for evolution of dusty galaxies and E/S0 galaxies, respectively, are developed to address two major questions concenring galaxy populations in deep infrared (IR) surveys: (1) Do normal late-type galaxies or starburst galaxies (including galaxies with obscured AGNs) dominate among sources in deep IR surveys? (2) How much do E/S0 galaxies contribute to the counts in deep mid-infrared (MIR: 3 -- 20$\mu m$) surveys? Among three new models for evolution of dusty galaxies, it is assumed in Model S1 that starburst galaxies are the dominant population, and in Model S2 that normal galaxies dominate. Model S3 is an intermediate model. Comparing the model predictions with a wide range of observational data collected from the literature, we find that none of these models can be ruled out, given the uncertainties of the data. We show that the most direct method to distinguish these models is to compare the predicted color distributions of IR galaxies with observations, which will soon be available from the SWIRE survey. The models for E/S0 galaxies follow a simple passive evolution approach. Among the three E/S0 models (E1, E2 and E3) investigated in this paper, Model E2 which is specified by a peak formation redshift $z_{peak}=2$, and an e-folding formation time scale $\omega=2$ Gyr, fits the data best. This suggests a synchronization between the evolution of E/S0 galaxies and of starburst galaxies, in the sense that the peak of the formation function of E/S0s ($z_{peak}=2$) is close to the peak of the evolution functions of starburst galaxies ($z_{peak} = 1.4$). New predictions for counts and confusion limits in the SIRTF bands are presented.

The soon-to-be-launched Space Infrared Telescope Facility (SIRTF) shall produce image data with an a-posteriori pointing knowledge of 1.4 arcsec (1 sigma radial) with a goal of 1.2 arcsec in the International Celestial Reference System (ICRS). To perform robust image coaddition, mosaic generation, extraction and position determination of faint sources, the pointing will need to be refined to better than a few-tenths of an arcsecond. We use a linear-sparse matrix solver to find a "global-minimization" of all relative image offsets in a mosaic from which refined pointings and orientations can be computed. This paper summarizes the pointing-refinement algorithm and presents the results of testing on simulated data.

Empirical 'backward' galaxy evolution models for IR-bright galaxies are constrained using multiband IR surveys. A new Monte-Carlo algorithm is developed for this task. It exploits a large library of realistic Spectral Energy Distributions (SEDs) of 837 local IR galaxies (IRAS 25$\mu m$ selected) from the UV (1000{\AA}) to the radio (20cm), including ISO-measured 3--13$\mu m$ unidentified broad features (UIBs). The basic assumption is that the local correlation between SEDs and Mid-Infrared (MIR) luminosities can be applied to earlier epochs of the Universe. Three populations of IR sources are considered in the evolution models. These include (1) starburst galaxies, (2) normal late-type galaxies, and (3) galaxies with AGN. A set of models so constructed are compared with data from the literature. Predictions for number counts, confusion limits, redshift distributions, and color-color diagrams are made for multiband surveys using the upcoming SIRTF satellite.

The Very Large Array (VLA) has been used in C-configuration to map an area \~0.3 square degrees at 1.4 GHz with 5sigma sensitivities of 0.305, 0.325, 0.380 and 0.450 mJy/beam over four equal subareas. Radio properties are presented for 62 detected sources. Deep optical imaging to Gunn r~25 mag using the Hale 5-m telescope covering ~0.21 square degrees is reported for a subset of 43 sources. This optical follow-up is much deeper than that of existing larger area radio surveys of similar radio sensitivity. Archival J, H and K-band photometry from the Two-Micron All Sky Survey (2MASS) is also presented. Using a robust likelihood-ratio technique, we optically identified 26 radio sources with probabilities >80%. Nine are uncertain/ambiguous detections and eight are empty fields. Comparisons with a stellar synthesis model that includes radio emission and dust reddening suggests that the near-infrared--to--optical emission in a small, bright sub-sample is reddened by `optically-thin' dust with absorption A_V ~ 2-2.5 mag, regardless of morphological type. This is consistent with other, more direct determinations of absorption. The radio--to--optical(--near-infrared) flux ratios of early-type galaxies require significant contamination in the radio by an active galactic nucleus (AGN), consistent with the current paradigm. Using our simple modeling approach, we also discuss a potential diagnostic for selecting Ultraluminous Infrared Galaxies (ULIGS) to z~1.6 from microJansky radio surveys.

We have used the star camera on the WIRE satellite to observe the K0 III star Alpha UMa, and we report the apparent detection of 10 oscillation modes. The lowest frequency mode is at 1.82 microhertz, and appears to be the fundamental mode. The mean spacing between the mode frequencies is 2.94 microhertz, which implies that all detected modes are radial. The mode frequencies are consistent with the physical parameters of a K0 III star, if we assume that only radial modes are excited. Mode amplitudes are 100 -- 400 micromagnitudes, which is consistent with the scaling relation of Kjeldsen & Beddinge (1995).

In this work we incorporate the newest ISO results on the mid-infrared spectral-energy-distributions (MIR SEDs) of galaxies into models for the number counts and redshift distributions of MIR surveys. A three-component model, with empirically determined MIR SED templates of (1) a cirrus/PDR component (2) a starburst component and (3) an AGN component, is developed for infrared (3--120\micron) SEDs of galaxies. The model includes a complete IRAS 25\micron selected sample of 1406 local galaxies ($z \leq 0.1$; Shupe et al. 1998a). Results based on these 1406 spectra show that the MIR emission features cause significant effects on the redshift dependence of the K-corrections for fluxes in the WIRE 25\micron band and ISOCAM 15\micron band. This in turn will affect deep counts and redshift distributions in these two bands, as shown by the predictions of two evolution models (a luminosity evolution model with $L\propto (1+z)^3$ and a density evolution model with $\rho\propto (1+z)^4$). The dips-and-bumps on curves of MIR number counts, caused by the emission features, should be useful indicators of evolution mode. The strong emission features at $\sim 6$--8\micron will help the detections of relatively high redshift ($z\sim 2$) galaxies in MIR surveys. On the other hand, determinations of the evolutionary rate based on the slope of source counts, and studies on the large scale structures using the redshift distribution of MIR sources, will have to treat the effects of the MIR emission features carefully. We have also estimated a 15\micron local luminosity function from the predicted 15\micron fluxes of the 1406 galaxies using the bivariate (15\micron vs. 25\micron luminosities) method. This luminosity function will improve our understanding of the ISOCAM 15\micron surveys.

We have established a model to systematically estimate the contribution of the mid-infrared emission features between 3 and 11.6 micron to the IRAS in-band fluxes, using the results of ISO PHT-S observation of 16 galaxies by Lu et al. (1997). The model is used to estimate more properly the k-corrections for calculating the restframe 12 and 25 micron fluxes and luminosities of IRAS galaxies. We have studied the 12-25 micron color-luminosity relation for a sample of galaxies selected at 25 microns. The color is found to correlate well with the 25 micron luminosity, the far-infrared and the blue luminosities. The relations with the mid-infrared luminosities are more sensitive to different populations of galaxies, while a single relation of the 12-25 micron color vs. the ratio of the far-infrared and the blue luminosities applies equally well to these different populations. The luminous and ultraluminous infrared galaxies have redder 12-25 micron colors than those of the quasars. These relations provide powerful tools to differentiate different populations of galaxies. We have also selected a sample of 668 galaxies from the IRAS Faint Source Survey flux-density limited at 200 mJy at 12 microns. A 12 micron local luminosity function is derived and, for the first time in the literature, effects of density variation in the local universe are considered and corrected in the calculation of the 12 micron luminosity function. It is also found that the 12 micron-selected sample are dominated by quasars and active galaxies, which therefore strongly affect the 12 micron luminosity function at high luminosities. The ultraluminous infrared galaxies are relatively rare at 12 micron comparing with a 25 micron sample.

The local luminosity function at 25 $\mu$m provides the basis for interpreting the results of deep mid-infrared surveys planned or in progress with space astrophysics missions including ISO, WIRE and SIRTF. We have selected a sample of 1458 galaxies from the IRAS Faint Source Survey with a flux density limit of 250 mJy at 25 $\mu$m. The local luminosity function is derived using both parametric and non-parametric maximum-likelihood techniques, and the classical $1/V_{max}$ estimator. Comparison of these results shows that the $1/V_{max}$ estimate of the luminosity function is significantly affected by the Local Supercluster. A maximum-likelihood fit to the radial density shows no systematic increase that would be caused by density evolution of the galaxy population. The density fit is used to correct the $1/V_{max}$ estimate. We also demonstrate the high quality and completeness of our sample by a variety of methods. The luminosity function derived from this sample is compared to previously published estimates, showing the prior estimates to have been strongly affected by the Local Supercluster. Our new luminosity function leads to lower estimates of mid-infrared backgrounds and number counts.