Researcher profile

F. Calura

F. Calura contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
43works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

43 published item(s)

preprint2026arXiv

A Pristine Star-Forming Complex at z=4.19

We report the discovery of a faint (M_1700 ~ -12.2), oxygen-deficient strongly-lensed ionizing source -- dubbed LAP2 (Lensed And Pristine 2) -- at a spectroscopic redshift of z=4.19. LAP2 appears to be isolated and lies very close to the caustic produced by the lensing galaxy cluster Abell 2744. It was observed with the James Webb Space Telescope (JWST) NIRSpec MSA in prism mode as part of the UNCOVER program. The NIRSpec spectrum reveals prominent Lya emission (7.1 sigma), clear Ha emission (6.2 sigma), tentative Hb emission (2.8 sigma), and no detectable [OIII]4959,5007 (~ 7 times fainter than Ha). The inferred [OIII] 2 sigma upper limit corresponds to an R3 index <0.85 (assuming the Ha/Hb ~ 2.86 case~B recombination ratio), which, under high-ionization conditions, implies a metallicity of Z < 6 X 10^(-3) Z_sun. The combination of faint ultraviolet luminosity, large rest-frame Ha equivalent width (~ 650 A), and extremely compact size (< 10 pc) suggests that LAP2 is being caught in an early, pristine formation phase consistent with an instantaneous-burst scenario, with an estimated stellar mass of at most a few 10^4 Msun. Deep VLT/MUSE observations further reveal copious Lya emission forming an arclet that straddles the critical line. LAP2 thus joins the rare class of extremely metal-poor star-forming complexes that JWST has started to unveil at redshifts 3 - 7, and provides a rare glimpse into a still very poorly explored faint-luminosity regime.

preprint2022arXiv

A new estimate of the cosmic star formation density from a radio-selected sample, and the contribution of $H$-dark galaxies at $z \geq 3$

The Star Formation Rate Density (SFRD) history of the Universe is well constrained up to redshift $z \sim 2$. At earlier cosmic epochs, the picture has been largely inferred from UV-selected galaxies (e.g. Lyman-break galaxies, LBGs). However, LBGs&#39; inferred SFRs strongly depend on the assumed dust extinction correction, which is not well-constrained at high-$z$, while observations in the radio domain are not affected by this issue. In this work we measure the SFRD from a 1.4 GHz-selected sample of $\sim$600 galaxies in the GOODS-N field up to redshift $\sim 3.5$. We take into account the contribution of Active Galactic Nuclei from the Infrared-Radio correlation. We measure the radio luminosity function, fitted with a modified Schechter function, and derive the SFRD. The cosmic SFRD shows a rise up to $z \sim 2$ and then an almost flat plateau up to $z \sim 3.5$. Our SFRD is in agreement with the ones from other FIR/radio surveys and a factor 2 higher than those from LBG samples. We also estimate that galaxies lacking a counterpart in the HST/WFC3 H-band ($H$-dark) make up $\sim 25\%$ of the $ϕ$-integrated SFRD relative to the full sample at z $\sim 3.2$, and up to $58\%$ relative to LBG samples.

preprint2022arXiv

Exploring the physical properties of lensed star-forming clumps at $2\lesssim z \lesssim6$

We study the physical properties (size, stellar mass, luminosity, star formation rate) and scaling relations for a sample of 166 star-forming clumps with redshift $z \sim 2-6.2$. They are magnified by the Hubble Frontier Field galaxy cluster MACS~J0416 and have robust lensing magnification ($2\lesssim μ\lesssim 82$) computed by using our high-precision lens model, based on 182 multiple images. Our sample extends by $\sim 3$ times the number of spectroscopically-confirmed lensed clumps at $z \gtrsim 2$. We identify clumps in ultraviolet continuum images and find that, whenever the effective spatial resolution (enhanced by gravitational lensing) increases, they fragment into smaller entities, likely reflecting the hierarchically-organized nature of star formation. Kpc-scale clumps, most commonly observed in field, are not found in our sample. The physical properties of our sample extend the parameter space typically probed by $z \gtrsim 1$ field observations and simulations, by populating the low mass (M$_\star \lesssim 10^7$ M$_\odot$), low star formation rate (SFR $\lesssim 0.5$ M$_\odot$ yr$^{-1}$), and small size (R$_\mathrm{eff} \lesssim 100$ pc) regime. The new domain probed by our study approaches the regime of compact stellar complexes and star clusters. In the mass-size plane, our sample spans the region between galaxies and globular clusters, with a few clumps in the region populated by young star clusters and globular-clusters. For the bulk of our sample, we measure star-formation rates which are higher than those observed locally in compact stellar systems, indicating different conditions for star formation at high redshift than in the local Universe.

preprint2022arXiv

Sub-parsec resolution cosmological simulations of star-forming clumps at high redshift with feedback of individual stars

We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z=6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified &#39;delayed cooling&#39; feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n~10^3-10^5 cm^{-3}). The sites where star formation ignites are characterised by maximum densities of the order of 10^5 cm^{-3} and gravitational pressures P/k>10^7 K/cm^3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z=6.14 is 3.4x10^7 M_sun, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of ~10^6 M_sun and half-mass sizes of ~100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one order of magnitude lower than the observed one. In the size-mass plane, our clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies.

preprint2022arXiv

The evolution of the Si IV content in the Universe from the epoch of reionization to cosmic noon

We investigate the abundance and distribution of metals in the high-redshift intergalactic medium and circum-galactic medium through the analysis of a sample of almost 600 SiIV absorption lines detected in high and intermediate resolution spectra of 147 quasars. The evolution of the number density of SiIV lines, the column density distribution function and the cosmic mass density are studied in the redshift interval 1.7 <= z <= 6.2 and for log N(SiIV) >= 12.5. All quantities show a rapid increase between z~6 and z< 5 and then an almost constant behaviour to z~2 in very good agreement with what is already observed for CIV absorption lines. The present results are challenging for numerical simulations: when simulations reproduce our SiIV results, they tend to underpredict the properties of CIV, and when the properties of CIV are reproduced, the number of strong SiIV lines (log N(SiIV) > 14) is overpredicted.

preprint2021arXiv

APOGEE DR16: a multi-zone chemical evolution model for the Galactic disc based on MCMC methods

