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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Water Ice Clouds in the Martian Atmosphere: A Comparison of Two Methods and Eras

Similar cloud features are seen in maps generated with each method with no obvious outliers. The temperature differencing method appears to possibly be somewhat more sensitive to weaker water ice signatures. We have also generated correlation plots comparing the two methods. At strong delta-T signals, the correlation between the two methods is quite good, and therefore extraction of opacities from earlier Viking data may be possible for these stronger detection levels. Weaker detections do not, however, show such a good correlation. We are currently analyzing why the correlation becomes poor at weak signal levels, though it may be due to the fact that the differencing method may be more sensitive to thin cloud hazes. Results of this ongoing analysis will be presented. A comparison of the Viking and Mars Global Surveyor (MGS) eras are also presented.

A S Hale↗

DOE OSTI2026

Identifying Anomalous DESI Galaxy Spectra with a Variational Autoencoder

The tens of millions of spectra being captured by the Dark Energy Spectroscopic Instrument (DESI) provide tremendous discovery potential. In this work we show how Machine Learning, in particular Variational Autoencoders (VAE), can detect anomalies in a sample of approximately 200,000 DESI spectra comprising galaxies, quasars and stars. We demonstrate that the VAE can compress the dimensionality of a spectrum by a factor of 100, while still retaining enough information to accurately reconstruct spectral features. We then detect anomalous spectra as those with high reconstruction error and those which are isolated in the VAE latent representation. The anomalies identified fall into two categories: spectra with artefacts and spectra with unique physical features. Awareness of the former can help to improve the DESI spectroscopic pipeline; whilst the latter can lead to the identification of new and unusual objects. To further curate the list of outliers, we use the Astronomaly package which employs Active Learning to provide personalised outlier recommendations for visual inspection. In this work we also explore the VAE latent space, finding that different object classes and subclasses are separated despite being unlabelled. We demonstrate the interpretability of this latent space by identifying tracks within it that correspond to various spectral characteristics. For example, we find tracks that correspond to increasing star formation and increase in broad emission lines along the Balmer series. In upcoming work we hope to apply the methods presented here to search for both systematics and astrophysically interesting objects in much larger datasets of DESI spectra.

