Showing 1–50 of 69 results for author: Lam, C

The Symbiotic X-ray Binary IGR J16194-2810: A Window on the Future Evolution of Wide Neutron Star Binaries From Gaia

Authors: Pranav Nagarajan, Kareem El-Badry, Casey Lam, Henrique Reggiani

Abstract: We present optical follow-up of IGR J16194-2810, a hard X-ray source discovered by the INTEGRAL mission. The optical counterpart is a $\sim500\,L_\odot$ red giant at a distance of $2.1$ kpc. We measured 17 radial velocities (RVs) of the giant over a period of $271$ days. Fitting these RVs with a Keplerian model, we find an orbital period of $P_{\rm orb} = 192.73 \pm 0.01$ days and a companion mass… ▽ More We present optical follow-up of IGR J16194-2810, a hard X-ray source discovered by the INTEGRAL mission. The optical counterpart is a $\sim500\,L_\odot$ red giant at a distance of $2.1$ kpc. We measured 17 radial velocities (RVs) of the giant over a period of $271$ days. Fitting these RVs with a Keplerian model, we find an orbital period of $P_{\rm orb} = 192.73 \pm 0.01$ days and a companion mass function $f(M_2) = 0.365 \pm 0.003 \,M_{\odot}$. We detect ellipsoidal variability with the same period in optical light curves from the ASAS-SN survey. Joint fitting of the RVs, light curves, and the broadband SED allows us to robustly constrain the masses of both components. We find a giant mass of $M_\star = 0.99^{+0.02}_{-0.03}\,M_{\odot}$ and a companion mass of $M_{2} = 1.23^{+0.05}_{-0.03}\,M_{\odot}$, implying that the companion is a neutron star (NS). We recover a $4.06$-hour period in the system's TESS light curve, which we tentatively associate with the NS spin period. The giant does not yet fill its Roche lobe, suggesting that current mass transfer is primarily via winds. MESA evolutionary models predict that the giant will overflow its Roche lobe in $5$-$10$ Myr, eventually forming a recycled pulsar + white dwarf binary with a $\sim 900$ day period. IGR J16194-2810 provides a window on the future evolution of wide NS + main sequence binaries recently discovered via Gaia astrometry. As with those systems, the binary's formation history is uncertain. Before the formation of the NS, it likely survived a common envelope episode with a donor-to-accretor mass ratio $\gtrsim 10$ and emerged in a wide orbit. The NS likely formed with a weak kick ($v_{\rm kick}\lesssim 50\,\rm km\,s^{-1}$), as stronger kicks would have disrupted the orbit. △ Less

Submitted 2 July, 2024; v1 submitted 27 May, 2024; originally announced May 2024.

Comments: 21 pages, 13 figures, Accepted to PASP

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Ages of "singles" versus "multis": Predictions for dynamical sculpting over Gyr in the Kepler Sample

Authors: Christopher Lam, Sarah Ballard

Abstract: The sample of host stars with multiple transiting planets has illuminated the orbital architectures of exoplanetary systems. These architectures may be shaped mostly by formation conditions, be continually sculpted by ongoing dynamical processes, or both. As more studies place planet occurrence within a galactic context, evidence has emerged for variable planet multiplicity over time. In this manu… ▽ More The sample of host stars with multiple transiting planets has illuminated the orbital architectures of exoplanetary systems. These architectures may be shaped mostly by formation conditions, be continually sculpted by ongoing dynamical processes, or both. As more studies place planet occurrence within a galactic context, evidence has emerged for variable planet multiplicity over time. In this manuscript, we investigate the use of transit multiplicity as a tool to constrain longer-timescale (>1 Gyr) dynamical sculpting. First, with a suite of injection-and-recovery tests, we quantify sensitivity to sculpting laws across different regimes. We employ a forward modeling framework in which we generate synthetic planetary systems, according to a prescribed sculpting speed and timescale, around the FGK dwarfs studied by the Kepler Mission. Some sculpting scenarios are hypothetically detectable in the Kepler sample, while others can be disfavored from Kepler transit statistics alone. Secondly, we apply our analysis to reverse-engineer the sculpting laws consistent with the true yield from Kepler. We confirm the present-day fraction of host stars containing dynamically cool "systems with tightly-packed inner planets" (STIPs) is 4-13%. A variety of Gyr-timescale sculpting laws successfully predict the transit multiplicity of the Kepler sample, but none of these laws succeeds in also producing a detectable trend with transit multiplicity and stellar age. An improvement to measured stellar age precision may help uncover such a sculpting law, but nevertheless reflects limitations in transit multiplicity as an observable. Therefore other phenomena, apart from Gyr-timescale dynamical sculpting, may be required to explain the Kepler yield. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Nuclear mass table in deformed relativistic Hartree-Bogoliubov theory in continuum, II: Even-$Z$ nuclei

Authors: DRHBc Mass Table Collaboration, Peng Guo, Xiaojie Cao, Kangmin Chen, Zhihui Chen, Myung-Ki Cheoun, Yong-Beom Choi, Pak Chung Lam, Wenmin Deng, Jianmin Dong, Pengxiang Du, Xiaokai Du, Kangda Duan, Xiaohua Fan, Wei Gao, Lisheng Geng, Eunja Ha, Xiao-Tao He, Jinniu Hu, Jingke Huang, Kun Huang, Yanan Huang, Zidan Huang, Kim Da Hyung, Hoi Yat Chan , et al. (58 additional authors not shown)

Abstract: The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-ne… ▽ More The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-neutron separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. A total of 4829 even-$Z$ nuclei are predicted to be bound, with an rms deviation of 1.477 MeV from the 1244 mass data. Good agreement with the available experimental odd-even mass differences, $α$ decay energies, and charge radii is also achieved. The description accuracy for nuclear masses and nucleon separation energies as well as the prediction for drip lines is compared with the results obtained from other relativistic and nonrelativistic density functional. The comparison shows that the DRHBc theory with PC-PK1 provides an excellent microscopic description for the masses of even-$Z$ nuclei. The systematics of the nucleon separation energies, odd-even mass differences, pairing energies, two-nucleon gaps, $α$ decay energies, rms radii, quadrupole deformations, potential energy curves, neutron density distributions, and neutron mean-field potentials are discussed. △ Less

Submitted 10 June, 2024; v1 submitted 5 February, 2024; originally announced February 2024.

Comments: 394 pages, 17 figures, 2 tables, published in Atomic Data and Nuclear Data Tables, data file in the TXT form is available for download under "Ancillary files"

Journal ref: Peng Guo, et. al. (DRHBc Mass Table Collaboration), Atomic Data and Nuclear Data Tables 158 (2024) 101661

Light-Curve Structure and Halpha Line Formation in the Tidal Disruption Event AT 2019azh

Authors: Sara Faris, Iair Arcavi, Lydia Makrygianni, Daichi Hiramatsu, Giacomo Terreran, Joseph Farah, D. Andrew Howell, Curtis McCully, Megan Newsome, Estefania Padilla Gonzalez, Craig Pellegrino, K. Azalee Bostroem, Wiam Abojanb, Marco C. Lam, Lina Tomasella, Thomas G. Brink, Alexei V. Filippenko, K. Decker French, Peter Clark, Or Graur, Giorgos Leloudas, Mariusz Gromadzki, Joseph P. Anderson, Matt Nicholl, Claudia P. Gutierrez , et al. (11 additional authors not shown)

