Massive galaxy clusters undergo strong evolution from z~1.6 to z~0.5, with overdense environments at high-z characterized by abundant dust-obscured star formation and stellar mass growth which rapidly give way to widespread quenching. Data spanning the near- to far-infrared (IR) spectrum can directly trace this transformation; however, such studies have largely been limited to the massive galaxy end of cluster populations. In this work, we present ``total light" stacking techniques spanning 3.4-500{\mu}m aimed at revealing the total cluster IR emission, including low mass members and potential intracluster dust. We detail our procedures for WISE, Spitzer, and Herschel imaging, including corrections to recover the total stacked emission in the case of high fractions of detected galaxies. We apply our stacking techniques to 232 well-studied massive (log M200/Msun~13.8) clusters across multiple z bins, recovering extended cluster emission at all wavelengths, typically at >5sigma. We measure the averaged near- to far-IR radial profiles and SEDs, quantifying the total stellar and dust content. The near-IR radial profiles are well described by an NFW model with a high (c~7) concentration parameter. Dust emission is similarly concentrated, albeit suppressed at small radii (r<0.2Mpc). The measured SEDs lack warm dust, consistent with the colder SEDs expected for low mass galaxies. We derive total stellar masses consistent with the theoretical Mhalo-M_star relation and specific-star formation rates that evolve strongly with redshift, echoing that of massive (log Mstar/Msun>10) cluster galaxies. Separating out the massive galaxy population reveals that the majority of cluster far-IR emission (~70-80%) is provided by the low mass constituents, which differs from field galaxies. This effect may be a combination of mass-dependent quenching and excess dust in low mass cluster galaxies.

We present a multi-wavelength investigation of the radio galaxy population in the galaxy cluster MOO J1506+5137 at $z$=1.09$\pm$0.03, which in previous work we identified as having multiple complex radio sources. The combined dataset used in this work includes data from the Low-Frequency Array Two-metre Sky Survey (LoTSS), NSF's Karl G. Jansky Very Large Array (VLA), the Robert C. Byrd Green Bank Telescope (GBT), the Spitzer Space Telescope, and the Dark Energy Camera Legacy Survey (DECaLS). We find that there are five radio sources which are all located within 500 kpc ($\sim$1$^{\prime}$) of the cluster center and have radio luminosities $P_{\mathrm{1.4GHz}}$ > 1.6$\times$10$^{24}$ W Hz$^{-1}$. The typical host galaxies are among the highest stellar mass galaxies in the cluster. The exceptional radio activity among the massive galaxy population appears to be linked to the dynamical state of the cluster. The galaxy distribution suggests an ongoing merger, with a subgroup found to the northwest of the main cluster. Further, two of the five sources are classified as bent-tail sources with one being a potential wide-angle tail (WAT)/hybrid morphology radio source (HyMoRS) indicating a dynamic environment. The cluster also lies in a region of the mass-richness plane occupied by other merging clusters in the Massive and Distant Clusters of WISE Survey (MaDCoWS). The data suggest that during the merger phase radio activity can be dramatically enhanced, which would contribute to the observed trend of increased radio activity in clusters with increasing redshift.

The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was implemented using the HACC simulation and analysis framework on the BG/Q system, Mira. The cosmological parameters are chosen to be consistent with the results from the Planck satellite. A range of analysis tools have been run in situ to enable a diverse set of science projects, and at the same time, to keep the resulting data amount manageable. Analysis outputs have been generated starting at redshift z~10 to allow for construction of synthetic galaxy catalogs using a semi-analytic modeling approach in post-processing. As part of our in situ analysis pipeline we employ a new method for tracking halo sub-structures, introducing the concept of subhalo cores. The production of multi-wavelength synthetic sky maps is facilitated by generating particle lightcones in situ, also beginning at z~10. We provide an overview of the simulation set-up and the generated data products; a first set of analysis results is presented. A subset of the data is publicly available.

We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA are among the most massive at $z$ = 2-4 (M$_{\star,med}$ = 8.5$ \times$ 10$^{10}$ M$_\odot$) and are either starburst or located in the upper part of the galaxy star-forming main sequence. A significant portion of our galaxy population ($\sim$ 40%), located at $z\sim$ 2.5-3, exhibits abnormally low gas fractions. The sizes of these galaxies, measured with ALMA, are compatible with the trend between $H$-band size and stellar mass observed for $z\sim2$ elliptical galaxies suggesting that they are building compact bulges. We show that there is a strong link between star formation surface density (at 1.1 mm) and gas depletion time: the more compact a galaxy's star-forming region is, the shorter its lifetime will be (without gas replenishment). The identified compact sources associated with relatively short depletion timescales ($\sim$100 Myr), are the ideal candidates to be the progenitors of compact elliptical galaxies at $z$ $\sim$ 2.

In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $\sigma$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $\mu$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$\sigma$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$\sigma$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues.

We construct an emulator for the halo mass function over group and cluster mass scales for a range of cosmologies, including the effects of dynamical dark energy and massive neutrinos. The emulator is based on the recently completed Mira-Titan Universe suite of cosmological $N$-body simulations. The main set of simulations spans 111 cosmological models with 2.1 Gpc boxes. We extract halo catalogs in the redshift range $z=[0.0, 2.0]$ and for masses $M_{200\mathrm{c}}\geq 10^{13}M_\odot/h$. The emulator covers an 8-dimensional hypercube spanned by {$\Omega_\mathrm{m}h^2$, $\Omega_\mathrm{b}h^2$, $\Omega_\nu h^2$, $\sigma_8$, $h$, $n_s$, $w_0$, $w_a$}; spatial flatness is assumed. We obtain smooth halo mass functions by fitting piecewise second-order polynomials to the halo catalogs and employ Gaussian process regression to construct the emulator while keeping track of the statistical noise in the input halo catalogs and uncertainties in the regression process. For redshifts $z\lesssim1$, the typical emulator precision is better than $2\%$ for $10^{13}-10^{14} M_\odot/h$ and $<10\%$ for $M\simeq 10^{15}M_\odot/h$. For comparison, fitting functions using the traditional universal form for the halo mass function can be biased at up to 30\% at $M\simeq 10^{14}M_\odot/h$ for $z=0$. Our emulator is publicly available at \url{https://github.com/SebastianBocquet/MiraTitanHMFemulator}.

Star formation depends critically on cooling mechanisms in the interstellar medium (ISM); however, thermal properties of gas in galaxies at the peak epoch of star formation (z ~ 2) remain poorly understood. A limiting factor in understanding the multiphase ISM is the lack of multiple tracers detected in the same galaxies, such as Polycyclic Aromatic Hydrocarbon (PAH) emission, a tracer of a critical photoelectric heating mechanism in interstellar gas, and [C II] 158{\mu}m fine-structure emission, a principal coolant. We present ALMA Band 9 observations targeting [C II] in six z ~ 2 star-forming galaxies with strong Spitzer IRS detections of PAH emission. All six galaxies are detected in dust continuum and marginally resolved. We compare the properties of PAH and [C II] emission, and constrain their relationship as a function of total infrared luminosity (LIR) and IR surface density. [C II] emission is detected in one galaxy at high signal-to-noise (34{\sigma}), and we place a secure upper limit on a second source. The rest of our sample are not detected in [C II] likely due to redshift uncertainties and narrow ALMA bandpass windows. Our results are consistent with the deficit in [C II]/LIR and PAH/LIR observed in the literature. However, the ratio of [C II] to PAH emission at z ~ 2 is possibly much lower than what is observed in nearby dusty star-forming galaxies. This could be the result of enhanced cooling via [O I] at high-z, hotter gas and dust temperatures, and/or a reduction in the photoelectric efficiency, in which the coupling between interstellar radiation and gas heating is diminished.

We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV slopes (as parametrised by beta), given their high values of the 'infrared excess' (IRX=L_IR/L_UV), relative to the average 'local' IRX-beta relation. We attribute this finding in part to the young ages of the underlying stellar populations but we find that the main reason behind the unusually blue UV slopes are the relatively shallow slopes of the corresponding dust attenuation curves. We show that, when stellar masses are being established via SED fitting, it is absolutely crucial to allow the attenuation curves to vary (rather than fixing it on Calzetti-like law), where we find that the inappropriate curves may underestimate the resulting stellar masses by a factor of ~2-3x on average. In addition, we find these LBGs to have relatively high specific star-formation rates (sSFRs), dominated by the dust component, as quantified via the fraction of obscured star formation ( f_obs = SFR_IR/SFR_(UV+IR)). We conclude that the ALMA-bright LBGs are, by selection, massive galaxies undergoing a burst of a star formation (large sSFRs, driven, for example, by secular or merger processes), with a likely geometrical disconnection of the dust and stars, responsible for producing shallow dust attenuation curves.

The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.

We present the results from a study with NSF's Karl G. Jansky Very Large Array (VLA) to determine the radio morphologies of extended radio sources and the properties of their host galaxies in 50 massive galaxy clusters at z~1. We find a majority of the radio morphologies to be Fanaroff-Riley (FR) type IIs. By analyzing the infrared counterparts of the radio sources, we find that ~40% of the host galaxies are the candidate brightest cluster galaxy (BCG) and ~83% are consistent with being one of the top six most massive galaxies in the cluster. We investigate the role of environmental factors on the radio-loud AGN population by examining correlations between environmental and radio-galaxy properties. We find that the highest stellar mass hosts ($M_{*} \gtrsim$ 4$\times 10^{11} M_{\odot}$) are confined to the cluster center and host compact jets. There is evidence for an increase in the size of the jets with cluster-centric radius, which may be attributed to the decreased ICM pressure confinement with increasing radius. Besides this correlation, there are no other significant correlations between the properties of the radio-AGN (luminosity, morphology, or size) and environmental properties (cluster richness and location within the cluster). The fact that there are more AGN in the cluster environment than the field at this epoch, combined with the lack of strong correlation between galaxy and environmental properties, argues that the cluster environment fosters radio activity but does not solely drive the evolution of these sources at this redshift.

We describe the first major public data release from cosmological simulations carried out with Argonne's HACC code. This initial release covers a range of datasets from large gravity-only simulations. The data products include halo information for multiple redshifts, down-sampled particles, and lightcone outputs. We provide data from two very large LCDM simulations as well as beyond-LCDM simulations spanning eleven w0-wa cosmologies. Our release platform uses Petrel, a research data service, located at the Argonne Leadership Computing Facility. Petrel offers fast data transfer mechanisms and authentication via Globus, enabling simple and efficient access to stored datasets. Easy browsing of the available data products is provided via a web portal that allows the user to navigate simulation products efficiently. The data hub will be extended by adding more types of data products and by enabling computational capabilities to allow direct interactions with simulation results.

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

The bulk of the star-formation rate density peak at cosmic noon was obscured by dust. How accurately we can assess the role of dust obscured star-formation is affected by inherent biases in our empirical methods -- both those that rely on direct dust emission and those that rely on the inferred dust attenuation of starlight. We use a library of hydrodynamic simulations with radiative transfer to explore these biases. We find that for IR luminous galaxies that are in rapidly quenching systems (e.g. post-coalescence) standard luminosity-to-SFR relations can strongly overestimate the true SFRs. We propose using the $L_{IR}/L_{1.6}$ color to both help identify such systems and provide more accurate SFRs. Conversely, we find that the diagnostic UVJ plot misidentifies a subset of dusty star-forming galaxies. This is due to variability in the effective attenuation curves including being much grayer in the optical-to-near-IR regime than the Calzetti starburst law. This is in agreement with recent observations of IR-selected galaxies at cosmic noon. Our results support the view that we need a panchromatic approach from the rest-frame UV through the IR and SED modeling that includes realistic SFHs and allows for variable attenuation curves if we want to fully account for dust obscured star-formation across the epochs of greatest galaxy build-up.

A challenging requirement posed by next-generation observations is a firm theoretical grasp of the impact of baryons on structure formation. Cosmological hydrodynamic simulations modeling gas physics are vital in this regard. A high degree of modeling flexibility exists in this space making it important to explore a range of methods in order to gauge the accuracy of simulation predictions. We present results from the first cosmological simulation using Conservative Reproducing Kernel Smoothed Particle Hydrodynamics (CRK-SPH). We employ two simulations: one evolved purely under gravity and the other with non-radiative hydrodynamics. Each contains 2x2304^3 cold dark matter plus baryon particles in an 800 Mpc/h box. We compare statistics to previous non-radiative simulations including power spectra, mass functions, baryon fractions, and concentration. We find self-similar radial profiles of gas temperature, entropy, and pressure and show that a simple analytic model recovers these results to better than 40% over two orders of magnitude in mass. We quantify the level of non-thermal pressure support in halos and demonstrate that hydrostatic mass estimates are biased low by 24% (10%) for halos of mass 10^15 (10^13) Msun/h. We compute angular power spectra for the thermal and kinematic Sunyaev-Zel'dovich effects and find good agreement with the low-l Planck measurements. Finally, artificial scattering between particles of unequal mass is shown to have a large impact on the gravity-only run and we highlight the importance of better understanding this issue in hydrodynamic applications. This is the first in a simulation campaign using CRK-SPH with future work including subresolution gas treatments.

In order to better understand how active galactic nuclei (AGN) effect the interstellar media of their host galaxies, we perform a meta-analysis of the CO emission for a sample of $z=0.01-4$ galaxies from the literature with existing CO detections and well-constrained AGN contributions to the infrared (67 galaxies). Using either Spitzer/IRS mid-IR spectroscopy or Spitzer+Herschel colors we determine the fraction of the infrared luminosity in each galaxy that can be attributed to heating by the AGN or stars. We calculate new average CO spectral line ratios (primarily from Carilli \& Walter 2013) to uniformly scale the higher-$J$ CO detections to the ground state and accurately determine our sample's molecular gas masses. We do not find significant differences in the gas depletion timescales/star formation efficiencies (SFEs) as a function of the mid-infrared AGN strength ($f_{\rm AGN}({\rm MIR})$ or $L_{\rm IR} ({\rm AGN})$), which indicates that the presence of an IR-bright AGN is not a sufficient sign-post of galaxy quenching. We also find that the dust-to-gas ratio is consistent for all sources, regardless of AGN emission, redshift, or $L_{\rm IR}$, indicating that dust is likely a reliable tracer of gas mass for massive dusty galaxies (albeit with a large degree of scatter). Lastly, if we classify galaxies as either AGN or star formation dominated, we do not find a robust statistically significant difference between their CO excitation.

We describe the Outer Rim cosmological simulation, one of the largest high-resolution N-body simulations performed to date, aimed at promoting science to be carried out with large-scale structure surveys. The simulation covers a volume of (4.225Gpc)^3 and evolves more than one trillion particles. It was executed on Mira, a BlueGene/Q system at the Argonne Leadership Computing Facility. We discuss some of the computational challenges posed by a system like Mira, a many-core supercomputer, and how the simulation code, HACC, has been designed to overcome these challenges. We have carried out a large range of analyses on the simulation data and we report on the results as well as the data products that have been generated. The full data set generated by the simulation totals more than 5PB of data, making data curation and data handling a large challenge in of itself. The simulation results have been used to generate synthetic catalogs for large-scale structure surveys, including DESI and eBOSS, as well as CMB experiments. A detailed catalog for the LSST DESC data challenges has been created as well. We publicly release some of the Outer Rim halo catalogs, downsampled particle information, and lightcone data.

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

Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.

Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

We illustrate the extraordinary potential of the (far-IR) Origins Survey Spectrometer (OSS) on board the Origins Space Telescope (OST) to address a variety of open issues on the co-evolution of galaxies and AGNs. We present predictions for blind surveys, each of 1000 h, with different mapped areas (a shallow survey covering an area of 10 deg$^{2}$ and a deep survey of 1 deg$^{2}$) and two different concepts of the OST/OSS: with a 5.9 m telescope (Concept 2, our reference configuration) and with a 9.1 m telescope (Concept 1, previous configuration). In 1000 h, surveys with the reference concept will detect from $\sim 1.9 \times 10^{6}$ to $\sim 8.7 \times 10^{6}$ lines from $\sim 4.8 \times 10^{5}$-$2.7 \times 10^{6}$ star-forming galaxies and from $\sim 1.4 \times 10^{4}$ to $\sim 3.8 \times 10^{4}$ lines from $\sim 1.3 \times 10^{4}$-$3.5 \times 10^{4}$ AGNs. The shallow survey will detect substantially more sources than the deep one; the advantage of the latter in pushing detections to lower luminosities/higher redshifts turns out to be quite limited. The OST/OSS will reach, in the same observing time, line fluxes more than one order of magnitude fainter than the SPICA/SMI and will cover a much broader redshift range. In particular it will detect tens of thousands of galaxies at $z \geq 5$, beyond the reach of that instrument. The polycyclic aromatic hydrocarbons lines are potentially bright enough to allow the detection of hundreds of thousands of star-forming galaxies up to $z \sim 8.5$, i.e. all the way through the re-ionization epoch. The proposed surveys will allow us to explore the galaxy-AGN co-evolution up to $z\sim 5.5-6$ with very good statistics. OST Concept 1 does not offer significant advantages for the scientific goals presented here.

We present new ALMA dust continuum observations of 101 $\log(\mathrm{M}_* / \mathrm{M}_\odot)$ > 9.5 galaxies in the COSMOS field to study the effect of environment on the interstellar medium at z ~ 0.7. At this redshift, our targets span a wide range of environments allowing for a diverse sample of galaxies with densities, $\Sigma$ = 0.16-10.5 Mpc$^{-2}$ (per $\Delta$ z = 0.024). Using the ALMA observations, we calculate the total ISM mass and look for depletion as a function of galaxy density in order to understand the quenching or triggering of star formation in galaxies in different environments. ISM mass is found to have a small dependence on environment, while the depletion timescale remains constant (~200 Myrs) across all environments. We find elevated ISM mass values at intermediate densities and lower values at high densities compared to low (field) densities. Our observed evolution in gas fraction with density in this single redshift slice is equivalent to the observed evolution with cosmic time over 2-3 Gyr. To explain the change in gas mass fraction seen in galaxies in intermediate and high densities, these results suggest environmental processes such as mergers and ram pressure stripping are likely playing a role in dense filamentary-cluster environments.

The presence of massive neutrinos affects the growth of large-scale structure in the universe, leaving a potentially observable imprint on the abundance and properties of massive dark matter-dominated halos. Cosmological surveys detect large numbers of these halos in the form of rich groups and clusters, using the information as an input to constraining the properties of dark energy. We use a suite of N-body simulations that include the effects of massive neutrinos as well as of dynamical dark energy to study the properties of the mass function. As in our previous work, we follow an approach valid at low neutrino mass, where the neutrino overdensities are assumed to be too small to act as a significant nonlinear source term for gravity. We study how well a universal form for the halo mass function describes our numerical results, finding that the use of an appropriate linear power spectrum within the formalism yields a good match to the simulation results, correctly accounting for the (neutrino mass-dependent) suppression of the mass function.

We present deep $J$ and $H$-band images in the extended Great Observatories Origins Deep Survey-North (GOODS-N) field covering an area of 0.22 $\rm{deg}^{2}$. The observations were taken using WIRCam on the 3.6-m Canada France Hawaii Telescope (CFHT). Together with the reprocessed $K_{\rm s}$-band image, the $5\sigma$ limiting AB magnitudes (in 2" diameter apertures) are 24.7, 24.2, and 24.4 AB mag in the $J$, $H$, and $K_{\rm s}$ bands, respectively. We also release a multi-band photometry and photometric redshift catalog containing 93598 sources. For non-X-ray sources, we obtained a photometric redshift accuracy $\sigma_{\mathrm{NMAD}}=0.036$ with an outlier fraction $\eta = 7.3\%$. For X-ray sources, which are mainly active galactic nuclei (AGNs), we cross-matched our catalog with the updated 2M-CDFN X-ray catalog from Xue et al. (2016) and found that 658 out of 683 X-ray sources have counterparts. $GALEX$ UV data are included in the photometric redshift computation for the X-ray sources to give $\sigma_{\mathrm{NMAD}} = 0.040$ with $\eta=10.5\%$. Our approach yields more accurate photometric redshift estimates compared to previous works in this field. In particular, by adopting AGN-galaxy hybrid templates, our approach delivers photometric redshifts for the X-ray counterparts with fewer outliers compared to the 3D-HST catalog, which fit these sources with galaxy-only templates.

We present the results from a pilot study with the Karl G. Jansky Very Large Array (JVLA) to determine the radio morphologies of extended radio sources and the properties of their host-galaxies in 10 massive galaxy clusters at z~1, an epoch in which clusters are assembling rapidly. These clusters are drawn from a parent sample of WISE-selected galaxy clusters that were cross-correlated with the VLA Faint Images of the Radio Sky at Twenty-Centimeters survey (FIRST) to identify extended radio sources within 1$^{\prime}$ of the cluster centers. Out of the ten targeted sources, six are FR II sources, one is an FR I source, and three sources have undetermined morphologies. Eight radio sources have associated Spitzer data, 75% presenting infrared counterparts. A majority of these counterparts are consistent with being massive galaxies. The angular extent of the FR sources exhibits a strong correlation with the cluster-centric radius, which warrants further investigation with a larger sample.

The Origins Space Telescope, one of four large Mission Concept studies sponsored by NASA for review in the 2020 US Astrophysics Decadal Survey, will open unprecedented discovery space in the infrared, unveiling our cosmic origins. We briefly describe in this article the key science themes and architecture for OST. With a sensitivity gain of up to a factor of 1,000 over any previous or planned mission, OST will open unprecedented discovery space, allow us to peer through an infrared window teeming with possibility. OST will fundamentally change our understanding of our cosmic origins - from the growth of galaxies and black holes, to uncovering the trail of water, to life signs in nearby Earth-size planets, and discoveries never imagined. Built to be highly adaptable, while addressing key science across many areas of astrophysics, OST will usher in a new era of infrared astronomy.

We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes

Profiles of dark matter-dominated halos at the group and cluster scales play an important role in modern cosmology. Using results from two very large cosmological $N$-body simulations, which increase the available volume at their mass resolution by roughly two orders of magnitude, we robustly determine the halo concentration-mass $(c-M)$ relation over a wide range of masses, employing multiple methods of concentration measurement. We characterize individual halo profiles, as well as stacked profiles, relevant for galaxy-galaxy lensing and next-generation cluster surveys; the redshift range covered is $0\leq z \leq 4$, with a minimum halo mass of $M_{200c}\sim2\times10^{11} M_\odot$. Despite the complexity of a proper description of a halo (environmental effects, merger history, nonsphericity, relaxation state), when the mass is scaled by the nonlinear mass scale $M_\star(z)$, we find that a simple non-power-law form for the $c-M/M_\star$ relation provides an excellent description of our simulation results across eight decades in $M/M_{\star}$ and for $0\leq z \leq 4$. Over the mass range covered, the $c-M$ relation has two asymptotic forms: an approximate power law below a mass threshold $M/M_\star\sim 500-1000$, transitioning to a constant value, $c_0\sim 3$ at higher masses. The relaxed halo fraction decreases with mass, transitioning to a constant value of $\sim 0.5$ above the same mass threshold. We compare Navarro-Frenk-White (NFW) and Einasto fits to stacked profiles in narrow mass bins at different redshifts; as expected, the Einasto profile provides a better description of the simulation results. At cluster scales at low redshift, however, both NFW and Einasto profiles are in very good agreement with the simulation results, consistent with recent weak lensing observations.

We present the results of a new study of the relationship between infrared excess (IRX), UV spectral slope (beta) and stellar mass at redshifts 2<z<3, based on a deep Atacama Large Millimeter Array (ALMA) 1.3-mm continuum mosaic of the Hubble Ultra Deep Field (HUDF). Excluding the most heavily-obscured sources, we use a stacking analysis to show that z~2.5 star-forming galaxies in the mass range 9.25 <= log(M/Msun) <= 10.75 are fully consistent with the IRX-beta relation expected for a relatively grey attenuation curve, similar to the commonly adopted Calzetti law. Based on a large, mass complete, sample of 2 <= z <= 3 star-forming galaxies drawn from multiple surveys, we proceed to derive a new empirical relationship between beta and stellar mass, making it possible to predict UV attenuation (A_1600) and IRX as a function of stellar mass, for any assumed attenuation law. Once again, we find that z~2.5 star-forming galaxies follow A_1600-mass and IRX-mass relations consistent with a relatively grey attenuation law, and find no compelling evidence that star-forming galaxies at this epoch follow a reddening law as steep as the Small Magellanic Cloud (SMC) extinction curve. In fact, we use a simple simulation to demonstrate that previous determinations of the IRX-beta relation may have been biased toward low values of IRX at red values of beta, mimicking the signature expected for an SMC-like dust law. We show that this provides a plausible mechanism for reconciling apparently contradictory results in the literature and that, based on typical measurement uncertainties, stellar mass provides a cleaner prediction of UV attenuation than beta. Although the situation at lower stellar masses remains uncertain, we conclude that for 2<z<3 star-forming galaxies with log(M/Msun) >= 9.75, both the IRX-beta and IRX-mass relations are well described by a Calzetti-like attenuation law.

The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs.

The total star formation budget of galaxies consists of the sum of the unobscured star formation, as observed in the rest-frame ultraviolet (UV), together with the obscured component that is absorbed and re-radiated by dust grains in the infrared. We explore how the fraction of obscured star formation depends on stellar mass for mass-complete samples of galaxies at $0<z<2.5$. We combine GALEX and WISE photometry for SDSS-selected galaxies with the 3D-HST treasury program and Spitzer/MIPS 24$\mu$m photometry in the well-studied 5 extragalactic CANDELS fields. We find a strong dependence of the fraction of obscured star formation ($f_{\mathrm{obscured}}$=SFR$_{\mathrm{IR}}$/SFR$_{\mathrm{UV+IR}}$) on stellar mass, with remarkably little evolution in this fraction with redshift out to $z$=2.5. 50\% of star formation is obscured for galaxies with log(M/M$_{\odot}$)=9.4; although unobscured star formation dominates the budget at lower masses, there exists a tail of low mass extremely obscured star-forming galaxies at $z>1$. For log(M/M$_{\odot}$)$>$10.5, $>$90\% of star formation is obscured at all redshifts. We also show that at fixed total SFR, $f_{\mathrm{obscured}}$ is lower at higher redshift. At fixed mass, high-redshift galaxies are observed to have more compact sizes and much higher star formation rates, gas fractions and hence surface densities (implying higher dust obscuration), yet we observe no redshift evolution in $f_{\mathrm{obscured}}$ with stellar mass. This poses a challenge to theoretical models to reproduce, where the observed compact sizes at high redshift seem in tension with lower dust obscuration.

Several recent studies have reported different intrinsic correlations between the AGN mid-IR luminosity ($L_{MIR}$) and the rest-frame 2-10 keV luminosity ($L_{X}$) for luminous quasars. To understand the origin of the difference in the observed $L_{X}-L_{MIR}$ relations, we study a sample of 3,247 spectroscopically confirmed type 1 AGNs collected from Bo\"{o}tes, XMM-COSMOS, XMM-XXL-North, and the SDSS quasars in the Swift/XRT footprint spanning over four orders of magnitude in luminosity. We carefully examine how different observational constraints impact the observed $L_{X}-L_{MIR}$ relations, including the inclusion of X-ray non-detected objects, possible X-ray absorption in type 1 AGNs, X-ray flux limits, and star formation contamination. We find that the primary factor driving the different $L_{X}-L_{MIR}$ relations reported in the literature is the X-ray flux limits for different studies. When taking these effects into account, we find that the X-ray luminosity and mid-IR luminosity (measured at rest-frame $6\mu m$, or $L_{6\mu m}$) of our sample of type 1 AGNs follow a bilinear relation in the log-log plane: $\log L_X =(0.84\pm0.03)\times\log L_{6\mu m}/10^{45}{\rm erg\;s^{-1}} + (44.60\pm0.01)$ for $L_{6\mu m} < 10^{44.79}{\rm erg\;s^{-1}} $, and $\log L_X = (0.40\pm0.03)\times\log L_{6\mu m}/10^{45}{\rm erg\;s^{-1}} +(44.51\pm0.01)$ for $L_{6\mu m} \geq 10^{44.79}{\rm erg\;s^{-1}} $. This suggests that the luminous type 1 quasars have a shallower $L_{X}-L_{MIR}$ correlation than the approximately linear relations found in local Seyfert galaxies. This result is consistent with previous studies reporting a luminosity-dependent $L_{X}-L_{MIR}$ relation, and implies that assuming a linear $L_{X}-L_{MIR}$ relation to infer the neutral gas column density for X-ray absorption might overestimate the column densities in luminous quasars.

The bulk of the stellar growth over cosmic time is dominated by IR luminous galaxies at cosmic noon (z=1-2), many of which harbor a hidden active galactic nucleus (AGN). We use state of the art infrared color diagnostics, combining Spitzer and Herschel observations, to separate dust-obscured AGN from dusty star forming galaxies (SFGs) in the CANDELS and COSMOS surveys. We calculate 24 micron counts of SFGs, AGN/star forming "Composites", and AGN. AGN and Composites dominate the counts above 0.8 mJy at 24 micron, and Composites form at least 25% of an IR sample even to faint detection limits. We develop methods to use the Mid-Infrared Instrument (MIRI) on JWST to identify dust-obscured AGN and Composite galaxies from z~1-2. With the sensitivity and spacing of MIRI filters, we will detect >4 times as many AGN hosts than with Spitzer/IRAC criteria. Any star formation rates based on the 7.7 micron PAH feature (likely to be applied to MIRI photometry) must be corrected for the contribution of the AGN, or the SFR will be overestimated by ~35% for cases where the AGN provides half the IR luminosity and ~50% when the AGN accounts for 90% of the luminosity. Finally, we demonstrate that our MIRI color technique can select AGN with an Eddington ratio of $\lambda_{\rm Edd}\sim0.01$ and will identify AGN hosts with a higher sSFR than X-ray techniques alone. JWST/MIRI will enable critical steps forward in identifying and understanding dust-obscured AGN and the link to their host galaxies.

At $z=1-3$, the formation of new stars is dominated by dusty galaxies whose far-IR emission indicates they contain colder dust than local galaxies of a similar luminosity. We explore the reasons for the evolving IR emission of similar galaxies over cosmic time using: 1) Local galaxies from GOALS $(L_{\rm IR}=10^{11}-10^{12}\,L_\odot)$; 2) Galaxies at $z\sim0.1-0.5$ from the 5MUSES ($L_{\rm IR}=10^{10}-10^{12}\,L_\odot$); 3) IR luminous galaxies spanning $z=0.5-3$ from GOODS and Spitzer xFLS ($L_{\rm IR}>10^{11}\,L_\odot$). All samples have Spitzer mid-IR spectra, and Herschel and ground-based submillimeter imaging covering the full IR spectral energy distribution, allowing us to robustly measure $L_{\rm IR}^{\rm\scriptscriptstyle SF}$, $T_{\rm dust}$, and $M_{\rm dust}$ for every galaxy. Despite similar infrared luminosities, $z>0.5$ dusty star forming galaxies have a factor of 5 higher dust masses and 5K colder temperatures. The increase in dust mass is linked with an increase in the gas fractions with redshift, and we do not observe a similar increase in stellar mass or star formation efficiency. $L_{160}^{\rm\scriptscriptstyle SF}/L_{70}^{\rm\scriptscriptstyle SF}$, a proxy for $T_{\rm dust}$, is strongly correlated with $L_{\rm IR}^{\rm\scriptscriptstyle SF}/M_{\rm dust}$ independently of redshift. We measure merger classification and galaxy size for a subsample, and there is no obvious correlation between these parameters and $L_{\rm IR}^{\rm \scriptscriptstyle SF}/M_{\rm dust}$ or $L_{160}^{\rm\scriptscriptstyle SF}/L_{70}^{\rm\scriptscriptstyle SF}$. In dusty star forming galaxies, the change in $L_{\rm IR}^{\rm\scriptscriptstyle SF}/M_{\rm dust}$ can fully account for the observed colder dust temperatures, suggesting that any change in the spatial extent of the interstellar medium is a second order effect.

