The galaxy-wide stellar initial mass function (gwIMF) of a galaxy in dependence of its metallicity and star formation rate (SFR) can be calculated by the integrated galactic IMF (IGIMF) theory. Lacchin et al. (2019) apply the IGIMF theory for the first time to study the chemical evolution of the ultra-faint dwarf (UFD) satellite galaxies and failed to reproduce the data. Here, we find that the IGIMF theory is naturally consistent with the data. We apply the time-evolving gwIMF calculated at each timestep. The number of type Ia supernova explosions per unit stellar mass formed is renormalised according to the gwIMF. The chemical evolution of Bo\"otes I, one of the best observed UFD, is calculated. Our calculation suggests a mildly bottom-light and top-light gwIMF for Bo\"otes I, and that this UFD has the same gas-consumption timescale as other dwarfs but was quenched about 0.1 Gyr after formation, being consistent with independent estimations and similar to Dragonfly 44. The recovered best fitting input parameters in this work are not covered in the work of Lacchin et al. (2019), creating the discrepancy between our conclusions. In addition, a detailed discussion of uncertainties is presented addressing how the results of chemical evolution models depend on applied assumptions. This study demonstrates the power of the IGIMF theory in understanding the star-formation in extreme environments and shows that UDFs are a promising pathway to constrain the variation of the low-mass stellar IMF.

Recent observational and theoretical studies indicate that the stellar initial mass function (IMF) varies systematically with the environment (star formation rate - SFR, metallicity). Although the exact dependence of the IMF on those properties is likely to change with improving observational constraints, the reported trend in the shape of the IMF appears robust. We present the first study aiming to evaluate the effect of the IMF variations on the measured cosmic SFR density (SFRD) as a function of metallicity and redshift, $f_{\rm SFR}$(Z,z). We also study the expected number and metallicity of white dwarf, neutron star and black hole progenitors under different IMF assumptions. Applying the empirically driven IMF variations described by the integrated galactic IMF (IGIMF) theory, we correct $f_{\rm SFR}$(Z,z) obtained by Chruslinska & Nelemans (2019) and find lower SFRD at high redshifts as well as a higher fraction of metal-poor stars being formed. In the local Universe, our calculation applying the IGIMF theory suggests more white dwarf and neutron star progenitors in comparison with the universal IMF scenario, while the number of black hole progenitors remains unaffected.

Several recent studies have shown that the Vela OB2 region hosts a complex constellation of sub-populations with ages in the range 10 to 50 Myr. Such populations might represent the best example of the outcome of clustered star formation in Giant Molecular clouds (GMC). We use Gaia DR2 data over an area of 40 deg radius around the open cluster Collinder 135 to extend the study of the stellar populations of the Vela OB2 region over an area of several hundreds of parsecs on sky. Detailed clustering algorithms combined with the exquisite astrometric quality of the GAIA catalogue allow us to detect a new cluster named BBJ 1 that shows the same age as NGC 2547 (30 to 35 Myr), but located at a distance of 260 pc from it. Deeper investigation of the region via clustering in 5D parameter space and in the colour-magnitude diagram allows us to detect a filamentary structure of stars that bridges the two clusters. Given the extent in space of such structure (260 pc) and the young age (~35 Myr), we exclude that such population originates by the same mechanism responsible to create tidal streams around older clusters. Even if we miss a complete picture of the 3D motion of the studied stellar structure because of the lack of accurate radial velocity measurements, we propose that such structure represent the detection of a 35 Myr old outcome of a mechanism of filamentary star formation in a GMC.

