Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://lnfm1.sai.msu.ru/~marat/articles/eng/rp_lam_Abubekerov_Gostev_2015_eng.pdf Äàòà èçìåíåíèÿ: Mon Jan 12 11:58:16 2015 Äàòà èíäåêñèðîâàíèÿ: Sun Apr 10 17:12:39 2016 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: þæíàÿ àòëàíòè÷åñêàÿ àíîìàëèÿ

ISSN 1063-7729, Astronomy Reports, 2015, Vol. 59, No. 1, pp. 1­11. c Pleiades Publishing, Ltd., 2015. Original Russian Text c M.K. Abubekerov, N.Yu. Gostev, A.M. Cherepashchuk, 2015, published in Astronomicheskii Zhurnal, 2015, Vol. 92, No. 1, pp. 3­14.

Limb-Darkening Anomalies in Stars Eclipsed by Exoplanets
M. K. Abubekerov* , N. Yu. Gostev, and A. M. Cherepashchuk
Sternberg Astronomical Institute, Moscow State University, Moscow, Russia
Received March 26, 2014; in final form, June 9, 2014

Abstract--The possibility of explaining discrepancies between observed and theoretical limb-darkening coefficients as being a manifestation of the extended, non-spherical atmosphere of an eclipsing exoplanet is analyzed, assuming the model stellar atmospheres used are correct. DOI: 10.1134/S1063772915010011

1. INTRODUCTION A detailed analysis of high-accuracy, multi-color light curves of the binary system with an exoplanet HD 209458 [1] demonstrated a discrepancy between the observed limb-darkening coefficients and theoretical values obtained from model stellar atmospheres. A more detailed study of the differences between observational and theoretical limb-darkening coefficients for HD 209458 showed that this difference can not be avoided, even when a fairly wide range of parameters of the model stellar atmosphere is considered [2]. Claret [2] referred to HD 209458 as a "challenge to the theory of stellar atmospheres." We attempted to identify a way of removing the discrepancy between the observed and theoretical limbdarkening coefficients in our earlier paper[3], where we showed that even the use of confidence intervals to estimate the "external" errors is not able to remove this discrepancy. L It is possible that the discrepancy is due to the use of a model that is too simple. For example, a model with two spherical stars is not able to fully take into account the interaction between a star and its planet. The importance of taking into account the effects of the interaction between the star and exoplanet has been confirmed by new observational data. For example, analysis of a transit near-UV light curve of the binary with an exoplanet WASP-12b showed a much deeper and wider eclipse compared to the photometric transit light curve in the visible [4]. This suggests deformation of the planetary magnetosphere and atmosphere under the action of the stellar wind [5, 6]. The WASP-12b system is composed of a late 1.35 M star with a surface temperature of 6300 K and a 1.4 MJ exoplanet located at a distance of 0.0229 AU (MJ is the mass of Jupiter). The orbital period of
*

E-mail: marat.abubekerov@gmail.com

the system is 1.09d . These parameters are fairly close to the parameters of binaries we have studied-- HD 209458 and HD 189733. The possible deformation of the magnetosphere and atmosphere of exoplanets due to the stellar winds of their central stars is also confirmed by three-dimensional gas-dynamical modeling of the behavior of the atmosphere of the hot Jupiter in WASP-12b, which is located close to its star [7]. The WASP-12b system resembles the binary systems HD 209458 and HD 189733. The results of the 3D gas-dynamical modeling [7, 8] indicate that the atmosphere of the exoplanet HD 209458b is perturbed by the stellar wind. Due to the comparatively high temperature of the planetary atmosphere, the outflow of matter from its surface is sufficiently strong for the planet to acquire a cometary tail [9, 10]. One of the reasons for this strong disturbance by the stellar wind of the central star is the weakness of the exoplanet's magnetic field, which is unable to withstand the stellar wind. Our previous studies [3, 11­14] treated the planet like a spherical body with a sharp edge. The new analyses of [4, 10] suggest that this idealized model is oversimplified. We thus cannot rule out that the discrepancy between the observed and theoretical limbdarkening coefficients of HD 209458 and HD 189733 is due to the application of model for the explanet atmospheres that is too idealized, rather than inadequacy of current models for thin stellar atmospheres. To investigate this question, we fixed the limbdarkening coefficients when fitting the eclipse data for HD 209458 and HD 189733 in our current study. We used values yielded by the latest theoretical computations of Claret [15] obtained using PHOENIX model for the limb-darkening coefficients of HD 209458 and HD 189733. Assuming spherical model components moving in circular orbits, we searched for the orbital inclinaton i, relative radius of the exoplanet rp , and
1


