Source author record

P. Andreani

P. Andreani 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

68works
7topics
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

68 published item(s)

preprint2022arXiv

The Bright Extragalactic ALMA Redshift Survey (BEARS) I: redshifts of bright gravitationally-lensed galaxies from the Herschel ATLAS

We present spectroscopic measurements for 71 galaxies associated with 62 of the brightest high-redshift submillimeter sources from the Southern fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), while targeting 85 sources which resolved into 142. We have obtained robust redshift measurements for all sources using the 12-m Array and an efficient tuning of ALMA to optimise its use as a redshift hunter, with 73 per cent of the sources having a robust redshift identification. Nine of these redshift identifications also rely on observations from the Atacama Compact Array. The spectroscopic redshifts span a range $1.41<z<4.53$ with a mean value of 2.75, and the CO emission line full-width at half-maxima range between $\rm 110\,km\,s^{-1} < FWHM < 1290\,km\,s^{-1}$ with a mean value of $\sim$ 500kms$^{-1}$, in line with other high-$z$ samples. The derived CO(1-0) luminosity is significantly elevated relative to line-width to CO(1-0) luminosity scaling relation, which is suggestive of lensing magnification across our sources. In fact, the distribution of magnification factors inferred from the CO equivalent widths is consistent with expectations from galaxy-galaxy lensing models, though there is a hint of an excess at large magnifications that may be attributable to the additional lensing optical depth from galaxy groups or clusters.

preprint2022arXiv

The molecular gas properties in local Seyfert 2 galaxies

We present a multi-wavelength study of the molecular gas properties of a sample of local Seyfert 2 galaxies to assess if, and to what extent, the presence of an active galactic nucleus (AGN) can affect the Interstellar Medium (ISM) properties in a sample of 33 local Seyfert 2 galaxies. We compare the molecular gas content (MH2), derived from new and archival low-J CO line measurements of a sample of AGN and a control sample of star-forming galaxies (SFGs). Both the AGN and the control sample are characterised in terms of host-galaxy properties (e.g., stellar and dust masses, Mstar and Mdust, respectively; and star formation rate, SFR). We also investigate the effect of AGN activity onto the emission of polycyclic aromatic hydrocarbon (PAH) molecules in the mid-infrared (MIR), a waveband where the dust-reprocessed emission from the obscured AGN contributes the most. The AGN hosted in less massive galaxies (i.e., log(Mstar/Msun)<10.5; log(Mdust/Msun)<7.5) show larger molecular gas content with respect to SFGs matched in stellar and dust mass. When comparing their depletion time(tdep~MH2/SFR), AGN show tdep~0.3-1.0 Gyr, similar to those observed in the control sample of SFGs. Seyfert 2 galaxies show fainter PAH luminosity at increasingly larger dominance of the nuclear activity in the MIR. We find no clear evidence for a systematic reduction of the molecular gas reservoir at galactic scale in Seyfert galaxies with respect to SFGs. This is in agreement with recent studies showing that molecular gas content only is reduced in regions of sub-kpc size, where the emission from the accreting supermassive black hole dominates. Nonetheless, we show that the impact of AGN activity on the ISM is clearly visible as suppression of the PAH luminosity.

preprint2022arXiv

The star-formation rates of QSOs

