The optimally-sampled galaxy-wide stellar initial mass function - Observational tests and the publicly available GalIMF code

Here we present a full description of the integrated galaxy-wide initial mass function (IGIMF) theory in terms of the optimal sampling and compare it with available observations. Optimal sampling is the method we use to discretize the IMF into stellar masses deterministically. Evidence has been indicating that nature may be closer to deterministic sampling as observations suggest a smaller scatter of various relevant observables than random sampling would give, which may result from a high level of self-regulation during the star formation process. The variation of the IGIMFs under various assumptions are documented. The results of the IGIMF theory are consistent with the empirical relation between the total mass of a star cluster and the mass of its most massive star, and the empirical relation between a galaxy's star formation rate (SFR) and the mass of its most massive cluster. Particularly, we note a natural agreement with the empirical relation between the IMF's power-law index and a galaxy's SFR. The IGIMF also results in a relation between the galaxy's SFR and the mass of its most massive star such that, if there were no binaries, galaxies with SFR $<10^{-4}$ M$_\odot$/yr should host no Type II supernova events. In addition, a specific list of initial stellar masses can be useful in numerical simulations of stellar systems. For the first time, we show optimally-sampled galaxy-wide IMFs (OSGIMF) which mimics the IGIMF with an additional serrated feature. Finally, A Python module, GalIMF, is provided allowing the calculation of the IGIMF and OSGIMF in dependence on the galaxy-wide SFR and metallicity.Comment: 15 pages, 15 figures, A&A, in press; paper remains unchanged (version1 equals version2); the GalIMF module is downloadable at githu

The possible role of stellar mergers for the formation of multiple stellar populations in globular clusters

Many possible scenarios for the formation of multiple stellar populations (MSP) in globular clusters (GCs) have been discussed so far, including the involvement of asymptotic giant branch stars, fast rotating main sequence stars, very massive main sequence stars and mass-transferring massive binaries based on stellar evolution modelling. But self-consistent, dynamical simulations of very young GCs are usually not considered. In this work, we perform direct $N$-body modelling such systems with total masses up to $3.2\times10^5$ M$_\odot$, taking into account the observationally constrained primordial binary properties, and discuss the stellar-mergers driven both by binary stellar evolution and dynamical evolution of GCs. The occurrence of stellar mergers is enhanced significantly in binary-rich clusters such that stars forming from the gas polluted by mergers-driven ejection/winds would appears as MSPs. We thus emphasize that stellar mergers can be an important process that connects MSP formation with star cluster dynamics, and that multiple MSP formation channels can naturally work together. The scenario studied here, also in view of a possible top-heavy IMF, may be particularly relevant for explaining the high mass fraction of MSPs (the mass budget problem) and the absence of MSPs in young and low-mass star clusters

The impact of the metallicity and star formation rate on the time-dependent galaxy-wide stellar initial mass function

The stellar initial mass function (IMF) is commonly assumed to be an invariant probability density distribution function of initial stellar masses being represented by the canonical IMF. As a consequence the galaxy-wide IMF (gwIMF), defined as the sum of the IMFs of all star forming regions, should also be invariant. Recent observational and theoretical results challenge the hypothesis that the gwIMF is invariant. In order to study the possible reasons for this variation we use the IMF determined in resolved star clusters and apply the IGIMF-theory to calculate a grid of gwIMF models for metallicities, -3<[Fe/H]<1, and galaxy-wide star formation rates, $10^{-5}$<SFR<$10^{5}\,\mathrm{M_{\odot}/yr}$. For a galaxy with metallicy [Fe/H]$1\,M_\odot$/yr, which is a common condition in the early Universe, we find that the gwIMF is top-heavy (more massive stars), when compared to the canonical IMF. For a SFR $< 1\,\mathrm{M_{\odot}/yr}$ the gwIMF becomes top-light regardless of the metallicity. For metallicities $\mathrm{[Fe/H]} > 0$ the gwIMF can become bottom-heavy regardless of the SFR. The IGIMF models predict that massive elliptical galaxies should have formed with a gwIMF that is top-heavy within the first few hundred Myr of the galaxy's life and that it evolves into a bottom-heavy gwIMF in the metal-enriched galactic center. We study the SFR$-$H$\alpha$ relation, its dependency on metallicity and the SFR, the correction factors to the Kennicutt SFR$_{\rm K}-$H$\alpha$ relation, and provide new fitting functions Late-type dwarf galaxies show significantly higher SFRs with respect to Kennicutt SFRs, while star forming massive galaxies have significantly lower SFRs than hitherto thought. This has implications for the gas-consumption time scales and for the main sequence of galaxies. The Leo P and ultra-faint dwarf galaxies are discussed explicitly. [abridged]Comment: Astronomy and Astrophysics (A&A) in press. 15 pages, 8 figure