2

ABUBEKEROV et al.

Table 1. Results of fitting the observed light curves of HD 209458 [16], using a quadratic limb-darkening law with fixed coefficients (errors are 1 , and were estimated using the differential-corrections method; the two rightmost columns list the reduced 2 and corresponding significance levels ) ° , A 3201 3750 4300 4849 5398 5802 6779 7755 8732 9708 rs 0.1131 0.1102 0.1112 0.1109 0.1117 0.1122 0.1123 0.1112 0.1119 0.1107 (rs ) 0.00053 0.00028 0.00018 0.00016 0.00016 0.00015 0.00014 0.00016 0.00020 0.00030 rp 0.01385 0.01323 0.01335 0.01326 0.01341 0.01356 0.01350 0.01336 0.01347 0.01339 (rp ) 0.000062 0.000030 0.000019 0.000018 0.000018 0.000016 0.000015 0.000019 0.000024 0.000035 rp /rs 0.12239 0.12003 0.12003 0.11958 0.12005 0.12092 0.12019 0.12010 0.12043 0.12096 (rp /rs ) 0.0046 0.0026 0.0016 0.0014 0.0014 0.0013 0.0012 0.0014 0.0018 0.0027 i 86.71 87.17 87.04 87.11 87.00 86.89 86.89 86.98 86.91 87.04 (i) 0.035 0.022 0.013 0.013 0.013 0.012 0.012 0.015 0.018 0.028
2 red

0.16 0.07 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00

1.07 1.10 1.19 1.30 1.20 1.46 1.46 1.32 1.20 1.23

The closeness of to 0 is related to the fixed coefficents in the quadratic limb-darkening law.

relative radius of the star rs as a function of the wavelength . The parameters i and, to an appreciable extent, rs should not depend on the wavelength. At the same time, if the exoplanet has a significant atmosphere, the radius rp may depend substantially on . If a wavelength dependence is also found for both i and rs , this implies inadequacy of the model used to describe the eclipsing binary, due to failure to take into account the exoplanet atmosphere, its cometary tail, and so on. 2. OBSERVATIONAL DATA AND FITTING We consider here the observed light curves of the binary systems HD 209458 and HD 189733, which have similar parameters. Despite this similarity, HD 189733 displays a satisfactory agreement between the observed and theoretical limb-darkening coefficients, while HD 209458 does not. As observational material for HD 209458, we used multi-color light curves obtained with the Hubble Space Telescope (HST) from May 3 to July 6, 2003 [16]. The central wavelengths of these transit light curves are = 3201, 3750, 4300, 4849, 5398, ° 5802, 6779, 7755, 8732, and 9708 A. Each light curve includes 500 individual brightness measurements for HD 209458. The rms error of the individual measurements varies from 1.79 â 10-4 to 6.09 â 10-4 (from 10-2 to 3 â 10-2 relative to the depth of the eclipse). More detailed information on these observational data and their reduction is provided in [16]. The observed transit light curves for HD 209458 are presented in Fig. 1.