Nicolaou, C. [University Coll. London] (ORCID:0000↗

NASA NTRS2018

Maximum Sample Temperature for Mars Sample Return: A Historical Perspective

Since the first Mars Sample Return (MSR) report published by the Jet Propulsion Laboratory (JPL) in 1974 [1], a series of panels, reports, and white papers have recognized the importance of sample temperature and offered an informed sample maximum temperature (henceforth SMT) limit for returning martian samples to Earth. The Mars Sample Handling and Requirements Panel (MSHARP, 1999) stated that "[t]he main issue in sample preservation is temperature" [2]. More recently, the Mars Exploration Program Analysis Group (MEPAG)'s "Science Priorities for Mars Sample Return" report (2008), declared that "[s]ignificant loss, particularly to biological studies, occurs if samples reach +50C for three hours", whereby "scientific objectives related to life goals could be seriously compromised" [3]. By contrast, the Mars 2020 mission has adopted a SMT of +60C as spelled out in Beaty et al., 2016 [4]. Samples will be collected and then deposited on the surface in sealed tubes for possible retrieval and return to Earth. Beaty et al. [4] calculates that the samples will experience maximum temperatures of ~+30 to +60C, depending on latitude. At present, there is no mission requirement for the measurement/data logging of sample temperature during this period. We will explore the history of martian SMTs, as they have been recorded since 1974 [1], effectively representing input across multiple generations of Mars scientists. Ten separate publications present SMTs for MSR samples [1-10]. One report [10] is for a mission concept specifically designed to exclude life detection investigations, and recommended an SMT of 50C. Another did not specify a temperature, recommending "Mars ambient temperature" [5]. Of the remaining eight, SMTs are given as: -30C [1], -20C [3], 60C [4], -73 to 41C depending on sample type [6], -40C [7], -43 to 13C depending on type [2,8], and -33C [9]. If we restrict the temperatures to samples highlighted in the Mars 2020 mission goals, i.e. organics-bearing and sedimentary rocks, then the average SMT is -28+/-39C (n=8). Applying a Dixon's Q Test at P=0.05 (two-tailed), the 60C SMT [4] fails with Q=0.602 versus Qcrit=0.526. Excluding the outlier produces an average SMT of -40+/-17C (n=7). Therefore, the average SMT expressed by the Mars science community over the past 44 years (two generations) is a sample temperature no greater than -40C. The difference in chemical reaction rates between this average SMT and Beaty et al [4] can be estimated using the Arrhenius equation. Assuming a generic chemical reaction with an activation energy of 50 kJ/mol and a pre-exponential factor invariant with temperature, this reaction will proceed 2300x faster at 60C than at -40C. To illustrate the effects of the increased reaction rate, consider 10 ppb of alanine in a Mars 2020 cache, and assume that it becomes unmeasurable if it degrades to 1 ppb, as per the Mars 2020 Organic Contamination Panel contamination limits [11]. If we illustrate the effect with an arbitrary degradation rate such that the alanine will become undetectable in ten years at -40C, then the same 10 ppb alanine degrades beyond detectability in only 38 days at 60C. Further research is required to quantify expected analyte losses in the cached samples due to thermal processing.

Fries, Marc↗

DOE OSTI2026

Probing the γ -Ray Emission Origin of Two Star-forming Galaxies NGC 2403 and NGC 3424 with the Fermi-LAT

Star-forming galaxies (SFGs) are a subclass of γ-ray emitters, and a correlation between their γ-ray luminosity (L γ ) and the total infrared (IR) luminosity (L IR ) has been established based on the Fermi Large Area Telescope (LAT) data. NGC 2403 and NGC 3424 have been reported as outliers in the L γ –L IR correlation with light curves showing significant variability, which contrasts with the temporally stable γ-ray emission in other SFGs, originating primarily from cosmic rays interacting with interstellar medium. In this study, we reanalyze the γ-ray emission in the directions of NGC 2403 and NGC 3424 using more than 16.5 yr Fermi-LAT data. NGC 3424 is found to be spatially coincident with the detected γ-ray source, while NGC 2403 is significantly offset from the nearest γ-ray source, suggesting an implausible association. We confirm the previously reported variability of both γ-ray sources and the significant deviation from the L γ –L IR correlation when assuming an association of both γ-ray sources with the two galaxies. Our findings lend further support to the interpretation that their γ-ray emission is driven primarily by alternative radiative processes—rather than by star formation activity—such as the ejecta of the Type IIP supernova SN 2004dj in NGC 2403 interacting with a surrounding high-density shell and an obscured active galactic nucleus in NGC 3424.

Liu, Linjie [Chinese Academy of Sciences (CAS), Ku↗

DOE OSTI2025

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↗

NASA NTRS2010

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.↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

Hydrometeor Identification from GMI Radiometer, Trained Using Polarimetric Radar

Daniel J. Cecil↗

NASA NTRS2018

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↗

DOE OSTI2020

Radar Derived Rainfall and Rain Gauge Measurements at SRS

Over the years rainfall data for the Savannah River Site has been obtained from ground level measurements made by rain gauges. These instruments have inherent errors or biases that can impact the measured rainfall totals but are assumed as ground truth for most climatological and weather applications. With the development of weather radar technologies, various methods to derive rainfall totals from radar reflectivity values have been developed and have continued to improve. The Z-R relationship, which uses an exponential relationship to estimate rainfall rate based on radar reflectivity values, provides estimates of rainfall amounts (Z-R Level III) for locations within the radar’s detection range. The recently developed Multi-Radar Multi-Sensor (MRMS) dataset combines reflectivity-based estimates using the Z-R relationship with a network of gauges and other rainfall estimates to produce a refined set of precipitation estimates for each grid point within its domain. Comparisons were done between gauge observations and radar estimates for various SRS locations to assess whether radar derived estimates are representative of rainfall measurements at the site. Results obtained show good agreement between radar derived amounts and ground measurements, with MRMS showing stronger correlations and lower spread than Z-R Level III estimates. Outliers and errors observed appear to be related to hydrometeor classification schemes. In general, MRMS proved to provide sufficiently representative estimates of daily rainfall amounts to replace existing SRS rain gauge measurements.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2019

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↗

NASA NTRS2022

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↗

DOE OSTI2022

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↗

DOE OSTI2024

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↗

NASA NTRS1996

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.↗

DOE OSTI2025

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↗

DOE OSTI2022

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↗

DOE OSTI2025

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↗