The alpha element to iron peak element ratio, for example [Mg/Fe], is a commonly applied indicator of the galaxy star formation timescale (SFT) since the two groups of elements are mainly produced by different types of supernovae that explode over different timescales. However, it is insufficient to consider only [Mg/Fe] when estimating the SFT. The [Mg/Fe] yield of a stellar population depends on its metallicity. Therefore, it is possible for galaxies with different SFTs and at the same time different total metallicity to have the same [Mg/Fe]. This effect has not been properly taken into consideration in previous studies. In this study, we assume the galaxy-wide stellar initial mass function (gwIMF) to be canonical and invariant. We demonstrate that our computation code reproduces the SFT estimations of previous studies where only the [Mg/Fe] observational constraint is applied. We then demonstrate that once both metallicity and [Mg/Fe] observations are considered, a more severe "downsizing relation" is required. This means that either low-mass ellipticals have longer SFTs (> 4 Gyr for galaxies with mass below $10^{10}$ M$_\odot$) or massive ellipticals have shorter SFTs ($\approx 200$ Myr for galaxies more massive than $10^{11}$ M$_\odot$) than previously thought. This modification increases the difficulty in reconciling such SFTs with other observational constraints. We show that applying different stellar yield modifications does not relieve this formation timescale problem. The quite unrealistically short SFT required by [Mg/Fe] and total metallicity would be prolonged if a variable stellar gwIMF were assumed. Since a systematically varying gwIMF has been suggested by various observations this could present a natural solution to this problem.

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 stellar initial mass function (IMF) is the key to understanding the matter cycle in the Universe. Edwin Salpeter's paper of 1955 founded this research field. Evidence today, however, challenges the initial mass function as an invariant probability distribution function.

Young star clusters (YSCs) with resolved stellar populations are well suited for studying star-cluster formation. In most cases, the (pre-main-sequence) stellar populations found in the YSCs are coeval with an intrinsic age spread of up to 1Myr. Such observations can be understood as the YSCs having formed in one burst, which star formation was truncated by stellar feedback. The recent discovery that the colour-magnitude diagram of the Orion Nebula Clusters (ONC) contains three well defined age-separated populations appears to shatter this model. The implication is that the ONC formed in three bursts, with star formation still on-going in the last burst. We present new observational results focusing on the three populations in the ONC using OmegaCAM photometry and Gaia DR2 measurements. We also describe a theoretical model which may explain these observations by an interplay between stellar feedback and cluster dynamics.

We report the discovery of the oldest stellar substructure in the Orion star forming region (OSFR), the Orion relic filament. The relic filament is physically associated with the OSFR as demonstrated by Gaia DR2 photometry and astrometry, as well as targeted radial velocity follow-up observations of a bright sub-sample of proper-motion selected candidate members. Gaia DR2 parallaxes place the Orion relic filament in the more distant part of the OSFR, approx. 430pc from the Sun. Given its age, velocity dispersion, spatial extent, and shape, it is not possible for the Orion relic filament to have formed as a single stellar cluster, even taking into account residual gas expulsion. The relic filament is also too young to be a tidal stream, since Galactic tides act on much longer time scales of order 100 Myr. It therefore appears likely that the structure formed from a molecular cloud filament similar to Orion A in the OSFR and retained its morphology despite decoupling from its natal gas. Hence, the Orion relic filament bears witness to the short-lived evolutionary phase between gas removal and dispersion due to shears and tides, and provides crucial new insights into how stars are formed in molecular clouds.

Growing evidence in recent years suggests a systematic variation of the stellar initial mass function (IMF), being top-heavy for starburst galaxies and possibly bottom-heavy for massive ellipticals. Galaxy chemical evolution simulations adopting an invariant canonical IMF face difficulty in simultaneously reproducing the metallicity and $\alpha$-enhancement of the massive elliptical galaxies. Applying a variable IMF that changes with time is a promising solution, however, it is non-trivial to couple a variable IMF theory with the existing galaxy evolution codes. Here we present the first open source simulation code which recalculates the galaxy-wide IMF at each time step according to the Integrated-Galactic-IMF (IGIMF) theory where the galaxy-wide IMF depends on the galactic star formation rate and metallicity. The resulting galaxy-wide IMF and metal abundance evolve with time. With this pilot work, we explore the effect of the IGIMF theory on galaxy chemical evolution in comparison with an invariant IMF.

