Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Outlier detection”

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.

144 records · Page 8

Model-free estimation of completeness, uncertainties, and outliers in atomistic machine learning using information theory

Abstract An accurate description of information is relevant for a range of problems in atomistic machine learning (ML), such as crafting training sets, performing uncertainty quantification (UQ), or extracting physical insights from large datasets. However, atomistic ML often relies on unsupervised learning or model predictions to analyze information contents from simulation or training data. Here, we introduce a theoretical framework that provides a rigorous, model-free tool to quantify information contents in atomistic simulations. We demonstrate that the information entropy of a distribution of atom-centered environments explains known heuristics in ML potential developments, from training set sizes to dataset optimality. Using this tool, we propose a model-free UQ method that reliably predicts epistemic uncertainty and detects out-of-distribution samples, including rare events in systems such as nucleation. This method provides a general tool for data-driven atomistic modeling and combines efforts in ML, simulations, and physical explainability.

36 MATERIALS SCIENCE↗

Material condition assessment with eddy current sensors

Eddy current sensors and sensor arrays are used for process quality and material condition assessment of conducting materials. In an embodiment, changes in spatially registered high resolution images taken before and after cold work processing reflect the quality of the process, such as intensity and coverage. These images also permit the suppression or removal of local outlier variations. Anisotropy in a material property, such as magnetic permeability or electrical conductivity, can be intentionally introduced and used to assess material condition resulting from an operation, such as a cold work or heat treatment. The anisotropy is determined by sensors that provide directional property measurements. The sensor directionality arises from constructs that use a linear conducting drive segment to impose the magnetic field in a test material. Maintaining the orientation of this drive segment, and associated sense elements, relative to a material edge provides enhanced sensitivity for crack detection at edges.

Goldfine, Neil J.↗

Empirical Hydrometeor Type Identification from GMI Brightness Temperature Measurements

Brightness temperatures (TB) from GMI and related sensors contain information about the types and amounts of precipitation particles in a column. Approaches like GPROF or machine learning are good for providing precipitation estimates and other quantified information, but it can be difficult to understand why / how a given set of TBs traces to a given retrieval solution. Outlier TBs that are not well-represented in a training sample can lead to dubious solutions. Our goal is to use the scattering signatures at multiple frequencies to determine whether hail is present, or otherwise graupel, or otherwise snow, or otherwise liquid rain without a substantial precipitation ice component. (A hierarchical approach is employed, so a prediction of hail can include the presence of all other categories, but a prediction of snow implies a lack of graupel or hail, and a prediction of liquid rain implies a lack of any detectable precipitation ice.)

Daniel J. Cecil↗

Hydrometeor Identification from GMI Radiometer, Trained Using Polarimetric Radar

Brightness temperatures (TB) from GMI and related sensors contain information about the types and amounts of precipitation particles in a column. Approaches like GPROF or machine learning are good for providing precipitation estimates and other quantified information, but it can be difficult to understand why / how a given set of TBs traces to a given retrieval solution. Outlier TBs that are not well-represented in a training sample can lead to dubious solutions. Our goal is to use the scattering signatures at multiple frequencies to determine whether hail is present, or otherwise graupel, or otherwise snow, or otherwise liquid rain without a substantial precipitation ice component. (A hierarchical approach is employed, so a prediction of hail can include the presence of all other categories, but a prediction of snow implies a lack of graupel or hail, and a prediction of liquid rain implies a lack of any detectable precipitation ice.)

Daniel J. Cecil↗

The NANOGrav 11-year Data Set: High-Precision Timing of 45 Millisecond Pulsars

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, and six high-timing-precision pulsars were observed weekly. All were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and (for binary pulsars) orbital parameters; time-variable dispersion delays; and parameters that quantify pulse-profile evolution with frequency. The timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of significant orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. A companion paper uses these data to constrain the strength of the gravitational-wave background

Arzoumanian, Zaven↗

GOES-16 and GOES-17 ABI INR Assessment

