arxiv: v2 [astro-ph.ep] 2 Oct 2018

Transcription

1 Exoplanet phase curves: observations and theory Vivien Parmentier and Ian J. M. Crossfield arxiv: v2 [astro-ph.ep] 2 Oct 2018 Abstract Phase curves are the best technique to probe the three dimensional structure of exoplanets atmospheres. In this chapter we first review current exoplanets phase curve observations and the particular challenges they face. We then describe the different physical mechanisms shaping the atmospheric phase curves of highly irradiated tidally locked exoplanets. Finally, we discuss the potential for future missions to further advance our understanding of these new worlds. Observing exoplanetary phase curves Origin and shape of a phase curve Planetary atmospheres are intrinsically three-dimensional objects, with both smalland large-scale variations of temperature, chemistry, and cloud coverage. This is even more important for the current population of characterizable exo-atmospheres, as most of them belong to planets in close-in orbits, probably tidally locked, with a large day/night temperature contrast induced by permanent radiative forcing. When ignored, the large spatial inhomogeneities of these planets can lead to a biased interpretation of transiting and secondary eclipse observations (Line and Parmentier 2016; Feng et al. 2016). Observing the phase curve of an exoplanet (i.e., the time-dependent change in the brightness of a planet as seen from Earth during one orbital period) is the most Vivien Parmentier University of Arizona, Lunar and Planetary Laboratory, 1629 E University Blvd, Tucson, AZ, USA vivien@lpl.arizona.edu Ian J. M. Crossfield Department of Physics, Massachusetts Institute of Technology, Cambridge, MA, USA iancross@mit.edu 1

2 2 Vivien Parmentier and Ian J. M. Crossfield straightforward way to probe the planet s longitudinal structure. The brightness of the planet is determined by the combined emitted and reflected light in the particular bandpass of the observations. For transiting planets, the shape and amplitude of the phase curve are determined by longitudinal inhomogeneities. Inhomogeneous illumination is always present, as over the course of one orbit, we see different hemispheres of the planet, ranging from its dayside (before and after the eclipse of the planet by the star) to its nightside. Inhomogeneities in temperatures can be probed by the thermal emission of the planet, with a nightside that is usually colder and thus dimmer than the dayside. Finally, inhomogeneous chemical composition and cloud coverage can be determined through the phase curve s wavelength dependence and shape. Short-period planets are good targets for atmospheric characterization through phase curve observations: their short orbit makes it easier to monitor the system during a full orbit of the planet. Moreover, planets in orbits 10 days are likely tidally locked (Guillot et al. 1996; Parmentier et al. 2015), meaning that their rotation period is their orbital period and that the longitude of the Earth-facing hemisphere can be related to the phase of the orbit. Finally, tidally locked planets in short orbit have a short radiative timescale, a weak Coriolis force, and an input of energy at large scale, making them likely to have planetary-scale atmospheric features that can be observed in the hemispherically averaged planetary flux (Showman and Guillot 2002). From phase curves to atmospheric properties For a steady-state atmosphere, a phase curve provides the longitudinal variation of the hemispherically averaged brightness of the planet. As such, only large-scale atmospheric structures can be inferred (see chapter by Cowan and Fuji). From current phase curve observations, three parameters can usually be retrieved (e.g. Cowan and Agol 2008; Demory et al. 2013; Knutson et al. 2009, 2012). For transiting planets, these parameters are the secondary eclipse depth, the phase of the maximum of the phase curve compared to the secondary eclipse (called the phase curve offset), and the phase curve relative amplitude, A F = (F p,max F p,min )/F p,max (see Fig. 1). The secondary eclipse depth gives the brightness of the dayside hemisphere (see chapter by Roi Alonso) and serves as a reference to calculate A F. The phase offset gives the longitude of the brightest hemisphere of the planet. A planet with a maximum brightness at the substellar point and a symmetrical brightness distribution will have a phase curve peaking during the secondary eclipse. A planet that is brighter east of the substellar point will have a phase curve peaking before secondary eclipse, leading, by definition, to a positive offset of the phase curve maximum. A planet that is brighter west of the substellar point will have a phase curve peaking after secondary eclipse, leading to a negative offset. Here we define east and west with respect to the rotation of the planet. For tidally locked planets,

3 Exoplanet phase curves: observations and theory 3 Relative Flux Eclipse Eclipse Transit Offset Relative Flux Min Max Eclipse depth Orbital Phase Fig. 1 Phase curve of HD189733b observed at 3.6µm with the Spitzer Space Telescope by Knutson et al. (2012). The bottom panel is a zoomed-in version of the top panel, with the quantities of interest annotated. the rotation and the revolution have the same direction and the link between phase curve offset and brightness distribution is independent of the observer s position. The phase curve relative amplitude A F provides insight into the brightness contrast between the brightest and the dimmest hemisphere. It goes from 0 (corresponding to no brightness variation) to 1 (when one hemisphere emits zero flux). Different flavors of this parameter can be used. As an example, Perez-Becker and Showman (2013) defined the day/night relative amplitude (F p,day F p,night )/F p,day ), whereas Komacek and Showman (2016) used the day/night brightness temperature contrast. The day/night contrast is useful when studying the energy balance of the atmosphere (Cowan et al. 2012; Schwartz and Cowan 2015), but its estimated value can depend on the measurement of the phase curve offset (Schwartz et al. 2017). The temperature contrast might seem a more intuitive quantity, however, for phase curves sparsely sampling the spectral space, the conversion from flux to temperatures is model dependent (Cowan et al. 2012). For massive planets, gravitational interactions between the planet and the star can lead to non-atmospheric signals in the phase curve such as the ellipsoidal variations and the beaming effect (Shporer 2017). Although techniques can be used to

