Source author record

G. Helou

G. Helou appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

108works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

108 published item(s)

preprint2020arXiv

Candidate Electromagnetic Counterpart to the Binary Black Hole Merger Gravitational Wave Event S190521g

We report the first plausible optical electromagnetic (EM) counterpart to a (candidate) binary black hole (BBH) merger. Detected by the Zwicky Transient Facility (ZTF), the EM flare is consistent with expectations for a kicked BBH merger in the accretion disk of an active galactic nucleus (AGN), and is unlikely ($<O(0.01\%$)) due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernovae and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing or a tidal disruption event, are ruled out or constrained to be $<O(0.1\%$). If the flare is associated with S190521g, we find plausible values of: total mass $ M_{\rm BBH} \sim 100 M_{\odot}$, kick velocity $v_k \sim 200\, {\rm km}\, {\rm s}^{-1}$ at $θ\sim 60^{\circ}$ in a disk with aspect ratio $H/a \sim 0.01$ (i.e., disk height $H$ at radius $a$) and gas density $ρ\sim 10^{-10}\, {\rm g}\, {\rm cm}^{-3}$. The merger could have occurred at a disk migration trap ($a \sim 700\, r_{g}$; $r_g \equiv G M_{\rm SMBH} / c^2$, where $M_{\rm SMBH}$ is the mass of the AGN supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a re-encountering with the disk in $\sim 1.6\, {\rm yr}\, (M_{\rm SMBH}/10^{8}M_{\odot})\, (a/10^{3}r_{g})^{3/2}$.

preprint2020arXiv

Joint Survey Processing of Euclid, Rubin and Roman: Final Report

The Euclid, Rubin/LSST and Roman (WFIRST) projects will undertake flagship optical/near-infrared surveys in the next decade. By mapping thousands of square degrees of sky and covering the electromagnetic spectrum between 0.3 and 2 microns with sub-arcsec resolution, these projects will detect several tens of billions of sources, enable a wide range of astrophysical investigations by the astronomical community and provide unprecedented constraints on the nature of dark energy and dark matter. The ultimate cosmological, astrophysical and time-domain science yield from these missions will require joint survey processing (JSP) functionality at the pixel level that is outside the scope of the individual survey projects. The JSP effort scoped here serves two high-level objectives: 1) provide precise concordance multi-wavelength images and catalogs over the entire sky area where these surveys overlap, which accounts for source confusion and mismatched isophotes, and 2) provide a science platform to analyze concordance images and catalogs to enable a wide range of astrophysical science goals to be formulated and addressed by the research community. For the cost of about 200WY, JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also allow for the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone.

preprint2020arXiv

The Koala: A Fast Blue Optical Transient with Luminous Radio Emission from a Starburst Dwarf Galaxy at $z=0.27$

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

preprint2019arXiv

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: The KINGFISH Sample

Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies.

preprint2019arXiv

Planck 2018 results. I. Overview and the cosmological legacy of Planck

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

preprint2019arXiv

Planck 2018 results. XII. Galactic astrophysics using polarized dust emission

We present 353 GHz full-sky maps of the polarization fraction $p$, angle $ψ$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $ψ$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas.

preprint2016arXiv

Planck 2015 results. XI. CMB power spectra, likelihoods, and robustness of parameters

This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlations of CMB data, using the hybrid approach employed previously: pixel-based at $\ell<30$ and a Gaussian approximation to the distribution of spectra at higher $\ell$. The main improvements are the use of more and better processed data and of Planck polarization data, and more detailed foreground and instrumental models, allowing further checks and enhanced immunity to systematics. Progress in foreground modelling enables a larger sky fraction. Improvements in processing and instrumental models further reduce uncertainties. For temperature, we perform an analysis of end-to-end instrumental simulations fed into the data processing pipeline; this does not reveal biases from residual instrumental systematics. The $Λ$CDM cosmological model continues to offer a very good fit to Planck data. The slope of primordial scalar fluctuations, $n_s$, is confirmed smaller than unity at more than 5σ from Planck alone. We further validate robustness against specific extensions to the baseline cosmology. E.g., the effective number of neutrino species remains compatible with the canonical value of 3.046. This first detailed analysis of Planck polarization concentrates on E modes. At low $\ell$ we use temperature at all frequencies and a subset of polarization. The frequency range improves CMB-foreground separation. Within the baseline model this requires a reionization optical depth $τ=0.078\pm0.019$, significantly lower than without high-frequency data for explicit dust monitoring. At high $\ell$ we detect residual errors in E, typically O($μ$K$^2$); we recommend temperature alone as the high-$\ell$ baseline. Nevertheless, Planck high-$\ell$ polarization allows a separate determination of $Λ$CDM parameters consistent with those from temperature alone.

preprint2016arXiv

Planck 2015 results. XIII. Cosmological parameters

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

preprint2016arXiv

Planck 2015 results. XIV. Dark energy and modified gravity

We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the cosmological constant scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of state, principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When estimating the density of DE at early times, we significantly improve present constraints. We then move to general parameterizations of the DE or MG perturbations that encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use baryonic acoustic oscillations, type-Ia supernovae and local measurements of the Hubble constant. We further show the impact of measurements of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. All results that include only background parameterizations are in agreement with LCDM. When testing models that also change perturbations (even when the background is fixed to LCDM), some tensions appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck TT+lowP when parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 sigma when external data sets are included. It however disappears when including CMB lensing.

preprint2016arXiv

Planck 2015 results. XIX. Constraints on primordial magnetic fields

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfvén waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

preprint2016arXiv

Planck 2015 results. XVII. Constraints on primordial non-Gaussianity

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

preprint2016arXiv

Planck 2015 results. XXV. Diffuse low-frequency Galactic foregrounds

(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...

preprint2016arXiv

Planck 2015 results. XXVI. The Second Planck Catalogue of Compact Sources

The Second Planck Catalogue of Compact Sources is a catalogue of sources detected in single-frequency maps from the full duration of the Planck mission and supersedes previous versions of the Planck compact source catalogues. It consists of compact sources, both Galactic and extragalactic, detected over the entire sky. Compact sources detected in the lower frequency channels are assigned to the PCCS2, while at higher frequencies they are assigned to one of two sub-catalogues, the PCCS2 or PCCS2E, depending on their location on the sky. The first of these catalogues covers most of the sky and allows the user to produce subsamples at higher reliabilities than the target 80% integral reliability of the catalogue. The PCCS2E contains sources detected in sky regions where the diffuse emission makes it difficult to quantify the reliability of the detections. Both the PCCS2 and PCCS2E include polarization measurements, in the form of polarized flux densities, or upper limits, and orientation angles for all seven polarization-sensitive Planck channels. The improved data-processing of the full-mission maps and their reduced noise levels allow us to increase the number of objects in the catalogue, improving its completeness for the target 80 % reliability as compared with the previous versions, the PCCS and ERCSC catalogues.

preprint2016arXiv

Planck intermediate results. XLIV. The structure of the Galactic magnetic field from dust polarization maps of the southern Galactic cap

We study the statistical properties of interstellar dust polarization at high Galactic latitude, using the Stokes parameter Planck maps at 353 GHz. Our aim is to advance the understanding of the magnetized interstellar medium (ISM), and to provide a model of the polarized dust foreground for cosmic microwave background component-separation procedures. Focusing on the southern Galactic cap, we examine the statistical distributions of the polarization fraction ($p$) and angle ($ψ$) to characterize the ordered and turbulent components of the Galactic magnetic field (GMF) in the solar neighbourhood. We relate patterns at large angular scales in polarization to the orientation of the mean (ordered) GMF towards Galactic coordinates $(l_0,b_0)=(70^\circ \pm 5^\circ,24^\circ \pm 5^\circ)$. The histogram of $p$ shows a wide dispersion up to 25 %. The histogram of $ψ$ has a standard deviation of $12^\circ$ about the regular pattern expected from the ordered GMF. We use these histograms to build a phenomenological model of the turbulent component of the GMF, assuming a uniform effective polarization fraction ($p_0$) of dust emission. To model the Stokes parameters, we approximate the integration along the line of sight (LOS) as a sum over a set of $N$ independent polarization layers, in each of which the turbulent component of the GMF is obtained from Gaussian realizations of a power-law power spectrum. We are able to reproduce the observed $p$ and $ψ$ distributions using: a $p_0$ value of (26 $\pm$ 3)%; a ratio of 0.9 $\pm$ 0.1 between the strengths of the turbulent and mean components of the GMF; and a small value of $N$. We relate the polarization layers to the density structure and to the correlation length of the GMF along the LOS.

preprint2016arXiv

Planck intermediate results. XLVI. Reduction of large-scale systematic effects in HFI polarization maps and estimation of the reionization optical depth

This paper describes the identification, modelling, and removal of previously unexplained systematic effects in the polarization data of the Planck High Frequency Instrument (HFI) on large angular scales, including new mapmaking and calibration procedures, new and more complete end-to-end simulations, and a set of robust internal consistency checks on the resulting maps. These maps, at 100, 143, 217, and 353 GHz, are early versions of those that will be released in final form later in 2016. The improvements allow us to determine the cosmic reionization optical depth $τ$ using, for the first time, the low-multipole $EE$ data from HFI, reducing significantly the central value and uncertainty, and hence the upper limit. Two different likelihood procedures are used to constrain $τ$ from two estimators of the CMB $E$- and $B$-mode angular power spectra at 100 and 143 GHz, after debiasing the spectra from a small remaining systematic contamination. These all give fully consistent results. A further consistency test is performed using cross-correlations derived from the Low Frequency Instrument maps of the Planck 2015 data release and the new HFI data. For this purpose, end-to-end analyses of systematic effects from the two instruments are used to demonstrate the near independence of their dominant systematic error residuals. The tightest result comes from the HFI-based $τ$ posterior distribution using the maximum likelihood power spectrum estimator from $EE$ data only, giving a value $0.055\pm 0.009$. In a companion paper these results are discussed in the context of the best-fit Planck $Λ$CDM cosmological model and recent models of reionization.

preprint2016arXiv

Planck intermediate results. XLVII. Planck constraints on reionization history

We investigate constraints on cosmic reionization extracted from the Planck cosmic microwave background (CMB) data. We combine the Planck CMB anisotropy data in temperature with the low-multipole polarization data to fit LCDM models with various parameterizations of the reionization history. We obtain a Thomson optical depth tau=0.058 +/- 0.012 for the commonly adopted instantaneous reionization model. This confirms, with only data from CMB anisotropies, the low value suggested by combining Planck 2015 results with other data sets and also reduces the uncertainties. We reconstruct the history of the ionization fraction using either a symmetric or an asymmetric model for the transition between the neutral and ionized phases. To determine better constraints on the duration of the reionization process, we also make use of measurements of the amplitude of the kinetic Sunyaev-Zeldovich (kSZ) effect using additional information from the high resolution Atacama Cosmology Telescope and South Pole Telescope experiments. The average redshift at which reionization occurs is found to lie between z=7.8 and 8.8, depending on the model of reionization adopted. Using kSZ constraints and a redshift-symmetric reionization model, we find an upper limit to the width of the reionization period of Dz < 2.8. In all cases, we find that the Universe is ionized at less than the 10% level at redshifts above z~10. This suggests that an early onset of reionization is strongly disfavoured by the Planck data. We show that this result also reduces the tension between CMB-based analyses and constraints from other astrophysical sources.

preprint2016arXiv

Planck intermediate results. XLVIII. Disentangling Galactic dust emission and cosmic infrared background anisotropies

Using the Planck 2015 data release (PR2) temperature maps, we separate Galactic thermal dust emission from cosmic infrared background (CIB) anisotropies. For this purpose, we implement a specifically tailored component-separation method, the so-called generalized needlet internal linear combination (GNILC) method, which uses spatial information (the angular power spectra) to disentangle the Galactic dust emission and CIB anisotropies. We produce significantly improved all-sky maps of Planck thermal dust emission, with reduced CIB contamination, at 353, 545, and 857 GHz. By reducing the CIB contamination of the thermal dust maps, we provide more accurate estimates of the local dust temperature and dust spectral index over the sky with reduced dispersion, especially at high Galactic latitudes above $b = \pm 20°$. We find that the dust temperature is $T = (19.4 \pm 1.3)$ K and the dust spectral index is $β= 1.6 \pm 0.1$ averaged over the whole sky, while $T = (19.4 \pm 1.5)$ K and $β= 1.6 \pm 0.2$ on 21 % of the sky at high latitudes. Moreover, subtracting the new CIB-removed thermal dust maps from the CMB-removed Planck maps gives access to the CIB anisotropies over 60 % of the sky at Galactic latitudes $|b| > 20°$. Because they are a significant improvement over previous Planck products, the GNILC maps are recommended for thermal dust science. The new CIB maps can be regarded as indirect tracers of the dark matter and they are recommended for exploring cross-correlations with lensing and large-scale structure optical surveys. The reconstructed GNILC thermal dust and CIB maps are delivered as Planck products.

preprint2015arXiv

A Joint Analysis of BICEP2/Keck Array and Planck Data

We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5°$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $μ$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies $\geq 150$ GHz to a lensed-$Λ$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$), using a prior on the frequency spectral behavior of polarized dust emission from previous \planck\ analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,σ$ significance.

preprint2015arXiv

Planck 2013 results. XXVII. Doppler boosting of the CMB: Eppur si muove

Our velocity relative to the rest frame of the cosmic microwave background (CMB) generates a dipole temperature anisotropy on the sky which has been well measured for more than 30 years, and has an accepted amplitude of v/c = 0.00123, or v = 369km/s. In addition to this signal generated by Doppler boosting of the CMB monopole, our motion also modulates and aberrates the CMB temperature fluctuations (as well as every other source of radiation at cosmological distances). This is an order 0.1% effect applied to fluctuations which are already one part in roughly one hundred thousand, so it is quite small. Nevertheless, it becomes detectable with the all-sky coverage, high angular resolution, and low noise levels of the Planck satellite. Here we report a first measurement of this velocity signature using the aberration and modulation effects on the CMB temperature anisotropies, finding a component in the known dipole direction, (l,b)=(264, 48) [deg], of 384km/s +- 78km/s (stat.) +- 115km/s (syst.). This is a significant confirmation of the expected velocity.

preprint2015arXiv

Planck 2015 results. I. Overview of products and scientific results

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

preprint2015arXiv

Planck 2015 results. IX. Diffuse component separation: CMB maps

We present foreground-reduced CMB maps derived from the full Planck data set in both temperature and polarization. Compared to the corresponding Planck 2013 temperature sky maps, the total data volume is larger by a factor of 3.2 for frequencies between 30 and 70 GHz, and by 1.9 for frequencies between 100 and 857 GHz. In addition, systematic errors in the forms of temperature-to-polarization leakage, analogue-to-digital conversion uncertainties, and very long time constant errors have been dramatically reduced, to the extent that the cosmological polarization signal may now be robustly recovered on angular scales $\ell\gtrsim40$. On the very largest scales, instrumental systematic residuals are still non-negligible compared to the expected cosmological signal, and modes with $\ell < 20$ are accordingly suppressed in the current polarization maps by high-pass filtering. As in 2013, four different CMB component separation algorithms are applied to these observations, providing a measure of stability with respect to algorithmic and modelling choices. The resulting polarization maps have rms instrumental noise ranging between 0.21 and 0.27$\,μ\textrm{K}$ averaged over 55 arcmin pixels, and between 4.5 and 6.1$\,μ\textrm{K}$ averaged over 3.4 arcmin pixels. The cosmological parameters derived from the analysis of temperature power spectra are in agreement at the $1σ$ level with the Planck 2015 likelihood. Unresolved mismatches between the noise properties of the data and simulations prevent a satisfactory description of the higher-order statistical properties of the polarization maps. Thus, the primary applications of these polarization maps are those that do not require massive simulations for accurate estimation of uncertainties, for instance estimation of cross-spectra and cross-correlations, or stacking analyses.

