Search

The Spitzer instability leads to the formation of a black hole subsystem (BHSub) at the centre of a star cluster providing energy to luminous stars (LSs) and increasing their rate of evaporation. When the self-depletion time of the BHSub exceeds the evaporation time of the LSs, a dark star cluster (DSC) will appear. Using the nbody7 code, we performed a comprehensive set of direct N -body simulations over a wide range of initial conditions to study the pure effect of the top-heaviness of the initial mass function (IMF) on the formation of the DSC phase. In the Galactic tidal field, top-heavy IMFs lead to the fast evaporation of LSs and the formation of DSCs. Therefore, DSCs can be present even in the outer region of the Milky Way (MW). To successfully transition to the DSC phase, the MW Globular Clusters (GCs) must possess an initial black hole (BH) mass fraction of |$\widetilde {\mathit {M}}_\mathrm{BH}(0)\gt 0.05$|⁠. For star clusters with |$\widetilde {\mathit {M}}_\mathrm{BH}(0)\gt 0.08$|⁠ , the DSC phase will be created for any given initial density of the cluster and Galactocentric distance. The duration of the cluster's lifetime spent in the DSC phase shows a negative (positive) correlation with the initial density, and Galactocentric distance of the star cluster if |$\widetilde {\mathit {M}}_\mathrm{BH}(0)\le 0.12$| (⁠|$\widetilde {\mathit {M}}_\mathrm{BH}(0)\ge 0.15$|⁠). Considering the canonical IMF, it is unlikely for any MW GCs to enter the DSC phase. We discuss the BH retention fraction in view of the observed properties of the GCs of the MW. [ABSTRACT FROM AUTHOR]

We present a comprehensive investigation on the mass function (MF) of a snake-like stellar structure in the solar neighbourhood, building on our previous discovery. To ensure the reliability of the data, we reselect the member stars of the Stellar 'Snake' in the latest Gaia Data Release 3 using the same approach as the initial series of articles. We also precisely measure the physical parameters of the clusters within the Stellar Snake. In light of the high completeness of the member stars in the cluster regions, we develop a simulated model colour–magnitude diagram-based inference method to derive the mass function, binary fraction, and mass-ratio distribution of the clusters in the Stellar Snake. Notably, despite their similar ages and metallicity, we discover systematic variations in the MFs along the elongation direction of the Snake in the mass range of 0.5 to 2.0 M⊙. The 'head' of the Snake conforms to a canonical initial mass function with a power-law slope of α ∼ −2.3. Extending towards the 'tail', the MF becomes more top-light, indicating a deficiency of massive stars within these clusters. This result provides evidence for the delayed formation of massive stars in the clusters. Such clues give support to the hypothesis that the Stellar Snake constitutes as a hierarchically primordial structure. [ABSTRACT FROM AUTHOR]

A dark star cluster (DSC) is a system in which the cluster potential is dominated by stellar remnants, such as black holes and neutron stars having larger masses than the long-lived low-mass stars. Due to mass segregation, these remnants are located in the central region of the cluster and form a dark core. We expect that at a few kpc from the Galactic Centre, the efficient evaporation of the lower-mass stars caused by the strong tidal force exposes the dark core, because the dynamical properties of the DSC are dominated by the remnants. Due to the invisibility of the remnants, finding a DSC by observation is challenging. In this project, we use N -body simulations to obtain models of DSCs and try to discern observables that signify a DSC. We consider four observables: the mass spectrum, the observational mass density profile, the observational velocity dispersion profile and the mass segregation. The models show that a DSC typically exhibits two distinct characteristics: for a given mass in stars and a given half-light radius, the expected velocity dispersion is underestimated when only visible stars are considered, and there is a lack of measurable mass segregation among the stars. These properties can be helpful for finding DSCs in observational data, such as the Gaia catalogue. [ABSTRACT FROM AUTHOR]

We present the observational mass–radius (M–R) relation for a sample of 47 magnetized white dwarfs (WDs) with the magnetic field strength (B) ranging from 1 to 773 MG, identified from the SDSS data release 7 (DR7). We derive their effective temperature, surface gravity (log g), luminosity, radius, and mass. While atmospheric parameters are derived using a Virtual Observatory Spectral Energy Distribution Analyzer (VOSA), the mass is derived using their location in the HR diagram in comparison with the evolutionary tracks of different masses. We implement this mass measurement instead of a more traditional method of deriving masses from log g , which is unreliable as is based on SED and generates errors from other physical parameters involved. The main disadvantage of this method is that we need to assume a core composition of WDs. As it is complicated to identify the exact composition of these WDs from low-resolution spectra, we use tracks for the masses 0.2 to 0.4 M⊙ assuming a He-core, 0.5 to 1.0 M⊙ assuming CO core, and above M⊙ assuming O–Ne–Mg core. We compare the observed M–R relation with those predicted by the finite temperature model by considering different B , which are well in agreement considering their relatively low-surface fields, ≲109 G. Currently, there is no direct observational detection of magnetized WDs with B > 109 G. We propose that our model can be further extrapolated to higher B , which may indicate the existence of super-Chandrasekhar mass (M > 1.4 M⊙) WDs at higher B. [ABSTRACT FROM AUTHOR]

The properties of the progenitors of gamma-ray bursts and of their environment are encoded in their luminosity function and cosmic formation rate. They are usually recovered from a flux-limited sample based on Lynden-Bell's |$c^{-}$| method. However, this method is based on the assumption that the luminosity is independent of the redshift. Observationally, if correlated, people use non-parametric |$\tau$| statistical method to remove this correlation through the transformation, |$L^{\prime }=L/g(z)$|⁠ , where z is the burst redshift, and |$g(z)=(1+z)^{k}$| parametrizes the underlying luminosity evolution. However, the application of this method to different observations could result in very different luminosity functions. By the means of Monte Carlo simulation, in this paper, we demonstrate that the origin of an observed correlation, measured by the |$\tau$| statistical method, is a complex combination of multiple factors when the underlying data are correlated. Thus, in this case, it is difficult to unbiasedly reconstruct the underlying population distribution from a truncated sample, unless the detailed information of the intrinsic correlation is accurately known in advance. In addition, we argue that an intrinsic correlation between luminosity function and formation rate is unlikely eliminated by a misconfigured transformation, and the |$g(z)$|⁠ , derived from a truncated sample with the |$\tau$| statistical method, does not necessarily represent its underlying luminosity evolution. [ABSTRACT FROM AUTHOR]

If primordial black holes (PBHs) constitute the dark matter (DM), stars forming in dark-matter dominated environments with low velocity dispersions, such as ultra-faint dwarf galaxies, may capture a black hole at birth. The capture probability is non-negligible for PBHs of masses around 1020 g, and increases with stellar mass. Moreover, infected stars are turned into virtually invisible black holes on cosmologically short time-scales. Hence, the number of observed massive main-sequence stars in ultra-faint dwarfs should be suppressed if the DM was made of asteroid-mass PBHs. This would impact the measured mass distribution of stars, making it top-light (i.e. depleted in the high-mass range). Using simulated data that mimic the present-day observational power of telescopes, we show that already existing measurements of the mass function of stars in local ultra-faint dwarfs could be used to constrain the fraction of DM composed of PBHs in the – currently unconstrained – mass range of 1019–1021 g. [ABSTRACT FROM AUTHOR]

We have studied a star-forming complex S193 using near-infrared (NIR) observations and other archival data covering optical to radio wavelengths. We identified stellar clusters in the complex using the NIR photometric data and estimated the membership and distance of the clusters. Using the mid-infrared (MIR) and far-infrared (FIR) images, the distribution of the dust emission around H  ii regions is traced in the complex. The Herschel column density and temperature maps analysis reveal 16 cold dust clumps in the complex. The H α image and 1.4 GHz radio continuum emission map are employed to study the ionized gas distribution and infer the spectral type and the dynamical age of each H  ii region/ionized clump in the complex. The 12CO(J  = 3 − 2) and 13CO(J  = 1 − 0) molecular line data hint at the presence of two velocity components around [−43, −46] and [−47, −50] km s−1, and their spatial distribution reveals two overlapping zones toward the complex. By investigating the immediate surroundings of the central cluster [BDS2003]57 and the pressure calculations, we suggest that the feedback from the massive stars seems responsible for the observed velocity gradient and might have triggered the formation of the central cluster [BDS2003]57. [ABSTRACT FROM AUTHOR]