4 4 Vivien Parmentier and Ian J. M. Crossfield isolate the atmospheric component of the phase curve, it inevitably leads to higher uncertainties in the derived atmospheric properties (Shporer and Hu 2015). Observational Challenges Phase curves are a uniquely challenging phenomenon to observe. The timescale and characteristic shape of phase curves make them much more difficult to observe than occultations. While the magnitude of a phase curve signal is comparable to that of the more commonly observed transits or eclipse, the phase curve timescale days is much longer than the hour timescale of occultations. Furthermore, phase curves contain only low-frequency components (Cowan and Agol 2008), while the high-frequency component is much more prominent in transits and eclipses. Occultation observations can be safely decorrelated even against high-order polynomials or complicated and/or periodic signals (e.g., Haynes et al. 2015), which is much more challenging for phase curves. In contrast, a single phase curve observation can be compromised by just a single monotonic trend. For example, it took five years to recognize that 40% of the originally published Spitzer/IRAC 8 µm phase curve of HD b was compromised by systematics (Knutson et al. 2007, 2012). Similarly, the Spitzer/MIPS 24 µm phase curve of HD b was irretrievably compromised by an instrumental drift lasting tens of hours (Crossfield et al. 2012). In such analyses, as with the case of the apparent ellipsoidal variation seen around WASP-12b at 4.5 µm (Cowan et al. 2012; Stevenson et al. 2014a), the question is often: how does one distinguish signal from systematic? The most straightforward solution is probably also the most effective: repeat the experiment, and observe multiple phase curves at the same wavelengths. This is the approach adopted in the first HST/WFC µm spectroscopic phase curve, of WASP-43b (Stevenson et al. 2014b). By observing three phase curves, a simultaneous analysis allowed the phase curve to be fit while separating out a long quadratic systematic in addition to the standard HST effects. While multiple observations with the same instrument do not guarantee repeatability (Stevenson et al. 2017), this approach at least ensures that questionable results are identified. Beyond observing multiple phase curves, a more fundamental requirement is to at least observe a full phase curve i.e., begin and end with a secondary eclipse whenever possible. Implicit in the above comments is the assumption of a continuously observed phase curve. The first phase curve observations employed a more parsimonious, snapshot approach using as few as five discrete-observing epochs (Harrington et al. 2006). This analysis was followed by a combined, snapshot plus continuous, program that measured a mid-infrared phase curve and confirmed the long-term stability of the MIPS 24 µm detector in most cases (Crossfield et al. 2010). Later, with the benefit of over a decade of Spitzer analyses, Krick et al. (2016) demonstrated the ability of Spitzer/IRAC to also obtain snapshot-mode phase curves. Caution may be advised since the phase curve amplitudes derived from the IRAC snapshot and from

5 Exoplanet phase curves: observations and theory 5 continuous observations (Wong et al. 2015) differ by 5σ. However, the snapshot technique is not obviously flawed even multiple continuous phase curves of the same target sometimes disagree by up to 4σ (Stevenson et al. 2017). Regardless, the snapshot technique is unlikely to dominate in the JWST era owing to the large overheads associated with even modest telescope slews. Wavelength dependence Observing exoplanet phase curves at different wavelengths is essential to obtain a complete view of their atmosphere. At optical wavelengths, where the brightness of the planet is often dominated by reflected light, the phase curve provides information on the longitudinal variation of the planet s albedo. At infrared wavelengths, where the brightness of the planet is often dominated by thermal emission, the phase curve provides information on the longitudinal variation of the planet s temperature and chemical composition. The spectral variation of thermal phase curves is shaped by molecular features. Inside a molecular absorption band, a phase curve probes low pressures (high altitudes) while probing deeper outside a molecular absorption band (Showman et al. 2009; Kataria et al. 2015). As a consequence, multiwavelength phase curves probe 2D (longitude, depth) thermal and chemical structure of the atmosphere (Knutson et al. 2009; Stevenson et al. 2014b, 2017). In general, there is no reason to believe that phase curves obtained in different bandpasses will be similar. When looking for trends between the amplitude or offset of phase curves and planetary parameters, it is therefore important to compare observations taken in the same bandpasses. In some bandpasses, both reflected light and thermal emission can contribute significantly. Fig. 2 shows the theoretical reflected and thermal phase curves of a hot Jupiter in the Kepler bandpass ( nm). The reflected phase curve is dominated by the presence of bright clouds west of the substellar point, whereas the thermal emission phase curve is dominated by the presence of a temperature maximum east of the substellar point. For this specific model (T eq = 1900K) and bandpass, thermal emission and reflected light compensate, and no offset is measured. For cooler planets, reflected light should dominate over thermal emission, whereas for hotter planets it should be the opposite, leading to the correlation between phase curve offset and equilibrium temperature for Kepler phase curves seen in the right panel of Fig. 2. Summary of observations to date Table 1 lists the current phase curve observations of transiting exoplanets, while Fig. 3 shows the observations for the thermal phase curves of hot Jupiters. In general, hot Jupiters thermal phase curves have a large amplitude (between 0.5 and 1)