As observational material for HD 189733, we used multi-color light curves obtained with the HST on May 22, May 26, and July 14, 2006 [17]. Ten light curves of eclipses in the ranges = 5500-6000, 6000-6500, 6500-7000, 7000-7500, 7500-8000, 8000-8500, 8500-9000, 9000-9500, 9500-10 000, ° 10 000- 10 500 A were obtained. More detailed information on these observational data can be found in [17]. In our analyses of these light curves, we adopted the central wavelengths = 5750, 6250, 6750, 7250, 7750, 8250, 8750, 9250, 9750, and ° 10 250 A. Each light curve contains 675 individual brightness measurements with uncertainties 5 â 10-4 of the normalized out-of-eclipse brightness. The observed transit light curves of HD 189733 are presented in Fig. 2. We fitted the transit light curves using a model with a spherical star eclipsed by a spherical exoplanet (the classical problem with two spherical stars), without taking into account the effects of the proximity of the two components. Such a model can easily be computed on modern computers, making it possible to solve a large number of model versions using a relatively small amount of computer time [11]. The model and algorithm used to fit the transit light curves are described in more detail in [3, 11­14]. Since the transit light curve of HD 189733 displays a disturbance due to the spottedness of the stellar surface [18] close to the minimum and on the right branch (Fig. 2), only the left branch of the light curve, up to orbital phase 178 , was fitted. It was shown in [13] that the left branch of the light curve is
ASTRONOMY REPORTS Vol. 59 No. 1 2015


STARS ECLIPSED BY EXOPLANETS

3

I 1.04

1.02

1.00

0.98

0.96

0.94

0.92 165 170 175 180 Phase angle, deg 185 190 195

Fig. 1. Observed light curves of HD 209458 [16] constructed for wavelengths 3201, 3750 4300, 4849 5398, 5802, 6779, 7755, ° 8732, and 9708 A (upper set of curves, with the wavelength increasing upward). The lower set of curves show the corresponding residuals. The solid curves show theoretical curves obtained for the model with a nonlinear (quadratic) limb-darkening law and the fixed coefficients from [15].

most reliable for the observational data we are using here. 3. RESULTS The results of our fitting of the transit light curves of HD 209458 are presented in Table 1, and are shown graphically in Figs. 3­6. The results of our fitting of the transit light curves of HD 189733 are presented in Table 2 and Figs. 7­10.
ASTRONOMY REPORTS Vol. 59 No. 1 2015

As can be seen from Table 1, the critical significance levels for the models derived from the transit ° light curves of HD 209458 with 4300 Aare close to zero. Thus, these models fit the observational data relatively poorly. This is most probably due to the fact that the limb-darkening coefficients were fixed. Note that this also affects the wavelength dependences of the derived parameters. The dependences presented


4

ABUBEKEROV et al.

I 1.04

1.02

1.00

0.98

0.96

0.94

0.92 172.5 175.0 177.5 180.0 182.5 Phase angle, deg 185.0 187.5 190.0

Fig. 2. Observed light curves of HD 189733 at = 5500-10 500 (upper set of curves, with the wavelength increasing upward). The lower set of curves shows the residuals for a quadratic limb-darkening law for the best-fit values of the parameters and the fixed darkening coefficients from [15].

in Figs. 3­6 show that the points are generally scattered irregularly. The spherical model [3] satisfies the observational data for HD 189733 significantly better. The wavelength dependences of the parameters shown in Figs. 7­10 are more regular, and we can distinguish a single trend in each case. Since the model derived

using the light curve for the central wavelength = ° 5750 A is rejected by the 2 criterion, the parameters red obtained from this light curve were not used in our overall analysis. We used linear fits and their rms values to find and analyze the wavelength dependences of the orbital
ASTRONOMY REPORTS Vol. 59 No. 1 2015


STARS ECLIPSED BY EXOPLANETS

5

i, deg 87.2 87.1 87.0 86.9 86.8 86.7 86.6 3000

4000

5000

6000

7000

8000

9000

10 000 , å