The analysis of the APOGEE DR16 data suggests the existence of a clear distinction between two sequences of disc stars at different Galactocentric distances in the [$α$/Fe] vs. [Fe/H] abundance ratio space: the so-called high-$α$ sequence, classically associated to an old population of stars in the thick disc, and the low-$α$ sequence, which mostly comprises relatively young stars in the thin disc. We perform a Bayesian analysis based on a Markov Chain Monte Carlo method to constrain a multi-zone two-infall chemical evolution model designed for regions at different Galactocentric distances using measured chemical abundances from the APOGEE DR16 sample. An inside-out formation of the Galaxy disc naturally emerges from the best fit of our two-infall chemical-evolution model to APOGEE-DR16: inner Galactic regions are assembled on shorter time-scales compared to the external ones. In the outer disc (with radii $R>6$ kpc), the chemical dilution due to a late accretion event of gas with primordial chemical composition is the main driver of the [Mg/Fe] vs. [Fe/H] abundance pattern in the low-$α$ sequence. In the inner disc, in the framework of the two-infall model, we confirm the presence of an enriched gas infall in the low-$α$ phase as suggested by chemo-dynamical models. Our Bayesian analysis of the recent APOGEE DR16 data suggests a significant delay time, ranging from $\sim$3.0 to 4.7 Gyr, between the first and second gas infall events for all the analyzed Galactocentric regions. Our results propose a clear interpretation of the [Mg/Fe] vs. [Fe/H] relations along the Galactic discs. The signatures of a delayed gas-rich merger which gives rise to a hiatus in the star formation history of the Galaxy are impressed in the [Mg/Fe] vs. [Fe/H] relation, determining how the low-$α$ stars are distributed in the abundance space at different Galactocentric distances.

preprint2021arXiv

High star cluster formation efficiency in the strongly lensed Sunburst Lyman-continuum galaxy at z=2.37

We investigate the strongly lensed (μx10-100) Lyman continuum (LyC) galaxy, dubbed Sunburst, at z=2.37, taking advantage of a new accurate model of the lens. A characterization of the intrinsic (delensed) properties of the galaxy yields a size of ~3 sq.kpc, a luminosity Muv=-20.3,and a stellar mass M~10^9 Msun;16% of the ultraviolet light is located in a 3 Myr old gravitationally-bound young massive star cluster (YMC) with an effective radius of Re~8 pc and a dynamical mass of ~10^7 Msun (similar to the stellar mass), from which LyC radiation is detected (λ< 912A). The inferred outflowing gas velocity (>300 km/s) exceeds the escape velocity of the star cluster. The resulting escape fraction of the ionizing radiation emerging from the Sunburst galaxy is >6-12%, whilst it is >46-93% if inferred from the YMC. 12 additional likely star clusters with 3<Re<20 pc are identified in the galaxy from which we derive a cluster formation efficiency Γ>~30%, which is consistent with the high Γderived in local galaxies experiencing extreme gas physical conditions. The presence of the YMC influences the morphology (nucleation), photometry (photometric jumps) and spectroscopic output (nebular emission) of the entire galaxy. The de-lensed LyC and UV (1600A) magnitudes of the YMC are ~30.6 and ~26.9, whilst the galaxy has m1600~24.8. A relatively large rest-frame equivalent width of EWrest(Hb+[OIII]4959-5007)~450A emerges from the galaxy with the YMC contributing to ~30%. If O-type stars are mainly forged in star clusters, then such engines were the key ionizing agents during reionization and the increasing occurrence of high EW lines (Hb+[OIII]) observed at z>6.5 might be an indirect signature of a high Γat reionization.Future facilities (like VLT/MAVIS or ELT), will probe bound clusters on moderately magnified (μ<5-10) galaxies across cosmic epochs up to reionization[ABRIDGED]

preprint2021arXiv

Second generation star formation in globular clusters of different masses

By means of three-dimensional hydrodynamical simulations, we investigate the formation of second generation (SG) stars in young globular clusters of different masses. We consider clusters with a first generation of asymptotic giant branch (AGB) stars with mass 10^5 and 10^6 Msun moving at constant velocity through a uniform gas with density 10^(-24) and 10^(-23) g cm^(-3). Our setup is designed to reproduce the encounter of a young cluster with a reservoir of dense gas, e. g. during its orbital motion in the host galaxy. In the low-density models, as a result of the cooling AGB ejecta which collect in the centre, weakly perturbed by the external ram pressure, a compact central He-rich SG stellar component is formed on a timescale which decreases with increasing initial cluster mass. Our high-density models are subject to stronger ram pressure, which prevents the accumulation of the most He-rich AGB ejecta in the cluster centre. As a result, the SG is more extended and less He-enhanced than in the low-density models. By combining our results with previous simulations, we are able to study relevant, cluster-related scaling relations across a dynamical range of two orders of magnitude in mass (from 10^5 Msun to 10^7 Msun). In agreement with current observationally-based estimates, we find positive correlations between the SG-to-total number ratio and maximum He enhancement in SG stars as a function of the initial cluster mass.

preprint2020arXiv

Candidate Population III stellar complex at z=6.629 in the MUSE Deep Lensed Field

We discovered a strongly lensed (μ>40) Lya emission at z=6.629 (S/N~18) in the MUSE Deep Lensed Field (MDLF) targeting the Hubble Frontier Field galaxy cluster MACS~J0416. Dedicated lensing simulations imply that the Lya emitting region necessarily crosses the caustic. The arc-like shape of the Lya extends 3 arcsec on the observed plane and is the result of two merged multiple images, each one with a de-lensed Lya luminosity L<~2.8 x 10^(40) erg/s arising from a confined region (< 150 pc effective radius). A spatially unresolved HST counterpart is barely detected at S/N~2 after stacking the near-infrared bands, corresponding to an observed(intrinsic) magnitude m_(1500)>~30.8(>~35.0). The inferred rest-frame Lya equivalent width is EWo > 1120 A if the IGM transmission is T(IGM)<0.5. The low luminosities and the extremely large Lya EWo match the case of a Population~III star complex made of several dozens stars (~ 10^4 Msun) which irradiate a HII region crossing the caustic. While the Lya and stellar continuum are among the faintest ever observed at this redshift, the continuum and the Lya emissions could be affected by differential magnification, possibly biasing the EWo estimate. The aforementioned tentative HST detection tend to favor a large EWo, making such a faint Pop~III candidate a key target for the James Webb Space Telescope and Extremely Large Telescopes.

preprint2020arXiv

Discovery of molecular gas fueling galaxy growth in a protocluster at z=1.7

Based on ALMA Band 3 observations of the CO(2-1) line transition, we report the discovery of three new gas-rich (M_H2 ~ 1.5-4.8 x 10^10 M_sun, SFRs in the range ~5-100 M_sun/yr) galaxies in an overdense region at z=1.7, that already contains eight spectroscopically confirmed members. This leads to a total of 11 confirmed overdensity members, within a projected distance of ~ 1.15 Mpc and in a redshift range of Dz = 0.012. Under simple assumptions, we estimate that the system has a total mass of >= 3-6 x 10^13 M_sun, and show that it will likely evolve into a >~ 10^14 M_sun cluster at z = 0. The overdensity includes a powerful Compton-thick Fanaroff-Riley type II (FRII) radio-galaxy, around which we discovered a large molecular gas reservoir (M_H2 ~ 2 x 10^11 M_sun). We fitted the FRII resolved CO emission with a 2-D Gaussian model with major (minor) axis of ~ 27 (~ 17) kpc, that is a factor of ~3 larger than the optical rest-frame emission. Under the assumption of a simple edge-on disk morphology, we find that the galaxy interstellar medium produces a column density towards the nucleus of ~ 5.5 x 10^23 cm^-2. Such a dense ISM may then contribute significantly to the total nuclear obscuration measured in the X-rays (N_(H,X) ~ 1.5 x 10^24 cm^-2) in addition to a small, pc-scale absorber around the central engine. The velocity map of this source unveils a rotational motion of the gas that is perpendicular to the radio-jets. The FRII is located at the center of the projected spatial distribution of the structure members, and its velocity offset from the peak of the redshift distribution is well within the structure&#39;s velocity dispersion. All this, coupled with the large amount of gas around the FRII, its stellar mass of ~ 3 x 10^11 M_sun, SFR of ~ 200-600 M_sun/yr, and powerful radio-to-X-ray emission, suggests that this source is the likely progenitor of the future brightest cluster galaxy.