The high-quality OmegaCAM photometry of the 3x3 deg around the Orion Nebula Cluster (ONC) in r, and i filters by Beccari et al.(2017) revealed three well-separated pre-main sequences in the color-magnitude diagram (CMD). The objects belonging to the individual sequences are concentrated towards the center of the ONC. The authors concluded that there are two competitive scenarios: a population of unresolved binaries and triples with an exotic mass ratio distribution, or three stellar populations with different ages. We use Gaia DR2 in combination with the photometric OmegaCAM catalog to test and confirm the presence of the putative three stellar populations. We also study multiple stellar systems in the ONC for the first time using Gaia DR2. We confirm that the second and third sequence members are more centrally concentrated towards the center of the ONC. In addition we find an indication that the parallax and proper motion distributions are different among the members of the stellar sequences. The age difference among stellar populations is estimated to be 1-2 Myr. We use Gaia measurements to identify and remove as many unresolved multiple system candidates as possible. Nevertheless we are still able to recover two well-separated sequences with evidence for the third one, supporting the existence of the three stellar populations. We were able to identify a substantial number of wide binary objects (separation between 1000-3000 au). This challenges previously inferred values that suggested no wide binary stars exist in the ONC. Our inferred wide-binary fraction is approx 5%. We confirm the three populations correspond to three separated episodes of star formation. Based on this result, we conclude that star formation is not happening in a single burst in this region. (abridged)

The so-called ultra-diffuse galaxy NGC~1052-DF2 was announced to be a galaxy lacking dark matter based on a spectroscopic study of its constituent globular clusters. Here we present the first spectroscopic analysis of the stellar body of this galaxy using the MUSE integral-field spectrograph at the (ESO) Very Large Telescope. The MUSE datacube simultaneously provides DF2's stellar velocity field and systemic velocities for seven globular clusters (GCs). We further discovered three planetary nebulae (PNe) that are likely part of this galaxy. While five of the clusters had velocities measured in the literature, we were able to confirm the membership of two more candidates through precise radial velocity measurements, which increases the measured specific frequency of GCs in DF2. The mean velocity of the diffuse stellar body, 1792.9$^{-1.8}_{+1.4}$~\kms, is consistent with the mean globular cluster velocity. We detect a weak but significant velocity gradient within the stellar body, with a kinematic axis close to the photometric major-axis, making it a prolate-like rotator. We estimate a velocity dispersion from the clusters and PNe of $\sigma_{\mathrm{int}}=10.6^{+3.9}_{-2.3}$~\kms. The velocity dispersion $\sigma_{\rm{DF2}\star}$(\re) for the stellar body within one effective radius is $10.8^{-4.0}_{+3.2}$~\kms. Considering various sources of systemic uncertainties this central value varies between 5 and 13~\kms, and we conservatively report a 95\% confidence upper limit to the dispersion within one \re\ of 21~\kms. We provide updated mass estimates based on these dispersions corresponding to the different distances to NGC~1052-DF2 that have been reported in the recent literature.

NGC 1052-DF2, an ultra diffuse galaxy (UDG), has been the subject of intense debate. Its alleged absence of dark matter, and the brightness and number excess of its globular clusters (GCs) at an initially assumed distance of 20Mpc, suggested a new formation channel for UDGs. We present the first systematic spectroscopic analysis of both the stellar body and the GCs (six of which were previously known, and one newly confirmed member) of this galaxy using MUSE@VLT. Even though NGC 1052-DF2 does not show any spatially extended emission lines we report the discovery of three planetary nebulae (PNe). We conduct full spectral fitting on the UDG and the stacked spectra of all GCs. The UDG's stellar population is old, 8.9$\pm$1.5 Gyr, metal-poor, with [M/H] = $-$1.07$\pm$0.12 with little or no $\alpha$-enrichment. The stacked spectrum of all GCs indicates a similar age of 8.9$\pm$1.8 Gyr, but lower metallicity, with [M/H] = $-$1.63$\pm$0.09, and similarly low $\alpha$-enrichment. There is no evidence for a variation of age and metallicity in the GC population with the available spectra. The significantly more metal-rich stellar body with respect to its associated GCs, the age of the population, its metallicity and alpha enrichment, are all in line with other dwarf galaxies. NGC 1052-DF2 thus falls on the same empirical mass-metallicity relation as other dwarfs, for the full distance range assumed in the literature. We find that both debated distance estimates (13 and 20 Mpc) are similarly likely, given the three discovered PNe.

