Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BLACK HOLES”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

R -mode Stability of GW190814's Secondary Component as a Supermassive and Superfast Pulsar

The nature of GW190814's secondary component m 2 of mass (2.50–2.67)M ⊙ in the mass gap between the currently known maximum mass of neutron stars and the minimum mass of black holes is currently under hot debate. Among the many possibilities proposed in the literature, the m 2 was suggested as a superfast pulsar while its r-mode stability against the run-away gravitational radiation through the Chandrasekhar-Friedman-Schutz mechanism is still unknown. Previously, Fortin et al. constructed a sample of 33 unified equations of state using the same nuclear interactions from the crust to the core consistently; from that sample we use those equations that fulfill all currently known astrophysical and nuclear physics constraints to compare the minimum frequency required for m 2 to rotationally sustain a mass greater than 2.50 M ⊙ with the critical frequency above which the r-mode instability occurs. Here, we use two extreme damping models assuming that the crust is either perfectly rigid or elastic. Using the stability of 19 observed low-mass X-ray binaries as an indication that the rigid crust damping of the r-mode dominates within the models studied, we find that m 2 is r-mode-stable while rotating with a frequency higher than 870.2 Hz (0.744 times its Kepler frequency of 1169.6 Hz) as long as its temperature is lower than about 3.9 × 10 7 K, further supporting the proposal that GW190814's secondary component is a supermassive and superfast pulsar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Fast Dissipation of Colliding Alfvén Waves in a Magnetically Dominated Plasma

Magnetic energy around compact objects often dominates over plasma rest mass, and its dissipation can power the object's luminosity. We describe a dissipation mechanism that works faster than magnetic reconnection. The mechanism involves two strong Alfvén waves with anti-aligned magnetic fields B 1 and B 2 that propagate in opposite directions along the background magnetic field B 0 and collide. The collision forms a thin current sheet perpendicular to B 0 , which absorbs the incoming waves. The current sheet is sustained by an electric field E breaking the magnetohydrodynamic condition E < B and accelerating particles to high energies. We demonstrate this mechanism with kinetic plasma simulations using a simple setup of two symmetric plane waves with amplitude A = B 1 / B 0 = B 2 /B 0 propagating in a uniform B 0 . The mechanism is activated when A > 1/2. It dissipates a large fraction of the wave energy, † = (2A - 1)/A 2 , reaching 100% when A = 1. The plane geometry allows one to see the dissipation process in a one-dimensional simulation. We also perform two-dimensional simulations, enabling spontaneous breaking of the plane symmetry by the tearing instability of the current sheet. At moderate A of main interest, the tearing instability is suppressed. Dissipation transitions to normal, slower, magnetic reconnection at A $\gg$ 1. The fast dissipation described in this paper may occur in various objects with perturbed magnetic fields, including magnetars, jets from accreting black holes, and pulsar wind nebulae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Collapsar R-process Yields Can Reproduce [Eu/Fe] Abundance Scatter in Metal-poor Stars

It is unclear if neutron star mergers can explain the observed r-process abundances of metal-poor stars. Collapsars, defined here as rotating massive stars whose collapse results in a rapidly accreting disk around a black hole that can launch jets, are a promising alternative. We find that we can produce a self-consistent model in which a population of collapsars with stochastic europium yields synthesizes all of the r-process material in metal-poor ([Fe/H] < - 2.5) stars. Our model reproduces the observed scatter and evolution of scatter of [Eu/Fe] abundances. We find that if collapsars are the dominant r-process site for metal-poor stars, r-process synthesis may be linked to supernovae that produce long γ-ray bursts. Our results also allow for the possibility that core-collapse supernovae beyond those that launch γ-ray bursts also produce r-process material (e.g., potentially a subset of Type Ic-BL supernovae). Furthermore, we identify collapsar jet properties (isotropic energy, engine luminosity, or engine time) that may trace r-process yield and verify that the amount of r-process yield produced per collapsar in our model ( ~ 0.07M⊙) is consistent with other independent estimates. In the future, achieving 0.05 dex precision on distribution scatter or a reliable selection function would further constrain our probe of r-process production. Our model would also hold for another prompt r-process site with a power-law yield, and work is needed to determine if, for example, fast-merging neutron stars can also explain abundance scatter.

Astronomy & Astrophysics↗
DOE OSTI2021

The Nature of Jets in Double-peaked Emission-line AGN in the KISSR Sample

