Showing 1–50 of 117 results for author: Lu, J R

Recommendations for Early Definition Science with the Nancy Grace Roman Space Telescope

Authors: Robyn E. Sanderson, Ryan Hickox, Christopher M. Hirata, Matthew J. Holman, Jessica R. Lu, Ashley Villar

Abstract: The Nancy Grace Roman Space Telescope (Roman), NASA's next flagship observatory, has significant mission time to be spent on surveys for general astrophysics in addition to its three core community surveys. We considered what types of observations outside the core surveys would most benefit from early definition, given 700 hours of mission time in the first two years of Roman's operation. We recom… ▽ More The Nancy Grace Roman Space Telescope (Roman), NASA's next flagship observatory, has significant mission time to be spent on surveys for general astrophysics in addition to its three core community surveys. We considered what types of observations outside the core surveys would most benefit from early definition, given 700 hours of mission time in the first two years of Roman's operation. We recommend that a survey of the Galactic plane be defined early, based on the broad range of stakeholders for such a survey, the added scientific value of a first pass to obtain a baseline for proper motions complementary to Gaia's, and the significant potential synergies with ground-based surveys, notably the Legacy Survey of Space and Time (LSST) on Rubin. We also found strong motivation to follow a community definition process for ultra-deep observations with Roman. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: Final report of the committee convened by NASA to assess community proposals for early definition science with Roman

Unveiling MOA-2007-BLG-192: An M Dwarf Hosting a Likely Super-Earth

Authors: Sean K. Terry, Jean-Philippe Beaulieu, David P. Bennett, Euan Hamdorf, Aparna Bhattacharya, Viveka Chaudhry, Andrew A. Cole, Naoki Koshimoto, Jay Anderson, Etienne Bachelet, Joshua W. Blackman, Ian A. Bond, Jessica R. Lu, Jean Baptiste Marquette, Clement Ranc, Natalia E. Rektsini, Kailash Sahu, Aikaterini Vandorou

Abstract: We present an analysis of high angular resolution images of the microlensing target MOA-2007-BLG-192 using Keck adaptive optics and the Hubble Space Telescope. The planetary host star is robustly detected as it separates from the background source star in nearly all of the Keck and Hubble data. The amplitude and direction of the lens-source separation allows us to break a degeneracy related to the… ▽ More We present an analysis of high angular resolution images of the microlensing target MOA-2007-BLG-192 using Keck adaptive optics and the Hubble Space Telescope. The planetary host star is robustly detected as it separates from the background source star in nearly all of the Keck and Hubble data. The amplitude and direction of the lens-source separation allows us to break a degeneracy related to the microlensing parallax and source radius crossing time. Thus, we are able to reduce the number of possible solutions by a factor of ${\sim}2$, demonstrating the power of high angular resolution follow-up imaging for events with sparse light curve coverage. Following Bennett et al. 2023, we apply constraints from the high resolution imaging on the light curve modeling to find host star and planet masses of $M_{\textrm{host}} = 0.28 \pm 0.04M_{\odot}$ and $m_p = 12.49^{+65.47}_{-8.03}M_{\oplus}$ at a distance from Earth of $D_L = 2.16 \pm 0.30\,$kpc. This work illustrates the necessity for the Nancy Grace Roman Galactic Exoplanet Survey (RGES) to use its own high resolution imaging to inform light curve modeling for microlensing planets that the mission discovers. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 26 pages, 10 figures, submitted to AJ

An Estimate of the Binary Star Fraction Among Young Stars at the Galactic Center: Possible Evidence of a Radial Dependence

Authors: Abhimat K. Gautam, Tuan Do, Andrea M. Ghez, Devin S. Chu, Matthew W. Hosek Jr., Shoko Sakai, Smadar Naoz, Mark R. Morris, Anna Ciurlo, Zoe Haggard, Jessica R. Lu

Abstract: We present the first estimate of the intrinsic binary fraction of young stars across the central $\approx$ 0.4 pc surrounding the supermassive black hole (SMBH) at the Milky Way Galactic center (GC). This experiment searched for photometric variability in 102 young stars, using 119 nights of 10"-wide adaptive optics imaging observations taken at Keck Observatory over 16 years in the K'- and H-band… ▽ More We present the first estimate of the intrinsic binary fraction of young stars across the central $\approx$ 0.4 pc surrounding the supermassive black hole (SMBH) at the Milky Way Galactic center (GC). This experiment searched for photometric variability in 102 young stars, using 119 nights of 10"-wide adaptive optics imaging observations taken at Keck Observatory over 16 years in the K'- and H-bands. We photometrically detected three binary stars, all of which are situated more than 1" (0.04 pc) from the SMBH and one of which, S2-36, is newly reported here with spectroscopic confirmation. To convert the observed binary fraction into an estimate of the underlying binary fraction, we determined experiment sensitivity through detailed light curve simulations, incorporating photometric effects of eclipses, irradiation, and tidal distortion in binaries. The simulations assumed a population of young binaries, with stellar ages (4 Myr) and masses matched to the most probable values measured for the GC young star population and underlying binary system parameters similar to those of local massive stars. The detections and simulations imply young, massive stars in the GC have a stellar binary fraction $\geq$ 71% (68% confidence), or $\geq$ 42% (95% confidence). This inferred GC young star binary fraction is consistent with that typically seen in young stellar populations in the solar neighborhood. Furthermore, our measured binary fraction is significantly higher than that recently reported by Chu et al. (2023) based on RV measurements of young stars <~1" of the SMBH. Constrained with these two studies, the probability that the same underlying young binary fraction extends across the entire region is <1.4%. This tension provides support for a radial dependence of the binary star fraction and, therefore, for the dynamical predictions of binary merger and evaporation events close to the SMBH. △ Less

Submitted 7 February, 2024; v1 submitted 23 January, 2024; originally announced January 2024.

Comments: 51 pages, 27 figures, 7 tables. Accepted for publication in The Astrophysical Journal. Abstract revised for final manuscript

Astrometric Microlensing by Primordial Black Holes with The Roman Space Telescope

Authors: James Fardeen, Peter McGill, Scott E. Perkins, William A. Dawson, Natasha S. Abrams, Jessica R. Lu, Ming-Feng Ho, Simeon Bird

Abstract: Primordial Black Holes (PBHs) could explain some fraction of dark matter and shed light on many areas of early-universe physics. Despite over half a century of research interest, a PBH population has so far eluded detection. The most competitive constraints on the fraction of dark matter comprised of PBHs ($f_{\rm DM}$) in the $(10^{-9}-10)M_{\odot}$ mass-ranges come from photometric microlensing… ▽ More Primordial Black Holes (PBHs) could explain some fraction of dark matter and shed light on many areas of early-universe physics. Despite over half a century of research interest, a PBH population has so far eluded detection. The most competitive constraints on the fraction of dark matter comprised of PBHs ($f_{\rm DM}$) in the $(10^{-9}-10)M_{\odot}$ mass-ranges come from photometric microlensing and bound $f_{\rm DM}\lesssim10^{-2}-10^{-1}$. With the advent of the Roman Space Telescope with its sub-milliarcsecond (mas) astrometric capabilities and its planned Galactic Bulge Time Domain Survey (GBTDS), detecting astrometric microlensing signatures will become routine. Compared with photometric microlensing, astrometric microlensing signals are sensitive to different lens masses-distance configurations and contains different information, making it a complimentary lensing probe. At sub-mas astrometric precision, astrometric microlensing signals are typically detectable at larger lens-source separations than photometric signals, suggesting a microlensing detection channel of pure astrometric events. We use a Galactic simulation to predict the number of detectable microlensing events during the GBTDS via this pure astrometric microlensing channel. Assuming an absolute astrometric precision floor for bright stars of 0.1 mas for the GBTDS, we find that the number of detectable events peaks at $\approx 10^{3} f_{\rm DM}$ for a population of $ 1 M_{\odot}$ PBHs and tapers to $\approx 10f_{\rm DM}$ and $\approx 100f_{\rm DM}$ at $10^{-4}M_{\odot}$ and $10^{3}M_{\odot}$, respectively. Accounting for the distinguishability of PBHs from Stellar lenses, we conclude the GBTDS will be sensitive to a PBH population at $f_{\rm DM}$ down to $\approx10^{-1}-10^{-3}$ for $(10^{-1}-10^{2})M_{\odot}$ likely yielding novel PBH constraints. △ Less

Submitted 8 March, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 21 pages, 11 figures, accepted to AAS Journals

The 3D Kinematics of the Orion Nebula Cluster II: Mass-dependent Kinematics of the Inner Cluster

Authors: Lingfeng Wei, Christopher A. Theissen, Quinn M. Konopacky, Jessica R. Lu, Chih-Chun Hsu, Dongwon Kim

Abstract: We present the kinematic anaylsis of $246$ stars within $4^\prime$ from the center of Orion Nebula Cluster (ONC), the closest massive star cluster with active star formation across the full mass range, which provides valuable insights in the the formation and evolution of star cluster on an individual-star basis. High-precision radial velocities and surface temperatures are retrieved from spectra… ▽ More We present the kinematic anaylsis of $246$ stars within $4^\prime$ from the center of Orion Nebula Cluster (ONC), the closest massive star cluster with active star formation across the full mass range, which provides valuable insights in the the formation and evolution of star cluster on an individual-star basis. High-precision radial velocities and surface temperatures are retrieved from spectra acquired by the NIRSPEC instrument used with adaptive optics (NIRSPAO) on the Keck II 10-m telescope. A three-dimensional kinematic map is then constructed by combining with the proper motions previously measured by the Hubble Space Telescope (HST) ACS/WFPC2/WFC3IR and Keck II NIRC2. The measured root-mean-squared velocity dispersion is $2.26\pm0.08~\mathrm{km}\,\mathrm{s}^{-1}$, significantly higher than the virial equilibrium's requirement of $1.73~\mathrm{km}\,\mathrm{s}^{-1}$, suggesting that the ONC core is supervirial, consistent with previous findings. Energy equipartition is not detected in the cluster. Most notably, the velocity of each star relative to its neighbors is found to be negatively correlated with stellar mass. Low-mass stars moving faster than their surrounding stars in a supervirial cluster suggests that the initial masses of forming stars may be related to their initial kinematic states. Additionally, a clockwise rotation preference is detected. A weak sign of inverse mass segregation is also identified among stars excluding the Trapezium stars, though it could be a sample bias. Finally, this study reports the discovery of four new candidate spectroscopic binary systems. △ Less

Submitted 6 March, 2024; v1 submitted 7 December, 2023; originally announced December 2023.

Comments: Published by ApJ; 28 pages, 20 figures, 3 tables

Journal ref: 2024 ApJ 962 174

Keck and Hubble Observations Show That MOA-2008-BLG-379Lb Is a Super-Jupiter Orbiting an M Dwarf

Authors: David P. Bennett, Aparna Bhattacharya, Jean-Philippe Beaulieu, Naoki Koshimoto, Joshua W. Blackman, Ian A. Bond, Clement Ranc, Natalia Rektsini, Sean K. Terry, Aikaterini Vandorou, Jessica R. Lu, Jean Baptiste Marquette, Greg Olmschenk, Daisuke Suzuki

Abstract: We present high angular resolution imaging that detects the MOA-2008-BLG-379L exoplanet host star using Keck adaptive optics and the Hubble Space Telescope. These observations reveal host star and planet masses of $M_{\rm host}=0.434\pm0.065 M_\odot$, and $m_p=2.44 \pm 0.49 M_{\rm Jupiter}$. They are located at a distance of $D_L=3.44\pm0.53\,$kpc, with a projected separation of $2.70\pm 0.42\,$AU… ▽ More We present high angular resolution imaging that detects the MOA-2008-BLG-379L exoplanet host star using Keck adaptive optics and the Hubble Space Telescope. These observations reveal host star and planet masses of $M_{\rm host}=0.434\pm0.065 M_\odot$, and $m_p=2.44 \pm 0.49 M_{\rm Jupiter}$. They are located at a distance of $D_L=3.44\pm0.53\,$kpc, with a projected separation of $2.70\pm 0.42\,$AU. These results contribute to our determination of exoplanet host star masses for the Suzuki et al. (2016) statistical sample, which will determine the dependence of the planet occurrence rate on the mass and distance of the host stars. We also present a detailed discussion of the image constrained modeling version of the eesunhong light curve modeling code that applies high angular resolution image constraints to the light curve modeling process. This code increases modeling efficiency by a large factor by excluding models that are inconsistent with the high angular resolution images. The analysis of this and other events from the Suzuki et al. (2016) statistical sample reveals the importance of including higher order effects, such as microlensing parallax and planetary orbital motion even when these features are not required to fit the light curve data. The inclusion of these effects may be needed to obtain accurate estimates of the uncertainty of other microlensing parameters that affect the inferred properties of exoplanet microlens systems. This will be important for the exoplanet microlensing survey of the Roman Space Telescope, which will use both light curve photometry and high angular resolution imaging to characterize planetary microlens systems. △ Less