We examine the far-IR properties of a sample of 5391 optically selected QSOs in the 0.5<z<2.65 redshift range down to log[nuLnu,2500 (erg/s)]>44.7, using SPIRE data from Herschel-ATLAS. We split the sample in a grid of 74 luminosity-redshift bins and compute the average optical-infrared spectral energy distribution (SED) in each bin. By normalising an intrinsic AGN template to the AGN optical power (at 5100A) we decompose the total infrared emission (L_IR; 8-1000um) into an AGN (L_IR,AGN) and star-forming component (L_IR,SF). We find that the AGN contribution to L_IR increases as a function of AGN power which manifests as a reduction of the `far-IR bump' in the average QSO SEDs. We note that L_IR,SF does not correlate with AGN power; the mean star formation rates (SFRs) of AGN host galaxies are a function of redshift only and they range from ~6 Msun/yr at z~0 to a plateau of <200 Msun/yr at z~2.6. Our results indicate that the accuracy of far-IR emission as a proxy for SFR decreases with increasing AGN luminosity. We show that, at any given redshift, observed trends between infrared luminosity (whether monochromatic or total) and AGN power (in the optical or X-rays) can be explained by a simple model which is the sum of two components: (A) the infrared emission from star-formation, uncorrelated with AGN power and (B) the infrared emission from AGN, directly proportional to AGN power in the optical or X-rays.

preprint2021arXiv

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).XII. Ionised gas emission in the inner regions of lenticular galaxies

As part of the VESTIGE survey, a blind narrow-band Ha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the CFHT, we discovered 8 massive lenticular galaxies with prominent ionised gas emission features in their inner (few kpc) regions. These features are either ionised gas filaments similar to those observed in cooling flows (2 gal), or thin discs with sizes 0.7<R(Ha)<2.0 kpc (6 gal), thus significantly smaller than those of the stellar disc. These discs have morphological properties similar to those of the dust seen in absorption in high-resolution HST images. Using a unique set of multifrequency data we show that while the gas located within these inner discs is photoionised by young stars, signaling ongoing star formation, the gas in the filamentary structures is shock-ionised. These discs have a star formation surface brightness similar to those observed in late-type galaxies. Because of their reduced size, however, these lenticular galaxies are located below the main sequence of unperturbed or cluster star-forming systems. By comparing the dust masses measured from absorption maps in optical images, from the Balmer decrement, or estimated by fitting the UV-to-far-IR spectral energy distribution of the target galaxies, we confirm that those derived from optical attenuation maps are heavily underestimated because of geometrical effects due to the relative distribution of the absorbing dust and the emitting stars. We have also shown that these galaxies have gas-to-dust ratios of G/D~80, and that the star formation within these discs follows the Schmidt relation, albeit with an efficiency reduced by a factor of ~ 2.5. Using our unique set of multifrequency data, we discuss the possible origin of the ionised gas in these objects, which suggests multiple and complex formation scenarios for massive lenticular galaxies in clusters.

preprint2020arXiv

A critical comparison of the Lomb-Scargle and the classical periodograms

The detection of signals hidden in noise is one of the oldest and common problems in astronomy. Various solutions have been proposed in the past such as the parametric approaches based on the least-squares fit of theoretical templates or the non-parametric techniques as the phase-folding method. Most of them, however, are suited only for signals with specific time evolution. For generic signals the spectral approach based on the periodogram is potentially the most effective. In astronomy the main problem in working with the periodogram is that often the sampling of the signals is irregular. This complicates its efficient computation (the fast Fourier transform cannot be directly used) but overall the determination of its statistical characteristics. The Lomb-Scargle periodogram (LSP) provides a solution to this last important issue, but its main drawback is the assumption of a very specific model of the datawhich is not correct for most of the practical applications. These issues are not always considered in literature with theoretical and practical consequences of no easy solution. Moreover, apart from pathological samplings, it is common believe that the LSP and the classical periodogram (CP) usually provide almost identical results. In general, this is true but here it is shown that there are situations where the LSP is less effective than the CP in the detection of signals in noise. There are no compelling reasons, therefore, to use the LSP instead of the CP which is directly connected to the correlation function of the observed signal with the sinusoidal functions at the various frequencies of interest.

preprint2020arXiv

A Virgo Environmental Survey Tracing Ionised Gas Emission. VESTIGE VIII. Bridging the cluster-ICM-galaxy evolution at small scales

We measure FIR emission from tails of stripped dust following the ionised and atomic gas components in galaxies undergoing ram pressure stripping. We study the dust-to-gas relative distribution and mass ratio in the stripped interstellar medium and relate them to those of the intra-cluster medium, thus linking the cluster-ICM-galaxy evolution at small-scales. The galaxy sample consists of three Scd Virgo galaxies with stellar masses in the range $10^9\lesssim \mathrm{M_{*}} \lesssim 10^{10}\, \mathrm{M_{\odot}}$, and within 1 Mpc from the cluster centre, namely NGC 4330, NGC 4522, and NGC 4654. Through the analysis of VESTIGE H$α$, $Herschel$ SPIRE far-infrared, and VIVA HI data, we trace the spatial distribution of the tails and infer the dust and gas masses from the measured far-infrared 250 $μ$m and HI flux densities. Dust-to-gas mass ratios (DGRs) in the tails are analysed as a function of the galaxy mass, metallicity, and dust temperature. Along the stripped component, the dust distribution closely follows the HI and H$α$ emitting gas, all extending beyond the optical disc. In these regions, the DGRs are $2.0\pm0.6\times10^{-3}$, $0.7\pm0.1\times10^{-3}$, and $0.4\pm0.03\times10^{-3}$, for NGC 4330, NGC 4522, and NGC 4654, respectively, i.e. up to a factor of 15 less than the values measured in the main body of nearby galaxies. We also find a negative trend in the DGR as a function of the metallicity that can be explained in terms of a dust component more centrally concentrated in more metal-rich systems. Together with the finding that the stripped dust is cold, $T_{d} \lesssim 25\, K$, our results support an outside-in stripping scenario of the galaxy interstellar medium. This study shows that ram pressure stripping is a key mechanism in the building up of the Virgo intra-cluster component injecting dust grains into the ICM, thus contributing to its metal enrichment.

preprint2020arXiv

Full of Orions: a 200-pc mapping of the interstellar medium in the redshift-3 lensed dusty star-forming galaxy SDP.81

We present a sub-kpc resolved study of the interstellar medium properties in SDP.81, a z=3.042 strongly gravitationally lensed dusty star-forming galaxy, based on high-resolution, multi-band ALMA observations of the FIR continuum, CO ladder and the [CII] line. Using a visibility-plane lens modelling code, we achieve a median source-plane resolution of ~200 pc. We use photon-dominated region (PDR) models to infer the physical conditions - far-UV field strength, density, and PDR surface temperature - of the star-forming gas on 200-pc scales, finding a FUV field strength of ~10^3-10^4 G0, gas density of ~10^5 cm^-3 and cloud surface temperatures up to 1500 K, similar to those in the Orion Trapezium region. The [CII] emission is significantly more extended than that FIR continuum: ~50 per cent of [CII] emission arises outside the FIR-bright region. The resolved [CII]/FIR ratio varies by almost 2 dex across the source, down to ~2x10^-4 in the star-forming clumps. The observed [CII]/FIR deficit trend is consistent with thermal saturation of the C+ fine-structure level occupancy at high gas temperatures. We make the source-plane reconstructions of all emission lines and continuum data publicly available.

preprint2020arXiv

SCUBA-2 overdensities associated with candidate protoclusters selected from Planck data

We measure the 850-$μ$m source densities of 46 candidate protoclusters selected from the Planck High-z catalogue (PHz) and the Planck Catalogue of Compact Sources (PCCS) that were followed up with Herschel-SPIRE and SCUBA-2. This paper aims to search for overdensities of 850-$μ$m sources in order to select the fields that are most likely to be genuine protoclusters. Of the 46 candidate protoclusters, 25 have significant overdensities ($>$5 times the field counts), 11 have intermediate overdensities (3--5 times the field counts) and 10 have no overdensity ($<$3 times the field counts) of 850-$μ$m sources. We find that the enhanced number densities are unlikely to be the result of sample variance. Compared with the number counts of another sample selected from Planck's compact source catalogues, this [PHz+PCCS]-selected sample has a higher fraction of candidate protoclusters with significant overdensities, though both samples show overdensities of 850-$μ$m sources above intermediate level. Based on the estimated star-formation rate densities (SFRDs), we suggest that both samples can efficiently select protoclusters with starbursting galaxies near the redshift at which the global field SFRD peaks ($2 < z < 3$). Based on the confirmation of overdensities found here, future follow-up observations on other PHz targets may greatly increase the number of genuine DSFG-rich clusters/protoclusters.

preprint2020arXiv

The GALEX Ultraviolet Virgo Cluster Survey (GUViCS) VIII. Diffuse dust in the Virgo intra-cluster space

We present the first detection of diffuse dust in the intra-cluster medium of the Virgo cluster out to $\sim$0.4 virial radii, and study the radial variation of its properties on a radial scale of the virial radius. Analysing near-UV - $i$ colours for a sample of $\sim12000$ background galaxies with redshifts $0.02 < z < 0.8$, we find significant colour reddening and relate it to variation in $E(B-V)$ values. The $E(B-V)$ mean profile shows a dust component characterised by an average reddening $E(B-V)\sim0.042 \pm 0.004$ mag within 1.5 degrees ($\sim0.3\, r_{vir}$) from the cluster centre. Assuming a Large Magellanic Cloud extinction law, we derive an average visual extinction $A_{V} = 0.14\pm 0.01$ for a total dust mass, $M_{d} = 2.5\pm0.2\times10^{9}M_{\odot}$, hence a dust-to-gas mass ratio $M_{d}/M_{g} = 3.0\pm 0.3 \times 10^{-4}$. Based on the upper limits on the flux density $\mathrm{I_{250μm} = 0.1\, MJy sr^{-1}} $ derived from $Herschel$ data, we estimate an upper limit for the dust temperature of $T_{d} \sim 10\, K$. However, similar densities can be obtained with dust at higher temperatures with lower emissivities. The Virgo cluster has diffuse dust in its intra-cluster medium characterised by different physical properties as those characterising the Milky Way dust. The diffuse dust in Virgo is transported into the cluster space through similar phenomena (stripping) as those building up the optical intra-cluster light, and it constitutes an additional cooling agent of the cluster gas.

preprint2019arXiv

A CO molecular gas wind 340 pc away from the Seyfert 2 nucleus in ESO 420-G13 probes an elusive radio jet

A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point at $440\, \rm{pc}$ where the jet hits a dense molecular core and shatters, dispersing the molecular gas into several clumps and filaments within the expansion cone. We also trace the jet in ionised gas within the inner $\lesssim 340\, \rm{pc}$ using the [NeII]$_{\rm 12.8 μm}$ line emission, where the molecular gas follows a circular rotation pattern. The wind outflow carries a mass of $\sim 8 \times 10^6\, \rm{M_\odot}$ at an average wind projected speed of $\sim 160\, \rm{km\,s^{-1}}$, which implies a mass outflow rate of $\sim 14\, \rm{M_\odot\,yr^{-1}}$. Based on the structure of the outflow and the budget of energy and momentum, we discard radiation pressure from the active nucleus, star formation, and supernovae as possible launching mechanisms. ESO 420-G13 is the second case after NGC 1377 where the presence of a previously unknown jet is revealed due to its interaction with the interstellar medium, suggesting that unknown jets in feeble radio nuclei might be more common than expected. Two possible jet-cloud configurations are discussed to explain the presence of an outflow at such distance from the AGN. The outflowing gas will likely not escape, thus a delay in the star formation rather than quenching is expected from this interaction, while the feedback effect would be confined within the central few hundred parsecs of the galaxy.

preprint2016arXiv

Far-infrared line spectra of active galaxies from the Herschel/PACS Spectrometer: the complete database

We present a coherent database of spectroscopic observations of far-IR fine-structure lines from the Herschel/PACS archive for a sample of 170 local AGN, plus a comparison sample of 20 starburst galaxies and 43 dwarf galaxies. Published Spitzer/IRS and Herschel/SPIRE line fluxes are included to extend our database to the full 10-600 $μm$ spectral range. The observations are compared to a set of CLOUDY photoionisation models to estimate the above physical quantities through different diagnostic diagrams. We confirm the presence of a stratification of gas density in the emission regions of the galaxies, which increases with the ionisation potential of the emission lines. The new [OIV]25.9$μm$/[OIII]88$μm$ vs [NeIII]15.6$μm$/[NeII]12.8$μm$ diagram is proposed as the best diagnostic to separate: $i)$ AGN activity from any kind of star formation; and $ii)$ low-metallicity dwarf galaxies from starburst galaxies. Current stellar atmosphere models fail to reproduce the observed [OIV]25.9$μm$/[OIII]88$μm$ ratios, which are much higher when compared to the predicted values. Finally, the ([NeIII]15.6$μm$ + [NeII]12.8$μm$)/([SIV]10.5$μm$ + [SIII]18.7$μm$) ratio is proposed as a promising metallicity tracer to be used in obscured objects, where optical lines fail to accurately measure the metallicity. The diagnostic power of mid- to far-infrared spectroscopy shown here for local galaxies will be of crucial importance to study galaxy evolution during the dust-obscured phase at the peak of the star formation and black-hole accretion activity ($1 < z < 4$). This study will be addressed by future deep spectroscopic surveys with present and forthcoming facilities such as JWST, ALMA, and SPICA.

preprint2016arXiv

Tracing black hole accretion with SED decomposition and IR lines: from local galaxies to the high-z Universe

We present new estimates of AGN accretion and star-formation luminosity in galaxies obtained for the local 12-$μ$m sample of Seyfert galaxies (12MGS), by performing a detailed broad-band spectral energy distribution (SED) decomposition including the emission of stars, dust heated by star formation and a possible AGN dusty torus. Thanks to the availability of data from the X-rays to the sub-millimetre, we constrain and test the contribution of the stellar, AGN and star-formation components to the SEDs. The availability of Spitzer-IRS low resolution mid-infrared (mid-IR) spectra is crucial to constrain the dusty torus component at its peak wavelengths. The results of SED-fitting are also tested against the available information in other bands: the reconstructed AGN bolometric luminosity is compared to those derived from X-rays and from the high excitation IR lines tracing AGN activity like [Ne V] and [O IV]. The IR luminosity due to star-formation (SF) and the intrinsic AGN bolometric luminosity are shown to be strongly related to the IR line luminosity. Variations of these relations with different AGN fractions are investigated, showing that the relation dispersions are mainly due to different AGN relative contribution within the galaxy. Extrapolating these local relations between line and SF or AGN luminosities to higher redshifts, by means of recent Herschel galaxy evolution results, we then obtain mid- and far-IR line luminosity functions useful to estimate how many star-forming galaxies and AGN we expect to detect in the different lines at different redshifts and luminosities with future IR facilities (e.g., JWST, SPICA).

preprint2015arXiv

ALMA imaging of SDP.81 - I. A pixelated reconstruction of the far-infrared continuum emission

We present a sub-50 pc-scale analysis of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large Millimetre/submillimetre Array (ALMA) science verification data. We model both the mass distribution of the gravitational lensing galaxy and the pixelated surface brightness distribution of the background source using a novel Bayesian technique that fits the data directly in visibility space. We find the 1 and 1.3 mm dust emission to be magnified by a factor of u_tot = 17.6+/-0.4, giving an intrinsic total star-formation rate of 315+/-60 M_sol/yr and a dust mass of 6.4+/-1.5*10^8 M_sol. The reconstructed dust emission is found to be non-uniform, but composed of multiple regions that are heated by both diffuse and strongly clumped star-formation. The highest surface brightness region is a ~1.9*0.7 kpc disk-like structure, whose small extent is consistent with a potential size-bias in gravitationally lensed starbursts. Although surrounded by extended star formation, with a density of 20-30+/-10 M_sol/yr/kpc^2, the disk contains three compact regions with densities that peak between 120-190+/-20 M_sol/yr/kpc^2. Such star-formation rate densities are below what is expected for Eddington-limited star-formation by a radiation pressure supported starburst. There is also a tentative variation in the spectral slope of the different star-forming regions, which is likely due to a change in the dust temperature and/or opacity across the source.

preprint2015arXiv

ALMA imaging of SDP.81 - II. A pixelated reconstruction of the CO emission lines

We present a sub-100 pc-scale analysis of the CO molecular gas emission and kinematics of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large Millimetre/submillimetre Array (ALMA) science verification data and a visibility-plane lens reconstruction technique. We find clear evidence for an excitation dependent structure in the unlensed molecular gas distribution, with emission in CO (5-4) being significantly more diffuse and structured than in CO (8-7). The intrinsic line luminosity ratio is r_8-7/5-4 = 0.30 +/- 0.04, which is consistent with other low-excitation starbursts at z ~ 3. An analysis of the velocity fields shows evidence for a star-forming disk with multiple velocity components that is consistent with a merger/post-coalescence merger scenario, and a dynamical mass of M(< 1.56 kpc) = 1.6 +/- 0.6 x 10^10 M_sol . Source reconstructions from ALMA and the Hubble Space Telescope show that the stellar component is offset from the molecular gas and dust components. Together with Karl G. Jansky Very Large Array CO (1-0) data, they provide corroborative evidence for a complex ~2 kpc-scale starburst that is embedded within a larger ~15 kpc structure.

preprint2015arXiv

ALMA Long Baseline Observations of the Strongly Lensed Submillimeter Galaxy HATLAS J090311.6+003906 at z=3.042

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

preprint2015arXiv

ALMA Observations of Asteroid 3 Juno at 60 Kilometer Resolution

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the asteroid 3 Juno obtained with an angular resolution of 0.042 arcseconds (60 km at 1.97 AU). The data were obtained over a single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. A sequence of ten consecutive images reveals continuous changes in the asteroid's profile and apparent shape, in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models from Inversion Techniques. We measure a geometric mean diameter of 259pm4 km, in good agreement with past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature is 215pm13 K, while the median over the whole surface is 197pm15 K. With the unprecedented resolution of ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the regolith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact crater. These results demonstrate ALMA's potential to resolve thermal emission from the surface of main belt asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics.

preprint2015arXiv

An Overview of the 2014 ALMA Long Baseline Campaign

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

preprint2015arXiv

Extragalactic sources in Cosmic Microwave Background maps

We discuss the potential of a next generation space-borne CMB experiment for studies of extragalactic sources with reference to COrE+, a project submitted to ESA in response to the M4 call. We consider three possible options for the telescope size: 1m, 1.5m and 2m (although the last option is probably impractical, given the M4 boundary conditions). The proposed instrument will be far more sensitive than Planck and will have a diffraction-limited angular resolution. These properties imply that even the 1m telescope option will perform substantially better than Planck for studies of extragalactic sources. The source detection limits as a function of frequency have been estimated by means of realistic simulations. The most significant improvements over Planck results are presented for each option. COrE+ will provide much larger samples of truly local star-forming galaxies, making possible analyses of the properties of galaxies (luminosity functions, dust mass functions, star formation rate functions, dust temperature distributions, etc.) across the Hubble sequence. Even more interestingly, COrE+ will detect, at |b|> 30 deg, thousands of strongly gravitationally lensed galaxies. Such large samples are of extraordinary astrophysical and cosmological value in many fields. Moreover, COrE+ high frequency maps will be optimally suited to pick up proto-clusters of dusty galaxies, i.e. to investigate the evolution of large scale structure at larger redshifts than can be reached by other means. Thanks to its high sensitivity COrE+ will also yield a spectacular advance in the blind detection of extragalactic sources in polarization. This will open a new window for studies of radio source polarization and of the global properties of magnetic fields in star forming galaxies and of their relationships with SFRs.

preprint2015arXiv

First Results from High Angular Resolution ALMA Observations Toward the HL Tau Region

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($α$), which ranges from $α\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

preprint2015arXiv

Star formation properties of sub-mJy radio sources

We investigate the star formation properties of ~800 sources detected in one of the deepest radio surveys at 1.4 GHz. Our sample spans a wide redshift range (~0.1 - 4) and about four orders of magnitude in star formation rate (SFR). It includes both star forming galaxies (SFGs) and active galactic nuclei (AGNs), further divided into radio-quiet and radio-loud objects. We compare the SFR derived from the far infrared luminosity, as traced by Herschel, with the SFR computed from their radio emission. We find that the radio power is a good SFR tracer not only for pure SFGs but also in the host galaxies of RQ AGNs, with no significant deviation with redshift or specific SFR. Moreover, we quantify the contribution of the starburst activity in the SFGs population and the occurrence of AGNs in sources with different level of star formation. Finally we discuss the possibility of using deep radio survey as a tool to study the cosmic star formation history.

preprint2014arXiv

A Multiwavelength Consensus on the Main Sequence of Star-Forming Galaxies at z~2

We compare various star formation rate (SFR) indicators for star-forming galaxies at $1.4<z<2.5$ in the COSMOS field. The main focus is on the SFRs from the far-IR (PACS-Herschel data) with those from the ultraviolet, for galaxies selected according to the BzK criterion. FIR-selected samples lead to a vastly different slope of the SFR-stellar mass ($M_*$) relation, compared to that of the dominant main sequence population as measured from the UV, since the FIR selection picks predominantly only a minority of outliers. However, there is overall agreement between the main sequences derived with the two SFR indicators, when stacking on the PACS maps the BzK-selected galaxies. The resulting logarithmic slope of the SFR-{$M_*$} relation is $\sim0.8-0.9$, in agreement with that derived from the dust-corrected UV-luminosity. Exploiting deeper 24$μ$m-Spitzer data we have characterized a sub-sample of galaxies with reddening and SFRs poorly constrained, as they are very faint in the $B$ band. The combination of Herschel with Spitzer data have allowed us to largely break the age/reddening degeneracy for these intriguing sources, by distinguishing whether a galaxy is very red in B-z because of being heavily dust reddened, or whether because star formation has been (or is being) quenched. Finally, we have compared our SFR(UV) to the SFRs derived by stacking the radio data and to those derived from the H$α$ luminosity of a sample of star-forming galaxies at $1.4<z<1.7$. The two sets of SFRs are broadly consistent as they are with the SFRs derived from the UV and by stacking the corresponding PACS data in various mass bins.

preprint2014arXiv

Karl G. Jansky Very Large Array observations of cold dust and molecular gas in starbursting quasar host galaxies at z~4.5

We present Karl G. Jansky Very Large Array (VLA) observations of 44 GHz continuum and CO J=2-1 line emission in BR1202-0725 at z=4.7 (a starburst galaxy and quasar pair) and BRI1335-0417 at z=4.4 (also hosting a quasar). With the full 8 GHz bandwidth capabilities of the upgraded VLA, we study the (rest-frame) 250 GHz thermal dust continuum emission for the first time along with the cold molecular gas traced by the Low-J CO line emission. The measured CO J=2-1 line luminosities of BR1202-0725 are L'(CO) = (8.7+/-0.8)x10^10 K km/s pc^2 and L'(CO) = (6.0+/-0.5)x10^10 K km/s pc^2 for the submm galaxy (SMG) and quasar, which are equal to previous measurements of the CO J=5-4 line luminosities implying thermalized line emission and we estimate a combined cold molecular gas mass of ~9x10^10 Msun. In BRI1335-0417 we measure L'(CO) = (7.3+/-0.6)x10^10 K km/s pc^2. We detect continuum emission in the SMG BR1202-0725 North (S(44GHz) = 51+/-6 microJy), while the quasar is detected with S(44GHz) = 24+/-6 microJy and in BRI1335-0417 we measure S(44GHz) = 40+/-7 microJy. Combining our continuum observations with previous data at (rest-frame) far-infrared and cm-wavelengths, we fit three component models in order to estimate the star-formation rates. This spectral energy distribution fitting suggests that the dominant contribution to the observed 44~GHz continuum is thermal dust emission, while either thermal free-free or synchrotron emission contributes less than 30%.

preprint2014arXiv

Relativistic cosmology number densities in void-Lemaitre-Tolman-Bondi models

The goal of this work is to compute the number density of far-IR selected galaxies in the comoving frame and along the past lightcone of observationally constrained LTB models and to compare those results with their standard model counterparts. We derived integral number densities and differential number densities using different cosmological distance definitions in the LTB models. Then, we computed selection functions and consistency functions for the luminosity functions in the combined fields of the Herschel/PACS evolutionary probe (PEP) survey in both standard and void cosmologies, from which we derived the observed values of the above-mentioned densities. We used the Kolmogorov-Smirnov statistics to study both the evolution of the consistency functions and its connection to the evolution of the comoving density of sources. Finally, we fitted the power-law behaviour of the densities along the observer's past lightcone. The analysis of the comoving number density shows that the increased flexibility of the LTB models is not enough to fit the observed redshift evolution of the number counts, if it is specialised to a recent best-fit giant void parametrisation. The results for the power-law fits of the densities along the observer's past lightcone show general agreement across both cosmological models studied here around a slope of -2.5 $\pm$ 0.1 for the integral number density on the luminosity-distance volumes. The differential number densities show much bigger slope discrepancies. We conclude that the differential number densities on the observer's past lightcone were still rendered dependent on the cosmological model by the flux limits of the PEP survey. In addition, we show that an intrinsic evolution of the sources must be assumed to fit the comoving number-density redshift evolution in the giant void parametrisation for the LTB models used in this work.

preprint2014arXiv

The bivariate K-band-submillimetre luminosity functions of the local HRS galaxy sample

We study the relationship between the K-band and the sub-millimetre (submm) emissions of nearby galaxies by computing the bivariate K-band-submm luminosity function (BLF) of the Herschel Reference Survey (HRS), a volume-limited sample observed in submm with Hersche/SPIRE. We derive the BLF from the K-band and submm cumulative distributions using a copula method. Using the BLF allows us to derive the relationship between the luminosities on more solid statistical ground. The analysis shows that over the whole HRS sample, no statistically meaningful conclusion can be derived for any relationship between the K-band and the submm luminosity. However, a very tight relationship between these luminosities is highlighted, by restricting our analysis to late-type galaxies. The luminosity function of late-type galaxies computed in the K-band and in the submm are dependent and the dependence is caused by the link, between the stellar mass and the cold dust mass, which has been already observed

preprint2013arXiv

An approach for the detection of point-sources in very high resolution microwave maps

This paper deals with the detection problem of extragalactic point-sources in multi-frequency, microwave sky maps that will be obtainable in future cosmic microwave background radiation (CMB) experiments with instruments capable of very high spatial resolution. With spatial resolutions that can be of order of 0.1-1.0 arcsec or better, the extragalactic point-sources will appear isolated. The same holds also for the compact structures due to the Sunyaev-Zeldovich (SZ) effect (both thermal and kinetic). This situation is different from the maps obtainable with instruments as WMAP or PLANCK where, because of the smaller spatial resolution (approximately 5-30 arcmin), the point-sources and the compact structures due to the SZ effect form a uniform noisy background (the "confusion noise"). Hence, the point-source detection techniques developed in the past are based on the assumption that all the emissions that contribute to the microwave background can be modeled with homogeneous and isotropic (often Gaussian) random fields and make use of the corresponding spatial power-spectra. In the case of very high resolution observations such an assumption cannot be adopted since it still holds only for the CMB. Here, we propose an approach based on the assumption that the diffuse emissions that contribute to the microwave background can be locally approximated by two-dimensional low order polynomials. In particular, two sets of numerical techniques are presented containing two different algorithms each. The performance of the algorithms is tested with numerical experiments that mimic the physical scenario expected for high Galactic latitude observations with the Atacama Large Millimeter/Submillimeter Array (ALMA).

preprint2013arXiv

Cosmological model dependence of the galaxy luminosity function: far-infrared results in the Lemaitre-Tolman-Bondi model

This is the first paper of a series aiming at investigating galaxy formation and evolution in the giant-void class of the Lemaitre-Tolman-Bondi (LTB) models that best fits current cosmological observations. Here we investigate the Luminosity Function (LF) methodology, and how its estimates would be affected by a change on the cosmological model assumed in its computation. Are the current observational constraints on the allowed Cosmology enough to yield robust LF results? We use the far-infrared source catalogues built on the observations performed with the Herschel/PACS instrument, and selected as part of the PACS evolutionary probe (PEP) survey. Schechter profiles are obtained in redshift bins up to z approximately 4, assuming comoving volumes in both the standard model, that is, Friedmann-Lemaitre-Robertson-Walker metric with a perfect fluid energy-momentum tensor, and non-homogeneous LTB dust models, parametrized to fit the current combination of results stemming from the observations of supernovae Ia, the cosmic microwave background, and baryonic acoustic oscillations. We find that the luminosity functions computed assuming both the standard model and LTB void models show in general good agreement. However, the faint-end slope in the void models shows a significant departure from the standard model up to redshift 0.4. We demonstrate that this result is not artificially caused by the used LF estimator which turns out to be robust under the differences in matter-energy density profiles of the models. The differences found in the LF slopes at the faint end are due to variation in the luminosities of the sources, which depend on the geometrical part of the model. It follows that either the standard model is over-estimating the number density of faint sources or the void models are under-estimating it.

preprint2013arXiv

Irregular time series in astronomy and the use of the Lomb-Scargle periodogram

Detection of a signal hidden by noise within a time series is an important problem in many astronomical searches, i.e. for light curves containing the contributions of periodic/semi-periodic components due to rotating objects and all other astrophysical time-dependent phenomena. One of the most popular tools for use in such studies is the "periodogram", whose use in an astronomical context is often not trivial. The "optimal" statistical properties of the periodogram are lost in the case of irregular sampling of signals, which is a common situation in astronomical experiments. Parts of these properties are recovered by the "Lomb-Scargle" (LS) technique, but at the price of theoretical difficulties, that can make its use unclear, and of algorithms that require the development of dedicated software if a fast implementation is necessary. Such problems would be irrelevant if the LS periodogram could be used to significantly improve the results obtained by approximated but simpler techniques. In this work we show that in many astronomical applications simpler techniques provide results similar to those obtainable with the LS periodogram. The meaning of the "Nyquist frequency" is also discussed in the case of irregular sampling.

preprint2013arXiv

Lyman break and UV-selected galaxies at z ~ 1: II. PACS-100um/160um FIR detections

We report the PACS-100um/160um detections of a sample of 42 GALEX-selected and FIR-detected Lyman break galaxies (LBGs) at z ~ 1 located in the COSMOS field and analyze their ultra-violet (UV) to far-infrared (FIR) properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger, more massive, dustier, redder in the UV continuum, and UV-brighter than PACS-undetected LBGs. PACS-detected LBGs at z ~ 1 are mostly disk-like galaxies and are located over the green-valley and red sequence of the color-magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star-formation rates (SFRs) than other PACS-detected UV-selected galaxies within their same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs are LIRGs. However, none of them are in the ULIRG regime, where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX-$β$ diagram, PACS-detected LBGs at z ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at their same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in age-dependent way: the overestimation factor is higher in younger galaxies.

preprint2013arXiv

The deepest Herschel-PACS far-infrared survey: number counts and infrared luminosity functions from combined PEP/GOODS-H observations

We present results from the deepest Herschel-PACS (Photodetector Array Camera and Spectrometer) far-infrared blank field extragalactic survey, obtained by combining observations of the GOODS (Great Observatories Origins Deep Survey) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and the construction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3-sigma depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 um, respectively, and resolve ~75% of the cosmic infrared background at 100um and 160um into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths and to determine the infrared luminosity function of galaxies down to LIR=10^11 Lsun at z~1 and LIR=10^12 Lsun at z~2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3-sigma) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts.

preprint2013arXiv

The evolution of the dust and gas content in galaxies

We use deep Herschel PACS and SPIRE observations in GOODSS, GOODSN and COSMOS to estimate the average dust mass (Mdust) of galaxies on a redshift-stellar mass (Mstar)-SFR grid. We study the scaling relations between Mdust, Mstar and SFR at z<=2.5. No clear evolution of Mdust is observed at fixed SFR and Mstar. We find a tight correlation between SFR and Mdust, likely a consequence of the Schmidt-Kennicutt (S-K) law. The Mstar-Mdust correlation observed by previous works flattens or sometimes disappears when fixing the SFR. Most of it likely derives from the combination of the Mdust-SFR and Mstar-SFR correlations. We then investigate the gas content as inferred by converting Mdust by assuming that the dust/gas ratio scales linearly with the gas metallicity. All galaxies in the sample follow, within uncertainties, the same SFR-Mgas relation (integrated S-K law), which broadly agrees with CO-based results for the bulk of the population, despite the completely different approaches. The majority of galaxies at z~2 form stars with an efficiency (SFE=SFR/Mgas) ~5 times higher than at z~0. It is not clear what fraction of such variation is an intrinsic redshift evolution and what fraction arises from selection effects. The gas fraction (fgas) decreases with Mstar and increases with SFR, and does not evolve with z at fixed Mstar and SFR. We explain these trends by introducing a universal relation between fgas, Mstar and SFR, non-evolving out to z~2.5. Galaxies move across this relation as their gas content evolves in time. We use the 3D fundamental fgas-Mstar-SFR relation and the redshift evolution of the Main Sequence to estimate the evolution of fgas in the average population of galaxies as a function of z and Mstar, and we find evidence a downsizing scenario.

preprint2013arXiv

The evolution of the dust temperatures of galaxies in the SFR$-M_{\ast}$ plane up to $z$$\,\thicksim\,$$2$

[Abridged] We study the evolution of the dust temperatures of galaxies in the SFR-M* plane up to z~2 using observations from the Herschel Space Observatory. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M*) and redshift estimates, we grid the SFR-M* parameter space in several redshift ranges and estimate the mean Tdust of each SFR-M*-z bin. Dust temperatures are inferred using the stacked far-infrared flux densities of our SFR-M*-z bins. At all redshifts, Tdust increases with infrared luminosities (LIR), specific SFRs (SSFR; i.e., SFR/M*) and distances with respect to the main sequence (MS) of the SFR-M* plane (i.e., D_SSFR_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). The Tdust-SSFR and Tdust-D_SSFR_MS correlations are statistically more significant than the Tdust-LIR one. While the slopes of these three correlations are redshift-independent, their normalizations evolve from z=0 and z~2. We convert these results into a recipe to derive Tdust from SFR, M* and z. The existence of a strong Tdust-D_SSFR_MS correlation provides us with information on the dust and gas content of galaxies. (i) The slope of the Tdust-D__SSFR_MS correlation can be explained by the increase of the star-formation efficiency (SFE; SFR/Mgas) with D_SSFR_MS as found locally by molecular gas studies. (ii) At fixed D_SSFR_MS, the constant Tdust observed in galaxies probing large ranges in SFR and M* can be explained by an increase or decrease of the number of star-forming regions with comparable SFE enclosed in them. (iii) At high redshift, the normalization towards hotter temperature of the Tdust-D_SSFR_MS correlation can be explained by the decrease of the metallicities of galaxies or by the increase of the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different.

preprint2013arXiv

The Herschel PEP/HerMES Luminosity Function. I: Probing the Evolution of PACS selected Galaxies to z~4

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

preprint2013arXiv

The Herschel-PEP survey: evidence for downsizing in the hosts of dusty star-forming systems

By making use of Herschel-PEP observations of the COSMOS and Extended Groth Strip fields, we have estimated the dependence of the clustering properties of FIR-selected sources on their 100um fluxes. Our analysis shows a tendency for the clustering strength to decrease with limiting fluxes: r0(S100um >8 mJy)~4.3 Mpc and r0(S100um >5 mJy)~5.8 Mpc. These values convert into minimum halo masses Mmin~10^{11.6} Msun for sources brighter than 8 mJy and Mmin~10^{12.4} Msun for S100um > 5 mJy galaxies. We show such an increase of the clustering strength to be due to an intervening population of z~2 sources, which are very strongly clustered and whose relative contribution, equal to about 10% of the total counts at S100um > 2 mJy, rapidly decreases for brighter flux cuts. By removing such a contribution, we find that z <~ 1 FIR galaxies have approximately the same clustering properties, irrespective of their flux level. The above results were then used to investigate the intrinsic dependence on cosmic epoch of the clustering strength of dusty star-forming galaxies between z~0 and z~2.5. In order to remove any bias in the selection process, the adopted sample only includes galaxies observed at the same rest-frame wavelength, lambda~60 um, which have comparable luminosities and therefore star-formation rates (SFR>~100 Msun/yr). Our analysis shows that the same amount of (intense) star forming activity takes place in extremely different environments at the different cosmological epochs. For z<~1 the hosts of such star forming systems are small, Mmin~10^{11} Msun, isolated galaxies. High (z~2) redshift star formation instead seems to uniquely take place in extremely massive/cluster-like halos, Mmin~10^{13.5} Msun, which are associated with the highest peaks of the density fluctuation field at those epochs. (abridged)

preprint2013arXiv

The Mean Star Formation Rate of X-ray selected Active Galaxies and its Evolution from z=2.5: Results from PEP-Herschel

(Abridged) We study relationships between the SFR and the nuclear properties of X-ray selected AGNs out to z=2.5, using far-IR data in three extragalactic deep fields as part of the PACS Evolutionary Probe (PEP) program. Guided by studies of intrinsic infra-red AGN SEDs, we show that the majority of the FIR emission in AGNs is produced by cold dust heated by star-formation. We uncover characteristic redshift-dependent trends between the mean FIR luminosity (L_fir) and accretion luminosity (L_agn) of AGNs. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L_fir and L_agn, but only among AGNs at low and moderate redshifts (z<1). We interpret this as a signature of the increasing importance of major-mergers in driving both the growth of super-massive black holes (SMBHs) and global star-formation in their hosts at high AGN luminosities. However, we also find that the enhancement of SFR in luminous AGNs weakens or disappears at high redshifts (z>1). This suggests that the role of mergers in SMBH-galaxy co-evolution is less important at these epochs. At all redshifts, we find essentially no relationship between L_fir and nuclear obscuration across five orders of magnitude in obscuring column density, suggesting that various different mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We explain our results within a scenario in which two different modes of SMBH fueling operate among low- and high-luminosity AGNs. We postulate, guided by emerging knowledge about the properties of high redshift galaxies, that the dominant mode of accretion among high-luminosity AGNs evolves with redshift.

preprint2012arXiv

A Herschel view of the far-infrared properties of submillimetre galaxies

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and thus might be evolving through a secular mode of star formation. [abridged]

preprint2012arXiv

Enhanced star formation rates in AGN hosts with respect to inactive galaxies from PEP-Herschel observations

We compare the average star formation (SF) activity in X-ray selected AGN hosts with mass-matched control inactive galaxies,including star forming and quiescent sources, at 0.5<z<2.5. Recent observations carried out by PACS, the 60-210um Herschel photometric camera, in GOODS-S, GOODS-N and COSMOS allow us to unbiasedly estimate the far-IR luminosity, and hence the SF properties, of the two samples. Accurate AGN host stellar masses are measured by decomposing their total emission into the stellar and nuclear components. We find a higher average SF activity in AGN hosts with respect to non-AGNs. The level of SF enhancement is modest (~0.26dex at ~3sigma) at low X-ray luminosities (Lx<~10^43.5erg/s) and more pronounced (0.56dex at >10sigma) for bright AGNs. However, when comparing to star forming galaxies only, AGN hosts are broadly consistent with the locus of their `main sequence'. We investigate the relative far-IR luminosity distributions of active and inactive galaxies, and find a higher fraction of PACS detected, hence normal and highly star forming systems among AGN hosts. Although different interpretations are possible, we explain our findings as a consequence of a twofold AGN growth path: faint AGNs evolve through secular processes, with instantaneous AGN accretion not tightly linked to the current total SF in the host, while luminous AGNs co-evolve with their hosts through periods of enhanced AGN activity and SF, possibly through major mergers. While an increased SF with respect to non-AGNs of similar mass is expected in the latter, we interpret the modest SF offsets measured in low-Lx AGN hosts as either a) generated by non-synchronous accretion and SF histories in a merger scenario or b) due to possible connections between instantaneous SF and accretion that can be induced by smaller scale (non-major merger) mechanisms. Far-IR luminosity distributions favour the latter scenario.

