Direct Estimate of Cirrus Noise in Herschel Hi-GAL Images

In Herschel images of the Galactic plane and many star forming regions, a major factor limiting our ability to extract faint compact sources is cirrus confusion noise, operationally defined as the "statistical error to be expected in photometric measurements due to confusion in a background of fluctuating surface brightness". The histogram of the flux densities of extracted sources shows a distinctive faint-end cutoff below which the catalog suffers from incompleteness and the flux densities become unreliable. This empirical cutoff should be closely related to the estimated cirrus noise and we show that this is the case. We compute the cirrus noise directly, both on Herschel images from which the bright sources have been removed and on simulated images of cirrus with statistically similar fluctuations. We connect these direct estimates with those from power spectrum analysis, which has been used extensively to predict the cirrus noise and provides insight into how it depends on various statistical properties and photometric operational parameters. We report multi-wavelength power spectra of diffuse Galactic dust emission from Hi-GAL observations at 70 to 500 microns within Galactic plane fields at l= 30 degrees and l= 59 degrees. We find that the exponent of the power spectrum is about -3. At 250 microns, the amplitude of the power spectrum increases roughly as the square of the median brightness of the map and so the expected cirrus noise scales linearly with the median brightness. Generally, the confusion noise will be a worse problem at longer wavelengths, because of the combination of lower angular resolution and the rising power spectrum of cirrus toward lower spatial frequencies, but the photometric signal to noise will also depend on the relative spectral energy distribution of the source compared to the cirrus.Comment: 4 pages (in journal), 3 figures, Astronomy and Astrophysics, accepted for publication 13 May 201

An Imprint of Molecular Cloud Magnetization in the Morphology of the Dust Polarized Emission