Fig. 3. Dependence of i on the photometric wavelength, obtained for HD 209458 by fitting the transit light curves using a quadratic limb-darkening law with the fixed coefficients [15]. The solid line shows the least-squares fit.

rp, deg 0.0140 0.0139 0.0138 0.0137 0.0136 0.0135 0.0134 0.0133 0.0132 0.0131 3000 4000 5000 6000 7000 8000 9000 10 000 , å

Fig. 4. Dependence of rp on the photometric wavelength, obtained for HD 209458 by fitting the transit light curves using a quadratic limb-darkening law with the fixed coefficients [15]. The solid line shows the dependence obtained for the entire range ° ° ° = 3201-9708 A, while the dashed lines show the least-squares fits for < 5000 Aand > 5000 A.

ASTRONOMY REPORTS

Vol. 59 No. 1 2015


6

ABUBEKEROV et al.

rp/rs 0.128 0.127 0.126 0.125 0.124 0.123 0.122 0.121 0.120 0.119 0.118 0.117 0.116 3000

4000

5000

6000

7000

8000

9000

10 000 , å

Fig. 5. Same as Fig. 4for theratio rp /rs .

rs, deg 0.1140 0.1135 0.1130 0.1125 0.1120 0.1115 0.1110 0.1105 0.1100 0.1095 0.1090 3000 4000 5000 6000 7000 8000 9000 10 000 , å

Fig. 6. Same as Fig. 4 for rs .

inclination i, exoplanet radius rp , and stellar radius r for HD 209458 and HD 189733.

s

on are more pronounced, although the scatter of the points on the corresponding trend is fairly large. Figure 4 shows that the scatter of points relative to the trend indicated by the rms deviations about the linear fit is fairly large. Deviations from the trend ° are particularly large for < 5000 A. This led us to analyze the trends for the dependence of rp on
ASTRONOMY REPORTS Vol. 59 No. 1 2015

4. DISCUSSION 4.1. HD 209458 For HD 209458, we note only a weak dependence of i on (Fig. 3). The dependences of rp , rs ,and rp /rs


STARS ECLIPSED BY EXOPLANETS

7

i, deg 86.2 86.1 86.0 85.9 85.8 85.7 85.6 85.5 6000

7000

8000

9000

10 000

11 000 , å

Fig. 7. Same as Fig. 3 for HD 189733.

rp, deg 0.0179 0.0178 0.0177 0.0176 0.0175 0.0174 0.0173 0.0172 6000

7000

8000

9000

10 000

11 000 , å

Fig. 8. Dependence of rp on the photometric wavelength, obtained for HD 189733 using a quadratic limb-darkening law with the fixed coefficients [15]. The solid line shows the least-squares fit.

° ° separately in the range < 5000 A and > 5000 A (Fig. 4). We proceeded in the same way in our analysis of the ratio of the exoplanet and stellar radii rp /rs and the stellar radius rs (Figs. 5 and 6). The wavelength dependence of the exoplanet radius rp for ° HD 209458b is stronger for < 5000 Athan for >
ASTRONOMY REPORTS Vol. 59 No. 1 2015

° 5000 A (Fig. 4), although the dependence changes its ° sign if the point corresponding to = 3201 A is not included. The noted regions in Fig. 4 differ qualitatively in the mean slopes of the fitted lines. As expected, the ° dependence is steeper for < 5000 A than for > ° 5000 A (Fig. 4). This break confirms that the effects


8

ABUBEKEROV et al.

rp/rs 0.162 0.161 0.160 0.159 0.158 0.157 0.156 0.155 0.154 0.153 0.152 6000 7000 8000 9000 10 000 11 000 , å

Fig. 9. Same as Fig. 8for theratio rp /rs .

rs, deg 0.1135 0.1130 0.1125 0.1120 0.1115 0.1110 0.1105 6000

7000

8000

9000

10 000

11 000 , å

Fig. 10. Same as Fig. 8 for rs .

