Search NASA⌕ Search

SEARCH · Search NASA

Results for “mass function”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Cosmic Reionization On Computers: The Galaxy-Halo Connection between $5 \leq z \leq10$

We explore the connection between the stellar component of galaxies and their host halos during the epoch of reionization (5 ≤ z ≤ 10) using the CROC (Cosmic Reionization on Computers) simulations. We compare simulated galaxies with observations and find that CROC underpredicts the abundance of luminous galaxies when compared to observed UV luminosity functions, and analogously the most massive galaxies when compared to observed stellar mass functions. We can trace the deficit of star formation to high redshifts, where the slope of the star formation rate to stellar mass relation is consistent with observations, but the normalization is systematically low. This results in a star formation rate density and stellar mass density that are systematically offset from observations. However, the less luminous or lower stellar mass objects have luminosities and stellar masses that agree fairly well with observational data. We explore the stellar-to-halo mass ratio (SHMR), a key quantity that is difficult to measure at high redshifts and that models do not consistently predict. In CROC, the SHMR decreases with redshift, a trend opposite to some abundance-matching studies. These discrepancies uncover where future effort should be focused in order to improve the fidelity of modeling cosmic reionization. We also compare the CROC galaxy bias with observational measurements using Lyman-break galaxy samples, finding reasonable consistency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A mathematical framework for ejecta cloud dynamics with application to source models and piezoelectric mass measurements

We present a mathematical framework for describing the dynamical evolution of an ejecta cloud generated by a generic ejecta source model. We consider a piezoelectric sensor fielded in the path of an ejecta cloud, for experimental configurations in which the ejecta are created at a singly shocked planar surface and fly ballistically through vacuum to the stationary sensor. To do so, we introduce the concept of a time- and velocity-dependent ejecta “areal mass function.” We derive expressions for the analytic (“true”) accumulated ejecta areal mass at the sensor and the measured (“inferred”) value obtained via the standard method for analyzing piezoelectric voltages. In this way, we derive an exact expression and upper bound for the error imposed upon a piezoelectric ejecta mass measurement (in a perfect system) by the assumption of instantaneous creation, which is commonly required for momentum diagnostic analyses. This error term is zero for truly instantaneous source models; otherwise, the standard piezoelectric analysis is guaranteed to overestimate the true mass. When combined with a piezoelectric dataset, this framework provides a unique solution for the ejecta particle velocity distribution, subject to the assumptions inherent in the data analysis. The framework also leads to strong boundary conditions that any ejecta source model must satisfy in order to be consistent with apparently global properties of piezoelectric measurements from a wide range of experiments. We demonstrate this methodology by applying it to the Richtmyer–Meshkov instability+self-similar velocity distribution ejecta source model currently under development at Los Alamos National Laboratory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The evolution of gamma-ray burst jet opening angle through cosmic time

ABSTRACT Jet opening angles of long gamma-ray bursts (lGRBs) appear to evolve in cosmic time, with lGRBs at higher redshifts being on average more narrowly beamed than those at lower redshifts. We examine the nature of this anticorrelation in the context of collimation by the progenitor stellar envelope. First, we show that the data indicate a strong correlation between gamma-ray luminosity and jet opening angle, and suggest this is a natural selection effect – only the most luminous GRBs are able to successfully launch jets with large opening angles. Then, by considering progenitor properties expected to evolve through cosmic time, we show that denser stars lead to more collimated jets; we argue that the apparent anticorrelation between opening angle and redshift can be accounted for if lGRB massive star progenitors at high redshifts have higher average density compared to those at lower redshifts. This may be viable for an evolving initial mass function (IMF) – under the assumption that average density scales directly with mass, this relationship is consistent with the form of the IMF mass evolution suggested in the literature. The jet angle–redshift anticorrelation may also be explained if the lGRB progenitor population is dominated by massive stars at high redshift, while lower redshift lGRBs allow for a greater diversity of progenitor systems (that may fail to collimate the jet as acutely). Overall, however, we find both the jet angle–redshift anticorrelation and jet angle–luminosity correlation are consistent with the conditions of jet launch through, and collimation by, the envelope of a massive star progenitor.

