Issue 
A&A
Volume 611, March 2018



Article Number  A63  
Number of page(s)  11  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201629676  
Published online  27 March 2018 
Lightcurve analysis of KOI 2700b: the second extrasolar planet with a cometlike tail
Astronomical Institute, Slovak Academy of Sciences,
059 60
Tatranská Lomnica, Slovakia
email: zgarai@ta3.sk
Received:
9
September
2016
Accepted:
3
December
2017
Context. The Kepler object KOI 2700b (KIC 8639908b) was discovered recently as the second exoplanet with a cometlike tail. It exhibits a distinctly asymmetric transit profile, likely indicative of the emission of dusty effluents and reminiscent of KIC 12557548b, the first exoplanet with a cometlike tail.
Aim. The scientific goal of this work is to verify the disintegratingplanet scenario of KOI 2700b by modeling its light curve and to put constraints on various tail and planet properties, as was done in the case of KIC 12557548b.
Methods. We obtained the phasefolded and binned transit light curve of KOI 2700b, which we subsequently iteratively modeled using the radiativetransfer code SHELLSPEC. We modeled the cometlike tail as part of a ring around the parent star and we also included the solid body of the planet in the model. During the modeling we applied selected species and dust particle sizes.
Results. We confirmed the disintegratingplanet scenario of KOI 2700b. Furthermore, via modeling, we derived some interesting features of KOI 2700b and its cometlike tail. It turns out that the orbital plane of the planet and its tail are not edgeon, but the orbital inclination angle is from the interval [85.1, 88.6] deg. In comparison with KIC 12557548b, KOI 2700b exhibits a relatively low dust density decreasing in its tail. We also derived the dust density at the beginning of the ring and the highest optical depth through the tail in front of the star, based on a tailmodel with a crosssection of 0.05 × 0.05 R_{⊙} at the beginning and 0.09 × 0.09 R_{⊙} at its end. Our results show that the dimension of the planet is R_{p}∕R_{s} ≤ 0.014 (R_{p} ≤ 0.871 R_{⊕}, or ≤5551 km). We also estimated the massloss rate from KOI 2700b, and we obtained Ṁ values from the interval [5.05 × 10^{7}, 4.41 × 10^{15}] g s^{−1}. On the other hand, we could not draw any satisfactory conclusions about the typical grain size in the dust tail.
Key words: planets and satellites: general / planetstar interactions / scattering
© ESO 2018
1 Introduction
Strong irradiation in closein exoplanet systems may cause mass loss from the planet (Burrows & Lunine 1995; Guillot et al. 1996), such as that detected, for example, in HD 209458b (VidalMadjar et al. 2003, 2004) and in HD 189733b (Lecavelier Des Etangs et al. 2010; Bourrier et al. 2013). Several theoretical studies have been devoted to this subject (Yelle 2004; Tian et al. 2005; Hubbard et al. 2007). In certain cases the mass loss from the exoplanet may cause formation of a cometlike tail. The hypothesis that a closein exoplanet may have a cometlike tail was suggested by Schneider et al. (1998) and was revisited by Mura et al. (2011). The transit light curve of dusty extrasolar comets was investigated, for example, by Lamers et al. (1997) and Lecavelier Des Etangs et al. (1999). The first good evidence for exocomet transits was presented by Rappaport et al. (2018).
The first exoplanet with a cometlike tail, KIC 12557548b, was discovered from Kepler longcadence data by Rappaport et al. (2012), and was found to be a closein exoplanet with an extremely short orbital period of P_{orb} ≃ 0.65356 days. Rappaport et al. (2012) suggested that the planet’s size is not larger than Mercury. The light curve of this planet was studied in more detail by Brogi et al. (2012), Budaj (2013), Kawahara et al. (2013), Croll et al. (2014), van Werkhoven et al. (2014), Bochinski et al. (2015), Schlawin et al. (2016) and van Lieshout et al. (2016). Brogi et al. (2012) and Budaj (2013) first validated the disintegratingplanet scenario using a model and both found that dust particles in the tail have typical radii of about 0.1–1 micron. PerezBecker & Chiang (2013) proposed a model of the atmospheric escape via the thermal wind that is only effective for planets which are less massive than Mercury. Gravity of the more massive planets would provide too deep a potential barrier for the wind. Garai et al. (2014) searched for cometlike tails in a sample of 20 closein exoplanet candidates with a period similar to KIC 12557548b from the Kepler mission, however, none of the exoplanet candidates showed signs of a cometlike tail. This result is in agreement with the model proposed by PerezBecker & Chiang (2013).
Recently, two more exoplanets have been discovered, KOI 2700b and K222b, whose transit shapes show evidence of a cometlike tail (Rappaport et al. 2014; SanchisOjeda et al. 2015). The Kepler object KOI 2700b (KIC 8639908b) was discovered by Rappaport et al. (2014) as the second exoplanet with a cometlike tail. It exhibits a distinctly asymmetric transit profile, likely indicative of the emission of dusty effluents and reminiscent of KIC 12557548b. The orbital period of KOI 2700b was determined by the discoverers using LombScargle analysis (P_{orb} = 0.910023(4) days) and χ^{2} minimization (P_{orb} = 0.910022(5) days). The host star has T_{eff} = 4435 K, M ≃ 0.63 M_{⊙} and R ≃ 0.57 R_{⊙}, comparable to the parameters ascribed to KIC 12557548. The upper limit of the planet radius R_{p} ~ 1.06 R_{⊕}, determined by the same authors, is in agreement with the theory of the thermal wind and planet evaporation (PerezBecker & Chiang 2013). The massloss rate, found by Rappaport et al. (2014), is ~6 × 10^{9} g s^{−1} or ~2 lunar masses per Gyr. The composition of the dust ejected by KIC 12557548b and KOI 2700b was proposed by van Lieshout et al. (2014). The observed tail lengths are consistent with dust grains composed of corundum (αAl_{2}O_{3}) or ironrich silicate minerals. Pure iron and carbonaceous compositions are not favored.
In this paper we aim to verify the disintegratingplanet scenario of KOI 2700b by modeling its light curve and put constraints on various tail and planet properties, as was done in the case of KIC 12557548b. The light curve of the object was analyzed similarly as in Budaj (2013), however, in comparison with Budaj (2013), in this work the lightcurve modeling process is more sophisticated. We first improved the orbital period and constructed the phasefolded and binned transit light curve of KOI 2700b using longcadence Kepler observationsfrom the quarters 1–17 (Sect. 2). Subsequently, we searched for longterm orbital period variations (Sect. 3). We then calculated phase functions and used the radiativetransfer code SHELLSPEC (Budaj & Richards 2004, 2010). We applied this code iteratively to model the observed light curve of KOI 2700b. Mie absorption and scattering on spherical dust grains with realistic dust opacities, phase functions, and finite radius of the source of the scattered light were taken into account (Sect. 4). The resulting parameters and model light curves are described and discussed in Sect. 5. Our findings are concluded in Sect. 6.
2 The light curve of KOI 2700b
2.1 Observations
We used the publicly available Kepler data from the quarters 1–17 in the form of Simple Aperture Photometry (SAP) fluxes. Presearch Data Conditioning Simple Aperture Photometry (PDCSAP) fluxes often deliver overcorrected light curves, in which the astrophysical signal is reduced or even canceled, hence we avoided using these data. A comparison between the SAP and PDCSAP Kepler phasefolded and binned transit light curves of KOI 2700b is depicted in Fig 1. Only the long cadence data were used to construct the light curve. These are 64 842 observations with an exposure time of about 30 min. Kepler observations were reduced in a similar manner to in Budaj (2013) and Garai et al. (2014).
Each quarter has a different flux level. Consequently, fluxes within each quarter were normalized to unity. We then improved the orbital period of the exoplanet found by Rappaport et al. (2014). For this purpose we used the method of phase dispersion minimization described in Sect. 3. First we phased the data with the orbital period of P_{orb } = 0.910023 days (Rappaport et al. 2014). We used the phase 0.5 for transits. Subsequently, data were cut into segments each covering oneorbital period. Each segment of data was fitted with a linear function. During the fitting procedure the part of the data between phases 0.45 and 0.65, covering the transit, was excluded from the fit. Consequently, the linear trend was removed from each chunk of data (including the transit data). This method can effectively remove the long term variability (mainly variability of the host star due to spots and rotation) while it does not introduce any nonlinear trend to the phased light curve. The final value of the orbital period was then found in these detrended data. We found the orbital period of P_{orb } = 0.9100259(15) days via thephase dispersion minimization method, which is in good agreement with the period presented by the discoverers. Finally, the data were phased with this new orbital period. Subsequently, we checked our result using the Fourier method (Deeming 1975) and found P_{orb } = 0.910042(7) days. This method confirmed the orbital period resulting from the phase dispersion minimization method. We estimated the error by performing Monte Carlo simulations generating and analyzing 100 artificial datasets.
To reduce the noise, the phased light curve was subjected to an averaging. We first applied the running window averaging technique. We used a window with the width in units of phase of 0.00075 (1min width, negligible in comparison with the exposure time), 0.005 (6.5min width, according to Rappaport et al. (2014)), 0.011 (15min width, about 1/2 of the exposure time), and 0.023 (30min width, comparable with the exposure time). The step was set as 1/10 window width in every case. Since the obtained averaged light curves showed systematic fluctuations at short timescales (about 20–30 min), easily visible, for example, in the light curve smoothed using the 15minwidth window (Fig. 2), we performed period analysis with the main aim toexplore its nature, using the method of phase dispersion minimization (see Sect. 3). For this purpose we used the available short cadence SAP Kepler detrended data, and searched for periodicities from 1 min to 1 h. Figure 3 shows that there are no significant frequencies in this interval, indicating that the observed systematic fluctuations of averaged light curves may be due to poor signal to noise ratio or to the exposure time. Subsequently, we applied averaging via data binning. Instead of stable window width, this technique uses a certain number of data points per bin. We examined the cases from 500 to 1500 data points per bin with a step of 50 data points per bin, and finally we selected the case of 1000 data points per bin, according to the minimum scatter and relatively good coverage of the phased light curve with averaged data points. The obtained light curve is depicted in Fig. 4, as well as in Fig. 1 (SAP flux). No systematic fluctuations are visible in this case and, if not mentioned otherwise, we used this averaged light curve in our analysis.
Fig. 1 Comparison between the SAP and PDCSAP Kepler phasefoldedand binned transit light curves of KOI 2700b. The illustrative vertical error bars on the left side represent the median uncertainties of data points depicted with the same color. 
Fig. 2 Averaged light curve of KOI 2700b, smoothed using the running window technique. The window width and step was 0.011 and 0.0011, respectively (in units of phase). The light curve shows systematic fluctuations at short timescales (20–30 min). See the text. 
Fig. 3 Results of the period analysis of KOI 2700b. We searched for periodicities from 1 min to 1 h using the method of phase dispersion minimization (see Sect. 3), but we did not find any significant frequency inthis interval. See the text. 
Fig. 4 Phasefolded and binned transit light curve of KOI 2700b. Red points are observations and green points represent the averaged data. We note that averaged data points have uncertainties that are too small to be discerned in this plot. Error bars of these data are well depicted, for example, in Fig 9. 
2.2 Comparison with the light curve of KIC 12557548b
The light curve of the exoplanet KOI 2700b is peculiar and very interesting. It exhibits a distinctly asymmetric transit profile with a sharp ingress followed by a short sharp egressand a long smooth egress; reminiscent of the light curve of the exoplanet KIC 12557548b. The observed light curve of the exoplanet KOI 2700b might show a relatively weak pretransit brightening (see for example, Fig. 1), however, we could not statistically identify this lightcurve feature, previously found at KIC 12557548b as significant. This latter exoplanet also exhibits a relatively weak posttransit brightening (Fig. 5). These brightenings are caused by the forward scattering on dust particles in the tail (Brogi et al. 2012; Budaj 2013). The light curves have one more similarity. The transit ingress begins at an orbital phase of about 0.45 and the transit egress ends at an orbital phase of about 0.65 (Fig. 5). The transit duration in phase units is 0.2 in both cases. The transit depth, however, is very different. KIC 12557548b has an average transit depth of about 0.005 in flux (i.e., 0.5%), while KOI 2700b has about ten times shallower average transit depth, 0.0004 in flux (i.e., 0.04%).
Fig. 5 Comparison between the SAP Kepler phasefolded and binned transit light curve of exoplanets KIC 12557548b and KOI 2700b. The illustrative vertical error bars on the left side represent the median uncertainties of data points depicted with the same color. 
2.3 Evolution of the light curve
The relatively shallow average transit depth of the exoplanet KOI 2700b (~ 0.04%) makes it difficult to study possible transittotransit variations in the transit depth, as was done in the case of KIC 12557548b. On the other hand, Rappaport et al. (2014) reported that the transit depth appears to decrease systematically with time during the Kepler mission. During this step we first dividedour detrended data into five approximately equal segments (in time) and each segment of data was subsequently phased with the orbital period of P_{orb } = 0.9100259 days (Sect. 2.1), found via the method of phase dispersion minimization (see Sect. 3). Since we used less data points per light curve (approximately 13 000 vs. 64 000), we applied 200 data points per bin (instead of 1000 data points per bin) during the data averaging. The resulting phased lightcurves are depicted in Fig. 6. Based on our analysis we can confirm the result concerning decreasing transit depth, found by Rappaport et al. (2014). On the other hand, we can also conclude that this depth decreasing is not monotonical. The last three light curves have approximately the same transit depth, which may indicate the possibility of a longterm periodic variability as well.
Fig. 6 Transit depth decreasing at KOI 2700b. The phasefolded and binned transit light curves of the exoplanet were constructed from 3 and 4 quarters (Q) of detrended data. The illustrative vertical error bar on the left side represents the median uncertainty of all data points. 
3 Search for longterm orbital period variations
PerezBecker & Chiang (2013) suggest that disintegrating planets undergo negligible orbital period variations due to the evaporation process. During this step we examined this prediction via searching for possible longterm changes of the orbital period of KOI 2700b. For this purpose we used the method of phase dispersion minimization (Stellingwerf 1978), application PDM2, version 4.13^{1}, which uses bin structure 50/2. The detrended data were used for this purpose (see Sect. 2.1). We assumed that the period changes linearly: (1) where β is a dimensionless value, but is often expressed in days/million years (d/Myr). The output from the analysis is a curve (Fig. 7), which shows the dependence of Θ _{min } as a function of β. The term Θ_{min} is a dimensionless statistical parameter (Stellingwerf 1978). The minimum value of Θ_{min} indicates the speed of the period change – β. We fitted a parabola to this curve and obtained β = 1.3 ± 3.2 d/Myr, which means that there is no significant evidence for the longterm orbital period change during the time span of the Kepler observations. This result is in agreement with the prediction, presented by PerezBecker & Chiang (2013). The error was estimated by means of Monte Carlo simulations. We generated and analyzed 100 artificial datasets with the same standard deviation as the original data.
Fig. 7 Search for a longterm orbital period change of KOI 2700b. The minimum of the parameter Θ _{min } corresponds to the periodchange speed β = 1.3 ± 3.2 d/Myr, which means that there is no significant evidence for the longterm orbital period change. 
4 Lightcurve modeling with the code SHELLSPEC
4.1 Calculation of the optical properties of the dust
According to van Lieshout et al. (2014), the cometlike tail of KOI 2700b may consist of dust grains composed of corundum (αAl_{2}O_{3}) or ironrich silicate minerals. Therefore during the modeling we applied selected species, which have similar chemical composition: γalumina, enstatite, forsterite, olivine (with 50% magnesium and 50% iron) and pyroxene (with 40% magnesium and 60% iron). Aluminum oxide (Al_{2 } O_{3}) occurs naturally in its crystalline polymorphic phase as the mineral corundum, which is αalumina. Aluminum oxide also exists in many other phases; each has a unique crystal structure and properties. In our model we applied cubic γAl_{2}O_{3}. Enstatite is the magnesiumrich endmember of the pyroxene silicate mineral series enstatite (MgSiO_{3 } ) – ferrosilite (FeSiO_{3 } ). Forsterite (Mg_{2 } SiO_{4}) is the magnesiumrich endmember of the olivine solid solution series. It is isomorphous with the ironrich endmember, fayalite. Besides forsterite, in our model we also applied olivine with 50% magnesium and 50% iron. The pyroxenes are a group of minerals with the general formula XY(Si,Al)_{2}O_{6}. The group has 20 endmembers and one of the endmembers is enstatite (see above). Pyroxenes are, however, usually isomorphous mixtures of some endmembers. Therefore, in our model we also applied a pyroxene with 40% magnesium and 60% iron.
Calculating the optical properties (i.e., opacities and phase functions, see below) of the dust grains composed of a mineral is relatively difficult and requires a lot of computing time. That is why Budaj et al. (2015) prepared tables of phase functions, opacities, albedos, equilibrium temperatures and radiative accelerations of dust grains in exoplanet systems. The tables cover the wavelength range of 0.2 to 500 microns and 21 particle radii from 0.01 to 100 microns for several species. Their assumptions include spherical grain shape, Deirmendjian particle size distribution (Deirmendjian 1964) and Mie theory. Budaj et al. (2015) used a widely available Mie scattering code CALLBHMIE^{2}, which calls iteratively the Mie scattering subroutine BHMIE (Bohren & Huffman 1983). The tables are freely available on the web^{3} and are applicable for our purposes in the code SHELLSPEC. From these tables we selected opacities and phase functions for γalumina, enstatite, forsterite, olivine (with 50% magnesium and 50% iron) and pyroxene (with 40% magnesium and 60% iron). In our model we applied the following particle radii: 0.01, 0.1 and 1.0 micron.
Dust can absorb the impinging radiation and convert it directly into heating of the grains. This process is called absorption or true absorption and it is quantified by the absorption opacity. Dust can also scatter radiation in a process called scattering without being heated. This process is quantified by the scattering opacity. The sum of the absorption opacity and the scattering opacity is the total opacity (or simply the opacity). Furthermore, scattering can be highly asymmetric, a property that is described by means of the phase function, which depends on the scattering angle (the deflection angle from the original direction of the impinging radiation). The most prominent feature is a strong forward scattering, when the scattering angle is nearly zero. Brightenings on the light curve of KIC 12557548b are also caused by the forward scattering on dust particles in the tail (Brogi et al. 2012; Budaj 2013). We also have to take into account the fact that from the viewpoint of dust particles the parent star has a nonnegligible angular dimension on the sky. Dust particles of KOI 2700b in the distance of a = 3.363 R_{⊙} from the parent star view this star as a disk with an angular diameter of about 19.4 deg. The true dimension of this star is R_{s} = 0.57 R_{⊙} (Rappaport et al. 2014). To take this effect into account we have to split the stellar disk into elementary surfaces and integrate the phase function over the disk. We calculated phase functions with a very fine step in the interval [0, 10] deg, because of the strong forward scattering and, consequently, the disk averaged phase function with a very fine step near the edge of the stellar disk. For this purpose we applied the software DISKAVER^{4}. It assumes a quadratic limb darkening of the stellar surface: (2) where I_{ν}(0) is intensity perpendicular to the surface of the source and Θ is angle between the line of sight and a normal to the surface. The quadratic limb darkening coefficients u_{1} and u_{2} were linearly interpolated based on the stellar parameters T_{eff} = 4433 K, logg = 4.721 (cgs) and Fe∕H = −0.2 (Rappaport et al. 2014) for the Kepler passband, using the online applet EXOFAST – Quadratic Limb Darkening^{5}, which is based on the IDLroutine QUADLD (Eastman et al. 2013), as u_{1} = 0.622 and u_{2} = 0.113. This software interpolates the Claret & Bloemen (2011) quadratic limb darkening tables. Calculations with the software DISKAVER were performed at a wavelength of 0.6 microns for consistency with the Kepler passband. An illustration of such phase functions, which take into account the finite dimension of the source of light, is depicted in the Fig. 8. Corresponding phase functions without taking into account the finite dimension of the source of light are also shown. These exhibit a strong peak near phase angle zero, which is the socalled forwardscattering peak. Larger particles and/or shorter wavelengths tend to have a stronger and narrower forwardscattering peak than smaller particles and/or longer wavelengths. Therefore, taking into account the finite dimension of the source of light is important mainly for larger particles and/or shorter wavelengths.
Fig. 8 Phase functions at 0.6 micron for different dust particle radii of γalumina. Lines without symbols are phase functions assuming a point source of light. Lines with symbols are phase functions that take into account the finite dimension of the stellar disk. The vertical line illustrates the angular dimension of the stellar disk as seen from KOI 2700b. 
4.2 Model and lightcurve modeling
The phased and averaged transit light curve of the exoplanet KOI 2700b (Fig. 4) was in principle modeled using the radiativetransfer code SHELLSPEC^{6} (Budaj & Richards 2004, 2010), version 39. This code calculates the light curves and spectra of interacting binaries or exoplanets immersed in the threedimensional (3D) circumstellar or circumplanetary environment. It solves simple radiative transfer along the line of sight and the scattered light is taken into account under the assumption that the medium is optically thin. A number of optional objects (such as a spot, disk, stream, ring, jet, shell) can be defined within the model, or it is possible to load a precalculated model from an extra file. Synthetic light curves or trailing spectrograms can be produced by changing our viewpoints on the 3D object. The software is written in Fortran77, however, some previous versions were transformed into the language Fortran90; see for example, Tkachenko et al. (2008); Šejnová et al. (2012).
For our purpose we used an optional object in the form of a ring, which we subsequently modified. From the transit shape we can conclude that transits are caused primarily by the cometlike tail, and not by the solid body of the planet. We can model such a cometlike tail as part of a ring with a nonnegligible thickness around a central star. Therefore during our calculations we assumed a spherical and limbdarkened central star (using the coefficients u_{1 } = 0.622 and u_{2 } = 0.113 from Eq. (2)), with a radius of R_{s } = 0.57 R_{⊙}, mass of M_{s } = 0.632 M_{⊙} and effective temperature of T_{eff} = 4433 K (Rappaport et al. 2014), located in the geometrical center of our ring. We modeled the cometlike tail as part of a ring with a radius of a = 3.363 R_{⊙}. Its geometrical crosssection is monotonically enlarging from the planet to the end of the ring, which is located at 60 deg behind the planet. At this point the ring is truncated. In our calculations we applied the same length of the tail as per Budaj (2013), because the transit duration in phase units is very similar in both cases (Fig. 5). The crosssection of the ring C and dust density along the ring ρ are allowed to change with the angle t [rad]: (3) or (4)
where ρ(0), C(0), and t(0) are the dust density, crosssection, and phase angle of view at the beginning of the ring and A1, A2 are the density exponents to model the dust destruction in the tail. We can see that there is a strong degeneracy between the crosssection and the dust density at a certain phase angle. In general, if we increase the crosssection and we want to obtain an appropriate model of the observed light curve, the result is that we need to decrease the dust density. There is a degeneracy relation Cρ = const. Since the dimension of the dust tail of KOI 2700b is unknown, we arbitrarily defined a geometrical crosssection of the tail, which was preferable for our computation process (in terms of the grid density, grid dimension, and computing time). Therefore, in our calculations we assumed a dust tail with a crosssection of 0.05 × 0.05 R_{⊙} at the beginning and 0.09 × 0.09 R_{⊙} at its end. These values are also in agreement with the escape velocity from a Mercurysized small planet (a few km s^{−1}). We note that in certain cases the chosen crosssection of the tail can affect the results about the grain size; for example, at higher inclinations (i < 78 deg) large amounts of dust will not transit the parent star and this can affect the identification of the typical particle radii of the tail (described in Sect. 5.2). However, as we described in Sect. 5.1, we obtained i in the range [85.1, 88.6] deg (the average value is i = 87.06 deg), that is, the latter situation is not the case here. The ring was located in the orbital plane of the planet.
On the other hand, based on the transit depth, we can assume that there is still ample room for a contribution to the light curve from the solid body of the planet. Therefore, we also included the solid body of the planet in the model, defined as a dark, nontransparent sphere and located at the beginning of the ring. The dimension of the planet was parametrized as R_{p}∕R_{s} (ratio of the radii). Although Rappaport et al. (2014) present the formal upper limit of the planet radius of R_{p} ≤ 1.06 R_{⊕} (2σ limit), which roughly corresponds to R_{p}∕R_{s} ≤ 0.017, based on the theory of the thermal wind and planet evaporation (PerezBecker & Chiang 2013) we can expect that the planet itself must be Mercurysized, or smaller than Mercury, which corresponds to R_{p}∕R_{s} ≤ 0.006.
In the SHELLSPEC code the central star with the defined object is located in a 3D grid. The code enables the user to look on the grid from different points of view and to calculate the corresponding flux. The flux is always calculated in the observer’s line of sight. The orbital inclination angle i corresponds the inclination of the intrinsic rotation axis of the model to the line of sight. At each point of view we calculated the final flux as (s + r)∕s = f, where s means modeling the flux from the parent star, s + r means modeling the parent star with the ring and planet, and f is the final and normalized flux from the system. In this way we also eliminated fluctuations due to the grid structure of objects in our model. To speed up the computation process, we reduced the number of points by selecting the neartransit part of the light curve. We calculated 33 artificial data points per synthetic light curve, between phases 0.25 and 0.75, in the same phases as per the observed light curve. This cropping of the synthetic light curve is justified, since we are interested only in the neartransit region. The synthetic light curves were subsequently convolved with a boxcar with a width of 30 min, simulating the integration time of the Kepler long cadence exposure. Convolved light curves were used for comparison with the observed light curve.
For the modeling process we generated and used an iterative procedure, which applied the SHELLSPEC code as a subroutine, and searched for the best fit parameters. Four free parameters were adjusted simultaneously during the fitting procedure – the orbital inclination angle i [deg], the dust density at the beginning of the ring ρ(0) [g cm^{−3 } ], the ratio of the radii R_{p}∕R_{s}, and the density exponent A2. Since model light curves calculated using A1 and A2 differed only slightly, we decided to use the density exponent A2 during the modeling, because it explains the dust destruction in the tail better than A1. One more free parameter – the transit midpoint phase shift of the synthetic light curve with respect to the observed light curve (Δ φ_{0 } ) – was adjusted only before the modeling process and then was kept fixed to its best value. This parameter reflects the unknown midtransit time of the planet. Every synthetic light curve was shifted in phase by Δ φ_{0 } = −0.235. The advantage of this treatment is that it saves computing time; on the other hand we cannot exclude the possibility that uncertainties on the resulting parameters are underestimated. We prepared a parameter space with a certain stepping, composed from the abovementioned four free parameters and limited based on the previous stepwise test calculations as follows. The parameter i: [90, 80] deg; the parameter ρ(0): [0.5000, 0.0001] × 10^{−15} g cm^{−3}; the parameter R_{p}∕R_{s}: [0.017, 0.001]; and the parameter A2: [−1.0, −25.0]. During the fitting procedure all combinations of free parameters were examined and the observed light curve was compared with the corresponding model light curve. A formal, quantitative goodnessoffit was measured via determination of reduced χ^{2} (). The best fit corresponds to the minimum value of . During the next iteration we reduced the range of the parameters and at the same time we used a finer stepping in the parameter space to fit the observed transit light curve better. As the first iteration we used the abovementioned range of the parameters and the stepping as follows. The parameter i: 2.0 deg; the parameter ρ(0): 0.0100 × 10^{−15} g cm^{−3}; the parameter R_{p}∕R_{s}: 0.004; and the parameter A2: 4.0. We selected the best value of a given parameter, found in the previous iteration, as a median () value of the new parameter range. As the second iteration we used the following range/stepping of the parameters: ± 2.0, stepping 1.0 deg; ± 0.0100 × 10^{−15}, stepping 0.0010 × 10^{−15} g cm^{−3}; ± 0.004, stepping 0.002; and ± 4.0, stepping 1.0 and during the last (third) iteration, similarly: ± 1.0, stepping 1.0 deg; , stepping 0.0001 × 10^{−15} g cm^{−3}; ± 0.002, stepping 0.001; and ± 1.0, stepping 0.5. To obtain final results we applied these three iterations. In this way we reduced the computing time.
To derive uncertainties on the resulting parameters we applied Monte Carlo simulations. First, we generated an artificial dataset with the same standard deviation as the observed data. Then, the artificial data were fitted with a model light curve and the best fit parameters were determined. The range of the parameters and the stepping corresponded to the range and stepping used in the last (i.e., in the 3rd) iteration (see the end of the previous paragraph). We note that the full procedure is not needed in this case, since the artificial data differ only slightly from the observed data. This process was repeated 100 times and, subsequently, to estimate the uncertainty in the given parameter the standard deviation was calculated.
We executed 15 joint fitting procedures – one for each combinationof species and dust particle size. When describing our main results in Sect. 5 we use the average value (without uncertainties) and the obtained range of the given parameter (including uncertainties). For individual values of parameters we refer the reader to Table 1.
Overview of the resulting parameter values.
5 Results and discussion
5.1 Inclination of the tail
The orbital inclination angle i is a very important system parameter. We first assumed that the orbital plane of the exoplanet and its tail have an inclination of i = 90 deg with respectto the plane of the sky. Our expectations, however, were not confirmed. During the modeling process it turned out that the transit of KOI 2700b is not edgeon, but based on our models, the orbital inclination angle is very probably close to the value of i = 87 deg. We obtained i in the range [85.1, 88.6] deg^{7} (see Table 1 for individual values). We obtained the value of i ≈ 88 deg in the cases of 0.1micron γalumina grains, 1micron forsterite grains and 0.01micron olivine grains. Similarly, we obtained a value of i ≈ 86 deg in the cases of 0.01micron γalumina grains and 1micron pyroxene grains. In other cases we obtained for this parameter the value of i ≈ 87 deg. The resulting parameter values are, however, in agreement within a 1σ limit (see Table 1).
5.2 Properties of the dust
A decrease in dust density was detected in the cometlike tail of KIC 12557548b by Budaj (2013). This feature of the tail is described by the density exponent A1 or A2 (see Eqs. (3) and (4)). As mentioned in Sect. 4.2, for our purposes we used the density exponent A2. It is one of the most important free parameters in our model. We first assumed that the dust density decreases along the tail of KOI 2700b in a similar manner to KIC 12557548b, with the exponent − 20 or − 25. Based on our fitting procedure, however, we obtained A2 in the range [−6.5, −8.2] (see Table 1 for individual values). The average value of the parameter is A2 = −7.13. In the case of 1micron γalumina grains the resulting parameter value was A2 ≈−8.0, and in the cases of 1micron enstatite grains and 1micron forsterite grains we obtained A2 ≈−7.5. For other dust particle sizes and species we obtained A2 ≈−7.0 as a resulting parameter value. We note that the range of values found for the parameter A2 is in good agreement with the exponential decay factor value (S), found by Rappaport et al. (2014). Ignoring the C(t) dependence of the dust density in the tail, the density exponent A2 can be converted to the exponential decay factor as S = −A2∕π; hence the range of A2 values ([−6.5, −8.2]; see Table 1 for individual values) corresponds to S of about [2.1, 2.6]. The exponential decay factor value, found by Rappaport et al. (2014), is S = 2.4 ± 1.0.
The density exponent A2 needs a certain density at the beginning of the ring ρ(0) to fit the transit depth properly at each chemical composition and dust particle size. In Table 1 we present the resulting values for this free parameter and also the resulting values for the derived parameter τ_{max}, which is the highest optical depth through the tail in front of the star, and which follows from ρ(0). These values were derived based on a tail model with a crosssection of 0.05× 0.05 R_{⊙} at the beginning and 0.09 × 0.09 R_{⊙} at its end. Table 1 shows that 0.01micron and 1micron grains have very similar individual values of ρ(0). On the other hand, the resulting values of τ_{max} are different in the case of these two groups. The difference is about 3.00× 10^{−3}. The third group, with 0.1micron grains, has individual values of ρ(0) that are about one order of magnitude smaller, but similar resulting values of τ_{max} to 0.01micron grains. The average values of the parameters are ρ(0)= 0.1667 × 10^{−15} g cm^{−3 } and τ_{max} = 5.24 × 10^{−3}, and we obtained these parameters in the range [0.0151, 0.3333] × 10^{−15} g cm^{−3 } and [4.28, 7.43] × 10^{−3}, respectively.
Based on modeling of the transit we cannot determine the chemical composition of the dust ejected by KOI 2700b. That is why we applied selected species, according to van Lieshout et al. (2014). Theoretically, the situation in the case of particle sizes is slightly better. We can estimate the typical particle size in the dust tail based on part of the light curve of KOI 2700b, where a pretransit brightening is expected or possible, as was done at KIC 12557548b by Budaj (2013). This brightening is caused by forward scattering on dust particles in the tail (Brogi et al. 2012; Budaj 2013) and the forward scattering is sensitive to the dust particle size (Fig. 8). Therefore, we can compare our model light curves with observations, especially in phases between 0.40 and 0.45, where the pretransit brightening was detected at KIC 12557548b (see Fig. 5), and search for the particle size, which satisfies the observed light curve of KOI 2700b better. Figures 9–11 show that this part of the model light curve is sensitive to the dust particle size, however, differences between the model pretransit brightenings are not enough to draw any satisfactory conclusions on the typical grain size in the dust tail. Since every model pretransit brightening satisfies the observations, we cannot select a model with appropriate grain size.
Fig. 9 Model light curves calculated for 0.01micron, 0.1micron, and 1micron grains of γalumina (left) and enstatite (right), compared with the observed light curve of KOI 2700b. 
Fig. 10 Model light curves calculated for 0.01micron, 0.1micron, and 1micron grains of forsterite (left) and olivine (right), compared with the observed light curve of KOI 2700b. 
5.3 The dimension of the planet solid body
The dimension of the solid body of KOI 2700b was a key point in our model. It was parametrized using the following free parameter: the ratio of the radii (R_{p}∕R_{s}). Based on the results of Rappaport et al. (2014), we selected the value of R_{p}∕R_{s} = 0.017 as an upper limit for this free parameter. However, as these authors pointed out, the corresponding value, R_{p} = 1.06 R_{⊕}, needs to be improved. Therefore, we strongly expected that the planet radius would be smaller, R_{p}∕R_{s} ~ 0.006, which corresponds to a Mercurysized planet. The parameter value R_{p}∕R_{s} ≤ 0.006 is also predicted by the theory of the thermal wind and planet evaporation (PerezBecker & Chiang 2013). The resulting parameter values confirmed our expectations. Table 1 shows that in the cases of 1micron γalumina, enstatite, and forsterite grains, we obtained R_{p}∕R_{s} ≈ 0.007, in the cases of 0.01micron enstatite and forsterite grains we obtained R_{p}∕R_{s} ≈ 0.008, and in other cases the resulting parameter value was 0.009. The average value of the parameter is R_{p}∕R_{s} = 0.008 (R_{p} = 0.497 R_{⊕}, or 3172 km). We obtained R_{p}∕R_{s} in the range [0.0, 0.014] (see Table 1 for individual values), hence we can consider the value of R_{p}∕R_{s} = 0.014, which corresponds to R_{p} = 0.871 R_{⊕} (or 5551 km), as an upper limit on the radius of the planet solid body of KOI 2700b. The upper limit of the planet radius is in agreement with the theory of the thermal wind and planet evaporation and is also derived with higher accuracy than the result found by Rappaport et al. (2014). A comparison of the observed light curve with the solidbody transit light curve of KOI 2700b is depicted in Fig. 12. For a comparison, we also plotted the transit light curve of an Earthsize and a Mercurysize planet. The light curves were modeled using the software JKTEBOP^{8} (Southworth et al. 2004) and its task No. 2. We applied i = 87 deg (according to the average of listed values in Table 1) and in the case of KOI 2700b the value of R_{p}∕R_{s} = 0.008 (the average value of the parameter) and R_{p}∕R_{s} = 0.014 (the upper limit of the planet radius).
Fig. 12 Phasefolded and binned transit light curve of KOI 2700b (observations) compared with the transit light curve of the solidbody of KOI 2700b (according to the average value and upper limit of the planet radius; see the text). Furthermore, the transit light curve of an Earthsize planet and a Mercurysize planet is also plotted. The orbital inclination angle is i = 87 deg. 
5.4 Massloss rate from the planet
We can estimate the massloss rate Ṁ from the planet assuming that the tail with the mass M_{tail} is replenished within a period, which is equal to the evaporation time scale of grains in the tail T_{evap}: (5)
The mass of the tail can be calculated by integration of Eq. (4) through the volume of the tail V as follows: (6)
Assuming that the geometrical crosssection of the tail is constant from the planet to the end of the ring, that is, C(t) = C(0), and substituting t − t(0) with L∕a, where L is the length of the tail and a is the radius of the ring in cm, we can simplify Eq. (6) to: (7) where the volume element of the tail can be expressed as dV = C(0)dL. After integration of Eq. (7) we get: (8)
The evaporation time scale of grains in the tail can be calculated using the equations and tables presented by van Lieshout et al. (2014). First, we can approximate the equilibrium vapor pressure as [dyn cm^{−2}], where T_{eqv} is the equilibrium temperature of the dust grain and 𝒜, ℬ are materialdependent sublimation parameters. The values of T_{eqv} can be interpolated from the tables of Budaj et al. (2015), and the values of 𝒜, ℬ are tabulated for some species in van Lieshout et al. (2014). Subsequently, we can approximate the massloss flux from the surface of the dust grain J [g cm^{−2} s^{−1}] as: (9) where α is the evaporation coefficient, μ is the molecular weight of the molecules that sublimate, m_{u} is the atomic mass unit and k_{B} is the Boltzmann constant in cgs units. The values of α and μ are also tabulated for some species in van Lieshout et al. (2014). Finally, we can calculate the evaporation time scale of grains in the tail as: (10)where r is the grain radius and ρ_{d} is the bulk density of dust species. The values of ρ_{d } are tabulated for some species in van Lieshout et al. (2014) or in Budaj et al. (2015). Using Eqs. (8) and (10) we can rewrite Eq. (5) to the following final form: (11)
For our purposes we used C(t) = C(0) = 0.05 × 0.05 R_{⊙}, the parameter values presented in Table 1 of this work, in Tables 2 and 3 of van Lieshout et al. (2014) and in the Tables of Budaj et al. (2015). We made three groups of estimations of the massloss rate from the planet: one for alumina, one for enstatite, and one for forsterite. We note that in the cases of enstatite and forsterite we interpolated the equilibrium temperature of pyroxene with 40% magnesium and 60% iron, and olivine with 50% magnesium and 50% iron, respectively. Pure enstatite or forsterite have too low a value of T_{eqv}, out of the temperature range for which 𝒜 and ℬ were determined. Our individual results are presented in Table 2. This procedure yields an average massloss rate from KOI 2700b of Ṁ = 3.03 × 10^{14} g s^{−1 } , and we obtained Ṁ in the range [5.05 × 10^{7}, 4.41 × 10^{15}] g s^{−1 } . This range of values is consistent with the massloss rate estimate of 6 × 10^{9}g s^{−1}, presented by Rappaport et al. (2014).
Estimation of the massloss rate from KOI 2700b.
5.5 Comparison with KIC 12557548b
Since the first exoplanet with a cometlike tail, KIC 12557548b, was analyzed similarly by Budaj (2013), it is very interesting to compare his results with the results obtained in this work. Unfortunately, we cannot compare the typical grain size of the tail, which we could not derive in the case of KOI 2700b, and the dimension of the planet solid body, which was not included in the model of KIC 12557548b, used by Budaj (2013). On the other hand, an upper limit on the radius of the planet solid body of KIC 12557548b was presented by Brogi et al. (2012) and by van Werkhoven et al. (2014), which we can use for our comparison. Similarly, the massloss rate from KIC 12557548b was calculated, for example, by Rappaport et al. (2012).
KOI 2700b and KIC 12557548b have very similar lightcurve shape with some peculiarities. Mainly, the light curve displays a sharp ingress, followed by a short sharp egress and a long smooth egress. The light curve of KIC 12557548b shows a relatively significant pretransit brightening as well, which in thelight curve of KOI 2700b is not significant or even present at all. A relatively weak posttransit brightening is visible only in the light curve of KIC 12557548b. The transit duration in phase units is also very similar. The transits of the planets are very probably not edgeon. It seems that KOI 2700b has an orbital inclination angle of i ≈ 87 deg, while KIC 12557548b has i ≈ 82 deg. On the other hand, the transit depth of the planets is very different, KOI 2700b has an average transit depth about ten times shallower. Due to the similar light curves, the modified ring model of the tail fits the observations in both cases well.
The tails are 60 deg long in both models, that is, in this work and in Budaj (2013). The cutoff at this distance from the planet is an assumption of the model. The dust density in the tail decreases in both cases with the distance from the planet; however, these density decreasings are not similar. The dust density in the tail of KIC 12557548b decreases rapidly along the tail, which is well defined by the density exponent A2 ≈−20 or A2 ≈−25. This may produce the week, but still observable posttransit brightening at KIC 12557548b. In comparison with this planet, KOI 2700b exhibits a relatively low dust density decreasing in its tail, which is defined by the density exponent A2 ≈−7. Therefore, in this case we cannot expect any posttransit brightening.
KOI 2700b and KIC 12557548b may have a planet solid body radius less than ~ 1 R_{⊕}. Since we obtained R_{p}∕R_{s} in the range [0.0, 0.014] (see Table 1 for individual values), we can consider the value of R_{p}∕R_{s} = 0.014, which corresponds to R_{p} = 0.871 R_{⊕}, or 5551 km, as an upper limit of the planet radius of KOI 2700b. In comparison with KOI 2700b, KIC 12557548b may have a slightly bigger radius as well, up to ~ 1 R_{⊕}.Brogi et al. (2012) present the value of R_{p} = 1.15 R_{⊕} and van Werkhoven et al. (2014) the value of R_{p} = 0.72 R_{⊕} as an upper limit for the radius of the planet solid body of KIC 12557548b, which corresponds to about 7300 km and 4600 km, respectively.
Although both planets are disintegrating and losing material, neither KOI 2700b nor KIC 12557548b show any evidence of a significant longterm orbital period variability. This is in agreement with the prediction, presented by PerezBecker & Chiang (2013), which suggests that disintegrating planets undergo negligible orbital period variations due to the evaporation process. A significant variability is observable in the lightcurve shape of the planets. In the case of KIC 12557548b there is a shortterm periodic variability in the transit core (Rappaport et al. 2012) and a longterm periodic variability in the egress (Budaj 2013). KOI 2700b exhibits only a longterm variability in the transit core, which may be periodic as well.
Based on our estimations, the massloss rate from KOI 2700b may be from 5.05 × 10^{7} to 4.41 × 10^{15} g s^{−1}. This interval of values is consistent with the upper limit of the massloss rate from KIC 12557548b, which was estimated by Rappaport et al. (2012), and which is 2 × 10^{11} g s^{−1}.
6 Conclusions
Our main scientific goal was to verify the disintegratingplanet scenario of KOI 2700b by modeling its light curve and put constraints on various tail and planet properties. The orbital period of the planet was improved based on the Kepler data from quarters 1–17. We obtained the value of P_{orb} = 0.9100259(15) days, which is in good agreement with the orbital period presented by the discoverers. We also searched for longterm orbital period changes of KOI 2700b, however, in agreement with PerezBecker & Chiang (2013), we found no significant evidence. We can confirm the decrease in transit depth of the exoplanet found by Rappaport et al. (2014). On the other hand, we cannot exclude the possibility that the transit depth varies periodically over a longterm time span longer than the available Kepler quarters.
We obtained the phasefolded and binned transit light curve of KOI 2700b, which we iteratively modeled using the radiativetransfer code SHELLSPEC. Mie absorption and scattering on spherical dust grains with realistic dust opacities, phase functions, and finite radius of the source of the scattered light were taken into account. During the modeling we applied selected species (γalumina, enstatite, forsterite, olivine with 50% magnesium and 50% iron and pyroxene with 40% magnesium and 60% iron) and dust particle sizes (0.01, 0.1 and 1 micron in radii). We executed 15 joint fitting procedures – one for each combinationof species and dust particle size. We modeled the cometlike tail as part of a ring around the parent star. The geometrical crosssection of the ring is monotonically enlarging from the planet to the end of the ring. Since this modified ring model of the tail satisfies the observations well, we confirm the disintegratingplanet scenario of KOI 2700b.
Furthermore, via modeling, we derived some interesting features of KOI 2700b and its cometlike tail. We first assumed that the orbital plane of the exoplanet and its tail had an inclination of i = 90 deg with respectto the plane of the sky, however it turns out that the orbital plane of the planet is not edgeon, but the orbital inclination angle is very probably close to the value of i = 87 deg. We obtained i in the range [85.1, 88.6] deg. KOI 2700b exhibits a relatively low dust density decreasing in its tail, which is defined by the density exponent A2 ≈−7. During our modeling process, we obtained this parameter in the range [−6.5, −8.2], which is in good agreement with Rappaport et al. (2014). We also derived the dust density at the beginning of the ring (ρ(0)) and the highest optical depth through the tail in front of the star (τ_{max}), based on a tail model with a crosssection of 0.05× 0.05 R_{⊙} at the beginning and 0.09 × 0.09 R_{⊙} at its end. The average values of the parameters are ρ(0) = 0.1667 × 10^{−15} g cm^{−3} and τ_{max} = 5.24 × 10^{−3}, and we obtained these parameters in the range [0.0151, 0.3333] × 10^{−15} g cm^{−3 } and [4.28, 7.43] × 10^{−3}, respectively. Since there is still ample room for a contribution to the light curve from the solid body of the planet, we also included it in the model. Our average value for the dimension of the planet is R_{p}∕R_{s} = 0.008 (R_{p } = 0.497 R_{⊕}, or 3172 km). Since we obtained this parameter in the range [0.0, 0.014], we can consider the value of R_{p } ∕R_{s} = 0.014, which corresponds to R_{p} = 0.871 R_{⊕} (or 5551 km), as an upper limit on the radius of the planet solid body of KOI 2700b. The upper limit of the planet radius is in agreement with the theory of the thermalwind and planet evaporation (PerezBecker & Chiang 2013) and is also derived with higher accuracy than the result found by Rappaport et al. (2014). We estimated the massloss rate from KOI 2700b, and we obtained Ṁ values from the interval [5.05 × 10^{7}, 4.41 × 10^{15}] g s^{−1 } , consistent with the estimate of Rappaport et al. (2014). Our average estimated massloss rate from KOI 2700b is Ṁ = 3.03 × 10^{14} g s^{−1 } . On the other hand, we could not draw any satisfactory conclusions on the typical grain size in the dust tail due to thesmall differences between the model pretransit brightenings at different particle radii. Finally we showed that the first exoplanet with a cometlike tail, KIC 12557548b, and KOI 2700b are similar in many respects.
Acknowledgements
I thank Dr. J. Budaj for the technical assistance, comments, and discussions. I also thank the anonymous referee for helpful comments and corrections. This work was supported by the Slovak Central Observatory Hurbanovo, by the VEGA grant of the Slovak Academy of Sciences No. 2/0031/18 and by the realization of the Project ITMS No. 26220120029, based on the Supporting Operational Research and Development Program financed from the European Regional Development Fund.
References
 Bochinski, J. J., Haswell, C. A., Marsh, T. R., Dhillon, V. S., & Littlefair, S. P. 2015, ApJ, 800, L21 [NASA ADS] [CrossRef] [Google Scholar]
 Bohren, C. F., & Huffman, D. R. 1983, Absorption and scattering of light by small particles (New York: John Wiley & Sons) [Google Scholar]
 Bourrier, V., Lecavelier des Etangs, A., Dupuy, H., et al. 2013, A&A, 551, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Brogi, M., Keller, C. U., de Juan Ovelar, M., et al. 2012, A&A, 545, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Budaj, J. 2013, A&A, 557, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Budaj, J., & Richards, M. T. 2004, Contributions Astron. Obs. Skalnate Pleso, 34, 167 [Google Scholar]
 Budaj, J., & Richards, M. T. 2010, in Binaries  Key to Comprehension of the Universe, eds. A. Prša & M. Zejda, ASP Conf. Ser., 435, 63 [NASA ADS] [Google Scholar]
 Budaj, J., Kocifaj, M., Salmeron, R., & Hubeny, I. 2015, MNRAS, 454, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Burrows, A., & Lunine, J. 1995, Nature, 378, 333 [NASA ADS] [CrossRef] [Google Scholar]
 Claret, A., & Bloemen, S. 2011, A&A, 529, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Croll, B., Rappaport, S., DeVore, J., et al. 2014, ApJ, 786, 100 [NASA ADS] [CrossRef] [Google Scholar]
 Deeming, T. J. 1975, Ap&SS, 36, 137 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Deirmendjian, D. 1964, Appl. Opt., 3, 187 [NASA ADS] [CrossRef] [Google Scholar]
 Eastman, J., Gaudi, B. S., & Agol, E. 2013, PASP, 125, 83 [NASA ADS] [CrossRef] [Google Scholar]
 Garai, Z., Zhou, G., Budaj, J., & Stellingwerf, R. F. 2014, Astron. Nachr., 335, 1018 [NASA ADS] [CrossRef] [Google Scholar]
 Guillot, T., Burrows, A., Hubbard, W. B., Lunine, J. I., & Saumon, D. 1996, ApJ, 459, L35 [NASA ADS] [CrossRef] [Google Scholar]
 Hubbard, W. B., Hattori, M. F., Burrows, A., Hubeny, I., & Sudarsky, D. 2007, Icarus, 187, 358 [NASA ADS] [CrossRef] [Google Scholar]
 Kawahara, H., Hirano, T., Kurosaki, K., Ito, Y., & Ikoma, M. 2013, ApJ, 776, L6 [NASA ADS] [CrossRef] [Google Scholar]
 Lamers, H. J. G. L. M., Lecavelier Des Etangs, A., & VidalMadjar, A. 1997, A&A, 328, 321 [NASA ADS] [Google Scholar]
 Lecavelier Des Etangs, A., VidalMadjar, A., & Ferlet, R. 1999, A&A, 343, 916 [NASA ADS] [Google Scholar]
 Lecavelier Des Etangs, A., Ehrenreich, D., VidalMadjar, A., et al. 2010, A&A,514, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mura, A., Wurz, P., Schneider, J., et al. 2011, Icarus, 211, 1 [NASA ADS] [CrossRef] [Google Scholar]
 PerezBecker, D., & Chiang, E. 2013, MNRAS, 433, 2294 [NASA ADS] [CrossRef] [Google Scholar]
 Rappaport, S., Levine, A., Chiang, E., et al. 2012, ApJ, 752, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Rappaport, S., Barclay, T., DeVore, J., et al. 2014, ApJ, 784, 40 [NASA ADS] [CrossRef] [Google Scholar]
 Rappaport, S., Vanderburg, A., Jacobs, T., et al. 2018, MNRAS, 474, 1453 [NASA ADS] [CrossRef] [Google Scholar]
 SanchisOjeda, R., Rappaport, S., Pallè, E., et al. 2015, ApJ, 812, 112 [NASA ADS] [CrossRef] [Google Scholar]
 Schlawin, E., Herter, T., Zhao, M., Teske, J. K., & Chen, H. 2016, ApJ, 826, 156 [NASA ADS] [CrossRef] [Google Scholar]
 Schneider, J., Rauer, H., Lasota, J. P., Bonazzola, S., & Chassefiere, E. 1998, in Brown Dwarfs and Extrasolar Planets, eds. R. Rebolo, E. L. Martin, & M. R. Zapatero Osorio, ASP Conf. Ser., 134, 241 [NASA ADS] [Google Scholar]
 Šejnová, K., Votruba, V., & Koubská, P. 2012, in Circumstellar Dynamics at High Resolution, eds. A. C. Carciofi, & T. Rivinius, ASP Conf. Ser., 464, 219 [NASA ADS] [Google Scholar]
 Southworth, J., Maxted, P. F. L., & Smalley, B. 2004, MNRAS, 351, 1277 [NASA ADS] [CrossRef] [Google Scholar]
 Stellingwerf, R. F. 1978, ApJ, 224, 953 [NASA ADS] [CrossRef] [Google Scholar]
 Tian, F., Toon, O. B., Pavlov, A. A., & De Sterck, H. 2005, ApJ, 621, 1049 [NASA ADS] [CrossRef] [Google Scholar]
 Tkachenko, A., Lehmann, H., Tsymbal, V., & Mkrtichian, D. 2008, Commun. Asteroseismol., 157, 377 [NASA ADS] [Google Scholar]
 van Lieshout, R., Min, M., & Dominik, C. 2014, A&A, 572, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 van Lieshout, R., Min, M., Dominik, C., et al. 2016, A&A, 596, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 van Werkhoven, T. I. M., Brogi, M., Snellen, I. A. G., & Keller, C. U.. 2014, A&A, 561, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 VidalMadjar, A., Lecavelier des Etangs, A., Désert, J.M., et al. 2003, Nature, 422, 143 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 VidalMadjar, A., Désert, J.M., Lecavelier des Etangs, A., et al. 2004, ApJ, 604, L69 [NASA ADS] [CrossRef] [Google Scholar]
 Yelle, R. V. 2004, Icarus, 170, 167 [NASA ADS] [CrossRef] [Google Scholar]