preprint2012arXiv

PACS-Herschel FIR detections of Lyman-alpha emitters at 2.0<z<3.5

In this work we analyze the physical properties of a sample of 56 spectroscopically selected star-forming (SF) Ly$α$ emitting galaxies at 2.0$\lesssim$z$\lesssim$3.5 using both a spectral energy distribution (SED) fitting procedure from rest-frame UV to mid-IR and direct 160$μ$m observations taken with the Photodetector Array Camera & Spectrometer (PACS) instrument onboard \emph{Herschel Space Observatory}. We define LAEs as those Ly$α$ emitting galaxies whose rest-frame Ly$α$ equivalent widths (Ly$α$ EW$_{rest-frame}$) are above 20Å, the typical threshold in narrow-band searches. Ly$α$ emitting galaxies with Ly$α$ EW$_{rest-frame}$ are called non-LAEs. As a result of an individual SED fitting for each object, we find that the studied sample of LAEs contains galaxies with ages mostly below 100Myr and a wide variety of dust attenuations, SFRs, and stellar masses. The heterogeneity in the physical properties is also seen in the morphology, ranging from bulge-like galaxies to highly clumpy systems. In this way, we find that LAEs at 2.0$\lesssim$z$\lesssim$3.5 are very diverse, and do not have a bimodal nature, as suggested in previous works. Furthermore, the main difference between LAEs and non-LAEs is their dust attenuation, because LAEs are not as dusty as non-LAEs. On the FIR side, four galaxies of the sample (two LAEs and two non-LAEs) have PACS-FIR counterparts. Their total IR luminosity place all of them in the ULIRG regime and are all dusty objects, with A$_{1200}$$\gtrsim$4mag. This is an indication from direct FIR measurements that dust and Ly$α$ emission are not mutually exclusive. This population of red and dusty LAEs is not seen at z$\sim$0.3, suggesting an evolution with redshift of the IR nature of galaxies selected via their Ly$α$ emission.

