We present a full-sky 100 μm map that is a reprocessed composite of the COBE/DIRBE and IRAS/ISSA maps, with the zodiacal foreground and confirmed point sources removed. Before using the ISSA maps, we remove the remaining artifacts from the IRAS scan pattern. Using the DIRBE 100 and 240 μm data, we have constructed a map of the dust temperature so that the 100 μm map may be converted to a map proportional to dust column density. The dust temperature varies from 17 to 21 K, which is modest but does modify the estimate of the dust column by a factor of 5. The result of these manipulations is a map with DIRBE quality calibration and IRAS resolution. A wealth of filamentary detail is apparent on many different scales at all Galactic latitudes. In high-latitude regions, the dust map correlates well with maps of H I emission, but deviations are coherent in the sky and are especially conspicuous in regions of saturation of H I emission toward denser clouds and of formation of H2 in molecular clouds. In contrast, high-velocity H I clouds are deficient in dust emission, as expected. To generate the full-sky dust maps, we must first remove zodiacal light contamination, as well as a possible cosmic infrared background (CIB). This is done via a regression analysis of the 100 μm DIRBE map against the Leiden-Dwingeloo map of H I emission, with corrections for the zodiacal light via a suitable expansion of the DIRBE 25 μm flux. This procedure removes virtually all traces of the zodiacal foreground. For the 100 μm map no significant CIB is detected. At longer wavelengths, where the zodiacal contamination is weaker, we detect the CIB at surprisingly high flux levels of 32 ± 13 nW m-2 sr-1 at 140 μm and of 17 ± 4 nW m-2 sr-1 at 240 μm (95% confidence). This integrated flux ~2 times that extrapolated from optical galaxies in the Hubble Deep Field. The primary use of these maps is likely to be as a new estimator of Galactic extinction. To calibrate our maps, we assume a standard reddening law and use the colors of elliptical galaxies to measure the reddening per unit flux density of 100 μm emission. We find consistent calibration using the B-R color distribution of a sample of the 106 brightest cluster ellipticals, as well as a sample of 384 ellipticals with B-V and Mg line strength measurements. For the latter sample, we use the correlation of intrinsic B-V versus Mg2 index to tighten the power of the test greatly. We demonstrate that the new maps are twice as accurate as the older Burstein-Heiles reddening estimates in regions of low and moderate reddening. The maps are expected to be significantly more accurate in regions of high reddening. These dust maps will also be useful for estimating millimeter emission that contaminates cosmic microwave background radiation experiments and for estimating soft X-ray absorption. We describe how to access our maps readily for general use.

/pdf/maps-of-dust-infrared-emission-for-use-in-estimation-of-207gvl50wo.pdf

Maps of Dust Infrared Emission for Use in Estimation of Reddening and Cosmic Microwave Background Radiation Foregrounds

We review recent determinations of the present-day mass function (PDMF) and initial mass function (IMF) in various components of the Galaxy—disk, spheroid, young, and globular clusters—and in conditions characteristic of early star formation. As a general feature, the IMF is found to depend weakly on the environment and to be well described by a power-law form forM , and a lognormal form below, except possibly for m!1 early star formation conditions. The disk IMF for single objects has a characteristic mass around M , m!0.08 c and a variance in logarithmic mass , whereas the IMF for multiple systems hasM , and . j!0.7 m!0.2 j!0.6 c The extension of the single MF into the brown dwarf regime is in good agreement with present estimates of L- and T-dwarf densities and yields a disk brown dwarf number density comparable to the stellar one, n!n! BD " pc !3 .T he IMF of young clusters is found to be consistent with the disk fi eld IMF, providing the same correction 0.1 for unresolved binaries, confirming the fact that young star clusters and disk field stars represent the same stellar population. Dynamical effects, yielding depletion of the lowest mass objects, are found to become consequential for ages!130 Myr. The spheroid IMF relies on much less robust grounds. The large metallicity spread in the local subdwarf photometric sample, in particular, remains puzzling. Recent observations suggest that there is a continuous kinematic shear between the thick-disk population, present in local samples, and the genuine spheroid one. This enables us to derive only an upper limit for the spheroid mass density and IMF. Within all the uncertainties, the latter is found to be similar to the one derived for globular clusters and is well represented also by a lognormal form with a characteristic mass slightly larger than for the disk, M , ,e xcluding as ignif icant population of m!0.2-0.3 c brown dwarfs in globular clusters and in the spheroid. The IMF characteristic of early star formation at large redshift remains undetermined, but different observational constraints suggest that it does not extend below!1M , .T hese results suggest a characteristic mass for star formation that decreases with time, from conditions prevailing at large redshift to conditions characteristic of the spheroid (or thick disk) to present-day conditions.Theseconclusions,however, remain speculative, given the large uncertainties in the spheroid and early star IMF determinations. These IMFs allow a reasonably robust determination of the Galactic present-day and initial stellar and brown dwarf contents. They also have important galactic implications beyond the Milky Way in yielding more accurate mass-to-light ratio determinations. The mass-to-light ratios obtained with the disk and the spheroid IMF yield values 1.8-1.4 times smaller than for a Salpeter IMF, respectively, in agreement with various recent dynamical determinations. This general IMF determination is examined in the context of star formation theory. None of the theories based on a Jeans-type mechanism, where fragmentation is due only to gravity, can fulfill all the observational constraints on star formation and predict a large number of substellar objects. On the other hand, recent numerical simulations of compressible turbulence, in particular in super-Alfvenic conditions, seem to reproduce both qualitatively and quantitatively the stellar and substellar IMF and thus provide an appealing theoretical foundation. In this picture, star formation is induced by the dissipation of large-scale turbulence to smaller scales through radiative MHD shocks, producing filamentary structures. These shocks produce local nonequilibrium structures with large density contrasts, which collapse eventually in gravitationally bound objects under the combined influence of turbulence and gravity. The concept of a single Jeans mass is replaced by a distribution of local Jeans masses, representative of the lognormal probability density function of the turbulent gas. Objects below the mean thermal Jeans mass still have a possibility to collapse, although with a decreasing probability.