Current version of the software package available at http://www.stellingwerf.com
The code is available at https://svn.ssec.wisc.edu/repos/geoffc/Mie/bhmief/callbhmie.f
The software is available at https://www.ta3.sk/budaj/dust/deirm/diskaver/
The software is available at https://www.ta3.sk/budaj/shellspec.html
The software is available at http://www.astro.keele.ac.uk/jkt/codes/jktebop.html
All Tables
All Figures
Fig. 1 Comparison between the SAP and PDCSAP Kepler phasefoldedand binned transit light curves of KOI 2700b. The illustrative vertical error bars on the left side represent the median uncertainties of data points depicted with the same color. 

In the text 
Fig. 2 Averaged light curve of KOI 2700b, smoothed using the running window technique. The window width and step was 0.011 and 0.0011, respectively (in units of phase). The light curve shows systematic fluctuations at short timescales (20–30 min). See the text. 

In the text 
Fig. 3 Results of the period analysis of KOI 2700b. We searched for periodicities from 1 min to 1 h using the method of phase dispersion minimization (see Sect. 3), but we did not find any significant frequency inthis interval. See the text. 

In the text 
Fig. 4 Phasefolded and binned transit light curve of KOI 2700b. Red points are observations and green points represent the averaged data. We note that averaged data points have uncertainties that are too small to be discerned in this plot. Error bars of these data are well depicted, for example, in Fig 9. 

