Searching for Exoplanets Around X-Ray Binaries with Accreting White Dwarfs, Neutron Stars, and Black Holes

TL;DR

This paper proposes using X-ray light curve analysis of X-ray binaries with white dwarfs, neutron stars, or black holes to detect exoplanets and potential extraterrestrial structures, expanding current search methods.

Contribution

It introduces a novel method for detecting planets and artificial structures around X-ray binaries through X-ray transit signatures, complementing existing optical and radio techniques.

Findings

X-ray transits can reveal planets in close orbits around X-ray binaries.

Transit durations range from minutes to hours depending on planet size and orbit.

The method can also potentially detect extraterrestrial Dyson spheres.

Abstract

We recommend that the search for exoplanets around binary stars be extended to include X-ray binaries in which the accretor is a white dwarf, neutron star, or black hole. We present a novel idea for detecting planets bound to such mass transfer binaries: we propose that the X-ray light curves of these binaries be inspected for signatures of transiting planets. X-ray transits may be the only way to detect planets around some systems, while providing a complementary approach to optical and/or radio observations in others. Any planets associated with X-ray binaries must be in stable orbits. We consider the range of allowable separations and find that orbital periods can be hours or longer, while transit durations extend upward from about a minute for Earth-radius planets in very close orbits, to hours for Jupiter-radius planets in wider orbits. The search for planets around mass transfer binaries could begin at once with existing X-ray observations of these systems. If and when a planet is detected around an X-ray binary, the size and mass of the planet may be readily measured, and it may also be possible to study the transmission and absorption of X-rays through its atmosphere. Finally, a noteworthy application of our proposal is that the same technique could be used to search for signals from extraterrestrial intelligence. If an advanced exocivilization placed a Dyson sphere or similar structure in orbit around the accretor of an X-ray binary in order to capture energy, such an artificial structure might cause detectable transits in the X-ray light curve.

Figures

Tables

Table 1: Variables used in this paper.
Symbol	Variable definition
$M_1\left(M_2\right)$	Mass of primary (donor) star
$R_1\left(R_2\right)$	Radius of primary (donor) star
$M_p,R_p,P_p$	Planet mass, radius, and period
$\mu$	Mass ratio: $\mu =M_2/\left(M_1+M_2\right)$
$q$	$q=M_2/M_1$
$k$	Planet-to-primary size ratio: $k=R_p/R_1$
$a_b$	Semi-major axis of binary
$a_c$	Critical semi-major axis of planetary orbit
$P$	Orbital period of planetary orbit
$e,b$	Eccentricity, impact parameter
$i_p\left(i_b\right)$	Inclination angle of planetary (binary) orbit
$\delta i$	Mutual inclination of planetary and binary orbits
$R_T$	Tidal disruption radius
$T\left(T_{\mathrm{tot}}\right)$	Full (total) duration time of transit
$\tau$	Ingress/egress time of transit
$\delta$	Transit depth
$P_{\mathrm{tran}}$	Transit probability

Equations

$$\begin{array}{rl}\frac{a_{c,\mathrm{out}}}{a_b}=& (1.60\pm 0.04)+(5.10\pm 0.05)e\\ & +(-2.22\pm 0.11)e^2+(4.12\pm 0.09)\mu \\ & +(-4.27\pm 0.17)e\mu +(-5.09\pm 0.11){\mu }^2\\ & +(4.61\pm 0.36)e^2{\mu }^2.\end{array}$$

$$\begin{array}{rl}\frac{a_{c,\mathrm{in}}}{a_b}=& (0.464\pm 0.006)+(-0.380\pm 0.010)\mu \\ & +(-0.631\pm 0.034)e+(0.586\pm 0.061)e\mu \\ & +(0.150\pm 0.041)e^2+(-0.198\pm 0.074)\mu e^2.\end{array}$$

$$\frac{a_b}{R_2}=\frac{0.6q^{2/3}+\ln (1+q^{1/3})}{0.49q^{2/3}},$$

$$T_{\mathrm{tot}}=\frac{P}{\pi }{\sin }^{-1}\left[\frac{R_1}{a_c}\frac{\sqrt{(1+k{)}^2-b^2}}{\sin i_p}\right],$$

$$T_{\mathrm{full}}=\frac{P}{\pi }{\sin }^{-1}\left[\frac{R_1}{a_c}\frac{\sqrt{(1-k{)}^2-b^2}}{\sin i_p}\right],$$

$$\tau =\frac{1}{2}(T_{\mathrm{tot}}-T_{\mathrm{full}}).$$

$$P^2=\frac{4{\pi }^2{a}_c^3}{G(M_{\star }+M_p)},$$

$$P_{\mathrm{cr}}=\frac{\min [\delta \mathrm{\Omega },2\pi ]}{2\pi },$$

$$\begin{array}{rl}{T}_{\mathrm{prec}}& =\frac{2\pi }{|\dot{\Omega }|}\\ & =P_p\frac{4}{3\cos \delta i}{\left(\frac{a_p}{a_b\mu }\right)}^2\frac{m_2}{m_1},\end{array}$$

$$\begin{array}{rl}d{l}_1& =t_{\mathrm{tran}}|v_p\cos \delta i-\frac{2v_1}{\pi }|+R_1+R_p\\ d{l}_2& =t_{\mathrm{tran}}|v_p\cos \delta i+\frac{2v_1}{\pi }|+2(R_1+R_p),\end{array}$$

$$P_{\mathrm{tran}}=P_{\mathrm{cr}}[1-(1-P_1{)}^n],$$