Double-peaked emission-line AGN (DPAGN) have been regarded as binary black hole candidates. We present here results from parsec-scale radio observations with the Very Long Baseline Array (VLBA) of five DPAGN belonging to the KISSR sample of emission-line galaxies. This work concludes our pilot study of nine type 2 Seyfert and LINER DPAGN from the KISSR sample. In the nine sources, dual compact cores are only detected in the “offset AGN,” KISSR 102. However, the overall incidence of jets in the eight sources that were detected with the VLBA is ≥ 60%. We find a difference in the “missing flux density” going from the Very Large Array to VLBA scales between Seyferts and LINERs, with LINERs showing less missing flux density on parsec-scales. Using the emission-line modeling code, MAPPINGS III, we find that the emission lines are likely to be influenced by jets in 5/9 sources. Furthermore, jet-medium interaction is the likely cause of the emission-line splitting observed in the SDSS spectra of these sources. Jets in radio-quiet AGN are therefore energetically capable of influencing their parsec- and kpcscale environments, which makes them agents of “radio AGN feedback” (similar to radio-loud AGN).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Igniting Weak Interactions in Neutron Star Postmerger Accretion Disks

The merger of two neutron stars or a neutron star and a black hole typically results in the formation of a postmerger accretion disk. Outflows from disks may dominate the overall ejecta from mergers and be a major source of r-process nuclei in our universe. We explore the parameter space of such disks and their outflows and r-process yields by performing 3D general-relativistic magnetohydrodynamic simulations with weak interactions and approximate neutrino transport. We discuss the mapping between the initial binary parameters and the parameter space of the resulting disks, chiefly characterized by their initial accretion rate. We demonstrate the existence of an ignition threshold for weak interactions at around ~10 -3 M⊙ s -1 for typical parameters by means of analytic calculations and numerical simulations. We find a degenerate, self-regulated, neutrino-cooled regime above the threshold and an advection-dominated regime below the threshold. Excess heating in the absence of neutrino cooling below the threshold leads to ≳60% of the initial disk mass being ejected in outflows, with typical velocities of ~(0.1–0.2)c, compared to ≲40% at ~(0.1–0.15)c above the threshold. While disks below the threshold show suppressed production of light r-process elements, disks above the threshold can produce the entire range of r-process elements, in good agreement with the observed solar system abundances. Disks below the ignition threshold may produce an overabundance of actinides seen in actinide-boost stars. As gravitational-wave detectors start to sample the neutron star merger parameter space, different disk realizations may be observable via their associated kilonova emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Equation-of-state Dependence of Gravitational Waves in Core-collapse Supernovae

Gravitational waves (GWs) provide unobscured insight into the birthplace of neutron stars and black holes in core-collapse supernovae (CCSNe). The nuclear equation of state (EOS) describing these dense environments is yet uncertain, and variations in its prescription affect the proto-neutron star (PNS) and the post-bounce dynamics in CCSN simulations, subsequently impacting the GW emission. In this study, we perform axisymmetric simulations of CCSNe with Skyrme-type EOSs to study how the GW signal and PNS convection zone are impacted by two experimentally accessible EOS parameters, (1) the effective mass of nucleons, m*, which is crucial in setting the thermal dependence of the EOS, and (2) the isoscalar incompressibility modulus, K sat . While K sat shows little impact, the peak frequency of the GWs has a strong effective mass dependence due to faster contraction of the PNS for higher values of m* owing to a decreased thermal pressure. These more compact PNSs also exhibit more neutrino heating, which drives earlier explosions and correlates with the GW amplitude via accretion plumes striking the PNS, exciting the oscillations. We investigate the spatial origin of the GWs and show the agreement between a frequency-radial distribution of the GW emission and a perturbation analysis. We do not rule out overshoot from below via PNS convection as another moderately strong excitation mechanism in our simulations. We also study the combined effect of effective mass and rotation. In all our simulations we find evidence for a power gap near ~1250 Hz; we investigate its origin and report its EOS dependence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. II. Connection with γ-Rays

Abstract It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multiwavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multiwavelength flares. Nevertheless, previous works have not directly evaluated γ -ray signatures from first-principles simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multiwavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and Fermi-LAT γ -ray bands as well as closely correlated optical and γ -ray flares. The swings in optical polarization angle are also accompanied by γ -ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self-Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fast γ -ray flares or orphan γ -ray flares under the leptonic scenario, complementary to the frequently considered minijet scenario. It may also imply neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Systematic Analysis of Stellar Populations in the Host Galaxies of Changing-look AGNs

