We report a candidate centi-parsec supermassive black hole binary (SMBHB) in the radio-quiet quasar SDSS J0159+0105 at z= 0.217. With a modified Lomb-Scargle code (GLSdeDRW) and auto-correlation ...analysis, we detect two significant (at P> 99%) periodic signals at ~741 day and ~1500 day from the 8.1 yr Catalina V-band light curve of this quasar. The period ratio, which is close to 1:2, is typical of a black hole binary system with a mass ratio of 0.05 < q< 0.8 according to recent numerical simulations. SDSS J0159+0105 has two SDSS spectroscopic observations separated by ~10 yr. There is a significant change in the broad Hbeta profile between the two epochs, which can be explained by a single broad-line region (BLR) around the binary system illuminated by the aforementioned mini-disks, or a stream of gas flowing from the circumbinary disk to one of the SMBHs. From the single BLR assumption and the orbital period t sub(orb)~ 1500 day, we estimate the total virial masses of M sub(SMBHB)~ 1.3 x 10 super(8)M sub(middot in circle), the average distances of BLR of ~0.04 pc (~50 lt-day, with + or -0.3 dex uncertainty), and an SMBHB separation of d= (0.01 pc)M super(1/3) sub(8.tot )(T sub(rest)/3.3 yr) super(2/3)~ 0.0.13 pc (15 lt-day). Based on analytical work, the postulated circumbinary disk has an inner radius of 2d= 0.026 pc (30 lt-day). SDSS J0159+0105 also displays unusual spectral energy distribution. The unique properties of SDSS J0159+0105 are consistent with it being a centi-parsec SMBHB.
With the coming data deluge from synoptic surveys, there is a need for frameworks that can quickly and automatically produce calibrated classification probabilities for newly observed variables based ...on small numbers of time-series measurements. In this paper, we introduce a methodology for variable-star classification, drawing from modern machine-learning techniques. We describe how to homogenize the information gleaned from light curves by selection and computation of real-numbered metrics (features), detail methods to robustly estimate periodic features, introduce tree-ensemble methods for accurate variable-star classification, and show how to rigorously evaluate a classifier using cross validation. On a 25-class data set of 1542 well-studied variable stars, we achieve a 22.8% error rate using the random forest (RF) classifier; this represents a 24% improvement over the best previous classifier on these data. This methodology is effective for identifying samples of specific science classes: for pulsational variables used in Milky Way tomography we obtain a discovery efficiency of 98.2% and for eclipsing systems we find an efficiency of 99.1%, both at 95% purity. The RF classifier is superior to other methods in terms of accuracy, speed, and relative immunity to irrelevant features; the RF can also be used to estimate the importance of each feature in classification. Additionally, we present the first astronomical use of hierarchical classification methods to incorporate a known class taxonomy in the classifier, which reduces the catastrophic error rate from 8% to 7.8%. Excluding low-amplitude sources, the overall error rate improves to 14%, with a catastrophic error rate of 3.5%.
Long-duration gamma-ray bursts (GRBs) are powerful cosmic explosions, signaling the death of massive stars. Among them, GRB 221009A is by far the brightest burst ever observed. Because of its ...enormous energy (
≈ 10
erg) and proximity (
≈ 0.15), GRB 221009A is an exceptionally rare event that pushes the limits of our theories. We present multiwavelength observations covering the first 3 months of its afterglow evolution. The x-ray brightness decays as a power law with slope ≈
, which is not consistent with standard predictions for jetted emission. We attribute this behavior to a shallow energy profile of the relativistic jet. A similar trend is observed in other energetic GRBs, suggesting that the most extreme explosions may be powered by structured jets launched by a common central engine.
We report the discovery by the Swift hard X-ray monitor of the transient source Swift J2058.4+0516 (Sw J2058+05). Our multi-wavelength follow-up campaign uncovered a long-lived (duration > ~ months), ...luminous X-ray (L sub(X,iso) approximate 3 x 10 super(47) erg s super(-1)) and radio (nuL sub(nu,iso) approximate 10 super(42) erg s super(-1)) counterpart. The associated optical emission, however, from which we measure a redshift of 1.1853, is relatively faint, and this is not due to a large amount of dust extinction in the host galaxy. Based on numerous similarities with the recently discovered GRB 110328A/Swift J164449.3+573451 (Sw J1644+57), we suggest that Sw J2058+05 may be the second member of a new class of relativistic outbursts resulting from the tidal disruption of a star by a supermassive black hole. If so, the relative rarity of these sources (compared with the expected rate of tidal disruptions) implies that either these outflows are extremely narrowly collimated (straighttheta < 1degrees) or only a small fraction of tidal disruptions generate relativistic ejecta. Analogous to the case of long-duration gamma-ray bursts and core-collapse supernovae, we speculate that rapid spin of the black hole may be a necessary condition to generate the relativistic component. Alternatively, if powered by gas accretion (i.e., an active galactic nucleus (AGN)), Sw J2058+05 would seem to represent a new mode of variability in these sources, as the observed properties appear largely inconsistent with known classes of AGNs capable of generating relativistic jets (blazars, narrow-line Seyfert 1 galaxies).
With growing data volumes from synoptic surveys, astronomers necessarily must become more abstracted from the discovery and introspection processes. Given the scarcity of follow-up resources, there ...is a particularly sharp onus on the frameworks that replace these human roles to provide accurate and well-calibrated probabilistic classification catalogs. Such catalogs inform the subsequent follow-up, allowing consumers to optimize the selection of specific sources for further study and permitting rigorous treatment of classification purities and efficiencies for population studies. Here, we describe a process to produce a probabilistic classification catalog of variability with machine learning from a multi-epoch photometric survey. In addition to producing accurate classifications, we show how to estimate calibrated class probabilities and motivate the importance of probability calibration. We also introduce a methodology for feature-based anomaly detection, which allows discovery of objects in the survey that do not fit within the predefined class taxonomy. Finally, we apply these methods to sources observed by the All-Sky Automated Survey (ASAS), and release the Machine-learned ASAS Classification Catalog (MACC), a 28 class probabilistic classification catalog of 50,124 ASAS sources in the ASAS Catalog of Variable Stars. We estimate that MACC achieves a sub-20% classification error rate and demonstrate that the class posterior probabilities are reasonably calibrated. MACC classifications compare favorably to the classifications of several previous domain-specific ASAS papers and to the ASAS Catalog of Variable Stars, which had classified only 24% of those sources into one of 12 science classes.
We calculate durations and spectral parameters for 218 Swift bursts detected by the BAT instrument between and including gamma-ray bursts (GRBs) 041220 and 070509, including 77 events with measured ...redshifts. Incorporating prior knowledge into the spectral fits, we are able to measure the characteristic unk spectral peak energy unk and the isotropic equivalent energy unk (1 10 super(4) keV) for all events. This complete and rather extensive catalog, analyzed with a unified methodology, allows us to address the persistence and origin of high-energy correlations suggested in pre-Swift observations. We find that the unk correlation is present in the Swift sample; however, the best-fit power-law relation is inconsistent with the best-fit pre-Swift relation at >5 a significance. It has a factor unk2 larger intrinsic scatter, after accounting for large errors on unk. A large fraction of the Swift events are hard and subluminous relative to (and inconsistent with) the pre-Swift relation, in agreement with indications from BATSE GRBs without redshift. Moreover, we determine an experimental threshold for the BAT detector and show how the unk correlation arises artificially due to partial correlation with the threshold. We show that pre-Swift correlations found by Amati et al., Yonetoku et al., and Firmani et al., and independently by others are likely unrelated to the physical properties of GRBs and are likely useless for tests of cosmology. Also, an explanation for these correlations in terms of a detector threshold provides a natural and quantitative explanation for why short-duration GRBs and events at low redshift tend to be outliers to the correlations.
We have recently witnessed the first multi-messenger detection of colliding neutron stars through gravitational waves (GWs) and electromagnetic (EM) waves (GW 170817) thanks to the joint efforts of ...LIGO/Virgo and Space/Ground-based telescopes. In this paper, we report on the RATIR follow-up observation strategies and show the results for the trigger G194575. This trigger is not of astrophysical interest; however, it is of great interest to the robust design of a follow-up engine to explore large sky-error regions. We discuss the development of an image-subtraction pipeline for the six-color, optical/NIR imaging camera RATIR. Considering a two-band (i and r) campaign in the fall of 2015, we find that the requirement of simultaneous detection in both bands leads to a factor ∼10 reduction in false alarm rate, which can be further reduced using additional bands. We also show that the performance of our proposed algorithm is robust to fluctuating observing conditions, maintaining a low false alarm rate with a modest decrease in system efficiency that can be overcome utilizing repeat visits. Expanding our pipeline to search for either optical or NIR detections (three or more bands), considering separately the optical riZ and NIR YJH bands, should result in a false alarm rate 1% and an efficiency 90%. RATIR's simultaneous optical/NIR observations are expected to yield about one candidate transient in the vast 100 deg2 LIGO error region for prioritized follow-up with larger aperture telescopes.
We show that the X-ray and Y-ray spectra of Swift GRBs and their afterglows are consistent with the emission characteristic of an expanding, relativistic fireball. The classical afterglow due to the ...impact of the fireball on the external medium is often not observed until 1 to several hours after the GRB. Focusing on GRBs 061121, 060614, and 060124, but generalizing to the full (>50 Ms XRT exposure) Swift sample up to and including GRB 061210, we show that the early emission in >90% of early afterglows has a characteristic vF sub(v) spectral energy E sub(peak), which likely evolves from the Y-rays through the soft X-ray band on timesoales of 10 super(2)-10 super(4) s after the GRB. The observed spectra are strongly curved when plotted with logarithmic axes and have often been incorrectly fitted in other studies with a time-varying soft X-ray absorption. The spectral evolution inferred from fitting instead models used to fit GRBs demonstrates a common evolution-a power-law hardness-intensity correlation and hard-to-soft evolution-for GRBs and the early X-ray afterglows and X-ray flares. Combined with studies of short-timescale variability, our findings indicate a central engine active for longer than previously suspected. The GRB spectra are observed to become very soft at late times due to an intrinsic spectral evolution and due to the surprising faintness of some afterglows. We discuss models for the early X-ray emission.
Several correlations among parameters derived from modeling the high-energy properties of gamma-ray bursts (GRBs) have been reported. We show that well known examples of these have common features ...indicative of strong contamination by selection effects. We focus here on the impact of detector threshold truncation on the spectral peak versus isotropic equivalent energy release (E pk-E iso) relation, extended to a large sample of 218 Swift and 56 HETE-2 GRBs with and without measured redshift. The existence of faint Swift events missing from pre-Swift surveys calls into question inferences based on pre-Swift surveys which must be subject to complicated incompleteness effects. We demonstrate a generalized method for treating data truncation in correlation analyses and apply this method to Swift and pre-Swift data. Also, we show that the E pk-E g ('Ghirlanda') correlation is effectively independent of the GRB redshifts, which suggests that its existence has little to do with intrinsic physics. We suggest that a physically based correlation, manifest observationally, must show significantly reduced scatter in the rest frame relative to the observer frame and must not persist if the assumed redshifts are scattered. As with the E pk-E g correlation, we find that the pre-Swift, bright GRB E pk-E iso correlation of Amati does not rigorously satisfy these conditions.
ABSTRACT
We present analytical and numerical models of the bright long GRB 210822A at z = 1.736. The intrinsic extreme brightness exhibited in the optical, which is very similar to other bright GRBs ...(e.g. GRBs 080319B, 130427A, 160625A 190114C, and 221009A), makes GRB 210822A an ideal case for studying the evolution of this particular kind of GRB. We use optical data from the RATIR instrument starting at T + 315.9 s, with publicly available optical data from other ground-based observatories, as well as Swift/UVOT, and X-ray data from the Swift/XRT instrument. The temporal profiles and spectral properties during the late stages align consistently with the conventional forward shock model, complemented by a reverse shock element that dominates optical emissions during the initial phases (T < 300 s). Furthermore, we observe a break at T = 80 000 s that we interpreted as evidence of a jet break, which constrains the opening angle to be about θj = (3–5) degrees. Finally, we apply a machine-learning technique to model the multiwavelength light curve of GRB 210822A using the afterglowpy library. We estimate the angle of sight θobs = (6.4 ± 0.1) × 10−1 degrees, the energy E0 = (7.9 ± 1.6) × 1053 erg, the electron index p = 2.54 ± 0.10, the thermal energy fraction in electrons ϵe = (4.63 ± 0.91) × 10−5 and in the magnetic field ϵB = (8.66 ± 1.01) × 10−6, the efficiency χ = 0.89 ± 0.01, and the density of the surrounding medium n0 = 0.85 ± 0.01 cm−3.