Abstract: AT 2019azh is a H+He tidal disruption event (TDE) with one of the most extensive ultraviolet and optical datasets available to date. We present our photometric and spectroscopic observations of this event starting several weeks before and out to approximately two years after g-band peak brightness and combine them with public photometric data. This extensive dataset robustly reveals a change in th… ▽ More AT 2019azh is a H+He tidal disruption event (TDE) with one of the most extensive ultraviolet and optical datasets available to date. We present our photometric and spectroscopic observations of this event starting several weeks before and out to approximately two years after g-band peak brightness and combine them with public photometric data. This extensive dataset robustly reveals a change in the light-curve slope and a bump in the rising light curve of a TDE for the first time, which may indicate more than one dominant emission mechanism contributing to the pre-peak light curve. We further confirm the relation seen in previous TDEs whereby the redder emission peaks later than the bluer emission. The post-peak bolometric light curve of AT 2019azh is better described by an exponential decline than by the canonical t^{-5/3} (and in fact any) power-law decline. We find a possible mid-infrared excess around peak optical luminosity, but cannot determine its origin. In addition, we provide the earliest measurements of the Halpha emission-line evolution and find no significant time delay between the peak of the V-band light curve and that of the Halpha luminosity. These results can be used to constrain future models of TDE line formation and emission mechanisms in general. More pre-peak 1-2 day cadence observations of TDEs are required to determine whether the characteristics observed here are common among TDEs. More importantly, detailed emission models are needed to fully exploit such observations for understanding the emission physics of TDEs. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: Submitted to ApJ

Microlensing Events in Five Years of Photometry from the Zwicky Transient Facility

Authors: Ruocheng Zhai, Antonio C. Rodriguez, Casey Y. Lam, Eric C. Bellm, Josiah Purdum, Frank J. Masci, Avery Wold

Abstract: Microlensing has a unique advantage for detecting dark objects in the Milky Way, such as free floating planets, neutron stars, and stellar-mass black holes. Most microlensing surveys focus towards the Galactic bulge, where higher stellar density leads to a higher event rate. However, microlensing events in the Galactic plane are closer, and take place over longer timescales. This enables a better… ▽ More Microlensing has a unique advantage for detecting dark objects in the Milky Way, such as free floating planets, neutron stars, and stellar-mass black holes. Most microlensing surveys focus towards the Galactic bulge, where higher stellar density leads to a higher event rate. However, microlensing events in the Galactic plane are closer, and take place over longer timescales. This enables a better measurement of the microlensing parallax, which serves as an independent constraint on the mass of the dark lens. In this work, we systematically searched for microlensing events in Zwicky Transient Facility (ZTF) Data Release 17 from 2018--2023 in the Galactic plane region $|b| < 20^\circ$. We find 124 high-confidence microlensing events and 54 possible events. In the event selection, we use the efficient \texttt{EventFinder} algorithm to detect microlensing signals, which could be used for large datasets such as future ZTF data releases or data from the Rubin Observatory Legacy Survey of Space and Time (LSST). With detection efficiencies of ZTF fields from catalog-level simulations, we calculate the mean Einstein timescale to be $\langle t_\mathrm{E}\rangle = 51.7 \pm 3.3$ days, smaller than previous results of the Galactic plane to within 1.5-$σ$. We calculate optical depths and event rates, which we interpret with caution due to the use of visual inspection in creating our final sample. With two years of additional ZTF data in DR17, we have more than doubled the amount of microlensing events (60) found in the three-year DR5 search and found events with longer Einstein timescales than before. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Comments: 9 figures, 3 tables. Submitted to ApJ

Natal Kicks from the Galactic Center and Implications on their Environment and the Roman Space Telescope

Authors: Carlos Jurado, Smadar Naoz, Casey Y. Lam, Bao-Minh Hoang

Abstract: Most galaxies, including the Milky Way, harbor a central supermassive black hole (SMBH) weighing millions to billions of solar masses. Surrounding these SMBHs are dense regions of stars and stellar remnants, such as neutron stars and black holes. Neutron stars and possibly black holes receive large natal kicks at birth on the order of hundreds of km s$^{-1}$. The natal kicks that occur in the vici… ▽ More Most galaxies, including the Milky Way, harbor a central supermassive black hole (SMBH) weighing millions to billions of solar masses. Surrounding these SMBHs are dense regions of stars and stellar remnants, such as neutron stars and black holes. Neutron stars and possibly black holes receive large natal kicks at birth on the order of hundreds of km s$^{-1}$. The natal kicks that occur in the vicinity of an SMBH may redistribute the orbital configuration of the compact objects and alter their underlying density distribution. We model the effects of natal kicks on a Galactic Center (GC) population of massive stars and stellar binaries with different initial density distributions. Using observational constraints from stellar orbits near the GC, we place an upper limit on the steepness of the initial stellar profile and find it to be core-like. In addition, we predict that $30-70 \%$ of compact objects become unbound from the SMBH due to their kicks and will migrate throughout the galaxy. Different black hole kick prescriptions lead to distinct spatial and kinematic distributions. We suggest that the Roman Space Telescope may be able to distinguish between these distributions and thus be able to differentiate natal kick mechanisms. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 18 pages, 11 Figures

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

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]

A hidden population of massive white dwarfs: two spotted K+WD binaries

Authors: D. M. Rowan, T. Jayasinghe, M. A. Tucker, C. Y. Lam, Todd A. Thompson, C. S. Kochanek, N. S. Abrams, B. J. Fulton, I. Ilyin, H. Isaacson, J. Lu, D. V. Martin, B. Nicholson

Abstract: The identification and characterization of massive ($\gtrsim 0.8~M_\odot$) white dwarfs is challenging in part due to their low luminosity. Here we present two candidate single-lined spectroscopic binaries, Gaia DR3 4014708864481651840 and 5811237403155163520, with K-dwarf primaries and optically dark companions. Both have orbital periods of $P\sim 0.45$ days and show rotational variability, ellip… ▽ More The identification and characterization of massive ($\gtrsim 0.8~M_\odot$) white dwarfs is challenging in part due to their low luminosity. Here we present two candidate single-lined spectroscopic binaries, Gaia DR3 4014708864481651840 and 5811237403155163520, with K-dwarf primaries and optically dark companions. Both have orbital periods of $P\sim 0.45$ days and show rotational variability, ellipsoidal modulations, and high-amplitude radial velocity variations. Using light curves from the Transiting Exoplanet Survey Satellite (TESS), radial velocities from ground-based spectrographs, and spectral energy distributions, we characterize these binaries to describe the nature of the unseen companion. We find that both systems are consistent with a massive white dwarf companion. Unlike simple ellipsoidal variables, star spots cause the light curve morphology to change between TESS sectors. We attempt to constrain the orbital inclination using PHOEBE binary light curve models, but degeneracies in the light curves of spotted stars prevent a precise determination. Finally, we search for similar objects using Gaia DR3 and TESS, and comment on these systems in the context of recently claimed compact object binaries. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 17 pages, 15 figures. Submitted to MNRAS

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

Rapidly Evolving Transients in Archival ZTF Public Alerts