The first two satellites of the US Geostationary Operational Environmental Satellite R-Series (GOES-R) were launched on November 19, 2016 and March 1, 2018 respectively. GOES-16 officially became GOES East on December 18, 2017, and the designation of GOES-17 as GOES West occurred on February 12 2019. The Advanced Baseline Imager (ABI) is the primary instrument on GOES-16 and GOES-17 for imaging Earth’s surface and atmosphere to significantly improve the detection and observation of severe environmental phenomena. The Image Navigation and Registration (INR) Performance Assessment Tool Set (IPATS) was developed to assess INR performance of GOES-R series ABI images. In this paper, we first describe the assessment of IPATS algorithm accuracy. Next, we present the relationship between view zenith angle (VZA) and the quality of the IPATS measurements. Lastly, we present GOES-16 and GOES-17 navigation (NAV) assessments results from flight data spanning from the start of INR assessment to June 2019. The results show a) IPATS “stair step” measurement error is less or equal to 0.06 ABI pixel with IPATS baseline configuration; b) VZA is an effective filter to exclude outliers of the measurements; and c) ABI INR for both satellites has improved over time as post-launch tests (PLT) were performed and corrections applied. This paper also shows that the post-launch INR tuning of GOES-17 was much shorter than GOES-16.

GOES-17↗

In Situ Measurements of Surface Texture with Virtual Environments Support Science-Driven Human Surface Operations on the Moon and Beyond

Visualization tools enabling real-time scientific analysis are important for supporting future astronaut operations on the lunar surface. Such tools can be built into virtual environments to support scientific investigations, as well as situational awareness, real-time decision making, and efficient communication between astronauts and ground and support systems. Understanding how these tools can be optimized for science is essential for upcoming Artemis missions. In this contribution, we discuss how measurements of surface texture at multiple length scales can greatly enhance in situ science on/of the Moon, and eventually Mars, asteroids, and beyond. Roughness measurements at various wavelengths directly support objectives defined in the Artemis Science Plan, including (O1) “understanding planetary processes,” (O2) “understanding volatile cycles,” and (O3) “interpreting the impact history of the Earth-Moon system” . Key scientific analyses enabled by texture measurements at different length scales include: ● Sub-centimeter scales: Texture measurements can help constrain lava flow crystallinity, lava rheology, emplacement flow dynamics, and cooling histories (O1). Measurements of lacunarity (voids in fractal fill space) can shed light on eruptive volatile content, residence time of migrating volatiles, and near-surface volume available for micro-cold trapping of volatiles (O1, O2). ● Centimeter–meter scales: Texture measurements can be used for the differentiation of individual lava flows, the reconstruction of local stratigraphies and emplacement sequences, characterization of post-emplacement surface modification processes (O1, O3). Derived roughness (polarization) metrics can be used in the detection of water ice and characterization of ice properties (e.g., purity, grade, depth, abundance). ● Hectometer–Kilometer scales: Texture measurements can be used to differentiate major geologic surface units and surface structures (O1), constrain the presence of abundant ground ices (O2), and analyze surface modification and estimate surface age (O3). Real-time measurements of surface texture across these multiple length scales will enable efficient sample identification and scientific investigations by future astronauts. To support these investigations and the objective classification of surface texture, virtual environments employed by astronauts should be able to instantaneously convert raw data into processed data (e.g., digital terrain and elevation models) and derived metrics (e.g., RMS, std, Hurst, CPR) and perform statistical analyses (e.g., PCA, outliers, correlation matrices). Such tools are being developed and tested by the Resource Exploration and Science of our Cosmic Environment (RESOURCE) team, a node of NASA’s Solar System Exploration Research Virtual Institute (SSERVI), and are an excellent example of the powerful synergies of human and robotic ground assets critical in the return of humans to the Moon.

Ariel N. Deutsch↗

Maximum Likelihood Spectrum Decomposition for Isotope Identification and Quantification

A spectral decomposition method has been implemented to identify and quantify isotopic source terms in high-resolution gamma-ray spectroscopy in static geometry and shielding scenarios. Monte Carlo simulations were used to build the response matrix of a shielded high-purity germanium detector monitoring an effluent stream with a Marinelli configuration. The decomposition technique was applied to a series of calibration spectra taken with the detector using a multi-nuclide standard. These results are compared with decay-corrected values from the calibration certificate. For most nuclei in the standard ( 241 Am, 109 Cd, 137 Cs, and 60 Co), the deviations from the certificate values were generally no more than 6% with a few outliers as high as 10%. Furthermore, for 57 Co, the radionuclide with the lowest activity, the deviations from the standard reached as high as 25%, driven by the meager statistics in the calibration spectra. In addition, a complete treatment of error propagation for the technique is presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