6 6 Vivien Parmentier and Ian J. M. Crossfield Table 1 Current exoplanets phase curve observations Planet T e eq(k) P (days) λ (µm) A F Offset ( ) Reference 55 CnC e ± ± 12 Demory et al. (2016) GJ436b g ± ± 26 Stevenson et al. (2012) HAT-P-7b ± 0.04 c 7.9 ± 0.3 a Esteves et al. (2015) ± ± 7.5 Wong et al. (2016) ± ± 7.5 Wong et al. (2016) HD149026b ± 0.34 n/a Knutson et al. (2009) HD189733b ± ± 4 Knutson et al. (2012) ± ± 5.5 Knutson et al. (2012) 8 > ± 2.7 Knutson et al. (2012) ± ± 8.1 Knutson et al. (2012) HD209458b ± ± 6 Zellem et al. (2014) υ And b f n/a 84.5 ± 6.3 Crossfield et al. (2010) WASP-12b ± 0.19 b 53 ± 7 d Cowan et al. (2012) ± 0.12 b 16 ± 4 Cowan et al. (2012) WASP-14b > ± 1.4 Wong et al. (2015) ± ± 1.4 Wong et al. (2015) WASP-18b > ±0.5 Arcangeli et al. (2017) ± ± 1 Maxted et al. (2013) b ± ± 2.3 Maxted et al. (2013) b WASP-19b ± ± 4 Wong et al. (2016) ± ± 3.6 Wong et al. (2016) WASP-43b ± ± 1.0 Stevenson et al. (2014b) ± ± 1.7 Stevenson et al. (2017) ± ± 1.8 Stevenson et al. (2017) WASP-103b ± ± 0.1 Kreidberg et al. (2018) ± ± 0.7 Kreidberg et al. (2018) ± ± 0.4 Kreidberg et al. (2018) Kepler-7b ± 0.38 c 25 ± 1.9 Esteves et al. (2015) Kepler-8b ± 0.45 c 24.8 ± 5.6 Esteves et al. (2015) Kepler-10b ± ± 6 Hu et al. (2015) Kepler-12b ± 0.26 c 44.3 ± 5.8 Esteves et al. (2015) Kepler-41b ± 0.48 c 31 ± 5 Esteves et al. (2015) Kepler-76b ± 0.05 c 8.3 ± 1.2 Esteves et al. (2015) For optical phase curves, only planets for which an offset is detected have been considered. The phase curves of highly eccentric planets are discussed in a later section. Only planets on circular orbits have been taken into account here, a discussion of eccentric planets can be found in the devoted section. Instruments used are Spitzer for the 3.6, 4.5, 8 and 24µm observations, HST/WFC3-IR for the 1.5µm observations, Kepler for the 0.6µm observations a The offset varies from orbit to orbit (Armstrong et al. 2016) b A F values were taken from Stevenson et al. (2017) c Uncertainty for phase curve maximum observed in the Kepler bandpass was assumed to be similar to the uncertainty of the secondary eclipse. d Using a different noise model the same authors get an offset of 0 ± 29. e We use T 4, a : semi-major axis, T star, R star : stellar effective temperature and radius. f This is a non transiting planet with a measured inclination of 24 ± 4 (Piskorz et al. 2017). g This planet has a non-zero but small eccentricity (e=0.1371) (see Stevenson et al. 2010). eq = T star 4 R2 star 4a 2

7 Exoplanet phase curves: observations and theory 7 Fp/F * [ppm] Total Reflected Emitted Phase ( ) Fig. 2 Left : Phase curve in the Kepler bandpass calculated from a hot Jupiter model with T eq = 1900K including silicate clouds. The phase curve (black) is a combination of the reflected phase curve (blue) and the thermal emission phase curve (red). The circles are brightness maps of the Earth-facing hemisphere. The bright strip on the east of the substellar point is reflected light due to clouds whereas the rest is thermal emission. Right : Observed offsets in the Kepler phase curve sample. Planets with a small equilibrium temperature, dominated by reflected light, have a negative offset, whereas planets with a high equilibrium temperature, dominated by thermal emission, have a positive offset. Figure from Parmentier et al. (2016). and have a positive offset, meaning that the brightest hemisphere is eastward of the substellar point. As seen in Fig. 3, there is no clear trend between the amplitude of thermal phase curves and planet equilibrium temperature. Earlier claims by Perez-Becker and Showman (2013) and Komacek and Showman (2016) were based on a smaller number of observations and by interpreting together observations taken in different bandpasses. In today s more complete dataset, no trend is seen, neither by looking at all bandpasses together nor by looking at them separately. A tentative trend, first proposed by Stevenson et al. (2017), is seen in the amplitude vs. planet rotation period plot: planets with a faster rotation rate might have a larger phase curve amplitude. Observations of planets with an orbital period between 1 and 2 days and, especially, larger than 3 days are currently being taken to confirm the presence of a correlation (Spitzer program 13038). The most striking feature of Fig. 3 might be in the phase curve offset vs. equilibrium temperature plot, with a lack of large offset for planets with equilibrium temperatures 1700 K. WASP-12b seems an outlier, but it s large offset might be the result of uncorrected instrument systematics (Cowan et al. 2012). The super- Earth 55 Cnc e does not fit in this trend defined by hot Jupiter observations. Its phase curve appears similar to cooler hot Jupiters, such as HD209458b, which has a large offset and a large phase curve amplitude at 4.5µm, pointing toward unique atmospheric or ground properties.