Authors: Wenxiong Li, Iair Arcavi, Ehud Nakar, Alexei V. Filippenko, Thomas G. Brink, WeiKang Zheng, Marco C. Lam, Ido Keinan, Seán J. Brennan, Noi Shitrit

Abstract: We search the archival Zwicky Transient Facility public survey for rapidly evolving transient (RET) candidates based on well-defined criteria between 2018 May and 2021 December. The search yielded 19 bona-fide RET candidates, corresponding to a discovery rate of $\sim 5.2$ events per year. Even with a Galactic latitude cut of $20^\circ$, 8 of the 19 events ($\sim 42$%) are Galactic, including one… ▽ More We search the archival Zwicky Transient Facility public survey for rapidly evolving transient (RET) candidates based on well-defined criteria between 2018 May and 2021 December. The search yielded 19 bona-fide RET candidates, corresponding to a discovery rate of $\sim 5.2$ events per year. Even with a Galactic latitude cut of $20^\circ$, 8 of the 19 events ($\sim 42$%) are Galactic, including one with a light-curve shape closely resembling that of the GW170817 kilonova (KN). An additional event is a nova in M31. Four out of the 19 events ($\sim 21$%) are confirmed extragalactic RETs (one confirmed here for the first time) and the origin of 6 additional events cannot be determined. We did not find any extragalactic events resembling the GW170817 KN, from which we obtain an upper limit on the volumetric rate of GW170817-like KNe of $R \le$ 2400 Gpc$^{-3}$ yr$^{-1}$ (95% confidence). These results can be used for quantifying contaminants to RET searches in transient alert streams, specifically when searching for kilonovae independently of gravitational-wave and gamma-ray-burst triggers. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: 16 pages, 7 figures, submitted to AAS Journals

AT 2021loi: A Bowen Fluorescence Flare with a Rebrightening Episode, Occurring in a Previously-Known AGN

Authors: Lydia Makrygianni, Benny Trakhtenbrot, Iair Arcavi, Claudio Ricci, Marco C. Lam, Assaf Horesh, Itai Sfaradi, K. Azalee Bostroem, Griffin Hosseinzadeh, D. Andrew Howell, Craig Pellegrino, Rob Fender, David A. Green, David R. A. Williams, Joe Bright

Abstract: AT 2021loi is an optical-ultraviolet transient located at the center of its host galaxy. Its spectral features identify it as a member of the ``Bowen Fluorescence Flare'' (BFF) class. The first member of this class was considered to be related to a tidal disruption event, but enhanced accretion onto an already active supermassive black hole was suggested as an alternative explanation. AT 2021loi,… ▽ More AT 2021loi is an optical-ultraviolet transient located at the center of its host galaxy. Its spectral features identify it as a member of the ``Bowen Fluorescence Flare'' (BFF) class. The first member of this class was considered to be related to a tidal disruption event, but enhanced accretion onto an already active supermassive black hole was suggested as an alternative explanation. AT 2021loi, having occurred in a previously-known unobscured AGN, strengthens the latter interpretation. Its light curve is similar to those of previous BFFs, showing a rebrightening approximately one year after the main peak (which was not explicitly identified, but might be the case, in all previous BFFs). An emission feature around 4680 A, seen in the pre-flare spectrum, strengthens by a factor of $\sim$2 around the optical peak of the flare, and is clearly seen as a double peaked feature then, suggesting a blend of NIII $λ4640$ with HeII $\lambda4686$ as its origin. The appearance of OIII $λ$3133 and possible NIII $λ\lambda4097,4103$ (blended with H$δ$) during the flare further support a Bowen Fluorescence classification. Here, we present ZTF, ATLAS, Keck, Las Cumbres Observatory, NEOWISE-R, $Swift$, AMI and VLA observations of AT 2021loi, making it one of the best observed BFFs to date. AT 2021loi thus provides some clarity on the nature of BFFs but also further demonstrates the diversity of nuclear transients. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: Submitted to ApJ. This version addresses comments from the referee

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

WDPhotTools -- A White Dwarf Photometric Toolkit in Python

Authors: M. C. Lam, K. W. Yuen, M. J. Green, W. Li

Abstract: From data collection to photometric fitting and analysis of white dwarfs to generating a white dwarf luminosity function requires numerous Astrophysical, Mathematical and Computational domain knowledge. The steep learning curve makes it difficult to enter the field, and often individuals have to reinvent the wheel to perform identical data reduction and analysis tasks. We have gathered a wide rang… ▽ More From data collection to photometric fitting and analysis of white dwarfs to generating a white dwarf luminosity function requires numerous Astrophysical, Mathematical and Computational domain knowledge. The steep learning curve makes it difficult to enter the field, and often individuals have to reinvent the wheel to perform identical data reduction and analysis tasks. We have gathered a wide range of publicly available white dwarf cooling models and synthetic photometry to provide a toolkit that allows (1) visualisation of various models, (2) photometric fitting of a white dwarf with or without distance and reddening, and (3) the computing of white dwarf luminosity functions with a choice of initial mass function, main sequence evolution model, star formation history, initial-final mass relation, and white dwarf cooling model. We have recomputed and compared the effective temperature of the white dwarfs from the Gaia EDR3 white dwarf catalogue. The two independent works show excellent agreement in the temperature solutions. △ Less

Submitted 22 October, 2022; v1 submitted 30 May, 2022; originally announced May 2022.

Comments: 18 pages, 12 figures, 2 tables. Accepted for publication in RAS Techniques and Instruments

R2-D2: Roman and Rubin -- From Data to Discovery

Authors: Suvi Gezari, Misty Bentz, Kishalay De, K. Decker French, Aaron Meisner, Michelle Ntampaka, Robert Jedicke, Ekta Patel, Daniel Perley, Robyn Sanderson, Christian Aganze, Igor Andreoni, Eric F. Bell, Edo Berger, Ian Dell'Antonio, Ryan Foley, Henry Hsieh, Mansi Kasliwal, Joel Kastner, Charles D. Kilpatrick, J. Davy Kirkpatrick, Casey Lam, Karen Meech, Dante Minniti, Ethan O. Nadler , et al. (6 additional authors not shown)

Abstract: The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to… ▽ More The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 29 pages, 12 figures, Table of Implementation Recommendations, Appendix of Community Science Pitches, AURA-commissioned whitepaper submitted to the Director of STScI (Ken Sembach) and the Director of NOIRLab (Pat McCarthy)

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

Identifying Blue Large Amplitude Pulsators from Gaia DR2 & ZTF DR3

Authors: Paul Ross McWhirter, Marco C. Lam

