SEARCH · Engineering Papers

Results for “BINARY STARS”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2020

The Pre-He White Dwarfs in Eclipsing Binaries. I. WASP 0131+28

We report the first BV light curves and high-resolution spectra of the post-mass transfer binary star WASP 0131+28 to study the absolute properties of extremely low-mass white dwarfs. From the observed spectra, the double-lined radial velocities were derived, and the effective temperature and rotational velocity of the brighter, more massive primary were found to be T{sub eff,1}=10,000±200 K and v{sub 1}sin i = 55 ± 10 km s{sup −1}, respectively. The combined analysis of the Transiting Exoplanet Survey Satellite (TESS) archive data and ours yielded the accurate fundamental parameters of the program target. The masses were derived to about 1.0% accuracy and the radii to 0.6%, or better. The secondary component's parameters of M{sub 2}=0.200±0.002 M{sub ⊙}, R{sub 2}=0.528±0.003 R{sub ⊙}, T{sub eff,2} = 11,186 ± 235 K, and L{sub 2}=3.9±0.3 L{sub ⊙} are in excellent agreement with the evolutionary sequence for a helium-core white dwarf of mass 0.203 M{sub ⊙}, and indicates that this star is halfway through the constant luminosity phase. The results presented in this article demonstrate that WASP 0131+28 is an EL CVn eclipsing binary in a thin disk, which is formed from the stable Roche-lobe overflow channel and composed of a main-sequence dwarf with a spectral type A0 and a pre-He white dwarf.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

High tide: a systematic search for ellipsoidal variables in ASAS-SN

ABSTRACT The majority of non-merging stellar mass black holes are discovered by observing high energy emission from accretion processes. Here, we pursue the large, but still mostly unstudied population of non-interacting black holes and neutron stars by searching for the tidally induced ellipsoidal variability of their stellar companions. We start from a sample of about 200 000 rotational variables, semiregular variables, and eclipsing binary stars from the All-Sky Automated Survey for Supernovae. We use a χ2 ratio test followed by visual inspection to identify 369 candidates for ellipsoidal variability. We also discuss how to combine the amplitude of the variability with mass and radius estimates for observed stars to calculate a minimum companion mass, identifying the most promising candidates for high mass companions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Birth of a Be star: an APOGEE search for Be stars forming through binary mass transfer

ABSTRACT Motivated by recent suggestions that many Be stars form through binary mass transfer, we searched the APOGEE survey for Be stars with bloated, stripped companions. From a well-defined parent sample of 297 Be stars, we identified one mass-transfer binary, HD 15124. The object consists of a main-sequence Be star ($M_{\rm Be}=5.3\pm 0.6\, {\rm M}_{\odot }$) with a low-mass ($M_{\rm donor}=0.92\pm 0.22\, {\rm M}_{\odot }$), subgiant companion on a 5.47-d orbit. The emission lines originate in an accretion disc caused by ongoing mass transfer, not from a decretion disc as in classical Be stars. Both stars have surface abundances bearing imprint of CNO processing in the donor’s core: the surface helium fraction is YHe ≈ 0.6, and the nitrogen-to-carbon ratio is 1000 times the solar value. The system’s properties are well-matched by binary evolution models in which mass transfer begins while a $3-5\, {\rm M}_{\odot }$ donor leaves the main sequence, with the originally less massive component becoming the Be star. These models predict that the system will soon become a detached Be + stripped star binary like HR 6819 and LB-1, with the stripped donor eventually contracting to become a core helium-burning sdO/B star. Discovery of one object in this short-lived (∼1 Myr) evolutionary phase implies the existence of many more that have already passed through it and are now Be + sdO/B binaries. We infer that $(10-60)\, {{\ \rm per\ cent}}$ of Be stars have stripped companions, most of which are $\sim 100\, \times$ fainter than the Be stars in the optical. Together with the dearth of main-sequence companions to Be stars and recent discovery of numerous Be + sdO/B binaries in the UV, our results imply that binarity plays an important role in the formation of Be stars.

El-Badry, Kareem (ORCID:0000000268711752)↗

DOE OSTI2023

