Extreme Mass-ratio Binaries Research Articles

Effective Field Theory for Extreme Mass Ratio Binaries.

We derive an effective field theory describing a pair of gravitationally interacting point particles in an expansion in their mass ratio, also known as the self-force (SF) expansion. The 0SF dynamics are trivially obtained to all orders in Newton's constant by the geodesic motion of the light body in a Schwarzschild background encoding the gravitational field of the heavy body. The corrections at 1SF and higher are generated by perturbations about this configuration-that is, the geodesic deviation of the light body and the fluctuation graviton-but crucially supplemented by an operator describing the recoil of the heavy body as it interacts with the smaller companion. Using this formalism we compute new results at third post-Minkowskian order for the conservative dynamics of a system of gravitationally interacting massive particles coupled to a set of additional scalar and vector fields.

Periodically modulated FRB as extreme mass ratio binaries

ABSTRACT The activity of at least one repeating fast radio burst (FRB) source is periodically modulated. If this modulation is the result of precession of the rotation axis and throat of an accretion disc around a black hole, driven by a companion that is also the source of accreted mass, then it may be possible to constrain the mass of the black hole. The dynamics is analogous to that of superorbital periods in ordinary mass-transfer binaries in which the accreting object may be a stellar-mass black hole, a neutron star or a white dwarf, but in the FRB source it may be an intermediate-mass black hole. In a semidetached (mass-transferring) binary, the orbital period is nearly proportional to the −1/2 power of the mean density of the mass-losing star and nearly independent of the mass of the primary, but the ratio of precessional to orbital periods scales approximately as the −2/3 power of the mass ratio for small mass ratios (massive accretors). Assuming a value for the secondary’s density and identifying the observed modulation period as a disc precession period would determine the mass ratio and the mass of the black hole. This model and magnetar-SNR (supernova remnant) models make distinguishable predictions of the evolution of the rotation measure that may soon be tested in FRB 121102.

Evidence for saturated and disrupted magnetic braking from samples of detached close binaries with M and K dwarfs

Context. Recent observations of close detached eclipsing M and K dwarf binaries have provided substantial support for magnetic saturation when stars rotate sufficiently fast, leading to a magnetic braking (MB) torque proportional to the spin of the star. Aims. We investigated here how strong MB torques need to be to reproduce the observationally inferred relative numbers of white dwarf plus M dwarf post-common-envelope binaries under the assumption of magnetic saturation. Methods. We carried out binary population simulations with the BSE code adopting empirically derived inter-correlated main-sequence binary distributions as initial binary populations and compared the simulation outcomes with observations. Results. We found that the dearth of extreme mass ratio binaries in the inter-correlated initial distributions is key to reproduce the large fraction of post-common-envelope binaries hosting low-mass M dwarfs (∼0.1 − 0.2 M⊙). In addition, orbital angular momentum loss rates due to MB should be high for M dwarfs with radiative cores and orders of magnitude smaller for fully convective stars to explain the observed dramatic change of the fraction of short-period binaries at the fully convective boundary. Conclusions. We conclude that saturated but disrupted, that is, dropping drastically at the fully convective boundary, MB can explain the observations of both close main-sequence binaries containing M and K dwarfs and post-common-envelope binaries. Whether a similar prescription can explain the spin down rates of single stars and of binaries containing more massive stars needs to be tested.

The critical mass ratio for W UMa-type contact binary systems

Contact binaries are close binary systems in which both components fill their inner Roche lobes so that the stars are in direct contact, and in potential mass and energy exchange. The most common such systems of low mass are the so-called W UMa-type. In the last few years, there has been a growing interest of the astronomical community in stellar mergers, primarily due to the detection of gravitational waves (mergers of black holes and neutron stars), but also because of an alternative model for the type Ia supernovae (merger of two white dwarfs), which are again particularly important in cosmology where they played a significant role in the discovery of dark energy and the accelerated expansion of the Universe. In that sense, contact systems of W UMa type with extremely low mass ratio are especially interesting because there are indications that, in their case too, stars can merge and possibly form fast-rotating stars such as FC Com stars and the blue-stragglers, and (luminous) red novae such as V1309 Sco. Namely, the previous theoretical research has shown that in the cases when the orbital angular momentum of the system is only about three times larger than the rotational angular momentum of the primary, a tidal Darwin's instability occurs, the components can no longer remain in synchronous rotation, orbit continue to shrink fast, and they finally merge into a single star. The above stability condition for contact systems can be linked to a specific critical mass ratio below which we expect a system to be unstable. We give an overview of this condition and show how it can be used to identify potential mergers. Finally, we discuss a number of known extreme mass ratio binaries from the literature and consider prospects for future research on this topic.

