skip to main content

Title: An improved test of the binary black hole hypothesis for quasars with double-peaked broad Balmer lines

Velocity offsets in the broad Balmer lines of quasars and their temporal variations serve as indirect evidence for bound supermassive black hole binaries (SBHBs) at sub-parsec separations. In this work, we test the SBHB hypothesis for 14 quasars with double-peaked broad emission lines using their long-term (14–41 yr) radial velocity curves. We improve on the previous work by (i) using elliptical instead of circular orbits for the SBHBs, (ii) adopting a statistical model for radial velocity jitter, (iii) employing a Markov chain Monte Carlo method to explore the orbital parameter space efficiently and build posterior distributions of physical parameters, and (iv) incorporating new observations. We determine empirically that jitter comprises approximately Gaussian distributed fluctuations about the smooth radial velocity curves that are larger than the measurement errors by factors of a few. We initially treat jitter by enlarging the effective error bars and then verify this approach via a variety of Gaussian process models for it. We find lower mass limits for the hypothesized SBHBs in the range 108–1011 M⊙. For seven objects, the SBHB scenario appears unlikely based on goodness-of-fit tests. For two additional objects, the minimum SBHB masses are unreasonably large (>1010 M⊙), strongly disfavouring the SBHB scenario. Using constraints on the orbital inclination angle (which requires some assumptions) makes the minimum masses of four more objects unreasonably large. We also cite physical and observational arguments against the SBHB hypothesis for nine objects. We conclude that the SBHB explanation is not the favoured explanation of double-peaked broad emission lines.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 1104-1126
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    CSS1603+19 is a cataclysmic variable (CV) with an orbital period of 81.96 min, near the minimal period of CVs. It is unusual in having a strong mid-infrared excess inconsistent with thermal emission from a brown dwarf companion. Here, we present time-resolved multiwavelength observations of this system. WISE photometry indicates that the mid-infrared excess displays a one-magnitude eclipsing-like variability during the orbit. We obtained near-infrared and optical spectroscopy using Gemini, MDM, and APO telescopes. Near-infrared spectra show possible cyclotron features indicating that the white dwarf has a magnetic field of about 5 MG. Optical and near-infrared spectra display double-peaked emission lines, with both components showing strong radial velocity variations during the orbital period and with the broad component leading the narrow component stably by about 0.2 of the orbital phase. We construct a physical model informed by existing observations of the system and determine that one component likely originates from the accretion column on to the magnetized white dwarf in synchronous rotation with the orbital motion and the other from the Roche overflow point. This allows us to constrain the masses of the binary components to be M1 > 0.24 M⊙ for the white dwarf accretor and M2 = 0.0644 ± 0.0074 M⊙ for the donor. We classify the system as an AM Herculis star, or a polar. It has likely completed its stint on the period gap, but has not yet gone through the period bounce.

    more » « less
  2. Context. Stellar evolution models are highly dependent on accurate mass estimates, especially for highly massive stars in the early stages of stellar evolution. The most direct method for obtaining model-independent stellar masses is derivation from the orbit of close binaries. Aims. Our aim was to derive the first astrometric plus radial velocity orbit solution for the single-lined spectroscopic binary star MWC 166 A, based on near-infrared interferometry over multiple epochs and ∼100 archival radial velocity measurements, and to derive fundamental stellar parameters from this orbit. A supplementary aim was to model the circumstellar activity in the system from K band spectral lines. Methods. The data used include interferometric observations from the VLTI instruments GRAVITY and PIONIER, as well as the MIRC-X instrument at the CHARA Array. We geometrically modelled the dust continuum to derive relative astrometry at 13 epochs, determine the orbital elements, and constrain individual stellar parameters at five different age estimates. We used the continuum models as a base to examine differential phases, visibilities, and closure phases over the Br γ and He  I emission lines in order to characterise the nature of the circumstellar emission. Results. Our orbit solution suggests a period of P  = 367.7 ± 0.1 d, approximately twice as long as found with previous radial velocity orbit fits. We derive a semi-major axis of 2.61 ± 0.04 au at d  = 990 ± 50 pc, an eccentricity of 0.498 ± 0.001, and an orbital inclination of 53.6 ± 0.3°. This allowed the component masses to be constrained to M 1  = 12.2 ± 2.2  M ⊙ and M 2  = 4.9 ± 0.5  M ⊙ . The line-emitting gas was found to be localised around the primary and is spatially resolved on scales of ∼11 stellar radii, where the spatial displacement between the line wings is consistent with a rotating disc. Conclusions. The large spatial extent and stable rotation axis orientation measured for the Br γ and He  I line emission are inconsistent with an origin in magnetospheric accretion or boundary-layer accretion, but indicate an ionised inner gas disc around this Herbig Be star. We observe line variability that could be explained either with generic line variability in a Herbig star disc or V/R variations in a decretion disc scenario. We have also constrained the age of the system, with relative flux ratios suggesting an age of ∼(7 ± 2)×10 5 yr, consistent with the system being composed of a main-sequence primary and a secondary still contracting towards the main-sequence stage. 
    more » « less
  3. Abstract

    In the first days of WZ Sge-type dwarf nova (DN) outbursts, the 2 : 1 resonance induces a spiral arm structure in the accretion disk, which is observed as early superhumps in optical light curves. We reports on our optical observations of an eclipsing WZ Sge-type DN PNV J00444033+4113068 during its 2021 superoutburst using the 3.8 m Seimei telescope and through the Variable Star Network collaboration. The eclipse analysis showed that its orbital period was 0.055425534(1) d. Our observations confirmed early superhumps with an amplitude of 0.7 mag, the largest amplitude among known WZ Sge-type DNe. More interestingly, its early superhumps became the reddest around their secondary minimum, whereas other WZ Sge-type DNe show the reddest color around the early superhump maximum. The spectrum around the peak of the outburst showed two double-peaked emission lines of He ii 4686 Å and Hα with a peak separation of ≥700 km s−1, supporting a very high-inclination system. With the early superhump mapping, the unique profile and color of the early superhump are successfully reproduced by an accretion disk with a vertically extended double arm structure. Therefore, a large amplitude and a unique color behavior of the early superhumps in PNV J00444033+4113068 can be explained by the 2 : 1 resonance model along with other WZ Sge-type DNe.

    more » « less
  4. Context. About a third of the hot subdwarfs of spectral type B (sdBs), which are mostly core-helium-burning objects on the extreme horizontal branch, are found in close binaries with cool, low-mass stellar, substellar, or white dwarf companions. They can show light variations due to different phenomena. Aims. Many hot subdwarfs now have space-based light curves with a high signal-to-noise ratio available. We used light curves from the Transiting Exoplanet Survey Satellite and the K2 space mission to look for more sdB binaries. Their light curves can be used to study the hot subdwarf primaries and their companions, and obtained orbital, atmospheric, and absolute parameters for those systems, when combined with other analysis methods. Methods. By classifying the light variations and combining these with the fit of the spectral energy distribution, the distance derived by the parallaxes obtained by Gaia , and the atmospheric parameters, mainly from the literature, we could derive the nature of the primaries and secondaries in 122 (75%) of the known sdB binaries and 82 newly found reflection effect systems. We derived absolute masses, radii, and luminosities for a total of 39 hot subdwarfs with cool, low-mass companions, as well 29 known and newly found sdBs with white dwarf companions. Results. The mass distribution of hot subdwarfs with cool, low-mass stellar and substellar companions, differs from those with white dwarf companions, implying they come from different populations. By comparing the period and minimum companion mass distributions, we find that the reflection effect systems all have M dwarf or brown dwarf companions, and that there seem to be several different populations of hot subdwarfs with white dwarf binaries – one with white dwarf minimum masses around 0.4  M ⊙ , one with longer periods and minimum companion masses up to 0.6  M ⊙ , and at the shortest period, another with white dwarf minimum masses around 0.8  M ⊙ . We also derive the first orbital period distribution for hot subdwarfs with cool, low-mass stellar or substellar systems selected from light variations instead of radial velocity variations. It shows a narrower period distribution, from 1.5 h to 35 h, compared to the distribution of hot subdwarfs with white dwarfs, which ranges from 1 h to 30 days. These period distributions can be used to constrain the previous common-envelope phase. 
    more » « less
  5. The hot nine-component system HD 93206, which contains a gravitationally bounded eclipsing Ac1+Ac2 binary ( P  = 5.9987 d) and a spectroscopic Aa1+Aa2 ( P  = 20.734 d) binary can provide important insights into the origin and evolution of massive stars. Using archival and new spectra, and a rich collection of ground-based and space photometric observations, we carried out a detailed study of this object. We provide a much improved description of both short orbits and a good estimate of the mutual period of both binaries of about 14 500 d (i.e. 40 years). For the first time, we detected weak lines of the fainter component of the 6.0 d eclipsing binary in the optical region of the spectrum, measured their radial velocities, and derived a mass ratio of M Ac2 / M Ac1  = 1.29, which is the opposite of what was estimated from the International Ultraviolet explorer (IUE) spectra. We confirm that the eclipsing subsystem Ac is semi-detached and is therefore in a phase of large-scale mass transfer between its components. The Roche-lobe filling and spectroscopically brighter component Ac1 is the less massive of the two and is eclipsed in the secondary minimum. We show that the bulk of the H α emission, so far believed to be associated with the eclipsing system, moves with the primary O9.7 I component Aa1 of the 20.73 d spectroscopic binary. However, the weak emission in the higher Balmer lines seems to be associated with the accretion disc around component Ac2. We demonstrate that accurate masses and other basic physical properties including the distance of this unique system can be obtained but require a more sophisticated modelling. A first step in this direction is presented in the accompanying Paper II (Brož et al.). 
    more » « less