Lloyd-Ronning, Nicole↗
DOE OSTI2021

SDSS-IV MaNGA: Stellar M/L gradients and the M/L-colour relation in galaxies

The stellar mass-to-light ratio gradient in SDSS r-band ∇(M * /L r ) of a galaxy depends on its mass assembly history, which is imprinted in its morphology and gradients of age, metallicity, and stellar initial mass function (IMF). Taking a MaNGA sample of 2051 galaxies with stellar masses ranging from 10 9 to 10 12 M ⊙ released in SDSS DR15, we focus on face-on galaxies, without merger and bar signatures, and investigate the dependence of the 2D ∇(M * /L r ) on other galaxy properties, including M * /L r -colour relationships by assuming a fixed Salpeter IMF as the mass normalization reference. The median gradient is ∇M * /L r ~ –0.1 (i.e. the M * /L r is larger at the centre) for massive galaxies, becomes flat around M * ~ 10 10 M ⊙ and change sign to ∇M * /L r ~ 0.1 at the lowest masses. The M * /L r inside a half-light radius increases with increasing galaxy stellar mass; in each mass bin, early-type galaxies have the highest value, while pure-disc late-type galaxies have the smallest. Correlation analyses suggest that the mass-weighted stellar age is the dominant parameter influencing the M * /L r profile, since a luminosity-weighted age is easily affected by star formation when the specific star formation rate (sSFR) inside the half-light radius is higher than 10 –3 Gyr –1 . With increased sSFR gradient, one can obtain a steeper negative ∇(M * /L r ). The scatter in the slopes of M * /L-colour relations increases with increasing sSFR, for example, the slope for post-starburst galaxies can be flattened to 0.45 from the global value 0.87 in the M * /L versus g – r diagram. Hence converting galaxy colours to M * /L should be done carefully, especially for those galaxies with young luminosity-weighted stellar ages, which can have quite different star formation histories.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Volume-complete Sample of M Dwarfs with Masses 0.1 ≤ M/M{sub ⊙} ≤ 0.3 within 15 Parsecs

M dwarfs with masses 0.1 ≤ M/M {sub ⊙} ≤ 0.3 are under increasing scrutiny because these fully convective stars pose interesting astrophysical questions regarding their magnetic activity and angular momentum history. They also afford the most accessible near future opportunity to study the atmospheres of terrestrial planets. Because they are intrinsically low in luminosity, the identification of the nearest examples of these M dwarfs is essential for progress. We present the volume-complete, all-sky list of 512 M dwarfs with masses 0.1 ≤ M/M {sub ⊙} ≤ 0.3 and with trigonometric distances placing them within 15 pc (π {sub trig} ≥ 66.67 mas) from which we have created a sample of 413 M dwarfs for spectroscopic study. We present the mass function for these 512 M dwarfs, which increases with decreasing stellar mass in linear mass space, but is flat in logarithmic mass space. As part of this sample, we present new V {sub J} R {sub KC} I {sub KC} photometry for 17 targets, measured as a result of the RECONS group’s long-term work at the CTIO/SMARTS 0.9 m telescope. We also note the details of targets that are known to be members of multiple systems and find a preliminary multiplicity rate of 21% ± 2% for the primary M dwarfs in our sample when considering known stellar and brown dwarf companions at all separations from their primaries. We further find that 43 ± 2% of all M dwarfs with masses 0.1 ≤ M/M {sub ⊙} ≤ 0.3 are found in multiple systems with primary stars of all masses within 15 pc.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

(U) A Theoretical Study of Asay Foil Trajectories (V.02)

