Search

The JWST North Ecliptic Pole (NEP) Time Domain Field (TDF) is a $>$14 arcmin diameter field optimized for multi-wavelength time-domain science with JWST. It has been observed across the electromagnetic spectrum both from the ground and from space, including with the Hubble Space Telescope (HST). As part of HST observations over 3 cycles (the "TREASUREHUNT" program), deep images were obtained with ACS/WFC in F435W and F606W that cover almost the entire JWST NEP TDF. Many of the individual pointings of these programs partially overlap, allowing an initial assessment of the potential of this field for time-domain science with HST and JWST. The cumulative area of overlapping pointings is ~88 arcmin$^2$, with time intervals between individual epochs that range between 1 day and 4$+$ years. To a depth of $m_{AB}$ $\simeq$ 29.5 mag (F606W), we present the discovery of 12 transients and 190 variable candidates. For the variable candidates, we demonstrate that Gaussian statistics are applicable, and estimate that ~80 are false positives. The majority of the transients will be supernovae, although at least two are likely quasars. Most variable candidates are AGN, where we find 0.42% of the general $z$ $<$ 6 field galaxy population to vary at the $~3\sigma$ level. Based on a 5-year timeframe, this translates into a random supernova areal density of up to ~0.07 transients per arcmin$^2$ (~245 deg$^{-2}$) per epoch, and a variable AGN areal density of ~1.25 variables per arcmin$^2$ (~4500 deg$^{-2}$) to these depths., Comment: 32 pages, 11 figures, 5 tables, 1 Appendix

Four ultra-luminous infrared galaxies (ULIRGs) observed with JWST/NIRcam in the Cosmos Evolution Early Release Science program offer an unbiased preview of the $z\approx2$ ULIRG population. The objects were originally selected at 24 $\mu$m and have strong polycyclic aromatic hydrocarbon emission features observed with Spitzer/IRS. The four objects have similar stellar masses of ${\sim}10^{11}$ M$_\odot$ but otherwise are quite diverse. One is an isolated disk galaxy, but it has an active nucleus as shown by X-ray observations and by a bright point-source nucleus. Two others are merging pairs with mass ratios of 6-7:1. One has active nuclei in both components, while the other has only one active nucleus: the one in the less-massive neighbor, not the ULIRG. The fourth object is clumpy and irregular and is probably a merger, but there is no sign of an active nucleus. The intrinsic spectral energy distributions for the four AGNs in these systems are typical of type-2 QSOs. This study is consistent with the idea that even if internal processes can produce large luminosities at $z\sim2$, galaxy merging may still be necessary for the most luminous objects. The diversity of these four initial examples suggests that large samples will be needed to understand the $z\approx2$ ULIRG population., Comment: 9 pages, 3 figures, 1 table, accepted for publication in ApJ. V2 updates author affiliations and acknowledgments, not scientific content

We report timing analysis of near-infrared (NIR), X-ray, and sub-millimeter (submm) data during a three-day coordinated campaign observing Sagittarius A*. Data were collected at 4.5 micron with the Spitzer Space Telescope, 2-8 keV with the Chandra X-ray Observatory, 3-70 keV with NuSTAR, 340 GHz with ALMA, and at 2.2 micron with the GRAVITY instrument on the Very Large Telescope Interferometer. Two dates show moderate variability with no significant lags between the submm and the infrared at 99% confidence. July 18 captured a moderately bright NIR flare (F_K ~ 15 mJy) simultaneous with an X-ray flare (F ~ 0.1 cts/s) that most likely preceded bright submm flux (F ~ 5.5 Jy) by about +34 (+14 -33) minutes at 99% confidence. The uncertainty in this lag is dominated by the fact that we did not observe the peak of the submm emission. A synchrotron source cooled through adiabatic expansion can describe a rise in the submm once the synchrotron-self-Compton NIR and X-ray peaks have faded. This model predicts high GHz and THz fluxes at the time of the NIR/X-ray peak and electron densities well above those implied from average accretion rates for Sgr A*. However, the higher electron density postulated in this scenario would be in agreement with the idea that 2019 was an extraordinary epoch with a heightened accretion rate. Since the NIR and X-ray peaks can also be fit by a non-thermal synchrotron source with lower electron densities, we cannot rule out an unrelated chance coincidence of this bright submm flare with the NIR/X-ray emission., Comment: Accepted for publication in The Astrophysical Journal

We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $\nu L_\nu\propto \nu^{0.5\pm0.2}$; the spectral slope observed in the X-ray band is $\nu L_\nu \propto \nu^{-0.7\pm0.5}$. We tested synchrotron and synchrotron self-Compton (SSC) scenarios. The observed near-infrared brightness and X-ray faintness, together with the observed spectral slopes, pose challenges for all models explored. We rule out a scenario in which the near-infrared emission is synchrotron emission and the X-ray emission is SSC. A one-zone model in which both the near-infrared and X-ray luminosity are produced by SSC and a model in which the luminosity stems from a cooled synchrotron spectrum can explain the flare. In order to describe the mean SED, both models require specific values of the maximum Lorentz factor $\gamma_{max}$, which however differ by roughly two orders of magnitude: the SSC model suggests that electrons are accelerated to $\gamma_{max}\sim 500$, while cooled synchrotron model requires acceleration up to $\gamma_{max}\sim5\times 10^{4}$. The SSC scenario requires electron densities of $10^{10}~\mathrm{cm^{-3}}$ much larger than typical ambient densities in the accretion flow, and thus require in an extraordinary accretion event. In contrast, assuming a source size of $1R_s$, the cooled synchrotron scenario can be realized with densities and magnetic fields comparable with the ambient accretion flow. For both models, the temporal evolution is regulated through the maximum acceleration factor $\gamma_{max}$, implying that sustained particle acceleration is required to explain at least a part of the temporal evolution of the flare., Comment: accepted for publication in Astronomy & Astrophysics; preview abstract shortened due to arXiv requirements

Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We have analyzed a comprehensive submillimeter (including new observations simultaneous with NIR monitoring), NIR, and 2-8 keV dataset. Submillimeter variations tend to lag those in the NIR by $\sim$30 minutes. An approximate Bayesian computation (ABC) fit to the X-ray first-order structure function shows significantly less power at short timescales in the X-rays than in the NIR. Less X-ray variability at short timescales combined with the observed NIR-X-ray correlations means the variability can be described as the result of two strictly correlated stochastic processes, the X-ray process being the low-pass-filtered version of the NIR process. The NIR--X-ray linkage suggests a simple radiative model: a compact, self-absorbed synchrotron sphere with high-frequency cutoff close to NIR frequencies plus a synchrotron self-Compton scattering component at higher frequencies. This model, with parameters fit to the submillimeter, NIR, and X-ray structure functions, reproduces the observed flux densities at all wavelengths, the statistical properties of all light curves, and the time lags between bands. The fit also gives reasonable values for physical parameters such as magnetic flux density $B\approx13$ G, source size $L \approx2.2R_{S}$, and high-energy electron density $n_{e}\approx4\times10^{7}$ cm$^{-3}$. An animation illustrates typical light curves, and we make public the parameter chain of our Bayesian analysis, the model implementation, and the visualization code., Comment: Accepted by ApJS, June 3, 2021. The journal version of figure 21 is animated. The animation can also be found here: https://doi.org/10.17617/1.kctx3s25. This version (version 2) has been revised according to the referee's suggestions and includes optimized figures and some changes and corrections of the text; no scientific conclusions have changed

