Issue |
A&A
Volume 570, October 2014
|
|
---|---|---|
Article Number | L5 | |
Number of page(s) | 6 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201424733 | |
Published online | 16 October 2014 |
Transit light curve and inner structure of close-in planets⋆
1
Departamento de Física, I3NUniversidade de Aveiro
Campus de Santiago
3810-193
Aveiro
Portugal
e-mail:
correia@ua.pt
2
ASD, IMCCE-CNRS UMR8028, Observatoire de Paris,
UPMC, 77 Av.
Denfert-Rochereau, 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 close-in 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 days1. These close-in 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., Ferraz-Mello 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 (fc) 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 (fa), 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 close-in planets at any point of its orbit, which is based in the equilibrium surface given by second-order 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 semi-principal 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
Y2, 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 Vp. 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/R2 is
the surface gravity, G is the gravitational constant, m is the mass of the planet,
and hf is the fluid second Love number for
radial displacement. For a homogeneous sphere hf = 5/2, but more
generally, it can can be obtained from the Darwin-Radau 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, close-in 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
r0, 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., Ferraz-Mello et al. 2008). Thus, on the planet’s surface, the non-spherical
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
semi-axes. However, there is a maximum value for q, corresponding to the Roche limit, rR (Chandrasekhar 1987)
(13)Adopting
the maximum value possible for hf = 5/2 (Eq. (8)) gives qmax ≈ 1/12 ≈
0.08, amax ≈ 1.25 b, and
cmax ≈ 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 close-in 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 limb-darkening 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
r0 =
(x0,y0,z0)
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 −
t0), 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 =
z0 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 limb-darkening, 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 close-in 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 semi-axes of the projected ellipse are given by A−1/2 = (cos2θ/a2 + sin2θ/b2)−1/2 and C− 1/2 = c (Eqs. (20)–(22)); that is, when a ≠ b, the long semi-axis 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 edge-on 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 Jupiter-like planet at the Roche limit).
4. Inner structure determination
Observational data (left side) and derived parameters (right side) for some close-in planets near the Roche limit (r0< 3 rR).
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 hf are also known from the observational data. It is thus possible to obtain an observational estimation for hf, 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 close-in 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 limb-darkening 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 non-uniform 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 Rs for the planet, hence a bulk density,
when coupled with the planetary mass determined from radial-velocity measurements. Most
studies assume a spherical planet in order to work out the volume
, thus
providing a spherical bulk density ρs =
m/Vs.
However, for close-in 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 Rs (Sect. 4). For low-mass planets such as Kepler-78b, 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 Rs 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 ≈
Rs, and adopting the solar system data,
hf =
1.5 for gaseous planets, and hf = 2 for rocky planets (Yoder 1995).
In Table 1 we compute the relative change in the bulk density, Δρ = 1 − ρ/ρs, for some close-in planets near the Roche limit (r0< 3 rR). 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 hf 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 r0. The largest possible correction in the density is ρmax ≈ 0.7 ρs, obtained from qmax (Eq. (13)) with i = 90°, which is 30% smaller than the spherical value. However, for realistic cases of planets near the Roche limit (r0 ≈ rR), 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 r0> 1.5 rR, 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 close-in planets, Δρ = 1 − ρ/ρs, as a function of the distance to the star r0. The red curve is relative to gaseous planets (hf = 1.5), while the blue one is relative to rocky planets (hf = 2.0). Both curves are obtained using expression (27) with i = 80°. The dots represent the planets listed in Table 1. For r0> 1.5 rR, 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 close-in 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 Kepler-78b has an Earth-like density (Howard et al. 2013; Pepe et al. 2013), but actually its real density is probably lower than 5 g/cm3 (Table 1). On the other hand, information on the internal structure differentiation can be obtained by determining the fluid Love number hf (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 hf (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 non-zero 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 −
t0), 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
semi-major axis a0 and eccentricity e0, 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 semi-axes and r0 = 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 (PEst-C/CTM/LA0025/2011), and from the PICS05998 France-Portugal program.
References
- Adams, E. R., López-Morales, 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]
- Ferraz-Mello, 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., Collier-Cameron, 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., Collier-Cameron, A., et al. 2011, ApJ, 730, L31 [NASA ADS] [CrossRef] [Google Scholar]
- Howard, A. W., Sanchis-Ojeda, 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., Collier-Cameron, 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 Jupiter-like planet orbiting at the Roche limit (r0 = rR = 2.51 R⊙) of a Sun-like 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 hf = 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 hf = 1.5 (Yoder 1995). |
All Tables
Observational data (left side) and derived parameters (right side) for some close-in planets near the Roche limit (r0< 3 rR).
All Figures
![]() |
Fig. 1 Difference between the transit light curves of an ellipsoidal and a spherical planet, for some real close-in 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 limb-darkening 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 close-in planets, Δρ = 1 − ρ/ρs, as a function of the distance to the star r0. The red curve is relative to gaseous planets (hf = 1.5), while the blue one is relative to rocky planets (hf = 2.0). Both curves are obtained using expression (27) with i = 80°. The dots represent the planets listed in Table 1. For r0> 1.5 rR, 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 Jupiter-like planet orbiting at the Roche limit (r0 = rR = 2.51 R⊙) of a Sun-like 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 hf = 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 hf = 1.5 (Yoder 1995). |
In the text |
Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.