Population Synthesis Models Research Articles

Surface chemical composition of single WNh stars

Context. Wolf–Rayet (WR) stars of the WNh category contain a significant fraction of hydrogen at their surface. They can be hydrogen-burning, very massive stars or stars in a post-main sequence phase of evolution. Also, WNh stars are sometimes not included in population synthesis models. Aims. We aim to better characterise the properties of single WNh stars in the Galaxy and the Magellanic Clouds. In particular, we want to constrain their surface chemistry beyond the hydrogen content by determining the helium, carbon, and nitrogen surface abundances. Methods. We perform a spectroscopic analysis of 22 single WNh stars. We fit their ultraviolet and/or optical spectra using synthetic spectra computed with the code CMFGEN. We determine the main stellar parameters (temperature, luminosity, mass-loss rates) and the surface H, He, C, and N mass fractions. We investigate the ability of current evolutionary models to reproduce all parameters at the same time. Results. We find that all WNh stars show the signatures of CNO-cycle material at their surface: they are carbon-depleted and nitrogen-rich. A clear trend of higher nitrogen content at higher metallicity is observed, as expected. The amount of hydrogen (X) varies significantly from one star to another, independently of luminosity. Values of X larger than 0.4 are not exceptional. The majority of Galactic WNh stars can be explained by evolutionary models, provided sufficient fine-tuning of the input parameters of evolutionary calculations. At lower metallicity, most stars escape predictions from evolutionary models. This has been noted in the literature but constraints on the surface nitrogen content exacerbate this severe issue. Conclusions. Our study highlights the need to refine the treatment of WR stars in both stellar evolution and population synthesis models.

Stellar Properties of Observed Stars Stripped in Binaries in the Magellanic Clouds

Massive stars (∼8–25 M ⊙) stripped of their hydrogen-rich envelopes via binary interaction are thought to be the main progenitors for merging neutron stars and stripped-envelope supernovae. We recently presented the discovery of the first set of such stripped stars in a companion paper. Here, we fit the spectra of 10 stars with new atmosphere models in order to constrain their stellar properties precisely. We find that the stellar properties align well with the theoretical expectations from binary evolution models for helium-core burning envelope-stripped stars. The fits confirm that the stars have high effective temperatures (T eff ∼ 50–100 kK), high surface gravities ( 5), and hydrogen-poor/helium-rich surfaces (X H,surf ∼ 0–0.4) while showing for the first time a range of bolometric luminosities (103–105 L ⊙), small radii (∼0.5–1 R ⊙), and low Eddington factors (Γ e ∼ 0.006–0.4). Using these properties, we derive intermediate current masses (∼1–8 M ⊙), which suggest that their progenitors were massive stars (∼5–25 M ⊙) and that a subset will reach core-collapse, leaving behind neutron stars or black holes. Using the model fits, we also estimate the emission rates of ionizing photons for these stars, which agree well with previous model expectations. Further, by computing models for a range of mass-loss rates, we find that the stellar winds are weaker than predicted by any existing scheme ( M ⊙ yr−1). The properties of this first sample of intermediate-mass helium stars suggest they both contain progenitors of type Ib and IIb supernovae, and provide important benchmarks for binary evolution and population synthesis models.

The impact of convective criteria on the properties of massive stars

Context. Libraries of stellar models computed with either the Ledoux or the Schwarzschild criterion to determine the sizes of convective regions are available in the literature. It is still not clear, however, which of these two criteria should be used, although many works have been devoted to that question in the past. Aims. In the framework of the evolution of single rotating stars, we study the differences between models computed with Ledoux and Schwarzschild criteria on the internal structure, evolutionary track in the Hertzsprung–Russell diagram (HRD), lifetimes, evolution of the surface abundances and velocities, and masses of the He and CO cores. We investigate the consequences on the nature of the supernova (SN) progenitors and the type of SN events, as well as on the stellar yields of light elements. We also study the impact on the outputs of population synthesis models. Methods. Models with initial masses between 7 and 120 M⊙ at solar metallicity (Z = 0.014) and with an initial rotation equal to 0 or 0.4 times the critical velocity at the zero-age main sequence were computed with either the Schwarzschild or the Ledoux criterion until the end of the C-burning phase. Results. Models with initial masses between 15 and 32 M⊙ computed with the Schwarzschild criterion show larger intermediate convective zones attached to the H-burning shell than models computed with the Ledoux criterion. Their CO cores and outer convective zones in the red supergiant (RSG) phase are also smaller. This impacts many outputs of stars during the core He-burning phase. Schwarzschild models have smaller CO cores and outer convective zones in the RSG phase, and their blue-to-red supergiant ratio is much higher than for Ledoux models. They also produce longer crossings of the Hertzsprung gap and favour blue loops. The upper luminosity of RSGs is little affected by the change in the convective criterion. The maximum luminosity of RSG progenitors for type II-P SN events is lowered from 5.2 to 4.95 when the Ledoux criterion is used instead of the Schwarzschild criterion in non-rotating models. The Schwarzschild criterion predicts longer-lasting, less nitrogen-enriched, and faster-rotating Cepheids. Rotational mixing tends to decrease the differences between Schwarzschild and Ledoux models. Conclusions. The results of this paper can be used as first guidelines to set up observational programs that may help to distinguish between these two model families.