preprint2012arXiv

Physical properties of Lyman-alpha emitters at $z\sim 0.3$ from UV-to-FIR measurements

The analysis of the physical properties of low-redshift Ly$α$ emitters (LAEs) can provide clues in the study of their high-redshift analogues. At $z \sim 0.3$, LAEs are bright enough to be detected over almost the entire electromagnetic spectrum and it is possible to carry out a more precise and complete study than at higher redshifts. In this study, we examine the UV and IR emission, dust attenuation, SFR and morphology of a sample of 23 GALEX-discovered star-forming (SF) LAEs at $z \sim 0.3$ with direct UV (GALEX), optical (ACS) and FIR (PACS and MIPS) data. Using the same UV and IR limiting luminosities, we find that LAEs at $z\sim 0.3$ tend to be less dusty, have slightly higher total SFRs, have bluer UV continuum slopes, and are much smaller than other galaxies that do not exhibit Ly$α$ emission in their spectrum (non-LAEs). These results suggest that at $z \sim 0.3$ Ly$α$ photons tend to escape from small galaxies with low dust attenuation. Regarding their morphology, LAEs belong to Irr/merger classes, unlike non-LAEs. Size and morphology represent the most noticeable difference between LAEs and non-LAEs at $z \sim 0.3$. Furthermore, the comparison of our results with those obtained at higher redshifts indicates that either the Ly$α$ technique picks up different kind of galaxies at different redshifts or that the physical properties of LAEs are evolving with redshift.