of absorption by the distorted exoplanet atmosphere are important in the blue part of the spectrum. According to the recent studies [7, 8], due to its weak magnetic field and interaction with the stellar wind from its central star, the hot Jupiter in the HD 209458b system is significantly non-spherical, and has a gas plume similar to a cometary tail. Thus,

the dependence shown in Fig. 4 fits the picture of the exoplanet HD 209458b as a cometary planet [7, 8]. Our conclusions concerning HD 209458 are in qualitative agreement with the conclusions derived by Fossatti et al. [4], who showed that the WASP-12b transit light curve in the UV differs significantly from the visual light curve, in terms of both its duration and the depth of the eclipse.
ASTRONOMY REPORTS Vol. 59 No. 1 2015


STARS ECLIPSED BY EXOPLANETS

9

­ ­ ­ ­ ­

x 1.5 1.4 1.3 1.2 1.1 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0.1 0.2 0.3 0.4 0.5 3000 4000 5000 6000 7000 8000 9000 10 000 11 000 12 000 13 000 , å

Fig. 11. Dependence of the observational [12] and theoretical [15] linear coefficient of the quadratic limb-darkening law for HD 209458 on the central photometric wavelength. The uncertainties in the limb-darkening coefficient x were obtained using the differential-corrections method. The errors shown are 1 .

­ ­ ­ ­ ­ ­

x 1.5 1.4 1.3 1.2 1.1 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 3000 4000 5000 6000 7000 8000 9000 10 000 11 000 12 000 13 000 , å

Fig. 12. Dependence of the observational [13] and theoretical [15] linear coefficients of the quadratic limb-darkening law for HD 189733 on the central photometric wavelengths. The uncertainties in the limb-darkening coefficient x were obtained using the differential-corrections method. The errors shown are 1 .

This suggests that the absorption in the visual ° at > 5000 A is lower than in the UV, so that the ° slope of the line representing the trend for > 5000 A in Fig. 4 is appreciably different from the slope for
ASTRONOMY REPORTS Vol. 59 No. 1 2015

° < 5000 A. The wavelength dependence of the exo° planet radius becomes less pronounced at > 5000 A (Fig. 4). For completness, Fig. 11 shows the wavelength dependence of the observational values of the linear


10

ABUBEKEROV et al.

Table 2. Results of fitting the observed light curves (left branch) of HD 189733 [18]. For remaining comments on the table, see Table 1. ° , A 5750 6250 6750 7250 7750 8250 8750 9250 9750 10 250
Values
2 red

rs 0.1102 0.1127 0.1121 0.1116 0.1118 0.1123 0.1123 0.1117 0.1110 0.1118

(rs ) 0.00028 0.00039 0.00024 0.00020 0.00020 0.00021 0.00021 0.00024 0.00032 0.00038

rp 0.01717 0.01775 0.01764 0.01750 0.01752 0.01762 0.01756 0.01743 0.01729 0.01746

(rp ) 0.000061 0.000092 0.000057 0.000048 0.000047 0.000049 0.000049 0.000055 0.000073 0.000090

rp /rs 0.1558 0.1576 0.1574 0.1568 0.1567 0.1569 0.1564 0.1561 0.1558 0.1561

(rp /rs ) 0.0025 0.0034 0.0021 0.0018 0.0017 0.0018 0.0018 0.0021 0.0028 0.0034

i 86.03 85.72 85.72 85.77 85.75 85.71 85.71 85.79 85.84 85.75


(i) 0.018 0.024 0.015 0.013 0.013 0.014 0.014 0.016 0.022 0.026



2 red

0.00 1.00 1.00 0.98 0.55 0.0005 0.81 0.61 0.054 1.00

5.56 0.58 0.52 0.69 1.02 1.80 0.88 0.99 1.38 0.62

< 1 may indicate correlations between the data points on the light curves.