We consider the trajectory of a generic Asay foil ejecta momentum diagnostic sensor, for scenarios where ejecta are produced at a planar surface and fly ballistically through a perfect vacuum to the sensor. To do so, we build upon a previously established mathematical formalism derived for the analytic study of stationary sensors (i.e., piezopins). First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source areal mass function (“source model”). This defines an integro-differential equation (IDE) for the foil trajectory. When ejecta production is instantaneous - as is generally assumed in momentum diagnostic data analyses - the IDE leads to an implicit and easily calculable closed-form solution for the foil trajectory in a perfect system, as long as the ejecta particle velocity distribution is twice-integrable. General properties of the instant-production trajectory solution indicate the existence of a boundary condition the particle velocity distribution must satisfy in order for the analytically predicted foil trajectories to be compatible with certain features commonly observed in foil data. This boundary condition is identical to one derived previously from a consideration of piezopin data. Armed with the analytic solution for instant production, we also consider various techniques used to extract time-dependent cumulative ejecta masses from foil trajectories, and derive an expression for the error imposed by using an approximated equation of motion. This analytic trajectory solution furthermore makes it possible to examine the common practice of presenting inferred cumulative ejecta masses as a function of a normalized implied velocity; we derive conditions under which this methodology is and is not meaningful. We also propose a strategy for extending the instant-production trajectory solution to time-dependent source functions.

42 ENGINEERING↗
DOE OSTI2021

Photometric and Kinematic Study of the Open Clusters SAI 44 and SAI 45

We carry out a detailed photometric and kinematic study of the poorly studied sparse open clusters SAI 44 and SAI 45 using ground-based BVR {sub c} I {sub c} data supplemented by archival data from Gaia eDR3 and Pan-STARRS. The stellar memberships are determined using a statistical method based on Gaia eDR3 kinematic data, and we found 204 members in SAI 44 while only 74 members are identified in SAI 45. The average distances to SAI 44 and SAI 45 are calculated to be 3670 ± 184 and 1668 ± 47 pc. The logarithmic age of the clusters are determined to be 8.82 ± 0.10 and 9.07 ± 0.10 yr for SAI 44 and SAI 45, respectively. The color–magnitude diagram of SAI 45 hosts an extended main-sequence turnoff (eMSTO). The apparent age spread is found to be similar to the apparent age spread predicted on the basis of the age spread and cluster age relation predicted by rotation models. This indicates that eMSTO is a stellar evolution rather than star formation phenomenon in SAI 45. We conclude that eMSTO in SAI 45 is mainly caused by the different rotation rates of stars as the SYCLIST synthetic population with different rotation rates was able to reproduce the observed eMSTO, and stars in the red part of the eMSTO were preferentially concentrated in the inner region, which again hints at different rotations being the reason for the extension in the upper MS. This finding supports the theory attributing the origin of eMSTO to the different rotations of eMSTO stars. The mass function slopes are obtained as -2.24 ± 0.66 and -2.58 ± 3.20 in the mass rages 2.426–0.990 M{sub ⊙} and 2.167–1.202 M{sub ⊙} for SAI 44 and SAI 45, respectively. SAI 44 exhibits the signature of mass segregation while we found weak evidence of mass segregation in SAI 45 possibly due to tidal stripping. The dynamical relaxation times of these clusters indicate that both clusters are in a dynamically relaxed state. Using the AD-diagram method, the apex coordinates are found to be ([[Formula]], -[[Formula]]) for SAI 44 and (-[[Formula]], -[[Formula]]) for SAI 45. The average space velocity components of the clusters SAI 44 and SAI 45 are calculated in units of km s{sup -1} as (-15.14 ± 3.90, -19.43 ± 4.41, -20.85 ± 4.57) and (28.13 ± 5.30, -9.78 ± 3.13, -19.59 ± 4.43), respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Comprehensive Analysis of NGC 2158 in the Gaia Era: Photometric Parameters, Apex, and Orbit