During star formation, dust plays the crucial role of coupling gas to stellar radiation fields, allowing radiation feedback to influence gas fragmentation and thus the stellar initial mass function (IMF). Variations in dust abundance therefore provide a potential avenue by which variation in galaxy metallicity might affect the IMF. In this paper, we present a series of radiation-magnetohydrodynamic simulations in which we vary the metallicity and thus the dust abundance from 1 per cent of solar to 3× solar, spanning the range from the lowest metallicity dwarfs to the most metal-rich early-type galaxies (ETGs) found in the local Universe. We design the simulations to keep all dimensionless parameters constant so that the interaction between feedback and star-forming environments of varying surface density and metallicity is the only factor capable of breaking the symmetry between the simulations and modifying the IMF, allowing us to isolate and understand the effects of each environmental parameter cleanly. We find that shifts in the IMF with varying metallicity at a fixed surface density are smaller than the shifts in varying surface density at a fixed surface metallicity. We also find that metallicity-induced IMF variations are too small to explain the mass-to-light ratio shifts seen in the ETGs. We therefore conclude that metallicity variations are much less important than variations in surface density in driving changes in the IMF and that the latter rather than the former are most likely responsible for the IMF variations found in ETGs. [ABSTRACT FROM AUTHOR]

The observed radial accelerations of 462 early-type galaxies (ETGs) at their half-mass radii are discussed. They are compared to the baryonic masses of the same galaxies, which are derived from theoretical expectations for their stellar populations and cover a range from |$\approx 10^4 \, {\rm M}_{\odot }$| to |$\approx 10^{11} \, {\rm M}_{\odot }$|⁠. Both quantities are plotted against each other, and it is tested whether they lie (within errors) along theoretical radial acceleration relations (RARs). We choose the Newtonian RAR and two Milgromian, or MONDian RARs. At low radial accelerations (corresponding to low masses), the Newtonian RAR fails without non-baryonic dark matter, but the two MONDian ones may work, provided moderate out-of-equilibrium dynamics in some of the low-mass ETGs. However all three RARs fail at high accelerations (corresponding to high masses) if all ETGs have formed their stellar populations with the canonical stellar initial mass function (IMF). A much better agreement with the observations can however be accomplished, if the theory of the integrated galaxy-wide stellar initial mass functions (IGIMFs) is used instead. This is because the IGIMF-theory predicts the formation of an overabundance of stellar remnants during the lifetime of the massive ETGs. Thus their baryonic masses today are higher than they would be if the ETGs had formed with a canonical IMF. Also the masses of the stellar-mass black holes should be rather high, which would mean that most of them probably formed when the massive ETGs were not as metal-enriched as they are today. The IGIMF-approach confirms downsizing. [ABSTRACT FROM AUTHOR]

When star clusters are formed at low star-formation rates (SFRs), their stellar initial mass function (IMF) can hardly be filled continuously with stars at each mass. This lack holds for massive stars and is verified observationally by the correlation between star-cluster mass and its most massive cluster star. Since galaxy evolution is strongly affected by massive stars, numerical models should account for this lack. Because a filled IMF is mostly applied even when only fractions of massive stars form, here we investigate, by means of 3D chemo-dynamical simulations of isolated dwarf galaxies, how deviations from a standard IMF in star clusters affect the evolution. We compare two different IMF recipes, a filled IMF and one truncated at a maximum mass at which a single complete star forms. Attention is given to energetic and chemical feedback by massive stars. Since their energy release is mass-dependent but steeper than the negative IMF slope, the energetic feedback retains a positive mass dependence, so that a filled IMF regulates star formation (SF) more strongly than truncated IMFs, though only stellar number fractions exist. The higher SFR of the truncated IMF in the simulation leads to more Type II supernovae (SNeII), driving galactic winds. Whether this results from the model-inherent larger SFR is questioned and therefore explored analytically. This shows the expected result for the Lyman continuum, but that the total SNII energy release is equal for both IMF modes, while the power is smaller for the truncated IMF. Reasonably, the different IMFs leave fingerprints in the abundance ratios of massive to intermediate-mass star elements. [ABSTRACT FROM AUTHOR]

Low-mass models of M-dwarfs that undergo the convective kissing instability fluctuate in luminosity and temperature resulting in a gap in the main sequence that is observed in the Gaia data. During this instability, the models have repeated periods of full convection where the material is mixed throughout the model. Stellar evolution models are performed using mesa with varying amounts of convective overshooting and semi-convection. We find that the amplitude and intensity of the instability is reduced with increasing amounts of overshooting but sustained when semi-convection is present. This is reflected in the loops in the evolutionary tracks in the Hertzsprung–Russell diagram. The surface abundances of 1H, 3He, 4He, 12C, 14N, and 16O increase or decrease over time due to the convective boundary, however the relative abundance changes are very small and not likely observable. The mass and magnitude values from the models are assigned to a synthetic population of stars from the mass–magnitude relation to create colour–magnitude diagrams, which reproduce the M-dwarf gap as a large indent into the blueward edge of the main sequence (MS). This is featured in the luminosity function as a small peak and dip. The width of the MS decreases over time along with the difference in width between the MS at masses higher and lower than the instability. The parallel offset and relative angle between the upper and lower parts of the MS also change with time along with the mass–magnitude relation. Potential age-dating methods for single stars and stellar populations are described. [ABSTRACT FROM AUTHOR]

We used machine learning techniques to predict period changes in variable stars using noisy and sparse time-series data, while inferring underlying physics and generalizing predictions about cycle-to-cycle variations. Our focus was on Mira variables, a well-known class of pulsating stars. Preprocessing data from Mira, R Andromedae, U Orionis, and Chi Cygni, obtained from the American Association of Variable Star Observers, we predicted luminosity magnitude uncertainty and classified pulsation states. Employing various classification and regression algorithms, along with feature engineering, we aimed to generalize predictions. We created a generalized data set with collective averaged data points, limiting our analysis to a common time duration. Linear regression models yielded no successful predictions, but decision tree and KNN regressors accurately predicted luminosity magnitude errors, indicative of variation over time. Feature engineering successfully aided regression and classification of pulsating star states. After hyper-parameter tuning using Bayesian neural networks, we achieved a classification accuracy of 0.8 and 0.94 for the KNN classifier, respectively, in classifying pulsation states of Mira variables. The regression model achieved an R2 score of 0.98. Our work provides a foundation for developing tools to analyze various pulsating star variables, including Cepheids, RR Lyrae, and Delta Scuti variables, as well as other astrophysical data. These techniques demonstrate impressive performance with time series data sets. [ABSTRACT FROM AUTHOR]

Many experimental quantities show a power-law distribution p (x) ∝  x −α. In astrophysics, examples are: size distribution of dust grains or luminosity function of galaxies. Such distributions are characterized by the exponent α and by the extremes x min x max where the distribution extends. There are no mathematical tools that derive the three unknowns at the same time. In general, one estimates a set of α corresponding to different guesses of x min x max. Then, the best set of values describing the observed data is selected a posteriori. In this paper, we present a tool that finds contextually the three parameters based on simple assumptions on how the observed values xi populate the unknown range between x min and x max for a given α. Our tool, freely downloadable, finds the best values through a non-linear least-squares fit. We compare our technique with the maximum likelihood estimators for power-law distributions, both truncated and not. Through simulated data, we show for each method the reliability of the computed parameters as a function of the number N of data in the sample. We then apply our method to observed data to derive: (i) the slope of the core mass function in the Perseus star-forming region, finding two power-law distributions: α = 2.576 between |$1.06\, \mathrm{M}_{\odot }$| and |$3.35\, \mathrm{M}_{\odot }$|⁠ , α = 3.39 between |$3.48\, \mathrm{M}_{\odot }$| and |$33.4\, \mathrm{M}_{\odot }$|⁠ ; (ii) the slope of the γ-ray spectrum of the blazar J0011.4+0057, extracted from the Fermi-LAT archive. For the latter case, we derive α = 2.89 between 1484 MeV and 28.7 GeV; then we derive the time-resolved slopes using subsets of 200 photons each. [ABSTRACT FROM AUTHOR]