coefficient x of the quadratic limb-darkening law obtained earlier in [12]. This figure also shows the latest theoretical values for the limb-darkening coefficients obtained by Claret [15]. Figure 11 shows that the theoretical and observational values are in agreement ° for short wavelengths ( < 5000 A), but disagree at ° > 5000 A. This all leads us to conclude that, although the ° long-wavelength range > 5000 A is less affected by absorption by the extended atmosphere of the exoplanet in the case of HD 209458 (including the cometary tail), the observational values of the linear darkening coefficient in the quadratic law do not agree ° with the theoretical values for > 5000 A. The reason for this discrepancy remains unclear.

according to a Rayleigh law. The fact that the stellar radius rs varies only slightly in this wavelength range (by 0.7%, see Fig. 10), suggests that the dependence of rp on we have found is fairly reliable. This dependence confirms the action of Rayleigh scattering in the exoplanet atmosphere [13, 17]. We carried out a comparative analysis of the previously obtained linear limb-darkening coefficients in the quadratic law x [13] with the latest theoretical limb-darkening coefficients computed by Claret [15]. The observational and theoretical limb-darkening coefficients are in good agreement for HD 189733 (see Fig. 12), apart from the limb-darkening coefficient x ° based on the transit light curve for = 5750 A, for 2 criterion which the model was rejected by the red (Table 2). 5. CONCLUSION We have analyzed multi-color transit light curves of the star­exoplanet binary systems HD 209458 and HD 189733 using the theoretical limb-darkening coefficients of [15] and a quadratic limb-darkening law. Supposing the model for the stellar atmosphere to be correct, we have checked whether the discrepancy between the observational [3, 11­14] and theoretical [2, 15] limb-darkening coefficients can be explained by the influence of an extended, non-spherical atmosphere around the exoplanet. In the case of HD 209458, the wavelength dependence of i, rs , and rp is most significant at short wavelengths, = ° 3201-5000 A, while these dependences are weaker ° at long wavelengths, = 5000-9708 A. This most
ASTRONOMY REPORTS Vol. 59 No. 1 2015

4.2. HD 189733
In the case of HD 189733, the model with two spherical bodies adequately describes the observed transit light curves according to the 2 criterion (apart from the transit light curve with the central ° wavelength = 5750 A, which was not included in our analysis of the results of the fitting). Orbital inclination of HD 189733 depends on the wavelength only weakly (Fig. 7). Unlike the exoplanet in the HD 209458 system, the wavelength dependence of the radius of the HD 209458 exoplanet is more regular. The radius of the exoplanet rp varies ° by 2.5% at wavelengths = 6000-10 000 A, with rp monotonically decreasing with increasing wavelength, providing evidence for the presence of an exoplanet atmosphere that scatters the stellar radiation


STARS ECLIPSED BY EXOPLANETS

11

likely indicates that the effects of the absorption of the stellar radiation by an extended, distorted atmosphere of the exoplanet are manifest at short wavelengths only, and are small at the longer wavelengths. In this connection, it would seem that the observational values of the limb-darkening coefficients should be most different from the theoretical val° ues at short wavelengths = 3201-5000 A [2, 15]. However, the situation is the opposite in the case of HD 209458: there is good agreement of the observational and theoretical limb-darkening coefficien at short wavelengths (Fig. 11), while a significant discrepancy is observed at longer wavelengths. The reason for this discrepancy is unclear, leading us to acknowledge that HD 209458 remains a "challenge for the theory of stellar atmospheres," as claimed by Claret [2]. In the case of HD 189733, the wavelength dependence of i is weak, and the relative variation of the stellar radius rs with wavelength is less than 1%, indicating that the effects of absorption of the stellar radiation by a distorted, extended exoplanet atmosphere are insignificant at long wave° lengths, = 6000-10 000 A, and a simple model for the binary system with two spherical bodies in circular orbits provides an adequate description of the available observational data. However, since the radius of the eclipsing exoplanet decreases with increasing wavelength by 2.5% at ° = 6000-10 000 A, this suggests that the exoplanet does have a quasi-atmosphere that scatters light according to a Rayleigh law. The observational limbdarkening coefficients [13] in the quadratic limbdarkening law are consistent with the theoretical ° values [2, 15] in the range = 6000-10 000 A in the case of HD 189733, confirming the validity of the modern theory of thin stellar atmospheres [15]. Further careful study of limb darkening in stars eclipsed by exoplanets are of considerable interest. ACKNOWLEDGMENTS This work was supported by a grant of the President of the Russian Federation for the State Support of Young Russian Phd Scientists (MK-893.2012.2) and the Russian Foundation for Basic Research (project 12-02-31466).

