Source author record

Alessandro Bressan

Alessandro Bressan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

52works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

52 published item(s)

preprint2023arXiv

The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come.

preprint2022arXiv

Dust evolution with MUPPI in cosmological volumes

We study the evolution of dust in a cosmological volume using a hydrodynamical simulation in which the dust production is coupled with the MUPPI (MUlti Phase Particle Integrator) sub-resolution model of star formation and feedback. As for the latter, we keep as reference the model setup calibrated previously to match the general properties of Milky Way like galaxies in zoom-in simulations. However, we suggest that an increase of the star formation efficiency with the local dust to gas ratio would better reproduce the observed evolution of the cosmic star formation density. Moreover, the paucity of quenched galaxies at low redshift demands a stronger role of AGN feedback. We tune the parameters ruling direct dust production from evolved stars and accretion in the inter stellar medium to get scaling relations involving dust, stellar mass and metallicity in good agreement with observations. In low mass galaxies the accretion process is inefficient. As a consequence, they remain poorer in silicate and small grains than higher mass ones. We reproduce reasonably well the few available data on the radial distribution of dust outside the galactic region, supporting the assumption that the dust and gas dynamics are well coupled at galactic scales.

preprint2022arXiv

Formation of black holes in the pair-instability mass gap: evolution of a post-collision star

The detection of GW190521 by the LIGO-Virgo collaboration has revealed the existence of black holes (BHs) in the pair-instability (PI) mass gap. Here, we investigate the formation of BHs in the PI mass gap via star -- star collisions in young stellar clusters. To avoid PI, the stellar-collision product must have a relatively small core and a massive envelope. We generate our initial conditions from the outputs of a hydro-dynamical simulation of the collision between a core helium burning star ($\sim 58$ M$_\odot$) and a main-sequence star ($\sim 42$ M$_\odot$). The hydro-dynamical simulation allows us to take into account the mass lost during the collision ($\sim 12$ M$_\odot$) and to build the chemical composition profile of the post-collision star. We then evolve the collision product with the stellar evolution codes PARSEC and MESA. We find that the post-collision star evolves through all the stellar burning phases until core collapse, avoiding PI. At the onset of core collapse, the post-collision product is a blue super-giant star. We estimate a total mass loss of about 1 M$_\odot$ during the post-collision evolution, due to stellar winds and shocks induced by neutrino emission in a failed supernova. The final BH mass is $\approx{87}$ M$_\odot$. Therefore, we confirm that the collision scenario is a suitable formation channel to populate the PI mass gap.

preprint2022arXiv

Simulating the Legacy Survey of Space and Time stellar content with TRILEGAL

We describe a large simulation of the stars to be observed by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The simulation is based on the TRILEGAL code, which resorts to large databases of stellar evolutionary tracks, synthetic spectra, and pulsation models, added to simple prescriptions for the stellar density and star formation histories of the main structures of the Galaxy, to generate mock stellar samples through a population synthesis approach. The main bodies of the Magellanic Clouds are also included. A complete simulation is provided for single stars, down to the $r=27.5$ mag depth of the co-added wide-fast-deep survey images. A second simulation is provided for a fraction of the binaries, including the interacting ones, as derived with the BinaPSE module of TRILEGAL. We illustrate the main properties and numbers derived from these simulations, including: comparisons with real star counts; the expected numbers of Cepheids, long-period variables and eclipsing binaries; the crowding limits as a function of seeing and filter; the star-to-galaxy ratios, etc. Complete catalogs are accessible through the NOIRLab Astro Data Lab, while the stellar density maps are incorporated in the LSST metrics analysis framework (MAF).

preprint2021arXiv

A Fresh Look at AGB stars in Galactic Open Clusters with Gaia: Impact on Stellar Models and the Initial-Final Mass Relation

Benefiting from the GAIA second and early third releases of photometric and astrometric data we examine the population of asymptotic giant branch (AGB) stars that appear in the fields of intermediate-age and young open star clusters. We identify 49 AGB star candidates, brighter than the tip of the red giant branch, with a good-to-high cluster membership probability. Among them we find 19 TP-AGB stars with known spectral type: 4 M stars, 3 MS/S stars and 12 C stars. By combining observations, stellar models, and radiative transfer calculations that include the effect of circumstellar dust, we characterize each star in terms of initial mass, luminosity, mass-loss rate, core mass, period and mode of pulsation. The information collected helps us shed light on the TP-AGB evolution at solar-like metallicity, placing constraints on the third dredge-up process, the initial masses of carbon stars, stellar winds, and the initial-final mass relation (IFMR). In particular, we find that two bright carbon stars, MSB 75 and BM IV 90, members of the clusters NGC 7789 and NGC 2660 (with similar ages of $\simeq 1.2-1.6$ Gyr and initial masses $ 2.1 \ge M_{\rm i}/M_{\odot} \ge 1.9$), have unusually high core masses, $M_{\rm c} \approx 0.67-0.7\,M_{\odot}$. These results support the findings of a recent work (Marigo et al. 2020) that identified a kink in the IFMR, which interrupts its monotonic trend just at the same initial masses. Finally, we investigate two competing scenarios to explain the $M_{\rm c}$ data: the role of stellar winds in single-star evolution, and binary interactions through the blue-straggler channel.

preprint2021arXiv

The Black Hole Mass Function Across Cosmic Times I. Stellar Black Holes and Light Seed Distribution

This is the first paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the intermediate to the (super)massive regime. In the present work we focus on stellar BHs and provide an ab-initio computation of their mass function across cosmic times. Specifically, we exploit the state-of-the-art stellar and binary evolutionary code \texttt{SEVN}, and couple its outputs with redshift-dependent galaxy statistics and empirical scaling relations involving galaxy metallicity, star-formation rate and stellar mass. The resulting relic mass function ${\rm d}N/{\rm d}V{\rm d}\log m_\bullet$ as a function of the BH mass $m_\bullet$ features a rather flat shape up to $m_\bullet\approx 50\, M_\odot$ and then a log-normal decline for larger masses, while its overall normalization at a given mass increases with decreasing redshift. We highlight the contribution to the local mass function from isolated stars evolving into BHs and from binary stellar systems ending up in single or binary BHs. We also include the distortion on the mass function induced by binary BH mergers, finding that it has a minor effect at the high-mass end. We estimate a local stellar BH relic mass density of $ρ_\bullet\approx 5\times 10^7\, M_\odot$ Mpc$^{-3}$, which exceeds by more than two orders of magnitude that in supermassive BHs; this translates into an energy density parameter $Ω_\bullet\approx 4\times 10^{-4}$, implying that the total mass in stellar BHs amounts to $\lesssim 1\%$ of the local baryonic matter. We show how our mass function for merging BH binaries compares with the recent estimates from gravitational wave observations by LIGO/Virgo, and discuss the possible implications for dynamical formation of BH binaries in dense environments like star clusters. [abridged]

preprint2020arXiv

Carbon star formation as seen through the non-monotonic initial-final mass relation

The initial-final mass relation (IFMR) links the birth mass of a star to the mass of the compact remnant left at its death. While the relevance of the IFMR across astrophysics is universally acknowledged, not all of its fine details have yet been resolved. A new analysis of a few carbon-oxygen white dwarfs in old open clusters of the Milky Way led us to identify a kink in the IFMR, located over a range of initial masses, $1.65 \lesssim M_{\rm i}/M_{\odot} \lesssim 2.10$. The kink's peak in WD mass of $\approx 0.70-0.75 \, M_{\odot}$ is produced by stars with $M_{\rm i} \simeq 1.8 - 1.9 \, M_{\odot}$, corresponding to ages of about $1.8 - 1.7 $ Gyr. Interestingly, this peak coincides with the initial mass limit between low-mass stars that develop a degenerate helium core after central hydrogen exhaustion, and intermediate-mass stars that avoid electron degeneracy. We interpret the IFMR kink as the signature of carbon star formation in the Milky Way. This finding is critical to constraining the evolution and chemical enrichment of low-mass stars, and their impact on the spectrophotometric properties of galaxies.

preprint2020arXiv

Constraining the thermally pulsing asymptotic giant branch phase with resolved stellar populations in the Large Magellanic Cloud