We analyse the structure of 15 protocluster forming regions in the Milky Way using their 1.3 mm continuum emission maps from the ALMA-IMF large program. The analysis of the cloud structure is performed using the delta-variance spectrum technique. The calculated spectra display a self-similar regime on small scales as well as the presence of a prominent bump on larger scales and whose physical size, L hub, falls in the range of ≈7000–60 000 au. These scales correspond to the sizes of the most compact clumps within the protocluster forming clouds. A significant correlation is found between L hub and the surface density of the free–free emission estimated from the integrated flux of the H41α recombination line |$\left(\Sigma _{\rm H41\alpha }^{\rm free \!-\! free}\right)$| as well as a significant anticorrelation between L hub and the ratio of the 1.3–3 mm continuum emission fluxes |$\left(S_{\rm 1.3 \, mm}^{\rm cloud}/S_{\rm 3 \, mm}^{\rm cloud}\right)$|⁠. Smaller values of |$\left(S_{\rm 1.3 \, mm}^{\rm cloud}/S_{\rm 3 \, mm}^{\rm cloud}\right)$| and larger values of |$\Sigma _{\rm H41\alpha }^{\rm free \!-\! free}$| correspond to more advanced evolutionary stages of the protocluster forming clumps. Hence, our results suggest that the sizes of the densest regions in the clouds are directly linked to their evolutionary stage and to their star formation activity with more evolved clouds having larger protocluster forming clumps. This is an indication that gravity plays a vital role in regulating the size, mass growth, and star formation activity of these clumps with ongoing gas accretion. [ABSTRACT FROM AUTHOR]

Context. The dominant formation channel of brown dwarfs (BDs) is not well constrained yet and a promising way to discriminate between scenarios is to test the environment-dependent efficiency in forming BDs. So far, the outcome of star formation, studied through the initial mass function, has been found to be very similar in all clusters that have been inspected. Aims. We aim to characterize the low-mass (sub)stellar population of the central portion (2.4 pc2) of the ∼2 Myr old cluster NGC 2244 using near-infrared spectroscopy. By studying this cluster, characterized by a low stellar density and numerous OB stars, we aim to explore the effect that OB stars may have on the production of BDs. Methods. We obtained near-infrared HK spectroscopy of 85 faint candidate members of NGC 2244. We derived the spectral type and extinction by comparison with spectral templates. We evaluated cluster membership using three gravity-sensitive spectral indices based on the shape of the H-band. Furthermore, we evaluated the infrared excess from Spitzer of all the candidate members of the cluster. Finally, we estimated the mass of all the candidate members of the cluster and derived the initial mass function, star-to-BD number ratio, and disk fraction. Results. The initial mass function is well represented by a power law (dN/dM ∝ M−α) below 0.4 M⊙, with a slope α = 0.7–1.1 depending on the fitted mass range. We calculated a star-to-BD number ratio of 2.2–2.8. We find the low-mass population of NGC 2244 to be consistent with nearby star-forming regions, although it is at the high end of BD production. We find BDs in NGC 2244 to be, on average, closer to OB stars than to low-mass stars, which could potentially be the first piece of evidence that OB stars affect the formation of BDs. We find a disk fraction of all the members with a spectral type later than K0 of 39 ± 9% which is lower than typical values found in nearby star-forming regions of similar ages. [ABSTRACT FROM AUTHOR]

We present the results of a set of radiation magnetohydrodynamic simulations of turbulent molecular clouds in which we vary the initial strength of the magnetic field within a range (1 ≲ μ ≲ 5) consistent with observations of local giant molecular clouds (GMCs). We find that as we increase the strength of the magnetic field, star formation transitions from unimodal (the baseline case, μ = 5, with a single burst of star formation and Salpeter IMF) to bimodal. This effect is clearest in the most strongly magnetized GMC (μ = 1): a first burst of star formation with duration, intensity, and IMF comparable to the baseline case, is followed by a second star formation episode in which only low-mass stars are formed. Overall, due to the second burst of star formation, the strongly magnetized case results in a longer star formation period and a higher efficiency of star formation. The second burst is produced by gas that is not expelled by radiative feedback, instead remains trapped in the GMC by the large-scale B -field, producing a nearly one-dimensional flow of gas along the field lines. The trapped gas has a turbulent and magnetic topology that differs from that of the first phase and strongly suppresses gas accretion onto protostellar cores, reducing their masses. We speculate that this star formation bimodality may be an important ingredient to understand the origin of multiple stellar populations observed in massive globular clusters. [ABSTRACT FROM AUTHOR]

The distribution of stars and stellar remnants (white dwarfs, neutron stars, and black holes) within globular clusters holds clues about their formation and long-term evolution, with important implications for their initial mass function (IMF) and the formation of black hole mergers. In this work, we present best-fitting multimass models for 37 Milky Way globular clusters, which were inferred from various data sets, including proper motions from Gaia EDR3 and HST, line-of-sight velocities from ground-based spectroscopy and deep stellar mass functions from HST. We use metallicity-dependent stellar evolution recipes to obtain present-day mass functions of stars and remnants from the IMF. By dynamically probing the present-day mass function of all objects in a cluster, including the mass distribution of remnants, these models allow us to explore in detail the stellar (initial) mass functions of a large sample of Milky Way GCs. We show that, while the low-mass mass function slopes are strongly dependent on the dynamical age of the clusters, the high-mass slope (α3; m > 1 M⊙) is not, indicating that the mass function in this regime has generally been less affected by dynamical mass loss. Examination of this high-mass mass function slope suggests an IMF in this mass regime consistent with a Salpeter IMF is required to reproduce the observations. This high-mass IMF is incompatible with a top-heavy IMF, as has been proposed recently. Finally, based on multimass model fits to our sample of Milky Way GCs, no significant correlation is found between the high-mass IMF slope and cluster metallicity. [ABSTRACT FROM AUTHOR]

The distribution of stars and stellar remnants (white dwarfs, neutron stars, and black holes) within globular clusters holds clues about their formation and long-term evolution, with important implications for their initial mass function (IMF) and the formation of black hole mergers. In this work, we present best-fitting multimass models for 37 Milky Way globular clusters, which were inferred from various data sets, including proper motions from Gaia EDR3 and HST, line-of-sight velocities from ground-based spectroscopy and deep stellar mass functions from HST. We use metallicity-dependent stellar evolution recipes to obtain present-day mass functions of stars and remnants from the IMF. By dynamically probing the present-day mass function of all objects in a cluster, including the mass distribution of remnants, these models allow us to explore in detail the stellar (initial) mass functions of a large sample of Milky Way GCs. We show that, while the low-mass mass function slopes are strongly dependent on the dynamical age of the clusters, the high-mass slope (α3; m > 1 M⊙) is not, indicating that the mass function in this regime has generally been less affected by dynamical mass loss. Examination of this high-mass mass function slope suggests an IMF in this mass regime consistent with a Salpeter IMF is required to reproduce the observations. This high-mass IMF is incompatible with a top-heavy IMF, as has been proposed recently. Finally, based on multimass model fits to our sample of Milky Way GCs, no significant correlation is found between the high-mass IMF slope and cluster metallicity. [ABSTRACT FROM AUTHOR]

