SEARCH · Engineering Papers

Results for “Transient sources”

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.

699 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Modelling populations of kilonovae

Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Kilonovae Across the Nuclear Physics Landscape: The Impact of Nuclear Physics Uncertainties on r-process-powered Emission

Merging neutron stars produce “kilonovae”—electromagnetic transients powered by the decay of unstable nuclei synthesized via rapid neutron capture (the r-process) in material that is gravitationally unbound during inspiral and coalescence. Kilonova emission, if accurately interpreted, can be used to characterize the masses and compositions of merger-driven outflows, helping to resolve a long-standing debate about the origins of r-process material in the Universe. We explore how the uncertain properties of nuclei involved in the r-process complicate the inference of outflow properties from kilonova observations. Using r-process simulations, we show how nuclear physics uncertainties impact predictions of radioactive heating and element synthesis. For a set of models that span a large range in both predicted heating and final abundances, we carry out detailed numerical calculations of decay product thermalization and radiation transport in a kilonova ejecta with a fixed mass and density profile. The light curves associated with our models exhibit great diversity in their luminosities, with peak brightness varying by more than an order of magnitude. We also find variability in the shape of the kilonova light curves and their color, which in some cases runs counter to the expectation that increasing levels of lanthanide and/or actinide enrichment will be correlated with longer, dimmer, redder emission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Enhanced production of 60 Fe in massive stars

Massive stars are a major source of chemical elements in the cosmos, ejecting freshly produced nuclei through winds and core-collapse supernova explosions into the interstellar medium. Among the material ejected, long-lived radioisotopes, such as 60 Fe (iron) and 26 Al (aluminum), offer unique signs of active nucleosynthesis in our galaxy. There is a long-standing discrepancy between the observed 60 Fe/ 26 Al ratio by γ-ray telescopes and predictions from supernova models. This discrepancy has been attributed to uncertainties in the nuclear reaction networks producing 60 Fe, and one reaction in particular, the neutron-capture on 59 Fe. Here we present experimental results that provide a strong constraint on this reaction. We use these results to show that the production of 60 Fe in massive stars is higher than previously thought, further increasing the discrepancy between observed and predicted 60 Fe/ 26 Al ratios. The persisting discrepancy can therefore not be attributed to nuclear uncertainties, and points to issues in massive-star models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Influence of substitution pattern on the dynamics of internal conversion and intersystem crossing in thiopyridone isomers

We report a combined experimental and theoretical investigation of the ultrafast internal conversion (IC) and intersystem crossing (ISC) dynamics of two thiopyridone (TP) isomers in solution. Our study used ultrafast transient X-ray absorption spectroscopy (XAS) at the sulfur K-edge, in conjunction with electronic excited state surface hopping molecular dynamics and simulations of the excited state XAS, to investigate the impact of the functional group substitution pattern and solvent on the dynamics of IC and ISC. The combination of the localized X-ray probe and the simulation results enables, in part, the differentiation between ππ* and nπ* character excited states, as well as singlet and triplet states. Access to nπ* character excitations has particular value since they often prove challenging to assess with optical spectroscopy. For 2-TP, the photoexcited S 2 (ππ*) state rapidly undergoes IC to the S 1 (nπ*) state below the instrument response time, followed by ISC to the T 1 (ππ*) state on a timescale of 600 fs in acetonitrile. For 4-TP, the timescale of S 2 to S 1 IC increases to 330 fs and the timescale of ISC increases to more than 10 ps. The differences between isomers are rationalized by considering the key role of the, nπ* intermediates in mediating the intersystem crossing of these systems. Varying the substitution pattern of the molecule can stabilize or destabilize these intermediates leading to the increase in ISC rate in the ortho isomer as compared to the para isomer, while changing the solvent from acetonitrile to water had minimal effect on the electronic excited state relaxation mechanism.

