Issue 
A&A
Volume 570, October 2014



Article Number  L5  
Number of page(s)  6  
Section  Letters  
DOI  https://doi.org/10.1051/00046361/201424733  
Published online  16 October 2014 
Transit light curve and inner structure of closein planets^{⋆}
^{1}
Departamento de Física, I3NUniversidade de Aveiro
Campus de Santiago
3810193
Aveiro
Portugal
email:
correia@ua.pt
^{2}
ASD, IMCCECNRS UMR8028, Observatoire de Paris,
UPMC, 77 Av.
DenfertRochereau, 75014
Paris,
France
Received:
1
August
2014
Accepted:
25
September
2014
Planets orbiting very close to their host stars have been found, some of them on the verge of tidal disruption. The ellipsoidal shape of these planets can significantly differ from a sphere, which modifies the transit light curves. Here we present an easy method for taking the effect of the tidal bulge into account in the transit photometric observations. We show that the differences in the light curve are greater than previously thought. When detectable, these differences provide an estimation of the fluid Love number, which is invaluable information on the internal structure of closein planets. We also derive a simple analytical expression to correct the bulk density of these bodies, that can be 20% smaller than current estimates obtained assuming a spherical radius.
Key words: planetary systems / techniques: photometric / planets and satellites: interiors
Appendix A is available in electronic form at http://www.aanda.org
© ESO, 2014
1. Introduction
About half of the more than 1000 known transiting exoplanets have orbital periods less than ten days^{1}. These closein planets undergo strong tidal effects raised by the parent star. One consequence is that their spins and orbits evolve until an equilibrium configuration is reached, corresponding to coplanarity, circularity, and synchronous rotation (e.g., Hut 1980; Correia 2009). Another consequence is that the shape of these planets differ from a spherical body, and it is approximated better by a triaxial ellipsoid (e.g., Chandrasekhar 1987; Correia & Rodríguez 2013). The asymmetry in the mass distribution increases with the proximity to the star, and it is particularly pronounced near the Roche limit (e.g., FerrazMello et al. 2008; Burton et al. 2014).
For simplicity, most observational works on transiting planets assume that its shape is spherical, so they determine an average radius. However, if there is enough precision in the data, it is possible to spot the polar oblateness (f_{c}) signature in the transit light curve (Seager & Hui 2002; Barnes & Fortney 2003; Ragozzine & Wolf 2009; Carter & Winn 2010), which gives us invaluable information on their internal structure. Previous studies have ignored the equatorial prolateness (f_{a}), which is actually more pronounced than the polar oblateness, since the long axis always points to the star, and thus should not be perceptible during the transit. Leconte et al. (2011) and Burton et al. (2014) have shown that failing to account for this distortion would result in a systematic underestimation of the planetary radius, hence of its bulk density. In addition, for planets near the Roche limit, the projected ellipsoid also depends on the inclination to the line of sight and on the rotation angle, which most likely modifies the transit light curve. Unfortunately, there is no comparison of transit light curves for ellipsoidal and spherical planets shown in the works by Leconte et al. (2011) and Burton et al. (2014).
Since a large number of transiting planets are being detected on the verge of tidal disruption (e.g., Valsecchi & Rasio 2014), it becomes important to understand the exact contribution of its shape to the photometric observations. The work by Leconte et al. (2011) is very complete, but as pointed by Burton et al. (2014), it requires complex internal structure models that are difficult to implement and to reconcile with observational parameters. Therefore, Burton et al. (2014) compute the tidal deformation by assuming surfaces of constant gravitational equipotential for the planet that are solely based on observable parameters. Nevertheless, this method is also not easy to implement, since it requires a numerical adjustment, such that the projected area of the model planet during the transit matches those given by the observations. Moreover, we lose the information relative to the internal structure, which is an important complement to the density determination.
In this Letter, we propose a simple analytical model for computing the projected area of closein planets at any point of its orbit, which is based in the equilibrium surface given by secondorder Love numbers. We thus obtain the transit light curve for these planets, which can be used to compare directly with the observations, and infer their internal structure and density.
2. Shape
The shape (or the figure) of a planet is usually described well by a reference ellipsoid (quadrupolar approximation). The standard equation of a triaxial ellipsoid centered on the origin of a Cartesian coordinate system (X,Y,Z) and aligned with the axes is (1)where a, b, and c are called the semiprincipal axes. For the Earth and the gaseous planets in the solar system, we have a ≃ b>c (oblate spheroids), but usually a>b>c for the main satellites, where the long axis a is directed to the central planet.
We let R = (X,Y,Z) be a generic point at the surface of the ellipsoid. Then, Eq. (1) can be rewritten as (2)where ^{T} denotes the transpose. Solving Eq. (1) for Y^{2}, the radial distance R =  R  can be expressed as (3)where (4)and (5)Assuming a ~ b ~ c, we can neglect terms in and in the previous expressions, and thus (6)The mass distribution inside the planet is a result of the self gravity, but also of the body’s deformation in response to any perturbing potential V_{p}. A very convenient way to define this deformation is through the Love number approach (e.g., Love 1911) in which the radial displacement ΔR is proportional to the equipotential perturbing surface (7)where g = Gm/R^{2} is the surface gravity, G is the gravitational constant, m is the mass of the planet, and h_{f} is the fluid second Love number for radial displacement. For a homogeneous sphere h_{f} = 5/2, but more generally, it can can be obtained from the DarwinRadau equation (e.g., Jeffreys 1976) (8)where ℐ is the mean moment of inertia, which depends on the internal mass differentiation.
Like the main satellites in the solar system, closein planets deform under the action of the centrifugal and tidal potentials. The first results from the planet’s rotation rate Ω about the c axis, while the second results from the differential attraction of a mass element by the nearby star with mass m_{∗}. For simplicity, we consider that the planet reached the final tidal equilibrium; that is, its orbit is circular with radius r_{0}, the spin axis is normal to the orbital plane (zero obliquity), and the rotation rate is synchronous with the orbital mean motion n, always pointing the long axis a to the star (e.g., FerrazMello et al. 2008). Thus, on the planet’s surface, the nonspherical contribution of the perturbing potential is given by (e.g., Correia & Rodríguez 2013) (9)Replacing the above perturbing potential in expression (7) and comparing with Eq. (6) for the surface of the ellipsoid, it becomes straightforward that (10)and (11)with (12)where . Since , the closer the planet is to the star, the greater is the difference between the ellipsoid semiaxes. However, there is a maximum value for q, corresponding to the Roche limit, r_{R} (Chandrasekhar 1987) (13)Adopting the maximum value possible for h_{f} = 5/2 (Eq. (8)) gives q_{max} ≈ 1/12 ≈ 0.08, a_{max} ≈ 1.25 b, and c_{max} ≈ 0.92 b. Thus, the approximation used to obtain expression (6) is still valid.
3. Transit light curve
Fig. 1 Difference between the transit light curves of an ellipsoidal and a spherical planet, for some real closein planets with different inclinations and distances to the star (Table 1). The red line shows the difference obtained with our model (Eq. (19)), while the blue line is obtained for an oblate planet with a = b ≠ c, both with stellar uniform flux. The dashed line also corresponds to the difference obtained with our model, but using a quadratic stellar limbdarkening correction I(μ) /I(1) = 1−0.4830(1 − μ) − 0.2023(1 − μ)^{2}, which corresponds to a stellar temperature of 5500 K (Sing 2010). 
We let (x,y,z) be the standard Cartesian coordinate system for exoplanets, centered on the star, where (x,y) are in the plane of the sky and z is along the line of sight. In this frame, a generic point on the planet’s surface r = (x,y,z) can be obtained as (14)where r_{0} = (x_{0},y_{0},z_{0}) is the position of the center of mass of the planet with respect to the star, and is a rotation matrix. When the planet orbits the star in a synchronous circular orbit with zero obliquity, we have (15)where θ is the rotation angle, i is the inclination angle between the normal to the orbit and the line of sight, and and are the standard 3 × 3 rotation matrix about the x and z axis, respectively (e.g., Murray & Correia 2010). In this particular tidal equilibrium configuration, for a given time t, we have θ = n(t − t_{0}), and thus . Moreover, the position of the planet in its orbit can also be easily obtained from as (16)From expression (14) we have (17)Replacing (17) in Eq. (2), we obtain the equation for the ellipsoidal surface of the planet in the new coordinate system (18)The projection of the ellipsoid in the (x,y) plane of the sky is obtained simply by setting z = z_{0} in previous equation. The result is an ellipse with general equation given by (19)with In the new coordinate system, the diminishing of the stellar flux is then given by the overlap between the planetary ellipse (Eq. (19)) and the stellar disk, defined by (23)where we assumed a spherical shape for the star with radius R_{∗}. There are many different methods of computing the overlap area. If the stellar flux is uniform, one easy way is to compute the segment areas (e.g., Eberly 2008; Hughes & Chraibi 2012). In more realistic problems (models that include limbdarkening, stellar rotation, stellar activity, etc.), one can use Monte Carlo integrations (e.g., Press et al. 1992; Carter & Winn 2010).
In Fig. 1 we plot the flux difference between the transit light curve of an ellipsoidal and a spherical planet for some real closein planets, chosen with different inclinations and distances to the star (Table 1). The projected ellipsoid is obtained with our model (Eq. (19)), while the spherical radius is chosen such that its circumference area is equal to the projected ellipse area just after the interior ingress (or just before the interior egress). We also show the flux difference for an oblate planet with a = b ≠ c, in order to compare with previous studies (Seager & Hui 2002; Barnes & Fortney 2003; Carter & Winn 2010).
There are essentially two main differences between the ellipsoidal and the spherical cases. One occurs at the ingress and egress phases, which corresponds to an oscillation in the flux difference. This feature is mainly due to the polar oblateness, and therefore it was already identified in previous studies (e.g., Seager & Hui 2002). The second feature, previously unnoticed, corresponds to a “bump” increase in the light curve difference during the whole transit, which is due to the rotation of the planet. For instance, for i = 90°, the semiaxes of the projected ellipse are given by A^{−1/2} = (cos^{2}θ/a^{2} + sin^{2}θ/b^{2})^{−1/2} and C^{− 1/2} = c (Eqs. (20)–(22)); that is, when a ≠ b, the long semiaxis A^{−1/2} is a function of the phase θ. Thus, the larger is the phase span during the transit, the higher the difference with respect to the spherical case. The maximum value is reached at the center of the transit (θ = 90°).
The effect of an ellipsoidal planet in the transit light curve is also maximized for edgeon orbits (i = 90°). For lower (or higher) inclinations, the planet spends less time in front of the star, so the “bump” increase due to the rotation becomes smaller. Nevertheless, for almost grazing orbits we can still observe significant differences in the light curve, although the two main features described above are no longer completely individualized (see Fig. A.1, for a Jupiterlike planet at the Roche limit).
4. Inner structure determination
Observational data (left side) and derived parameters (right side) for some closein planets near the Roche limit (r_{0}< 3 r_{R}).
For those systems where the differences in the light curve can be detected, one can use an ellipsoid instead of a sphere to fit the photometric observations (Eq. (19)). In addition to the equatorial radius b, the projected ellipsoid only depends on q (Eqs. (10)–(12)). Then, when we adjust our model to the observational data, q is the only supplementary parameter to fit, which accounts for all the observed differences in the transit light curve. In the expression of q (Eq. (12)), all parameters but h_{f} are also known from the observational data. It is thus possible to obtain an observational estimation for h_{f}, which represents an important additional constraint for the inner structure differentiation (Eq. (8)).
Seager & Hui (2002) originally showed that the oblateness of transiting exoplanets can be constrained from the variations in the light curve during the ingress and egress phases. However, Barnes & Fortney (2003) noticed that these oscillations are short in time and therefore difficult to observe, although some relatively weak constraints on the Love number were derived for the system HD 189733 (Carter & Winn 2010).
In Fig. 1 we observe that the “bump” increase due to the rotation of the prolate planet is present during the whole transit, not only at ingress and egress phases. As a consequence, the flux differences can be significantly larger than previously thought, increasing our chances of observing the shape of closein planets. For giant planets near the Roche limit, the flux differences can reach 10^{4}, which is above the intrinsic stellar photometric variability on transit timescales, estimated to be ~10^{5} (Borucki et al. 1997). For longer distances to the star, the flux difference can drop below the 10^{5} threshold, so the shape is more difficult to determine. For rocky planets, the flux difference is also smaller than 10^{5}. However, since the Roche limit is closer to the star for these planets, the “bump” becomes more prominent, increasing again the chances of detection.
In Fig. 1 we also observe that the flux difference is about the same for a stellar uniform flux, or when limbdarkening is considered. The reason is that the differences in the shape of the planet are so tiny with respect to the size of the star that locally the flux can be considered approximately uniform. Therefore, for realistic light curves obtained with nonuniform stellar flux, one can still obtain the correction introduced by the axial asymmetry using a uniform flux model, which is much easier to compute.
5. Density determination
Transiting systems provide a radius R_{s} for the planet, hence a bulk density, when coupled with the planetary mass determined from radialvelocity measurements. Most studies assume a spherical planet in order to work out the volume , thus providing a spherical bulk density ρ_{s} = m/V_{s}. However, for closein planets, the true volume of the ellipsoid is given by V = 4πabc/ 3, which gives for the bulk density (Eqs. (10), (11)) (24)When there is enough precision in the data, b and q can be directly adjusted from the transit light curve, instead of R_{s} (Sect. 4). For lowmass planets such as Kepler78b, it may be difficult to spot the tiny differences in the light curve with respect to a spherical planet (Fig. 1). In those cases, one can use the spherical model, where R_{s} is related to the projected ellipsoid (Eq. (19)) through (25)At the center of the transit (θ = 90°), we thus have (Eqs. (20)–(22)) (26)which gives for the bulk density (Eq. (24)) (27)Here, q cannot be adjusted from the observations, but it can be estimated from expression (12) using b ≈ R_{s}, and adopting the solar system data, h_{f} = 1.5 for gaseous planets, and h_{f} = 2 for rocky planets (Yoder 1995).
In Table 1 we compute the relative change in the bulk density, Δρ = 1 − ρ/ρ_{s}, for some closein planets near the Roche limit (r_{0}< 3 r_{R}). These corrections agree with those obtained numerically by Burton et al. (2014). The small differences observed between the two methods result from the assumption of different inner structure models. Indeed, an exact match could be obtained using slightly different h_{f} values in our model.
In Fig. 2 we show the relative change in the bulk density as a function of the distance to the star r_{0}. The largest possible correction in the density is ρ_{max} ≈ 0.7 ρ_{s}, obtained from q_{max} (Eq. (13)) with i = 90°, which is 30% smaller than the spherical value. However, for realistic cases of planets near the Roche limit (r_{0} ≈ r_{R}), we get values lower than 20% for rocky planets and lower than 15% for gaseous planets. As in Burton et al. (2014), we also observe that the relative change in the density decreases rapidly with the distance to the star. This correction is less than 5% for planets with r_{0}> 1.5 r_{R}, which is equivalent to the measurement error on the currently published bulk densities (e.g., Hellier et al. 2011). Therefore, with the current observational precision, a correction in the density determination is only justified for planets that are close to the Roche limit.
Fig. 2 Relative change in the bulk density of closein planets, Δρ = 1 − ρ/ρ_{s}, as a function of the distance to the star r_{0}. The red curve is relative to gaseous planets (h_{f} = 1.5), while the blue one is relative to rocky planets (h_{f} = 2.0). Both curves are obtained using expression (27) with i = 80°. The dots represent the planets listed in Table 1. For r_{0}> 1.5 r_{R}, the change is negligible when compared to the present observational errors in the density determination. 
6. Conclusions
We have presented here a simple model for correcting the transit light curve of closein planets, which is easy to implement and which allows us to obtain important constraints for the internal structure of these planets. On one hand, it provides a better determination for the bulk density (Eq. (27)), which is lower than previously assumed (Fig. 2). For instance, it has been announced that Kepler78b has an Earthlike density (Howard et al. 2013; Pepe et al. 2013), but actually its real density is probably lower than 5 g/cm^{3} (Table 1). On the other hand, information on the internal structure differentiation can be obtained by determining the fluid Love number h_{f} (Eq. (8)), even when the correction in the density is negligible. If there is enough precision in the transit light curve (<10^{4}), one can adjust an ellipsoid instead of a sphere to the data (Fig. 1). Then, in addition to the equatorial radius b, one can also quantify the planet asymmetry q, and thus obtain an observational estimation for h_{f} (Eqs. (10)–(12)).
There are other methods that allow indirect determination of the Love number (e.g., Batygin et al. 2009; Ragozzine & Wolf 2009; Mardling 2010), but they require that the global dynamics of the system is known (the presence of planetary companions, precise eccentricities and inclinations, etc.), which is very unlikely to achieve at present. Our method does not require any additional knowledge on the system, apart from the variations in the transit light curve. It thus provides a direct determination of the Love number, so is therefore much more reliable.
In our model we have assumed that the planet is in a tidal final equilibrium configuration. However, our method can be generalized to planets in eccentric orbits, with any rotation rate and nonzero obliquity. For that purpose, a point at the surface of the planet (Eq. (14)) is given by a new rotation matrix (28)where ε is the obliquity (the angle between the equator and the orbital plane), and ϕ is the precession angle. Also, we have now
θ = Ω(t − t_{0}), since Ω ≠ n. Therefore, the position of the planet on its orbit (Eq. (16)) also needs to be corrected through (29)where ω is the argument of the pericentre and v is the true anomaly. For elliptical orbits with semimajor axis a_{0} and eccentricity e_{0}, we also have .
Finally, our model can also be used to correct the light curve of eclipsing close binary stars, where the secondary replaces the planet. For the primary, instead of using Eq. (23), we need to use an analog to Eq. (19), where a, b and c are replaced by the primary semiaxes and r_{0} = 0. However, in this case the approximation for the disturbing potential (Eq. (9)) may require additional corrective terms (e.g., Budaj 2011; Burton et al. 2014).
Acknowledgments
The author acknowledges support from Fundação para a Ciência e a Tecnologia, Portugal (PEstC/CTM/LA0025/2011), and from the PICS05998 FrancePortugal program.
References
 Adams, E. R., LópezMorales, M., Elliot, J. L., et al. 2011, ApJ, 741, 102 [NASA ADS] [CrossRef] [Google Scholar]
 Bakos, G. Á., Hartman, J., Torres, G., et al. 2011, ApJ, 742, 116 [NASA ADS] [CrossRef] [Google Scholar]
 Barge, P., Baglin, A., Auvergne, M., et al. 2008, A&A, 482, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barnes, J. W., & Fortney, J. J. 2003, ApJ, 588, 545 [NASA ADS] [CrossRef] [Google Scholar]
 Batygin, K., Bodenheimer, P., & Laughlin, G. 2009, ApJ, 704, L49 [NASA ADS] [CrossRef] [Google Scholar]
 Borucki, W. J., Koch, D. G., Dunham, E. W., & Jenkins, J. M. 1997, in Planets Beyond the Solar System and the Next Generation of Space Missions, ed. D. Soderblom, ASP Conf. Ser., 119, 153 [Google Scholar]
 Budaj, J. 2011, AJ, 141, 59 [NASA ADS] [CrossRef] [Google Scholar]
 Burton, J. R., Watson, C. A., Fitzsimmons, A., et al. 2014, ApJ, 789, 113 [NASA ADS] [CrossRef] [Google Scholar]
 Carter, J. A., & Winn, J. N. 2010, ApJ, 709, 1219 [NASA ADS] [CrossRef] [Google Scholar]
 Chan, T., Ingemyr, M., Winn, J. N., et al. 2011, AJ, 141, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S. 1987, Ellipsoidal figures of equilibrium (New York: Dover) [Google Scholar]
 Correia, A. C. M. 2009, ApJ, 704, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Correia, A. C. M., & Rodríguez, A. 2013, ApJ, 767, 128 [NASA ADS] [CrossRef] [Google Scholar]
 Dumusque, X., Bonomo, A. S., Haywood, R. D., et al. 2014, ApJ, 789, 154 [NASA ADS] [CrossRef] [Google Scholar]
 Eberly, D. 2008, http://www.geometrictools.com/ [Google Scholar]
 Enoch, B., Anderson, D. R., Barros, S. C. C., et al. 2011, AJ, 142, 86 [NASA ADS] [CrossRef] [Google Scholar]
 FerrazMello, S., Rodríguez, A., & Hussmann, H. 2008, Celest. Mech. Dyn. Astron., 101, 171 [Google Scholar]
 Gillon, M., Smalley, B., Hebb, L., et al. 2009, A&A, 496, 259 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gillon, M., Demory, B.O., Benneke, B., et al. 2012a, A&A, 539, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gillon, M., Triaud, A. H. M. J., Fortney, J. J., et al. 2012b, A&A, 542, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gillon, M., Anderson, D. R., CollierCameron, A., et al. 2014, A&A, 562, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hatzes, A. P., Fridlund, M., Nachmani, G., et al. 2011, ApJ, 743, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Hébrard, G., Collier Cameron, A., Brown, D. J. A., et al. 2013, A&A, 549, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hellier, C., Anderson, D. R., CollierCameron, A., et al. 2011, ApJ, 730, L31 [NASA ADS] [CrossRef] [Google Scholar]
 Howard, A. W., SanchisOjeda, R., Marcy, G. W., et al. 2013, Nature, 503, 381 [NASA ADS] [CrossRef] [Google Scholar]
 Hughes, G. B., & Chraibi, M. 2012, Computing and Visualization in Science, 15, 291 [CrossRef] [Google Scholar]
 Hut, P. 1980, A&A, 92, 167 [NASA ADS] [Google Scholar]
 Jeffreys, H. 1976, The earth. Its origin, history and physical constitution (CUP) [Google Scholar]
 Leconte, J., Lai, D., & Chabrier, G. 2011, A&A, 528, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Love, A. E. H. 1911, Some Problems of Geodynamics (CUP) [Google Scholar]
 Mardling, R. A. 2010, MNRAS, 407, 1048 [NASA ADS] [CrossRef] [Google Scholar]
 Murray, C. D., & Correia, A. C. M. 2010, in Exoplanets (University of Arizona Press), 15 [Google Scholar]
 Pepe, F., Cameron, A. C., Latham, D. W., et al. 2013, Nature, 503, 377 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN. The art of scientific computing, 2nd edn. (Cambridge: University Press) [Google Scholar]
 Ragozzine, D., & Wolf, A. S. 2009, ApJ, 698, 1778 [NASA ADS] [CrossRef] [Google Scholar]
 Seager, S., & Hui, L. 2002, ApJ, 574, 1004 [NASA ADS] [CrossRef] [Google Scholar]
 Sing, D. K. 2010, A&A, 510, A21 [Google Scholar]
 Smalley, B., Anderson, D. R., CollierCameron, A., et al. 2012, A&A, 547, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Southworth, J., Hinse, T. C., Dominik, M., et al. 2009, ApJ, 707, 167 [NASA ADS] [CrossRef] [Google Scholar]
 Valsecchi, F., & Rasio, F. A. 2014, ApJ, 787, L9 [NASA ADS] [CrossRef] [Google Scholar]
 Yoder, C. F. 1995, in Global Earth Physics: A Handbook of Physical Constants (Washington D.C: American Geophysical Union), 1 [Google Scholar]