We introduce a new cosmic emulator for the matter power spectrum covering eight cosmological parameters. Targeted at optical surveys, the emulator provides accurate predictions out to a wavenumber k~5/Mpc and redshift z<=2. Besides covering the standard set of LCDM parameters, massive neutrinos and a dynamical dark energy of state are included. The emulator is built on a sample set of 36 cosmological models, carefully chosen to provide accurate predictions over the wide and large parameter space. For each model, we have performed a high-resolution simulation, augmented with sixteen medium-resolution simulations and TimeRG perturbation theory results to provide accurate coverage of a wide k-range; the dataset generated as part of this project is more than 1.2Pbyte. With the current set of simulated models, we achieve an accuracy of approximately 4%. Because the sampling approach used here has established convergence and error-control properties, follow-on results with more than a hundred cosmological models will soon achieve ~1% accuracy. We compare our approach with other prediction schemes that are based on halo model ideas and remapping approaches. The new emulator code is publicly available.

ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample has known far-infrared luminosities and, hence, star formation rates (SFRs), and stellar masses (M$_{\rm *}$) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, M$_{\rm *}$, and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFR$_MS$. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolutionary and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function we estimate the accretion rates of gas required to maintain continuity of the MS evolution ($>100$ \msun yr$^{-1}$ at z $>$ 2.5). Simple power-law dependences are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and AGN activity at early epochs.

We directly detect dust emission in an optically-detected, multiply-imaged galaxy lensed by the Frontier Fields cluster MACSJ0717.5+3745. We detect two images of the same galaxy at 1.1mm with the AzTEC camera on the Large Millimeter Telescope leaving no ambiguity in the counterpart identification. This galaxy, MACS071_Az9, is at z>4 and the strong lensing model (mu=7.5) allows us to calculate an intrinsic IR luminosity of 9.7e10 Lsun and an obscured star formation rate of 14.6 +/- 4.5 Msun/yr. The unobscured star formation rate from the UV is only 4.1 +/- 0.3 Msun/yr which means the total star formation rate (18.7 +/- 4.5 Msun/yr) is dominated (75-80%) by the obscured component. With an intrinsic stellar mass of only 6.9e9Msun, MACS0717_Az9 is one of only a handful of z>4 galaxies at these lower masses that is detected in dust emission. This galaxy lies close to the estimated star formation sequence at this epoch. However, it does not lie on the dust obscuration relation (IRX-beta) for local starburst galaxies and is instead consistent with the Small Magellanic Cloud (SMC) attenuation law. This remarkable lower mass galaxy showing signs of both low metallicity and high dust content may challenge our picture of dust production in the early Universe.

The James Webb Space Telescope's Medium Resolution Spectrometer (MRS), will offer nearly 2 orders of magnitude improvement in sensitivity and >3X improvement in spectral resolution over our previous space-based mid-IR spectrometer, the Spitzer IRS. In this paper, we make predictions for spectroscopic pointed observations and serendipitous detections with the MRS. Specifically, pointed observations of Herschel sources require only a few minutes on source integration for detections of several star-forming and active galactic nucleus lines, out to z$=$3 and beyond. But the same data will also include tens of serendipitous 0$\lesssim$z$\lesssim$4 galaxies per field with infrared luminosities ranging $\sim10^6-10^{13}$L$_{\odot}$. In particular, for the first time and for free we will be able to explore the $L_{IR}<10^{9}L_{\odot}$ regime out to $z\sim3$. We estimate that with $\sim$100 such fields, statistics of these detections will be sufficient to constrain the evolution of the low-$L$ end of the infrared luminosity function, and hence the star formation rate function. The above conclusions hold for a wide range in potential low-$L$ end of the IR luminosity function, and accounting for the PAH deficit in low-$L$, low-metallicity galaxies.

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts (<z> = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4.

We present $\simeq$0$.\!\!^{\prime\prime}4$-resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at $z = 1.3-3.0$. These galaxies are selected from sensitive, blank-field surveys of the $2' \times 2'$ Hubble Ultra-Deep Field at $\lambda = 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA). They have star-formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at $z \sim 2$. Morphological classification performed on spatially-resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGN), thereby representing a diversity of $z \sim 2$ SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of $4.2 \pm 1.8$ kpc. This provides direct evidence for galaxy-wide star formation in distant, blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGN, which are found to be co-spatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies whose SFRs are $3-8$ times larger, providing a constraint on the characteristic SFR ($\sim300$ M$_{\odot}$yr$^{-1}$) above which a significant population of more compact star-forming galaxies appears to emerge.

A key question in extragalactic studies is the determination of the relative roles of stars and AGN in powering dusty galaxies at $z\sim$1-3 where the bulk of star-formation and AGN activity took place. In Paper I, we present a sample of $336$ 24$\mu$m-selected (Ultra)Luminous Infrared Galaxies, (U)LIRGs, at $z \sim 0.3$-$2.8$, where we focus on determining the AGN contribution to the IR luminosity. Here, we use hydrodynamic simulations with dust radiative transfer of isolated and merging galaxies, to investigate how well the simulations reproduce our empirical IR AGN fraction estimates and determine how IR AGN fractions relate to the UV-mm AGN fraction. We find that: 1) IR AGN fraction estimates based on simulations are in qualitative agreement with the empirical values when host reprocessing of the AGN light is considered; 2) for star-forming galaxy-AGN composites our empirical methods may be underestimating the role of AGN, as our simulations imply $>$50% AGN fractions, $\sim$3$\times$ higher than previous estimates; 3) 6% of our empirically classified "SFG" have AGN fractions $\gtrsim$50%. While this is a small percentage of SFGs, if confirmed, would imply the true number density of AGN may be underestimated; 4) this comparison depends on the adopted AGN template -- those that neglect the contribution of warm dust lower the empirical fractions by up to 2$\times$; and 5) the IR AGN fraction is only a good proxy for the intrinsic UV-mm AGN fraction when the extinction is high ($A_V\gtrsim 1$ or up to and including coalescence in a merger).

This paper describes the beginning of the Far-Infrared Surveyor mission study for NASA's Astrophysics Decadal 2020. We describe the scope of the study, and the open process approach of the Science and Technology Definition Team. We are currently developing the science cases and provide some preliminary highlights here. We note key areas for technological innovation and improvements necessary to make a Far-Infrared Surveyor mission a reality.

ALMA Cycle 2 observations of the long wavelength dust emission in 145 star-forming galaxies are used to probe the evolution of star-forming ISM. We also develop the physical basis and empirical calibration (with 72 low-z and z ~ 2 galaxies) for using the dust continuum as a quantitative probe of interstellar medium (ISM) masses. The galaxies with highest star formation rates (SFRs) at <z> = 2.2 and 4.4 have gas masses up to 100 times that of the Milky Way and gas mass fractions reaching 50 to 80%, i.e. gas masses 1 - 4 times their stellar masses. We find a single high-z star formation law: SFR = 35 M_ mol^0.89 x (1+z)_{z=2}^0.95 x (sSFR)_{MS}^0.23 \msun yr^-1 -- an approximately linear dependence on the ISM mass and an increased star formation efficiency per unit gas mass at higher redshift. Galaxies above the Main Sequence (MS) have larger gas masses but are converting their ISM into stars on a timescale only slightly shorter than those on the MS -- thus these 'starbursts' are largely the result of having greatly increased gas masses rather than and increased efficiency for converting gas to stars. At z $> 1$, the entire population of star-forming galaxies has $\sim$ 2 - 5 times shorter gas depletion times than low-z galaxies. These shorter depletion times indicate a different mode of star formation in the early universe -- most likely dynamically driven by compressive, high-dispersion gas motions -- a natural consequence of the high gas accretion rates.

We present a detailed, multi-wavelength study of star formation (SF) and AGN activity in 11 near-infrared (IR) selected, spectroscopically confirmed, massive ($\gtrsim10^{14}\,\rm{M_{\odot}}$) galaxy clusters at $1<z<1.75$. Using new, deep $Herschel$/PACS imaging, we characterize the optical to far-IR spectral energy distributions (SEDs) for IR-luminous cluster galaxies, finding that they can, on average, be well described by field galaxy templates. Identification and decomposition of AGN through SED fittings allows us to include the contribution to cluster SF from AGN host galaxies. We quantify the star-forming fraction, dust-obscured SF rates (SFRs), and specific-SFRs for cluster galaxies as a function of cluster-centric radius and redshift. In good agreement with previous studies, we find that SF in cluster galaxies at $z\gtrsim1.4$ is largely consistent with field galaxies at similar epochs, indicating an era before significant quenching in the cluster cores ($r<0.5\,$Mpc). This is followed by a transition to lower SF activity as environmental quenching dominates by $z\sim1$. Enhanced SFRs are found in lower mass ($10.1< \log \rm{M_{\star}}/\rm{M_{\odot}}<10.8$) cluster galaxies. We find significant variation in SF from cluster-to-cluster within our uniformly selected sample, indicating that caution should be taken when evaluating individual clusters. We examine AGN in clusters from $z=0.5-2$, finding an excess AGN fraction at $z\gtrsim1$, suggesting environmental triggering of AGN during this epoch. We argue that our results $-$ a transition from field-like to quenched SF, enhanced SF in lower mass galaxies in the cluster cores, and excess AGN $-$ are consistent with a co-evolution between SF and AGN in clusters and an increased merger rate in massive haloes at high redshift.