Reliable models of the thermally pulsing asymptotic giant branch (TP-AGB) phase are of critical importance across astrophysics, including our interpretation of the spectral energy distribution of galaxies, cosmic dust production, and enrichment of the interstellar medium. With the aim of improving sets of stellar isochrones that include a detailed description of the TP-AGB phase, we extend our recent calibration of the AGB population in the Small Magellanic Cloud (SMC) to the more metal rich Large Magellanic Cloud (LMC). We model the LMC stellar populations with the TRILEGAL code, using the spatially-resolved star formation history derived from the VISTA survey. We characterize the efficiency of the third dredge-up by matching the star counts and the $K_{\rm s}$-band luminosity functions of the AGB stars identified in the LMC. In line with previous findings, we confirm that, compared to the SMC, the third dredge-up in AGB stars of the LMC is somewhat less efficient, as a consequence of the higher metallicity. The predicted range of initial mass of C-rich stars is between $M_{\rm i} \approx 1.7 - 3~\mathrm{M}_{\odot}$ at $Z_{\rm i} = 0.008$. We show how the inclusion of new opacity data in the carbon star spectra will improve the performance of our models. We discuss the predicted lifetimes, integrated luminosities and mass-loss rate distributions of the calibrated models. The results of our calibration are included in updated stellar isochrones publicly available.

preprint2020arXiv

Formation of GW190521 from stellar evolution: the impact of the hydrogen-rich envelope, dredge-up and $^{12}$C($α$, $γ$)$^{16}$O rate on the pair-instability black hole mass gap

Pair-instability (PI) is expected to open a gap in the mass spectrum of black holes (BHs) between $\approx{}40-65$ M$_\odot$ and $\approx{}120$ M$_\odot$. The existence of the mass gap is currently being challenged by the detection of GW190521, with a primary component mass of $85^{+21}_{-14}$ M$_{\odot}$. Here, we investigate the main uncertainties on the PI mass gap: the $^{12}$C($α$, $γ$)$^{16}$O reaction rate and the H-rich envelope collapse. With the standard $^{12}$C($α$, $γ$)$^{16}$O rate, the lower edge of the mass gap can be 70 M$_\odot$ if we allow for the collapse of the residual H-rich envelope at metallicity $Z\leq{}0.0003$. Adopting the uncertainties given by the STARLIB database, for models computed with the $^{12}$C($α$, $γ$)$^{16}$O rate $-1\, σ$, we find that the PI mass gap ranges between $\approx{}80$ M$_\odot$ and $\approx{}150$ M$_\odot$. Stars with $M_{\rm ZAMS}>110$ M$_\odot$ may experience a deep dredge-up episode during the core helium-burning phase, that extracts matter from the core enriching the envelope. As a consequence of the He-core mass reduction, a star with $M_{\rm ZAMS} =160$ M$_\odot$ may avoid the PI and produce a BH of 150 M$_\odot$. In the $-2\,{}σ$ case, the PI mass gap ranges from 92 M$_\odot$ to 110 M$_\odot$. Finally, in models computed with $^{12}$C($α$, $γ$)$^{16}$O $-3\,{}σ$, the mass gap is completely removed by the dredge-up effect. The onset of this dredge-up is particularly sensitive to the assumed model for convection and mixing. The combined effect of H-rich envelope collapse and low $^{12}$C($α$, $γ$)$^{16}$O rate can lead to the formation of BHs with masses consistent with the primary component of GW190521.

preprint2020arXiv

Impact of the Rotation and Compactness of Progenitors on the Mass of Black Holes

We investigate the impact of stellar rotation on the formation of black holes (BHs), by means of our population-synthesis code SEVN. Rotation affects the mass function of BHs in several ways. In massive metal-poor stars, fast rotation reduces the minimum zero-age main sequence (ZAMS) mass for a star to undergo pair instability and pulsational pair instability. Moreover, stellar winds are enhanced by rotation, peeling-off the entire hydrogen envelope. As a consequence of these two effects, the maximum BH mass we expect from the collapse a rotating metal-poor star is only $\sim{}45$ M$_\odot$, while the maximum mass of a BH born from a non-rotating star is $\sim{}60$ M$_\odot$. Furthermore, stellar rotation reduces the minimum ZAMS mass for a star to collapse into a BH from $\sim{}18-25$ M$_\odot$ to $\sim{}13-18$ M$_\odot$. Finally, we have investigated the impact of different core-collapse supernova (CCSN) prescriptions on our results. While the threshold value of compactness for direct collapse and the fallback efficiency strongly affect the minimum ZAMS mass for a star to collapse into a BH, the fraction of hydrogen envelope that can be accreted onto the final BH is the most important ingredient to determine the maximum BH mass. Our results confirm that the interplay between stellar rotation, CCSNe and pair instability plays a major role in shaping the BH mass spectrum.

preprint2020arXiv

PHAT XX. AGB stars and other cool giants in M31 star clusters

The presence of AGB stars in clusters provides key constraints for stellar models, as has been demonstrated with historical data from the Magellanic Clouds. In this work, we look for candidate AGB stars in M31 star clusters from the Panchromatic Hubble Andromeda Treasury (PHAT) survey. Our photometric criteria selects stars brighter than the tip of the red giant branch, which includes the bulk of the thermally-pulsing AGB stars as well as early-AGB stars and other luminous cool giants expected in young stellar populations (e.g. massive red supergiants, and intermediate-mass red helium-burning stars). The AGB stars can be differentiated, a posteriori, using the ages already estimated for our cluster sample. 937 candidates are found within the cluster aperture radii, half (450) of which are very likely cluster members. Cross-matching with additional databases reveals two carbon stars and ten secure variables among them. The field-corrected age distribution reveals the presence of young supergiants peaking at ages smaller than 100 Myr, followed by a long tail of AGB stars extending up to the oldest possible ages. This long tail reveals the general decrease in the numbers of AGB stars from initial values of 50e-6/Msun at 100 Myr down to 5e-6/Msun at 10 Gyr. Theoretical models of near-solar metallicity reproduce this general trend, although with localized discrepancies over some age intervals, whose origin is not yet identified. The entire catalogue is released together with finding charts to facilitate follow-up studies.

preprint2019arXiv

Impact of AGN feedback on galaxies and their multiphase ISM across cosmic time

We present simulations of galaxy formation, based on the GADGET-3 code, in which a sub-resolution model for star formation and stellar feedback is interfaced with a new model for AGN feedback. Our sub-resolution model describes a multiphase ISM, accounting for hot and cold gas within the same resolution element: we exploit this feature to investigate the impact of coupling AGN feedback energy to the different phases of the ISM over cosmic time. Our fiducial model considers that AGN feedback energy coupling is driven by the covering factors of the hot and cold phases. We perform a suite of cosmological hydrodynamical simulations of disc galaxies ($M_{\rm halo, \, DM} \simeq 2 \cdot 10^{12}$ M$_{\odot}$, at $z=0$), to investigate: $(i)$ the effect of different ways of coupling AGN feedback energy to the multiphase ISM; $(ii)$ the impact of different prescriptions for gas accretion (i.e. only cold gas, both cold and hot gas, with the additional possibility of limiting gas accretion from cold gas with high angular momentum); $(iii)$ how different models of gas accretion and coupling of AGN feedback energy affect the coevolution of supermassive BHs and their host galaxy. We find that at least a share of the AGN feedback energy has to couple with the diffuse gas, in order to avoid an excessive growth of the BH mass. When the BH only accretes cold gas, it experiences a growth that is faster than in the case in which both cold and hot gas are accreted. If the accretion of cold gas with high angular momentum is reduced, the BH mass growth is delayed, the BH mass at $z=0$ is reduced by up to an order of magnitude, and the BH is prevented from accreting below $z \lesssim 2$, when the galaxy disc forms.

preprint2016arXiv

Constraining dust properties in Circumstellar Envelopes of C-stars in the Small Magellanic Cloud: optical constants and grain size of Carbon dust

We present a new approach aimed at constraining the typical size and optical properties of carbon dust grains in Circumstellar envelopes (CSEs) of carbon-rich stars (C-stars) in the Small Magellanic Cloud (SMC). To achieve this goal, we apply our recent dust growth description, coupled with a radiative transfer code to the CSEs of C-stars evolving along the TP-AGB, for which we compute spectra and colors. Then we compare our modeled colors in the near- and mid-infrared (NIR and MIR) bands with the observed ones, testing different assumptions in our dust scheme and employing several data sets of optical constants for carbon dust available in the literature. Different assumptions adopted in our dust scheme change the typical size of the carbon grains produced. We constrain carbon dust properties by selecting the combination of grain size and optical constants which best reproduces several colors in the NIR and MIR at the same time. The different choices of optical properties and grain size lead to differences in the NIR and MIR colors greater than two magnitudes in some cases. We conclude that the complete set of observed NIR and MIR colors are best reproduced by small grains, with sizes between $\sim$0.035 and $\sim$0.12~$μ$m, rather than by large grains between $\sim0.2$ and $0.7$~$μ$m. The inability of large grains to reproduce NIR and MIR colors seems independent of the adopted optical data set. We also find a possible trend of the grain size with mass-loss and/or carbon excess in the CSEs of these stars.

