Star Formation Rate Tracers Research Articles

An uncertainty principle for star formation – III. The characteristic emission time-scales of star formation rate tracers

ABSTRACT We recently presented a new statistical method to constrain the physics of star formation and feedback on the cloud scale by reconstructing the underlying evolutionary timeline. However, by itself this new method only recovers the relative durations of different evolutionary phases. To enable observational applications, it therefore requires knowledge of an absolute ‘reference time-scale’ to convert relative time-scales into absolute values. The logical choice for this reference time-scale is the duration over which the star formation rate (SFR) tracer is visible because it can be characterized using stellar population synthesis (SPS) models. In this paper, we calibrate this reference time-scale using synthetic emission maps of several SFR tracers, generated by combining the output from a hydrodynamical disc galaxy simulation with the SPS model slug2. We apply our statistical method to obtain self-consistent measurements of each tracer’s reference time-scale. These include H α and 12 ultraviolet (UV) filters (from GALEX, Swift, and HST), which cover a wavelength range 150–350 nm. At solar metallicity, the measured reference time-scales of H α are ${4.32^{+0.09}_{-0.23}}$ Myr with continuum subtraction, and 6–16 Myr without, where the time-scale increases with filter width. For the UV filters we find 17–33 Myr, nearly monotonically increasing with wavelength. The characteristic time-scale decreases towards higher metallicities, as well as to lower star formation rate surface densities, owing to stellar initial mass function sampling effects. We provide fitting functions for the reference time-scale as a function of metallicity, filter width, or wavelength, to enable observational applications of our statistical method across a wide variety of galaxies.

An uncertainty principle for star formation – V. The influence of dust extinction on star formation rate tracer lifetimes and the inferred molecular cloud lifecycle

ABSTRACT Recent observational studies aiming to quantify the molecular cloud lifecycle require the use of known ‘reference time-scales’ to turn the relative durations of different phases of the star formation process into absolute time-scales. We previously constrained the characteristic emission time-scales of different star formation rate (SFR) tracers, as a function of the SFR surface density and metallicity. However, we omitted the effects of dust extinction. Here, we extend our suite of SFR tracer emission time-scales by accounting for extinction, using synthetic emission maps of a high-resolution hydrodynamical simulation of an isolated, Milky Way-like disc galaxy. The stellar feedback included in the simulation is inefficient compared to observations, implying that it represents a limiting case in which the duration of embedded star formation (and the corresponding effect of extinction) is overestimated. Across our experiments, we find that extinction mostly decreases the SFR tracer emission time-scale, changing the time-scales by factors of 0.04–1.74, depending on the gas column density. UV filters are more strongly affected than H α filters. We provide the limiting correction factors as a function of the gas column density and flux sensitivity limit for a wide variety of SFR tracers. Applying these factors to observational characterizations of the molecular cloud lifecycle produces changes that broadly fall within the quoted uncertainties, except at high kpc-scale gas surface densities ($\Sigma _{\rm g}\gtrsim 20~{\mathrm{M_{\odot }\, pc^{-2}}}$). Under those conditions, correcting for extinction may decrease the measured molecular cloud lifetimes and feedback time-scales, which further strengthens previous conclusions that molecular clouds live for a dynamical time and are dispersed by early, pre-supernova feedback.

A panchromatic spatially resolved analysis of nearby galaxies – I. Sub-kpc-scale main sequence in grand-design spirals

ABSTRACT We analyse the spatially resolved relation between stellar mass (M⋆) and star formation rate (SFR) in disc galaxies (i.e. the main sequence, MS). The studied sample includes eight nearby face-on grand-design spirals, e.g. the descendant of high-redshift, rotationally supported star-forming galaxies. We exploit photometric information over 23 bands, from the UV to the far-IR, from the publicly available DustPedia data base to build spatially resolved maps of stellar mass and SFRs on sub-galactic scales of 0.5–1.5 kpc, by performing a spectral energy distribution fitting procedure that accounts for both the observed and obscured star formation processes, over a wide range of internal galaxy environments (bulges, spiral arms, and outskirts). With more than 30 000 physical cells, we have derived a definition of the local spatially resolved MS per unit area for discs, log (ΣSFR) = 0.82log (Σ*) − 8.69. This is consistent with the bulk of recent results based on optical IFU, using the H α line emission as an SFR tracer. Our work extends the analysis at lower sensitivities in both M⋆ and SFR surface densities, up to a factor of ∼10. The self-consistency of the MS relation over different spatial scales, from sub-galactic to galactic, as well as with a rescaled correlation obtained for high-redshift galaxies, clearly proves its universality.