Submitted 7 May, 2024; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: 42 pages, with 4 figures, accepted for publication in the Astronomical Journal

The JWST Galactic Center Survey -- A White Paper

Authors: Rainer Schoedel, Steve Longmore, Jonny Henshaw, Adam Ginsburg, John Bally, Anja Feldmeier, Matt Hosek, Francisco Nogueras Lara, Anna Ciurlo, Mélanie Chevance, J. M. Diederik Kruijssen, Ralf Klessen, Gabriele Ponti, Pau Amaro-Seoane, Konstantina Anastasopoulou, Jay Anderson, Maria Arias, Ashley T. Barnes, Cara Battersby, Giuseppe Bono, Lucía Bravo Ferres, Aaron Bryant, Miguel Cano Gonzáalez, Santi Cassisi, Leonardo Chaves-Velasquez , et al. (85 additional authors not shown)

Abstract: The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of… ▽ More The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\,pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value. △ Less

Submitted 14 March, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: This White Paper will be updated when required (e.g. new authors joining, editing of content). Most recent update: 24 Oct 2023

Disentangling the Black Hole Mass Spectrum with Photometric Microlensing Surveys

Authors: Scott Ellis Perkins, Peter McGill, William Dawson, Natasha S. Abrams, Casey Y. Lam, Ming-Feng Ho, Jessica R. Lu, Simeon Bird, Kerianne Pruett, Nathan Golovich, George Chapline

Abstract: From the formation mechanisms of stars and compact objects to nuclear physics, modern astronomy frequently leverages surveys to understand populations of objects to answer fundamental questions. The population of dark and isolated compact objects in the Galaxy contains critical information related to many of these topics, but is only practically accessible via gravitational microlensing. However,… ▽ More From the formation mechanisms of stars and compact objects to nuclear physics, modern astronomy frequently leverages surveys to understand populations of objects to answer fundamental questions. The population of dark and isolated compact objects in the Galaxy contains critical information related to many of these topics, but is only practically accessible via gravitational microlensing. However, photometric microlensing observables are degenerate for different types of lenses, and one can seldom classify an event as involving either a compact object or stellar lens on its own. To address this difficulty, we apply a Bayesian framework that treats lens type probabilistically and jointly with a lens population model. This method allows lens population characteristics to be inferred despite intrinsic uncertainty in the lens-class of any single event. We investigate this method's effectiveness on a simulated ground-based photometric survey in the context of characterizing a hypothetical population of primordial black holes (PBHs) with an average mass of $30 M_{\odot}$. On simulated data, our method outperforms current black hole (BH) lens identification pipelines and characterizes different subpopulations of lenses while jointly constraining the PBH contribution to dark matter to ${\approx}25$\%. Key to robust inference, our method can marginalize over population model uncertainty. We find the lower mass cutoff for stellar origin BHs, a key observable in understanding the BH mass gap, particularly difficult to infer in our simulations. This work lays the foundation for cutting-edge PBH abundance constraints to be extracted from current photometric microlensing surveys. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 31 pages, 18 figures, submitted to AAS

Microlensing Discovery and Characterization Efficiency in the Vera C. Rubin Legacy Survey of Space and Time

Authors: Natasha S. Abrams, Markus P. G. Hundertmark, Somayeh Khakpash, Rachel A. Street, R. Lynne Jones, Jessica R. Lu, Etienne Bachelet, Yiannis Tsapras, Marc Moniez, Tristan Blaineauu, Rosanne Di Stefano, Martin Makler, Anibal Varela, Markus Rabus

Abstract: The Vera C. Rubin Legacy Survey of Space and Time will discover thousands of microlensing events across the Milky Way Galaxy, allowing for the study of populations of exoplanets, stars, and compact objects. It will reach deeper limiting magnitudes over a wider area than any previous survey. We evaluate numerous survey strategies simulated in the Rubin Operation Simulations (OpSims) to assess the d… ▽ More The Vera C. Rubin Legacy Survey of Space and Time will discover thousands of microlensing events across the Milky Way Galaxy, allowing for the study of populations of exoplanets, stars, and compact objects. It will reach deeper limiting magnitudes over a wider area than any previous survey. We evaluate numerous survey strategies simulated in the Rubin Operation Simulations (OpSims) to assess the discovery and characterization efficiencies of microlensing events. We have implemented three metrics in the Rubin Metric Analysis Framework: a discovery metric and two characterization metrics, where one estimates how well the lightcurve is covered and the other quantifies how precisely event parameters can be determined. We also assess the characterizability of microlensing parallax, critical for detection of free-floating black hole lenses, in a representative bulge and disk field. We find that, given Rubin's baseline cadence, the discovery and characterization efficiency will be higher for longer duration and larger parallax events. Microlensing discovery efficiency is dominated by observing footprint, where more time spent looking at regions of high stellar density including the Galactic bulge, Galactic plane, and Magellanic clouds, leads to higher discovery and characterization rates. However, if the observations are stretched over too wide an area, including low-priority areas of the Galactic plane with fewer stars and higher extinction, event characterization suffers by > 10%, which could impact exoplanet, binary star, and compact object events alike. We find that some rolling strategies (where Rubin focuses on a fraction of the sky in alternating years) in the Galactic bulge can lead to a 15-20% decrease in microlensing parallax characterization, so rolling strategies should be chosen carefully to minimize losses. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 21 pages, 15 figures, submitted to ApJS Rubin Survey Strategy Edition

Near-Infrared Flux Distribution of Sgr A* from 2005-2022: Evidence for an Enhanced Accretion Episode in 2019

Authors: Grant C. Weldon, Tuan Do, Gunther Witzel, Andrea M. Ghez, Abhimat K. Gautam, Eric E. Becklin, Mark R. Morris, Gregory D. Martinez, Shoko Sakai, Jessica R. Lu, Keith Matthews, Matthew W. Hosek Jr., Zoë Haggard

Abstract: Sgr A* is the variable electromagnetic source associated with accretion onto the Galactic center supermassive black hole. While the near-infrared (NIR) variability of Sgr A* was shown to be consistent over two decades, unprecedented activity in 2019 challenges existing statistical models. We investigate the origin of this activity by re-calibrating and re-analyzing all of our Keck Observatory Sgr… ▽ More Sgr A* is the variable electromagnetic source associated with accretion onto the Galactic center supermassive black hole. While the near-infrared (NIR) variability of Sgr A* was shown to be consistent over two decades, unprecedented activity in 2019 challenges existing statistical models. We investigate the origin of this activity by re-calibrating and re-analyzing all of our Keck Observatory Sgr A* imaging observations from 2005-2022. We present light curves from 69 observation epochs using the NIRC2 imager at 2.12 $μ$m with laser guide star adaptive optics. These observations reveal that the mean luminosity of Sgr A* increased by a factor of $\sim$3 in 2019, and the 2019 light curves had higher variance than in all time periods we examined. We find that the 2020-2022 flux distribution is statistically consistent with the historical sample and model predictions, but with fewer bright measurements above 0.6 mJy at the $\sim$2$σ$ level. Since 2019, we have observed a maximum $K_s$ (2.2 $μ$m) flux of 0.9 mJy, compared to the highest pre-2019 flux of 2.0 mJy and highest 2019 flux of 5.6 mJy. Our results suggest that the 2019 activity was caused by a temporary accretion increase onto Sgr A*, possibly due to delayed accretion of tidally-stripped gas from the gaseous object G2 in 2014. We also examine faint Sgr A* fluxes over a long time baseline to search for a quasi-steady quiescent state. We find that Sgr A* displays flux variations over a factor of $\sim$500, with no evidence for a quiescent state in the NIR. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: Accepted to ApJ Letters, 26 pages, 15 figures

A re-analysis of the isolated black hole candidate OGLE-2011-BLG-0462/MOA-2011-BLG-191

Authors: Casey Y. Lam, Jessica R. Lu

Abstract: There are expected to be $\sim 10^8$ isolated black holes (BHs) in the Milky Way. OGLE-2011-BLG-0462/MOA-2011-BLG-191 (OB110462) is the only such BH with a mass measurement to date. However, its mass is disputed: Lam et al. (2022a,b) measured a lower mass of $1.6 - 4.4 M_\odot$, while Sahu et al. (2022); Mróz et al. (2022) measured a higher mass of $5.8 - 8.7 M_\odot$. We re-analyze OB110462, incl… ▽ More There are expected to be $\sim 10^8$ isolated black holes (BHs) in the Milky Way. OGLE-2011-BLG-0462/MOA-2011-BLG-191 (OB110462) is the only such BH with a mass measurement to date. However, its mass is disputed: Lam et al. (2022a,b) measured a lower mass of $1.6 - 4.4 M_\odot$, while Sahu et al. (2022); Mróz et al. (2022) measured a higher mass of $5.8 - 8.7 M_\odot$. We re-analyze OB110462, including new data from the Hubble Space Telescope (HST) and re-reduced Optical Gravitational Lensing Experiment (OGLE) photometry. We also re-reduce and re-analyze the HST dataset with newly available software. We find significantly different ($\sim 1$ mas) HST astrometry than Lam et al. (2022a,b) in the de-magnified epochs due to the amount of positional bias induced by a bright star $\sim$0.4 arcsec from OB110462. After modeling the updated photometric and astrometric datasets, we find the lens of OB110462 is a $6.0^{+1.2}_{-1.0} M_\odot$ BH. Future observations with the Nancy Grace Roman Space Telescope, which will have an astrometric precision comparable or better to HST but a field of view $100\times$ larger, will be able to measure hundreds of isolated BH masses via microlensing. This will enable the measurement of the BH mass distribution and improve understanding of massive stellar evolution and BH formation channels. △ Less

Submitted 8 August, 2023; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 23 pages, 18 figures, 8 tables. Accepted for publication in ApJ on 2 Aug 2023 [Same as v1, just fixed typo in email address]

Maximizing science return by coordinating the survey strategies of Roman with Rubin, and other major facilities

Authors: R. A. Street, S. Gough-Kelly, C. Lam, A. Varela, M. Makler, E. Bachelet, J. R. Lu, N. Abrams, A. Pusack, S. Terry, R. Di~Stefano, Y. Tsapras, M. P. G. Hundertmark, R. J. J. Grand, T. Daylan, J. Sobeck

Abstract: [Abridged] The Nancy Grace Roman Space Telescope will be one of several flagship survey facilities operating over the next decade starting $\sim$2025. The deep near-IR imaging that Roman will deliver will be highly complementary to the capabilities of other survey telescopes that will operate contemporaneously, particularly those that can provide data at different wavelengths and messengers, or di… ▽ More [Abridged] The Nancy Grace Roman Space Telescope will be one of several flagship survey facilities operating over the next decade starting $\sim$2025. The deep near-IR imaging that Roman will deliver will be highly complementary to the capabilities of other survey telescopes that will operate contemporaneously, particularly those that can provide data at different wavelengths and messengers, or different time intervals. Combining data from multiple facilities can provide important astrophysical insights, provided the data acquisition is carefully scheduled, and careful plans are made for appropriate joint data analyses. In this White Paper, we discuss the broad range of science that would be enabled by coordinating Roman observations of the Galactic Bulge with those of the Vera C. Rubin Observatory. Specifically, we discuss how Roman's characterization of lensing events caused by exoplanets, stellar systems and stellar remnants can be enhanced by data from Rubin. The same data will also be highly advantageous for the determination of stellar properties, and for distinguishing exoplanetary transits. It will enable more accurate period-color-luminosity relationships to be measured for RR~Lyrae throughout the Milky Way Bulge and Bar, probing galactic structure and dynamics. But we stress that this is only a sample of the full potential and advocate for a more complete study to be made as a joint effort between these major projects. We note that we do not suggest any changes beyond the established Science Requirements for the RGBTDS, in terms of survey footprint or filter selection. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: Submitted to the 2023 call for White Papers regarding the Roman Mission Core Community Surveys

Roman CCS White Paper: Characterizing the Galactic population of isolated black holes