Spectroscopic Analysis of Milky Way Outer Halo Satellites: Aquarius II and Boötes II

In this paper, we present a chemical and kinematic analysis of two ultrafaint dwarf galaxies (UFDs), Aquarius II (Aqu II) and Boötes II (Boo II), using Magellan/IMACS spectroscopy. We present the largest sample of member stars for Boo II (12), and the largest sample of red giant branch members with metallicity measurements for Aqu II (eight). In both UFDs, over 80% of targets selected based on Gaia proper motions turned out to be spectroscopic members. In order to maximize the accuracy of stellar kinematic measurements, we remove the identified binary stars and RR Lyrae variables. For Aqu II, we measure a systemic velocity of -65.3 ± 1.8 km s-1 and a metallicity of [Fe/H] = $-2.57^{+0.17}_{-0.17}$. When compared with previous measurements, these values display a ~6 km s -1 difference in radial velocity and a decrease of 0.27 dex in metallicity. Similarly for Boo II, we measure a systemic velocity of $-130.4^{+1.4}_{-1.1}$ km s -1 , more than 10 km s -1 different from the literature, a metallicity almost 1 dex smaller at [Fe/H] = $-2.71^{+0.11}_{-0.10}$, and a velocity dispersion 3 times smaller at σ v hel = $2.9^{+1.6}_{-1.2}$km s -1 . Additionally, we derive systemic proper-motion parameters and model the orbits of both UFDs. Finally, we highlight the extremely dark-matter-dominated nature of Aqu II and compute the J-factor for both galaxies to aid searches of dark matter annihilation. Despite the small size and close proximity of Boo II, it is an intermediate target for the indirect detection of dark matter annihilation due to its low-velocity dispersion and corresponding low dark matter density.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Statistical estimates of the binary properties of rotational variables

ABSTRACT We present a model to estimate the average primary masses, companion mass ranges, the inclination limit for recognizing a rotational variable, and the primary mass spreads for populations of binary stars. The model fits a population’s binary mass function distribution and allows for a probability that some mass functions are incorrectly estimated. Using tests with synthetic data, we assess the model’s sensitivity to each parameter, finding that we are most sensitive to the average primary mass and the minimum companion mass, with less sensitivity to the inclination limit and little to no sensitivity to the primary mass spread. We apply the model to five populations of binary spotted rotational variables identified in ASAS-SN, computing their binary mass functions using RV data from APOGEE. Their average primary mass estimates are consistent with our expectations based on their CMD locations ($\sim 0.75 \, {\rm M}_{\odot }$ for lower main sequence primaries and $\sim 0.9$–$1.2 \, {\rm M}_{\odot }$ for RS CVn and sub-subgiants). Their companion mass range estimates allow companion masses down to $M_2/M_1\simeq 0.1$, although the main sequence population may have a higher minimum mass fraction ($\sim 0.4$). We see weak evidence of an inclination limit $\gtrsim 50^{\circ }$ for the main sequence and sub-subgiant groups and no evidence of an inclination limit in the other groups. No groups show strong evidence for a preferred primary mass spread. We conclude by demonstrating that the approach will provide significantly better estimates of the primary mass and the minimum mass ratio and reasonable sensitivity to the inclination limit with 10 times as many systems.

Phillips, Anya (ORCID:000900051914974X)↗

DOE OSTI2021

Flaring and Spot Activities on the Semi-detached Binary System KIC 06852488

The continuous and unbroken time-series photometry data of Kepler and the Transiting Exoplanet Survey Satellite (TESS) provide a good opportunity to study the continuous variations in the light curve of binary stars. This paper reports the discoveries of the correlation between the two maxima on the light curve of KIC 06852488 and the correlation between the O’Connell effect and the O – C curve. The Max I and Max II vary with a same cycle length of ~2000 days and a 180° phase difference, and the variation of the Max II coincides with the O – C curve of the primary light minima. After analyzing the Kepler and TESS light curves, it is detected that this binary is a semi-detached system with a mass ratio of 0.4622(±0.0036). The secondary component is filling its critical Roche lobe, while the primary component is only filling 12.25(±0.10)% of the lobe volume that contributes more than 85% luminosity of the total system. The variation of the O’Connell effect could be explained by an evolving hot spot on the primary component and an evolving cool spot on the secondary component. It is found that their positions are symmetrical with the inner Lagrange L1 point. After subtracting the binary brightness changes, six optical flares are detected that come from the cool secondary. We calculated the energy of flares and found that all of the flares are superflares. Both the flares and an evolving dark spot on the lobe-filling secondary component make the system a very interesting source to investigate the effects of binary interaction on the magnetic activities.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Decaying Orbit of the Hot Jupiter WASP-12b: Confirmation with TESS Observations