We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively matched-aperture photometry across all bands, based upon object size information available in the catalog and Moffat profile point spread function fits. We test this catalog with a simple machine learning-based photometric redshift algorithm based upon Random Forest regression, and find that the corrected aperture photometry leads to significant improvement in photo-z accuracy compared to the original SExtractor catalogs from CFHTLS and Subaru. The deep ugrizY photometry and spectroscopic redshifts are well-suited for empirical tests of photometric redshift algorithms for LSST. The resulting catalogs are publicly available. We include a basic summary of the strategy of the DEEP3 Galaxy Redshift Survey to accompany the recent public release of DEEP3 data., Comment: 20 pages, 9 figures; published in MNRAS; for associated catalogs, see http://d-scholarship.pitt.edu/36064/

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

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

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science., Comment: White Paper in response to LSST Call for Observing Strategy Input

We present multi-wavelength global star formation rate (SFR) estimates for 326 galaxies from the Star Formation Reference Survey (SFRS) in order to determine the mutual scatter and range of validity of different indicators. The widely used empirical SFR recipes based on 1.4 GHz continuum, 8.0 $\mu$m polycyclic aromatic hydrocarbons (PAH), and a combination of far-infrared (FIR) plus ultraviolet (UV) emission are mutually consistent with scatter of $\raise{-0.8ex}\stackrel{\textstyle <}{\sim }$0.3 dex. The scatter is even smaller, $\raise{-0.8ex}\stackrel{\textstyle <}{\sim }$0.24 dex, in the intermediate luminosity range 9.3}5\times10^9$ L$_{\odot}$ or with implied far-UV extinction <1 mag, the UV spectral slope gives extinction corrections with 0.22~dex uncertainty., Comment: To be published in Monthly Notices of the Royal Astronomical Society. Data tables available at https://www.cfa.harvard.edu/irac/SFRS/ . V2 has minor clarifications from V1; no changes in data or results

Sgr A*, the supermassive black hole (SMBH) at the center of our Milky Way Galaxy, is known to be a variable source of X-ray, near-infrared (NIR), and submillimeter (submm) radiation and therefore a prime candidate to study the electromagnetic radiation generated by mass accretion flow onto a black hole and/or a related jet. Disentangling the power source and emission mechanisms of this variability is a central challenge to our understanding of accretion flows around SMBHs. Simultaneous multiwavelength observations of the flux variations and their time correlations can play an important role in obtaining a better understanding of possible emission mechanisms and their origin. This paper presents observations of two flares that both apparently violate the previously established patterns in the relative timing of submm/NIR/X-ray flares from Sgr A*. One of these events provides the first evidence of coeval structure between NIR and submm flux increases, while the second event is the first example of the sequence of submm/X-ray/NIR flux increases all occurring within ~1 hr. Each of these two events appears to upend assumptions that have been the basis of some analytic models of flaring in Sgr A*. However, it cannot be ruled out that these events, even though unusual, were just coincidental. These observations demonstrate that we do not fully understand the origin of the multiwavelength variability of Sgr A*, and show that there is a continued and important need for long-term, coordinated, and precise multiwavelength observations of Sgr A* to characterize the full range of variability behavior., Comment: 9 pages, 3 figures, to appear in the Astrophysical Journal

Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $\mu$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $\mu$m data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 $\mu$m data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current Fast Fourier Transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is $\tau_{b} = 243^{+82}_{-57}$ minutes ($90\%$ credible interval). The PSD has no detectable features on timescales down to 8.5 minutes ($95\%$ credible level), which is the ISCO orbital frequency for a dimensionless spin parameter $a = 0.92$. One light curve measured simultaneously at 2.12 and 4.5 $\mu$m during a low flux-density phase gave a spectral index $\alpha_s = 1.6 \pm 0.1$ ($F_\nu \propto \nu^{-\alpha_s}$). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined datasets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a non-thermal electron population from below 20 GHz through the NIR., Comment: Accepted for publication in ApJ on May 30, 2018. A machine readable version of the light curve data is included in the journal's online publication. Version 2 includes proof corrections

We present Spitzer Space Telescope 3.6 and 4.5 micron observations of the binary neutron star merger GW170817 at 43, 74, and 264 days post-merger. Using the final observation as a template, we uncover a source at the position of GW170817 at 4.5 micron with a brightness of 22.9+/-0.3 AB mag at 43 days and 23.8+/-0.3 AB mag at 74 days (the uncertainty is dominated by systematics from the image subtraction); no obvious source is detected at 3.6 micron to a 3-sigma limit of >23.3 AB mag in both epochs. The measured brightness is dimmer by a factor of about 2-3 times compared to our previously published kilonova model, which is based on UV, optical, and near-IR data at <30 days. However, the observed fading rate and color (m_{3.6}-m_{4.5}> 0 AB mag) are consistent with our model. We suggest that the discrepancy is likely due to a transition to the nebular phase, or a reduced thermalization efficiency at such late time. Using the Spitzer data as a guide, we briefly discuss the prospects of observing future binary neutron star mergers with Spitzer (in LIGO/Virgo Observing Run 3) and the James Webb Space Telescope (in LIGO/Virgo Observing Run 4 and beyond)., Comment: 6 pages, 2 figures, submitted to ApJL

VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection., Comment: 19 pages, 9 figures, accepted for publication in MNRAS

We study a large galaxy sample from the Spitzer Matching Survey of the UltraVISTA ultra-deep Stripes (SMUVS) to search for sources with enhanced 3.6 micron fluxes indicative of strong Halpha emission at z=3.9-4.9. We find that the percentage of "Halpha excess" sources reaches 37-40% for galaxies with stellar masses log10(M*/Msun) ~ 9-10, and decreases to <20% at log10(M*/Msun) ~ 10.7. At higher stellar masses, however, the trend reverses, although this is likely due to AGN contamination. We derive star formation rates (SFR) and specific SFR (sSFR) from the inferred Halpha equivalent widths (EW) of our "Halpha excess" galaxies. We show, for the first time, that the "Halpha excess" galaxies clearly have a bimodal distribution on the SFR-M* plane: they lie on the main sequence of star formation (with log10(sSFR/yr^{-1})<-8.05) or in a starburst cloud (with log10(sSFR/yr^{-1}) >-7.60). The latter contains ~15% of all the objects in our sample and accounts for >50% of the cosmic SFR density at z=3.9-4.9, for which we derive a robust lower limit of 0.066 Msun yr^{-1} Mpc^{-3}. Finally, we identify an unusual >50sigma overdensity of z=3.9-4.9 galaxies within a 0.20 x 0.20 sq. arcmin region. We conclude that the SMUVS unique combination of area and depth at mid-IR wavelengths provides an unprecedented level of statistics and dynamic range which are fundamental to reveal new aspects of galaxy evolution in the young Universe., Comment: 18 pages, 11 figures, 1 table. Re-submitted to the ApJ, after addressing referee report. Main changes with respect to v1: a new section and a new appendix have been added to investigate further the origin and robustness of the sSFR bimodality. No conclusion changed