Absolute parameters of four W UMa stars with extreme low mass ratios

Photometric observations of four extreme low mass ratio W UMa binaries with wide total eclipses, and their absolute parameters are presented. The main results are as follows: (i) ASAS J015429+2042.8 is of A subtype while V755 Cep, VESPA V22 and GSC 02800-01387 are of W subtype; (ii) The mass ratios are close to the lower mass-ratio limit and can be classified as Extreme Low Mass Ratio Binaries (ELMRB). (iii) The fill-out factors of ASAS J015429+2042.8 and VESPA V22 are 0.16 and 0.21 respectively, the one of GSC 02800-01387 is 0.63 while the one of V755 Cep is 0.83; (iv) The smaller in size (and less massive) components of all targets are located above the MS band of the evolutional diagrams while the bigger (and more massive) components are located around ZAMS; (v) Our ELMRB systems fall into the stability zone with Jspin/Jorb lower than the Darwin’s instability ∼1/3.

Divergences in gravitational-wave emission and absorption from extreme mass ratio binaries

A powerful technique to calculate gravitational radiation from binary systems involves a perturbative expansion: if the masses of the two bodies are very different, the ``small'' body is treated as a point particle of mass ${m}_{p}$ moving in the gravitational field generated by the large mass $M$, and one keeps only linear terms in the small mass ratio ${m}_{p}/M$. This technique usually yields finite answers, which are often in good agreement with fully nonlinear numerical relativity results, even when extrapolated to nearly comparable mass ratios. Here we study two situations in which the point-particle approximation yields a divergent result: the instantaneous flux emitted by a small body as it orbits the light ring of a black hole, and the total energy absorbed by the horizon when a small body plunges into a black hole. By integrating the Teukolsky (or Zerilli/Regge-Wheeler) equations in the frequency and time domains we show that both of these quantities diverge. We find that these divergences are an artifact of the point-particle idealization, and are able to interpret and regularize this behavior by introducing a finite size for the point particle. These divergences do not play a role in black-hole imaging, e.g., by the Event Horizon Telescope.

A GPU-Accelerated Mixed-Precision WENO Method for Extremal Black Hole and Gravitational Wave Physics Computations

We develop and use a novel mixed-precision weighted essentially non-oscillatory (WENO) method for solving the Teukolsky equation, which arises when modeling perturbations of Kerr black holes. We show that WENO methods outperform higher-order finite-difference methods, standard in the discretization of the Teukolsky equation, due to the need to add dissipation for stability purposes in the latter. In particular, as the WENO scheme uses no additional dissipation, it is well suited for scenarios requiring long-time evolution such as the study of price tails and gravitational wave emission from extreme mass ratio binaries. In the mixed-precision approach, the expensive computation of the WENO weights is performed in reduced floating-point precision that results in a significant speedup factor of \(\approx 3.3\). In addition, we use state-of-the-art Nvidia general-purpose graphics processing units and cluster parallelism to further accelerate the WENO computations. Our optimized WENO solver can be used to quickly generate accurate results of significance in the field of black hole and gravitational wave physics. We apply our solver to study the behavior of the Aretakis charge—a conserved quantity, that if detected by a gravitational wave observatory like LIGO/Virgo would prove the existence of extremal black holes.

AT 2020iko: A WZ Sge-type Dwarf Nova Candidate with an Anomalous Precursor Event

Abstract The ongoing Zwicky Transient Facility (ZTF) survey is generating a massive alert rate from a variety of optical transients and variable stars, which are being filtered down to subsets meeting user-specified criteria by broker systems such as the Arizona-NOIRLab Temporal Analysis and Response to Events System (ANTARES). In a beta implementation of the algorithm of Soraisam et al. on ANTARES, we flagged AT 2020iko from the ZTF real-time alert stream as an anomalous source. This source is located close to a red extended Sloan Digital Sky Survey source. In the first few epochs of detection, it exhibited a V-shaped brightness profile, preceded by nondetections both in ZTF and in the All-Sky Automated Survey for Supernovae extending to 2014. Its full light curve shows a precursor event, followed by a main superoutburst and at least two rebrightenings. A low-resolution spectrum of this source points to a dwarf nova (DN) nature. Although some of the features of AT 2020iko indicate an SU UMa-type DN, its large amplitude, presence of rebrightenings, and inferred supercycle period of ≥6 yr are in favor of AT 2020iko being a new WZ Sge-type DN candidate, a subset of rare DNe consisting of extreme mass-ratio (<0.1) binaries with an orbital period around the period minimum. The precusor event of AT 2020iko brightened by 6.5 mag, while its decay spanned 3–5 mag. We speculate this superoutburst is associated with a less expanded accretion disk than in typical superoutbursts in WZ Sge systems, with the large depth of the precursor decay implying an extremely small mass ratio. To the best of our knowledge, such a precursor event has not been recorded for any DN. This result serves to demonstrate the efficacy of our real-time anomaly search algorithm.