Gaia DR2 view of the Lupus V-VI clouds: the candidate diskless young stellar objects are mainly background contaminants

Extensive surveys of star-forming regions with Spitzer have revealed populations of disk-bearing young stellar objects. These have provided crucial constraints, such as the timescale of dispersal of protoplanetary disks, obtained by carefully combining infrared data with spectroscopic or X-ray data. While observations in various regions agree with the general trend of decreasing disk fraction with age, the Lupus V and VI regions appeared to have been at odds, having an extremely low disk fraction. Here we show, using the recent Gaia data release 2 (DR2), that these extremely low disk fractions are actually due to a very high contamination by background giants. Out of the 83 candidate young stellar objects (YSOs) in these clouds observed by Gaia, only five have distances of 150 pc, similar to YSOs in the other Lupus clouds, and have similar proper motions to other members in this star-forming complex. Of these five targets, four have optically thick (Class II) disks. On the one hand, this result resolves the conundrum of the puzzling low disk fraction in these clouds, while, on the other hand, it further clarifies the need to confirm the Spitzer selected diskless population with other tracers, especially in regions at low galactic latitude like Lupus V and VI. The use of Gaia astrometry is now an independent and reliable way to further assess the membership of candidate YSOs in these, and potentially other, star-forming regions.Comment: Accepted for publication on Astronomy&Astrophysics Letter

On the degree of stochastic asymmetry in the tidal tails of star clusters

Context: Tidal tails of star clusters are commonly understood to be populated symmetrically. Recently, the analysis of Gaia data revealed large asymmetries between the leading and trailing tidal tail arms of the four open star clusters Hyades, Praesepe, Coma Berenices and NGC 752. Aims: As the evaporation of stars from star clusters into the tidal tails is a stochastic process, the degree of stochastic asymmetry is quantified in this work. Methods: For each star cluster 1000 configurations of test particles are integrated in the combined potential of a Plummer sphere and the Galactic tidal field over the life time of the particular star cluster. For each of the four star clusters the distribution function of the stochastic asymmetry is determined and compared with the observed asymmetry. Results: The probabilities for a stochastic origin of the observed asymmetry of the four star clusters are: Praesepe ~1.7 sigma, Coma Berenices ~2.4 sigma, Hyades ~6.7 sigma, NGC 752 ~1.6 sigma. Conclusions: In the case of Praesepe, Coma Berenices and NGC 752 the observed asymmetry can be interpreted as a stochastic evaporation event. However, for the formation of the asymmetric tidal tails of the Hyades additional dynamical processes beyond a pure statistical evaporation effect are required.Comment: accepted for publication by A&

Gaia DR2 view of the Lupus V-VI clouds: The candidate diskless young stellar objects are mainly background contaminants

