Issue 
A&A
Volume 555, July 2013



Article Number  A51  
Number of page(s)  15  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201321661  
Published online  01 July 2013 
The effect of rotation and tidal heating on the thermal lightcurves of super Mercuries
^{1} Univ. Bordeaux, LAB, UMR 5804, 33270 Floirac, France
email: selsis@obs.ubordeaux1.fr; hersant@obs.ubordeaux1.fr; asophie.maurin@gmail.com; rayray.sean@gmail.com
^{2} CNRS, LAB, UMR 5804, 33270 Floirac, France
^{3} Laboratoire de Météorologie Dynamique, Institut Pierre Simon Laplace, 75231 Paris Cedex 5, France
email: jeremy.leconte@lmd.jussieu.fr
^{4} UNSCNRS – Observatoire de la Côte d’Azur, Laboratoire Lagrange, BP 4229, 06304 Nice Cedex 4, France
email: marco.delbo@oca.eu
Received: 5 April 2013
Accepted: 16 May 2013
Shortperiod (<50 day), lowmass (< 10 M_{⊕}) exoplanets are abundant, and the few of them whose radius and mass have been measured already reveal a diversity in composition. Some of these exoplanets are found on eccentric orbits and are subjected to strong tides that affect their rotation and result in significant tidal heating. Within this population, some planets are likely to be depleted in volatiles and have no atmosphere. We modeled the thermal emission of these super Mercuries to study the signatures of rotation and tidal dissipation on their infrared lightcurve. We computed the timedependent temperature map on the surface and in the subsurface of the planet and the resulting diskintegrated emission spectrum received by a distant observer for any observation geometry. We calculated the illumination of the planetary surface for any Keplerian orbit and rotation. We included the internal tidal heat flow, vertical heat diffusion in the subsurface and generated synthetic lightcurves. We show that the different rotation periods predicted by tidal models (spinorbit resonances, pseudosynchronization) produce different photometric signatures, which are observable provided that the thermal inertia of the surface is high, as for solid or melted rocks (but not regolith). Tidal dissipation can also directly affect the lightcurves and make the inference of the rotation more difficult or easier depending on the existence of hot spots on the surface. Infrared lightcurve measurement with the James Webb Space Telescope and EChO can be used to infer exoplanets’ rotation periods and dissipation rates and thus to test tidal models. This data will also constrain the nature of the (sub)surface by constraining the thermal inertia.
Key words: techniques: photometric / planetary systems
© ESO, 2013
1. Introduction
Exoplanet Doppler searches have revealed a high frequency of exoplanets in the 1−10 M_{⊕} mass range. Around F, G, K stars, Mayor et al. (2011) have found a debiased occurrence rate of f = 0.41 ± 0.16 for these masses and orbital periods shorter than 50 days. More than one quarter of all FGK stars host this type of hot superEarth, as they are sometimes called. For the same types of stars, orbital periods and the 1 − 1.8 R_{⊕} range (corresponding to 1−10 M_{⊕} for rocky planets), Petigura et al. (2013) find an occurrence f ~ 17% from Kepler candidates, a result that significantly departs from the estimates derived from radial velocity detections. From the HARPS survey of M stars and for the same mass range, Bonfils et al. (2013) inferred a frequency of and for orbital period within 1 − 10 d and 10−100 d, respectively. Statistics of Kepler candidates for cool stars imply a similar frequency: 0.51 planet per star for the 0.5 − 1.4 R_{⊕} range and periods shorter than 50 days (Dressing & Charbonneau 2013).
Despite the abundance of Kepler transiting candidates with radii below 2 R_{⊕}, only a few planets have wellmeasured masses below 10 M_{⊕}. These lowmass planets exhibit a variety of densities that goes from lower than 1 g/cm^{3} for Kepler 11 c, e, f (Lissauer et al. 2013) up to 8.8 ± 2.9 g/cm^{3} for Kepler 10 b (Batalha et al. 2011). Within this mass range, we can expect to find volatilerich planets with dense atmospheres but also airless bodies that have either accreted from dry, refractory material or lost their volatile component due to the strong irradiation. At 0.04 AU from a Sunlike star, a 5 M_{⊕} planet that initially contained 10% of its mass in H_{2}O (a proportion 100 times higher than that of Earth) could lose its whole volatile content in 500 Myr (Selsis et al. 2007a).
Some known lowmass exoplanets on eccentric orbits.
Fig. 1
Number of rotations per orbital period for a pseudosynchronized planet as a function of eccentricity. The filled circle shows the pseudosynchronization (ω_{eq} / n = 2.38) for the default value of eccentricity used in this study (e = 0.45). 
Although most of shortperiod lowmass exoplanets have circular orbits, some have significant eccentricities. Their orbit could have been excited by other planets in the system, continuously or at some moment in the past, or circularization might be slow due to a low dissipation. Table 1 gives a list of known exoplanets found by radial velocity surveys with 4.5 < msini < 12 M_{⊕} and 0.1 < e < 0.45.
For planets on eccentric orbits and subjected to strong tides the rotation period is not known a priori. The planet can be captured into a spinorbit resonance like Mercury, which rotates three times every two orbits (resonance 3:2). It can also have reached pseudosynchronization, a theoretical state that minimizes the dissipation (Hut 1981; Levrard et al. 2007). This equilibrium rotation (ω_{equ} / n, see formula (9) and Fig. 1) is always faster than synchronization (1:1 resonance). The existence of the pseudosynchronous state has recently been put in question for the case of solid planets by Makarov & Efroimsky (2013), who claim that only spinorbit resonances should be expected. The planet can, however, be captured into different spinorbit resonances depending on its initial rotation, the evolution of its eccentricity, its internal structure and rheology, its oblateness (which depends on its rotation rate) and interactions with other planets. In general, the order of the most likely spinorbit resonances is expected to increase with the eccentricity (Makarov et al. 2012). Measuring the rotation period of eccentric exoplanets would therefore provide valuable information for tidal theories as well as useful constraints on the nature and history of the planets. In the present study, we explore how rotation can be measured or constrained thanks to the infrared photometric variations of the planet.
In addition to transits and eclipses when they occur, the photometric variation of a spatiallyunresolved planetary system contains two superimposed signals: the intrinsic variability of the star, and the change in apparent brightness of the planet(s). One obvious planetary modulation is associated with the orbital phase. This phase modulation (or phase curve), which is obvious in reflected light, also exists at thermal wavelengths provided that a sufficient brightness temperature contrast exists between the day and night sides at the wavelength of the observation. Phase curves have been observed for hot Jupiters both in the infrared (Knutson et al. 2007) and in visible light (Welsh et al. 2010), for both transiting and nontransiting planets (Crossfield et al. 2010), and for eccentric planets (Lewis et al. 2013; Cubillos et al. 2013). The reflected lightcurve of Kepler 10 b, a very shortperiod transiting “terrestrial” planet, was obtained by Batalha et al. (2011) despite a relative amplitude of only ~5 × 10^{6} (5 ppm), and the dayside emission of a similar exoplanet, 55 Cnc e was measured with Spitzer (Demory et al. 2012).
Thermal and reflected lightcurves have been modeled for multiplanet systems (Kane et al. 2011; Kane & Gelino 2013). Atmosphere models have been used to compute multiband phase curves, assuming radiative equilibrium (Barman et al. 2005), a timedependent radiativeconvective model (Iro et al. 2005) and more recently 3D general circulation models (Showman et al. 2009; Selsis et al. 2011; Menou 2012). The photometric signature of eccentric exoplanets has been studied by several authors (Iro & Deming 2010; Kane & Gelino 2011a; Cowan & Agol 2011; Lewis et al. 2013). Maurin et al. (2012) modeled the spectral phase curves of synchronously rotating super Mercuries and showed that the inclination, albedo and radius of the planet could be inferred for outoftransit configurations by multiwavelength observations. Lightcurves of multiplanet systems have been modeled by Kane et al. (2011) and Kane & Gelino (2013). Bolmont et al. (2013) investigated the influence of tidal heating on the thermal emission and phase curve of 55 Cnc e.
Thermal lightcurves have been proposed as a technique to characterize habitable Earthlike planets and to detect the presence of a dense atmosphere with future instruments (Selsis 2004). Several aspects of characterizing habitable worlds using orbital photometry have been studied: the effect of the planet’s obliquity (Gaidos & Williams 2004) and eccentricity (Cowan et al. 2012), the presence of a moon (Selsis 2004; Moskovitz et al. 2009), and the rotation period in the visible (Ford et al. 2001; Pallé et al. 2008) and infrared (GómezLeal et al. 2012).
Observing thermal phase curves of planets smaller than 2 R_{⊕} will require future space observatories like James Webb Space Telescope (JWST) and EChO (Tinetti et al. 2012). Today, it may be achievable with Spitzer only for a very hot and very nearby object like 55 Cnc e. EChO aims at performing high precision spectrophotometry with an ability to detect relative photometric variations better than 10 ppm with days to weeks stability and a broad spectral coverage (0.4 − 16 μm). Observing phase curves to characterize lowmass exoplanets and their atmosphere is an important part of the science of EChO.
Here we study the photometric signature of SuperMercuries, big rocky exoplanets with no atmospheres on eccentric orbits. In Sect. 2, we describe the insolation, thermophysical and tidal models used to produce thermal lightcurves. In Sect. 3, we present the results and describe the influence of rotation, eccentricity, observing geometry, thermal inertia and tidal dissipation on the lightcurves. Our approach and results are discussed in Sect. 4 and the main conclusions of the study are given in Sect. 5.
2. Model
2.1. Illumination of the planetary surface
Our first step is to compute the timedependent orbital distance r and zenith angle θ at any point of the planetary surface. We assume a Keplerian orbit and a rotation defined by a period, obliquity and argument of periastron. Although our model can be applied to any orientation of the rotation axis, we do not consider oblique planets in this study and the (prograde) rotation is thus only defined by its period.
The bolometric stellar flux φ_{⋆} absorbed by the planetary surface at a given longitude and latitude and at an orbital distance r is given by (1)where σ is the StefanBoltzmann constant, T_{eff} and R_{⋆} are the effective temperature and radius of the star, A is the surface bolometric albedo (assumed to be constant over the planet and thus equal to the Bond albedo), μ the cosine of the zenith angle θ and ϵ the bolometric emissivity. Both r and μ are timedependent. The distance r is obtained by solving Kepler’s equation by NewtonRaphson iterations using the cubic approximation by Mikkola (1987) as a first guess (IDL keplereq routine). For a null obliquity the expression for μ is given by (2)where lon_{⋆}(t) is the longitude of the substellar point, which can be calculated as follows for a null obliquity and a prograde rotation: (3)where ν is the true anomaly of the planet, ω the angular spin velocity and C is an arbitrary constant depending on the position of the planet at t = 0 and the origin of the longitude coordinate. If ν(t = 0) and C are set to 0° then the planet is initially at periastron with the substellar point at 0° longitude.
2.2. Internal heat flow from tidal dissipation
In this study, we neglect the internal heat flux due to the decay of radiogenic species or the release of the energy of accretion/differentiation, which should only affect the planetary thermal emission at a measurable level for very young objects. We include the surface heat flow resulting from tides given by where Ė_{tides} is the rate of tidal dissipation. To calculate the dissipation rate, we use the constant timelag model described in Leconte et al. (2010). For a planet with a mass M_{p} orbiting a star with mass M_{⋆} at a semimajor axis a and with a null obliquity, the rate of tidal dissipation is given by (4)with where e is the eccentricity, n is the orbital mean motion, ω is the rotation rate of the planet. The quantity K_{p} is defined as (8)where G is the gravitational constant, R_{p} and M_{p} are the radius and mass of the planet, and M_{⋆} is the mass of the star. Δt_{p} and k_{2,p} are the timelag and potential Love number of degree 2 for the planet, for which we use the Earth quantities (Δt_{p} = 630 s and k_{2,p} = 0.3) as the default dissipation (Neron de Surgy & Laskar 1997).
In this study, we assume that eccentric planets are either captured into a spinorbit resonance or have reached the equilibrium rotation rate (or pseudosynchronization) that minimizes Ė_{tides}. This equilibrium rotation rate is a function of e (see Fig. 1): (9)Figure 2 shows the eccentricityorbital period domain in which three observation criteria are fulfilled. First, tidal dissipation increases the planet thermal emission by at least 10%. Second, the flux from the planet exceeds 10^{5} times that of the star at 10 μm. Third, the eccentricity damping time e / ė is longer than 1 Gyr. This damping timescale is calculated according to Leconte et al. (2010, Eq. (6)) assuming no forcing from planetary companions but perturbations from massive or eccentric planets in the same system could maintain a high eccentricity for longer Beust et al. (2012); Bolmont et al. (2013). This criterion is therefore conservative and the observability domain could be extended to higher eccentricities in multiple systems. Figure 2 also shows how this observability domain varies with the mass of the central star, the radius of the planet and the dissipation factor (which is poorly constrained for exoplanets). In these graphs the thermal emission of the planet is calculated assuming an isothermal surface heated by the fully redistributed stellar illumination (averaged over one orbital period) and the tidallyproduced internal heat flow.
Fig. 2
Observability of the 10 μm emission excess due to tidal heating in a periodeccentricity diagram. a) Dashed contours indicate the circularization timescale (e / ė), solid contours give the planet/star contrast ratio in ppm, and the dotted contours show the ratio between the emission due to tidal heating and the emission without tidal heating. The contours are calculated for a 0.2 M_{⊙} star, A = 0.2, an isothermal surface (stellar and internal flux fully redistributed), a 2 R_{⊕} planet and the dissipation parameters of the Earth. The colored area represents the observability domain for the tidal heating, where the configuration lasts for more than 1 Gyr, produce a planet/star contrast higher than 10 ppm and results in tidal heating contributing to more than 10% of the emission. Panels b), c) and d) respectively show the influence of the stellar mass, the planetary radius and the dissipation factor on this tidal heating observability domain. 
2.3. Subsurface heat diffusion
We compute the surface and subsurface timedependent temperature taking into account the internal heat flow and the forcing by the absorbed stellar light by modeling the vertical diffusion of heat.
The 1D heat equation can be written: (10)with z the depth, ρ the density of the layer, κ the thermal conductivity and c the heat capacity at constant volume (usually noted c_{V}). Maurin (2012) investigated the effect of using a realistic temperaturedependency for c(T) and κ(T) compared with using constant values for c and κ calculated for the planet’s equilibrium temperature T_{eq}. Maurin (2012) conducted a broad exploration of the parameter space of T_{eq}, rotation rate, and surface materials, and found maximum differences of only a few K locally in the surface temperature with a negligible impact on the diskintegrated thermal emission. Therefore, and although our model can include the temperature dependency of κ and c, we assume κ and c to be constant. As in most heat diffusion models used in planetary science (Spencer et al. 1989), Eq. (10) simplifies to: (11)
2.3.1. Boundary conditions
At a given longitudelatitude point of the planet, the upperboundary conditions is determined by the absorbed stellar light and the thermal emission is given by (12)where z_{max} corresponds to the surface and ϵ is the surface emissivity (assumed to be equal to 1 in this study). The lowerboundary condition is given by the internal heat flow: (13)where z_{min} is the bottom of the modeled vertical column and Φ_{i} is computed from the tidal dissipation. We assume that z_{min} is large enough that the variations of surface irradiation do not propagate down to z_{min} such that temperature eventually reaches a steadystate.
2.3.2. Dimensionless equations
We introduce a timescale τ (without initially specifying it) and the typical length over which diffusion acts in a duration τ(14)These quantities can then be used to make the change of variables (15)In the dimensionless variables, Eq. (11) takes the simple form: (16)while Eqs. (12) and (13) become Once the illumination flux φ_{⋆}(t) and the internal flux Φ_{i} are known, the surface temperature T(z,t) is determined by a single physical parameter, the thermal inertia Γ: (19)
2.3.3. Numerical scheme
To calculate the temperature structure we discretize Eq. (10) assuming layers of constant thickness. To ensure conservation, the temperatures are defined at the center of each layer and fluxes are calculated at each interface. The system is then solved using a CrankNicolson scheme, with Neumann boundary conditions computed from Eqs. (12) and (13). One version of our model uses a variable thickness (small near the surface, larger close to the lower boundary) without significantly reducing the required number of layers.
In order to set the thickness of the layers and the total depth of the column, we need to evaluate the characteristic length scales of the problem. If the surface insolation is periodic with a pulsation ω_{⋆}, we can use τ = 1 / ω_{⋆} and L = l_{s}, where l_{s} is the thermal skin depth as defined by Spencer et al. (1989): (20)For an eccentric orbit, the insolation at a given point of the surface is not periodic. The length of the day – or duration between two consecutive star rises – varies throughout the orbit, with the exception of some spinorbit resonances. In addition, the zenith angle does not vary periodically (except, again, for some spinorbit resonances). We consider two timescales: the mean diurnal cycle and the orbital period, and their associated length. The illumination modulation with the shortest periodicity determines the value of the thermal skin and imposes the vertical spatial resolution. If the orbital motion induces the shortest periodicity then the thermal skin l_{s} is calculated assuming ω_{⋆} = n, where n is the mean motion. If the averaged length of the day (calculated assuming a mean orbital motion) is shorter than the orbital period then the thermal skin l_{s} is calculated using ω_{⋆} = ω_{rot} − n, where ω_{rot} is the pulsation associated with the sidereal rotation. The modulation with the longest period sets the depth of the column to be modeled. In practice, however, we use only one single (the smallest) value for the thermal skin l_{s} and typically place our deepest layer at 10 to 20l_{s}, using 64 to 256 layers, depending on the case.
This method has the advantage of simplicity (compared for instance with a timedependent adaptativegrid) and effectiveness. However, one should keep in mind that using a resolution based on an average value for ω_{⋆} and l_{s} may prevent the model from resolving heat waves with very short wavelengths that are produced during a small fraction of the orbit, for instance near periastron for highly eccentric orbits. This is one reason why we systematically verify that the spatial resolution and total depth we use are sufficient, by making sure that the results remain unchanged after an increase of the number of layers and/or the total depth.
Temperature profiles are typically calculated for 684 surface points (19 latitudes, 36 longitudes), enough to produce smooth lightcurves. A smaller number of points is used when testing the sensitivity to initial conditions, the required total depth, number of layers, and number of orbits. As the use of the model is not limited by computation time we can afford to use a small and constant time step (typically 0.2% of the orbital period), which is the same for all points on the planetary surface.
Default parameters of the study (unless other values are specified).
2.3.4. Initial temperature profile
For the result to be independent of the arbitrary initial (sub)surface temperature profile, a minimum number of orbits must be simulated. To keep this number as small as possible it is important to start from initial conditions that are as close as possible to an average thermal state. For each point on the planet, we calculate , the mean stellar energy absorbed over all orbits of the simulation. In the absence of an internal heat flux we start with an isothermal profile with temperature given by . If there is an internal heat flux Φ_{i} then the initial temperature profile is: We note that the real timeaveraged temperature and internal energy differ from this initial profile, even in the simplest case of radiative equilibrium and no internal heat flux. This is a consequence of the RogersHölders inequality:
Because we calculate a temperature from an averaged flux, this initial condition is always too hot but the heat excess is not known a priori. For large variations of T(t) this excess can be significant, especially if T(t) approaches 0 K. For this reason, it is important to check that simulations have been run over a sufficient number of orbits for the modeled subsurface to release this excess heat. The deeper the modeled column and the larger the thermal inertia the longer it takes to transport the initial excess heat to the surface and radiate it to space. The total depth and thermal inertia thus have an important effect on the required number of orbits. In most of our simulations 12 orbits were sufficient but some cases required up to 36 orbits. This number also depends on the actual duration of an orbit.
Fig. 3
Evolution of the rotation period of a planet starting from the 3:1 and 1:2 spinorbit resonances. All the parameters of the system, except e = 0, are set to the default values of Table 2. 
2.4. Lightcurves for a distant observer
For the configurations and at the wavelengths considered in this work, the contribution from reflected light remains negligible compared with the thermal emission. We do not include the contribution of the reflected light in the formula, although it is calculated by the model. The surface of the planet is divided into a longitudelatitude grid. Each cell j of the grid has a surface temperature T_{j} and an area S_{j}. The flux spectral density received by a distant observer at a distance d is
where α_{j} is the angle between the normal to the cell and the direction toward the observer and I_{j} is the specific intensity of the cell given by
where B_{λ} is the Planck function, ϵ_{λ} is the surface emissivity. This formula assumes a Lambertian (isotropic) distribution of intensity. We thus neglect the beaming (increase of intensity at low phase angles) that results from rugosity and craters. It would be interesting in a future study to include an anisotropic emission to quantify the photometric effect of thermal beaming. In practice, cosα_{j} is calculated as the scalar product of the direction vectors attached to the center of the planet and pointing towards the observer and the center of the j cell. Only locations visible to the observer (cosα_{j} > 0) contribute to φ_{P,λ}.
The spectral energy distribution of the stellar irradiation is not needed to compute the planetary emission, which is sensitive only to the bolometric luminosity and the Bond albedo. However, it is convenient to express the flux received by the planet (in a given spectral band) in terms of the planet/star contrast ratio. To do this we approximate the stellar radiation by a blackbody at the effective temperature of the star. When considering a specific stellar system or a comparison with real observations it is more appropriate to use a more realistic spectral density obtained from a model and/or observations, as we did in a separate study dedicated to 55 Cnc e (Bolmont et al. 2013).
Fig. 4
lightcurves and temperatures maps (as seen by the observer) for a circular orbit. a) These lightcurves are calculated for the default parameters (except for the eccentricity, which is null) and for different rotation periods. Tidal dissipation is not included. Only panel a) includes the case ω / n = 1 / 2. b) The lightcurves for ω / n = 3 / 2 and 3 are superimposed with that of ω / n = 2 because their thermal inertia has been modified in order to keep the quantity constant. c) same as b) but with the tidal dissipation included. d) same as c) but the half of the tidal dissipation is released in a localized hot spot. The temperature maps correspond to panel a). As in all following figures, these maps are displayed for the same orbit and duration as the lightcurves above and are separated by a constant time interval: one 24^{th} of the full duration (which can be 1 or 2 orbits depending on the figure). By default, time 0 corresponds to a phase angle of 90° (and periastron when e ≠ 0). 
3. Results
3.1. Circular orbits
We first consider a planet on a circular orbit and the effect that a nonsynchronized rotation has on the thermal lightcurve and on tidal heating signatures. In reality, planets that are sufficiently close to their star to exhibit an observable emission (above 10^{5} times the stellar flux) are subjected to strong tidal interactions that synchronize very rapidly the rotation with the orbital motion. We considered five planetary rotation periods: ω / n = 0.5, 1, 1.5, 2 and 3. These configurations would not realistically be maintained for long given the effects of tidal dissipation. Rather, any nonsynchronous states should actually evolve extremely rapidly towards synchronization, as illustrated in Fig. 3. Indeed, the planet initially rotating 3 times per orbit slows to 2 rotations per orbit in 10 000 years and is synchronized in less than 100 000 years. It is therefore unlikely that such a configuration would be observed. Nonetheless, we use the simple case of planets on circular orbits with varying rotation states to illustrate how eccentricity breaks degeneracies and gives rise to a rotational signature.
Figure 4 shows the phase curves of planets on circular orbits with a range of spin states. All other parameters were set to their default values listed in Table 2. Tidal dissipation is not included. For a circular orbit, the synchronized case is equivalent to a null thermal inertia as all points on the planetary surface receive a constant illumination. In non synchronized cases, the surface thermal inertia has two effects. First, it damps the amplitude of the lightcurve by lowering the zonal temperature gradient. Second, it induces a phase shift of the lightcurve, corresponding to a delay of its maximum (westward shift of the hottest point) if the rotation is faster than synchronization, and to an advance (eastward shift of the hottest point) if rotation is slower than synchronization, which is the case for the 1:2 spinorbit resonance. A projected inclination of 0° (a faceon orbit) would produce a flat curve (unlike in the case of an eccentric orbit), as would an infinite thermal inertia for any observer position even if the planet is not synchronized.
Panel b of Fig. 4 illustrates the degeneracy between rotation rate, thermal inertia and albedo. This is well known for solar system small bodies: in the absence of internal heat flow, the lightcurve is indeed controlled by a single thermal parameter proportional to (Spencer et al. 1989; Lagerros 1996). By modifying the thermal inertia of the ω / n = 3 / 2 and 3 cases to match the same as the ω / n = 2 case, the same lightcurves are obtained.
When including tidal dissipation and assuming the same dissipation factor k_{2,p}Δt_{p} the different rotation rates produce different internal heat fluxes. The heat flux increases with the departure from the synchronous rotation. Tidal heating increases the planetary flux and lowers the amplitude of the lightcurve by attenuating the day/night temperature contrast, as shown in panel c of Fig. 4. Although the previous degeneracy is broken, inferring the rotation rate requires a priori knowledge of the dissipation. If, however, the tidal dissipation is not released uniformly at the surface but mainly through one or a few large hot spots, like on Io (Veeder et al. 2012), then the planet behaves like a lighthouse and its lightcurve exhibits a modulation at the rotation period. In the phase curves in panel d of Fig. 4, half of the tidal dissipation is released through one hot spot centered on the equator. This region of increased internal heat flow is modeled with a Gaussian profile and a halfwidth of 45°. Changing the number, location, intensity and size of hot spots would of course affect the amplitude of the rotation modulation.
3.2. Eccentric orbits
The degeneracy between thermal inertia and rotation rate is broken for eccentric planets. The combination of the change of insolation during the orbital motion and the rotation of the planet produces a nonhomogeneous zonal heating: two meridians do not receive in general the same insolation when they face the star. Due to the surface thermal inertia, the resulting zonal temperature differences produce a photometric variation associated with the rotation that cannot, a priori, be obtained with a different set of rotation rates and thermal inertia. It is easier to understand the phenomenon for a very eccentric orbit. In this case the periastron passage is very brief and produces an intense heating that is limited to a small area of the planet, creating a hot region. After this periastron branding the rotating planet shines as an IR lighthouse, producing a rotation modulation that decays as the heated area cools down. At the next periastron passage, it is either another area of the planet that is heated or the same one for some spinorbit resonances (1:2, 2:1, 3:1, ...). In all cases, the photometric signature is periodic for a uniform planet because the shift between a heated region and the previous one remains constant.
Fig. 5
Effect of the thermal inertia on a pseudosynchronized planet (ω / n = 2.38) with the default parameters of Table 2. No tidal heat flux is included. The lightcurves and the surface temperature maps correspond to the same observation geometry and cover the two same orbits. On the temperature maps, the equator and the reference meridian (longitude 0° in our coordinate system) are indicated. 
3.2.1. Influence of the thermal inertia
The surface thermal inertia must be high enough for the rotation to produce an observable signature. This is why our default thermal inertia is 3000 SI (i.e. Js^{−1/2} K^{1} m^{2}). This value is typical of rocks but much higher than that of solar system airless bodies: 50 SI for the Moon, 80 SI for Mercury (Spencer et al. 1989). For asteroids it varies from a few tens of SI to a few hundreds SI (Delbo’ & Tanga 2009) and exceptionally reaches 1000 SI (Shepard et al. 2008). Low thermal inertia are due to the presence of a regolith, while bare, solid rocks have a much higher thermal inertia. On the surface of Mars, the thermal inertia varies from about 50 SI in areas covered with fine dust and very little bedrock, to 2500 SI for bedrock (Putzig & Mellon 2007). On Mars the thin atmosphere affects the effective thermal inertia of porous or dusty soils, as the heat transport by the gas particle increases the thermal conductivity of the medium compared to the heat transfer in vacuum.
Figure 5 illustrates the effect of thermal inertia on planetary phase curves, using the default values from Table 2. Given the expected precision of 10 ppm attainable by EChO and JWST for quiet stars, Fig. 5 shows that 1000 SI is roughly the minimum thermal inertia needed to have a measurable deviation from the lightcurves obtained with a null thermal inertia. Of course, there are no direct constraints on the (sub)surfaces of exoplanets of a few M_{⊕}; these could be partially made of materials other than regolith such as rocks and metals, solid or melted. A high surface thermal inertia cannot be ruled out. A layer of liquid at the surface or in the subsurface can result in a high thermal inertia. For instance, the 25 m average mixed layer at the surface of Earth’s oceans has an equivalent thermal inertia of ~25 000 SI (Forget, priv. comm.). A melted (sub)surface is consistent with the high rates of tidal dissipation and resulting internal heat flow associated with eccentric shortperiod orbit (Barnes et al. 2010). A melted surface implies an atmosphere, which would be inconsistent with our airless assumption. However, vapor pressures of refractory materials are low and the resulting atmospheres could have negligible opacities and negligible influence on the emission and the surface temperature. For instance, a pressure below 1.5 Pa was calculated by Léger et al. (2011) above a substellar magma ocean on CoRoT 7b. The highest value used for thermal inertia in this this graph, 10 × Γ_{0}, is probably unrealistic as it exceeds that of pure iron but it illustrates the effect of increasing Γ: for an infinite thermal inertia the surface behaves as if it were exposed to the timeaveraged insolation and produces a flat lightcurve.
Fig. 6
Same system as above but for Γ_{0} and for different “ecliptic” longitudes of the observer. 
Fig. 7
Influence of the inclination on the lightcurves and apparent temperature maps. (Default parameters, equilibrium rotation: ω / n = 2.38, tidal heat flux not included). For the inclination of 90° the secondary eclipses appear on the lightcurves and the moments of transit are indicated. 
If observations of a given planetary system were to show a significant departure from the Γ = 0 lightcurve, it would represent an extremely valuable target. An elevated thermal inertia – with a minimum value determined by the sensitivity of the measurement (see Fig. 5) – would provide a direct constraint on the nature of the planetary surface. In addition, the rotation period can be inferred or at least constrained for planets with large thermal inertia. This of course requires that the existence of an atmosphere could be discarded by other means; that will be the subject of a future study but is already discussed in Maurin et al. (2012). We do not address here the case of a planetary surface with an inhomogeneous thermal inertia, for instance with large highinertia spots on an overall lowinertia surface, or viceversa.
Fig. 8
Top: lightcurves for different eccentrics and rotation rates. The solid gray lines indicate the equilibrium rotation (pseudosynchronization). Bottom: surface temperature maps for e = 0.5. Parameters other than eccentricity are set to default values and tidal heat flux is not included. 
3.2.2. Influence of the observation geometry
For a given inclination (or subobserver latitude), Fig. 6 shows that the shape of the lightcurve depends strongly on the ecliptic^{1} longitude of the observer. It matters whether one observes the day side of the planet near the periastron or the apoastron (even for a null thermal inertia), or whether one observes the quarter that precedes or follows the periastron. For example, when comparing the lightcurves seen by observers 1 and 3 in Fig. 6, the much larger amplitude seen by observer 3 comes from the fact that observer 3 sees the planet’s day side during periastron whereas observer 1 does not. A precise knowledge of the orbital ephemeris from radial velocity or intertransit duration measurements is therefore required to further analyze the lightcurves in terms of rotation and thermal inertia.
Our default inclination in this study is 90° (transit configuration) but the lightcurve variations are in theory observable for any inclination (see Kane & Gelino 2011b, for the effect on reflected lightcurves). Unlike for circular orbits, a photometric modulation is observable at null inclination due to the variation of the planet reflection and emission with orbital distance. Fig. 7 shows how the observed thermal lightcurve changes with the inclination. The amplitude of the lightcurve variations is of course the largest for i = 90°, an inclination that also provides constraints on the radius and the absolute planetary flux from to primary transit and secondary eclipse observations, respectively. However, the i = 90° and i = 60° lightcurves differ by less than 10%. Half of randomly oriented orbits present an inclination between 60 and 90° and similarly strong observable photometric signatures.
Note that we represented the secondary eclipse (the planet disappearing behind the stellar disk) only in Fig. 7 although all the other figures are also done for i = 90°. Our code can produce lightcurves that include the secondary eclipses but we chose not to show them as we focus on the orbital photometry. Primary transits do not affect the planet’s emission but do somewhat increase the planet/star contrast ratio due to the apparent dimming of the star. For our default parameters, however, there is only a 0.3% increase in contrast ratio.
3.2.3. Interplay between eccentricity and rotation
Figure 8 shows lightcurves for different eccentricities and rotation rates. Several spinorbit resonances are considered: 1:2, 1:1 (synchronization), 3:2, 2:1 as well as pseudosynchronization. First, one can notice that, unlike circular orbits, eccentric orbits do not exhibit the same lightcurve for a null thermal inertia and synchronization. The surface of an eccentric synchronized planet does not receive a constant illumination because of changes in orbital distance and socalled optical librations from the variations of the orbital motion angular velocity. The position of the substellar point oscillates and the permanently dark area of the planet is smaller than a full hemisphere. This dark region can be seen as an “eye of Sauron” (Tolkien 1954) in the temperature map of Fig. 8 calculated for ω / n = 1 and on Fig. A.1. For a null obliquity, this region covers a fraction s of the planetary surface given by the following series expansion (23)The exact analytical calculation of this fractional area s is given in Appendix A. For eccentricities larger than ~0.724, all longitudes receive starlight at some position on the orbit. For our default eccentricity of 0.45, only 20.6% of the surface of a synchronized planet never receives starlight. Selsis et al. (2007b) also gives a series expansion for s, but which is valid only for eccentricities lower than about 0.3.
The effect of the rotation rate on the lightcurve increases with the eccentricity. For e = 0.7 and pseudosynchronization (ω / n = 6.5), one can clearly see the effect of the periastron branding at time = 0 and the damped rotation modulation after the periastron passages that produces 6 local maxima in the orbital phase curve.
It is also important to note that, in particular for high eccentricities, the curves obtained for different rotations but the same thermal inertia converge to similar values near the periastron passage. Again, this is easiest to understand for high eccentricities: near periastron, the variation of zenith angle on a point of the surface of the planet is dominated by the angular velocity of the (fast) orbital motion over the angular velocity of the rotation. As a consequence, the lightcurve close to periastron provides a constraint on the thermal inertia. Closer to apoastron, it is the rotation angular velocity that dominates the modulation of surface irradiation and the lightcurve is sensitive to both thermal inertia and rotation. Therefore, rotation and thermal inertia should in theory be inferred, at least for a transiting planet with known inclination and radius, and a reference for the absolute stellar flux obtained during the secondary eclipse.
Fig. 9
Effect of tidal heating on pseudosynchronized planets. All parameters are set to default values except the eccentricity. The dashed line shows the lightcurve at radiative equilibrium (Γ = 0) and without tidal heating. The solid line is obtained when including the subsurface heat transport but no tidal heating. The dotted line is calculated with heat diffusion and a uniform internal tidal heat flow. To produce the dasheddotted line, half of the tidal heat is released through a hot spot with a 45° radius and centered on the equator (see text). Bottom: Surface temperatures maps corresponding to the lightcurves calculated for e = 0.5. 
3.2.4. Influence of the stellar spectral type at constant T_{eq}
For a null surface thermal inertia and a given ω / n ratio the shape and intensity of the thermal lightcurve does not change with the spectral type if we scale the orbital distance to keep the equilibrium temperature of the planet constant. A nonnull thermal inertia makes the thermal lightcurve sensitive to the orbital period, as illustrated in Fig. 10 (top). The irradiation of the surface changes more slowly for a longer orbital period, giving more time for the surface to adapt to the change of illumination. The more massive the star, the longer the orbital period, and the closer the lightcurve to the case of radiative equilibrium (Γ = 0).
For a given equilibrium temperature, lowmass host stars provide better candidates to attempt a measurement of the rotation period. The signature of rotation is stronger for short orbital periods. In addition, Fig. 10 (bottom) shows that the infrared planet/star contrast ratio strongly decreases with increasing stellar mass due to the massluminosity relationship.
3.2.5. Influence of tidal dissipation
We now include tidal heating computed with the constant timelag model as described in Sect. 2.2. It is important to note that the terrestrial value we use for the dissipation factor is arbitrary and could differ by orders of magnitude from that of an exoplanet. In particular, the very strong dissipation rates and heat flows considered here would dramatically alter the internal structure and, thus, the dissipation properties of the planet. A more consistent coupling between the dissipation factor and the dissipation rate could be included in the future, provided some assumptions on the composition of the planet.
As in the circular case described in Sect. 3.1, we assume that the tidal heat is either released uniformly over the planetary surface or half is released uniformly and the other half through a hot spot. This hot spot is arbitrarily centered on the equator with a gaussian profile of heat flow and a halfwidth of 45°. In our model the properties, location and size of the hot spot(s) can be changed but the purpose here is not to present an exhaustive zoology of all the lightcurves that can be produced by different distributions of surface heat flows. We want to illustrate the fact that a nonuniform heat flow produces a rotation modulation that could be observable for strong tidal heating.
Figure 9 shows the resulting lightcurves with and without the tidal heating for different eccentricities. For the default parameters and e = 0.2 the combination of rotation and thermal inertia produce some moderate damping of the curve but tides have a negligible effect. For e = 0.3 the signature of rotation appears as well as a slight increase of the planet brightness. For e = 0.4 and 0.5, the tidal heat flux imparts a strong signal on the phase curves. In the case of a uniform release of tidal heat, the signature of rotation is attenuated, in particular at e = 0.5 as the emission is dominated by the tidal flux. In the case of the hot spot, the rotation modulation is strong and does not decay over an orbit as periastron branding does. In the case of a hot spot, only spinorbit resonances produce a periodic signal (with a period that can be one or more orbital periods). In Fig. 9 the curves are calculated for pseudosynchronization. This results for e = 0.5 in ω / n = 2.84, and the nonperiodicity of the curve can be seen, for instance by looking at successive maxima, which intensity change depending whether the hot spot and periastron heating are superimposed or shifted. Assuming a single strong hot spot obviously maximizes the rotation modulation. Multiple or weaker hot spots would naturally produce a weaker signature.
Fig. 10
Lightcurves calculated for a planet receiving an average insolation of 1021 Wm^{2} (T_{eq} = 366 K with A = 0.1) but orbiting different stars. All parameters are set to default values except the stellar mass and orbital period. The flux from the planet received at 10 pc is given in the top panel, while the bottom one shows the planet/star contrast ratio. Corresponding orbital periods P and semimajor axes a are indicated. 
Fig. 11
lightcurves for low thermal inertia halfdark halfbright planet. Dashed lines show the lightcurves for a planet with the default parameters but a low thermal inertia (100 SI) and albedo values between 0.1 and 0.9. The solid curve is obtained if the surface of the planet is divided into two hemispheres (separated by two meridians) with highly different albedos: A = 0.1 and 0.9. On the top panel, the system is in a 2:1 spinorbit resonance. The lightcurve is periodic but its shape depends on the arbitrary position of the hemispheres. On the bottom panel, the system is in pseudosynchronization and the signal is no longer periodic. 
4. Discussions
4.1. Retrieval of planet parameters
The retrieval of the rotation period from lightcurves obtained with future space IR telescopes requires a realistic description of the instrument, sources of noise and stellar variability, and will be the purpose of future studies. Although such work has not yet been done, we did perform χ^{2} minimizations to test the nondegeneracy between ω and Γ. We considered a transiting 2 R_{⊕} planet on a 20 day orbit around a 0.3 M_{⊙} star with an eccentricity of 0.3 (with these parameters the contribution of tidal heating to the thermal emission can safely be neglected). Assuming the radius, inclination and orbital period of the planet are known, which is the case for a transiting planet, the 3 remaining parameters controlling the lightcurves are A, ω and Γ. We produced lightcurves for 2 consecutive orbits for five 2μmwidth bands covering the 6−16 μm range on a grid of A, ω and Γ. We then simulated an observation with EChO (Tinetti et al. 2012) for a set of values: A_{p}, ω_{p} and Γ_{p} assuming a distance of 10 pc, a detector throughput of 30%, and a full coverage of the 2 orbits. Considering only the stellar photon noise, we randomly added a gaussian noise to every exposure (>24 per orbit as the result is independent of this number above this value) and produced 500 noisy observations. We computed the individual χ^{2} between each noisy lightcurve and those from the grid and produced a 3D map of the mean χ^{2} in the AωΓ space. We then identified the n− σ contours around the actual A_{p}, ω_{p} and Γ_{p} values.
The values we set for the albedo and rotation period are A_{p} = 0.3 and ω_{p} / n = 2 and we considered two values for the thermal inertia: Γ = 1500 SI (Γ_{0} / 2) and Γ = 3000 SI (Γ_{0}). As noted in Sect. 3.2.1 the specific signature of rotation increases with the thermal inertia and we thus expect a better retrieval with the highest value of Γ. With Γ = 1500 SI we found the following 2σ errors: , SI and A = 0.3 ± 0.05. With Γ = 3000 SI, the error decreases as expected: , Γ = 3000 ± 800 SI and A = 0.3 ± 0.05 (error on the albedo is unchanged). Again, these results are more a practical confirmation that the effects of rotation and thermal inertia, which are indistinguishable for a circular orbit, are broken apart by eccentricity, rather than a robust prediction of what can be expected from real observations. It is however promising that the spin might be constrained even with EChO (a 1.1 m telescope) while JWST will be able to obtain a much higher signaltonoise ratio with a collecting area 25 times larger. For outoftransit exoplanets, the radius and the inclination would also have to be constrained and the retrieval of the five parameters has not been tested yet.
4.2. Rotation signature at low thermal inertia from albedo spots
As shown in Sect. 3.2.1, a strong thermal inertia is needed to impart a measurable signature of planetary rotation on the phase curve. This is true for a planetary surface with uniform properties but not necessarily for a planet with a very heterogeneous albedo. An extreme case is illustrated in Fig. 11. This figure shows that a nonsynchronized planet with low surface thermal inertia but with dark and bright hemispheres (separated by two meridians) produces a rotation modulation that superimposes with the orbital response. A more realistic albedo surface distribution with less contrast and smaller patterns would, of course, result in a weak modulation with an amplitude unlikely above 10 ppm.
4.3. Anisotropic and spectral radiative properties of the surface
Our model assumes an emissivity of 1 independent of the wavelength. In reality, planetary surfaces are not blackbodies and their emissivity (although usually close to unity) vary with wavelength, depending on the surface material. There is for instance a well known absorption band around 10 μm for silicates. These spectral properties are used to determine the composition of planetary surfaces in the solar system. The emissivity of solar system regoliths have been studied in details (e.g. Lim et al. 2005) but regolith have a low thermal inertia and would not exhibit a measurable signature of rotation. The emissity of high thermal inertia refractory surfaces (bedrocks, magma) is still to be studied. Hu et al. (2012) produced synthetic spectra of planetary surfaces of various compositions that could be expected for highly irradiated exoplanets, but significant emissity spectral features are associated with dusty/porous the surface (Vernazza et al. 2012). Assuming that a high thermal inertia is inconsistent with a dusty/porous surface, the spectral signatures of highinertia surfaces is not expected to exhibit strong spectral features.
The spectral dependency of the albedo is not important as long as we do not include the reflected light, as we are only sensitive to the surface albedo averaged over the wavelengths of the insolation. The surface albedo may however depend on the incidence angle, as it does on the Moon and generally on low albedo regoliths. This is why the full Moon appears flatter than a Lambertian sphere. We performed some tests comparing lightcurves with a LommelSeeliger bidirectional reflectance distribution function (BRDF) for the reflected light (Fairbairn 2005) and with an equivalent constant Bond albedo. We found negligible effect on the thermal lightcurve due to the fact that areas under grazing illumination have a low temperature and do not contribute significantly to the diskintegrated emission. The effect on the temperature is very small and the latitudedependence of temperature on the starlit hemisphere of the Moon (for which we can neglect the very low thermal inertia) follows fairly well a T_{SS}μ^{1 / 4} law (Lawson et al. 2000), where T_{SS} is the subsolar point.
Fig. 12
Influence of the obliquity β. Parameters are set to their default values and ω = ω_{eq} = 2.38n. Panels a) and b): lightcurves for different obliquities and two orientations of the spin axis, without tidal heating. Panels c) and d): Same as b) but with tidal heating included, uniformly in c) and with half of the tidal dissipation released through one hot spot (Gaussian, half width: 45°) centered on the equator in d). 
The thermal emission of solar system planetary surfaces is known to be slightly anisotropic, with a beaming toward the direction of illumination (similar to the opposition surge known for reflected light). This is due to macroscopic rugosity and craters. For low zenith angles, a crater floor is not affected by the shadows of its edges but, on the contrary, the crater floor receives the thermal emission from the edges in addition to the stellar illumination, resulting in a higher brightness temperature for an observer that is roughly in the direction of the star. The effect of rugosity can be understood if the surface is modeled as a layer of spheres. Each illuminated sphere has as hot substellar point (assuming a low thermal inertia) and emits more in the direction of the illumination. As the opposition effect, this beaming is the strongest for very low phase angles. For non transiting exoplanets and transiting exoplanets, the phase angle for a distant observer cannot be lower than the angular radius of the star, which limits the effect of the beaming. This is in particularly true for very closein exoplanets for which the angular stellar radius can be very high (~15° for Kepler 10 b and 55 Cnc e). New thermal modeling of asteroids shows that the beaming effect is not only important at low phase angles, and that the most important effect is coming not from the subsolar point but from the rough surface elements near the terminator (Rozitis & Green 2011). While we neglect these effects here, it might be important to include them in a forthcoming study.
Taking into account spectral and anisotropic effects of realistic planetary surfaces can be done by computing BRDFs with Hapke functions and associated parameters (Hapke 1981). This may be done in the future, at least to assess the impact on the retrieval of the rotation period but it would imply adding several parameters, as unknown as the nature of exoplanetary surfaces. In addition, Hapke functions have been developed to reproduce the BRDF of regolith and may not be adapted for different types of surfaces that could be found on exoplanets.
4.4. Obliquity
To this point we have assumed that obliquity tides have damped any initial obliquity to zero. For isolated planets with a planet/star IR contrast ratio above 10 ppm, obliquity damping time is indeed very efficient: for our default parameters, an initial obliquity of 80° decreases to 5° in just 20 000 yrs. Dissipation factors must be 5 orders of magnitude lower than Earth’s to preserve a high initial obliquity for more than 1 Gyr, although this situation cannot be ruled out as uncertainties on the dissipation span many orders of magnitude. In the presence of planetary companions, however, the planet can be captured into a Cassini state, a resonance between spin precession and orbital precession, with any nonzero obliquity (Colombo 1966; Peale 1969). Our model is designed to include an obliquity but this adds to the parameters list the two angles that define the orientation of the spin axis: the obliquity and, for instance, the periastron longitude (which is the angle between the spring equinox and the periastron passage).
Figure 12 shows a limited exploration to test the effect of obliquity on the phase curve. Without tidal dissipation, the effect on the lightcurve is small (<1%) if the obliquity is smaller than 10° (panels a) and b) of Fig. 12). The viewing and orbital geometries play an important role and can introduce new degeneracies, in particular between thermal inertia and obliquity. Indeed, in panel a) of Fig. 12 the obliquity produces, for this particular observation geometry, a damping and phase shift of the phase curve that mimic a thermal inertia.
Panel graph (c) of Fig. 12 shows that, when tidal heating is included with our default dissipation factor, the planet emission increases with the obliquity. However, unless the dissipation can be inferred by independent measurement, this cannot be used directly to constrain the obliquity. Panel graph (d) of the same figure shows the thermal lightcurve obtained when 50% of the tidal heating is released through a hot spot (with arbitrary profile, position and size). Cases like (c) and (d), in which both the obliquity and the tidal dissipation are high, should not last for long and are unlikely to be observed: the high dissipation rate would either have eroded the initial obliquity or prevented/swept away possible Cassini states (Fabrycky et al. 2007).
5. Conclusions
In this paper we studied the effect of a planet’s rotation on its infrared photometric lightcurve. We also included the effect of tidal dissipation within the planet, the strength of which is a strong function of the planet’s orbital eccentricity (see Eq. (4)). We generated a wide range of synthetic IR phasecurves to determine the most important parameters.
For an airless planet on a circular orbit, the effect of surface thermal inertia and rotation on the thermal lightcurve cannot be distinguished. This degeneracy is broken for an eccentric orbit by an effect we call periastron branding – the peak of heating at periastron produces a longitudinal gradient of temperature, which in turns result in a rotation modulation – allowing for the rotation and the surface thermal inertia to be measured or, at least, constrained.
For an eccentric planet captured into the 1:1 spinorbit resonance we provide an analytical expression as well as a series expansion for the fractional area that never receives starlight as a function of eccentricity (Appendix A below).
The higher the eccentrity and surface thermal inertia, the better the constraint on the rotation period. Below a value of about 1000 SI, rotation and thermal inertia have a negligible effect on the thermal lightcurve. Constraining these parameters requires a super Mercury with a super thermal inertia. This implies solid or melted, nonporous surface different from the regoliths found on airless bodies in the solar system.
Tidal heating can affect the radiative budget of the planet and its thermal lightcurve. If released homogeneously over the surface of the planet, tidal heat increases the overall infrared emission and damps the thermal lightcurve variations associated with the insolation and rotation. If there is a longitudinal asymmetry in the surface tidal heat flux (like on Io) then the rotation signature can be enhanced. With hot spot tidal heating, a planet behaves as an IR lighthouse regardless of the thermal inertia. Albedo or thermal inertia spots can produce a similar “lighthouse” effect.
A photometric precision of 10 ppm, stable over at least two orbital periods, is needed to retrieve the rotation period, in the best case scenario. This is compatible with the planned characteristics of EChO (Tinetti et al. 2012). Assuming photonnoise limited observations and a perfectly stable star, preliminary retrieval calculations are promising for EChO and JWST. However, we note that these calculations were made for synthetic eccentric and transiting exoplanets whereas the known candidates with significant eccentricities are outoftransit (see Table 1). Stellar variability and instrument stability (Crossfield et al. 2012, e.g.) are expected to be the main obstacle.
Measuring the rotation period of tidallyevolved eccentric planets would be extremely important to test tidal models in conditions that cannot be studied in the solar system. The rotation period can have dramatic consequences on climate and habitability and the first potentially habitable planets discovered, GJ581d (Selsis et al. 2007b; Makarov et al. 2012; Wordsworth et al. 2011) and GJ667Cc (Delfosse et al. 2013) are both eccentric and subjected to strong tides. The study of super Mercuries thermal lightcurves could provide a better understanding of tidal interactions and, indirectly, on habitability.
Acknowledgments
We thank the anonymous referee for helpful comments. F.S. acknowledges support from the European Research Council (Starting Grant 209622: E_{3}ARTHs). J.L. acknowledges support from the DIM ACAV. F.H. and S.N.R. abide.
References
 Barman, T. S., Hauschildt, P. H., & Allard, F. 2005, ApJ, 632, 1132 [NASA ADS] [CrossRef] [Google Scholar]
 Barnes, R., Raymond, S. N., Greenberg, R., Jackson, B., & Kaib, N. A. 2010, ApJ, 709, L95 [NASA ADS] [CrossRef] [Google Scholar]
 Batalha, N. M., Borucki, W. J., Bryson, S. T., et al. 2011, ApJ, 729, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Beust, H., Bonfils, X., Montagnier, G., Delfosse, X., & Forveille, T. 2012, A&A, 545, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bolmont, E., Selsis, F., Raymond, S. N., et al. 2013, A&A, in press, DOI: 10.1051/00046361/201220837 [Google Scholar]
 Bonfils, X., Mayor, M., Delfosse, X., et al. 2007, A&A, 474, 293 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bonfils, X., Delfosse, X., Udry, S., et al. 2013, A&A, 549, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bouchy, F., Mayor, M., Lovis, C., et al. 2009, A&A, 496, 527 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Colombo, G. 1966, AJ, 71, 891 [NASA ADS] [CrossRef] [Google Scholar]
 Cowan, N. B., & Agol, E. 2011, ApJ, 726, 82 [Google Scholar]
 Cowan, N. B., Voigt, A., & Abbot, D. S. 2012, ApJ, 757, 80 [NASA ADS] [CrossRef] [Google Scholar]
 Crossfield, I. J. M., Hansen, B. M. S., Harrington, J., et al. 2010, ApJ, 723, 1436 [NASA ADS] [CrossRef] [Google Scholar]
 Crossfield, I. J. M., Knutson, H., Fortney, J., et al. 2012, ApJ, 752, 81 [NASA ADS] [CrossRef] [Google Scholar]
 Cubillos, P., Harrington, J., Madhusudhan, N., et al. 2013, ApJ, 768, 42 [NASA ADS] [CrossRef] [Google Scholar]
 Delbo’, M., & Tanga, P. 2009, Planet. Space Sci., 57, 259 [NASA ADS] [CrossRef] [Google Scholar]
 Delfosse, X., Bonfils, X., Forveille, T., et al. 2013, A&A, 553, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Demory, B.O., Gillon, M., Seager, S., et al. 2012, ApJ, 751, L28 [Google Scholar]
 Dressing, C. D., & Charbonneau, D. 2013, ApJ, 767, 95 [NASA ADS] [CrossRef] [Google Scholar]
 Fabrycky, D. C., Johnson, E. T., & Goodman, J. 2007, ApJ, 665, 754 [NASA ADS] [CrossRef] [Google Scholar]
 Fairbairn, M. B. 2005, JRASC, 99, 92 [NASA ADS] [Google Scholar]
 Fischer, D. A., Gaidos, E., Howard, A. W., et al. 2012, ApJ, 745, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Ford, E. B., Seager, S., & Turner, E. L. 2001, Nature, 412, 885 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Forveille, T., Bonfils, X., Delfosse, X., et al. 2011, A&A, submitted [arXiv:1109.2505] [Google Scholar]
 Gaidos, E., & Williams, D. M. 2004, New Astron., 10, 67 [NASA ADS] [CrossRef] [Google Scholar]
 GómezLeal, I., Pallé, E., & Selsis, F. 2012, ApJ, 752, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Hapke, B. 1981, J. Geophys. Res., 86, 4571 [NASA ADS] [Google Scholar]
 Hu, R., Ehlmann, B. L., & Seager, S. 2012, ApJ, 752, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Hut, P. 1981, A&A, 99, 126 [NASA ADS] [Google Scholar]
 Iro, N., & Deming, L. D. 2010, ApJ, 712, 218 [NASA ADS] [CrossRef] [Google Scholar]
 Iro, N., Bézard, B., & Guillot, T. 2005, A&A, 436, 719 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kane, S. R., & Gelino, D. M. 2011a, ApJ, 741, 52 [NASA ADS] [CrossRef] [Google Scholar]
 Kane, S. R., & Gelino, D. M. 2011b, ApJ, 729, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Kane, S. R., & Gelino, D. M. 2013, ApJ, 762, 129 [NASA ADS] [CrossRef] [Google Scholar]
 Kane, S. R., Gelino, D. M., Ciardi, D. R., Dragomir, D., & von Braun, K. 2011, ApJ, 740, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Knutson, H. A., Charbonneau, D., Allen, L. E., et al. 2007, Nature, 447, 183 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Lagerros, J. S. V. 1996, A&A, 310, 1011 [NASA ADS] [Google Scholar]
 Lawson, S. L., Jakosky, B. M., Park, H., & Mellon, M. T. 2000, J. Geophys. Res., 105, 4273 [NASA ADS] [CrossRef] [Google Scholar]
 Leconte, J., Chabrier, G., Baraffe, I., & Levrard, B. 2010, A&A, 516, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Léger, A., Grasset, O., Fegley, B., et al. 2011, Icarus, 213, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Levrard, B., Correia, A. C. M., Chabrier, G., et al. 2007, A&A, 462, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lewis, N. K., Knutson, H. A., Showman, A. P., et al. 2013, ApJ, 766, 95 [NASA ADS] [CrossRef] [Google Scholar]
 Lim, L. F., McConnochie, T. H., Bell, J. F., & Hayward, T. L. 2005, Icarus, 173, 385 [Google Scholar]
 Lissauer, J. J., JontofHutter, D., Rowe, J. F., et al. 2013 [arXiv:1303.0227] [Google Scholar]
 Makarov, V. V., & Efroimsky, M. 2013, ApJ, 764, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Makarov, V. V., Berghea, C., & Efroimsky, M. 2012, ApJ, 761, 83 [NASA ADS] [CrossRef] [Google Scholar]
 Maurin, A. S. 2012, Ph.D. Thesis, Université de Bordeaux I [Google Scholar]
 Maurin, A. S., Selsis, F., Hersant, F., & Belu, A. 2012, A&A, 538, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mayor, M., Marmier, M., Lovis, C., et al. 2011, A&A, submitted [arXiv:1109.2497] [Google Scholar]
 Menou, K. 2012, ApJ, 744, L16 [Google Scholar]
 Mikkola, S. 1987, Celest. Mech., 40, 329 [Google Scholar]
 Moskovitz, N. A., Gaidos, E., & Williams, D. M. 2009, Astrobiology, 9, 269 [NASA ADS] [CrossRef] [Google Scholar]
 Neron de Surgy, O., & Laskar, J. 1997, A&A, 318, 975 [NASA ADS] [Google Scholar]
 Pallé, E., Ford, E. B., Seager, S., MontañésRodríguez, P., & Vazquez, M. 2008, ApJ, 676, 1319 [NASA ADS] [CrossRef] [Google Scholar]
 Peale, S. J. 1969, AJ, 74, 483 [NASA ADS] [CrossRef] [Google Scholar]
 Pepe, F., Correia, A. C. M., Mayor, M., et al. 2007, A&A, 462, 769 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Petigura, E. A., Marcy, G. W., & Howard, A. W. 2013 [arXiv:1304.0460] [Google Scholar]
 Putzig, N. E., & Mellon, M. T. 2007, Icarus, 191, 68 [NASA ADS] [CrossRef] [Google Scholar]
 Rivera, E. J., Laughlin, G., Butler, R. P., et al. 2010, ApJ, 719, 890 [Google Scholar]
 Rozitis, B., & Green, S. F. 2011, MNRAS, 415, 2042 [NASA ADS] [CrossRef] [Google Scholar]
 Selsis, F. 2004, in Extrasolar Planets: Today and Tomorrow, eds. J. Beaulieu, A. Lecavelier , & C. Terquem, ASP Conf. Ser., 321, 170 [Google Scholar]
 Selsis, F., Chazelas, B., Bordé, P., et al. 2007a, Icarus, 191, 453 [NASA ADS] [CrossRef] [Google Scholar]
 Selsis, F., Kasting, J. F., Levrard, B., et al. 2007b, A&A, 476, 1373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Selsis, F., Wordsworth, R. D., & Forget, F. 2011, A&A, 532, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Shepard, M. K., Clark, B. E., Nolan, M. C., et al. 2008, Icarus, 193, 20 [NASA ADS] [CrossRef] [Google Scholar]
 Showman, A. P., Fortney, J. J., Lian, Y., et al. 2009, ApJ, 699, 564 [NASA ADS] [CrossRef] [Google Scholar]
 Spencer, J. R., Lebofsky, L. A., & Sykes, M. V. 1989, Icarus, 78, 337 [NASA ADS] [CrossRef] [Google Scholar]
 Tinetti, G., Beaulieu, J. P., Henning, T., et al. 2012, Exp. Astron., 34, 311 [NASA ADS] [CrossRef] [Google Scholar]
 Tolkien, J. R. R. 1954, The Fellowship of the Ring (George Allen & Unwin ) [Google Scholar]
 Veeder, G. J., Davies, A. G., Matson, D. L., et al. 2012, Icarus, 219, 701 [NASA ADS] [CrossRef] [Google Scholar]
 Vernazza, P., Delbo, M., King, P. L., et al. 2012, Icarus, 221, 1162 [NASA ADS] [CrossRef] [Google Scholar]
 Welsh, W. F., Orosz, J. A., Seager, S., et al. 2010, ApJ, 713, L145 [NASA ADS] [CrossRef] [Google Scholar]
 Wordsworth, R. D., Forget, F., Selsis, F., et al. 2011, ApJ, 733, L48 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Analytical calculation of the permanently dark area for synchronous planets on eccentric orbits