preprint2020arXiv

Dust and gas content of high-redshift galaxies hosting obscured AGN in the CDF-S

Obscured AGN represent a significant fraction of the entire AGN population, especially at high redshift (~70% at z=3--5). They are often characterized by the presence of large gas and dust reservoirs that are thought to sustain and possibly obscure vigorous star formation processes that make these objects shine at far-IR and sub-mm wavelengths. We exploit ALMA Cycle 4 observations of the continuum (~2.1mm) and high-J CO emission of a sample of six X-ray selected SMGs hosting an obscured AGN at z_spec>2.5 in the 7 Ms CDF-S. We measured the masses and sizes of the dust and molecular gas and we derived the gas density and column density on the basis of a uniform sphere geometry. Finally, we compared the measured column densities with those derived from the Chandra X-ray spectra. We detected both the continuum and line emission for three sources for which we measured both the flux density and size. For the undetected sources, we derived an upper limit on the flux density. We found that the detected galaxies are rich in gas and dust (molecular gas mass in the range <0.5 - 2.7 x 10^10 M_sun for α_CO=0.8 and up to ~2 x 10^11~M_sun for α_CO=6.5, and dust mass <0.9 - 4.9 x 10^8 M_sun) and compact (gas major axis 2.1-3.0 kpc, dust major axis 1.4-2.7 kpc). The column densities associated with the ISM are on the order of 10^(23-24) cm-2, which is comparable with those derived from the X-ray spectra. For the detected sources we also derived dynamical masses in the range 0.8 - 3.7 x 10^10 M_sun. We conclude that the ISM of high redshift galaxies can substantially contribute to nuclear obscuration up to the Compton-thick (>10^24 cm-2) regime. In addition, we found that all the detected sources show a velocity gradient reminding one rotating system, even though two of them show peculiar features in their morphology that can be associated with a chaotic, possibly merging, structure.

preprint2020arXiv

Galactic Archaeology with asteroseismic ages. II. Confirmation of a delayed gas infall using Bayesian analysis based on MCMC methods

With the wealth of information from large surveys and observational campaigns in the contemporary era, it is critical to properly exploit the data to constrain the parameters of Galactic chemical evolution models and quantify the associated uncertainties. We aim at constraining the two-infall chemical evolution models for the solar annulus using the measured chemical abundance ratios and seismically inferred age of stars in the APOKASC sample. In the revised two-infall chemical evolution models by Spitoni et al. (2019), a significant delay of $\sim 4.3$ Gyr has been invoked between the two episodes of gas accretion. In this work, we wish to test its robustness and statistically confirm/quantify the delay. For the first time, a Bayesian framework based on Markov Chain Monte Carlo methods has been used for fitting the two-infall chemical evolution models to the data. In addition to fitting the data for stars in the APOKASC sample, our best fit models also reproduce other important observational constraints of the chemical evolution of the disk: i) present day stellar surface density; ii) present day supernova and star formation rates; iii) the metallicity distribution function; and iv) solar abundance values. We found a significant delay between the two gas accretion episodes for various models explored with different values for the star formation efficiencies. The values for the delay lie in the range $4.5-5.5$ Gyr. The results suggest that the APOKASC sample carries the signature of delayed gas-rich merger, with the dilution as main process determining the shape of low-$α$ stars in the abundance ratios space.

preprint2020arXiv

The deep Chandra survey in the SDSS J1030+0524 field

We present the X-ray source catalog for the 479 ks Chandra exposure of the SDSS J1030+0524 field, that is centered on a region that shows the best evidence to date of an overdensity around a z > 6 quasar, and also includes a galaxy overdensity around a Compton-thick Fanaroff-Riley type II radio galaxy at z = 1.7. Using wavdetect for initial source detection and ACIS Extract for source photometry and significance assessment, we create preliminary catalogs of sources that are detected in the full, soft, and hard bands, respectively. We produce X-ray simulations that mirror our Chandra observation to filter our preliminary catalogs and get a completeness level of > 91% and a reliability level of 95% in each band. The catalogs in the three bands are then matched into a final main catalog of 256 unique sources. Among them, 244, 193, and 208 are detected in the full, soft, and hard bands, respectively. This makes J1030 field the fifth deepest extragalactic X-ray survey to date. The field is part of the Multiwavelength Survey by Yale-Chile (MUSYC), and is also covered by optical imaging data from the Large Binocular Camera (LBC) at the Large Binocular Telescope, near-IR imaging data from the Canada France Hawaii Telescope WIRCam, and Spitzer IRAC. Thanks to its dense multi-wavelength coverage, J1030 represents a legacy field for the study of large-scale structures around distant accreting supermassive black holes. Using a likelihood ratio analysis, we associate multi-band counterparts for 252 (98.4%) of the 256 Chandra sources, with an estimated reliability of 95%. Finally, we compute the cumulative number of sources in each X-ray band, finding that they are in general agreement with the results from the Chandra Deep Fields.

preprint2020arXiv

The MUSE Deep Lensed Field on the Hubble Frontier Field MACS~J0416

Context: A census of faint and tiny star forming complexes at high redshift is key to improving our understanding of reionizing sources, galaxy growth and the formation of globular clusters. Aims: We present the MUSE Deep Lensed Field (MDLF) program. Methods: We describe Deep MUSE observations of 17.1 hours integration on a single pointing over the Hubble Frontier Field galaxy cluster MACS~J0416. Results: We confirm spectroscopic redshifts for all 136 multiple images of 48 source galaxies at 0.9<z<6.2. Within those galaxies, we securely identify 182 multiple images of 66 galaxy components that we use to constrain our lens model. We identify 116 clumps belonging to background high-z galaxies; the majority of them are multiple images and span magnitude, size and redshift intervals of [-18,-10], [~400-3] parsec and 1<z<6.6, respectively, with the most magnified ones probing possible single gravitationally bound star clusters. The depth of the MDLF combined with lensing magnification lead us to reach a detection limit for unresolved emission lines of a few 10$^{-20}$ erg/s/cm2, after correction for lensing magnification. Ultraviolet high-ionization metal lines (and HeII1640) are detected with S/N>10 for individual objects down to de-lensed magnitude 28-30 suggesting that they are common in such faint sources. Conclusions:Deep MUSE observations, in combination with existing HST imaging, allowed us to:(1) confirm redshifts for extremely faint high-z sources;(2) peer into their internal clumps (down to 100-200 pc scale);(3) in some cases break down such clumps into bound star clusters (<20 pc scale);(4) double the number of constraints for the lens model,reaching an unprecedented set of 182 bona-fide multiple images and confirming up to 213 galaxy cluster members. These results demonstrate the power that JWST and future ELTs will have when combined to study gravitational telescopes.[abridged]