We present an investigation of NGC 2158 using Gaia DR2 data. We identified 800 most likely cluster members with membership probability higher than 90%. The mean proper motions of this object are determined as (μ {sub x} = −0.203 ± 0.003, μ {sub y} = −1.99 ± 0.004) mas yr{sup −1}. The limiting radius, log(age), and distance of the cluster are obtained as 23.5 arcmin, 9.38 ± 0.04 Gyr, and 4.69 ± 0.22 kpc, respectively. The overall mass function slope (0.93 ± 0.14) is flatter than the Salpeter value (1.35) within the mass range of 1.17–1.44 M {sub ⊙}. This cluster also shows the mass-segregation effect and our study demonstrates that NGC 2158 is a dynamically relaxed open cluster. Using the apex diagram (AD), the apex coordinates of the cluster are obtained in different ways and are examined using the (μ {sub U},μ {sub T}) diagram. The best value of the apex coordinates is determined as A = 8724 ± 160, D = −3661 ± 530. We also determined the orbit of the cluster and found that NGC 2158 moves almost in the solar antapex direction. The resulting spatial velocity of NGC 2158 is 51 km s{sup −1}. A significant oscillation along the Z-coordinate up to 529 pc is detected. Various scenarios regarding the origin of this cluster are also discussed.

42 ENGINEERING↗
DOE OSTI2021

Wide Dust Gaps in Protoplanetary Disks Induced by Eccentric Planets: A Mass-eccentricity Degeneracy

The tidal perturbation of embedded protoplanets on their natal disks has been widely attributed to be the cause of gap-ring structures in submillimeter images of protoplanetary disks around T Tauri stars. Numerical simulations of this process have been used to propose scaling of characteristic dust-gap width/gap-ring distance with respect to planet mass. Applying such scaling to analyze observed gap samples yields a continuous mass distribution for a rich population of hypothetical planets in the range of several Earth to Jupiter masses. In contrast, the conventional core-accretion scenario of planet formation predicts a bimodal mass function due to (1) the onset of runaway gas accretion above ~20 Earth masses and (2) suppression of accretion induced by gap opening. Here, we examine the dust disk response to the tidal perturbation of eccentric planets as a possible resolution of this paradox. Based on simulated gas and dust distributions, we show the gap-ring separation of Neptune-mass planets with small eccentricities might become comparable to that induced by Saturn-mass planets on circular orbits. Furthermore, this degeneracy may obliterate the discrepancy between the theoretical bimodal mass distribution and the observed continuous gap width distribution. Despite damping due to planet–disk interaction, modest eccentricity may be sustained either in the outer regions of relatively thick disks or through resonant excitation among multiple super Earths. Moreover, the ring-like dust distribution induced by planets with small eccentricities is axisymmetric even in low viscosity environments, consistent with the paucity of vortices in Atacama Large Millimeter/submillimeter Array images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Photometry and astrometry with JWST – III. A NIRCam-Gaia DR3 analysis of the open cluster NGC 2506

ABSTRACT In the third paper of this series aimed at developing the tools for analysing resolved stellar populations using the cameras on board of the James Webb Space Telescope (JWST), we present a detailed multiband study of the 2 Gyr Galactic open cluster NGC 2506. We employ public calibration data sets collected in multiple filters to: (i) derive improved effective Point Spread Functions (ePSFs) for 10 NIRCam filters; (ii) extract high-precision photometry and astrometry for stars in the cluster, approaching the main sequence (MS) lower mass of ∼0.1 M⊙; and (iii) take advantage of the synergy between JWST and Gaia DR3 to perform a comprehensive analysis of the cluster’s global and local properties. We derived a MS binary fraction of ∼57.5 per cent, extending the Gaia limit (∼0.8 M⊙) to lower masses (∼0.4 M⊙) with JWST. We conducted a study on the mass functions (MFs) of NGC 2506, mapping the mass segregation with Gaia data, and extending MFs to lower masses with the JWST field. We also combined information on the derived MFs to infer an estimate of the cluster present-day total mass. Lastly, we investigated the presence of white dwarfs (WDs) and identified a strong candidate. However, to firmly establish its cluster membership, as well as that of four other WD candidates and of the majority of faint low-mass MS stars, further JWST equally deep observations will be required. We make publicly available catalogues, atlases, and the improved ePSFs.