8 8 Vivien Parmentier and Ian J. M. Crossfield Fig. 3 Amplitude (top) and offset (bottom) of exoplanet thermal phase curves as a function of equilibrium temperature (left) and orbital period (right). Different symbols and colors represent different wavelengths. Values predicted at 4.5µm by the Jupiter radius, solar composition, chemical equilibrium, cloudless and drag-free global circulation model of Parmentier et al. (2016) are overplotted. The model systematically underpredicts the amplitude and overpredicts the offset of the phase curve. For the 0.6µm case, only planets with a phase curve dominated by thermal emission are considered (i.e. planets with T eq > 2000K). For WASP-12b, the two possible offsets derived from the same phase curve are linked with a dashed curve (Cowan et al. 2012). Colored points are hot Jupiters, gray points are super-earths and mini-neptunes. The data is from Table 1. The physics of exoplanet phase curves Thermal structure Atmospheric regimes In response to the large day-to-night irradiation gradient, tidally locked planets develop a strong atmospheric circulation that advects energy from the dayside to the nightside and reduces the day/night temperature contrast compared to the local radiative equilibrium case. However, energy losses by radiation or by the damping of

9 Exoplanet phase curves: observations and theory 9 the winds can prevent a fully efficient redistribution of energy. This balance between energy transport and energy losses can be understood in terms of timescales (Showman and Guillot 2002; Komacek and Showman 2016). The radiative timescale is the characteristic time it takes for a parcel of gas to lose its energy by radiation. The wave timescale is the time it takes for a gravity wave to travel horizontally over one planetary radius. The drag timescale is the time it takes for the waves or the winds to lose a significant part of their kinetic energy. Depending on the balance between these three timescales, the atmospheric circulation on tidally locked planets is expected to follow at least two different regimes (see also chapter by N. Lewis). When the radiative timescale and the drag timescales are similar or longer than the wave timescale, the presence of an eastward, equatorial jet is expected. The resulting temperature map is shifted eastward compared to a local radiative equilibrium case, and the hottest spot of the atmosphere lays eastward of the substellar point. When the radiative timescale or the drag timescale is smaller than the wave timescale, the atmospheric circulation is characterized by a day-to-night flow with a temperature map symmetric around the substellar point and a somewhat larger hemispheric temperature contrast (e.g. Showman and Polvani 2011; Tsai et al. 2014). Phenomenological models vs. global circulation models The energy balance of a planet can be determined by measuring the thermal emission from two opposite hemispheres. When considering the day- and nightside hemispheres, one can understand the energy balance in terms of how much of the incoming energy is reflected to space (i.e., the Bond albedo) and how much of the absorbed stellar light is transferred to the nightside (i.e., the redistribution factor). In order to have a complete energy balance, the full dayside and nightside spectra should be used. However, current observations only cover a few bandpasses, and the effective temperature of each hemisphere must be estimated by extrapolating the measurements with assumptions that are only exact if each hemisphere emits like a blackbody. These approximations should work for the planet s dayside (Cowan and Agol 2011a) but might fail when estimating the nightside temperatures (Schwartz and Cowan 2015) where the vertical temperature gradient is larger. The latest results from Schwartz et al. (2017) show a large dispersion in both redistribution efficiency and geometric albedo between planets, with no apparent trend between equilibrium temperature and day/night temperature contrast, echoing the lack of trends seen in the top left panel of Fig. 3. Schwartz et al. (2017) also conclude that many hot Jupiters have Bond albedos larger than 0.3, in apparent contradiction with their measured low geometric albedos in the Kepler bandpass (Heng and Demory 2013). This apparent discrepancy could be the result of an asymmetric scattering function (e.g., Dyudina et al. 2005), a lower geometric albedo in the Kepler bandpass than outside of it (Crossfield 2015), or an intrinsic bias in the method. Linking the phase curve amplitude and the phase offset is a more complex task that requires a model of the longitudinal distribution of the temperature and the

10 10 Vivien Parmentier and Ian J. M. Crossfield opacities. One (and two) -dimensional models of the longitudinal (and latitudinal) variation of the temperature have been calculated by taking into account the competing effects of longitudinal advection of energy and radiative losses (Cowan and Agol 2011b; Hu et al. 2015; Zhang and Showman 2017). They lead to a temperature map that is determined by one parameter: the ratio of the advective and the radiative timescale. Consequently, they always predict a correlation between the phase curve offset and the phase curve amplitude. These models, however, lack vertical transport, which is proven to be an important factor setting the day/night temperature contrast (Komacek and Showman 2016). The phase curve offset and the phase curve amplitude might therefore be set by different mechanisms, which could be the reason why they are not correlated in the observations (Crossfield 2015) nor in the more complex, three-dimensional models of tidally locked exoplanets (Komacek et al. 2017; Parmentier et al. 2016). Whereas these phenomenological models are useful to retrieve parameters such as the radiative timescale from phase curve observations, they can t say much about the mechanisms setting these parameters. Global circulation models, solving for the hydrodynamics, the radiative transfer, and/or the magnetic effects in three dimensions, have been used to quantitatively link the observed day/night contrast and phase curve offset to planetary parameters. As seen in Fig. 3, cloudless, dragless, solar composition models systematically underpredict the amplitude and overpredict the offset of thermal phase curves (Showman et al. 2009; Kataria et al. 2015). Other mechanisms, such as metallicity (Kataria et al. 2015), cloud composition (Oreshenko et al. 2016; Parmentier et al. 2016), rotation period (Showman et al. 2009; Showman and Polvani 2011), disequilibrium chemistry (Cooper and Showman 2006), or magnetic field strength (Rogers 2017) must play an important role. Radiative timescale The radiative timescale is expected to vary linearly with pressure and with the inverse of the cube of the temperature (Iro et al. 2005; Showman et al. 2008). Since the pressures probed by atmospheric sensing span several orders of magnitude, the pressure dependence of the radiative timescale has a major consequence on the observations. Phase curves obtained at wavelengths probing deep in the atmosphere are expected to probe layers with a large radiative timescale and thus a small amplitude and a large offset, whereas phase curves obtained at wavelengths probing shallower layers (e.g., inside molecular bands) should have a larger amplitude and smaller offset. This is illustrated in the cloudless case of Fig. 6. The temperature dependence of the radiative timescale is also expected to impact the observations. With all things being equal, hotter planets are expected to cool more efficiently, leading to a larger day/night temperature contrast and thus a larger phase curve amplitude (Perez-Becker and Showman 2013; Komacek and Showman 2016). As seen in Fig. 3, the lack of a clear trend in the amplitude vs. equilibrium

