skip to main content


Title: An efficient method for simulating light curves of cosmological microlensing and caustic crossing events
ABSTRACT

A new window to observing individual stars and other small sources at cosmological distances was opened recently, with the detection of several caustic-crossing events in galaxy cluster fields. Many more such events are expected soon from dedicated campaigns with the Hubble Space Telescope and the James Webb Space Telescope. These events can not only teach us about the lensed sources themselves, such as individual high-redshift stars, star clusters, or accretion discs, but through their light curves they also hold information about the point-mass function of the lens, and thus, potentially, the composition of dark matter. We present here a simple method for simulating light curves of such events, i.e. the change in apparent magnitude of the source as it sweeps over the net of caustics generated by microlenses embedded around the critical region of the lens. The method is recursive and so any reasonably sized small source can be accommodated, down to sub-solar scales, in principle. We compare the method, which we dub Adaptive Boundary Method, with other common methods such as simple inverse ray shooting, and demonstrate that it is significantly more efficient and accurate in the small-source and high-magnification regime of interest. A python version of the code is made publicly available in an open-source fashion for simulating future events.

 
more » « less
NSF-PAR ID:
10368137
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
514
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 2545-2560
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present the analysis of five black hole candidates identified from gravitational microlensing surveys. Hubble Space Telescope astrometric data and densely sampled light curves from ground-based microlensing surveys are fit with a single-source, single-lens microlensing model in order to measure the mass and luminosity of each lens and determine if it is a black hole. One of the five targets (OGLE-2011-BLG-0462/MOA-2011-BLG-191 or OB110462 for short) shows a significant >1 mas coherent astrometric shift, little to no lens flux, and has an inferred lens mass of 1.6–4.4M. This makes OB110462 the first definitive discovery of a compact object through astrometric microlensing and it is most likely either a neutron star or a low-mass black hole. This compact-object lens is relatively nearby (0.70–1.92 kpc) and has a slow transverse motion of <30 km s−1. OB110462 shows significant tension between models well fit to photometry versus astrometry, making it currently difficult to distinguish between a neutron star and a black hole. Additional observations and modeling with more complex system geometries, such as binary sources, are needed to resolve the puzzling nature of this object. For the remaining four candidates, the lens masses are <2M, and they are unlikely to be black holes; two of the four are likely white dwarfs or neutron stars. We compare the full sample of five candidates to theoretical expectations on the number of black holes in the Milky Way (∼108) and find reasonable agreement given the small sample size.

     
    more » « less
  2. Abstract

    The Transiting Exoplanet Survey Satellite (TESS) has an exceptionally large plate scale of 21″ px−1, causing most TESS light curves to record the blended light of multiple stars. This creates a danger of misattributing variability observed by TESS to the wrong source, which would invalidate any analysis. We developed a method that can localize the origin of variability on the sky to better than one fifth of a pixel. Given measured frequencies of variability (e.g., from periodogram analysis), we show that the best-fit sinusoid amplitudes to raw light curves extracted from each pixel are distributed in the same way as light from the variable source. The primary assumption of this method is that other nearby stars are not variable at the same frequencies. Essentially, we are using the high frequency resolution of TESS to overcome limitations from its low spatial resolution. We have implemented our method in an open-source Python package,TESS_localize(github.com/Higgins00/TESS-Localize), that determines the location of a variable source on the sky and the most likely Gaia source given TESS pixel data and a set of observed frequencies of variability. Our method utilizes models of the TESS pixel response function, and we characterize systematics in the residuals of fitting these models to data. We find that even stars more than three pixels outside a photometric aperture can produce significant contaminant signals in the extracted light curves. Given the ubiquity of source blending in TESS light curves, verifying the source of observed variability should be a standard step in TESS analyses.

     
    more » « less
  3. Abstract

    Given its large plate scale of 21″ pixel−1, analyses of data from the Transiting Exoplanet Survey Satellite (TESS) space telescope must be wary of source confusion from blended light curves, which creates the potential to attribute observed photometric variability to the wrong astrophysical source. We explore the impact of light curve contamination on the detection of fast yellow pulsating supergiant (FYPS) stars as a case study to demonstrate the importance of confirming the source of detected signals in the TESS pixel data. While some of the FYPS signals have already been attributed to contamination from nearby eclipsing binaries, others are suggested to be intrinsic to the supergiant stars. In this work, we carry out a detailed analysis of the TESS pixel data to fit the source locations of the dominant signals reported for 17 FYPS stars with the Python packageTESS_localize. We are able to reproduce the detections of these signals for 14 of these sources, obtaining consistent source locations for four. Three of these originate from contaminants, while the signal reported for BZ Tuc is likely a spurious frequency introduced to the light curve of this 127 day Cepheid by the data processing pipeline. Other signals are not significant enough to be localized with our methods, or have long periods that are difficult to analyze given other TESS systematics. Since no localizable signals hold up as intrinsic pulsation frequencies of the supergiant targets, we argue that unambiguous detection of pulsational variability should be obtained before FYPS are considered a new class of pulsator.

     
    more » « less
  4. Abstract

    We measure the Einstein radius of the single-lens microlensing event KMT-2022-BLG-2397 to beθE= 24.8 ± 3.6μas, placing it at the upper shore of the Einstein Desert, 9 ≲θE/μas ≲ 25, between free-floating planets (FFPs) and bulge brown dwarfs (BDs). In contrast to the six BD (25 ≲θE≲ 50) events presented by Gould et al. (2022), which all had giant-star source stars, KMT-2022-BLG-2397 has a dwarf-star source, with angular radiusθast∼ 0.9μas. This prompts us to study the relative utility of dwarf and giant sources for characterizing FFPs and BDs from finite-source point-lens (FSPL) microlensing events. We find “dwarfs” (including main-sequence stars and subgiants) are likely to yield twice as manyθEmeasurements for BDs and a comparable (but more difficult to quantify) improvement for FFPs. We show that neither current nor planned experiments will yield complete mass measurements of isolated bulge BDs, nor will any other planned experiment yield as manyθEmeasurements for these objects as the Korea Microlensing Telescope (KMT). Thus, the currently anticipated 10 yr KMT survey will remain the best way to study bulge BDs for several decades to come.

     
    more » « less
  5. ABSTRACT

    As we observe a rapidly growing number of astrophysical transients, we learn more about the diverse host galaxy environments in which they occur. Host galaxy information can be used to purify samples of cosmological Type Ia supernovae, uncover the progenitor systems of individual classes, and facilitate low-latency follow-up of rare and peculiar explosions. In this work, we develop a novel data-driven methodology to simulate the time-domain sky that includes detailed modelling of the probability density function for multiple transient classes conditioned on host galaxy magnitudes, colours, star formation rates, and masses. We have designed these simulations to optimize photometric classification and analysis in upcoming large synoptic surveys. We integrate host galaxy information into the snana simulation framework to construct the simulated catalogue of optical transients and correlated hosts (SCOTCH, a publicly available catalogue of 5-million idealized transient light curves in LSST passbands and their host galaxy properties over the redshift range 0 < z < 3. This catalogue includes supernovae, tidal disruption events, kilonovae, and active galactic nuclei. Each light curve consists of true top-of-the-galaxy magnitudes sampled with high (≲2 d) cadence. In conjunction with SCOTCH, we also release an associated set of tutorials and transient-specific libraries to enable simulations of arbitrary space- and ground-based surveys. Our methodology is being used to test critical science infrastructure in advance of surveys by the Vera C. Rubin Observatory and the Nancy G. Roman Space Telescope.

     
    more » « less