Authors: Casey Y. Lam, Natasha Abrams, Jeff Andrews, Etienne Bachelet, Arash Bahramian, David Bennett, Valerio Bozza, Floor Broekgaarden, Sukanya Chakrabarti, William Dawson, Kareem El-Badry, Maya Fishbach, Giacomo Fragione, Scott Gaudi, Abhimat Gautam, Ryosuke Hirai, Daniel Holz, Matthew Hosek Jr., Macy Huston, Tharindu Jayasinghe, Samson Johnson, Daisuke Kawata, Naoki Koshimoto, Jessica R. Lu, Ilya Mandel , et al. (12 additional authors not shown)

Abstract: Although there are estimated to be 100 million isolated black holes (BHs) in the Milky Way, only one has been found so far, resulting in significant uncertainty about their properties. The Galactic Bulge Time Domain Survey provides the only opportunity in the coming decades to grow this catalog by order(s) of magnitude. This can be achieved if 1) Roman's astrometric potential is fully realized in… ▽ More Although there are estimated to be 100 million isolated black holes (BHs) in the Milky Way, only one has been found so far, resulting in significant uncertainty about their properties. The Galactic Bulge Time Domain Survey provides the only opportunity in the coming decades to grow this catalog by order(s) of magnitude. This can be achieved if 1) Roman's astrometric potential is fully realized in the observation strategy and software pipelines, 2) Roman's observational gaps of the Bulge are minimized, and 3) observations with ground-based facilities are taken of the Bulge to fill in gaps during non-Bulge seasons. A large sample of isolated BHs will enable a broad range of astrophysical questions to be answered, such as massive stellar evolution, origin of gravitational wave sources, supernova physics, and the growth of supermassive BHs, maximizing Roman's scientific return. △ Less

Submitted 29 June, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 20 pages. Submitted in response to Nancy Grace Roman Space Telescope white paper call: https://roman.gsfc.nasa.gov/science/ccs_white_papers.html. v2 fixes a typo in Figure 5 axis label (days --> year)

The Galactic Center with Roman

Authors: Sean K. Terry, Matthew W. Hosek Jr., Jessica R. Lu, Casey Lam, Natasha Abrams, Arash Bahramian, Richard Barry, Jean-Phillipe Beaulieu, Aparna Bhattacharya, Devin Chu, Anna Ciurlo, Will Clarkson, Tuan Do, Kareem El-Badry, Ryan Felton, Matthew Freeman, Abhimat Gautam, Andrea Ghez, Daniel Huber, Jason Hunt, Macy Huston, Tharindu Jayasinghe, Naoki Koshimoto, Madeline Lucey, Florian Peißker , et al. (9 additional authors not shown)

Abstract: We advocate for a Galactic center (GC) field to be added to the Galactic Bulge Time Domain Survey (GBTDS). The new field would yield high-cadence photometric and astrometric measurements of an unprecedented ${\sim}$3.3 million stars toward the GC. This would enable a wide range of science cases, such as finding star-compact object binaries that may ultimately merge as LISA-detectable gravitational… ▽ More We advocate for a Galactic center (GC) field to be added to the Galactic Bulge Time Domain Survey (GBTDS). The new field would yield high-cadence photometric and astrometric measurements of an unprecedented ${\sim}$3.3 million stars toward the GC. This would enable a wide range of science cases, such as finding star-compact object binaries that may ultimately merge as LISA-detectable gravitational wave sources, constraining the mass function of stars and compact objects in different environments, detecting populations of microlensing and transiting exoplanets, studying stellar flares and variability in young and old stars, and monitoring accretion onto the central supermassive black hole. In addition, high-precision proper motions and parallaxes would open a new window into the large-scale dynamics of stellar populations at the GC, yielding insights into the formation and evolution of galactic nuclei and their co-evolution with the growth of the supermassive black hole. We discuss the possible trade-offs between the notional GBTDS and the addition of a GC field with either an optimal or minimal cadence. Ultimately, the addition of a GC field to the GBTDS would dramatically increase the science return of Roman and provide a legacy dataset to study the mid-plane and innermost regions of our Galaxy. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 19 pages, 3 figures. Submitted to the NASA Roman Core Community Surveys White Paper Call

Evidence of a Decreased Binary Fraction for Massive Stars Within 20 Milliparsecs of the Supermassive Black Hole at the Galactic Center

Authors: Devin S. Chu, Tuan Do, Andrea Ghez, Abhimat K. Gautam, Anna Ciurlo, Kelly Kosmo O'Neil, Matthew W. Hosek Jr., Aurélien Hees, Smadar Naoz, Shoko Sakai, Jessica R. Lu, Zhuo Chen, Rory O. Bentley, Eric E. Becklin, Keith Matthews

Abstract: We present the results of the first systematic search for spectroscopic binaries within the central 2 x 3 arcsec$^2$ around the supermassive black hole at the center of the Milky Way galaxy. This survey is based primarily on over a decade of adaptive optics-fed integral-field spectroscopy (R$\sim$4000), obtained as part of the Galactic Center Orbits Initiative at Keck Observatory, and has a limiti… ▽ More We present the results of the first systematic search for spectroscopic binaries within the central 2 x 3 arcsec$^2$ around the supermassive black hole at the center of the Milky Way galaxy. This survey is based primarily on over a decade of adaptive optics-fed integral-field spectroscopy (R$\sim$4000), obtained as part of the Galactic Center Orbits Initiative at Keck Observatory, and has a limiting $K$'-band magnitude of 15.8, which is at least 4 magnitudes deeper than previous spectroscopic searches for binaries at larger radii within the central nuclear star cluster. From this primary dataset, over 600 new radial velocities are extracted and reported, increasing by a factor of 3 the number of such measurements. We find no significant periodic signals in our sample of 28 stars, of which 16 are massive, young (main-sequence B) stars and 12 are low-mass, old (M and K giant) stars. Using Monte Carlo simulations, we derive upper limits on the intrinsic binary star fraction for the young star population at 47% (at 95% confidence) located $\sim$20 mpc from the black hole. The young star binary fraction is significantly lower than that observed in the field (70%). This result is consistent with a scenario in which the central supermassive black hole drives nearby stellar binaries to merge or be disrupted and may have important implications for the production of gravitational waves and hypervelocity stars. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: Accepted by ApJ. 22 pages, 14 figures

Smartphone screens as astrometric calibrators

Authors: Aidan Walk, Charles-Antoine Claveau, Michael Bottom, Mark Chun, Shane Jacobson, Maxwell Service, Jessica R. Lu

Abstract: Geometric optical distortion is a significant contributor to the astrometric error budget in large telescopes using adaptive optics. To increase astrometric precision, optical distortion calibration is necessary. We investigate using smartphone OLED screens as astrometric calibrators. Smartphones are low cost, have stable illumination, and can be quickly reconfigured to probe different spatial fre… ▽ More Geometric optical distortion is a significant contributor to the astrometric error budget in large telescopes using adaptive optics. To increase astrometric precision, optical distortion calibration is necessary. We investigate using smartphone OLED screens as astrometric calibrators. Smartphones are low cost, have stable illumination, and can be quickly reconfigured to probe different spatial frequencies of an optical system's geometric distortion. In this work, we characterize the astrometric accuracy of a Samsung S20 smartphone, with a view towards providing large format, flexible astrometric calibrators for the next generation of astronomical instruments. We find the placement error of the pixels to be 189 nm +/- 15 nm RMS. At this level of error, milliarcsecond astrometric accuracy can be obtained on modern astronomical instruments. △ Less

Submitted 6 July, 2023; v1 submitted 8 March, 2023; originally announced March 2023.

Comments: 15 pages, 11 figures; accepted, Journal of Astronomical Instrumentation

Stellar Populations in the Central 0.5 pc of Our Galaxy III: The Dynamical Sub-structures

Authors: Siyao Jia, Ningyuan Xu, Jessica R. Lu, D. S Chu, K. Kosmo O'Neil, W. B. Drechsler, M. W. Hosek Jr., S. Sakai, T. Do, A. Ciurlo, A. K. Gautam, A. M. Ghez, E. Becklin, M. R. Morris, R. O. Bentley

Abstract: We measure the 3D kinematic structures of the young stars within the central 0.5 parsec of our Galactic Center using the 10 m telescopes of the W.~M.~Keck Observatory over a time span of 25 years. Using high-precision measurements of positions on the sky, and proper motions and radial velocities from new observations and the literature, we constrain the orbital parameters for each young star. Our… ▽ More We measure the 3D kinematic structures of the young stars within the central 0.5 parsec of our Galactic Center using the 10 m telescopes of the W.~M.~Keck Observatory over a time span of 25 years. Using high-precision measurements of positions on the sky, and proper motions and radial velocities from new observations and the literature, we constrain the orbital parameters for each young star. Our results show two statistically significant sub-structures: a clockwise stellar disk with 18 candidate stars, as has been proposed before, but with an improved disk membership; a second, almost edge-on plane of 10 candidate stars oriented East-West on the sky that includes at least one IRS 13 star. We estimate the eccentricity distribution of each sub-structure and find that the clockwise disk has <$e$> = 0.39 and the edge-on plane has <$e$> = 0.68. We also perform simulations of each disk/plane with incompleteness and spatially-variable extinction to search for asymmetry. Our results show that the clockwise stellar disk is consistent with a uniform azimuthal distribution within the disk. The edge-on plane has an asymmetry that cannot be explained by variable extinction or incompleteness in the field. The orientation, asymmetric stellar distribution, and high eccentricity of the edge-on plane members suggest that this structure may be a stream associated with the IRS 13 group. The complex dynamical structure of the young nuclear cluster indicates that the star formation process involved complex gas structures and dynamics and is inconsistent with a single massive gaseous disk. △ Less

Submitted 3 February, 2023; originally announced February 2023.

Comments: 41 pages, 26 figures, 13 tables, 2 appendices. Accepted for publication in ApJ

The Swansong of the Galactic Center Source X7: An Extreme Example of Tidal Evolution near the Supermassive Black Hole

Authors: Anna Ciurlo, Randall D. Campbell, Mark R. Morris, Tuan Do, Andrea M. Ghez, Eric E. Becklin, Rory O. Bentley, Devin S. Chu, Abhimat K. Gautam, Yash A. Gursahani, Aurelien Hees, Kelly Kosmo O'Neil, Jessica R. Lu, Gregory D. Martinez, Smadar Naoz, Shoko Sakai, Rainer Schoedel

Abstract: We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-γ line-emission and Lp (3.8 μm) imaging data, we provide the first estimate of its orbital p… ▽ More We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-γ line-emission and Lp (3.8 μm) imaging data, we provide the first estimate of its orbital parameters and quantitative characterization of the evolution of its morphology and mass. We find that the leading edge of X7 appears to be on a mildly eccentric (e~0.3), relatively short-period (170 years) orbit and is headed towards periapse passage, estimated to occur in ~2036. Furthermore, our kinematic measurements rule out the earlier suggestion that X7 is associated with the stellar source S0-73 or with any other point source that has overlapped with X7 during our monitoring period. Over the course of our observations, X7 has (1) become more elongated, with a current length-to-width ratio of 9, (2) maintained a very consistent long-axis orientation (position angle of 50 deg), (3) inverted its radial velocity differential from tip to tail from -50 to +80 km/sec, and (4) sustained its total brightness (12.8 Lp magnitudes at the leading edge) and color temperature (425 K), which suggest a constant mass of ~50 MEarth. We present a simple model showing that these results are compatible with the expected effect of tidal forces exerted on it by the central black hole and we propose that X7 is the gas and dust recently ejected from a grazing collision in a binary system. △ Less

Submitted 16 January, 2023; originally announced January 2023.

Comments: accepted by ApJ, 20 pages, 16 figures

The Star Formation History of the Milky Way's Nuclear Star Cluster

Authors: Zhuo Chen, Tuan Do, Andrea M. Ghez, Matthew Hosek Jr., Anja Feldmeier-Krause, Devin Chu, Rory Bentley, Jessica R. Lu, Mark R. Morris