preprint2016arXiv

Dynamics of tidally captured planets in the Galactic Center

Recent observations suggest ongoing planet formation in the innermost parsec of our Galaxy. The super-massive black hole (SMBH) might strip planets or planetary embryos from their parent star, bringing them close enough to be tidally disrupted. We investigate the chance of planet tidal captures by running three-body encounters of SMBH-star-planet systems with a high-accuracy regularized code. We show that tidally captured planets have orbits close to those of their parent star. We conclude that the final periapsis distance of the captured planet from the SMBH will be much larger than 200 AU, unless its parent star was already on a highly eccentric orbit.

preprint2016arXiv

Dynamics of tidally captured planets in the Galactic Center

Recent observations suggest ongoing planet formation in the innermost parsec of the Galactic center (GC). The super-massive black hole (SMBH) might strip planets or planetary embryos from their parent star, bringing them close enough to be tidally disrupted. Photoevaporation by the ultraviolet field of young stars, combined with ongoing tidal disruption, could enhance the near-infrared luminosity of such starless planets, making their detection possible even with current facilities. In this paper, we investigate the chance of planet tidal captures by means of high-accuracy N-body simulations exploiting Mikkola's algorithmic regularization. We consider both planets lying in the clockwise (CW) disk and planets initially bound to the S-stars. We show that tidally captured planets remain on orbits close to those of their parent star. Moreover, the semi-major axis of the planet orbit can be predicted by simple analytic assumptions in the case of prograde orbits. We find that starless planets that were initially bound to CW disk stars have mild eccentricities and tend to remain in the CW disk. However, we speculate that angular momentum diffusion and scattering with other young stars in the CW disk might bring starless planets on low-angular momentum orbits. In contrast, planets initially bound to S-stars are captured by the SMBH on highly eccentric orbits, matching the orbital properties of the G1 and G2 clouds. Our predictions apply not only to planets but also to low-mass stars initially bound to the S-stars and tidally captured by the SMBH.

preprint2016arXiv

Evolution of Thermally Pulsing Asymptotic Giant Branch Stars V: Constraining the Mass Loss and Lifetimes of Intermediate Mass, Low Metallicity AGB Stars

Thermally-Pulsing Asymptotic Giant Branch (TP-AGB) stars are relatively short lived (less than a few Myr), yet their cool effective temperatures, high luminosities, efficient mass-loss and dust production can dramatically effect the chemical enrichment histories and the spectral energy distributions of their host galaxies. The ability to accurately model TP-AGB stars is critical to the interpretation of the integrated light of distant galaxies, especially in redder wavelengths. We continue previous efforts to constrain the evolution and lifetimes of TP-AGB stars by modeling their underlying stellar populations. Using Hubble Space Telescope (HST) optical and near-infrared photometry taken of 12 fields of 10 nearby galaxies imaged via the ACS Nearby Galaxy Survey Treasury and the near-infrared HST/SNAP follow-up campaign, we compare the model and observed TP-AGB luminosity functions as well as the number ratio of TP-AGB to red giant branch stars. We confirm the best-fitting mass-loss prescription, introduced by Rosenfield et al. 2014, in which two different wind regimes are active during the TP-AGB, significantly improves models of many galaxies that show evidence of recent star formation. This study extends previous efforts to constrain TP-AGB lifetimes to metallicities ranging -1.59 < [Fe/H] < -0.56 and initial TP-AGB masses up to ~ 4 Msun, which include TP-AGB stars that undergo hot-bottom burning.

preprint2016arXiv

Synthetic photometry for M and K giants and stellar evolution: hydrostatic dust-free model atmospheres and chemical abundances

Based on a grid of hydrostatic spherical COMARCS models for cool stars we have calculated observable properties of these objects, which will be mainly used in combination with stellar evolution tracks and population synthesis tools. The high resolution opacity sampling and low resolution convolved spectra as well as bolometric corrections for a large number of filter systems are made electronically available. We exploit those data to study the effect of mass, C/O ratio and nitrogen abundance on the photometry of K and M giants. Depending on effective temperature, surface gravity and the chosen wavelength ranges variations of the investigated parameters cause very weak to moderate and, in the case of C/O values close to one, even strong shifts of the colours. For the usage with stellar evolution calculations they will be treated as correction factors applied to the results of an interpolation in the main quantities. When we compare the synthetic photometry to observed relations and to data from the Galactic Bulge, we find in general a good agreement. Deviations appear for the coolest giants showing pulsations, mass loss and dust shells, which cannot be described by hydrostatic models.

preprint2015arXiv

Envelope Overshooting in Low Metallicity Intermediate- and High-mass Stars: a test with the Sagittarius Dwarf Irregular Galaxy

We check the performance of the {\sl\,PARSEC} tracks in reproducing the blue loops of intermediate age and young stellar populations at very low metallicity. We compute new evolutionary {\sl\,PARSEC} tracks of intermediate- and high-mass stars from 2\Msun to 350\Msun with enhanced envelope overshooting (EO), EO=2\HP and 4\HP, for very low metallicity, Z=0.0005. The input physics, including the mass-loss rate, has been described in {\sl\,PARSEC}~V1.2 version. By comparing the synthetic color-magnitude diagrams (CMDs) obtained from the different sets of models with envelope overshooting EO=0.7\HP (the standard {\sl\,PARSEC} tracks), 2\HP and 4\HP, with deep observations of the Sagittarius dwarf irregular galaxy (SagDIG), we find an overshooting scale EO=2\HP to best reproduce the observed loops. This result is consistent with that obtained by \citet{Tang_etal14} for Z in the range 0.001-0.004. We also discuss the dependence of the blue loop extension on the adopted instability criterion and find that, contrary to what stated in literature, the Schwarzschild criterion, instead of the Ledoux criterion, favours the development of blue loops. Other factors that could affect the CMD comparisons such as differential internal extinction or the presence of binary systems are found to have negligible effects on the results. We thus confirm that, in presence of core overshooting during the H-burning phase, a large envelope overshooting is needed to reproduce the main features of the central He-burning phase of intermediate- and high-mass stars.

preprint2015arXiv

Lithium evolution in metal-poor stars: from Pre-Main Sequence to the Spite plateau

Lithium abundance derived in metal-poor main sequence stars is about three times lower than the value of primordial Li predicted by the standard Big Bang nucleosynthesis when the baryon density is taken from the CMB or the deuterium measurements. This disagreement is generally referred as the lithium problem. We here reconsider the stellar Li evolution from the pre-main sequence to the end of the main sequence phase by introducing the effects of convective overshooting and residual mass accretion. We show that $^7$Li could be significantly depleted by convective overshooting in the pre-main sequence phase and then partially restored in the stellar atmosphere by a tail of matter accretion which follows the Li depletion phase and that could be regulated by EUV photo-evaporation. By considering the conventional nuclear burning and microscopic diffusion along the main sequence we can reproduce the Spite plateau for stars with initial mass $m_0=0.62 - 0.80 M_{\odot}$, and the Li declining branch for lower mass dwarfs, e.g, $m_0=0.57 - 0.60 M_{\odot}$, for a wide range of metallicities (Z=0.00001 to Z=0.0005), starting from an initial Li abundance $A({\rm Li}) =2.72$. This environmental Li evolution model also offers the possibility to interpret the decrease of Li abundance in extremely metal-poor stars, the Li disparities in spectroscopic binaries and the low Li abundance in planet hosting stars.

preprint2015arXiv

OI and CaII observations in intermediate redshift quasars