11 Exoplanet phase curves: observations and theory 11 temperature plot indicates that other mechanisms must contribute to the shape of exoplanets phase curves (Komacek et al. 2017). Drag, composition and mean molecular weight The thermal structure of the planet can be affected when either the radiative timescale, the drag timescale, or the wave timescales are changed. Atmospheric drag is often used and parametrized as Rayleigh drag in models to represent different types of physical mechanisms such as shocks (Perna et al. 2012; Heng 2012), hydrodynamic instabilities (Fromang et al. 2016) or Lorentz forces (Perna et al. 2010; Batygin et al. 2013). Its main effect is to weaken the superrotating jet, leading to a more pronounced substellar to anti-stellar flow and a weaker equatorial eastward jet. As seen in Fig. 4, small Rayleigh drag (τ drag 10 5 s) has a small effect on the phase curve amplitude but a large effect on the phase offset. By shifting the circulation pattern from a jet-dominated to a day-to-night flow, a small drag does not change much the heat transport to the nightside but reduces the longitudinal asymmetry of the temperature and therefore the phase curve offset. A larger drag (e.g., τ drag s) can increase the phase curve amplitude significantly, but lead to an extremely small phase curve offset. It is important to keep in mind that the physical mechanism responsible for the drag might not be homogeneous in the atmosphere or might not even act as a Rayleigh drag at all. For example, Lorenz forces should have a fundamentally different effect on the flow than Rayleigh drag, leading to a potential inversion of the equatorial jet speed (Rogers and Komacek 2014; Rogers 2017), whereas hydrodynamic instabilities and shocks should drag the flow where large velocity or temperature gradients are present (Fromang et al. 2016; Heng 2012). Atmospheric composition can also play an important role in shaping the phase curve of exoplanets. By increasing the abundance of metals in the atmosphere, metallicity increases the opacities, leading to two consequences. First the atmospheric layers probed are shifted to smaller pressures, where the radiative timescale is smaller leading to a large phase curve amplitude and a smaller phase curve offset. Second the enhanced opacities lead to a warmer atmosphere and thus to a smaller radiative timescale. Overall, increasing the atmospheric metallicity should increase the amplitude and decrease the offset of the phase curve. For WASP-43b, global circulation models with a metallicity of five times solar provide a better, but not satisfactory, match to the observations compared to models assuming a solar metallicity (Kataria et al. 2015). In planets with equilibrium temperatures larger than 1900K, titanium and vanadium oxide are expected to be present in the dayside atmosphere and create a strong thermal inversion (Hubeny et al. 2003; Fortney et al. 2008; Parmentier et al. 2015, 2016). The detection of such an inverted temperature profile has recently been claimed through the presence of molecular emission features in the dayside of both WASP-33b (Haynes et al. 2015) and WASP-121b (Evans et al. 2017). Thermal in-

12 12 Vivien Parmentier and Ian J. M. Crossfield versions are expected to disappear in the planet s nightside where these molecular features are expected to become absorption features. In small planets, where the atmospheric composition can be very diverse, the mean molecular weight of the atmosphere is expected to impact both the radiative timescale and the wave timescale. Zhang and Showman (2017) show that the day/night temperature contrast is expected to increase, and the eastward shift of the hot spot is expected to decrease, when the mean molecular weight is increased. In real planets, atmospheres with different mean molecular weight will have a very different composition, leading to variations of several orders of magnitude in the opacities and potentially a comparable or stronger effect than the effect of the mean molecular weight alone. Fig. 4 Observed thermal phase curve amplitude (upper panel) and offset of the maximum (bottom panel) compared with semi-gray global circulation models from (Komacek et al. 2017) using different strength of drag. The color and shape of the points are the same as in Fig. 3.