We describe a morphological imprint of magnetization found when considering the relative orientation of the magnetic field direction with respect to the density structures in simulated turbulent molecular clouds. This imprint was found using the Histogram of Relative Orientations (HRO): a new technique that utilizes the gradient to characterize the directionality of density and column density structures on multiple scales. We present results of the HRO analysis in three models of molecular clouds in which the initial magnetic field strength is varied, but an identical initial turbulent velocity field is introduced, which subsequently decays. The HRO analysis was applied to the simulated data cubes and mock-observations of the simulations produced by integrating the data cube along particular lines of sight. In the 3D analysis we describe the relative orientation of the magnetic field $\mathbf{B}$ with respect to the density structures, showing that: 1.The magnetic field shows a preferential orientation parallel to most of the density structures in the three simulated cubes. 2.The relative orientation changes from parallel to perpendicular in regions with density over a critical density $n_{T}$ in the highest magnetization case. 3.The change of relative orientation is largest for the highest magnetization and decreases in lower magnetization cases. This change in the relative orientation is also present in the projected maps. In conjunction with simulations HROs can be used to establish a link between the observed morphology in polarization maps and the physics included in simulations of molecular clouds.Comment: (16 pages, 11 figures, submitted to ApJ 05MAR2013, accepted 07JUL2013

Structure formation in a colliding flow: The Herschel view of the Draco nebula

The Draco nebula is a high Galactic latitude interstellar cloud likely to have been formed by the collision of a Galactic halo cloud entering the disk of the Milky Way. Such conditions are ideal to study the formation of cold and dense gas in colliding flows of warm gas. We present Herschel-SPIRE observations that reveal the fragmented structure of the interface between the infalling cloud and the Galactic layer. This front is characterized by a Rayleigh-Taylor instability structure. From the determination of the typical length of the periodic structure (2.2 pc) we estimated the gas kinematic viscosity and the turbulence dissipation scale (0.1 pc) that is compatible with that expected if ambipolar diffusion is the main mechanism of energy dissipation in the WNM. The small-scale structures of the nebula are typical of that seen in some molecular clouds. The gas density has a log-normal distribution with an average value of $10^3$ cm$^{-3}$. The size of the structures is 0.1-0.2 pc but this estimate is limited by the resolution of the observations. The mass ranges from 0.2 to 20 M$_{\odot}$ and the distribution of the more massive clumps follows a power law $dN/d\log(M) \sim M^{-1.4}$. We identify a mass-size relation with the same exponent as that found in GMCs ($M\sim L^{2.3}$) but only 15% of the mass of the cloud is in gravitationally bound structures. We conclude that the increase of pressure in the collision is strong enough to trigger the WNM-CNM transition caused by the interplay between turbulence and thermal instability as self-gravity is not dominating the dynamics.Comment: 16 pages, A&A, in pres

Magnetic field morphology in nearby molecular clouds as revealed by starlight and submillimetre polarization

Within four nearby (d < 160 pc) molecular clouds, we statistically evaluate the structure of the interstellar magnetic field, projected on the plane of the sky and integrated along the line of sight, as inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz and from the optical and NIR polarization of background starlight. We compare the dispersion of the field orientation directly in vicinities with an area equivalent to that subtended by the Planck effective beam at 353 GHz (10') and using the second-order structure functions of the field orientation angles. We find that the average dispersion of the starlight-inferred field orientations within 10'-diameter vicinities is less than 20 deg, and that at these scales the mean field orientation is on average within 5 deg of that inferred from the submillimetre polarization observations in the considered regions. We also find that the dispersion of starlight polarization orientations and the polarization fractions within these vicinities are well reproduced by a Gaussian model of the turbulent structure of the magnetic field, in agreement with the findings reported by the Planck collaboration at scales greater than 10' and for comparable column densities. At scales greater than 10', we find differences of up to 14.7 deg between the second-order structure functions obtained from starlight and submillimetre polarization observations in the same positions in the plane of the sky, but comparison with a Gaussian model of the turbulent structure of the magnetic field indicates that these differences are small and are consistent with the difference in angular resolution between both techniques.Comment: 15 pages, 10 figures, submitted to A&

BLAST05: Power Spectra of Bright Galactic Cirrus at Submillimeter Wavelengths

We report multi-wavelength power spectra of diffuse Galactic dust emission from BLAST observations at 250, 350, and 500 microns in Galactic Plane fields in Cygnus X and Aquila. These submillimeter power spectra statistically quantify the self-similar structure observable over a broad range of scales and can be used to assess the cirrus noise which limits the detection of faint point sources. The advent of submillimeter surveys with the Herschel Space Observatory makes the wavelength dependence a matter of interest. We show that the observed relative amplitudes of the power spectra can be related through a spectral energy distribution (SED). Fitting a simple modified black body to this SED, we find the dust temperature in Cygnus X to be 19.9 +/- 1.3 K and in the Aquila region 16.9 +/- 0.7 K. Our empirical estimates provide important new insight into the substantial cirrus noise that will be encountered in forthcoming observations.Comment: Submitted to the Astrophysical Journal. Maps and other data are available at http://blastexperiment.info

Component separation methods for the Planck mission

The Planck satellite will map the full sky at nine frequencies from 30 to 857 GHz. The CMB intensity and polarization that are its prime targets are contaminated by foreground emission. The goal of this paper is to compare proposed methods for separating CMB from foregrounds based on their different spectral and spatial characteristics, and to separate the foregrounds into components of different physical origin. A component separation challenge has been organized, based on a set of realistically complex simulations of sky emission. Several methods including those based on internal template subtraction, maximum entropy method, parametric method, spatial and harmonic cross correlation methods, and independent component analysis have been tested. Different methods proved to be effective in cleaning the CMB maps from foreground contamination, in reconstructing maps of diffuse Galactic emissions, and in detecting point sources and thermal Sunyaev-Zeldovich signals. The power spectrum of the residuals is, on the largest scales, four orders of magnitude lower than that of the input Galaxy power spectrum at the foreground minimum. The CMB power spectrum was accurately recovered up to the sixth acoustic peak. The point source detection limit reaches 100 mJy, and about 2300 clusters are detected via the thermal SZ effect on two thirds of the sky. We have found that no single method performs best for all scientific objectives. We foresee that the final component separation pipeline for Planck will involve a combination of methods and iterations between processing steps targeted at different objectives such as diffuse component separation, spectral estimation and compact source extraction.Comment: Matches version accepted by A&A. A version with high resolution figures is available at http://people.sissa.it/~leach/compsepcomp.pd

Radio to infrared spectra of late-type galaxies with Planck and WMAP data

We use the Planck Early Release Compact Source Catalogue combined with WMAP and other archival measurements to construct continuum spectra of three nearby dusty star-forming galaxies: Messier 82, NGC 253 and NGC 4945. We carry out a least-squares fit to the spectra using a combination of simple synchrotron, free-free and thermal dust models, and look for evidence of anomalous microwave emission (AME). We find that the radio spectra of all three galaxies are consistent with steep spectrum synchrotron emission, with a significant amount of free-free emission required to explain the Planck and WMAP data points in the frequency range 30-150 GHz. This brings the star-formation rate based on free-free emission into better agreement with that from the non-thermal emission. We place limits on the presence of AME in these galaxies, finding that it is lower than expectations based on the ratio of far infrared to AME from the Galaxy. Nevertheless, the shape of the spectrum of NGC 4945 hints at the presence of AME with a peak around 30 GHz. Future Planck data will let us look more closely at these galaxies, as well as to extend the analysis to many more galaxies.Comment: 5 pages, 1 figure (6 panels), 1 table. Submitted to MNRAS letter

The pre-launch Planck Sky Model: a model of sky emission at submillimetre to centimetre wavelengths

We present the Planck Sky Model (PSM), a parametric model for the generation of all-sky, few arcminute resolution maps of sky emission at submillimetre to centimetre wavelengths, in both intensity and polarisation. Several options are implemented to model the cosmic microwave background, Galactic diffuse emission (synchrotron, free-free, thermal and spinning dust, CO lines), Galactic H-II regions, extragalactic radio sources, dusty galaxies, and thermal and kinetic Sunyaev-Zeldovich signals from clusters of galaxies. Each component is simulated by means of educated interpolations/extrapolations of data sets available at the time of the launch of the Planck mission, complemented by state-of-the-art models of the emission. Distinctive features of the simulations are: spatially varying spectral properties of synchrotron and dust; different spectral parameters for each point source; modeling of the clustering properties of extragalactic sources and of the power spectrum of fluctuations in the cosmic infrared background. The PSM enables the production of random realizations of the sky emission, constrained to match observational data within their uncertainties, and is implemented in a software package that is regularly updated with incoming information from observations. The model is expected to serve as a useful tool for optimizing planned microwave and sub-millimetre surveys and to test data processing and analysis pipelines. It is, in particular, used for the development and validation of data analysis pipelines within the planck collaboration. A version of the software that can be used for simulating the observations for a variety of experiments is made available on a dedicated website.Comment: 35 pages, 31 figure

Evidence for dust evolution within the Taurus Complex from Spitzer images

We present Spitzer images of the Taurus Complex (TC) and take advantage of the sensitivity and spatial resolution of the observations to characterize the diffuse IR emission across the cloud. This work highlights evidence of dust evolution within the translucent sections of the archetype reference for studies of quiescent molecular clouds. We combine Spitzer 160 um and IRAS 100 um observations to produce a dust temperature map and a far-IR dust opacity map at 5' resolution. The average dust temperature is about 14.5K with a dispersion of +/-1K across the cloud. The far-IR dust opacity is a factor 2 larger than the average value for the diffuse ISM. This opacity increase and the attenuation of the radiation field (RF) both contribute to account for the lower emission temperature of the large grains. The structure of the TC significantly changes in the mid-IR images that trace emission from PAHs and VSGs. We focus our analysis of the mid-IR emission to a range of ecliptic latitudes where the zodiacal light residuals are small. Within this cloud area, there are no 8 and 24 um counterparts to the brightest 160 um emission features. Conversely, the 8 and 24 um images reveal filamentary structure that is strikingly inconspicuous in the 160 um and extinction maps. The IR colors vary over sub-parsec distances across this filamentary structure. We compare the observed colors with model calculations quantifying the impact of the RF intensity and the abundance of stochastically heated particles on the dust SED. To match the range of observed colors, we have to invoke variations by a factor of a few of both the interstellar RF and the abundance of PAHs and VSGs. We conclude that within this filamentary structure a significant fraction of the dust mass cycles in and out the small size end of the dust size distribution.Comment: 43 pages, 13 figures, accepted for publication in Ap

ISOGAL: A deep survey of the obscured inner Milky Way with ISO at 7 and 15 micron and with DENIS in the near-infrared

The ISOGAL project is an infrared survey of specific regions sampling the Galactic Plane selected to provide information on Galactic structure,stellar populations,stellar mass-loss and the recent star formation history of the inner disk and Bulge of the Galaxy. ISOGAL combines 7 and 15 micron ISOCAM observations - with a resolution of 6'' at worst - with DENIS IJKs data to determine the nature of the sources and theinterstellar extinction. We have observed about 16 square degrees with a sensitivity approaching 10-20mJy, detecting ~10^5 sources,mostly AGB stars,red giants and young stars. The main features of the ISOGAL survey and the observations are summarized in this paper,together with a brief discussion of data processing and quality. The primary ISOGAL products are described briefly (a full description is given in Schuller et al. 2003, astro-ph/0304309): viz. the images and theISOGAL-DENIS five-wavelength point source catalogue. The main scientific results already derived or in progress are summarized. These include astrometrically calibrated 7 and 15um images,determining structures of resolved sources; identification and properties of interstellar dark clouds; quantification of the infrared extinction law and source dereddening; analysis of red giant and (especially) AGB stellar populations in the central Bulge,determining luminosity,presence of circumstellar dust and mass--loss rate,and source classification,supplemented in some cases by ISO/CVF spectroscopy; detection of young stellar objects of diverse types,especially in the inner Bulge with information about the present and recent star formation rate; identification of foreground sources with mid-IR excess. These results are the subject of about 25 refereed papers published or in preparation.Comment: A&A in press. 19 pages,10 Ps figures; problems with figures fixe