We present an unprecedented spectroscopic survey of the CaII triplet + OI for a sample of 14 luminous ($-$26 $\gtrsim$ M$_V$ $\gtrsim$ $-$29), intermediate redshift (0.85 $\lesssim$ $z$ $\lesssim$ 1.65) quasars. The ISAAC spectrometer at ESO VLT allowed us to cover the CaII NIR spectral region redshifted into the H and K windows. We describe in detail our data analysis which enabled us to detect CaII triplet emission in all 14 sources (with the possible exception of HE0048-2804) and to retrieve accurate line widths and fluxes of the triplet and OI $λ$8446. The new measurements show trends consistent with previous lower $z$ observations, indicating that CaII and optical FeII emission are probably closely related. The ratio between the CaII triplet and the optical FeII blend at $λ$4570 $Å$ is apparently systematically larger in our intermediate redshift sample relative to a low-$z$ control sample. Even if this result needs a larger sample for adequate interpretation, higher CaII/optical FeII should be associated with recent episodes of star formation in the intermediate redshift quasars and, at least in part, explain an apparent correlation of CaII triplet equivalent width with $z$ and $L$. The CaII triplet measures yield significant constraints on the emitting region density and ionization parameter, implying CaII triplet emission from log n$_H$ $\gtrsim$ 11 [cm$^{-3}$] and ionization parameter log $U$ $\lesssim$ 1.5. Line width and intensity ratios suggest properties consistent with emission from the outer part of a high density broad line region (a line emitting accretion disk?).

preprint2015arXiv

On the Interpretation of Sub-Giant Branch Morphologies of Intermediate-Age Star Clusters with Extended Main Sequence Turnoffs

Recent high-quality photometry of many star clusters in the Magellanic Clouds with ages of 1$\,-\,$2 Gyr revealed main sequence turnoffs (MSTOs) that are significantly wider than can be accounted for by a simple stellar population (SSP). Such extended MSTOs (eMSTOs) are often interpreted in terms of an age spread of several $10^8$ yr, challenging the traditional view of star clusters as being formed in a single star formation episode. Li et al. and Bastian & Niederhofer recently investigated the sub-giant branches (SGBs) of NGC 1651, NGC 1806, and NGC 1846, three star clusters in the Large Magellanic Cloud (LMC) that exhibit an eMSTO. They argued that the SGB of these star clusters can be explained only by a SSP. We study these and two other similar star clusters in the LMC, using extensive simulations of SSPs including unresolved binaries. We find that the shapes of the cross-SGB profiles of all star clusters in our sample are in fact consistent with their cross-MSTO profiles when the latter are interpreted as age distributions. Conversely, SGB morphologies of star clusters with eMSTOs are found to be inconsistent with those of simulated SSPs. Finally, we create PARSEC isochrones from tracks featuring a grid of convective overshoot levels and a very fine grid of stellar masses. A comparison of the observed photometry with these isochrones shows that the morphology of the red clump (RC) of such star clusters is also consistent with that implied by their MSTO in the age spread scenario. We conclude that the SGB and RC morphologies of star clusters featuring eMSTOs are consistent with the scenario in which the eMSTOs are caused by a distribution of stellar ages.

preprint2015arXiv

PARSEC evolutionary tracks of massive stars up to $350 M_\odot$ at metallicities 0.0001$\leq Z \leq$0.04

We complement the PARSEC data base of stellar evolutionary tracks with new models of massive stars, from the pre-main sequence phase to the central carbon ignition. We consider a broad range of metallicities, 0.0001$\leq Z \leq$0.04 and initial masses up to $M_{\rm ini}=350\,M_\odot$. The main difference with respect to our previous models of massive stars is the adoption of a recent formalism accounting for the mass-loss enhancement when the ratio of the stellar to the Eddington luminosity, $Γ_e$, approaches unity. With this new formalism, the models are able to reproduce the Humphreys-Davidson limit observed in the Galactic and Large Magellanic Cloud colour-magnitude diagrams, without an ad hoc mass-loss enhancement. We also follow the predictions of recent wind models indicating that the metallicity dependence of the mass-loss rates becomes shallower when $Γ_e$ approaches unity. We thus find that the more massive stars may suffer from substantial mass-loss even at low metallicity. We also predict that the Humphreys-Davidson limit should become brighter at decreasing metallicity. We supplement the evolutionary tracks with new tables of theoretical bolometric corrections, useful to compare tracks and isochrones with the observations. For this purpose, we homogenize existing stellar atmosphere libraries of hot and cool stars (PoWR, ATLAS9 and Phoenix) and we add, where needed, new atmosphere models computed with WM-basic. The mass, age and metallicity grids are fully adequate to perform detailed investigations of the properties of very young stellar systems, both in local and distant galaxies. The new tracks supersede the previous old Padova models of massive stars.

preprint2015arXiv

Predictions for Ultra-Deep Radio Counts of Star-Forming Galaxies

We have worked out predictions for the radio counts of star-forming galaxies down to nJy levels, along with redshift distributions down to the detection limits of the phase 1 Square Kilometer Array MID telescope (SKA1-MID) and of its precursors. Such predictions were obtained by coupling epoch dependent star formation rate (SFR) functions with relations between SFR and radio (synchrotron and free-free) emission. The SFR functions were derived taking into account both the dust obscured and the unobscured star-formation, by combining far-infrared (FIR), ultra-violet (UV) and H_alpha luminosity functions up to high redshifts. We have also revisited the South Pole Telescope (SPT) counts of dusty galaxies at 95\,GHz performing a detailed analysis of the Spectral Energy Distributions (SEDs). Our results show that the deepest SKA1-MID surveys will detect high-z galaxies with SFRs two orders of magnitude lower compared to Herschel surveys. The highest redshift tails of the distributions at the detection limits of planned SKA1-MID surveys comprise a substantial fraction of strongly lensed galaxies. We predict that a survey down to 0.25 microJy at 1.4 GHz will detect about 1200 strongly lensed galaxies per square degree, at redshifts of up to 10. For about 30% of them the SKA1-MID will detect at least 2 images. The SKA1-MID will thus provide a comprehensive view of the star formation history throughout the re-ionization epoch, unaffected by dust extinction. We have also provided specific predictions for the EMU/ASKAP and MIGHTEE/MeerKAT surveys.

preprint2015arXiv

The AGB bump: a calibrator for the core mixing

The efficiency of convection in stars affects many aspects of their evolution and remains one of the key-open questions in stellar modelling. In particular, the size of the mixed core in core-He-burning low-mass stars is still uncertain and impacts the lifetime of this evolutionary phase and, e.g., the C/O profile in white dwarfs. One of the known observables related to the Horizontal Branch (HB) and Asymptotic Giant Branch (AGB) evolution is the AGB bump. Its luminosity depends on the position in mass of the helium-burning shell at its first ignition, that is affected by the extension of the central mixed region. In this preliminary work we show how various assumptions on near-core mixing and on the thermal stratification in the overshooting region affect the luminosity of the AGB bump, as well as the period spacing of gravity modes in core-He-burning models.

preprint2015arXiv

The influence of dense gas rings on the dynamics of a stellar disk in the Galactic center

The Galactic center hosts several hundred early-type stars, about 20% of which lie in the so-called clockwise disk, while the remaining 80% do not belong to any disks. The circumnuclear ring (CNR), a ring of molecular gas that orbits the supermassive black hole (SMBH) with a radius of 1.5 pc, has been claimed to induce precession and Kozai-Lidov oscillations onto the orbits of stars in the innermost parsec. We investigate the perturbations exerted by a gas ring on a nearly-Keplerian stellar disk orbiting a SMBH by means of combined direct N-body and smoothed particle hydrodynamics simulations. We simulate the formation of gas rings through the infall and disruption of a molecular gas cloud, adopting different inclinations between the infalling gas cloud and the stellar disk. We find that a CNR-like ring is not efficient in affecting the stellar disk on a timescale of 3 Myr. In contrast, a gas ring in the innermost 0.5 pc induces precession of the longitude of the ascending node Omega, significantly affecting the stellar disk inclination. Furthermore, the combined effect of two-body relaxation and Omega-precession drives the stellar disk dismembering, displacing the stars from the disk. The impact of precession on the star orbits is stronger when the stellar disk and the inner gas ring are nearly coplanar. We speculate that the warm gas in the inner cavity might have played a major role in the evolution of the clockwise disk.

preprint2015arXiv

The mass spectrum of compact remnants from the PARSEC stellar evolution tracks