13 Exoplanet phase curves: observations and theory 13 Opacity structure Chemistry In chemical equilibrium, temperature inhomogeneities often pair up with chemical composition inhomogeneities. For a given atomic composition, equilibrium favors different molecules at different pressures and temperatures. A day/night temperature contrast on a tidally locked planet could result in a day/night chemical gradient. Such a large-scale change in the chemistry should trigger a change in the opacities and affect the offset and the amplitude of the phase curve. Horizontal and vertical advection are expected to further complicate this picture. When the chemical reaction timescale is long compared to the horizontal and the vertical mixing timescale, chemical reactions cannot happen fast enough during the transport of gas from one side to the other side of the planet and quenching happens. The atmospheric circulation can then drive the chemical abundances out of their local equilibrium state and erase any chemical gradient expected from the local equilibrium (Cooper and Showman 2006; Visscher and Moses 2011). The exact chemical composition of the atmosphere is determined by a combination of the chemical equilibrium abundances of the deep atmosphere (vertical quenching) and of the hot dayside (horizontal quenching) (Agúndez et al. 2014). Chemical quenching is expected to affect the opacities in very specific bands, leading to a peculiar signature in the day/night contrast vs. wavelength relationship (Steinrueck et al. 2017). In very hot planets (T eq 2500K) molecules such as water should be thermally dissociated at the dayside photosphere but not in the nightside, whereas more strongly bond molecules, such as CO, should be present at all longitudes (see Fig. 5). When molecules dissociate, continuum opacities from hydrogen ions are expected to become dominant, leading to a lack of spectral features and a blackbody-like thermal emission in the dayside (Kreidberg et al. 2018; Arcangeli et al. 2017). Clouds The large temperature variations in the atmosphere of tidally locked planets in closein orbits are responsible for the prevalence of large longitudinal inhomogeneities in their cloud coverage as inferred from the observations (Demory et al. 2013; Shporer and Hu 2015). In each planet, some species can condense in the cold nightside of the planet but must be in gaseous form in the hot dayside (Parmentier et al. 2013). Since the temperature map of most hot Jupiters is not symmetric, but shifted eastward, the western part of the dayside is usually cold and clouds can also form there. The reflected phase curve can then be dominated by the bright reflective part of the dayside atmosphere where cloud are present, leading to a phase curve peaking after secondary eclipse, the opposite sign than for thermal phase curve (see Hu et al. (2015); Garcia Munoz and Isaak (2015); Webber et al. (2015) and Fig 2). The longitudinal distribution of a given type of cloud in the dayside is determined to first order by the thermal structure of the planet and the cloud-specific condensation

14 14 Vivien Parmentier and Ian J. M. Crossfield α=-180 α=-120 α=-60 α=0 α=60 α=120 α=180 Temp Temperature CO H 2 O H log 10 (Abundance) Fig. 5 Temperature (top) and abundances of water, CO and H at the 1.4µm photosphere of WASP-103b as predicted from global circulation models of Parmentier et al. (2018). The different columns are for different orbital phases, α = 0 being the secondary eclipse. The water abundance drops by two orders of magnitude in the dayside due to thermal dissociation (Kreidberg et al. 2018). The abundances are calculated assuming local chemical equilibrium (Visscher et al. 2006). temperature (Lee et al. 2016). If the thermal structure is known, either by observing a thermal phase curve or by modeling the temperature distribution a priori, the cloud map derived from the reflected light curve can be used to constrain the cloud chemical composition (Oreshenko et al. 2016; Parmentier et al. 2016). Clouds should also affect the thermal emission of the planet. The presence of nightside clouds produces a large opacity gradient between the day and the night that can suppress the thermal emission from the cloudy regions by raising the photosphere to low pressures. As a consequence, clouds are expected to increase the phase curve amplitude and decrease the phase curve offset of hot Jupiters (see Fig. 6), even if they are present only on the planet s nightside. When clouds are present, the brightest hemisphere is not necessarily the hottest one, and the shift of the maximum of thermal phase curves does not track the shift of the hottest point of the planet anymore (Parmentier et al. 2017). Similarly, the day/night flux contrast cannot be easily converted to a day/night horizontal temperature contrast, a redistribution parameter or an advective timescale. Any inference of the thermal structure and atmospheric composition of an atmosphere based on the phase offset and the day/night contrast of a thermal phase curve can be highly biased by the presence of clouds. An optical phase curve, probing the planet albedo as a function of longitude seems a necessary complement to any thermal phase curve.

15 Exoplanet phase curves: observations and theory 15 Fig. 6 Relative amplitude (left) and phase offset (right) of currently known exoplanet phase curves as a function of wavelength. The color of the points represents the equilibrium temperature of the planet. Predictions from three different global circulation models of a hot Jupiter with T eq = 1400K have been overplotted (Parmentier et al. 2017): an example of a cloudless model (orange), a model with a partial cloud coverage (green, assuming the presence of MnS clouds) and an example of a fully cloudy model (red, assuming the presence of MgSiO 3 ). Time variability Atmospheres can be time variable, leading to a time dependence of the phase curve amplitude and offset between orbits. Given the large number of observations needed to observe one phase curve and thus the time needed to constrain differences between phase curves, atmospheric variability has not yet been widely characterized. The most notable exception is the planet HAT-P-7b, a hot Jupiter that has been observed continuously by the Kepler spacecraft over more than 700 orbital periods. As shown by Armstrong et al. (2016), the shift of the maximum of the phase curve may vary from positive to negative values over a timescale of 50 orbits. Two scenarios have been proposed to explain this behavior. In the first scenario, a global change in the dayside temperature of the planet, without necessarily a change in the position of the hot spot, would lead to a variation in the cloud coverage, changing the relative ratio of thermal and reflected contribution to the light curve (Armstrong et al. 2016). In the second scenario, proposed by Rogers (2017), the entire phase curve of HAT-P-7b is dominated by thermal emission. The coupling between the ionized atmosphere and the magnetic field leads to an oscillation in the strength and direction of the equatorial jet and the center of the hottest hemisphere itself oscillates from east to west of the substellar point. Observations of variability at infrared wavelength should be able to disentangle between the two scenarios. In the first case, no variation is expected at infrared wavelengths, where reflected light is never dominant. In the second scenario, the change of the hottest hemisphere longitude should affect similarly the Kepler phase curve and the infrared phase curve.