“Changing-look” active galactic nuclei (CL-AGNs) are a newly discovered class of AGNs that show the appearance (or disappearance) of broad emission lines within short timescales (months to years),and are often associated with dramatic changes in their continuum emissions. They provide us with an unprecedented chance to directly investigate the host galaxy properties with minimal contamination from the luminous central engine during the turn-off state, which is difficult for normal luminous AGNs. In this work, for the first time, we systematically characterize the stellar populations and star formation histories of host galaxies for 26 turn-off CL-AGNs using the stellar population synthesis code STARLIGHT. We find that the stellar populations of CL-AGNs are similar to those of normal AGNs, except that the intermediate stellar populations contribute more fractions. We estimate their stellar velocity dispersions (σ * ) and black hole masses (M BH , vir ), and find that CL-AGNs also follow the overall M BH –σ * relationship. We also confirm the previous claims that CL-AGNs tend to be biased toward lower Eddington ratios, and that their extreme variabilities are more likely due to the intrinsic changes of the accretion rates. In addition, CL-AGNs with recent star formations tend to have higher Eddington ratios. Compared with previous studies, our analysis suggests that there may be a correlation between CL-AGN host galaxy properties and their CL phenomena.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Radio and X-Ray Observations of the Luminous Fast Blue Optical Transient AT 2020xnd

We present deep X-ray and radio observations of the fast blue optical transient (FBOT) AT 2020xnd/ZTF 20acigmel at z = 0.2433 from 13 days to 269 days after explosion. AT 2020xnd belongs to the category of optically luminous FBOTs with similarities to the archetypal event AT 2018cow. AT 2020xnd shows luminous radio emission reaching L ν ≈ 8 × 10 29 erg s -1 Hz -1 at 20 GHz and 75 days post-explosion, accompanied by luminous and rapidly fading soft X-ray emission peaking at L X ≈ 6 × 10 42 erg s -1 . Interpreting the radio emission in the context of synchrotron radiation from the explosion's shock interaction with the environment, we find that AT 2020xnd launched a high-velocity outflow (v ~ 0.1c–0.2c) propagating into a dense circumstellar medium (effective $\dot{M}\approx {10}^{-3}\,{M}_{\odot }$ yr -1 for an assumed wind velocity of v w = 1000 km s -1 ). Similar to AT 2018cow, the detected X-ray emission is in excess compared to the extrapolated synchrotron spectrum and constitutes a different emission component, possibly powered by accretion onto a newly formed black hole or neutron star. These properties make AT 2020xnd a high-redshift analog to AT 2018cow, and establish AT 2020xnd as the fourth member of the class of optically luminous FBOTs with luminous multiwavelength counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Properties of High-redshift Gamma-Ray Bursts

Abstract The immense power of gamma-ray bursts (GRBs) makes them ideal probes of the early universe. By using absorption lines in the afterglows of high-redshift GRBs, astronomers can study the evolution of metals in the early universe. With an understanding of the nature of GRB progenitors, the rate and properties of GRBs observed at high redshift can probe the star formation history and the initial mass function of stars at high redshift. This paper presents a detailed study of the dependence on metallicity and mass of the properties of long-duration GRBs under the black hole accretion disk paradigm to predict the evolution of these properties with redshift. These models are calibrated on the current GRB observations and then used to make predictions for new observations and new missions (e.g., the proposed Gamow mission) studying high-redshift GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Lick AGN Monitoring Project 2016: Dynamical Modeling of Velocity-resolved H β Lags in Luminous Seyfert Galaxies

We have modeled the velocity-resolved reverberation response of the Hβ broad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitoring Project 2016 sample, drawing inferences on the geometry and structure of the low-ionization broad-line region (BLR) and the mass of the central supermassive black hole. Overall, we find that the Hβ BLR is generally a thick disk viewed at low to moderate inclination angles. We combine our sample with prior studies and investigate line-profile shape dependence, such as ${\mathrm{log}}_{10}(\mathrm{FWHM}/\sigma )$, on BLR structure and kinematics and search for any BLR luminosity-dependent trends. We find marginal evidence for an anticorrelation between the profile shape of the broad Hβ emission line and the Eddington ratio, when using the rms spectrum. However, we do not find any luminosity-dependent trends, and conclude that AGNs have diverse BLR structure and kinematics, consistent with the hypothesis of transient AGN/BLR conditions rather than systematic trends.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Near-infrared Look at AGN Feedback in Bulgeless Galaxies