For a synchronous planet with a null obliquity on an eccentric orbit, the fraction of the surface that receives starlight at some point of the orbit is greater than 1/2 because the varying orbital motion will cause a libration of the star as seen from the planet. In other words, the longitude of the substellar point (A.1)varies with time (see Sect. 2.1). The amplitude of the libration is given by extrema of lon_{⋆}(t) that we calculate hereafter.
Fig. A.1
The eye of Sauron. Temperature maps showing the shrinking of the permanently dark area with increasing eccentricities for a synchronously rotating planet. Except for e, default parameters are used. 
Because lon_{⋆}(t) is a smooth function, we search for the values of t that cancel its derivative. However, we cannot use the time variable, which requires to solve Kepler’s equation. The idea is thus to use the true anomaly as our variable (the change of variable is smooth) and to search for ν_{ext} such that (A.2)where the operator stands for the temporal derivative. For a Keplerian orbit, (A.3)where n is the orbital mean motion. Substituting Eq. (A.3) in Eq. (A.2), recognizing that ω = n for a synchronous rotation, and keeping only the real solutions, we find (A.4)which defines the locations (ν_{ext} ≡ ν_{+} and ν_{−}) of our two extrema over one orbit. To find the values of the extrema themselves, we have to integrate (A.5)where ν_{0} is an arbitrary choice for the reference of the true anomaly that changes C but not the physical results. Again using Eq. (A.3), we find that (A.6)Finally, the fraction of the surface that never receives daylight is given by (A.7)which can be computed using Eqs. (A.4) and (A.6) When this fraction becomes negative (around e_{0} ≈ 0.724), it just means that the anti stellar point at the periastron passage receives light once between periastron and apastron, and once between apastron and periastron. Over the range of eccentricities where s is meaningful ([0,e_{0}]), the first terms of expansion of this function, (A.8)reproduce the latter with an absolute error below 0.002.
This calculation assumes a pointlike star and should be corrected to account for the angular extension of the star, which reduces the value of s and the critical eccentricty at which this area disappears.
All Tables
All Figures
Fig. 1
Number of rotations per orbital period for a pseudosynchronized planet as a function of eccentricity. The filled circle shows the pseudosynchronization (ω_{eq} / n = 2.38) for the default value of eccentricity used in this study (e = 0.45). 