16 16 Vivien Parmentier and Ian J. M. Crossfield Eccentric planets Highly eccentric planets are laboratories to understand the response of a planet atmosphere to a time-dependent forcing. Quantities that are degenerate in the phase curves of tidally locked planets, such as the radiative timescale and the advective timescale, can potentially be measured separately in the eccentric case (Lewis et al. 2013; de Wit et al. 2016). When passing through periastron, the planet is inhomogeneously and transiently heated. Then, depending of the rotation period of the planet (Hut 1981), the heated hemisphere is expected to face Earth periodically after periastron, leading observable oscillations (or ringing ) in the phase curve (Langton and Laughlin 2008; Kataria et al. 2013). While the period of these oscillations are related to the rotation period of the planet, their amplitude is linked to the radiative timescale of the atmosphere (Cowan and Agol 2011b). So far, the thermal phase curves of two highly eccentric planets have been published, while the phase curve of a third one (XO-3b) has been observed but is not yet published. No oscillations were found in the thermal phase curves of HAT-P-2b and HD80606b, which was interpreted by the presence of a short ( 2 6h) radiative timescale at the photosphere (Lewis et al. 2013; de Wit et al. 2016). Such a short radiative timescale is reproduced by global circulation models of HAT-P-2b (Langton and Laughlin 2008; Lewis et al. 2013) but not by models of HD80606b, where the cooler temperatures predict a longer timescale. Any opacity source, such as clouds, could raise the photosphere of HD80606b to lower pressures and account for this short radiative timescale (Lewis et al. 2017). Planned Hubble Space Telescope observations should measure the albedo of HD80606b and its variation around periastron. Non-transiting planets Phase curve observations can be used to characterize the atmospheres of nontransiting planets on short orbits. The number of non-transiting planets is far greater than the number of transiting ones, so non-transiting planets can potentially provide a large sample of bright planets to study (Millholland and Laughlin 2017). Unfortunately, the phase curve relative amplitude cannot be determined directly for nontransiting, even when the inclination of the system is known, as both the radius of the planet and the stellar absolute flux cannot be measured directly (Crossfield et al. 2010). The phase curve offset, however, can be directly determined from the observation, given a good enough knowledge of the planet ephemeris. So far, υ And b is the only non-transiting planet with a known inclination (Piskorz et al. 2017) that has a measured phase curve offset (Crossfield et al. 2010). Measuring a phase curve offset for planets over a wide range of inclinations could be used to infer the latitudinal thermal structure of exoplanets atmospheres.

17 Exoplanet phase curves: observations and theory 17 Future prospects with JWST The James Webb Space Telescope (JWST) will provide valuable information on giant tidally locked planets through phase curve observations. By observing both optical infrared wavelengths, NIRISS should be able, in one spectroscopic phase curve observation, to provide both the cloud map and the temperature map of hot Jupiters and show the change of sign of the phase curve offset with wavelength expected for a partially cloudy atmosphere (see Fig. 6). NIRCAM and NIRSPEC will provide valuable measurements of molecular abundances and their variation with longitude. These near-infrared wavelengths, however, are not best suited to observe the cooler nightside of tidally locked planets where most observations measured a close-to-zero flux and could not detect any molecular features. These cooler nightsides should, however, be brighter at longer wavelengths, where MIRI can observe. CHEOPS, TESS, and PLATO will significantly increase the number of optical phase curves of hot Jupiters, giving important constraint on the cloud distribution on exoplanets (see Shporer 2017, for a complete review). Moreover, these missions will find numerous planets, particularly in the Neptune size range, that current and future observatories (HST, JWST) will be able to characterize. As such, numerous phase curves of Neptune and sub-neptune size planets could be obtained during the next decade, increasing drastically the current sample of two (55 CnC e and GJ436b, see Figure 3). FINESSE and ARIEL, currently into the selection process, would provide hundreds of exoplanets spectroscopic phase curves covering a wide range of planet mass, radius, irradiation and rotation period. Such a survey mission is needed to understand the diversity of climates on extrasolar worlds. JWST will probably be used to search for the phase curves of temperate rocky planets orbiting M dwarfs (Kreidberg and Loeb 2016), such as the Trappist-1 planets, Proxima centauri b or LHS 1140b. Theoretically, the phase curve amplitude of these systems can be used to determine whether the planet possesses a fluid that can efficiently transport energy (i.e., an atmosphere or an ocean), since the amplitude of the phase curve of a bare rock would be larger than any other cases. In practice, these observations will be extremely challenging as they require an exquisite stability of the telescope and the star over a period of days to tens of days. Future and more precise observatories, such as the Origins Space Telescope in the infrared or LUVOIR in the optical, will be more likely to obtain phase curve measurements of temperate planets. Interpreting further the shape of these temperate rocky planet phase curves will be complex, as the number of unconstrained parameters that can influence the phase curve of a rocky planet is larger than the ones influencing a hot Jupiter, for which we already face many challenges. Among them are the large range of possible atmospheric composition (Leconte et al. 2015a), the possibility of asynchronous rotation (Leconte et al. 2015b), the range of possible atmospheric mass (Koll and Abbot 2015), the thermal inertia of the ground (Selsis et al. 2013), the peculiar behavior of clouds (Yang et al. 2013), the effects of non-dilute atmospheres (Pierrehumbert

18 18 Vivien Parmentier and Ian J. M. Crossfield and Ding 2016) or the intertwined signal of different planets (including undetected ones) on the same system (Kane and Gelino 2013). Conclusion Exoplanet phase curves are fundamental probes of the three-dimensional structure of exoplanets in close-in orbits. The current sample shows that every exoplanet studied so far is unique, making the search for trends between phase curve shape and planetary parameters harder. State-of-the-art models are unable to quantitatively explain current phase curve observations. Whether clouds, drag, atmospheric variability, chemistry, or instrumental systematics are behind the unexpectedly large amplitude and small offsets of hot Jupiters, thermal phase curves should be determined in the near future through the wavelength dependence of the phase curve shape as will be observed by coming spacecrafts. Numerous observations of cooler and smaller planets are expected in the next decade and will probably challenge our current understanding of planetary atmospheres. Cross-References Observing Exoplanets with the Spitzer Space Telescope Characterization of Exoplanets: Secondary Eclipses Mapping Exoplanets Exoplanet Atmosphere Observations from Transmission Spectroscopy and Other Planet-Star Combined Light Observational Techniques Atmospheric Circulation for Exoplanet Atmospheres Radiative Transfer for Exoplanet Atmospheres Acknowledgements We thanks Kevin Stevenson for useful feedback on the manuscript. V.P. acknowledges support from the Sagan Postdoctoral Fellowship through the NASA Exoplanet Science Institute. References Agúndez M, Parmentier V, Venot O, Hersant F Selsis F (2014) Pseudo 2d chemical model of hotjupiter atmospheres: application to hd b and hd b. A&A 564:A73, URL http: //dx.doi.org/ / / Arcangeli J, Desert JM, Parmentier V et al. (2017) Climate of a massive gas giant: Phase-curve of wasp-18b with hst. in prep