preprint2015arXiv

Are z>2 Herschel galaxies proto-spheroids?

We present a backward approach for the interpretation of the evolution of the near-infrared and the far-infrared luminosity functions across the redshift range 0<z<3. In our method, late-type galaxies are treated by means of a parametric phenomenological method based on PEP/HerMES data up to z~4, whereas spheroids are described by means of a physically motivated backward model. The spectral evolution of spheroids is modelled by means of a single-mass model, associated to a present-day elliptical with K-band luminosity comparable to the break of the local early-type luminosity function. The formation of proto-spheroids is assumed to occurr across the redshift range 1< z < 5. The key parameter is represented by the redshift z_0.5 at which half proto-spheroids are already formed. A statistical study indicates for this parameter values between z_0.5=1.5 and z_0.5=3. We assume as fiducial value z_0.5~2, and show that this assumption allows us to describe accourately the redshift distributions and the source counts. By assuming z_0.5 ~ 2 at the far-IR flux limit of the PEP-COSMOS survey, the PEP-selected sources observed at z>2 can be explained as progenitors of local spheroids caught during their formation. We also test the effects of mass downsizing by dividing the spheroids into three populations of different present-day stellar masses. The results obtained in this case confirm the validity of our approach, i.e. that the bulk of proto-spheroids can be modelled by means of a single model which describes the evolution of galaxies at the break of the present-day early type K-band LF.

preprint2015arXiv

Feedback from massive stars and gas expulsion from proto-Globular Clusters

Globular clusters are considerably more complex structures than previously thought, harbouring at least two stellar generations which present clearly distinct chemical abundances. Scenarios explaining the abundance patterns in globular clusters mostly assume that originally the clusters had to be much more massive than today, and that the second generation of stars originates from the gas shed by stars of the first generation (FG). The lack of metallicity spread in most globular clusters further requires that the supernova-enriched gas ejected by the FG is completely lost within ~30 Myr, a hypothesis never tested by means of three-dimensional hydrodynamic simulations. In this paper, we use 3D hydrodynamic simulations including stellar feedback from winds and supernovae, radiative cooling and self-gravity to study whether a realistic distribution of OB associations in a massive proto-GC of initial mass M_tot ~ 10^7 M_sun is sufficient to expel its entire gas content. Our numerical experiment shows that the coherence of different associations plays a fundamental role: as the bubbles interact, distort and merge, they carve narrow tunnels which reach deeper and deeper towards the innermost cluster regions, and through which the gas is able to escape. Our results indicate that after 3 Myr, the feedback from stellar winds is responsible for the removal of ~40% of the pristine gas, and that after 14 Myr, ~ 99% of the initial gas mass has been removed.

preprint2015arXiv

Star Formation in Herschel&#39;s Monsters versus Semi-Analytic Models

We present a direct comparison between the observed star formation rate functions (SFRF) and the state-of-the-art predictions of semi-analytic models (SAM) of galaxy formation and evolution. We use the PACS Evolutionary Probe Survey (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) data-sets in the COSMOS and GOODS-South fields, combined with broad-band photometry from UV to sub-mm, to obtain total (IR+UV) instantaneous star formation rates (SFRs) for individual Herschel galaxies up to z~4, subtracted of possible active galactic nucleus (AGN) contamination. The comparison with model predictions shows that SAMs broadly reproduce the observed SFRFs up to z~2, when the observational errors on the SFR are taken into account. However, all the models seem to under-predict the bright-end of the SFRF at z>2. The cause of this underprediction could lie in an improper modelling of several model ingredients, like too strong (AGN or stellar) feedback in the brighter objects or too low fall-back of gas, caused by weak feedback and outflows at earlier epochs.

preprint2014arXiv

ALMA reveals a warm and compact starburst around a heavily obscured supermassive black hole at z=4.75

We report ALMA Cycle 0 observations at 1.3mm of LESS J033229.4-275619 (XID403), an Ultraluminous Infrared Galaxy at $z=4.75$ in the Chandra Deep Field South hosting a Compton-thick QSO. The source is not resolved in our data at a resolution of $\sim$0.75 arcsec, placing an upper-limit of 2.5 kpc to the half-light radius of the continuum emission from heated-dust. After deconvolving for the beam size, however, we found a $\sim3σ$ indication of an intrinsic source size of $0.27\pm0.08$ arcsec (Gaussian FWHM), which would correspond to $r_{half}\sim0.9\pm0.3$ kpc. We build the far-IR SED of XID403 by combining datapoints from both ALMA and Herschel and fit it with a modified blackbody spectrum. For the first time, we measure the dust temperature $T_d=58.5\pm5.3$ K in this system, which is comparable to what has been observed in other high-z submillimeter galaxies. The measured star formation rate is SFR=$1020\pm150$ $M_{\odot}$ yr$^{-1}$, in agreement with previous estimates at lower S/N. Based on the measured SFR and source size, we constrain the SFR surface density to be $Σ_{SFR}>26\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ ($\sim200\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ for $r_{half}\sim0.9$ kpc). The compactness of this starburst is comparable to what has been observed in other local and high-z starburst galaxies. If the gas mass measured from previous [CII] and CO(2-1) observations at low resolution is confined within the same dust region, assuming $r_{half}\sim0.9\pm0.3$ kpc, this would produce a column density of $N_H\sim0.3-1.1\times10^{24}$cm$^{-2}$ towards the central SMBH, similar to the column density of $\approx1.4\times10^{24}$cm$^{-2}$ measured from the X-rays. Then, in principle, if both gas and dust were confined on sub-kpc scales, this would be sufficient to produce the observed X-ray column density without any need of a pc-scale absorber [abridged].

preprint2014arXiv

Chemical evolution models: GRB host identification and cosmic dust predictions

The nature of some GRB host galaxies has been investigated by means of chemical evolution models of galaxies of different morphological type following the evolution of the abundances of H, He, C, N, O, $α$-elements, Ni, Fe, Zn, and including also the evolution of dust. By comparing predictions with abundance data, we were able to constrain nature and age of GRB hosts. We also computed a theoretical cosmic dust rate, including stellar dust production, accretion and destruction, under the hypotheses of pure luminosity evolution and strong number density evolution of galaxies. We suggest that one of the three GRB hosts is a massive proto-spheroid catched during its formation, while for the other two the situation is more uncertain, although one could perhaps be a spheroid and the other a spiral galaxy. We estimated the chemical ages of the host galaxies which vary from 15 to 320 Myr. Concerning the cosmic effective dust production rate in an unitary volume of the Universe, our results show that in the case of pure luminosity evolution there is a first peak between redshift $z=8$ and $9$ and another at $z\sim 5$, whereas in the case of strong number density evolution it increases slightly from $z=10$ to $z\sim 2$ and then it decreases down to $z=0$. Finally, we found tha the total cosmic dust mass density at the present time is: $Ω_{dust} \sim 3.5\cdot 10^{-5}$in the case of pure luminosity evolution and $Ω_{dust} \sim 7\cdot 10^{-5}$ in the case of number density evolution.

