skip to main content


Title: A Reanalysis of Public Galactic Bulge Gravitational Microlensing Events from OGLE-III and -IV
Abstract Modern surveys of gravitational microlensing events have progressed to detecting thousands per year, and surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is the microlensing Einstein radius crossing time ( t E ) distribution. However, systematics in individual light curves as well as oversimplistic modeling can lead to biased results. To address this, we developed a model to simultaneously handle the microlensing parallax due to Earth's motion, systematic instrumental effects, and unlensed stellar variability with a Gaussian process model. We used light curves for nearly 10,000 OGLE-III and -IV Milky Way bulge microlensing events and fit each with our model. We also developed a forward model approach to infer the t E distribution by forward modeling from the data rather than using point estimates from individual events. We find that modeling the variability in the baseline removes a source of significant bias in individual events, and the previous analyses overestimated the number of t E > 100 day events due to their oversimplistic model ignoring parallax effects. We use our fits to identify the hundreds filling a regime in the microlensing parameter space that are 50% pure of black holes. Finally, we have released the largest-ever catalog of Markov Chain Monte Carlo parameter estimates for microlensing events.  more » « less
Award ID(s):
1909641
NSF-PAR ID:
10357843
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
The Astrophysical Journal Supplement Series
Volume:
260
Issue:
1
ISSN:
0067-0049
Page Range / eLocation ID:
2
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We present a phenomenological forward Monte Carlo model for forecasting the population of active galactic nuclei (AGNs) in dwarf galaxies observable via their optical variability. Our model accounts for expected changes in the spectral energy distribution of AGNs in the intermediate-mass black hole (IMBH) mass range and uses observational constraints on optical variability as a function of black hole (BH) mass to generate mock light curves. Adopting several different models for the BH occupation function, including one for off-nuclear IMBHs, we quantify differences in the predicted local AGN mass and luminosity functions in dwarf galaxies. As a result, we are able to model the fraction of variable AGNs as a function of important galaxy host properties, such as host galaxy stellar mass, in the presence of selection effects. We find that our adopted occupation fractions for the ‘heavy’ and ‘light’ initial BH seeding scenarios can be distinguished with variability at the 2–3σ level for galaxy host stellar masses below ∼108M⊙ with data from the upcoming Vera C. Rubin Observatory. We also demonstrate the prevalence of a selection bias whereby recovered IMBH masses fall, on average, above the predicted value from the local host galaxy–BH mass scaling relation with the strength of this bias dependent on the survey sensitivity. Our methodology can be used more broadly to calibrate AGN demographic studies in synoptic surveys. Finally, we show that a targeted ∼ hourly cadence program over a few nights with the Rubin Observatory can provide strong constraints on IMBH masses given their expected rapid variability time-scales.

     
    more » « less
  2. Abstract

    Uncertainty in the initial–final mass relation (IFMR) has long been a problem in understanding the final stages of massive star evolution. One of the major challenges of constraining the IFMR is the difficulty of measuring the mass of nonluminous remnant objects (i.e., neutron stars and black holes). Gravitational-wave detectors have opened the possibility of finding large numbers of compact objects in other galaxies, but all in merging binary systems. Gravitational lensing experiments using astrometry and photometry are capable of finding compact objects, both isolated and in binaries, in the Milky Way. In this work we improve the Population Synthesis for Compact object Lensing Events (PopSyCLE)microlensing simulation code in order to explore the possibility of constraining the IFMR using the Milky Way microlensing population. We predict that the Roman Space Telescope’s microlensing survey will likely be able to distinguish different IFMRs based on the differences at the long end of the Einstein crossing time distribution and the small end of the microlensing parallax distribution, assuming the small (πE≲ 0.02) microlensing parallaxes characteristic of black hole lenses are able to be measured accurately. We emphasize that future microlensing surveys need to be capable of characterizing events with small microlensing parallaxes in order to place the most meaningful constraints on the IFMR.

     
    more » « less
  3. ABSTRACT We present a sample of 706, z < 1.5 active galactic nuclei (AGNs) selected from optical photometric variability in three of the Dark Energy Survey (DES) deep fields (E2, C3, and X3) over an area of 4.64 deg2. We construct light curves using difference imaging aperture photometry for resolved sources and non-difference imaging PSF photometry for unresolved sources, respectively, and characterize the variability significance. Our DES light curves have a mean cadence of 7 d, a 6-yr baseline, and a single-epoch imaging depth of up to g ∼ 24.5. Using spectral energy distribution (SED) fitting, we find 26 out of total 706 variable galaxies are consistent with dwarf galaxies with a reliable stellar mass estimate ($M_{\ast }\lt 10^{9.5}\, {\rm M}_\odot$; median photometric redshift of 0.9). We were able to constrain rapid characteristic variability time-scales (∼ weeks) using the DES light curves in 15 dwarf AGN candidates (a subset of our variable AGN candidates) at a median photometric redshift of 0.4. This rapid variability is consistent with their low black hole (BH) masses. We confirm the low-mass AGN nature of one source with a high S/N optical spectrum. We publish our catalogue, optical light curves, and supplementary data, such as X-ray properties and optical spectra, when available. We measure a variable AGN fraction versus stellar mass and compare to results from a forward model. This work demonstrates the feasibility of optical variability to identify AGNs with lower BH masses in deep fields, which may be more ‘pristine’ analogues of supermassive BH seeds. 
    more » « less
  4. Aims. We investigate the data collected by the high-cadence microlensing surveys during the 2022 season in search of planetary signals appearing in the light curves of microlensing events. From this search, we find that the lensing event MOA-2022-BLG-249 exhibits a brief positive anomaly that lasted for about one day, with a maximum deviation of ~0.2 mag from a single-source, single-lens model. Methods. We analyzed the light curve under the two interpretations of the anomaly: one originated by a low-mass companion to the lens (planetary model) and the other originated by a faint companion to the source (binary-source model). Results. We find that the anomaly is better explained by the planetary model than the binary-source model. We identified two solutions rooted in the inner-outer degeneracy and for both of them, the estimated planet-to-host mass ratio, q ~ 8 × 10 −5 , is very small. With the constraints provided by the microlens parallax and the lower limit on the Einstein radius, as well as the blend-flux constraint, we find that the lens is a planetary system, in which a super-Earth planet, with a mass of (4.83 ± 1.44) Μ ⊕ , orbits a low-mass host star, with a mass of (0.18 ± 0.05) M ⊙ , lying in the Galactic disk at a distance of (2.00 ± 0.42) kpc. The planet detection demonstrates the elevated microlensing sensitivity of the current high-cadence lensing surveys to low-mass planets. 
    more » « less
  5. Abstract

    Stellar variability is driven by a multitude of internal physical processes that depend on fundamental stellar properties. These properties are our bridge to reconciling stellar observations with stellar physics and to understand the distribution of stellar populations within the context of galaxy formation. Numerous ongoing and upcoming missions are charting brightness fluctuations of stars over time, which encode information about physical processes such as the rotation period, evolutionary state (such as effective temperature and surface gravity), and mass (via asteroseismic parameters). Here, we explore how well we can predict these stellar properties, across different evolutionary states, using only photometric time-series data. To do this, we implement a convolutional neural network, and with data-driven modeling we predict stellar properties from light curves of various baselines and cadences. Based on a single quarter of Kepler data, we recover the stellar properties, including the surface gravity for red giant stars (with an uncertainty of ≲0.06 dex) and rotation period for main-sequence stars (with an uncertainty of ≲5.2 days, and unbiased from ≈5 to 40 days). Shortening the Kepler data to a 27 days Transiting Exoplanet Survey Satellite–like baseline, we recover the stellar properties with a small decrease in precision, ∼0.07 for loggand ∼5.5 days forProt, unbiased from ≈5 to 35 days. Our flexible data-driven approach leverages the full information content of the data, requires minimal or no feature engineering, and can be generalized to other surveys and data sets. This has the potential to provide stellar property estimates for many millions of stars in current and future surveys.

     
    more » « less