Using ground-based gravitational-wave detectors, we probe the mass function of intermediate-mass black holes (IMBHs) wherein we also include BHs in the upper mass gap at ∼60–130 M ⊙. Employing the projected sensitivity of the upcoming LIGO and Virgo fourth observing run (O4), we perform Bayesian analysis on quasi-circular nonprecessing, spinning IMBH binaries (IMBHBs) with total masses 50–500 M ⊙, mass ratios 1.25, 4, and 10, and dimensionless spins up to 0.95, and estimate the precision with which the source-frame parameters can be measured. We find that, at 2σ, the mass of the heavier component of IMBHBs can be constrained with an uncertainty of ∼10%–40% at a signal-to-noise ratio of 20. Focusing on the stellar-mass gap with new tabulations of the 12C(α, γ)16O reaction rate and its uncertainties, we evolve massive helium core stars using MESA to establish the lower and upper edges of the mass gap as ≃ 59 − 13 + 34 M ⊙ and ≃ 139 − 14 + 30 M ⊙ respectively, where the error bars give the mass range that follows from the ±3σ uncertainty in the 12C(α, γ)16O nuclear reaction rate. We find that high resolution of the tabulated reaction rate and fine temporal resolution are necessary to resolve the peak of the BH mass spectrum. We then study IMBHBs with components lying in the mass gap and show that the O4 run will be able to robustly identify most such systems. Finally, we reanalyze GW190521 with a state-of-the-art aligned-spin waveform model, finding that the primary mass lies in the mass gap with 90% credibility.

The gravitationally lensed star WHL 0137–LS, nicknamed Earendel, was identified with a photometric redshift z phot = 6.2 ± 0.1 based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera images of Earendel in eight filters spanning 0.8–5.0 μm. In these higher-resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to μ > 4000 and restricting the source plane radius further to r < 0.02 pc, or ∼4000 au. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system and support the previous photometric redshift estimate. Fitting grids of stellar spectra to our photometry yields a stellar temperature of T eff ≃ 13,000–16,000 K, assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from log ( L ) = 5.8 to 6.6 L ⊙, which is in the range where one expects luminous blue variable stars. Follow-up observations, including JWST NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe.

We explore the evolution of a select grid of solar metallicity stellar models from their pre-main sequence phase to near their final fates in a neutrino Hertzsprung-Russell diagram, where the neutrino luminosity replaces the traditional photon luminosity. Using a calibrated MESA solar model for the solar neutrino luminosity (L ν, = 0.02398 · L γ, = 9.1795 × 10 31 erg s −1 ) as a normalization, we identify 0.3 MeV electron neutrino emission from helium burning during the helium flash (peak L ν /L ν, 10 4 , flux Φ ν,He flash 170 (10 pc/d) 2 cm −2 s −1 for a star located at a distance of d parsec, timescale 3 days) and the thermal pulse (peak L ν /L ν, 10 9 , flux Φ ν,TP 1.7×10 7 (10 pc/d) 2 cm −2 s −1 , timescale 0.1 yr) phases of evolution in low mass stars as potential probes for stellar neutrino astronomy. We also delineate the contribution of neutrinos from nuclear reactions and thermal processes to the total neutrino loss along the stellar tracks in a neutrino Hertzsprung-Russell diagram. We find, broadly but with exceptions, that neutrinos from nuclear reactions dominate whenever hydrogen and helium burn, and that neutrinos from thermal processes dominate otherwise.

We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M ⊙ carbon–oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of 22Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6−7 significant figures of precision. Usually white dwarf models without 22Ne are fit to the observed periods and other properties. The rms residuals to the ≃150−400 s low-order g-mode periods are typically in the range of σ rms ≲ 0.3 s, for a fit precision of σ rms/P ≲ 0.3%. We find average relative period shifts of ΔP/P ≃ ±0.5% for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of ΔP/P depending on the specific g-mode, abundance of 22Ne, effective temperature, and the mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when 22Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of 22Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including 22Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model-fitting processes.

Gravitational-wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernovae (PISNe). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections of merging BHs. We use MESA to evolve single, nonrotating, massive helium cores with a metallicity of Z = 10−5, until they either collapse to form a BH or explode as a PISN, without leaving a compact remnant. We calculate the boundaries of the lower BH mass gap for S-factors in the range S(300 keV) = (77,203) keV b, corresponding to the ±3σ uncertainty in our high-resolution tabulated 12C(α,γ)16O reaction rate probability distribution function. We extensively test temporal and spatial resolutions for resolving the theoretical peak of the BH mass spectrum across the BH mass gap. We explore the convergence with respect to convective mixing and nuclear burning, finding that significant time resolution is needed to achieve convergence. We also test adopting a minimum diffusion coefficient to help lower-resolution models reach convergence. We establish a new lower edge of the upper mass gap as M lower ≃ 60 − 14 + 32 M ⊙ from the ±3σ uncertainty in the 12C(α, γ)16O rate. We explore the effect of a larger 3α rate on the lower edge of the upper mass gap, finding M lower ≃ 69 − 18 + 34 M ⊙. We compare our results with BHs reported in the Gravitational-Wave Transient Catalog.

Abstract.Los Alamos National Laboratory has calculated a new generation of radiative opacities (OPLIB data using the ATOMIC code) for elements with atomic number Z=1-30 with improved physics input, updated atomic data, and finer temperature grid to replace the Los Alamos LED-COP opacities released in the year 2000. We calculate the evolution of standard solar models including these new opacities, and compare with models evolved using the Lawrence Livermore National Laboratory OPAL opacities (Iglesias & Rogers 1996). We use the solar abundance mixture of Asplund et al. 2009. The Los Alamos ATOMIC opacities (Colgan et al. 2013a(Colgan et al. , 2013b have steeper opacity derivatives than those of OPAL for temperatures and densities of the solar interior radiative zone. We compare the calculated nonadiabatic solar oscillation frequencies and solar interior sound speed to observed frequencies and helioseismic inferences. The calculated sound-speed profiles are similar for models evolved using either the updated Iben evolution code (see Guzik & Mussack 2010), or the MESA evolution code (Paxton et al. 2015). The LANL ATOMIC opacities partially mitigate the 'solar abundance problem'.

We seek signatures of the current experimental 12C α , γ 16 O reaction rate probability distribution function in the pulsation periods of carbon–oxygen white dwarf (WD) models. We find that adiabatic g-modes trapped by the interior carbon-rich layer offer potentially useful signatures of this reaction rate probability distribution function. Probing the carbon-rich region is relevant because it forms during the evolution of low-mass stars under radiative helium-burning conditions, mitigating the impact of convective mixing processes. We make direct quantitative connections between the pulsation periods of the identified trapped g-modes in variable WD models and the current experimental 12C α , γ 16 O reaction rate probability distribution function. We find an average spread in relative period shifts of ΔP/P ≃ ±2% for the identified trapped g-modes over the ±3σ uncertainty in the 12C α , γ 16 O reaction rate probability distribution function—across the effective temperature range of observed DAV and DBV WDs and for different WD masses, helium shell masses, and hydrogen shell masses. The g-mode pulsation periods of observed WDs are typically given to six to seven significant figures of precision. This suggests that an astrophysical constraint on the 12C α , γ 16 O reaction rate could, in principle, be extractable from the period spectrum of observed variable WDs.

No abstract