/pdf/galactic-stellar-and-substellar-initial-mass-function-1fntlgtk9q.pdf

Galactic stellar and substellar initial mass function

This paper presents the first cosmological results based on Planck measurements of the cosmic microwave background (CMB) temperature and lensing-potential power spectra. We find that the Planck spectra at high multipoles (l ≳ 40) are extremely well described by the standard spatially-flat six-parameter ΛCDM cosmology with a power-law spectrum of adiabatic scalar perturbations. Within the context of this cosmology, the Planck data determine the cosmological parameters to high precision: the angular size of the sound horizon at recombination, the physical densities of baryons and cold dark matter, and the scalar spectral index are estimated to be θ∗ = (1.04147 ± 0.00062) × 10-2, Ωbh2 = 0.02205 ± 0.00028, Ωch2 = 0.1199 ± 0.0027, and ns = 0.9603 ± 0.0073, respectively(note that in this abstract we quote 68% errors on measured parameters and 95% upper limits on other parameters). For this cosmology, we find a low value of the Hubble constant, H0 = (67.3 ± 1.2) km s-1 Mpc-1, and a high value of the matter density parameter, Ωm = 0.315 ± 0.017. These values are in tension with recent direct measurements of H0 and the magnitude-redshift relation for Type Ia supernovae, but are in excellent agreement with geometrical constraints from baryon acoustic oscillation (BAO) surveys. Including curvature, we find that the Universe is consistent with spatial flatness to percent level precision using Planck CMB data alone. We use high-resolution CMB data together with Planck to provide greater control on extragalactic foreground components in an investigation of extensions to the six-parameter ΛCDM model. We present selected results from a large grid of cosmological models, using a range of additional astrophysical data sets in addition to Planck and high-resolution CMB data. None of these models are favoured over the standard six-parameter ΛCDM cosmology. The deviation of the scalar spectral index from unity isinsensitive to the addition of tensor modes and to changes in the matter content of the Universe. We find an upper limit of r0.002< 0.11 on the tensor-to-scalar ratio. There is no evidence for additional neutrino-like relativistic particles beyond the three families of neutrinos in the standard model. Using BAO and CMB data, we find Neff = 3.30 ± 0.27 for the effective number of relativistic degrees of freedom, and an upper limit of 0.23 eV for the sum of neutrino masses. Our results are in excellent agreement with big bang nucleosynthesis and the standard value of Neff = 3.046. We find no evidence for dynamical dark energy; using BAO and CMB data, the dark energy equation of state parameter is constrained to be w = -1.13-0.10+0.13. We also use the Planck data to set limits on a possible variation of the fine-structure constant, dark matter annihilation and primordial magnetic fields. Despite the success of the six-parameter ΛCDM model in describing the Planck data at high multipoles, we note that this cosmology does not provide a good fit to the temperature power spectrum at low multipoles. The unusual shape of the spectrum in the multipole range 20 ≲ l ≲ 40 was seen previously in the WMAP data and is a real feature of the primordial CMB anisotropies. The poor fit to the spectrum at low multipoles is not of decisive significance, but is an “anomaly” in an otherwise self-consistent analysis of the Planck temperature data.