HST FUV Spectroscopy of Super Star Cluster A in the Green Pea Analog Mrk 71: Revealing the Presence of Very Massive Stars

Mrk 71 is a low-metallicity (Z = 0.16 Z ☉) starburst region in the local dwarf galaxy NGC 2366, hosting two super star clusters (SSCs A and B), and it is recognized as a Green Pea (GP) analog with SSC A responsible for the GP properties. We present STIS and FOS far-ultraviolet (FUV) spectra of the embedded SSC Mrk 71-A obtained with the Hubble Space Telescope. The STIS FUV spectrum shows the characteristic features of very massive stars (VMS; masses >100 M ⊙) and we derive an age of 1 ± 1 Myr by comparison with the Charlot & Bruzual suite of spectral population synthesis models with upper mass limits of 300 and 600 M ⊙. We compare the STIS spectrum with all known SSC spectra exhibiting VMS signatures: NGC 5253-5, R136a, NGC 3125-A1, and the z = 2.37 Sunburst cluster. We find that the cluster mass-loss rates and wind velocities, as characterized by the C iv P Cygni profiles and the He ii emission line strengths, are very similar over Z = 0.16–0.4 Z☉. This agrees with predictions that the optically thick winds of VMS will be enhanced near the Eddington limit and show little metallicity dependence. We find very strong damped Lyα absorption with N(H i) =1022.2 cm−2 associated with Mrk 71-A. We discuss the natal environment of this young SSC in terms of radiatively driven winds, catastrophic cooling, and recent models where the cluster is surrounded by highly pressurized clouds with large neutral columns.

Mapping the core of the Tarantula Nebula with VLT-MUSE – III. A template for metal-poor starburst regions in the visual and far-ultraviolet

ABSTRACT We present the integrated VLT-MUSE spectrum of the central 2 × 2 arcmin2 (30 × 30 pc2) of NGC 2070, the dominant giant H ii region of the Tarantula Nebula in the Large Magellanic Cloud, together with an empirical far-ultraviolet spectrum constructed via LMC template stars from the ULLYSES survey and Hubble Tarantula Treasury Project UV photometry. NGC 2070 provides a unique opportunity to compare results from individual stellar populations (e.g. VLT FLAMES Tarantula Survey) in a metal-poor starburst region to the integrated results from the population synthesis tools Starburst99, Charlot & Bruzual, and BPASS. The metallicity of NGC 2070 inferred from standard nebular strong line calibrations is ∼0.4 ± 0.1 dex lower than obtained from direct methods. The Hα inferred age of 4.2 Myr from Starburst99 is close to the median age of OB stars within the region, although individual stars span a broad range of 1–7 Myr. The inferred stellar mass is close to that obtained for the rich star cluster R136 within NGC 2070, although this contributes only 21 per cent to the integrated far-UV continuum. He ii λ1640 emission is dominated by classical WR stars and main sequence WNh + Of/WN stars. Around 18 per cent of the NGC 2070 far UV continuum flux arises from very massive stars with ≥100 M⊙, including several very luminous Of supergiants. None of the predicted population synthesis models at low metallicities are able to successfully reproduce the far-UV spectrum of NGC 2070. We attribute issues to the treatment of mass-loss in very massive stars, the lack of contemporary empirical metal-poor templates, plus WR stars produced via binary evolution.

Considering the Single and Binary Origins of the Type IIP SN 2017eaw