Abstract: We report the first star formation history study of the Milky Way's nuclear star cluster (NSC) that includes observational constraints from a large sample of stellar metallicity measurements. These metallicity measurements were obtained from recent surveys from Gemini and VLT of 770 late-type stars within the central 1.5 pc. These metallicity measurements, along with photometry and spectroscopical… ▽ More We report the first star formation history study of the Milky Way's nuclear star cluster (NSC) that includes observational constraints from a large sample of stellar metallicity measurements. These metallicity measurements were obtained from recent surveys from Gemini and VLT of 770 late-type stars within the central 1.5 pc. These metallicity measurements, along with photometry and spectroscopically derived temperatures, are forward modeled with a Bayesian inference approach. Including metallicity measurements improves the overall fit quality, as the low-temperature red giants that were previously difficult to constrain are now accounted for, and the best fit favors a two-component model. The dominant component contains 93%$\pm$3% of the mass, is metal-rich ($\overline{[M/H]}\sim$0.45), and has an age of 5$^{+3}_{-2}$ Gyr, which is $\sim$3 Gyr younger than earlier studies with fixed (solar) metallicity; this younger age challenges co-evolutionary models in which the NSC and supermassive black holes formed simultaneously at early times. The minor population component has low metallicity ($\overline{[M/H]}\sim$ -1.1) and contains $\sim$7% of the stellar mass. The age of the minor component is uncertain (0.1 - 5 Gyr old). Using the estimated parameters, we infer the following NSC stellar remnant population (with $\sim$18% uncertainty): 1.5$\times$10$^5$ neutron stars, 2.5$\times$10$^5$ stellar mass black holes (BHs) and 2.2$\times$10$^4$ BH-BH binaries. These predictions result in 2-4 times fewer neutron stars compared to earlier predictions that assume solar metallicity, introducing a possible new path to understand the so-called "missing pulsar problem". Finally, we present updated predictions for the BH-BH merger rates (0.01-3 Gpc$^{-3}$yr$^{-1}$). △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: 37 pages, 23 figures, Accepted to ApJ

Primordial Black Hole Dark Matter Simulations Using PopSyCLE

Authors: K. Pruett, W. Dawson, M. S. Medford, C. Lam, J. R. Lu, S. Perkins, P. McGill, N. Golovich, G. Chapline

Abstract: Primordial black holes (PBHs), theorized to have originated in the early universe, are speculated to be a viable form of dark matter. If they exist, they should be detectable through photometric and astrometric signals resulting from gravitational microlensing of stars in the Milky Way. Population Synthesis for Compact-object Lensing Events, or PopSyCLE, is a simulation code that enables users to… ▽ More Primordial black holes (PBHs), theorized to have originated in the early universe, are speculated to be a viable form of dark matter. If they exist, they should be detectable through photometric and astrometric signals resulting from gravitational microlensing of stars in the Milky Way. Population Synthesis for Compact-object Lensing Events, or PopSyCLE, is a simulation code that enables users to simulate microlensing surveys, and is the first of its kind to include both photometric and astrometric microlensing effects, which are important for potential PBH detection and characterization. To estimate the number of observable PBH microlensing events we modify PopSyCLE to include a dark matter halo consisting of PBHs. We detail our PBH population model, and demonstrate our PopSyCLE + PBH results through simulations of the OGLE-IV and Roman microlensing surveys. We provide a proof-of-concept analysis for adding PBHs into PopSyCLE, and thus include many simplifying assumptions, such as $f_{\text{DM}}$, the fraction of dark matter composed of PBHs, and $\bar{m}_{\text{PBH}}$, mean PBH mass. Assuming $\bar{m}_{\text{PBH}}=30$ $M_{\odot}$, we find $\sim$ 3.6$f_{\text{DM}}$ times as many PBH microlensing events than stellar evolved black hole events, a PBH average peak Einstein crossing time of $\sim$ 91.5 days, estimate on order of $10^2f_{\text{DM}}$ PBH events within the 8 year OGLE-IV results, and estimate Roman to detect $\sim$ 1,000$f_{\text{DM}}$ PBH microlensing events throughout its planned microlensing survey. △ Less

Submitted 8 April, 2024; v1 submitted 12 November, 2022; originally announced November 2022.

The Impact of Initial-Final Mass Relations on Black Hole Microlensing

Authors: Sam Rose, Casey Y. Lam, Jessica R. Lu, Michael Medford, Matthew W. Hosek Jr., Natasha S. Abrams, Emily Ramey, Sergiy S. Vasylyev

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 non-luminous remnant objects (i.e. neutron stars and black holes). Gravitational wave detectors have opened the possibility of finding large numbers of compact objec… ▽ More 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 non-luminous 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 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 \lesssim 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. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: 24 pages, 17 figures Accepted to ApJ

A Search for Predicted Astrometric Microlensing Events by Nearby Brown Dwarfs

Authors: Judah Luberto, Emily C. Martin, Peter McGill, Alexie Leauthaud, Andrew J. Skemer, Jessica R. Lu

Abstract: Gravitational microlensing has the potential to provide direct gravitational masses of single, free-floating brown dwarfs, independent of evolutionary and atmospheric models. The proper motions and parallaxes of nearby brown dwarfs can be used to predict close future alignments with distant background stars that cause a microlensing event. Targeted astrometric follow up of the predicted microlensi… ▽ More Gravitational microlensing has the potential to provide direct gravitational masses of single, free-floating brown dwarfs, independent of evolutionary and atmospheric models. The proper motions and parallaxes of nearby brown dwarfs can be used to predict close future alignments with distant background stars that cause a microlensing event. Targeted astrometric follow up of the predicted microlensing events permits the brown dwarf's mass to be measured. Predicted microlensing events are typically found via searching for a peak threshold signal using an estimate of the lens mass. We develop a novel method that finds predicted events that instead will lead to a target lens mass precision. The main advantage of our method is that it does not require a lens mass estimate. We use this method to search for predicted astrometric microlensing events occurring between 2014 - 2032 using a catalog of 1225 low mass star and brown dwarf lenses in the Solar Neighborhood of spectral type M6 or later and a background source catalog from DECaLS Data Release 9. The background source catalog extends to $g = $ 23.95, providing a more dense catalog compared to Gaia. Our search did not reveal any upcoming microlensing events. We estimate the rate of astrometric microlensing event for brown dwarfs in the Legacy Survey and find it to be low $\sim10^{-5}$yr$^{-1}$. We recommend carrying out targeted searches for brown dwarfs in front of the Galactic Bulge and Plane to find astrometric microlensing events that will allow the masses of single, free-floating brown dwarfs to be measured. △ Less

Submitted 25 October, 2022; originally announced October 2022.

Comments: 12 pages, 6 figures, Accepted to AJ

AIROPA II: Modeling Instrumental Aberrations for Off-Axis Point Spread Functions in Adaptive Optics

Authors: Anna Ciurlo, Paolo Turri, Gunther Witzel, Jessica R. Lu, Tuan Do, Breann N. Sitarski, Michael P. Fitzgerald, Andrea M. Ghez, Carlos Alvarez, Sean K. Terry, Greg Doppmann, James E. Lyke, Sam Ragland, Randall Campbell, Keith Matthews

Abstract: Images obtained with single-conjugate adaptive optics (AO) show spatial variation of the point spread function (PSF) due to both atmospheric anisoplanatism and instrumental aberrations. The poor knowledge of the PSF across the field of view strongly impacts the ability to take full advantage of AO capabilities. The AIROPA project aims to model these PSF variations for the NIRC2 imager at the Keck… ▽ More Images obtained with single-conjugate adaptive optics (AO) show spatial variation of the point spread function (PSF) due to both atmospheric anisoplanatism and instrumental aberrations. The poor knowledge of the PSF across the field of view strongly impacts the ability to take full advantage of AO capabilities. The AIROPA project aims to model these PSF variations for the NIRC2 imager at the Keck Observatory. Here, we present the characterization of the instrumental phase aberrations over the entire NIRC2 field of view and we present a new metric for quantifying the quality of the calibration, the fraction of variance unexplained (FVU). We used phase diversity measurements obtained on an artificial light source to characterize the variation of the aberrations across the field of view and their evolution with time. We find that there is a daily variation of the wavefront error (RMS of the residuals is 94~nm) common to the whole detector, but the differential aberrations across the field of view are very stable (RMS of the residuals between different epochs is 59~nm). This means that instrumental calibrations need to be monitored often only at the center of the detector, and the much more time-consuming variations across the field of view can be characterized less frequently (most likely when hardware upgrades happen). Furthermore, we tested AIROPA's instrumental model through real data of the fiber images on the detector. We find that modeling the PSF variations across the field of view improves the FVU metric by 60\% and reduces the detection of fake sources by 70\%. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Journal ref: Journal of Astronomical Telescopes, Instruments, and Systems, Vol. 8, Issue 3, 038007 (September 2022)

AIROPA IV: Validating Point Spread Function Reconstruction on Various Science Cases

Authors: Sean K. Terry, Jessica R. Lu, Paolo Turri, Anna Ciurlo, Abhimat Gautam, Tuan Do, Michael P. Fitzgerald, Andrea Ghez, Matthew Hosek Jr., Gunther Witzel

Abstract: We present an analysis of six independent on-sky datasets taken with the Keck-II/NIRC2 instrument. Using the off-axis point spread function (PSF) reconstruction software AIROPA, we extract stellar astrometry, photometry, and other fitting metrics in order to characterize the performance of this package. We test the effectiveness of AIROPA to reconstruct the PSF across the field of view in varying… ▽ More We present an analysis of six independent on-sky datasets taken with the Keck-II/NIRC2 instrument. Using the off-axis point spread function (PSF) reconstruction software AIROPA, we extract stellar astrometry, photometry, and other fitting metrics in order to characterize the performance of this package. We test the effectiveness of AIROPA to reconstruct the PSF across the field of view in varying atmospheric conditions, number and location of PSF reference stars, stellar crowding and telescope position angle (PA). We compare the astrometric precision and fitting residuals between a static PSF model and a spatially varying PSF model that incorporates instrumental aberrations and atmospheric turbulence during exposures. Most of the fitting residuals we measure show little to no improvement in the variable-PSF mode over the single-PSF mode. For one of the data sets, we find photometric performance is significantly improved (by ${\sim}10\times$) by measuring the trend seen in photometry as a function of off-axis location. For nearly all other metrics we find comparable astrometric and photometric precision across both PSF modes, with a ${\sim}13$% smaller astrometric uncertainty in variable-PSF mode in the best case. We largely confirm that the spatially variable PSF does not significantly improve the astrometric and other PSF fitting residuals over the static PSF for on-sky observations. We attribute this to unaccounted instrumental aberrations that are not characterized through afternoon adaptive optics (AO) bench calibrations. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: 23 pages, 13 figures, submitted to JATIS

Measuring the Orbits of the Arches and Quintuplet Clusters using HST and Gaia: Exploring Scenarios for Star Formation Near the Galactic Center

Authors: Matthew W. Hosek Jr., Tuan Do, Jessica R. Lu, Mark R. Morris, Andrea M. Ghez, Gregory D. Martinez, Jay Anderson

Abstract: We present new absolute proper motion measurements for the Arches and Quintuplet clusters, two young massive star clusters near the Galactic center. Using multi-epoch HST observations, we construct proper motion catalogs for the Arches ($\sim$35,000 stars) and Quintuplet ($\sim$40,000 stars) fields in ICRF coordinates established using stars in common with the Gaia EDR3 catalog. The bulk proper mo… ▽ More We present new absolute proper motion measurements for the Arches and Quintuplet clusters, two young massive star clusters near the Galactic center. Using multi-epoch HST observations, we construct proper motion catalogs for the Arches ($\sim$35,000 stars) and Quintuplet ($\sim$40,000 stars) fields in ICRF coordinates established using stars in common with the Gaia EDR3 catalog. The bulk proper motions of the clusters are measured to be ($μ_{α*}$, $μ_δ$) = (-0.80 $\pm$ 0.032, -1.89 $\pm$ 0.021) mas/yr for the Arches and ($μ_{α*}$, $μ_δ$) = (-0.96 $\pm$ 0.032, -2.29 $\pm$ 0.023) mas/yr for the Quintuplet, achieving $\sim$5x higher precision than past measurements. We place the first constraints on the properties of the cluster orbits that incorporate the uncertainty in their current line-of-sight distances. The clusters will not approach closer than $\sim$25 pc to SgrA*, making it unlikely that they will inspiral into the Nuclear Star Cluster within their lifetime. Further, the cluster orbits are not consistent with being circular; the average value of r$_{apo}$ / r$_{peri}$ is $\sim$1.9 (equivalent to eccentricity of $\sim$0.31) for both clusters. Lastly, we find that the clusters do not share a common orbit, challenging one proposed formation scenario in which the clusters formed from molecular clouds on the open stream orbit derived by Kruijssen et al. (2015). Meanwhile, our constraints on the birth location and velocity of the clusters offer mild support for a scenario in which the clusters formed via collisions between gas clouds on the x1 and x2 bar orbit families. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Comments: Accepted for publication in ApJ. 38 pages, 25 figures. Proper motion catalogs included in ancillary materials

