2041-8205/853/2/L27

Previously, planets have been detected only in the Milky Way galaxy. Here, we show that quasar microlensing provides a means to probe extragalactic planets in the lens galaxy, by studying the microlensing properties of emission close to the event horizon of the supermassive black hole of the background quasar, using the current generation telescopes. We show that a population of unbound planets between stars with masses ranging from Moon to Jupiter masses is needed to explain the frequent Fe K α line energy shifts observed in the gravitationally lensed quasar RXJ 1131–1231 at a lens redshift of z = 0.295 or 3.8 billion lt-yr away. We constrain the planet mass-fraction to be larger than 0.0001 of the halo mass, which is equivalent to 2000 objects ranging from Moon to Jupiter mass per main-sequence star.

Export citation and abstract BibTeX RIS

1. Introduction

Over the past two decades, it has been established that planets are ubiquitous in the Milky Way galaxy (e.g., Wolszczan & Frail 1992; Mayor & Queloz 1995; Udry & Santos 2007; Lissauer et al. 2014; Winn & Fabrycky 2015). Extrapolating to the extragalactic regime, it is natural to hypothesize that planets are common in external galaxies as well. However, we lack the observational techniques to test this hypothesis because compared to their Galactic brethren, extragalactic planets are much farther away and much more difficult to separate from the host stars/galaxies. Just as gravitational microlensing provides a unique tool to detect planets in the Galaxy (e.g., Mao & Paczynski 1991; Gould & Loeb 1992; Gaudi 2012), it can also provide the capability to detect planets in extragalactic galaxies by combining microlensing and a galaxy-scale gravitational lens. We are interested in quasar–galaxy strong lensing systems, where a background quasar is gravitationally lensed by a foreground galaxy and multiple images of the quasar form (Walsh et al. 1979). Light from these quasar images crosses different locations of the foreground galaxy and is further lensed by nearby stars in the region in the lens galaxy. This effect is called quasar microlensing (Wambsganss 2006; Kochanek et al. 2007) and has been used extensively to measure the structure of the quasar accretion disk around the supermassive black hole (SMBH) at the center (e.g., Kochanek 2004; Dai et al. 2010; Chen et al. 2011, 2012; Mosquera et al. 2013; Chartas et al. 2017; Guerras et al. 2017) and the properties of mass distributions in the lens galaxy (e.g., Morgan et al. 2008; Bate et al. 2011; Blackburne et al. 2014). As we probe smaller and smaller emission regions of the accretion disk close to the event horizon of the SMBH, the gravitational fields of planets in the lensing galaxy start to contribute to the overall gravitational lensing effect, providing us with an opportunity to probe planets in extragalactic galaxies. An important length scale in gravitational lensing is the Einstein ring size. For a point mass, the Einstein ring in the source plane is where D ol , D os , and D ls are the angular diameter distances between the observer, lens, and source, respectively. For a typical lens redshift of z l = 0.5 and source redshift of z s = 1.5, the Einstein ring for an Earth mass object is If an emission region is smaller or comparable to this Einstein ring size, the emission will be significantly affected by the microlensing effect, where the microlensing flux magnification can be a factor of several or higher. The Schwarzschild radius of a typical 108 black hole is which is comparable to the Einstein ring size of Earth-size planets. Therefore, emission close to the Schwarzschild radius of the SMBH in the central engine of the source quasar will be affected by planets in the lensing galaxy. In this Letter, we show that the lensing effects from planets can explain some of the observational data of the gravitational lens RXJ 1131–1231. Throughout the Letter, we assume a flat cosmology with the cosmological parameters , , and .

2. Data

