Stellar Population Synthesis Models Research Articles

The delay time distribution of Type-Ia supernovae in galaxy clusters: the impact of extended star-formation histories

ABSTRACT The delay time distribution (DTD) of Type-Ia supernovae (SNe Ia) is important for understanding chemical evolution, SN Ia progenitors, and SN Ia physics. Past estimates of the DTD in galaxy clusters have been deduced from SN Ia rates measured in cluster samples observed at various redshifts, corresponding to different time intervals after a presumed initial brief burst of star formation. A recent analysis of a cluster sample at z = 1.13–1.75 confirmed indications from previous studies of lower redshift clusters, that the DTD has a power-law form, DTD(t) = R1(t/Gyr)α, with amplitude R1, at delay $t=1\,\rm Gyr$, several times higher than measured in field-galaxy environments. This implied that SNe Ia are somehow produced in larger numbers by the stellar populations in clusters. This conclusion, however, could have been affected by the implicit assumption that the stars were formed in a single brief starburst at high z. Here, we re-derive the DTD from the cluster SN Ia data, but relax the single-burst assumption. Instead, we allow for a range of star-formation histories and dust extinctions for each cluster. Via MCMC modelling, we simultaneously fit, using stellar population synthesis models and DTD models, the integrated galaxy-light photometry in several bands, and the SN Ia numbers discovered in each cluster. With these more-realistic assumptions, we find a best-fitting DTD with power-law index $\alpha =-1.09_{-0.12}^{+0.15}$, and amplitude $R_1=0.41_{-0.10}^{+0.12}\times 10^{-12}\,{\rm yr}^{-1}\, {\rm M}_\odot ^{-1}$. We confirm a cluster-environment DTD with a larger amplitude than the field-galaxy DTD, by a factor ∼2–3 (at 3.8σ). Cluster and field DTDs have consistent slopes of α ≈ −1.1.

On the Thermally Pulsing Asymptotic Giant Branch Contribution to the Light of Nearby Disk Galaxies

Abstract The study of the luminosity contribution from thermally pulsing asymptotic giant branch (TP-AGB) stars to the stellar populations of galaxies is crucial to determine their physical parameters (e.g., stellar mass and age). We use a sample of 84 nearby disk galaxies to explore diverse stellar population synthesis models with different luminosity contributions from TP-AGB stars. We fit the models to optical and near-infrared (NIR) photometry, on a pixel-by-pixel basis. The statistics of the fits show a preference for a low-luminosity contribution (i.e., high mass-to-light ratio in the NIR) from TP-AGB stars. Nevertheless, for 30%–40% of the pixels in our sample a high-luminosity contribution (hence low mass-to-light ratio in the NIR) from TP-AGB stars is favored. According to our findings, the mean TP-AGB star luminosity contribution in nearby disk galaxies may vary with Hubble type. This may be a consequence of the variation of the TP-AGB mass-loss rate with metallicity, if metal-poor stars begin losing mass earlier than metal-rich stars, because of a pre-dust wind that precedes the dust-driven wind.

Powderday: Dust Radiative Transfer for Galaxy Simulations