ALMA Cycle 2 observations of the long wavelength dust emission in 180 star-forming (SF) galaxies are used to investigate the evolution of ISM masses at z = 1 to 6.4. The ISM masses exhibit strong increases from z = 0 to $\rm <z>$ = 1.15 and further to $\rm <z>$ = 2.2 and 4.8, particularly amongst galaxies above the SF galaxy main sequence (MS). The galaxies with highest SFRs at $\rm <z>$ = 2.2 and 4.8 have gas masses 100 times that of the Milky Way and gas mass fractions reaching 50 to 80\%, i.e. gas masses 1 - 4$\times$ their stellar masses. For the full sample of galaxies, we find a single, very simple SF law: $\rm SFR \propto M_{\rm ISM}^{0.9}$, i.e. a `linear' dependence on the ISM mass -- on and above the MS. Thus, the galaxies above the MS are converting their larger ISM masses into stars on a timescale similar to those on the MS. At z $> 1$, the entire population of star-forming galaxies has $\sim$5 - 10$\times$ shorter gas depletion times ($\sim0.2$ Gyr) than galaxies at low redshift. These {\bf shorter depletion times are due to a different, dominant mode of SF in the early universe} -- dynamically driven by compressive, high dispersion gas motions and/or galaxy interactions. The dispersive gas motions are a natural consequence of the extraordinarily high gas accretion rates which must occur to maintain the prodigious SF.

Large-scale structure surveys in the coming years will measure the redshift-space power spectrum to unprecedented accuracy, allowing for powerful new tests of the LambdaCDM picture as well as measurements of particle physics parameters such as the neutrino masses. We extend the Time-RG perturbative framework to redshift space, computing the power spectrum P_s(k,mu) in massive neutrino cosmologies with time-dependent dark energy equations of state w(z). Time-RG is uniquely capable of incorporating scale-dependent growth into the P_s(k,mu) computation, which is important for massive neutrinos as well as modified gravity models. Although changes to w(z) and the neutrino mass fraction both affect the late-time scale-dependence of the non-linear power spectrum, we find that the two effects depend differently on the line-of-sight angle mu. Finally, we use the HACC N-body code to quantify errors in the perturbative calculations. For a LambdaCDM model at redshift z=1, our procedure predicts the monopole~(quadrupole) to 1% accuracy up to a wave number 0.19h/Mpc (0.28h/Mpc), compared to 0.08h/Mpc (0.07h/Mpc) for the Kaiser approximation and 0.19h/Mpc (0.16h/Mpc) for the current state-of-the-art perturbation scheme. Our calculation agrees with the simulated redshift-space power spectrum even for neutrino masses above the current bound, and for rapidly-evolving dark energy equations of state, |dw/dz| ~ 1. Along with this article, we make our redshift-space Time-RG implementation publicly available as the code redTime.

We characterize infrared spectral energy distributions of 343 (Ultra) Luminous Infrared Galaxies from $z=0.3-2.8$. We diagnose the presence of an AGN by decomposing individual Spitzer mid-IR spectroscopy into emission from star-formation and an AGN-powered continuum; we classify sources as star-forming galaxies (SFGs), AGN, or composites. Composites comprise 30% of our sample and are prevalent at faint and bright $S_{24}$, making them an important source of IR AGN emission. We combine spectroscopy with multiwavelength photometry, including Herschel imaging, to create three libraries of publicly available templates (2-1000 $\mu$m). We fit the far-IR emission using a two temperature modified blackbody to measure cold and warm dust temperatures ($T_c$ and $T_w$). We find that $T_c$ does not depend on mid-IR classification, while $T_w$ shows a notable increase as the AGN grows more luminous. We measure a quadratic relationship between mid-IR AGN emission and total AGN contribution to $L_{\rm IR}$. AGN, composites, and SFGs separate in $S_8/S_{3.6}$ and $S_{250}/S_{24}$, providing a useful diagnostic for estimating relative amounts of these sources. We estimate that >40% of IR selected samples host an AGN, even at faint selection thresholds ($S_{24}$>100 $\mu$Jy). Our decomposition technique and color diagnostics are relevant given upcoming observations with the James Webb Space Telescope.

Recent observations with the GISMO (Goddard-IRAM Superconducting 2 Millimeter Observer) 2 mm camera revealed a detection 8 arcsec away from the lensed galaxy MACS1149-JD1 at z=9.6. Within the 17.5 arcsec FWHM GISMO beam, this detection is consistent with the position of the high-redshift galaxy and therefore, if confirmed, this object could be claimed to be the youngest galaxy producing significant quantities of dust. We present higher resolution (8.5 arcsec) observations of this system taken with the AzTEC 1.1 mm camera mounted on the Large Millimeter Telescope Alfonso Serrano. Dust continuum emission at the position of MACS1149-JD1 is not detected with an r.m.s. of 0.17 mJy/beam. However, we find a detection ~ 11 arcsec away from MACS1149-JD1, still within the GISMO beam which is consistent with an association to the GISMO source. Combining the AzTEC and GISMO photometry, together with Herschel ancillary data, we derive a z_phot= 0.7-1.6 for the dusty galaxy. We conclude therefore that the GISMO and AzTEC detections are not associated with MACS1149-JD1. From the non-detection of MACS1149-JD1 we derive the following (3 \sigma) upper limits corrected for gravitational lensing magnification and for cosmic microwave background effects: dust mass < 1.6 x 10^7 M_sun, IR luminosity < 8 x 10^10 L_sun, star formation rate < 14 M_sun/yr, and UV attenuation < 2.7 mag. These limits are comparable to those derived for other high-redshift galaxies from deep Atacama Large Millimeter/submillimeter Array (ALMA) observations.

Measuring redshifted CO line emission is an unambiguous method for obtaining an accurate redshift and total cold gas content of optically faint, dusty starburst systems. Here, we report the first successful spectroscopic redshift determination of AzTEC J095942.9+022938 ("COSMOS AzTEC-1"), the brightest 1.1mm continuum source found in the AzTEC/JCMT survey (Scott et al. 2008), through a clear detection of the redshifted CO (4-3) and CO (5-4) lines using the Redshift Search Receiver on the Large Millimeter Telescope. The CO redshift of $z=4.3420\pm0.0004$ is confirmed by the detection of the redshifted 158 micron [C II] line using the Submillimeter Array. The new redshift and Herschel photometry yield $L_{FIR}=(1.1\pm0.1)\times 10^{13} L_\odot$ and $SFR = 1300\, M_\odot$ yr$^{-1}$. Its molecular gas mass derived using the ULIRG conversion factor is $1.4\pm0.2 \times 10^{11} M_\odot$ while the total ISM mass derived from the 1.1mm dust continuum is $3.7\pm0.7 \times 10^{11} M_\odot$ assuming dust temperature of 35 K. Our dynamical mass analysis suggests that the compact gas disk ($r\approx 1.1$ kpc, inferred from dust continuum and SED analysis) has to be nearly face-on, providing a natural explanation for the uncommonly bright, compact stellar light seen by the HST. The [C II] line luminosity $L_{[C~II]} = 7.8\pm1.1 \times 10^9 L_\odot$ is remarkably high, but it is only 0.04 per cent of the total IR luminosity. AzTEC COSMOS-1 and other high redshift sources with a spatially resolved size extend the tight trend seen between [C II]/FIR ratio and $\Sigma_{FIR}$ among IR-bright galaxies reported by Diaz-Santos et al. (2013) by more than an order of magnitude, supporting the explanation that the higher intensity of the IR radiation field is responsible for the "[C II] deficiency" seen among luminous starburst galaxies.

We have discovered an optically rich galaxy cluster at z=1.7089 with star formation occurring in close proximity to the central galaxy. The system, SpARCS104922.6+564032.5, was detected within the Spitzer Adaptation of the red-sequence Cluster Survey, (SpARCS), and confirmed through Keck-MOSFIRE spectroscopy. The rest-frame optical richness of Ngal(500kpc) = 30+/-8 implies a total halo mass, within 500kpc, of ~3.8+/-1.2 x 10^14 Msun, comparable to other clusters at or above this redshift. There is a wealth of ancillary data available, including Canada-France-Hawaii Telescope optical, UKIRT-K, Spitzer-IRAC/MIPS, and Herschel-SPIRE. This work adds submillimeter imaging with the SCUBA2 camera on the James Clerk Maxwell Telescope and near-infrared imaging with the Hubble Space Telescope (HST). The mid/far-infrared (M/FIR) data detect an Ultra-luminous Infrared Galaxy spatially coincident with the central galaxy, with LIR = 6.2+/-0.9 x 10^12 Lsun. The detection of polycyclic aromatic hydrocarbons (PAHs) at z=1.7 in a Spitzer-IRS spectrum of the source implies the FIR luminosity is dominated by star formation (an Active Galactic Nucleus contribution of 20%) with a rate of ~860+/-30 Msun/yr. The optical source corresponding to the IR emission is likely a chain of of > 10 individual clumps arranged as "beads on a string" over a linear scale of 66 kpc. Its morphology and proximity to the Brightest Cluster Galaxy imply a gas-rich interaction at the center of the cluster triggered the star formation. This system indicates that wet mergers may be an important process in forming the stellar mass of BCGs at early times.

The halo occupation distribution (HOD) approach has proven to be an effective method for modeling galaxy clustering and bias. In this approach, galaxies of a given type are probabilistically assigned to individual halos in N-body simulations. In this paper, we present a fast emulator for predicting the fully nonlinear galaxy power spectrum over a range of freely specifiable HOD modeling parameters. The emulator is constructed using results from 100 HOD models run on a large LCDM N-body simulation, with Gaussian Process interpolation applied to a PCA-based representation of the galaxy power spectrum. The total error is currently ~3% (~2% in the simulation and ~1% in the emulation process) from z=1 to z=0, over the considered parameter range. We use the emulator to investigate parametric dependencies in the HOD model, as well as the behavior of galaxy bias as a function of HOD parameters. The emulator is publicly available at http://www.hep.anl.gov/cosmology/CosmicEmu/emu.html.

Ground and space-based sky surveys enable powerful cosmological probes based on measurements of galaxy properties and the distribution of galaxies in the Universe. These probes include weak lensing, baryon acoustic oscillations, abundance of galaxy clusters, and redshift space distortions; they are essential to improving our knowledge of the nature of dark energy. On the theory and modeling front, large-scale simulations of cosmic structure formation play an important role in interpreting the observations and in the challenging task of extracting cosmological physics at the needed precision. These simulations must cover a parameter range beyond the standard six cosmological parameters and need to be run at high mass and force resolution. One key simulation-based task is the generation of accurate theoretical predictions for observables, via the method of emulation. Using a new sampling technique, we explore an 8-dimensional parameter space including massive neutrinos and a variable dark energy equation of state. We construct trial emulators using two surrogate models (the linear power spectrum and an approximate halo mass function). The new sampling method allows us to build precision emulators from just 26 cosmological models and to increase the emulator accuracy by adding new sets of simulations in a prescribed way. This allows emulator fidelity to be systematically improved as new observational data becomes available and higher accuracy is required. Finally, using one LCDM cosmology as an example, we study the demands imposed on a simulation campaign to achieve the required statistics and accuracy when building emulators for dark energy investigations.

Do spatial distributions of dust grains in galaxies have typical forms, as do spatial distributions of stars? We investigate whether or not the distributions resemble uniform foreground screens, as commonly assumed by the high-redshift galaxy community. We use rest-frame infrared, ultraviolet, and H$\alpha$ line luminosities of dust-poor and dusty galaxies at z ~ 0 and z ~ 1 to compare measured H$\alpha$ escape fractions with those predicted by the Calzetti attenuation formula. The predictions, based on UV escape fractions, overestimate the measured H$\alpha$ escape fractions for all samples. The interpretation of this result for dust-poor z ~ 0 galaxies is that regions with ionizing stars have more dust than regions with nonionizing UV-emitting stars. Dust distributions for these galaxies are nonuniform. The interpretation of the overestimates for dusty galaxies at both redshifts is less clear. If the attenuation formula is inapplicable to these galaxies, perhaps the disagreements are unphysical; perhaps dust distributions in these galaxies are uniform. If the attenuation formula does apply, then dusty galaxies have nonuniform dust distributions; the distributions are more uniform than they are in dust-poor galaxies. A broad range of H$\alpha$ escape fractions at a given UV escape fraction for z ~ 1 dusty galaxies, if real, indicates diverse dust morphologies and the implausibility of the screen assumption.

We propose a new method of pushing $Herschel$ to its faintest detection limits using universal trends in the redshift evolution of the far infrared over 24$\mu$m colours in the well-sampled GOODS-North field. An extension to other fields with less multi-wavelength information is presented. This method is applied here to raise the contribution of individually detected $Herschel$ sources to the cosmic infrared background (CIRB) by a factor 5 close to its peak at 250$\mu$m and more than 3 in the 350$\mu$m and 500$\mu$m bands. We produce realistic mock $Herschel$ images of the deep PACS and SPIRE images of the GOODS-North field from the GOODS-$Herschel$ Key Program and use them to quantify the confusion noise at the position of individual sources, i.e., estimate a "local confusion noise". Two methods are used to identify sources with reliable photometric accuracy extracted using 24$\mu$m prior positions. The clean index (CI), previously defined but validated here with simulations, which measures the presence of bright 24$\mu$m neighbours and the photometric accuracy index (PAI) directly extracted from the mock $Herschel$ images. After correction for completeness, thanks to our mock $Herschel$ images, individually detected sources make up as much as 54% and 60% of the CIRB in the PACS bands down to 1.1 mJy at 100$\mu$m and 2.2 mJy at 160$\mu$m and 55, 33, and 13% of the CIRB in the SPIRE bands down to 2.5, 5, and 9 mJy at 250$\mu$m, 350$\mu$m, and 500$\mu$m, respectively. The latter depths improve the detection limits of $Herschel$ by factors of 5 at 250$\mu$m, and 3 at 350$\mu$m and 500$\mu$m as compared to the standard confusion limit. Interestingly, the dominant contributors to the CIRB in all $Herschel$ bands appear to be distant siblings of the Milky Way ($z$$\sim$0.96 for $\lambda$$<$300$\mu$m) with a stellar mass of $M_{\star}$$\sim$9$\times$10$^{10}$M$_{\odot}$.

We present a measurement of the star formation properties of a uniform sample of mid-IR selected, unobscured and obscured quasars (QSO1s and QSO2s) in the Bo\"otes survey region. We use an spectral energy distribution (SED) analysis for photometric data spanning optical to far-IR wavelengths to decompose AGN and host galaxy components. We find that when compared to a matched sample of QSO1s, the QSO2s have higher far-IR detection fractions, far-IR fluxes and infrared star formation luminosities ($L_{\rm IR}^{\rm SF}$) by a factor of $\sim2$. Correspondingly, we show that the AGN obscured fraction rises from 0.3 to 0.7 between $4-40\times10^{11}L_\odot$. We also find evidence associating the absorption in the X-ray emission with the presence of far-IR emitting dust. Overall, these results are consistent with galaxy evolution models in which quasar obscurations can be associated with a dust-enshrouded starburst galaxies.

We measure the star formation rates (SFRs) of massive ($M_{\star}>10^{10.1}M_{\odot}$) early-type galaxies (ETGs) in a sample of 11 high-redshift ($1.0 < z < 1.5$) galaxy clusters drawn from the IRAC Shallow Cluster Survey (ISCS). We identify ETGs visually from Hubble Space Telescope imaging and select likely cluster members as having either an appropriate spectroscopic redshift or red sequence color. Mid-infrared SFRs are measured using Spitzer 24 $\mu$m data for isolated cluster galaxies for which contamination by neighbors, and active galactic nuclei, can be ruled out. Cluster ETGs show enhanced specific star formation rates (sSFRs) compared to cluster galaxies in the local Universe, but have sSFRs more than four times lower than that of field ETGs at $1 < z < 1.5$. Relative to the late-type cluster population, isolated ETGs show substantially quenched mean SFRs, yet still contribute 12% of the overall star formation activity measured in $1 < z < 1.5$ clusters. We find that new ETGs are likely being formed in ISCS clusters; the fraction of cluster galaxies identified as ETGs increases from 34% to 56% from $z \sim 1.5 \rightarrow 1.25$. While the fraction of cluster ETGs that are highly star-forming ($\textrm{SFR}\geq26\ M_{\odot}$ yr$^{-1}$) drops from 27% to 10% over the same period, their sSFRs are roughly constant. All these factors taken together suggest that, particularly at $z\gtrsim1.25$, the events that created these distant cluster ETGs$-$likely mergers, at least among the most massive$-$were both recent and gas-rich.

Modeling large-scale sky survey observations is a key driver for the continuing development of high resolution, large-volume, cosmological simulations. We report the first results from the 'Q Continuum' cosmological N-body simulation run carried out on the GPU-accelerated supercomputer Titan. The simulation encompasses a volume of (1300 Mpc)^3 and evolves more than half a trillion particles, leading to a particle mass resolution of ~1.5 X 10^8 M_sun. At this mass resolution, the Q Continuum run is currently the largest cosmology simulation available. It enables the construction of detailed synthetic sky catalogs, encompassing different modeling methodologies, including semi-analytic modeling and sub-halo abundance matching in a large, cosmological volume. Here we describe the simulation and outputs in detail and present first results for a range of cosmological statistics, such as mass power spectra, halo mass functions, and halo mass-concentration relations for different epochs. We also provide details on challenges connected to running a simulation on almost 90% of Titan, one of the fastest supercomputers in the world, including our usage of Titan's GPU accelerators.

The molecular gas, H$_2$, that fuels star formation in galaxies is difficult to observe directly. As such, the ratio of $L_{\rm IR}$ to $L^\prime_{\rm CO}$ is an observational estimation of the star formation rate compared with the amount of molecular gas available to form stars, which is related to the star formation efficiency and the inverse of the gas consumption timescale. We test what effect an IR luminous AGN has on the ratio $L_{\rm IR}/L^\prime_{\rm CO}$ in a sample of 24 intermediate redshift galaxies from the 5 mJy Unbiased Spitzer Extragalactic Survey (5MUSES). We obtain new CO(1-0) observations with the Redshift Search Receiver on the Large Millimeter Telescope. We diagnose the presence and strength of an AGN using Spitzer IRS spectroscopy. We find that removing the AGN contribution to $L_{\rm IR}^{\rm tot}$ results in a mean $L_{\rm IR}^{\rm SF}/L^\prime_{\rm CO}$ for our entire sample consistent with the mean $L_{\rm IR}/L^\prime_{\rm CO}$ derived for a large sample of star forming galaxies from $z\sim0-3$. We also include in our comparison the relative amount of polycyclic aromatic hydrocarbon emission for our sample and a literature sample of local and high redshift Ultra Luminous Infrared Galaxies and find a consistent trend between $L_{6.2}/L_{\rm IR}^{\rm SF}$ and $L_{\rm IR}^{\rm SF}/L^\prime_{\rm CO}$, such that small dust grain emission decreases with increasing $L_{\rm IR}^{\rm SF}/L^\prime_{\rm CO}$ for both local and high redshift dusty galaxies.

Current and future surveys of large-scale cosmic structure are associated with a massive and complex datastream to study, characterize, and ultimately understand the physics behind the two major components of the 'Dark Universe', dark energy and dark matter. In addition, the surveys also probe primordial perturbations and carry out fundamental measurements, such as determining the sum of neutrino masses. Large-scale simulations of structure formation in the Universe play a critical role in the interpretation of the data and extraction of the physics of interest. Just as survey instruments continue to grow in size and complexity, so do the supercomputers that enable these simulations. Here we report on HACC (Hardware/Hybrid Accelerated Cosmology Code), a recently developed and evolving cosmology N-body code framework, designed to run efficiently on diverse computing architectures and to scale to millions of cores and beyond. HACC can run on all current supercomputer architectures and supports a variety of programming models and algorithms. It has been demonstrated at scale on Cell- and GPU-accelerated systems, standard multi-core node clusters, and Blue Gene systems. HACC's design allows for ease of portability, and at the same time, high levels of sustained performance on the fastest supercomputers available. We present a description of the design philosophy of HACC, the underlying algorithms and code structure, and outline implementation details for several specific architectures. We show selected accuracy and performance results from some of the largest high resolution cosmological simulations so far performed, including benchmarks evolving more than 3.6 trillion particles.

Over the next decade, cosmological measurements of the large-scale structure of the Universe will be sensitive to the combined effects of dynamical dark energy and massive neutrinos. The matter power spectrum is a key repository of this information. We extend higher-order perturbative methods for computing the power spectrum to investigate these effects over quasi-linear scales. Through comparison with N-body simulations we establish the regime of validity of a Time-Renormalization Group (Time-RG) perturbative treatment that includes dynamical dark energy and massive neutrinos. We also quantify the accuracy of Standard (SPT), Renormalized (RPT) and Lagrangian Resummation (LPT) perturbation theories without massive neutrinos. We find that an approximation that neglects neutrino clustering as a source for nonlinear matter clustering predicts the Baryon Acoustic Oscillation (BAO) peak position to 0.25% accuracy for redshifts 1 < z < 3, justifying the use of LPT for BAO reconstruction in upcoming surveys. We release a modified version of the public Copter code which includes the additional physics discussed in the paper.

We report the detection of 6.2um polycyclic aromatic hydrocarbon (PAH) and rest-frame 4-7um continuum emission in the z=4.055 submillimeter galaxy GN20, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. This represents the first detection of PAH emission at z>4. The strength of the PAH emission feature is consistent with a very high star formation rate of ~1600Msun/yr. We find that this intense starburst powers at least ~1/3 of the faint underlying 6um continuum emission, with an additional, significant (and perhaps dominant) contribution due to a power-law-like hot dust source, which we interpret to likely be a faint, dust-obscured active galactic nucleus (AGN). The inferred 6um AGN continuum luminosity is consistent with a sensitive upper limit on the hard X-ray emission as measured by the Chandra X-Ray Observatory if the previously undetected AGN is Compton-thick. This is in agreement with the finding at optical/infrared wavelengths that the galaxy and its nucleus are heavily dust-obscured. Despite the strong power-law component enhancing the mid-infrared continuum emission, the intense starburst associated with the photon-dominated regions that give rise to the PAH emission appears to dominate the total energy output in the infrared. GN20 is one of the most luminous starburst galaxies known at any redshift, embedded in a rich protocluster of star-forming galaxies. This investigation provides an improved understanding of the energy sources that power such exceptional systems, which represent the extreme end of massive galaxy formation at early cosmic times.

The use of submm dust continuum emission to probe the mass of interstellar dust and gas in galaxies is empirically calibrated using samples of local star forming galaxies, Planck observations of the Milky Way and high redshift submm galaxies (SMGs). All of these objects suggest a similar calibration, strongly supporting the view that the Rayleigh-Jeans (RJ) tail of the dust emission can be used as an accurate and very fast probe of the ISM in galaxies. We present ALMA Cycle 0 observations of the Band 7 (350 GHz) dust emission in 107 galaxies from z = 0.2 to 2.5. Three samples of galaxies with a total of 101 galaxies were stellar mass-selected from COSMOS to have $M_* \simeq10^{11}$\msun: 37 at z$\sim0.4$, 33 at z$\sim0.9$ and 31 at z$=2$. A fourth sample with 6 IR luminous galaxies at z = 2 was observed for comparison with the purely mass-selected samples. From the fluxes detected in the stacked images for each sample, we find that the ISM content has decreased a factor $\sim 6$ from $1 - 2 \times 10^{10}$\msun at both z = 2 and 0.9 down to $\sim 2 \times 10^9$\msun at z = 0.4. The IR luminous sample at z = 2 shows a further $\sim 4$ times increase in M$_{ISM}$ compared to the equivalent non-IR bright sample at the same redshift. The gas mass fractions are $\sim 2\pm0.5, 12\pm3, 14\pm2 ~\rm{and} ~53\pm3$ $%$ for the four subsamples (z = 0.4, 0.9, 2 and IR bright galaxies).

We analyze the star formation properties of 16 infrared-selected, spectroscopically confirmed galaxy clusters at $1 < z < 1.5$ from the Spitzer/IRAC Shallow Cluster Survey (ISCS). We present new spectroscopic confirmation for six of these high-redshift clusters, five of which are at $z>1.35$. Using infrared luminosities measured with deep Spitzer/MIPS observations at 24 $\mu$m, along with robust optical+IRAC photometric redshifts and SED-fitted stellar masses, we present the dust-obscured star-forming fractions, star formation rates and specific star formation rates in these clusters as functions of redshift and projected clustercentric radius. We find that $z\sim 1.4$ represents a transition redshift for the ISCS sample, with clear evidence of an unquenched era of cluster star formation at earlier times. Beyond this redshift the fraction of star-forming cluster members increases monotonically toward the cluster centers. Indeed, the specific star formation rate in the cores of these distant clusters is consistent with field values at similar redshifts, indicating that at $z>1.4$ environment-dependent quenching had not yet been established in ISCS clusters. Combining these observations with complementary studies showing a rapid increase in the AGN fraction, a stochastic star formation history, and a major merging episode at the same epoch in this cluster sample, we suggest that the starburst activity is likely merger-driven and that the subsequent quenching is due to feedback from merger-fueled AGN. The totality of the evidence suggests we are witnessing the final quenching period that brings an end to the era of star formation in galaxy clusters and initiates the era of passive evolution.

We compare the star formation (SF) activity in cluster galaxies to the field from z=0.3-1.5 using $Herschel$ SPIRE 250$\mu$m imaging. We utilize 274 clusters from the IRAC Shallow Cluster Survey (ISCS) selected as rest-frame near-infrared overdensities over the 9 square degree Bootes field . This analysis allows us to quantify the evolution of SF in clusters over a long redshift baseline without bias against active cluster systems. Using a stacking analysis, we determine the average star formation rates (SFRs) and specific-SFRs (SSFR=SFR/M$_{\star}$) of stellar mass-limited (M>1.3x10$^{10}$ M$_{\odot}$), statistical samples of cluster and field galaxies, probing both the star forming and quiescent populations. We find a clear indication that the average SF in cluster galaxies is evolving more rapidly than in the field, with field SF levels at z>1.2 in the cluster cores (r<0.5 Mpc), in good agreement with previous ISCS studies. By quantifying the SF in cluster and field galaxies as an exponential function of cosmic time, we determine that cluster galaxies are evolving ~2 times faster than the field. Additionally, we see enhanced SF above the field level at z~1.4 in the cluster outskirts (r>0.5 Mpc). These general trends in the cluster cores and outskirts are driven by the lower mass galaxies in our sample. Blue cluster galaxies have systematically lower SSFRs than blue field galaxies, but otherwise show no strong differential evolution with respect to the field over our redshift range. This suggests that the cluster environment is both suppressing the star formation in blue galaxies on long time-scales and rapidly transitioning some fraction of blue galaxies to the quiescent galaxy population on short time-scales. We argue that our results are consistent with both strangulation and ram pressure stripping acting in these clusters, with merger activity occurring in the cluster outskirts.

We explore the relationship between gas, dust and star formation in a sample of 12 ultra-luminous infrared galaxies (ULIRGs) at high redshift compared to a similar sample of local galaxies. We present new CO observations and/or Spitzer mid-IR spectroscopy for 6 70 micron selected galaxies at z~1 in order to quantify the properties of the molecular gas reservoir, the contribution of an active galactic nuclei (AGN) to the mid-IR luminosity and the star formation efficiency (SFE=LIR/L'CO). The mid-IR spectra show strong polycyclic aromatic hydrocarbon (PAH) emission and our spectral decomposition suggests that the AGN makes a minimal contribution (<25%) to the mid-IR luminosity. The 70 micron selected ULIRGs which we find to be spectroscopic close pairs, are observed to have high SFE, similar to local ULIRGs and high redshift submillimeter galaxies, consistent with enhanced IR luminosity due to an ongoing major merger. Combined with existing observations of local and high redshift ULIRGs, we further compare the PAH, IR and CO luminosities. We show that the ratio LPAH6.2/LIR decreases with increasing IR luminosity for both local and high redshift galaxies but the trend for high redshift galaxies is shifted to higher IR luminosities; the average LPAH6.2/LIR ratio at a given LIR is ~3 times higher at high redshift. When we normalize by the molecular gas, we find this trend to be uniform for galaxies at all redshifts and that the molecular gas is correlated with the PAH dust emission.The similar trends seen in the [CII] to molecular gas ratios in other studies suggests that PAH emission, like [CII], continues to be a good tracer of photodissociation regions even at high redshift. Together the CO, PAH and far-IR fine structure lines should be useful for constraining the interstellar medium conditions in high redshift galaxies.

We present a measurement of the average supermassive black hole accretion rate (BHAR) as a function of star formation rate (SFR) for galaxies in the redshift range 0.25<z<0.8. We study a sample of 1,767 far-IR selected star-forming galaxies in the 9 deg^2 Bo\"otes multiwavelength survey field. The SFR is estimated using 250 micron observations from the Herschel Space Observatory, for which the contribution from the AGN is minimal. In this sample, 121 AGNs are directly identified using X-ray or mid-IR selection criteria. We combined these detected AGNs and an X-ray stacking analysis for undetected sources to study the average BHAR for all of the star-forming galaxies in our sample. We find an almost linear relation between the average BHAR (in M_sun/year) and the SFR (in M_sun/year) for galaxies across a wide SFR range 0.85<log SFR<2.56 : log BHAR=(-3.72\pm0.52)+(1.05\pm0.33) log SFR. This global correlation between SFR and average BHAR is consistent with a simple picture in which SFR and AGN activity are tightly linked over galaxy evolution timescales.

Current X-ray surveys have proved to be essential tools in order to identify and study AGNs across cosmic time. However, there is evidence that the most heavily obscured AGNs are largely missing even in the deepest surveys. The search for these obscured AGNs is one of the most outstanding issues of extragalactic astronomy, since they are expected to make a major contribution to the high energy peak of the X-ray background (XRB) and might constitute a particularly active and dusty phase of black hole and galaxy evolution. Using a newly developed SED fitting technique to decompose the AGN from the star-formation emission in the infrared band we identify a sample of 17 IR bright quasars at z=1-3. For the majority of these sources the X-ray spectrum is well characterised by an absorbed power law model, revealing that ~25% of the sources are unabsorbed; the remainder of these IR quasars have moderate-to-high absorption (N_H>10^23 cm^-2), and ~25% are not detected in the X-rays and are likely to be Compton thick (N_H>10^24 cm^-2). We therefore find a much higher fraction of obscured than unobscured quasars, indicating that a large fraction of the luminous black hole accretion was very heavily obscured at z~2.

We present the design and performance of the multi-object fiber spectrographs for the Sloan Digital Sky Survey (SDSS) and their upgrade for the Baryon Oscillation Spectroscopic Survey (BOSS). Originally commissioned in Fall 1999 on the 2.5-m aperture Sloan Telescope at Apache Point Observatory, the spectrographs produced more than 1.5 million spectra for the SDSS and SDSS-II surveys, enabling a wide variety of Galactic and extra-galactic science including the first observation of baryon acoustic oscillations in 2005. The spectrographs were upgraded in 2009 and are currently in use for BOSS, the flagship survey of the third-generation SDSS-III project. BOSS will measure redshifts of 1.35 million massive galaxies to redshift 0.7 and Lyman-alpha absorption of 160,000 high redshift quasars over 10,000 square degrees of sky, making percent level measurements of the absolute cosmic distance scale of the Universe and placing tight constraints on the equation of state of dark energy. The twin multi-object fiber spectrographs utilize a simple optical layout with reflective collimators, gratings, all-refractive cameras, and state-of-the-art CCD detectors to produce hundreds of spectra simultaneously in two channels over a bandpass covering the near ultraviolet to the near infrared, with a resolving power R = \lambda/FWHM ~ 2000. Building on proven heritage, the spectrographs were upgraded for BOSS with volume-phase holographic gratings and modern CCD detectors, improving the peak throughput by nearly a factor of two, extending the bandpass to cover 360 < \lambda < 1000 nm, and increasing the number of fibers from 640 to 1000 per exposure. In this paper we describe the original SDSS spectrograph design and the upgrades implemented for BOSS, and document the predicted and measured performances.

We present a large-scale study of the X-ray properties and near-IR-to-radio SEDs of submillimetre galaxies (SMGs) detected at 1.1mm with the AzTEC instrument across a ~1.2 square degree area of the sky. Combining deep 2-4 Ms Chandra data with Spitzer IRAC/MIPS and VLA data within the GOODS-N/S and COSMOS fields, we find evidence for AGN activity in ~14 percent of 271 AzTEC SMGs, ~28 percent considering only the two GOODS fields. Through X-ray spectral modeling and SED fitting using Monte Carlo Markov Chain techniques to Siebenmorgen et al. (2004) (AGN) and Efstathiou et al. (2000) (starburst) templates, we find that while star formation dominates the IR emission, with SFRs ~100-1000 M_sun/yr, the X-ray emission for most sources is almost exclusively from obscured AGNs, with column densities in excess of 10^23 cm^-2. Only for ~6 percent of our sources do we find an X-ray-derived SFR consistent with NIR-to-radio SED derived SFRs. Inclusion of the X-ray luminosities as a prior to the NIR-to-radio SED effectively sets the AGN luminosity and SFR, preventing significant contribution from the AGN template. Our SED modeling further shows that the AGN and starburst templates typically lack the required 1.1 mm emission necessary to match observations, arguing for an extended, cool dust component. The cross correlation function between the full samples of X-ray sources and SMGs in these fields does not indicate a strong correlation between the two populations at large scales, suggesting that SMGs and AGNs do not necessarily trace the same underlying large scale structure. Combined with the remaining X-ray-dim SMGs, this suggests that sub-mm bright sources may evolve along multiple tracks, with X-ray-detected SMGs representing transitionary objects between periods of high star formation and AGN activity while X-ray-faint SMGs represent a brief starburst phase of more normal galaxies.

We present a new submm/mm galaxy counterpart identification technique which builds on the use of Spitzer IRAC colors as discriminators between likely counterparts and the general IRAC galaxy population. Using 102 radio- and SMA-confirmed counterparts to AzTEC sources across three fields (GOODS-N, GOODS-S, and COSMOS), we develop a non-parametric IRAC color-color characteristic density distribution (CDD), which, when combined with positional uncertainty information via likelihood ratios, allows us to rank all potential IRAC counterparts around SMGs and calculate the significance of each ranking via the reliability factor. We report all robust and tentative radio counterparts to SMGs, the first such list available for AzTEC/COSMOS, as well as the highest ranked IRAC counterparts for all AzTEC SMGs in these fields as determined by our technique. We demonstrate that the technique is free of radio bias and thus applicable regardless of radio detections. For observations made with a moderate beamsize (~18"), this technique identifies ~85 per cent of SMG counterparts. For much larger beamsizes (>30"), we report identification rates of 33-49 per cent. Using simulations, we demonstrate that this technique is an improvement over using positional information alone for observations with facilities such as AzTEC on the LMT and SCUBA-2 on JCMT.

We have compiled a large sample of 151 high redshift (z=0.5-4) galaxies selected at 24 microns (S24>100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-infrared spectrum into contributions from star formation and activity in the galactic nuclei. In addition, we have a wealth of photometric data from Spitzer IRAC/MIPS and Herschel PACS/SPIRE. We explore how effective different infrared color combinations are at separating our mid-IR spectroscopically determined active galactic nuclei from our star forming galaxies. We look in depth at existing IRAC color diagnostics, and we explore new color-color diagnostics combining mid-IR, far-IR, and near-IR photometry, since these combinations provide the most detail about the shape of a source's IR spectrum. An added benefit of using a color that combines far-IR and mid-IR photometry is that it is indicative of the power source driving the IR luminosity. For our data set, the optimal color selections are S250/S24 vs. S8.0/S3.6 and S100/S24 vs. S8.0/S3.6; both diagnostics have ~10% contamination rate in the regions occupied primarily by star forming galaxies and active galactic nuclei, respectively. Based on the low contamination rate, these two new IR color-color diagnostics are ideal for estimating both the mid-IR power source of a galaxy when spectroscopy is unavailable and the dominant power source contributing to the IR luminosity. In the absence of far-IR data, we present color diagnostics using the WISE mid-IR bands which can efficiently select out high z (z~2) star forming galaxies.

Remarkable observational advances have established a compelling cross-validated model of the Universe. Yet, two key pillars of this model -- dark matter and dark energy -- remain mysterious. Sky surveys that map billions of galaxies to explore the `Dark Universe', demand a corresponding extreme-scale simulation capability; the HACC (Hybrid/Hardware Accelerated Cosmology Code) framework has been designed to deliver this level of performance now, and into the future. With its novel algorithmic structure, HACC allows flexible tuning across diverse architectures, including accelerated and multi-core systems. On the IBM BG/Q, HACC attains unprecedented scalable performance -- currently 13.94 PFlops at 69.2% of peak and 90% parallel efficiency on 1,572,864 cores with an equal number of MPI ranks, and a concurrency of 6.3 million. This level of performance was achieved at extreme problem sizes, including a benchmark run with more than 3.6 trillion particles, significantly larger than any cosmological simulation yet performed.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