Abstract: Blue Large Amplitude Pulsators (BLAPs) are hot, subluminous stars undergoing rapid variability with periods of under 60 mins. They have been linked with the early stages of pre-white dwarfs and hot subdwarfs. They are a rare class of variable star due to their evolutionary history within interacting binary systems and the short timescales relative to their lifetime in which they are pulsationally… ▽ More Blue Large Amplitude Pulsators (BLAPs) are hot, subluminous stars undergoing rapid variability with periods of under 60 mins. They have been linked with the early stages of pre-white dwarfs and hot subdwarfs. They are a rare class of variable star due to their evolutionary history within interacting binary systems and the short timescales relative to their lifetime in which they are pulsationally unstable. All currently known BLAPs are relatively faint (15-19 mag) and are located in the Galactic plane. These stars have intrinsically blue colours but the large interstellar extinction in the Galactic plane prevents them from swift identification using colour-based selection criteria. In this paper, we correct the Gaia $G$-band apparent magnitude and $G_{\mathrm{BP}}-G_{\mathrm{RP}}$ colours of 89.6 million sources brighter than 19 mag in the Galactic plane with good quality photometry combined with supplementary all-sky data totalling 162.3 million sources. Selecting sources with colours consistent with the known population of BLAPs and performing a cross-match with the Zwicky Transient Facility (ZTF) DR3, we identify 98 short period candidate variables. Manual inspection of the period-folded light curves reveals 22 candidate BLAPs. Of these targets, 6 are consistent with the observed periods and light curves of the known BLAPs, 10 are within the theoretical period range of BLAPs and 6 are candidate high-gravity BLAPs. We present follow-up spectra of 21 of these candidate sources and propose to classify 1 of them as a BLAP, and tentatively assign an additional 8 of them as BLAPs for future population studies. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: 11 pages, 5 figures, 2 tables, accepted to MNRAS

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

Single-lens mass measurement in the high-magnification microlensing event Gaia19bld located in the Galactic disc

Authors: K. A. Rybicki, Ł. Wyrzykowski, E. Bachelet, A. Cassan, P. Zieliński, A. Gould, S. Calchi Novati, J. C. Yee, Y. -H. Ryu, M. Gromadzki, P. Mikołajczyk, N. Ihanec, K. Kruszyńska, F. -J. Hambsch, S. Zoła, S. J. Fossey, S. Awiphan, N. Nakharutai, F. Lewis, F. Olivares E., S. Hodgkin, A. Delgado, E. Breedt, D. L. Harrison, M. vanLeeuwen , et al. (44 additional authors not shown)

Abstract: We present the photometric analysis of Gaia19bld, a high-magnification ($A\approx60$) microlensing event located in the southern Galactic plane, which exhibited finite source and microlensing parallax effects. Due to a prompt detection by the Gaia satellite and the very high brightness of $I = 9.05~$mag at the peak, it was possible to collect a complete and unique set of multi-channel follow-up ob… ▽ More We present the photometric analysis of Gaia19bld, a high-magnification ($A\approx60$) microlensing event located in the southern Galactic plane, which exhibited finite source and microlensing parallax effects. Due to a prompt detection by the Gaia satellite and the very high brightness of $I = 9.05~$mag at the peak, it was possible to collect a complete and unique set of multi-channel follow-up observations, which allowed us to determine all parameters vital for the characterisation of the lens and the source in the microlensing event. Gaia19bld was discovered by the Gaia satellite and was subsequently intensively followed up with a network of ground-based observatories and the Spitzer Space Telescope. We collected multiple high-resolution spectra with Very Large Telescope (VLT)/X-Shooter to characterise the source star. The event was also observed with VLT Interferometer (VLTI)/PIONIER during the peak. Here we focus on the photometric observations and model the light curve composed of data from Gaia, Spitzer, and multiple optical, ground-based observatories. We find the best-fitting solution with parallax and finite source effects. We derived the limit on the luminosity of the lens based on the blended light model and spectroscopic distance. We compute the mass of the lens to be $1.13 \pm 0.03~M_{\odot}$ and derive its distance to be $5.52^{+0.35}_{-0.64}~\mathrm{kpc}$. The lens is likely a main sequence star, however its true nature has yet to be verified by future high-resolution observations. Our results are consistent with interferometric measurements of the angular Einstein radius, emphasising that interferometry can be a new channel for determining the masses of objects that would otherwise remain undetectable, including stellar-mass black holes. △ Less

Submitted 2 December, 2021; originally announced December 2021.

Comments: accepted to Astronomy&Astrophysics

Automated SpectroPhotometric Image REDuction (ASPIRED)

Authors: Marco C. Lam, Robert J. Smith, Iair Arcavi, Iain A. Steele, Josh Veitch-Michaelis, Lukasz Wyrzykowski

Abstract: We provide a suite of public open-source spectral data-reduction software to rapidly obtain scientific products from all forms of long-slit-like spectroscopic observations. Automated SpectroPhotometric REDuction (ASPIRED) is a Python-based spectral data-reduction toolkit. It is designed to be a general toolkit with high flexibility for users to refine and optimize their data-reduction routines for… ▽ More We provide a suite of public open-source spectral data-reduction software to rapidly obtain scientific products from all forms of long-slit-like spectroscopic observations. Automated SpectroPhotometric REDuction (ASPIRED) is a Python-based spectral data-reduction toolkit. It is designed to be a general toolkit with high flexibility for users to refine and optimize their data-reduction routines for the individual characteristics of their instruments. The default configuration is suitable for low-resolution long-slit spectrometers and provides a quick-look quality output. However, for repeatable science-ready reduced spectral data, some moderate one-time effort is necessary to modify the configuration. Fine-tuning and additional (pre)processing may be required to extend the reduction to systems with more complex setups. It is important to emphasize that although only a few parameters need updating, ensuring their correctness and suitability for generalization to the instrument can take time due to factors such as instrument stability. We compare some example spectra reduced with ASPIRED to published data processed with iraf-based and STARLINK-based pipelines, and find no loss in the quality of the final product. The Python-based, iraf-free ASPIRED can significantly ease the effort of an astronomer in constructing their own data-reduction workflow, enabling simpler solutions to data-reduction automation. This availability of near-real-time science-ready data will allow adaptive observing strategies, particularly important in, but not limited to, time-domain astronomy. △ Less

Submitted 14 June, 2023; v1 submitted 3 November, 2021; originally announced November 2021.

Comments: 23 pages, 10 figures, 1 table, 1 appendix

Journal ref: AJ, 166, 13 (2023)

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

An update on the development of ASPIRED

Authors: Marco C Lam, Robert J Smith, Josh Veitch-Michaelis, Iain A Steele, Paul R McWhirter

Abstract: We are reporting the updates in version 0.2.0 of the Automated SpectroPhotometric REDuction (ASPIRED) pipeline, designed for common use on different instruments. The default settings support many typical long-slit spectrometer configurations, whilst it also offers a flexible set of functions for users to refine and tailor-make their automated pipelines to an instrument's individual characteristics… ▽ More We are reporting the updates in version 0.2.0 of the Automated SpectroPhotometric REDuction (ASPIRED) pipeline, designed for common use on different instruments. The default settings support many typical long-slit spectrometer configurations, whilst it also offers a flexible set of functions for users to refine and tailor-make their automated pipelines to an instrument's individual characteristics. Such automation provides near real-time data reduction to allow adaptive observing strategies, which is particularly important in the Time Domain Astronomy. Over the course of last year, significant improvement was made in the internal data handling as well as data I/O, accuracy and repeatability in the wavelength calibration. △ Less

Submitted 7 December, 2020; originally announced December 2020.

Comments: 4 pages, 2 figures, conference proceeding article for the ADASS XXX to be published in ASP Conference Series

Automating Inference of Binary Microlensing Events with Neural Density Estimation

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