The mass spectrum of stellar-mass black holes (BHs) is highly uncertain. Dynamical mass measurements are available only for few ($\sim{}10$) BHs in X-ray binaries, while theoretical models strongly depend on the hydrodynamics of supernova (SN) explosions and on the evolution of massive stars. In this paper, we present and discuss the mass spectrum of compact remnants that we obtained with SEVN, a new public population-synthesis code, which couples the PARSEC stellar evolution tracks with up-to-date recipes for SN explosion (depending on the Carbon-Oxygen mass of the progenitor, on the compactness of the stellar core at pre-SN stage, and on a recent two-parameter criterion based on the dimensionless entropy per nucleon at pre-SN stage). SEVN can be used both as a stand-alone code and in combination with direct-summation N-body codes (Starlab, HiGPUs). The PARSEC stellar evolution tracks currently implemented in SEVN predict significantly larger values of the Carbon-Oxygen core mass with respect to previous models. For most of the SN recipes we adopt, this implies substantially larger BH masses at low metallicity ($\leq{}2\times{}10^{-3}$), than other population-synthesis codes. The maximum BH mass found with SEVN is $\sim{}$25, 60 and 130 M$_{\odot}$ at metallicity $Z =2 \times{} 10^{-2}$ , $2 \times{}10^{-3}$ and $2\times{} 10^{-4}$ , respectively. Mass loss by stellar winds plays a major role in determining the mass of BHs for very massive stars ($\geq{}90$ M$_\odot{}$), while the remnant mass spectrum depends mostly on the adopted SN recipe for lower progenitor masses. We discuss the implications of our results for the transition between NS and BH mass, and for the expected number of massive BHs (with mass $>25$ M$_\odot{}$) as a function of metallicity.

preprint2015arXiv

The VMC Survey - XIV. First results on the look-back time star-formation rate tomography of the Small Magellanic Cloud

We analyse deep images from the VISTA survey of the Magellanic Clouds in the YJKs filters, covering 14 sqrdeg (10 tiles), split into 120 subregions, and comprising the main body and Wing of the Small Magellanic Cloud (SMC). We apply a colour--magnitude diagram reconstruction method that returns their best-fitting star formation rate SFR(t), age-metallicity relation (AMR), distance and mean reddening, together with 68% confidence intervals. The distance data can be approximated by a plane tilted in the East-West direction with a mean inclination of 39 deg, although deviations of up to 3 kpc suggest a distorted and warped disk. After assigning to every observed star a probability of belonging to a given age-metallicity interval, we build high-resolution population maps. These dramatically reveal the flocculent nature of the young star-forming regions and the nearly smooth features traced by older stellar generations. They document the formation of the SMC Wing at ages <0.2 Gyr and the peak of star formation in the SMC Bar at 40 Myr. We clearly detect periods of enhanced star formation at 1.5 Gyr and 5 Gyr. The former is possibly related to a new feature found in the AMR, which suggests ingestion of metal-poor gas at ages slightly larger than 1 Gyr. The latter constitutes a major period of stellar mass formation. We confirm that the SFR(t) was moderately low at even older ages.

preprint2015arXiv

Uncertainties on near-core mixing in red-clump stars: effects on the period spacing and on the luminosity of the AGB bump

Low-mass stars in the He-core-burning phase (HeCB) play a major role in stellar, galactic, and extragalactic astrophysics. The ability to predict accurately the properties of these stars, however, depends on our understanding of convection, which remains one of the key open questions in stellar modelling. We argue that the combination of the luminosity of the AGB bump (AGBb) and the period spacing of gravity modes (DP) during the HeCB phase, provides us with a decisive test to discriminate between competing models of these stars. We use the MESA, BaSTI, and PARSEC stellar evolution codes to model a typical giant star observed by Kepler. We explore how various near-core-mixing scenarios affect the predictions of the above-mentioned constraints, and we find that DP depends strongly on the prescription adopted. Moreover we show that the detailed behaviour of DP shows the signature of sharp variations in the Brunt-Vaisala frequency, which could potentially give additional information about near-core features. We find evidence for the AGBb among Kepler targets, and a first comparison with observations shows that, even if standard models are able to reproduce the luminosity distribution, no standard model can account for satisfactorily the period spacing of HeCB stars. Our analysis allows us to outline a candidate model to describe simultaneously the two observed distributions: a model with a moderate overshooting region characterized by an adiabatic thermal stratification. This prescription will be tested in the future on cluster stars, to limit possible observational biases.

preprint2014arXiv

A physical model for the evolving UV luminosity function of high redshift galaxies and their contribution to the cosmic reionization

[Abridged] We present a physical model for the evolution of the ultraviolet (UV) luminosity function (LF) of high-z galaxies taking into account in a self-consistent way their chemical evolution and the associated evolution of dust extinction. The model yields good fits of the UV and Lyman-alpha LFs at z>~2. The weak evolution of both LFs between z=2 and z=6 is explained as the combined effect of the negative evolution of the halo mass function, of the increase with redshift of the star formation efficiency, and of dust extinction. The slope of the faint end of the UV LF is found to steepen with increasing redshift, implying that low luminosity galaxies increasingly dominate the contribution to the UV background at higher and higher redshifts. The observed range of UV luminosities at high-z implies a minimum halo mass capable of hosting active star formation M_crit <~ 10^9.8 M_odot, consistent with the constraints from hydrodynamical simulations. From fits of Lyman-alpha LFs plus data on the luminosity dependence of extinction and from the measured ratios of non-ionizing UV to Lyman-continuum flux density for samples of z=~3 Lyman break galaxies and Lyman-alpha emitters, we derive a simple relationship between the escape fraction of ionizing photons and the star formation rate, impling larger escape fraction for less massive galaxies. Galaxies already represented in the UV LF (M_UV <~ -18) can keep the universe fully ionized up to z=~6, consistent with (uncertain) data pointing to a rapid drop of the ionization degree above z~6. On the other side, the electron scattering optical depth, tau_es, inferred from CMB experiments favor an ionization degree close to unity up to z=~9-10. Consistency with CMB data can be achieved if M_crit =~ 10^8.5 M_odot, implying that the UV LFs extend to M_UV =~ -13, although the corresponding tau_es is still on the low side of CMB-based estimates.

preprint2014arXiv

Constraining Mass-Loss & Lifetimes of Low Mass, Low Metallicity AGB Stars

The evolution and lifetimes of thermally pulsating asymptotic giant branch (TP-AGB) stars suffer from significant uncertainties. We present a detailed framework for constraining model luminosity functions of TP-AGB stars using resolved stellar populations. We show an example of this method that compares various TP-AGB mass-loss prescriptions that differ in their treatments of mass loss before the onset of dust-driven winds (pre-dust). We find that models with more efficient pre-dust driven mass loss produce results consistent with observations, as opposed to more canonical mass-loss models. Efficient pre-dust driven mass-loss predicts for [Fe/H] < -1.2, lower mass TP-AGB stars (M < 1 Msun) must have lifetimes less than about 1.2 Myr.

preprint2014arXiv

Dust production from sub-solar to super-solar metallicity in Thermally Pulsing Asymptotic Giant Branch Stars

We discuss the dust chemistry and growth in the circumstellar envelopes (CSEs) of Thermally Pulsing Asymptotic Giant Branch (TP-AGB) star models computed with the COLIBRI code, at varying initial mass and metallicity (Z=0.001, 0.008, 0.02, 0.04, 0.06). A relevant result of our analysis deals with the silicate production in M-stars. We show that, in order to reproduce the observed trend between terminal velocities and mass-loss rates in Galactic M-giants, one has to significantly reduce the efficiency of chemisputtering by H2 molecules, usually considered as the most effective dust destruction mechanism. This indication is also in agreement with the most recent laboratory results, which show that silicates may condense already at T=1400 K, instead than at Tcond=1000 K, as obtained by models that include chemisputtering. From the analysis of the total dust ejecta, we find that the total dust-to-gas ejecta of intermediate-mass stars are much less dependent on metallicity than usually assumed. In a broader context, our results are suitable to study the dust enrichment of the interstellar medium provided by TP-AGB stars in both nearby and high redshift galaxies.

preprint2014arXiv

Evolution of Thermally Pulsing Asymptotic Giant Branch Stars IV. Constraining Mass-Loss & Lifetimes of Low Mass, Low Metallicity AGB Stars