preprint2015arXiv

Planck 2015 results. X. Diffuse component separation: Foreground maps

Planck has mapped the microwave sky in nine frequency bands between 30 and 857 GHz in temperature and seven bands between 30 and 353 GHz in polarization. In this paper we consider the problem of diffuse astrophysical component separation, and process these maps within a Bayesian framework to derive a consistent set of full-sky astrophysical component maps. For the temperature analysis, we combine the Planck observations with the 9-year WMAP sky maps and the Haslam et al. 408 MHz map to derive a joint model of CMB, synchrotron, free-free, spinning dust, CO, line emission in the 94 and 100 GHz channels, and thermal dust emission. Full-sky maps are provided with angular resolutions varying between 7.5 arcmin and 1 deg. Global parameters (monopoles, dipoles, relative calibration, and bandpass errors) are fitted jointly with the sky model, and best-fit values are tabulated. For polarization, the model includes CMB, synchrotron, and thermal dust emission. These models provide excellent fits to the observed data, with rms temperature residuals smaller than 4 uK over 93% of the sky for all Planck frequencies up to 353 GHz, and fractional errors smaller than 1% in the remaining 7% of the sky. The main limitations of the temperature model at the lower frequencies are degeneracies among the spinning dust, free-free, and synchrotron components; additional observations from external low-frequency experiments will be essential to break these. The main limitations of the temperature model at the higher frequencies are uncertainties in the 545 and 857 GHz calibration and zero-points. For polarization, the main outstanding issues are instrumental systematics in the 100-353 GHz bands on large angular scales in the form of temperature-to-polarization leakage, uncertainties in the analog-to-digital conversion, and very long time constant corrections, all of which are expected to improve in the near future.

preprint2015arXiv

Planck 2015 results. XXIII. The thermal Sunyaev-Zeldovich effect--cosmic infrared background correlation

We use Planck data to detect the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) effect and the infrared emission from the galaxies that make up the the cosmic infrared background (CIB). We first perform a stacking analysis towards Planck-confirmed galaxy clusters. We detect infrared emission produced by dusty galaxies inside these clusters and demonstrate that the infrared emission is about 50% more extended than the tSZ effect. Modelling the emission with a Navarro--Frenk--White profile, we find that the radial profile concentration parameter is $c_{500} = 1.00^{+0.18}_{-0.15}$. This indicates that infrared galaxies in the outskirts of clusters have higher infrared flux than cluster-core galaxies. We also study the cross-correlation between tSZ and CIB anisotropies, following three alternative approaches based on power spectrum analyses: (i) using a catalogue of confirmed clusters detected in Planck data; (ii) using an all-sky tSZ map built from Planck frequency maps; and (iii) using cross-spectra between Planck frequency maps. With the three different methods, we detect the tSZ-CIB cross-power spectrum at significance levels of (i) 6 $σ$, (ii) 3 $σ$, and (iii) 4 $σ$. We model the tSZ-CIB cross-correlation signature and compare predictions with the measurements. The amplitude of the cross-correlation relative to the fiducial model is $A_{\rm tSZ-CIB}= 1.2\pm0.3$. This result is consistent with predictions for the tSZ-CIB cross-correlation assuming the best-fit cosmological model from Planck 2015 results along with the tSZ and CIB scaling relations.

preprint2015arXiv

Planck 2015 results. XXVIII. The Planck Catalogue of Galactic Cold Clumps

We present the Planck Catalogue of Galactic Cold Clumps (PGCC), an all-sky catalogue of Galactic cold clump candidates detected by Planck. This catalogue is the full version of the Early Cold Core (ECC) catalogue, which was made available in 2011 with the Early Release Compact Source Catalogue (ERCSC) and contained 915 high S/N sources. It is based on the Planck 48 months mission data that are currently being released to the astronomical community. The PGCC catalogue is an observational catalogue consisting exclusively of Galactic cold sources. The three highest Planck bands (857, 545, 353 GHz) have been combined with IRAS data at 3 THz to perform a multi-frequency detection of sources colder than their local environment. After rejection of possible extragalactic contaminants, the PGCC catalogue contains 13188 Galactic sources spread across the whole sky, i.e., from the Galactic plane to high latitudes, following the spatial distribution of the main molecular cloud complexes. The median temperature of PGCC sources lies between 13 and 14.5 K, depending on the quality of the flux density measurements, with a temperature ranging from 5.8 to 20 K after removing sources with the 1% largest temperature estimates. Using seven independent methods, reliable distance estimates have been obtained for 5574 sources, which allows us to derive their physical properties such as their mass, physical size, mean density and luminosity. The PGCC sources are located mainly in the solar neighbourhood, up to a distance of 10.5 kpc towards the Galactic centre, and range from low-mass cores to large molecular clouds. Because of this diversity and because the PGCC catalogue contains sources in very different environments, the catalogue is useful to investigate the evolution from molecular clouds to cores. Finally, the catalogue also includes 54 additional sources located in the SMC and LMC.

preprint2015arXiv

Planck intermediate results. XL. The Sunyaev-Zeldovich signal from the Virgo cluster

The Virgo cluster is the largest Sunyaev-Zeldovich (SZ) source in the sky, both in terms of angular size and total integrated flux. Planck's wide angular scale and frequency coverage, together with its high sensitivity, allow a detailed study of this large object through the SZ effect. Virgo is well resolved by Planck, showing an elongated structure, which correlates well with the morphology observed from X-rays, but extends beyond the observed X-ray signal. We find a good agreement between the SZ signal (or Compton paranmeter, y_c) observed by Planck and the expected signal inferred from X-ray observations and simple analytical models. Due to its proximity to us, the gas beyond the virial radius can be studied with unprecedented sensitivity by integrating the SZ signal over tens of square degrees. We study the signal in the outskirts of Virgo and compare it with analytical models and a constrained simulation of the environment of Virgo. Planck data suggest that significant amounts of low-density plasma surround Virgo out to twice the virial radius. We find the SZ signal in the outskirts of Virgo to be consistent with a simple model that extrapolates the inferred pressure at lower radii while assuming that the temperature stays in the keV range beyond the virial radius. The observed signal is also consistent with simulations and points to a shallow pressure profile in the outskirts of the cluster. This reservoir of gas at large radii can be linked with the hottest phase of the elusive warm/hot intergalactic medium. Taking the lack of symmetry of Virgo into account, we find that a prolate model is favoured by the combination of SZ and X-ray data, in agreement with predictions.

preprint2015arXiv

Planck intermediate results. XXI. Comparison of polarized thermal emission from Galactic dust at 353 GHz with interstellar polarization in the visible

The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust, bringing new constraints on the properties of dust. The dust grains that emit the radiation seen by Planck in the submillimetre also extinguish and polarize starlight in the visible. Comparison of the polarization of the emission and of the interstellar polarization on selected lines of sight probed by stars provides unique new diagnostics of the emission and light scattering properties of dust. Using ancillary catalogues of interstellar polarization and extinction of starlight, we obtain the degree of polarization, $p_V$ , and the optical depth in the V band to the star, $τ_V$. Toward these stars we measure the submillimetre polarized intensity, $P_S$, and total intensity, $I_S$, in the Planck 353 GHz channel. For those lines of sight through the diffuse interstellar medium with comparable values of the estimated column density and polarization directions close to orthogonal, we correlate properties in the submillimetre and visible to find two ratios, $R_{S/V} = (P_S/I_S)/(p_V/τ_V)$ and $R_{P/p} = P_S/p_V$ , the latter focusing directly on the polarization properties of the aligned grain population alone. We find $R_{S/V}$ = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and $R_{P/p}$ = 5.4 MJy sr$^{-1}$, with uncertainties 0.2 and 0.3 MJy sr$^{-1}$, respectively. Our estimate of $R_{S/V}$ is compatible with predictions based on a range of polarizing dust models that have been developed for the diffuse interstellar medium. However, our estimate of $R_{P/p}$ is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, $R_{P/p}$, indicates that changes to the optical properties in the models of the aligned grain population are required.

preprint2015arXiv

Planck intermediate results. XXII. Frequency dependence of thermal emission from Galactic dust in intensity and polarization

Planck has mapped the intensity and polarization of the sky at microwave frequencies with unprecedented sensitivity. We make use of the Planck 353 GHz I, Q, and U Stokes maps as dust templates, and cross-correlate them with the Planck and WMAP data at 12 frequencies from 23 to 353 GHz, over circular patches with 10 degree radius. The cross-correlation analysis is performed for both intensity and polarization data in a consistent manner. We use a mask that focuses our analysis on the diffuse interstellar medium at intermediate Galactic latitudes. We determine the spectral indices of dust emission in intensity and polarization between 100 and 353 GHz, for each sky-patch. The mean values, $1.59\pm0.02$ for polarization and $1.51\pm0.01$ for intensity, for a mean dust temperature of 19.6 K, are close, but significantly different ($3.6\,σ$). We determine the mean spectral energy distribution (SED) of the microwave emission, correlated with the 353 GHz dust templates, by averaging the results of the correlation over all sky-patches. We find that the mean SED increases for decreasing frequencies at $ν< 60$ GHz, for both intensity and polarization. The rise of the polarization SED towards low frequencies may be accounted for by a synchrotron component correlated with dust, with no need for any polarization of the anomalous microwave emission. We use a spectral model to separate the synchrotron and dust polarization and to characterize the spectral dependence of the dust polarization fraction. The polarization fraction ($p$) of the dust emission decreases by $(21\pm6)$ % from 353 to 70 GHz. The decrease of $p$ could indicate differences in polarization efficiency among components of interstellar dust (e.g., carbon versus silicate grains). Our observational results provide inputs to quantify and optimize the separation between Galactic and cosmological polarization.

preprint2015arXiv

Planck intermediate results. XXVII. High-redshift infrared galaxy overdensity candidates and lensed sources discovered by Planck and confirmed by Herschel-SPIRE

[Abridged] We use the Planck all-sky submm and mm maps to search for rare sources distinguished by extreme brightness, a few hundreds of mJy, and their potential for being situated at high redshift. These "cold" Planck sources, selected using the High Frequency Instrument (HFI) directly from the maps and from the Planck Catalogue of Compact Sources (PCCS), all satisfy the criterion of having their rest-frame far-infrared peak redshifted to the frequency range 353 and 857 GHz. This colour-selection favours galaxies in the redshift range z=2-4, which we consider as cold peaks in the cosmic infrared background (CIB). We perform a dedicated Herschel-SPIRE follow-up of 234 such Planck targets, finding a significant excess of red 350 and 500um sources, in comparison to reference SPIRE fields. About 94% of the SPIRE sources in the Planck fields are consistent with being overdensities of galaxies peaking at 350um. About 3% are candidate lensed systems, all 12 of which have secure spectroscopic confirmations, placing them at redshifts z>2.2. The galaxy overdensities are detected with high significance, half of the sample showing statistical significance above 10sigma. The SPIRE photometric redshifts of galaxies in overdensities suggest a peak at z~2. Under the Td=35K assumption, we derive an infrared (IR) luminosity for each SPIRE source of about 4x10^12 Lsun, yielding star formation rates of typically 700 Msun.yr^-1. If the observed overdensities are actual gravitationally-bound structures, the total total star formation rates reaches 7x10^3 Msun.yr^-1. Taken together, these sources show the signatures of high-z (z>$) protoclusters of intensively star-forming galaxies. All these observations confirm the uniqueness of our sample and demonstrate the ability of the all-sky Planck-HFI cold sources to select populations of cosmological and astrophysical interest for structure formation studies.

preprint2015arXiv

Planck intermediate results. XXXIX. The Planck list of high-redshift source candidates

The Planck mission, thanks to its large frequency range and all-sky coverage, has a unique potential for systematically detecting the brightest, and rarest, submillimetre sources on the sky, including distant objects in the high-redshift Universe traced by their dust emission. A novel method, based on a component-separation procedure using a combination of Planck and IRAS data, has been applied to select the most luminous cold submm sources with spectral energy distributions peaking between 353 and 857GHz at 5' resolution. A total of 2151 Planck high-z source candidates (the PHZ) have been detected in the cleanest 26% of the sky, with flux density at 545GHz above 500mJy. Embedded in the cosmic infrared background close to the confusion limit, these high-z candidates exhibit colder colours than their surroundings, consistent with redshifts z>2, assuming a dust temperature of 35K and a spectral index of 1.5. First follow-up observations obtained from optical to submm have confirmed that this list consists of two distinct populations. A small fraction (around 3%) of the sources have been identified as strongly gravitationally lensed star-forming galaxies, which are amongst the brightest submm lensed objects (with flux density at 545GHz ranging from 350mJy up to 1Jy) at redshift 2 to 4. However, the vast majority of the PHZ sources appear as overdensities of dusty star-forming galaxies, having colours consistent with z>2, and may be considered as proto-cluster candidates. The PHZ provides an original sample, complementary to the Planck Sunyaev-Zeldovich Catalogue; by extending the population of the virialized massive galaxy clusters to a population of sources at z>1.5, the PHZ may contain the progenitors of today's clusters. Hence the PHZ opens a new window on the study of the early ages of structure formation, and the understanding of the intensively star-forming phase at high-z.

preprint2015arXiv

Planck intermediate results. XXXV. Probing the role of the magnetic field in the formation of structure in molecular clouds

Within ten nearby (d < 450 pc) Gould Belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky, inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz, and the gas column density structures, quantified by the gradient of the column density, $N_H$. The selected regions, covering several degrees in size, are analyzed at an effective angular resolution of 10' FWHM, thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from $N_H \approx 10^{21}$ to $10^{23}$ cm$^{-2}$, and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analyzed in bins of column density using the novel statistical tool called "Histogram of Relative Orientations". Throughout this study, we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region, i.e., we assume a constant dust grain alignment efficiency, independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing $N_H$, from preferentially parallel or having no preferred orientation to preferentially perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfvénic or sub-Alfvénic turbulence, implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck. We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.

preprint2015arXiv

Planck intermediate results. XXXVIII. E- and B-modes of dust polarization from the magnetized filamentary structure of the interstellar medium

The quest for a B-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the 353 GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We filter the intensity and polarization maps to isolate filaments in the range of angular scales where the power asymmetry between E-modes and B-modes is observed. Using the Smoothed Hessian Major Axis Filament Finder, we identify 259 filaments at high Galactic latitude, with lengths larger or equal to 2° (corresponding to 3.5 pc in length for a typical distance of 100 pc). These filaments show a preferred orientation parallel to the magnetic field projected onto the plane of the sky, derived from their polarization angles. We present mean maps of the filaments in Stokes I, Q, U, E, and B, computed by stacking individual images rotated to align the orientations of the filaments. Combining the stacked images and the histogram of relative orientations, we estimate the mean polarization fraction of the filaments to be 11 %. Furthermore, we show that the correlation between the filaments and the magnetic field orientations may account for the E and B asymmetry and the $C_{\ell}^{TE}/C_{\ell}^{EE}$ ratio, reported in the power spectra analysis of the Planck 353 GHz polarization maps. Future models of the dust foreground for CMB polarization studies will need to take into account the observed correlation between the dust polarization and the structure of interstellar matter.

preprint2015arXiv

The Palomar Transient Factory: High Quality Realtime Data Processing in a Cost-Constrained Environment

The Palomar Transient Factory (PTF) is a synoptic sky survey in operation since 2009. PTF utilizes a 7.1 square degree camera on the Palomar 48-inch Schmidt telescope to survey the sky primarily at a single wavelength (R-band) at a rate of 1000-3000 square degrees a night. The data are used to detect and study transient and moving objects such as gamma ray bursts, supernovae and asteroids, as well as variable phenomena such as quasars and Galactic stars. The data processing system at IPAC handles realtime processing and detection of transients, solar system object processing, high photometric precision processing and light curve generation, and long-term archiving and curation. This was developed under an extremely limited budget profile in an unusually agile development environment. Here we discuss the mechanics of this system and our overall development approach. Although a significant scientific installation in of itself, PTF also serves as the prototype for our next generation project, the Zwicky Transient Facility (ZTF). Beginning operations in 2017, ZTF will feature a 50 square degree camera which will enable scanning of the entire northern visible sky every night. ZTF in turn will serve as a stepping stone to the Large Synoptic Survey Telescope (LSST), a major NSF facility scheduled to begin operations in the early 2020s.