While it is generally believed that supermassive black holes (SMBHs) lie in most galaxies with bulges, few SMBHs have been confirmed in bulgeless galaxies. Identifying such a population could provide important insights to the BH seed population and secular BH growth. To this end, we obtained near-infrared (NIR) spectroscopic observations of a sample of low-redshift bulgeless galaxies with mid-infrared colors suggestive of active galactic nuclei (AGNs). We find additional evidence of AGN activity (such as coronal lines and broad permitted lines) in 69% (9/13) of the sample, demonstrating that mid-infrared selection is a powerful tool to detect AGNs. More than half of the galaxies with confirmed AGN activity show fast outflows in [O III ] in the optical and/or [Si vi] in the NIR, with the latter generally having much faster velocities that are also correlated to their spatial extent. We are also able to obtain virial BH masses for some targets and find they fall within the scatter of other late-type galaxies in the M BH –M stellar relation. The fact that they lack a significant bulge component indicates that secular processes, likely independent of major mergers, grew these BHs to supermassive sizes. Finally, we analyze the rotational gas kinematics and find two notable exceptions: two AGN hosts with outflows that appear to be rotating faster than expected. There is an indication that these two galaxies have stellar masses significantly lower than expected from their dark matter halo masses. This, combined with the observed AGN activity and strong gas outflows, may be evidence of the effects of AGN feedback.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Constraints on the Cosmic Expansion History from GWTC–3

We use 47 gravitational wave sources from the Third LIGO–Virgo–Kamioka Gravitational Wave Detector Gravitational Wave Transient Catalog (GWTC–3) to estimate the Hubble parameter H(z), including its current value, the Hubble constant H0. Each gravitational wave (GW) signal provides the luminosity distance to the source, and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog. Using the binary black hole (BBH) redshifted masses, we simultaneously infer the source mass distribution and H(z). The source mass distribution displays a peak around 34 M ⊙ , followed by a drop-off. Assuming this mass scale does not evolve with the redshift results in a H(z) measurement, yielding H 0 = $68$$^{+12}_{–8}$km s –1 Mpc –1 (68% credible interval) when combined with the H 0 measurement from GW170817 and its electromagnetic counterpart. This represents an improvement of 17% with respect to the H 0 estimate from GWTC–1. The second method associates each GW event with its probable host galaxy in the catalog GLADE+, statistically marginalizing over the redshifts of each event's potential hosts. Assuming a fixed BBH population, we estimate a value of H 0 = $68$$^{+8}_{–6}$km s –1 Mpc –1 with the galaxy catalog method, an improvement of 42% with respect to our GWTC–1 result and 20% with respect to recent H 0 studies using GWTC–2 events. However, we show that this result is strongly impacted by assumptions about the BBH source mass distribution; the only event which is not strongly impacted by such assumptions (and is thus informative about H 0 ) is the well-localized event GW190814.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Resolving the Inner Parsec of the Blazar J1924–2914 with the Event Horizon Telescope

The blazar J1924–2914 is a primary Event Horizon Telescope (EHT) calibrator for the Galactic center’s black hole Sagittarius A*. Here we present the first total and linearly polarized intensity images of this source obtained with the unprecedented 20 μas resolution of the EHT. J1924–2914 is a very compact flat-spectrum radio source with strong optical variability and polarization. In April 2017 the source was observed quasi-simultaneously with the EHT (April 5–11), the Global Millimeter VLBI Array (April 3), and the Very Long Baseline Array (April 28), giving a novel view of the source at four observing frequencies, 230, 86, 8.7, and 2.3 GHz. These observations probe jet properties from the subparsec to 100 pc scales. We combine the multifrequency images of J1924–2914 to study the source morphology. We find that the jet exhibits a characteristic bending, with a gradual clockwise rotation of the jet projected position angle of about 90° between 2.3 and 230 GHz. Linearly polarized intensity images of J1924–2914 with the extremely fine resolution of the EHT provide evidence for ordered toroidal magnetic fields in the blazar compact core.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Limits on X-Ray Polarization at the Core of Centaurus A as Observed with the Imaging X-Ray Polarimetry Explorer