Abstract: Automated inference of binary microlensing events with traditional sampling-based algorithms such as MCMC has been hampered by the slowness of the physical forward model and the pathological likelihood surface. Current analysis of such events requires both expert knowledge and large-scale grid searches to locate the approximate solution as a prerequisite to MCMC posterior sampling. As the next gen… ▽ More Automated inference of binary microlensing events with traditional sampling-based algorithms such as MCMC has been hampered by the slowness of the physical forward model and the pathological likelihood surface. Current analysis of such events requires both expert knowledge and large-scale grid searches to locate the approximate solution as a prerequisite to MCMC posterior sampling. As the next generation, space-based microlensing survey with the Roman Space Observatory is expected to yield thousands of binary microlensing events, a new scalable and automated approach is desired. Here, we present an automated inference method based on neural density estimation (NDE). We show that the NDE trained on simulated Roman data not only produces fast, accurate, and precise posteriors but also captures expected posterior degeneracies. A hybrid NDE-MCMC framework can further be applied to produce the exact posterior. △ Less

Submitted 15 February, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: 7 pages, 1 figure. This article is superseded by arXiv:2102.05673. Accepted to the ML4PS workshop at NeurIPS 2020

Two New Nova Shells associated with V4362 Sagittarii and DO Aquilae

Authors: E. J. Harvey, M. P. Redman, P. Boumis, S. Akras, K. Fitzgerald, S. Dulaimi, S. C. Williams, M. J. Darnley, M. C. Lam, 1 M. Kopsacheilli, S. Derlopa

Abstract: A classical nova is an eruption on the surface of a white dwarf in an accreting binary system. The material ejected from the white dwarf surface generally forms an axisymmetric shell. The shaping mechanisms of nova shells are probes of the processes that take place at energy scales between planetary nebulae and supernova remnants. We report on the discovery of nova shells surrounding the post-nova… ▽ More A classical nova is an eruption on the surface of a white dwarf in an accreting binary system. The material ejected from the white dwarf surface generally forms an axisymmetric shell. The shaping mechanisms of nova shells are probes of the processes that take place at energy scales between planetary nebulae and supernova remnants. We report on the discovery of nova shells surrounding the post-nova systems V4362 Sagittarii (1994) and more limited observations of DO Aquilae (1925). Distance measurements of 0.5p/m1.4 kpc for V4362 Sgr and 6.7 p/m 3.5 kpc -0.2 for DO Aql are found based on the expansion parallax method. The growth rates are measured to be 0.07``/year for DO Aql and 0.32``/year for V4362 Sgr. A preliminary investigation into the ionisation structure of the nova shell associated with V4362 Sgr is presented. The observed ionisation structure of nova shells depends strongly on their morphology and the orientation of the central component towards the observer. X-ray, IR and UV observations as well as optical integral field unit spectroscopy are required to better understand these interesting objects. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 18 pages, 14 Figures, 9 Tables, accepted MNRAS

MSC Class: 85-11

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

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

Confirmation of Monoperiodicity Above 20 Seconds for Two Blue Large-Amplitude Pulsators

Authors: Paul Ross McWhirter, Marco C. Lam, Iain A. Steele

Abstract: Blue Large-Amplitude Pulsators (BLAPs) are a new class of pulsating variable star. They are located close to the hot subdwarf branch in the Hertzsprung-Russell diagram and have spectral classes of late O or early B. Stellar evolution models indicate that these stars are likely radially pulsating, driven by iron group opacity in their interiors. A number of variable stars with a similar driving mec… ▽ More Blue Large-Amplitude Pulsators (BLAPs) are a new class of pulsating variable star. They are located close to the hot subdwarf branch in the Hertzsprung-Russell diagram and have spectral classes of late O or early B. Stellar evolution models indicate that these stars are likely radially pulsating, driven by iron group opacity in their interiors. A number of variable stars with a similar driving mechanism exist near the hot subdwarf branch with multi-periodic oscillations caused by either pressure (p) or gravity (g) modes. No multi-periodic signals were detected in the OGLE discovery light curves since it would be difficult to detect short period signals associated with higher-order p modes with the OGLE cadence. Using the RISE instrument on the Liverpool Telescope, we produced high cadence light curves of two BLAPs, OGLE-BLAP-009 ($m_{\mathrm{v}}=15.65$ mag) and OGLE-BLAP-014 ($m_{\mathrm{v}}=16.79$ mag) using a $720$ nm longpass filter. Frequency analysis of these light curves identify a primary oscillation with a period of $31.935\pm0.0098$ mins and an amplitude from a Fourier series fit of $0.236$ mag for BLAP-009. The analysis of BLAP-014 identifies a period of $33.625\pm0.0214$ mins and an amplitude of $0.225$ mag. Analysis of the residual light curves reveals no additional short period variability down to an amplitude of $15.20\pm0.26$ mmag for BLAP-009 and $58.60\pm3.44$ mmag for BLAP-014 for minimum periods of $20$ s and $60$ s respectively. These results further confirm that the BLAPs are monoperiodic. △ Less

Submitted 5 June, 2020; originally announced June 2020.

Comments: 10 pages, 9 figures, 4 tables, accepted to MNRAS

First discovery of an ultra-cool white dwarf benchmark in common proper motion with an M dwarf

Authors: M. C. Lam, N. C. Hambly, N. Lodieu, S. Blouin, E. J. Harvey, R. J. Smith, M. C. Galvez-Ortiz, Z. H. Zhang

Abstract: Ultra-cool white dwarfs are among the oldest stellar remnants in the Universe. Their efficient gravitational settling and low effective temperatures are responsible for the smooth spectra they exhibit. For that reason, it is not possible to derive their radial velocities or to find the chemistry of the progenitors. The best that can be done is to infer such properties from associated sources, whic… ▽ More Ultra-cool white dwarfs are among the oldest stellar remnants in the Universe. Their efficient gravitational settling and low effective temperatures are responsible for the smooth spectra they exhibit. For that reason, it is not possible to derive their radial velocities or to find the chemistry of the progenitors. The best that can be done is to infer such properties from associated sources, which are coeval. The simplest form of such a system is a common proper motion pair where one star is an evolved degenerate and the other a main-sequence star. In this work, we present the discovery of the first of such a system, the M dwarf LHS 6328 and the ultra-cool white dwarf PSO J1801+625, from the Pan-STARRS 1 3$π$ survey and the $Gaia$ Data Release 2. Follow-up spectra were collected covering a usable wavelength range of 3750-24500 Å. Their spectra show that the white dwarf has an effective temperature of 3550K and surface gravity of $\log g = 7.45 \pm 0.13$ or $\log g = 7.49 \pm 0.13$ for a CO or He core, respectively, when compared against synthetic spectra of ultra-cool white dwarf atmosphere models. The system has slightly subsolar metallicity with -0.25 < [Fe/H] < 0.0, and a spatial velocity of (U, V, W) = ($-114.26 \pm 0.24$, $222.94 \pm 0.60$, $10.25 \pm 0.34$) km s $^{-1}$ the first radial velocity and metallicity measurements of an ultra-cool white dwarf. This makes it the first and only benchmark of its kind to date. △ Less

Submitted 7 April, 2020; v1 submitted 25 February, 2020; originally announced February 2020.

Comments: 11 pages, 4 figures, 5 tables. Some typographical changes

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 493, Issue 4, April 2020, Pages 6001-6010

The First Habitable Zone Earth-sized Planet from TESS. I: Validation of the TOI-700 System