The effect of diffuse background on the spatially-resolved Schmidt relation in nearby spiral galaxies

Context. The global Schmidt law of star formation provides a power-law relation between the surface densities of star-formation rate (SFR) and gas, and successfully explains plausible scenarios of galaxy formation and evolution. However, star formation being a multi-scale process, requires spatially-resolved analysis for a better understanding of the physics of star formation. Aims. It has been shown that the removal of a diffuse background from SFR tracers, such as Hα, far-ultraviolet (FUV), infrared, leads to an increase in the slope of the sub-galactic Schmidt relation. We reinvestigate the local Schmidt relations in nine nearby spiral galaxies taking into account the effect of inclusion and removal of diffuse background in SFR tracers as well as in the atomic gas. Methods. We used multiwavelength data obtained as part of the Spitzer Infrared Nearby Galaxies Survey, Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel, The H I Nearby Galaxy Survey, and HERA CO-Line Extragalactic Survey. Making use of a novel split of the overall light distribution as a function of spatial scale, we subtracted the diffuse background in the SFR tracers as well as the atomic gas. Using aperture photometry, we study the Schmidt relations on background subtracted and unsubtracted data at physical scales varying between 0.5–2 kpc. Results. The fraction of diffuse background varies from galaxy to galaxy and accounts to ∼34% in Hα, ∼43% in FUV, ∼37% in 24 μm, and ∼75% in H I on average. We find that the inclusion of diffuse background in SFR tracers leads to a linear molecular gas Schmidt relation and a bimodal total gas Schmidt relation. However, the removal of diffuse background in SFR tracers leads to a super-linear molecular gas Schmidt relation. A further removal of the diffuse background from atomic gas results in a slope ∼1.4 ± 0.1, which agrees with dynamical models of star formation accounting for flaring effects in the outer regions of galaxies.

A LOFAR-IRAS cross-match study: the far-infrared radio correlation and the 150 MHz luminosity as a star-formation rate tracer

Aims. We aim to study the far-infrared radio correlation (FIRC) at 150 MHz in the local Universe (at a median redshift ⟨z⟩∼0.05) and improve the use of the rest-frame 150 MHz luminosity, L150, as a star-formation rate (SFR) tracer, which is unaffected by dust extinction. Methods. We cross-match the 60 μm selected Revised IRAS Faint Source Survey Redshift (RIFSCz) catalogue and the 150 MHz selected LOFAR value-added source catalogue in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Spring Field. We estimate L150 for the cross-matched sources and compare it with the total infrared (IR) luminosity, LIR, and various SFR tracers. Results. We find a tight linear correlation between log L150 and log LIR for star-forming galaxies, with a slope of 1.37. The median qIR value (defined as the logarithm of the LIR to L150 ratio) and its rms scatter of our main sample are 2.14 and 0.34, respectively. We also find that log L150 correlates tightly with the logarithm of SFR derived from three different tracers, i.e., SFRHα based on the Hα line luminosity, SFR60 based on the rest-frame 60 μm luminosity and SFRIR based on LIR, with a scatter of 0.3 dex. Our best-fit relations between L150 and these SFR tracers are, log L150 (L⊙) = 1.35(±0.06) × log SFRHα (M⊙ yr−1) + 3.20(±0.06), log L150 (L⊙) = 1.31(±0.05) × log SFR60 (M⊙ yr−1) + 3.14(±0.06), and log L150 (L⊙) = 1.37 (±0.05) × log SFRIR (M⊙ yr−1) + 3.09(±0.05), which show excellent agreement with each other.

Resolved Molecular Gas and Star Formation Properties of the Strongly Lensed z = 2.26 Galaxy SDSS J0901+1814