preprint2014arXiv

Chemodynamics of a Simulated Disc Galaxy: Initial Mass Functions and SNIa Progenitors

We trace the formation and advection of several elements within a cosmological adaptive mesh refinement simulation of an L* galaxy. We use nine realisations of the same initial conditions with different stellar Initial Mass Functions (IMFs), mass limits for type-II and type-Ia supernovae (SNII, SNIa) and stellar lifetimes to constrain these sub-grid phenomena. Our code includes self-gravity, hydrodynamics, star formation, radiative cooling and feedback from multiple sources within a cosmological framework. Under our assumptions of nucleosynthesis we find that SNII with progenitor masses of up to 100 Msun are required to match low metallicity gas oxygen abundances. Tardy SNIa are necessary to reproduce the classical chemical evolution knee in [O/Fe]-[Fe/H]: more prompt SNIa delayed time distributions do not reproduce this feature. Within our framework of hydrodynamical mixing of metals and galaxy mergers we find that chemical evolution is sensitive to the shape of the IMF and that there exists a degeneracy with the mass range of SNII. We look at the abundance plane and present the properties of different regions of the plot, noting the distinct chemical properties of satellites and a series of nested discs that have greater velocity dispersions, are more alpha-rich and metal poor with age.

preprint2014arXiv

VLT/X-shooter spectroscopy of the GRB 120327A afterglow

We present a study of the environment of the Swift long gamma-ray burst GRB 120327A at z ~2.8 through optical spectroscopy of its afterglow. We analyzed medium-resolution, multi-epoch spectroscopic observations (~7000 - 12000, corresponding to ~ 15 - 23 km/s, S/N = 15- 30 and wavelength range 3000-25000AA) of the optical afterglow of GRB 120327A, taken with X-shooter at the VLT 2.13 and 27.65 hr after the GRB trigger. The first epoch spectrum shows that the ISM in the GRB host galaxy at z = 2.8145 is extremely rich in absorption features, with three components contributing to the line profiles. The hydrogen column density associated with GRB 120327A has log NH / cm^(-2) = 22.01 +/- 0.09, and the metallicity of the host galaxy is in the range [X/H] = -1.3 to -1.1. In addition to the ground state lines, we detect absorption features associated with excited states of CII, OI, SiII, FeII, and NiII, which we used to derive information on the distance between the host absorbing gas and the site of the GRB explosion. The variability of the FeI\lambda2396 excited line between the two epochs proves that these features are excited by the GRB UV flux. Moreover, the distance of component I is found to be dI=200+100-60 pc, while component II is located closer to the GRB, at dII=100+40-30 pc. These values are among the lowest found in GRBs. Component III does not show excited transitions, so it should be located farther away from the GRB. The presence of H2 molecules is firmly established, with a molecular fraction in the range f=4 X 10^(-7) - 10^(-4). This particularly low value can be attributed to the small dust content. This represents the third positive detection of molecules in a GRB environment.

preprint2013arXiv

A fast and accurate method to compute the mass return from multiple stellar populations

The mass returned to the ambient medium by aging stellar populations over cosmological times sums up to a significant fraction (20% - 30% or more) of their initial mass. This continuous mass injection plays a fundamental role in phenomena such as galaxy formation and evolution, fueling of supermassive black holes in galaxies and the consequent (negative and positive) feedback phenomena, and the origin of multiple stellar populations in globular clusters. In numerical simulations the calculation of the mass return can be time consuming, since it requires at each time step the evaluation of a convolution integral over the whole star formation history, so the computational time increases quadratically with the number of time-steps. The situation can be especially critical in hydrodynamical simulations, where different grid points are characterized by different star formation histories, and the gas cooling and heating times are shorter by orders of magnitude than the characteristic stellar lifetimes. In this paper we present a fast and accurate method to compute the mass return from stellar populations undergoing arbitrarily complicated star formation histories. At each time-step the mass return is calculated from its value at the previous time, and the star formation rate over the last time-step only. Therefore in the new scheme there is no need to store the whole star formation history, and the computational time increases linearly with the number of time-steps.

preprint2013arXiv

The [alpha/Fe] ratios of very metal-poor stars within the IGIMF theory

The aim of this paper is to quantify the amplitude of the predicted plateau in [alpha/Fe] ratios associated with the most metal-poor stars of a galaxy. We assume that the initial mass function in galaxies is steeper if the star formation rate (SFR) is low -- as per the integrated galactic initial mass function (IGIMF) theory. A variant of the theory, in which the IGIMF depends upon the metallicity of the parent galaxy, is also considered. The IGIMF theory predicts low [alpha/Fe] plateaus in dwarf galaxies, characterised by small SFRs. The [alpha/Fe] plateau is up to 0.7dex lower than the corresponding plateau of the Milky Way. For a universal IMF one should expect instead that the [alpha/Fe] plateau is the same for all the galaxies, irrespective of their masses or SFRs. Assuming a strong dependence of the IMF on the metallicity of the parent galaxy, dwarf galaxies can show values of the [alpha/Fe] plateau similar to those of the Milky Way, and almost independent on the SFR. The [Mg/Fe] ratios of the most metal-poor stars in dwarf galaxies satellites of the Milky Way can be reproduced either if we consider metallicity-dependent IMFs or if the early SFRs of these galaxies were larger than we presently think. Present and future observations of dwarf galaxies can help disentangle between these different IGIMF formulations.

preprint2013arXiv

The dust content of QSO hosts at high redshift

Infrared observations of high-z quasar (QSO) hosts indicate the presence of large masses of dust in the early universe. When combined with other observables, such as neutral gas masses and star formation rates, the dust content of z~6 QSO hosts may help constraining their star formation history. We have collected a database of 58 sources from the literature discovered by various surveys and observed in the FIR. We have interpreted the available data by means of chemical evolution models for forming proto-spheroids, investigating the role of the major parameters regulating star formation and dust production. For a few systems, given the derived small dynamical masses, the observed dust content can be explained only assuming a top-heavy initial mass function, an enhanced star formation efficiency and an increased rate of dust accretion. However, the possibility that, for some systems, the dynamical mass has been underestimated cannot be excluded. If this were the case, the dust mass can be accounted for by standard model assumptions. We provide predictions regarding the abundance of the descendants of QSO hosts; albeit rare, such systems should be present and detectable by future deep surveys such as Euclid already at z>4.

preprint2013arXiv

The two regimes of the cosmic sSFR evolution are due to spheroids and discs