We present here a new spectral energy distribution (SED) fitting approach that we adopt to select radio-excess sources amongst distant star-forming galaxies in the GOODS-Herschel (North) field and to reveal the presence of hidden, highly obscured AGN. Through extensive SED analysis of 458 galaxies with radio 1.4 GHz and mid-IR 24 um detections using some of the deepest Chandra X-ray, Spitzer and Herschel infrared, and VLA radio data available to date, we have robustly identified a sample of 51 radio-excess AGN (~1300 deg^-2) out to redshift z~3. These radio-excess AGN have a significantly lower far-IR/radio ratio (q<1.68) than the typical relation observed for star-forming galaxies (q~2.2). We find that ~45% of these radio-excess sources have a dominant AGN component in the mid-IR band, while for the remainders the excess radio emission is the only indicator of AGN activity. The fraction of radio-excess AGN increases with X-ray luminosity reaching ~60% at Lx~10^44-10^45 erg/s, making these sources an important part of the total AGN population. However, almost half (24/51) of these radio-excess AGN are not detected in the deep Chandra X-ray data, suggesting that some of these sources might be heavily obscured. We also find that the specific star formation rates (sSFRs) of the radio-excess AGN are on average lower that those observed for X-ray selected AGN hosts, indicating that our sources are forming stars more slowly than typical AGN hosts, and possibly their star formation is progressively quenching.

We explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared (0.3-1000 microns) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We have compiled a large sample of 151 galaxies selected at 24 microns (S24 > 100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-IR spectrum into contributions from star formation and AGN activity. A significant portion (~25%) of our sample is dominated by an AGN in the mid-IR. Based on the mid-IR classification, we divide our full sample into four sub-samples: z~1 star-forming (SF) sources; z~2 SF sources; AGN with clear 9.7 micron silicate absorption; and AGN with featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, including deep Herschel imaging at 100, 160, 250, 350, and 500 microns, we use 95 galaxies with complete spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked difference in the warmer dust temperatures. We calculate the average effective temperature of the dust in each sub-sample and find a significant (~20 K) difference between the SF and AGN systems. We compare our composite SEDs to local templates and find that local templates do not accurately reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR luminous galaxies contain significantly more cool dust than their local counterparts. We find that a full suite of photometry spanning the IR peak is necessary to accurately account for the dominant dust temperature components in high redshift IR luminous galaxies.

Dusty galaxies at z ~ 2 span a wide range of relative brightness between rest-frame mid-infrared (8um) and ultraviolet wavelengths. We attempt to determine the physical mechanism responsible for this diversity. Dust-obscured galaxies (DOGs), which have rest-frame mid-IR to UV flux density ratios > 1000, might be abnormally bright in the mid-IR, perhaps due to prominent AGN and/or PAH emission, or abnormally faint in the UV. We use far-infrared data from the GOODS-Herschel survey to show that most DOGs with 10^12 L_Sun < L_IR < 10^13 L_Sun are not abnormally bright in the mid-IR when compared to other dusty galaxies with similar IR (8--1000um) luminosities. We observe a relation between the median IR to UV luminosity ratios and the median UV continuum power-law indices for these galaxies, and we find that only 24% have specific star formation rates which indicate the dominance of compact star-forming regions. This circumstantial evidence supports the idea that the UV- and IR-emitting regions in these galaxies are spatially coincident, which implies a connection between the abnormal UV faintness of DOGs and dust obscuration. We conclude that the range in rest-frame mid-IR to UV flux density ratios spanned by dusty galaxies at z ~ 2 is due to differing amounts of UV obscuration. Of galaxies with these IR luminosities, DOGs are the most obscured. We attribute differences in UV obscuration to either: 1) differences in the degree of alignment between the spatial distributions of dust and massive stars, or 2) differences in the total dust content.

We report the Herschel SPIRE detection of dust emission arising from UV-luminous (L>~L*) star-forming galaxies at 3.3<z<4.3. Our sample of 1,913 Lyman Break Galaxy (LBG) candidates is selected over an area of 5.3 deg2 in the Bootes Field of the NOAO Deep Wide-Field Survey. This is one of the largest samples of UV-luminous galaxies at this epoch and enables an investigation of the bright end of the galaxy luminosity function. We divide our sample into three luminosity bins and stack the Herschel SPIRE data to measure the average spectral energy distribution (SED) of LBGs at far-infrared (FIR) wavelengths. We find that these galaxies have average IR luminosities of (3-5)x10^{11}L_sun and 60-70% of their star-formation obscured by dust. The FIR SEDs peak at lambda>100um, suggesting dust temperatures (T_d=27-30K) significantly colder than that of local galaxies of comparable IR luminosities. The observed IR-to-UV luminosity ratio (IRX=L_IR/L_UV) is low (~3-4) compared with that observed for z~2 LBGs (IRX(z~2)~7.1+/-1.1). The correlation between the slope of the UV continuum and IRX for galaxies in the two lower luminosity bins suggests dust properties similar to those of local starburst galaxies. However, the galaxies in the highest luminosity bin appear to deviate from the local relation, suggesting that their dust properties may differ from those of their lower-luminosity and low-redshift counterparts. We speculate that the most UV luminous galaxies at this epoch are being observed in a short-lived and young evolutionary phase.

Using deep 100-160 micron observations in GOODS-S from the GOODS-H survey, combined with HST/WFC3 NIR imaging from CANDELS, we present the first morphological analysis of a complete, FIR selected sample of 52 ULIRGs at z~2. We also make use of a comparison sample of galaxies without Herschel detections but with the same z and magnitude distribution. Our visual classifications of these two samples indicate that the fraction of objects with disk and spheroid morphologies is roughly the same but that there are significantly more mergers, interactions, and irregular galaxies among the ULIRGs. The combination of disk and irregular/interacting morphologies suggests that early stage interactions and minor mergers could play an important role in ULIRGs at z~2. We compare these fractions with those of a z~1 sample across a wide luminosity range and find that the fraction of disks decreases systematically with L_IR while the fraction of mergers and interactions increases, as has been observed locally. At comparable luminosities, the fraction of ULIRGs with various morphological classifications is similar at z~2 and z~1. We investigate the position of the ULIRGs, along with 70 LIRGs, on the specific star formation rate versus redshift plane, and find 52 systems to be starbursts (lie more than a factor of 3 above the main sequence relation). The morphologies of starbursts are dominated by interacting and merging systems (50%). If irregular disks are included as potential minor mergers, then we find that up to 73% of starbursts are involved in a merger or interaction at some level. Although the final coalescence of a major merger may not be required for the high luminosities of ULIRGs at z~2 as is the case locally, the large fraction of interactions at all stages and potential minor mergers suggest that the high star formation rates of ULIRGs are still largely externally triggered at z~2.