RXJ 1131–1231 is a quadruple lens systems with source and lens redshifts of z s = 0.658 and z l = 0.295, respectively (Sluse et al. 2003). The central black hole mass is measured to be (Dai et al. 2010), and thus the hole has a gravitational radius cm. The Einstein ring radius in the source plane for objects in the lensing galaxy is cm or cm. Thus, the X-ray emission especially the reflection components within ~10 r g around black holes can be significantly affected by the lensing effect of planet-size objects. During the past decade, RXJ 1131–1231 has been observed by the Chandra X-ray Observatory for 38 times, and we have detected microlensing signatures in the reflection component of the X-ray emission in the source (Chartas et al. 2009, 2012, 2017), as illustrated in Figure 1 (left). First, the Fe Kα line in the reflection component is measured to have large blueshifted or redshifted peak energies compared to the rest-frame peak energy of with ranging from 0.5 to 1.3, and in some cases double lines are detected. Figure 1 (right) shows the distribution of the measured energy shifts in RXJ 1131–1231 (Chartas et al. 2017). Although the relativistic Fe Kα line can peak at a range of energies from ~5 to 8 keV, depending on the emissivity profile, black hole spin, and observer's viewing angle (e.g., Brenneman & Reynolds 2006), the line is never observed to vary to this large range with g changing by a factor of two for a single object. For example, in the frequently observed Seyferts, such as MCG−6–30–15, NGC 4151, and MCG−05–23–16, the peaks of the Fe Kα lines are measured to be constant (Kara et al. 2014; Beuchert et al. 2017; Wang et al. 2017). Second, the line energy variations are detected with a very high frequency in the 38 total observations in all four images. Other non-microlensing interpretations of these line shifts including emission from patches of an inhomogeneous disk or ionized accretion disk, or intrinsic absorption or occultation, are unlikely (Section 5 of Chartas et al. 2017). We calculate the number of epochs where a line energy shift is observed, based on Table 5 of Chartas et al. (2017), using the following steps. We first identify the peaks of the g distributions for the four images, g = 1 for images A and C and g = 0.85 for images B and D (Figure 1, right), and then for each detected line we calculate the significance of the line energy shift from the peak based on its measured energy and uncertainties, i.e., . When two lines are detected for the same image in the same observation, we only use the stronger line in the analysis. We find 13, 16, 11, and 5 epochs with shifted lines of more than energy shifts, corresponding to rates of 34%, 42%, 29%, and 13% for a total of 38 observations for images A to D, respectively. For energy shifts, the corresponding numbers of epochs are 10, 13, 8, and 4 with rates of 26%, 34%, 21%, and 11% for images A to D. Figure 1. (Left) Simulated redshift/blueshift distribution of an emission region, with a power-law emissivity index of n = 5.8, close to a rapidly spinning supermassive black hole with a = 0.9 at the quasar center, viewed at an inclination angle of 70°. The color bar below indicates the color-coded energy shift values, , produced by the general and special relativistic effects by the black hole. The green line shows a microlensing caustic (high magnification regions and discontinuities of the magnification map) on the emission region, preferentially magnifying a portion of the emission region with a specific subset of g values, which will result in a line energy shift. (Right) Distributions of the Fe Kα line energy shifts measured in the four quasar images of RXJ 1131–1231 (Chartas et al. 2017). We choose a bin size of 0.05 for illustration purposes, and the choice of bin size affects little the line shift rates estimated in this Letter. Download figure: Standard image High-resolution image Export PowerPoint slide The combination of these unique features has never been observed in a non-lensed AGN, while two other lensed quasars have similar line energy shifts reported, but with less monitoring epochs (Chartas et al. 2017). The line energy shift will occur when a microlensing caustic lands on the line emission region as illustrated in Figure 1 (left) and preferentially magnifies a portion of the line emission with different observed energies due to the special and general relativistic effects (Popović et al. 2006; Krawczynski & Chartas 2017). However, explaining the high frequency of line energy shifts imposes an additional difficulty—the high frequency implies a large number of additional microlenses.

3. Microlensing Analysis