We present Hubble WFC3/IR slitless grism spectra of a remarkably bright $z\gtrsim10$ galaxy candidate, GN-z11, identified initially from CANDELS/GOODS-N imaging data. A significant spectroscopic continuum break is detected at $\lambda=1.47\pm0.01~\mu$m. The new grism data, combined with the photometric data, rule out all plausible lower redshift solutions for this source. The only viable solution is that this continuum break is the Ly$\alpha$ break redshifted to ${z_\mathrm{grism}=11.09^{+0.08}_{-0.12}}$, just $\sim$400 Myr after the Big Bang. This observation extends the current spectroscopic frontier by 150 Myr to well before the Planck (instantaneous) cosmic reionization peak at z~8.8, demonstrating that galaxy build-up was well underway early in the reionization epoch at z>10. GN-z11 is remarkably and unexpectedly luminous for a galaxy at such an early time: its UV luminosity is 3x larger than L* measured at z~6-8. The Spitzer IRAC detections up to 4.5 $\mu$m of this galaxy are consistent with a stellar mass of ${\sim10^{9}~M_\odot}$. This spectroscopic redshift measurement suggests that the James Webb Space Telescope (JWST) will be able to similarly and easily confirm such sources at z>10 and characterize their physical properties through detailed spectroscopy. Furthermore, WFIRST, with its wide-field near-IR imaging, would find large numbers of similar galaxies and contribute greatly to JWST's spectroscopy, if it is launched early enough to overlap with JWST., Comment: 12 pages, 11 figures, accepted for publication in ApJ

We analyze the spatial distribution of dusty young stellar objects (YSOs) identified in the Spitzer Survey of the Orion Molecular clouds, augmenting these data with Chandra X-ray observations to correct for incompleteness in dense clustered regions. We also devise a scheme to correct for spatially varying incompleteness when X-ray data are not available. The local surface densities of the YSOs range from 1 pc$^{-2}$ to over 10,000 pc$^{-2}$, with protostars tending to be in higher density regions. This range of densities is similar to other surveyed molecular clouds with clusters, but broader than clouds without clusters. By identifying clusters and groups as continuous regions with surface densities $\ge10$ pc$^{-2}$, we find that 59% of the YSOs are in the largest cluster, the Orion Nebular Cluster (ONC), while 13% of the YSOs are found in a distributed population. A lower fraction of protostars in the distributed population is evidence that it is somewhat older than the groups and clusters. An examination of the structural properties of the clusters and groups show that the peak surface densities of the clusters increase approximately linearly with the number of members. Furthermore, all clusters with more than 70 members exhibit asymmetric and/or highly elongated structures. The ONC becomes azimuthally symmetric in the inner 0.1 pc, suggesting that the cluster is only $\sim 2$ Myr in age. We find the star formation efficiency (SFE) of the Orion B cloud is unusually low, and that the SFEs of individual groups and clusters are an order of magnitude higher than those of the clouds. Finally, we discuss the relationship between the young low mass stars in the Orion clouds and the Orion OB 1 association, and we determine upper limits to the fraction of disks that may be affected by UV radiation from OB stars or by dynamical interactions in dense, clustered regions., Comment: Accepted to AJ. Higher resolution and two column version of paper as well as photometry data available at http://astro1.physics.utoledo.edu/~megeath/Orion/The_Spitzer_Orion_Survey.html

We present mid-IR (19 - 37 microns) imaging observations of S106 from SOFIA/FORCAST, complemented with IR observations from Spitzer/IRAC (3.6 - 8.0 microns), IRTF/MIRLIN (11.3 and 12.5 microns), and Herschel/PACS (70 and 160 microns). We use these observations, observations in the literature, and radiation transfer modeling to study the heating and composition of the warm (~ 100 K) dust in the region. The dust is heated radiatively by the source S106 IR, with little contributions from grain-electron collisions and Ly-alpha radiation. The dust luminosity is >~ (9.02 +/- 1.01) x 10^4 L_sun, consistent with heating by a mid- to late-type O star. We find a temperature gradient (~ 75 - 107 K) in the lobes, which is consistent with a dusty equatorial geometry around S106 IR. Furthermore, the SOFIA observations resolve several cool (~ 65 - 70 K) lanes and pockets of warmer (~ 75 - 90 K) dust in the ionization shadow, indicating that the environment is fragmented. We model the dust mass as a composition of amorphous silicates, amorphous carbon, big grains, very small grains, and PAHs. We present the relative abundances of each grain component for several locations in S106., Comment: 18 pages, 11 figures, 4 tables, accepted by ApJ

We present Spitzer/Infrared Spectrograph (IRS) 5-21 micron spectroscopic maps towards 12 regions in the Andromeda galaxy (M31). These regions include the nucleus, bulge, an active region in the star-forming ring, and 9 other regions chosen to cover a range of mid-to-far-infrared colours. In line with previous results, PAH feature ratios (6.2 micron and 7.7 micron features compared to the 11.2 micron feature) measured from our extracted M31 spectra, except the nucleus, strongly correlate. The equivalent widths of the main PAH features, as a function of metallicity and radiation hardness, are consistent with those observed for other nearby spiral and starburst galaxies. Reprocessed data from the ISOCAM instrument on the Infrared Space Observatory agree with the IRS data; early reports of suppressed 6-8 micron features and enhanced 11.3 micron feature intensity and FWHM apparently resulted from background-subtraction problems. The nucleus does not show any PAH emission but does show strong silicate emission at 9.7 micron. Furthermore, different spectral features (11.3 micron PAH emission, silicate emission and [NeIII] 15.5 micron line emission) have distinct spatial distributions in the nuclear region: the silicate emission is strongest towards the stellar nucleus, while the PAH emission peaks 15 arcsec north of the nucleus. The PAH feature ratios at this position are atypical with strong emission at 11.2 microns and 15-20 microns but weak emission at 6--8 microns. The nucleus itself is dominated by stellar light giving rise to a strong blue continuum and silicate emission., Comment: 14 pages, 13 figures; MNRAS in press

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