We present ~1" resolution (~2 kpc in the source plane) observations of the CO(1-0), CO(3-2), Halpha, and [N II] lines in the strongly-lensed z=2.26 star-forming galaxy SDSS J0901+1814. We use these observations to constrain the lensing potential of a foreground group of galaxies, and our source-plane reconstructions indicate that SDSS J0901+1814 is a nearly face-on (i~30 degrees) massive disk with r_{1/2}>~4 kpc for its molecular gas. Using our new magnification factors (mu_tot~30), we find that SDSS J0901+1814 has a star formation rate (SFR) of 268^{+63}_{-61} M_sun/yr, M_gas=(1.6^{+0.3}_{-0.2})x10^11x(alpha_CO/4.6) M_sun, and M_star=(9.5^{+3.8}_{-2.8})x10^10 M_sun, which places it on the star-forming galaxy "main sequence." We use our matched high-angular resolution gas and SFR tracers (CO and Halpha, respectively) to perform a spatially resolved (pixel-by-pixel) analysis of SDSS J0901+1814 in terms of the Schmidt-Kennicutt relation. After correcting for the large fraction of obscured star formation (SFR_Halpha/SFR_TIR=0.054^{+0.015}_{-0.014}), we find SDSS J0901+1814 is offset from "normal" star-forming galaxies to higher star formation efficiencies independent of assumptions for the CO-to-H_2 conversion factor. Our mean best-fit index for the Schmidt-Kennicutt relation for SDSS J0901+1814, evaluated with different CO lines and smoothing levels, is n=1.54+/-0.13; however, the index may be affected by gravitational lensing, and we find n=1.24+/-0.02 when analyzing the source-plane reconstructions. While the Schmidt-Kennicutt index largely appears unaffected by which of the two CO transitions we use to trace the molecular gas, the source-plane reconstructions and dynamical modeling suggest that the CO(1-0) emission is more spatially extended than the CO(3-2) emission.

Dense gas and star formation in individual Giant Molecular Clouds in M31

Studies both of entire galaxies and of local Galactic star formation indicate a dependency of a molecular cloud's star formation rate (SFR) on its dense gas mass. In external galaxies, such measurements are derived from HCN(1-0) observations, usually encompassing many Giant Molecular Clouds (GMCs) at once. The Andromeda galaxy (M31) is a unique laboratory to study the relation of the SFR and HCN emission down to GMC scales at solar-like metallicities. In this work, we correlate our composite SFR determinations with archival HCN, HCO$^+$, and CO observations, resulting in a sample of nine reasonably representative GMCs. We find that, at the scale of individual clouds, it is important to take into account both obscured and unobscured star formation to determine the SFR. When correlated against the dense-gas mass from HCN, we find that the SFR is low, in spite of these refinements. We nevertheless retrieve an SFR - dense-gas mass correlation, confirming that these SFR tracers are still meaningful on GMC scales. The correlation improves markedly when we consider the HCN/CO ratio instead of HCN by itself. This nominally indicates a dependency of the SFR on the dense-gas fraction, in contradiction to local studies. However, we hypothesize that this partly reflects the limited dynamic range in dense-gas mass, and partly that the ratio of single-pointing HCN and CO measurements may be less prone to systematics like sidelobes. In this case, the HCN/CO ratio would importantly be a better empirical measure of the dense-gas content itself.

Populating H2 and CO in galaxy simulation with dust evolution

There are two major theoretical issues for the star formation law (the relation between the surface densities of molecular gas and star formation rate on a galaxy scale): (i) At low metallicity, it is not obvious that star-forming regions are rich in H$_2$ because the H$_2$ formation rate depends on the dust abundance; and (ii) whether or not CO really traces H$_2$ is uncertain, especially at low metallicity. To clarify these issues, we use a hydrodynamic simulation of an isolated disc galaxy with a spatial resolution of a few tens parsecs. The evolution of dust abundance and grain size distribution is treated consistently with the metal enrichment and the physical state of the interstellar medium. We compute the H$_2$ and CO abundances using a subgrid post-processing model based on the dust abundance and the dissociating radiation field calculated in the simulation. We find that when the metallicity is $\lesssim 0.4$ Z$_\odot$ ($t<1$ Gyr), H$_2$ is not a good tracer of star formation rate because H$_2$-rich regions are limited to dense compact regions. At $Z\gtrsim 0.8$ Z$_\odot$, a tight star formation law is established for both H$_2$ and CO. At old ($t \sim 10$ Gyr) ages, we also find that adopting the so-called MRN grain size distribution with an appropriate dust-to-metal ratio over the entire disc gives reasonable estimates for the H$_2$ and CO abundances. For CO, improving the spatial resolution of the simulation is important while the H$_2$ abundance is not sensitive to sub-resolution structures at $Z\gtrsim 0.4$ Z$_\odot$.