AIROPA III: Testing Simulated and On-Sky Data

Authors: Paolo Turri, Jessica R. Lu, Gunther Witzel, Anna Ciurlo, Tuan Do, Andrea M. Ghez, Matthew C. Britton, Sam Ragland, Sean K. Terry

Abstract: Adaptive optics images from the W. M. Keck Observatory have delivered numerous influential scientific results, including detection of multi-system asteroids, the supermassive black hole at the center of the Milky Way, and directly imaged exoplanets. Specifically, the precise and accurate astrometry these images yield was used to measure the mass of the supermassive black hole using orbits of the s… ▽ More Adaptive optics images from the W. M. Keck Observatory have delivered numerous influential scientific results, including detection of multi-system asteroids, the supermassive black hole at the center of the Milky Way, and directly imaged exoplanets. Specifically, the precise and accurate astrometry these images yield was used to measure the mass of the supermassive black hole using orbits of the surrounding star cluster. Despite these successes, one of the major obstacles to improved astrometric measurements is the spatial and temporal variability of the point-spread function delivered by the instruments. AIROPA is a software package for the astrometric and photometric analysis of adaptive optics images using point-spread function fitting together with the technique of point-spread function reconstruction. In adaptive optics point-spread function reconstruction, the knowledge of the instrument performance and of the atmospheric turbulence is used to predict the long-exposure point-spread function of an observation. In this paper we present the results of our tests using AIROPA on both simulated and on-sky images of the Galactic Center. We find that our method is very reliable in accounting for the static aberrations internal to the instrument, but it does not improve significantly the accuracy on sky, possibly due to uncalibrated telescope aberrations. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Searching for a Hypervelocity White Dwarf Companion: A Proper Motion Survey of SN 1006

Authors: Joshua V. Shields, Wolfgang Kerzendorf, Matthew W. Hosek Jr., Ken J. Shen, Armin Rest, Tuan Do, Jessica R. Lu, Andrew G. Fullard, Giovanni Strampelli, Alfredo Zenteno

Abstract: Type Ia Supernovae (SNe Ia) are securely understood to come from the thermonuclear explosion of a white dwarf as a result of binary interaction, but the nature of that binary interaction and the secondary object is uncertain. Recently, a double white dwarf model known as the dynamically driven double-degenerate double-detonation (D6) model has become a promising explanation for these events. One r… ▽ More Type Ia Supernovae (SNe Ia) are securely understood to come from the thermonuclear explosion of a white dwarf as a result of binary interaction, but the nature of that binary interaction and the secondary object is uncertain. Recently, a double white dwarf model known as the dynamically driven double-degenerate double-detonation (D6) model has become a promising explanation for these events. One realization of this scenario predicts that the companion may survive the explosion and reside within the remnant as a fast moving ($V_{peculiar} >1000$ km s$^{-1}$), overluminous ($L > 0.1 L_\odot$) white dwarf. Recently, three objects which appear to have these unusual properties have been discovered in the Gaia survey. We obtained photometric observations of the SN Ia remnant SN 1006 with the Dark Energy Camera over four years to attempt to discover a similar star. We present a deep, high precision astrometric proper motion survey of the interior stellar population of the remnant. We rule out the existence of a high proper motion object consistent with our tested realization of the D6 scenario ($V_{transverse} > 600$ km s$^{-1}$ with $m_r < 21$ corresponding to an intrinsic luminosity of $L > 0.0176 L_\odot$). We conclude that such a star does not exist within the remnant, or is hidden from detection by either strong localized dust or the unlikely possibility of ejection from the binary system near parallel to the line of sight. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: 15 pages, 10 figures

Journal ref: https://iopscience.iop.org/article/10.3847/2041-8213/ac7950

Adaptive Optics Imaging Breaks the Central Caustic Cusp Approach Degeneracy in High Magnification Microlensing Events

Authors: Sean K. Terry, David P. Bennett, Aparna Bhattacharya, Naoki Koshimoto, Jean-Phillipe Beaulieu, Joshua W. Blackman, Ian A. Bond, Andrew A. Cole, Jessica R. Lu, Jean Baptiste Marquette, Clément Ranc, Natalia Rektsini, Aikaterini Vandorou

Abstract: We report new results for the gravitational microlensing target OGLE-2011-BLG-0950 from adaptive optics (AO) images using the Keck observatory. The original analysis by Choi et al. 2012 reports degenerate solutions between planetary and stellar binary lens systems. This is due to a degeneracy in high magnification events where the shape of the light curve peak can be explained by a source approach… ▽ More We report new results for the gravitational microlensing target OGLE-2011-BLG-0950 from adaptive optics (AO) images using the Keck observatory. The original analysis by Choi et al. 2012 reports degenerate solutions between planetary and stellar binary lens systems. This is due to a degeneracy in high magnification events where the shape of the light curve peak can be explained by a source approach to two different cusp geometries with different source radius crossing times. This particular case is the most important type of degeneracy for exoplanet demographics, because the distinction between a planetary mass or stellar binary companion has direct consequences for microlensing exoplanet statistics. The 8 and 10-year baselines between the event and the Keck observations allow us to directly measure a relative proper motion of $4.20\pm 0.21\,$mas/yr, which confirms the detection of the lens star system and directly rules out the planetary companion models that predict a ${\sim}4 \times$ smaller relative proper motion. The combination of the lens brightness and close stellar binary light curve parameters yield primary and secondary star masses of $M_{A} = 1.12^{+0.06}_{-0.04}M_\odot$ and $M_{B} = 0.47^{+0.04}_{-0.03}M_\odot$ at a distance of $D_L = 6.70^{+0.55}_{-0.30}\,$kpc, and a primary-secondary projected separation of $0.39^{+0.05}_{-0.04}\,$AU. Since this degeneracy is likely to be common, the high resolution imaging method described here will be used to disentangle the central caustic cusp approach degeneracy for events observed by the \textit{Roman} exoplanet microlensing survey using the \textit{Roman} images taken near the beginning or end of the survey. △ Less

Submitted 1 November, 2022; v1 submitted 7 June, 2022; originally announced June 2022.

Comments: Revised version, 19 pages, 8 figures. AJ, 164, 217

An isolated mass gap black hole or neutron star detected with astrometric microlensing

Authors: Casey Y. Lam, Jessica R. Lu, Andrzej Udalski, Ian Bond, David P. Bennett, Jan Skowron, Przemek Mroz, Radek Poleski, Takahiro Sumi, Michal K. Szymanski, Szymon Kozlowski, Pawel Pietrukowicz, Igor Soszynski, Krzysztof Ulaczyk, Lukasz Wyrzykowski, Shota Miyazaki, Daisuke Suzuki, Naoki Koshimoto, Nicholas J. Rattenbury, Matthew W. Hosek Jr., Fumio Abe, Richard Barry, Aparna Bhattacharya, Akihiko Fukui, Hirosane Fujii , et al. (20 additional authors not shown)