Extensive surveys of star-forming regions with Spitzer have revealed populations of disk-bearing young stellar objects. These have provided crucial constraints, such as the timescale of dispersal of protoplanetary disks, obtained by carefully combining infrared data with spectroscopic or X-ray data. While observations in various regions agree with the general trend of decreasing disk fraction with age, the Lupus V and VI regions appeared to have been at odds, having an extremely low disk fraction. Here we show, using the recent Gaia data release 2 (DR2), that these extremely low disk fractions are actually due to a very high contamination by background giants. Out of the 83 candidate young stellar objects (YSOs) in these clouds observed by Gaia, only five have distances of ∼150 pc, similar to YSOs in the other Lupus clouds, and have similar proper motions to other members in this star-forming complex. Of these five targets, four have optically thick (Class II) disks. On the one hand, this result resolves the conundrum of the puzzling low disk fraction in these clouds, while, on the other hand, it further clarifies the need to confirm the Spitzer selected diskless population with other tracers, especially in regions at low galactic latitude like Lupus V and VI. The use of Gaia astrometry is now an independent and reliable way to further assess the membership of candidate YSOs in these, and potentially other, star-forming regions

Scientific Text Analysis with Robots applied to observatory proposals

To test the potential disruptive effect of Artificial Intelligence (AI) transformers (e.g., ChatGPT) and their associated Large Language Models on the time allocation process, both in proposal reviewing and grading, an experiment has been set-up at ESO for the P112 Call for Proposals. The experiment aims at raising awareness in the ESO community and build valuable knowledge by identifying what future steps ESO and other observatories might need to take to stay up to date with current technologies. We present here the results of the experiment, which may further be used to inform decision-makers regarding the use of AI in the proposal review process. We find that the ChatGPT-adjusted proposals tend to receive lower grades compared to the original proposals. Moreover, ChatGPT 3.5 can generally not be trusted in providing correct scientific references, while the most recent version makes a better, but far from perfect, job. We also studied how ChatGPT deals with assessing proposals. It does an apparent remarkable job at providing a summary of ESO proposals, although it doesn\u27t do so good to identify weaknesses. When looking at how it evaluates proposals, however, it appears that ChatGPT systematically gives a higher mark than humans, and tends to prefer proposals written by itself.Proceedings for the SPIE Astronomical Telescopes and Instrumentation 202

Revisiting the globular clusters of NGC 1052-DF2

The ultra-diffuse galaxy (UDG) NGC 1052-DF2 has captured the interest of astronomers ever since the low velocity dispersion measured from ten globular clusters (GCs) suggested a low dark matter fraction. Also, its GC system was found to be unusually bright, with a GC luminosity function peak at least one magnitude brighter than expected for a galaxy at a distance of 20 Mpc. In this work we present an updated view of the GC system of NGC 1052-DF2. We analysed archival MUSE data of NGC 1052-DF2 to confirm the membership of four additional GCs based on their radial velocities, thereby raising the number of spectroscopically confirmed GCs to 16. We measured the ages and metallicities of 11 individual GCs, finding them to be old (> 9 Gyr) and with a range of metallicities from [M/H] = - 0.7 to - 1.8 dex. The majority of GCs are found to be more metal-poor than the host galaxy, with some metal-rich GCs sharing the metallicity of the host ([M/H] = - 1.09- 0.07+0.09 dex). The host galaxy shows a flat age and metallicity gradient out to 1 Re. Using a distance measurement based on the internal GC velocity dispersions (D = 16.2 Mpc), we derived photometric GC masses and find that the peak of the GC mass function compares well with that of the Milky Way. From updated GC velocities, we estimated the GC system velocity dispersion of NGC 1052-DF2 with a simple kinematic model and find σGCS = 14.86- 2.83+3.89 km s- 1. However, this value is reduced to σGCS = 8.63- 2.14+2.88 km s- 1 when the GC that has the highest relative velocity based on a low S/N spectrum is considered an interloper. We discuss the possible origin of NGC 1502-DF2, taking the lower distance, spread in GC metallicities, flat stellar population profiles, and dynamical mass estimate into consideration