Characterizing pre-stellar cores in star-forming regions is an important step towards the validation of theoretical models of star formation. Due to their sub-arcsecond resolution, ALMA observations can potentially provide samples of pre-stellar cores up to distances of a few kpc, where regions of massive star formation can be targeted. However, the extraction of real cores from dust-continuum observations of turbulent star-forming clouds is affected by complex projection effects. In this work, we study the problem of core extraction both in the idealized case of column-density maps and in the more realistic case of synthetic 1.3 mm ALMA observations. The analysis is carried out on 12 regions of high column density from our 250 pc simulation. We find that derived core masses are highly unreliable, with only a weak correlation between the masses of cores selected in the synthetic ALMA maps and those of the corresponding three-dimensional cores. The fraction of real three-dimensional cores detected in the synthetic maps increases monotonically with mass and remains always below 50 per cent. Above |$\sim 1\, {\rm M}_{\odot }$| , the core mass function derived from the column-density maps is steeper than that of the three-dimensional cores, while the core mass function from the synthetic ALMA maps has a slope closer to that of the real three-dimensional cores. Because of the mass uncertainties, proper guidance from realistic simulations is essential if ALMA observations of protoclusters at kpc distances are to be used to test star formation models. [ABSTRACT FROM AUTHOR]

Low-mass star formation regions are unlikely to fully populate their initial mass functions (IMFs), leading to a deficit of massive stars. In binary stellar populations, the full range of binary separations and mass ratios will also be underpopulated. To explore the effects of stochastic sampling in the integrated light of stellar clusters, we calculate models at a broad range of cluster masses, from 102 to 107 M⊙, using a binary stellar population synthesis code. For clusters with stellar masses less than 105 M⊙, observable quantities show substantial scatter and their mean properties reflect the expected deficit of massive stars. In common with previous work, we find that purely stochastic sampling of the IMF appears to underestimate the mass of the most massive star in known clusters. However, even with this constraint, the majority of clusters likely inject sufficient kinetic energy to clear their birth clusters of gas. For quantities that directly measure the impact of the most massive stars, such as N ion, ξion, and βUV, uncertainties due to stochastic sampling dominate over those from the IMF shape or distribution of binary parameters, while stochastic sampling has a negligible effect on the stellar continuum luminosity density. [ABSTRACT FROM AUTHOR]

We derive the local stellar formation history from the Gaia -defined 40 pc white dwarf sample. This is currently the largest volume-complete sample of white dwarfs for which spectroscopy is available, allowing for classification of the chemical abundances at the photosphere, and subsequently accurate determination of the atmospheric parameters. We create a population synthesis model and show that a uniform stellar formation history for the last ≈10.5 Gyr provides a satisfactory fit to the observed distribution of absolute Gaia G magnitudes. To test the robustness of our derivation, we vary various assumptions in the population synthesis model, including the initial mass function, initial-to-final mass relation, kinematic evolution, binary fraction, and white dwarf cooling time-scales. From these tests, we conclude that the assumptions in our model have an insignificant effect on the derived relative stellar formation rate as a function of look-back time. However, the onset of stellar formation (age of Galactic disc) is sensitive to a variety of input parameters including the white dwarf cooling models. Our derived stellar formation history gives a much better fit to the absolute Gaia G magnitudes than most previous studies. [ABSTRACT FROM AUTHOR]

We have determined stellar mass functions of 120 Milky Way globular clusters and massive Large Magellanic Cloud/Small Magellanic Cloud star clusters based on a comparison of archival Hubble Space Telescope photometry with a large grid of direct N -body simulations. We find a strong correlation of the global mass function slopes of star clusters with both their internal relaxation times and their lifetimes. Once dynamical effects are being accounted for, the mass functions of most star clusters are compatible with an initial mass function described by a broken power-law distribution N (m) ∼ m α with break masses at 0.4 and 1.0 M⊙ and mass function slopes of αLow = −0.3 for stars with masses m < 0.4 M⊙, αHigh = −2.30 for stars with m > 1.0 M⊙, and αMed = −1.65 for intermediate-mass stars. Alternatively, a lognormal mass function with a characteristic mass log  M C = −0.36 and width σC = 0.28 for low-mass stars and a power-law mass function for stars with m > 1 M⊙ also fit our data. We do not find a significant environmental dependence of the initial mass function on cluster mass, density, global velocity dispersion, or metallicity. Our results lead to a larger fraction of high-mass stars in globular clusters compared to canonical Kroupa/Chabrier mass functions, increasing the efficiency of self-enrichment in clusters and helping to alleviate the mass budget problem of multiple stellar populations in globular clusters. By comparing our results with direct N -body simulations, we finally find that only simulations in which most black holes are ejected by natal birth kicks correctly reproduce the observed correlations. [ABSTRACT FROM AUTHOR]

Molecular hydrogen allows cooling in primordial gas, facilitating its collapse into Population III stars within primordial halos. Lyman–Werner (LW) radiation from these stars can escape the halo and delay further star formation by destroying H2 in other halos. As cosmological simulations show that increasing the background LW field strength increases the average halo mass required for star formation, we perform follow-up simulations of selected halos to investigate the knock-on effects this has on the Population III IMF. We follow 5 halos for each of the J 21  = 0, 0.01, and 0.1 LW field strengths, resolving the pre-stellar core density of 10−6 g cm−3 (1018 cm−3) before inserting sink particles and following the fragmentation behaviour for hundreds of years further. We find that the mass accreted onto sinks by the end of the simulations is proportional to the mass within the ∼10−2 pc molecular core, which is not correlated to the initial mass of the halo. As such, the IMFs for masses above the brown dwarf limit show little dependence on the LW strength, although they do show variance in the number of low-mass clumps formed. As the range of background LW field strengths tested here covers the most likely values from literature, we conclude that the IMF for so-called Pop III.2 stars is not significantly different from the initial population of Pop III.1 stars. The primordial IMF therefore likely remains unchanged until the formation of the next generation of Population II stars. [ABSTRACT FROM AUTHOR]

Context. Star-forming regions are excellent benchmarks for testing and validating theories of star formation and stellar evolution. The Perseus star-forming region, being one of the youngest (< 10 Myr), closest (280−320 pc), and most studied in the literature, is a fundamental benchmark. Aims. We aim to study the membership, phase-space structure, mass, and energy (kinetic plus potential) distribution of the Perseus star-forming region using public catalogues (Gaia, APOGEE, 2MASS, and Pan-STARRS). Methods. We used Bayesian methodologies that account for extinction to identify the Perseus physical groups in the phase-space, retrieve their candidate members, derive their properties (age, mass, 3D positions, 3D velocities, and energy), and attempt to reconstruct their origin. Results. We identify 1052 candidate members in seven physical groups (one of them new) with ages between 3 and 10 Myr, dynamical super-virial states, and large fractions of energetically unbounded stars. Their mass distributions are broadly compatible with that of Chabrier for masses ≳0.1 M⊙ and do not show hints of over-abundance of low-mass stars in NGC 1333 with respect to IC 348. These groups' ages, spatial structure, and kinematics are compatible with at least three generations of stars. Future work is still needed to clarify if the formation of the youngest was triggered by the oldest. Conclusions. The exquisite Gaia data complemented with public archives and mined with comprehensive Bayesian methodologies allow us to identify 31% more members than previous studies, discover a new physical group (Gorgophone: 7 Myr, 191 members, and 145 M⊙), and confirm that the spatial, kinematic, and energy distributions of these groups support the hierarchical star formation scenario. [ABSTRACT FROM AUTHOR]

We present a spectrophotometric analysis of 2880 cool white dwarfs within 100 pc of the Sun and cooler than T eff ∼ 10 000 K, with grizy Pan-STARRS photometry and Gaia trigonometric parallaxes available. We also supplement our data sets with near-infrared JHK photometry, when available, which is shown to be essential for interpreting the coolest white dwarfs in our sample. We perform a detailed analysis of each individual object using state-of-the-art model atmospheres appropriate for each spectral type, including DA, DC, DQ, DZ, He-rich DA, and the so-called IR-faint white dwarfs. We discuss the temperature and mass distributions of each subsample, as well as revisit the spectral evolution of cool white dwarfs. We find little evidence in our sample for the transformation of a significant fraction of DA stars into He-atmosphere white dwarfs through the process of convective mixing between T eff = 10 000 and ∼6500 K, although the situation changes drastically in the range T eff = 6500–5500 K where the fraction of He-atmosphere white dwarfs reaches ∼45 per cent. However, we also provide strong evidence that at even cooler temperatures (T eff ≲ 5200 K), most DC white dwarfs have H atmospheres. We discuss a possible mechanism to account for this sudden transformation from He- to H-atmosphere white dwarfs involving the onset of crystallization and the occurrence of magnetism. Finally, we also argue that DQ, DZ, and DC white dwarfs may form a more homogeneous population than previously believed. [ABSTRACT FROM AUTHOR]