We perform a microlensing analysis focusing on the occurrence of line energy shifts, i.e., the frequency of caustic crossing events for the Fe Kα emission region of RXJ 1131–1231. To simplify the calculations, we assume that the line energy shift is occurring when discontinuities (caustics) in the microlensing magnification maps land on the source region, such that a portion of the disk with different g factors is magnified differently. If the source is located inside of a highly magnified but smooth region in the magnification map, the line shape will not be significantly distorted because all the emission regions with different g factors are magnified similarly. The microlensing magnification maps are generated with the inverse polygon mapping algorithm (Mediavilla et al. 2011). We focus on the three brighter images A, B, and C in this analysis because the line shift frequencies measured in the faintest image D is subject more to selection effects due to the sensitivity limits of the observations. The macro model parameters are adopted from the most-likely model of Dai et al. (2010), with the the global convergence and shear as , 0.53, and 0.55 and , 0.41, and 0.39, respectively, for images A, B, and C. The stellar population is modeled with a broken power-law mass function in the range of 0.05 to 2 with the break at 0.5 , the low and high power-law indices are and , and the normalizations are set such that the surface mass fraction in stars , 0.10, and 0.10 for images A, B, and C. This stellar population and global convergence and shear parameters are typical of the models used in previous quasar microlensing calculations for this system (e.g., Dai et al. 2010). The brown dwarf population is approximately modeled within the 0.05–0.08 mass range, which provides a consistent fractional contribution to the total stellar/brown dwarf population, compared to the recent constraint of Sumi et al. (2011) with the mass range between 0.01 to 0.08 and a power-law slope of 0.5. Figure 2 (left) shows the simulated microlensing magnification pattern for image A of RXJ 1131–1231, and we can see that the density of microlensing caustics is too low to explain the ~30% rate of line energy shifts observed in RXJ 1131–1231 for a compact source. Figure 2. (Left) Microlensing magnification map of RXJ 1131–1231A with only stars and a dimension of (8000 r g )2. A random track (3740 r g ) for a compact source moving across the map with a 10 year duration is overplotted. The probability for a caustic to land on the source region is only a few percent, significantly below the observed rate of ~30%. (Right) Magnification map, with a dimension of (400 r g )2, with the additional planet population and a planet mass-fraction of for RXJ 1131–1231A. The caustic density is much higher with the additional planets. Download figure: Standard image High-resolution image Export PowerPoint slide We next add smaller bodies, planets, to the lens population. Since planets bound to stars will not significantly alter the magnification map because they are located far inside the Einstein ring of their parent stars (Section 4), we focus on free floating planets (Sumi et al. 2011; Strigari et al. 2012). We model the mass distribution of the floating planets ranging from Moon mass to Jupiter mass by a power-law model, , with the index fixed at t = 2 (Strigari et al. 2012), and the normalization , the fraction of planet surface mass density to the total surface mass density, is left as a free parameter. We generate a range of magnification maps with ranging from 0.0001 to 0.001, equivalent to 103–104 planets in the Moon to Jupiter mass range per main-sequence star, where the upper boundary is selected to match the recent limit of floating Jupiters in the Milky Way (Mróz et al. 2017). Since it is computationally expensive to calculate large magnification maps with a huge number of lenses, we constrain the size of the maps to be with each pixel 0.375 r g (e.g., Figure 2, right), and for each set of parameters, we generate 30 random maps to sample the large-scale variation of the magnification pattern. Figure 2 shows one example of the magnification map with the additional planets. Discontinuities in the magnification maps are found by convolving the map with a sharpening kernel, where n is the source size in pixels, the central n × n pixels of the kernel have the value X, and the remaining background pixels have the value −1. X is set to be , such that when in the magnification map, the convolution will result in a positive value, and the remaining pixel will be zero or negative. The adaptive source size n takes into account of the finite source size effect, which smooths the magnification pattern. We also produce convolutions with n × n constant kernels and compare with the sharpened magnification patterns, and Figure 3 shows two examples. We can see that the sharpened maps capture the main characteristics (peak and troughs) of the maps from convolutions with constant kernels. We then calculate the ratio of positive values over all valid pixels in the sharpened map to estimate the model probability of magnification discontinuity landing on the source region. We use the 30 maps for each set of input parameters to estimate the variance of the model probabilities. We set a prior size limit of r g for the Fe Kα emission region because first the large energy shifts can only be possible if the emission region is close to the black hole, where the general and special relativistic effects are large, and second, several studies suggest that the emissivity profile of the Fe Kα region is steep and more compact than the X-ray continuum, which is measured to be ~10 r g (e.g., Fabian & Vaughan 2003; Chen et al. 2012). Figure 3. Comparisons between the convolution of the magnification map with a constant kernel (left in each panel) and with the sharpening kernels (right in each panel) defined in this Letter, showing that the sharpening kernel has successfully captured the characteristics of the magnification maps for different source sizes. The left panel shows the comparison between the convolutions for a 4 pixel source size, and the right panel shows the comparison for a 16 pixel source size. Download figure: Standard image High-resolution image Export PowerPoint slide Figure 4 compares the microlensing model predictions with the observed rates of Fe Kα line energy shifts in images A, B, and C in RXJ 1131–1231. The model predictions depend on two parameters, the size of the Fe Kα emission region and the planet surface mass fractions. We can clearly see that microlensing models with only stars are significantly ruled out for any source size considered, and the statistical significances are more than 4.5, 6.4, and 5.6σ for images A, B, and C, respectively. With the additional lenses from planets, the microlensing model predictions increase and match the observed rates for large source sizes considered ~10 r g . We set the confidence limit on the lower limit of the planet mass-fraction, using the 3σ observed rates for the three images. The model with a planet surface density is excluded by 0.84, 1.73 and 1.11σ in images A, B, and C, respectively. Combining the probabilities, we find that the model with is excluded by 99.9%. Thus, the combined probabilities favor models with planet mass-fractions larger than . Figure 4. Model probabilities of observing an Fe Kα line energy shift in images A, B, and C of RXJ 1131–1231 as a function of source size. Different symbols show the model probabilities for different planet mass-fractions, , and the error bars show the model uncertainties. The dotted and dashed lines show the observed line energy shift rates of 1σ and 3σ confidence levels. The model with lenses composed of only stars is ruled out by more than 4.5, 6.4, and 5.6σ for images A, B, and C, respectively, for all source sizes. Download figure: Standard image High-resolution image Export PowerPoint slide