Theory suggests that the orbits of some close-in giant planets should decay due to tidal interactions with their host stars. To date, WASP-12b is the only hot Jupiter reported to have a decaying orbit, at a rate of 29 ± 2 ms yr{sup −1}. We analyzed data from NASA’s Transiting Exoplanet Survey Satellite (TESS) to verify that WASP-12b’s orbit is indeed changing. We find that the TESS transit and occultation data are consistent with a decaying orbit with an updated period of 1.091420090 ± 0.000000041 days and a decay rate of 32.53 ± 1.62 ms yr{sup −1}. We find an orbital decay timescale of τ=P/| P-dot |=2.90±0.14 Myr. If the observed decay results from tidal dissipation, the modified tidal quality factor is Q{sub ⋆}{sup ′}=1.39±0.15×10{sup 5}, which falls at the lower end of values derived for binary star systems and hot Jupiters. Our result highlights the power of space-based photometry for investigating the orbital evolution of short-period exoplanets.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

HATS-47b, HATS-48Ab, HATS-49b, and HATS-72b: Four Warm Giant Planets Transiting K Dwarfs

We report the discovery of four transiting giant planets around K dwarfs. The planets HATS-47b, HATS-48Ab, HATS-49b, and HATS-72b have masses of 0.369{sub −0.021}{sup +0.031} M{sub J}, 0.243{sub −0.030}{sup +0.022} M{sub J}, 0.353{sub −0.027}{sup +0.038} M{sub J}, and 0.1254±0.0039 M{sub J}, respectively, and radii of 1.117±0.014 R{sub J}, 0.800±0.015 R{sub J}, 0.765±0.013 R{sub J}, and 0.7224±0.0032 R{sub J}, respectively. The planets orbit close to their host stars with orbital periods of 3.9228 days, 3.1317 days, 4.1480 days, and 7.3279 days, respectively. The hosts are main-sequence K dwarfs with masses of 0.674{sub −0.012}{sup +0.016} M{sub ⊙}, 0.7279 ±0.0066 M{sub ⊙}, 0.7133±0.0075 M{sub ⊙}, and 0.7311±0.0028, and with V-band magnitudes of V=14.829±0.010, 14.35±0.11, 14.998±0.040 and 12.469±0.010. The super-Neptune HATS-72b (a.k.a. WASP-191b and TOI 294.01) was independently identified as a transiting planet candidate by the HATSouth, WASP, and TESS surveys, and we present a combined analysis of all of the data gathered by each of these projects (and their follow-up programs). An exceptionally precise mass is measured for HATS-72b thanks to high-precision radial velocity (RV) measurements obtained with VLT/ESPRESSO, FEROS, HARPS, and Magellan/PFS. We also incorporate TESS observations of the warm Saturn–hosting systems HATS-47 (a.k.a. TOI 1073.01), HATS-48A, and HATS-49. HATS-47 was independently identified as a candidate by the TESS team, while the other two systems were not previously identified from the TESS data. The RV orbital variations are measured for these systems using Magellan/PFS. HATS-48A has a resolved [[Formula]] neighbor in Gaia DR2, which is a common-proper-motion binary star companion to HATS-48A with a mass of 0.22 M{sub ⊙} and a current projected physical separation of ∼1400 au.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Excitation and Excavation of the Claws of the Southern Crab