preprint2012arXiv

The AGN content in luminous IR galaxies at z\sim2 from a global SED analysis including Herschel data

We use Herschel-PACS far-infrared data, combined with previous multi-band information and mid-IR spectra, to properly account for the presence of an active nucleus and constrain its energetic contribution in luminous infrared (IR) sources at z\sim2. The sample is composed of 24 sources in the GOODS-South field, with typical IR luminosity of 10^{12} Lo. Data from the 4 Ms Chandra X-ray imaging in this field are also used to identify and characterize AGN emission. We reproduce the observed spectral energy distribution (SED), decomposed into a host-galaxy and an AGN component. A smooth-torus model for circum-nuclear dust is used to account for the direct and re-processed contribution from the AGN. We confirm that galaxies with typical L_{8-1000um}\sim10^{12}Lo at z\sim2 are powered predominantly by star-formation. An AGN component is present in nine objects (\sim35% of the sample) at the 3sigma confidence level, but its contribution to the 8-1000 um emission accounts for only \sim5% of the energy budget. The AGN contribution rises to \sim23% over the 5-30 um range (in agreement with Spitzer IRS results) and to \sim60% over the narrow 2-6 um range. The presence of an AGN is confirmed by X-ray data for 3 (out of nine) sources, with X-ray spectral analysis indicating the presence of significant absorption, i.e. NH\sim10^{23}-10^{24} cm^{-2}. An additional source shows indications of obscured AGN emission from X-ray data. The comparison between the mid-IR--derived X-ray luminosities and those obtained from X-ray data suggests that obscuration is likely present also in the remaining six sources that harbour an AGN according to the SED-fitting analysis.

preprint2012arXiv

Two modified ILC methods to detect point sources in Cosmic Microwave Background maps

We propose two detection techniques that take advantage of a small sky area approximation and are based on modifications of the "internal linear combination" (ILC) method, an approach widely used in Cosmology for the separation of the various components that contribute to the microwave background. The main advantage of the proposed approach, especially in handling multi-frequency maps of the same region, is that it does not require the "a priori" knowledge of the spatial power-spectrum of either the CMB and/or the Galactic foreground. Hence, it is more robust, easier and more intuitive to use. The performance of the proposed algorithms is tested with numerical experiments that mimic the physical scenario expected for high Galactic latitude observations with the Atacama Large Millimeter/submillimeter Array (ALMA).

preprint2011arXiv

Building the cosmic infrared background brick by brick with Herschel/PEP

The cosmic infrared background (CIB) includes roughly half of the energy radiated by all galaxies at all wavelengths across cosmic time, as observed at the present epoch. The PACS Evolutionary Probe (PEP) survey is exploited here to study the CIB and its redshift differential, at 70, 100 and 160 micron, where the background peaks. Combining PACS observations of the GOODS-S, GOODS-N, Lockman Hole and COSMOS areas, we define number counts spanning over more than two orders of magnitude in flux: from ~1 mJy to few hundreds mJy. Stacking of 24 micron sources and P(D) statistics extend the analysis down to ~0.2 mJy. Taking advantage of the wealth of ancillary data in PEP fields, differential number counts and CIB are studied up to z=5. Based on these counts, we discuss the effects of confusion on PACS blank field observations and provide confusion limits for the three bands considered. The total CIB surface brightness emitted above PEP 3 sigma flux limits is 4.52 +/- 1.18, 8.35 +/- 0.95 and 9.49 +/- 0.59 [nW/m2/sr] at 70, 100, and 160 micron, respectively. These values correspond to 58 +/- 7% and 74 +/- 5% of the COBE/DIRBE CIB direct measurements at 100 and 160 micron. Employing the P(D) analysis, these fractions increase to ~65% and ~89%. More than half of the resolved CIB was emitted at redshift z<=1. The 50%-light redshifts lie at z=0.58, 0.67 and 0.73 at the three PACS wavelengths. The distribution moves towards earlier epochs at longer wavelengths: while the 70 micron CIB is mainly produced by z<=1.0 objects, the contribution of z>1.0 sources reaches 50% at 160 micron. Most of the CIB resolved in the three PACS bands was emitted by galaxies with infrared luminosities in the range 1e11-1e12 L(sun).