In the text 
Fig. 2
Observability of the 10 μm emission excess due to tidal heating in a periodeccentricity diagram. a) Dashed contours indicate the circularization timescale (e / ė), solid contours give the planet/star contrast ratio in ppm, and the dotted contours show the ratio between the emission due to tidal heating and the emission without tidal heating. The contours are calculated for a 0.2 M_{⊙} star, A = 0.2, an isothermal surface (stellar and internal flux fully redistributed), a 2 R_{⊕} planet and the dissipation parameters of the Earth. The colored area represents the observability domain for the tidal heating, where the configuration lasts for more than 1 Gyr, produce a planet/star contrast higher than 10 ppm and results in tidal heating contributing to more than 10% of the emission. Panels b), c) and d) respectively show the influence of the stellar mass, the planetary radius and the dissipation factor on this tidal heating observability domain. 

In the text 
Fig. 3
Evolution of the rotation period of a planet starting from the 3:1 and 1:2 spinorbit resonances. All the parameters of the system, except e = 0, are set to the default values of Table 2. 

In the text 
Fig. 4
lightcurves and temperatures maps (as seen by the observer) for a circular orbit. a) These lightcurves are calculated for the default parameters (except for the eccentricity, which is null) and for different rotation periods. Tidal dissipation is not included. Only panel a) includes the case ω / n = 1 / 2. b) The lightcurves for ω / n = 3 / 2 and 3 are superimposed with that of ω / n = 2 because their thermal inertia has been modified in order to keep the quantity constant. c) same as b) but with the tidal dissipation included. d) same as c) but the half of the tidal dissipation is released in a localized hot spot. The temperature maps correspond to panel a). As in all following figures, these maps are displayed for the same orbit and duration as the lightcurves above and are separated by a constant time interval: one 24^{th} of the full duration (which can be 1 or 2 orbits depending on the figure). By default, time 0 corresponds to a phase angle of 90° (and periastron when e ≠ 0). 

