In this work we investigate the far-UV to NIR shape of the dust attenuation curve of a sample of IR selected dust obscured (U)LIRGs at z$\sim$2. The spectral energy distributions (SEDs) are fitted with CIGALE, a physically-motivated spectral synthesis model based on energy balance. Its flexibility allows us to test a wide range of different analytical prescriptions for the dust attenuation curve, including the well-known Calzetti and Charlot & Fall curves, and modified versions of them. The attenuation curves computed under the assumption of our reference double power-law model are in very good agreement with those derived, in previous works, with radiative transfer (RT) SED fitting. We investigate the position of our galaxies in the IRX-$\beta$ diagram and find this to be consistent with grayer slopes, on average, in the UV. We also find evidence for a flattening of the attenuation curve in the NIR with respect to more classical Calzetti-like recipes. This larger NIR attenuation yields larger derived stellar masses from SED fitting, by a median factor of $\sim$ 1.4 and up to a factor $\sim$10 for the most extreme cases. The star formation rate appears instead to be more dependent on the total amount of attenuation in the galaxy. Our analysis highlights the need for a flexible attenuation curve when reproducing the physical properties of a large variety of objects.

We explore from a statistical point of view the far-infrared (far-IR) and sub-millimeter (sub-mm) properties of a large sample of LBGs (22,000) at z~3 in the COSMOS field. The large number of galaxies allows us to split it in several bins as a function of UV luminosity, UV slope, and stellar mass to better sample their variety. We perform stacking analysis in PACS (100 and 160 um), SPIRE (250, 350 and 500 um) and AzTEC (1.1 mm) images. Our stacking procedure corrects the biases induced by galaxy clustering and incompleteness of our input catalogue in dense regions. We obtain the full IR spectral energy distributions (SED) of subsamples of LBGs and derive the mean IR luminosity as a function of UV luminosity, UV slope, and stellar mass. The average IRX is roughly constant over the UV luminosity range, with a mean of 7.9 (1.8 mag). However, it is correlated with UV slope, and stellar mass. We investigate using a statistically-controlled stacking analysis as a function of (stellar mass, UV slope) the dispersion of the IRX-UVslope and IRX-M* plane. Our results enable us to study the average relation between star-formation rate (SFR) and stellar mass, and we show that our LBG sample lies on the main sequence of star formation at z~3.

We describe the search for Lyman-break galaxies (LBGs) near the sub-millimeter bright starburst galaxy HFLS3 at $z$$=$6.34 and a study on the environment of this massive galaxy during the end of reionization.We performed two independent selections of LBGs on images obtained with the \textit{Gran Telescopio Canarias} (GTC) and the \textit{Hubble Space Telescope} (HST) by combining non-detections in bands blueward of the Lyman-break and color selection. A total of 10 objects fulfilling the LBG selection criteria at $z$$>$5.5 were selected over the 4.54 and 55.5 arcmin$^2$ covered by our HST and GTC images, respectively. The photometric redshift, UV luminosity, and the star-formation rate of these sources were estimated with models of their spectral energy distribution. These $z$$\sim$6 candidates have physical properties and number densities in agreement with previous results. The UV luminosity function at $z$$\sim$6 and a Voronoi tessellation analysis of this field shows no strong evidence for an overdensity of relatively bright objects (m$_{F105W}$$<$25.9) associated with \textit{HFLS3}. However, the over-density parameter deduced from this field and the surface density of objects can not excluded definitively the LBG over-density hypothesis. Moreover we identified three faint objects at less than three arcseconds from \textit{HFLS3} with color consistent with those expected for $z$$\sim$6 galaxies. Deeper data are needed to confirm their redshifts and to study their association with \textit{HFLS3} and the galaxy merger that may be responsible for the massive starburst.

We complement our previous analysis of a sample of z~1-2 luminous and ultra-luminous infrared galaxies ((U)LIRGs), by adding deep VLA radio observations at 1.4 GHz to a large data-set from the far-UV to the sub-mm, including Spitzer and Herschel data. Given the relatively small number of (U)LIRGs in our sample with high S/N radio data, and to extend our study to a different family of galaxies, we also include 6 well sampled near IR-selected BzK galaxies at z~1.5. From our analysis based on the radiative transfer spectral synthesis code GRASIL, we find that, while the IR luminosity may be a biased tracer of the star formation rate (SFR) depending on the age of stars dominating the dust heating, the inclusion of the radio flux offers significantly tighter constraints on SFR. Our predicted SFRs are in good agreement with the estimates based on rest-frame radio luminosity and the Bell (2003) calibration. The extensive spectro-photometric coverage of our sample allows us to set important constraints on the SF history of individual objects. For essentially all galaxies we find evidence for a rather continuous SFR and a peak epoch of SF preceding that of the observation by a few Gyrs. This seems to correspond to a formation redshift of z~5-6. We finally show that our physical analysis may affect the interpretation of the SFR-M* diagram, by possibly shifting, with respect to previous works, the position of the most dust obscured objects to higher M* and lower SFRs.