The deep, wide-area (~800-900 arcmin**2) near-infrared/WFC3/IR + Spitzer/IRAC observations over the CANDELS fields have been a remarkable resource for constraining the bright end of high redshift UV luminosity functions (LFs). However, the lack of HST 1.05-micron observations over the CANDELS fields has made it difficult to identify z~9-10 sources robustly, since such data are needed to confirm the presence of an abrupt Lyman break at 1.2 microns. We report here on the successful identification of many such z~9-10 sources from a new HST program (z9-CANDELS) that targets the highest-probability z~9-10 galaxy candidates with observations at 1.05 microns, to search for a robust Lyman-break at 1.2 microns. The potential z~9-10 candidates are preselected from the full HST, Spitzer/IRAC S-CANDELS observations, and the deepest-available ground-based optical+near-infrared observations. We identified 15 credible z~9-10 galaxies over the CANDELS fields. Nine of these galaxies lie at z~9 and 5 are new identifications. Our targeted follow-up strategy has proven to be very efficient in making use of scarce HST time to secure a reliable sample of z~9-10 galaxies. Through extensive simulations, we replicate the selection process for our sample (both the preselection and follow-up) and use it to improve current estimates for the volume density of bright z~9 and z~10 galaxies. The volume densities we find are 5(-2)(+3)x and 8(-3)(+9)x lower, respectively, than found at z~8. When compared with the best-fit evolution (i.e., dlog_{10} rho(UV)/dz=-0.29+/-0.02) in the UV luminosities densities from z~8 to z~4 integrated to 0.3L*(z=3) (-20 mag), these luminosity densities are 2.6(-0.9)(+1.5)x and 2.2(-1.1)(+2.0)x lower, respectively, than the extrapolated trends. Our new results are broadly consistent with the "accelerated evolution" scenario at z>8, as seen in many theoretical models., Comment: 23 pages, 15 figures, 7 tables, updated to match the version in press, including some minor textual corrections identified at the proof stage

We present the results of deep \chandra\ imaging of the central region of the Extended Groth Strip, the AEGIS-X Deep (AEGIS-XD) survey. When combined with previous \chandra\ observations of a wider area of the strip, AEGIS-X Wide (AEGIS-XW; Laird et~al. 2009), these provide data to a nominal exposure depth of 800ks in the three central ACIS-I fields, a region of approximately $0.29$~deg$^{2}$. This is currently the third deepest X-ray survey in existence, a factor $\sim 2-3$ shallower than the Chandra Deep Fields (CDFs) but over an area $\sim 3$ times greater than each CDF. We present a catalogue of 937 point sources detected in the deep \chandra\ observations. We present identifications of our X-ray sources from deep ground-based, Spitzer, GALEX and HST imaging. Using a likelihood ratio analysis, we associate multi band counterparts for 929/937 of our X-ray sources, with an estimated 95~\% reliability, making the identification completeness approximately 94~\% in a statistical sense. Reliable spectroscopic redshifts for 353 of our X-ray sources are provided predominantly from Keck (DEEP2/3) and MMT Hectospec, so the current spectroscopic completeness is $\sim 38$~per cent. For the remainder of the X-ray sources, we compute photometric redshifts based on multi-band photometry in up to 35 bands from the UV to mid-IR. Particular attention is given to the fact that the vast majority the X-ray sources are AGN and require hybrid templates. Our photometric redshifts have mean accuracy of $\sigma=0.04$ and an outlier fraction of approximately 5\%, reaching $\sigma=0.03$ with less than 4\% outliers in the area covered by CANDELS . The X-ray, multi-wavelength photometry and redshift catalogues are made publicly available., Comment: 35 pages, 15 figures, accepted for publication in ApJS

We present Herschel (PACS and SPIRE) far-infrared (FIR) photometry of a complete sample of z>1 3CR sources, from the Herschel GT project The Herschel Legacy of distant radio-loud AGN (PI: Barthel). Combining these with existing Spitzer photometric data, we perform an infrared (IR) spectral energy distribution (SED) analysis of these landmark objects in extragalactic research to study the star formation in the hosts of some of the brightest active galactic nuclei (AGN) known at any epoch. Accounting for the contribution from an AGN-powered warm dust component to the IR SED, about 40% of our objects undergo episodes of prodigious, ULIRG-strength star formation, with rates of hundreds of solar masses per year, coeval with the growth of the central supermassive black hole. Median SEDs imply that the quasar and radio galaxy hosts have similar FIR properties, in agreement with the orientation-based unification for radio-loud AGN. The star-forming properties of the AGN hosts are similar to those of the general population of equally massive non-AGN galaxies at comparable redshifts, thus there is no strong evidence of universal quenching of star formation (negative feedback) within this sample. Massive galaxies at high redshift may be forming stars prodigiously, regardless of whether their supermassive black holes are accreting or not., Comment: 30 pages, 13 figures, 4 tables. Accepted for publication in A&A

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.54. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts., Comment: 26 pages, 23 figures. Accepted for publication on A&A

We present the first detection from the Spitzer Space Telescope of 4.5 micron variability from Sgr A*, the emitting source associated with the Milky Way's central black hole. The >23 hour continuous light curve was obtained with the IRAC instrument in 2013 December. The result characterizes the variability of Sgr A* prior to the closest approach of the G2 object, a putative infalling gas cloud that orbits close to Sgr A*. The high stellar density at the location of Sgr A* produces a background of ~250 mJy at 4.5 microns in each pixel with a large pixel-to-pixel gradient, but the light curve for the highly variable Sgr A* source was successfully measured by modeling and removing the variations due to pointing wobble. The observed flux densities range from the noise level of ~0.7 mJy rms in a 6.4-s measurement to ~10 mJy. Emission was seen above the noise level ~34% of the time. The light curve characteristics, including the flux density distribution and structure function, are consistent with those previously derived at shorter infrared wavelengths. We see no evidence in the light curve for activity attributable to the G2 interaction at the observing epoch, ~100 days before the expected G2 periapsis passage. The IRAC light curve is more than a factor of two longer than any previous infrared observation, improving constraints on the timescale of the break in the power spectral distribution of Sgr A* flux densities. The data favor the longer of the two previously published values for the timescale., Comment: 13 pages, 10 figures, 2 tables, accepted for publication in the ApJ

We report on observations of near-Earth asteroid 2011 MD with the Spitzer Space Telescope. We have spent 19.9 h of observing time with channel 2 (4.5 {\mu}m) of the Infrared Array Camera and detected the target within the 2{\sigma} positional uncertainty ellipse. Using an asteroid thermophysical model and a model of nongravitational forces acting upon the object we constrain the physical properties of 2011 MD, based on the measured flux density and available astrometry data. We estimate 2011 MD to be 6 (+4/-2) m in diameter with a geometric albedo of 0.3 (+0.4/-0.2) (uncertainties are 1{\sigma}). We find the asteroid's most probable bulk density to be 1.1 (+0.7/-0.5) g cm^{-3}, which implies a total mass of (50-350) t and a macroporosity of >=65%, assuming a material bulk density typical of non-primitive meteorite materials. A high degree of macroporosity suggests 2011 MD to be a rubble-pile asteroid, the rotation of which is more likely to be retrograde than prograde., Comment: 20 pages, 4 figures