Dissipation in extreme mass-ratio binaries with a spinning secondary

We present the gravitational-wave flux balance law in an extreme mass-ratio binary with a spinning secondary. This law relates the flux of energy (angular momentum) radiated to null infinity and through the event horizon to the local change in the secondary's orbital energy (angular momentum) for generic (non-resonant) bound orbits in Kerr spacetime. As an explicit example we compute these quantities for a spin-aligned body moving on a circular orbit around a Schwarzschild black hole. We perform this calculation both analytically, via a high-order post-Newtonian expansion, and numerically in two different gauges. Using these results we demonstrate explicitly that our new balance law holds.

Gravitational burst radiation from pulsars in the Galactic centre and stellar clusters

ABSTRACT Pulsars (PSRs) orbiting intermediate or supermassive black holes at the centre of galaxies and globular clusters are known as Extreme Mass Ratio Binaries (EMRBs) and have been identified as precision probes of strong-field GR. For appropriate orbital parameters, some of these systems may also emit gravitational radiation in a ‘burst-like’ pattern. The observation of this burst radiation in conjunction with the electromagnetic radio timing signal would allow for multimessenger astronomy in strong-field gravitational regimes. In this work we investigate gravitational radiation from these PSR-EMRBs, calculating the waveforms and SNRs and explore the influence of this GW on the pulsar radio signal. We find that for typical PSR-EMRBs, gravitational burst radiation should be detectable from both the Galactic centre and the centre of stellar clusters, and that this radiation will not meaningfully affect the pulsar timing signal, allowing PSR-EMRB to remain ‘clean’ test-beds of strong-field GR.

Polar planets around highly eccentric binaries are the most stable

ABSTRACT We study the orbital stability of a non-zero mass, close-in circular orbit planet around an eccentric orbit binary for various initial values of the binary eccentricity, binary mass fraction, planet mass, planet semimajor axis, and planet inclination by means of numerical simulations that cover 5 × 104 binary orbits. For small binary eccentricity, the stable orbits that extend closest to the binary (most stable orbits) are nearly retrograde and circulating. For high binary eccentricity, the most stable orbits are highly inclined and librate near the so-called generalized polar orbit which is a stationary orbit that is fixed in the frame of the binary orbit. For more extreme mass ratio binaries, there is a greater variation in the size of the stability region (defined by initial orbital radius and inclination) with planet mass and initial inclination, especially for low binary eccentricity. For low binary eccentricity, inclined planet orbits may be unstable even at large orbital radii (separation ${\gt}5 \, a_{\rm b}$). The escape time for an unstable planet is generally shorter around an equal mass binary compared with an unequal mass binary. Our results have implications for circumbinary planet formation and evolution and will be helpful for understanding future circumbinary planet observations.

Repeated ringing of black holes: Quasinormal bursts from highly eccentric, extreme mass-ratio binaries

Recent studies of scalar and gravitational waveforms from high-eccentricity, extreme mass-ratio black hole binaries show the presence of quasi-normal bursts (QNB) -- lingering high frequency decaying oscillations (also known as ``wiggles'') -- soon after each periapsis passage. One puzzle associated with these QNB is that in the case of a nearly-extreme rotating central black hole the frequency of the QNB has been found to be in a range which is lower than the corresponding range of relevant quasi-normal modes. We reproduce these results using a different approach and perform a detailed analysis to find evidence for the resolution of the puzzle and for the origin of the QNB. We find that the QNB frequency as measured at future null infinity evolves in (retarded) time and approaches the dominant quasi-normal frequency exponentially in time. We also show that the QNB amplitude decays inversely in (retarded) time. We discuss the time dependence of both the QNB waveform frequency and its amplitude and argue that this behavior arises as a result of the excitation of many quasi-normal overtones and the summation thereof.

Pulsar timing in extreme mass ratio binaries: a general relativistic approach