preprint2015arXiv

The Type IIb SN 2011dh - Two years of observations and modelling of the lightcurves

We present optical and near-infrared (NIR) photometry and spectroscopy as well as modelling of the lightcurves of the Type IIb supernova (SN) 2011dh. Our extensive dataset spans two years, and complemented with Spitzer mid-infrared (MIR) data, we use it to build an optical-to-MIR bolometric lightcurve between days 3 and 732. To model the bolometric lightcurve before day 400 we use a grid of hydrodynamical SN models and a bolometric correction determined with steady-state NLTE modelling. Using this method we find a helium core mass of 3.1 (+0.7-0.4) solar masses for SN 2011dh, consistent within error bars with previous results. We present bolometric and broad-band lightcurves between days 100 and 500 for the Jerkstrand et al. (2014) steady-state NLTE models, and the preferred 12 solar masses (initial mass) model shows a good overall agreement with the observed lightcurves. We find an excess in the K and the MIR bands developing between days 100 and 250, during which an increase in the optical decline rate is also observed. A local origin of the excess is suggested by the depth of the HeI 2.058 micron absorption. Steady-state NLTE models with a modest dust opacity in the core, turned on during this period, reproduce the observed behaviour, but an additional excess in the Spitzer 4.5 micron band remains. CO first-overtone band emission is detected at day 206, and possibly at day 89, and assuming the additional excess to be dominated by CO fundamental band emission, we find fundamental to first-overtone band ratios considerably higher than observed in SN 1987A. This paper concludes our extensive observational and modelling work on SN 2011dh. The results from hydrodynamical modelling, steady-state NLTE modelling, and stellar evolutionary progenitor analysis are all consistent, and suggest an initial mass of 12 solar masses for the progenitor.

preprint2015arXiv

The Whole is Greater than the Sum of the Parts: Optimizing the Joint Science Return from LSST, Euclid and WFIRST

The focus of this report is on the opportunities enabled by the combination of LSST, Euclid and WFIRST, the optical surveys that will be an essential part of the next decade's astronomy. The sum of these surveys has the potential to be significantly greater than the contributions of the individual parts. As is detailed in this report, the combination of these surveys should give us multi-wavelength high-resolution images of galaxies and broadband data covering much of the stellar energy spectrum. These stellar and galactic data have the potential of yielding new insights into topics ranging from the formation history of the Milky Way to the mass of the neutrino. However, enabling the astronomy community to fully exploit this multi-instrument data set is a challenging technical task: for much of the science, we will need to combine the photometry across multiple wavelengths with varying spectral and spatial resolution. We identify some of the key science enabled by the combined surveys and the key technical challenges in achieving the synergies.

preprint2014arXiv

[CII] 158 $μ$m Emission as a Star Formation Tracer

The [CII] 157.74 $μ$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $Σ_{\rm [CII]}-Σ_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $Σ_{\rm[CII]}-Σ_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $Σ_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies.

preprint2014arXiv

Cool dust heating and temperature mixing in nearby star-forming galaxies

Physical conditions of the interstellar medium in galaxies are closely linked to the ambient radiation field and the heating of dust grains. In order to characterize dust properties in galaxies over a wide range of physical conditions, we present here the radial surface brightness profiles of the entire sample of 61 galaxies from Key Insights into Nearby Galaxies: Far-Infrared Survey with Herschel (KINGFISH). The main goal of our work is the characterization of the grain emissivities, dust temperatures, and interstellar radiation fields responsible for heating the dust. After fitting the dust and stellar radial profiles with exponential functions, we fit the far-infrared spectral energy distribution (SED) in each annular region with single-temperature modified black bodies using both variable (MBBV) and fixed (MBBF) emissivity indices beta, as well as with physically motivated dust models. Results show that while most SED parameters decrease with radius, the emissivity index beta also decreases with radius in some galaxies, but in others is increasing, or rising in the inner regions and falling in the outer ones. Despite the fixed grain emissivity (average beta~ 2.1) of the physically-motivated models, they are well able to accommodate flat spectral slopes with beta<= 1. We find that flatter slopes (beta<= 1.5) are associated with cooler temperatures, contrary to what would be expected from the usual Tdust-beta degeneracy. This trend is related to variations in Umin since beta and Umin are very closely linked over the entire range in Umin sampled by the KINGFISH galaxies: low Umin is associated with flat beta<=1. Both these results strongly suggest that the low apparent βvalues (flat slopes) in MBBV fits are caused by temperature mixing along the line-of-sight, rather than by intrinsic variations in grain properties. Abstract truncated for arXiv.

preprint2014arXiv

Planck 2013 results. I. Overview of products and scientific results

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

preprint2014arXiv

Planck 2013 results. VI. High Frequency Instrument data processing

We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (hereafter HFI), which we performed to produce six temperature maps from the first 473 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9.7 to 4.6 arcmin. The detector noise per (effective) beam solid angle is respectively, 10, 6, 12 and 39 microKelvin in HFI four lowest frequency channel (100--353 GHz) and 13 and 14 kJy/sr for the 545 and 857 GHz channels. Using the 143 GHz channel as a reference, these two high frequency channels are intercalibrated within 5% and the 353 GHz relative calibration is at the percent level. The 100 and 217 GHz channels, which together with the 143 GHz channel determine the high-multipole part of the CMB power spectrum (50 < l <2500), are intercalibrated at better than 0.2 %.

preprint2014arXiv

Planck 2013 results. VIII. HFI photometric calibration and mapmaking

This paper describes the processing applied to the HFI cleaned time-ordered data to produce photometrically calibrated maps. HFI observes the sky over a broad range of frequencies, from 100 to 857 GHz. To get the best accuracy on the calibration on such a large range, two different photometric calibration schemes have to be used. The 545 and 857 \GHz\ data are calibrated using Uranus and Neptune flux density measurements, compared with models of their atmospheric emissions to calibrate the data. The lower frequencies (below 353 GHz) are calibrated using the cosmological microwave background dipole.One of the components of this anisotropy results from the orbital motion of the satellite in the Solar System, and is therefore time-variable. Photometric calibration is thus tightly linked to mapmaking, which also addresses low frequency noise removal. The 2013 released HFI data show some evidence for apparent gain variations of the HFI bolometers' detection chain. These variations were identified by comparing observations taken more than one year apart in the same configuration. We developed an effective correction to limit its effect on calibration, and assess its accuracy. We present several methods used to estimate the precision of the photometric calibration. We distinguish relative (from one detector to another, or from one frequency to another) and absolute uncertainties. In both cases, we found that these uncertainties range from a few $10^{-3}$ to several per cents from 100 to 857 GHz. We describe the pipeline producing the maps from the HFI timelines, based on the photometric calibration parameters and we detail the scheme used to a posteriori set the zero level of the maps. We also briefly discuss the cross-calibration between HFI and the SPIRE instrument on board Herschel. We finally summarize the basic characteristics of the set of the HFI maps from the 2013 Planck data release.

preprint2014arXiv

Planck 2013 results. XI. All-sky model of thermal dust emission

This paper presents an all-sky model of dust emission from the Planck 857, 545 and 353 GHz, and IRAS 100 micron data. Using a modified black-body fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a tight representation of the data at 5 arc min. It shows variations of the order of 30 % compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions. An increase of the dust opacity at 353 GHz, tau_353/N_H, from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T_obs, that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at HI column densities lower than 10^20 cm^-2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anti-correlation between tau_353/N_H and T_obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant. The implication is that in the diffuse high-latitude ISM tau_353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of tau_353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B-V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data.

preprint2014arXiv

Planck 2013 results. XIV. Zodiacal emission

The Planck satellite provides a set of all-sky maps at nine frequencies from 30 GHz to 857 GHz. Planets, minor bodies, and diffuse interplanetary dust emission (IPD) are all observed. The IPD can be separated from Galactic and other emissions because Planck views a given point on the celestial sphere multiple times, through different columns of IPD. We use the Planck data to investigate the behaviour of zodiacal emission over the whole sky at sub-millimetre and millimetre wavelengths. We fit the Planck data to find the emissivities of the various components of the COBE zodiacal model -- a diffuse cloud, three asteroidal dust bands, a circumsolar ring, and an Earth-trailing feature. The emissivity of the diffuse cloud decreases with increasing wavelength, as expected from earlier analyses. The emissivities of the dust bands, however, decrease less rapidly, indicating that the properties of the grains in the bands are different from those in the diffuse cloud. We fit the small amount of Galactic emission seen through the telescope's far sidelobes, and place limits on possible contamination of the CMB results from both zodiacal and far-sidelobe emission. When necessary, the results are used in the Planck pipeline to make maps with zodiacal emission and far sidelobes removed. We show that the zodiacal correction to the CMB maps is small compared to the Planck CMB temperature power spectrum and give a list of flux densities for small Solar System bodies.

preprint2014arXiv

Planck 2013 results. XXIII. Isotropy and statistics of the CMB

The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3 sigma). However, we find little evidence for non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance of the CMB signal. A power asymmetry is now found to persist to scales corresponding to about l=600, and can be described in the low-l regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence for a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking.

preprint2014arXiv

Planck intermediate results. XIX. An overview of the polarized thermal emission from Galactic dust

This paper presents the large-scale polarized sky as seen by Planck HFI at 353 GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse large-scale maps of dust polarization fraction and polarization direction, while taking account of noise bias and possible systematic effects. We find that the maximum observed dust polarization fraction is high (pmax > 18%), in particular in some of the intermediate dust column density (AV < 1mag) regions. There is a systematic decrease in the dust polarization fraction with increasing dust column density, and we interpret the features of this correlation in light of both radiative grain alignment predictions and fluctuations in the magnetic field orientation. We also characterize the spatial structure of the polarization angle using the angle dispersion function and find that, in nearby fields at intermediate latitudes, the polarization angle is ordered over extended areas that are separated by filamentary structures, which appear as interfaces where the magnetic field sky projection rotates abruptly without apparent variations in the dust column density. The polarization fraction is found to be anti-correlated with the dispersion of the polarization angle, implying that the variations are likely due to fluctuations in the 3D magnetic field orientation along the line of sight sampling the diffuse interstellar medium.We also compare the dust emission with the polarized synchrotron emission measured with the Planck LFI, with low-frequency radio data, and with Faraday rotation measurements of extragalactic sources. The two polarized components are globally similar in structure along the plane and notably in the Fan and North Polar Spur regions. A detailed comparison of these three tracers shows, however, that dust and cosmic rays generally sample different parts of the line of sight and confirms that much of the variation observed in the Planck data is due to the 3D structure of the magnetic field.

preprint2014arXiv

Planck intermediate results. XVII. Emission of dust in the diffuse interstellar medium from the far-infrared to microwave frequencies

The dust-HI correlation is used to characterize the emission properties of dust in the diffuse interstellar medium. We cross-correlate sky maps from Planck, WMAP, and DIRBE, at 17 frequencies from 23 to 3000 GHz, with the Parkes survey of the 21-cm line emission of neutral atomic hydrogen, over a contiguous area of 7500 deg$^2$ centred on the southern Galactic pole. Our analysis yields four specific results. (1) The dust temperature is observed to be anti-correlated with the dust emissivity and opacity. We interpret this result as evidence for dust evolution within the diffuse ISM. The mean dust opacity is measured to be $(7.1 \pm 0.6) 10^{-27} cm^2/H \times (ν/353\, GHz)^{1.53\pm0.03}$ for $100 < ν<353$GHz. (2) We map the spectral index of dust emission at millimetre wavelengths, which is remarkably constant at $β_{mm} = 1.51\pm 0.13$. We compare it with the far infrared spectral index beta_FIR derived from greybody fits at higher frequencies, and find a systematic difference, $β_{mm}-β_{FIR} = -0.15$, which suggests that the dust SED flattens at $ν< 353\,$GHz. (3) We present spectral fits of the microwave emission correlated with HI from 23 to 353 GHz, which separate dust and anomalous microwave emission. The flattening of the dust SED can be accounted for with an additional component with a blackbody spectrum, which accounts for $(26 \pm 6)$% of the dust emission at 100 GHz and could represent magnetic dipole emission. Alternatively, it could account for an increasing contribution of carbon dust, or a flattening of the emissivity of amorphous silicates, at millimetre wavelengths. These interpretations make different predictions for the dust polarization SED. (4) We identify a Galactic contribution to the residuals of the dust-HI correlation, which we model with variations of the dust emissivity on angular scales smaller than that of our correlation analysis.

preprint2014arXiv

Planck intermediate results. XX. Comparison of polarized thermal emission from Galactic dust with simulations of MHD turbulence

Polarized emission observed by Planck HFI at 353 GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions $p$ and angles $ψ$. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction $p_\mathrm{max}$ decreases with column density $N_\mathrm{H}$ in the more opaque fields with $N_\mathrm{H} > 10^{21}\,\mathrm{cm}^{-2}$; and (iii) the polarization fraction along a given line of sight is correlated with the local spatial coherence of the polarization angle. These observations are compared to polarized emission maps computed in simulations of anisotropic magnetohydrodynamical (MHD) turbulence in which we assume a uniform intrinsic polarization fraction of the dust grains. We find that an estimate of this parameter may be recovered from the maximum polarization fraction $p_\mathrm{max}$ in diffuse regions where the magnetic field is ordered on large scales and perpendicular to the line of sight. This emphasizes the impact of anisotropies of the magnetic field on the emerging polarization signal. The decrease of the polarization fraction with column density in nearby molecular clouds is well reproduced in the simulations, indicating that it is essentially due to the turbulent structure of the magnetic field: an accumulation of variously polarized structures along the line of sight leads to such an anti-correlation. In the simulations, polarization fractions are also found to anti-correlate with the angle dispersion function $\mathcal{S}$. [abridged]

preprint2014arXiv

Planck intermediate results. XXX. The angular power spectrum of polarized dust emission at intermediate and high Galactic latitudes

The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_\ell^{EE,BB}$ over the range $40<\ell<600$ well away from the Galactic plane. These will bring new insights into interstellar dust physics and a precise determination of the level of contamination for CMB polarization experiments. We show that statistical properties of the emission can be characterized over large fractions of the sky using $C_\ell$. For the dust, they are well described by power laws in $\ell$ with exponents $α^{EE,BB}=-2.42\pm0.02$. The amplitudes of the polarization $C_\ell$ vary with the average brightness in a way similar to the intensity ones. The dust polarization frequency dependence is consistent with modified blackbody emission with $β_d=1.59$ and $T_d=19.6$K. We find a systematic ratio between the amplitudes of the Galactic $B$- and $E$-modes of 0.5. We show that even in the faintest dust-emitting regions there are no "clean" windows where primordial CMB $B$-mode polarization could be measured without subtraction of dust emission. Finally, we investigate the level of dust polarization in the BICEP2 experiment field. Extrapolation of the Planck 353GHz data to 150GHz gives a dust power $\ell(\ell+1)C_\ell^{BB}/(2π)$ of $1.32\times10^{-2}μ$K$_{CMB}^2$ over the $40<\ell<120$ range; the statistical uncertainty is $\pm0.29$ and there is an additional uncertainty (+0.28,-0.24) from the extrapolation, both in the same units. This is the same magnitude as reported by BICEP2 over this $\ell$ range, which highlights the need for assessment of the polarized dust signal even in the cleanest windows of the sky.

preprint2013arXiv

A Detailed Study of the Radio--FIR Correlation in NGC6946 with Herschel-PACS/SPIRE from KINGFISH