e RPCA : Robust Principal Component Analysis for Exponential Family Distributions

Abstract Robust principal component analysis (RPCA) is a widely used method for recovering low‐rank structure from data matrices corrupted by significant and sparse outliers. These corruptions may arise from occlusions, malicious tampering, or other causes for anomalies, and the joint identification of such corruptions with low‐rank background is critical for process monitoring and diagnosis. However, existing RPCA methods and their extensions largely do not account for the underlying probabilistic distribution for the data matrices, which in many applications are known and can be highly non‐Gaussian. We thus propose a new method called RPCA for exponential family distributions (), which can perform the desired decomposition into low‐rank and sparse matrices when such a distribution falls within the exponential family. We present a novel alternating direction method of multiplier optimization algorithm for efficient decomposition, under either its natural or canonical parametrization. The effectiveness of is then demonstrated in two applications: the first for steel sheet defect detection and the second for crime activity monitoring in the Atlanta metropolitan area.

Zheng, Xiaojun↗

SeaWiFS Technical Report Series: SeaWiFS Calibration and Validation Quality Control Procedures - Volume 38

This document provides five brief reports that address several quality control procedures under the auspices of the Calibration and Validation Element (CVE) within the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Project. Chapter 1 describes analyses of the 32 sensor engineering telemetry streams. Anomalies in any of the values may impact sensor performance in direct or indirect ways. The analyses are primarily examinations of parameter time series combined with statistical methods such as auto- and cross-correlation functions. Chapter 2 describes how the various onboard (solar and lunar) and vicarious (in situ) calibration data will be analyzed to quantify sensor degradation, if present. The analyses also include methods for detecting the influence of charged particles on sensor performance such as might be expected in the South Atlantic Anomaly (SAA). Chapter 3 discusses the quality control of the ancillary environmental data that are routinely received from other agencies or projects which are used in the atmospheric correction algorithm (total ozone, surface wind velocity, and surface pressure; surface relative humidity is also obtained, but is not used in the initial operational algorithm). Chapter 4 explains the procedures for screening level-, level-2, and level-3 products. These quality control operations incorporate both automated and interactive procedures which check for file format errors (all levels), navigation offsets (level-1), mask and flag performance (level-2), and product anomalies (all levels). Finally, Chapter 5 discusses the match-up data set development for comparing SeaWiFS level-2 derived products with in situ observations, as well as the subsequent outlier analyses that will be used for evaluating error sources.

Hooker, Stanford B.↗

Comparison of removal and spatial mark‐resight models for estimating wild pig density

Density estimation is critical to effectively manage invasive species and elucidate areas of highest concern. For wild pigs (Sus scrofa), the ability to estimate density is complicated because of their variable home range sizes and social structure. Common methods for estimating density (e.g., mark-recapture) may be unsuitable in management applications because additional data needs to be collected before and after management. Removal models offer a suitable alternative to estimate density changes following management and can be applied broadly across areas where management of wild pigs is ongoing. We collected wild pig removal and camera trap data from 25 private properties ranging in size from approximately 0.5 km 2 to 95 km 2 across 3 ecoregions in South Carolina, USA, from 2020–2023. We compared factors affecting consistency and precision of property-level density estimates between removal and spatial mark-resight (SMR) models. In general, excluding 1 large outlier, density estimates from removal models were between 0.60 and 15.85 wild pigs/km 2 (median = 5.34) with a median coefficient of variation (CV) of 0.76 and 95% confidence intervals for the CV between 0.70 and 0.94. Similarly, excluding 1 large outlier, density estimates from SMR were between 0.22 and 30.97 wild pigs/km 2 (median = 5.48) with a median CV of 0.39 and 95% confidence intervals for the CV between 0.38 and 1.20. We found the precision of removal models was affected primarily by the number of wild pigs dispatched in the removal period (3 months) and the ecoregion in which they were removed. None of the covariates, including the number of recaptures (a corresponding measure of sample size), influenced precision of the SMR models, although recaptures did influence the density estimates. At the individual property level, density estimates from our 2 estimators were dissimilar from each other in approximately 80% of instances, although none of the covariates we examined influenced dissimilarity. Our results provide unique insight into how sample size affects density estimates using 2 common methods and into novel SMR models that incorporate both marked and unmarked detections. In addition, the density estimates in this study can be used as a reference for wild pig densities in common land cover types throughout the southeastern United States.