We report the statistical properties of stars and brown dwarfs obtained from three radiation hydrodynamical simulations of star cluster formation with metallicities of 1, 1/10, and 1/100 of the solar value. The star-forming clouds are subjected to cosmic microwave background radiation that is appropriate for star formation at a redshift z  = 5. The results from the three calculations are compared to each other, and to similar previously published calculations that had levels of background radiation appropriate for present-day (z  = 0) star formation. Each of the calculations treats dust and gas temperatures separately and includes a thermochemical model of the diffuse interstellar medium. We find that whereas the stellar mass distribution is insensitive to the metallicity for present-day star formation, at z  = 5 the characteristic stellar mass increases with increasing metallicity and the mass distribution has a deficit of brown dwarfs and low-mass stars at solar metallicity compared to the Galactic initial mass function. We also find that the multiplicity of M-dwarfs decreases with increasing metallicity at z  = 5. These effects are a result of metal-rich gas being unable to cool to as low temperatures at z  = 5 compared to at z  = 0 due to the hotter cosmic microwave background radiation, which inhibits fragmentation at high densities. [ABSTRACT FROM AUTHOR]

Star formation models predict that the metal-poor initial mass function (IMF) can be substantially different from that observed in the metal-rich Milky Way. This changeover occurs because metal-poor gas clouds cool inefficiently due to their lower abundance of metals and dust. However, predictions for the metal-poor IMF to date rely on assuming solar-scaled abundances, i.e. [X/O] = 0 at all [O/H]. There is now growing evidence that elements such as C and O that dominate metal line cooling in the ISM do not follow solar scaling at low metallicities. In this work, we extend models that predict the variation in the characteristic (or the peak) IMF mass as a function of metallicity using [C/O] ratios derived from observations of metal-poor Galactic stars and of H  ii regions in dwarf galaxies. These data show [C/O] < 0 at subsolar [O/H], which leads to a substantially different metal-poor IMF in the metallicity range where C  i and C  ii cooling dominate ISM thermodynamics, resulting in an increase in the characteristic mass by a factor as large as 7. An important consequence of this difference is a shift in the location of the transition from a top- to a bottom-heavy IMF upwards by 0.5–1 dex in metallicity. Our findings indicate that the IMF is very sensitive to the assumptions around solar-scaled ISM compositions in metal-poor systems (e.g. dwarf galaxies, the Galactic halo, and metal-poor stars) that are a key focus of JWST. [ABSTRACT FROM AUTHOR]

We explore the assumption, widely used in many astrophysical calculations, that the stellar initial mass function (IMF) is universal across all galaxies. By considering both a canonical broken-power-law IMF and a non-universal IMF, we are able to compare the effect of different IMFs on multiple observables and derived quantities in astrophysics. Specifically, we consider a non-universal IMF that varies as a function of the local star formation rate, and explore the effects on the star formation rate density (SFRD), the extragalactic background light, the supernova (both core-collapse and thermonuclear) rates, and the diffuse supernova neutrino background. Our most interesting result is that our adopted varying IMF leads to much greater uncertainty on the SFRD at |$z \approx 2-4$| than is usually assumed. Indeed, we find an SFRD (inferred using observed galaxy luminosity distributions) that is a factor of |$\gtrsim 3$| lower than canonical results obtained using a universal IMF. Secondly, the non-universal IMF we explore implies a reduction in the supernova core-collapse rate of a factor of |$\sim 2$|⁠ , compared against a universal IMF. The other potential tracers are only slightly affected by changes to the properties of the IMF. We find that currently available data do not provide a clear preference for universal or non-universal IMF. However, improvements to measurements of the star formation rate and core-collapse supernova rate at redshifts |$z \gtrsim 2$| may offer the best prospects for discernment. [ABSTRACT FROM AUTHOR]

We use a series of magnetohydrodynamic simulations including both radiative and protostellar outflow feedback to study environmental variation of the initial mass function (IMF). The simulations represent a carefully-controlled experiment whereby we keep all dimensionless parameters of the flow constant except for those related to feedback. We show that radiation feedback suppresses the formation of lower mass objects more effectively as the surface density increases, but this only partially compensates for the decreasing Jeans mass in denser environments. Similarly, we find that protostellar outflows are more effective at suppressing the formation of massive stars in higher surface density environments. The combined effect of these two trends is towards an IMF with a lower characteristic mass and a narrower overall mass range in high surface density environments. We discuss the implications for these findings for the interpretation of observational evidence of IMF variation in early type galaxies. [ABSTRACT FROM AUTHOR]

We determine the mass functions (MFs) and the dynamical parameters of 15 nearby open clusters (OCs), using the unprecedented data set of the Gaia Early Data Release 3. We select the members of each cluster by combining the photometric (colour and magnitude) and astrometric (parallax and proper motions) parameters of stars, minimizing the contamination from Galactic field interlopers. By comparing the observed distribution of stars along the cluster main sequence with the best-fitting synthetic population, we find the present-day MF and the binary fraction of the OCs, along with their dynamical parameters like mass, half-mass radius, and half-mass relaxation time. We found that the global present-day MF of OCs are consistent with a single power-law function, F (m) ∝  m α, with slopes −3 < α < −0.6 including both subsolar, 0.2 < m /M⊙ < 1, and supersolar mass regimes. A significant correlation between the MF slope and the ratio of age to half-mass relaxation time is evidenced, similarly to the same conclusion already observed among Galactic globular clusters. However, OCs evolve along different tracks in comparison with the globular clusters, possibly indicating primordial differences in their initial mass function (IMF). The comparison with Monte Carlo simulations suggests that all the analysed OCs could be born with an IMF with slope αIMF < −2.3. We also show that the less evolved OCs have an MF consistent with that of the solar neighbourhood, indicating a possible connection between the dissolution of OCs and the formation of the Galactic disc. [ABSTRACT FROM AUTHOR]

Population III (Pop III) stars form in groups due to the fragmentation of primordial gas. While uniform magnetic fields have been shown to support against fragmentation in present-day star formation, it is unclear whether realistic k 3/2 primordial fields can have the same effect. We bypass the issues associated with simulating the turbulent dynamo by introducing a saturated magnetic field at equipartition with the velocity field when the central densities reach 10−13 g cm−3. We test a range of sink particle creation densities from 10−10 to 10−8 g cm−3. Within the range tested, the fields did not suppress fragmentation of the gas and hence could not prevent the degree of fragmentation from increasing with increased resolution. The number of sink particles formed and total mass in sink particles was unaffected by the magnetic field across all seed fields and resolutions. The magnetic pressure remained sub-dominant to the gas pressure except in the highest density regions of the simulation box, where it became equal to but never exceeded gas pressure. Our results suggest that the inclusion of magnetic fields in numerical simulations of Pop III star formation is largely unimportant. [ABSTRACT FROM AUTHOR]

One of the key mysteries of star formation is the origin of the stellar initial mass function (IMF). The IMF is observed to be nearly universal in the Milky Way and its satellites, and significant variations are only inferred in extreme environments, such as the cores of massive elliptical galaxies and the Central Molecular Zone. In this work, we present simulations from the STARFORGE project that are the first cloud-scale radiation-magnetohydrodynamic simulations that follow individual stars and include all relevant physical processes. The simulations include detailed gas thermodynamics, as well as stellar feedback in the form of protostellar jets, stellar radiation, winds, and supernovae. In this work, we focus on how stellar radiation, winds, and supernovae impact star-forming clouds. Radiative feedback plays a major role in quenching star formation and disrupting the cloud; however, the IMF peak is predominantly set by protostellar jet physics. We find that the effect of stellar winds is minor, and supernovae 'occur too late' to affect the IMF or quench star formation. We also investigate the effects of initial conditions on the IMF. We find that the IMF is insensitive to the initial turbulence, cloud mass, and cloud surface density, even though these parameters significantly shape the star formation history of the cloud, including the final star formation efficiency. Meanwhile, the characteristic stellar mass depends weakly on metallicity and the interstellar radiation field, which essentially set the average gas temperature. Finally, while turbulent driving and the level of magnetization strongly influence the star formation history, they only influence the high-mass slope of the IMF. [ABSTRACT FROM AUTHOR]