preprint2011arXiv

FIR measurements of Ly-$α$ emitters at z$\lesssim$1.0: dust attenuation from PACS-\emph{Herschel}

One remaining open question regarding the physical properties of Ly$α$ emitters (LAEs) is their dust content and its evolution with redshift. The variety of results is large and with those reported by now is difficult to establish clear relations between dust, other fundamental parameters of galaxies (star-formation rate, metallicity or age) and redshift. In this Letter, we report \emph{Herschel} PACS-100$μ$m, PACS-160$μ$m and \emph{Spitzer} MIPS-24$μ$m detections of a sample of spectroscopically GALEX-selected LAEs at z$\sim$0.3 and $\sim$1.0. Five out of ten and one out of two LAEs are detected in, at least, one PACS band at z$\sim$0.3 and $\sim$1.0, respectively. These measurements have a great importance given that they allow us to quantify, for the first time, the dust content of LAEs from direct FIR observations. MIPS-24$μ$m detections allow us to determine IR properties of the PACS-undetected LAEs. We obtain that mid-IR/FIR detected star-forming (SF) LAEs at z$\sim$0.3 have dust content within 0.75$\lesssim$ $A_{1200Å}$ $\lesssim$2.0, with a median value of A$_{1200\textrmÅ}$$\sim$1.1. This range broadens out to 0.75$\lesssim$ $A_{1200Å}$ $\lesssim$2.5 when considering those LAEs at z$\sim$1.0. Only one SF LAE is undetected both in MIPS-24$μ$m and PACS, with $A_{1200Å}$ $\lesssim$0.75. These results seem to be larger than those reported for high-redshift LAEs and, therefore, although an evolutionary trend is not clearly seen, it could point out that low-redshift LAEs are dustier than high-redshift ones. However, the diverse methods used could introduce a systematic offset in the results.

preprint2011arXiv

PACS Evolutionary Probe (PEP) - A Herschel Key Program

Deep far-infrared photometric surveys studying galaxy evolution and the nature of the cosmic infrared background are a key strength of the Herschel mission. We describe the scientific motivation for the PACS Evolutionary Probe (PEP) guaranteed time key program and its role in the complement of Herschel surveys, and the field selection which includes popular multiwavelength fields such as GOODS, COSMOS, Lockman Hole, ECDFS, EGS. We provide an account of the observing strategies and data reduction methods used. An overview of first science results illustrates the potential of PEP in providing calorimetric star formation rates for high redshift galaxy populations, thus testing and superseeding previous extrapolations from other wavelengths, and enabling a wide range of galaxy evolution studies.

preprint2011arXiv

The evolution of the star formation activity per halo mass up to redshift ~ 1.6 as seen by Herschel

Star formation in massive galaxies is quenched at some point during hierarchical mass assembly. To understand where and when the quenching processes takes place, we study the evolution of the total star formation rate per unit total halo mass (Σ(SFR/M)) in three different mass scales: low mass halos (field galaxies), groups, and clusters, up to a redshift ~1.6. We use deep far-infrared PACS data at 100 and 160 um to accurately estimate the total star formation rate of the Luminous Infrared Galaxy population of 9 clusters with mass ~10^{15} M_{\odot}, and 9 groups/poor clusters with mass ~ 5 x 10^{13} M_{\odot}. Estimates of the field Σ(SFR/M) are derived from the literature, by dividing the star formation rate density by the mean comoving matter density of the universe. The field Σ(SFR/M) increases with redshift up to z~1 and it is constant thereafter. The evolution of the Σ(SFR/M)-z relation in galaxy systems is much faster than in the field. Up to redshift z~0.2, the field has a higher Σ(SFR/M) than galaxy groups and galaxy clusters. At higher redshifts, galaxy groups and the field have similar Σ(SFR/M), while massive clusters have significantly lower Σ(SFR/M) than both groups and the field. There is a hint of a reversal of the SFR activity vs. environment at z~1.6, where the group Σ(SFR/M) lies above the field Σ(SFR/M)-z relation. We discuss possible interpretations of our results in terms of the processes of downsizing, and star-formation quenching.

preprint2011arXiv

The impact of evolving infrared spectral energy distributions of galaxies on star formation rate estimates

We combine Herschel-PACS data from the PEP program with Spitzer 24 um and 16 um photometry and ultra deep IRS mid-infrared spectra, to measure the mid- to far-infrared spectral energy distribution (SED) of 0.7<z<2.5 normal star forming galaxies around the main sequence (the redshift-dependent relation of star formation rate and stellar mass). Our deep data confirm from individual far-infrared detections that z~2 star formation rates are overestimated if based on 24 um fluxes and SED templates that are calibrated via local trends with luminosity. Galaxies with similar ratios of rest-frame nuLnu(8) to 8-1000 um infrared luminosity (LIR) tend to lie along lines of constant offset from the main sequence. We explore the relation between SED shape and offset in specific star formation rate (SSFR) from the redshift-dependent main sequence. Main sequence galaxies tend to have a similar nuLnu(8)/LIR regardless of LIR and redshift, up to z~2.5, and nuLnu(8)/LIR decreases with increasing offset above the main sequence in a consistent way at the studied redshifts. We provide a redshift-independent calibration of SED templates in the range of 8--60 um as a function of log(SSFR) offset from the main sequence. Redshift dependency enters only through the evolution of the main sequence with time. Ultra deep IRS spectra match these SED trends well and verify that they are mostly due to a change in ratio of PAH to LIR rather than continua of hidden AGN. Alternatively, we discuss the dependence of nuLnu(8)/LIR on LIR. Same nuLnu(8)/LIR is reached at increasingly higher LIR at higher redshift, with shifts relative to local by 0.5 and 0.8 dex in log(LIR) at redshifts z~1 and z~2. Corresponding SED template calibrations are provided for use if no stellar masses are in hand. For most of those z~2 star forming galaxies that also host an AGN, the mid-infrared is dominated by the star forming component.

preprint2011arXiv

The lesser role of starbursts for star formation at z=2

Two main modes of star formation are know to control the growth of galaxies: a relatively steady one in disk-like galaxies, defining a tight star formation rate (SFR)-stellar mass sequence, and a starburst mode in outliers to such a sequence which is generally interpreted as driven by merging. Such starburst galaxies are rare but have much higher SFRs, and it is of interest to establish the relative importance of these two modes. PACS/Herschel observations over the whole COSMOS and GOODS-South fields, in conjunction with previous optical/near-IR data, have allowed us to accurately quantify for the first time the relative contribution of the two modes to the global SFR density in the redshift interval 1.5<z<2.5, i.e., at the cosmic peak of the star formation activity. The logarithmic distributions of galaxy SFRs at fixed stellar mass are well described by Gaussians, with starburst galaxies representing only a relatively minor deviation that becomes apparent for SFRs more than 4 times higher than on the main sequence. Such starburst galaxies represent only 2% of mass-selected star forming galaxies and account for only 10% of the cosmic SFR density at z~2. Only when limited to SFR>1000M(sun)/yr, off-sequence sources significantly contribute to the SFR density (46+/-20%). We conclude that merger-driven starbursts play a relatively minor role for the formation of stars in galaxies, whereas they may represent a critical phase towards the quenching of star formation and morphological transformation in galaxies.

preprint2011arXiv

The PEP survey: clustering of infrared-selected galaxies and structure formation at z~2 in the GOODS South

ABRIDGED-This paper presents the first direct estimate of the 3D clustering properties of far-infrared sources up to z~3. This has been possible thanks to the Pacs Evolutionary Probe (PEP) survey of the GOODS South field performed with the PACS instrument onboard the Herschel Satellite. An analysis of the two-point correlation function over the whole redshift range spanned by the data reports for the correlation length, r_0~6.3 Mpc and r_0~6.7 Mpc, respectively at 100um and 160um, corresponding to dark matter halo masses M>~10^{12.4} M_sun. Objects at z~2 instead seem to be more strongly clustered, with r_0~19 Mpc and r_0~17 Mpc in the two considered PACS channels. This dramatic increase of the correlation length between z~1 and z~2 is connected with the presence of a wide, M>~10^{14} M_sun, filamentary structure which includes more than 50% of the sources detected at z~2. An investigation of the properties of such sources indicates the possibility for boosted star-forming activity in those which reside within the overdense environment with respect of more isolated galaxies found in the same redshift range. Lastly, we also present our results on the evolution of the relationship between luminous and dark matter in star-forming galaxies between z~1 and z~2. We find that the increase of (average) stellar mass in galaxies <M*> between z~1 and z~2 is about a factor 10 lower than that of the dark matter haloes hosting such objects (<M*>[z~1]/<M*>[z~2] ~ 0.4 vs M_{halo}[z~1]/M_{halo}[z~2] ~ 0.04). Our findings agree with the evolutionary picture of downsizing whereby massive galaxies at z~2 were more actively forming stars than their z~1 counterparts, while at the same time contained a lower fraction of their mass in the form of luminous matter.

preprint2010arXiv

Dissecting the cosmic infra-red background with Herschel/PEP

The constituents of the cosmic IR background (CIB) are studied at its peak wavelengths (100 and 160 um) by exploiting Herschel/PACS observations of the GOODS-N, Lockman Hole, and COSMOS fields in the PACS Evolutionary Probe (PEP) guaranteed-time survey. The GOODS-N data reach 3 sigma depths of ~3.0 mJy at 100 um and ~5.7 mJy at 160 um. At these levels, source densities are 40 and 18 beams/source, respectively, thus hitting the confusion limit at 160 um. Differential number counts extend from a few mJy up to 100-200 mJy, and are approximated as a double power law, with the break lying between 5 and 10 mJy. The available ancillary information allows us to split number counts into redshift bins. At z<=0.5 we isolate a class of luminous sources (L(IR)~1e11 Lsun), whose SEDs resemble late-spiral galaxies, peaking at ~130 um restframe and significantly colder than what is expected on the basis of pre-Herschel models. By integrating number counts over the whole covered flux range, we obtain a surface brightness of 6.36 +/- 1.67 and 6.58 +/-1.62 [nW m^-2 sr^-1] at 100 and 160 um, resolving ~45% and ~52% of the CIB, respectively. When stacking 24 um sources, the inferred CIB lies within 1.1 sigma and 0.5 sigma from direct measurements in the two bands, and fractions increase to 50% and 75%.Most of this resolved CIB fraction was radiated at z<=1.0, with 160 um sources found at higher redshift than 100 um ones.