We report on Spitzer Space Telescope IRAC observations of near-Earth object (NEO) 2009 BD that were carried out in support of the NASA Asteroid Robotic Retrieval Mission (ARRM) concept. We did not detect 2009 BD in 25 hrs of integration at 4.5 micron. Based on an upper-limit flux density determination from our data, we present a probabilistic derivation of the physical properties of this object. The analysis is based on the combination of a thermophysical model with an orbital model accounting for the non-gravitational forces acting upon the body. We find two physically possible solutions. The first solution shows 2009 BD as a 2.9+/-0.3 m diameter rocky body (rho = 2.9+/-0.5 g cm-3) with an extremely high albedo of 0.85(+0.20/-0.10) that is covered with regolith-like material, causing it to exhibit a low thermal inertia (Gamma = 30(+20/-10) SI units). The second solution suggests 2009 BD to be a 4+/-1 m diameter asteroid with pV = 0.45(+0.35/-0.15) that consists of a collection of individual bare rock slabs (Gamma = 2000+/-1000 SI units, rho = 1.7(+0.7/-0.4) g cm-3). We are unable to rule out either solution based on physical reasoning. 2009 BD is the smallest asteroid for which physical properties have been constrained, in this case using an indirect method and based on a detection limit, providing unique information on the physical properties of objects in the size range smaller than 10 m., Comment: 28 pages, 8 figures, accepted for publication in ApJ

We present a study of the infrared properties for a sample of seven spectroscopically confirmed submillimeter galaxies at $z>$4.0. By combining ground-based near-infrared, Spitzer IRAC and MIPS, Herschel SPIRE, and ground-based submillimeter/millimeter photometry, we construct their Spectral Energy Distributions (SED) and a composite model to fit the SEDs. The model includes a stellar emission component at $\lambda_{\rm rest} <$ 3.5$ \mu$m; a hot dust component peaking at $\lambda_{rest} \sim$ 5$\,\mu$m; and cold dust component which becomes significant for $\lambda_{\rm rest} >$ 50$\,\mu$m. Six objects in the sample are detected at 250 and 350$ \mu$m. The dust temperatures for the sources in this sample are in the range of 40$-$80 K, and their $L_{\rm FIR}$ $\sim$ 10$^{13}$ L$_{\odot}$ qualifies them as Hyper$-$Luminous Infrared Galaxies (HyperLIRGs). The mean FIR-radio index for this sample is around $< q > = 2.2$ indicating no radio excess in their radio emission. Most sources in the sample have 24$ \mu$m detections corresponding to a rest-frame 4.5$ \mu$m luminosity of Log$_{10}$(L$_{4.5}$ / L$_{\odot}$) = 11 $\sim$ 11.5. Their L$_{\rm 4.5}$/$L_{\rm FIR}$ ratios are very similar to those of starburst dominated submillimeter galaxies at $z \sim$ 2. The $L_{\rm CO}-L_{\rm FIR}$ relation for this sample is consistent with that determined for local ULIRGs and SMGs at $z \sim$ 2. We conclude that submillimeter galaxies at $z >$ 4 are hotter and more luminous in the FIR, but otherwise very similar to those at $z \sim$ 2. None of these sources show any sign of the strong QSO phase being triggered., Comment: 49 pages, 15 figures. Submitted to ApJ

Out of several dozen z > 7 candidate galaxies observed spectroscopically, only five have been confirmed via Lyman-alpha emission, at z=7.008, 7.045, 7.109, 7.213 and 7.215. The small fraction of confirmed galaxies may indicate that the neutral fraction in the intergalactic medium (IGM) rises quickly at z > 6.5, as Lyman-alpha is resonantly scattered by neutral gas. However, the small samples and limited depth of previous observations makes these conclusions tentative. Here we report the results of a deep near-infrared spectroscopic survey of 43 z > 6.5 galaxies. We detect only a single galaxy, confirming that some process is making Lyman-alpha difficult to detect. The detected emission line at 1.0343 um is likely to be Lyman-alpha emission, placing this galaxy at a redshift z = 7.51, an epoch 700 million years after the Big Bang. This galaxy's colors are consistent with significant metal content, implying that galaxies become enriched rapidly. We measure a surprisingly high star formation rate of 330 Msol/yr, more than a factor of 100 greater than seen in the Milky Way. Such a galaxy is unexpected in a survey of our size, suggesting that the early universe may harbor more intense sites of star-formation than expected., Comment: To be published in Nature on 24 October. Under press embargo until 1pm EDT on 23 October. 12 pages, 8 figures, 4 tables (including supplementary information)

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

[abridged] We present the discovery of four surprisingly bright (H_160 ~ 26 - 27 mag AB) galaxy candidates at z~9-10 in the complete HST CANDELS WFC3/IR GOODS-N imaging data, doubling the number of z~10 galaxy candidates that are known, just ~500 Myr after the Big Bang. Two similarly bright sources are also detected in a systematic re-analysis of the GOODS-S data set. Three of the four galaxies in GOODS-N are significantly detected at 4.5-6.2sigma in the very deep Spitzer/IRAC 4.5 micron data, as is one of the GOODS-S candidates. Furthermore, the brightest of our candidates (at z=10.2+-0.4) is robustly detected also at 3.6 micron (6.9sigma), revealing a flat UV spectral energy distribution with a slope beta=-2.0+-0.2, consistent with demonstrated trends with luminosity at high redshift. The abundance of such luminous candidates suggests that the luminosity function evolves more significantly in phi_* than in L_* at z>~8 with a higher number density of bright sources than previously expected. Despite the discovery of these luminous candidates, the cosmic star formation rate density for galaxies with SFR >0.7 M_sun/yr shows an order-of-magnitude increase in only 170 Myr from z ~ 10 to z ~ 8, consistent with previous results. Based on the IRAC detections, we derive galaxy stellar masses at z~10, finding that these luminous objects are typically 10^9 M_sun. The cosmic stellar mass density at z~10 is log10 rho_* = 4.7^+0.5_-0.8 M_sun Mpc^-3 for galaxies brighter than M_UV~-18. The remarkable brightness, and hence luminosity, of these z~9-10 candidates highlights the opportunity for deep spectroscopy to determine their redshift and nature, demonstrates the value of additional search fields covering a wider area to understand star-formation in the very early universe, and highlights the opportunities for JWST to map the buildup of galaxies at redshifts much earlier than z~10., Comment: 20 pages, 12 figures, changed to match resubmitted version to ApJ

We present a UV-to-mid infrared multi-wavelength catalog in the CANDELS/GOODS-S field, combining the newly obtained CANDELS HST/WFC3 F105W, F125W, and F160W data with existing public data. The catalog is based on source detection in the WFC3 F160W band. The F160W mosaic includes the data from CANDELS deep and wide observations as well as previous ERS and HUDF09 programs. The mosaic reaches a 5$\sigma$ limiting depth (within an aperture of radius 0.17 arcsec) of 27.4, 28.2, and 29.7 AB for CANDELS wide, deep, and HUDF regions, respectively. The catalog contains 34930 sources with the representative 50% completeness reaching 25.9, 26.6, and 28.1 AB in the F160W band for the three regions. In addition to WFC3 bands, the catalog also includes data from UV (U-band from both CTIO/MOSAIC and VLT/VIMOS), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), and infrared (HST/WFC3 F098M, VLT/ISAAC Ks, VLT/HAWK-I Ks, and Spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) observations. The catalog is validated via stellar colors, comparison with other published catalogs, zeropoint offsets determined from the best-fit templates of the spectral energy distribution of spectroscopically observed objects, and the accuracy of photometric redshifts. The catalog is able to detect unreddened star-forming (passive) galaxies with stellar mass of 10^{10}M_\odot at a 50% completeness level to z$\sim$3.4 (2.8), 4.6 (3.2), and 7.0 (4.2) in the three regions. As an example of application, the catalog is used to select both star-forming and passive galaxies at z$\sim$2--4 via the Balmer break. It is also used to study the color--magnitude diagram of galaxies at 0