Current population synthesis modeling suggests that 30%–50% of Type II supernovae originate from binary progenitors; however, the identification of a binary progenitor is challenging. One indicator of a binary progenitor is that the surrounding stellar population is too old to contain a massive single star. Measurements of the progenitor mass of SN 2017eaw are starkly divided between observations made temporally close to core collapse, which show a progenitor mass of 13–15 M ⊙ (final helium-core mass MHe,core=4.4–6.0M⊙ —which is a more informative property than initial mass) and those from the stellar population surrounding the SN, which find M ≤ 10.8 M ⊙ ( MHe,core≤3.4M⊙ ). In this paper, we reanalyze the surrounding stellar population with improved astrometry and photometry, finding a median age of 16.8−1.0+3.2 Myr for all stars younger than 50 Myr ( MHe,core=4.7M⊙ ) and 85.9−6.5+3.2 Myr for stars younger than 150 Myr. 16.8 Myr is now consistent with the helium-core mass range derived from the temporally near-explosion observations for single stars. Applying the combined constraints to population synthesis models, we determine that the probability of the progenitor of SN 2017eaw being an initially single star is 65% compared to 35% for prior binary interaction. 85.9 Myr is inconsistent with any formation scenarios. We demonstrate that combining progenitor age constraints with helium-core mass estimates from red supergiant SED modeling, late-time spectra, and indirectly from light-curve modeling can help to differentiate single and binary progenitor scenarios and provide a framework for the application of this technique to future observations.

Linking the Metallicity Enrichment History to the Star Formation History: An SFH-regulated Chemical Evolution Model and Its Implications for the Gas Cycling Process

The metallicity enrichment history (MEH) of a galaxy is determined by its star formation history (SFH) and the gas cycling process. In this paper, we construct a chemical evolution model that is regulated by the SFH of the system. In this SFH-regulated model, the evolution of all other variables, including the MEH, can be determined by the SFH. We test this model on six locally isolated dwarf galaxies covering three dwarf types that were observed by the Local Cosmology from Isolated Dwarfs (LCID) project. The SFHs and MEHs of these LCID galaxies are relatively reliable because they have been measured from deep color–magnitude diagrams that reach the main-sequence turnoff stars with good photometric accuracy. With simple assumptions of the star formation law and the mass-dependent outflows, our SFH-regulated model successfully reproduces the MEHs of all six LCID galaxies from their SFHs, with only one free parameter, the wind efficiency η ∼ 1.0, for all six galaxies. This model provides a physically motivated link that directly connects the SFH and MEH of a galaxy, which will be useful to accommodate into the state-of-the-art stellar population synthesis models to help relieve the nuisance of the heavy degeneracy between the ages and metallicities of the stellar populations.

Relation between spectral indices and binary fractions in globular clusters

Context. We study the relation between the known binary fraction and spectral absorption feature index to judge whether (and potentially which) spectral absorption feature indices are suitable for determining the binary fraction. Aims. The determination of the binary fraction is important in studies of binary star formation, evolutionary population synthesis models, and other works. The number of binary stars is difficult to determine for nearly all stellar systems because the individual stars are need to be resolved photometrically or spectroscopically. By comparison, their integrated spectra or spectral absorption feature indices are relatively easy to obtain. Methods. We used Galactic globular clusters (GCs) as our sample since they have relatively accurate binary fraction measurements and spectroscopic observations along the radial direction. When studying the relation between binary fractions and the spectral absorption feature index, we used three types of binary fractions: one with a mass ratio of q > 0.5 (f(q > 0.5)) and two with a total binary fraction (the methods of counting (f(tot)mf) and fitting (f(tot)mc)), calculated and obtained the equivalent widths or magnitudes of 46 spectral absorption feature indices at three spectral resolutions (FWHMLick/IDS, 5, and 15 Å). Since the regions for the binary fraction measurements (photometric) and the spectroscopic observations are different, we used the method of constructing the radial binary-fraction profile to get the binary fractions corresponding to the regions in the spectroscopic observations. The construction of the radial binary-fraction profile was obtained by using the python curve_fit module to fit the measured and analytic binary fraction values. The analytic value was expressed by taking advantage of the King surface-density profile and the assumed forms with respect to the radial binary-fraction profile (linear, quadratic, exponential, and Gaussian). Results. We find that the low-resolution (15 Å) spectrum is not suitable for this study and the binary fraction type would affect the results: f(q > 0.5) and f(tot)mc exhibit better correlations with the spectral absorption feature index than f(tot)mf and the difference in metallicity would significantly affect the above relationship. Finally, to eliminate the effects of metallicity, age, and dynamical evolution, we only used those GCs with multiple spectra observed among different regions. We find that OIII-1, OIII-2, HγF, HδF, HγA, HδA, Hβ, Ca4455, C24668, and TiO1 indices have strong correlations with binary fraction. The two OIII indices are the most sensitive to the binary fraction, followed by four Balmer indices – the two narrower central bandpass Balmer indices (∼20 Å, F-definition) are more sensitive than the wider two (∼40 Å, A-definition) and, lastly, the Ca4455, C24668, and TiO1 indices. Using the binary fraction-sensitive spectral absorption feature indices in combination with the age- and metallicity-sensitive spectral absorption feature indices, we can determine the ages or metallicities first. Then we can go on to obtain the binary fractions for those stellar systems in which the individual stars cannot be resolved (e.g., dense or distant stellar systems). Furthermore, we are then able to carry out further studies of binary star formation and improve evolutionary population synthesis models.

A resolved study of the inner regions of nearby galaxies with an excess of young massive stars: missing link in the AGN–starburst connection?

ABSTRACT We have selected galaxies with very high levels of H α emission (EQW(H α) >700 Å.) in their central regions from the final data release of the MaNGA survey. Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range 9.5 < log M*/(M⊙) < 11.5. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populations of young massive stars, finding that selection in the 2D space of extinction-corrected H α EQW and [S iii]/[S ii] line ratio produces the best results. By comparing stacked spectra covering these regions with stacked spectra covering normal starburst regions with 100 Å<EQW(H α) <200 Å, we obtain the following main results: (1) Clear signatures of excess Wolf–Rayet stars are found in half of the H α excess regions, (2) galaxy regions containing excess Wolf–Rayet stars are more often associated with the presence of high-ionization emission lines characteristic of accreting black holes. Excess [Ne iii] is detected in 4 out of 8 of the WR regions and there are tentative [Fe x] detections in 2 galaxies. (3) Regions of the galaxy with excess Wolf–Rayet stars are located where the interstellar medium has larger ionized gas turbulent velocities and higher neutral gas overdensities. We make a first attempt to constrain changes in the high-mass end of the stellar initial mass function (IMF) using the HR-pyPopStar evolutionary population synthesis models that include high-wavelength-resolution theoretical atmosphere libraries for Wolf–Rayet stars.

Inferring the presence of very massive stars in local star-forming regions

Context. Very massive stars (VMS) have masses in excess of 100 M⊙ and are rare. However, owing to their powerful winds, very high luminosity, and efficient nucleosynthesis, they are key players of star-forming regions. In particular, their strong ionizing fluxes impact the surrounding interstellar medium. Aims. We aim at detecting VMS in local star-forming region from the imprint they leave on the integrated UV and optical light. Methods. We analyzed a sample of 27 star-forming regions and galaxies in the local Universe. We selected sources with a metallicity close to 12+log(O/H) = 8.3, which is typical of the Large Magellanic Cloud. We defined empirical criteria to distinguish sources dominated by VMS and Wolf-Rayet stars (WR), using template spectra of VMS- and WR-dominated regions. We subsequently built population synthesis models with an updated treatment of VMS. We investigated the successes and failures of these new models in accounting for the UV-optical spectroscopy of our sample sources. Results. We show that the UV range alone is not sufficient to distinguish between VMS- and WR-dominated sources because their spectra are almost identical in this range. The region of the WR bumps in the optical breaks the degeneracy. In particular, the morphology of the blue bump at 4640–4686 Å is a key diagnostic. Beyond the prototypical R136 region, which contains VMS, we identify two galaxies showing clear signatures of VMS. In two other galaxies or regions the presence of VMS can be suspected, as already discussed in the literature. The stellar population is clearly dominated by WR stars in seven other sources. The most recent BPASS population synthesis models can neither account for the strong He II 1640 emission, nor for the shape of the blue bump in VMS- and WR-dominated sources. Our models that include VMS more realistically reproduce the UV-optical spectra of VMS-dominated sources. Conclusions. We conclude that VMS are present in some local star-forming regions, but that separating them from WR-dominated populations requires optical spectroscopy with a high signal-to-noise ratio. A high equivalent width of He II 1640 is not a sufficient condition for identifying VMS. Populations synthesis models need to take VMS into account by incorporating not only evolutionary tracks, but also dedicated spectral libraries. Finally, we stress that the treatment of WR stars needs to be improved as well.

Abundances and Transients from Neutron Star–White Dwarf Mergers

We systematically investigate the mergers of neutron star–white dwarf binaries from beginning to end, with a focus on the properties of the inflows and outflows in accretion disks and their electromagnetic emissions. Using population synthesis models, we determine a subset of these binaries in which the white dwarf companion undergoes unstable mass transfer and complete tidal disruption, forming a large accretion disk around the neutron star. The material evolves according to a one-dimensional advection-dominated accretion-disk model with nuclear burning, neutrino emissions, and disk-surface wind ejection. The extreme dynamics of the entire process have proven difficult to analyze, and thus currently, the properties are poorly understood. The outflows from the mergers are iron- and nickel-rich, giving rise to optical and infrared emissions powered by the decay of the radioactive iron-type isotopes, calculated via the SuperNu light-curve code. We find these systems capable of powering bright, yet short-lived, optical transients with the potential to power gamma-ray bursts.

Initial mass function variability from the integrated light of diverse stellar systems

ABSTRACT We present a uniform analysis of the stellar initial mass function (IMF) from integrated light spectroscopy of 15 compact stellar systems (11 globular clusters in M31 and 4 ultra compact dwarfs in the Virgo cluster, UCDs) and two brightest Coma cluster galaxies (BCGs), covering a wide range of metallicities (−1.7 < [Fe/H] < 0.01) and velocity dispersions (7.4 km s−1 <σ < 275 km s−1). The S/N ∼100 Å−1 Keck LRIS spectra are fitted over the range 4000 < λ/Å < 10 000 with flexible full-spectrum stellar population synthesis models. We use the models to fit simultaneously for ages, metallicities, and individual elemental abundances of the population, allowing us to decouple abundance variations from variations in IMF slope. We show that compact stellar systems do not follow the same trends with physical parameters that have been found for early-type galaxies. Most globular clusters in our sample have an IMF consistent with that of the Milky Way, over a wide range of [Fe/H] and [Mg/Fe]. There is more diversity among the UCDs, with some showing evidence for a bottom-heavy IMF, but with no clear correlation with metallicity, abundance, or velocity dispersion. The two Coma BCGs have similar velocity dispersion and metallicity, but we find the IMF of NGC 4874 is consistent with that of the Milky Way while NGC 4889 presents evidence for a significantly bottom-heavy IMF. For this sample, the IMF appears to vary between objects in a way that is not explained by a single metallicity-dependent prescription.

Astrophysical uncertainties in the gravitational-wave background from stellar-mass compact binary mergers

ABSTRACT We investigate the stochastic gravitational-wave background (SGWB) produced by merging binary black holes (BBHs) and binary neutron stars (BNSs) in the frequency ranges of Laser Interferometer Gravitational-Wave Observatory (LIGO)/Virgo/Kagra and Laser Interferometer Space Antenna (LISA). We develop three analytical models, which are calibrated to the measured local merger rates, and complement them with three population synthesis models based on the cosmic code. We discuss the uncertainties, focusing on the impact of the BBH mass distribution, the effect of the metallicity of the progenitor stars, and the time delay distribution between star formation and compact binary merger. We also explore the effect of uncertainties in binary stellar evolution on the background. For BBHs, our analytical models predict ΩGW in the range [4 × 10−10 to 1 × 10−9] (25 Hz) and [1 × 10−12 to 4 × 10−12] (3 mHz), and between [2 × 10−10 to 2 × 10−9] (25 Hz) and [7 × 10−13 to 7 × 10−12] (3 mHz) for our population synthesis models. This background is unlikely to be detected during the LIGO/Virgo/Kagra O4 run, but could be detectable with LISA. We predict about 10 BBH and no BNS mergers that could be individually detectable by LISA for a period of observation of 4 yr. Our study provides new insights into the population of compact binaries and the main sources of uncertainty in the astrophysical SGWB.

The galaxy cluster AC114 – II. Stellar populations and the mass–metallicity relation

ABSTRACT We investigate the mass–metallicity relation (MZR)for galaxies in the Abell cluster AC114 from 7 h of VIMOS/MR data collected at the European Southern Observatory-Very Large Telescope in 2009. The dynamical analysis completed in our previous paper allowed us to select cluster members, whose spectra are here analysed with stellar population synthesis models. Active and passive galaxies are identified based on the presence/absence of the [O ii] λ3727, [O iii] λλ4959, 5007 and/or Hβ emission lines, depending on the galaxy redshift. We find that active galaxies have lower average masses than passive ones, and have lower average metallicities. The MZR of the cluster is found to be steeper than that for galaxies in the local universe.

Full spectrum fitting with photometry in ppxf: stellar population versus dynamical masses, non-parametric star formation history and metallicity for 3200 LEGA-C galaxies at redshift z ≈ 0.8

ABSTRACT I introduce some improvements to the ppxf method, which measures the stellar and gas kinematics, star formation history (SFH) and chemical composition of galaxies. I describe the new optimization algorithm that ppxf uses and the changes I made to fit both spectra and photometry simultaneously. I apply the updated ppxf method to a sample of 3200 galaxies at redshift 0.6 < z < 1 (median z = 0.76, stellar mass $M_\ast \gtrsim 3\times 10^{10}$ M⊙), using spectroscopy from the LEGA-C survey (DR3) and 28-bands photometry from two different sources. I compare the masses from new JAM dynamical models with the ppxf stellar population M* and show the latter are more reliable than previous estimates. I use three different stellar population synthesis (SPS) models in ppxf and both photometric sources. I confirm the main trend of the galaxies’ global ages and metallicity [M/H] with stellar velocity dispersion σ* (or central density), but I also find that [M/H] depends on age at fixed σ*. The SFHs reveal a sharp transition from star formation to quenching for galaxies with $\lg (\sigma _\ast /\mathrm{km}\, s^{-1})\gtrsim 2.3$ ($\sigma _\ast \gtrsim 200$$\mathrm{km}\, s^{-1}$), or average mass density within 1 kpc $\lg (\Sigma _1^{\rm JAM}/\mathrm{\mathrm{M}_{\odot }kpc^{-2}})\gtrsim 9.9$ ($\Sigma _1^{\rm JAM}\gtrsim 7.9\times 10^9\, \mathrm{\mathrm{M}_{\odot }\ kpc^{-2}}$), or with $[M/H]\gtrsim -0.1$, or with Sersic index $\lg n_{\rm Ser}\gtrsim 0.5$ ($n_{\rm Ser}\gtrsim 3.2$). However, the transition is smoother as a function of M*. These results are consistent for two SPS models and both photometric sources, but they differ significantly from the third SPS model, which demonstrates the importance of comparing model assumptions.

An integration modeling framework for individual-scale daily mobility estimation

This paper presents a novel microsimulation framework (BayABM) that integrates population synthesis, Bayesian Networks (BNs), and activity-based models (activity-BM) to model individual heterogeneity in mobility. The population synthesis simulates a synthetic population with socioeconomic and household-level attributes. The Bayesian Network is adopted to learn and compute the probability distributions of the mobility attributes and infer the mobility that most likely arises at a given time. A hierarchical location model with a random forest is developed to refine the destination choice based on the probability obtained from BN. A total of six mobility attributes are simulated, i.e., travel frequency, travel purpose, travel destination location, travel mode, activity duration, and travel duration. The strengths of the proposed framework are: (1) it can model the individual mobility dynamics; (2) it can generate individual mobility for a large-scale population; and (3) it has broader coverage with better representativeness and less strict data requirement. The proposed model is tested with a Miami-Dade County, Florida case study. A total of 7,253,672 travels for 2,680,607 persons and 860,380 households across the county are generated from a public travel diary survey. Results show that individuals at 20–50 age and 0–20 age travel most frequently, with a total of 3,413,710 and 2,074,677 activities per day. The travel peak of a day occurs at 11:00 and 18:00. The results demonstrate the efficacy of the enhanced model, especially in large-scale applications.

The Demographics, Stellar Populations, and Star Formation Histories of Fast Radio Burst Host Galaxies: Implications for the Progenitors

We present a comprehensive catalog of observations and stellar population properties for 23 highly secure host galaxies of fast radio bursts (FRBs). Our sample comprises 6 repeating FRBs and 17 apparent nonrepeaters. We present 82 new photometric and 8 new spectroscopic observations of these hosts. Using stellar population synthesis modeling and employing nonparametric star formation histories (SFHs), we find that FRB hosts have a median stellar mass of ≈109.9 M ⊙, mass-weighted age ≈5.1 Gyr, and ongoing star formation rate ≈1.3 M ⊙ yr−1 but span wide ranges in all properties. Classifying the hosts by degree of star formation, we find that 87% (20 of 23 hosts) are star-forming, two are transitioning, and one is quiescent. The majority trace the star-forming main sequence of galaxies, but at least three FRBs in our sample originate in less-active environments (two nonrepeaters and one repeater). Across all modeled properties, we find no statistically significant distinction between the hosts of repeaters and nonrepeaters. However, the hosts of repeating FRBs generally extend to lower stellar masses, and the hosts of nonrepeaters arise in more optically luminous galaxies. While four of the galaxies with the clearest and most prolonged rises in their SFHs all host repeating FRBs, demonstrating heightened star formation activity in the last ≲100 Myr, one nonrepeating host shows this SFH as well. Our results support progenitor models with short delay channels (i.e., magnetars formed via core-collapse supernova) for most FRBs, but the presence of some FRBs in less-active environments suggests a fraction form through more delayed channels.

The stability of unevenly spaced planetary systems

Studying the orbital stability of multi-planet systems is essential to understand planet formation, estimate the stable time of an observed planetary system, and advance population synthesis models. Although previous studies have primarily focused on ideal systems characterized by uniform orbital separations, in reality a diverse range of orbital separations exists among planets within the same system. This study focuses on investigating the dynamical stability of systems with non-uniform separation. We considered a system with 10 planets with masses of 10−7 solar masses around a central star with a mass of 1 solar mass. We performed more than 100,000 runs of N-body simulations with different parameters. Results demonstrate that reducing merely one pair of planetary spacing leads to an order of magnitude shorter orbital crossing times that could be formulated based on the Keplerian periods of the closest separation pair. Furthermore, the first collisions are found to be closely associated with the first encounter pair that is likely to be the closest separation pair initially. We conclude that when estimating the orbital crossing time and colliding pairs in a realistic situation, updating the formula derived for evenly spaced systems would be necessary.

Lowest accreting protoplanetary discs consistent with X-ray photoevaporation driving their final dispersal

ABSTRACT Photoevaporation from high-energy stellar radiation has been thought to drive the dispersal of protoplanetary discs. Different theoretical models have been proposed, but their predictions diverge in terms of the rate and modality at which discs lose their mass, with significant implications for the formation and evolution of planets. In this paper, we use disc population synthesis models to interpret recent observations of the lowest accreting protoplanetary discs, comparing predictions from EUV-driven, FUV-driven, and X-ray-driven photoevaporation models. We show that the recent observational data of stars with low accretion rates (low accretors) point to X-ray photoevaporation as the preferred mechanism driving the final stages of protoplanetary disc dispersal. We also show that the distribution of accretion rates predicted by the X-ray photoevaporation model is consistent with observations, while other dispersal models tested here are clearly ruled out.

Relation between photometric and parameter errors of star clusters

ABSTRACT Many observations supply only photometry data with large uncertainties. The study of star clusters based on such data is important for astrophysical studies, although photometric uncertainty affects the accuracy of results. In order to estimate the parameter uncertainties that are caused by photometric errors, this work quantifies the influence of photometric error on the uncertainties of some basic parameters when colour-magnitude diagram (CMD) fitting is used to determine these parameters. The observed data are derived from Hubble Space Telescope (HST) observations with the WFPC2 in four bands because of the large coverage of photometric error. The photometric errors that are estimated by Artificial Star Test (AST) are taken, and a binary star stellar population synthesis model and the Powerful CMD code are adopted to determine the parameters of star clusters. Two popular types of stellar population models, i.e. binary star simple stellar population (bsSSP) and single star simple stellar population (ssSSP), are used. The effects of photometric errors on the uncertainties of distance modulus (m − M), colour excess (E), metallicity (Z), and age (t), i.e. Δm − M, ΔE, ΔZ, and Δt, are studied via 19 globular clusters. The results show that Δm − M, ΔE, ΔZ, and Δt have positive correlations with photometric error. Some linear fitting formulae are given to make the results easy to use. The results can be used to estimate the errors of parameters that are caused by photometric error when determining the parameters by CMD fitting, for star clusters with relatively large photometric uncertainties (>∼0.01 mag).