In the text 
Fig. 5
Effect of the thermal inertia on a pseudosynchronized planet (ω / n = 2.38) with the default parameters of Table 2. No tidal heat flux is included. The lightcurves and the surface temperature maps correspond to the same observation geometry and cover the two same orbits. On the temperature maps, the equator and the reference meridian (longitude 0° in our coordinate system) are indicated. 

In the text 
Fig. 6
Same system as above but for Γ_{0} and for different “ecliptic” longitudes of the observer. 

In the text 
Fig. 7
Influence of the inclination on the lightcurves and apparent temperature maps. (Default parameters, equilibrium rotation: ω / n = 2.38, tidal heat flux not included). For the inclination of 90° the secondary eclipses appear on the lightcurves and the moments of transit are indicated. 

In the text 
Fig. 8
Top: lightcurves for different eccentrics and rotation rates. The solid gray lines indicate the equilibrium rotation (pseudosynchronization). Bottom: surface temperature maps for e = 0.5. Parameters other than eccentricity are set to default values and tidal heat flux is not included. 

In the text 
Fig. 9
Effect of tidal heating on pseudosynchronized planets. All parameters are set to default values except the eccentricity. The dashed line shows the lightcurve at radiative equilibrium (Γ = 0) and without tidal heating. The solid line is obtained when including the subsurface heat transport but no tidal heating. The dotted line is calculated with heat diffusion and a uniform internal tidal heat flow. To produce the dasheddotted line, half of the tidal heat is released through a hot spot with a 45° radius and centered on the equator (see text). Bottom: Surface temperatures maps corresponding to the lightcurves calculated for e = 0.5. 