Abstract: We present the analysis of five black hole candidates identified from gravitational microlensing surveys. Hubble Space Telescope astrometric data and densely sampled lightcurves 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 (OGL… ▽ More We present the analysis of five black hole candidates identified from gravitational microlensing surveys. Hubble Space Telescope astrometric data and densely sampled lightcurves 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.4 $M_\odot$. 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. OB110462 shows significant tension between models well-fit to photometry vs. 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 $<2 M_\odot$ 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 ($\sim 10^8$) and find reasonable agreement given the small sample size. △ Less

Submitted 31 May, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: Accepted to ApJ Letters, with corresponding ApJ Supplement. 10 page Letter (6 figures, 2 tables) + 51 page Supplement (27 figures, 20 tables, 9 appendices). Some minor updates from the refereeing process, but no change to main conclusions

60 Microlensing Events from the Three Years of Zwicky Transient Facility Phase One

Authors: Michael S. Medford, Natasha S. Abrams, Jessica R. Lu, Peter Nugent, Casey Y. Lam

Abstract: Microlensing events have historically been discovered throughout the Galactic bulge and plane by surveys designed solely for that purpose. We conduct the first multi-year search for microlensing events on the Zwicky Transient Facility (ZTF), an all-sky optical synoptic survey that observes the entire visible Northern sky every few nights. We discover 60 high quality microlensing events in the thre… ▽ More Microlensing events have historically been discovered throughout the Galactic bulge and plane by surveys designed solely for that purpose. We conduct the first multi-year search for microlensing events on the Zwicky Transient Facility (ZTF), an all-sky optical synoptic survey that observes the entire visible Northern sky every few nights. We discover 60 high quality microlensing events in the three years of ZTF-I using the bulk lightcurves in the ZTF Public Data Release 5. 19 of our events are found outside of the Galactic plane ($|b| \geq 15^\circ$), nearly doubling the number of previously discovered events in the stellar halo from surveys pointed toward the Magellanic Clouds and the Andromeda Galaxy. We also record 1,558 ongoing candidate events as potential microlensing that can continue to be observed by ZTF-II for identification. The scalable and computationally efficient methods developed in this work can be applied to future synoptic surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope, as they attempt to find microlensing events in even larger and deeper datasets. △ Less

Submitted 8 February, 2023; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: 33 pages, 15 figures, 7 tables, Accepted to ApJ, microlensing catalogs available at https://zenodo.org/record/7277537

Developing Adaptive Secondary Mirror Concepts for the APF and W.M. Keck Observatory Based on HVR Technology

Authors: Philip M. Hinz, Rachel Bowens-Rubin, Christoph Baranec, Kevin Bundy, Mark Chun, Daren Dillon, Brad Holden, Wouter Jonker, Molly Kosiarek, Renate Kupke, Stefan Kuiper, Olivier Lai, Jessica R. Lu, Matthew Maniscalco, Matthew Radovan, Sam Ragland, Stephanie Sallum, Andrew Skemer, Peter Wizinowich

Abstract: An Adaptive secondary mirror (ASM) allows for the integration of adaptive optics (AO) into the telescope itself. Adaptive secondary mirrors, based on hybrid variable reluctance (HVR) actuator technology, developed by TNO, provide a promising path to telescope-integrated AO. HVR actuators have the advantage of allowing mirrors that are stiffer, more power efficient, and potentially less complex tha… ▽ More An Adaptive secondary mirror (ASM) allows for the integration of adaptive optics (AO) into the telescope itself. Adaptive secondary mirrors, based on hybrid variable reluctance (HVR) actuator technology, developed by TNO, provide a promising path to telescope-integrated AO. HVR actuators have the advantage of allowing mirrors that are stiffer, more power efficient, and potentially less complex than similar, voice-coil based ASM's. We are exploring the application of this technology via a laboratory testbed that will validate the technical approach. In parallel, we are developing conceptual designs for ASMs at several telescopes including the Automated Planet Finder Telescope (APF) and for Keck Observatory. An ASM for APF has the potential to double the light through the slit for radial velocity measurements, and dramatically improved the image stability. An ASM for WMKO enables ground layer AO correction and lower background infrared AO observations, and provides for more flexible deployment of instruments via the ability to adjust the location of the Cassegrain focus. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: 16 pages, Proceedings of SPIE

Trigonometric parallaxes of two T dwarfs with Keck and ShaneAO astrometry

Authors: Jayke Nguyen, S. Mark Ammons, Kaitlin Dennison, E. Victor Garcia, Jessica R. Lu, Stephen McMillan, Maissa Salama

Abstract: We present trigonometric parallax and proper motion measurements for two T-type brown dwarfs. We derive our measurements from infrared laser guide star adaptive optics observations spanning five years from the ShaneAO/SHARCS and NIRC2/medium-cam instruments on the Shane and Keck telescopes, respectively. To improve our astrometric precision, we measure and apply a distortion correction to our fiel… ▽ More We present trigonometric parallax and proper motion measurements for two T-type brown dwarfs. We derive our measurements from infrared laser guide star adaptive optics observations spanning five years from the ShaneAO/SHARCS and NIRC2/medium-cam instruments on the Shane and Keck telescopes, respectively. To improve our astrometric precision, we measure and apply a distortion correction to our fields for both instruments. We also transform the Keck and ShaneAO astrometric reference frames onto the ICRS using five-parameter parallax and proper motion solutions for background reference stars from Gaia DR2. Fitting for parallax and proper motion, we measure parallaxes of $73.5\pm9.2$ mas and $70.1\pm6.7$ mas for WISEJ19010703+47181688 (WISE1901) and WISEJ21543294+59421370 (WISE2154), respectively. We utilize Monte Carlo methods to estimate the error in our sparse field methods, taking into account overfitting and differential atmospheric refraction. Comparing to previous measurements in the literature, all of our parallax and proper motion values fall within $2σ$ of the published measurements, and 4 of 6 measurements are within $1σ$. These data are among the first parallax measurements of these T dwarfs and serve as precise measurements for calibrating stellar formation models. These two objects are the first results of an ongoing survey of T dwarfs with Keck/NIRC2 and the Shane Adaptive Optics system at Lick Observatory. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: 18 Pages, 19 Figures. Accepted by PASP

Report number: LLNL-JRNL-821003

The 3-D Kinematics of the Orion Nebula Cluster: NIRSPEC-AO Radial Velocities of the Core Population

Authors: Christopher A. Theissen, Quinn M. Konopacky, Jessica R. Lu, Dongwon Kim, Stella Y. Zhang, Chih-Chun Hsu, Laurie Chu, Lingfeng Wei

Abstract: The kinematics and dynamics of stellar and substellar populations within young, still-forming clusters provides valuable information for constraining theories of formation mechanisms. Using Keck II NIRSPEC+AO data, we have measured radial velocities for 56 low-mass sources within 4' of the core of the Orion Nebula Cluster (ONC). We also re-measure radial velocities for 172 sources observed with SD… ▽ More The kinematics and dynamics of stellar and substellar populations within young, still-forming clusters provides valuable information for constraining theories of formation mechanisms. Using Keck II NIRSPEC+AO data, we have measured radial velocities for 56 low-mass sources within 4' of the core of the Orion Nebula Cluster (ONC). We also re-measure radial velocities for 172 sources observed with SDSS/APOGEE. These data are combined with proper motions measured using $HST$ ACS/WFPC2/WFC3IR and Keck II NIRC2, creating a sample of 135 sources with all three velocity components. The velocities measured are consistent with a normal distribution in all three components. We measure intrinsic velocity dispersions of ($σ_{v_α}$, $σ_{v_δ}$, $σ_{v_r}$) = ($1.64\pm0.12$, $2.03\pm0.13$, $2.56^{+0.16}_{-0.17}$) km s$^{-1}$. Our computed intrinsic velocity dispersion profiles are consistent with the dynamical equilibrium models from Da Rio et al. (2014) in the tangential direction, but not in the line of sight direction, possibly indicating that the core of the ONC is not yet virialized, and may require a non-spherical potential to explain the observed velocity dispersion profiles. We also observe a slight elongation along the north-south direction following the filament, which has been well studied in previous literature, and an elongation in the line of sight to tangential velocity direction. These 3-D kinematics will help in the development of realistic models of the formation and early evolution of massive clusters. △ Less

Submitted 20 October, 2021; v1 submitted 12 May, 2021; originally announced May 2021.

Comments: Accepted to ApJ; 33 pages, 5 tables; 15 figures

On the Possibility of Stellar Lenses in the Black Hole Candidate Microlensing Events MACHO-96-BLG-5 and MACHO-98-BLG-6

Authors: Fatima N. Abdurrahman, Haynes F. Stephens, Jessica R. Lu

Abstract: Though stellar-mass black holes (BHs) are likely abundant in the Milky Way (N=10^8-10^9), only ~20 have been detected to date, all in accreting binary systems (Casares 2006). Gravitational microlensing is a proposed technique to search for isolated BHs, which to date have not been detected. Two microlensing events, MACHO-1996-BLG-5 (M96-B5) and MACHO-1998-BLG-6 (M98-B6), initially observed near th… ▽ More Though stellar-mass black holes (BHs) are likely abundant in the Milky Way (N=10^8-10^9), only ~20 have been detected to date, all in accreting binary systems (Casares 2006). Gravitational microlensing is a proposed technique to search for isolated BHs, which to date have not been detected. Two microlensing events, MACHO-1996-BLG-5 (M96-B5) and MACHO-1998-BLG-6 (M98-B6), initially observed near the lens-source minimum angular separation in 1996 and 1998, respectively, have long Einstein crossing times (>300 days), identifying the lenses as candidate black holes. Twenty years have elapsed since the time of lens-source closest approach for each of these events, indicating that if the lens and source are both luminous, and if their relative proper motion is sufficiently large, the two components should be spatially resolvable. We attempt to eliminate the possibility of a stellar lens for these events by: (1) using Keck near-infrared adaptive optics images to search for a potentially now-resolved, luminous lens; and (2) examining multi-band photometry of the source to search for flux contributions from a potentially unresolved, luminous lens. We combine detection limits from NIRC2 images with light curve data to eliminate all non-BH lenses for relative lens-source proper motions above 0.81 mas/yr for M96-B5 and 2.48 mas/yr for M98-B6. Further, we use WFPC2 broadband images to eliminate the possibility of stellar lenses at any proper motion. We present the narrow range of non-BH possibilities allowed by our varied analyses. Finally, we suggest future observations that would constrain the remaining parameter space with the methods developed in this work. △ Less

Submitted 17 March, 2021; originally announced March 2021.

Comments: Accepted to ApJ on March 13, 2021

Real-Time Likelihood-Free Inference of Roman Binary Microlensing Events with Amortized Neural Posterior Estimation

Authors: Keming Zhang, Joshua S. Bloom, B. Scott Gaudi, Francois Lanusse, Casey Lam, Jessica R. Lu

Abstract: Fast and automated inference of binary-lens, single-source (2L1S) microlensing events with sampling-based Bayesian algorithms (e.g., Markov Chain Monte Carlo; MCMC) is challenged on two fronts: high computational cost of likelihood evaluations with microlensing simulation codes, and a pathological parameter space where the negative-log-likelihood surface can contain a multitude of local minima tha… ▽ More Fast and automated inference of binary-lens, single-source (2L1S) microlensing events with sampling-based Bayesian algorithms (e.g., Markov Chain Monte Carlo; MCMC) is challenged on two fronts: high computational cost of likelihood evaluations with microlensing simulation codes, and a pathological parameter space where the negative-log-likelihood surface can contain a multitude of local minima that are narrow and deep. Analysis of 2L1S events usually involves grid searches over some parameters to locate approximate solutions as a prerequisite to posterior sampling, an expensive process that often requires human-in-the-loop domain expertise. As the next-generation, space-based microlensing survey with the Roman Space Telescope is expected to yield thousands of binary microlensing events, a new fast and automated method is desirable. Here, we present a likelihood-free inference (LFI) approach named amortized neural posterior estimation, where a neural density estimator (NDE) learns a surrogate posterior $\hat{p}(θ|x)$ as an observation-parametrized conditional probability distribution, from pre-computed simulations over the full prior space. Trained on 291,012 simulated Roman-like 2L1S simulations, the NDE produces accurate and precise posteriors within seconds for any observation within the prior support without requiring a domain expert in the loop, thus allowing for real-time and automated inference. We show that the NDE also captures expected posterior degeneracies. The NDE posterior could then be refined into the exact posterior with a downstream MCMC sampler with minimal burn-in steps. △ Less

Submitted 30 March, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 15 pages, 8 figures, 3 tables. Accepted for publication in the Astronomical Journal. This article supersedes arXiv:2010.04156

MOA-2009-BLG-319Lb: A Sub-Saturn Planet Inside the Predicted Mass Desert

Authors: Sean K. Terry, Aparna Bhattacharya, David P. Bennett, Jean-Phillipe Beaulieu, Naoki Koshimoto, Joshua W. Blackman, Ian A. Bond, Andrew A. Cole, Calen B. Henderson, Jessica R. Lu, Jean Baptiste Marquette, Clement Ranc, Aikaterini Vandorou

Abstract: We present an adaptive optics (AO) analysis of images from the Keck-II telescope NIRC2 instrument of the planetary microlensing event MOA-2009-BLG-319. The $\sim$10 year baseline between the event and the Keck observations allows the planetary host star to be detected at a separation of $66.5\pm 1.7\,$mas from the source star, consistent with the light curve model prediction. The combination of th… ▽ More We present an adaptive optics (AO) analysis of images from the Keck-II telescope NIRC2 instrument of the planetary microlensing event MOA-2009-BLG-319. The $\sim$10 year baseline between the event and the Keck observations allows the planetary host star to be detected at a separation of $66.5\pm 1.7\,$mas from the source star, consistent with the light curve model prediction. The combination of the host star brightness and light curve parameters yield host star and planet masses of M_host = 0.514 $\pm$ 0.063M_Sun and m_p = 66.0 $\pm$ 8.1M_Earth at a distance of $D_L = 7.0 \pm 0.7\,$kpc. The star-planet projected separation is $2.03 \pm 0.21\,$AU. The planet-star mass ratio of this system, $q = (3.857 \pm 0.029)\times 10^{-4}$, places it in the predicted "planet desert" at $10^{-4} < q < 4\times 10^{-4}$ according to the runaway gas accretion scenario of the core accretion theory. Seven of the 30 planets in the Suzuki et al. (2016) sample fall in this mass ratio range, and this is the third with a measured host mass. All three of these host stars have masses of 0.5 $\leq$ M_host/M_Sun $\leq$ 0.7, which implies that this predicted mass ratio gap is filled with planets that have host stars within a factor of two of 1M_Sun. This suggests that runaway gas accretion does not play a major role in determining giant planet masses for stars somewhat less massive than the Sun. Our analysis has been accomplished with a modified DAOPHOT code that has been designed to measure the brightness and positions of closely blended stars. This will aid in the development of the primary method that the Nancy Grace Roman Space Telescope mission will use to determine the masses of microlens planets and their hosts. △ Less

Submitted 6 January, 2021; v1 submitted 17 September, 2020; originally announced September 2020.

Comments: 16 pages, 6 figures, AJ in press

Journal ref: AJ 161 54 (2021)

A Reanalysis of Public Galactic Bulge Gravitational Microlensing Events from OGLE-III and IV

Authors: Nathan Golovich, William A. Dawson, Fran Bartolić, Casey Y. Lam, Jessica R. Lu, Michael S. Medford, Michael D. Schneider, George Chapline, Edward F. Schlafly, Alex Drlica-Wagner, Kerianne Pruett

Abstract: Modern surveys of gravitational microlensing events have progressed to detecting thousands per year. 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 studying the microlensing Einstein radius crossing time distribution ($t_E$). However, systematics in individual light… ▽ More Modern surveys of gravitational microlensing events have progressed to detecting thousands per year. 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 studying the microlensing Einstein radius crossing time distribution ($t_E$). However, systematics in individual light curves as well as over-simplistic 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 timescale 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 previous analyses over-estimated the number of long timescale ($t_E>100$ days) events due to their over simplistic models ignoring parallax effects and stellar variability. We use our fits to identify hundreds of events that are likely black holes. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: Submitted version, in review, 33 pages, 18 figures, MCMC posterior samples available by publisher after acceptance

MOA-2007-BLG-400 A Super-Jupiter Mass Planet Orbiting a Galactic BulgeK-dwarf Revealed by Keck Adaptive Optics Imaging

Authors: Aparna Bhattacharya, David P. Bennett, Jean P. Beaulieu, Ian A. Bond, Naoki Koshimoto, Jessica R. Lu, Joshua W. Blackman, Aikaterini Vandorou, Sean K. Terry, Virginie Batista, Jean B. Marquette, Andrew A. Cole, Akihiko Fukui, Calen B. Henderson

Abstract: We present Keck/NIRC2 adaptive optics imaging of planetary microlensing event MOA-2007-BLG-400 that resolves the lens star system from the source. We find that the MOA-2007-BLG-400L planetary system consists of a $1.71\pm 0.27 M_{\rm Jup}$ planet orbiting a $0.69\pm 0.04M_{\odot}$ K-dwarf host star at a distance of $6.89\pm 0.77\,$kpc from the Sun. So, this planetary system probably resides in the… ▽ More We present Keck/NIRC2 adaptive optics imaging of planetary microlensing event MOA-2007-BLG-400 that resolves the lens star system from the source. We find that the MOA-2007-BLG-400L planetary system consists of a $1.71\pm 0.27 M_{\rm Jup}$ planet orbiting a $0.69\pm 0.04M_{\odot}$ K-dwarf host star at a distance of $6.89\pm 0.77\,$kpc from the Sun. So, this planetary system probably resides in the Galactic bulge. The planet-host star projected separation is only weakly constrained due to the close-wide light curve degeneracy; the 2$σ$ projected separation range is 0.6--$7.2\,$AU. This host mass is at the top end of the range of masses predicted by a standard Bayesian analysis that assumes that all stars have an equal chance of hosting a star of the observed mass ratio. This and the similar result for event MOA-2013-BLG-220 suggests that more massive stars may be more likely to host planets with a mass ratio in the $0.002 < q < 0.004$ range that orbit beyond the snow line. These results also indicate the importance of host star mass measurements for exoplanets found by microlensing. The microlensing survey imaging data from NASA's Nancy Grace Roman Space Telescope (formerly WFIRST) mission will be doing mass measurements like this for a huge number of planetary events. This host lens is the highest contrast lens-source detected in microlensing mass measurement analysis (the lens being 10$\times$ fainter than the source). We present an improved method of calculating photometry and astrometry uncertainties based on the Jackknife method, which produces more accurate errors that are $\sim$$2.5 \times$ larger than previous estimates. △ Less

Submitted 4 September, 2020; originally announced September 2020.

Comments: 29 pages, 4 figures, AJ submitted

SPISEA: A Python-Based Simple Stellar Population Synthesis Code for Star Clusters

Authors: M. W. Hosek Jr, J. R. Lu, C. Y. Lam, A. K. Gautam, K. E. Lockhart, D. Kim, S. Jia

Abstract: We present SPISEA (Stellar Population Interface for Stellar Evolution and Atmospheres), an open-source Python package that simulates simple stellar populations. The strength of SPISEA is its modular interface which offers the user control of 13 input properties including (but not limited to) the Initial Mass Function, stellar multiplicity, extinction law, and the metallicity-dependent stellar evol… ▽ More We present SPISEA (Stellar Population Interface for Stellar Evolution and Atmospheres), an open-source Python package that simulates simple stellar populations. The strength of SPISEA is its modular interface which offers the user control of 13 input properties including (but not limited to) the Initial Mass Function, stellar multiplicity, extinction law, and the metallicity-dependent stellar evolution and atmosphere model grids used. The user also has control over the Initial-Final Mass Relation in order to produce compact stellar remnants (black holes, neutron stars, and white dwarfs). We demonstrate several outputs produced by the code, including color-magnitude diagrams, HR-diagrams, luminosity functions, and mass functions. SPISEA is object-oriented and extensible, and we welcome contributions from the community. The code and documentation are available on GitHub and ReadtheDocs, respectively. △ Less

Submitted 10 July, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 18 pages, 7 figures. Version update: Software package renamed from "PyPopStar" to "SPISEA" due to a naming conflict. Also added a few references to intro. Otherwise, paper is the same

Gravitational Microlensing Event Statistics for the Zwicky Transient Facility

Authors: Michael S. Medford, Jessica R. Lu, William A. Dawson, Casey Y. Lam, Nathan R. Golovich, Edward F. Schlafly, Peter Nugent

Abstract: Microlensing surveys have discovered thousands of events with almost all events discovered within the Galactic bulge or toward the Magellanic clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky every few nights including the Galactic plane. ZTF observes $\sim10^9$ stars in g-band and r… ▽ More Microlensing surveys have discovered thousands of events with almost all events discovered within the Galactic bulge or toward the Magellanic clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky every few nights including the Galactic plane. ZTF observes $\sim10^9$ stars in g-band and r-band and can significantly contribute to the observed microlensing population. We predict that ZTF will observe $\sim$1100 microlensing events in three years of observing within $10^\circ$ degrees latitude of the Galactic plane, with $\sim$500 events in the outer Galaxy ($\ell \geq 10^\circ$). This yield increases to $\sim$1400 ($\sim$800) events by combining every three ZTF exposures, $\sim$1800 ($\sim$900) events if ZTF observes for a total of five years, and $\sim$2400 ($\sim$1300) events for a five year survey with post-processing image stacking. Using the microlensing modeling software PopSyCLE, we compare the microlensing populations in the Galactic bulge and the outer Galaxy. We also present an analysis of the microlensing event ZTF18abhxjmj to demonstrate how to leverage these population statistics in event modeling. ZTF will constrain Galactic structure, stellar populations, and primordial black holes through photometric microlensing. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Comments: 19 pages, 13 figures, 5 tables, accepted to ApJ (6/4/2020), microlensing simulation catalogs available at https://portal.nersc.gov/project/uLens/Galactic_Microlensing_Distributions

Journal ref: The Astrophysical Journal, 2020, Volume 897, Number 2, 144

Search for a Variation of the Fine Structure around the Supermassive Black Hole in Our Galactic Center

Authors: A. Hees, T. Do, B. M. Roberts, A. M. Ghez, S. Nishiyama, R. O. Bentley, A. K. Gautam, S. Jia, T. Kara, J. R. Lu, H. Saida, S. Sakai, M. Takahashi, Y. Takamori

Abstract: Searching for space-time variations of the constants of Nature is a promising way to search for new physics beyond General Relativity and the standard model motivated by unification theories and models of dark matter and dark energy. We propose a new way to search for a variation of the fine-structure constant using measurements of late-type evolved giant stars from the S-star cluster orbiting the… ▽ More Searching for space-time variations of the constants of Nature is a promising way to search for new physics beyond General Relativity and the standard model motivated by unification theories and models of dark matter and dark energy. We propose a new way to search for a variation of the fine-structure constant using measurements of late-type evolved giant stars from the S-star cluster orbiting the supermassive black hole in our Galactic Center. A measurement of the difference between distinct absorption lines (with different sensitivity to the fine structure constant) from a star leads to a direct estimate of a variation of the fine structure constant between the star's location and Earth. Using spectroscopic measurements of 5 stars, we obtain a constraint on the relative variation of the fine structure constant below $10^{-5}$. This is the first time a varying constant of Nature is searched for around a black hole and in a high gravitational potential. This analysis shows new ways the monitoring of stars in the Galactic Center can be used to probe fundamental physics. △ Less

Submitted 26 February, 2020; originally announced February 2020.

Comments: 7 pages + 10 pages appendix, 3 figures, version accepted for publication

Journal ref: Phys. Rev. Letters 124, 081101, 2020

PopSyCLE: A New Population Synthesis Code for Compact Object Microlensing Events

Authors: Casey Y. Lam, Jessica R. Lu, Matthew W. Hosek, Jr., William A. Dawson, Nathan R. Golovich

Abstract: We present a new Milky Way microlensing simulation code, dubbed PopSyCLE (Population Synthesis for Compact object Lensing Events). PopSyCLE is the first resolved microlensing simulation to include a compact object distribution derived from numerical supernovae explosion models and both astrometric and photometric microlensing effects. We demonstrate the capabilities of PopSyCLE by investigating th… ▽ More We present a new Milky Way microlensing simulation code, dubbed PopSyCLE (Population Synthesis for Compact object Lensing Events). PopSyCLE is the first resolved microlensing simulation to include a compact object distribution derived from numerical supernovae explosion models and both astrometric and photometric microlensing effects. We demonstrate the capabilities of PopSyCLE by investigating the optimal way to find black holes (BHs) with microlensing. Candidate BHs have typically been selected from wide-field photometric microlensing surveys, such as OGLE, by selecting events with long Einstein crossing times ($t_E>120$ days). These events can be selected at closest approach and monitored astrometrically in order to constrain the mass of each lens; PopSyCLE predicts a BH detection rate of ~40% for such a program. We find that the detection rate can be enhanced to ~85% by selecting events with both $t_E>120$ days and a microlensing parallax of $π_E<0.08$. Unfortunately, such a selection criterion cannot be applied during the event as $π_E$ requires both pre- and post-peak photometry. However, historical microlensing events from photometric surveys can be revisited using this new selection criteria in order to statistically constrain the abundance of BHs in the Milky Way. The future WFIRST microlensing survey provides both precise photometry and astrometry and will yield individual masses of $\mathcal{O}(100-1000) black holes, which is at least an order of magnitude more than is possible with individual candidate follow-up with current facilities. The resulting sample of BH masses from WFIRST will begin to constrain the shape of the black hole present-day mass function, BH multiplicity, and BH kick velocity distributions. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: 33 pages, 22 figures, accepted to ApJ

Keck Observations Confirm a Super-Jupiter Planet Orbiting M-dwarf OGLE-2005-BLG-071L

Authors: David P. Bennett, Aparna Bhattacharya, Jean-Philippe Beaulieu, Joshua W. Blackman, Aikaterini Vandorou, Sean K. Terry, Andrew A. Cole, Calen B. Henderson, Naoki Koshimoto, Jessica R. Lu, Jean Baptiste Marquette, Clement Ranc, Andrzej Udalski

Abstract: We present adaptive optics imaging from the NIRC2 instrument on the Keck-2 telescope that resolves the exoplanet host (and lens) star as it separates from the brighter source star. These observations yield the $K$-band brightness of the lens and planetary host star, as well as the lens-source relative proper motion, $μ_{\rm rel,H}$. in the heliocentric reference frame. The $μ_{\rm rel,H}$ measurem… ▽ More We present adaptive optics imaging from the NIRC2 instrument on the Keck-2 telescope that resolves the exoplanet host (and lens) star as it separates from the brighter source star. These observations yield the $K$-band brightness of the lens and planetary host star, as well as the lens-source relative proper motion, $μ_{\rm rel,H}$. in the heliocentric reference frame. The $μ_{\rm rel,H}$ measurement allows determination of the microlensing parallax vector, $π_E$, which had only a single component determined by the microlensing light curve. The combined measurements of $μ_{\rm rel,H}$ and $K_L$ provide the masses of the host stat, $M_{\rm host} = 0.426\pm 0.037 M_\odot$, and planet, $m_p = 3.27 \pm 0.32 M_{\rm Jup}$ with a projected separation of $3.4\pm 0.5\,$AU. This confirms the tentative conclusion of a previous paper (Dong et al. 2009) that this super-Jupiter mass planet, OGLE-2005-BLG-071Lb, orbits an M-dwarf. Such planets are predicted to be rare by the core accretion theory and have been difficult to find with other methods, but there are two such planets with firm mass measurements from microlensing, and an additional 11 planetary microlens events with host mass estimates $< 0.5M_\odot$ and planet mass estimates $> 2$ Jupiter masses that could be confirmed by high angular follow-up observations. We also point out that OGLE-2005-BLG-071L has separated far enough from its host star that it should be possible to measure the host star metallicity withspectra from a high angular resolution telescope such as Keck, the VLT, the Hubble Space Telescope or the James Webb Space Telescope. △ Less

Submitted 12 January, 2020; v1 submitted 10 September, 2019; originally announced September 2019.

Comments: accepted by AJ

Geometric Distortion Calibration with Photo-lithographic Pinhole Masks for High-Precision Astrometry

Authors: Maxwell Service, Jessica R. Lu, Mark Chun, Ryuiji Suzuki, Matthias Schoeck, Jenny Atwood, David Andersen, Glen Herriot

Abstract: Adaptive optics (AO) systems deliver high-resolution images that may be ideal for precisely measuring positions of stars (i.e. astrometry) if the system has stable and well-calibrated geometric optical distortions. A calibration unit, equipped with back-illuminated pinhole mask, can be utilized to measure instrumental optical distortions. AO systems on the largest ground-based telescopes, such as… ▽ More Adaptive optics (AO) systems deliver high-resolution images that may be ideal for precisely measuring positions of stars (i.e. astrometry) if the system has stable and well-calibrated geometric optical distortions. A calibration unit, equipped with back-illuminated pinhole mask, can be utilized to measure instrumental optical distortions. AO systems on the largest ground-based telescopes, such as the W. M. Keck Observatory and the Thirty Meter Telescope require pinhole positions known to 20 nm to achieve an astrometric precision of 0.001 of a resolution element. We characterize a photo-lithographic pinhole mask and explore the systematic errors that result from different experimental setups. We characterized the nonlinear geometric distortion of a simple imaging system using the mask; and we measured 857 nm RMS of optical distortion with a final residual of 39 nm (equivalent to 20 μas for TMT). We use a sixth order bivariate Legendre polynomial to model the optical distortion and allow the reference positions of the individual pinholes to vary. The nonlinear deviations in the pinhole pattern with respect to the manufacturing design of a square pattern are 47.2 nm +/- 4.5 nm (random) +/- 10.8 nm (systematic) over an area of 1788 mm$^2$. These deviations reflect the additional error induced when assuming the pinhole mask is manufactured perfectly square. We also find that ordered mask distortions are significantly more difficult to characterize than random mask distortions as the ordered distortions can alias into optical camera distortion. Future design simulations for astrometric calibration units should include ordered mask distortions. We conclude that photo-lithographic pinhole masks are >10 times better than the pinhole masks deployed in first generation AO systems and are sufficient to meet the distortion calibration requirements for the upcoming thirty meter class telescopes. △ Less

Submitted 13 August, 2019; originally announced August 2019.

Comments: 20 pages, 13 figures

Relativistic redshift of the star S0-2 orbiting the Galactic center supermassive black hole

Authors: Tuan Do, Aurelien Hees, Andrea Ghez, Gregory D. Martinez, Devin S. Chu, Siyao Jia, Shoko Sakai, Jessica R. Lu, Abhimat K. Gautam, Kelly Kosmo O'Neil, Eric E. Becklin, Mark R. Morris, Keith Matthews, Shogo Nishiyama, Randy Campbell, Samantha Chappell, Zhuo Chen, Anna Ciurlo, Arezu Dehghanfar, Eulalia Gallego-Cano, Wolfgang E. Kerzendorf, James E. Lyke, Smadar Naoz, Hiromi Saida, Rainer Schödel , et al. (4 additional authors not shown)

Abstract: General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach… ▽ More General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach to the black hole. We detect the combination of special relativistic- and gravitational-redshift, quantified using a redshift parameter, $Υ$. Our result, $Υ=0.88 \pm 0.17$, is consistent with General Relativity ($Υ=1$) and excludes a Newtonian model ($Υ=0$ ) with a statistical significance of 5 $σ$. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 78 pages, 21 figures, accepted to Science, astrometry and radial velocity measurements, and posterior chains available here, http://www.astro.ucla.edu/~ghezgroup/gc/gr/data.zip

An adaptive scheduling tool to optimize measurements to reach a scientific objective: methodology and application to the measurements of stellar orbits in the Galactic Center

Authors: A. Hees, A. Dehghanfar, T. Do, A. M. Ghez, G. D. Martinez, R. Campbell, J. R. Lu

Abstract: In various fields of physics and astronomy, access to experimental facilities or to telescopes is becoming more and more competitive and limited. It becomes therefore important to optimize the type of measurements and their scheduling to reach a given scientific objective and to increase the chances of success of a scientific project. In this communication, extending the work of Ford (2008) and of… ▽ More In various fields of physics and astronomy, access to experimental facilities or to telescopes is becoming more and more competitive and limited. It becomes therefore important to optimize the type of measurements and their scheduling to reach a given scientific objective and to increase the chances of success of a scientific project. In this communication, extending the work of Ford (2008) and of Loredo et al. (2012), we present an efficient adaptive scheduling tool aimed at prioritzing measurements in order to reach a scientific goal. The algorithm, based on the Fisher matrix, can be applied to a wide class of measurements. We present this algorithm in detail and discuss some practicalities such as systematic errors or measurements losses due to contigencies (such as weather, experimental failure, ...). As an illustration, we consider measurements of the short-period star S0-2 in our Galactic Center. We show that the radial velocity measurements at the two turning points of the radial velocity curve are more powerful for detecting the gravitational redshift than measurements at the maximal relativistic signal. We also explicitly present the methodology that was used to plan measurements in order to detect the relativistic redshift considering systematics and possible measurements losses. For the future, we identify the astrometric turning points to be highly sensitive to the relativistic advance of the periastron. Finally, we also identify measurements particularly sensitive to the distance to our Galactic Center: the radial velocities around periastron and the astrometric measurements just before closest approach and at the maximal right ascension astrometric turning point. △ Less

Submitted 7 June, 2019; originally announced June 2019.

Comments: 12 pages, 6 figures, submitted

Journal ref: ApJ 880, 87, 2019

The Quintuplet Cluster: Extended Structure and Tidal Radius

Authors: Nicholas Z. Rui, Matthew W. Hosek Jr., Jessica R. Lu, William I. Clarkson, Jay Anderson, Mark R. Morris, Andrea M. Ghez

Abstract: The Quintuplet star cluster is one of only three known young ($<10$ Myr) massive (M $>10^4$ M$_\odot$) clusters within $\sim100$ pc of the Galactic Center. In order to explore star cluster formation and evolution in this extreme environment, we analyze the Quintuplet's dynamical structure. Using the HST WFC3-IR instrument, we take astrometric and photometric observations of the Quintuplet covering… ▽ More The Quintuplet star cluster is one of only three known young ($<10$ Myr) massive (M $>10^4$ M$_\odot$) clusters within $\sim100$ pc of the Galactic Center. In order to explore star cluster formation and evolution in this extreme environment, we analyze the Quintuplet's dynamical structure. Using the HST WFC3-IR instrument, we take astrometric and photometric observations of the Quintuplet covering a $120''\times120''$ field-of-view, which is $19$ times larger than those of previous proper motion studies of the Quintuplet. We generate a catalog of the Quintuplet region with multi-band, near-infrared photometry, proper motions, and cluster membership probabilities for $10,543$ stars. We present the radial density profile of $715$ candidate Quintuplet cluster members with $M\gtrsim4.7$ M$_\odot$ out to $3.2$ pc from the cluster center. A $3σ$ lower limit of $3$ pc is placed on the tidal radius, indicating the lack of a tidal truncation within this radius range. Only weak evidence for mass segregation is found, in contrast to the strong mass segregation found in the Arches cluster, a second and slightly younger massive cluster near the Galactic Center. It is possible that tidal stripping hampers a mass segregation signature, though we find no evidence of spatial asymmetry. Assuming that the Arches and Quintuplet formed with comparable extent, our measurement of the Quintuplet's comparatively large core radius of $0.62^{+0.10}_{-0.10}$ pc provides strong empirical evidence that young massive clusters in the Galactic Center dissolve on a several Myr timescale. △ Less

Submitted 4 April, 2019; originally announced April 2019.

Comments: 25 pages (21-page main text, 4-page appendix), 18 figures, submitted to ApJ

Astro2020: From Stars to Compact Objects: The Initial-Final Mass Relation

Authors: Jessica R. Lu, Casey Lam, Will Dawson, B. Scott Gaudi, Nathan Golovich, Michael Medford, Fatima Abdurrahman, Rachael L. Beaton

Abstract: One of the key phases of stellar evolution that remains poorly understood is stellar death. We lack a predictive model for how a star of a given mass explodes and what kind of remnant it leaves behind (i.e. the initial-final mass relation, IFMR). Progress has been limited due to the difficulty in finding and weighing black holes and neutron stars in large numbers. Technological advances that allow… ▽ More One of the key phases of stellar evolution that remains poorly understood is stellar death. We lack a predictive model for how a star of a given mass explodes and what kind of remnant it leaves behind (i.e. the initial-final mass relation, IFMR). Progress has been limited due to the difficulty in finding and weighing black holes and neutron stars in large numbers. Technological advances that allow for sub-milliarcsecond astrometry in crowded fields have opened a new window for finding black holes and neutron stars: astrometric gravitational lensing. Finding and weighing a sample of compact objects with astrometric microlensing will allow us to place some of the first constraints on the present-day mass function of isolated black holes and neutron stars, their multiplicity, and their kick velocities. All of these are fundamental inputs into understanding the death phase of stellar evolution, improving supernovae models, and interpreting LIGO detections in an astrophysical context. To achieve these goals, we require large area surveys, such as the WFIRST exoplanet microlensing survey, to photometrically identify long-duration (>120 day), un-blended microlensing events as candidate compact objects. We also require high-precision astrometric follow-up monitoring using extremely large telescopes, equipped with adaptive optics, such as TMT and GMT. △ Less

Submitted 3 April, 2019; originally announced April 2019.

Comments: Astro2020 Science White Paper

Primordial Black Hole Microlensing: The Einstein Crossing Time Distribution

Authors: Jessica R. Lu, Casey Y. Lam, Michael Medford, William Dawson, Nathan Golovich

Abstract: Gravitational microlensing is one of the few means of finding primordial black holes (PBHs), if they exist. Recent LIGO detections of 30 Msun black holes have re-invigorated the search for PBHs in the 10-100 Msun mass regime. Unfortunately, individual PBH microlensing events cannot easily be distinguished from stellar lensing events from photometry alone. However, the distribution of microlensing… ▽ More Gravitational microlensing is one of the few means of finding primordial black holes (PBHs), if they exist. Recent LIGO detections of 30 Msun black holes have re-invigorated the search for PBHs in the 10-100 Msun mass regime. Unfortunately, individual PBH microlensing events cannot easily be distinguished from stellar lensing events from photometry alone. However, the distribution of microlensing timescales (tE, the Einstein radius crossing time) can be analyzed in a statistical sense using models of the Milky Way with and without PBHs. While previous works have presented both theoretical models and observational constrains for PBHs (e.g. Calcino et al. 2018; Niikura et al. 2019), surprisingly, they rarely show the observed quantity -- the tE distribution -- for different abundances of PBHs relative to the total dark matter mass (fPBH). We present a simple calculation of how the tE distribution changes between models with and without PBHs. △ Less

Submitted 3 April, 2019; originally announced April 2019.

Comments: RNAAS

Journal ref: 2019, RNAAS, Vol 3, 58

"Auxiliary" Science with the WFIRST Microlensing Survey

Authors: B. Scott Gaudi, Rachel Akeson, Jay Anderson, Etienne Bachelet, David P. Bennett, Aparna Bhattacharya, Valerio Bozza, Sebastiano Calchi Novati, Calen B. Henderson, Samson A. Johnson, Jeffrey Kruk, Jessica R. Lu, Shude Mao, Benjamin T. Montet, David M. Nataf, Matthew T. Penny, Radoslaw Poleski, Clément Ranc, Kailash Sahu, Yossi Shvartzvald, David N. Spergel, Daisuke Suzuki, Keivan G. Stassun, Rachel A. Street

Abstract: The Wide Field Infrared Survey Telescope (WFIRST) will monitor $\sim 2$ deg$^2$ toward the Galactic bulge in a wide ($\sim 1-2~μ$m) W149 filter at 15-minute cadence with exposure times of $\sim$50s for 6 seasons of 72 days each, for a total $\sim$41,000 exposures taken over $\sim$432 days, spread over the 5-year prime mission. This will be one of the deepest exposures of the sky ever taken, reachi… ▽ More The Wide Field Infrared Survey Telescope (WFIRST) will monitor $\sim 2$ deg$^2$ toward the Galactic bulge in a wide ($\sim 1-2~μ$m) W149 filter at 15-minute cadence with exposure times of $\sim$50s for 6 seasons of 72 days each, for a total $\sim$41,000 exposures taken over $\sim$432 days, spread over the 5-year prime mission. This will be one of the deepest exposures of the sky ever taken, reaching a photon-noise photometric precision of 0.01 mag per exposure and collecting a total of $\sim 10^9$ photons over the course of the survey for a W149$_{\rm AB}\sim 21$ star. Of order $4 \times 10^7$ stars will be monitored with W149$_{\rm AB}$<21, and 10$^8$ stars with W145$_{\rm AB}$<23. The WFIRST microlensing survey will detect $\sim$54,000 microlensing events, of which roughly 1% ($\sim$500) will be due to isolated black holes, and $\sim$3% ($\sim$1600) will be due to isolated neutron stars. It will be sensitive to (effectively) isolated compact objects with masses as low as the mass of Pluto, thereby enabling a measurement of the compact object mass function over 10 orders of magnitude. Assuming photon-noise limited precision, it will detect $\sim 10^5$ transiting planets with sizes as small as $\sim 2~R_\oplus$, perform asteroseismology of $\sim 10^6$ giant stars, measure the proper motions to $\sim 0.3\%$ and parallaxes to $\sim 10\%$ for the $\sim 6 \times 10^6$ disk and bulge stars in the survey area, and directly detect $\sim 5 \times 10^3$ Trans-Neptunian objects (TNOs) with diameters down to $\sim 10$ km, as well as detect $\sim 10^3$ occulations of stars by TNOs during the survey. All of this science will completely serendipitous, i.e., it will not require modifications of the WFIRST optimal microlensing survey design. Allowing for some minor deviation from the optimal design, such as monitoring the Galactic center, would enable an even broader range of transformational science. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: 9 pages, 3 figures, submitted to the Astro2020 Science White Paper call. This white paper draws heavily from Chapter 2.5 of Spergel et al. 2015 (arXiv:1503.03757), which was written by the lead author of this white paper