Authors: Emily A. Gilbert, Thomas Barclay, Joshua E. Schlieder, Elisa V. Quintana, Benjamin J. Hord, Veselin B. Kostov, Eric D. Lopez, Jason F. Rowe, Kelsey Hoffman, Lucianne M. Walkowicz, Michele L. Silverstein, Joseph E. Rodriguez, Andrew Vanderburg, Gabrielle Suissa, Vladimir S. Airapetian, Matthew S. Clement, Sean N. Raymond, Andrew W. Mann, Ethan Kruse, Jack J. Lissauer, Knicole D. Colón, Ravi kumar Kopparapu, Laura Kreidberg, Sebastian Zieba, Karen A. Collins , et al. (70 additional authors not shown)

Abstract: We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follo… ▽ More We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system. △ Less

Submitted 10 July, 2020; v1 submitted 3 January, 2020; originally announced January 2020.

Comments: Accepted for publication in AJ (30 pages, 13 figures, 4 tables, 2 appendices)

Automated SpectroPhotometric Image REDuction (ASPIRED)

Authors: Marco C Lam, Robert J Smith, Iain A Steele

Abstract: We aim to provide a suite of publicly available spectral data reduction software to facilitate rapid scientific outcomes from time-domain observations. For time resolved observations, an automated pipeline frees astronomers from performance of the routine data analysis tasks to concentrate on interpretation, planning future observations and communication with international collaborators. The proje… ▽ More We aim to provide a suite of publicly available spectral data reduction software to facilitate rapid scientific outcomes from time-domain observations. For time resolved observations, an automated pipeline frees astronomers from performance of the routine data analysis tasks to concentrate on interpretation, planning future observations and communication with international collaborators. The project consists of two parts: data processing (ASPIRED) and a graphical user interface (gASPIRED).ASPIREDis written in python 3, and was intentionally developed as a self-consistent reduction pipeline with its own diagnostics and error handling. The pipeline can reduce2D spectral data from raw image to wavelength and flux calibrated 1D spectrum automatically without any user input.gASPIREDis a cross-platform software developedwithElectronon a single code base. It brings interactivity to the software with a well-maintained and user-friendly environment. △ Less

Submitted 12 December, 2019; originally announced December 2019.

Comments: 4 pages, 1 figure, 29th Astronomical Data Analysis Software & Systems (ADASS) conference

RASCAL: Towards automated spectral wavelength calibration

Authors: Josh Veitch-Michaelis, Marco C Lam

Abstract: Wavelength calibration is a routine and critical part of any spectral work-flow, but many astronomers still resort to matching detected peaks and emission lines by hand. We present RASCAL (RANSAC Assisted Spectral CALibration), a python library for automated wavelength calibration of astronomical spectrographs. RASCAL implements recent state-of-the-art methods for wavelength calibration and requir… ▽ More Wavelength calibration is a routine and critical part of any spectral work-flow, but many astronomers still resort to matching detected peaks and emission lines by hand. We present RASCAL (RANSAC Assisted Spectral CALibration), a python library for automated wavelength calibration of astronomical spectrographs. RASCAL implements recent state-of-the-art methods for wavelength calibration and requires minimal input from a user. In this paper we discuss the implementation of the library and apply it to real-world calibration spectra. △ Less

Submitted 12 December, 2019; originally announced December 2019.

Comments: 4 pages, 2 figures, 29th Astronomical Data Analysis Software & Systems (ADASS) conference

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

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

The L 98-59 System: Three Transiting, Terrestrial-Sized Planets Orbiting a Nearby M-dwarf

Authors: Veselin B. Kostov, Joshua E. Schlieder, Thomas Barclay, Elisa V. Quintana, Knicole D. Colon, Jonathan Brande, Karen A. Collins, Adina D. Feinstein, Samuel Hadden, Stephen R. Kane, Laura Kreidberg, Ethan Kruse, Christopher Lam, Elisabeth Matthews, Benjamin T. Montet, Francisco J. Pozuelos, Keivan G. Stassun, Jennifer G. Winters, George Ricker, Roland Vanderspek, David Latham, Sara Seager, Joshua Winn, Jon M. Jenkins, Dennis Afanasev , et al. (90 additional authors not shown)

Abstract: We report the Transiting Exoplanet Survey Satellite (TESS) discovery of three terrestrial-sized planets transiting L 98-59 (TOI-175, TIC 307210830) -- a bright M dwarf at a distance of 10.6 pc. Using the Gaia-measured distance and broad-band photometry we find that the host star is an M3 dwarf. Combined with the TESS transits from three sectors, the corresponding stellar parameters yield planet ra… ▽ More We report the Transiting Exoplanet Survey Satellite (TESS) discovery of three terrestrial-sized planets transiting L 98-59 (TOI-175, TIC 307210830) -- a bright M dwarf at a distance of 10.6 pc. Using the Gaia-measured distance and broad-band photometry we find that the host star is an M3 dwarf. Combined with the TESS transits from three sectors, the corresponding stellar parameters yield planet radii ranging from 0.8REarth to 1.6REarth. All three planets have short orbital periods, ranging from 2.25 to 7.45 days with the outer pair just wide of a 2:1 period resonance. Diagnostic tests produced by the TESS Data Validation Report and the vetting package DAVE rule out common false positive sources. These analyses, along with dedicated follow-up and the multiplicity of the system, lend confidence that the observed signals are caused by planets transiting L 98-59 and are not associated with other sources in the field. The L 98-59 system is interesting for a number of reasons: the host star is bright (V = 11.7 mag, K = 7.1 mag) and the planets are prime targets for further follow-up observations including precision radial-velocity mass measurements and future transit spectroscopy with the James Webb Space Telescope; the near resonant configuration makes the system a laboratory to study planetary system dynamical evolution; and three planets of relatively similar size in the same system present an opportunity to study terrestrial planets where other variables (age, metallicity, etc.) can be held constant. L 98-59 will be observed in 4 more TESS sectors, which will provide a wealth of information on the three currently known planets and have the potential to reveal additional planets in the system. △ Less

Submitted 28 May, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Comments: 27 pages, 22 figures, AJ accepted

Dark Matter Science in the Era of LSST

Authors: Keith Bechtol, Alex Drlica-Wagner, Kevork N. Abazajian, Muntazir Abidi, Susmita Adhikari, Yacine Ali-Haïmoud, James Annis, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Carlo Baccigalupi, Arka Banerjee, Nilanjan Banik, Charles Bennett, Florian Beutler, Simeon Bird, Simon Birrer, Rahul Biswas, Andrea Biviano, Jonathan Blazek, Kimberly K. Boddy, Ana Bonaca, Julian Borrill, Sownak Bose, Jo Bovy , et al. (155 additional authors not shown)

Abstract: Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We… ▽ More Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

Probing the Fundamental Nature of Dark Matter with the Large Synoptic Survey Telescope

Authors: Alex Drlica-Wagner, Yao-Yuan Mao, Susmita Adhikari, Robert Armstrong, Arka Banerjee, Nilanjan Banik, Keith Bechtol, Simeon Bird, Kimberly K. Boddy, Ana Bonaca, Jo Bovy, Matthew R. Buckley, Esra Bulbul, Chihway Chang, George Chapline, Johann Cohen-Tanugi, Alessandro Cuoco, Francis-Yan Cyr-Racine, William A. Dawson, Ana Díaz Rivero, Cora Dvorkin, Denis Erkal, Christopher D. Fassnacht, Juan García-Bellido, Maurizio Giannotti , et al. (75 additional authors not shown)