The evolution and lifetimes of thermally pulsating asymptotic giant branch (TP-AGB) stars suffer from significant uncertainties. In this work, we analyze the numbers and luminosity functions of TP-AGB stars in six quiescent, low metallicity ([Fe/H] $\lesssim -0.86$) galaxies taken from the ANGST sample, using HST photometry in both optical and near-infrared filters. The galaxies contain over 1000 TP-AGB stars (at least 60 per field). We compare the observed TP-AGB luminosity functions and relative numbers of TP-AGB and RGB stars, to models generated from different suites of TP-AGB evolutionary tracks after adopting star formation histories (SFH) derived from the HST deep optical observations. We test various mass-loss prescriptions that differ in their treatments of mass-loss before the onset of dust-driven winds (pre-dust). These comparisons confirm that pre-dust mass-loss is important, since models that neglect pre-dust mass-loss fail to explain the observed TP-AGB/RGB ratio or the luminosity functions. In contrast, models with more efficient pre-dust mass-loss produce results consistent with observations. We find that for [Fe/H]$\lesssim-0.86$, lower mass TP-AGB stars ($M\lesssim1M_\odot$) must have lifetimes of $\sim0.5$ Myr and higher masses ($M\lesssim 3M_\odot$) must have lifetimes $\lesssim 1.2$ Myr. In addition, assuming our best-fitting mass-loss prescription, we show that the third dredge up has no significant effect on TP-AGB lifetimes in this mass and metallicity range.

preprint2014arXiv

Extended Main Sequence Turnoffs in Intermediate-Age Star Clusters: A Correlation Between Turnoff Width and Early Escape Velocity

We present color-magnitude diagram analysis of deep Hubble Space Telescope imaging of a mass-limited sample of 18 intermediate-age (1 - 2 Gyr old) star clusters in the Magellanic Clouds, including 8 clusters for which new data was obtained. We find that ${\it all}$ star clusters in our sample feature extended main sequence turnoff (eMSTO) regions that are wider than can be accounted for by a simple stellar population (including unresolved binary stars). FWHM widths of the MSTOs indicate age spreads of 200-550 Myr. We evaluate dynamical evolution of clusters with and without initial mass segregation. Our main results are: (1) the fraction of red clump (RC) stars in secondary RCs in eMSTO clusters scales with the fraction of MSTO stars having pseudo-ages $\leq 1.35$ Gyr; (2) the width of the pseudo-age distributions of eMSTO clusters is correlated with their central escape velocity $v_{\rm esc}$, both currently and at an age of 10 Myr. We find that these two results are unlikely to be reproduced by the effects of interactive binary stars or a range of stellar rotation velocities. We therefore argue that the eMSTO phenomenon is mainly caused by extended star formation within the clusters; (3) we find that $v_{\rm esc} \geq 15$ km/s out to ages of at least 100 Myr for ${\it all}$ clusters featuring eMSTOs, while $v_{\rm esc} \leq 12$ km/s at all ages for two lower-mass clusters in the same age range that do ${\it not}$ show eMSTOs. We argue that eMSTOs only occur for clusters whose early escape velocities are higher than the wind velocities of stars that provide material from which second-generation stars can form. The threshold of 12-15 km/s is consistent with wind velocities of intermediate-mass AGB stars in the literature.

preprint2014arXiv

Improving PARSEC models for very low mass stars

Many stellar models present difficulties in reproducing basic observational relations of very low mass stars (VLMS), including the mass--radius relation and the optical colour--magnitudes of cool dwarfs. Here, we improve PARSEC models on these points. We implement the T--tau relations from PHOENIX BT-Settl model atmospheres as the outer boundary conditions in the PARSEC code, finding that this change alone reduces the discrepancy in the mass--radius relation from 8 to 5 per cent. We compare the models with multi--band photometry of clusters Praesepe and M67, showing that the use of T--tau relations clearly improves the description of the optical colours and magnitudes. But anyway, using both Kurucz and PHOENIX model spectra, model colours are still systematically fainter and bluer than the observations. We then apply a shift to the above T--tau relations, increasing from 0 at T_eff = 4730 K to ~14% at T_eff = 3160 K, to reproduce the observed mass--radius radius relation of dwarf stars. Taking this experiment as a calibration of the T--tau relations, we can reproduce the optical and near infrared CMDs of low mass stars in the old metal--poor globular clusters NGC6397 and 47Tuc, and in the intermediate--age and young solar--metallicity open clusters M67 and Praesepe. Thus, we extend PARSEC models using this calibration, providing VLMS models more suitable for the lower main sequence stars over a wide range of metallicities and wavelengths. Both sets of models are available on PARSEC webpage.

preprint2014arXiv

New PARSEC evolutionary tracks of massive stars at low metallicity: testing canonical stellar evolution in nearby star forming dwarf galaxies

We extend the {\sl\,PARSEC} library of stellar evolutionary tracks by computing new models of massive stars, from 14\Msun to 350\Msun. The input physics is the same used in the {\sl\,PARSEC}~V1.1 version, but for the mass-loss rate which is included by considering the most recent updates in literature. We focus on low metallicity, $Z$=0.001 and $Z$=0.004, for which the metal poor dwarf irregular star forming galaxies, Sextans A, WLM and NCG6822, provide simple but powerful workbenches. The models reproduce fairly well the observed CMDs but the stellar colour distributions indicate that the predicted blue loop is not hot enough in models with canonical extent of overshooting. In the framework of a mild extended mixing during central hydrogen burning, the only way to reconcile the discrepancy is to enhance the overshooting at the base of the convective envelope (EO) during the first dredge-UP. The mixing scales required to reproduce the observed loops, EO=2\HP or EO=4\HP, are definitely larger than those derived from, e.g., the observed location of the RGB bump in low mass stars. This effect, if confirmed, would imply a strong dependence of the mixing scale below the formal Schwarzschild border, on the stellar mass or luminosity. Reproducing the features of the observed CMDs with standard values of envelope overshooting would require a metallicity significantly lower than the values measured in these galaxies. Other quantities, such as the star formation rate and the initial mass function, are only slightly sensitive to this effect. Future investigations will consider other metallicities and different mixing schemes.

preprint2014arXiv

The impact of metallicity-dependent mass loss versus dynamical heating on the early evolution of star clusters

We have run direct N-body simulations to investigate the impact of stellar evolution and dynamics on the structural properties of young massive (3x10^4 solar masses) star clusters (SCs) with different metallicities (Z=1, 0.1, 0.01 solar metallicity). Metallicity drives the mass loss by stellar winds and supernovae (SNe), with SCs losing more mass at high metallicity. We have simulated three sets of initial conditions, with different initial relaxation timescale. We find that the evolution of the half-mass radius of SCs depends on how fast two-body relaxation is with respect to the lifetime of massive stars. If core collapse is slow in comparison with stellar evolution, then mass loss by stellar winds and SNe is the dominant mechanism driving SC evolution, and metal-rich SCs expand more than metal-poor ones. In contrast, if core collapse occurs on a comparable timescale with respect to the lifetime of massive stars, then SC evolution depends on the interplay between mass loss and three-body encounters: dynamical heating by three-body encounters (mass loss by stellar winds and SNe) is the dominant process driving the expansion of the core in metal-poor (metal-rich) SCs. As a consequence, the half-mass radius of metal-poor SCs expands more than that of metal-rich ones. We also find core radius oscillations, which grow in number and amplitude as metallicity decreases.

preprint2014arXiv

Uncertainties in stellar evolution models: convective overshoot

In spite of the great effort made in the last decades to improve our understanding of stellar evolution, significant uncertainties remain due to our poor knowledge of some complex physical processes that require an empirical calibration, such as the efficiency of the interior mixing related to convective overshoot. Here we review the impact of convective overshoot on the evolution of stars during the main Hydrogen and Helium burning phases.

preprint2013arXiv

An extended main sequence turn-off in the Small Magellanic Cloud star cluster NGC411

Based on new observations with the Wide Field Camera 3 onboard the Hubble Space Telescope, we report the discovery of an extended main sequence turn-off (eMSTO) in the intermediate-age star cluster NGC411. This is the second case of an eMSTO being identified in a star cluster belonging to the Small Magellanic Cloud (SMC), after NGC419. Despite the present masses of these two SMC clusters differ by a factor of 4, the comparison between their colour--magnitude diagrams (CMD) shows striking similarities, especially regarding the shape of their eMSTOs. The loci of main CMD features are so similar that they can be well described, in a first approximation, by the same mean metallicity, distance and extinction. NGC411, however, presents merely a trace of secondary red clump as opposed to its prominent manifestation in NGC419. This could be due either to the small number statistics in NGC411, or by the star formation in NGC419 having continued for 60 Myr longer than in NGC411. Under the assumption that the eMSTOs are caused by different generations of stars at increasing age, both clusters are nearly coeval in their first episodes of star formation. The initial period of star formation, however, is slightly more marked in NGC419 than in NGC411. We discuss these findings in the context of possible scenarios for the origin of eMSTOs.