(abridged) The stellar initial mass function (IMF) can be conveniently represented as a canonical two-part power law function and is largely invariant for star formation regions evident in the Local Group of galaxies. The IMF is a hilfskonstrukt. It is a mathematical formulation of an idealised population of stars formed together in one star formation event. The nature of the IMF (is it a probability density or an optimal sampling distribution function?) is raised. Binary stars, if unresolved, have a very significant influence at low stellar masses. Especially important is to take care of the changing binary fraction as a result of stellar-dynamical evolution of the embedded clusters which spawn the field populations of galaxies, given that the binary fraction at birth is very high and independent of primary-star mass. The high multiplicity fraction amongst massive stars leads to a substantial fraction of these being ejected out of their birth clusters and to massive stars merging. This explains the top-lightness of the IMF in star clusters in M31. In close binaries also the masses of the components can be changed due to mass transfer. A large amount of evidence points to the IMF becoming top-heavy with decreasing metallicity and above a star-formation-rate density of about 0.1 Msun/(pc^3 yr) of the cluster-forming cloud core. This is also indicated through the observed supernova rates in star-bursting galaxies. At the same time, the IMF may be bottom light at low metallicity and bottom-heavy at high metallicity, possibly accounting for the results on elliptical galaxies and ultra-faint dwarf galaxies, respectively.

Recently, three stellar sequences separated in age by about 1 Myr were discovered in the Orion Nebula Cluster (ONC; Beccari et al. 2017). Kroupa et al. (2018) suggest that such small dense subpopulations eject all their OB stars via the decay of unstable few-body systems such that the gas can recombine and form new stars. This explains the multisequence phenomenon without introducing an extra mechanism into star formation theory. In this work, we apply the recently updated primordial binary distribution model (Belloni et al. 2017; implemented here in a new version of MCLUSTER) and perform a large set of direct N-body simulations to investigate the feasibility of this dynamical scenario. Our results suggest that if 3-4 OB stars in the ONC formed primordially mass-segregated in the cluster centre with a maximum separation of about 0.003 pc, all OB stars have a high chance (50-70 per cent) to escape from the centre and do not come back within 1 Myr and the dynamical ejection scenario is a viable channel to form short-age-interval multipopulation sequences as observed in the ONC. This is also consistent with self-regulated star formation.

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]$<0$ and SFR$\,> 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]

Using Gaia DR2 data, combined with OmegaCAM ground-based optical photometry from the AD-HOC survey, and detailed Radial Velocity measurements from ESO-Gaia, we analyse in detail a 10x5 deg region around the Wolf-Rayet star $\gamma^2$ Vel, including the previously known clusters Gamma Vel and NGC2547. Using clustering analysis that considers positions, proper motions and parallax, we discover 6 clusters or associations -- 4 of which appear new. Analysis of the colour-magnitude diagram for these clusters show that 4 of them formed coevally from the same molecular clouds 10 Myr ago, while NGC 2547 formed together with a newly discovered cluster 30 Myr ago. This study shows the incredible wealth of data provided by Gaia for the study of young stellar clusters.

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.