60 APPLIED LIFE SCIENCES↗

Implications of the correlation between bulge-to-total baryonic mass ratio and the number of satellites for SAGA galaxies

We searched for correlations between the number of satellites and fundamental galactic properties for the Milky Way-like host galaxies in order to better understand their diverse satellite populations. We specifically aim to understand why galaxies that are very similar in stellar mass content, star formation rate, and local environment have very different numbers of satellites. Deep and extensive spectroscopic observations are needed to characterize the complete satellite luminosity function beyond the Local Group. One such endeavor is an ongoing Satellites of Galactic Analogs (SAGA) spectroscopic survey that has completed spectroscopic observations of 36 Milky Way-like galaxies within their virial radii down to the luminosity of Leo I dwarf galaxy. We correlated the number of satellites of SAGA galaxies with several fundamental properties of their hosts – including total specific angular momentum, which is considered to be well preserved throughout galaxy lifetime – in an attempt to identify the main driver of their diverse satellite populations. We aim to reveal some intrinsic galactic property decisive in making more or less satellites irrespective of baryonic mass or the environment in which galaxies reside. We modeled Spitzer Heritage Archive images of SAGA host galaxies at 3.6 and 4.5 microns with GALFIT code to obtain their stellar masses. We also searched the Extragalactic Database for information on their gas content and rotation velocities. Empirical correlations, like the baryonic Tully–Fisher relation and the stellar mass–size relation were used to exclude outliers. All the available galactic properties from the literature along with measured stellar masses were correlated with the number of satellites and no significant correlation was found. However, when we considered the “expected” number of satellites based on the correlation between the baryonic bulge-to-total ratio and the number of satellites confirmed for several nearby galaxies then strong correlations emerge between this number and (1) the mass of the bulge, and (2) the total specific angular momentum. The first correlation is positive, implying that galaxies with more massive bulges have more satellites, as already confirmed. Furthermore, the second correlation with the angular momentum is negative, meaning that, the smaller the angular momentum, the greater the number of expected satellites. This would imply that either satellites cannot form if galaxy angular momentum is too high, or that satellites form inside-out, so that angular momentum is being transferred to the outer parts of the galaxies. However, deeper spectroscopic observations are needed to confirm these findings, because they rely on the expected rather than detected number of satellites. There was a luminosity limit to the SAGA survey equivalent to the luminosity of Leo I dwarf satellite of the Milky Way galaxy (the SAGA limit). In particular, correlations found in this work are very susceptible to the total number of satellites of the NGC 4158 galaxy. This galaxy is predicted to have many more satellites than detected up to the SAGA limit.

79 ASTRONOMY AND ASTROPHYSICS↗

On Finding Black Holes in Photometric Microlensing Surveys

There are expected to be millions of isolated black holes in the galaxy resulting from the deaths of massive stars. Measuring the abundance and properties of this remnant population would shed light on the end stages of stellar evolution and the evolution paths of black hole systems. Detecting isolated black holes is currently only possible via gravitational microlensing, which has so far yielded one definitive detection. The difficulty in finding microlensing black holes lies in having to choose a small subset of events, based on characteristics of their light curves, to allocate expensive and scarce follow-up resources to confirm the identity of the lens. Current methods either rely on simple cuts in parameter space without using the full distribution information or are only effective on small subsets of events. In this paper, we present a new lens classification method. The classifier takes in posterior constraints on light-curve parameters and combines them with a Galactic simulation to estimate the lens class probability. This method is flexible and can be used with any set of microlensing light-curve parameters, making it applicable to large samples of events. We make this classification framework available via the popclass Python package. We apply the classifier to ~10,000 microlensing events from the Optical Gravitational Lensing Experiment survey and find 23 high-probability black hole candidates. Our classifier also suggests that the only known isolated black hole is an observational outlier, according to current Galactic models, and the allocation of astrometric follow-up on this event was a high-risk strategy.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing dark matter with strong gravitational lensing through an effective density slope