We show that the Southern Crab (aka Hen2–104) presents an auspicious opportunity to study the form and speed of the invisible winds that excavate and shock the lobes of various types of bipolar nebulae associated with close and highly evolved binary stars. A deep three-color image overlay of Hen2–104 reveals that the ionization state of its lobe edges, or “claws,” increases steadily from singly to doubly ionized values with increasing wall latitude. This “reverse” ionization pattern is unique among planetary nebulae (and similar objects) and incompatible with UV photoionization from a central source. We show that the most self-consistent explanation for the ionization pattern is shock ionization by a fast (~600 km s -1 ) “tapered” stellar wind in which the speed and momentum flux of the wind increase with equatorial latitude. We present a hydrodynamic simulation that places the latitude-dependent form, the knotty walls, and the reverse ionization of the outer lobes of Hen2–104 into a unified context.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Visual Orbits of Spectroscopic Binaries with the CHARA Array. III. HD 8374 and HD 24546

We present the visual orbits of two long-period spectroscopic binary stars, HD 8374 and HD 24546, using interferometric observations acquired with the CHARA Array and the Palomar Testbed Interferometer. We also obtained new radial velocities from echelle spectra using the APO 3.5 m and Fairborn 2.0 m telescopes. By combining the visual and spectroscopic observations, we solve for the full, three-dimensional orbits and determine the stellar masses and distances to within 3% uncertainty. We then estimate the effective temperature and radius of each component star through Doppler tomography and spectral energy distribution analyses, in order to compare the observed stellar parameters to the predictions of stellar evolution models. For HD 8374, we find masses of M {sub 1} = 1.636 ± 0.050M {sub ⊙} and M {sub 2} = 1.587 ± 0.049M {sub ⊙}, radii of R {sub 1} = 1.84 ± 0.05R {sub ⊙} and R {sub 2} = 1.66 ± 0.12R {sub ⊙}, temperatures of T{sub eff1}=7280±110 K and T{sub eff2}=7280±120 K, and an estimated age of 1.0 Gyr. For HD 24546, we find masses of M {sub 1} = 1.434 ± 0.014M {sub ⊙} and M {sub 2} = 1.409 ± 0.014M {sub ⊙}, radii of R {sub 1} = 1.67 ± 0.06R {sub ⊙} and R {sub 2} = 1.60 ± 0.10R {sub ⊙}, temperatures of T{sub eff1}=6790±120 K and T{sub eff2}=6770±90 K, and an estimated age of 1.4 Gyr. HD 24546 is therefore too old to be a member of the Hyades cluster, despite its physical proximity to the group.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Tidal Circularization of Binaries by Resonance Locking. I. The Importance of the Pre-main Sequence

Although tidal dissipation in binary stars has been studied for over a century, theoretical predictions have yet to match the observed properties of binary populations. This work quantitatively examines the recent proposal of tidal circularization by resonance locking, where tidal dissipation arises from resonances between the star’s natural oscillation frequencies and harmonics of the orbital frequency, and where resonances are “locked” for an extended period of time due to concurrent stellar evolution. We focus on tidal resonances with axi-symmetric gravity modes, and examine binaries with primary masses from one to two solar masses. We find that orbital evolution via resonance locking occurs primarily during the star’s pre-main-sequence phase, with the main-sequence phase contributing negligibly. Resonance locking, ignoring nonlinearity, can circularize binaries with peri-center distances out to ∼10 stellar radii, corresponding to circular periods of ∼4–6 days. However, we find resonantly excited gravity modes will become nonlinear in stellar cores, which prevents them from reaching their full, linear amplitudes. We estimate that such a “saturated resonance lock” reduces the circularization period by about a third, but resonance locking remains much more effective than the cumulative actions of equilibrium tides. In a companion paper, we examine recent binary data to compare against theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Carbon-enhanced metal-poor star candidates from BP/RP spectra in Gaia DR3