The LOFAR window on star-forming galaxies and AGNs – curved radio SEDs and IR–radio correlation at 0&lt;z&lt;2.5

We present a study of the low-frequency radio properties of star forming (SF) galaxies and active galactic nuclei (AGN) up to redshift $z=2.5$. The new spectral window probed by the Low Frequency Array (LOFAR) allows us to reconstruct the radio continuum emission from 150 MHz to 1.4 GHz to an unprecedented depth for a radio-selected sample of $1542$ galaxies in $\sim 7~ \rm{deg}^2$ of the LOFAR Bo\"otes field. Using the extensive multi-wavelength dataset available in Bo\"otes and detailed modelling of the FIR to UV spectral energy distribution (SED), we are able to separate the star-formation (N=758) and the AGN (N=784) dominated populations. We study the shape of the radio SEDs and their evolution across cosmic time and find significant differences in the spectral curvature between the SF galaxy and AGN populations. While the radio spectra of SF galaxies exhibit a weak but statistically significant flattening, AGN SEDs show a clear trend to become steeper towards lower frequencies. No evolution of the spectral curvature as a function of redshift is found for SF galaxies or AGN. We investigate the redshift evolution of the infrared-radio correlation (IRC) for SF galaxies and find that the ratio of total infrared to 1.4 GHz radio luminosities decreases with increasing redshift: $ q_{\rm 1.4GHz} = (2.45 \pm 0.04) \times (1+z)^{-0.15 \pm 0.03} $. Similarly, $q_{\rm 150MHz}$ shows a redshift evolution following $ q_{\rm 150GHz} = (1.72 \pm 0.04) \times (1+z)^{-0.22 \pm 0.05}$. Calibration of the 150 MHz radio luminosity as a star formation rate tracer suggests that a single power-law extrapolation from $q_{\rm 1.4GHz}$ is not an accurate approximation at all redshifts.

Observations of M31 and M33 with the Fermi Large Area Telescope: A Galactic Center Excess in Andromeda?

Abstract The Fermi Large Area Telescope (LAT) has opened the way for comparative studies of cosmic rays (CRs) and high-energy objects in the Milky Way (MW) and in other, external, star-forming galaxies. Using 2 yr of observations with the Fermi LAT, Local Group galaxy M31 was detected as a marginally extended gamma-ray source, while only an upper limit has been derived for the other nearby galaxy M33. We revisited the gamma-ray emission in the direction of M31 and M33 using more than 7 yr of LAT Pass 8 data in the energy range 0.1 – 100 GeV , presenting detailed morphological and spectral analyses. M33 remains undetected, and we computed an upper limit of 2.0 × 10 − 12 erg cm − 2 s − 1 on the 0.1 – 100 GeV energy flux (95% confidence level). This revised upper limit remains consistent with the observed correlation between gamma-ray luminosity and star formation rate tracers and implies an average CR density in M33 that is at most half of that of the MW. M31 is detected with a significance of nearly 10 σ . Its spectrum is consistent with a power law with photon index Γ = 2.4 ± 0.1 stat + syst and a 0.1 – 100 GeV energy flux of ( 5.6 ± 0.6 stat + syst ) × 10 − 12 erg cm − 2 s − 1 . M31 is detected to be extended with a 4 σ significance. The spatial distribution of the emission is consistent with a uniform-brightness disk with a radius of 0.°4 and no offset from the center of the galaxy, but nonuniform intensity distributions cannot be excluded. The flux from M31 appears confined to the inner regions of the galaxy and does not fill the disk of the galaxy or extend far from it. The gamma-ray signal is not correlated with regions rich in gas or star formation activity, which suggests that the emission is not interstellar in origin, unless the energetic particles radiating in gamma rays do not originate in recent star formation. Alternative and nonexclusive interpretations are that the emission results from a population of millisecond pulsars dispersed in the bulge and disk of M31 by disrupted globular clusters or from the decay or annihilation of dark matter particles, similar to what has been proposed to account for the so-called Galactic center excess found in Fermi-LAT observations of the MW.