4. Discussion

We have shown that quasar microlensing can probe planets, especially the unbound ones, in extragalactic galaxies, by studying the microlensing behavior of emission very close to the innermost stable circular orbit of the supermassive black hole of the source quasar. For bound planets, they do not contribute much to the overall magnification pattern in this study. The Einstein ring size on the lens plane is proportional to . While the ratio of is of the same order for Galactic and extragalactic microlensing, the lens distance here is at a cosmological distance. Thus, the bound planets are located quite inside the Einstein ring of the parent star and do not significantly change the magnification patterns. Similarly, the star multiplicity factor is unimportant here because the Einstein ring is much larger in the extragalactic situation, and multiple star systems can be treated by a single stars. The unbound planet population is very difficult to constrain, even in the Milky Way galaxy. Using the density of caustics in the magnification pattern, we are able to constrain the surface mass density of planets with respect to the total mass as and the planet to star mass ratio . This planet to star mass ratio is equivalent to objects per main-sequence star in the mass range between Moon and Jupiter, or objects in Mars to Jupiter range including 0.08 Jupiters. This constraint is consistent with the upper end of the theoretical estimate of ~105 between Moon and Jupiter (Strigari et al. 2012), and the recently observed constraint, 0.25 Jupiters per main-sequence star, in the Milky Way (Mróz et al. 2017). It is possible that a population of distant but bound planets (Sumi et al. 2011) can contribute to a significant fraction of the planet population, which we defer to future investigations. Because of the much larger Einstein ring size for extragalactic microlensing, we expect that two models, the unbound and the distant but bound planets, can be better distinguished in the extragalactic regime.