ABSTRACT Carbon-enhanced metal-poor (CEMP) stars comprise almost a third of stars with [Fe/H] < −2, although their origins are still poorly understood. It is highly likely that one sub-class (CEMP-s stars) is tied to mass-transfer events in binary stars, while another sub-class (CEMP-no stars) are enriched by the nucleosynthetic yields of the first generations of stars. Previous studies of CEMP stars have primarily concentrated on the Galactic halo, but more recently they have also been detected in the thick disc and bulge components of the Milky Way. Gaia DR3 has provided an unprecedented sample of over 200 million low-resolution (R ≈ 50) spectra from the BP and RP photometers. Training on the CEMP catalogue from the SDSS/SEGUE database, we use XGBoost to identify the largest all-sky sample of CEMP candidate stars to date. In total, we find 58 872 CEMP star candidates, with an estimated contamination rate of 12 per cent. When comparing to literature high-resolution catalogues, we positively identify 60–68 per cent of the CEMP stars in the data, validating our results and indicating a high completeness rate. Our final catalogue of CEMP candidates spans from the inner to outer Milky Way, with distances as close as r ∼ 0.8 kpc from the Galactic centre, and as far as r > 30 kpc. Future higher resolution spectroscopic follow-up of these candidates will provide validations of their classification and enable investigations of the frequency of CEMP-s and CEMP-no stars throughout the Galaxy, to further constrain the nature of their progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

XRISM view of a stellar flare: High-resolution Fe K spectra of HR 1099, an RS CVn-type star

A high-resolution X-ray spectroscopic observation was made of the RS CVn-type binary star HR 1099 using the Resolve instrument onboard XRISM for its calibration purposes. During the $\sim$400 ks telescope time covering 1.5 binary orbit, a flare lasting for $\sim$100 ks was observed with a released X-ray radiation energy of ${\sim }10^{34}$ erg, making it the first stellar flare ever observed with an X-ray microcalorimeter spectrometer. The flare peak count rate is 6.4 times higher than that in quiescence and is distinguished clearly in time thanks to the long telescope time. Many emission lines were detected in the 1.7–10 keV range both in the flare and quiescent phases. Using the high spectral resolution of Resolve in the Fe K band (6.5–7.0 keV), we resolved the inner-shell lines of Fe XIX-XXVI as well as the outer-shell lines of Fe XXV-XXVI . These lines have peaks in the contribution functions at different temperatures over a wide range, allowing us to construct the differential emission measure (DEM) distribution over the electron temperature of 1–10 keV (roughly 10–100 MK) based only on Fe lines, thus without an assumption of the elemental abundance. The reconstructed DEM has a bimodal distribution, and only the hotter component increased during the flare. The elemental abundance was derived based on the DEM distribution thus constructed. A significant abundance increase was observed during the flare for Ca and Fe, which are some of the elements with the lowest first ionization potential among those analyzed, but not for Si, S, and Ar. This behavior is seen in some giant solar flares and the present result is a clear example in stellar flares.

Astronomy and AstroPhysics↗

DOE OSTI2022

The SDSS-HET Survey of Kepler Eclipsing Binaries. A Sample of Four Benchmark Binaries

The purpose of this work is to extend a sample of accurately modeled, benchmark-grade eclipsing binaries (EBs) with accurately determined masses and radii. We select four "well-behaved" Kepler binaries, KIC 2306740, KIC 4076952, KIC 5193386 and KIC 5288543, each with at least eight double-lined spectra from the Apache Point Observatory Galactic Evolution Experiment instrument that is part of the Sloan Digital Sky Surveys III and IV, and from the Hobby–Eberly High Resolution Spectrograph. We obtain masses and radii with uncertainties of 2.5% or less for all four systems. Three of these systems have orbital periods longer than 9 days, and thus populate an undersampled region of the parameter space for extremely well-characterized detached EBs. We compare the derived masses and radii against MESA MIST isochrones to determine the ages of the systems. All systems were found to be coeval, showing that the results are consistent across MESA MIST and PHOEBE .

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Two-epoch Orbit Estimation for Wide Binaries Resolved in Hipparcos and Gaia