preprint2010arXiv

Evolution of Dust Temperature of Galaxies through Cosmic Time as seen by Herschel

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

preprint2010arXiv

Far-infrared properties of submillimeter and optically faint radio galaxies

We use deep observations obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel space observatory to study the far-infrared (FIR) properties of submm and optically faint radio galaxies (SMGs and OFRGs). From literature we compiled a sample of 35 securely identified SMGs and nine OFRGs located in the GOODS-N and the A2218 fields. This sample is cross-matched with our PACS 100 um and 160 um multi-wavelength catalogs. About half of the galaxies in our sample are detected with PACS. The dust temperatures and the infrared luminosities of our galaxies are derived by fitting their PACS and SCUBA 850 um (only the upper limits for the OFRGs) flux densities with a single modified (beta=1.5) black body function. The median dust temperature of our SMG sample is T=36+/-8K while for our OFRG sample it is T=47+/-3K. For both samples, median dust temperatures derived from Herschel data agree well with previous estimates. In particular, Chapman et al. (2005) found a dust temperature of T=36+/-7K for a large sample of SMGs assuming the validity of the FIR/radio correlation. The agreement between our studies confirms that the local FIR/radio correlation effectively holds at high redshift even though we find <q>=2.17+/-0.19, a slightly lower value than that observed in local systems. The median IR luminosities of SMGs and OFRGs are 4.6*10^12 Lsun and 2.6*10^12 Lsun, respectively. We note that for both samples the IR luminosity estimates from the radio part of the spectral energy distribution are accurate, while estimates from the mid-IR are considerably (x3) more uncertain. Our observations confirm the remarkably high luminosities of SMGs and thus imply median star-formation rates of 960Msun yr^-1 for SMGs with S(850um)>5mJy and 460Msun yr^-1 for SMGs with S(850um)>2mJy, assuming a Chabrier IMF and no dominant AGN contribution to the far-infrared luminosity.

preprint2010arXiv

Herschel deep far-infrared counts through the Abell 2218 cluster-lens

Gravitational lensing by massive galaxy clusters allows study of the population of intrinsically faint infrared galaxies that lie below the sensitivity and confusion limits of current infrared and submillimeter telescopes. We present ultra-deep PACS 100 and 160 microns observations toward the cluster lens Abell 2218, to penetrate the Herschel confusion limit. We derive source counts down to a flux density of 1 mJy at 100 microns and 2 mJy at 160 microns, aided by strong gravitational lensing. At these levels, source densities are 20 and 10 beams/source in the two bands, approaching source density confusion at 160 microns. The slope of the counts below the turnover of the Euclidean-normalized differential curve is constrained in both bands and is consistent with most of the recent backwards evolutionary models. By integrating number counts over the flux range accessed by Abell 2218 lensing (0.94-35 mJy at 100 microns and 1.47-35 mJy at 160 microns, we retrieve a cosmic infrared background (CIB) surface brightness of ~8.0 and ~9.9 nW m^-2 sr^-1, in the respective bands. These values correspond to 55% (+/- 24%) and 77% (+/- 31%) of DIRBE direct measurements. Combining Abell 2218 results with wider/shallower fields, these figures increase to 62% (+/- 25%) and 88% (+/- 32%) CIB total fractions, resolved at 100 and 160 microns, disregarding the high uncertainties of DIRBE absolute values.

preprint2010arXiv

Herschel FIR counterparts of selected Ly-alpha emitters at z~2.2. Fast evolution since z~3 or missed obscured AGNs?

Ly-alpha emitters (LAEs) are seen everywhere in the redshift domain from local to z~7. Far-infrared (FIR) counterparts of LAEs at different epochs could provide direct clues on dust content, extinction, and spectral energy distribution (SED) for these galaxies. We search for FIR counterparts of LAEs that are optically detected in the GOODS-North field at redshift z~2.2 using data from the Herschel Space Telescope with the Photodetector Array Camera and Spectrometer (PACS). The LAE candidates were isolated via color-magnitude diagram using the medium-band photometry from the ALHAMBRA Survey, ancillary data on GOODS-North, and stellar population models. According to the fitting of these spectral synthesis models and FIR/optical diagnostics, most of them seem to be obscured galaxies whose spectra are AGN-dominated. From the analysis of the optical data, we have observed a fraction of AGN or composite over source total number of ~0.75 in the LAE population at z~2.2, which is marginally consistent with the fraction previously observed at z=2.25 and even at low redshift (0.2<z<0.45), but significantly different from the one observed at redshift ~3, which could be compatible either with a scenario of rapid change in the AGN fraction between the epochs involved or with a non detection of obscured AGN in other z=2-3 LAE samples due to lack of deep FIR observations. We found three robust FIR (PACS) counterparts at z~2.2 in GOODS-North. This demonstrates the possibility of finding dust emission in LAEs even at higher redshifts.

preprint2010arXiv

Herschel PEP: The star-formation rates of 1.5<z<2.5 massive galaxies

The star formation rate (SFR) is a key parameter in the study of galaxy evolution. The accuracy of SFR measurements at z~2 has been questioned following a disagreement between observations and theoretical models. The latter predict SFRs at this redshift that are typically a factor 4 or more lower than the measurements. We present star-formation rates based on calorimetric measurements of the far-infrared (FIR) luminosities for massive 1.5<z<2.5, normal star-forming galaxies (SFGs), which do not depend on extinction corrections and/or extrapolations of spectral energy distributions. The measurements are based on observations in GOODS-N with the Photodetector Array Camera & Spectrometer (PACS) onboard Herschel, as part of the PACS Evolutionary Probe (PEP) project, that resolve for the first time individual SFGs at these redshifts at FIR wavelengths. We compare FIR-based SFRs to the more commonly used 24 micron and UV SFRs. We find that SFRs from 24 micron alone are higher by a factor of ~4-7.5 than the true SFRs. This overestimation depends on luminosity: gradually increasing for log L(24um)>12.2 L_sun. The SFGs and AGNs tend to exhibit the same 24 micron excess. The UV SFRs are in closer agreement with the FIR-based SFRs. Using a Calzetti UV extinction correction results in a mean excess of up to 0.3 dex and a scatter of 0.35 dex from the FIR SFRs. The previous UV SFRs are thus confirmed and the mean excess, while narrowing the gap, is insufficient to explain the discrepancy between the observed SFRs and simulation predictions.

preprint2010arXiv

Herschel unveils a puzzling uniformity of distant dusty galaxies

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

preprint2010arXiv

Herschel-ATLAS: blazars in the SDP field

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$μ$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $γ$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

preprint2010arXiv

PEP: first Herschel probe of dusty galaxy evolution up to z~3

We exploit the deepest existing far-infrared (FIR) data obtained so far by Herschel at 100 and 160 um in the GOODS-N, as part of the PACS Evolutionary Probe (PEP) survey, to derive for the first time the evolution of the rest-frame 60-um, 90-um, and total IR luminosity functions (LFs) of galaxies and AGNs from z=0 to unprecedented high redshifts (z~2-3). The PEP LFs were computed using the 1/Vmax method. The FIR sources were classified by means of a detailed broad- band SED-fitting analysis and spectral characterisation. Based on the best-fit model results, k-correction and total IR (8-1000 um) luminosity were obtained for each source. LFs (monochromatic and total) were then derived for various IR populations separately in different redshift bins and compared to backward evolution model predictions. We detect strong evolution in the LF to at least z~2. Objects with SEDs similar to local spiral galaxies are the major contributors to the star formation density (SFD) at z< 0.3, then, as redshift increases, moderate SF galaxies - most likely containing a low-luminosity AGN - start dominating up to z ~= 1.5. At >1.5 the SFD is dominated by the contributions of starburst galaxies. In agreement with previous findings, the comoving IR LD derived from our data evolves approximately as (1 + z)^(3.8+/-0.3) up to z~1, there being some evidence of flattening up to z~2.

preprint2010arXiv

Star formation in AGN hosts in GOODS-N

Sensitive Herschel far-infrared observations can break degeneracies that were inherent to previous studies of star formation in high-z AGN hosts. Combining PACS 100 and 160um observations of the GOODS-N field with 2Msec Chandra data, we detect ~20% of X-ray AGN individually at >3sig. The host far-infrared luminosity of AGN with L2-10~10^43erg/s increases with redshift by an order of magnitude from z=0 to z~1. In contrast, there is little dependence of far-infrared luminosity on AGN luminosity, for L2-10<~10^44erg/s AGN at z>~1. We do not find a dependence of far-infrared luminosity on X-ray obscuring column, for our sample which is dominated by L2-10<10^44erg/s AGN. In conjunction with properties of local and luminous high-z AGN, we interpret these results as reflecting the interplay between two paths of AGN/host coevolution. A correlation of AGN luminosity and host star formation is traced locally over a wide range of luminosities and also extends to luminous high z AGN. This correlation reflects an evolutionary connection, likely via merging. For lower AGN luminosities, star formation is similar to that in non-active massive galaxies and shows little dependence on AGN luminosity. The level of this secular, non-merger driven star formation increasingly dominates over the correlation at increasing redshift.

preprint2010arXiv

The dust content of high-z submillimeter galaxies revealed by Herschel

We use deep observations taken with the Photodetector Array Camera and Spectrometer (PACS), on board the Herschel satellite as part of the PACS evolutionary probe (PEP) guaranteed project along with submm ground-based observations to measure the dust mass of a sample of high-z submillimeter galaxies (SMGs). We investigate their dust content relative to their stellar and gas masses, and compare them with local star-forming galaxies. High-z SMGs are dust rich, i.e. they have higher dust-to-stellar mass ratios compared to local spiral galaxies (by a factor of 30) and also compared to local ultraluminous infrared galaxies (ULIRGs, by a factor of 6). This indicates that the large masses of gas typically hosted in SMGs have already been highly enriched with metals and dust. Indeed, for those SMGs whose gas mass is measured, we infer dust-to-gas ratios similar or higher than local spirals and ULIRGs. However, similarly to other strongly star-forming galaxies in the local Universe and at high-z, SMGs are characterized by gas metalicities lower (by a factor of a few) than local spirals, as inferred from their optical nebular lines, which are generally ascribed to infall of metal-poor gas. This is in contrast with the large dust content inferred from the far-IR and submm data. In short, the metalicity inferred from the dust mass is much higher (by more than an order of magnitude) than that inferred from the optical nebular lines. We discuss the possible explanations of this discrepancy and the possible implications for the investigation of the metalicity evolution at high-z.

preprint2010arXiv

The far-infrared/radio correlation as probed by Herschel

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

preprint2010arXiv

The first Herschel view of the mass-SFR link in high-z galaxies

We exploit deep observations of the GOODS-N field taken with PACS, on board of Herschel, as part of the PEP guaranteed time, to study the link between star formation and stellar mass in galaxies to z~2. Starting from a stellar mass-selected sample of ~4500 galaxies with mag[4.5mu]<23 (AB), we identify ~350 objects with a PACS detection at 100 or 160mu and ~1500 with only Spitzer 24 mu counterpart. Stellar masses and total IR luminosities (LIR) are estimated by fitting the SEDs. Consistently with other Herschel results, we find that LIR based only on 24 mu data is overestimated by a median factor ~1.8 at z~2, whereas it is underestimated (with our approach) up to a factor ~1.6 at 0.5<z<1.0. We then exploit this calibration to correct LIR based on the MIPS fluxes. These results clearly show how Herschel is fundamental to constrain LIR, and hence the SFR, of high redshift galaxies. Using the galaxies detected with PACS (and/or MIPS), we investigate the existence and evolution of the relations between the star formation rate (SFR), the specific star formation rate (SSFR=SFR/mass) and the stellar mass. Moreover, in order to avoid selection effects, we also repeat this study through a stacking analysis on the PACS images to fully exploit the far-IR information also for the Herschel and Spitzer undetected subsamples. We find that the SSFR-mass relation steepens with redshift, being almost flat at z<1.0 and reaching a slope of alpha=-0.50^(+0.13)_(-0.16) at z~2, at odds with recent works based on radio-stacking analysis at the same redshift. The mean SSFR of galaxies increases with redshift, by a factor ~15 for massive M>10^(11) (M_sun) galaxies from z=0 to z=2, and seems to flatten at z>1.5 in this mass range. Moreover, the most massive galaxies have the lowest SSFR at any z, implying that they have formed their stars earlier and more rapidly than their low mass counterparts.

preprint2010arXiv

The Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory

The Photodetector Array Camera and Spectrometer (PACS) is one of the three science instruments on ESA's far infrared and submillimetre observatory. It employs two Ge:Ga photoconductor arrays (stressed and unstressed) with 16x25 pixels, each, and two filled silicon bolometer arrays with 16x32 and 32x64 pixels, respectively, to perform integral-field spectroscopy and imaging photometry in the 60-210μ m wavelength regime. In photometry mode, it simultaneously images two bands, 60-85μ m or 85-125μ\m and 125-210μ m, over a field of view of ~1.75'x3.5', with close to Nyquist beam sampling in each band. In spectroscopy mode, it images a field of 47"x47", resolved into 5x5 pixels, with an instantaneous spectral coverage of ~1500km/s and a spectral resolution of ~175km/s. We summarise the design of the instrument, describe observing modes, calibration, and data analysis methods, and present our current assessment of the in-orbit performance of the instrument based on the Performance Verification tests. PACS is fully operational, and the achieved performance is close to or better than the pre-launch predictions.

preprint2010arXiv

Unevenly-sampled signals: a general formalism of the Lomb-Scargle periodogram

The periodogram is a popular tool that tests whether a signal consists only of noise or if it also includes other components. The main issue of this method is to define a critical detection threshold that allows identification of a component other than noise, when a peak in the periodogram exceeds it. In the case of signals sampled on a regular time grid, determination of such a threshold is relatively simple. When the sampling is uneven, however, things are more complicated. The most popular solution in this case is to use the "Lomb-Scargle" periodogram, but this method can be used only when the noise is the realization of a zero-mean, white (i.e. flat-spectrum) random process. In this paper, we present a general formalism based on matrix algebra, which permits analysis of the statistical properties of a periodogram independently of the characteristics of noise (e.g. colored and/or non-stationary), as well as the characteristics of sampling.

preprint2009arXiv

MAMBO Observations at 240GHz of optically obscured Spitzer sources: source clumps and radio activity at high redshift

Optically very faint R>25.5 sources detected by the Spitzer Space Telescope at 24um represent a very interesting population at redshift z [1.5-3]. They exhibit strong clustering properties, implying that they are hosted by very massive halos, and their mid-IR emission could be powered by either dust-enshrouded star-formation and/or by an obscured AGN. We report observations carried out with the MAMBO array at the IRAM 30m antenna on Pico Veleta of a candidate protocluster with five optically obscured sources selected from the 24um Spitzer sample of the First Look Survey. Interestingly, these sources appear to lie on a high density filament aligned with the two radio jets of an AGN. Four out of five of the observed sources were detected. We combine these measurements with optical, infrared and radio observations to probe}the nature of the candidate protocluster members. Our preliminary conclusions can be summarized as follows: the Spectral Energy Distributions of all sources include both AGN and starburst contributions; the AGN contribution to the bolometric luminosities ranges between 14 and 26% of the total. Such a contribution is enough for the AGN to dominate the emission at 5.8, 8 and 24um, while the stellar component, inferred from SED fitting, prevails at 1.25mm and at lambda<4.5um. The present analysis suggests a coherent interplay at high-z between extended radio activity and the development of filamentary large-scale structures.