REFERENCES
1. J. Southworth, Mon. Not. R. Astron. Soc. 386, 1644 (2008). 2. A. Claret, Astron. Astrophys. 506, 1335 (2009). 3. M. K. Abubekerov, N. Yu. Gostev, and A. M. Cherepashchuk, Astron. Rep. 52, 99 (2008). 4. L. Fossati, C. A. Haswell, C. S. Froning, L. Hebb, S. Holmes, U. Kolb, Ch. Helling, A. Carter, P. Wheatley, A. C. Cameron, B. Loeillet, D. Pollacco, R. Street, H. C. Stempels, E. Simpson, S. Udry, Y. C. Joshi, R. G. West, I. Skillen, and D. Wilson, Astrophys. J. Lett. 714, L222 (2010). 5. J. Llama, K. Wood, M. Jardine, A. A. Vidotto, Ch. Helling, L. Fossati, and C. A. Haswell, Mon. Not. R. Astron. Soc. 416, L41 (2011). 6. A. A. Vidotto, M. Jardine, and Ch. Helling, Mon. Not. R. Astron. Soc. 414, 1573 (2011). 7. D. Bisikalo, P. Kaygorodov, D. Ionov, V. Shematovich, H. Lammer, and L. Fossati, Astropys. J. 764, 19 (2013). 8. D. V. Bisikalo, P. V. Kaygorodov, D. E. Ionov, and V. I. Shematovich, Astron. Rep. 57, 715 (2013); arXiv:1311.4441v2 [astro-ph.EP] (2013). 9. E. M. Schneiter, P. F. Velazquez, A. Esquivel, and A. C. Raga, Astrophys. J. Lett. 671, L57 (2007). 10. V. Bourrier and A. Lecavelier des Etangs, Astron. Astrophys. 557, A124 (2013). 11. M. K. Abubekerov, N. Yu. Gostev, and A. M. Cherepashchuk, Astron. Rep. 53, 722 (2009). 12. M. K. Abubekerov, N. Yu. Gostev, and A. M. Cherepashchuk, Astron. Rep. 54, 1005 (2010). 13. M. K. Abubekerov, N. Yu. Gostev, and A. M. Cherepashchuk, Astron. Rep. 55, 1051 (2011). 14. M. K. Abubekerov and N. Yu. Gostev, Mon. Not. R. Astron. Soc. 432, 2216 (2013). 15. A. Claret, Astron. Astrophys. 552, A16 (2013). 16. H. A. Knutson, D. Charbonneau, R. W. Noyes, T. M. Brown, and R. L. Gilliland, Astrophys. J. 655, 564 (2007). 17. F. Pont, R. L. Gilliland, C. Moutou D. Charbonneau, F. Bouchy, T. M. Brown, M. Mayor, D. Queloz, N. Santos, and S. Udry, Astron. Astrophys. 476, 1347 (2007). 18. F. Pont, H. Knutson, R. L. Gilliland, C. Moutou, and D. Charbonneau, Mon. Not. R. Astron. Soc. 385, 109 (2008).

Translated by L. Yungel'son

ASTRONOMY REPORTS

Vol. 59 No. 1 2015