The Hipparcos catalog and its Double and Multiple System Annex (DMSA) lists 4099 components with individual proper motions and coordinates on the epoch 1991.25. Many of these long-period binary stars are also present in Gaia Data Release 2 (DR2). Using the available relative positions and proper motions separated by 24.25 yr, the equations of relative orbital motion can be solved for the two epoch eccentric anomalies, orbital period, and eccentricity. This method employs elimination of the linear Thiele–Innes unknowns and nonlinear optimization of the remaining condition equations. The quality of these solutions is compromised by the insufficient condition and modest precision of the Hipparcos astrometric data, as revealed by Monte Carlo simulations with artificially perturbed data points. The presence of multiple systems and optical pairs can also perturb the results. Limited experiments with artificial data indicate that useful estimates can be obtained with a 25 yr epoch difference for wide binaries with orbital periods up to ~500 yr. The prospects of this method dramatically improve with the proposed next-generation space astrometry missions such as Gaia-NIR and Theia, especially when additional conditions are included from astrometric or spectroscopic measurements. An ancillary catalog of cross-identification and astrometric information for 1295 double-star pairs cross-matched in Gaia DR2 and Hipparcos is also published.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Estimating outflow masses and velocities in merger simulations: Impact of r -process heating and neutrino cooling

Here, the determination of the mass, composition, and geometry of matter outflows in black hole-neutron star and neutron star-neutron star binaries is crucial to current efforts to model kilonovae and to understand the role of neutron star merger in r-process nucleosynthesis. In this manuscript, we review the simple criteria currently used in merger simulations to determine whether matter is unbound and what the asymptotic velocity of ejected material will be. We then show that properly accounting for both heating and cooling during r-process nucleosynthesis is important to accurately predict the mass and kinetic energy of the outflows. These processes are also likely to be crucial to predict the fall-back timescale of any bound ejecta. We derive a model for the asymptotic velocity of unbound matter and binding energy of bound matter that accounts for both of these effects and that can easily be implemented in merger simulations. We show, however, that the detailed velocity distribution and geometry of the outflows can currently only be captured by full three-dimensional fluid simulations of the outflows, as nonlocal effect ignored by the simple criteria used in merger simulations cannot be safely neglected when modeling these effects. Finally, we propose the introduction of simple source terms in the fluid equations to approximately account for heating/cooling from r-process nucleosynthesis in future seconds-long three-dimensional simulations of merger remnants, without the explicit inclusion of out-of-nuclear statistical equilibrium reactions in the simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Stellar kinematics of dwarf galaxies from multi-epoch spectroscopy: application to Triangulum II

ABSTRACT We present new MMT/Hectochelle spectroscopic measurements for 257 stars observed along the line of sight to the ultrafaint dwarf galaxy Triangulum II (Tri II). Combining results from previous Keck/DEIMOS spectroscopy, we obtain a sample that includes 16 likely members of Tri II, with up to 10 independent redshift measurements per star. To this multi-epoch kinematic data set, we apply methodology that we develop in order to infer binary orbital parameters from sparsely sampled radial velocity curves with as few as two epochs. For a previously identified (spatially unresolved) binary system in Tri II, we infer an orbital solution with period $296.0_{-3.3}^{+3.8} \rm ~ d$, semimajor axis $1.12^{+0.41}_{-0.24}\rm ~au$, and systemic velocity $-380.0 \pm 1.7 \rm ~km ~s^{-1}$ that we then use in the analysis of Tri II’s internal kinematics. Despite this improvement in the modelling of binary star systems, the current data remain insufficient to resolve the velocity dispersion of Tri II. We instead find a 95 per cent confidence upper limit of $\sigma _{v} \lesssim 3.4 \rm ~km~s^{-1}$.

Buttry, Rachel (ORCID:0000000155017609)↗

DOE OSTI2023

Jetted and Turbulent Stellar Deaths: New LVK-detectable Gravitational-wave Sources

Upcoming LIGO–Virgo–KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events that come from black hole and neutron star binary mergers. Detection of noninspiral GW sources is also anticipated. We report the discovery of a new class of noninspiral GW sources—the end states of massive stars—that can produce the brightest simulated stochastic GW burst signal in the LVK bands known to date, and could be detectable in LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble—cocoon—inside of the star. We simulate such a system using state-of-the-art 3D general relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, ~10–100 Hz, over a characteristic jet activity timescale ~10–100 s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third-generation GW detectors with an estimated rate of ~10 events yr –1 . The detection rate drops to ~1% of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multimessenger events.

79 ASTRONOMY AND ASTROPHYSICS↗