Nardiello, D. (ORCID:0000000311493659)↗
DOE OSTI2020

An Investigation of Poorly Studied Open Cluster NGC 4337 Using Multicolor Photometric and Gaia DR2 Astrometric Data

We present a comprehensive analysis (photometric and kinematical) of the poorly studied open cluster NGC 4337 using 2MASS, WISE, APASS, and Gaia DR2 databases. By determining the membership probabilities of stars, we identified the 624 most probable members with membership probability higher than 50% by using proper motion and parallax data taken from Gaia DR2. The mean proper motion of the cluster is obtained as μ{sub x}=−8.83±0.01 and μ {sub y} = 1.49 ± 0.006 mas yr{sup −1}. We find the normal interstellar extinction toward the cluster region. The radial distribution of members provides a cluster radius of 7.′75 (5.63 pc). The estimated age of 1600 ± 180 Myr indicates that NGC 4337 is an old open cluster with a bunch of red giant stars. The overall mass function slope for main-sequence stars is found as 1.46±0.18 within the mass range 0.75−2.0 M{sub ⊙}, which is in fair agreement with Salpeter’s value (x = 1.35) within uncertainty. The present study demonstrates that NGC 4337 is a dynamically relaxed open cluster. Using the Galactic potential model, Galactic orbits are obtained for NGC 4337. We found that this object follows a circular path around the Galactic center. Under the kinematical analysis, we compute the apex coordinates (A, D) by using two methods: (i) the classical convergent point method and (ii) the AD-diagram method. The obtained coordinates are (A {sub conv}, D {sub conv}) = (96.°27 ± 0.°10, 13.°14 ± 0.°27) and (A {sub ◦}, D {sub ◦}) = (100.°282 ± 0.°10, 9.°577 ± 0.°323) respectively. We also computed the Velocity Ellipsoid Parameters, matrix elements (μ {sub ij}), direction cosines (l {sub j}, m {sub j}, n {sub j}), and the Galactic longitude of the vertex (l {sub 2}).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

COOL-LAMPS. II. Characterizing the Size and Star Formation History of a Bright Strongly Lensed Early-type Galaxy at Redshift 1.02

We present COOL J1323+0343, an early-type galaxy at z = 1.0153 ± 0.0006, strongly lensed by a cluster of galaxies at z = 0.353 ± 0.001. This object was originally imaged by DECaLS and noted as a gravitational lens by COOL-LAMPS, a collaboration initiated to find strong-lensing systems in recent public optical imaging data, and confirmed with follow-up data. With ground-based grzH imaging and optical spectroscopy from the Las Campanas Observatory and the Nordic Optical Telescope, we derive a stellar mass, metallicity, and star formation history from stellar-population synthesis modeling. The lens modeling implies a total magnification, summed over the three images in the arc, of μ ~ 113. The stellar mass in the source plane is M * ~ 10.64 M ⊙ and the 1σ upper limit on the star formation rate (SFR) in the source plane is SFR ~ 3.75 × 10 -2 M ⊙ yr -1 (log sSFR = -12.1 yr -1 ) in the youngest two age bins (0–100 Myr), closest to the epoch of observation. Our measurements place COOL J1323+0343 below the characteristic mass of the stellar mass function, making it an especially compelling target that could help clarify how intermediate-mass quiescent galaxies evolve. We reconstruct COOL J1323+0343 in the source plane and fit its light profile. This object is below the expected size evolution of an early-type galaxy at this mass with an effective radius r e ~ 0.5 kpc. This extraordinarily magnified and bright lensed early-type galaxy offers an exciting opportunity to study the morphology and star formation history of an intermediate-mass early-type galaxy in detail at z ~ 1.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

SIDM Concerto: Compilation and Data Release of Self-interacting Dark Matter Zoom-in Simulations