/pdf/planck-2013-results-xvi-cosmological-parameters-12gx9fa74h.pdf

Planck 2013 results. XVI. Cosmological parameters

We present the first results based on Planck measurements of the CMB temperature and lensing-potential power spectra. The Planck spectra at high multipoles are extremely well described by the standard spatially-flat six-parameter LCDM cosmology. In this model Planck data determine the cosmological parameters to high precision. We find a low value of the Hubble constant, H0=67.3+/-1.2 km/s/Mpc and a high value of the matter density parameter, Omega_m=0.315+/-0.017 (+/-1 sigma errors) in excellent agreement with constraints from baryon acoustic oscillation (BAO) surveys. Including curvature, we find that the Universe is consistent with spatial flatness to percent-level precision using Planck CMB data alone. We present results from an analysis of extensions to the standard cosmology, using astrophysical data sets in addition to Planck and high-resolution CMB data. None of these models are favoured significantly over standard LCDM. The deviation of the scalar spectral index from unity is insensitive to the addition of tensor modes and to changes in the matter content of the Universe. We find a 95% upper limit of r<0.11 on the tensor-to-scalar ratio. There is no evidence for additional neutrino-like relativistic particles. Using BAO and CMB data, we find N_eff=3.30+/-0.27 for the effective number of relativistic degrees of freedom, and an upper limit of 0.23 eV for the summed neutrino mass. Our results are in excellent agreement with big bang nucleosynthesis and the standard value of N_eff=3.046. We find no evidence for dynamical dark energy. Despite the success of the standard LCDM model, this cosmology does not provide a good fit to the CMB power spectrum at low multipoles, as noted previously by the WMAP team. While not of decisive significance, this is an anomaly in an otherwise self-consistent analysis of the Planck temperature data.

We present far-infrared (FIR) photometry at 150 and 205 micron(s) of eight low-redshift starburst galaxies obtained with the Infrared Space Observatory (ISO) ISOPHOT. Five of the eight galaxies are detected in both wave bands, and these data are used, in conjunction with IRAS archival photometry, to model the dust emission at lambda approximately greater than 40 microns. The FIR spectral energy distributions (SEDs) are best fitted by a combination of two modified Planck functions, with T approx. 40 - 55 K (warm dust) and T approx. 20-23 K (cool dust) and with a dust emissivity index epsilon = 2. The cool dust can be a major contributor to the FIR emission of starburst galaxies, representing up to 60% of the total flux. This component is heated not only by the general interstellar radiation field, but also by the starburst itself. The cool dust mass is up to approx. 150 times larger than the warm dust mass, bringing the gas-to-dust ratios of the starbursts in our sample close to Milky Way values, once resealed for the appropriate metallicity. The ratio between the total dust FIR emission in the range 1-1000 microns and the IRAS FIR emission in the range 40 - 120 microns is approx. 1.75, with small variations from galaxy to galaxy. This ratio is about 40% larger than previously inferred from data at millimeter wavelengths. Although the galaxies in our sample are generally classified as "UV bright," for four of them the UV energy emerging shortward of 0.2 microns is less than 15% of the FIR energy. On average, about 30% of the bolometric flux is coming out in the UV-to-near-IR wavelength range; the rest is emitted in the FIR. Energy balance calculations show that the FIR emission predicted by the dust reddening of the UV-to-near-IR stellar emission is within a factor of approx. 2 of the observed value in individual galaxies and within 20% when averaged over a large sample. If our sample of local starbursts is representative of high-redshift (z approx. greater than 1), UV - bright star-forming galaxies, these galaxies' FIR emission will be generally undetected in submillimeter surveys, unless: (1) their bolometric luminosity is comparable to or larger than that of ultraluminous FIR galaxies and (2) their FIR SED contains a cool dust component.

/pdf/the-dust-content-and-opacity-of-actively-star-forming-3rvrlvzkyq.pdf

The Dust Content and Opacity of Actively Star-Forming Galaxies