This paper aims at explaining the two phases in the observed specific star formation rate (sSFR), namely the high (>3/Gyr) values at z>2 and the smooth decrease since z=2. In order to do this, we compare to observations the specific star formation rate evolution predicted by well calibrated models of chemical evolution for elliptical and spiral galaxies, using the additional constraints on the mean stellar ages of these galaxies (at a given mass). We can conclude that the two phases of the sSFR evolution across cosmic time are due to different populations of galaxies. At z>2 the contribution comes from spheroids: the progenitors of present-day massive ellipticals (which feature the highest sSFR) as well as halos and bulges in spirals (which contribute with average and lower-than-average sSFR). In each single galaxy the sSFR decreases rapidly and the star formation stops in <1 Gyr. However the combination of different generations of ellipticals in formation might result in an apparent lack of strong evolution of the sSFR (averaged over a population) at high redshift. The z<2 decrease is due to the slow evolution of the gas fraction in discs, modulated by the gas accretion history and regulated by the Schmidt law. The Milky Way makes no exception to this behaviour.

preprint2012arXiv

Comparison of star formation rates from Halpha and infrared luminosities as seen by Herschel

We empirically test the relation between the SFR(LIR) derived from the infrared luminosity, LIR, and the SFR(Ha) derived from the Ha emission line luminosity using simple conversion relations. We use a sample of 474 galaxies at z = 0.06 - 0.46 with both Ha detection (from 20k zCOSMOS survey) and new far-IR Herschel data (100 and 160 μm). We derive SFR(Ha) from the Ha extinction corrected emission line luminosity. We find a very clear trend between E(B - V) and LIR that allows to estimate extinction values for each galaxy even if the Ha emission line measurement is not reliable. We calculate the LIR by integrating from 8 up to 1000 μm the SED that is best fitting our data. We compare SFR(Ha) with the SFR(LIR). We find a very good agreement between the two SFR estimates, with a slope of m = 1.01 \pm 0.03 in the SFR(LIR) vs SFR(Ha) diagram, a normalization constant of a = -0.08 \pm 0.03 and a dispersion of sigma = 0.28 dex.We study the effect of some intrinsic properties of the galaxies in the SFR(LIR)-SFR(Ha) relation, such as the redshift, the mass, the SSFR or the metallicity. The metallicity is the parameter that affects most the SFR comparison. The mean ratio of the two SFR estimators log[SFR(LIR)/SFR(Ha)] varies by approx. 0.6 dex from metal-poor to metal-rich galaxies (8.1 < log(O/H) + 12 < 9.2). This effect is consistent with the prediction of a theoretical model for the dust evolution in spiral galaxies. Considering different morphological types, we find a very good agreement between the two SFR indicators for the Sa, Sb and Sc morphologically classified galaxies, both in slope and normalization. For the Sd, irregular sample (Sd/Irr), the formal best-fit slope becomes much steeper (m = 1.62 \pm 0.43), but it is still consistent with 1 at the 1.5 sigma level, because of the reduced statistics of this sub-sample.

preprint2012arXiv

Cosmic star formation rate: a theoretical approach

The cosmic star formation rate (CSFR), is an important clue to investigate the history of the assembly and evolution of galaxies. Here, we develop a method to study the CSFR from a purely theoretical point of view. Starting from detailed models of chemical evolution, we obtain the histories of star formation of galaxies of different morphological types. These histories are then used to determine the luminosity functions of the same galaxies by means of a spectro-photometric code. We obtain the CSFR under different hypothesis. First, we study the hypothesis of a pure luminosity evolution scenario, in which all galaxies are supposed to form at the same redshift and then evolve only in luminosity. Then we consider scenarios in which the number density or the slope of the LFs are assumed to vary with redshift. After comparison with available data we conclude that a pure luminosity evolution does not provide a good fit to the data, especially at very high redshift, although many uncertainties are still present in the data. On the other hand, a variation in the number density of ellipticals and spirals as a function of redshift can provide a better fit to the observed CSFR. We also explore cases of variable slope of the LFs with redshift and variations of number density and slope at the same time. We cannot find any of those cases which can improve the fit to the data respect to the solely number density variation. Finally, we compute the evolution of the average cosmic metallicity in galaxies with redshift.

preprint2012arXiv

Metallicity Gradients in Disks: Do Galaxies Form Inside-Out?

We examine radial and vertical metallicity gradients using a suite of disk galaxy simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient and reconcile differences between extant models and observations within the `inside-out&#39; disk growth paradigm. A sample of 25 disks is used, consisting of 19 from our RaDES (Ramses Disk Environment Study) sample, realised with the adaptive mesh refinement code RAMSES. Four disks are selected from the MUGS (McMaster Unbiased Galaxy Simulations) sample, generated with the smoothed particle hydrodynamics (SPH) code GASOLINE, alongside disks from Rahimi et al. (GCD+) and Kobayashi & Nakasato (GRAPE-SPH). Two chemical evolution models of inside-out disk growth were employed to contrast the temporal evolution of their radial gradients with those of the simulations. We find that systematic differences exist between the predicted evolution of radial abundance gradients in the RaDES and chemical evolution models, compared with the MUGS sample; specifically, the MUGS simulations are systematically steeper at high-redshift, and present much more rapid evolution in their gradients. We find that the majority of the models predict radial gradients today which are consistent with those observed in late-type disks, but they evolve to this self-similarity in different fashions, despite each adhering to classical `inside-out&#39; growth. We find that radial dependence of the efficiency with which stars form as a function of time drives the differences seen in the gradients; systematic differences in the sub-grid physics between the various codes are responsible for setting these gradients. Recent, albeit limited, data at redshift z=1.5 are consistent with the steeper gradients seen in our SPH sample, suggesting a modest revision of the classical chemical evolution models may be required.

preprint2012arXiv

The Distribution of Metals in Cosmological Hydrodynamical Simulations of Dwarf Disk Galaxies

We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass -- halo mass relationship. Each realisation employed identical initial conditions and assembly histories, but differed from one another in their underlying baryonic physics prescriptions, including (a) the efficiency with which each supernova energy couples to the ISM, (b) the impact of feedback associated with massive star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling of Type II SNe remnants, (d) the treatment of metal diffusion, and (e) varying the IMF. Our analysis focusses on the resulting stellar metallicity distribution functions (MDFs) in each simulated (analogous) `solar neighbourhood&#39; and central `bulge&#39; region. We compare the simulated MDFs&#39; skewness, kurtosis, and dispersion (inter-quartile, inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empirical solar neighbourhood MDF and Local Group dwarfs. We find that the MDFs of the simulated discs are more negatively skewed, with higher kurtosis, than those observed locally. We can trace this difference to the simulations&#39; tight and correlated age-metallicity relations (compared with that of the Milky Way), suggesting that these relations within `dwarf&#39; discs might be steeper than in L* discs and/or the degree of stellar orbital re-distribution and migration inferred locally has not been captured in their entirety, at the resolution of our simulations. The important role of metal diffusion in ameliorating the over-production of extremely metal-poor stars is highlighted.

preprint2012arXiv

The Lyman alpha forest flux probability distribution at z>3