We present SIDM Concerto: 14 cosmological zoom-in simulations in cold dark matter (CDM) and self-interacting dark matter (SIDM) models based on the Symphony and Milky Way-est suites. SIDM Concerto includes one Large Magellanic Cloud– (LMC-) mass system (host mass ∼10 11 M ⊙ ), two Milky Way (MW) analogs (∼10 12 M ⊙ ), two group-mass hosts (∼10 13 M ⊙ ), and one low-mass cluster (∼10 14 M ⊙ ). Each host contains ≈2 × 10 7 particles and is run in CDM and one or more strong, velocity-dependent SIDM models. Our analysis of SIDM (sub)halo populations over seven subhalo mass decades reveals that (1) the fraction of core-collapsed isolated halos and subhalos peaks at a maximum circular velocity corresponding to the transition of the SIDM cross section from a v −4 to v 0 scaling; (2) SIDM subhalo mass functions are suppressed by ≈50% relative to CDM in LMC, MW, and group-mass hosts but are consistent with CDM in the low-mass cluster host; (3) subhalos’ inner density profile slopes, which are more diverse in SIDM than in CDM, are sensitive to both the amplitude and shape of the SIDM cross section. Our simulations provide a benchmark for testing SIDM predictions with astrophysical observations of field and satellite galaxies, strong lensing systems, and stellar streams. Data products are publicly available at doi:10.5281/zenodo.14933624.

dark matter↗
DOE OSTI2022

Fuzzy dark matter and the Dark Energy Survey Year 1 data

ABSTRACT Gravitational weak lensing by dark matter haloes leads to a measurable imprint in the shear correlation function of galaxies. Fuzzy dark matter (FDM), composed of ultralight axion-like particles of mass m ∼ 10−22 eV, suppresses the matter power spectrum and shear correlation with respect to standard cold dark matter. We model the effect of FDM on cosmic shear using the optimized halo model HMCode, accounting for additional suppression of the mass function and halo concentration in FDM as observed in N-body simulations. We combine Dark Energy Survey Year 1 (DES-Y1) data with the Planck cosmic microwave background anisotropies to search for shear correlation suppression caused by FDM. We find no evidence of suppression compared to the preferred cold dark matter model, and thus set a new lower limit to the FDM particle mass. Using a log-flat prior and marginalizing over uncertainties related to the non-linear model of FDM, we find a new, independent 95 per cent C.L. lower limit log10m > −23 combining Planck and DES-Y1 shear, an improvement of almost two orders of magnitude on the mass bound relative to CMB-only constraints. Our analysis is largely independent of baryonic modelling, and of previous limits to FDM covering this mass range. Our analysis highlights the most important aspects of the FDM non-linear model for future investigation. The limit to FDM from weak lensing could be improved by up to three orders of magnitude with $\mathcal {O}(0.1)$ arcmin cosmic shear angular resolution, if FDM and baryonic feedback can be simultaneously modelled to high precision in the halo model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Deep Study of an Intermediate-age Open Cluster SAI 35 (Juchert 20) Using Ground-based Imaging and Gaia EDR3 Astrometry

We present a CCD UBVI photometric study of poorly studied intermediate-age open cluster SAI 35 (Juchert 20) for the first time. To accomplish this study, we also used LAMOST DR5, Two Micron All Sky Survey, and Gaia EDR3 databases. We identified 214 most probable cluster members with membership probability higher than 50%. The mean proper motion of the cluster is found as μ{sub α}cos δ=1.10±0.01 and μ {sub δ} = −1.66 ± 0.01 mas yr{sup −1}. We find the normal interstellar extinction law using the various two-color diagrams. The age, distance, reddening, and radial velocity of the cluster are estimated to be 360 ± 40 Myr, 2.9 ± 0.15 kpc, 0.72 ± 0.05 mag, and −91.62 ± 6.39 km s{sup −1}, respectively. The overall mass function slope for main-sequence stars is found to be 1.49 ± 0.16 within the mass range 1.1–3.1 M {sub ⊙}, which is in agreement with Salpeter’s value within uncertainty. The present study demonstrates that SAI 35 is dynamically relaxed. Galactic orbital parameters are determined using Galactic potential models. We found that this object follows a circular path around the Galactic center.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Inferring Type II-P Supernova Progenitor Masses from Plateau Luminosities

