skip to main content

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

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
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 2545-2560
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(, 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 We measured the precise masses of the host and planet in the OGLE-2003-BLG-235 system, when the lens and source were resolving, with 2018 Keck high resolution images. This measurement is in agreement with the observation taken in 2005 with the Hubble Space Telescope (HST). In the 2005 data, the lens and sources were not resolved and the measurement was made using color-dependent centroid shift only. The Nancy Grace Roman Space Telescope will measure masses using data typically taken within 3–4 yr of the peak of the event, which is a much shorter baseline when compared to most of the mass measurements to date. Hence, the color-dependent centroid shift will be one of the primary methods of mass measurements for the Roman telescope. Yet, mass measurements of only two events (OGLE-2003-BLG-235 and OGLE-2005-BLG-071) have been done using the color-dependent centroid shift method so far. The accuracy of the measurements using this method are neither completely known nor well studied. The agreement of the Keck and HST results, as shown in this paper, is very important because this agreement confirms the accuracy of the mass measurements determined at a small lens-source separation using the color-dependent centroid shift method. It also shows that with >100 high resolution images, the Roman telescope will be able to use color-dependent centroid shift at a 3–4 yr time baseline and produce mass measurements. We find that OGLE-2003-BLG-235 is a planetary system that consists of a 2.34 ± 0.43 M Jup planet orbiting a 0.56 ± 0.06 M ⊙ K-dwarf host star at a distance of 5.26 ± 0.71 kpc from the Sun. 
    more » « less