preprint2009arXiv

The Herschel ATLAS

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

preprint2008arXiv

A Statistical Analysis of the "Internal Linear Combination" Method in Problems of Signal Separation as in CMB Observations

AIMS: The separation of foreground contamination from cosmic microwave background (CMB) observations is one of the most challenging and important problem of digital signal processing in Cosmology. In literature, various techniques have been presented, but no general consensus about their real performances and properties has been reached. This is due to the characteristics of these techniques that have been studied essentially through numerical simulations based on semi-empirical models of the CMB and the Galactic foregrounds. Such models often have different level of sophistication and/or are based on different physical assumptions (e.g., the number of the Galactic components and the level of the noise). Hence, a reliable comparison is difficult. What actually is missing is a statistical analysis of the properties of the proposed methodologies. Here, we consider the "Internal Linear Combination" method (ILC) which, among the separation techniques, requires the smallest number of "a priori" assumptions. This feature is of particular interest in the context of the CMB polarization measurements at small angular scales where the lack of knowledge of the polarized backgrounds represents a serious limit. METHODS: The statistical characteristics of ILC are examined through an analytical approach and the basic conditions are fixed in a way to work satisfactorily. RESULTS: ILC provides satisfactory results only under rather restrictive conditions. This is a critical fact to take into consideration in planning the future ground-based observations (e.g., with ALMA) where, contrary to the satellite experiments, there is the possibility to have a certain control of the experimental conditions.

preprint2006arXiv

Stochastic modeling of kHz QPO light curves

The Kluzniak & Abramowicz model explains high frequency, double peak, "3:2" QPOs observed in neutron star and black hole sources in terms of a non-linear parametric resonance between radial and vertical epicyclic oscillations of an almost Keplerian accretion disk. The 3:2 ratio of epicyclic frequencies occurs only in strong gravity. Rebusco (2004) and Horak (2004) studied the model analytically: they proved that a small forcing may indeed excite the parametric 3:2 resonance, but they have not explained the physical nature of the forcing. Here we integrate their equations numerically, dropping the ad hoc forcing, and adding instead a stochastic term to mimic the action of the very complex processes that occur in disks as, for example, MRI turbulence. We demonstrate that the presence of the stochastic term triggers the resonance in epicyclic oscillations of nearly Keplerian disks, and influences their pattern.

preprint1996arXiv

1.25 mm Observations of a Complete Sample of IRAS Galaxies: (II) Dust Properties

We present 1.25 mm continuum data for a southern galaxy sample selected from the IRAS PSC and complete to S_60=2 Jy. Two thirds of the galaxies have been detected and significant limits on the remaining objects have been set. We find, on a statistical basis, indications that the dust emission in these galaxies is somewhat more centrally concentrated than that of the optical light, possibly tracing a higher metal content in the inner galactic regions. This result also allows to estimate the aperture corrections to the millimetric data. The latter, together with IRAS photometric data, have been used to compare the broad-band FIR/mm spectra with a simple dust model. According to their far-IR/mm spectrum, the sample galaxies show a dichothomy: almost half of the objects, those displaying bright 25-60um fluxes ascribed to warm dust residing in 'starburst' regions, are characterized by higher values of the bolometric (optical + FIR) luminosity, of the dust-to-gas mass ratio, of the dust optical depths and of the overall extinction. A complementary class of objects dominated by cold dust ('cirrus') shows opposite trends. Because of the favourable observational setup, selection wavelength and completeness, we believe these data provide an exhaustive and unbiased view of dust properties in spiral galaxies.

preprint1996arXiv

MM Observations of IRAS Galaxies: Dust Properties, Luminosity Functions and Contributions to the Sub-MM Background

We have studied the FIR/{\it mm} spectrum of IR galaxies by combining IRAS photometry with new {\it mm} data on a complete southern IRAS galaxy sample. The observed spectra and a dust model emphasize a dicothomy in the galaxy population: half of the objects with a lot of warm dust are characterized by higher values of the bolometric (UV-FIR) luminosity, of the dust-to-gas mass ratio, of the dust optical depths and extinction, while those dominated by cold ({\it cirrus}) dust show opposite trends. From these data we derive the {\it mm} luminosity function of galaxies and estimate their contribution to the sub-{\it mm} background (BKG).