Abstract Connecting observations of core-collapse supernova explosions to the properties of their massive star progenitors is a long-sought, and challenging, goal of supernova science. Recently, Barker et al. presented bolometric light curves for a landscape of progenitors from spherically symmetric neutrino-driven core-collapse supernova (CCSN) simulations using an effective model. They find a tight relationship between the plateau luminosity of the Type II-P CCSN light curve and the terminal iron-core mass of the progenitor. Remarkably, this allows us to constrain progenitor properties with photometry alone. We analyze a large observational sample of Type II-P CCSN light curves and estimate a distribution of iron-core masses using the relationship of Barker et al. The inferred distribution matches extremely well with the distribution of iron-core masses from stellar evolutionary models and namely, contains high-mass iron cores that suggest contributions from very massive progenitors in the observational data. We use this distribution of iron-core masses to infer minimum and maximum masses of progenitors in the observational data. Using Bayesian inference methods to locate optimal initial mass function parameters, we find M min = 9.8 − 0.27 + 0.37 and M max = 24.0 − 1.9 + 3.9 solar masses for the observational data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Filling the black hole mass gap: Avoiding pair instability in massive stars through addition of nonnuclear energy

In standard stellar evolution, stars with masses ranging from approximately 150 to 240M ⊙ are expected to evolve to a pair instability supernova with no black hole (BH) remnant. This evolutionary behavior leads to a predicted gap in the black hole mass function from approximately 50 to 140M ⊙ . Yet the LIGO and Virgo Collaborations recently discovered black holes of masses 66M ⊙ and 85M ⊙ in the gravitational wave event GW190521. Here, we propose a new method to populate the BH mass gap. If an energy source is added throughout the star in addition to nuclear fusion, it is possible for the altered evolution to avoid the complete destruction of a pair instability supernova, and instead a BH remnant is left behind. An example of an extra energy source is dark matter annihilation within the star, but our results hold more generally. We show this phenomenon by exploring the effect of adding an energy source independent of temperature and density to a 180M ⊙ star, using the MESA one-dimensional stellar evolution software. If ~50% of the star’s energy is due to this new source, the star is capable of avoiding the pair instability entirely and evolving towards a core-collapse supernova and ultimately a BH remnant with mass ~ 120M ⊙ .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Discovery of a Gamma-Ray Black Widow Pulsar by GPU-accelerated Einstein@Home

We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653-0158. The associated Fermi Large Area Telescope gamma-ray source 4FGL J1653.6-0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations—whose periodic brightness modulations suggested an orbit—no radio pulsations had been found in many searches. The pulsar was discovered by directly searching the gamma-ray data using the GPU-accelerated Einstein@Home distributed volunteer computing system. The multidimensional parameter space was bounded by positional and orbital constraints obtained from the optical counterpart. More sensitive analyses of archival and new radio data using knowledge of the pulsar timing solution yield very stringent upper limits on radio emission. Any radio emission is thus either exceptionally weak, or eclipsed for a large fraction of the time. The pulsar has one of the three lowest inferred surface magnetic-field strengths of any known pulsar with B surf ≈ 4 × 10 7 G. The resulting mass function, combined with models of the companion star's optical light curve and spectra, suggests a pulsar mass gsim 2 M ⊙ . The companion is lightweight with mass ~0.01 M ⊙ , and the orbital period is the shortest known for any rotation-powered binary pulsar. This discovery demonstrates the Fermi Large Area Telescope's potential to discover extreme pulsars that would otherwise remain undetected.

79 ASTRONOMY AND ASTROPHYSICS↗