We report deep ALMA observations complemented with associated HST imaging for a luminous (m_uv=25) galaxy, `Himiko', at a redshift z=6.595. The galaxy is remarkable for its high star formation rate, 100 Mo/yr, securely estimated from our deep HST and Spitzer photometry, and the absence of any evidence for strong AGN activity or gravitational lensing magnification. Our ALMA observations probe an order of magnitude deeper than previous IRAM observations, yet fail to detect a 1.2mm dust continuum, indicating a flux <52uJy comparable with or weaker than that of local dwarf irregulars with much lower star formation rates. We likewise provide a strong upper limit for the flux of [CII] 158um, L([CII]) < 5.4x10^7 Lo, a diagnostic of the hot interstellar gas often described as a valuable probe for early galaxies. In fact, our observations indicate Himiko lies off the local L([CII]) - star formation rate scaling relation by a factor of more than 30. Both aspects of our ALMA observations suggest Himiko is an unique object with a very low dust content and perhaps nearly primordial interstellar gas. Our HST images provide unique insight into the morphology of this remarkable source, highlighting an extremely blue core of activity and two less extreme associated clumps. Himiko is undergoing a triple major merger event whose extensive ionized nebula of Lyman alpha emitting gas, discovered in our earlier work with Subaru, is powered by star formation and the dense circum-galactic gas. We are likely witnessing an early massive galaxy during a key period of its mass assembly close to the end of the reionization era., Comment: 13 pages, 10 figures. Accepted for publication in ApJ. Revised [CII] detection limit, added Figure 8, conclusions unchanged

Chandra observations of a complete, flux-limited sample of 38 high-redshift (10) indicating obscuration (log N_H ~ 22-24 cm^-2). These properties and the correlation between obscuration and radio core-fraction are consistent with orientation-dependent obscuration as in Unification models. About half the NLRGs have soft X-ray hardness ratios and/or high [OIII] emission line to X-ray luminosity ratio suggesting obscuration by Compton thick (CT) material so that scattered nuclear or extended X-ray emission dominates (as in NGC1068). The ratios of unobscured to Compton-thin (10^{22}< N_H(int) < 1.5 x 10^{24} cm^-2) to CT (N_H(int) > 1.5 x 10^{24} cm^-2) is 2.5:1.4:1 in this high luminosity, radio-selected sample. The obscured fraction is 0.5, higher than is typically reported for AGN at comparable luminosities from multi-wavelength surveys (0.1-0.3). Assuming random nuclear orientation, the unobscured half-opening angle of the disk/wind/torus structure is ~ 60deg and the obscuring material covers 30deg of which ~ 12deg is Compton thick. The multi-wavelength properties reveal that many NLRGs have intrinsic absorption 10-1000x higher than indicated by their X-ray hardness ratios, and their true L_X values are ~10--100x larger than the hardness-ratio absorption corrections would indicate., Comment: 42 pages, 11 figures, 3 tables; ApJ, in press

We present the multiwavelength - ultraviolet to mid-infrared - catalog of the UKIDSS Ultra-Deep Survey (UDS) field observed as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Based on publicly available data, the catalog includes: the CANDELS data from the Hubble Space Telescope (near-infrared WFC3 F125W and F160W data and visible ACS F606W and F814W data), u-band data from CFHT/Megacam, B, V, Rc, i' and z' band data from Subaru/Suprime-Cam, Y and Ks band data from VLT/HAWK-I, J, H and K bands data from UKIDSS (Data Release 8), and Spitzer/IRAC data (3.6, 4.5 from SEDS, 5.8 and 8.0um from SpUDS). The present catalog is F160W-selected and contains 35932 sources over an area of 201.7 square arcmin and includes radio and X-ray detected sources and spectroscopic redshifts available for 210 sources., Comment: 22 pages, 20 figures, ApJ in press. The CANDELS UDS catalog (in ascii & fits format) and associated ReadMe file are made available on the CANDELS website (http://candels.ucolick.org/) as well as on the MAST CANDELS page (http://archive.stsci.edu/prepds/candels/)

Star-forming galaxies at redshifts z>6 are likely responsible for the reionization of the universe, and it is important to study the nature of these galaxies. We present three candidates for z~7 Lyman-break galaxies (LBGs) from a 155 arcmin^2 area in the CANDELS/COSMOS field imaged by the deep FourStar Galaxy Evolution (zFourGE) survey. The FourStar medium-band filters provide the equivalent of R~10 spectroscopy, which cleanly distinguishes between z~7 LBGs and brown dwarf stars. The distinction between stars and galaxies based on an object's angular size can become unreliable even when using HST imaging; there exists at least one very compact z~7 candidate (FWHM~0.5-1 kpc) that is indistinguishable from a point source. The medium-band filters provide narrower redshift distributions compared with broad-band-derived redshifts. The UV luminosity function derived using the three z~7 candidates is consistent with previous studies, suggesting an evolution at the bright end (MUV -21.6 mag) from z~7 to z~5. Fitting the galaxies' spectral energy distributions, we predict Lyman-alpha equivalent widths for the two brightest LBGs, and find that the presence of a Lyman-alpha line affects the medium-band flux thereby changing the constraints on stellar masses and UV spectral slopes. This illustrates the limitations of deriving LBG properties using only broad-band photometry. The derived specific star-formation rates for the bright LBGs are ~13 per Gyr, slightly higher than the lower-luminosity LBGs, implying that the star-formation rate increases with stellar mass for these galaxies., Comment: Accepted in ApJ

We present a multi-wavelength study of a 3.6 $\mu$m-selected galaxy sample in the Extended Groth strip. The sample is complete for galaxies with stellar mass $>10^{9.5}$ \Msun and redshift $0.4}3\sigma$ errors) is about 3.5%. A new method of validation based on pair statistics confirms the estimate of standard deviation even for galaxies lacking spectroscopic redshifts. Basic galaxy properties measured include rest-frame $U-B$ colors, $B$- and $K$-band absolute magnitudes, and stellar masses. We divide the sample into quiescent and star-forming galaxies according to their rest-frame $U-B$ colors and 24 to 3.6 \micron\ flux density ratios and derive rest $K$-band luminosity functions and stellar mass functions for quiescent, star forming, and all galaxies. The results show that massive, quiescent galaxies were in place by $z\approx1$, but lower mass galaxies generally ceased their star formation at later epochs., Comment: 64 pages, 24 figures

In order to understand the nature of the sources producing the recently uncovered cosmic infrared background (CIB) fluctuations, we study cross-correlations between the fluctuations in the source-subtracted CIB from Spitzer/IRAC data and the unresolved cosmic X-ray background from deep Chandra observations. Our study uses data from the EGS/AEGIS field, where both data sets cover an 8'x45' region of the sky. Our measurement is the cross-power spectrum between the IR and X-ray data. The cross-power signal between the IRAC maps at 3.6 um and 4.5 um and the Chandra [0.5-2] keV data has been detected, at angular scales 20'', with an overall significance of 3.8 sigma and 5.6 sigma, respectively. At the same time we find no evidence of significant cross-correlations at the harder Chandra bands. The cross-correlation signal is produced by individual IR sources with 3.6 um and 4.5 um magnitudes m_AB 25-26 and [0.5-2] keV X-ray fluxes 7e-17 cgs.We determine that at least 15%-25% of the large scale power of the CIB fluctuations is correlated with the spatial power spectrum of the X-ray fluctuations. If this correlation is attributed to emission from accretion processes at both IR and X-ray wavelengths, this implies a much higher fraction of accreting black holes than among the known populations. We discuss the various possible origins for the cross-power signal and show that neither local foregrounds nor the known remaining normal galaxies and active galactic nuclei can reproduce the measurements. These observational results are an important new constraint on theoretical modeling of the near-IR CIB fluctuations., Comment: 28 pages, 9 Figures, Accepted by ApJ

We present a survey of the Orion A and B molecular clouds undertaken with the IRAC and MIPS instruments onboard Spitzer. In total, five distinct fields were mapped covering 9 sq. degrees in five mid-IR bands spanning 3-24 microns. The survey includes the Orion Nebula Cluster, the Lynds 1641, 1630 and 1622 dark clouds, and the NGC 2023, 2024, 2068 and 2071 nebulae. These data are merged with the 2MASS point source catalog to generate a catalog of eight band photometry. We identify 3479 dusty young stellar objects (YSOs) in the Orion molecular clouds by searching for point sources with mid-IR colors indicative of reprocessed light from dusty disks or infalling envelopes. The YSOs are subsequently classified on the basis of their mid-IR colors and their spatial distributions are presented. We classify 2991 of the YSOs as pre-main sequence stars with disks and 488 as likely protostars. Most of the sources were observed with IRAC in 2-3 epochs over 6 months; we search for variability between the epochs by looking for correlated variability in the 3.6 and 4.5 micron bands. We find that 50% of the dusty YSOs show variability. The variations are typically small (0.2 mag.) with the protostars showing a higher incidence of variability and larger variations. The observed correlations between the 3.6, 4.5, 5.8 and 8 micron variability suggests that we are observing variations in the heating of the inner disk due to changes in the accretion luminosity or rotating accretion hot spots., Comment: Accepted to the Astronomical Journal. A full resolution version of the paper and the point source catalog can be found at http://astro1.physics.utoledo.edu/~megeath/Orion/The_Spitzer_Orion_Survey.html

We report the results of a study exploring the stellar populations of 13 luminous (L>L*), spectroscopically confirmed, galaxies in the redshift interval 5.5= 300 Myr, the degeneracies introduced by dust extinction mean that only two of these objects actually require a >300 Myr old stellar population to reproduce the observed photometry. Moreover, when considering only smoothly-varying star-formation histories, we observe a clear tension between the data and models such that a galaxy SED template with an old age is often chosen in order to try and fit objects with blue UV-slopes but red UV-to-optical colours. To break this tension we explore SED fitting with two-component models (burst plus on-going star-formation) and allow for nebular emission. On average, the inclusion of nebular emission leads to lower stellar-mass estimates (median offset 0.18 dex), moderately higher specific star-formation rates, and allows for a wider range of plausible stellar ages. However, based on our SED modelling, we find no strong evidence for extremely young ages in our sample (<50 Myr). Finally, considering all of the different star-formation histories explored, we find that the median best-fitting ages are of the order 200-300 Myr and that the objects with the tightest constraints indicate ages in the range 50-200 Myr (Abridged)., Comment: 23 pages, 12 figures, for submission to MNRAS

We present a study on Spectral Energy Distributions, Morphologies, and star formation for an IRAC-selected extremely red object sample in the GOODS Chandra Deep Field-South. This work was enabled by new HST/WFC3 near-IR imaging from the CANDELS survey as well as the deepest available X-ray data from Chandra 4 Ms observations. This sample consists of 133 objects with the 3.6um limiting magnitude of [3.6] = 21.5, and is approximately complete for galaxies with M >10^{11}M_sun at 1.5 < z < 2.5. We classify this sample into two types, quiescent and star-forming galaxies, in the observed infrared color-color ([3.6]-[24] vs K-[3.6]) diagram. The further morphological study of this sample show a consistent result with the observed color classification. The classified quiescent galaxies are bulge dominated and star-forming galaxies in the sample have disk or irregular morphologies. Our observed infrared color classification is also consistent with the rest-frame color (U-V vs V-J) classification. We also found that quiescent and star-forming galaxies are well separated in the nonparametric morphology parameter (Gini vs M_{20}) diagram measuring their concentration and clumpiness: quiescent galaxies have Gini coefficient higher than 0.58 and star forming galaxies have Gini coefficient lower that 0.58. We argue that the star formation quenching process must lead to or be accompanied by the increasing galaxy concentration. One prominent morphological feature of this sample is that disks are commonly seen in this massive galaxy sample at 1.5 < z < 2.5: 30% of quiescent galaxies and 70% of star forming galaxies with M >10^{11}M_sun have disks in their rest-frame optical morphologies. The prevalence of these extended, relatively undisturbed disks challenges the merging scenario as the main mode of massive galaxy formation., Comment: 36 pages,16 figures,Accepted by ApJ

We present Spitzer IRAC and MIPS observations of the star-forming region containing intermediate-mass young stellar object (YSO) AFGL 490. We supplement these data with near-IR 2MASS photometry and with deep SQIID observations off the central high extinction region. We have more than doubled the known membership of this region to 57 Class I and 303 Class II YSOs via the combined 1-24 um photometric catalog derived from these data. We construct and analyze the minimum spanning tree of their projected positions, isolating one locally over-dense cluster core containing 219 YSOs (60.8% of the region's members). We find this cluster core to be larger yet less dense than similarly analyzed clusters. Although the structure of this cluster core appears irregular, we demonstrate that the parsec-scale surface densities of both YSOs and gas are correlated with a power law slope of 2.8, as found for other similarly analyzed nearby molecular clouds. We also explore the mass segregation implications of AFGL 490's offset from the center of its core, finding that it has no apparent preferential central position relative to the low-mass members., Comment: 44 pages, 13 figures, accepted to ApJ

We extend previous measurements of cosmic infrared background (CIB) fluctuations to ~ 1 deg using new data from the Spitzer Extended Deep Survey. Two fields, with depths of ~12 hr/pixel over 3 epochs, are analyzed at 3.6 and 4.5 mic. Maps of the fields were assembled using a self-calibration method uniquely suitable for probing faint diffuse backgrounds. Resolved sources were removed from the maps to a magnitude limit of AB mag ~ 25, as indicated by the level of the remaining shot noise. The maps were then Fourier-transformed and their power spectra were evaluated. Instrumental noise was estimated from the time-differenced data, and subtracting this isolates the spatial fluctuations of the actual sky. The power spectra of the source-subtracted fields remain identical (within the observational uncertainties) for the three epochs indicating that zodiacal light contributes negligibly to the fluctuations. Comparing to 8 mic power spectra shows that Galactic cirrus cannot account for the fluctuations. The signal appears isotropically distributed on the sky as required for an extragalactic origin. The CIB fluctuations continue to diverge to > 10 times those of known galaxy populations on angular scales out to < 1 deg. The low shot noise levels remaining in the diffuse maps indicate that the large scale fluctuations arise from the spatial clustering of faint sources well below the confusion noise. The spatial spectrum of these fluctuations is in reasonable agreement with an origin in populations clustered according to the standard cosmological model (LCDM) at epochs coinciding with the first stars era., Comment: ApJ, to be published

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

We found the first direct evidence that the Cepheid class namesake, delta Cephei, is currently losing mass. These observations are based on data obtained with the Spitzer Space Telescope in the infrared, and with the Very Large Array in the radio. We found that delta Cephei is associated with a vast circumstellar structure, reminiscent of a bow shock. This structure is created as the wind from the star interacts with the local interstellar medium. We measure an outflow velocity of ~35.5 km/s and a mass loss rate of ~1E-7 - 1E-6 Mo/yr. The very low dust content of the outflow suggests that the wind is possibly pulsation-driven, rather than dust-driven as common for other classes of evolved stars., Comment: Proceeding submitted for the 20th Stellar Pulsation Conference (Impact of New Instruments and Insights in Stellar Pulsation), Granada, Spain, September 2011

We report detection of four IRAC sources in the GOODS-South field with an extremely red color of H$-$[3.6]$>$4.5. The four sources are not detected in the deep HST WFC3 H-band image with H$_{limit}$=28.3 mag. We find that only 3 types of SED templates can produce such a red H$-$[3.6] color: a very dusty SED with the Calzetti extinction of A$_V$=16 mag at z=0.8; a very dusty SED with the SMC extinction of A$_V$=8 mag at z=2.0$\sim$2.2; and an 1Gyr SSP with A$_V \sim$0.8 at z=5.7. We argue that these sources are unlikely dusty galaxies at z$\leq$2.2 based on absent strong MIPS 24$\mu$m emission. The old stellar population model at z$>$4.5 remains a possible solution for the 4 sources. At z$>$4.5, these sources have stellar masses of Log(M$_*$/M$_{\odot}$)=10.6$\sim$11.2. One source, ERS-1, is also a type-II X-ray QSO with L$_{2-8keV}$=1.6$\times 10^{44}$ erg s$^{-1}$. One of the four sources is an X-ray QSO and another one is a HyperLIRG, suggesting a galaxy-merging scenario for the formation of these massive galaxies at high redshifts., Comment: 15 pages, 3 figures, accepted by ApJ Lett

Chandra X-ray observations of the high redshift (z =1.532) radio-loud quasar 3C270.1 in 2008 February show the nucleus to have a power-law spectrum, Gamma = 1.66 +/- 0.08, typical of a radio-loud quasar, and a marginally-detected Fe Kalpha emission line. The data also reveal extended X-ray emission, about half of which is associated with the radio emission from this source. The southern emission is co-spatial with the radio lobe and peaks at the position of the double radio hotspot. Modeling this hotspot including Spitzer upper limits rules out synchrotron emission from a single power-law population of electrons, favoring inverse-Compton emission with a field of ~11nT, roughly a third of the equipartition value. The northern emission is concentrated close to the location of a 40 deg. bend where the radio jet is presumed to encounter external material. It can be explained by inverse Compton emission involving Cosmic Microwave Background photons with a field of ~3nT, roughly a factor of nine below the equipartition value. The remaining, more diffuse X-ray emission is harder (HR=-0.09 +/- 0.22). With only 22.8+/-5.6 counts, the spectral form cannot be constrained. Assuming thermal emission with a temperature of 4 keV yields an estimate for the luminosity of 1.8E44 erg/s, consistent with the luminosity-temperature relation of lower-redshift clusters. However deeper Chandra X-ray observations are required to delineate the spatial distribution, and better constrain the spectrum of the diffuse emission to verify that we have detected X-ray emission from a high-redshift cluster., Comment: 20 pages, 7 figures, accepted by ApJ

Star formation is arguably the most important physical process in the cosmos. It is a fundamental driver of galaxy evolution and the ultimate source of most of the energy emitted by galaxies. A correct interpretation of star formation rate (SFR) measures is therefore essential to our understanding of galaxy formation and evolution. Unfortunately, however, no single SFR estimator is universally available or even applicable in all circumstances: the numerous galaxies found in deep surveys are often too faint (or too distant) to yield significant detections with most standard SFR measures, and until now there have been no global, multi-band observations of nearby galaxies that span all the conditions under which star-formation is taking place. To address this need in a systematic way, we have undertaken a multi-band survey of all types of star-forming galaxies in the local Universe. This project, the Star Formation Reference Survey (SFRS), is based on a statistically valid sample of 369 nearby galaxies that span all existing combinations of dust temperature, SFR, and specific SFR. Furthermore, because the SFRS is blind with respect to AGN fraction and environment it serves as a means to assess the influence of these factors on SFR. Our panchromatic global flux measurements (including GALEX FUV+NUV, SDSS ugriz, 2MASS JHKs, Spitzer 3-8{\mu}m, and others) furnish uniform SFR measures and the context in which their reliability can be assessed. This paper describes the SFRS survey strategy, defines the sample, and presents the multi-band photometry collected to date., Comment: 48 pages, 12 figures, 10 tables. Accepted by PASP. This version edited to correct references and typographical errors

We present the discovery of an infrared nebula around the Cepheid prototype delta Cephei and its hot companion HD 213307. Large scale (~2.1x10^4 AU) nebulosity is detected at 5.8, 8.0, 24 and 70 um. Surrounding the two stars, the 5.8 and 8.0 um emission is largely attributable to Polycyclic Aromatic Hydrocarbon (PAH) emission swept from the ISM by a wind originating from delta Cephei and/or its companion. Stochastically heated small dust grains are the most likely source of the 24 and 70 um extended emission. The 70 um emission, in particular, resembles a bow shock aligned in the direction of the proper motion of delta Cephei. This discovery supports the hypothesis that delta Cephei may be currently losing mass, at a rate in the range ~ 5x10^-9 to 6x10^-8 Mo/yr., Comment: Published on The Astrophysical Journal