In the text 
Fig. 5 Comparison between the SAP Kepler phasefolded and binned transit light curve of exoplanets KIC 12557548b and KOI 2700b. The illustrative vertical error bars on the left side represent the median uncertainties of data points depicted with the same color. 

In the text 
Fig. 6 Transit depth decreasing at KOI 2700b. The phasefolded and binned transit light curves of the exoplanet were constructed from 3 and 4 quarters (Q) of detrended data. The illustrative vertical error bar on the left side represents the median uncertainty of all data points. 

In the text 
Fig. 7 Search for a longterm orbital period change of KOI 2700b. The minimum of the parameter Θ _{min } corresponds to the periodchange speed β = 1.3 ± 3.2 d/Myr, which means that there is no significant evidence for the longterm orbital period change. 

In the text 
Fig. 8 Phase functions at 0.6 micron for different dust particle radii of γalumina. Lines without symbols are phase functions assuming a point source of light. Lines with symbols are phase functions that take into account the finite dimension of the stellar disk. The vertical line illustrates the angular dimension of the stellar disk as seen from KOI 2700b. 

In the text 
Fig. 9 Model light curves calculated for 0.01micron, 0.1micron, and 1micron grains of γalumina (left) and enstatite (right), compared with the observed light curve of KOI 2700b. 

In the text 
Fig. 11 As in Fig. 10, but for pyroxene grains. 

In the text 
Fig. 10 Model light curves calculated for 0.01micron, 0.1micron, and 1micron grains of forsterite (left) and olivine (right), compared with the observed light curve of KOI 2700b. 

In the text 
Fig. 12 Phasefolded and binned transit light curve of KOI 2700b (observations) compared with the transit light curve of the solidbody of KOI 2700b (according to the average value and upper limit of the planet radius; see the text). Furthermore, the transit light curve of an Earthsize planet and a Mercurysize planet is also plotted. The orbital inclination angle is i = 87 deg. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.