Online material
Appendix A: Additional figure
Fig. A.1 Transit light curves for a Jupiterlike planet orbiting at the Roche limit (r_{0} = r_{R} = 2.51 R_{⊙}) of a Sunlike star with uniform flux (in red). For comparison, we also show the light curve of a spherical planet (dashed line), and the difference between the two curves (bottom figure). The dots correspond to the relative positions shown in the top figure. The shape of the planet in the top figure is obtained with an exaggerated h_{f} = 6, so that we can better spot the modifications in the projected shape, but the transit light curves are obtained with the realistic Jupiter’s value h_{f} = 1.5 (Yoder 1995). 
All Tables
Observational data (left side) and derived parameters (right side) for some closein planets near the Roche limit (r_{0}< 3 r_{R}).
All Figures
Fig. 1 Difference between the transit light curves of an ellipsoidal and a spherical planet, for some real closein planets with different inclinations and distances to the star (Table 1). The red line shows the difference obtained with our model (Eq. (19)), while the blue line is obtained for an oblate planet with a = b ≠ c, both with stellar uniform flux. The dashed line also corresponds to the difference obtained with our model, but using a quadratic stellar limbdarkening correction I(μ) /I(1) = 1−0.4830(1 − μ) − 0.2023(1 − μ)^{2}, which corresponds to a stellar temperature of 5500 K (Sing 2010). 

In the text 
Fig. 2 Relative change in the bulk density of closein planets, Δρ = 1 − ρ/ρ_{s}, as a function of the distance to the star r_{0}. The red curve is relative to gaseous planets (h_{f} = 1.5), while the blue one is relative to rocky planets (h_{f} = 2.0). Both curves are obtained using expression (27) with i = 80°. The dots represent the planets listed in Table 1. For r_{0}> 1.5 r_{R}, the change is negligible when compared to the present observational errors in the density determination. 

In the text 
Fig. A.1 Transit light curves for a Jupiterlike planet orbiting at the Roche limit (r_{0} = r_{R} = 2.51 R_{⊙}) of a Sunlike star with uniform flux (in red). For comparison, we also show the light curve of a spherical planet (dashed line), and the difference between the two curves (bottom figure). The dots correspond to the relative positions shown in the top figure. The shape of the planet in the top figure is obtained with an exaggerated h_{f} = 6, so that we can better spot the modifications in the projected shape, but the transit light curves are obtained with the realistic Jupiter’s value h_{f} = 1.5 (Yoder 1995). 

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.