19 Exoplanet phase curves: observations and theory 19 Armstrong DJ, de Mooij E, Barstow J et al. (2016) Variability in the atmosphere of the hot giant planet hat-p-7 b. Nature Astronomy 1:0004 EP, URL s Batygin K, Stanley S Stevenson DJ (2013) Magnetically Controlled Circulation on Hot Extrasolar Planets. ApJ776:53 Cooper CS Showman AP (2006) Dynamics and Disequilibrium Carbon Chemistry in Hot Jupiter Atmospheres, with Application to HD b. ApJ649: Cowan NB Agol (2011a) The Statistics of Albedo and Heat Recirculation on Hot Exoplanets. ApJ729:54 Cowan NB Agol E (2008) Inverting Phase Functions to Map Exoplanets. ApJ678:L129 L132 Cowan NB Agol E (2011b) A Model for Thermal Phase Variations of Circular and Eccentric Exoplanets. ApJ726:82 Cowan NB, Machalek P, Croll B et al. (2012) Thermal Phase Variations of WASP-12b: Defying Predictions. ApJ747:82 Crossfield IJM (2015) Observations of Exoplanet Atmospheres. PASP127: Crossfield IJM, Hansen BMS, Harrington J et al. (2010) A New 24 µm Phase Curve for υ Andromedae b. ApJ723: Crossfield IJM, Knutson H, Fortney J et al. (2012) Spitzer/MIPS 24 µm Observations of HD b: Three Eclipses, Two and a Half Transits, and a Phase Curve Corrupted by Instrumental Sensitivity Variations. ApJ752:81 de Wit J, Lewis NK, Langton J et al. (2016) Direct Measure of Radiative and Dynamical Properties of an Exoplanet Atmosphere. ApJ820:L33 Demory BO, de Wit J, Lewis N et al. (2013) Inference of Inhomogeneous Clouds in an Exoplanet Atmosphere. ApJ776:L25 Demory BO, Gillon M, de Wit J et al. (2016) A map of the large day-night temperature gradient of a super-earth exoplanet. Nature532: Dyudina UA, Sackett PD, Bayliss DDR et al. (2005) Phase Light Curves for Extrasolar Jupiters and Saturns. ApJ618: Esteves LJ, De Mooij EJW Jayawardhana R (2015) Changing Phases of Alien Worlds: Probing Atmospheres of Kepler Planets with High-precision Photometry. ApJ804:150 Evans TM, Sing DK, Kataria T et al. (2017) An ultrahot gas-giant exoplanet with a stratosphere. Nature548:58 61 Feng YK, Line MR, Fortney JJ et al. (2016) The Impact of Non-uniform Thermal Structure on the Interpretation of Exoplanet Emission Spectra. ApJ829:52 Fortney JJ, Lodders K, Marley MS Freedman RS (2008) A Unified Theory for the Atmospheres of the Hot and Very Hot Jupiters: Two Classes of Irradiated Atmospheres. ApJ678: Fromang S, Leconte J Heng K (2016) Shear-driven instabilities and shocks in the atmospheres of hot Jupiters. A&A591:A144 Garcia Munoz A Isaak KG (2015) Probing exoplanet clouds with optical phase curves. ArXiv e-prints Guillot T, Burrows A, Hubbard WB, Lunine JI Saumon D (1996) Giant Planets at Small Orbital Distances. ApJ459:L35 Harrington J, Hansen BM, Luszcz SH et al. (2006) The Phase-Dependent Infrared Brightness of the Extrasolar Planet υ Andromedae b. Science 314: Haynes K, Mandell AM, Madhusudhan N, Deming D Knutson H (2015) Spectroscopic Evidence for a Temperature Inversion in the Dayside Atmosphere of Hot Jupiter WASP-33b. ApJ806:146 Heng K (2012) On the Existence of Shocks in Irradiated Exoplanetary Atmospheres. ApJ761:L1 Heng K Demory BO (2013) Understanding Trends Associated with Clouds in Irradiated Exoplanets. ApJ777:100 Hu R, Demory BO, Seager S, Lewis N Showman AP (2015) A Semi-analytical Model of Visiblewavelength Phase Curves of Exoplanets and Applications to Kepler- 7 b and Kepler- 10 b. ApJ802:51 Hubeny I, Burrows A Sudarsky D (2003) A Possible Bifurcation in Atmospheres of Strongly Irradiated Stars and Planets. ApJ594: