Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STELLAR MASS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

An Educational Guide for 2D Stellar Structure Calculations of Rapidly Rotating Stars using the ESTER code

The Evolution STEllaire en Rotation (ESTER) code is the first 2D stellar structure code to be made open-source and freely available to the astronomy and astrophysics community. An important and novel advancement of this code is that it can reproduce the distorted shape and observable signatures (e.g., gravity darkening) of rapidly rotating stars. ESTER also calculates the steady-state large-scale flows within the star, namely their differential rotation and associated meridional circulation. In this report, we explore and document the physics implemented within version 1.1.0rc2 of the ESTER code, in a way that complements published descriptions. We illustrate this physics by plotting how stellar structure parameters vary through stellar interiors at a range of latitudes and at different angular velocities. We investigate how the thin convective envelopes of intermediate mass stars vary with latitude when rapidly rotating, becoming deeper and thicker near the equator. Simple comparisons of ESTER model predictions (e.g., central temperature and density, luminosity) with the output from the Modules for Experiments in Stellar Astrophysics (MESA) code [Paxton et al., 2010] shows generally good agreement. Additional comparisons provide important benchmarking and verification for ESTER as a comparatively young code. Finally, we provide a guide for installing and running the code on our local university cluster, aimed at helping students to begin work.

79 ASTRONOMY AND ASTROPHYSICS↗

Stellar multiplicity and stellar rotation: insights from APOGEE

ABSTRACT We measure rotational broadening in spectra taken by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey to characterize the relationship between stellar multiplicity and rotation. We create a sample of 2786 giants and 24 496 dwarfs with stellar parameters and multiple radial velocities from the APOGEE pipeline, projected rotation speeds vsin i determined from our own pipeline, and distances, masses, and ages measured by Sanders & Das. We use the statistical distribution of the maximum shift in the radial velocities, ΔRVmax, as a proxy for the close binary fraction to explore the interplay between stellar evolution, rotation, and multiplicity. Assuming that the minimum orbital period allowed is the critical period for Roche Lobe overflow and rotational synchronization, we calculate theoretical upper limits on expected vsin i and ΔRVmax values. These expectations agree with the positive correlation between the maximum ΔRVmax and vsin i values observed in our sample as a function of log(g). We find that the fast rotators in our sample have a high occurrence of short-period [log(P/d) ≲ 4] companions. We also find that old, rapidly rotating main-sequence stars have larger completeness-corrected close binary fractions than their younger peers. Furthermore, rapidly rotating stars with large ΔRVmax consistently show differences of 1–10 Gyr between the predicted gyrochronological and measured isochronal ages. These results point towards a link between rapid rotation and close binarity through tidal interactions. We conclude that stellar rotation is strongly correlated with stellar multiplicity in the field, and caution should be taken in the application of gyrochronology relations to cool stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Galaxy formation with wave/fuzzy dark matter: The core-halo structure and the solitonic imprint

Dark matter-dominated cores have long been claimed for the well-studied local group dwarf galaxies. More recently, extended stellar halos have been uncovered around several of these dwarfs through deeper imaging and spectroscopy. Such core-halo structures (inner flat core and a characteristic r −3 asymptotic outer halo profile) are not a feature of conventional cold dark matter (CDM). In contrast, smooth and prominent dark matter cores are predicted for wave/fuzzy dark matter (ψDM). The question arises as to what extent the visible stellar profiles should reflect this dark matter core structure. Here we compare cosmological hydrodynamical simulations of CDM, “WDM” (model used as a proxy for ψDM) & ψDM, aiming to predict the stellar profiles for these three DM scenarios. We show that cores surrounded by extended halos are distinguishable for ψDM, where the stellar density is enhanced in the core due to the presence of the relatively dense soliton. Our analysis demonstrates that, in our simulations, a distinctive core-halo structure does not appear in the case of CDM in the DM, gas, or stars. Whereas we do find a core-halo transition for DM, gas, and stars for ψDM, and the scale of this transition is in line with the predicted core radius set by the soliton scale anticipated for the adopted boson mass of 2.5×10 −22 eV. The presence of a core-halo structure in the stellar profile for Galaxy 1 for ψDM is visible for the most massive and the first galaxy to form in the simulation. Clearly, further simulations are needed to establish how strict this possible relationship is between the DM and stellar core-halo profile as a potential observational discriminator. Furthermore, we observe the anticipated asymmetry for ψDM due to the soliton's motion (jumping and random walk), a distinctive characteristic not found in the symmetric distributions of stars in the warm and CDM models.

dark matter↗

The role of the turbulence driving mode for the initial mass function

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating, and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ∼1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF.

79 ASTRONOMY AND ASTROPHYSICS↗

GRB 171205A: Hypernova and Newborn Neutron Star

GRB 171205A is a low-luminosity, long-duration gamma-ray burst (GRB) associated with SN 2017iuk, a broad-line type Ic supernova (SN). It is consistent with having been formed in the core collapse of a widely separated binary, which we have called the binary-driven hypernova of type III. The core collapse of the CO star forms a newborn NS (νNS) and the SN explosion. Fallback accretion transfers mass and angular momentum to the νNS, here assumed to be born non-rotating. The accretion energy injected into the expanding stellar layers powers the prompt emission. The multiwavelength power-law afterglow is explained by the synchrotron radiation of electrons in the SN ejecta, powered by energy injected by the spinning νNS. We calculate the amount of mass and angular momentum gained by the νNS, as well as the νNS rotational evolution. The νNS spins up to a period of 47 ms, then releases its rotational energy powering the synchrotron emission of the afterglow. The paucity of the νNS spin explains the low-luminosity characteristic and that the optical emission of the SN from the nickel radioactive decay outshines the optical emission from the synchrotron radiation. From the νNS evolution, we infer that the SN explosion had to occur at most 7.36 h before the GRB trigger. Therefore, for the first time, the analysis of the GRB data leads to the time of occurrence of the CO core collapse leading to the SN explosion and the electromagnetic emission of the GRB event.

79 ASTRONOMY AND ASTROPHYSICS↗

Mass measurements show slowdown of rapid proton capture process at waiting-point nucleus 64Ge

Abstract X-ray bursts are among the brightest stellar objects frequently observed in the sky by space-based telescopes. A type-I X-ray burst is understood as a violent thermonuclear explosion on the surface of a neutron star, accreting matter from a companion star in a binary system. The bursts are powered by a nuclear reaction sequence known as the rapid proton capture process (rp process), which involves hundreds of exotic neutron-deficient nuclides. At so-called waiting-point nuclides, the process stalls until a slower β + decay enables a bypass. One of the handful of rp process waiting-point nuclides is 64 Ge, which plays a decisive role in matter flow and therefore the produced X-ray flux. Here we report precision measurements of the masses of 63 Ge, 64,65 As and 66,67 Se—the relevant nuclear masses around the waiting-point 64 Ge—and use them as inputs for X-ray burst model calculations. We obtain the X-ray burst light curve to constrain the neutron-star compactness, and suggest that the distance to the X-ray burster GS 1826–24 needs to be increased by about 6.5% to match astronomical observations. The nucleosynthesis results affect the thermal structure of accreting neutron stars, which will subsequently modify the calculations of associated observables.

Physics↗

Carbon- and Oxygen-rich stars in MaStar: identification and classification

ABSTRACT Carbon- and Oxygen-rich stars populating the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase of stellar evolution are relevant contributors to the spectra of ∼1 Gyr old populations. Atmosphere models for these types are uncertain, due to complex molecules and mass-loss effects. Empirical spectra are then crucial, but samples are small due to the short (∼3 Myr) TP-AGB lifetime. Here we exploit the vastness of the MaNGA Stellar library MaStar (∼60 000 spectra) to identify C and O-rich type stars. We define an optical colour selection with cuts of (g − r) > 2 and (g − i) < 1.55(g − r) − 0.07, calibrated with known C- and O-rich spectra. This identifies C- and O-rich stars along clean, separated sequences. An analogue selection is found in V, R, I bands. Our equation identifies C-rich and O-rich spectra with predictive performance metric F1-scores of 0.72 and 0.74 (over 1), respectively. We finally identify 41 C- and 87 O-rich type AGB stars in MaStar, 5 and 49 of which do not have a SIMBAD counterpart. We also detect a sample of non-AGB, dwarf C-stars. We further design a fitting procedure to classify the spectra into broad spectral types, by using as fitting templates empirical C- and O-rich spectra. We find remarkably good fits for the majority of candidates and categorize them into C- and O-rich bins following existing classifications, which correlate to effective temperature. Our selection models can be applied to large photometric surveys (e.g. Euclid, Rubin). The classified spectra will facilitate future evolutionary population synthesis models.

Hill, Lewis (ORCID:0000000180528080)↗

Stellar population analysis of MaNGA early-type galaxies: IMF dependence and systematic effects

We study systematics associated with estimating simple stellar population (SSP) parameters – age, metallicity [M/H], α-enhancement [α/Fe], and initial mass function (IMF) shape – and associated M * /L gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs), and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) having to subtract emission before estimating absorption line strengths; (iii) the decision to fit the whole spectrum or just a few absorption lines; (iv) SSP model differences (e.g. isochrones, enrichment, IMF). The MILES+Padova SSP models, fit to the H β , 〈Fe〉, TiO 2SDSS , and [MgFe] Lick indices in the stacks, indicate that out to the half-light radius R e : (a) ages are younger and [α/Fe] values are lower in the central regions but the opposite is true of [M/H]; (b) the IMF is more bottom-heavy in the center, but is close to Kroupa beyond about R e /2; (c) this makes M * /L about 2 × larger in the central regions than beyond R e /2. While the models of Conroy et al. return similar [M/H] and [α/Fe] profiles, the age and (hence) M * /L profiles can differ significantly even for solar abundances and a Kroupa IMF; different responses to non-solar abundances and IMF parametrization further compound these differences. There are clear (model independent) differences between E-SRs, E-FRs, and S0s: younger ages and less enhanced [α/Fe] values suggest that E-FRs and S0s are not SSPs, but relaxing this assumption is unlikely to change their inferred M * /L gradients significantly.

79 ASTRONOMY AND ASTROPHYSICS↗

New Fe 59 Stellar Decay Rate with Implications for the Fe 60 Radioactivity in Massive Stars

The discrepancy between observations from γ-ray astronomy of the 60 Fe/ 26 Al γ-ray flux ratio and recent calculations is an unresolved puzzle in nuclear astrophysics. The stellar β-decay rate of 59 Fe is one of the major nuclear uncertainties impeding us from a precise prediction. The important Gamow-Teller strengths from the low-lying states in 59 Fe to the 59 Co ground state are measured for the first time using the exclusive measurement of the 59 Co(t,3He+γ) 59 Fe charge-exchange reaction. The new stellar decay rate of 59 Fe is a factor of 3.5±1.1 larger than the currently adopted rate at T=1.2 GK. Stellar evolution calculations show that the 60 Fe production yield of an 18 solar mass star is decreased significantly by 40% when using the new rate. Our result eliminates one of the major nuclear uncertainties in the predicted yield of 60 Fe and alleviates the existing discrepancy of the 60 Fe/ 26 Al ratio.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Constraints on Metastable Helium in the Atmospheres of WASP-69b and WASP-52b with Ultranarrowband Photometry

Infrared observations of metastable 2{sup 3}S helium absorption with ground- and space-based spectroscopy are rapidly maturing, as this species is a unique probe of exoplanet atmospheres. Specifically, the transit depth in the triplet feature (with vacuum wavelengths near 1083.3 nm) can be used to constrain the temperature and mass-loss rate of an exoplanet’s upper atmosphere. Here, we present a new photometric technique to measure metastable 2{sup 3}S helium absorption using an ultranarrowband filter (FWHM 0.635 nm) coupled to a beam-shaping diffuser installed in the Wide-field Infrared Camera on the 200 inch Hale Telescope at Palomar Observatory. We use telluric OH lines and a helium arc lamp to characterize refractive effects through the filter and to confirm our understanding of the filter transmission profile. We benchmark our new technique by observing a transit of WASP-69b and detect an excess absorption of 0.498% ± 0.045% (11.1σ), consistent with previous measurements after considering our bandpass. We then use this method to study the inflated gas giant WASP-52b and place a 95th percentile upper limit on excess absorption in our helium bandpass of 0.47%. Using an atmospheric escape model, we constrain the mass-loss rate for WASP-69b to be 5.25{sub −0.46}{sup +0.65}×10{sup −4} M{sub J} Gyr{sup −1} (3.32{sub −0.56}{sup +0.67}×10{sup −3} M{sub J} Gyr{sup −1}) at 7000 K (12,000 K). Additionally, we set an upper limit on the mass-loss rate of WASP-52b at these temperatures of 2.1×10{sup −4} M{sub J} Gyr{sup −1} (2.1×10{sup −3} M{sub J} Gyr{sup −1}). These results show that ultranarrowband photometry can reliably quantify absorption in the metastable helium feature.

79 ASTRONOMY AND ASTROPHYSICS↗

The Ubiquity of AGN Winds in Seyfert 1 Galaxies at Low Redshift

The optical spectra of 3896 Seyfert 1 (Sy1) galaxies detected with WISE at z < 0.4 were analyzed for evidence of outflows. In 37% of the Sy1s in our sample, the outflows appear as broad, blue-shifted, spectrally resolved components of the [O iii]λ5007 Å emission line, with a mean maximum velocity V {sub max} ∼ 1014 km s{sup −1} that is consistent with AGN winds. For each Sy1, we deduced the black hole (BH) mass, bolometric luminosity, Eddington ratio, and power-law index of the continuum, which we compared with the star formation rate (SFR) and host morphology. Having separated our sample in two spectroscopic subgroups—Sy1s with only broad Balmer lines (Sy1B) and those with both narrow and broad (Sy1N) lines—and distinguishing those that show as outflow (Sy1Bw and Sy1Nw), we report the following differences: (1) the BH mass is systematically higher and the power-law steeper in the Sy1B–Sy1Bw than in the Sy1N–Sy1Nw; (2) V {sub max} is higher in the Sy1Bw than in the Sy1Nw, correlated in both groups with the BH mass and bolometric luminosity; (3) the Eddington ratio and SFR are higher in the Sy1 with outflows; and (4) the specific star formation rates (sSFRs) of the Sy1s are normal for their morphology and mass, typical of early-type spiral galaxies in the green valley, far from the quenched regime. From these results, we conclude that AGN winds in Sy1s are triggered by higher accretion rates and probably radiatively launched, and there is no clear evidence of an effect on the star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

The Mass-loss History of the Red Hypergiant VY CMa

Imaging and spectroscopy of the knots, clumps, and extended arcs in the complex ejecta of VY CMa confirm a record of high mass-loss events over the past few hundred years. Hubble Space Telescope/Space Telescope Imaging Spectrograph spectroscopy of numerous small knots close to the star allow us to measure their radial velocities from the strong K i emission and determine their separate motions, spatial orientations, and time since ejecta. Their ages concentrate around 70, 120, 200, and 250 yr ago. A K i emission knot only 50 mas from the star ejected as recently as 1985–1995 may coincide with an H{sub 2}O maser. Comparison with VY CMa’s historic light curve from 1800 to the present shows several knots with ejection times that correspond with extended periods of variability and deep minima. The similarity of this correspondence in VY CMa with the remarkable recent dimming of Betelgeuse and an outflow of gas is apparent. The evidence for similar outflows from the surface of a more typical red supergiant suggests that discrete ejections are more common and surface or convective activity is a major source of mass loss for red supergiants.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

LBT Reveals Large Dust Particles and a High Mass-loss Rate for K2-22 b

The disintegrating planet candidate K2-22 b shows periodic and stochastic transits best explained by an escaping debris cloud. However, the mechanism that creates the debris cloud is unknown. The grain size of the debris as well as its sublimation rate can be helpful in understanding the environment that disintegrates the planet. Here, we present simultaneous photometry with the g band at 0.48 μm and K {sub S} band at 2.1 μm using the Large Binocular Telescope. During an event with very low dust activity, we put a new upper limit on the size of the planet of 0.71 R {sub ⊕} or 4500 km. We also detected a medium depth transit that can be used to constrain the dust particle sizes. We find that the median particle size must be larger than about 0.5–1.0 μm, depending on the composition of the debris. This leads to a high mass-loss rate of about 3 × 10{sup 8} kg s{sup −1}, which is consistent with hydrodynamic escape models. If they are produced by some alternate mechanism such as explosive volcanism, it would require extraordinary geological activity. Combining our upper limits on the planet size with the high mass-loss rate, we find a lifetime of the planet of less than 370 Myr. This drops to just 21 Myr when adopting the 0.02 M {sub ⊕} mass predicted from hydrodynamical models.

47 OTHER INSTRUMENTATION↗

Nondetection of Helium in the Upper Atmospheres of TRAPPIST-1b, e, and f

We obtained high-resolution spectra of the ultracool M-dwarf TRAPPIST-1 during the transit of its planet “b” using two high-dispersion near-infrared spectrographs, the Infrared Doppler (IRD) instrument on the Subaru 8.2m telescope, and the Habitable Zone Planet Finder (HPF) instrument on the 10 m Hobby–Eberly Telescope. These spectroscopic observations are complemented by a photometric transit observation for planet “b” using the APO/ARCTIC, which assisted us in capturing the correct transit times for our transit spectroscopy. Using the data obtained by the new IRD and HPF observations, as well as the prior transit observations of planets “b,” “e” and “f” from IRD, we attempt to constrain the atmospheric escape of the planet using the He i triplet 10830 Å absorption line. We do not detect evidence for any primordial extended H-He atmospheres in all three planets. To limit any planet-related absorption, we place an upper limit on the equivalent widths of <7.754 mÅ for planet “b,” <10.458 mÅ for planet “e,” <4.143 mÅ for planet “f” at 95% confidence from the IRD data, and <3.467 mÅ for planet “b” at 95% confidence from HPF data. Using these limits along with a solar-like composition isothermal Parker wind model, we attempt to constrain the mass-loss rates for the three planets. For TRAPPIST-1b, our models exclude the highest possible energy-limited rate for a wind temperature <5000 K. This nondetection of extended atmospheres with low mean-molecular weights in all three planets aids in further constraining their atmospheric composition by steering the focus toward the search of high-molecular-weight species in their atmospheres.

74 ATOMIC AND MOLECULAR PHYSICS↗

The Initial–Final Mass Relation for Hydrogen-deficient White Dwarfs

The initial–final mass relation represents the total mass lost by a star during the entirety of its evolution from the zero age main sequence to the white-dwarf cooling track. The semiempirical initial–final mass relation (IFMR) is largely based on observations of DA white dwarfs, the most common spectral type of white dwarf and the simplest atmosphere to model. We present a first derivation of the semiempirical IFMR for hydrogen-deficient (non-DA) white dwarfs in open star clusters. We identify a possible discrepancy between the DA and non-DA IFMRs, with non-DA white dwarfs ≈0.07 M {sub ⊙} less massive at a given initial mass. Such a discrepancy is unexpected based on theoretical models of non-DA formation and observations of field white dwarf mass distributions. If real, the discrepancy is likely due to enhanced mass loss during the final thermal pulse and renewed post-AGB evolution of the star. However, we are dubious that the mass discrepancy is physical and instead is due to the small sample size, to systematic issues in model atmospheres of non-DAs, and to the uncertain evolutionary history of Procyon B (spectral type DQZ). A significantly larger sample size is needed to test these assertions. In addition, we also present Monte Carlo models of the correlated errors for DA and non-DA white dwarfs in the initial–final mass plane. We find the uncertainties in initial–final mass determinations for individual white dwarfs can be significantly asymmetric, but the recovered functional form of the IFMR is grossly unaffected by the correlated errors.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Anticorrelation between Duration and Redshift in Gamma-Ray Bursts

Abstract For gamma-ray bursts (GRBs) with durations greater than two seconds (so-called long GRBs), the intrinsic prompt gamma-ray emission appears, on average, to last longer for bursts at lower redshifts. We explore the nature of this duration–redshift anticorrelation, describing systems and conditions in which this cosmological evolution could arise. In particular, we explore its dependence on the metallicity of a massive star progenitor, because we can securely count on the average stellar metallicity to increase with decreasing redshift. Although stars with higher metallicity/lower redshift lose mass and angular momentum through line-driven winds, in some cases these stars are able to form more extended accretion disks when they collapse, potentially leading to longer-duration GRBs. We also examine how this duration–redshift trend may show up in interacting binary models composed of a massive star and compact object companion, recently suggested to be the progenitors of radio-bright GRBs. Under certain conditions, mass loss and equation-of-state effects from massive stars with higher metallicity and lower redshift can decrease the binary separation. This can then lead to spin-up of the massive star and allow for a longer-duration GRB upon the massive star’s collapse. Finally, the duration–redshift trend may also be supported by a relatively larger population of small-separation binaries born in situ at low redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way

We present the first Event Horizon Telescope (EHT) observations of Sagittarius A* (Sgr A*), the Galactic center source associated with a supermassive black hole. These observations were conducted in 2017 using a global interferometric array of eight telescopes operating at a wavelength of λ = 1.3 mm. The EHT data resolve a compact emission region with intrahour variability. A variety of imaging and modeling analyses all support an image that is dominated by a bright, thick ring with a diameter of 51.8 ± 2.3 μas (68% credible interval). The ring has modest azimuthal brightness asymmetry and a comparatively dim interior. Using a large suite of numerical simulations, we demonstrate that the EHT images of Sgr A* are consistent with the expected appearance of a Kerr black hole with mass ∼4 × 10$^{6}$ M $_{⊙}$, which is inferred to exist at this location based on previous infrared observations of individual stellar orbits, as well as maser proper-motion studies. Our model comparisons disfavor scenarios where the black hole is viewed at high inclination (i > 50°), as well as nonspinning black holes and those with retrograde accretion disks. Our results provide direct evidence for the presence of a supermassive black hole at the center of the Milky Way, and for the first time we connect the predictions from dynamical measurements of stellar orbits on scales of 10$^{3}$–10$^{5}$ gravitational radii to event-horizon-scale images and variability. Furthermore, a comparison with the EHT results for the supermassive black hole M87* shows consistency with the predictions of general relativity spanning over three orders of magnitude in central mass.

79 ASTRONOMY AND ASTROPHYSICS↗

Predicting the Galactic population of free-floating planets from realistic initial conditions

ABSTRACT We present the first prediction for the mass distribution function of Galactic free-floating planets (FFPs) that aims to accurately include the relative contributions of multiple formation pathways and stellar populations. We derive our predicted distribution from dedicated simulations of planet birth, growth, migration, and ejection around circumbinary systems and extend these results to also include the contributions from single and wide binary systems. Our resulting FFP mass distribution shows several distinct features, including a strong peak at ${\sim} 8\, {\rm M}_{\oplus }$ arising from the transition between pebble and gas accretion regimes and a trough at ${\sim} 1\, {\rm M}_{\oplus }$ due to the shift in the dominant ejection process from planet–planet scattering to ejection through interactions with stars in circumbinary systems. We find that interactions with the central binary in close circumbinary systems are likely the dominant progenitor for FFPs more massive than Earth, leading to a steep power-law dependence in mass that agrees well with existing observations. In contrast, we find planet–planet scattering events in single and wide binary systems likely produce the majority of planets at Mars mass and below, resulting in a shallower power-law dependence. Our results suggest that existing extrapolations into the sub-terrestrial mass range may significantly overestimate the true FFP abundance. The features we predict in the mass distribution of FFPs will be detectable by upcoming space-based microlensing surveys and, if observed, will provide key insight into the origins of FFPs and the environments in which they form.

Coleman, Gavin A. L. (ORCID:0000000151118963)↗