ABSTRACT Many dark matter (DM) models that are consistent with current cosmological data show differences in the predicted (sub)halo mass function, especially at sub-galactic scales, where observations are challenging due to the inefficiency of star formation. Strong gravitational lensing has been shown to be a useful tool for detecting dark low-mass (sub)haloes through perturbations in lensing arcs, therefore allowing the testing of different DM scenarios. However, measuring the total mass of a perturber from strong lensing data is challenging. Overestimating or underestimating perturber masses can lead to incorrect inferences about the nature of DM. In this paper, we argue that inferring an effective slope of the DM density profile, which is the power-law slope of perturbers at intermediate radii, where we expect the perturber to have the largest observable effect, is a promising way to circumvent these challenges. Using N-body simulations, we show that (sub)halo populations under different DM scenarios differ in their effective density slope distributions. Using realistic mocks of Hubble Space Telescope observations of strong lensing images, we show that the effective density slope of perturbers can be robustly measured with high enough accuracy to discern between different models. We also present our measurement of the effective density slope $\gamma =1.96\substack{+0.12 \\ -0.12}$ for the perturber in JVAS B1938+666, which is a 2σ outlier of the cold DM scenario. More measurements of this kind are needed to draw robust conclusions about the nature of DM.

79 ASTRONOMY AND ASTROPHYSICS↗

Meta-Study of Particulate Detection Losses on Radioactive Air Sample Filters

Several mathematical relationships between air sample filter mass loading and the correlated analytical self-absorption factor were developed using data from other published research in this meta-study. Gross-alpha and -beta applications are addressed for this research. As filter media becomes loaded with particulate matter, there is potential for measurement losses due to self-absorption by mass loading. Components contributing to absorption include particulate dust, radioactive particulates, and filter material. Standards indicate a correction factor should be used when the penetration of radioactive material into the collection media or self-absorption of radiation by the material collected would reduce the detection rate by more than 5%. Previously, losses due to self-absorption have been reported up to 100% over a range up to ~10 mg∙cm -2 mass loading. These absorption losses then can be used to determine a correction factor for sample results. For low mass loadings (e.g., ≤0.1 mg∙cm -2 ) corrections factors in the 0.85 - 1 range have been recommended and used, while at higher mass loadings nearer to 10 mg∙cm -2 correction factors closer to 0 (representing near 100% losses) are used. Based on data from published studies, the different methods for relating percent loss due to self-absorption to mass loading include linear, exponential, quadratic, and trinomial derived functions. Where applicable, both forced zero and non-forced zero results were evaluated. From the derived functions evaluated, the trinomial function provided the best fit. Once the sample filter mass loading is known, the trinomial function can be applied to estimate losses and the corresponding self-absorption factor. When applied to routine operating conditions for radiological facility stacks monitored at the Pacific Northwest National Laboratory for an average sample filter mass loading of 0.09 ± 0.12 (2σ) mg∙cm -2 (excluding negative values and outliers) and a range from 0 - 0.24 mg∙cm -2 , the estimated trinomial function nominal self-absorption losses are less than 5% at 0.09 mg∙cm -2 and less than 10% at 0.24 mg∙cm -2 . The trinomial function is one method that may be used to adjust the activity results of an air sample when the sample-specific mass loading is determined. The application of no correction factor when the ANSI/HPS N13.1-2021 guidance of a 5% threshold for loss is not reached with typical stack sample mass loadings may be reasonable in high-efficiency particulate air filtered systems. For simplicity, it would be conservative in assigning the self-absorption correction factor at the 5% threshold (i.e., 0.95) for general uses but in cases of heavy mass loading to calculate the factor.

air sampling↗

HD 222925: A New Opportunity to Explore the Astrophysical and Nuclear Conditions of r-process Sites

Abstract With the most trans-iron elements detected of any star outside the solar system, HD 222925 represents the most complete chemical inventory among metal-poor stars enhanced with elements made by the rapid neutron capture (“ r ”) process. As such, HD 222925 may be a new “template” for the observational r -process, where before the (much higher-metallicity) solar r -process residuals were used. In this work, we test under which conditions a single site accounts for the entire elemental r -process abundance pattern of HD 222925. We found that several of our tests—with the single exception of the black hole–neutron star merger case—challenge the single-site assumption by producing an ejecta distribution that is highly constrained, in disagreement with simulation predictions. However, we found that ejecta distributions that are more in line with simulations can be obtained under the condition that the nuclear data near the second r -process peak are changed. Therefore, for HD 222925 to be a canonical r -process template likely as a product of a single astrophysical source, the nuclear data need to be reevaluated. The new elemental abundance pattern of HD 222925—including the abundances obtained from space-based, ultraviolet (UV) data—call for a deeper understanding of both astrophysical r -process sites and nuclear data. Similar UV observations of additional r -process–enhanced stars will be required to determine whether the elemental abundance pattern of HD 222925 is indeed a canonical template (or an outlier) for the r -process at low metallicity.