We derive the distribution of the synchrotron spectral index across NGC6946 and investigate the correlation between the radio continuum (synchrotron) and far-infrared (FIR) emission using the KINGFISH Herschel PACS and SPIRE data. The radio--FIR correlation is studied as a function of star formation rate, magnetic field strength, radiation field strength, and the total gas surface brightness. The synchrotron emission follows both star-forming regions and the so-called magnetic arms present in the inter-arm regions. The synchrotron spectral index is steepest along the magnetic arms ($α_n \sim 1$), while it is flat in places of giant H{\sc ii} regions and in the center of the galaxy ($α_n \sim 0.6-0.7$). The map of $α_n$ provides an observational evidence for aging and energy loss of cosmic ray electrons propagating in the disk of the galaxy. Variations in the synchrotron--FIR correlation across the galaxy are shown to be a function of both star formation and magnetic fields. We find that the synchrotron emission correlates better with cold rather than with warm dust emission, when the interstellar radiation field is the main heating source of dust. The synchrotron--FIR correlation suggests a coupling between the magnetic field and the gas density. NGC6946 shows a power-law behavior between the total (turbulent) magnetic field strength B and the star formation rate surface density $Σ_{\rm SFR}$ with an index of 0.14\,(0.16)$\pm$0.01. This indicates an efficient production of the turbulent magnetic field with the increasing gas turbulence expected in actively star forming regions. The scale-by-scale analysis of the synchrotron--FIR correlation indicates that the ISM affects the propagation of old/diffused cosmic ray electrons, resulting in a diffusion coefficient of $D_0=4.6\times 10^{28}$\,cm$^2$\,s$^{-1}$ for 2.2\,GeV CREs.

preprint2013arXiv

Calibration of the total infrared luminosity of nearby galaxies from Spitzer and Herschel bands

We present new empirical calibrations to estimate resolved and integrated total infrared luminosities from Spitzer and Herschel bands used as monochromatic or combined tracers. We base our calibrations on resolved elements of nearby galaxies (3 to 30 Mpc) observed with Herschel. We perform a resolved SED modelling of these objects using the Draine and Li (2007) dust models and investigate the influence of the addition of SPIRE measurements in the estimation of LTIR. We find that using data up to 250 um leads to local LTIR values consistent with those obtained with a complete coverage (up to 500 um) within 10 per cent for most of our resolved elements. We then study the distribution of energy in the resolved SEDs of our galaxies. The bulk of energy (30-50 per cent) is contained in the (70-160 um) band. The (24-70 um) fraction decreases with increasing metallicity. The (160-1100 um) submillimeter band can account for up to 25 per cent of the LTIR in metal-rich galaxies. We investigate the correlation between TIR surface brightnesses/luminosities and monochromatic Spitzer and Herschel surface brightnesses/luminosities. The three PACS bands can be used as reliable monochromatic estimators of the LTIR, the 100 um band being the most reliable monochromatic tracer. There is also a strong correlation between the SPIRE 250um and LTIR, although with more scatter than for the PACS relations. We also study the ability of our monochromatic relations to reproduce integrated LTIR of nearby galaxies as well as LTIR of z=1-3 sources. Finally, we provide calibration coefficients that can be used to derive TIR surface brightnesses/luminosities from a combination of Spitzer and Herschel surface brightnesses/fluxes and analyse the associated uncertainties.

preprint2013arXiv

Mid- to far infrared properties of star-forming galaxies and active galactic nuclei

We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative contribution of the AGN to the infrared energy budget of the sources. We found that the total infrared emission of galaxies with weak 6.2um PAH emission (EW<0.2um) is dominated by AGN activity, while for galaxies with EW>0.2um more than 50% of the LIR arises from star formation. We also found that for galaxies detected in the 250-500um Herschel bands an AGN has a statistically insignificant effect on the temperature of the cold dust and the far-IR colours of the host galaxy, which are primarily shaped by star formation activity. For star-forming galaxies we reveal an anti-correlation between the LIR-to-rest-frame 8um luminosity ratio, IR8 = LIR\L8, and the strength of PAH features. We found that this anti-correlation is primarily driven by variations in the PAHs emission, and not by variations in the 5-15um mid-IR continuum emission. Using the [NeIII]/[NeII] line ratio as a tracer of the hardness of the radiation field, we confirm that galaxies with harder radiation fields tend to exhibit weaker PAH features, and found that they have higher IR8 values and higher dust-mass-weighted luminosities (LIR/Md), the latter being a proxy for the dust temperature (Td). We argue that these trends originate either from variations in the environment of the star-forming regions or are caused by variations in the age of the starburst. Finally, we provide scaling relations that will allow estimating LIR, based on single-band observations with the mid-infrared instrument, on board the upcoming JWST.

preprint2013arXiv

Planck 2013 results. XII. Component separation

Planck has produced detailed all-sky observations over nine frequency bands between 30 and 857 GHz. These observations allow robust reconstruction of the primordial cosmic microwave background (CMB) temperature fluctuations over nearly the full sky, as well as new constraints on Galactic foregrounds. This paper describes the component separation framework adopted by Planck. We test four foreground-cleaned CMB maps derived using qualitatively different component separation algorithms. The quality of our reconstructions is evaluated through detailed simulations and internal comparisons, and shown through various tests to be internally consistent and robust for CMB power spectrum and cosmological parameter estimation up to l = 2000. The parameter constraints on LambdaCDM cosmologies derived from these maps are consistent with those presented in the cross-spectrum based Planck likelihood analysis. We choose two of the CMB maps for specific scientific goals. We also present maps and frequency spectra of the Galactic low-frequency, CO, and thermal dust emission. The component maps are found to provide a faithful representation of the sky, as evaluated by simulations. For the low-frequency component, the spectral index varies widely over the sky, ranging from about beta = -4 to -2. Considering both morphology and prior knowledge of the low frequency components, the index map allows us to associate a steep spectral index (beta < -3.2) with strong anomalous microwave emission, corresponding to a spinning dust spectrum peaking below 20 GHz, a flat index of beta > -2.3 with strong free-free emission, and intermediate values with synchrotron emission.

preprint2013arXiv

Planck 2013 results. XV. CMB power spectra and likelihood

We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the error budget at higher l's. For l < 50, our likelihood exploits all Planck frequency channels from 30 to 353 GHz through a physically motivated Bayesian component separation technique. At l >= 50, we employ a correlated Gaussian likelihood approximation based on angular cross-spectra derived from the 100, 143 and 217 GHz channels. We validate our likelihood through an extensive suite of consistency tests, and assess the impact of residual foreground and instrumental uncertainties on cosmological parameters. We find good internal agreement among the high-l cross-spectra with residuals of a few uK^2 at l <= 1000. We compare our results with foreground-cleaned CMB maps, and with cross-spectra derived from the 70 GHz Planck map, and find broad agreement in terms of spectrum residuals and cosmological parameters. The best-fit LCDM cosmology is in excellent agreement with preliminary Planck polarisation spectra. The standard LCDM cosmology is well constrained by Planck by l <= 1500. For example, we report a 5.4 sigma deviation from n_s /= 1. Considering various extensions beyond the standard model, we find no indication of significant departures from the LCDM framework. Finally, we report a tension between the best-fit LCDM model and the low-l spectrum in the form of a power deficit of 5-10% at l <~ 40, significant at 2.5-3 sigma. We do not elaborate further on its cosmological implications, but note that this is our most puzzling finding in an otherwise remarkably consistent dataset. (Abridged)

preprint2013arXiv

Planck 2013 results. XXX. Cosmic infrared background measurements and implications for star formation

We present new measurements of CIB anisotropies using Planck. Combining HFI data with IRAS, the angular auto- and cross frequency power spectrum is measured from 143 to 3000 GHz, and the auto-bispectrum from 217 to 545 GHz. The total areas used to compute the CIB power spectrum and bispectrum are about 2240 and 4400 deg^2, respectively. After careful removal of the contaminants, and a complete study of systematics, the CIB power spectrum and bispectrum are measured with unprecedented signal to noise ratio from angular multipoles ell~150 to 2500, and ell~130 to 1100, respectively. Two approaches are developed for modelling CIB power spectrum anisotropies. The first approach takes advantage of the unique measurements by Planck at large angular scales, and models only the linear part of the power spectrum, with a mean bias of dark matter halos hosting dusty galaxies at a given redshift weighted by their contribution to the emissivities. The second approach is based on a model that associates star-forming galaxies with dark matter halos and their subhalos, using a parametrized relation between the dust-processed infrared luminosity and (sub-)halo mass. The two approaches simultaneously fit all auto- and cross- power spectra very well. We find that the star formation history is well constrained up to z~2. However, at higher redshift, the accuracy of the star formation history measurement is strongly degraded by the uncertainty in the spectral energy distribution of CIB galaxies. We also find that CIB galaxies have warmer temperatures as redshift increases. The CIB bispectrum is steeper than that expected from the power spectrum, although well fitted by a power law; this gives some information about the contribution of massive halos to the CIB bispectrum.

preprint2013arXiv

Star Formation Rates in Resolved Galaxies: Calibrations with Near and Far Infrared Data for NGC5055 and NGC6946

We use the near--infrared Brγhydrogen recombination line as a reference star formation rate (SFR) indicator to test the validity and establish the calibration of the {\it Herschel} PACS 70 μm emission as a SFR tracer for sub--galactic regions in external galaxies. Brγoffers the double advantage of directly tracing ionizing photons and of being relatively insensitive to the effects of dust attenuation. For our first experiment, we use archival CFHT Brγand Ks images of two nearby galaxies: NGC\,5055 and NGC\,6946, which are also part of the {\it Herschel} program KINGFISH (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel). We use the extinction corrected Brγemission to derive the SFR(70) calibration for H{\sc ii} regions in these two galaxies. A comparison of the SFR(70) calibrations at different spatial scales, from 200 pc to the size of the whole galaxy, reveals that about 50% of the total 70μm emission is due to dust heated by stellar populations that are unrelated to the current star formation. We use a simple model to qualitatively relate the increase of the SFR(70) calibration coefficient with decreasing region size to the star formation timescale. We provide a calibration for an unbiased SFR indicator that combines the observed Hαwith the 70 μm emission, also for use in H{\sc ii} regions. We briefly analyze the PACS 100 and 160 μm maps and find that longer wavelengths are not as good SFR indicators as 70μm, in agreement with previous results. We find that the calibrations show about 50% difference between the two galaxies, possibly due to effects of inclination.

preprint2012arXiv

A Study of Heating and Cooling of the ISM in NGC 1097 with Herschel-PACS and Spitzer-IRS

NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared (FIR) [CII]158um and [OI]63um lines and mid-infrared H2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([CII]158um+[OI]63um)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7um PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [CII]158μm/PAH(5.5 - 14um) is found. PAHs in the ring are responsible for a factor of two more [CII]158um and [OI]63um emission per unit mass than PAHs in the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high intensity photodissociation regions (PDRs), in which case G0 ~ 10^2.3 and nH ~ 10^3.5 cm^-3 in the ring. For these values of G0 and nH PDR models cannot reproduce the observed H2 emission. Much of the the H2 emission in the starburst ring could come from warm regions in the diffuse ISM that are heated by turbulent dissipation or shocks.

preprint2012arXiv

Characterizing Cosmic-Ray Propagation in Massive Star-Forming Regions: The Case of 30 Doradus and the Large Magellanic Cloud

Using infrared, radio, and gamma-ray data,we investigate the propagation characteristics of cosmic-ray (CR) electrons and nuclei in the 30 Doradus (30\,Dor) star-forming region in the Large Magellanic Cloud (LMC) using a phenomenological model based on the radio-far-infrared correlation within galaxies. Employing a correlation analysis, we derive an average propagation length of \sim 100-140 pc for \sim 3 GeV CR electrons resident in 30 Dor from consideration of the radio and infrared data. Assuming that the observed gamma-ray emission towards 30 Dor is associated with the star-forming region, and applying the same methodology to the infrared and gamma-ray data, we estimate a \sim 20 GeV propagation length of 200-320 pc for the CR nuclei. This is approximately twice as large as for\sim 3 GeV CR electrons, corresponding to a spatial diffusion coefficient that is \sim 4 times higher, scaling as (R/GV)^δ with δ\approx 0.7-0.8 depending on the smearing kernel used in the correlation analysis. This value is in agreement with the results found by extending the correlation analysis to include \sim 70 GeV CR nuclei traced by the 3-10 GeV gamma-ray data (δ\approx 0.66+/-0.23). Using the mean age of the stellar populations in 30 Dor and the results from our correlation analysis, we estimate a diffusion coefficient D_{R} \approx 0.9-1.0 \times10^27 (R/GV)^{0.7} cm^2 s^-1. We compare the values of the CR electron propagation length and surface brightness for 30 Dor and the LMC as a whole with those of entire disk galaxies. We find that the trend of decreasing average CR propagation distance with increasing disk-averaged star formation activity holds for the LMC, and extends down to single star-forming regions, at least for the case of 30 Dor.

preprint2012arXiv

Detection of Powerful Mid-IR H2 Emission in the Bridge between the Taffy Galaxies

We report the detection of strong, resolved emission from warm H2 in the Taffy galaxies and bridge. Relative to the continuum and faint PAH emission, the H2 emission is the strongest in the connecting bridge, approaching L(H2)/L(PAH8μm) = 0.1 between the two galaxies, where the purely rotational lines of H2 dominate the mid-infrared spectrum in a way very reminiscent of the group-wide shock in the interacting group Stephan's Quintet. The surface brightness in the 0-0 S(0) and S(1) H2 lines in the bridge is more than twice that observed at the center of the Stephan's Quintet shock. We observe a warm H2 mass of 4.2 \times 108 M\odot in the bridge, but taking into account the unobserved bridge area, the total warm mass is likely to be twice this value. We use excitation diagrams to characterize the warm molecular gas, finding an average surface mass of 5 \times 106 M\odot kpc-2 and typical excitation temperatures of 150-175 K. H2 emission is also seen in the galaxy disks, although there the emission is more consistent with normal star forming galaxies. We investigate several possible heating mechanisms for the bridge gas, but favor the conversion of kinetic energy from the head-on collision via turbulence and shocks as the main heating source. Since the cooling time for the warm H2 is short (5000 yr), shocks must be permeating the molecular gas in bridge region in order to continue heating the H2.

preprint2012arXiv

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: NGC 628 and NGC 6946