Stars form in dense, clustered environments, where feedback from newly formed stars eventually ejects the gas, terminating star formation and leaving behind one or more star clusters. Using the STARFORGE simulations, it is possible to simulate this process in its entirety within a molecular cloud, while explicitly evolving the gas radiation and magnetic fields and following the formation of individual, low-mass stars. We find that individual star-formation sites merge to form ever larger structures, while still accreting gas. Thus clusters are assembled through a series of mergers. During the cluster assembly process, a small fraction of stars are ejected from their clusters; we find no significant difference between the mass distribution of the ejected stellar population and that of stars inside clusters. The star-formation sites that are the building blocks of clusters start out mass segregated with one or a few massive stars at their centre. As they merge the newly formed clusters maintain this feature, causing them to have mass-segregated substructures without themselves being centrally condensed. The merged clusters relax to a centrally condensed mass-segregated configuration through dynamical interactions between their members, but this process does not finish before feedback expels the remaining gas from the cluster. In the simulated runs, the gas-free clusters then become unbound and breakup. We find that turbulent driving and a periodic cloud geometry can significantly reduce clustering and prevent gas expulsion. Meanwhile, the initial surface density and level of turbulence have little qualitative effect on cluster evolution, despite the significantly different star formation histories. [ABSTRACT FROM AUTHOR]

Mapping cosmic dawn with 21-cm tomography offers an exciting new window into the era of primordial star formation. However, self-consistent implementation of both the process of star formation and the related 21-cm signal is challenging, due to the multiscale nature of the problem. In this study, we develop a flexible semi-analytical model to follow the formation of the first stars and the process of gradual transition from primordial to metal-enriched star formation. For this transition, we use different scenarios with varying time-delays (or recovery times) between the first supernovae and the formation of the second generation of stars. We use recovery times between 10 and 100 Myr and find that these delays have a strong impact on the redshift at which the transition to metal-enriched star formation occurs. We then explore the effect of this transition on the 21-cm signal and find that the recovery time has a distinctive imprint in the signal. Together with an improved understanding of how this time-delay relates to the properties of Population III stars, future 21-cm observations can give independent constraints on the earliest epoch of star formation. As the transition away from the primordial star formation is expected to occur at high redshifts, here we ignore the impact of X-ray and ionizing radiation, focusing on the effect of Ly α photons on the 21-cm signal. [ABSTRACT FROM AUTHOR]

In this work, we restudy the dependence of luminosity function and event rates for different gamma-ray burst samples on the criteria of sample selection and threshold effect. To compare with many previous studies, we have chosen two samples including 88 and 118 long bursts with known redshift and peak flux over 2.6 ph cm−2 s−1, from which 79 bursts are picked out to constitute our complete sample. It is found that the evolution of luminosity with redshift can be expressed by L ∝(1 + z) k with a diverse k relied more on the sample selection. Interestingly, the cumulative distributions of either non-evolving luminosities or redshifts are found to be also determined by the sample selection rather than the instrumental sensitivity. Nevertheless, the non-evolving luminosities of our samples are similarly distributed with a comparable break luminosity of L 0 ∼ 1051 erg s−1. Importantly, we verify with a K-S test that three cases of event rates for the two burst samples evolve with redshift similarly except a small discrepancy due to sampling differences at low-redshift of z < 1, in which all event rates show an excess of Gaussian profile instead of monotonous decline no matter whether the sample is complete. Most importantly, it is found that the burst rates violate the star formation rate at low redshift, while both of them are good in agreement with each other in the higher-redshift regions as many authors discovered previously. Therefore, we predict that two types of long bursts are favored in terms of their associations with both the star formation and the cosmic metallicity. [ABSTRACT FROM AUTHOR]

Recent observations of active galactic nuclei (AGNs) have shown a high Fe  ii /Mg  ii line-flux ratio in their broad-line regions, nearly independent of redshift up to z ≳ 6. The high flux ratio requires rapid production of iron in galactic nuclei to reach an abundance ratio of [Fe/Mg] ≳ 0.2 as high as those observed in matured galaxies in the local universe. We propose a possible explanation of rapid iron enrichment in AGNs by massive star formation that follows a top-heavy initial mass function (IMF) with a power-law index of Γ larger than the canonical value of Γ = −2.35 for a Salpeter IMF. Taking into account metal production channels from different types of SNe, we find that the high value of [Fe/Mg] ≳ 0.2 requires the IMF to be characterized with Γ ≳ −1 (Γ ≳ 0) and a high-mass cutoff at M max ≃ 100–150 M⊙ (M max ≳ 250 M⊙). Given the conditions, core-collapse SNe with M * ≳ 70 M⊙ and pair-instability SNe give a major contribution for iron enrichment. Such top-heavy stellar IMFs would be a natural consequence from mass growth of stars formed in dense AGN discs under Bondi-like gas accretion that is regulated by feedback at M * ≳ 10 M⊙. The massive stellar population formed in AGN discs also leave stellar-mass black hole remnants, whose mergers associated with gravitational-wave emission account for at most 10  per cent of the merger rate inferred from LIGO/Virgo observations to simultaneously explain the high [Fe/Mg] ratio with metal ejection. [ABSTRACT FROM AUTHOR]

We present a new SOFIA / FORCAST mid-infrared survey of luminous protostars and crowded star-forming environments in Cygnus X, the nearest million-solar mass molecular cloud complex. We derive bolometric luminosities for over 1000 sources in the region with these new data in combination with extant Spitzer and UKIDSS photometry, with 63 new luminous protostar candidates identified by way of the high-quality SOFIA / FORCAST data. By including FORCAST data, we construct protostellar luminosity functions (PLFs) with improved completeness at the high luminosity end. The PLFs are well described by a power-law function with an index of ∼−0.5. Based on the Herschel temperature and column density measurements, we find no obvious dependence of the PLFs on the local gas temperature, but PLFs in regions of high stellar density or gas column density exhibit some excess at higher luminosities. Through the comparison between our observed PLFs and existing accretion models, both the turbulent core and the competitive accretion models are consistent with our results, while the isothermal sphere model is disfavoured. The implications of these results on the star formation process are discussed. [ABSTRACT FROM AUTHOR]

We determine the distribution of cooling ages of massive Gaia EDR3 white dwarfs identified with over 90 per cent probability within 200 pc and with mass in the range 0.95–1.25 M⊙. Using three sets of publicly available models, we consider sub-samples of these white dwarfs sorted into three equally spaced mass bins. Under the assumption of a constant white dwarf formation rate, we find an excess of white dwarfs, both along the Q branch and below it, corresponding respectively to stars that are in the process of freezing and those that are completely frozen. We compare the cooling age distributions for each of these bins to the recently determined time-varying star formation rate of Gaia DR2 main sequence stars. For white dwarfs in the two lightest mass bins, spanning the mass range 0.95–1.15 M⊙, we find that the cumulative cooling age distribution is statistically consistent with the expectation from the star formation rate. For white dwarfs in the heaviest mass bin, 1.15–1.25 M⊙, we find that their cumulative distribution is inconsistent with the star formation rate for all of the models considered; instead, we find that their cooling age distribution is well fitted by a linear combination of the distribution expected for single stellar evolution products and the distribution expected for double white dwarf merger products when approximately 40–50 per cent of the 1.15–1.25  M ⊙ white dwarfs that formed over the past 4 Gyr are produced through double white dwarf mergers. [ABSTRACT FROM AUTHOR]