Long-duration gamma-ray bursts (LGRBs) are produced by the collapse of very massive stars. Due to the short lifetime of their progenitors, LGRBs pinpoint star-forming galaxies. We present here a multi-band search for the host galaxy of the long dark GRB 050219A within the enhanced Swift/XRT error circle. We used spectroscopic observations acquired with VLT/X-shooter to determine the redshift and star-formation rate of the putative host galaxy. We compared the results with the optical/IR spectral energy distribution obtained with different facilities. Surprisingly, the host galaxy is a old and quiescent early-type galaxy at z = 0.211 characterised by an unprecedentedly low specific star-formation rate. It is the first LGRB host to be also an early-type post-starburst galaxy. This is further evidence that GRBs can explode in all kind of galaxies, with the only requirement being an episode of star-formation.

In this work we explore the impact of the presence of an active galactic nucleus (AGN) on the mid- and far-infrared (IR) properties of galaxies as well as the effects of simultaneous AGN and starburst activity in these same galaxies. To do this we apply a multi-component, multi-band spectral synthesis technique to a sample of 250 micron selected galaxies of the Herschel Multi-tiered Extragalactic Survey, with IRS spectra available for all galaxies. Our results confirm that the inclusion of the IRS spectra plays a crucial role in the spectral analysis of galaxies with an AGN component improving the selection of the best-fit hot dust model (torus). We find a correlation between the obscured star formation rate (SFR) derived from the IR luminosity of the starburst component, SFR_IR and SFR_PAH, derived from the luminosity of the PAH features, L_PAH, with SFR_FIR taking higher values than SFR_PAH. The correlation is different for AGN- and starburst-dominated objects. The ratio of L_PAH to that of the starburst component, L_PAH/L_SB, is almost constant for AGN-dominated objects but decreases with increasing L_SB for starburst-dominated objects. SFR_FIR increases with the accretion luminosity, L_acc, with the increase less prominent for the very brightest, unobscured AGN-dominated sources. We find no correlation between the masses of the hot and cold dust components. We interpret this as a non-constant fraction of gas driven by the gravitational effects to the AGN while the starburst is ongoing. We also find no evidence of the AGN affecting the temperature of the cold dust component, though this conclusion is mostly based on objects with a non-dominant AGN component. We conclude that our findings do not provide evidence that the presence of AGN affects the star formation process in the host galaxy, but rather that the two phenomena occur simultaneously over a wide range of luminosities.

We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $\mu$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $\Delta$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $\Delta$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

We computed a comprehensive set of theoretical ultraviolet spectra of hot, massive stars with the radiation-hydrodynamics code WM-Basic. This model atmosphere and spectral synthesis code is optimized for computing the strong P Cygni-type lines originating in the winds of hot stars, which are the strongest features in the ultraviolet spectral region. The computed set is suitable as a spectral library for inclusion in evolutionary synthesis models of star clusters and star-forming galaxies. The chosen stellar parameters cover the upper left Hertzsprung-Russell diagram at L >~ 10^2.75 Lsun and T_eff >~ 20,000 K. The adopted elemental abundances are 0.05 Zsun, 0.2 Zsun, 0.4 Zsun, Zsun, and 2 Zsun. The spectra cover the wavelength range from 900 to 3000 {\AA} and have a resolution of 0.4 {\AA}. We compared the theoretical spectra to data of individual hot stars in the Galaxy and the Magellanic Clouds obtained with the International Ultraviolet Explorer (IUE) and Far Ultraviolet Spectroscopic Explorer (FUSE) satellites and found very good agreement. We built a library with the set of spectra and implemented it into the evolutionary synthesis code Starburst99 where it complements and extends the existing empirical library towards lower chemical abundances. Comparison of population synthesis models at solar and near-solar composition demonstrates consistency between synthetic spectra generated with either library. We discuss the potential of the new library for the interpretation of the rest-frame ultraviolet spectra of star-forming galaxies. Properties that can be addressed with the models include ages, initial mass function, and heavy-element abundance. The library can be obtained both individually or as part of the Starburst99 package.

It has recently been shown that galaxy formation models within the LambdaCDM cosmology predict that, compared to the observed population, small galaxies (with stellar masses < 10^{11} M_sun) form too early, are too passive since z ~ 3 and host too old stellar populations at z=0. We then expect an overproduction of small galaxies at z > 4 that should be visible as an excess of faint Lyman-break galaxies. To check whether this excess is present, we use the MORGANA galaxy formation model and GRASIL spectro-photometric + radiative transfer code to generate mock catalogues of deep fields observed with HST-ACS. We add observational noise and the effect of Lyman-alpha emission, and perform color-color selections to identify Lyman-break galaxies. The resulting mock candidates have plausible properties that closely resemble those of observed galaxies. We are able to reproduce the evolution of the bright tail of the luminosity function of Lyman-break galaxies (with a possible underestimate of the number of the brightest i-dropouts), but uncertainties and degeneracies in dust absorption parameters do not allow to give strong constraints to the model. Besides, our model shows a clear excess with respect to observations of faint Lyman-break galaxies, especially of z_{850} ~ 27 V-dropouts at z ~ 5. We quantify the properties of these "excess" galaxies and discuss the implications: these galaxies are hosted in dark matter halos with circular velocities in excess of 100 km s^{-1}, and their suppression may require a deep re-thinking of stellar feedback processes taking place in galaxy formation.