We characterize the dust in NGC628 and NGC6946, two nearby spiral galaxies in the KINGFISH sample. With data from 3.6um to 500um, dust models are strongly constrained. Using the Draine & Li (2007) dust model, (amorphous silicate and carbonaceous grains), for each pixel in each galaxy we estimate (1) dust mass surface density, (2) dust mass fraction contributed by polycyclic aromatic hydrocarbons (PAH)s, (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in regions with high starlight intensity. We obtain maps for the dust properties, which trace the spiral structure of the galaxies. The dust models successfully reproduce the observed global and resolved spectral energy distributions (SEDs). The overall dust/H mass ratio is estimated to be 0.0082+/-0.0017 for NGC628, and 0.0063+/-0.0009 for NGC6946, consistent with what is expected for galaxies of near-solar metallicity. Our derived dust masses are larger (by up to a factor 3) than estimates based on single-temperature modified blackbody fits. We show that the SED fits are significantly improved if the starlight intensity distribution includes a (single intensity) "delta function" component. We find no evidence for significant masses of cold dust T<12K. Discrepancies between PACS and MIPS photometry in both low and high surface brightness areas result in large uncertainties when the modeling is done at PACS resolutions, in which case SPIRE, MIPS70 and MIPS160 data cannot be used. We recommend against attempting to model dust at the angular resolution of PACS.

preprint2012arXiv

Planck Early Results. II. The thermal performance of Planck

The performance of the Planck instruments in space is enabled by their low operating temperatures, 20K for LFI and 0.1K for HFI, achieved through a combination of passive radiative cooling and three active mechanical coolers. The scientific requirement for very broad frequency coverage led to two detector technologies with widely different temperature and cooling needs. Active coolers could satisfy these needs; a helium cryostat, as used by previous cryogenic space missions (IRAS, COBE, ISO, Spitzer, AKARI), could not. Radiative cooling is provided by three V-groove radiators and a large telescope baffle. The active coolers are a hydrogen sorption cooler (<20K), a 4He Joule-Thomson cooler (4.7K), and a 3He-4He dilution cooler (1.4K and 0.1K). The flight system was at ambient temperature at launch and cooled in space to operating conditions. The HFI bolometer plate reached 93mK on 3 July 2009, 50 days after launch. The solar panel always faces the Sun, shadowing the rest of Planck, andoperates at a mean temperature of 384K. At the other end of the spacecraft, the telescope baffle operates at 42.3K and the telescope primary mirror operates at 35.9K. The temperatures of key parts of the instruments are stabilized by both active and passive methods. Temperature fluctuations are driven by changes in the distance from the Sun, sorption cooler cycling and fluctuations in gas-liquid flow, and fluctuations in cosmic ray flux on the dilution and bolometer plates. These fluctuations do not compromise the science data.

preprint2012arXiv

Planck Intermediate Results. IX. Detection of the Galactic haze with Planck

Using precise full-sky observations from Planck, and applying several methods of component separation, we identify and characterize the emission from the Galactic "haze" at microwave wavelengths. The haze is a distinct component of diffuse Galactic emission, roughly centered on the Galactic centre, and extends to |b| ~35 deg in Galactic latitude and |l| ~15 deg in longitude. By combining the Planck data with observations from the WMAP we are able to determine the spectrum of this emission to high accuracy, unhindered by the large systematic biases present in previous analyses. The derived spectrum is consistent with power-law emission with a spectral index of -2.55 +/- 0.05, thus excluding free-free emission as the source and instead favouring hard-spectrum synchrotron radiation from an electron population with a spectrum (number density per energy) dN/dE ~ E^-2.1. At Galactic latitudes |b|<30 deg, the microwave haze morphology is consistent with that of the Fermi gamma-ray "haze" or "bubbles," indicating that we have a multi-wavelength view of a distinct component of our Galaxy. Given both the very hard spectrum and the extended nature of the emission, it is highly unlikely that the haze electrons result from supernova shocks in the Galactic disk. Instead, a new mechanism for cosmic-ray acceleration in the centre of our Galaxy is implied.

preprint2012arXiv

Resolving the Far-IR Line Deficit: Photoelectric Heating and Far-IR Line Cooling in NGC 1097 and NGC 4559

The physical state of interstellar gas and dust is dependent on the processes which heat and cool this medium. To probe heating and cooling of the ISM over a large range of infrared surface brightness, on sub-kiloparsec scales, we employ line maps of [C \ii] 158 $μ$m, [O \one] 63 $μ$m, and [N \ii] 122 $μ$m in NGC 1097 and NGC 4559, obtained with the PACS spectrometer onboard {\it Herschel}. We matched new observations to existing Spitzer-IRS data that trace the total emission of polycyclic aromatic hydrocarbons (PAHs). We confirm at small scales in these galaxies that the canonical measure of photoelectric heating efficiency, ([C \ii] + [O \one])/TIR, decreases as the far-infrared color, $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m), increases. In contrast, the ratio of far-infrared (far-IR) cooling to total PAH emission, ([C \ii] + [O \one])/PAH, is a near constant $\sim$6% over a wide range of far-infrared color, 0.5 \textless\ $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m) $\lesssim$ 0.95. In the warmest regions, where $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m) $\gtrsim$ 0.95, the ratio ([C \ii] + [O \one])/PAH drops rapidly to 4%. We derived representative values of the local UV radiation density, $G_0$, and the gas density, $n_H$, by comparing our observations to models of photodissociation regions. The ratio $G_0/n_H$, derived from fine-structure lines, is found to correlate with the mean dust-weighted starlight intensity, $<U>$ derived from models of the IR SED. Emission from regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small grains are susceptible to ionization effects. We note that there is a shift in the 7.7 / 11.3 $μ$m PAH ratio in regions that exhibit a deficit in ([C \ii] + [O \one])/PAH, suggesting that small grains are ionized in these environments.

preprint2012arXiv

The Palomar Transient Factory photometric catalog 1.0

We construct a photometrically calibrated catalog of non-variable sources from the Palomar Transient Factory (PTF) observations. The first version of this catalog presented here, the PTF photometric catalog 1.0, contains calibrated R_PTF-filter magnitudes for about 21 million sources brighter than magnitude 19, over an area of about 11233 deg^2. The magnitudes are provided in the PTF photometric system, and the color of a source is required in order to convert these magnitudes into other magnitude systems. We estimate that the magnitudes in this catalog have typical accuracy of about 0.02 mag with respect to magnitudes from the Sloan Digital Sky Survey. The median repeatability of our catalog's magnitudes for stars between 15 and 16 mag, is about 0.01 mag, and it is better than 0.03 mag for 95% of the sources in this magnitude range. The main goal of this catalog is to provide reference magnitudes for photometric calibration of visible light observations. Subsequent versions of this catalog, which will be published incrementally online, will be extended to a larger sky area and will also include g_PTF-filter magnitudes, as well as variability and proper motion information.

preprint2012arXiv

The Star Formation in Radio Survey: GBT 33 GHz Observations of Nearby Galaxy Nuclei and Extranuclear Star-Forming Regions

We present 33\,GHz photometry of 103 galaxy nuclei and extranuclear star-forming complexes taken with the Green Bank Telescope (GBT) as part of the Star Formation in Radio Survey (SFRS). Among the sources without evidence for an AGN, and also having lower frequency radio data, we find a median thermal fraction at 33GHz of ~76% with a dispersion of ~24%. For all sources resolved on scales <0.5kpc, the thermal fraction is even larger, being >90%. This suggests that the rest-frame 33GHz emission provides a sensitive measure of the ionizing photon rate from young star-forming regions, thus making it a robust star formation rate indicator. Taking the 33GHz star formation rates as a reference, we investigate other empirical calibrations relying on different combinations of warm 24μm dust, total infrared (IR; 8-1000μm), Hα line, and far-UV continuum emission. The recipes derived here generally agree with others found in the literature, albeit with a large dispersion that most likely stems from a combination of effects. Comparing the 33GHz to total IR flux ratios as a function of the radio spectral index, measured between 1.7 and 33GHz, we find that the ratio increases as the radio spectral index flattens which does not appear to be a distance effect. Consequently, the ratio of non-thermal to total IR emission appears relatively constant, suggesting only moderate variations in the cosmic-ray electron injection spectrum and ratio of synchrotron to total cooling processes among star-forming complexes. Assuming that this trend solely arises from an increase in the thermal fraction sets a maximum on the scatter of the non-thermal spectral indices among the star-forming regions of σ_α^{NT} < 0.13.

preprint2012arXiv

Ultraluminous Star-Forming Galaxies and Extremely Luminous Warm Molecular Hydrogen Emission at z=2.16 in the PKS 1138-26 Radio Galaxy Protocluster

A deep Spitzer Infrared Spectrograph map of the PKS 1138-26 galaxy protocluster reveals ultraluminous PAH emission from obscured star formation in three protocluster galaxies, including Halpha-emitter (HAE) 229, HAE 131, and the central Spiderweb Galaxy. Star formation rates of 500-1100 Msun/yr are estimated from the 7.7 micron PAH feature. At such prodigious formation rates, the galaxy stellar masses will double in 0.6-1.1 Gyr. We are viewing the peak epoch of star formation for these protocluster galaxies. However, it appears that extinction of Halpha is much greater (factor of 40) in the two ULIRG HAEs compared to the Spiderweb. This may be attributed to different spatial distributions of star formation--nuclear star formation in the HAEs versus extended star formation in accreting satellite galaxies in the Spiderweb. We find extremely luminous mid-IR rotational line emission from warm molecular hydrogen in the Spiderweb Galaxy, with L(H2 0-0 S(3))= 1.4E44 erg/s (3.7E10 Lsun), 20 times more luminous than any previously known H2 emission galaxy (MOHEG). Depending on temperature, this corresponds to a very large mass of >9E6-2E9 Msun of T>300 K molecular gas, heated by the PKS 1138-26 radio jet, acting to quench nuclear star formation. There is >8 times more warm H2 at these temperatures in the Spiderweb than what has been seen in low-redshift (z<0.2) radio galaxies, indicating that the Spiderweb may have a larger reservoir of molecular gas than more evolved radio galaxies. This is the highest redshift galaxy yet in which warm molecular hydrogen has been directly detected.

preprint2012arXiv

Wide-Field InfraRed Survey Telescope (WFIRST) Final Report

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.

preprint2011arXiv

Calibrating Extinction-Free Star Formation Rate Diagnostics with 33GHz Free-Free Emission in NGC6946

Abridged: Using free-free emission measured in the Ka-band (26-40GHz) for 10 star-forming regions in the nearby galaxy NGC6946, including its starbursting nucleus, we compare a number of SFR diagnostics that are typically considered to be unaffected by interstellar extinction: i.e., non-thermal radio (i.e., 1.4GHz), total infrared (IR; 8-1000um), and warm dust (i.e., 24um) emission, along with the hybrid (obscured + unobscured) indicators of Hα+24um and UV+IR. The 33GHz free-free emission is assumed to provide the most accurate measure of the current SFR. Among the extranuclear star-forming regions, the 24um, Hα+24um and UV+IR SFR calibrations are in good agreement with the 33GHz free-free SFRs. However, each of the SFR calibrations relying on some form of dust emission overestimate the nuclear SFR by a factor of ~2. This is more likely the result of excess dust heating through an accumulation of non-ionizing stars associated with an extended episode of star formation in the nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the non-thermal radio continuum yield values which only agree with the free-free SFRs for the nucleus, and underestimate the SFRs from the extranuclear star-forming regions by a factor of ~2. This result likely arises from the CR electrons decaying within the starburst region with negligible escape compared to the young extranuclear star-forming regions. Finally, we find that the SFRs estimated using the total 33GHz emission agree well with the free-free SFRs due to the large thermal fractions present at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33GHz observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift without the need of ancillary radio data to account for the non-thermal emission.

preprint2011arXiv

Herschel Far-Infrared and Sub-millimeter Photometry for the KINGFISH Sample of Nearby Galaxies