We present a measurement of the Lyman alpha flux probability distribution function (PDF) measured from a set of eight high resolution quasar spectra with emission redshifts at 3.3 < z < 3.8. We carefully study the effect of metal absorption lines on the shape of the PDF. Metals have a larger impact on the PDF measurements at lower redshift, where there are fewer Lyman alpha absorption lines. This may be explained by an increase in the number of metal lines which are blended with Lyman alpha absorption lines toward higher redshift, but may also be due to the presence of fewer metals in the intergalactic medium with increasing lookback time. We also provide a new measurement of the redshift evolution of the effective optical depth, tau_eff, at 2.8 < z < 3.6, and find no evidence for a deviation from a power law evolution in the log(tau_eff)-log(1+z) plane. The flux PDF measurements are furthermore of interest for studies of the thermal state of the intergalactic medium (IGM) at z ~ 3 . By comparing the PDF to state-of-the-art cosmological hydrodynamical simulations, we place constraints on the temperature of the IGM and compare our results with previous measurements of the PDF at lower redshift. At redshift z=3, our new PDF measurements are consistent with an isothermal temperature-density relation, T=T_0 Delta^{gamma-1}, with a temperature at the mean density of T_0 = 19250 +/- 4800 K and a slope gamma=0.90+/-0.21 (1 sigma uncertainties). In comparison, joint constraints with previous PDF measurements at z<3 favour an inverted (gamma<1) temperature-density relation with T_0=17900 +/- 3500 K and gamma=0.70 +/- 0.12, in broad agreement with previous analyses.

preprint2012arXiv

The stellar metallicity distribution of disc galaxies and bulges in cosmological simulations

By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born &#34;in situ&#34; relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.

preprint2011arXiv

The Chemical and Dynamical Evolution of Isolated Dwarf Galaxies

Using a suite of simulations (Governato et al. 2010) which successfully produce bulgeless (dwarf) disk galaxies, we provide an analysis of their associated cold interstellar media (ISM) and stellar chemical abundance patterns. A preliminary comparison with observations is undertaken, in order to assess whether the properties of the cold gas and chemistry of the stellar components are recovered successfully. To this end, we have extracted the radial and vertical gas density profiles, neutral hydrogen velocity dispersion, and the power spectrum of structure within the ISM. We complement this analysis of the cold gas with a brief examination of the simulations&#39; metallicity distribution functions and the distribution of alpha-elements-to-iron.

preprint2011arXiv

The Cold Gas Content of Bulgeless Dwarf Galaxies

We present an analysis of the neutral hydrogen (HI) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al. (2010), the high resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the HI distribution and kinematics of this simulated bulgeless disk with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g., velocity ellipsoid, turbulence), and the power spectrum of structure within the cold interstellar medium from the simulations. The highest resolution dwarf, when using a high density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the SMC. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within pc-scale molecular clouds.

preprint2010arXiv

Dynamical properties of AMAZE and LSD galaxies from gas kinematics and the Tully-Fisher relation at z~3

We present a SINFONI integral field kinematical study of 33 galaxies at z~3 from the AMAZE and LSD projects which are aimed at studying metallicity and dynamics of high-redshift galaxies. The number of galaxies analyzed in this paper constitutes a significant improvement compared to existing data in the literature and this is the first time that a dynamical analysis is obtained for a relatively large sample of galaxies at z~3. 11 galaxies show ordered rotational motions (~30% of the sample), in these cases we estimate dynamical masses by modeling the gas kinematics with rotating disks and exponential mass distributions. We find dynamical masses in the range 2 \times 10^9 M\odot - 2 \times 10^11 M\odot with a mean value of ~ 2 \times 10^10 M\odot. By comparing observed gas velocity dispersion with that expected from models, we find that most rotating objects are dynamically &#34;hot&#34;, with intrinsic velocity dispersions of the order of ~90 km s-1. The median value of the ratio between the maximum disk rotational velocity and the intrinsic velocity dispersion for the rotating objects is 1.6, much lower than observed in local galaxies value (~10) and slightly lower than the z~2 value (2 - 4). Finally we use the maximum rotational velocity from our modeling to build a baryonic Tully-Fisher relation at z~3. Our measurements indicate that z~3 galaxies have lower stellar masses (by a factor of ten on average) compared to local galaxies with the same dynamical mass. However, the large observed scatter suggests that the Tully-Fisher relation is not yet &#34;in place&#34; at these early cosmic ages, possibly due to the young age of galaxies. A smaller dispersion of the Tuly-Fisher relation is obtained by taking into account the velocity dispersion with the use of the S_0.5 indicator, suggesting that turbulent motions might have an important dynamical role.

preprint2010arXiv

Effects of the integrated galactic IMF on the chemical evolution of the solar neighbourhood

The initial mass function determines the fraction of stars of different intial mass born per stellar generation. In this paper, we test the effects of the integrated galactic initial mass function (IGIMF) on the chemical evolution of the solar neighbourhood. The IGIMF (Weidner & Kroupa 2005) is computed from the combination of the stellar intial mass function (IMF), i.e. the mass function of single star clusters, and the embedded cluster mass function, i.e. a power law with index beta. By taking into account also the fact that the maximum achievable stellar mass is a function of the total mass of the cluster, the IGIMF becomes a time-varying IMF which depends on the star formation rate. We applied this formalism to a chemical evolution model for the solar neighbourhood and compared the results obtained by assuming three possible values for beta with the results obtained by means of a standard, well-tested, constant IMF. In general, a lower absolute value of beta implies a flatter IGIMF, hence a larger number of massive stars and larger metal ejection rates. This translates into higher type Ia and II supernova rates, higher mass ejection rates from massive stars and a larger amount of gas available for star formation, coupled with lower present-day stellar mass densities. (abridged) We also discuss the importance of the present day stellar mass function (PDMF) in providing a way to disentangle among various assumptions for beta. Our results indicate that the model adopting the IGIMF computed with beta ~2 should be considered the best since it allows us to reproduce the observed PDMF and to account for most of the chemical evolution constraints considered in this work.

preprint2010arXiv

Spectral catalogue of bright gamma-ray bursts detected with the BeppoSAX/GRBM

The emission process responsible for the so-called &#34;prompt&#34; emission of gamma-ray bursts is still unknown. A number of empirical models fitting the typical spectrum still lack a satisfactory interpretation. A few GRB spectral catalogues derived from past and present experiments are known in the literature and allow to tackle the issue of spectral properties of gamma-ray bursts on a statistical ground. We extracted and studied the time-integrated photon spectra of the 200 brightest GRBs observed with the Gamma-Ray Burst Monitor which flew aboard the BeppoSAX mission (1996-2002) to provide an independent statistical characterisation of GRB spectra. The spectra were fit with three models: a simple power-law, a cut-off power law or a Band function. The typical photon spectrum of a bright GRB consists of a low-energy index around 1.0 and a peak energy of the nuFnu spectrum E_p~240 keV in agreement with previous results on a sample of bright CGRO/BATSE bursts. Spectra of ~35% of GRBs can be fit with a power-law with a photon index around 2, indicative of peak energies either close to or outside the GRBM energy boundaries. We confirm the correlation between E_p and fluence, with a logarithmic dispersion of 0.13 around the power-law with index 0.21+-0.06. The low-energy and peak energy distributions are not yet explained in the current literature. The capability of measuring time-resolved spectra over a broadband energy range, ensuring precise measurements of parameters such as E_p, will be crucial for future experiments (abridged).