Abstract We present powderday (available at https://github.com/dnarayanan/powderday), a flexible, fast, open-source dust radiative transfer package designed to interface with both idealized and cosmological galaxy formation simulations. powderday builds on fsps stellar population synthesis models, and hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, allowing significant flexibility in the assumed stellar physics and nebular line emission. The dust content follows either simple observationally motivated prescriptions (i.e., constant dust-to-metals ratios, or dust-to-gas ratios that vary with metallicity), direct modeling from galaxy formation simulations that include dust physics, as well as a novel approach that includes the dust content via learning-based algorithms from the simba cosmological galaxy formation simulation. Active galactic nuclei (AGNs) can additionally be included via a range of prescriptions. The output of these models are broadband (912 Å–1 mm) spectral energy distributions (SEDs), as well as filter-convolved monochromatic images. powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models and seamlessly interfaces with gizmo, arepo, gasoline, changa, and enzo. We demonstrate the capabilities of the code via three applications: a model for the star formation rate–infrared luminosity relation in galaxies (including the impact of AGNs), the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs, and the impact of galaxy inclination angle on dust attenuation laws.

Spatially-resolved stellar population properties of the M51–NGC 5195 system from multi-wavelength photometric data

Using multi-band photometric images of M51 and its companion NGC 5195 from ultraviolet to optical and infrared, we investigate spatially resolved stellar population properties of this interacting system with stellar population synthesis models. The observed infrared excess (IRX) is used to constrain dust extinction. Stellar mass is also inferred from the model fitting. By fitting observed spectral energy distributions (SEDs) with synthetical ones, we derive two-dimensional distributions of stellar age, metallicity, dust extinction and stellar mass. In M51, two grand-designed spiral arms extending from the bulge show young age, rich metallicity and abundant dust. The inter-arm regions are filled with older, metal-poorer and less dusty stellar populations. Except for the spiral arm extending from M51 into NGC 5195, the stellar population properties of NGC 5195 are quite featureless. NGC 5195 is much older than M51, and its core is very dusty with AV up to 1.67 mag and dense in stellar mass surface density. The close encounters might drive the dust in the spiral arm of M51 into the center of NGC 5195.

COOL-LAMPS. I. An Extraordinarily Bright Lensed Galaxy at Redshift 5.04*

Abstract We report the discovery of COOL J1241+2219, a strongly lensed galaxy at redshift z = 5.043 ± 0.002 with observed magnitude z AB = 20.47, lensed by a moderate-mass galaxy cluster at z = 1.001 ± 0.001. COOL J1241+2219 is the brightest lensed galaxy currently known at optical and near-infrared wavelengths at z ≳ 5; it is ∼5 times brighter than the prior record-holder lensed galaxy, and several magnitudes brighter than the brightest unlensed galaxies known at these redshifts. It was discovered as part of COOL-LAMPS, a collaboration initiated to find strongly lensed systems in recent public optical imaging data. We characterize the lensed galaxy, as well as the central galaxy of the lensing cluster using ground-based grizJH imaging and optical spectroscopy. We report model-based magnitudes, and derive stellar masses, dust content, metallicity, and star-formation rates via stellar-population synthesis modeling. Our lens mass modeling, based on ground-based imaging, implies a median source magnification of ∼30, which puts the stellar mass and star-formation rate (in the youngest age bin, closest to the epoch of observation) at logM * = 10.11 − 0.26 + 0.21 and SFR = 27 − 9 + 13 M ⊙ yr−1, respectively. We constrain a star-formation history for COOL J1241+2219 consistent with constant star formation across ∼1 Gyr of cosmic time, and that places this galaxy on the high-mass end of the star-forming main sequence. COOL J1241+2219 is two to four times more luminous than a galaxy with the characteristic UV luminosity at these redshifts. The UV continuum slope β = −2.2 ± 0.2 places this galaxy on the blue side of the observed distribution of galaxies at z = 5, although the lack of Lyα emission indicates dust sufficient to suppress this emission.

Sengi: A small, fast, interactive viewer for spectral outputs from stellar population synthesis models

We present Sengi, https://christopherlovell.github.io/sengi , an online tool for viewing the spectral outputs of stellar population synthesis (SPS) codes. Typical SPS codes require significant disk space or computing resources to produce spectra for simple stellar populations with arbitrary parameters. This makes it difficult to present their results in an interactive, web-friendly format. Sengi uses Non-negative Matrix Factorisation (NMF) and bilinear interpolation to estimate output spectra for arbitrary values of stellar age and metallicity. The reduced disk requirements and computational expense allows the result to be served as a client-based Javascript application. In this paper we present the method for generating grids of spectra, fitting those grids with NMF, bilinear interpolation across the fitted coefficients, and finally provide estimates of the prediction and interpolation errors.

Comparisons of Different Fitting Methods for the Physical Parameters of a Star Cluster Sample of M33 with Spectroscopy and Photometry

Star clusters are good tracers for formation and evolution of galaxies. We compared different fitting methods by using spectra (or by combining photometry) to determine the physical parameters. We choose a sample of 17 star clusters in M33, which previously lacked spectroscopic observations. The low-resolution spectra were taken with the Xinglong 2.16-m reflector of NAOC. The photometry used in the fitting includes $\rm u_{SC}$ and $\rm v_{SAGE}$ bands from the SAGE survey, as well as the published $UBVRI$ and $ugriz$ photometry. We firstly derived ages and metallicities with the {\sc ULySS} (Vazdekis et al. and {\sc pegase-hr}) SSP model and the Bruzual \& Charlot (2003) (BC03) stellar population synthesis models for the full-spectrum fitting. The fitting results of both the BC03 and {\sc ULySS} models seem consistent with those of previous works as well. Then we add the SAGE $\rm u_{SC}$ and $\rm v_{SAGE}$ photometry in the spectroscopic fitting with the BC03 models. It seems the results become much better, especially for the Padova 2000+Chabrier IMF set. Finally we add more photometry data, $UBVRI$ and $ugriz$, in the fitting and we found that the results do not improve significantly. Therefore, we conclude that the photometry is useful for improving the fitting results, especially for the blue bands ($\lambda <4000$ {\AA}), e.g., $\rm u_{SC}$ and $\rm v_{SAGE}$ band. At last, we discuss the "UV-excess" for the star clusters and we find five star clusters have UV-excess, based on the $GALEX$ FUV, NUV photometry.

The origin and impact of Wolf-Rayet-type mass loss

AbstractClassical Wolf-Rayet (WR) stars mark an important stage in the late evolution of massive stars. As hydrogen-poor massive stars, these objects have lost their outer layers, while still losing further mass through strong winds indicated by their prominent emission line spectra. Wolf-Rayet stars have been detected in a variety of different galaxies. Their strong winds are a major ingredient of stellar evolution and population synthesis models. Yet, a coherent theoretical picture of their strong mass-loss is only starting to emerge. In particular, the occurrence of WR stars as a function of metallicity (Z) is still far from being understood.To uncover the nature of the complex and dense winds of Wolf-Rayet stars, we employ a new generation of model atmospheres including a consistent solution of the wind hydrodynamics in an expanding non-LTE situation. With this technique, we can dissect the ingredients driving the wind and predict the resulting mass-loss for hydrogen-depleted massive stars. Our modelling efforts reveal a complex picture with strong, non-linear dependencies on the luminosity-to-mass ratio and Z with a steep, but not totally abrupt onset for WR-type winds in helium stars. With our findings, we provide a theoretical motivation for a population of helium stars at low Z, which cannot be detected via WR-type spectral features. Our study of massive He-star atmosphere models yields the very first mass-loss recipe derived from first principles in this regime. Implementing our first findings in stellar evolution models, we demonstrate how traditional approaches tend to overpredict WR-type mass loss in the young Universe.

Wolf–Rayet stars in the Antennae unveiled by MUSE

ABSTRACT We present the analysis of archival Very Large Telescope (VLT) Multi-Unit Spectroscopic Explorer (MUSE) observations of the interacting galaxies NGC 4038/39 (a.k.a. the Antennae) at a distance of 18.1 Mpc. Up to 38 young star-forming complexes with evident contribution from Wolf–Rayet (WR) stars are unveiled. We use publicly available templates of Galactic WR stars in conjunction with available photometric extinction measurements to quantify and classify the WR population in each star-forming region, on the basis of its nearly Solar oxygen abundance. The total estimated number of WR stars in the Antennae is 4053 ± 84, of which there are 2021 ± 60 WNL and 2032 ± 59 WC-types. Our analysis suggests a global WC to WN-type ratio of 1.01 ± 0.04, which is consistent with the predictions of the single star evolutionary scenario in the most recent bpass stellar population synthesis models.

The MOSDEF Survey: The First Direct Measurements of the Nebular Dust Attenuation Curve at High Redshift*

Abstract We use a sample of 532 star-forming galaxies at redshifts z ≃ 1.4–2.6 with deep rest-frame optical spectra from the MOSFIRE Deep Evolution Field (MOSDEF) survey to place the first constraints on the nebular attenuation curve at high redshift. Based on the first five low-order Balmer emission lines detected in the composite spectra of these galaxies (Hα through ), we derive a nebular attenuation curve that is similar in shape to that of the Galactic extinction curve, suggesting that the dust covering fraction and absorption/scattering properties along the lines of sight to massive stars at high redshift are similar to those of the average Milky Way sight line. The curve derived here implies nebular reddening values that are, on average, systematically larger than those derived for the stellar continuum. In the context of stellar population synthesis models that include the effects of stellar multiplicity, the difference in reddening of the nebular lines and stellar continuum may imply molecular cloud crossing timescales that are a factor of longer than those inferred for local molecular clouds, star formation rates that are constant or increasing with time such that newly formed and dustier OB associations always dominate the ionizing flux, and/or that the dust responsible for reddening the nebular emission may be associated with nonmolecular (i.e., ionized and neutral) phases of the interstellar medium. Our analysis points to a variety of investigations of the nebular attenuation curve that will be enabled with the next generation of ground- and space-based facilities.

No missing photons for reionization: moderate ionizing photon escape fractions from the FIRE-2 simulations

ABSTRACT We present the escape fraction of hydrogen ionizing photons (fesc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z ≥ 5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 haloes in Mvir ∼ 108–$10^{12}\, M_{\odot }$ (M* ∼ 104–$10^{10}\, M_{\odot }$) at z = 5–12. We find the sample average 〈fesc〉increases with halo mass for Mvir ∼ 108–$10^{9.5}\, M_{\odot }$, becomes nearly constant for 109.5–$10^{11}\, M_{\odot }$, and decreases at ${\gtrsim}10^{11}\, M_{\odot }$. Equivalently, 〈fesc〉 increases with stellar mass up to $M_{\ast }\sim 10^8\, M_{\odot }$ and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate 〈fesc〉 ∼ 0.2 for galaxies at $M_{\ast }\sim 10^8\, M_{\odot }$. Nearly half of the escaped ionizing photons come from stars 1–3 Myr old and the rest from stars 3–10 Myr old. Binaries only have a modest effect, boosting 〈fesc〉 by ∼25–35 per cent and the number of escaped photons by 60–80 per cent. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of 〈fesc〉 at the high-mass end is due to dust attenuation, while at the low-mass end, 〈fesc〉 decreases owing to inefficient star formation and hence feedback. At fixed mass, 〈fesc〉 tends to increase with redshift. Although the absolute 〈fesc〉does not fully converge with resolution in our simulations, the mass- and redshift-dependence of 〈fesc〉 is likely robust. Our simulations produce sufficient ionizing photons for cosmic reionization.

An upgraded interpolator of the radial basis function network for spectral calculation based on empirical stellar spectral library

Stellar population synthesis is an important method in galaxy and star-cluster studies. In stellar population synthesis models, a stellar spectral library is necessary for the integrated spectra of the stellar population. Usually, the stellar spectral library is utilized for the transformation between the stellar atmospheric parameters and stellar spectrum. The empirical stellar spectral library has distinct advantages over the theoretical library. However, for the empirical spectral library, the distribution of stars is irregular in the stellar atmospheric parameter space, which makes the traditional interpolator difficult to get accurate results. In this work, we will provide an improved radial basis function interpolator which is implemented to obtain the interpolated stellar spectra based on the empirical stellar spectral library. For this interpolator, we use the relation between the standard variance σ in the Gaussian radial basis function and the density distribution of stars in the stellar atmospheric parameter space to give the prior constraint on this σ. Moreover, we also consider the anisotropic radial basis function by the advantage of the local dispersion of stars in the stellar atmospheric parameter space. Furthermore, we use the empirical stellar spectral library MILES to test this interpolator. On the whole, the interpolator has a good performance except for the edge of the low-temperature region. At last, we compare this interpolator with our work in 2018, and the interpolation result shows an obvious improvement. Users can apply this interpolator to retrieve the interpolated spectra based on the stellar spectral library quickly and easily.

Delayed Photons from Binary Evolution Help Reionize the Universe

Abstract High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (f esc) of hydrogen-ionizing photons often assume stellar radiation based on single-stellar population synthesis models. However, strong evidence suggests the binary fraction of massive stars is ≳70%. Moreover, simulations so far have yielded values of f esc falling only on the lower end of the ∼10%–20% range, the amount presumed necessary to reionize the universe. Analyzing a high-resolution (4 pc) cosmological radiation-hydrodynamic simulation, we study how f esc changes when we include two different products of binary stellar evolution—stars stripped of their hydrogen envelopes and massive blue stragglers. Both produce significant amounts of ionizing photons 10–200 Myr after each starburst. We find the relative importance of these photons to be amplified with respect to escaped ionizing photons, because peaks in star formation rates (SFRs) and f esc are often out of phase by this 10–200 Myr. Additionally, low-mass, bursty galaxies emit Lyman continuum radiation primarily from binary products when SFRs are low. Observations of these galaxies by the James Webb Space Telescope could provide crucial information on the evolution of binary stars as a function of redshift. Overall, including stripped stars and massive blue stragglers increases our photon-weighted mean escape fraction ( ) by ∼13% and ∼10%, respectively, resulting in . Our results emphasize that using updated stellar population synthesis models with binary stellar evolution provides a more sound physical basis for stellar reionization.

An uncertainty principle for star formation – III. The characteristic emission time-scales of star formation rate tracers

ABSTRACT We recently presented a new statistical method to constrain the physics of star formation and feedback on the cloud scale by reconstructing the underlying evolutionary timeline. However, by itself this new method only recovers the relative durations of different evolutionary phases. To enable observational applications, it therefore requires knowledge of an absolute ‘reference time-scale’ to convert relative time-scales into absolute values. The logical choice for this reference time-scale is the duration over which the star formation rate (SFR) tracer is visible because it can be characterized using stellar population synthesis (SPS) models. In this paper, we calibrate this reference time-scale using synthetic emission maps of several SFR tracers, generated by combining the output from a hydrodynamical disc galaxy simulation with the SPS model slug2. We apply our statistical method to obtain self-consistent measurements of each tracer’s reference time-scale. These include H α and 12 ultraviolet (UV) filters (from GALEX, Swift, and HST), which cover a wavelength range 150–350 nm. At solar metallicity, the measured reference time-scales of H α are ${4.32^{+0.09}_{-0.23}}$ Myr with continuum subtraction, and 6–16 Myr without, where the time-scale increases with filter width. For the UV filters we find 17–33 Myr, nearly monotonically increasing with wavelength. The characteristic time-scale decreases towards higher metallicities, as well as to lower star formation rate surface densities, owing to stellar initial mass function sampling effects. We provide fitting functions for the reference time-scale as a function of metallicity, filter width, or wavelength, to enable observational applications of our statistical method across a wide variety of galaxies.

Setting the Stage for Cosmic Chronometers. II. Impact of Stellar Population Synthesis Models Systematics and Full Covariance Matrix

Abstract The evolution of differential ages of passive galaxies at different redshifts (cosmic chronometers) has been proved to be a method potentially able to constrain the Hubble parameter in a cosmology-independent way, but the systematic uncertainties must be carefully evaluated. In this paper, we compute the contribution to the full covariance matrix of systematic uncertainties due to the choice of initial mass function, stellar library, and metallicity, exploring a variety of stellar population synthesis models. Through simulations in the redshift range 0 &lt; z &lt; 1.5, we find that the choice of the stellar population synthesis model dominates the total error budget on H(z), with contributions at a level of ∼4.5%, discarding the most discordant model. The contribution due to the choice of initial mass function is &lt;0.5%, while that due to the stellar library is ∼6.6%, on average. We also assess the impact of an uncertainty in the stellar metallicity determination, finding that an error of ∼10% (5%) on the stellar metallicity propagates to a 9% (4%) error on H(z). These results are used to provide the combined contribution of these systematic effects on the error budget. For current H(z) measurements, where the uncertainties due to metallicity and star formation history were already included, we show that, using the more modern stellar libraries, the additional systematic uncertainty is between 5.4% (at z = 0.2) and 2.3% (at z = 1.5). To reach the goal of keeping the systematic error budget below the 1% level, we discuss the efforts needed to obtain higher resolution and signal-to-noise spectra and improvements in the modeling of stellar population synthesis.

SPECULATOR: Emulating Stellar Population Synthesis for Fast and Accurate Galaxy Spectra and Photometry

We present SPECULATOR - a fast, accurate, and flexible framework for emulating stellar population synthesis (SPS) models for predicting galaxy spectra and photometry. For emulating spectra, we use principal component analysis to construct a set of basis functions, and neural networks to learn the basis coefficients as a function of the SPS model parameters. For photometry, we parameterize the magnitudes (for the filters of interest) as a function of SPS parameters by a neural network. The resulting emulators are able to predict spectra and photometry under both simple and complicated SPS model parameterizations to percent-level accuracy, giving a factor of $10^3$-$10^4$ speed up over direct SPS computation. They have readily-computable derivatives, making them amenable to gradient-based inference and optimization methods. The emulators are also straightforward to call from a GPU, giving an additional order-of-magnitude speed-up. Rapid SPS computations delivered by emulation offers a massive reduction in the computational resources required to infer the physical properties of galaxies from observed spectra or photometry and simulate galaxy populations under SPS models, whilst maintaining the accuracy required for a range of applications.

Beyond Gaia: Asteroseismic Distances of M Giants Using Ground-based Transient Surveys

Abstract Evolved stars near the tip of the red giant branch show solar-like oscillations with periods spanning hours to months and amplitudes ranging from ∼1 mmag to ∼100 mmag. The systematic detection of the resulting photometric variations with ground-based telescopes would enable the application of asteroseismology to a much larger and more distant sample of stars than is currently accessible with space-based telescopes such as Kepler or the ongoing Transiting Exoplanet Survey Satellite mission. We present an asteroseismic analysis of 493 M giants using data from two ground-based surveys: the Asteroid Terrestrial-impact Last Alert System (ATLAS) and the All-Sky Automated Survey for Supernovae (ASAS-SN). By comparing the extracted frequencies with constraints from Kepler, the Sloan Digital Sky Survey Apache Point Observatory Galaxy Evolution Experiment, and Gaia we demonstrate that ground-based transient surveys allow accurate distance measurements to oscillating M giants with a precision of ∼15%. Using stellar population synthesis models we predict that ATLAS and ASAS-SN can provide asteroseismic distances to ∼2 × 106 galactic M giants out to typical distances of 20–50 kpc, vastly improving the reach of Gaia and providing critical constraints for Galactic archeology and galactic dynamics.

Bridging between the Integrated and Resolved Main Sequence of Star Formation

Abstract The position of galaxies on the stellar mass, star formation rate (SFR) plane with respect to the star-forming main sequence at each redshift is a convenient way to infer where the galaxy is in its evolution compared to the rest of the population. We use Hubble Space Telescope high-resolution images in the GOODS-S field from the the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) and fit multiwavelength lights in resolution elements of galaxies with stellar population synthesis models. We then construct resolved kpc-scale stellar mass, SFR surface density curves for galaxies at z ∼ 1. Fitting these resolved main sequence curves with Schechter functions, we parameterize and explain the multiwavelength structure of galaxies with three variables: ϕ*, α, and M*. For quenched galaxies below the main sequence, we find an average high-mass slope (α) of the resolved main sequence curves to be ∼−0.4. The scatter of this slope is higher among the lower mass star-forming galaxies and those above the main sequence compared to quenched galaxies, due to lack of an evolved bulge. Our findings agree well with an inside-out quenching of star formation. We find that the knee of the Schechter fits (M*) for galaxies below the main sequence occurs at lower stellar mass surface densities compared to star-forming galaxies, which hints at how far quenching has proceeded outward.

A Comparison of Rotating and Binary Stellar Evolution Models: Effects on Massive Star Populations

Abstract Both rotation and interactions with binary companions can significantly affect massive star evolution, altering interior and surface abundances, mass loss rates and mechanisms, observed temperatures and luminosities, and their ultimate core-collapse fates. The Geneva and the Binary Population and Spectral Synthesis (BPASS) stellar evolution codes include detailed treatments of the effects of rotation and binary evolution, respectively, and can illustrate the impact of these phenomena on massive stars and stellar populations. However, a direct comparison of these two widely used codes is vital if we hope to use their predictions for interpreting observations. In particular, rotating and binary models will predict different young stellar populations, impacting the outputs of stellar population synthesis (SPS) and the resulting interpretation of large samples of massive stars based on commonly used tools such as star count ratios. Here we compare the Geneva and BPASS evolutionary models, using an interpolated SPS scheme introduced in our previous work and a novel Bayesian framework to present the first in-depth direct comparison of massive stellar populations produced from single, rotating, and binary nonrotating evolution models. We calculate both models’ predicted values of star count ratios and compare the results to observations of massive stars in Westerlund 1, h + χ Persei, and both Magellanic Clouds. We also consider the limitations of both the observations and the models, and how to quantitatively include observational completeness limits in SPS models. We demonstrate that the methods presented here, when combined with robust stellar evolutionary models, offer a potential means of estimating the physical properties of massive stars in large stellar populations.

Dynamics of DiskMass Survey galaxies in refracted gravity

We aim to verify whether refracted gravity (RG) is capable of describing the dynamics of disk galaxies without resorting to the presence of dark matter. RG is a classical theory of gravity in which the standard Poisson equation is modified with the introduction of the gravitational permittivity, which is a universal monotonic function of the local mass density. We used the rotation curves and the radial profiles of the stellar velocity dispersion perpendicular to the galactic disks of 30 disk galaxies from the DiskMass Survey (DMS) to determine the gravitational permittivity. RG describes the rotation curves and the vertical velocity dispersions by requiring galaxy mass-to-light ratios that are in agreement with stellar population synthesis models, and disk thicknesses that are in agreement with observations, once observational biases are taken into account. Our results rely on setting the three free parameters of the gravitational permittivity for each individual galaxy. However, we show that the differences of these parameters from galaxy to galaxy can, in principle, be ascribed to statistical fluctuations. We adopted an approximate procedure to estimate a single set of parameters that may properly describe the kinematics of the entire sample and suggest that the gravitational permittivity is indeed a universal function. Finally, we showed that the RG models of the individual rotation curves can only partly describe the radial acceleration relation (RAR) between the observed centripetal acceleration derived from the rotation curve and the Newtonian gravitational acceleration originating from the baryonic mass distribution. Evidently, the RG models underestimate the observed accelerations by 0.1 to 0.3 dex at low Newtonian accelerations. An additional problem that ought to be considered is the strong correlation, at much more than 5σ, between the residuals of the RAR models and three radially-dependent properties of the galaxies, whereas the DMS data show a considerably less significant correlation, at more than 4σ, for only two of these quantities. These correlations might be the source of the non-null intrinsic scatter of the RG models: this non-null scatter is at odds with the observed intrinsic scatter of other galaxy samples different from DMS, which is consistent with zero. Further investigations are required to assess whether these discrepancies in the RAR originate from the DMS sample, which might not be ideal for deriving the RAR, or whether they are genuine failures of RG.