Abstract: Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We d… ▽ More Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ . △ Less

Submitted 24 April, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 96 pages, 22 figures, 1 table

Report number: FERMILAB-PUB-19-048-A-AE

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) Spectrographs

Authors: J. C. Wilson, F. R. Hearty, M. F. Skrutskie, S. R. Majewski, J. A. Holtzman, D. Eisenstein, J. Gunn, B. Blank, C. Henderson, S. Smee, M. Nelson, D. Nidever, J. Arns, R. Barkhouser, J. Barr, S. Beland, M. A. Bershady, M. R. Blanton, S. Brunner, A. Burton, L. Carey, M. Carr, J. P. Colque, J. Crane, G. J. Damke , et al. (64 additional authors not shown)

Abstract: We describe the design and performance of the near-infrared (1.51--1.70 micron), fiber-fed, multi-object (300 fibers), high resolution (R = lambda/delta lambda ~ 22,500) spectrograph built for the Apache Point Observatory Galactic Evolution Experiment (APOGEE). APOGEE is a survey of ~ 10^5 red giant stars that systematically sampled all Milky Way populations (bulge, disk, and halo) to study the Ga… ▽ More We describe the design and performance of the near-infrared (1.51--1.70 micron), fiber-fed, multi-object (300 fibers), high resolution (R = lambda/delta lambda ~ 22,500) spectrograph built for the Apache Point Observatory Galactic Evolution Experiment (APOGEE). APOGEE is a survey of ~ 10^5 red giant stars that systematically sampled all Milky Way populations (bulge, disk, and halo) to study the Galaxy's chemical and kinematical history. It was part of the Sloan Digital Sky Survey III (SDSS-III) from 2011 -- 2014 using the 2.5 m Sloan Foundation Telescope at Apache Point Observatory, New Mexico. The APOGEE-2 survey is now using the spectrograph as part of SDSS-IV, as well as a second spectrograph, a close copy of the first, operating at the 2.5 m du Pont Telescope at Las Campanas Observatory in Chile. Although several fiber-fed, multi-object, high resolution spectrographs have been built for visual wavelength spectroscopy, the APOGEE spectrograph is one of the first such instruments built for observations in the near-infrared. The instrument's successful development was enabled by several key innovations, including a "gang connector" to allow simultaneous connections of 300 fibers; hermetically sealed feedthroughs to allow fibers to pass through the cryostat wall continuously; the first cryogenically deployed mosaic volume phase holographic grating; and a large refractive camera that includes mono-crystalline silicon and fused silica elements with diameters as large as ~ 400 mm. This paper contains a comprehensive description of all aspects of the instrument including the fiber system, optics and opto-mechanics, detector arrays, mechanics and cryogenics, instrument control, calibration system, optical performance and stability, lessons learned, and design changes for the second instrument. △ Less

Submitted 3 February, 2019; originally announced February 2019.

Comments: 81 pages, 67 figures, PASP, accepted

The White Dwarf Luminosity Functions from the Pan-STARRS 1 3π Steradian Survey

Authors: Marco C. Lam, Nigel C. Hambly, Nicholas Rowell, Kenneth C. Chambers, Bertrand Goldman, Klaus W. Hodapp, Nick Kaiser, Rolf-Peter Kudritzki, Eugene A. Magnier, John L. Tonry, Richard J. Wainscoat, Christopher Waters

Abstract: A large sample of white dwarfs is selected by both proper motion and colours from the Pan-STARRS 1 3π Steradian Survey Processing Version 2 to construct the White Dwarf Luminosity Functions of the discs and halo in the solar neighbourhood. Four-parameter astrometric solutions were recomputed from the epoch data. The generalised maximum volume method is then used to calculate the density of the pop… ▽ More A large sample of white dwarfs is selected by both proper motion and colours from the Pan-STARRS 1 3π Steradian Survey Processing Version 2 to construct the White Dwarf Luminosity Functions of the discs and halo in the solar neighbourhood. Four-parameter astrometric solutions were recomputed from the epoch data. The generalised maximum volume method is then used to calculate the density of the populations. After removal of crowded areas near the Galactic plane and centre, the final sky area used by this work is 7.833 sr, which is 83% of the 3π sky and 62% of the whole sky. By dividing the sky using Voronoi tessellation, photometric and astrometric uncertainties are recomputed at each step of the integration to improve the accuracy of the maximum volume. Interstellar reddening is considered throughout the work. We find a disc-to-halo white dwarf ratio of about 100. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 18 pages, 13 figures, 2 tables in text, 4 tables in Appendix A, supplementary online materials are available on MNRAS. Accepted for publication on 3 Oct 2018

Laying the Groundwork for the Development of the Data Archive of the New Robotic Telescope

Authors: Marco C. Lam, Bovornpratch Vijarnwannaluk, Pathompong Butpan, Christopher M. Copperwheat, Andrzej S. Piascik, Utane Sawangwit, Robert J. Smith, Iain A. Steele

Abstract: The Liverpool Telescope has been in fully autonomous operation since 2004. The supporting data archive facility has largely been untouched. The data provision service has not been an issue although some modernisation of the system is desirable. This project is timely. Not only does it suit the upgrade of the current LT data archive, it is in line with the design phase of the New Robotic Telescope… ▽ More The Liverpool Telescope has been in fully autonomous operation since 2004. The supporting data archive facility has largely been untouched. The data provision service has not been an issue although some modernisation of the system is desirable. This project is timely. Not only does it suit the upgrade of the current LT data archive, it is in line with the design phase of the New Robotic Telescope which will be online in the early-2020s; and with the development of a new data archive facility for a range of telescopes at the National Astronomical Research Institute of Thailand. The Newton Fund enabled us to collaborate in designing a new versatile generic system that serves all purposes. In the end, we conclude that a single system would not meet the needs of all parties and only adopt similar front-ends while the back-ends are bespoke to our respective systems and data-flows. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 11 pages, 3 tables

Journal ref: Proc. SPIE 10707, Software and Cyberinfrastructure for Astronomy V, 1070721 (6 July 2018)

A machine learns to predict the stability of circumbinary planets

Authors: Christopher Lam, David Kipping

Abstract: Long-period circumbinary planets appear to be as common as those orbiting single stars and have been found to frequently have orbital radii just beyond the critical distance for dynamical stability. Assessing the stability is typically done either through N-body simulations or using the stability criterion first considered by Dvorak and later developed by Holman and Wiegert: a second-order polynom… ▽ More Long-period circumbinary planets appear to be as common as those orbiting single stars and have been found to frequently have orbital radii just beyond the critical distance for dynamical stability. Assessing the stability is typically done either through N-body simulations or using the stability criterion first considered by Dvorak and later developed by Holman and Wiegert: a second-order polynomial calibrated to broadly match numerical simulations. However, the polynomial is unable to capture islands of instability introduced by mean motion resonances, causing the accuracy of the criterion to approach that of a random coin-toss when close to the boundary. We show how a deep neural network (DNN) trained on N-body simulations generated with REBOUND is able to significantly improve stability predictions for circumbinary planets on initially coplanar, circular orbits. Specifically, we find that the accuracy of our DNN never drops below 86%, even when tightly surrounding the boundary of instability, and is fast enough to be practical for on-the-fly calls during likelihood evaluations typical of modern Bayesian inference. Our binary classifier DNN is made publicly available at https://github.com/CoolWorlds/orbital-stability △ Less