preprint2013arXiv

Evolution of Thermally Pulsing Asymptotic Giant Branch Stars I. The COLIBRI Code

We present the COLIBRI code for computing the evolution of stars along the TP-AGB phase. Compared to purely synthetic TP-AGB codes, COLIBRI relaxes a significant part of their analytic formalism in favour of a detailed physics applied to a complete envelope model, in which the stellar structure equations are integrated from the atmosphere down to the bottom of the hydrogen-burning shell. This allows to predict self-consistently: (i) the effective temperature, and more generally the convective envelope and atmosphere structures, correctly coupled to the changes in the surface chemical abundances and gas opacities; (ii) sphericity effects in the atmospheres; (iii) the core mass-luminosity relation and its break-down due to hot bottom burning (HBB) in the most massive AGB stars, (iv) the HBB nucleosynthesis via the solution of a complete nuclear network (pp chains, and the CNO, NeNa, MgAl cycles), including also the production of 7Li via the Cameron-Fowler beryllium transport mechanism; (v) the chemical composition of the pulse-driven convective zone; (vi) the onset and quenching of the third dredge-up, with a suitable temperature criterion. At the same time COLIBRI pioneers new techniques in the treatment of the physics of stellar interiors. It is the first evolutionary code ever to use accurate on-the-fly computation of the equation of state for roughly 800 atoms, ions, molecules, and of the Rosseland mean opacities throughout the deep envelope. Another distinguishing aspect of COLIBRI is its high computational speed. This feature is necessary for calibrating the uncertain parameters and processes that characterize the TP-AGB phase, a step of paramount importance for producing reliable stellar population synthesis models of galaxies up to high redshift. (abridged)

preprint2013arXiv

Evolution of Thermally Pulsing Asymptotic Giant Branch Stars II. Dust production at varying metallicity

We present the dust ejecta of the new stellar models for the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase computed with the COLIBRI code. We use a formalism of dust growth coupled with a stationary wind for both M and C-stars. In the original version of this formalism, the most efficient destruction process of silicate dust in M-giants is chemisputtering by H2 molecules. For these stars we find that dust grains can only form at relatively large radial distances (r~5 R*), where they cannot be efficiently accelerated, in agreement with other investigations. In the light of recent laboratory results, we also consider the alternative case that the condensation temperature of silicates is determined only by the competition between growth and free evaporation processes (i.e. no chemisputtering). With this latter approach we obtain dust condensation temperatures that are significantly higher (up to Tcond~1400 K) than those found when chemisputtering is included (Tcond~900 K), and in better agreement with condensation experiments. As a consequence, silicate grains can remain stable in inner regions of the circumstellar envelopes (r~2 R*), where they can rapidly grow and can be efficiently accelerated. With this modification, our models nicely reproduce the observed trend between terminal velocities and mass loss rates of Galactic M-giants. For C-stars the formalism is based on the homogeneous growth scheme where the key role is played by the carbon over oxygen excess. The models reproduce fairly well the terminal velocities of Galactic stars and there is no need to invoke changes in the standard assumptions. At decreasing metallicity the carbon excess becomes more pronounced and the efficiency of dust formation increases. This trend could be in tension with recent observational evidence in favour of a decreasing efficiency, at decreasing metallicity.

preprint2013arXiv

Evolution of Thermally Pulsing Asymptotic Giant Branch Stars III. Dust production at supersolar metallicities

We extend the formalism presented in our recent calculations of dust ejecta from the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase, to the case of super-solar metallicity stars. The TP-AGB evolutionary models are computed with the COLIBRI code. We adopt our preferred scheme for dust growth. For M-giants, we neglect chemisputtering by H$_2$ molecules and, for C-stars we assume a homogeneous growth scheme which is primarily controlled by the carbon over oxygen excess. At super-solar metallicities, dust forms more efficiently and silicates tend to condense significantly closer to the photosphere (r~1.5 R$_*$) - and thus at higher temperatures and densities - than at solar and sub-solar metallicities (r~2-3 R$_*$). In such conditions, the hypothesis of thermal decoupling between gas and dust becomes questionable, while dust heating due to collisions plays an important role. The heating mechanism delays dust condensation to slightly outer regions in the circumstellar envelope. We find that the same mechanism is not significant at solar and sub-solar metallicities. The main dust products at super-solar metallicities are silicates. We calculate the total dust ejecta and dust-to-gas ejecta, for various values of the stellar initial masses and initial metallicities Z=0.04, 0.06. Merging these new calculations with those for lower metallicities it turns out that, contrary to what often assumed, the total dust-to-gas ejecta of intermediate-mass stars exhibit only a weak dependence on the initial metal content.

preprint2013arXiv

Low ionization lines in high luminosity quasars: The calcium triplet

In order to investigate where and how low ionization lines are emitted in quasars we are studying a new collection of spectra of the CaII triplet at $λ$8498, $λ$8542, $λ$8662 observed with the Very Large Telescope (VLT) using the Infrared Spectrometer And Array Camera (ISAAC). Our sample involves luminous quasars at intermediate redshift for which CaII observations are almost nonexistent. We fit the CaII triplet and the OI $λ$8446 line using the H$β$ profile as a model. We derive constraints on the line emitting region from the relative strength of the CaII triplet, OI $λ$8446 and H$β$.

preprint2013arXiv

Red Giant evolution and specific problems

In spite of the great effort made in the last decades to improve our understanding of stellar evolution, significant uncertainties still remain due to our poor knowledge of some complex physical processes that still require an empirical calibration, such as the efficiency of convective heat transport and interior mixing. Here we will review the impact of these uncertainties on the evolution of red giant stars.

preprint2013arXiv

The insidious boosting of TP-AGB stars in intermediate-age Magellanic Cloud clusters

(Abridged) In the recent controversy about the role of TP-AGB stars in evolutionary population synthesis (EPS) models of galaxies, one particular aspect is puzzling: TP-AGB models aimed at reproducing the lifetimes and integrated fluxes of the TP-AGB phase in Magellanic Cloud (MC) clusters, when incorporated into EPS models, are found to overestimate the TP-AGB contribution in resolved star counts and integrated spectra of galaxies. In this paper, we call attention to a particular evolutionary aspect that in all probability is the main cause of this conundrum. As soon as stellar populations intercept the ages at which RGB stars first appear, a sudden change in the lifetime of the core He-burning phase causes a temporary boost in the production rate of subsequent evolutionary phases, including the TP-AGB. For a timespan of about 0.1 Gyr, triple TP-AGB branches develop at slightly different initial masses, causing their frequency and contribution to the integrated luminosity of the stellar population to increase by a factor of 2. The boost occurs just in the proximity of the expected peak in the TP-AGB lifetimes, and for ages of 1.6 Gyr. Coincidently, this relatively narrow age interval happens to contain the few very massive MC clusters that host most of the TP-AGB stars used to constrain stellar evolution and EPS models. This concomitance makes the AGB-boosting particularly insidious in the context of present EPS models. The effect brings about three main consequences. (1) Present estimates of the TP-AGB contribution to the integrated light of galaxies derived from MC clusters, are biased towards too large values. (2) The relative TP-AGB contribution of single-burst populations falling in this critical age range cannot be accurately derived by the fuel consumption theorem. (3) A careful revision of AGB star populations in intermediate-age MC clusters is urgently demanded.

preprint2013arXiv

The Star Formation History of Redshift z~2 Galaxies: The Role of The Infrared Prior