Type Ia supernovae, SNeIa, are used as standard candles in cosmology to determine the distances of the galaxies harbouring them. We show that the luminosity of an SNIa depends on its distance from us when physical constants (the speed of light c , the gravitational constant G , and the Planck constant h) are permitted to evolve. It is because the Chandrasekhar mass of the white dwarf that explodes to create SNIa depends on the values of the constants at the epoch the SNIa is formed. We show that the SNeIa luminosities could be about four times higher in the past than they are now. Thus, the luminosity distance estimation of the earliest SNeIa could be off by up to a factor of 2. Cosmological parameters, determined with this correction applied to the redshift versus distance modulus data base (Pantheon SNeIa), are not very different from those from the standard ΛCDM model without this correction, except for the dark-energy density and the curvature energy density; the latter increases at the cost of the former. Variations of the constants are given by |$\dot{G}/G = \ 3.90 \ ({ \pm 0.04}) \times {10^{ - 10}}\ {\rm y{r^{ - 1}}}$| and |$\dot{c}/c = \dot{h}/h\ = \ 1.30\ ({ \pm 0.01}) \times {10^{ - 10}}\ {\rm y{r^{ - 1}}}$| at present. These variations are valid only when |$G,\ c,\ $| and |$\ h$| are permitted to vary concurrently rather than individually. [ABSTRACT FROM AUTHOR]

A tidal disruption event (TDE) occurs when a star is destroyed by the strong tidal shear of a massive black hole (MBH). The accumulation of TDE observations over the last years has revealed that post-starburst galaxies are significantly overrepresented in the sample of TDE hosts. Here we address the post-starburst preference by investigating the decline of TDE rates in a Milky-Way like nuclear stellar cluster featuring either a monochromatic (1 |$\, \mathrm{M}{\odot {}}$|⁠) or a complete, evolved stellar mass function. In the former case, the decline of TDE rates with time is very mild, and generally up to a factor of a few in 10 Gyr. Conversely, if a complete mass function is considered, a strong TDE burst over the first 0.1–1 Gyr is followed by a considerable rate drop, by at least an order of magnitude over 10 Gyr. The decline starts after a mass segregation time-scale, and it is more pronounced assuming a more top-heavy initial mass function and/or an initially denser nucleus. Our results thus suggest that the post-starburst preference can be accounted for in realistic systems featuring a complete stellar mass function, even in moderately dense galactic nuclei. Overall, our findings support the idea that starbursting galactic nuclei are characterized by a top-heavy initial mass function; we speculate that accounting for this can reconcile the discrepancy between observed and theoretically predicted TDE rates even in quiescent galaxies. [ABSTRACT FROM AUTHOR]

We present a pilot library of synthetic NUV, U, B, V , and I photometry of star clusters with stochastically sampled IMFs and ionized gas for initial masses, Mi  = 103, 104, and 105 M⊙; t  = 1, 3, 4, and 8 Myr; Z  = 0.014 and Z  = 0.002; and log(U S) = −2 and −3. We compare the library with predictions from deterministic models and observations of isolated low-mass (<104 M⊙) star clusters with co-spatial compact H  ii  regions. The clusters are located in NGC 7793, one of the nearest galaxies observed as part of the HST  LEGUS and Hα-LEGUS surveys. (1) For model magnitudes that only account for the stars: (a) the residual |deterministic mag - median stochastic mag| can be ≥0.5 mag, even for Mi  = 105 M⊙; and (b) the largest spread in stochastic magnitudes occurs when Wolf–Rayet stars are present. (2) For Mi  = 105 M⊙: (a) the median stochastic mag with gas can be >1.0 mag more luminous than the median stochastic magnitude without gas; and (b) nebular emission lines can contribute with |$\gt 50{{\ \rm per\ cent}}$| and |$\gt 30{{\ \rm per\ cent}}$| to the total emission in the V and I bands, respectively. (3) Age-dating OB-star clusters via deterministic tracks in the U - B versus V - I plane is highly uncertain at Z  = 0.014 for Mi ∼ 103 M⊙ and Z  = 0.002 for Mi ∼ 103–105 M⊙. (4) For low-mass clusters, the V -band extinction derived with stochastic models significantly depends on the value of log(U S). (5) The youngest clusters tend to have higher extinction. (6) The majority of clusters have multi-peaked age PDFs. (7) Finally, we discuss the importance of characterizing the true variance in the number of stars per mass bin in nature. [ABSTRACT FROM AUTHOR]

The evolution of white dwarfs is essentially a gravothermal process of cooling in which the basic ingredients for predicting their evolution are well identified, although not always well understood. There are two independent ways to test the cooling rate. One is the luminosity function of the white dwarf population, and another is the secular drift of the period of pulsation of those individuals that experience variations. Both scenarios are sensitive to the cooling or heating time scales, for which reason, the inclusion of any additional source or sink of energy will modify these properties and will allow to set bounds to these perturbations. These studies also require complete and statistical significant samples for which current large data surveys are providing an unprecedented wealth of information. In this paper we review how these techniques are applied to several cases like the secular drift of the Newton gravitational constant, neutrino magnetic moments, axions and weakly interacting massive particles (WIMPS).

As star-forming clouds collapse, the gas within them fragments to ever-smaller masses. Naively one might expect this process to continue down to the smallest mass that is able to radiate away its binding energy on a dynamical time-scale, the opacity limit for fragmentation, at ~0.01 M⊙. However, the observed peak of the initial mass function (IMF) lies a factor of 20-30 higher in mass, suggesting that some other mechanism halts fragmentation before the opacity limit is reached. In this paper we analyse radiation-magnetohydrodynamic simulations of star cluster formation in typical Milky Way environments in order to determine what physical process limits fragmentation in them. We examine the regions in the vicinity of stars that form in the simulations to determine the amounts of mass that are prevented from fragmenting by thermal and magnetic pressure. We show that, on small scales, thermal pressure enhanced by stellar radiation heating is the dominant mechanism limiting the ability of the gas to further fragment. In the brown dwarf mass regime, ~0.01 M⊙, the typical object that forms in the simulations is surrounded by gas whose mass is several times its own that is unable to escape or fragment, and instead is likely to accrete. This mechanism explains why ~0.01 M⊙ objects are rare: unless an outside agent intervenes (e.g. a shock strips away the gas around them), they will grow by accreting the warmed gas around them. In contrast, by the time stars grow to masses of ~0.2 M⊙, the mass of heated gas is only tens of percent of the central star mass, too small to alter its final mass by a large factor. This naturally explains why the IMF peak is at ~0.2 M⊙.

We use Gaia DR2 to measure the initial mass function (IMF) of stars within 250 pc and masses in the range of 0.2 < m /M⊙ < 1.0, separated according to kinematics and metallicity, as determined from Gaia transverse velocity, v T, and location on the Hertzsprung–Russell diagram (HRD). The predominant thin-disc population (v T < 40 km s−1) has an IMF similar to traditional (e.g. Kroupa) stellar IMFs, with star numbers per mass interval d N /d m described by a broken power law, m −α, and index |$\alpha _\textrm {high}=2.03^{+0.14}_{-0.05}$| above m ∼ 0.5, shallowing to |$\alpha _\textrm {low}=1.34^{+0.11}_{-0.22}$| at m ≲ 0.5. Thick-disc stars (60 km s−1 < v T < 150 km s−1) and stars belonging to the 'high-metallicity' or 'red-sequence' halo (v T > 100 km s−1 or v T > 200 km s−1, and located above the isochrone on the HRD with a metallicity [M/H] > −0.6) have a somewhat steeper high-mass slope, |$\alpha _\textrm {high}=2.35^{+0.97}_{-0.19}$| (and a similar low-mass slope |$\alpha _\textrm {low}=1.14^{+0.42}_{-0.50}$|⁠). Halo stars from the 'blue sequence', which are characterized by low metallicity ([M/H] < −0.6), however, have a distinct, bottom-heavy IMF, well described by a single power law with |$\alpha =1.82^{+0.17}_{-0.14}$| over most of the mass range probed. The IMF of the low-metallicity halo is reminiscent of the Salpeter-like IMF that has been measured in massive early-type galaxies, a stellar population that, like Milky Way halo stars, has a high ratio of α elements to iron, [α/Fe]. Blue-sequence stars are likely the debris from accretion by the Milky Way, ∼10 Gyr ago, of the Gaia-Enceladus dwarf galaxy, or similar events. These results hint at a distinct mode of star formation common to two ancient stellar populations – elliptical galaxies and galaxies possibly accreted early-on by ours. [ABSTRACT FROM AUTHOR]