We present optical photometry of Hubble Space Telescope (HST) ACS/WFC data of the resolved stellar populations in the outer disc of the dwarf irregular galaxy DDO 154. The photometry reveals that young main sequence stars are almost absent from the outermost HI disc. Instead, most are clustered near the main stellar component of the galaxy. We constrain the stellar initial mass function (IMF) by comparing the luminosity function of the main sequence stars to simulated stellar populations assuming a constant star formation rate over the dynamical timescale. The best-fitting IMF is deficient in high mass stars compared to a canonical Kroupa IMF, with a best-fit slope $\alpha = -2.45$ and upper mass limit $M_U = 16\ M_{\odot}$. This top-light IMF is consistent with predictions of the Integrated Galaxy-wide IMF theory. Combining the HST images with HI data from The HI Nearby Galaxy Survey Treasury (THINGS) we determine the star formation law (SFL) in the outer disc. The fit has a power law exponent $N = 2.92 \pm0.22$ and zero point $A=4.47 \pm 0.65 \times 10^{-7} \ M_{\odot} \ \text{yr}^{-1} \ \text{kpc}^{-2}$. This is depressed compared to the Kennicutt-Schmidt Star Formation Law, but consistent with weak star formation observed in diffuse HI environments. Extrapolating the SFL over the outer disc implies that there could be significant star formation occurring that is not detectable in H$\alpha$. Last, we determine the Toomre stability parameter $Q$ of the outer disc of DDO 154 using the THINGS HI rotation curve and velocity dispersion map. 72% of the HI in our field has $Q\leq 4$ and this incorporates 96% of the observed MS stars. Hence 28% of the HI in the field is largely dormant.

(abridged) A scenario for the formation of multiple co-eval populations separated in age by about 1 Myr in very young clusters (VYCs, ages less than 10 Myr) and with masses in the range 600-20000 Msun is outlined. It rests upon a converging inflow of molecular gas building up a first population of pre-main sequence stars. The associated just-formed O stars ionise the inflow and suppress star formation in the embedded cluster. However, they typically eject each other out of the embedded cluster within 10^6 yr, that is before the molecular cloud filament can be ionised entirely. The inflow of molecular gas can then resume forming a second population. This sequence of events can be repeated multiply. This model is applied to the Orion Nebula Cluster (ONC), in which three well-separated pre-main sequences in the color-magnitude diagram of the cluster have recently been discovered. The mass-inflow history is constrained using this model and the number of OB stars ejected from each population are estimated for verification using Gaia data. As a consequence of the proposed model, the three runaway O star systems, AE Aur, mu Col and iota Ori, are considered as significant observational evidence for stellar-dynamical ejections of massive stars from the oldest population in the ONC.

The stellar initial mas function (IMF) has been described as being invariant, bottom heavy or top-heavy in extremely dense star burst conditions. To provide usable observable diagnostic we calculate redshift dependent spectral energy distributions of stellar populations in extreme star burst clusters which are likely to have been the precursors of present day massive globular clusters (GCs) and of ultra compact dwarf galaxies (UCDs). The retention fraction of stellar remnants is taken into account to asses the mass to light ratios of the ageing star-burst. Their redshift dependent photometric properties are calculated as predictions for James Webb Space Telescope (JWST) observations. While the present day GCs and UCDs are largely degenerate concerning bottom-heavy or top-heavy IMFs, a metallicity- and density-dependent top-heavy IMF implies the most massive UCDs, at ages <100 Myr, to appear as objects with quasar-like luminosities with a 0.1-10% variability on a monthly time scale due to core collapse supernovae.

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.

Recent observations have revealed that at least several old globular clusters (GCs) in the Galaxy have discrete distributions of stars along the Mg-Al anti-correlation. In order to discuss this recent observation, we construct a new one-zone GC formation model in which the maximum stellar mass (m_max) in the initial mass function (IMF) of stars in a forming GC depends on the star formation rate (SFR), as deduced from independent observations. We investigate the star formation histories of forming GCs. The principal results are as follows. About 30 Myr after the formation of the first generation (1G) of stars within a particular GC, new stars can be formed from ejecta from asymptotic giant branch (AGB) stars of 1G. However, the formation of this second generation (2G) of stars can last only for [10-20] Myr, because the most massive SNe of 2G expel all of the remaining gas. The third generation (3G) of stars are then formed from AGB ejecta ~ 30$ Myr after the truncation of 2G star formation. This cycle of star formation followed by its truncation by SNe can continue until all AGB ejecta is removed from the GC by some physical process. Thus, it is inevitable that GCs have discrete multiple stellar populations in the [Mg/Fe]-[Al/Fe] diagram. Our model predicts that low-mass GCs are unlikely to have discrete multiple stellar populations, and young massive clusters may not have massive OB stars owing to low m_max ([20-30] M_sun) during the secondary star formation.