preprint2010arXiv

The chemical evolution of elliptical galaxies with stellar and QSO dust production

The presence of dust strongly affects the way we see galaxies and also the chemical abundances we measure in gas. It is therefore important to study he chemical evolution of galaxies by taking into account dust evolution. We aim at performing a detailed study of abundance ratios of high redshift objects and their dust properties. We focus on Lyman-Break galaxies (LBGs) and Quasar (QSO) hosts and likely progenitors of low- and high-mass present-day elliptical galaxies, respectively. We have adopted a chemical evolution model for elliptical galaxies taking account the dust production from low and intermediate mass stars, supernovae Ia, supernovae II, QSOs and both dust destruction and accretion processes. By means of such a model we have followed the chemical evolution of ellipticals of different baryonic masses. Our model complies with chemical downsizing. We made predictions for the abundance ratios versus metallicity trends for models of differing masses that can be used to constrain the star formation rate, initial mass function and dust mass in observed galaxies. We predict the existence of a high redshift dust mass-stellar mass relationship. We have found a good agreement with the properties of LBGs if we assume that they formed at redshift z=2-4. In particular, a non-negligible amount of dust is needed to explain the observed abundance pattern. We studied the QSO SDSS J114816, one of the most distant QSO ever observed (z=6.4), and we have been able to reproduce the amount of dust measured in this object. The dust is clearly due to the production from supernovae and the most massive AGB stars as well as from the grain growth in the interstellar medium. The QSO dust is likely to dominate only in the very central regions of the galaxies and during the early development of the galactic wind.

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 Origin of the Mass-Metallicity relation: an analytical approach

The existence of a mass-metallicity (MZ) relation in star forming galaxies at all redshift has been recently established. We aim at studying some possible physical mechanisms contributing to the MZ relation by adopting analytical solutions of chemical evolution models including infall and outflow. We explore the hypotheses of a variable galactic wind rate, infall rate and yield per stellar generation (i.e. a variation in the IMF), as possible causes for the MZ relation. By means of analytical models we compute the expected O abundance for galaxies of a given total baryonic mass and gas mass.The stellar mass is derived observationally and the gas mass is derived by inverting the Kennicutt law of star formation, once the star formation rate is known. Then we test how the parameters describing the outflow, infall and IMF should vary to reproduce the MZ relation, and we exclude the cases where such a variation leads to unrealistic situations. We find that a galactic wind rate increasing with decreasing galactic mass or a variable IMF are both viable solutions for the MZ relation. A variable infall rate instead is not acceptable. It is difficult to disentangle among the outflow and IMF solutions only by considering the MZ relation, and other observational constraints should be taken into account to select a specific solution. For example, a variable efficiency of star formation increasing with galactic mass can also reproduce the MZ relation and explain the downsizing in star formation suggested for ellipticals. The best solution could be a variable efficiency of star formation coupled with galactic winds, which are indeed observed in low mass galaxies.

preprint2009arXiv

The evolution of the mass-metallicity relation in galaxies of different morphological types

By means of chemical evolution models for ellipticals, spirals and irregular galaxies, we aim at investigating the physical meaning and the redshift evolution of the mass-metallicity relation as well as how this relation is connected with galaxy morphology. {abridged} We assume that galaxy morphologies do not change with cosmic time. We present a method to account for a spread in the epochs of galaxy formation and to refine the galactic mass grid. (abridged) We compare our predictions to observational results obtained for galaxies between redshifts 0.07 and 3.5. We reproduce the mass-metallicity (MZ) relation mainly by means of an increasing efficiency of star formation with mass in galaxies of all morphological types, without any need to invokegalactic outflows favoring the loss of metals in the less massive galaxies. Our predictions can help constraining the slope and the zero point of the observed local MZ relation, both affected by uncertainties related to the use of different metallicity calibrations. We show how, by considering the MZ, the O/H vs star formation rate (SFR), and the SFR vs galactic mass diagrams at various redshifts, it is possible to constrain the morphology of the galaxies producing these relations. Our results indicate that the galaxies observed at z=3.5 should be mainly proto-ellipticals, whereas at z=2.2 the observed galaxies consist of a morphological mix of proto-spirals and proto-ellipticals. At lower redshifts, the observed MZ relation is well reproduced by considering both spirals and irregulars. (abridged)

preprint2009arXiv

The Mass-Metallicity relation in galaxies of different morphological types

By means of chemical evolution models of different morphological types, we study the mass-metallicity (MZ) relation and its evolution with redshift. Our aim is to understand the role of galaxies of different morphological types in the MZ relation at various redshift. One major result is that at high redshift, the majority of the galaxies falling on the MZ plot are apparently proto-ellipticals. Finally, we show some preliminary results of a study of the MZ relation in a framework of hierarchical galaxy formation.

preprint2009arXiv

Theoretical cosmic Type Ia supernova rates

The aim of this work is the computation of the cosmic Type Ia supernova rates at very high redshifts (z>2). We adopt various progenitor models in order to predict the number of explosions in different scenarios for galaxy formation and to check whether it is possible to select the best delay time distribution model, on the basis of the available observations of Type Ia supernovae. We also computed the Type Ia supernova rate in typical elliptical galaxies of different initial luminous masses and the total amount of iron produced by Type Ia supernovae in each case. It emerges that: it is not easy to select the best delay time distribution scenario from the observational data and this is because the cosmic star formation rate dominates over the distribution function of the delay times; the monolithic collapse scenario predicts an increasing trend of the SN Ia rate at high redshifts whereas the predicted rate in the hierarchical scheme drops dramatically at high redshift; for the elliptical galaxies we note that the predicted maximum of the Type Ia supernova rate depends on the initial galactic mass. The maximum occurs earlier (at about 0.3 Gyr) in the most massive ellipticals, as a consequence of downsizing in star formation. We find that different delay time distributions predict different relations between the Type Ia supernova rate per unit mass at the present time and the color of the parent galaxies and that bluer ellipticals present higher supernova Type Ia rates at the present time.

preprint2008arXiv

AMAZE. I. The evolution of the mass-metallicity relation at z>3

We present initial results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z>3 by means of deep near-IR spectroscopy. Gas metallicities are measured, for an initial sample of nine star forming galaxies at z~3.5, by means of optical nebular lines redshifted into the near-IR. Stellar masses are accurately determined by using Spitzer-IRAC data, which sample the rest-frame near-IR stellar light in these distant galaxies. When compared with previous surveys, the mass-metallicity relation inferred at z~3.5 shows an evolution much stronger than observed at lower redshifts. The evolution is prominent even in massive galaxies, indicating that z~3 is an epoch of major action in terms of star formation and metal enrichment also for massive systems. There are also indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.5 is difficult to reconcile with the predictions of some hierarchical evolutionary models. Such discrepancies suggest that at z>3 galaxies are assembled mostly with relatively un-evolved sub-units, i.e. small galaxies with low star formation efficiency. The bulk of the star formation and metallicity evolution probably occurs once small galaxies are already assembled into bigger systems.