79 ASTRONOMY AND ASTROPHYSICS↗

A preliminary investigation into the feasibility of laser ablation U–Pb isotope ratio measurement via all‐faraday cup detection with 10 11 ‐ and 10 13 ‐Ω amplifiers on the Neoma multicollector‐inductively coupled plasma‐mass spectrometer

Rationale Signal detection for uranium–lead (U–Pb) dating of zircon is typically performed via ion counters. Here, we develop a preliminary understanding of the strengths and limitations of faraday‐cup‐based detection. Methods A suite of zircon reference materials and the NIST‐610 glass were sampled using laser ablation followed by U–Pb isotope ratio measurement on a Neoma multicollector‐inductively coupled plasma‐mass spectrometer. Results We were able to produce geologically accurate 207 Pb/ 206 Pb, 206 Pb/ 238 U, and 207 Pb/ 235 U ratios for the NIST‐610 glass and the zircon standards, with ages ranging from ~2.5 Ga to ~337 Ma (TanBrown A, Oracle, 91550, Mud Tank, Temora, and Plešovice). Two of the younger zircon standards examined (94‐35, ~55.6 Ma, and Fish Canyon, 28.6 Ma) yielded accurate 206 Pb/ 238 U but not 207 Pb/ 235 U or 207 Pb/ 206 Pb ratios, whereas the youngest zircon standard (Penglai, ~4.4 Ma) failed for all three ratios of interest. The accuracy and precision of the all‐faraday method are directly tied to signal intensity, with reliable data capable of being produced even when both isotopes in a ratio have signals below ~0.001 V (equivalent to ~62 500 cps on an ion counter). Conclusion The all‐faraday cup multicollection method provides sufficient sensitivity to obtain geologically meaningful U–Pb data, with possible advantages being that laser pit depth‐dependent changes in the observed interelemental fractionation behavior may be easier to correct using a static collector configuration compared to when the ion beam is swept across a single detector while also removing the need for an interdetector‐type calibration. Further work is needed to refine the all‐faraday cup method (e.g., application of background subtraction and common Pb corrections, outlier removal, and interelement as well as down‐hole fractionation corrections), but our initial results demonstrate that the faraday detector method has sufficient sensitivity to warrant further study.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Population Properties of Compact Objects from the Second LIGO-Virgo Gravitational-Wave Transient Catalog

We report on the population of 47 compact binary mergers detected with a false-alarm rate of <1yr^(−1) in the second LIGO–Virgo Gravitational-Wave Transient Catalog. We observe several characteristics of the merging binary black hole (BBH) population not discernible until now. First, the primary mass spectrum contains structure beyond a power law with a sharp high-mass cutoff; it is more consistent with a broken power law with a break at 39.7^(+20.3)(-9.1)M or a power law with a Gaussian feature peaking at 33.1^(+4.0)(-5.6)M (90% credible interval). While the primary mass distribution must extend to ~65M or beyond, only 2.9^(+3.5)(-1.7)M% of systems have primary masses greater than 45M. Second, we find that a fraction of BBH systems have component spins misaligned with the orbital angular momentum, giving rise to precession of the orbital plane. Moreover, 12% to 44% of BBH systems have spins tilted by more than 90°, giving rise to a negative effective inspiral spin parameter, χeff. Under the assumption that such systems can only be formed by dynamical interactions, we infer that between 25% and 93% of BBHs with nonvanishing |χeff| > 0.01 are dynamically assembled. Third, we estimate merger rates, finding RBBH = 23.9^(+14.3)(-8.6) Gpc^(-3) yr^(-1) for BBHs and RBNS = 320^(+490)(-240) Gpc^(-3) yr^(-1) for binary neutron stars. We find that the BBH rate likely increases with redshift (85% credibility) but not faster than the star formation rate (86% credibility). Additionally, we examine recent exceptional events in the context of our population models, finding that the asymmetric masses of GW190412 and the high component masses of GW190521 are consistent with our models, but the low secondary mass of GW190814 makes it an outlier.

R Abbott↗