We build a sample of 298 spectroscopically-confirmed galaxies at redshift z~2, selected in the z-band from the GOODS-MUSIC catalog. By exploiting the rest frame 8 um luminosity as a proxy of the star formation rate (SFR) we check the accuracy of the standard SED-fitting technique, finding it is not accurate enough to provide reliable estimates of the galaxy physical parameters. We then develop a new SED-fitting method that includes the IR luminosity as a prior and a generalized Calzetti law with a variable RV . Then we exploit such a new method to re-analyze our galaxy sample, and to robustly determine SFRs, stellar masses and ages. We find that there is a general trend of increasing attenuation with the SFR. Moreover, we find that the SFRs range between a few to 1000 solar mass per year, the masses from one billion to 400 billion solar masses, while the ages from a few tens of Myr to more than 1 Gyr. We discuss how individual age easurements of highly attenuated objects indicate that dust must form within a few tens of Myr and be copious already at ~100 Myr. In addition, we find that low luminous galaxies harbor, on average, significantly older stellar populations and are also less massive than brighter ones; we discuss how these findings and the well known 'downsizing' scenario are consistent in a framework where less massive galaxies form first, but their star formation lasts longer. Finally, we find that the near-IR attenuation is not scarce for luminous objects, contrary to what is customarily assumed; we discuss how this affects the interpretation of the observed mass-to-light ratios.

preprint2013arXiv

The star formation history of the Large Magellanic Cloud star clusters NGC1846 and NGC1783

NGC1846 and NGC1783 are two massive star clusters in the Large Magellanic Cloud, hosting both an extended main sequence turn-off and a dual clump of red giants. They present similar masses but differ mainly in angular size. Starting from their high-quality ACS data in the F435W, F555W and F814W filters, and updated sets of stellar evolutionary tracks, we derive their star formation rates as a function of age, SFR(t), by means of the classical method of CMD reconstruction which is usually applied to nearby galaxies. The method confirms the extended periods of star formation derived from previous analysis of the same data. When the analysis is performed for a finer resolution in age, we find clear evidence for a 50-Myr long hiatus between the oldest peak in the SFR(t), and a second prolonged period of star formation, in both clusters. For the more compact cluster NGC1846, there seems to be no significant difference between the SFR(t) in the cluster centre and in an annulus with radii between 20 and 60 arcsec (from 4.8 to 15.4 pc). The same does not occur in the more extended NGC1783 cluster, where the outer ring (between 33 and 107 arcsec, from 8.0 to 25.9 pc) is found to be slightly younger than the centre. We also explore the best-fitting slope of the present-day mass function and binary fraction for the different cluster regions, finding hints of a varying mass function between centre and outer ring in NGC1783. These findings are discussed within the present scenarios for the formation of clusters with multiple turn-offs.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury I: Bright UV Stars in the Bulge of M31

As part of the Panchromatic Hubble Andromeda Treasury (PHAT) multi-cycle program, we observed a 12' \times 6.5' area of the bulge of M31 with the WFC3/UVIS filters F275W and F336W. From these data we have assembled a sample of \sim4000 UV-bright, old stars, vastly larger than previously available. We use updated Padova stellar evolutionary tracks to classify these hot stars into three classes: Post-AGB stars (P-AGB), Post-Early AGB (PE-AGB) stars and AGB-manqué stars. P-AGB stars are the end result of the asymptotic giant branch (AGB) phase and are expected in a wide range of stellar populations, whereas PE-AGB and AGB-manqué (together referred to as the hot post-horizontal branch; HP-HB) stars are the result of insufficient envelope masses to allow a full AGB phase, and are expected to be particularly prominent at high helium or α abundances when the mass loss on the RGB is high. Our data support previous claims that most UV-bright sources in the bulge are likely hot (extreme) horizontal branch stars (EHB) and their progeny. We construct the first radial profiles of these stellar populations, and show that they are highly centrally concentrated, even more so than the integrated UV or optical light. However, we find that this UV-bright population does not dominate the total UV luminosity at any radius, as we are detecting only the progeny of the EHB stars that are the likely source of the UVX. We calculate that only a few percent of MS stars in the central bulge can have gone through the HP-HB phase and that this percentage decreases strongly with distance from the center. We also find that the surface density of hot UV-bright stars has the same radial variation as that of low-mass X-ray binaries. We discuss age, metallicity, and abundance variations as possible explanations for the observed radial variation in the UV-bright population.

preprint2011arXiv

Tracing rejuvenation events in nearby S0 galaxies

With the aim of characterizing rejuvenation processes in early-type galaxies, we analyzed five barred S0 galaxies showing prominent outer ring in ultraviolet (UV) imaging. We analyzed GALEX far- (FUV) and near- (NUV) UV and optical data using stellar population models and estimated the age and the stellar mass of the entire galaxies and of the UV-bright ring structures. Outer rings consist of young (<200 Myr old) stellar populations, accounting for up to 70% of the FUV flux but containing only a few % of the total stellar mass. Integrated photometry of the whole galaxies places four of these objects on the green valley, indicating a globally evolving nature. We suggest such galaxy evolution is likely driven by bar induced instabilities, i.e. inner secular evolution, that conveys gas to the nucleus and to the outer rings. At the same time, HI observations of NGC 1533 and NGC 2962 suggest external gas re-fueling can play a role in the rejuvenation processes of such galaxies.

preprint2008arXiv

Evolution of asymptotic giant branch stars II. Optical to far-infrared isochrones with improved TP-AGB models

We present a large set of theoretical isochrones, whose distinctive features mostly reside on the greatly improved treatment of the thermally pulsing asymptotic giant branch (TP-AGB) phase. Essentially, we have coupled the TP-AGB tracks described in Paper I, at their stages of pre-flash quiescent H-shell burning, with the evolutionary tracks for the previous evolutionary phases from Girardi et al. (2000). Theoretical isochrones for any intermediate value of age and metallicity are then derived by interpolation in the grids. We take care that the isochrones keep, to a good level of detail, the several peculiarities present in these TP-AGB tracks. Theoretical isochrones are then converted to about 20 different photometric systems -- including traditional ground-based systems, and those of recent major wide-field surveys such as SDSS, OGLE, DENIS, 2MASS, UKIDSS, etc., -- by means of synthetic photometry applied to an updated library of stellar spectra, suitably extended to include C-type stars. Finally, we correct the predicted photometry by the effect of circumstellar dust during the mass-losing stages of the AGB evolution, which allows us to improve the results for the optical-to-infrared systems, and to simulate mid- and far-IR systems such as those of Spitzer and AKARI. Access to the data is provided both via a web repository of static tables (http://stev.oapd.inaf.it/dustyAGB07 and CDS), and via an interactive web interface (http://stev.oapd.inaf.it/cmd) that provides tables for any intermediate value of age and metallicity, for several photometric systems, and for different choices of dust properties.

preprint2001arXiv

Dust and Nebular Emission in Star Forming Galaxies

Star forming galaxies exhibit a variety of physical conditions, from quiescent normal spirals to the most powerful dusty starbursts. In order to study these complex systems, we need a suitable tool to analyze the information coming from observations at all wavelengths. We present a new spectro-photometric model which considers in a consistent way starlight as reprocessed by gas and dust. We discuss preliminary results to interpret some observed properties of VLIRGs.

preprint1999arXiv

The third dredge-up and the carbon star luminosity functions in the Magellanic Clouds

We investigate the formation of carbon stars as a function of the stellar mass and parent metallicity. Theoretical modelling is based on an improved scheme for treating the third dredge-up in synthetic calculations of thermally pulsing asymptotic giant branch (TP-AGB) stars. In this approach, the usual criterion (based on a constant minimum core mass for the occurrence of dredge-up, M_c^min) is replaced by one on the minimum temperature at the base of the convective envelope, T_b^dred, at the stage of the post-flash luminosity maximum. Envelope integrations then allow determination of M_c^min as a function of stellar mass, metallicity, and pulse strength (see Wood 1981), thus inferring if and when dredge-up first occurs. Moreover, the final possible shut down of the process is predicted. Extensive grids of TP-AGB models were computed using this scheme. We present and discuss the calibration of the two dredge-up parameters (lambda and T_b^dred) aimed at reproducing the carbon star luminosity function (CSLF) in the LMC. It turns out that the faint tail is almost insensitive to the history of star formation rate (SFR) in the parent galaxy (it is essentially determined by T_b^dred), in contrast to the bright wing which may be more affected by the details of the recent SFR. Once the faint end is reproduced, the peak location is a stringent calibrator of lambda. The best fit to the observed CSLF in the LMC is obtained with Z=0.008, lambda=0.50, log(T_b^dred)=6.4, and a constant SFR up to 5x10^8 yr ago. A good fit to the CSLF in the SMC is then easily derived from the Z=0.004 models, with a single choice of parameters, and a constant SFR over the entire significant age interval. The results are consistent with the theoretical expectation that the third dredge-up is more efficient at lower Zs.