In the text 
Fig. 10
Lightcurves calculated for a planet receiving an average insolation of 1021 Wm^{2} (T_{eq} = 366 K with A = 0.1) but orbiting different stars. All parameters are set to default values except the stellar mass and orbital period. The flux from the planet received at 10 pc is given in the top panel, while the bottom one shows the planet/star contrast ratio. Corresponding orbital periods P and semimajor axes a are indicated. 

In the text 
Fig. 11
lightcurves for low thermal inertia halfdark halfbright planet. Dashed lines show the lightcurves for a planet with the default parameters but a low thermal inertia (100 SI) and albedo values between 0.1 and 0.9. The solid curve is obtained if the surface of the planet is divided into two hemispheres (separated by two meridians) with highly different albedos: A = 0.1 and 0.9. On the top panel, the system is in a 2:1 spinorbit resonance. The lightcurve is periodic but its shape depends on the arbitrary position of the hemispheres. On the bottom panel, the system is in pseudosynchronization and the signal is no longer periodic. 

In the text 
Fig. 12
Influence of the obliquity β. Parameters are set to their default values and ω = ω_{eq} = 2.38n. Panels a) and b): lightcurves for different obliquities and two orientations of the spin axis, without tidal heating. Panels c) and d): Same as b) but with tidal heating included, uniformly in c) and with half of the tidal dissipation released through one hot spot (Gaussian, half width: 45°) centered on the equator in d). 

In the text 
Fig. A.1
The eye of Sauron. Temperature maps showing the shrinking of the permanently dark area with increasing eccentricities for a synchronously rotating planet. Except for e, default parameters are used. 

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.