Using a sample of BzK-selected galaxies at z~2 identified from the CFHT/WIRCAM near-infrared survey of GOODS-North, we discuss the relation between star formation rate (SFR), specific star formation rate (SSFR), and stellar mass (M_{*}), and the clustering of galaxies as a function of these parameters. For star-forming galaxies (sBzKs), the UV-based SFR, corrected for extinction, scales with the stellar mass as SFR ~ M_{*}^{alpha} with alpha = 0.74+/-0.20 down to M_{*} ~ 10^{9} M_{solar}, indicating a weak dependence on the stellar mass of the SSFR. We also measure the angular correlation function and hence infer the correlation length for sBzK galaxies as a function of M_{*}, SFR, and SSFR, as well as K-band apparent magnitude. We show that passive galaxies (pBzKs) are more strongly clustered than sBzK galaxies at a given stellar mass, mirroring the color-density relation seen at lower redshifts. We also find that the correlation length of sBzK galaxies ranges from 4 to 20 h^{-1}Mpc, being a strong function of M_{K}, M_{*}, and SFR. On the other hand, the clustering dependence on SSFR changes abruptly at 2x10^{-9} yr^{-1}, which is the typical value for "main sequence" star-forming galaxies at z~2. We show that the correlation length reaches a minimum at this characteristic value, and is larger for galaxies with both smaller and larger SSFRs; a dichotomy that is only marginally implied from the predictions of the semi-analytical models. Our results suggest that there are two types of environmental effects at work at z~2. Stronger clustering for relatively quiescent galaxies implies that the environment has started to play a role in quenching star formation. At the same time, stronger clustering for galaxies with elevated SSFRs ("starbursts") might be attributed to an increased efficiency for galaxy interactions and mergers in dense environments.

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase.

We present observations of CO J=2-1 line emission in infrared-luminous cluster galaxies at z~1 using the IRAM Plateau de Bure Interferometer. Our two primary targets are optically faint, dust-obscured galaxies (DOGs) found to lie within 2 Mpc of the centers of two massive (>10^14 Msun) galaxy clusters. CO line emission is not detected in either DOG. We calculate 3-sigma upper limits to the CO J=2-1 line luminosities, L'_CO < 6.08x10^9 and < 6.63x10^9 K km/s pc^2. Assuming a CO-to-H_2 conversion factor derived for ultraluminous infrared galaxies in the local Universe, this translates to limits on the cold molecular gas mass of M_H_2 < 4.86x10^9 Msun and M_H_2 < 5.30x10^9 Msun. Both DOGs exhibit mid-infrared continuum emission that follows a power-law, suggesting that an AGN contributes to the dust heating. As such, estimates of the star formation efficiencies in these DOGs are uncertain. A third cluster member with an infrared luminosity, L_IR < 7.4x10^11 Lsun, is serendipitously detected in CO J=2-1 line emission in the field of one of the DOGs located roughly two virial radii away from the cluster center. The optical spectrum of this object suggests that it is likely an obscured AGN, and the measured CO line luminosity is L'_CO = (1.94 +/- 0.35)x10^10 K km/s pc^2, which leads to an estimated cold molecular gas mass M_H_2 = (1.55+/-0.28)x10^10 Msun. A significant reservoir of molecular gas in a z~1 galaxy located away from the cluster center demonstrates that the fuel can exist to drive an increase in star-formation and AGN activity at the outskirts of high-redshift clusters.

Dust-obscured galaxies (DOGs) are a subset of high-redshift (z \approx 2) optically-faint ultra-luminous infrared galaxies (ULIRGs, e.g. L_{IR} > 10^{12} Lsun). We present new far-infrared photometry, at 250, 350, and 500 um (observed-frame), from the Herschel Space Telescope for a large sample of 113 DOGs with spectroscopically measured redshifts. Approximately 60% of the sample are detected in the far-IR, confirming their high IR luminosities, which range from 10^{11.6} Lsun < L_{IR} (8-1000 um) <10^{13.6} Lsun. 90% of the Herschel detected DOGs in this sample are ULIRGs and 30% have L_{IR} > 10^{13} Lsun. The rest-frame near-IR (1 - 3 um) SEDs of the Herschel detected DOGs are predictors of their SEDs at longer wavelengths. DOGs with "power-law" SEDs in the rest-frame near-IR show observed-frame 250/24 um flux density ratios similar to the QSO-like local ULIRG, Mrk 231. DOGs with a stellar "bump" in their rest-frame near-IR show observed-frame 250/24 um flux density ratios similar to local star-bursting ULIRGs like NGC 6240. For the Herschel detected DOGs, accurate estimates (within \approx 25%) of total IR luminosity can be predicted from their rest-frame mid-IR data alone (e.g. from Spitzer observed-frame 24 um luminosities). Herschel detected DOGs tend to have a high ratio of infrared luminosity to rest-frame 8 um luminosity (the IR8= L_{IR}(8-1000 um)/v L_{v}(8 um) parameter of Elbaz et al. 2011). Instead of lying on the z=1-2 "infrared main-sequence" of star forming galaxies (like typical LIRGs and ULIRGs at those epochs) the DOGs, especially large fractions of the bump sources, tend to lie in the starburst sequence. While, Herschel detected DOGs are similar to scaled up versions of local ULIRGs in terms of 250/24 um flux density ratio, and IR8, they tend to have cooler far-IR dust temperatures (20-40 K for DOGs vs. 40-50 K for local ULIRGs). Abridged.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

To optimally identify candidates of the Compton-thick (CT) active galactic nuclei (AGN) that contribute to the unresolved X-ray background in infrared surveys, a tracer of column density is desirable in addition to an AGN indicator. In a recent study, we aimed to test whether the 9.7 micron silicate absorption feature can be used for this purpose when seen at high optical depths. We found that the extreme criterion of optical thickness at 9.7 micron is efficient in identifying CT objects among local AGN. Having identified six of the nine CT AGN in the 12 micron sample with Spitzer and X-ray spectra, we expanded this analysis at intermediate/high z, using all GOODS and FLS sources with Spitzer and X-ray observations. We found 12 sources with tau_9.7>1 that host an AGN between 0.8<z<2.7. Four of them are likely to be CT according to their low X-ray to 6 micron luminosity ratio. Surveys with complete coverage in both mid-infrared spectra and X-ray data can provide large populations of such sources, as at least 5-9% of all infrared bright galaxies in the GOODS and FLS samples are tau_9.7>1 AGN.

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

Using extremely deep PACS 100- and 160um Herschel data from the GOODS-Herschel program, we identify 21 infrared bright galaxies previously missed in the deepest 24um surveys performed by MIPS. These MIPS dropouts are predominantly found in two redshift bins, centred at z ~0.4 and ~1.3. Their S_100/S_24 flux density ratios are similar to those of local LIRGs and ULIRGs, whose silicate absorption features at 18um (at z ~ 0.4) and 9.7um (at z ~ 1.3) are shifted into the 24um MIPS band at these redshifts. The high-z sub-sample consists of 11 infrared luminous sources, accounting for ~2% of the whole GOODS-Herschel sample and putting strong upper limits on the fraction of LIRGs/ULIRGs at 1.0<z<1.7 that are missed by the 24um surveys. We find that a S_100/S_24 > 43 colour cut selects galaxies with a redshift distribution similar to that of the MIPS dropouts and when combined with a second colour cut, S_16/S_8 > 4, isolates sources at 1.0 < z < 1.7. We show that these sources have elevated specific star formation rates (sSFR) compared to main sequence galaxies at these redshifts and are likely to be compact starbursts with moderate/strong 9.7um silicate absorption features in their mid-IR spectra. Herschel data reveal that their infrared luminosities extrapolated from the 24um flux density are underestimated, on average, by a factor of ~3. These silicate break galaxies account for 16% (8%) of the ULIRG (LIRG) population in the GOODS fields, indicating a lower limit in their space density of 2.0 \times 10^(-5) Mpc^(-3). Finally, we provide estimates of the fraction of z < 2 MIPS dropout sources as a function of the 24-, 100-, 160-, 250- and 350um sensitivity limits, and conclude that previous predictions of a population of silicate break galaxies missed by the major 24um extragalactic surveys have been overestimated.

Spitzer/IRS has revealed many sources with very deep Si features at 9.7micron (tau>1). We set out to investigate whether a strong Si absorption feature is a good indicator for the presence of a heavily obscured AGN. We compile X-ray spectroscopic observations available in the literature on the optically-thick,tau(9.7)>1 sources from the IRAS Seyfert sample. We find that the majority of the high-tau optically confirmed Seyferts (6/9) in this sample are probably CT. Thus we provide direct evidence for a connection between mid-IR optically-thick galaxies and CT AGN, with the success rate being close to 70% in the local Universe. This is at least comparable, if not better, than other rates obtained with photometric information in the mid to far-IR, or even mid-IR to Xray. However, this technique cannot provide complete CT AGN samples,ie there are many CT AGN which do not show significant Si absorption, with the most notable example being N1068. Having assessed the validity of the high 9.7micron technique locally, we attempt to construct a sample of candidate CT AGN at higher redshifts. We compile a sample of 7 high-tau sources in the GOODS and 5 in the Spitzer FLS. All these have been selected to have no PAH features EW(6.2)<0.3 in order to maximize the probability that they are AGN. 6 out of 7 sources in the GOODS have been detected in X-rays, while for the five FLS sources only X-ray flux upper limits are available. The high X-ray luminosities of the detected GOODS sources corroborates that these are AGN. For FLS, ancillary optical spectroscopy reveals hidden nuclei in two more sources. SED fitting can support the presence of an AGN in the vast majority of sources. We cannot derive useful X-ray spectroscopy constraints on whether these are CT. However, the low LX/L6 ratios, suggest that at least 4 out of the 6 detected sources in GOODS may be associated with CT AGN.

