skip to main content


Search for: All records

Creators/Authors contains: "Graham, Matthew J."

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract

    Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar’s brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described by a stochastic differential equation (SDE). Upcoming wide-field telescopes such as the Rubin Observatory Legacy Survey of Space and Time (LSST) are expected to observe tens of millions of AGN in multiple filters over a ten year period, so there is a need for efficient and automated modeling techniques that can handle the large volume of data. Latent SDEs are machine learning models well suited for modeling quasar variability, as they can explicitly capture the underlying stochastic dynamics. In this work, we adapt latent SDEs to jointly reconstruct multivariate quasar light curves and infer their physical properties such as the black hole mass, inclination angle, and temperature slope. Our model is trained on realistic simulations of LSST ten year quasar light curves, and we demonstrate its ability to reconstruct quasar light curves even in the presence of long seasonal gaps and irregular sampling across different bands, outperforming a multioutput Gaussian process regression baseline. Our method has the potential to provide a deeper understanding of the physical properties of quasars and is applicable to a wide range of other multivariate time series with missing data and irregular sampling.

     
    more » « less
  2. ABSTRACT

    The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6σ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ∼10 per cent of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares.

     
    more » « less
  3. ABSTRACT

    Cataclysmic variables (CVs) that have evolved past the period minimum during their lifetimes are predicted to be systems with a brown dwarf donor. While population synthesis models predict that around 40–70 per cent of the Galactic CVs are post-period minimum systems referred to as ‘period bouncers’, only a few dozen confirmed systems are known. We report the study and characterization of a new eclipsing CV, SRGeJ041130.3+685350 (SRGeJ0411), discovered from a joint SRG/eROSITA and ZTF programme. The optical spectrum of SRGeJ0411 shows prominent hydrogen and helium emission lines, typical for CVs. We obtained optical high-speed photometry to confirm the eclipse of SRGeJ0411 and determine the orbital period to be Porb ≈ 97.530 min. The spectral energy distribution suggests that the donor has an effective temperature of ≲ 1800 K. We constrain the donor mass with the period–density relationship for Roche lobe-filling stars and find that Mdonor ≲ 0.04 M⊙. The binary parameters are consistent with evolutionary models for post-period minimum CVs, suggesting that SRGeJ0411 is a new period bouncer. The optical emission lines of SRGeJ0411 are single-peaked despite the system being eclipsing, which is typically only seen due to stream-fed accretion in polars. X-ray spectroscopy hints that the white dwarf in SRGeJ0411 could be magnetic, but verifying the magnetic nature of SRGeJ0411 requires further investigation. The lack of optical outbursts has made SRGeJ0411 elusive in previous surveys, and joint X-ray and optical surveys highlight the potential for discovering similar systems in the near future.

     
    more » « less
  4. Abstract

    SX Phoenicis (SXP) variables are short-period pulsating stars that exhibit a period–luminosity (PL) relation. We derived thegri-band PL and extinction-free period–Wesenheit (PW) relations, as well as the period-color and reddening-free period-Q-index relations for 47 SXP stars located in 21 globular clusters, using the optical light curves taken from Zwicky Transient Facility. These empirical relations were derived for the first time in thegrifilters except for theg-band PL relation. We used ourgi-band PL and PW relations to derive a distance modulus to Crater II dwarf spheroidal which hosts one SXP variable. Assuming that the fundamental and first-overtone pulsation mode for the SXP variable in Crater II, we found distance moduli of 20.03 ± 0.23 mag and 20.37 ± 0.24 mag, respectively, using the PW relation, where the latter is in excellent agreement with independent RR Lyrae based distance to Crater II dwarf galaxy.

     
    more » « less
  5. ABSTRACT

    The scarce optical variability studies in spectrally classified Type 2 active galactic nuclei (AGNs) have led to the discovery of anomalous objects that are incompatible with the simplest unified models (UMs). This paper focuses on the exploration of different variability features that allow to distinguish between obscured, Type 2 AGNs and the variable, unobscured Type 1s. We analyse systematically the Zwicky Transient Facility, 2.5-yr-long light curves of ∼15 000 AGNs from the Sloan Digital Sky Survey Data Release 16, which are generally considered Type 2s due to the absence of strong broad emission lines (BELs). Consistent with the expectations from the UM, the variability features are distributed differently for distinct populations, with spectrally classified weak Type 1s showing one order of magnitude larger variances than the Type 2s. We find that the parameters given by the damped random walk model lead to broader H α equivalent width for objects with τg > 16 d and long-term structure function SF∞, g > 0.07 mag. By limiting the variability features, we find that ∼11 per cent of Type 2 sources show evidence for optical variations. A detailed spectral analysis of the most variable sources (∼1 per cent of the Type 2 sample) leads to the discovery of misclassified Type 1s with weak BELs and changing-state candidates. This work presents one of the largest systematic investigations of Type 2 AGN optical variability to date, in preparation for future large photometric surveys.

     
    more » « less
  6. Abstract We report multiwavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka “Barbie”) identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak-observed pseudo-bolometric luminosity of log( L max / [ erg s − 1 ] ) = 45.7 from slowly fading ztf- g and ztf- r light curves spanning over 1000 observer-frame days. The host galaxy is not detected in archival Pan-STARRS observations ( g > 23.3 mag), implying a lower limit to the outburst amplitude of more than 5 mag relative to the quiescent host galaxy. Optical spectra exhibit strong emission lines with narrow cores from the H Balmer series and ultraviolet semi-forbidden lines of Si iii ] λ 1892, C iii ] λ 1909, and  C ii ] λ 2325. Typical nebular lines in Active Galactic Nucleus (AGN) spectra from ions such as [O ii ] and [O iii ] are not detected. These spectral features, along with the smooth light curve that is unlike most AGN flaring activity and the luminosity that exceeds any observed or theorized supernova, lead us to conclude that AT 2021lwx is most likely an extreme tidal disruption event (TDE). Modeling of ZTF photometry with MOSFiT suggests that the TDE was between a ≈14 M ⊙ star and a supermassive black hole of mass M BH ∼ 10 8 M ⊙ . Continued monitoring of the still-evolving light curve along with deep imaging of the field once AT 2021lwx has faded can test this hypothesis and potentially detect the host galaxy. 
    more » « less
    Free, publicly-accessible full text available May 1, 2024
  7. Abstract We present the first gri -band period–luminosity (PL) and period–Wesenheit (PW) relations for 37 Type II Cepheids (TIICs) located in 18 globular clusters based on photometric data from the Zwicky Transient Facility. We also updated BVIJHK -band absolute magnitudes for 58 TIICs in 24 globular clusters using the latest homogeneous distances to the globular clusters. The slopes of g / r / i - and B / V / I -band PL relations are found to be statistically consistent when using the same sample of distance and reddening. We employed the calibration of ri -band PL/PW relations in globular clusters to estimate a distance to M31 based on a sample of ∼270 TIICs from the PAndromeda project. The distance modulus to M31, obtained using calibrated ri -band PW relation, agrees well with the recent determination based on classical Cepheids. However, distance moduli derived using the calibrated r - and i -band PL relations are systematically smaller by ∼0.2 mag, suggesting there are possible additional systematic errors on the PL relations. Finally, we also derive the period–color (PC) relations and for the first time the period–Q-index (PQ) relations, where the Q -index is reddening free, for our sample of TIICs. The PC relations based on ( r − i ) and near-infrared colors and the PQ relations are found to be relatively independent of the pulsation periods. 
    more » « less
  8. Abstract Magnetic cataclysmic variables (CVs) are luminous Galactic X-ray sources, which have been difficult to find in purely optical surveys due to their lack of outburst behavior. The eROSITA telescope on board the Spektr-RG mission is conducting an all-sky X-ray survey and recently released the public eROSITA Final Equatorial Depth Survey (eFEDS) catalog. We crossmatched the eFEDS catalog with photometry from the Zwicky Transient Facility and discovered two new magnetic CVs. We obtained high-cadence optical photometry and phase-resolved spectroscopy for each magnetic CV candidate and found them both to be polars. Among the newly discovered magnetic CVs is eFEDS J085037.2+044359/ZTFJ0850+0443, an eclipsing polar with orbital period P orb = 1.72 hr and WD mass M WD = 0.81 ± 0.08 M ⊙ . We suggest that eFEDS J085037.2+044359/ZTFJ0850+0443 is a low magnetic field strength polar, with B WD ≲ 10 MG. We also discovered a non-eclipsing polar, eFEDS J092614.1+010558/ZTFJ0926+0105, with orbital period P orb = 1.47 hr and magnetic field strength B WD = 36–42 MG. 
    more » « less
  9. Abstract

    We present the firstgri-band period–luminosity (PL) and period–Wesenheit (PW) relations for the fundamental mode anomalous Cepheids. These PL and PW relations were derived from a combined sample of five anomalous Cepheids in globular cluster M92 and the Large Magellanic Cloud, both of which have distance accurate to ∼1% available from literature. Ourg-band PL relation is similar to theB-band PL relation as reported in previous study. We applied our PL and PW relations to anomalous Cepheids discovered in dwarf galaxy Crater II, and found a larger but consistent distance modulus than the recent measurements based on RR Lyrae. Our calibrations ofgri-band PL and PW relations, even though less precise due to small number of anomalous Cepheids, will be useful for distance measurements to dwarf galaxies.

     
    more » « less
  10. Abstract

    One of the open questions following the discovery of GW170817 is whether neutron star (NS) mergers are the only astrophysical sites capable of producingr-process elements. Simulations have shown that 0.01–0.1Mofr-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both NS mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature ofr-process nucleosynthesis in the binary NS merger GW170817 was its long-lasting near-infrared (NIR) emission, thus motivating a systematic photometric study of the light curves of broad-lined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL—including 18 observed with the Zwicky Transient Facility and 7 from the literature—in the optical/NIR bands to determine what quantity ofr-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account forr-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on ther-process mass for these SNe. We also perform independent light curve fits to models without ther-process. We find that ther-process-free models are a better fit to the light curves of the objects in our sample. Thus, we find no compelling evidence ofr-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities ofr-process ejecta mass or indicate whether all collapsars are completely devoid ofr-process nucleosynthesis.

     
    more » « less