GAMA/H-ATLAS: common star formation rate indicators and their dependence on galaxy physical parameters

We compare common star-formation rate (SFR) indicators in the local Universe in the GAMA equatorial fields (around 160 sq. deg.), using ultraviolet (UV) photometry from GALEX, far-infrared (FIR) and sub-millimetre (sub-mm) photometry from H-ATLAS, and Halpha spectroscopy from the GAMA survey. With a high-quality sample of 745 galaxies (median redshift 0.08), we consider three SFR tracers: UV luminosity corrected for dust attenuation using the UV spectral slope beta (SFRUV,corr), Halpha line luminosity corrected for dust using the Balmer decrement (BD) (SFRHalpha,corr), and the combination of UV and IR emission (SFRUV+IR). We demonstrate that SFRUV,corr can be reconciled with the other two tracers after applying attenuation corrections by calibrating IRX (i.e. the IR to UV luminosity ratio) and attenuation in the Halpha (derived from BD) against beta. However, beta on its own is very unlikely to be a reliable attenuation indicator. We find that attenuation correction factors depend on parameters such as stellar mass, z and dust temperature (Tdust), but not on Halpha equivalent width (EW) or Sersic index. Due to the large scatter in the IRX vs beta correlation, when compared to SFRUV+IR, the beta-corrected SFRUV,corr exhibits systematic deviations as a function of IRX, BD and Tdust.

Molecular gas and star formation in the tidal dwarf galaxy VCC 2062

The physical mechanisms driving star formation (SF) in galaxies are still not fully understood. Tidal dwarf galaxies (TDGs), made of gas ejected during galaxy interactions, seem to be devoid of dark matter and have a near-solar metallicity. The latter makes it possible to study molecular gas and its link to SF using standard tracers (CO, dust) in a peculiar environment. We present a detailed study of a nearby TDG in the Virgo Cluster, VCC 2062, using new high-resolution CO(1--0) data from the Plateau de Bure, deep optical imaging from the Next Generation Virgo Cluster Survey (NGVS), and complementary multiwavelength data. Until now, there was some doubt whether VCC 2062 was a true TDG, but the new deep optical images from the NGVS reveal a stellar bridge between VCC 2062 and its parent galaxy, NGC 4694, which is clear proof of its tidal origin. Several high-resolution tracers (\halpha, UV, 8~\mi, and 24~\mi) of the star formation rate (SFR) are compared to the molecular gas distribution as traced by the CO(1-0). Coupled with the SFR tracers, the NGVS data are used with the CIGALE code to model the stellar populations throughout VCC 2062, yielding a declining SFR in the recent past, consistent with the low \halpha/UV ratio, and a high burst strength. HI emission covers VCC 2062, whereas the CO is concentrated near the HI maxima. The CO peaks correspond to two very distinct regions: one with moderate SF to the NE and one with only slightly weaker CO emission but with nearly no SF. Even where SF is clearly present, the SFR is below the value expected from the surface density of the molecular and the total gas as compared to spiral galaxies and other TDGs. After discussing different possible explanations, we conclude that the low surface brightness is a crucial parameter to understand the low SFR.

Spatially-resolved SFR in nearby disk galaxies using IFS data