Submitted 8 February, 2018; v1 submitted 11 January, 2018; originally announced January 2018.

Comments: Accepted by MNRAS. Deep neural network and implementation can be found at https://github.com/CoolWorlds/orbital-stability

Laboratory demonstration of a cryogenic deformable mirror for wavefront correction of space-borne infrared telescopes

Authors: Aoi Takahashi, Keigo Enya, Kanae Haze, Hirokazu Kataza, Takayuki Kotani, Hideo Matsuhara, Tomohiro Kamiya, Tomoyasu Yamamuro, Paul Bierden, Steven Cornelissen, Charlie Lam, Michael Feinberg

Abstract: This paper demonstrates a cryogenic deformable mirror (DM) with 1,020 actuators based on micro-electrical mechanical systems (MEMS) technology. Cryogenic space-borne infrared telescopes can experience a wavefront error due to a figure error of their mirror surface, which makes the imaging performance worse. For on-orbit wavefront correction as one solution, we developed a MEMS-processed electro-st… ▽ More This paper demonstrates a cryogenic deformable mirror (DM) with 1,020 actuators based on micro-electrical mechanical systems (MEMS) technology. Cryogenic space-borne infrared telescopes can experience a wavefront error due to a figure error of their mirror surface, which makes the imaging performance worse. For on-orbit wavefront correction as one solution, we developed a MEMS-processed electro-static DM with a special surrounding structure for use under the cryogenic temperature. We conducted a laboratory demonstration of its operation in three cooling cycles between 5 K and 295 K. Using a laser interferometer, we detected the deformation corresponding to the applied voltages under the cryogenic temperature for the first time. The relationship between voltages and displacements was qualitatively expressed by the quadratic function, which is assumed based on the principle of electro-static DMs. We also found that it had a high operating repeatability of a few nm RMS and no significant hysteresis. Using the measured values of repeatability, we simulated the improvement of PSF by wavefront correction with our DM. These results show that our developed DM is effective in improving imaging performance and PSF contrast of space-borne infrared telescopes. △ Less

Submitted 29 September, 2017; originally announced October 2017.

Comments: 29 pages, 23 figures

Journal ref: Applied Optics, Vol.56, Issue 23, pp. 6694-6708 (2017)

Maximizing Survey Volume for Large-Area Multi-Epoch Surveys with Voronoi Tessellation

Authors: Marco C. Lam

Abstract: The survey volume of a proper motion-limited sample is typically much smaller than a magnitude-limited sample. This is because of the noisy astrometric measurements from detectors that are not dedicated for astrometric missions. In order to apply an empirical completeness correction, existing works limit the survey depth to the shallower parts of the sky that hamper the maximum potential of a surv… ▽ More The survey volume of a proper motion-limited sample is typically much smaller than a magnitude-limited sample. This is because of the noisy astrometric measurements from detectors that are not dedicated for astrometric missions. In order to apply an empirical completeness correction, existing works limit the survey depth to the shallower parts of the sky that hamper the maximum potential of a survey. The number of epoch of measurement is a discrete quantity that cannot be interpolated across the projected plane of observation, so that the survey properties change in discrete steps across the sky. This work proposes a method to dissect the survey into small parts with Voronoi tessellation using candidate objects as generating points, such that each part defines a `mini-survey' that has its own properties. Coupling with a maximum volume density estimator, the new method is demonstrated to be unbiased and recovered {\sim}20% more objects than the existing method in a mock catalogue of a white dwarf-only solar neighbourhood with Pan--STARRS 1-like characteristics. Towards the end of this work, we demonstrate one way to increase the tessellation resolution with artificial generating points, which would be useful for analysis of rare objects with small number counts. △ Less

Submitted 27 April, 2017; originally announced April 2017.

Comments: 10 pages, 2 tables, 8 figures, accepted by MNRAS on 10 April 2017

White Dwarf Luminosity Functions from the Pan{\ndash}STARRS1 3π Survey

Authors: Marco C. Lam

Abstract: White dwarfs are among the most common objects in the stellar halo; however, due to their low luminosity and low number density compared to the stars in the discs of the Milky Way, they are scarce in the observable volume. Hence, they are still poorly understood one hundred years after their discovery as relatively few have been observed. They are crucial to the understanding of several fundamenta… ▽ More White dwarfs are among the most common objects in the stellar halo; however, due to their low luminosity and low number density compared to the stars in the discs of the Milky Way, they are scarce in the observable volume. Hence, they are still poorly understood one hundred years after their discovery as relatively few have been observed. They are crucial to the understanding of several fundamental properties of the Galaxy {\mdash} the geometry, kinematics and star formation history, as well as to the study of the end-stage of stellar evolution for low- and intermediate-mass stars. White dwarfs were traditionally identified by their ultraviolet (UV) excess, however, if they have cooled for a long time, they become so faint in that part of the spectrum that they cannot be seen by the most sensitive modern detectors. Proper motion was then used as a means to identify white dwarf candidates, due to their relatively large space motions compared to other objects with the same colour. The use of proper motion as a selection criterion has proven effective and has yielded large samples of candidates with the SuperCOSMOS Sky Survey and Sloan Digital Sky Survey. In this work I will further increase the sample size with the Panchromatic Synoptic Telescope And Rapid Response System 1 (Pan{\ndash}STARRS1). (cont.) △ Less

Submitted 7 February, 2017; originally announced February 2017.

Comments: 173 pages, Thesis, University of Edinburgh

Get ready for Gaia: cool white dwarfs in common proper motion with Tycho stars

Authors: Nigel Hambly, Nick Rowell, Marco Lam

Abstract: We discuss the Gaia Data Release 1 (September 2016) and preliminary work on maximising the benefit for cool white dwarf (WD) science in advance of the full parallax catalogue which will appear around one year later in DR2. The Tycho catalogue is used in conjunction with the all-sky ground based astrometric/photometric SuperCOSMOS Sky Survey in order to identify candidate faint common proper motion… ▽ More We discuss the Gaia Data Release 1 (September 2016) and preliminary work on maximising the benefit for cool white dwarf (WD) science in advance of the full parallax catalogue which will appear around one year later in DR2. The Tycho catalogue is used in conjunction with the all-sky ground based astrometric/photometric SuperCOSMOS Sky Survey in order to identify candidate faint common proper motion objects to the Tycho stars. Gaia DR1 is supplemented by the Tycho-Gaia Astrometric Solution catalogue containing some 2 million parallaxes with Hipparcos-like precision for Tycho stars. While hotter, brighter WDs are present in Tycho, cooler examples are much rarer (if present at all) and CPM offers one method to infer precision distances for a statistically useful sample of these very faint WDs. △ Less

Submitted 7 February, 2017; originally announced February 2017.

Comments: 4 pages, 2 figures, EuroWD16: 20th European White Dwarf Workshop, to be published in PASP Vol. 509