The detection of a pulsar (PSR) in a tight, relativistic orbit around a supermassive or intermediate mass black hole - such as those in the Galactic centre or in the centre of Globular clusters - would allow for precision tests of general relativity (GR) in the strong-field, non-linear regime. We present a framework for calculating the theoretical time-frequency signal from a PSR in such an Extreme Mass Ratio Binary (EMRB). This framework is entirely relativistic with no weak-field approximations and so able to account for all higher-order strong-field gravitational effects, relativistic spin dynamics, the convolution with astrophysical effects and the combined impact on the PSR timing signal. Specifically we calculate both the spacetime path of the pulsar radio signal and the complex orbital and spin dynamics of a spinning pulsar around a Kerr black hole, accounting for spacetime curvature and frame dragging, relativistic and gravitational time delay, gravitational light bending, temporal and spatial dispersion induced by the presence of plasma along the line of sight and relativistic aberration. This then allows for a consistent time-frequency solution to be generated. Such a framework is key for assessing the use of PSR as probes of strong field GR, helping to inform the detection of an EMRB system hosting a PSR and, most essentially, for providing an accurate theoretical basis to then compare with observations to test fundamental physics.

Spin dynamics of a millisecond pulsar orbiting closely around a massive black hole

We investigate the spin dynamics of a millisecond pulsar (MSP) in a tightly bounded orbit around a massive black hole. These binaries are progenitors of the extreme-mass-ratio-inspirals (EMRIs) and intermediate-mass-ratio-inspirals (IMRIs) gravitational wave events. The Mathisson-Papapetrou-Dixon (MPD) formulation is used to determine the orbital motion and spin modulation and evolution. We show that the MSP will not be confined in a planar Keplerian orbit and its spin will exhibit precession and nutation induced by spin-orbit coupling and spin-curvature interaction. These spin and orbital behaviours will manifest observationally in the temporal variations in the MSP's pulsed emission and, with certain geometries, in the self-occultation of the pulsar's emitting poles. Radio pulsar timing observations will be able to detect such signatures. These extreme-mass-ratio binaries (EMRBs) and intermediate-mass-ratio binaries (IMRBs) are also strong gravitational wave sources. Combining radio pulsar timing and gravitational wave observations will allow us to determine the dynamics of these systems in high precision and hence the subtle behaviours of spinning masses in strong gravity.

Tidal disruption by extreme mass ratio binaries and application to ASASSN-15lh

Tidal disruption events (TDEs) observed in massive galaxies with inferred central black hole masses $M_h > 10^8 \ M_\odot$ are presumptive candidates for TDEs by lower mass secondaries in binary systems. We use hydrodynamic simulations to quantify the characteristics of such TDEs, focusing on extreme mass ratio binaries and mpc separations where the debris stream samples the binary potential. The simulations are initialised with disruption trajectories from 3-body integrations of stars with parabolic orbits with respect to the binary center of mass. The most common outcome is found to be the formation of an unbound debris stream, with either weak late-time accretion or no accretion at all. A substantial fraction of streams remain bound, however, and these commonly yield structured fallback rate curves that exhibit multiple peaks or sharp drops. We apply our results to the superluminous supernova candidate ASASSN-15lh and show that its features, including its anomalous rebrightening at $\sim 100$ days after detection, are consistent with the tidal disruption of a star by a supermassive black hole in a binary system.

THE FIRST PHOTOMETRIC INVESTIGATION OF THE NEGLECTED W-UMa-TYPE BINARY STAR V532 MONOCEROTIS

ABSTRACT The first two sets of complete charge-coupled device light curves in the B and V bands of the short-period binary system, V532 Mon, are presented. The light curves are analyzed with spot models using the Wilson–Devinney code. V532 Mon is found to be an A-subtype intermediate-contact binary with a degree of contact factor of , a mass ratio of , and a cool spot on the primary component. The period changes of the system are investigated by combining newly determined times of light minimum with others published in the literature. We found that the general trend of the observed-calculated curve shows a downward parabolic variation that corresponds to a long-term decrease in the orbital period at a rate of −1. The long-term decrease of the orbital period can be explained by mass transfer from the more-massive component to the less massive one. The long-time period decrease, the intermediate-contact configuration, and the astrophysical parameters of the binary system suggest that V532 Mon will evolve into high fill-out, extreme mass-ratio overcontact binaries.

Analytical high-order post-Newtonian expansions for spinning extreme mass ratio binaries