AbstractExploring the spatial distribution of the star formation rate (SFR) in nearby galaxies is essential to understand their evolution through cosmic time. With this aim in mind, we use a representative sample that contains a variety of morphological types, the CALIFA Integral Field Spectroscopy (IFS) sample. Previous to this work, we have verified that our extinction-corrected Hα measurements successfully reproduce the values derived from other SFR tracers such as Hαobs + IR or UVobs + IR (Catalán-Torrecilla et al. 2015).Now, we go one step further applying 2-dimensional photometric decompositions (Méndez-Abreu et al. (2008), Méndez-Abreu et al. (2014)) over these datacubes. This method allows us to obtain the amount of SFR in the central part (bulge or nuclear source), the bar and the disk, separately. First, we determine the light coming from each component as the ratio between the luminosity in every component (bulge, bar or disk) and the total luminosity of the galaxy. Then, for each galaxy we multiply the IFS datacubes by these previous factors to recover the luminosity in each component. Finally, we derive the spectrum associated to each galaxy component integrating the spatial information in the weighted datacube using an elliptical aperture covering the whole galaxy.2D photometric decomposition applied over 3D datacubes will give us a more detailed understanding of the role that disks play in more massive galaxies. Knowing if the disks in more massive SF galaxies have on average a lower or higher level of star formation activity and how these results are affected by the presence of nuclear bars are still open questions that we can now solve. We describe the behavior of these components in the SFR vs. stellar mass diagram. In particular, we highlight the role of the disks and their contribution to both the integrated SFR for the whole galaxy and the SFR in the disk at different stellar masses in the SFR vs. stellar mass diagram together with their relative position to the star forming Main Sequence.

AGN feedback and star formation in young and old radio galaxies

AbstractPowerful radio galaxies show evidence of ongoing active galactic nuclei (AGN) feedback, mainly in the form of fast, massive outflows. Yet, it is not clear how these outflows affect the star formation of their hosts. We investigated the different manifestations of AGN feedback in the evolved 3C 293 radio source, and in the young, reactivated 3C 236. Both sources harbor young star‐forming regions and fast outflows of H I and ionized gas. Due to the different evolution stages of these sources, they are expected to be at different phases of the AGNISM interaction process. Using high spatial resolution observations of CO lines taken with the IRAM Plateau de Bure interferometer, we studied the distribution and kinematics of the molecular gas in these sources, and compared the results with the dust and star‐formation images of the host. We searched for signatures of outflow motions in the CO kinematics, and derived the star formation rate (SFR) and star formation efficiency (SFE) of the host with all available SFR tracers. Based on the gas mass derived from our observations and on the SFR estimates, we compared the star‐formation efficiency of 3C 236 and 3C 293 against a sample of powerful radio galaxies. Our results suggest that the apparently low SFE of evolved radio galaxies is not necessarily due to AGN feedback, but may be caused by an underestimation of the SFR and/or an overestimation of the molecular gas densities in radio galaxies. (© 2016 WILEY‐VCH Verlag GmbH &amp; Co. KGaA, Weinheim)

Star formation laws in the Andromeda galaxy: gas, stars, metals and the surface density of star formation

We use hierarchical Bayesian regression analysis to investigate star formation laws in the Andromeda galaxy (M31) in both local (30, 155, and 750pc) and global cases. We study and compare the well-known Kennicutt-Schmidt law, the extended Schmidt law and the metallicity/star formation correlation. Using a combination of H$\alpha$ and 24 $\mu$m emission, a combination of far-ultraviolet and 24$\mu$m, and the total infrared emission, we estimate the total star formation rate (SFR) in M31 to be between $0.35\pm 0.04$M$_{\odot}$yr$^{-1}$ and $0.4\pm 0.04$M$_{\odot}$yr$^{-1}$. We produce a stellar mass surface density map using IRAC 3.6$\mu$m emission and measured the total stellar mass to be $6.9 \times 10^{10}$M$_{\odot}$. For the Kennicutt-Schmidt law in M31, we find the power-law index $N$ to be between 0.49 and 1.18, for all the laws, the power-law index varies more with changing gas tracer than with SFR tracer. The power-law index also changes with distance from the centre of the galaxy. We also applied the commonly-used ordinary least squares fitting method and showed that using different fitting methods leads to different power-law indices. There is a correlation between the surface density of SFR and the stellar mass surface density, which confirms that the Kennicutt-Schmidt law needs to be extended to consider the other physical properties of galaxies. We found a weak correlation between metallicity, the SFR and the stellar mass surface density.

Gamma-Ray Bursts: A Radio Perspective