We present measurements of the polarization of X-rays in the 2–8 keV band from the nucleus of the radio galaxy Centaurus A (Cen A), using a 100 ks observation from the Imaging X-ray Polarimetry Explorer (IXPE). Nearly simultaneous observations of Cen A were also taken with the Swift, NuSTAR, and INTEGRAL observatories. No statistically significant degree of polarization is detected with IXPE. These observations have a minimum detectable polarization at 99% confidence (MDP 99 ) of 6.5% using a weighted, spectral model-independent calculation in the 2–8 keV band. The polarization angle ψ is consequently unconstrained. Spectral fitting across three orders of magnitude in X-ray energy (0.3–400 keV) demonstrates that the SED of Cen A is well described by a simple power law with moderate intrinsic absorption (N H ~ 10 23 cm -2 ) and a Fe Kα emission line, although a second unabsorbed power law is required to account for the observed spectrum at energies below 2 keV. This spectrum suggests that the reprocessing material responsible for this emission line is optically thin and distant from the central black hole. Our upper limits on the X-ray polarization are consistent with the predictions of Compton scattering, although the specific seed photon population responsible for the production of the X-rays cannot be identified. The low polarization degree, variability in the core emission, and the relative lack of variability in the Fe Kα emission line support a picture where electrons are accelerated in a region of highly disordered magnetic fields surrounding the innermost jet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Angular-momentum Transport in Proto-neutron Stars and the Fate of Neutron Star Merger Remnants

Both the core collapse of rotating massive stars, and the coalescence of neutron star (NS) binaries result in the formation of a hot, differentially rotating NS remnant. The timescales over which differential rotation is removed by internal angular-momentum transport processes (viscosity) have key implications for the remnant’s long-term stability and the NS equation of state (EOS). Guided by a nonrotating model of a cooling proto-NS, we estimate the dominant sources of viscosity using an externally imposed angular-velocity profile Ω(r). Although the magneto-rotational instability provides the dominant source of effective viscosity at large radii, convection and/or the Tayler–Spruit dynamo dominate in the core of merger remnants where dΩ/dr ≥ 0. Furthermore, the viscous timescale in the remnant core is sufficiently short that solid-body rotation will be enforced faster than matter is accreted from rotationally supported outer layers. Guided by these results, we develop a toy model for how the merger remnant core grows in mass and angular momentum due to accretion. We find that merger remnants with sufficiently massive and slowly rotating initial cores may collapse to black holes via envelope accretion, even when the total remnant mass is less than the usually considered threshold ≈1.2 M TOV for forming a stable solid-body rotating NS remnant (where M TOV is the maximum nonrotating NS mass supported by the EOS). This qualitatively new picture of the post-merger remnant evolution and stability criterion has important implications for the expected electromagnetic counterparts from binary NS mergers and for multimessenger constraints on the NS EOS

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Standard Siren Measurement of the Hubble Constant Using Gravitational-wave Events from the First Three LIGO/Virgo Observing Runs and the DESI Legacy Survey

We present a new constraint on the Hubble constant H 0 using a sample of well-localized gravitational-wave (GW) events detected during the first three LIGO/Virgo observing runs as dark standard sirens. In the case of dark standard sirens, a unique host galaxy is not identified, and the redshift information comes from the distribution of potential host galaxies. From the third LIGO/Virgo observing run detections, we add the asymmetric-mass binary black hole GW190412 and the high-confidence GW candidates S191204r, S200129m, and S200311bg to the sample of dark standard sirens analyzed in Palmese et al. Our sample contains the top 20% (based on localization) GW events and candidates to date with significant coverage by the Dark Energy Spectroscopic Instrument Legacy Survey. We combine the H 0 posterior for eight dark siren events, finding H 0 = $79.8^{+19.1}_{-12.8}$ km s -1 Mpc -1 (68% highest density interval) for a prior in H 0 uniform between [20, 140] km s -1 Mpc -1 . This result shows that a combination of eight well-localized dark sirens combined with an appropriate galaxy catalog is able to provide an H 0 constraint that is competitive (~20% versus 18% precision) with a single bright standard siren analysis (i.e., assuming the electromagnetic counterpart) using GW170817. When combining the posterior with that from GW170817, we obtain H 0 = $72.77^{+11.0}_{-7.55}$ km s -1 Mpc -1 . This result is broadly consistent with recent H 0 estimates from both the cosmic microwave background and supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