$$P_1=\{\begin{array}{ll}1& \text{if}\frac{d{l}_1+d{l}_2}{2}>2a_b\mu ,\\ \frac{d{l}_1+d{l}_2}{4a_b\mu }& \text{otherwise},\end{array}$$

$$C=\frac{dc}{dt}\,\Bigg{\{}{\rm min}\Bigg{[}\Bigg{(}\frac{{\cal A}_{\rm cov}}{\pi R_{\rm em}^{2}}\Bigg{)}^{2},1\Bigg{]}\Bigg{\}}\,\times\Delta\,t,$$

$$\Delta\,t=17~{}{\rm min}\,\Bigg{(}\frac{R_{\rm eff}}{R_{J}}\Bigg{)}\,\Bigg{(}\frac{a_{b}}{10\,R_{\odot}}\,\frac{M_{\odot}}{M_{\star}}\Bigg{)}^{\frac{1}{2}},$$

Full text

1. Introduction

X-ray binaries (XRBs) are gravitationally interacting binaries in which a stellar remnant accretes matter from its companion. The term XRB is often used for binaries in which the accretor is a neutron star or black hole, and here we will apply it as well to binaries in which the accretor is a white dwarf. We include both cataclysmic variables (CVs), in which the white dwarf accretes from a close low-mass companion, and symbiotics, in which the donor is a giant. We consider the possibility that some XRBs host planetary systems, and we present a method to discover planets orbiting within or around such interacting binaries.

1.1. Planets in Mass-Transfer Binaries: Feasibility

Planets in binaries may be common. At the time of this writing, roughly 70 examples are known of planets orbiting one member of a binary (e.g., Butler et al., 1997; Cochran et al., 1997; Roell et al., 2012), and more than a dozen cases are known in which the planet orbits both members (e.g., Thorsett et al., 1999; Doyle et al., 2011; Hinse et al., 2015; Orosz et al., 2012a, b). Following the terminology of Dvorak (1986), these two types of systems are S-type and P-type, respectively. The first known P-type system consists of a planet orbiting a pulsar and a white dwarf and was detected from observed changes in the travel time of the radio pulses (Thorsett et al., 1999). This system passed through an epoch as an X-ray binary, during which the giant progenitor of today’s white dwarf transferred mass to the neutron star, spinning it up Here we consider how to search for planets in active mass-transfer systems.

Planets have also been discovered in P-type orbits around binaries which will become CVs during their next evolutionary phase (e.g., HW Vir, Lee et al. 2009; HW Vir, Beuermann et al. 2012; DP Leo, Qian et al. 2010; NN Ser, Beuermann et al. 2010). These binaries are referred to as pre-CVs or post-common-envelope binaries. Pre-CVs contain a white dwarf which is the remnant of a giant that filled its Roche lobe under circumstances inconsistent with dynamically stable mass transfer. The giant’s envelope was transformed into a common envelope within which the white dwarf and its stellar companion spiraled closer to each other, but did not merge. Eventually, the white dwarf’s companion will fill its Roche lobe, and the binary will become a CV. If the planet is retained, it will be in orbit with an XRB.

Since planets have been found in P-type orbits both prior to and after the binary is an XRB, it is clearly important to consider planets during the XRB phase. Planets in S-type orbits around accreting compact objects or companion donor stars may also be present, especially in wider binaries, where mass transfer may proceed through winds. Those planets not destroyed by mass transfer may transit the compact object and portions of its accretion disk, potentially causing detectable X-ray eclipses. The same considerations apply, whether the accretor is a white dwarf, neutron star, or black hole.

1.2. The Possibility of Detection Through X-ray Transits

X-ray light curves provide ideal opportunities to search for transits, because the sizes of the X-ray emitting regions in many XRBs are so small that the passage of a planet produces a sharply-defined total eclipse. For example, X-ray emission from black holes can be dominated by a soft ($\sim$ keV) component emanating primarily from the inner portion of an accretion disk, which may have an extent of only tens or hundreds of kilometers (e.g., Remillard & McClintock, 2006). X-rays from neutron stars may also emerge from very small regions. In certain states, emission from a white dwarf emanates from an area comparable to its surface area, so that an Earth-size planet could cause a deep eclipse, and somewhat larger planets could produce total eclipses. This would apply to a white dwarf in a supersoft state (e.g., van den Heuvel et al., 1992), or post-nova, and also to some magnetic CVs (e.g., Nucita et al., 2009).

The question of detectability hinges on whether the number of photons missing due to the eclipse is large enough that the deficit and its duration can be reliably measured. If reliable detection cannot be made during a single transit, then the signal could possibly be identified through repeated transits, if the orbital period is short enough. The challenges are similar to those that have been successfully met by analyses of data from the Kepler mission, which has identified Earth-size and even smaller planets orbiting Sun-like stars. The X-ray data analyses we advocate would take place in a complementary part of the discovery parameter space: the X-ray transits are deeper, in some cases total eclipses; the background is small or nonexistent, eliminating problems of blending with light from nearby stars; and the numbers of X-ray photons are small. The paucity of X-ray photons is a new challenge which provides a lower bound on the baseline count rate relative to the eclipse time scale and the total exposure time (see §4). Nevertheless, if planets orbit within or around XRBs, then detectability is possible for many systems.

1.3. Complementary Methods

X-ray binaries are also copious emitters of optical light. In cases in which the optical emission, like the X-ray emission, emanates from a small region, optical transits may also be detected. HT Cas and OY Car are examples of CVs in which both optical and X-ray eclipses have been studied (Wood et al., 1989; Horne et al., 1991; Wheatley & West, 2003; Thorstensen et al., 2008; Nucita et al., 2009). In many XRBs, however, the sources of optical emission can be complex and extended. The donor star and outer regions of the accretion disk can both contribute to or dominate the optical emission, diluting the effects of stellar or planetary transits. In these cases, X-rays may provide the best route to planet detection. Moreover, the detection of transits in two widely separated wavebands increases confidence in the interpretation that the dips are related to the passage of a solid body. In addition, different wavebands reveal complementary aspects of the emission during mass transfer. Periodic pulsations, like those emanating from millisecond pulsars, can reveal the presence of planets through radio timing measurements. As with optical observations, combining information across wavebands is an exciting prospect.

1.4. Plan of this Paper

We demonstrate that planets may produce detectable X-ray eclipses. We consider two extremes of planet masses and sizes, from Earth- to Jupiter-sized. Our main goals are (1) to examine the parameter space of stable orbits for planets in XRBs; (2) to study the expected light curve characteristics; and (3) to consider which XRBs would likely be most profitable to study with our approach.

In Section 2, we present the equations for investigating the stability limits of planets orbiting XRBs, with a specific focus on CVs. We provide examples of theoretical transiting light curves in Section 3. We address the issue of detectability and consider the properties of XRBs that are the most promising targets for light curve analysis in Section 4. In Section 5 we summarize our conclusions and discuss the implications of our work.

2. Orbital Stability

Our goal in this section is to explore a realistic range within the parameter space of stable orbits for planets in XRBs. For simplicity, we restrict our considerations to single planets in circular, coplanar orbits. An illustration of a planet orbiting and transiting the compact accretor of an XRB is shown in Figure 1.

Stability in triple systems generally requires a hierarchical structure, with the outer orbit being at least a few times larger than the inner orbit. Thus, for P-type (or, circumbinary) orbits, there is a critical value, $a_{c,{\rm out}}$, for the semimajor axis of the planetary orbit in terms of the separation $a_{b}$ between the components of the binary. Holman & Wiegert (1999) determined an expression for $a_{c,{\rm out}}$ by conducting numerical simulations of test particles moving in the gravitational field of binaries:

[TABLE]

Here, $e$ is the orbital eccentricity, which will generally be zero in the mass transfer binaries we consider. The binary mass ratio is $\mu=M_{2}/(M_{1}+M_{2})$, where $M_{1}$ is the mass of primary, and $M_{2}$ is the mass of its stellar companion. In our case, without loss of generality, we define the primary to be the compact accretor. This equation applies over the range $0.1\leq\mu\leq 0.9$ (Holman & Wiegert, 1999). For $a>a_{c,\rm}$ the orbit is stable.

It is also possible for the planet to be located within the binary orbit. In this case, Holman & Wiegert (1999) find that the critical orbital radius of an S-type planet, $a_{c,\rm in},$ is

[TABLE]

Equation (1) is valid to within 3–6% over the range $0.0\leq e\leq 0.7$, and equation (2) is valid to within 4–11% for $0.0\leq e\leq 0.8$. The variables used in the above equations and the rest of this paper are summarized in Table 1.

The specific values we derive for the innermost and outermost stable orbits are not expected to represent exact results for all planets in binaries, because some of the input assumptions may not be satisfied. For example, the orbits may not be coplanar, there may be a distant third star, multiple planets, or the process of mass transfer may alter the results somewhat. Nevertheless, the general result should be robust: the planetary and stellar orbits must be hierarchical, with the semimajor axis of a stable planetary orbit differing from the semimajor axis of the binary orbit by a factor of a few.

When mass transfer occurs because the donor fills its Roche lobe, the binary’s orbital separation is set by the condition that the radius of the donor, $R_{2}$, equals the size of its Roche lobe, $R_{L}$. Using the function derived by Eggleton (1983), we have:

[TABLE]

where $q=M_{2}/M_{1}$.

In Figure 2, we display the binary separation of XRBs as a function of $\mu$, for a fixed accretor mass of $M_{1}=0.8~{}M_{\odot}$, corresponding to a white dwarf. We consider a case in which a main-sequence donor, with $R_{2}=M_{2}(R_{\odot}/M_{\odot})$, fills its Roche lobe. The orbits are circular ($e=0$). The plot also shows the critical semimajor axes for a single planet. We show both outer ($a_{c,\rm out}$; dotted line) and inner ($a_{c,\rm in}$; dashed line) critical orbits. Overplotted are data points from Zorotovic et al. (2011) of binary separations for real CVs. There are no CVs on the upper right portion of the curve because mass transfer becomes unstable when the donor is much more massive than the accretor. The CVs depart from our curve on the lower left because the low-mass donors in these systems are degenerate and have larger radii than would be given by the extrapolation of the radius-mass relation for main-sequence stars.

We note that the orbital separation can be significantly larger in mass transfer binaries. The value of $a_{b}$ is generally $2$-$3$ times the radius of the donor, when the donor fills its Roche lobe. Thus, when the donor is a subgiant or giant, $a_{b}$ is larger than in the main-sequence case we have shown. In addition, mass transfer is often effected through winds. The winds can be driven by radiation from the accretor or else may be generated by a highly evolved donor. In these cases, the separation is larger than it would be for Roche-lobe-filling systems.

Another factor that must be taken into account is the value of the mass ratio, since large values are associated with mass transfer that is unstable on dynamical time scales. For mass transfer mediated by winds, there is not a fixed analytic limit on the mass ratio, since larger mass ratios are allowed if much of the mass lost from the donor exits the system carrying only a small amount of angular momentum per unit mass. For mass transfer that occurs through the L1 Lagrange point, the limit depends on the evolutionary state of the donor, the fraction of the mass it loses that is ejected from the system, and the amount of angular momentum carried by ejected mass. The critical mass ratio for stability may have values somewhat smaller than unity, although values as large as $4$ can apply to main-sequence donors (e.g., Ivanova, 2015).

A planet that ventures too close to a white dwarf or other compact object would not be able to survive disruption by tidal forces. The tidal disruption radius (the Roche limit) is $R_{T}=1.16~{}R_{\odot}(\rho_{p}/\rho_{\earth})^{-1/3}(M_{1}/0.6~{}M_{\odot})^{1/3}$, where $\rho_{p,\earth}$ are the mass densities of the planet and of Earth, respectively. Tidal forces present a challenge for the survival of inner orbit (S-orbit) planets.

In Figure 3, we show the orbital parameters for several Roche-lobe filling binaries. Each panel shows the critical inner semimajor axis, $a_{c,\rm in}$, as a function of donor mass. We show plots for four different values of the primary star’s mass, $M_{1}=0.8$, 1.4, 10, and 50 $M_{\odot}$. The horizontal dashed line marks the tidal radius for an Earth-density planet, ($\rho_{\earth}\approx 5.5~{}\rm{g~{}cm}^{-3};~{}R_{T,\earth}$, and a Jupiter-density planet, ($\rho_{\rm J}\approx 0.24\rho_{\earth};~{}R_{T,\rm J}$). Thus, for instance, a Jupiter-density planet in an XRB with a primary of mass $0.8~{}M_{\odot}$ (top panel) would be ripped apart for binary orbital separations less than $R_{T,\rm J}\approx 2.1R_{\odot}$. An Earth-density planet in such a system could not survive orbital separations less than $R_{T,\earth}\approx 1.3R_{\odot}$.

These plots illustrate several key points. First, as the formula for $R_{T}$ indicates, survival is more likely for planets of higher density. Second, for accretors with masses in the range expected for white dwarfs and neutron stars, we do not expect to find stable S-type orbits if the donor star fills its Roche lobe. Thus, neither CVs nor low-mass X-ray binaries with neutron-star accretors are likely to have planets in stable S-type orbits. It is only for accretors of higher mass that we can find stable S-type orbits for binaries in which the donor is less massive than the accretor, ensuring the dynamical stability of mass transfer. In fact, the discovery of a planet in an S-type orbit may place lower bounds on the accretor’s mass, something very useful in the study of black holes.

The above results do not rule out S-type orbits in mass transfer binaries. First, when the donor is larger, either a subgiant or giant, the orbital separation is larger even when the donor fills its Roche lobe. Second, when mass transfer occurs through winds, the orbit can be even larger. When the stellar orbits are wide, planets can orbit around either or both stars.

3. Light Curves

To model the transit light curves, we need to assess the duration and the ingress/egress times of transits for a set of orbital parameters. A cartoon illustration of a transit in Figure 1 shows that, for a non-grazing eclipse, the stellar and planetary disks touch at four contact times $t_{\rm I}$–$t_{\rm IV}$. (For a grazing eclipse, the second and third contacts do not take place.) The transit duration $T$ and ingress/egress duration $\tau$ depend on the planet period, $P$, the radius of the primary star, $R_{1}$, the radius of the planet, $R_{p}$, the impact parameter, $b$, and the orbital inclination of the planet projected onto the plane of the sky, $i_{p}$. For a circular orbit, the total duration, $T_{\rm tot}\equiv t_{\rm IV}-t_{\rm I}$, and full duration, $T_{\rm full}\equiv t_{\rm III}-t_{\rm II}$, are given by (e.g., Winn 2010),

[TABLE]

[TABLE]

where $k=R_{p}/R_{1}$. The ingress and egress durations, $\tau_{\rm ing}=t_{\rm II}-t_{\rm I}$ and $\tau_{\rm egr}=t_{\rm IV}-t_{\rm III}$ are usually unequal for eccentric orbits, though the difference is minor in practice. Here, since we consider only circular orbits, we take $\tau\equiv\tau_{\rm ing}=\tau_{\rm egr}$, and thus

[TABLE]

In many XRBs, the source of X-rays is a small region. It could for example be the inner accretion disk of a black hole or of a neutron star. In such cases the duration of the eclipse is determined by the amount of time it takes for the planet to cross this very small region. In the case where the primary is a white dwarf of mass $M_{1}$, the size of the X-ray-emitting region can be close to the size of the white dwarf (e.g., Nucita et al., 2009), $R_{1}$, which is determined by the white dwarf mass-radius relation: $R_{1}\approx 0.01(M_{1}/M_{\odot})^{-1/3}R_{\odot}$.

The orbital period of the planet is simply derived using Kepler’s third law,

[TABLE]

where $M_{p}$ is the planet mass, and $a_{c}$ is determined from equations (1) or (2) for outer or inner orbits, respectively. For inner orbits around the primary star, $M_{\star}$ would refer to the mass of the primary star alone. For outer orbits, however, the planet circles the combined mass of the primary and donor star, so $M_{\star}=M_{1}+M_{2}$. We consider Earth-mass and Jupiter-mass planets. Figure 4 shows the planetary period for outer and inner orbits as a function of $\mu$ for a CV in which the primary is a 0.8 $M_{\odot}$ white dwarf. Because $\mu$ varies and $M_{1}$ is fixed, this is equivalent to plotting the period for a range of donor masses, $M_{2}=\mu M_{1}(1-\mu)$. Since $M_{\star}\gg M_{p}$, the periods of Earth- and Jupiter-mass planets are essentially equal. Recall that we have derived the period for values of the critical semimajor axis, within which planetary orbits are unstable. But orbits greater than $a_{c}$ are also permissible, and therefore the periods we calculate should be considered lower limits for a given $M_{\star}$. In Figure 4, we overplot the periods of potential planets in circumbinary (CB) orbits around the Zorotovic et al. (2011) sample of CVs. Using the measured masses of the binary members to define $\mu$, we calculate the periods assuming the planets are in stable orbits at the critical semimajor axis, $a_{c,\rm out}$. These calculations suggest that the periods of CB planets around XRBs may be as short as a few to tens of hours.

For completeness, we plot the period for inner orbits in Figure 4 as well, though based on the analysis of the tidal disruption radius in §2, such orbits may be physically implausible in CVs and other low-mass XRBs.

Figure 5 plots the full duration and ingress/egress transit times versus $\mu$ for both Earth- and Jupiter-mass planets orbiting outside the binary. Again, we show the scenario where the primary mass is fixed to that of a 0.8 $M_{\odot}$ white dwarf. And we focus our attention on simplified cases where the planet transits directly across the equator of the primary, i.e., $i_{p}=90\arcdeg$ and $b=0$. The full duration and ingress/egress times displayed in Figure 5 are lower limits for the scenarios considered; that is, these are the minimum transit times expected for CB systems in which the planets are near the stability limit. If they exist, planets in XRB systems may certainly populate orbits beyond the stability limit. Yet in the currently known CB systems discovered with Kepler—if they can provide us with any direction—the planets are mostly close to the stability limit (Welsh et al., 2015; Li et al., 2016). The general trend is for both transit and ingress/egress durations to increase with binary mass ratio, $\mu$. This is because as either the primary or donor mass increases, the orbit at which a planet is stable is pushed further outward, corresponding to an increase in the orbital period.

Small planets are expected to have short transit durations, on the order of a few seconds, and comparatively longer ingress/egress times, because their sizes are comparable to the sizes of white dwarfs ($\sim R_{\earth}$). Due to their similar sizes, a small planet will spend relatively little time fully covering the stellar disk and a longer time crossing the edge of the stellar disk (ingress/egress). The inverse is true for larger, Jupiter-like planets, which have long transit durations, on the order of several minutes, and comparatively shorter ingress/egress times. Jupiter, for instance, is about 8.5 times larger than the 0.01 $R_{\odot}$ white dwarf in the HT Cas binary observed by Horne et al. (1991) and Nucita et al. (2009). Because white dwarfs (and other compact objects) are so small, the ingress/egress times of both small and large planets should be comparable.

Simulated transit light curves for Earth- and Jupiter-mass planets in CB (P-type) orbits around CVs with different mass ratios $\mu$ are shown in Figures 6 and 7. The orbits are circular and have no inclination. In Figure 6, the white dwarf mass varies ($M_{1}=0.6$, 0.8, 1, and 1.2 $M_{\odot}$) in each of the panels, while the donor mass remains fixed at $M_{2}=0.1~{}M_{\odot}$. In Figure 7, the white dwarf mass is fixed at $M_{1}=0.8$ $M_{\odot}$ in each of the panels, while the donor mass varies ($M_{2}=0.05$, 0.1, 0.5, and 1 $M_{\odot}$). In each plot in Figures 6 and 7, the critical semimajor axis of the planet, (which depends on the mass ratio of the binary), sets the orbital period and therefore the transit duration.

4. Discussion

In this section we give examples of planetary transit light curves based on simple theoretical considerations. Our theoretical examples illustrate key features, but are not designed to capture the full complexity of real X-ray light curves, which may be affected by geometrical effects and time variability. Fortunately, there are many published examples of X-ray eclipses by stars, detected in spite of all of the complicating factors (e.g. van Teeseling, 1997; Ramsay et al., 2001; Porquet et al., 2005; Cherepashchuk et al., 2009; Nucita et al., 2009; Ponti et al., 2017; Kennedy et al., 2017). In most cases, a dip in the X-ray light curve that is shorter in duration and/or more shallow than the stellar eclipse would also be detectable. These systems therefore provide examples in which planetary transits can either be discovered, or else the size and orbit of any transiting planet can be constrained if no signal is detected.

4.1. Example light curves

In this section, we discuss the signatures we might expect in the light curves of CB planets in mass-transfer binary systems. In Figure 6 we plot light curves for a set of CVs in which the mass of the donor stays fixed at 0.1 $M_{\odot}$, and the white dwarf mass increases from 0.6 $M_{\odot}$ to near the Chandrasekhar limit at 1.2 $M_{\odot}$. As the white dwarf mass increases, $\mu$ decreases, and therefore the transit duration becomes increasingly short. The transit depth, $\delta=(R_{p}/R_{1})^{2}$, for Earth-like planets gets deeper as $M_{1}$ increases, while Jupiter-like planets entirely eclipse the white dwarf for all mass ratios. The former effect is due to the decrease in white dwarf size with increasing mass, according to the relationship $R_{1}\approx 0.01(M_{1}/M_{\odot})^{-1/3}R_{\odot}$. As the size of the white dwarf decreases—assuming the size of the X-ray emitting region in the CV is comparable in size to the white dwarf—an Earth-like (or larger) planet is capable of blocking out increasingly more of the flux, giving rise to deeper transit depths. For $M_{1}=0.6$, an Earth-like planet would block $\delta\approx 0.4=40\%$ of the X-ray flux for a few seconds, while for $M_{1}=1.2$, it would block roughly 95% of the flux for less time.

If we now consider CVs in which the white dwarf mass is fixed at a typical value of 0.8 $M_{\odot}$ and the donor mass increases, we find that the transit durations become increasingly longer, most noticeable for large, Jupiter-like planets (Figure 7). For a CV in which $M_{1}=0.8~{}M_{\odot}$ and $M_{2}=0.05~{}M_{\odot}$, a Jupiter-like planet would eclipse the white dwarf for roughly 4 minutes. If instead, $M_{2}=M_{\odot}$, a Jupiter would cover the white dwarf for roughly 9 minutes. This trend is due to the increasing $a_{c}$, which increases with the mass ratio of the binary. For smaller planets, the transit depth $\delta\approx 0.68$ remains fixed for all mass ratios; that is, roughly 68% of the X-ray flux would be blocked out by an Earth-like planet orbiting a CV with a 0.8 $M_{\odot}$ white dwarf.

We note that for small, Earth-like planets, the duration over which the transit is deepest can be very short, lasting only a few seconds (see Figure 5). This might raise doubts as to the possibility of detecting such small planets around XRBs. However, it should be kept in mind that the ingress/egress durations of small planets are longer, around $\sim 30$–45 seconds. Thus, the entire time scale of such a transit ($T+\tau_{\rm ing}+\tau_{\rm egr}$) will be $\gtrsim 1$ minute, an interval over which there could plausibly be detectable effects.

4.2. Probability of Transits

To estimate transit probabilities for circumbinary planets orbiting XRB systems, we follow Li et al. (2016), who provide analytical expressions for calculating transit probabilities involving eclipsing stellar binaries. The transit criterion can be separated into two parts. First, we calculate the geometric condition called transitability, in which the planet orbit intersects the stellar binary orbit. Transitability does not guarantee actual transits on every crossing of the planet past the binary orbit, due to the relative motion of the three objects (Martin & Triaud, 2015). Thus, our second step is to determine the planet transit probability, given that it intersects the binary orbit. We provide the following summary of the probability calculations and refer the reader to Li et al. (2016) for the full details.

Both the mutual inclination between the planetary and binary orbits, $\delta i$, and orbital precession can affect the transitability. Precession is described by time variation in the codependent quantities $\delta i$ and $\Omega$, the longitude of ascending node of the planet. Li et al. (2016) derive the probability for a planet to cross the binary orbit, $P_{\rm cr}$, as

[TABLE]

where $\delta\Omega=\Delta\Omega+\dot{\Omega}T_{\rm obs}$ is the total range of ascending node longitudes for the primary star, $\Omega$, during the total observation time, $T_{\rm obs}$. The total range of $\Omega$ over which the planet can cross the binary orbit is $\Delta\Omega$, which depends on the true anomaly of the planet, the planetary and binary semimajor axes, the size of the star being transited, the mutual inclination angle $\delta i$, and the inclination angle of the binary orbit, $i_{b}$ (equations $10$-$12$ in Li et al., 2016). We slightly modify equations (11) and (12) in Li et al. to include the size of the transiting planet, $R_{p}$. Li et al. do not include $R_{p}$ since, for the scenarios they consider, $R_{p}\ll R_{1}$. However, for XRB systems, the transiting planet may have a size comparable to or larger than the size of the X-ray emitting region and cannot be ignored. Thus, in Li et al.’s equations (11) and (12), we replace every occurrence of $R_{1}$ with $(R_{1}+R_{p})$ (Gongjie Li, priv. communication).

The time rate of change of $\Omega$ due to precession is related to the orbital precession timescale of the planetary orbit by (Schneider, 1994)

[TABLE]

where $P_{p}$ is the planetary period, and $a_{p}$ is the semi-major axis of the planet. In Figure 8, we plot $T_{\rm prec}$ as a function of $\mu$ for three scenarios, $m_{1}=0.8$, $1.4$, and $10$ $M_{\odot}$, for the coplanar case $\delta i=0\arcdeg$. Since the stellar masses used to calculate $P_{p}$ are much greater than the planet masses, the periods of Earth- and Jupiter-mass planets are indistinguishable in the plot. The precession time is given in terms of an observing time of $T_{\rm obs}=100$ ks (see below), and the figure demonstrates that for a typical observing window, the precession periods of the considered circumbinary systems are comparatively long, $T_{\rm prec}>T_{\rm obs}$. This confirms that precession should be taken into account for calculations of the XRB transit probability.

Next, we calculate the probability of transit, $P_{\rm tran}$, given that the planet and binary orbits intersect. Li et al. (2016) show that $P_{\rm tran}$ depends on the ratio of the relative displacement of the planet and star with respect to the observer, $dl$, and the projected width of the stellar orbit. The relative displacements are

[TABLE]

where $v_{p,1}$ are the orbital velocities of the planet and primary star, and $t_{\rm tran}=\rm{min}[(\pi(R_{1}+P_{p}))/(2v_{p}\sin\delta i),T_{\rm obs}]$ is the time it takes for the planet to cross the stellar orbit. When the planet and primary star are both closer to the observer with respect to the center of mass, the relative displacement is $dl=dl_{1}$. When the star is on the other side of the center of mass, the planet and star move in opposite directions, and $dl=dl_{2}$.

With these definitions, Li et al. (2016) derive the probability for the planet to transit the primary star at least once in a finite observing window, of length $T_{\rm obs}$, as

[TABLE]

where

[TABLE]

and where $2a_{b}\mu$ is the projected width of the binary orbit. The exponent in equation (11), $n$, is the number of times the planet crosses the binary orbit during a finite observing window. The number of crossings is given by a conditional expression that depends on $T_{\rm obs}$, the planetary orbital period $P_{p}$, and precession (equation 18 in Li et al.).

Figure 9 displays the transit probabilities as a function of mass ratio, $\mu$, for a planet orbiting three different binary types: $m_{1}=0.8$ $M_{\odot}$ and $R_{1}=0.01~{}R_{\odot}$ (solid lines); $m_{1}=1.4$ $M_{\odot}$ and $R_{1}=10$ km (dashed lines); and $m_{1}=10$ $M_{\odot}$ and $R_{1}=100$ km (dotted lines). We assume a total observing time of $T_{\rm obs}=100$ ks. In the left panel, probabilities are calculated for coplanar systems, i.e., $\delta i=0\arcdeg$. In the right panel, we assume a mutual inclination of $\delta i=2\arcdeg$. This value is motivated by observations of circumbinary planets detected by Kepler, which exhibit a mean value of $1\fdg 7$ (Martin & Triaud, 2015; Li et al., 2016). For both panels, we assume an orbital inclination for the binary of $i_{b}=80\arcdeg$. Modeling of HT Cas suggests that the systems has an inclination of $81\arcdeg$ (Horne et al., 1991).

Figure 9 shows that for fixed values of $\mu$, $P_{\rm tran}$ increases as the mass of the primary star and size of the X-ray emitting region decrease. Thus, systems representing CVs, with $m_{1}=0.8$ $M_{\odot}$, have the highest transit probabilities, while systems representing $10$ $M_{\odot}$ blacks holes and their companions have the lowest. When mutual inclination effects are included, the relative behavior of $P_{\rm tran}$ is more complex. In this example, increasing $\delta i$ by only $2\arcdeg$ causes $P_{\rm tran}$ to decrease for most realistic values of the mass ratio, $0.05\lesssim\mu\lesssim 0.6$. However, increasing $\delta i$ even more would lead to an increase in the transit probability. This is due to the competing effects between $T_{\rm obs}$, $dl_{1,2}$, and $\Delta\Omega$, all of which have different functional dependencies on $\delta i$. As an example, if HT Cas ($\mu=0.19$, $i_{B}\approx 80\arcdeg$) has a circumbinary Earth- or Jupiter-like planet at the stability limit, the transit probability would be $\lesssim 20\%$.

In Figure 10, we display plots of the transit probability, assuming that the binary orbital plane is completely edge-on with respect to the observer, i.e., $i_{b}=90\arcdeg$. In these highly optimistic scenarios, values of $P_{\rm tran}$ are much higher than those for the examples depicted in Figure 9, where $i_{b}=80\arcdeg$. Again, we see that a small increase in the mutual inclination between the binary and planetary orbital planes tends to cause a decrease in $P_{\rm tran}$. In both Figures 9 and 10, this decrease in $P_{\rm tran}$ is especially pronounced for Earth-like planets (black curves).

In general, our calculations show that the transit probabilities for potential circumbinary planets in XRB systems are favorable for detection. Taking the right panel in Figure 9 to represent the most conservative and plausible scenarios ($\delta i=2\arcdeg$ and $i_{ib}=80\arcdeg$), planets orbiting CVs may have transit probabilities in the range $2.5\%\lesssim P_{\rm tran}\lesssim 40\%$. Planets in neutron star and black hole systems may have transit probabilities of $0.9\%\lesssim P_{\rm tran}\lesssim 20\%$ and $0.1\%\lesssim P_{\rm tran}\lesssim 4\%$, respectively, depending on the mass ratio of the system.

As indicated in Figures 9-10 and in equation (8), the transit probability depends on the total observation time, $T_{\rm obs}$. We note that $T_{\rm obs}$ need not be continuous. Rather, we are concerned with the total duration of all of the observations of a given X-ray source for which archived data is available. The archived resources are vast, including, from the present generation of X-ray missions alone, hundreds of mega-seconds of data. The distribution of exposure times within a given year provides a general guideline. For example, observing programs approved by Chandra during Cycle 20111The Chandra web page with statistics for proposals: http://cxc.harvard.edu/target_lists/cycle19/cycle19_peer_results_ stats.html. included 171 targets with exposures of 16 ks or less, with 120 of these having exposure times between 8 ks and 16 ks. In addition, there were 160 targets with exposures having durations ranging from 16 ks to roughly 1 Ms. Many systems of special interest have been targeted for long observations, favorable for the discovery of transits. These include some individual X-ray binaries, such as IC 10 X-1, an eclipsing accreting black hole in a Local Group dwarf galaxy. They also include globular cluster such as M15, and external galaxies located within 10 Mpc, such as M51, M82, and M101.

Finally, we comment on our choices for $R_{1}$, in equation (8) and Figures 9–10, necessary to calculate $P_{\rm tran}$. The sizes of X-ray emitting regions depend on the nature of the compact accretor and also on the details of the accretion and emission processes. For instance, for a white dwarf that has undergone a nova explosion, soft X-ray emission emanates from the surface of the star, and the effective radius of emission is comparable to its radius, typically on the order of $0.01\,R_{\odot}$. Several luminous supersoft X-ray sources (SSSs) in the Magellanic Clouds and Milky Way are thought to be quasisteady-burning white dwarfs, and their X-rays are emitted from an effective area comparable to that of a white dwarf. The SSS class also includes white dwarfs post-nova (Henze, 2016). Neutron stars can emit from their surfaces or, in some states, from the regions around the magnetic poles (Wang et al., 2017). Emission can also come from a somewhat larger accretion disk, which is hottest near the inner radius. Urquhart & Soria (2016) have used long exposures taken of M51 by both Chandra and XMM-Newton to discover X-ray transits of two separate X-ray sources, likely to contain black hole accretors. These particular systems can exhibit emission from an inner disk having dimensions equal to that of the last stable orbit, 10s to 100s of kilometers for stellar mass black holes, and scaling up for black holes of larger mass (Soria et al., 2017).

In sum, if Earth- and Jupiter-size planets exist in circumbinary XRB systems, and if their orbits are close to the minimum values set by stability requirements, then the transit probabilities show promise for detection. A planet orbiting a single star has a transit probability of $(R_{1}+R_{p})/a_{p}$, and it is often the case that $R_{1}\ll a_{p}$, corresponding to a small value for the transit probability. In an XRB, however, the semimajor axis of a planet may be so small that—despite the emitting region in the XRB also being potentially very small (see §1.2)—the ratio determining the transitability probability (equation 8), and therefore $P_{\rm tran}$, is comparatively much larger. Other factors, such as inclination and precession, may also favor high transit probabilities.

We reiterate that our calculations assume that planets are near the stability limit. If planets happen to be in wider orbits, this would decrease the transit probability. Although there would be slight modifications in our calculations, the transit probability would essentially fall off as the inverse of the planetary separation, $a_{p}$ (Gongjie Li, priv. communication). Thus, in our example above for HT Cas, if a planet was in an orbit at a radius $N$ times the critically stable orbit, the transit probability would be $\lesssim(1/N)\times 20\%$.

In the next section, we discuss conditions that influence the detectability of transiting planets.

4.3. Detectability and Identifying Candidate Systems

For any given XRB with transiting planets, the identification of transits depends on (1) the baseline, i.e., the out-of-eclipse background-subtracted count rate of the source, $dc/dt$; (2) the duration, $\Delta\,t$, of the occultation; and (3) the depth of the occultation. The depth of the occultation depends on the relative surface areas of the source and planet, as well as on the inclination of their orbit. Let ${\cal A}(t)$ be the area of the source which is covered by the planet at some time $t$ after the first contact has been made and before the last contact is made. The average value of the area covered is ${\cal A}_{\rm cov}=\langle{\cal A}(t)\rangle$.

The number of counts “missing” during the transit is then

[TABLE]

where $R_{\rm em}$ is the size of the X-ray-emitting region.

Exact expressions for $\Delta\,t$, which include the geometrical effects, are given in § 2. Here we express $\Delta\,t$ simply as $2\,R_{\rm eff}/v$, where $v$ is the orbital speed during the eclipse, and $R_{\rm eff}$ is an effective radius.

[TABLE]

where $M_{\star}$ is the total mass enclosed by the orbit, and $R_{J}$ is the radius of Jupiter. In the limit of a small source and/or a large planet, $R_{\rm eff}$ may be taken to be the radius of the planet, $R_{p}$. More generally, $R_{\rm eff}=f\times(R_{p}+R_{\rm em})$, where $f$ is a number (generally $f\lesssim 1$) which takes into account orientation and ingress/egress effects.

Equation (13) shows that, for a given count rate, more counts are missed during a transit for larger orbital separation and smaller total mass. These same factors yield longer orbital periods, however, producing fewer opportunities to detect additional transits that could validate the interpretation that a dip in the count rate was associated with a transit. Thus, for any given system, there is an optimal compromise between larger $a_{b}$ and smaller $M_{\star}$, which promotes the detection of a single passage, and smaller $a_{b}$ and larger $M_{\star}$, which may make it possible to detect multiple transits.

To set scales relevant for the analysis of archived data, we consider the distribution of background-subtracted count rates obtained during observations with the Chandra X-ray observatory (Wang et al., 2016). There were 772,000 X-ray states were found, and roughly 1000 (6000, 21,000) of them exhibited count rates larger than 1 s*-1* (0.1 s*-1*, 0.01 s*-1*).222 Wang et al. (2016) analyzed all of the Chandra data collected during the years $\sim$ 1999-2014. Many sources were detected multiple times, so that the total number of distinct states observed is about $3$ times larger than the total number of distinct sources. Many of the X-ray sources producing these states are XRBs. With the XMM-Newton observatory, count rates of the same sources would generally be several times larger, depending on the observational set up (Snowden, 2002). However, the number of background counts per unit time would also be larger.

For the highest count rates, a total eclipse lasting a minute would be readily detected; for the lowest of these count rates, the passage might have to last for tens of minutes to assure the reliable detection of a single planetary passage. Repeated transits could be detected, allowing better analysis of the dips, if the total exposure time, generally comprised of a sum of exposure times from separate observations, encompasses multiple orbits. Some regions of the sky contain specific individual XRBs of high interest (e.g., black hole candidates), and have been imaged many times. In addition, dozens of galaxies have had total exposure times in the range of 1-2 days, with some sources (e.g., in M31) having had more coverage. It therefore seems well within the realm of possibility to detect planetary transits.

4.4. Physical Systems

Already-discovered eclipsing XRBs can be used to guide searches through archival data for planetary transits for a number of reasons. First, these systems provide specific examples of a compact source being obscured, in that each exhibits a measured eclipse depth and duration. By studying the signal-to-noise and other characteristics of the eclipse, we can determine how small and how fast-moving a planet orbiting the X-ray source could be and still be detected. Such feasibility studies will also inform analyses of other XRBs which do not exhibit stellar eclipses.

Second, because stellar and planetary orbits may be coplanar, an XRB in which there is a stellar eclipse may be more likely to also exhibit planetary eclipses, should that system host planets. Thus, XRBs with stellar eclipses may be the most promising systems in which to search for planetary transits. The period of planetary transits would be different from the binary period; the stability issues we discuss here would provide useful guides.

Finally, a stellar eclipse establishes that the X-ray source is compact. For example, the long-lasting X-ray eclipse in the black-hole/Wolf-Rayet binary X1 in the dwarf galaxy IC 10, illustrates that some black hole binaries do have compact regions emitting distinctive X-rays. Similarly, several eclipsing white dwarfs are known, where X-rays emanate from near the surface of the white dwarf (e.g., Nucita et al., 2009). These indicate that it could be productive to search for planets orbiting white dwarfs that have recently experienced novae. Post-nova, many systems pass through a supersoft phase in which X-rays with energies $<1$ keV are emitted from near the white dwarf’s surface. Such searches may find planets in P-type orbits around close-binary novae and in S-type orbits in symbiotic novae, which generally have wide-orbit giant donors.

Equation (13) allows an estimate of the count rate necessary to allow transits to be identified. Count rates that are large enough for this purpose can be provided by nearby XRBs or by distant yet bright XRBs. Consider a source with X-ray luminosity equal to $10^{39}$ erg s*-1* at $5$ Mpc; this corresponds to an ultraluminous XRB in an external galaxy. If the X-ray spectrum is a power law of index $1.7$, XMM-Newton would collect $0.1$ counts per second (cps).333HEASARC PIMMS predicts $0.1$ cps with XMM PN THIN (0.4–10 keV) for an on-axis observation and $N_{H}=10^{21}$ cm*-2*. (https://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3pimms/w3pimms.pl.) This can be sufficient for the discovery of transiting planets. The same flux would be obtained from a source of $10^{33}$ erg s*-1* at $5$ kpc; this could correspond to a quiescent low-mass X-ray binary (LMXB) in a Galactic globular cluster, while the flux from a bright globular cluster LMXB would be $1000$ times larger. Alternatively, a very dim source with X-ray luminosity of $10^{29}$ erg s*-1* at $50$ pc, would provide the same flux.

As discussed in §4.2, the probability $P_{\rm tran}$ that the orientation is favorable for transit detection declines with increasing orbital separation. If an observation capable of detecting transits lasts over a time interval of one or more orbital periods, the probability remains equal to $P_{\rm tran}$. Even if a sequence of shorter observations are spread out over a long interval, the probability is roughly equal to $P_{\rm tran}$ as long as the total observation time is longer than an orbital period. On the other hand, when the orbital period is longer than the total observing time, then the likelihood that the time interval of an observation will include a transit depends on the fraction of the time the system is in transit, the transit duty cycle: $\sim R/(\pi\,a)$.

The question of how the duty cycle affects detectability primarily hinges on the typical time duration of the X-ray observations relative to the orbital period. Many X-ray sources that are promising targets for transit searches have been observed for hundreds of kiloseconds, including various regions in M31 (e.g., Vulic et al., 2016). The same is true of roughly 10 nearby galaxies, in which the observing footprint covers most of the galaxy, and with each galaxy typically containing dozens of bright sources (e.g., Jordán et al., 2004; Sell et al., 2011). Finally, Galactic globular clusters and many individual sources in our galaxy and the Magellanic Clouds have been targets of long-term monitoring by several different X-ray observatories.

These exposures are long enough to ensure that transits occurring in systems with orbital periods up to tens of hours have a probability $P_{\rm tran}$ of being observed. Planets in wider orbits could produce detectable transits—and, in fact, the transit duration would be longer, enhancing the signal. But the probability would be reduced as described above.

We don’t know the orbital period distribution of planets in X-ray binaries. However, among known exoplanets discovered via the transit method, roughly one in eight have orbital periods shorter than two days.444See exoplanets.eu We can gain insight about possible X-ray binary planetary systems from observations and calculations of the evolution of planetary orbits during stellar evolution. Planets in XRB systems would likely have been part of the planetary system of one of the binary’s original component stars. Dynamical interactions during the binary’s evolution would have ejected some planets and placed others into closer orbits. Recently, low-mass objects have been discovered orbiting WDs. For example, eclipses of a WD (SDSS J141126+200911) at optical wavelengths identify an object interpreted as a brown dwarf in a 121.73 minute orbit, producing eclipses 125 s in duration (Beuermann et al., 2013). Such discoveries focus attention on the question of how such systems evolve. Calculations of the tidal interactions between evolving stars and their low-mass companions, and calculations of the possible engulfment of the low-mass companions find that there should be a period gap in the post-evolution orbits, with brown dwarfs able to occupy orbits with periods smaller than $0.1$ day (Nordhaus et al., 2010). Planets may come to be in close orbits around WDs if they survive engulfment, or if their initial orbits were wide enough to avoid engulfment, but they were subsequently scattered into close orbits (Nordhaus & Spiegel, 2013).

4.5. Search for extraterrestrial intelligence

An interesting application of our proposal to search for planets around XRBs is that the same technique could be used in the search for extraterrestrial intelligence (SETI), since if advanced civilizations exist, they may produce anomalies in astronomical systems. One possible anomaly could be caused by placing around a star a Dyson sphere, a surface designed to capture a star’s energy and convert the starlight into far-infrared emission (Dyson, 1960). It has also been suggested that smaller, multiple structures could be placed in orbit around a star (e.g., Harrop & Schulze-Makuch, 2010). Should they exist, some such artificial structures might caused detectable transits.

Such artificial structures are often discussed in conjunction with the possible energy needs of an advanced civilization, and because any convincing discovery of them would provide evidence of exocivilizations. Here we make a possible connection to XRBs, which are (1) highly luminous and therefore potentially desirable as sources of energy, and also (2) distinctive and rare enough to be targeted for X-ray observations by emerging technological civilizations.

We therefore suggest that searching for transits around XRBs could—in addition to or instead of leading to evidence of natural planets—lead to evidence of artificial structures. While it is not clear that we could distinguish between natural and artificial structures, one could attempt to answer such questions as: Are the planetary orbits stable? Does the light curve, especially during ingress and egress, suggest an unusual shape?

We also make a particular connection to globular clusters, which have been suggested as environments within which advanced civilizations may be able to thrive because of (1) the absence of young stars and the attendant risks, such as supernovae; (2) the stability of orbits in the habitable zones of the low-mass stars comprising the primary cluster populations; and (3) the relatively high stellar density, which could facilitate interstellar travel, should technologically advanced civilizations inhabit globular clusters (Di Stefano & Ray, 2016). XRBs are rare, but are more common in globular clusters than in the rest of the Galaxy (Clark, 1975; Di Stefano & Rappaport, 1994). Furthermore, they could be significant sources of energy in an environment dominated by M dwarfs, while at the same time lying within a relatively small distance of any point in the cluster, due to the limited (generally $<10$ pc) size of globular clusters.

5. Summary

We have presented a novel idea for searching for exoplanets around mass transfer X-ray binaries. For demonstration of concept, we have focused on CVs, in which the typical mass of the white dwarf is $\sim 0.8$ $M_{\odot}$ and the donor is a main-sequence star or a brown dwarf. If a planet is associated with such a binary, it is most likely to have a circumbinary (P-type) orbit. While, for a mass ratio of $0.5$, the closest possible orbit would exhibit transits of just under a minute for an Earth-size planet and about $10$ times longer for a Jupiter-size planet, the transit durations would be larger in wider orbits. We have shown how these quantities scale for accretors of higher mass (neutron stars and black holes) and for a range of mass ratios. If the binary has a wider orbit, then either the compact object or the donor star, may support a planet in an S-type orbit.

We have demonstrated that a number of factors may favor high transit probabilities form potential circumbinary XRB planets. If planets are in orbits near their critical semi-major axes, and if inclination and precession are important, then transit probabilities may range from roughly 0.1% to 40%, depending on the mass ratio of the binary and size of the X-ray emitting region. Given our probability calculations, we have also shown that transits of planets orbiting XRBs, including symbiotic binaries, are detectable with present-day technology. Future observations with missions such as the proposed Lynx X-ray mission,555https://wwwastro.msfc.nasa.gov/lynx/ which may have $50$ times the sensitivity of Chandra, will dramatically extend the opportunities for transit detection. The discovery of a planetary transit in X-ray binaries, while not assured, would be of major scientific importance. Any well-designed archival study will either discover planets or else will place limits on planetary orbits. If and when the first planet is detected around a mass transfer binary, follow-up observations could be conducted with other telescopes to constrain its atmospheric properties and to search for biosignatures.

From Nia Imara: I dedicate this paper to my grandmother, Carrie Antonia Lee, and to Lewis Bundy. I thank Matthew Holman, Gongjie Li, and Jason Wright for discussions and comments that helped improve this paper. This work is supported by the John Harvard Distinguished Scientist Fellowship.