Gamma-ray bursts (GRBs) are extremely energetic events at cosmological distances. They provide unique laboratory to investigate fundamental physical processes under extreme conditions. Due to extreme luminosities, GRBs are detectable at very high redshifts and potential tracers of cosmic star formation rate at early epoch. While the launch ofSwiftandFermihas increased our understanding of GRBs tremendously, many new questions have opened up. Radio observations of GRBs uniquely probe the energetics and environments of the explosion. However, currently only 30% of the bursts are detected in radio bands. Radio observations with upcoming sensitive telescopes will potentially increase the sample size significantly and allow one to follow the individual bursts for a much longer duration and be able to answer some of the important issues related to true calorimetry, reverse shock emission, and environments around the massive stars exploding as GRBs in the early Universe.

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.

Determining the fraction of dust heating from young and old stellar populations with 3D dust radiative transfer

AbstractA major difficulty hampering the accuracy of UV/optical star formation rate tracers is the effect of interstellar dust, absorbing and scattering light produced by both young and old stellar populations (SPs). Although empirically calibrated corrections or energy balance SED fitting are often used for fast de-reddening of galaxy stellar emission, eventually only radiative transfer calculations can provide self-consistent predictions of galaxy model spectra, taking into account important factors such as galaxy inclination, different morphological components, non-local heating of the dust and scattered radiation. In addition, dust radiative transfer can be used to determine the fraction of monochromatic dust emission powered by either young or old SPs. This calculation needs to take into account the different response of the dust grains to the UV and optical radiation field, depending on the grain size and composition. We determined the dust heating fractions, on both global and local scales, for a high-resolution galaxy model by using our 3D ray-tracing dust radiative transfer code “DART-Ray”. We show the results obtained using this method and discuss the consequences for star formation rate indicators.

Are LGRBs biased tracers of star formation? Clues from the host galaxies of the Swift/BAT6 complete sample of LGRBs

AbstractLong gamma-ray bursts (LGRBs) are associated to the deaths of massive stars and could thus be used as a potentially powerful tool to trace cosmic star formation. However the conditions needed to produce a LGRBs may introduce a bias in the LGRB rate versus star formation rate (SFR) relation (called LGRB efficiency hereafter).We have undertaken a study of the properties of the host galaxies of the BAT6 complete sample of LGRB to improve our knowledge on the LGRB efficiency, its redshift evolution, and the factor affecting it. This is the base to properly use LGRBs as SFR tracers.We show that at z &lt; 1 LGRBs are not direct SFR tracers because they tend to avoid high-metallicity galaxies. The use of the BAT6 complete sample keeps this result from being affected by possible biases that could have influenced past results based on incomplete samples. The preference for low (but not extremely low) metallicities can be a consequence of the particular conditions needed for the progenitor star to produce a GRB.

AN UNEXPECTEDLY LOW-REDSHIFT EXCESS OF SWIFT GAMMA-RAY BURST RATE

Gamma-ray bursts (GRBs) are the most violent explosions in the Universe and can be used to explore the properties of high-redshift universe. It is believed that the long GRBs are associated with the deaths of massive stars. So it is possible to use GRBs to investigate the star formation rate (SFR). In this paper, we use Lynden-Bell's $c^-$ method to study the luminosity function and rate of \emph{Swift} long GRBs without any assumptions. We find that the luminosity of GRBs evolves with redshift as $L(z)\propto g(z)=(1+z)^k$ with $k=2.43_{-0.38}^{+0.41}$. After correcting the redshift evolution through $L_0(z)=L(z)/g(z)$, the luminosity function can be expressed as $\psi(L_0)\propto L_0^{-0.14\pm0.02}$ for dim GRBs and $\psi(L_0)\propto L_0^{-0.70\pm0.03}$ for bright GRBs, with the break point $L_{0}^{b}=1.43\times10^{51}~{\rm erg~s^{-1}}$. We also find that the formation rate of GRBs is almost constant at $z<1.0$ for the first time, which is remarkably different from the SFR. At $z>1.0$, the formation rate of GRB is consistent with the SFR. Our results are dramatically different from previous studies. Some possible reasons for this low-redshift excess are discussed. We also test the robustness of our results with Monte Carlo simulations. The distributions of mock data (i.e., luminosity-redshift distribution, luminosity function, cumulative distribution and $\log N-\log S$ distribution) are in good agreement with the observations. Besides, we also find that there are remarkable difference between the mock data and the observations if long GRB are unbiased tracers of SFR at $z<1.0$.