We present an analytic computation of Detweiler's redshift invariant for a point mass in a circular orbit around a Kerr black hole, giving results up to 8.5 post-Newtonian order while making no assumptions on the magnitude of the spin of the black hole. Our calculation is based on the functional series method of Mano, Suzuki and Takasugi, and employs a rigorous mode-sum regularization prescription based on the Detweiler-Whiting singular-regular decomposition. The approximations used in our approach are minimal; we use the standard self-force expansion to linear order in the mass ratio, and the standard post-Newtonian expansion in the separation of the binary. A key advantage of this approach is that it produces expressions that include contributions at all orders in the spin of the Kerr black hole. While this work applies the method to the specific case of Detweiler's redshift invariant, it can be readily extended to other gauge invariant quantities and to higher post-Newtonian orders.

Analytical high-order post-Newtonian expansions for extreme mass ratio binaries

We present analytic computations of gauge invariant quantities for a point mass in a circular orbit around a Schwarzschild black hole, giving results up to 15.5 post-Newtonian order in this paper and up to 21.5 post-Newtonian order in an online repository. Our calculation is based on the functional series method of Mano, Suzuki and Takasugi (MST) and a recent series of results by Bini and Damour. We develop an optimised method for generating post-Newtonian expansions of the MST series, enabling significantly faster computations. We also clarify the structure of the expansions for large values of $\ell$, and in doing so develop an efficient new method for generating the MST renormalised angular momentum, $\nu$.

A NEW CLASS OF NASCENT ECLIPSING BINARIES WITH EXTREME MASS RATIOS

Early B-type main-sequence (MS) stars (M$_1$ = 5-16 M$_{\odot}$) with closely orbiting low-mass stellar companions (q = M$_2$/M$_1$ < 0.25) can evolve to produce Type Ia supernovae, low-mass X-ray binaries, and millisecond pulsars. However, the formation mechanism and intrinsic frequency of such close extreme mass-ratio binaries have been debated, especially considering none have hitherto been detected. Utilizing observations of the Large Magellanic Cloud galaxy conducted by the Optical Gravitational Lensing Experiment, we have discovered a new class of eclipsing binaries in which a luminous B-type MS star irradiates a closely orbiting low-mass pre-MS companion that has not yet fully formed. The primordial pre-MS companions have large radii and discernibly reflect much of the light they intercept from the B-type MS primaries ($\Delta$I$_{\rm refl}$ = 0.02-0.14 mag). For the 18 definitive MS + pre-MS eclipsing binaries in our sample with good model fits to the observed light curves, we measure short orbital periods P = 3.0-8.5 days, young ages $\tau$ = 0.6-8 Myr, and small secondary masses M$_2$ = 0.8-2.4 M$_{\odot}$ (q = 0.07-0.36). The majority of these nascent eclipsing binaries are still associated with stellar nurseries, e.g. the system with the deepest eclipse $\Delta$I$_1$ = 2.8 mag and youngest age $\tau$ = 0.6$\pm$0.4 Myr is embedded in the bright H II region 30 Doradus. After correcting for selection effects, we find that (2.0$\pm$0.6)% of B-type MS stars have companions with short orbital periods P = 3.0-8.5 days and extreme mass ratios q = 0.06-0.25. This is $\approx$10 times greater than that observed for solar-type MS primaries. We discuss how these new eclipsing binaries provide invaluable insights, diagnostics, and challenges for the formation and evolution of stars, binaries, and H II regions.

Light-curve analysis and orbital period change of the extreme mass-ratio overcontact binary AW CrB

Extreme mass-ratio contact binaries with a high degree of overcontact may be in the late evolutionary stages of the contact phase. Detailed photometric analyses and orbital period studies of those systems can provide invaluable information for the coalescence scenario of close binary systems, as recently observed in V1309 Sco. In this paper the first light curve analysis and period study for the totally eclipsing contact binary AW CrB is presented. The VRcIc CCD photometric light curves are analysed by means of the WD code. The asymmetry of the light curves is modelled by a cool star spot on the primary component. It is shown that AW CrB has a high degree of overcontact f=75% and an extreme mass-ratio of q=0.10, placing it among the few contact binaries with the lowest known mass-ratios. The mean density of the primary component suggest that it has evolved away from the ZAMS. Eighty-four times of minimum light are reported, derived from observations available in public archives, the literature, and the observations presented in this paper. The orbital period study shows a continuously increasing period, at a rate of dP/dt=3.58x10^-7 d yr^-1, suggesting a mass transfer from the secondary component to the primary one, implying a further decreasing mass-ratio.