Garratt, Douglas [SLAC National Accelerator Labora↗

DOE OSTI2022

The role of supernova convection for the lower mass gap in the isolated binary formation of gravitational wave sources

ABSTRACT Understanding the astrophysical phenomena involving compact objects requires an insight about the engine behind the core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of the growing populations of the detected black hole–black hole, black hole–neutron star, and neutron star–neutron star mergers. The time-scale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. We adopt new formulas for the relation between the pre-SN star properties and their remnants and check how they impact the population of double compact object (DCO) mergers formed via the isolated binary evolution. The new formulas give one the ability to test a wide spectrum of assumptions on the convection growth time. In particular, the different variants allow for a smooth transition between having a deep mass gap and a remnant mass distribution filled by massive neutron stars and low-mass black holes. We present the distribution of masses, mass ratios, and the local merger rate densities of DCO for the different variants of new formulas and test them together with different approaches to other highly uncertain processes. We find that the mass distribution of DCO mergers is sensitive to the adopted assumption on the SN convection growth time-scale up to $m_1+m_2 \lesssim 35 \, \mathrm{M}_{\odot }$. Between the two extreme tested variants the probability of compact object formation within the mass gap may differ by up to approximately two orders of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Pursuing Heteroleptic Ligand Design Principles for Photoactive Fe Complexes with Ultrafast X-ray Emission and Variable-Temperature Optical Spectroscopies

Understanding the key parameters that govern the photophysical and photochemical properties of transition metal complexes is essential for the development of efficient photosensitizers for photocatalytic applications. Achieving this objective necessitates clear and detailed investigations of their electronic excited states, for which time-resolved metal Kβ X-ray emission spectroscopy (XES) has proven highly effective. Here, we present a time-resolved Fe Kβ XES study of a heteroleptic Fe(II) polypyridyl carbene complex, [Fe(phen) 2 (C 4 H 10 N 4 )] 2+ (1; phen = 1,10-phenanthroline), utilizing both the valence-to-core and Kβ mainline spectral regions, complemented by variable-temperature transient optical absorption (VT-TA) spectroscopy. Detailed analysis of the time-resolved Kβ XES data, supported by density functional theory (DFT) calculations and an Eyring analysis of the VT-TA data, reveals parallel excited state relaxation dynamics that support an assignment of the long-lived excited state to a triplet metal-centered state. Placing these results in the context of prior studies of heteroleptic Fe(II) polypyridyl cyanide complexes motivated a series of DFT calculations to investigate the effects of ligand structural flexibility and arrangement. These calculations reinforce the experimentally derived conclusion that constraining structural flexibility with multidentate ligands significantly impacts the excited state relaxation dynamics. Furthermore, our study emphasizes that the arrangement of strong field ligands in heteroleptic complexes substantially affects the energy of Jahn–Teller active triplet metal-centered states in low-spin d 6 metal complexes. Together, these findings provide synthetic design principles for extending metal-to-ligand charge transfer excited state lifetimes of heteroleptic Fe complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Purification and Concentration of 99 Mo from a UREX raffinate

The medical isotope 99m Tc ( t ½ = 6.0 h) accounts for over 80% of isotopes used in diagnostic nuclear medicine today and is normally derived from its transient equilibrium parent 99 Mo ( t ½ = 66 h). To reduce U.S. dependence on 99 Mo derived from HEU by foreign suppliers, the NNSA/M3 program—under the American Medical Isotope Production Act—has been tasked with facilitating the work of domestic 99 Mo suppliers that do not utilize HEU. Superconducting electron linear accelerators with high-Z converter targets can generate bremsstrahlung photons and neutron fluxes that can induce photonuclear reactions and LEU fission. After sufficient production intervals, targets can be rotated out and processed while another batch is irradiated. The process flow chemistry is shown in Figure 1. The irradiated U 3 O 8 targets are retrieved and dissolved in HNO 3 ; the volatile fission products are expelled and captured during this step. The HNO 3 liquor bearing the uranyl, fission and activation products is injected into the UREX suite of liquid-liquid extraction banks. The uranyl is partitioned in typical fashion using tri-$\textit{n}$-butyl phosphate (TBP) in a hydrocarbon diluent. The raffinate—comprising mainly 99 Mo and fission products—is then injected into the MoLLE (Molybdenum Liquid-Liquid Extraction) flowsheet, where the Mo is selectively extracted by an organophosphorous acid extractant such as di-(2-ethylhexyl phosphoric acid) (HDEHP). Trace amounts of Nb, Zr, Np, I, and Te are co-extracted. The Mo is stripped using acetohydroxamic acid (AHA) and fed onto an anion exchange column. Following a series of hydroxide, HCl, and oxalic acid wash steps, the Mo can be recovered in NaOH/NaCl. A manuscript describing this process chemistry in more detail was recently published.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2024

Hydroxide Exchange Membrane Carbon Capture (HEMCC) Using Nickel Hydroxide Batteries and Flow-through Membranes

Proposed is an electrochemical nickel hydroxide based hydroxide exchange membrane carbon capture (HEMCC) device for Direct Air Capture (DAC) of CO2. DAC has been identified as one of the key net negative carbon technologies to achieve net zero carbon emissions. Net negative carbon technologies are required to offset continued emissions from dilute CO2 sources such as agriculture and construction. The majority of current DAC technologies at scale (>1 KT∙yr-1) are adsorbent based technologies with significant energy cost. The traditional DAC energy cost is primarily driven by the temperature swing required to regenerate the sorbent and has been shown to be 1.8 MWh·ton-1 at the system level. Electrochemical pH swing devices are a growing research area for carbon capture devices with the goal of lowering the energy cost required for DAC. A pH gradient is built by generating OH- at the cathode and consuming OH- at the anode. An acid-base equilibrium with CO2 allows for the capture of CO2 at the cathode and release at the anode. This extends from other electrochemical CO2 capture devices based on pKa shifts of an electrochemically active organic species allowing for the capture and release of CO2. Electrochemical CO2 capture is considered promising based on potentially low energy costs to capture CO2 in comparison with current temperature swing adsorption technologies. This work explores Ni(OH)2 electrodes to produce the pH gradient for CO2 capture and release. At the cathode NiOOH is reduced to Ni(OH)2 while at the anode Ni(OH)2 is oxidized to NiOOH. The symmetrical electrodes allow for a low voltage requirement; the thermodynamic potential difference of standard electrochemical reactions is zero. Most of the voltage observed is to produce the pH gradient with the remainder driving the polarization of the electrodes. There is a resistance component as well, but this is small in comparison due to the low current densities used in the device, nominally 2 mA·cm-1. Two similar devices are presented, a traditional MEA (membrane electrode assembly) and a flow-through MEA. The traditional MEA separates the two Ni(OH)2 electrodes with an 80μm Piperion® membrane. While the flow-through membrane separates the electrodes with a three piece membrane consisting of two 80μm Piperion® membranes with a porous membrane between them. In the traditional MEA system air is passed over the cathode for capture, while the flow-through MEA the air is passed through the porous membrane isolated from the electrodes. The traditional MEA has been used to establish a baseline performance of the device and has been shown to capture CO2 at an energy cost of 1 MWh·ton-1 at the device level. An understanding has been built around the components of that energy cost including the relationship of flux to current density, effect of a regeneration process, transient battery behavior, and gas losses coinciding with changing the polarization of the batteries. The flow-through MEA looks to address of transient battery behavior and gas losses. It allows for denser, higher capacity electrodes, which can lean on traditional Ni-MH battery technology used in alkaline batteries used today. The higher capacities, limit the transient battery effect on flux in the device. Gas losses are addressed by having a continuous inlet air stream to the device and continuous outlet product.

Buchen, James↗

DOE OSTI2024

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage

Industrial process heating (IPH) accounts for approximately 70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately 1500 TWht (approximately 60%) of IPH demand is in the temperature range of 100-300. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. TES development has historically been dominated by technologies suitable for deployment with concentrating solar power (CSP). State-of-the-art thermal storage deployed commercially with power tower CSP plants uses a 60%/40% NaNO3/KNO3 molten salt and operates between temperatures of approximately 280 degrees Celsius and 570 degrees Celsius using a two-tank configuration. However, these nitrate salts are unsuitable for operation outside of this temperature range due to a high freezing point of approximately 220 degrees Celsius, and limits on high-temperature salt stability and corrosion resistance of containment alloys. Other materials being investigated for TES include those based on: (1) sensible energy storage (various molten salt compositions, inert solid particles, rocks or pebble beds, sulfur, water, concrete, graphite, etc.), (2) latent energy storage in materials that undergo solid-liquid phase change at relevant temperatures (organic materials for low-temperature applications, inorganic salts and/or metals for high-temperature applications), or (3) thermochemical energy storage (hydrides, hydroxides, carbonates, metal oxides, etc.). The application temperature and challenges pertaining to storage material and/or containment cost, energy density, long-term thermal and cyclic stability, and charge/discharge heat transfer effectiveness drive material selection for a given IPH or electricity generation application. Sulfur is a cheap commodity at $80/ton compared to $1100 - 1300/ton for conventional salts. When using a metric of storage cost per kWh, sulfur costs around 2-3 $/kWh. Previous sulfur TES development focused on high temperature (>600 degrees) concentrated solar power applications with sulfur encapsulated in pipes and flow of gaseous HTF (air) in the shell side. However, for lower-temperature IPH applications in the range of approximately 100-300 degrees Celsius Element 16 adopted a compact and scalable TES design with molten sulfur in the shell and HTF pipes submerged in the molten sulfur bath. The low-cost molten sulfur TES for dispatchable IPH has deployment potential for broad applications. The spatial and temporal evolution of the HTF and sulfur temperature is critical to the TES system performance, and thus detailed modeling can improve understanding of the performance and facilitate design improvements. Using high performance computing and computational fluid dynamics (CFD) a low-cost molten sulfur thermal energy storage (TES) system for industrial process heating (IPH) applications was developed. The unique challenges in CFD modeling of sulfur TES are the sharp property changes of sulfur relevant to the working temperatures. Above 159, liquid sulfur undergoes polymerization, and the viscosity of sulfur rapidly increases by several orders of magnitude between 159 degrees Celsius and 188 degrees Celsius, followed by a decrease in viscosity beyond 188 degrees Celsius due to thermal bound dissociation. In addition, various concentrations of H2S impurities can also modify sulfur viscosity. This numerical challenge is especially relevant to transient simulation of the sulfur TES charging and discharging processes as the extreme property variations limit the applicability of traditional heat transfer correlations. Transient CFD simulations including the temperature-dependent sulfur properties and geometric complexity of the TES design were used to predict the effect of natural convection during charging and discharging on the heat transfer process, sulfur temperature uniformity, charge/discharge rates, and performance of the storage devices. The CFD model was validated with experimental results for a full charge and discharge cycle. The work will show 3D and 2D simulation comparisons aimed to facilitate rapid design iterations and a machine learning based design optimization approach.

CFD↗

DOE OSTI2025

Excited State Covalency, Dynamics, and Photochemistry of Square Planar Ni-Thiolate Complexes Revealed by Ultrafast X-ray Absorption

Highly covalent Ni bis(dithiolene) and related complexes provide an ideal platform for investigating the effects of metal–ligand orbital hybridization on excited state character and dynamics. In particular, we focus on the ligand field excited states that dominate the photophysics of first-row transition metal complexes. Here, we investigate if they can be significantly delocalized off the metal center, possibly yielding photochemical reactivity more similar to charge transfer excited states than metal-centered ligand field excited states. Here [Ni(mpo) 2 ] (mpo = 2-mercaptopyridine- N-oxide) provides a representative example for the larger chemical class and is an active electro- and photo-catalyst for proton reduction. A detailed characterization of the excited state electronic structure, dynamics, and photochemistry of [Ni(mpo) 2 ] is presented based on ultrafast transient X- ray absorption spectroscopy at the Ni and S 1s core absorption K-edges. By comparing the ultrafast Ni K-edge absorption to ab initio calculations, we identify an excited state relaxation mechanism where an initial ligand-to-metal charge transfer excitation results in both excited state electron transfer (generating a catalytically relevant reduced photoproduct [Ni(mpo) 2 ] − ) and relaxation to a pseudo-tetrahedral triplet ligand field excited state. From the ultrafast S K-edge absorption, the ligand field excited state is found to be highly delocalized onto the thiolate ligands and a tetrahedral structural distortion is shown to substantially influence the degree of delocalization. The results identify a significant structural coordinate to target when aiming to control the excited state covalency in square planar complexes.

Lim, Hyeongtaek [SLAC National Accelerator Laborat↗

DOE OSTI2026

Exponential crystallization in corals

Corals form their reef-building aragonite (CaCO 3 ) skeletons via transient precursor phases yet understanding of the dynamics of these early-stage transformations remains incomplete. Using time-independent myriad mapping (MM) at 50 nm resolution, we map five mineral phases near the skeleton surface of Stylophora pistillata corals grown in varying seawater pH. All precursors, crystalline and amorphous, exhibit a consistent exponential decay from the growth front, with a shared decay length of 0.7 ± 0.1 μm, independent of time, phase, or pH. This spatial decay, paired with the constant growth rate of the skeleton, reveals a decay time of 5.1 ± 0.5 minutes. The dominant precursor is not amorphous but crystalline: calcium carbonate hemihydrate (CCHH, CaCO₃·½H₂O). These results suggest that exponential crystallization kinetics govern coral biomineralization and may be a widespread feature in biogenic, geologic, and synthetic systems—traceable long after initial mineral deposition.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Liquid-like spin dynamics in a hybrid Heisenberg-Ising antiferromagnet

Two-dimensional (2D) materials and their heterostructures enable unconventional electronic properties and functionalities not accessible in their bulk counterparts. This approach is now being extended to magnetic materials to engineer their spin structures and magnetic fields produced by them. However, spin dynamics of 2D magnetic heterostructures remain largely unexplored. Here, we demonstrate that heterointerfacing Heisenberg square-lattice antiferromagnet (AF) Sr 2 IrO 4 with its bilayer variant Ising AF Sr 3 Ir 2 O 7 in a superlattice leads to liquid-like spin dynamics in the former, characterized by slow recovery of the AF order after its transient suppression by an optical pump, and complete absence of spin waves except in an immediate vicinity of the ordering wavevector. Instead, the spin excitation spectra are dominated by isotropic continua, which in previous works have been interpreted as fractional spin excitations, or spinons, that extends to unprecedentedly low energies. Thus, our results provide a pathway to frustrated magnetism in square lattices by heterointerfacing two distinct types of AFs.

RIXS↗

DOE OSTI2025

The La Silla Schmidt Southern Survey

We present the La Silla Schmidt Southern Survey (LS4), a new wide-field, time-domain survey to be conducted with the 1 m ESO Schmidt telescope. The 268 megapixel LS4 camera mosaics 32 2k × 4k fully depleted CCDs, providing a ∼20 deg 2 field of view with 1″ pixel −1 resolution. The LS4 camera will have excellent performance at longer wavelengths: in a standard 45 s exposure the expected 5σ limiting magnitudes in g, i, z are ∼21.5, ∼20.9, and ∼20.3 mag (AB), respectively. The telescope design requires a novel filter holder that fixes different bandpasses over each quadrant of the detector. Two quadrants will have i band, while the other two will be g and z band with color information obtained by dithering targets across the different quadrants. The majority (90%) of the observing time will be used to conduct a public survey that monitors the extragalactic sky at both moderate (3 days) and high (1 day) cadence, as well as focused observations within the Galactic plane and bulge. Alerts from the public survey will be broadcast to the community via established alert brokers. LS4 will run concurrently with the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST). The combination of LS4+LSST will enable detailed holistic monitoring of many nearby transients: high-cadence LS4 observations will resolve the initial rise and peak of the light curve while less-frequent but deeper observations by LSST will characterize the years before and after explosion. Here, we summarize the primary science objectives of LS4 including microlensing events in the Galaxy, extragalactic transients powered by massive black holes or stellar explosions, the search for electromagnetic counterparts to multi-messenger events, and supernova cosmology.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Developing the S 32 ( p , d ) S * 31 ( p ) ( γ ) reaction to probe the P 30 ( p , γ ) S 31 reaction rate in classical novae

Background: The 30 P(p,γ) 31 S reaction rate is one of the largest remaining sources of uncertainty in the final abundances of nuclei created in a classical nova involving a ONe white dwarf. The reaction rate directly influences silicon isotopic ratios, which are used as identifiers of presolar grains with nova origins. Additionally, the uncertainty in the 30 P(p,γ) 31 S reaction rate has been found to limit the use of nova nuclear thermometers based on observations of elemental ratios in nova ejecta. Purpose: Reduce uncertainties in the nuclear data for proton-unbound states in 31 S, which act as resonances for the 30 P(p,γ) 31 S reaction at classical nova temperatures, and develop a technique for high efficiency, high-resolution reaction-decay coincidence measurements. Methods: Here, the 32 S(p,d) 31 S reaction was used to populate the states of interest in 31 S. The experiment was performed at the Texas A&M Cyclotron Institute using the LLNL Hyperion array for the detection of charged particles and γ rays. A downstream silicon telescope was used to select reaction deuterons, and a single upstream silicon detector was used to measure protons emitted in the decay of unbound 31 S levels. Results: Several states in 31 S above the proton separation energy were observed to have been populated. Decay protons from the resonant states in 31 S were identified as events in the upstream silicon detectors that came in coincidence with deuterons in the downstream telescope. Protons emitted from these states were measured and branching ratios extracted. Conclusions: While no new reaction rate is derived, spin-parity assignments for several higher-lying proton unbound states have been confirmed. Measured p 0 branching ratios for these levels have been compared to previous measurements with good agreement, and in some cases provided a reduction in uncertainty. The previously identified T = 3/2 state may have been incorrectly assigned a large p 0 branching ratio in a previous measurement. The technique of measuring reaction-decay coincidences with a particle-gamma setup appears promising.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Population of Merging Compact Binaries Inferred Using Gravitational Waves through GWTC-3

We report on the population properties of compact binary mergers inferred from gravitational-wave observations of these systems during the first three LIGO-Virgo observing runs. The Gravitational-Wave Transient Catalog 3 (GWTC-3) contains signals consistent with three classes of binary mergers: binary black hole, binary neutron star, and neutron star–black hole mergers. We infer the binary neutron star merger rate to be between 10 and 1700 Gpc − 3 yr − 1 and the neutron star–black hole merger rate to be between 7.8 and 140 Gpc − 3 yr − 1 , assuming a constant rate density in the comoving frame and taking the union of 90% credible intervals for methods used in this work. We infer the binary black hole merger rate, allowing for evolution with redshift, to be between 17.9 and 44 Gpc − 3 yr − 1 at a fiducial redshift ( z = 0.2 ). The rate of binary black hole mergers is observed to increase with redshift at a rate proportional to ( 1 + z ) κ with κ = 2. 9 − 1.8 + 1.7 for z ≲ 1 . Using both binary neutron star and neutron star–black hole binaries, we obtain a broad, relatively flat neutron star mass distribution extending from 1.2 − 0.2 + 0.1 to 2.0 − 0.3 + 0.3 M ⊙ . We confidently determine that the merger rate as a function of mass sharply declines after the expected maximum neutron star mass, but cannot yet confirm or rule out the existence of a lower mass gap between neutron stars and black holes. We also find the binary black hole mass distribution has localized over- and underdensities relative to a power-law distribution, with peaks emerging at chirp masses of 8.3 − 0.5 + 0.3 and 27.9 − 1.8 + 1.9 M ⊙ . While we continue to find that the mass distribution of a binary’s more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above approximately 60 M ⊙ , which would indicate the presence of a upper mass gap. Observed black hole spins are small, with half of spin magnitudes below χ i ≈ 0.25 . While the majority of spins are preferentially aligned with the orbital angular momentum, we infer evidence of antialigned spins among the binary population. We observe an increase in spin magnitude for systems with more unequal-mass ratio. We also observe evidence of misalignment of spins relative to the orbital angular momentum. Published by the American Physical Society 2023

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