Dusty hyperluminous galaxies in the early universe provide unique environments for studying the role of massive stars in the formation and destruction of dust. At redshifts above 6, when the universe was less than 1 Gyr old, dust could have only condensed in the explosive ejecta of Type II supernovae (SNe), since most of the progenitors of the AGB stars, the major alternative source of interstellar dust, did not have time to evolve off the main sequence. We present analytical models for the evolution of the gas, dust, and metals in high redshift galaxies, with a special application to SDSS J1148+5251, a hyperluminous quasar at z = 6.4. We show that an average SN must condense at least 1 Msun of dust to account for the mass of dust in this object, when grain destruction by supernova remnants (SNRs) is taken into account. This required yield is in excess of ~0.05 Msun, the largest mass of dust inferred from infrared observations of Cas A. If the yield of Cas A is typical, then other processes, such as accretion onto preexisting grains in molecular clouds is needed to produce the mass of dust in J1148+5251. For such process to be effective, SNR must significantly increase, presumably by non-evaporative grain-grain collisions during the late stages of their evolution, the number of nucleation centers onto which refractory elements can condense in molecular clouds.

The origin and evolution of dust in high-redshift galaxies

The research for this project was supported by NASA 16-ATP16-0004. We thank the referee for pointing out the detection of the [C II] emission in Arp 220 East, and for suggesting the discussion on the applicability of this line as a SFR tracer in heavily obscured systems. We also acknowledge helpful discussions with Jackie Fischer, Jason Glenn, and Nick Scoville.

Dust Masses, Compositions, and Luminosities in the Nuclear Disks and the Diffuse Circumnuclear Medium of Arp 220

The presence of an infrared (IR) continuum in the spectrum of SN 1987A was discovered in the early stages of its evolution by various photometric and spectral observations in the 5 — 20 µm wavelength regime [1,2]. This IR continuum component represented an IR excess over that expected from a Planck spectrum fitted to the SN emission at optical wavelengths. Subsequent observations with the Kuiper Airborne Observatory (KAO) showed this continuum to be quite flat out to 100 μm with a value of ≈ 5 – 10 Jy around day 265 after core collapse [3,4]. AITKEN et al. [5], were first to interpret this IR excess as free-free emission from the expanding ionized gas. As expected from an expanding cooling gas, the emission from this component has been continuously declining since its first discovery, its 10 μm emissivity reaching a minimum of about 1 Jy around day 400 [6]. Subsequent observations by ROCHE et al. [6] revealed an upturn in the 10 µm continuum light curve of the SN, indicating either the delayed arrival of thermal emission from pre-existing circumstellar/interstellar dust heated by the UV-visual output of the SN (i.e. an IR echo), or the onset of dust formation in the SN ejecta.

Infrared Emission from SN 1987A: Light Echoes or Dust Formation?

We present the main design features for the GISMO-2 bolometer camera, which we build for background-limited operation at the IRAM 30 m telescope on Pico Veleta, Spain. GISMO-2 will operate simultaneously in the 1 and 2 mm atmospherical windows. The 1 mm channel uses a 32×40 TES-based backshort under grid (BUG) bolometer array, the 2 mm channel operates with a 16×16 BUG array. The camera utilizes almost the entire full field of view provided by the telescope. The optical design of GISMO-2 was strongly influenced by our experience with the GISMO 2 mm bolometer camera, which is successfully operating at the 30 m telescope. GISMO is accessible to the astronomical community through the regular IRAM call for proposals.

/pdf/design-and-expected-performance-of-gismo-2-a-two-color-45h0al3x42.pdf

Design and Expected Performance of GISMO-2, a Two Color Millimeter Camera for the IRAM 30 m Telescope

SAFIRE is a versatile imaging Fabry-Perot spectrograph covering 145 to 655 microns, with spectral resolving powers ranging over 5 - 10,000. Selected as a `PI' instrument for the airborne Stratospheric Observatory for Infrared Astronomy (SOFIA). SAFIRE will apply 2D pop-up bolometer arrays to provide background-limited imaging spectrometry. Superconducting transition edge bolometers and SQUID multiplexers are being developed for these detectors. SAFIRE is expected to be a `First Light' instrument, usable during the initial SOFIA operations. Although a PI instrument rather than a `Facility Class' science instrument, it will be highly integrated with the standard SOFIA planning, observation, and data analysis tools.

/pdf/submillimeter-and-far-infrared-experiment-safire-a-pi-class-c4ds90z016.pdf

Eli Dwek

Papers

The origin and evolution of dust in high-redshift galaxies

Dust Masses, Compositions, and Luminosities in the Nuclear Disks and the Diffuse Circumnuclear Medium of Arp 220

Infrared Emission from SN 1987A: Light Echoes or Dust Formation?

Design and Expected Performance of GISMO-2, a Two Color Millimeter Camera for the IRAM 30 m Telescope

Submillimeter and far-infrared experiment (SAFIRE): a PI class instrument for SOFIA