Next-generation telescopes such as the James Webb Space Telescope (JWST) and the Nancy Grace Roman Space Telescope (NGRST) will enable us to study the first billion years of our Universe in unprecedented detail. In this work, we use the astraeus (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework, which couples galaxy formation and reionization (for a wide range of reionization feedback models), to estimate the cosmic variance expected in the UV Luminosity Function (UV LF) and the stellar mass function in JWST surveys. We find that different reionization scenarios play a minor role in the cosmic variance. Most of the cosmic variance is completely driven by the underlying density field and increases above |$100{{\ \rm per\ cent}}$| for |${\it M}_{\rm UV}\sim -17.5 ~ (-20)$| at z  = 12 (6) for the JADES-deep survey (the deep JWST Advanced Extragalactic Survey with an area of 46 arcmin2); the cosmic variance decreases with an increasing survey area roughly independently of redshift. We find that the faint-end (⁠|${\it M}_{\rm UV}\gtrsim -17$|⁠) slope of the Lyman Break Galaxies (LBGs) UV LF becomes increasingly shallower with increasing reionization feedback and show how JWST observations will be able to distinguish between different models of reionization feedback at z > 9, even accounting for cosmic variance. We also show the environments (in terms of density and ionization fields) of LBGs during the EoR, finding that the underlying overdensity and ionization fraction scale positively with the UV luminosity. Finally, we also provide a public software tool to allow interested readers to compute cosmic variance for different redshifts and survey areas. [ABSTRACT FROM AUTHOR]

Statistical methods are frequently built upon assumptions that limit their applicability to certain problems and conditions. Failure to recognize these limitations can lead to conclusions that may be inaccurate or biased. An example of such methods is the non-parametric Efron–Petrosian test statistic used in the studies of truncated data. We argue and show how the inappropriate use of this statistical method can lead to biased conclusions when the assumptions under which the method is valid do not hold. We do so by reinvestigating the evidence recently provided by multiple independent reports on the evolution of the luminosity/energetics distribution of cosmological Long-duration Gamma-Ray Bursts (LGRBs) with redshift. We show that the effects of detection threshold have been likely significantly underestimated in the majority of previous studies. This underestimation of detection threshold leads to severely incomplete LGRB samples that exhibit strong apparent luminosity-redshift or energetics-redshift correlations. We further confirm our findings by performing extensive Monte Carlo simulations of the cosmic rates and the luminosity/energy distributions of LGRBs and their detection process. [ABSTRACT FROM AUTHOR]

The surface brightness fluctuation (SBF) technique is one of the distance measurement methods that has been applied on the low surface brightness (LSB) galaxy NGC 1052-DF2 yielding a distance of about 20 Mpc implying it to be a dark matter deficient galaxy. We assume the number of giant stars above a given luminosity threshold to represent the SBF magnitude. The SBF magnitude depends on the distance, but this is degenerate with the star formation history (SFH). Using a stellar population synthesis model, we calculate the number of giant stars for stellar populations with different galaxy-wide stellar initial mass functions (gwIMFs), ages, metallicities, and SFHs. If the gwIMF is the invariant canonical IMF, the 1σ (3σ) uncertainty in colour allows a distance as low as 12 Mpc (8 Mpc). If instead the true underlying gwIMF is the integrated galaxy-wide IMF (IGIMF) then overestimating distances for low-mass galaxies would be a natural result, allowing NGC 1052-DF2 to have a distance of 11 Mpc within the 1σ colour uncertainty. Finally, we show that our main conclusion on the existence of a bias in the SBF distance estimation is not much affected by changing the luminosity lower limit for counting giant stars. [ABSTRACT FROM AUTHOR]

In the star-formation process, the vital impact of environmental factors such as feedback from massive stars and stellar density on the form of the initial mass function (IMF) at the low-mass end is yet to be understood. Hence a systematic highly sensitive observational analysis of a sample of regions under diverse environmental conditions is essential. We analyse the IMF of eight young clusters (<5 Myr), namely IC 1848-West, IC 1848-East, NGC 1893, NGC 2244, NGC 2362, NGC 6611, Stock 8, and Cygnus OB2, which are located at the Galactocentric distance (R g) range ∼6–12 kpc along with the nearby cluster IC 348 using deep near-IR photometry and Gaia DR2. These clusters are embedded in massive stellar environments of radiation strength |$\log(L_\mathrm{FUV}/\mathrm{L}_{\odot })\, \sim 2.6$| –6.8, |$\log(L_\mathrm{EUV})\, \sim$|  42.2–50.85 photon s−1, with stellar density in the range of ∼170–1220 star pc−2. After structural analysis and field decontamination we obtain an unbiased uniformly sensitive sample of pre-main-sequence members of the clusters down to the brown-dwarf regime. The lognormal fit to the IMF of nine clusters gives the mean characteristic mass (m c) and σ of 0.32 ± 0.02 M⊙ and 0.47 ± 0.02, respectively. We compare the IMF with that of low- and high-mass clusters across the Milky Way. We also check for any systematic variation with respect to the radiation field strength and the stellar density as well with R g. We conclude that there is no strong evidence for an environmental effect in the underlying form of the IMF of these clusters. [ABSTRACT FROM AUTHOR]

The non-detection of zero-metallicity stars in ultra-faint dwarf galaxies (UFDs) can be used to constrain the initial mass function (IMF) of the first (PopIII) stars by means of a statistical comparison between available data and predictions from chemical evolution models. To this end we develop a model that follows the formation of isolated UFDs, calibrated to best reproduce the available data for the best studied system: Boötes I. Our statistical approach shows that UFDs are the best suitable systems to study the implications of the persisting non-detection of zero-metallicity stars on the PopIII IMF, i.e. its shape, the minimum mass (m min), and the characteristic mass (m ch). We show that accounting for the incomplete sampling of the IMF is essential to compute the expected number of long-lived PopIII stars in inefficiently star-forming UFDs. By simulating the colour–magnitude diagram of Boötes I, and thus take into account the mass-range of the observed stars, we can obtain even tighter constrains on m min. By exploiting the 96 stars with measured metallicities (⁠|$\rm i \lt 19$|⁠) in the UFDs represented by our model, we demonstrate that |$m_{\mathrm{ ch}} \gt 1\: \rm {M_{\odot }}$| or |$m_{\mathrm{ min}} \gt 0.8 \:\rm {M_{\odot }}$| at |$99\%$| confidence level. This means that a present-day IMF for PopIII stars is excluded by our model, and a top-heavy PopIII IMF is strongly favoured. We can limit |$m_{\mathrm{ min}} \gt 0.8\: \rm {M_{\odot }}$| independent of the PopIII IMF shape by targeting the four UFDs Boötes I, Hercules, Leo IV, and Eridanus II with future generation instruments, such as ELT/MOSAIC (⁠|$\rm i \lt 25$|⁠), which can provide samples of >10 000 stars. [ABSTRACT FROM AUTHOR]

We aim at identifying very low-mass isolated planetary-mass member candidates in the nearest OB association to the Sun, Upper Scorpius (USco) (145 pc; 5–10 Myr), to constrain the form and shape of the luminosity function and mass spectrum in this regime. We conducted a deep multiband (Y  = 21.2, J  = 20.5, Z  = 22.0 mag) photometric survey of 6 deg2 in the central region of USco. We extend the current sequence of astrometric and spectroscopic members by about two magnitudes in Y and one magnitude in J , reaching potentially T-type free-floating members in the association with predicted masses below 5 Jupiter masses, well into the planetary-mass regime. We extracted a sample of 57 candidates in this area and present infrared spectroscopy confirming two of them as young L-type members with characteristic spectral features of 10-Myr-old brown dwarfs. Among the 57 candidates, we highlight 10 new candidates fainter than the coolest members previously confirmed spectroscopically. We do not see any obvious sign of decrease in the mass spectrum of the association, suggesting that star processes can form substellar objects with masses down to 4–5 Jupiter masses. [ABSTRACT FROM AUTHOR]