We measure the angular two-point correlation function of sub-millimeter galaxies (SMGs) from 1.1-millimeter imaging of the COSMOS field with the AzTEC camera and ASTE 10-meter telescope. These data yields one of the largest contiguous samples of SMGs to date, covering an area of 0.72 degrees^2 down to a 1.26 mJy/beam (1-sigma) limit, including 189 (328) sources with S/N greater than 3.5 (3). We can only set upper limits to the correlation length r_0, modeling the correlation function as a power-law with pre-assigned slope. Assuming existing redshift distributions, we derive 68.3% confidence level upper limits of r_0 < 6-8 h^-1 Mpc at 3.7 mJy, and r_0 < 11-12 h^-1 Mpc at 4.2 mJy. Although consistent with most previous estimates, these upper limits imply that the real r_0 is likely smaller. This casts doubts on the robustness of claims that SMGs are characterized by significantly stronger spatial clustering, (and thus larger mass), than differently selected galaxies at high-redshift. Using Monte Carlo simulations we show that even strongly clustered distributions of galaxies can appear unclustered when sampled with limited sensitivity and coarse angular resolution common to current sub-millimeter surveys. The simulations, however, also show that unclustered distributions can appear strongly clustered under these circumstances. From the simulations, we predict that at our survey depth, a mapped area of two degrees^2 is needed to reconstruct the correlation function, assuming smaller beam sizes of future surveys (e.g. the Large Millimeter Telescope's 6" beam size). At present, robust measures of the clustering strength of bright SMGs appear to be below the reach of most observations.

(Abridged) We report on an accounting of the star formation and accretion driven energetics of 24um detected sources in GOODS North. For sources having infrared (IR; 8-1000um) luminosities >3x10^12 L_sun when derived by fitting local SEDs to 24um photometry alone, we find these IR luminosity estimates to be a factor of ~4 times larger than those estimated when the SED fitting includes additional 16 and 70um data (and in some cases mid-infrared spectroscopy and 850um data). This discrepancy arises from the fact that high luminosity sources at z>>0 appear to have far- to mid-infrared ratios, as well as aromatic feature equivalent widths, typical of lower luminosity galaxies in the local Universe. Using our improved estimates for IR luminosity and AGN contributions, we investigate the evolution of the IR luminosity density versus redshift arising from star formation and AGN processes alone. We find that, within the uncertainties, the total star formation driven IR luminosity density is constant between 1.15 < z < 2.35, although our results suggest a slightly larger value at z>2. AGN appear to account for <18% of the total IR luminosity density integrated between 0< z < 2.35, contributing <25% at each epoch. LIRG appear to dominate the star formation rate (SFR) density along with normal star-forming galaxies (L_IR < 10^11 L_sun) between 0.6 < z < 1.15. Once beyond z >2, the contribution from ultraluminous infrared galaxies ULIRGs becomes comparable with that of LIRGs. Using our improved IR luminosity estimates, we find existing calibrations for UV extinction corrections based on measurements of the UV spectral slope typically overcorrect UV luminosities by a factor of ~2, on average, for our sample of 24um-selected sources; accordingly we have derived a new UV extinction correction more appropriate for our sample.

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and shallow mapping project at 850 microns, designed to find rare objects, both Galactic and extragalactic. Two distinct sets of exploratory observations were undertaken, and used to test the SASSy approach and data reduction pipeline. The first was a 0.5 by 0.5 degrees map around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but no other convincing sources are present in the map. Comparison with other galaxies with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. The second observations cover 1 square degree around the W5-E HII region. As well as diffuse structure in the map, a filtering approach was able to extract 27 compact sources with signal-to-noise greater than 6. By matching with data at other wavelengths we can see that the SCUBA-2 data can be used to discriminate the colder cores. Together these observations show that the SASSy project will be able to meet its original goals of detecting new bright sources which will be ideal for follow-up observations with other facilities.

The advent of large cosmological sky surveys - ushering in the era of precision cosmology - has been accompanied by ever larger cosmological simulations. The analysis of these simulations, which currently encompass tens of billions of particles and up to trillion particles in the near future, is often as daunting as carrying out the simulations in the first place. Therefore, the development of very efficient analysis tools combining qualitative and quantitative capabilities is a matter of some urgency. In this paper we introduce new analysis features implemented within ParaView, a parallel, open-source visualization toolkit, to analyze large N-body simulations. The new features include particle readers and a very efficient halo finder which identifies friends-of-friends halos and determines common halo properties. In combination with many other functionalities already existing within ParaView, such as histogram routines or interfaces to Python, this enhanced version enables fast, interactive, and convenient analyses of large cosmological simulations. In addition, development paths are available for future extensions.

We present a study of the cosmic infrared background, which is a measure of the dust obscured activity in all galaxies in the Universe. We venture to isolate the galaxies responsible for the background at 1mm; with spectroscopic and photometric redshifts we constrain the redshift distribution of these galaxies. We create a deep 1.16mm map (sigma ~ 0.5mJy) by combining the AzTEC 1.1mm and MAMBO 1.2mm datasets in GOODS-N. This combined map contains 41 secure detections, 13 of which are new. By averaging the 1.16mm flux densities of individually undetected galaxies with 24um flux densities > 25uJy, we resolve 31--45 per cent of the 1.16mm background. Repeating our analysis on the SCUBA 850um map, we resolve a higher percentage (40--64 per cent) of the 850um background. A majority of the background resolved (attributed to individual galaxies) at both wavelengths comes from galaxies at z > 1.3. If the ratio of the resolved submillimeter to millimeter background is applied to a reasonable scenario for the origins of the unresolved submillimeter background, 60--88 per cent of the total 1.16mm background comes from galaxies at z > 1.3.

Submillimetre continuum radiation allows us to probe cold objects, particularly the earliest, dusty phases of star formation, high-redshift galaxies and circumstellar disks. The submillimetre window gives a unique view of the physical and dynamical conditions in the neutral and molecular interstellar medium. In the next decade a combination of wide-field surveys with single-dish telescopes and targeted follow-up with ALMA and other facilities should enable rapid progress in answering questions about the origins of planetary systems, stars and galaxies.

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

We present deep {\it Spitzer} mid-infrared spectroscopy, along with 16, 24, 70, and 850\,$\micron$\ photometry, for 22 galaxies located in the Great Observatories Origins Deep Survey-North (GOODS-N) field. The sample spans a redshift range of $0.6\la z \la 2.6$, 24~$\mu$m flux densities between $\sim$0.2$-$1.2 mJy, and consists of submillimeter galaxies (SMGs), X-ray or optically selected active galactic nuclei (AGN), and optically faint ($z_{AB}>25$\,mag) sources. We find that infrared (IR; $8-1000~\micron$) luminosities derived by fitting local spectral energy distributions (SEDs) with 24~$\micron$ photometry alone are well matched to those when additional mid-infrared spectroscopic and longer wavelength photometric data is used for galaxies having $z\la1.4$ and 24~$\micron$-derived IR luminosities typically $\la 3\times 10^{12}~L_{\sun}$. However, for galaxies in the redshift range between $1.4\la z \la 2.6$, typically having 24~$\micron$-derived IR luminosities $\ga 3\times 10^{12}~L_{\sun}$, IR luminosities are overestimated by an average factor of $\sim$5 when SED fitting with 24~$\micron$ photometry alone. This result arises partly due to the fact that high redshift galaxies exhibit aromatic feature equivalent widths that are large compared to local galaxies of similar luminosities. Through a spectral decomposition of mid-infrared spectroscopic data, we are able to isolate the fraction of IR luminosity arising from an AGN as opposed to star formation activity. This fraction is only able to account for $\sim$30\% of the total IR luminosity among the entire sample.

We explore a technique for identifying the highest redshift (z>4) sources in Herschel/SPIRE and BLAST submillimeter surveys by localizing the position of the far-infrared dust peak. Just as Spitzer/IRAC was used to identify stellar `bump' sources, the far-IR peak is also a redshift indicator; although, the latter also depends on the average dust temperature. We demonstrate the wide range of allowable redshifts for a reasonable range of dust temperatures and show that it is impossible to constraint the redshift of individual objects using solely the position of the far-IR peak. By fitting spectral energy distribution models to simulated Herschel/SPIRE photometry we show the utility of radio and/or far-infrared data in breaking this degeneracy. With prior knowledge of the dust temperature distribution it is possible to obtain statistical samples of high redshift submillimeter galaxy candidates. We apply this technique to the BLAST survey of ECDFS to constrain the number of dusty galaxies at z>4. We find 8 +/- 2 galaxies with flux density ratios of S500>S350; this sets an upper limit of 17 +/- 4 deg-2 if we assume all are at z>4. This is <35% of all 500 micron-selected galaxies down to S500>45 mJy (LIR>2e13Lsun for z>4). Modeling with conventional temperature and redshift distributions estimates the percentage of these 500 micron peak galaxies at z>4 to be between 10-85%. Our results are consistent with other estimates of the number density of very high redshift submillimeter galaxies and follows the decline in the star formation rate density at z>4.

We provide measurements of the integrated galaxy light at 70, 160, 250, 350 and 500 micron using deep far-infrared and submillimeter data from space (Spitzer) and balloon platform (BLAST) extragalactic surveys. We use the technique of stacking at the positions of 24 micron sources, to supplement the fraction of the integrated galaxy light that is directly resolved through direct detections. We demonstrate that the integrated galaxy light even through stacking, falls short by factors of 2-3 in resolving the extragalactic far-infrared background. We also show that previous estimates of the integrated galaxy light (IGL) through stacking, have been biased towards high values. This is primarily due to multiple counting of the far-infrared/submillimeter flux from 24 micron sources which are clustered within the large point spread function of a brighter far-infrared source. Using models for the evolution of the luminosity function at z<1.2 which are constrained by observations at 24 micron and 70 micron, and which are consistent with the results from the stacking analysis, we find that galaxies at z<1.2, account for ~95-55% of the extragalactic far-infrared background in the ~70-500 micron range respectively. This places strong upper limits on the fraction of dust obscured star-formation at z>1, which are remarkably, below the values derived from the extinction corrected ultraviolet luminosities of galaxies. We use the results to make predictions for the nature of galaxies that extragalactic surveys with Herschel Space Observatory will reveal. Finally, from our constraints on the far-infrared IGL, we provide evidence for the existence of ice mantle dust, orbiting the sun at a distance of ~40 AU, which is contributing intensity to both the near- and far-infrared background. [ABRIDGED]

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

We present a set of ultra-large particle-mesh simulations of the LyA forest targeted at understanding the imprint of baryon acoustic oscillations (BAO) in the inter-galactic medium. We use 9 dark matter only simulations which can, for the first time, simultaneously resolve the Jeans scale of the intergalactic gas while covering the large volumes required to adequately sample the acoustic feature. Mock absorption spectra are generated using the fluctuating Gunn-Peterson approximation which have approximately correct flux probability density functions (PDFs) and small-scale power spectra. On larger scales there is clear evidence in the redshift space correlation function for an acoustic feature, which matches a linear theory template with constant bias. These spectra, which we make publicly available, can be used to test pipelines, plan future experiments and model various physical effects. As an illustration we discuss the basic properties of the acoustic signal in the forest, the scaling of errors with noise and source number density, modified statistics to treat mean flux evolution and misestimation, and non-gravitational sources such as fluctuations in the photo-ionizing background and temperature fluctuations due to HeII reionization.

It is widely believed that ultraluminous infrared (IR) galaxies and active galactic nuclei (AGN) activity are triggered by galaxy interactions and merging, with the peak of activity occurring at z~2, where submillimetre galaxies are thousands of times more numerous than local ULIRGs. In this evolutionary picture, submillimetre galaxies (SMGs) would host an AGN, which would eventually grow a black hole (BH) strong enough to blow off all of the gas and dust leaving an optically luminous QSO. To probe this evolutionary sequence we have focussed on the 'missing link' sources, which demonstrate both strong starburst (SB) and AGN signatures, in order to determine if the SB is the main power source even in SMGs when we have evidence that an AGN is present from their IRAC colours. The best way to determine if a dominant AGN is present is to look for their signatures in the mid-infrared with the Spitzer IRS, since often even deep X-ray observations miss identifying the presence of AGN in heavily dust-obscured SMGs. We present the results of our audit of the energy balance between star-formation and AGN within this special sub-population of SMGs -- where the BH has grown appreciably to begin heating the dust emission.

We have employed the GMRT and the VLA to map the Lockman Hole. At 610 and 1,400 MHz, we reach noise levels of 15 and 6 uJy/beam, respectively, with well-matched resolutions (~5"). At this depth we obtained reliable detections for about half of the known submm galaxies (SMGs) in the field. For radio-identified SMGs, which are typically at z ~ 2, we measure a mean radio spectral index of alpha = -0.75 +/- 0.06 and standard deviation of 0.29, between rest-frame ~1.8 and ~4.2 GHz. The slope of their continuum emission is indistinguishable from that of local star-forming galaxies and suggests that extended optically-thin synchrotron emission dominates the radio output of SMGs. Cooling effects by synchrotron emission and Inverse Compton scattering off the CMB do not seem to affect their radio SEDs. For those SMGs judged by Spitzer mid-IR colours and spectroscopy to host obscured AGN, we find a clear deviation from the rest of the sample - they typically have steeper radio spectral indices, alpha ~< -1.0. These findings suggest these mid-IR-/AGN-selected SMGs may have an intrinsically different injection mechanism for relativistic particles, or they might reside in denser environments. This work provides a reliable spectral template for the estimation of far-IR/radio photometric redshifts, and will enable accurate statistical K-corrections for the large samples of SMGs expected with SCUBA-2 and Herschel.

The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample represents the final set of galaxies observed using the original SDSS target selection criteria. We analyse the clustering of galaxies within this sample, including both the Luminous Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift Survey (2dFGRS) data. Baryon Acoustic Oscillations are observed in power spectra measured for different slices in redshift; this allows us to constrain the distance--redshift relation at multiple epochs. We achieve a distance measure at redshift z=0.275, of r_s(z_d)/D_V(0.275)=0.1390+/-0.0037 (2.7% accuracy), where r_s(z_d) is the comoving sound horizon at the baryon drag epoch, D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3), D_A(z) is the angular diameter distance and H(z) is the Hubble parameter. We find an almost independent constraint on the ratio of distances D_V(0.35)/D_V(0.2)=1.736+/-0.065, which is consistent at the 1.1sigma level with the best fit Lambda-CDM model obtained when combining our z=0.275 distance constraint with the WMAP 5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, but the discrepancy is now of lower significance, a change caused by a revised error analysis and a change in the methodology adopted, as well as the addition of more data. Using WMAP5 constraints on Omega_bh^2 and Omega_ch^2, and combining our BAO distance measurements with those from the Union Supernova sample, places a tight constraint on Omega_m=0.286+/-0.018 and H_0 = 68.2+/-2.2km/s/Mpc that is robust to allowing curvature and non-Lambda dark energy. This result is independent of the behaviour of dark energy at redshifts greater than those probed by the BAO and supernova measurements. (abridged)

We present an analysis of the X-ray properties of 35 submm galaxies(SMGs) in the CDF-N. Using a sample of robust 850micron-selected galaxies, with accurate positions from Spitzer and/or radio counterparts, we find 16 objects (45+/-8%) with significant X-ray detections in the 2Ms Chandra data. 6 of these SMGs (~17+/-6%) have measured X-ray luminosities or upper limits consistent with those expected based on the far-infrared (FIR) or radio-derived star formation rate (SFR), and hence with the X-rays coming solely from star formation. In another 7 sources (20+/-7%) a dominant AGN contribution to the X-ray emission is required, while in 3 more it is unclear whether stellar process or accretion are responsible. Stacking of the X-ray undetected SMGs reveals a highly significant detection. If due to star formation, this corresponds to an average X-ray derived SFR of ~150 Msun/yr. We deduce that the AGN fraction in SMGs based on X-ray observations is 20-29 (+/-7) %, which is towards the lower limit of previous estimates. Spectral analysis shows that in general the SMGs are not heavily obscured in the X-ray but most of the SMGs classfied as AGN show absorption with N_H in excess of 10^22 cm^-2. Of the secure AGN, the bolometric luminosity appears to be dominated by the AGN in only 3 cases. In ~85% of the SMGs, the X-ray spectrum effectively rules out an AGN contribution that dominates the bolometric emission, even if the AGN is Compton thick. The evidence therefore suggests that intense star formation accounts for both the FIR and X-ray emission in most SMGs. We argue that, rather than having an especially high AGN fraction or duty cycle, SMGs have a high X-ray detection rate at very faint fluxes partly because of their high star formation rates and, in rarer cases, because the submm emission is from an AGN.(abridged)

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

Using the IRAM Plateau-de-Bure Interferometer (PdBI) we have searched for the upper fine structure line of neutral carbon (CI(2->1), nu_rest = 809 GHz) and 12CO(J=7->6) (nu_rest=806GHz) towards the submillimetre galaxies (SMGs) GN20 (SMM J123711.9+622212, z = 4.055) and GN20.2 (SMM J123708.8+622202, z = 4.051). The far-infrared (FIR) continuum is detected at 8sigma significance in GN20, with a flux density of S_1.8mm = 1.9+-0.2mJy, while no continuum is detected in GN20.2. Both sources are statistically undetected in both CI and CO(7-6) lines; we derive line luminosity limits for both CI and CO of L'<~2x10^10K km s^-1 pc^2. Assuming carbon excitation temperatures of T_ex = 30K (the galaxies' measured dust temperatures), we infer CI mass limits of M_CI < 5.4x10^6 Msun (GN20) and M_CI < 6.8x10^6 Msun (GN20.2). The derived CI abundance limits are <1.8x10^-5 for GN20 and <3.8x10^-5 for GN20.2 implying that the systems have Milky Way level neutral carbon enrichment (X[CI]/X[H2]) or lower, similar to high-redshift carbon-detected systems (at 5x10^-5) but about 50 times less than the neutral carbon enrichment of local starburst galaxies. Observations of GN20 and GN20.2 in high-resolution MERLIN+VLA radio maps of GOODS-N are used to further constrain the sizes and locations of active regions. We conclude that the physical gas properties of young rapidly evolving systems like GN20 and GN20.2 are likely significantly different than starburst/ULIRG environments in the local Universe yet similar to z~2 SMGs. Unless gravitationally amplified examples can be found, observations of galaxies like GN20 will require the order of magnitude increase in sensitivity of the Atacama Large Millimetre Array (ALMA) to constrain their CI and high-J CO content, despite the fact that they are the brightest systems at z~4.

We present the power spectrum of the reconstructed halo density field derived from a sample of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Release (DR7). The halo power spectrum has a direct connection to the underlying dark matter power for k <= 0.2 h/Mpc, well into the quasi-linear regime. This enables us to use a factor of ~8 more modes in the cosmological analysis than an analysis with kmax = 0.1 h/Mpc, as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 2006). The observed halo power spectrum for 0.02 < k < 0.2 h/Mpc is well-fit by our model: chi^2 = 39.6 for 40 degrees of freedom for the best fit LCDM model. We find \Omega_m h^2 * (n_s/0.96)^0.13 = 0.141^{+0.009}_{-0.012} for a power law primordial power spectrum with spectral index n_s and \Omega_b h^2 = 0.02265 fixed, consistent with CMB measurements. The halo power spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to an effective distance to z=0.35: r_s/D_V(0.35) = 0.1097^{+0.0039}_{-0.0042}. Combining the halo power spectrum measurement with the WMAP 5 year results, for the flat LCDM model we find \Omega_m = 0.289 +/- 0.019 and H_0 = 69.4 +/- 1.6 km/s/Mpc. Allowing for massive neutrinos in LCDM, we find \sum m_{\nu} < 0.62 eV at the 95% confidence level. If we instead consider the effective number of relativistic species Neff as a free parameter, we find Neff = 4.8^{+1.8}_{-1.7}. Combining also with the Kowalski et al. (2008) supernova sample, we find \Omega_{tot} = 1.011 +/- 0.009 and w = -0.99 +/- 0.11 for an open cosmology with constant dark energy equation of state w.

We present results from a multi-wavelength study of 29 sources (false detection probabilities <5%) from a survey of the Great Observatories Origins Deep Survey-North field at 1.1mm using the AzTEC camera. Comparing with existing 850um SCUBA studies in the field, we examine differences in the source populations selected at the two wavelengths. The AzTEC observations uniformly cover the entire survey field to a 1-sigma depth of ~1mJy. Searching deep 1.4GHz VLA, and Spitzer 3--24um catalogues, we identify robust counterparts for 21 1.1mm sources, and tentative associations for the remaining objects. The redshift distribution of AzTEC sources is inferred from available spectroscopic and photometric redshifts. We find a median redshift of z=2.7, somewhat higher than z=2.0 for 850um-selected sources in the same field, and our lowest redshift identification lies at a spectroscopic redshift z=1.1460. We measure the 850um to 1.1mm colour of our sources and do not find evidence for `850um dropouts', which can be explained by the low-SNR of the observations. We also combine these observed colours with spectroscopic redshifts to derive the range of dust temperatures T, and dust emissivity indices $\beta$ for the sample, concluding that existing estimates T~30K and $\beta$~1.75 are consistent with these new data.

(abridged) We present deep Spitzer mid-infrared spectroscopy, along with 16, 24, 70, and 850um photometry, for 22 galaxies located in GOODS-N. The sample spans a redshift range of 0.6 < z < 2.6, 24um flux densities between ~0.2-1.2 mJy, and consists of SMGs, AGN, and optically faint (z_AB > 25) sources. We find that IR luminosities derived by fitting local SEDs with 24um photometry alone are well matched to those when additional mid-infrared spectroscopic and longer wavelength photometric data is used for galaxies having z < 1.4 and 24um-derived IR luminosities typically > 3x10^12 L_sun. However, for galaxies in the redshift range between 1.4 < z < 2.6, typically having 24um-derived IR luminosities > 3x10^12 L_sun, IR luminosities are overestimated by an average factor of ~5 when SED fitting with 24um photometry alone. This result arises partly due to the fact that high redshift galaxies exhibit PAH EQWs that are large compared to local galaxies of similar luminosities. Using improved estimates for the IR luminosities of these sources, we investigate whether their IR emission is found to be in excess relative to that expected based on extinction corrected UV SFRs, possibly suggesting the presence of an obscured AGN. Through a spectral decomposition of IRS spectroscopic data, we are able to isolate the fraction of IR luminosity arising from an AGN as opposed to star formation activity. This fraction is only ~30% of the total IR luminosity among the entire sample, on average. Of the sources identified as having mid-infrared excesses, half are accounted for by using proper bolometric corrections while half show the presence of obscured AGN. We do not find evidence for evolution in the FIR-radio correlation over this redshift range, although the SMGs have IR/radio ratios which are, on average, ~3 times lower than the nominal value.

We present the serendipitous discovery of z=4.05 molecular gas CO emission lines with the IRAM Plateau de Bure interferometer coincident with GN20 and GN20.2, two luminous submillimeter galaxies (SMGs) in the Great Observatories Origins Deep Survey North field (GOODS-N). These are among the most distant submillimeter-selected galaxies reliably identified through CO emission and also some of the most luminous known. In terms of CO to bolometric luminosities, stellar mass and star formation rates (SFRs), these newly discovered z>4 SMGs are similar to z~1.5-3 SMGs studied to date. These z~4 SMGs have much higher specific SFRs than typical B-band dropout Lyman break galaxies at the same redshift. The stellar mass-SFR correlation for normal galaxies does not seem to evolve much further, between z~2 and z~4. A significant z=4.05 spectroscopic redshift spike is observed in GOODS-N, and a strong spatial overdensity of B-band dropouts and IRAC selected z>3.5 galaxies appears to be centered on the GN20 and GN20.2 galaxies. This suggests a proto-cluster structure with total mass ~10^14 Msun. Using photometry at mid-IR, submm and radio wavelengths, we show that reliable photometric redshifts (Dz/(1+z)~0.1) can be derived for SMGs over 1<z<4. This new photometric redshift technique has been used to provide a first estimate of the space density of 3.5<z<6 hyper-luminous starburst galaxies, and to show that they contribute substantially to the SFR density at early epochs. Many of these high-redshift starbursts will be within reach of Herschel. We find that the radio to mid-IR flux density ratio can be used to select z>3.5 starbursts, regardless of their submm/mm emission [abridged].

We present the angular autocorrelation function of 2603 dust-obscured galaxies (DOGs) in the Bootes field of the NOAO Deep Wide-Field Survey. DOGs are red, obscured galaxies, defined as having R-[24] \ge 14 (F_24/F_R \ga 1000). Spectroscopy indicates that they are located at 1.5 \la z \la 2.5. We find strong clustering, with r_0 = 7.40^{+1.27}_{-0.84} Mpc/h for the full F_24 > 0.3 mJy sample. The clustering and space density of the DOGs are consistent with those of submillimeter galaxies, suggestive of a connection between these populations. We find evidence for luminosity-dependent clustering, with the correlation length increasing to r_0 = 12.97^{+4.26}_{-2.64} Mpc/h for brighter (F_24 > 0.6 mJy) DOGs. Bright DOGs also reside in richer environments than fainter ones, suggesting these subsamples may not be drawn from the same parent population. The clustering amplitudes imply average halo masses of log M = 12.2^{+0.3}_{-0.2} Msun for the full DOG sample, rising to log M = 13.0^{+0.4}_{-0.3} Msun for brighter DOGs. In a biased structure formation scenario, the full DOG sample will, on average, evolve into ~ 3 L* present-day galaxies, whereas the most luminous DOGs may evolve into brightest cluster galaxies.

We use deep far-IR, submm, radio and X-ray imaging and mid-IR spectroscopy to explore the nature of a sample of Spitzer-selected dust-obscured galaxies (DOGs) in GOODS-N. A sample of 79 galaxies satisfy the criteria R-[24]>14 (Vega) down to S24>100 microJy. Twelve of these galaxies have IRS spectra available which we use to measure redshifts and classify these objects as being dominated by star formation or active galactic nuclei (AGN) activity in the mid-IR. The IRS spectra and Spitzer photometric redshifts confirm that the DOGs lie in a tight redshift distribution around z~2. Based on mid-IR colors, 80% of DOGs are likely dominated by star formation; the stacked X-ray emission from this sub-sample of DOGs is also consistent with star formation. Since only a small number of DOGs are individually detected at far-IR and submm wavelengths, we use a stacking analysis to determine the average flux from these objects and plot a composite IR (8-1000 microns) spectral energy distribution (SED). The average luminosity of these star forming DOGs is LIR~1e12 Lsun. We compare the average star forming DOG to the average bright (S850>5 mJy) submillimeter galaxy (SMG); the S24>100 microJy DOGs are 3 times more numerous but 8 times less luminous in the IR. The far-IR SED shape of DOGs is similar to that of SMGs (average dust temperature of around 30 K) but DOGs have a higher mid-IR to far-IR flux ratio. The average star formation-dominated DOG has a star formation rate of 200 Msun/yr which, given their space density, amounts to a contribution of 0.01 Msun/yr/Mpc3 (or 5-10%) to the star formation rate density at z~2. We use the composite SED to predict the average flux of DOGs in future Herschel/PACS 100 micron and SCUBA-2 450 micron surveys and show that the majority of them will be detected.

We seek to improve estimates of the power spectrum covariance matrix from a limited number of simulations by employing a novel statistical technique known as shrinkage estimation. The shrinkage technique optimally combines an empirical estimate of the covariance with a model (the target) to minimize the total mean squared error compared to the true underlying covariance. We test this technique on N-body simulations and evaluate its performance by estimating cosmological parameters. Using a simple diagonal target, we show that the shrinkage estimator significantly outperforms both the empirical covariance and the target individually when using a small number of simulations. We find that reducing noise in the covariance estimate is essential for properly estimating the values of cosmological parameters as well as their confidence intervals. We extend our method to the jackknife covariance estimator and again find significant improvement, though simulations give better results. Even for thousands of simulations we still find evidence that our method improves estimation of the covariance matrix. Because our method is simple, requires negligible additional numerical effort, and produces superior results, we always advocate shrinkage estimation for the covariance of the power spectrum and other large-scale structure measurements when purely theoretical modeling of the covariance is insufficient.

Many models that seek to explain the origin of the unresolved X-ray background predict that Compton-thick Active Galactic Nuclei (AGNs) are ubiquitious at high redshift. However, few distant Compton-thick AGNs have been reliably identified to date. Here we present Spitzer-IRS spectroscopy and 3.6-70um photometry of a z=2.2 optically identified AGN (HDF-oMD49) that is formally undetected in the 2Ms Chandra Deep Field-North (CDF-N) survey. The Spitzer-IRS spectrum and spectral energy distribution of this object is AGN dominated, and a comparison of the energetics at X-ray wavelengths to those derived from mid-infrared (mid-IR) and optical spectroscopy shows that the AGN is intrinsically luminous (L_X~3x10^44 erg/s) but heavily absorbed by Compton-thick material (N_H>>10^24 cm^{-2}); i.e., this object is a Compton-thick quasar. Adopting the same approach that we applied to HDF-oMD49, we found a further six objects at z~2-2.5 in the literature that are also X-ray weak/undetected but have evidence for AGN activity from optical and/or mid-IR spectroscopy, and show that all of these sources are also Compton-thick quasars with L_X>10^44 erg/s. On the basis of the definition of Daddi etal. (2007), these Compton-thick quasars would be classified as mid-IR excess galaxies, and our study provides the first spectroscopic confirmation of Compton-thick AGN activity in a subsample of these z~2 mid-IR bright galaxies. Using the four objects that lie in the CDF-N field, we estimate the space-density of Compton-thick quasars [Phi~(0.7-2.5)x10^-5 Mpc^-3 for L_X>10^44 erg/s objects at z~2-2.5] and show that Compton-thick accretion is as ubiquitious as unobscured accretion in the distant Universe.

We have conducted a deep and uniform 1.1 mm survey of the GOODS-N field with AzTEC on the James Clerk Maxwell Telescope (JCMT). Here we present the first results from this survey including maps, the source catalogue, and 1.1 mm number-counts. The results presented here were obtained from a 245 sq-arcmin region with near uniform coverage to a depth of 0.96-1.16 mJy/beam. Our robust catalogue contains 28 source candidates detected with S/N >= 3.75, only 1-2 of which are expected to be spurious detections. Of these source candidates, 8 are also detected by SCUBA at 850 um in regions where there is good overlap between the two surveys. The major advantage of our survey over that with SCUBA is the uniformity of coverage. We calculate number counts using two different techniques: the first using a frequentist parameter estimation, and the second using a Bayesian method. The two sets of results are in good agreement. We find that the 1.1 mm differential number counts are well described in the 2-6 mJy range by the functional form dN/dS = N' (S'/S) exp(-S/S') with fitted parameters S' = 1.25 +/-0.38 mJy and dN/dS = 300 +/- 90 per mJy per sq-deg at 3 mJy.

We present a 1200-micron image of the Great Observatories Origin Deep Survey North (GOODS-N) field, obtained with the Max Planck Millimeter Bolometer array (MAMBO) on the IRAM 30-m telescope. The survey covers a contiguous area of 287 square arcmin to a near-uniform noise level of ~0.7mJy/beam. After Bayesian flux deboosting, a total of 30 sources are recovered (>=3.5sigma). An optimal combination of our 1200-micron data and an existing 850-micron image from the Submillimetre Common-User Bolometer Array (SCUBA) yielded 33 sources (>=4sigma). We combine our GOODS-N sample with those obtained in the Lockman Hole and ELAIS-N2 fields (Scott et al. 2002; Greve et al. 2004) in order to explore the degree of overlap between 1200-micron- and 850-micron-selected galaxies (hereafter SMGs), finding no significant difference between their 850-micron to 1200-micron flux density distributions. However, a noise-weighted stacking analysis yields a significant detection of the 1200-micron-blank SCUBA sources, whereas no significant 850-micron signal is found for the 850-micron-blank MAMBO sources. The hypothesis that the 850/1200-micron flux density distribution of SCUBA sources is also representative of the MAMBO population is rejected at the ~4sigma level, via Monte Carlo simulations. Therefore, although the populations overlap, galaxies selected at 850 and 1200micron are different, and there is compelling evidence for a significant 1200-micron-detected population which is not recovered at 850micron. These are submm drop-outs (SDOs), with S_850/S_1200 = 0.7-1.7, requiring very cold dust or unusual spectral energy distributions (T_d ~ 10K; beta ~ 1), unless SDOs reside beyond the redshift range observed for radio-identified SMGs, i.e. at z > 4.

We report the detection of CO(2-1) emission from the z=1.2 ultraluminous infrared galaxy (ULIRG) GOODS J123634.53+621241.3 (also known as the sub-millimeter galaxy GN26). These observations represent the first discovery of high-redshift CO emission using the new Combined Array for Research in Millimeter-Wave Astronomy (CARMA). Of all high-redshift (z>1) galaxies within the GOODS-North field, this source has the largest far-infrared (FIR) flux observed in the Spitzer 70um and 160um bands. The CO redshift confirms the optical identification of the source, and the bright CO(2-1) line suggests the presence of a large molecular gas reservoir of about 7x10^10 M(sun). The infrared-to-CO luminosity ratio of L(IR)/L'(CO) = 80+/-30 L(sun) (K Km/s pc^2)^-1 is slightly smaller than the average ratio found in local ULIRGs and high-redshift sub-millimeter galaxies. The short star-formation time scale of about 70 Myr is consistent with a starburst associated with the merger event and is much shorter than the time scales for spiral galaxies and estimates made for high-redshift galaxies selected on the basis of their B-z and z-K colors.

We present evidence for the existence of an IRAC excess in the spectral energy distribution (SED) of 5 galaxies at 0.6<z<0.9 and 1 galaxy at z=1.7. These 6 galaxies, located in the Great Observatories Origins Deep Survey field (GOODS-N), are star forming since they present strong 6.2, 7.7, and 11.3 um polycyclic aromatic hydrocarbon (PAH) lines in their Spitzer IRS mid-infrared spectra. We use a library of templates computed with PEGASE.2 to fit their multiwavelength photometry and derive their stellar continuum. Subtraction of the stellar continuum enables us to detect in 5 galaxies a significant excess in the IRAC band pass where the 3.3 um PAH is expected. We then assess if the physical origin of the IRAC excess is due to an obscured active galactic nucleus (AGN) or warm dust emission. For one galaxy evidence of an obscured AGN is found, while the remaining four do not exhibit any significant AGN activity. Possible contamination by warm dust continuum of unknown origin as found in the Galactic diffuse emission is discussed. The properties of such a continuum would have to be different from the local Universe to explain the measured IRAC excess, but we cannot definitively rule out this possibility until its origin is understood. Assuming that the IRAC excess is dominated by the 3.3 um PAH feature, we find good agreement with the observed 11.3 um PAH line flux arising from the same C-H bending and stretching modes, consistent with model expectations. Finally, the IRAC excess appears to be correlated with the star-formation rate in the galaxies. Hence it could provide a powerful diagnostic for measuring dusty star formation in z>3 galaxies once the mid-infrared spectroscopic capabilities of the James Webb Space Telescope become available.

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

We present 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.

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

We explore the epoch dependence of number density and star-formation rate for submillimetre galaxies (SMGs) found at 850 um. The study uses a sample of 38 SMG in the GOODS-N field, for which cross-waveband identifications have been obtained for 35/38 members together with redshift measurements or estimates. A maximum-likelihood analysis is employed, along with the `single-source-survey' technique. We find a diminution in both space density and star formation rate at z > 3, closely mimicking the redshift cut-offs found for QSOs selected in different wavebands. The diminution in redshift is particularly marked, at a significance level too small to measure. The data further suggest, at a significance level of about 0.001, that two separately-evolving populations may be present, with distinct luminosity functions. These results parallel the different evolutionary behaviours of LIRGs and ULIRGs, and represent another manifestation of `cosmic down-sizing', suggesting that differential evolution extends to the most extreme star-forming galaxies.

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 microns has been carried out to map the regions around several 850 micron-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared luminosities and hence star-formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350 micron flux densities for 24 850 micron sources, seven of which are detected at >2.5-sigma within a 10 arcsec search radius of the 850 micron positions. When results from the literature are included the total number of 350 micron flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-infrared (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich (Mdust~9x10^8 Msun), luminous (Lfir~2x10^12 Lsun), star-forming galaxies with intrinsic dust temperatures of ~35 K and star-formation rates of ~400 Msun/yr. We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350 micron source counts, N350(>25mJy)~200-500 deg^-2.

We present deep mid-infrared spectroscopy with the Spitzer Space Telescope of 13 submillimeter galaxies (SMGs) in the Great Observatories Origins Deep Survey North (GOODS-N) field. We find strong polycyclic aromatic hydrocarbon (PAH) emission in all of our targets which allows us to measure mid-IR spectroscopic redshifts and place constraints on the contribution from star formation and active galactic nuclei (AGN) activity to the mid-IR emission. In the high signal-to-noise ratio composite spectrum of SMGs, we find that the hot dust continuum from an AGN contributes at most 30% of the mid-IR luminosity. We find that the mid-IR properties of SMGs are distinct from those of 24 micron -selected ULIRGs at z~2; the former are predominantly dominated by star formation while the latter are a more heterogeneous sample with many showing significant AGN activity. We fit the IRS spectrum and the mid-IR to radio photometry of SMGs with template spectral energy distributions to determine the best estimate of the total infrared luminosity from star formation. While many SMGs contain an AGN as evinced by their X-ray properties, our multi-wavelength analysis shows that the total infrared luminosity, LIR, in SMGs, is dominated by star formation and not AGN activity. We find that high redshift SMGs lie on the relation between LIR and LPAH6.2 (or LPAH7.7 or LPAH11.3) that has been established for local starburst galaxies. This suggests that PAH luminosity can be used as a proxy for the star formation rate in SMGs. SMGs are consistent with being a short lived cool phase in a massive merger. Although there appears to be an AGN present in many SMGs, it does not appear to have become strong enough to heat the dust and dominate the mid- or far-infrared emission.

We introduce a method to constrain general cosmological models using Baryon Acoustic Oscillation (BAO) distance measurements from galaxy samples covering different redshift ranges, and apply this method to analyse samples drawn from the SDSS and 2dFGRS. BAO are detected in the clustering of the combined 2dFGRS and SDSS main galaxy samples, and measure the distance--redshift relation at z=0.2. BAO in the clustering of the SDSS luminous red galaxies measure the distance--redshift relation at z=0.35. The observed scale of the BAO calculated from these samples and from the combined sample are jointly analysed using estimates of the correlated errors, to constrain the form of the distance measure D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3). Here D_A is the angular diameter distance, and H(z) is the Hubble parameter. This gives r_s/D_V(0.2)=0.1980+/-0.0058 and r_s/D_V(0.35)=0.1094+/-0.0033 (1sigma errors), with correlation coefficient of 0.39, where r_s is the comoving sound horizon scale at recombination. Matching the BAO to have the same measured scale at all redshifts then gives D_V(0.35)/D_V(0.2)=1.812+/-0.060. The recovered ratio is roughly consistent with that predicted by the higher redshift SNLS supernovae data for Lambda cosmologies, but does require slightly stronger cosmological acceleration at low redshift. If we force the cosmological model to be flat with constant w, then we find Om_m=0.249+/-0.018 and w=-1.004+/-0.089 after combining with the SNLS data, and including the WMAP measurement of the apparent acoustic horizon angle in the CMB.

We have studied the submillimetre (submm) properties of the following classes of near-infrared (NIR)-selected massive galaxies at high redshifts: BzK-selected star-forming galaxies (BzKs); distant red galaxies (DRGs); and extremely red objects (EROs). We used the SCUBA HAlf Degree Extragalactic Survey (SHADES), the largest uniform submm survey to date. Partial overlap of SIRIUS/NIR images and SHADES in SXDF has allowed us to identify 4 submm-bright NIR-selected galaxies, which are detected in the mid-infrared, 24 micron, and the radio, 1.4 GHz. We find that all of our submm-bright NIR-selected galaxies satisfy the BzK selection criteria, except for one galaxy whose B-z and z-K colours are however close to the BzK colour boundary. Two of the submm-bright NIR-selected galaxies satisfy all of the selection criteria we considered, i.e. they belong to the BzK-DRG-ERO overlapping population, or `extremely red' BzKs. Although these extremely red BzKs are rare (0.25 arcmin^{-2}), up to 20 % of this population could be submm galaxies. This fraction is significantly higher than that found for other galaxy populations studied here. Via a stacking analysis, we have detected the 850 micron flux of submm-faint BzKs and EROs in our SCUBA maps. While the contribution of z~2 BzKs to the submm background is about 10--15 % and similar to that from EROs typically at z~1, BzKs have a higher fraction (~30 %) of submm flux in resolved sources compared with EROs and submm sources as a whole. From the SED fitting analysis for both submm-bright and submm-faint BzKs, we found no clear signature that submm-bright BzKs are experiencing a specifically luminous evolutionary phase, compared with submm-faint BzKs. An alternative explanation might be that submm-bright BzKs are more massive than submm-faint ones.

Approximately 20-30% of 1.4<z<2.5 galaxies with K<22 (Vega) detected with Spitzer MIPS at 24um show excess mid-IR emission relative to that expected based on the rates of star formation measured from other multiwavelength data.These galaxies also display some near-IR excess in Spitzer IRAC data, with a spectral energy distribution peaking longward of 1.6um in the rest frame, indicating the presence of warm-dust emission usually absent in star forming galaxies. Stacking Chandra data for the mid-IR excess galaxies yields a significant hard X-ray detection at rest-frame energies >6.2 keV. The stacked X-ray spectrum rises steeply at >10 keV, suggesting that these sources host Compton-thick Active Galactic Nuclei (AGNs) with column densities N_H~10^{24} cm^-2 and an average, unobscured X-ray luminosity L_{2-8 keV}~(1-4)x10^43 erg/s. Their sky density(~3200 deg^-2) and space density (~2.6x10^-4 Mpc^-3) are twice those of X-ray detected AGNs at z~2, and much larger than those of previously-known Compton thick sources at similar redshifts. The mid-IR excess galaxies are part of the long sought-after population of distant heavily obscured AGNs predicted by synthesis models of the X-ray background. The fraction of mid-IR excess objects increases with galaxy mass, reaching ~50-60% for M~10^11 M_sun, an effect likely connected with downsizing in galaxy formation. The ratio of theinferred black hole growth rate from these Compton-thick sources to the global star formation rate at z=2 is similar to the mass ratio of black holes to stars inlocal spheroids, implying concurrent growth of both within the precursors oftoday's massive galaxies.

Examining a sample of massive galaxies at 1.4<z<2.5 with K_{Vega}<22 from the Great Observatories Origins Deep Survey, we compare photometry from Spitzer at mid- and far-IR, to submillimeter, radio and rest-frame ultraviolet wavelengths, to test the agreement between different tracers of star formation rates (SFRs) and to explore the implications for galaxy assembly. For z~2 galaxies with moderate luminosities(L_{8um}<10^{11}L_sun), we find that the SFR can be estimated consistently from the multiwavelength data based on local luminosity correlations. However,20--30% of massive galaxies, and nearly all those with L_{8um}>10^{11}L_sun, show a mid-IR excess which is likely due to the presence of obscured active nuclei, as shown in a companion paper. There is a tight and roughly linear correlation between stellar mass and SFR for 24um-detected galaxies. For a given mass, the SFR at z=2 was larger by a factor of ~4 and ~30 relative to that in star forming galaxies at z=1 and z=0, respectively. Typical ultraluminous infrared galaxies (ULIRGs) at z=2 are relatively 'transparent' to ultraviolet light, and their activity is long lived (~400 Myr), unlike that in local ULIRGs and high redshift submillimeter-selected galaxies. ULIRGs are the common mode of star formation in massive galaxies at z=2, and the high duty cycle suggests that major mergers are not the dominant trigger for this activity.Current galaxy formation models underpredict the normalization of the mass-SFR correlation by about a factor of 4, and the space density of ULIRGs by an orderof magnitude, but give better agreement for z>1.4 quiescent galaxies.

We present the redshift distribution of the SHADES galaxy population based on the rest-frame radio-mm-FIR colours of 120 robustly detected 850um sources in the Lockman Hole East (LH) and Subaru XMM-Newton Deep Field (SXDF). The redshift distribution derived from the full SED information is shown to be narrower than that determined from the radio-submm spectral index, as more photometric bands contribute to a higher redshift accuracy. The redshift distribution of sources derived from at least two photometric bands peaks at z ~ 2.4 and has a near-Gaussian distribution, with 50 per cent (interquartile range) of sources at z=1.8-3.1. We find a statistically-significant difference between the measured redshift distributions in the two fields; the SXDF peaking at a slightly lower redshift (median z ~ 2.2) than the LH (median z ~ 2.7), which we attribute to the noise-properties of the radio observations. We demonstrate however that there could also be field-to-field variations that are consistent with the measured differences in the redshift distributions, and hence, that the incomplete area observed by SHADES with SCUBA, despite being the largest sub-mm survey to date, may still be too small to fully characterize the bright sub-mm galaxy population. Finally we present a brief comparison with the predicted, or assumed, redshift distributions of sub-mm galaxy formation and evolution models, and we derive the contribution of these SHADES sources and the general sub-mm galaxy population to the star formation-rate density at different epochs.

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

The sub-millimetre wavelength regime is perhaps the most poorly explored over large areas of the sky, despite the considerable effort that has been expended in making deep maps over small regions. As a consequence the properties of the sub-millimetre sky as a whole, and of rare bright objects in particular, remains largely unknown. Here we describe a forthcoming survey (the SCUBA-2 ``All-Sky'' Survey, or SASSy) designed to address this issue by making a large-area map of approximately one-fifth of the sky visible from the JCMT (4800 square degrees) down to a 1 sigma noise level of 30 mJy/beam. This map forms the pilot for a much larger survey, which will potentially map the remaining sky visible from the JCMT, with the region also visible to ALMA as a priority. SASSy has been awarded 500 hours for the 4800 square degree pilot phase and will commence after the commissioning of SCUBA-2, expected in early 2008.

The study of the dominant population of high redshift IR-luminous galaxies (10^11 - 10^12 Lsun at 1<z<3), requires observation of sources at the ~0.1 mJy level in the mid-IR. We present the deepest spectra taken to date with the Infrared Spectrograph (IRS) on the Spitzer Space Telescope. We targeted two faint (f24~0.15 mJy) sources in the Southern GOODS field at z=1.09 and z=2.69. Spectra of the lower redshift target were taken in the observed-frame 8--21 micron range, while the spectrum of the higher redshift target covered 21--37 microns. We also present the spectra of two secondary sources within the slit. We detect strong PAH emission in all four targets, and compare the spectra to those of local galaxies observed by the IRS. The z=1.09 source appears to be a typical, star-formation dominated IR-luminous galaxy, while the z=2.69 source is a composite source with strong star formation and a prominent AGN. The IRAC colors of this source show no evidence of rest-frame near-infrared stellar photospheric emission. We demonstrate that an AGN which contributes only a small (~10%) fraction of the bolometric luminosity can produce enough hot dust emission to overwhelm the near-IR photospheric emission from stars. Such sources would be excluded from photometric surveys which rely on the near-IR bump to identify starbursts, leading to an underestimate of the star formation rate density.

We present 70 micron properties of submillimeter galaxies (SMGs) in the Great Observatories Origins Deep Survey (GOODS) North field. Out of thirty submillimeter galaxies (S_850 > 2 mJy) in the central GOODS-N region, we find two with secure 70 micron detections. These are the first 70 micron detections of SMGs. One of the matched SMGs is at z ~ 0.5 and has S_70/S_850 and S_70/S_24 ratios consistent with a cool galaxy. The second SMG (z = 1.2) has infrared-submm colors which indicate it is more actively forming stars. We examine the average 70 micron properties of the SMGs by performing a stacking analysis, which also allows us to estimate that S_850 > 2 mJy SMGs contribute 9 +- 3% of the 70 micron background light. The S_850/S_70 colors of the SMG population as a whole is best fit by cool galaxies, and because of the redshifting effects these constraints are mainly on the lower z sub-sample. We fit Spectral Energy Distributions (SEDs) to the far-infrared data points of the two detected SMGs and the average low redshift SMG (z_{median}= 1.4). We find that the average low-z SMG has a cooler dust temperature than local ultraluminous infrared galaxies (ULIRGs) of similar luminosity and an SED which is best fit by scaled up versions of normal spiral galaxies. The average low-z SMG is found to have a typical dust temperature T = 21 -- 33 K and infrared luminosity L_{8-1000 micron} = 8.0 \times 10^11 L_sun. We estimate the AGN contribution to the total infrared luminosity of low-z SMGs is less than 23%.

Submillimetre (sub-mm) galaxies have very high infrared (IR) luminosities and are thousands of times more numerous at z~2 than local ultra-luminous IR galaxies. They therefore represent a key phase in galaxy evolution which can be missed in optical surveys. Determining their contribution to the global star formation rate requires dissecting their IR emission into contributions from starbursts (SB) and active galactic nuclei (AGN). There are several examples of AGN systems which masquerade as SBs in either the IR or X-ray, and SBs can often look like AGN in some wavebands. A combination of SB and AGN emission is not unreasonable, given models of merger-driven evolution. To assess in detail what powers the intense IR luminosity of sub-mm galaxies it is important to obtain a complete multi-wavelength picture. Mid-IR spectroscopy is a particularly good probe of where the intense IR luminosity is coming from. We present the first results from a program to obtain Spitzer IRS spectroscopy of a sample of high redshift galaxies in the GOODS-N field, a large fraction of which are sub-mm galaxies. This field is already home to the deepest X-ray, optical, IR and radio data. We piece together the sub-mm data with the Spitzer photometry and IRS spectra to provide a well sampled IR spectral energy distribution (SED) of sub-mm galaxies and determine the contribution to the bolometric luminosity from the AGN and SB components.

We measure the cosmological matter density by observing the positions of baryon acoustic oscillations in the clustering of galaxies in the Sloan Digital Sky Survey (SDSS). We jointly analyse the main galaxies and LRGs in the SDSS DR5 sample, using over half a million galaxies in total. The oscillations are detected with 99.74% confidence (3.0sigma assuming Gaussianity) compared to a smooth power spectrum. When combined with the observed scale of the peaks within the CMB, we find a best-fit value of Omega_m=0.256+0.029-0.024 (68% confidence interval), for a flat Lambda cosmology when marginalising over the Hubble parameter and the baryon density. This value of the matter density is derived from the locations of the baryon oscillations in the galaxy power spectrum and in the CMB, and does not include any information from the overall shape of the power spectra. This is an extremely clean cosmological measurement as the physics of the baryon acoustic oscillation production is well understood, and the positions of the oscillations are expected to be independent of systematics such as galaxy bias.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

We present a Fourier analysis of the clustering of galaxies in the combined Main galaxy and Luminous Red Galaxy (LRG) Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) sample. The aim of our analysis is to consider how well we can measure the cosmological matter density using the signature of the horizon at matter-radiation equality embedded in the large-scale power spectrum. The new data constrains the power spectrum on scales 100--600h^-1Mpc with significantly higher precision than previous analyses of just the SDSS Main galaxies, due to our larger sample and the inclusion of the LRGs. This improvement means that we can now reveal a discrepancy between the shape of the measured power and linear CDM models on scales 0.01<k<0.15hMpc^-1, with linear model fits favouring a lower matter density (Omega_m=0.22+/-0.04) on scales 0.01<k<0.06hMpc^-1 and a higher matter density (Omega_m=0.32+/-0.01) when smaller scales are included, assuming a flat LCDM model with h=0.73 and n_s=0.96. This discrepancy could be explained by scale-dependent bias and, by analysing subsamples of galaxies, we find that the ratio of small-scale to large-scale power increases with galaxy luminosity, so all of the SDSS galaxies cannot trace the same power spectrum shape over 0.01<k<0.2hMpc^-1. However, the data are insufficient to clearly show a luminosity-dependent change in the largest scale at which a significant increase in clustering is observed, although they do not rule out such an effect. Significant scale-dependent galaxy bias on large-scales, which changes with the r-band luminosity of the galaxies, could potentially explain differences in our Omega_m estimates and differences previously observed between 2dFGRS and SDSS power spectra and the resulting parameter constraints.

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

We present SEDs, Spitzer colours, and IR luminosities for 850 micron selected galaxies in the GOODS-N field. Using the deep Spitzer Legacy images and new data and reductions of the VLA-HDF radio data, we find statistically secure counterparts for 60 per cent (21/35) of our submm sample, and identify tentative counterparts for another 12 objects. This is the largest sample of submm galaxies with statistically secure counterparts detected in the radio and with Spitzer. We find that in most cases the 850 micron emission is dominated by a single 24 micron source. A composite rest-frame SED shows that the submm sources peak at longer wavelengths than those of local ULIRGs of the same luminosity and therefore appear to be cooler. The SEDs of submm galaxies are also different from those of their high redshift neighbours, the near-IR selected BzK galaxies, whose mid-IR to radio SEDs are more like those of local ULIRGs. Using 24 micron. 850 micron and 1.4 GHz observations, we fit templates that span the mid-IR through radio to derived the integrated IR luminosity of the submm galaxies. By themselves, 24 micron and radio fluxes are able to predict LIR reasonably well because they are relatively insensitive to temperature. However, the submm flux by itself consistently overpredicts LIR when using spectral templates which obey the local ULIRG temperature-luminosity relation. The shorter Spitzer wavelengths sample the stellar bump at the redshifts of the submm sources, and we find that the Spitzer photometry alone provides a model independent estimate of the redshift. Using X-ray and mid-IR data, only 5 per cent of our secure counterparts show strong evidence for an active galactic nucleus dominating the IR luminosity.

How does the submm galaxy population detected by SCUBA fit into galaxy evolution? How do these rare starbursting systems, which contribute significantly to high redshift star-formation, relate to other galaxy populations? Deep radio observations have been most useful for studying these systems, but still leave a significant fraction of the population unidentified. Now with the deep IRAC and MIPS images from the GOODS Spitzer Legacy program and a re-analysis of the deep radio data, we are able to identify counterparts for a large fraction of SCUBA galaxies in GOODS-N. All of these counterparts are detected by Spitzer. Given the vast multi-wavelength data in this field, we can study the spectral energy distributions (SEDs) of these systems and determine what is fueling their intense infrared luminosities. A rest-frame composite optical-to-radio SED for all spectroscopically identified submm sources shows that the average SCUBA galaxy is consistent with models of ultraluminous starburst galaxies, although cooler than those observed locally. Because of this, the submm flux alone consistently overestimates IR when using spectral templates which obey the local ULIRG temperature-luminosity relation. The wide range of 24/850 micron flux ratios as a function of redshift indicates the presence of strong mid-IR features, to be confirmed with deep IRS spectroscopy. The IRAC colours of the submm systems provide useful redshift constraints, since, at these redshifts, IRAC samples the stellar bump. The Spitzer photometry of this large sample of submm galaxies has allowed us to put constraints on many of the outstanding issues in submm astronomy.

We report on a preliminary study concerning the origin of radio emission within radio galaxies at L(1.4GHz)>1E24 W/Hz in the GOODS-N field. In the local universe, Condon et al. (2002) and Yun et al. (2001) have shown that in galaxies with radio luminosities greater than 1E23 W/Hz the majority of the radio emission originates from a `monster' i.e., an AGN. Using the Chandra 2Msec X-ray image centered on the GOODS-N field and a reprocessed VLA HDF A-array data plus newly acquired VLA B-array data (rms=5.3microJy), we find that radio galaxies (with spectroscopic redshifts; all have z>1) with L(1.4GHz)>1E24 W/Hz typically have an X-ray detection rate of 72% (60% emit hard X-rays suggesting an AGN origin for the radio emission) in contrast to 25% for radio galaxies with L < 1E23 W/Hz. The ACS images of these L(1.4 GHz) > 1E24 W/Hz galaxies typically show compact rather than extended galaxy morphology which is generally found for the less luminous radio emitting galaxies but a few appear to be ongoing galaxy mergers. We also present SED fitting for these luminous radio galaxies including Spitzer IRAC & MIPS 24um photometry and 60% show distinct power-law SED indicative of an AGN. Initial results tell us that the X-ray emitting radio galaxy population are generally not submm sources but the few (~10%) that are SCUBA sources appear to be the small AGN population found by Pope et al. and others.

We present measurements of the normalised redshift-space three-point correlation function (Q_z) of galaxies from the Sloan Digital Sky Survey (SDSS) main galaxy sample. We have applied our "npt" algorithm to both a volume-limited (36738 galaxies) and magnitude-limited sample (134741 galaxies) of SDSS galaxies, and find consistent results between the two samples, thus confirming the weak luminosity dependence of Q_z recently seen by other authors. We compare our results to other Q_z measurements in the literature and find it to be consistent within the full jack-knife error estimates. However, we find these errors are significantly increased by the presence of the ``Sloan Great Wall'' (at z ~ 0.08) within these two SDSS datasets, which changes the 3-point correlation function (3PCF) by 70% on large scales (s>=10h^-1 Mpc). If we exclude this supercluster, our observed Q_z is in better agreement with that obtained from the 2dFGRS by other authors, thus demonstrating the sensitivity of these higher-order correlation functions to large-scale structures in the Universe. This analysis highlights that the SDSS datasets used here are not ``fair samples'' of the Universe for the estimation of higher-order clustering statistics and larger volumes are required. We study the shape-dependence of Q_z(s,q,theta) as one expects this measurement to depend on scale if the large scale structure in the Universe has grown via gravitational instability from Gaussian initial conditions. On small scales (s <= 6h^-1 Mpc), we see some evidence for shape-dependence in Q_z, but at present our measurements are consistent with a constant within the errors (Q_z ~ 0.75 +/- 0.05). On scales >10h^-1 Mpc, we see considerable shape-dependence in Q_z.

We determine the number counts and z=0-5 luminosity function for a well-defined, homogeneous sample of quasars from the Sloan Digital Sky Survey (SDSS). We conservatively define the most uniform statistical sample possible, consisting of 15,343 quasars within an effective area of 1622 deg^2 that was derived from a parent sample of 46,420 spectroscopically confirmed broad-line quasars in the 5282 deg^2 of imaging data from SDSS Data Release Three. The sample extends from i=15 to i=19.1 at z<3 and to i=20.2 for z>3. The number counts and luminosity function agree well with the results of the 2dF QSO Survey, but the SDSS data probe to much higher redshifts than does the 2dF sample. The number density of luminous quasars peaks between redshifts 2 and 3, although uncertainties in the selection function in this range do not allow us to determine the peak redshift more precisely. Our best fit model has a flatter bright end slope at high redshift than at low redshift. For z<2.4 the data are best fit by a redshift-independent slope of beta = -3.1 (Phi(L) propto L^beta). Above z=2.4 the slope flattens with redshift to beta=-2.37 at z=5. This slope change, which is significant at a >5-sigma level, must be accounted for in models of the evolution of accretion onto supermassive black holes.

Bright sub-mm galaxies are expected to arise in massive highly-biased haloes, and hence exhibit strong clustering. We argue that a valuable tool for measuring these clustering properties is the cross-correlation of sub-mm galaxies with faint optically-selected sources. We analyze populations of SCUBA-detected and optical galaxies in the GOODS-N survey area. Using optical/IR photometric-redshift information, we search for correlations induced by two separate effects: (1) cosmic magnification of background sub-mm sources by foreground dark matter haloes traced by optical galaxies at lower redshifts; and (2) galaxy clustering due to sub-mm and optical sources tracing the same population of haloes where their redshift distributions overlap. Regarding cosmic magnification, we find no detectable correlation. Our null result is consistent with a theoretical model for the cosmic magnification, and we show that a dramatic increase in the number of sub-mm sources will be required to measure the effect reliably. Regarding clustering, we find evidence at the 3.5-sigma level for a cross-correlation between sub-mm and optical galaxies analyzed in identical photometric redshift slices. The data hint that the sub-mm sources have an enhanced bias parameter compared to the optically-selected population (with a significance of 2-sigma). The next generation of deep sub-mm surveys can potentially perform an accurate measurement of each of these cross-correlations, adding a new set of diagnostics for understanding the development of massive structure in the Universe.

We present high resolution 890 micron images of two 20 mJy submillimeter galaxies, SMMJ123711+622212 and MIPS~J142824.0+352619, obtained using the Submillimeter Array (SMA). Using submillimeter interferometric observations with an angular resolution of 2.5", the coordinates of these high redshift sources are determined with an accuracy of 0.2". The new SMA data on SMMJ123711+622212 reveal an unresolved submm source offset to the east by 0.8" from an optical galaxy found in deep HST images, suggesting either a large galaxy with a dusty central region, or an interacting galaxy system. The SMA image of hyper-luminous (LFIR = 3.2 x 10^{13} Lsun) source MIPS~J142824.0+352619 provides a firm upper limit to the source size of < 1.2. This constraint provides evidence that the foreground lens is only weakly affecting the observed high FIR luminosity.

We present Spitzer/MIPS 24 micron observations of 50 supernova host galaxies at 0.1<z<1.7 in the Great Observatories Origins Deep Survey (GOODS) fields. We also discuss the detection of SN host galaxies in SCUBA/850 micron observations of GOODS-N and Spitzer/Infrared Spectrograph (IRS) 16 micron observations of GOODS-S. About 60% of the host galaxies of both Type Ia and core-collapse supernovae are detected at 24 microns, a detection rate which is a factor of 1.5 higher than the field galaxy population. Among the 24 micron detected hosts, 80% have far-infrared luminosities that are comparable to or greater than the optical luminosity indicating the presence of substantial amounts of dust in the hosts. The median bolometric luminosity of the Type Ia SN hosts is \~10^10.5 L_sun, very similar to that of core-collapse SN hosts. Using the high resolution Hubble/ACS data, we have studied the variation of rest-frame optical/ultraviolet colors within the 24 micron detected galaxies at z<1 to understand the origin of the dust emission. The 24 micron detected galaxies have average colors which are redder by ~0.1 mag than the 24 micron undetected hosts while the latter show greater scatter in internal colors. This suggests that a smooth distribution of dust is responsible for the observed mid- and far-infrared emission. 70% of supernovae that have been detected in the GOODS fields are located within the half-light radius of the hosts where dust obscuration effects are significant. Although the dust emission that we detect cannot be translated into a line of sight A_V, we suggest that the factor of 2-3 larger scatter in the peak B-V colors that is seen in the high-z Type Ia supernova sample relative to the low-z supernovae might be partially due to the dust that we detect in the hosts.

We investigate the multi-wavelength emission of BzK selected star forming galaxies at z~2 in the Great Observatories Origins Deep Survey (GOODS) North region. Most (82%) of the sources are individually detected at 24um in the Spitzer MIPS imaging, and one fourth (26%) in the VLA radio data. Significant detections of the individually undetected objects are obtained through stacking in the radio, submm and X-ray domains. The typical star forming galaxy with stellar mass ~10^{11}Mo at z=2 is an Ultra-luminous Infrared Galaxy (ULIRG), with L_IR ~ 1-2x10^{12}Lo and star formation rate SFR 200-300Mo/yr, implying a comoving density of ULIRGs at z=2 at least 3 orders of magnitude above the local one. SFRs derived from the reddening corrected UV luminosities agree well, on average, with the longer wavelength estimates. The high 24um detection rate suggests a relatively large duty cycle for the BzK star forming phase, consistently with the available independent measurements of the space density of passively evolving galaxies at z>1.4. If the IMF at z=2 is similar to the local one, and in particular is not a top-heavy IMF, this suggests that a substantial fraction of the high mass tail (>10^{11}Mo) of the galaxy stellar mass function was completed by z~1.4.

We present a new sub-mm Super-map in the HDF-North region (GOODS-North field), containing 40 statistically robust sources at 850 microns. This map contains additional data, and several new sources, including one of the brightest blank-sky extragalactic sub-mm sources ever detected. We have used the ACS HST images and ground-based near-IR observations from GOODS to develop a systematic approach for counterpart identification. 72 per cent of our sources with optical coverage have a unique optical counterpart using our new techniques for counterpart identification, and an additional 18 per cent have more than one possibility that meet our criteria in the ACS images. We have found a much higher ERO rate than other sub-mm surveys, due to the increased depth in the optical images. The median photometric redshift (and quartile range), from optical and near-infrared data, is 1.7 (1.3-2.5) for the radio-detected sub-mm sources, and rises to 2.3 (1.3-2.7) for the radio-undetected sub-sample. We find interesting correlations between the 850 micron flux and both the i magnitude and the photometric redshift, from which there appears to be an absence of high redshift faint counterparts to the lower flux density SCUBA sources. While the quantitative morphologies span a range of values, in general the sub-mm galaxies show larger sizes and a higher degree of asymmetry than other galaxy populations at the same redshifts.

We present the large-scale correlation function measured from a spectroscopic sample of 46,748 luminous red galaxies from the Sloan Digital Sky Survey. The survey region covers 0.72 h^{-3} Gpc^3 over 3816 square degrees and 0.16<z<0.47, making it the best sample yet for the study of large-scale structure. We find a well-detected peak in the correlation function at 100h^{-1} Mpc separation that is an excellent match to the predicted shape and location of the imprint of the recombination-epoch acoustic oscillations on the low-redshift clustering of matter. This detection demonstrates the linear growth of structure by gravitational instability between z=1000 and the present and confirms a firm prediction of the standard cosmological theory. The acoustic peak provides a standard ruler by which we can measure the ratio of the distances to z=0.35 and z=1089 to 4% fractional accuracy and the absolute distance to z=0.35 to 5% accuracy. From the overall shape of the correlation function, we measure the matter density Omega_mh^2 to 8% and find agreement with the value from cosmic microwave background (CMB) anisotropies. Independent of the constraints provided by the CMB acoustic scale, we find Omega_m = 0.273 +- 0.025 + 0.123 (1+w_0) + 0.137 Omega_K. Including the CMB acoustic scale, we find that the spatial curvature is Omega_K=-0.010+-0.009 if the dark energy is a cosmological constant. More generally, our results provide a measurement of cosmological distance, and hence an argument for dark energy, based on a geometric method with the same simple physics as the microwave background anisotropies. The standard cosmological model convincingly passes these new and robust tests of its fundamental properties.

We investigate the photometric redshift distribution, colours and morphologies of a large sample of submillimetre sources in the GOODS-North region. With the depth achieved by the deep ACS images, optical counterparts have been found for all the radio-detected sub-mm sources, and we have used several techniques to identify counterparts to the radio-undetected sources. The sub-mm counterparts have a median photometric redshift, from optical and near-IR photometry, of 1.9. Our sample contains a much higher fraction of EROs than other sub-mm surveys, and larger angular sizes than the average galaxy at all redshifts. The sub-mm galaxy counterparts are consistent with massive star-forming galaxies at high redshift.

We present radio, optical and X-ray detected counterparts to the sub-mm sources found using SCUBA in the Hubble Deep Field North region (GOODS-N). A new counterpart identification statistic is developed to identify properties of galaxies detected at other wavelengths that can be used to aid counterpart identification. We discriminate between criteria that can be used to pre-select sub-mm bright objects, and those that identify the counterpart to a known sub-mm object. Optically faint galaxies detected in the deepest 1.4 GHz radio continuum maps are the only effective way of pre-selecting SCUBA galaxies, and radio sources are the best way to identify counterparts to known sub-mm detections. Looking at radio spectral indices, only the steeper sources (indicative of star formation) are detected in the sub-mm. Although we find several X-ray identifications, we show that deep Chandra images do not contribute to counterpart identifications, since in all cases they are already detected in the more easily obtained VLA radio maps. We also find find no evidence for clustering between Chandra and SCUBA sources in this field. For a known SCUBA position, the reddest source tends to be the correct association, although we can find no cut on colour, magnitude, or clustering property that efficiently pre-selects for SCUBA sources. 15 micron ISO sources are statistically detected by SCUBA, but the limiting mid-IR flux is not low enough to provide useful constraints. We present postage stamp strips for each SCUBA detection in separate bands from X-ray to radio, providing direct visual evidence that approximately half of the sub-mm sources in this field remain unidentified, despite an abundance of deep multi-wavelength data.

We present the two--point correlation function (2PCF) of narrow-line active galactic nuclei (AGN) selected within the First Data Release of the Sloan Digital Sky Survey. Using a sample of 13605 AGN in the redshift range 0.055 < z < 0.2, we find that the AGN auto--correlation function is consistent with the observed galaxy auto--correlation function on scales 0.2h^{-1}Mpc to >100h^{-1}Mpc. The AGN hosts trace an intermediate population of galaxies and are not detected in either the bluest (youngest) disk--dominated galaxies or many of the reddest (oldest) galaxies. We show that the AGN 2PCF is dependent on the luminosity of the narrow [OIII] emission line (L_{[OIII]}), with low L_{[OIII]} AGN having a higher clustering amplitude than high L_{[OIII]} AGN. This is consistent with lower activity AGN residing in more massive galaxies than higher activity AGN, and L_{[OIII]} providing a good indicator of the fueling rate. Using a model relating halo mass to black hole mass in cosmological simulations, we show that AGN hosted by ~ 10^{12} M_{odot} dark matter halos have a 2PCF that matches that of the observed sample. This mass scale implies a mean black hole mass for the sample of M_{BH} ~ 10^8 M_{odot}.

We present estimates of cosmological parameters from the application of the Karhunen-Loeve transform to the analysis of the 3D power spectrum of density fluctuations using Sloan Digital Sky Survey galaxy redshifts. We use Omega_m*h and f_b = Omega_b/Omega_m to describe the shape of the power spectrum, sigma8 for the (linearly extrapolated) normalization, and beta to parametrize linear theory redshift space distortions. On scales k < 0.16 h/Mpc, our maximum likelihood values are Omega_m*h = 0.264 +/-0.043, f_b = 0.286 +/- 0.065, sigma8 = 0.966 +/- 0.048, and beta = 0.45 +/- 0.12. When we take a prior on Omega_b from WMAP, we find Omega_m*h = 0.207 +/- 0.030, which is in excellent agreement with WMAP and 2dF. This indicates that we have reasonably measured the gross shape of the power spectrum but we have difficulty breaking the degeneracy between Omega_m*h and f_b because the baryon oscillations are not resolved in the current spectroscopic survey window function.

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

Eigenmode analysis is one of the most promising methods of analyzing large data sets in ongoing and near-future galaxy surveys. In such analyses, a fast evaluation of the correlation matrix in arbitrary cosmological models is crucial. The observational effects, including peculiar velocity distortions in redshift space, light-cone effects, selection effects, and effects of the complex shape of the survey geometry, should be taken into account in the analysis. In the framework of the linear theory of gravitational instability, we provide the methodology to quickly compute the correlation matrix. Our methods are not restricted to shallow redshift surveys, arbitrarily deep samples can be dealt with as well. Therefore, our methods are useful in constraining the geometry of the universe and the dark energy component, as well as the power spectrum of galaxies, since ongoing and near-future galaxy surveys probe the universe at intermediate to deep redshifts, z ~ 0.2--5. In addition to the detailed methods to compute the correlation matrix in 3-dimensional redshift surveys, methods to calculate the matrix in 2-dimensional projected samples are also provided. Prospects of applying our methods to likelihood estimation of the cosmological parameters are discussed.

We measure the large-scale real-space power spectrum P(k) using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 square degrees with mean redshift z~0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.02 h/Mpc < k < 0.3h/Mpc. We pay particular attention to modeling, quantifying and correcting for potential systematic errors, nonlinear redshift distortions and the artificial red-tilt caused by luminosity-dependent bias. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k<0.1h/Mpc, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the WMAP satellite. As a simple characterization of the data, our measurements are well fit by a flat scale-invariant adiabatic cosmological model with h Omega_m =0.201+/- 0.017 and L* galaxy sigma_8=0.89 +/- 0.02 when fixing the baryon fraction Omega_b/Omega_m=0.17 and the Hubble parameter h=0.72; cosmological interpretation is given in a companion paper.

We present measurements of parameters of the 3-dimensional power spectrum of galaxy clustering from 222 square degrees of early imaging data in the Sloan Digital Sky Survey. The projected galaxy distribution on the sky is expanded over a set of Karhunen-Loeve eigenfunctions, which optimize the signal-to-noise ratio in our analysis. A maximum likelihood analysis is used to estimate parameters that set the shape and amplitude of the 3-dimensional power spectrum. Our best estimates are Gamma=0.188 +/- 0.04 and sigma_8L = 0.915 +/- 0.06 (statistical errors only), for a flat Universe with a cosmological constant. We demonstrate that our measurements contain signal from scales at or beyond the peak of the 3D power spectrum. We discuss how the results scale with systematic uncertainties, like the radial selection function. We find that the central values satisfy the analytically estimated scaling relation. We have also explored the effects of evolutionary corrections, various truncations of the KL basis, seeing, sample size and limiting magnitude. We find that the impact of most of these uncertainties stay within the 2-sigma uncertainties of our fiducial result.

On 26 May 1999, one of the Sloan Digital Sky Survey (SDSS) fiber-fed spectrographs saw astronomical first light. This was followed by the first spectroscopic commissioning run during the dark period of June 1999. We present here the first hour of extra-galactic spectroscopy taken during these early commissioning stages: an observation of the Coma cluster of galaxies. Our data samples the Southern part of this cluster, out to a radius of 1.5degrees and thus fully covers the NGC 4839 group. We outline in this paper the main characteristics of the SDSS spectroscopic systems and provide redshifts and spectral classifications for 196 Coma galaxies, of which 45 redshifts are new. For the 151 galaxies in common with the literature, we find excellent agreement between our redshift determinations and the published values. As part of our analysis, we have investigated four different spectral classification algorithms: spectral line strengths, a principal component decomposition, a wavelet analysis and the fitting of spectral synthesis models to the data. We find that a significant fraction (25%) of our observed Coma galaxies show signs of recent star-formation activity and that the velocity dispersion of these active galaxies (emission-line and post-starburst galaxies) is 30% larger than the absorption-line galaxies. We also find no active galaxies within the central (projected) 200 h-1 Kpc of the cluster. The spatial distribution of our Coma active galaxies is consistent with that found at higher redshift for the CNOC1 cluster survey. Beyond the core region, the fraction of bright active galaxies appears to rise slowly out to the virial radius and are randomly distributed within the cluster with no apparent correlation with the potential merger of the NGC 4839 group. [ABRIDGED]

Abridged] We present here new results on the space density of rich, optically-selected, clusters of galaxies at low redshift (z<0.15). These results are based on the application of the matched filter cluster-finding algorithm to 1067 deg2 of the Edinburgh/Durham Southern Galaxy Catalogue (EDSGC). This is the first major application of this methodology at low redshift and in total, we have detected 2109 clusters above a richness cut-off of R_m >= 50 (or Lambda_cl >= 10). We have used extensive Monte Carlo simulations to define the detection thresholds for our algorithm, to measure the effective area of the EDCCII and to determine our spurious detection rate. These simulations have shown that our detection efficiency is strongly correlated with the presence of large-scale structure in the EDSGC data which we believe is due to the assumption of a flat, uniform background. Using these simulations, we are able to compute the space density of clusters in this new survey. We find 83.5 10-6 h-3 Mpc-3 for 100<= R_m < 200 (Lambda_cl~20) systems, 10.1 10-6 h-3 Mpc-3 for 200 <= R_m < 400 (Lambda_cl~40) systems and 2.3 10-6 h-3 Mpc-3 for R_m > 400 (Lambda_cl > 80) systems. These three richness bands roughly correspond to Abell Richness Classes 0, 1 and >=2 respectfully. These new measurements of the local space density of clusters are in agreement with those found at higher redshift (0.2 < z_est < 0.6) in the Palomar Distant Cluster Survey (PDCS; Postman et al. 1996 & Holden et al. 1999) thus confirming that the space density of clusters remains nearly constant out to z~0.6 in agreement with previous optical and X-ray measurements of the space density of clusters.