New far-infrared and sub-millimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially-integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500um emission shows evidence for a sub-millimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photo-dissociation regions is found to be (21+/-4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.

preprint2011arXiv

KINGFISH -- Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel: Survey Description and Image Atlas

The KINGFISH project (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) is an imaging and spectroscopic survey of 61 nearby (d < 30 Mpc) galaxies, chosen to cover a wide range of galaxy properties and local interstellar medium (ISM) environments found in the nearby Universe. Its broad goals are to characterize the ISM of present-day galaxies, the heating and cooling of their gaseous and dust components, and to better understand the physical processes linking star formation and the ISM. KINGFISH is a direct descendant of the Spitzer Infrared Nearby Galaxies Survey (SINGS), which produced complete Spitzer imaging and spectroscopic mapping and a comprehensive set of multi-wavelength ancillary observations for the sample. The Herschel imaging consists of complete maps for the galaxies at 70, 100, 160, 250, 350, and 500 microns. The spectal line imaging of the principal atomic ISM cooling lines ([OI]63um, [OIII]88um, [NII]122,205um, and [CII]158um) covers the subregions in the centers and disks that already have been mapped in the mid-infrared with Spitzer. The KINGFISH and SINGS multi-wavelength datasets combined provide panchromatic mapping of the galaxies sufficient to resolve individual star-forming regions, and tracing the important heating and cooling channels of the ISM, across a wide range of local extragalactic ISM environments. This paper summarizes the scientific strategy for KINGFISH, the properties of the galaxy sample, the observing strategy, and data processing and products. It also presents a combined Spitzer and Herschel image atlas for the KINGFISH galaxies, covering the wavelength range 3.6 -- 500 microns. All imaging and spectroscopy data products will be released to the Herschel user generated product archives.

preprint2011arXiv

Planck Early Results XVIII: The power spectrum of cosmic infrared background anisotropies

Using Planck maps of six regions of low Galactic dust emission with a total area of about 140 square degrees, we determine the angular power spectra of cosmic infrared background (CIB) anisotropies from multipole l = 200 to l = 2000 at 217, 353, 545 and 857 GHz. We use 21-cm observations of HI as a tracer of thermal dust emission to reduce the already low level of Galactic dust emission and use the 143 GHz Planck maps in these fields to clean out cosmic microwave background anisotropies. Both of these cleaning processes are necessary to avoid significant contamination of the CIB signal. We measure correlated CIB structure across frequencies. As expected, the correlation decreases with increasing frequency separation, because the contribution of high-redshift galaxies to CIB anisotropies increases with wavelengths. We find no significant difference between the frequency spectrum of the CIB anisotropies and the CIB mean, with Delta I/I=15% from 217 to 857 GHz. In terms of clustering properties, the Planck data alone rule out the linear scale- and redshift-independent bias model. Non-linear corrections are significant. Consequently, we develop an alternative model that couples a dusty galaxy, parametric evolution model with a simple halo-model approach. It provides an excellent fit to the measured anisotropy angular power spectra and suggests that a different halo occupation distribution is required at each frequency, which is consistent with our expectation that each frequency is dominated by contributions from different redshifts. In our best-fit model, half of the anisotropy power at l=2000 comes from redshifts z<0.8 at 857 GHz and z<1.5 at 545 GHz, while about 90% come from redshifts z>2 at 353 and 217 GHz, respectively.

preprint2011arXiv

Planck Early Results XXIV: Dust in the diffuse interstellar medium and the Galactic halo

This paper presents the first results of comparison of Planck along with IRAS data with Green Bank Telescope 21-cm observations in 14 fields covering more than 800 square degrees at high Galactic latitude. Galactic dust emission for fields with average HI column density lower than 2 x 10^20 cm^-2 is well correlated with 21-cm emission. The residual emission in these fields, once the HI-correlated emission is removed, is consistent with the expected statistical properties of the cosmic infrared background fluctuations. Fields with larger column densities show significant excess dust emission compared to the HI column density. Regions of excess lie in organized structures that suggest the presence of hydrogen in molecular form, though they are not always correlated with CO emission. Dust emission from intermediate-velocity clouds is detected with high significance. Its spectral properties are consistent with, compared to the local ISM values, significantly hotter dust (T~20 K), lower sub-millimeter dust opacity, and a relative abundance of very small grains to large grains about four times higher. These results are compatible with expectations for clouds that are part of the Galactic fountain in which there is dust shattering and fragmentation. Correlated dust emission in HVCs is not detected; the average of the 99.9% confidence upper limits to the emissivity is 0.15 times the local ISM value at 857 and 3000 GHz, in accordance with gas phase evidence for lower metallicity and depletion in these clouds. Unexpected anti-correlated variations of the dust temperature and emission cross-section per H atom are identified in the local ISM and IVCs, a trend that continues into molecular environments. This suggests that dust growth through aggregation, seen in molecular clouds, is active much earlier in the cloud condensation and star formation processes.

preprint2011arXiv

Planck Early Results: All sky temperature and dust optical depth from Planck and IRAS: Constraints on the "dark gas" in our galaxy

We construct an all-sky map of the apparent temperature and optical depth of thermal dust emission using the Planck-HFI and IRAS data. The optical depth maps are correlated to tracers of the atomic and molecular gas. The correlation is linear in the lowest column density regions at high galactic latitudes. At high NH, the correlation is consistent with that of the lowest NH. In the intermediate NH range, we observe departure from linearity, with the dust optical depth in excess to the correlation. We attribute this excess emission to thermal emission by dust associated with a Dark-Gas phase, undetected in the available HI and CO measurements. We show the 2D spatial distribution of the Dark-Gas in the solar neighborhood and show that it extends around known molecular regions traced by CO. The average dust emissivity in the HI phase in the solar neighborhood follows roughly a power law distribution with beta = 1.8 all the way down to 3 mm, although the SED flattens slightly in the millimetre. The threshold for the existence of the Dark-Gas is found at NH = (8.0\pm 0.58) 10^{20} Hcm-2. Assuming the same dust emissivity at high frequencies for the dust in the atomic and molecular phases leads to an average XCO = (2.54\pm0.13) 10^{20} H2cm-2/(K km s-1). The mass of Dark-Gas is found to be 28% of the atomic gas and 118% of the CO emitting gas in the solar neighborhood. A possible origin for the Dark-Gas is the existence of a dark molecular phase, where H2 survives photodissociation but CO does not. The observed transition for the onset of this phase in the solar neighborhood (AV = 0.4 mag) appears consistent with recent theoretical predictions. We also discuss the possibility that up to half of the Dark-Gas could be in atomic form, due to optical depth effects in the HI measurements.

preprint2011arXiv

Planck early results: first assessment of the High Frequency Instrument in-flight performance

The Planck High Frequency Instrument (HFI) is designed to measure the temperature and polarization anisotropies of the Cosmic Microwave Background and galactic foregrounds in six wide bands centered at 100, 143, 217, 353, 545 and 857 GHz at an angular resolution of 10' (100 GHz), 7' (143 GHz), and 5' (217 GHz and higher). HFI has been operating flawlessly since launch on 14 May 2009. The bolometers cooled to 100 mK as planned. The settings of the readout electronics, such as the bolometer bias current, that optimize HFI's noise performance on orbit are nearly the same as the ones chosen during ground testing. Observations of Mars, Jupiter, and Saturn verified both the optical system and the time response of the detection chains. The optical beams are close to predictions from physical optics modeling. The time response of the detection chains is close to pre-launch measurements. The detectors suffer from an unexpected high flux of cosmic rays related to low solar activity. Due to the redundancy of Planck's observations strategy, the removal of a few percent of data contaminated by glitches does not affect significantly the sensitivity. The cosmic rays heat up significantly the bolometer plate and the modulation on periods of days to months of the heat load creates a common drift of all bolometer signals which do not affect the scientific capabilities. Only the high energy cosmic rays showers induce inhomogeneous heating which is a probable source of low frequency noise.

preprint2011arXiv

Planck Early Results: Origin of the submm excess dust emission in the Magellanic Clouds

The integrated Spectral Energy Distributions of the Large and Small Magellanic Cloud appear significantly flatter than expected from dust models based on their FIR and radio emission. The origin of this millimetre excess is still unexplained, and is here investigated using the Planck data. The background CMB contribution is subtracted using an ILC method performed locally around the galaxies. The foreground emission from the Milky Way is subtracted. After subtraction, the emission of both galaxies correlates closely with the gas emission of the LMC and SMC. The millimetre excess in the LMC can be explained by CMB fluctuations, but a significant excess is still present in the SMC SED. The Planck and IRIS data at 100 micron are combined to produce thermal dust temperature and optical depth maps of the two galaxies. The LMC temperature map shows the presence of a warm inner arm already found with the Spitzer data, but also shows the existence of a previously unidentified cold outer arm. Several cold regions are found along this arm, some of which are associated with known molecular clouds. The average emissivity spectral index is found to be consistent with beta=1.5 and beta=1.2 below 500 microns for the LMC and SMC respectively, significantly flatter than the values observed in the Milky Way. Furthermore, there is evidence in the SMC for a further flattening of the SED in the sub-mm. The spatial distribution of the millimetre dust excess in the SMC follows the gas and thermal dust distribution. Different models are explored in order to fit the dust emission in the SMC. It is concluded that the millimetre excess is unlikely to be caused by very cold dust emission and that it could be due to a combination of spinning dust. emission and thermal dust emission by more amorphous dust grains than those present in our Galaxy.

preprint2011arXiv

Planck Early Results: The Galactic Cold Core Population revealed by the first all-sky survey

We present the statistical properties of the first version of the Cold Core Catalogue of Planck Objects (C3PO), in terms of their spatial distribution, temperature, distance, mass, and morphology. We also describe the statistics of the Early Cold Core Catalogue (ECC, delivered with the Early Release Compact Source Catalogue, ERCSC) that is the subset of the 915 most reliable detections of the complete catalogue. We have used the CoCoCoDeT algorithm to extract 10783 cold sources. Temperature and dust emission spectral index β values are derived using the fluxes in the IRAS 100 \mum band and the three highest frequency Planck bands. Temperature spans from 7K to 17K, and peaks around 13K. Data are not consistent with a constant value of β over the all temperature range. β ranges from 1.4 to 2.8 with a mean value around 2.1, and several possible scenarios are possible, including β(T) and the effect of multiple T components folded into the measurements. For one third of the objects the distances are obtained. Most of the detections are within 2 kpc in the Solar neighbourhood, but a few are at distances greater than 4 kpc. The cores are distributed from the deep Galactic plane, despite the confusion, to high latitudes (>30$^{\circle}$). The associated mass estimates range from 1 to $10^5$ solar masses. Using their physical properties these cold sources are shown to be cold clumps, defined as the intermediate cold sub-structures between clouds and cores. These cold clumps are not isolated but mostly organized in filaments associated with molecular clouds. The Cold Core Catalogue of Planck Objects (C3PO) is the first unbiased all-sky catalogue of cold objects. It gives an unprecedented statistical view to the properties of these potential pre-stellar clumps and offers a unique possibility for their classification in terms of their intrinsic properties and environment.

preprint2011arXiv

Planck Early Results: The Planck mission

The European Space Agency's Planck satellite was launched on 14 May 2009, and has been surveying the sky stably and continuously since 13 August 2009. Its performance is well in line with expectations, and it will continue to gather scientific data until the end of its cryogenic lifetime. We give an overview of the history of Planck in its first year of operations, and describe some of the key performance aspects of the satellite. This paper is part of a package submitted in conjunction with Planck's Early Release Compact Source Catalogue, the first data product based on Planck to be released publicly. The package describes the scientific performance of the Planck payload, and presents results on a variety of astrophysical topics related to the sources included in the Catalogue, as well as selected topics on diffuse emission.

preprint2011arXiv

Planck Early Results: The Planck View of Nearby Galaxies

The all-sky coverage of the Planck Early Release Compact Source Catalogue (ERCSC) provides an unsurpassed survey of galaxies at submillimetre (submm) wavelengths, representing a major improvement in the numbers of galaxies detected, as well as the range of far-IR/submm wavelengths over which they have been observed. We here present the first results on the properties of nearby galaxies using these data. We match the ERCSC catalogue to IRAS-detected galaxies in the Imperial IRAS Faint Source Redshift Catalogue (IIFSCz), so that we can measure the spectral energy distributions (SEDs) of these objects from 60 to 850 microns. This produces a list of 1717 galaxies with reliable associations between Planck and IRAS, from which we select a subset of 468 for SED studies, namely those with strong detections in the three highest frequency Planck bands and no evidence of cirrus contamination. The SEDs are fitted using parametric dust models to determine the range of dust temperatures and emissivities. We find evidence for colder dust than has previously been found in external galaxies, with T<20K. Such cold temperatures are found using both the standard single temperature dust model with variable emissivity beta, or a two dust temperature model with beta fixed at 2. We also compare our results to studies of distant submm galaxies (SMGs) which have been claimed to contain cooler dust than their local counterparts. We find that including our sample of 468 galaxies significantly reduces the distinction between the two populations. Fits to SEDs of selected objects using more sophisticated templates derived from radiative transfer models confirm the presence of the colder dust found through parameteric fitting. We thus conclude that cold (T<20K) dust is a significant and largely unexplored component of many nearby galaxies.

preprint2011arXiv

Planck Early Results: The submillimetre properties of a sample of Galactic cold clumps

(abridged) We perform a detailed investigation of sources from the Cold Cores Catalogue of Planck Objects (C3PO). Our goal is to probe the reliability of the detections, validate the separation between warm and cold dust emission components, provide the first glimpse at the nature, internal morphology and physical characterictics of the Planck-detected sources. We focus on a sub-sample of ten sources from the C3PO list, selected to sample different environments, from high latitude cirrus to nearby (150pc) and remote (2kpc) molecular complexes. We present Planck surface brightness maps and derive the dust temperature, emissivity spectral index, and column densities of the fields. With the help of higher resolution Herschel and AKARI continuum observations and molecular line data, we investigate the morphology of the sources and the properties of the substructures at scales below the Planck beam size.

preprint2011arXiv

Planck early results. VI. The High Frequency Instrument data processing

We describe the processing of the 336 billion raw data samples from the High Frequency Instrument (HFI) which we performed to produce six temperature maps from the first 295 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545 and 857 GHz with an angular resolution ranging from 9.9 to 4.4^2. The white noise level is around 1.5 μK degree or less in the 3 main CMB channels (100--217GHz). The photometric accuracy is better than 2% at frequencies between 100 and 353 GHz and around 7% at the two highest frequencies. The maps created by the HFI Data Processing Centre reach our goals in terms of sensitivity, resolution, and photometric accuracy. They are already sufficiently accurate and well-characterised to allow scientific analyses which are presented in an accompanying series of early papers. At this stage, HFI data appears to be of high quality and we expect that with further refinements of the data processing we should be able to achieve, or exceed, the science goals of the Planck project.

preprint2011arXiv

Planck Early Results. VII. The Early Release Compact Source Catalog

(Abridged) A brief description of the methodology of construction, contents and usage of the Planck Early Release Compact Source Catalogue (ERCSC), including the Early Cold Cores (ECC) and the Early Sunyaev-Zeldovich (ESZ) cluster catalogue is provided. The catalogue is based on data that consist of mapping the entire sky once and 60% of the sky a second time by Planck, thereby comprising the first high sensitivity radio/submillimetre observations of the entire sky. A Monte-Carlo algorithm based on the injection and extraction of artificial sources into the Planck maps was implemented to select reliable sources among all extracted candidates such that the cumulative reliability of the catalogue is >=90%. The 10sigma photometric flux density limit of the catalogue at |b|>30 deg is 0.49, 1.0, 0.67, 0.5, 0.33, 0.28, 0.25, 0.47 and 0.82 Jy at each of the nine frequencies between 30 and 857 GHz. Sources which are up to a factor of ~2 fainter than this limit, and which are present in "clean" regions of the Galaxy where the sky background due to emission from the interstellar medium is low, are included in the ERCSC if they meet the high reliability criterion. The Planck ERCSC sources have known associations to stars with dust shells, stellar cores, radio galaxies, blazars, infrared luminous galaxies and Galactic interstellar medium features. A significant fraction of unclassified sources are also present in the catalogs. In addition, two early release catalogs that contain 915 cold molecular cloud core candidates and 189 SZ cluster candidates that have been generated using multi-frequency algorithms are presented. The entire source list, with more than 15000 unique sources, is ripe for follow-up characterisation with Herschel, ATCA, VLA, SOFIA, ALMA and other ground-based observing facilities.

preprint2011arXiv

Planck early results. XX. New light on anomalous microwave emission from spinning dust grains

Anomalous microwave emission (AME) has been observed by numerous experiments in the frequency range ~10-60 GHz. Using Planck maps and multi-frequency ancillary data, we have constructed spectra for two known AME regions: the Perseus and Rho Ophiuchi molecular clouds. The spectra are well fitted by a combination of free-free radiation, cosmic microwave background, thermal dust, and electric dipole radiation from small spinning dust grains. The spinning dust spectra are the most precisely measured to date, and show the high frequency side clearly for the first time. The spectra have a peak in the range 20-40 GHz and are detected at high significances of 17.1 sigma for Perseus and 8.4 sigma for Rho Ophiuchi. In Perseus, spinning dust in the dense molecular gas can account for most of the AME; the low density atomic gas appears to play a minor role. In Rho Ophiuchi, the ~30 GHz peak is dominated by dense molecular gas, but there is an indication of an extended tail at frequencies 50-100 GHz, which can be accounted for by irradiated low density atomic gas. The dust parameters are consistent with those derived from other measurements. We have also searched the Planck map at 28.5 GHz for candidate AME regions, by subtracting a simple model of the synchrotron, free-free, and thermal dust. We present spectra for two of the candidates; S140 and S235 are bright HII regions that show evidence for AME, and are well fitted by spinning dust models.

preprint2011arXiv

The Displaced Dusty ISM of NGC 3077: Tidal Stripping in the M 81 Triplet

We present the detection of extended (~30 kpc^2) dust emission in the tidal \hi\ arm near NGC 3077 (member of the M\,81 triplet) using SPIRE on board Herschel. Dust emission in the tidal arm is typically detected where the \hi\ column densities are >10^21 cm^-2. The SPIRE band ratios show that the dust in the tidal arm is significantly colder (~13 K) than in NGC 3077 itself (~31 K), consistent with the lower radiation field in the tidal arm. The total dust mass in the tidal arm is ~1.8 x 10^6 M_sun (assuming beta=2), i.e. substantially larger than the dust mass associated with NGC 3077 (~2 x 10^5 M_sun). Where dust is detected, the dust-to-gas ratio is 6+/-3 x 10^-3, consistent within the uncertainties with what is found in NGC 3077 and nearby spiral galaxies with Galactic metallicities. The faint HII regions in the tidal arm can not be responsible for the detected enriched material and are not the main source of the dust heating in the tidal arm. We conclude that the interstellar medium (atomic HI, molecules and dust) in this tidal feature was pre-enriched and stripped off NGC 3077 during its recent interaction (~3 x 10^8 yr ago) with M 82 and M 81. This implies that interaction can efficiently remove heavy elements and enriched material (dust, molecular gas) from galaxies. As interactions were more frequent at large lookback times, it is conceivable that they could substantially contribute (along with galactic outflows) to the enrichment of the intergalactic medium.

preprint2011arXiv

The Palomar Transient Factory photometric calibration

The Palomar Transient Factory (PTF) provides multiple epoch imaging for a large fraction of the sky. Here we describe the photometric calibration of the PTF data products that relates the PTF magnitudes to other mag systems. The calibration process utilizes SDSS r~16 mag point source objects as photometric standards. During photometric conditions, this allows us to solve for the extinction coefficients and color terms, and to estimate the camera illumination correction. This also enables the calibration of fields that are outside the SDSS footprint. We test the precision and repeatability of the PTF photometric calibration. Given that PTF is observing in a single filter each night, we define a PTF calibrated magnitude system for R-band and g-band. We show that, in this system, 59% (47%) of the photometrically calibrated PTF R-band (g-band) data achieve a photometric precision of 0.02-0.04 mag. Given the objects' color, the PTF magnitude system can be converted to other systems. A night-by-night comparison of the calibrated magnitudes of individual stars observed on multiple nights shows that they are consistent to a level of ~0.02 mag. Most of the data that were taken under non-photometric conditions can be calibrated relative to other epochs of the same sky footprint obtained during photometric conditions. We provide a guide describing the use of the PTF photometric calibration data products, as well as the transformations between the PTF magnitude system and the SDSS and Johnson-Cousins systems. (abridged).

preprint2010arXiv

Far-Infrared Line Imaging of the Starburst Ring in NGC 1097 with the Herschel/PACS Spectrometer

NGC 1097 is a nearby SBb galaxy with a Seyfert nucleus and a bright starburst ring. We study the physical properties of the interstellar medium (ISM) in the ring using spatially resolved far-infrared spectral maps of the circumnuclear starburst ring of NGC 1097, obtained with the PACS spectrometer on board the Herschel Space Telescope. In particular, we map the important ISM cooling and diagnostic emission lines of [OI] 63 $μ$m, [OIII] 88 $μ$m, [NII] 122 $μ$m, [CII] 158 $μ$m and [NII] 205 $μ$m. We observe that in the [OI] 63 $μ$m, [OIII] 88 $μ$m, and [NII] 122 $μ$m line maps, the emission is enhanced in clumps along the NE part of the ring. We observe evidence of rapid rotation in the circumnuclear ring, with a rotation velocity of ~220$ km s$^{-1}$ (inclination uncorrected) measured in all lines. The [OI] 63 $μ$m/[CII] 158 $μ$m ratio varies smoothly throughout the central region, and is enhanced on the northeastern part of the ring, which may indicate a stronger radiation field. This enhancement coincides with peaks in the [OI] 63 $μ$m and [OIII] 88 $μ$m maps. Variations of the [NII] 122 $μ$m/[NII] 205 $μ$m ratio correspond to a range in the ionized gas density between 150 and 400 cm$^{-3}$.

preprint2010arXiv

Mapping far-IR emission from the central kiloparsec of NGC 1097

Using photometry of NGC 1097 from the Herschel PACS (Photodetector Array Camera and Spectrometer) instrument, we study the resolved properties of thermal dust continuum emission from a circumnuclear starburst ring with a radius ~ 900 pc. These observations are the first to resolve the structure of a circumnuclear ring at wavelengths that probe the peak (i.e. lambda ~ 100 micron) of the dust spectral energy distribution. The ring dominates the far-infrared (far-IR) emission from the galaxy - the high angular resolution of PACS allows us to isolate the ring's contribution and we find it is responsible for 75, 60 and 55% of the total flux of NGC 1097 at 70, 100 and 160 micron, respectively. We compare the far-IR structure of the ring to what is seen at other wavelengths and identify a sequence of far-IR bright knots that correspond to those seen in radio and mid-IR images. The mid- and far-IR band ratios in the ring vary by less than +/- 20% azimuthally, indicating modest variation in the radiation field heating the dust on ~ 600 pc scales. We explore various explanations for the azimuthal uniformity in the far-IR colors of the ring including a lack of well-defined age gradients in the young stellar cluster population, a dominant contribution to the far-IR emission from dust heated by older (> 10 Myr) stars and/or a quick smoothing of local enhancements in dust temperature due to the short orbital period of the ring. Finally, we improve previous limits on the far-IR flux from the inner ~ 600 pc of NGC 1097 by an order of magnitude, providing a better estimate of the total bolometric emission arising from the active galactic nucleus and its associated central starburst.

preprint2010arXiv

No Evolution in the IR-Radio Relation for IR-Luminous Galaxies at z<2 in the COSMOS Field

Previous observational studies of the infrared (IR)-radio relation out to high redshift employed any detectable star forming systems at a given redshift within the restricted area of cosmological survey fields. Consequently, the evolution inferred relies on a comparison between the average IR/radio properties of (i) very IR-luminous high-z sources and (ii) more heterogeneous low(er)-z samples that often lack the strongest IR emitters. In this report we consider populations of objects with comparable luminosities over the last 10 Gyr by taking advantage of deep IR (esp. Spitzer 24 micron) and VLA 1.4 GHz observations of the COSMOS field. Consistent with recent model predictions, both Ultra Luminous Infrared Galaxies (ULIRGs) and galaxies on the bright end of the evolving IR luminosity function do not display any change in their average IR/radio ratios out to z~2 when corrected for bias. Uncorrected data suggested ~0.3 dex of positive evolution.

preprint2010arXiv

On the black hole mass - velocity dispersion relation for type-1 and type-2 AGN

We present results from infrared spectroscopic projects that aim to test the relation between the mass of a black hole, M_BH, and the velocity dispersion of the stars in its host-galaxy bulge. We demonstrate that near-infrared, high-resolution spectroscopy assisted by adaptive optics is key in populating the high-luminosity end of the relation. We show that the velocity dispersions of mid-infrared, high-ionization lines originating from gas in the narrow-line region of the active galactic nucleus follow the same relation. This result provides a way of inferring M_BH estimates for the cosmologically significant population of obscured, type-2 AGN that can be applicable to data from spectrographs on the next generation infrared telescopes.

preprint2009arXiv

A Search for Dust Emission in the Leo Intergalactic Cloud

We present a search for infrared dust emission associated with the Leo cloud, a large intergalactic cloud in the M96 group. Mid-infrared and far-infrared images were obtained with IRAC and MIPS on the Spitzer Space Telescope. Our analysis of these maps is done at each wavelength relative to the HI spatial distribution. We observe a probable detection at 8 microns and a marginal detection at 24 microns associated with the highest HI column densities in the cloud. At 70 and 160 microns, upper limits on the dust emission are deduced. The level of the detection is low so that the possibility of a fortuitous cirrus clump or of an overdensity of extragalactic sources along the line of sight can not be excluded. If this detection is confirmed, the quantities of dust inferred imply a dust to gas ratio in the intergalactic cloud up to a few times solar but no less than 1/20 solar. A confirmed detection would therefore exclude the possibility that the intergalactic cloud has a primordial origin. Instead, this large intergalactic cloud could therefore have been formed through interactions between galaxies in the group.

preprint2009arXiv

A Spitzer Unbiased Ultradeep Spectroscopic Survey

We carried out an unbiased, spectroscopic survey using the low-resolution module of the infrared spectrograph (IRS) on board Spitzer targeting two 2.6 square arcminute regions in the GOODS-North field. IRS was used in spectral mapping mode with 5 hours of effective integration time per pixel. One region was covered between 14 and 21 microns and the other between 20 and 35 microns. We extracted spectra for 45 sources. About 84% of the sources have reported detections by GOODS at 24 microns, with a median F_nu(24um) ~ 100 uJy. All but one source are detected in all four IRAC bands, 3.6 to 8 microns. We use a new cross-correlation technique to measure redshifts and estimate IRS spectral types; this was successful for ~60% of the spectra. Fourteen sources show significant PAH emission, four mostly SiO absorption, eight present mixed spectral signatures (low PAH and/or SiO) and two show a single line in emission. For the remaining 17, no spectral features were detected. Redshifts range from z ~ 0.2 to z ~ 2.2, with a median of 1. IR Luminosities are roughly estimated from 24 microns flux densities, and have median values of 2.2 x 10^{11} L_{\odot} and 7.5 x 10^{11} L_{\odot} at z ~ 1 and z ~ 2 respectively. This sample has fewer AGN than previous faint samples observed with IRS, which we attribute to the fainter luminosities reached here.

preprint2009arXiv

Spitzer 70 and 160-micron Observations of the COSMOS Field

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

preprint2009arXiv

The Detection of Anomalous Dust Emission in the Nearby Galaxy NGC 6946

We report on the Ka-band (26-40 GHz) emission properties for 10 star-forming regions in the nearby galaxy NGC 6946. From a radio spectral decomposition, we find that the 33 GHz flux densities are typically dominated by thermal (free-free) radiation. However, we also detect excess Ka-band emission for an outer-disk star-forming region relative to what is expected given existing radio, submillimeter, and infrared data. Among the 10 targeted regions, measurable excess emission at 33 GHz is detected for half of them, but in only one region is the excess found to be statistically significant ($\approx7σ$). We interpret this as the first likely detection of so called `anomalous' dust emission outside of the Milky Way. We find that models explaining this feature as the result of dipole emission from rapidly rotating ultrasmall grains are able to reproduce the observations for reasonable interstellar medium conditions. While these results suggest that the use of Ka-band data as a measure of star formation activity in external galaxies may be complicated by the presence of anomalous dust, it is unclear how significant a factor this will be for globally integrated measurements as the excess emission accounts for $\la$10% of the total Ka-band flux density from all 10 regions.

preprint2009arXiv

The Palomar Transient Factory: System Overview, Performance and First Results

The Palomar Transient Factory (PTF) is a fully-automated, wide-field survey aimed at a systematic exploration of the optical transient sky. The transient survey is performed using a new 8.1 square degree camera installed on the 48-inch Samuel Oschin telescope at Palomar Observatory; colors and light curves for detected transients are obtained with the automated Palomar 60-inch telescope. PTF uses eighty percent of the 1.2-m and fifty percent of the 1.5-m telescope time. With an exposure of 60-s the survey reaches a depth of approximately 21.3 in g' and 20.6 in R (5 sigma, median seeing). Four major experiments are planned for the five-year project: 1) a 5-day cadence supernova search; 2) a rapid transient search with cadences between 90 seconds and 1 day; 3) a search for eclipsing binaries and transiting planets in Orion; and 4) a 3-pi sr deep H-alpha survey. PTF provides automatic, realtime transient classification and follow up, as well as a database including every source detected in each frame. This paper summarizes the PTF project, including several months of on-sky performance tests of the new survey camera, the observing plans and the data reduction strategy. We conclude by detailing the first 51 PTF optical transient detections, found in commissioning data.

preprint2008arXiv

Environmental Effects in Clusters: Modified Far-Infrared--Radio Relations within Virgo Cluster Galaxies

(abridged) We present a study on the effects of the intracluster medium (ICM) on the interstellar medium (ISM) of 10 Virgo cluster spiral galaxies using {\it Spitzer} far-infrared (FIR) and VLA radio continuum imaging. Relying on the FIR-radio correlation within normal galaxies, we use our infrared data to create model radio maps which we compare to the observed radio images. For 6 of our sample galaxies we find regions along their outer edges that are highly deficient in the radio compared with our models. We believe these observations are the signatures of ICM ram pressure. For NGC 4522 we find the radio deficit region to lie just exterior to a region of high radio polarization and flat radio spectral index, although the total 20 cm radio continuum in this region does not appear strongly enhanced. These characteristics seem consistent for other galaxies with radio polarization data in the literature. The strength of the radio deficit is inversely correlated with the time since peak pressure as inferred from stellar population studies and gas stripping simulations, suggesting the strength of the radio deficit is good indicator of the strength of the current ram pressure. We also find that galaxies having {\it local} radio {\it deficits} appear to have {\it enhanced global} radio fluxes. Our preferred physical picture is that the observed radio deficit regions arise from the ICM wind sweeping away cosmic-ray (CR) electrons and the associated magnetic field, thereby creating synchrotron tails as observed for some of our galaxies. We propose that CR particles are also re-accelerated by ICM-driven shocklets behind the observed radio deficit regions which in turn enhances the remaining radio disk brightness.

preprint2008arXiv

The relations among 8, 24, and 160 micron dust emission within nearby spiral galaxies

We investigate the relations among the stellar continuum-subtracted 8 micron polycyclic aromatic hydrocarbon (PAH 8 micron) emission, 24 micron hot dust emission, and 160 micron cold dust emission in fifteen nearby face-on spiral galaxies in the Spitzer Infrared Nearby Galaxies Survey sample. The relation between PAH 8 and 24 micron emission measured in ~2 kpc regions is found to exhibit a significant amount of scatter, and strong spatial variations are observed in the (PAH 8 micron)/24 micron surface brightness ratio. In particular, the (PAH 8 micron)/24 micron surface brightness ratio is observed to be high in the diffuse interstellar medium and low in bright star-forming regions and other locations with high 24 micron surface brightness. PAH 8 micron emission is found to be well-correlated with 160 micron emission on spatial scales of ~2 kpc, and the (PAH 8 micron)/160 micron surface brightness ratio is generally observed to increase as the 160 micron surface brightness increases. These results suggest that the PAHs are associated with the diffuse, cold dust that produces most of the 160 micron emission in these galaxies, and the variations in the (PAH 8 micron)/160 micron ratio may generally be indicative of either the intensity or the spectrum of the interstellar radiation field that is heating both the PAHs and the diffuse interstellar dust.

preprint2008arXiv

The Spitzer Infrared Nearby Galaxies Survey: A High-Resolution Spectroscopy Anthology

High resolution mid-infrared spectra are presented for 155 nuclear and extranuclear regions from the Spitzer Infrared Nearby Galaxies Survey (SINGS). The fluxes for nine atomic forbidden and three molecular hydrogen mid-infrared emission lines are also provided, along with upper limits in key lines for infrared-faint targets. The SINGS sample shows a wide range in the ratio of [SIII]18.71um/[SIII]33.48um, but the average ratio of the ensemble indicates a typical interstellar electron density of 300-400 cm^{-3} on ~23"x15" scales and 500-600 cm^{-3} using ~11"x9" apertures, independent of whether the region probed is a star-forming nuclear, a star-forming extranuclear, or an AGN environment. Evidence is provided that variations in gas-phase metallicity play an important role in driving variations in radiation field hardness, as indicated by [NeIII]15.56um/[NeII]12.81um, for regions powered by star formation. Conversely, the radiation hardness for galaxy nuclei powered by accretion around a massive black hole is independent of metal abundance. Furthermore, for metal-rich environments AGN are distinguishable from star-forming regions by significantly larger [NeIII]15.56um/[NeII]12.81um ratios. Finally, [FeII]25.99um/[NeII]12.81um versus [SiII]34.82um/[SIII]33.48um also provides an empirical method for discerning AGN from normal star-forming sources. However, similar to [NeIII]15.56um/[NeII]12.81um, these mid-infrared line ratios lose their AGN/star-formation diagnostic powers for very low metallicity star-forming systems with hard radiation fields.

preprint2007arXiv

Dust and Atomic Gas in Dwarf Irregular Galaxies of the M81 Group: The SINGS and THINGS view

[abridged] We present observations of the dust and atomic gas phase in seven dwarf irregular galaxies of the M81 group from the SINGS and THINGS surveys. The Spitzer observations provide a first glimpse of the nature of the non-atomic ISM in these metal-poor (Z~0.1 Z_sun), quiescent (SFR~0.001-0.1 M_sun/yr) dwarf galaxies. Dust emission is detected in five out of the seven targets. Most detected dust emission is restricted to HI column densities >1x10^21 cm^-2. Spitzer spectroscopy of two regions in the brightest galaxies (IC 2574 and Holmberg II) show distinctly different spectral shapes. The spectrum of IC 2574 shows aromatic features that are less luminous (relative to the FIR luminosity) compared to an average SINGS spiral galaxy by a factor of \~7 . The aromatic features in Holmberg~II (which has only a slightly lower gas-phase metallicity) are fainter than in IC 2574 by an order of magnitude. This result emphazises that the strength of the aromatic features is not a simple linear function of metallicity. We estimate dust masses of ~10^4-10^6 M_sun for the M81 dwarf galaxies, resulting in an average dust--to--gas ratio (M_dust/M_HI) of ~3x10^-4 (1.5x10^-3 if only the HI that is associated with dust emission is considered); this is an order of magnitude lower than the typical value derived for the SINGS spirals. The dwarf galaxies are underluminous per unit star formation rate at 70um as compared to the more massive galaxies in SINGS by a factor of ~2. However, the average 70um/160um ratio in the sample dwarf galaxies is higher than what is found in the other galaxies of the SINGS sample. This can be explained by a combination of a lower dust content in conjunction with a higher dust temperature in the dwarfs.

preprint2007arXiv

The infrared luminosity function of galaxies at redshifts z=1 and z~2 in the GOODS fields

We present the rest-frame 8 micron luminosity function (LF) at redshifts z=1 and ~2, computed from Spitzer 24 micron-selected galaxies in the GOODS fields over an area of 291 sq. arcmin. Using classification criteria based on X-ray data and IRAC colours, we identify the AGN in our sample. The rest-frame 8 micron LF for star-forming galaxies at redshifts z=1 and ~2 have the same shape as at z~0, but with a strong positive luminosity evolution. The number density of star-forming galaxies with log_{10}(nu L_nu(8 micron))>11 increases by a factor >250 from redshift z~0 to 1, and is basically the same at z=1 and ~2. The resulting rest-frame 8 micron luminosity densities associated with star formation at z=1 and ~2 are more than four and two times larger than at z~0, respectively. We also compute the total rest-frame 8 micron LF for star-forming galaxies and AGN at z~2 and show that AGN dominate its bright end, which is well-described by a power-law. Using a new calibration based on Spitzer star-forming galaxies at 0<z<0.6 and validated at higher redshifts through stacking analysis, we compute the bolometric infrared (IR) LF for star-forming galaxies at z=1 and ~2. We find that the respective bolometric IR luminosity densities are (1.2+/-0.2) x 10^9 and (6.6^{+1.2}_{-1.0}) x 10^8 L_sun Mpc^{-3}, in agreement with previous studies within the error bars. At z~2, around 90% of the IR luminosity density associated with star formation is produced by luminous and ultraluminous IR galaxies (LIRG and ULIRG), with the two populations contributing in roughly similar amounts. Finally, we discuss the consistency of our findings with other existing observational results on galaxy evolution.

preprint2007arXiv

The Role of Galaxy Interactions and Mergers in Star Formation at z<1.3: Mid-Infrared Properties in the Spitzer First Look Survey

By combining the 0.12 square degree F814W Hubble Space Telescope (HST) and Spitzer MIPS 24 micron imaging in the First Look Survey (FLS), we investigate the properties of interacting and merging Mid-Infrared bright and faint sources at 0.2< z <1.3. We find a marginally significant increase in the pair fraction for MIPS 24 micron detected, optically selected close pairs, pair fraction=0.25+/-0.10 at z~1, in contrast to 0.11+/-0.08 at z~0.4, while galaxies below our 24 micron MIPS detection limit show a pair fraction consistent with zero at all redshifts. Additionally, 24 micron detected galaxies with fluxes >0.1mJy are on average five times more likely to be in a close galaxy pair between 0.2< z <1.3 than galaxies below this flux limit. Using the 24 micron flux to derive the total Far-IR luminosity we find that paired galaxies (early stage mergers) are responsible for 27% +/-9% of the IR luminosity density resulting from star formation at z~1 while morphologically classified (late stage) mergers make up 34%+/-11%. This implies that 61%+/-14% of the infrared luminosity density and in turn ~40% of the star formation rate density at z~1 can be attributed to galaxies at some stage of a major merger or interaction. We argue that, close pairs/mergers in a LIRG/ULIRG phase become increasingly important contributers to the IR luminosity and star formation rate density of the Universe at z>0.7.

preprint2007arXiv

The Spitzer Legacy Survey of the HST-ACS 2 sq. deg. COSMOS Field: survey strategy and first analysis

The Spitzer-COSMOS survey (S-COSMOS) is a Legacy program (Cycles 2+3) designed to carry out a uniform deep survey of the full 2 sq deg COSMOS field in all seven Spitzer bands (3.6, 4.5, 5.6, 8.0, 24.0, 70.0, 160.0 u). This paper describes the survey parameters, mapping strategy, data reduction procedures, achieved sensitivities to date, and the complete data set for future reference. We show that the observed infrared backgrounds in the S-COSMOS field are within 10% of the predicted background levels. The fluctuations in the background at 24u have been measured and do not show any significant contribution from cirrus, as expected. In addition, we report on the number of asteroid detections in the low galactic latitude COSMOS field. We use the Cycle 2 S-COSMOS data to determine preliminary number counts, and compare our results with those from previous Spitzer Legacy surveys (e.g. SWIRE, GOODS). The results from this "first analysis" confirm that the S-COSMOS survey will have sufficient sensitivity with IRAC to detect ~ L* disks and spheroids out to z ~ 3, and with MIPS to detect ultraluminous starbursts and AGN out to z ~3 at 24u and out to z ~1.5-2 at 70u and 160u.

preprint2007arXiv

Warm molecular hydrogen in the Spitzer SINGS galaxy sample

(simplified) Results on the properties of warm H2 in 57 normal galaxies are derived from H2 rotational transitions, obtained as part of SINGS. This study extends previous extragalactic surveys of H2, the most abundant constituent of the molecular ISM, to more common systems (L_FIR = e7 to 6e10 L_sun) of all morphological and nuclear types. The S(1) transition is securely detected in the nuclear regions of 86% of SINGS galaxies with stellar masses above 10^9.5 M_sun. The derived column densities of warm H2 (T > ~100 K), even though averaged over kiloparsec-scale areas, are commensurate with those of resolved PDRs; the median of the sample is 3e20 cm-2. They amount to between 1% and >30% of the total H2. The power emitted in the sum of the S(0) to S(2) transitions is on average 30% of the [SiII] line power, and ~4e-4 of the total infrared power (TIR) within the same area for star-forming galaxies, which is consistent with excitation in PDRs. The fact that H2 emission scales tightly with PAH emission, even though the average radiation field intensity varies by a factor ten, can also be understood if both tracers originate predominantly in PDRs, either dense or diffuse. A large fraction of the 25 LINER/Sy targets, however, strongly depart from the rest of the sample, in having warmer H2 in the excited states, and an excess of H2 emission with respect to PAHs, TIR and [SiII]. We propose a threshold in H2 to PAH power ratios, allowing the identification of low-luminosity AGNs by an excess H2 excitation. A dominant contribution from shock heating is favored in these objects. Finally, we detect, in nearly half the star-forming targets, non-equilibrium ortho to para ratios, consistent with FUV pumping combined with incomplete ortho-para thermalization by collisions, or possibly non-equilibrium PDR fronts advancing into cold gas.

preprint2006arXiv

Spitzer and JCMT Observations of the Active Galactic Nucleus in the Sombrero Galaxy (NGC 4594)

We present Spitzer 3.6-160 micron images, Spitzer mid-infrared spectra, and JCMT SCUBA 850 micron images of the Sombrero Galaxy (NGC 4594), an Sa galaxy with a 10^9 M_solar low luminosity active galactic nucleus (AGN). The brightest infrared sources in the galaxy are the nucleus and the dust ring. The spectral energy distribution of the AGN demonstrates that, while the environment around the AGN is a prominent source of mid-infrared emission, it is a relatively weak source of far-infrared emission, as had been inferred for AGN in previous research. The weak nuclear 160 micron emission and the negligible polycyclic aromatic hydrocarbon emission from the nucleus also implies that the nucleus is a site of only weak star formation activity and the nucleus contains relatively little cool interstellar gas needed to fuel such activity. We propose that this galaxy may be representative of a subset of low ionization nuclear emission region galaxies that are in a quiescent AGN phase because of the lack of gas needed to fuel circumnuclear star formation and Seyfert-like AGN activity. Surprisingly, the AGN is the predominant source of 850 micron emission. We examine the possible emission mechanisms that could give rise to the 850 micron emission and find that neither thermal dust emission, CO line emission, bremsstrahlung emission, nor the synchrotron emission observed at radio wavelengths can adequately explain the measured 850 micron flux density by themselves. The remaining possibilities for the source of the 850 micron emission include a combination of known emission mechanisms, synchrotron emission that is self-absorbed at wavelengths longer than 850 microns, or unidentified spectral lines in the 850 micron band.

preprint2006arXiv

The Effect of Star Formation on the Far-Infrared--Radio Correlation within Galaxies

Using data obtained for twelve galaxies as part of the {\it Spitzer} Infrared Nearby Galaxies Survey (SINGS) and the Westerbork Synthesis Radio Telescope (WSRT)-SINGS radio continuum survey, we study how star formation activity affects the far-infrared (FIR)--radio correlation {\it within} galaxies by testing a phenomenological model, which describes the radio image as a smeared version of the FIR image. The physical basis of this description is that cosmic-ray (CR) electrons will diffuse measurably farther than the mean free path of dust-heating photons before decaying by synchrotron radiation. This description works well in general. Galaxies with higher infrared surface brightnesses have best-fit smoothing scale-lengths of a few hundred parsecs, substantially shorter than those for lower surface brightness galaxies. We interpret this result to suggest that galaxies with higher disk averaged star formation rates have had a recent episode of enhanced star formation and are characterized by a higher fraction of young CR electrons that have traveled only a few hundred parsecs from their acceleration sites in supernova remnants compared to galaxies with lower star formation activity.

preprint2006arXiv

The Mid-Infrared Spectrum of Star-Forming Galaxies: Global Properties of PAH Emission

We present a sample of low-resolution 5-38um Spitzer IRS spectra of the inner few square kiloparsecs of 59 nearby galaxies spanning a large range of star formation properties. A robust method for decomposing mid-infrared galaxy spectra is described, and used to explore the behavior of PAH emission and the prevalence of silicate dust extinction. Evidence for silicate extinction is found in ~1/8 of the sample, at strengths which indicate most normal galaxies undergo A_V < ~3 magnitudes averaged over their centers. The contribution of PAH emission to the total infrared power is found to peak near 10% and extend up to ~20%, and is suppressed at metallicities Z < ~Z_sun/4, as well as in low-luminosity AGN environments. Strong inter-band PAH feature strength variations (2-5x) are observed, with the presence of a weak AGN and, to a lesser degree, increasing metallicity shifting power to the longer wavelength bands. A peculiar PAH emission spectrum with markedly diminished 5-8um features arises among the sample solely in systems with relatively hard radiation fields harboring low-luminosity AGN. The AGN may modify the emitting grain distribution and provide the direct excitation source of the unusual PAH emission, which cautions against using absolute PAH strength to estimate star formation rates in systems harboring active nuclei. Alternatively, the low star formation intensity often associated with weak AGN may affect the spectrum. The effect of variations in the mid-infrared spectrum on broadband infrared surveys is modeled, and points to more than a factor of two uncertainty in results which assume a fixed PAH emission spectrum, for redshifts z=0-2.5.

preprint2005arXiv

Spitzer 70 and 160-micron Observations of the Extragalactic First Look Survey

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

preprint2005arXiv

Star Formation Rates and Extinction Properties of IR-Luminous Galaxies in the Spitzer First Look Survey

We investigate the instantaneous star formation rates (SFR) and extinction properties for a large (N=274), near-infrared (2.2 micron) + mid-infrared (24 micron) selected sample of normal to ultra-luminous infrared galaxies (ULIRGs) [10^9<L_{IR}/L_{sun}<10^12.5] with <z>~0.8 in the Spitzer Extragalactic First Look Survey. We combine Spitzer MIPS 24-micron observations with high-resolution, optical Keck Deimos spectroscopy to derive optical emission-line and infrared star formation rates (SFR_{opt} & SFR_{IR}, respectively). Direct comparison of these SFR diagnostics reveals that our sample exhibits a wide range of extinction (1.0<Av<4.0 mag). This is after removing spectroscopic and IRAC color-selected AGN candidates that account for 12% of the sample. Objects with SFRs of a few solar masses per year have Av values consistent with those of normal spirals (Av~1.0 mag). By contrast, LIRGs at z>1, which make up a large fraction of our sample, have SFR~100 M_{sun}/yr and <Av>~2.5 mag. This translates to a 97% mean attenuation of the [OII] forbidden line doublet, with the most extreme sources having as much as 99.7% of their [OII] line flux extinguished by dust. Based on a SFR_{IR}/SFR_{opt} diagnostic, we derive an IR-luminosity-dependent Av^{IR} function [Av^{IR}=0.75*log(L_{IR}/L_{sun})-6.35 mag] that we use to extinction correct our emission line luminosities. Application of this correction results in a correlation between SFR_{IR} and SFR_{opt} that has a dispersion of 0.2 dex. Investigation of the Av dependence on redshift reveals that for a fixed L_{IR}, there is no significant Av evolution. The mean attenuation of our sample is intermediate between that of local optical/UV- and radio-selected samples and has a marginally stronger L_{IR} dependence.

preprint2004arXiv

IRAC Observations of M81

IRAC images of M81 show three distinct morphological constituents: a smooth distribution of evolved stars with bulge, disk, and spiral arm components; a clumpy distribution of dust emission tracing the spiral arms; and a pointlike nuclear source. The bulge stellar colors are consistent with M-type giants, and the disk colors are consistent with a slightly younger population. The dust emission generally follows the blue and ultraviolet emission, but there are large areas that have dust emission without ultraviolet and smaller areas with ultraviolet but little dust emission. The former are presumably caused by extinction, and the latter may be due to cavities in the gas and dust created by supernova explosions. The nucleus appears fainter at 8 um than expected from ground-based 10 um observations made four years ago.

preprint2004arXiv

Mid-Infrared IRS Spectroscopy of NGC 7331: A First Look at the SINGS Legacy

The nearby spiral galaxy NGC 7331 was spectrally mapped from 5-38um using all modules of Spitzer's IRS spectrograph. A strong new dust emission feature, presumed due to PAHs, was discovered at 17.1um. The feature's intensity is nearly half that of the ubiquitous 11.3um band. The 7-14um spectral maps revealed significant variation in the 7.7 and 11.3um PAH features between the stellar ring and nucleus. Weak [OIV] 25.9um line emission was found to be centrally concentrated in the nucleus, with an observed strength over 10% of the combined neon line flux, indicating an AGN or unusually active massive star photo-ionization. Two [SIII] lines fix the characteristic electron density in the HII regions at n_e < ~200 cm^-3. Three detected H_2 rotational lines, tracing warm molecular gas, together with the observed IR continuum, are difficult to match with standard PDR models. Either additional PDR heating or shocks are required to simultaneously match lines and continuum.

preprint2004arXiv

Spatially Resolved Ultraviolet, H-alpha, Infrared, and Radio Star Formation in M81

We present Multiband Imaging Photometer for Spitzer (MIPS) observations of M81 at 24, 70, and 160 micron. The grand design nature of M81 is clearly seen, showing two well resolved spiral arms containing numerous bright star forming regions. The MIPS images reveal a significant amount of cold dust associated with the spiral arms. We investigate the variation of the ultraviolet (UV), H-alpha, and infrared (IR) luminosities and star formation rate (SFR) indicators across the face of M81 using the MIPS images and archival UV and H-alpha images. For regions in M81, we find that UV and H-alpha SFRs (uncorrected for dust attenuation) are always lower than the IR SFR. The cause of this behavior is dust attenuation and/or using SFR calibrations appropriate for entire galaxies, not regions in galaxies. The characteristics of the dust attenuation for the regions indicate the dust grains and/or geometry are different from those in starburst galaxies. The behavior of the infrared-radio correlation in M81 is seen to vary from the global average, with variations correlated with the morphology of M81.

preprint2004arXiv

The Anatomy of Star Formation in NGC 300

The Spitzer Space Telescope was used to study the mid-infrared to far-infrared properties of NGC 300, and to compare dust emission to Halpha to elucidate the heating of the ISM and the star formation cycle at scales < 100 pc. The new data allow us to discern clear differences in the spatial distribution of 8 micron dust emission with respect to 24 micron dust and to HII regions traced by the Halpha light. The 8 micron emission highlights the rims of HII regions, and the 24 micron emission is more strongly peaked in star forming regions than at 8 microns. We confirm the existence and approximate amplitude of interstellar dust emission at 4.5 microns, detected statistically in Infrared Space Observatory (ISO) data, and conclude it arises in star forming regions. When averaging over regions larger than ~ 1 kpc, the ratio of Halpha to Aromatic Feature emission in NGC 300 is consistent with the values observed in disks of spiral galaxies. The mid-to-far-infrared spectral energy distribution of dust emission is generally consistent with pre-Spitzer models.

preprint2003arXiv

SINGS: The SIRTF Nearby Galaxies Survey

The SIRTF Nearby Galaxy Survey is a comprehensive infrared imaging and spectroscopic survey of 75 nearby galaxies. Its primary goal is to characterize the infrared emission of galaxies and their principal infrared-emitting components, across a broad range of galaxy properties and star formation environments. SINGS will provide new insights into the physical processes connecting star formation to the interstellar medium properties of galaxies, and provide a vital foundation for understanding infrared observations of the distant universe and ultraluminous and active galaxies. The galaxy sample and observing strategy have been designed to maximize the scientific and archival value of the data set for the SIRTF user community at large. The SIRTF images and spectra will be supplemented by a comprehensive multi-wavelength library of ancillary and complementary observations, including radio continuum, HI, CO, submillimeter, BVRIJHK, H-alpha, Paschen-alpha, ultraviolet, and X-ray data. This paper describes the main astrophysical issues to be addressed by SINGS, the galaxy sample and the observing strategy, and the SIRTF and other ancillary data products.