1 Introduction

Hot-Jupiter Inflation due to Deep Energy Deposition


Some extrasolar giant planets in close orbits—“hot Jupiters”—exhibit larger radii than that of a passively cooling planet. The extreme irradiation these hot Jupiters receive from their close in stars creates a thick isothermal layer in their envelopes, which slows down their convective cooling, allowing them to retain their inflated size for longer. This is yet insufficient to explain the observed sizes of the most inflated planets. Some models invoke an additional power source, deposited deep in the planet’s envelope. Here we present an analytical model for the cooling of such irradiated, and internally heated gas giants. We show that a power source , deposited at an optical depth , creates an exterior convective region, between optical depths and , beyond which a thicker isothermal layer exists, which in extreme cases may extend to the center of the planet. This convective layer, which occurs only for , further delays the cooling of the planet. Such a planet is equivalent to a planet irradiated with , where is an effective power-law index describing the radiative energy density as function of the optical depth for a convective planet . Our simple analytical model reproduces the main trends found in previous numerical works, and provides an intuitive understanding. We derive scaling laws for the cooling rate of the planet, its central temperature, and radius. These scaling laws can be used to estimate the effects of tidal or Ohmic dissipation, wind shocks, or any other mechanism involving energy deposition, on sizes of hot Jupiters.

Subject headings:
planetary systems — planets and satellites: general

1. Introduction

Hot Jupiters are a class of giant planets that orbit their parent stars much closer than Jupiter orbits the sun, therefore exhibiting high surface temperatures (see, e.g, Mayor & Queloz, 1995). The observed radii of many of these exoplanets are greater than Jupiter’s radius , with the largest planets reaching radii (Baraffe et al., 2010; Anderson et al., 2011; Chan et al., 2011; Hartman et al., 2011; Spiegel & Burrows, 2013). This discovery challenges theoretical evolution models, which predict that gas giants cool and contract to smaller radii, closer to , at their inferred age, which is often greater than 1Gyr (Burrows et al., 2000, 2007; Laughlin et al., 2005). To some extent, the larger radii are the result of the intense stellar irradiation, which induces a deep radiative, nearly isothermal, layer at the outer edge of the otherwise fully convective planets (Guillot et al., 1996; Arras & Bildsten, 2006). This radiative layer slows down the evolutionary cooling of the planet significantly, in comparison with isolated planets (Burrows et al., 2000; Chabrier et al., 2004). The slower cooling results in a higher bulk entropy at a given age, and therefore, a larger radius (see, e.g, Arras & Bildsten, 2006; Spiegel & Burrows, 2012; Marleau & Cumming, 2014).

Yet, the radii of at least some of the inflated hot Jupiters exceed theoretical evolutionary models, even with stellar irradiation taken into account (Baraffe et al., 2003; Burrows et al., 2007; Liu et al., 2008). To settle this discrepancy, a number of explanations have been suggested (see Baraffe et al., 2010; Fortney & Nettelmann, 2010; Spiegel & Burrows, 2013; Baraffe et al., 2014, for comprehensive reviews). The different explanations may be divided into several categories. While some works study the effects of extra power, deposited at some depth inside the planet, other focus on enhanced atmospheric opacities (Burrows et al., 2007), suppression of convective heat loss (double-diffusive convection, see Chabrier & Baraffe, 2007; Leconte & Chabrier, 2012), turbulent mixing (Youdin & Mitchell, 2010), or more consistent coupling of the day side and night side of the planet, and taking into account the three-dimensional nature of the problem (Hansen, 2008; Guillot, 2010; Spiegel & Burrows, 2010, 2013; Budaj et al., 2012; Dobbs-Dixon et al., 2012; Perez-Becker & Showman, 2013; Mayne et al., 2014; Rauscher & Showman, 2014).

In this work we focus on the effects of extra power sources deep in the planet. Possible heat sources include tidal dissipation due to orbital eccentricity (Bodenheimer et al., 2001, 2003; Gu et al., 2003; Winn & Holman, 2005; Jackson et al., 2008; Liu et al., 2008; Ibgui & Burrows, 2009; Miller et al., 2009; Ibgui et al., 2010, 2011; Leconte et al., 2010), “thermal tides” (Arras & Socrates, 2009a, b, 2010; Socrates, 2013), Ohmic heating (Batygin & Stevenson, 2010; Perna et al., 2010, 2012; Batygin et al., 2011; Huang & Cumming, 2012; Rauscher & Menou, 2013; Wu & Lithwick, 2013; Rogers & Showman, 2014), and conversion of a portion of the absorbed stellar flux into kinetic energy in the atmosphere, which is then dissipated at a greater depth (Guillot & Showman, 2002; Showman & Guillot, 2002). While some of these mechanisms (e.g. Ohmic heating, kinetic energy dissipation) are driven by the stellar irradiation, others (e.g. tidal dissipation) are not. If, as some studies suggest, excess planet inflation is correlated with the stellar irradiation (Laughlin et al., 2011; Schneider et al., 2011), then the former class of heat deposition mechanisms is favored.

We do not attempt to determine which power sources are more plausible. Alternatively, in the current work we study the effects of a general power source on the cooling, and therefore the radius, of a hot Jupiter. Specifically, we are interested in the dependence of the planet’s cooling history on the amount of additional power, and on the depth at which this power is deposited inside the planet. Recently, Spiegel & Burrows (2013) have conducted such a study numerically (see also Guillot & Showman, 2002; Baraffe et al., 2003). In this work we present a simple analytical model for the cooling of insolated gas giants in the presence of heat deposition somewhere in their envelope. Despite our simplifications, this model describes the essential processes, and provides an intuitive explanation to the main trends of the scenario.

The outline of the paper is as follows. In Section 2 we present our model for irradiated planets, without additional heat sources. The effects of additional deposited heat are studied in Sections 3 and 4. In Section 5 we compare our results to the calculations of Spiegel & Burrows (2013), and the main conclusions are summarized in Section 6.

2. Irradiated Planets

In this section we present our simplified “toy model” for irradiated gas giants, without additional heat sources. The model is one dimensional, and does not differentiate between the day and night sides of the planet (see Spiegel & Burrows, 2013). Many aspects of the deposition-free model are essentially similar to a previous analytical work by Arras & Bildsten (2006), and are also implemented by Youdin & Mitchell (2010).

Isolated gas giants, such as Jupiter, are almost fully convective, with convection beginning at an optical depth (Trafton, 1964, 1967; Hubbard, 1968; Stevenson & Salpeter, 1977); see also Guillot et al. (1994) for a reexamination. Insolated planets, on the other hand, are expected to develop a deep radiative envelope (Guillot et al., 1996).

2.1. Convective Interior Structure

The analysis of irradiated planets is simpler when considering their profiles in the plane, with


denoting the optical depth at radius , from the planet surface, at radius , and is the radiative energy density. The density, temperature, and opacity are denoted by , , and , respectively, and is the radiation constant. Our main simplification is to assume that this profile can be approximated by a power law in the convective interior


with denoting the central radiation energy density, determined by the central temperature .

This scaling can be derived in the case of a power-law opacity , as assumed by Arras & Bildsten (2006) and Youdin & Mitchell (2010). The opacity is combined with the adiabatic temperature and pressure relation in the convective region (see, e.g., Hubbard, 1977; Saumon et al., 1996) , and the resulting polytropic density profile at the edge of the planet ( denotes the central density; see, e.g., Matzner & McKee, 1999). We neglect super-adiabatic corrections to the convective temperature profile (Trafton, 1967; Chabrier & Baraffe, 2007). These power laws yield


and , with denoting the estimate for the central opacity, if the power-law opacity could be continued to the center.

The Rosseland opacity (we assume a gray approximation) in the range of interest is depicted in Figure 1 of Arras & Bildsten (2006). The molecular dominated opacity below is fit by . Above , the opacity rises steeply, with (), due to , whose population increases primarily as a consequence of metal ionization (see, e.g., Kippenhahn & Weigert, 1994). At even higher temperatures, above , the vanishes, due to lack of neutral hydrogen, and the opacity is given by the Kramers bound-free opacity law, with , (see, e.g., Kippenhahn & Weigert, 1994; Padmanabhan, 2000).

As explained in Sections 2.3 and 2.4, the internal luminosity of a planet is determined by the location of the radiative-convective boundary. The optical depth (and pressure level) of this boundary is given by the opacity power-law which is suitable for the (almost isothermal) radiative-zone temperature, and does not depend on the deeper opacity structure. For planets without additional heating, this temperature is typically below , and we therefore get , by taking (Peebles, 1964; Hubbard, 1977; Saumon et al., 1996). Throughout the majority of this work we adopt the approximation of small inflations relative to the zero-temperature radius (see Section 4), implying roughly constant central densities and pressures. Consequently, a constant is suitable for a Jupiter-size planet. Inspection of numerical evolution models (see Table 1 in Burrows et al., 1997) reveals that the ratio between central and photosphere temperatures, which is given by setting in Equation (2), , is , and changes only by a factor of 2 during the cooling of a Jupiter mass isolated planet, consistent with our estimate.

In Section 3 we show that heat deposition raises the radiative-convective boundary temperature above , requiring a more general form, which takes into account the strong dependence of the opacity on the temperature


In this case we get from Equation (3), which is suitable for temperatures above . As we show in Section 2.3, the convective profile must satisfy .

2.2. Radiative Envelope Structure

Proximity to a star dictates an equilibrium temperature on the planet surface, which we take as the photosphere, at . Our model is spherical so we assume an even distribution of the absorbed heat over the planet surface, for which (see, e.g., Guillot et al., 1996)


with and denoting the stellar temperature and radius, respectively, and , are the planet’s orbital distance and albedo, respectively. This equilibrium temperature defines an energy density , and luminosity , where is the Stefan-Boltzmann constant.

The radiation energy is related to the internal luminosity of the planet through the diffusion approximation


where is the speed of light. Close to the surface, is constant (in space), since it is the result of cooling of the entire planet, and the heat capacity of the outer layers is negligible. Therefore, we can integrate Equation (6) and find the radiative profile


where the distinction between the boundary conditions at and is negligible, since we are generally interested in . Bodenheimer et al. (2003) use a similar boundary condition (with at ), while Arras & Bildsten (2006) and Youdin & Mitchell (2010) define a temperature for their boundary condition, which is correlated with , therefore leading to similar conclusions. More elaborate boundary conditions, which incorporate a frequency-dependent opacity, and require the solution of the radiative transfer equation, or the use of a two-stream approximation (Barman et al., 2001; Hubeny et al., 2003; Hansen, 2008; Guillot, 2010; Guillot & Havel, 2011), are beyond the scope of this work, and their effect on the temperature deep inside the radiative zone (and therefore on the cooling rate) can be modeled by a somewhat different effective . We remark that an atmosphere which is more transparent in the optical regime (although this does not seem to be the case for hot Jupiters; see, e.g., Spiegel & Burrows, 2013) may induce a greenhouse effect, which can be treated with our model as heat deposited at , where the depth for optical light is unity.

Planar symmetry, and the focus on the outer layers of the planet throughout this work, are justified both self-consistently in Section 2.4, and by Guillot et al. (1996), who find that during the 10Gyr evolution of a hot Jupiter, the radiative region penetrates only a few percent of the radius.

2.3. Radiative-Convective Transition

According to the Schwarzschild criterion, when the radiative temperature gradient becomes larger than the adiabatic gradient, convective instability develops. Differentiating Equations (2) and (7) shows that convection sets in when


We see from Equation (8) that at the profile is radiative. We also see that to ensure convection at , we must have .

Figure 1.— Transition between convective (solid black line) and radiative (dashed blue line) profiles. The convective profile follows Equation (2), while the radiative profile is according to Equation (7). The combined profile follows the radiative profile for , and the convective profile for . Two additional extreme cases are plotted for the same values of and (dotted blue lines). A very high value corresponds to a fully radiative planet, while a very low value of leads to , and therefore a fully convective planet. Typical values of the optical depth and temperature are given.

The transition between convective and radiative profiles is illustrated schematically in Figure 1. For a planet with given interior (, ) and exterior () boundary conditions, the convective profile is given by Equation (2), while the radiative profile is given by Equation (7). Since the radiative profile is linear in , the luminosity and the radiative-convective transition point are calculated by finding the tangent to the adiabatic profile. For , the radiative gradient is smaller than the adiabatic, so the profile is radiative. For , the radiative gradient is larger than the adiabatic, so the profile is convective. Quantitatively, the tangency point is found by equating Equations (2) and (7) and their derivatives. The solution for is


which corresponds to a radiation energy density


implying that the radiative zone is isothermal to within a factor of . Equation (9) indicates that the radiative region thickens with increasing solar irradiation and with decreasing central temperature (see also Guillot et al., 1996; Burrows et al., 2000; Arras & Bildsten, 2006; Youdin & Mitchell, 2010).

2.4. Internal Luminosity and Cooling

Equations (8) and (9) connect the internal luminosity to the central temperature


Note that Equation (11) breaks down when . In this case, the planet is considered to be fully convective, and its luminosity is determined by the radiation energy density at the photosphere . Substituting in Equation (2) we find


and luminosity


Comparing Equations (9) and (12) shows that the condition is equivalent to the condition . In other words, as long as the central temperature of the planet is high enough to ensure that a convective profile will reach without dropping below , the planet remains fully convective. However, as the planet cools, the central temperature decreases, and the convective profile would have predicted temperatures below . This situation is unphysical, and an approximately isothermal radiative zone develops.

To summarize, using Equations (11) and (13), the internal luminosity of the planet, as a function of its central radiation energy density is given by


where the transition is when . Since , Equation (14) indicates that the proximity to a star, which is manifested for , diminishes the internal luminosity and slows down the planet’s cooling (see also Guillot et al., 1996; Burrows et al., 2000; Arras & Bildsten, 2006; Youdin & Mitchell, 2010; Spiegel & Burrows, 2013).

As a planet cools, it makes the transition from an isolated planet, which is not affected by stellar irradiation, and matches the high limit of Equation (14), to an insolated planet, which matches the low limit. The age, at which this transition takes place, may be estimated by writing an evolution equation for the central temperature of an isolated planet


where is Boltzmann’s constant, is the proton mass, is the mass of the planet, and denotes time. Equation (15) is valid for the cooling degenerate phase of the planet’s evolution, which follows a rapid non-degenerate contraction phase (see Guillot, 2005; Arras & Bildsten, 2006). Substituting the luminosity from Equation (14), we get


Assuming that the planet cools down from high temperatures, Equation (16) predicts the cooling time to a temperature , which marks the transition to the insolated regime (at this stage the photosphere reaches a temperature ):


were the numerical value is estimated for (Jupiter mass) and . This estimate roughly fits Jupiter, which is a few Gyr old (Guillot, 2005), has an equilibrium temperature , and is close to the transition, since (Hanel et al., 1981). Strongly irradiated planets orbit their parent star closer than Jupiter by a factor of 100, thus having, using Equation (5), , and therefore reach the insolated regime and develop a radiative zone after less than years (consistent with Guillot et al., 1996).

In the following, insolated stage of planetary evolution, Equation (15) is combined with the low limit of Equation (14) to produce a cooling equation


with the numerical values estimated for , . This result is consistent with Guillot et al. (1996). At an age of 8Gyr, the radiative region encompasses a few percent in radius, , consistent with Guillot et al. (1996), and self-consistently justifying our focus on the outer layers of the planet.

A third stage in a cooling planet’s life, equilibrium, is reached when . In this case, according to Equation (9), the planet becomes fully radiative, and the luminosity vanishes (since the temperature gradient vanishes). This condition determines the planet’s equilibrium central temperature, and therefore its equilibrium radius. More precisely, the planet becomes fully radiative when drops below . It is easy to see, from Equation (18), that without additional heat deposition, planets never reach equilibrium (the timescales are too long).

The three stages of planetary evolution, described in this section, are displayed in Figure 1, for planets with the same central temperature and optical depth, but different equilibrium temperatures.

3. Deposition of Additional Power

In this section we discuss how the internal luminosity is affected by an additional heat source, which deposits power at some depth inside the planet. Our goal is to find , with denoting the deposited luminosity, and denoting the optical depth where the power is deposited. In Section 2.4, and specifically in Equation (14), we found , which we abbreviate hereafter as . We focus on planets which have developed a substantial radiative zone, , so is determined by the low limit of Equation (14). The motivation is that strongly irradiated planets develop such a radiative zone in less than years, as explained in Section 2.4 (see also Guillot et al., 1996).

Without heat sources, the radiative zone is described by Equation (7). The addition of a deposited power at introduces a jump in the profile slope. Explicitly, the radiative profile is altered:


with denoting the total luminosity for . The radiative profile, with deposited heat included, is displayed in Figure 2. As seen in Figure 2, the deposited heat effectively changes to some higher . The internal luminosity is found, as in Section 2.3, by calculating the tangent from the point to the convective profile, which is equivalent to calculating the tangent from the point . Therefore, the results of Section 2.3 are reproduced, but with instead of . Since, from Equation (11), , the internal luminosity decreases. This result is also evident graphically in Figure 2: the slope of the radiative tangent decreases, as well as the deeper penetration of the radiative zone into the planet interior.

3.1. Secondary Convective Region

In Figure 2 we assumed that the outer layers remain radiative. However, if the total power, , which has to be evicted from the planet, is high enough, then it cannot be delivered radiatively, and a secondary convective instability develops. According to Equations (8) and (19), with replacing as the evicted power, convection appears at an optical depth


where the radiation energy density reaches . If we focus on intense deposition , then we may approximate , and using Equation (20), .

Figure 2.— Modification of the radiative profile due to the deposition of heat at optical depth . The convective profile (solid black line), and the radiative profile without a heat source (dashed blue line) are similar to Figure 1. The modified radiative profile (dot-dashed red line) is described by Equation (19). The imaginary continuation of the profile with a slope down to (dotted red line) defines . As in Figure 1, the combined profile follows the radiative profile up to the tangency point, and then it follows the convective profile.
Figure 3.— Emergence of a secondary convective region from to (solid green line). The convective (solid black line) and radiative (dashed blue line) profiles without a heat source are similar to Figures 1 and 2. The radiative profile with heat deposited at (upper dot-dashed red line) is similar to Figure 2, and is inadequate in this case. Due to the secondary convective region, the radiative profile is altered, and it follows the lower dot-dashed red line, up to the tangency point with the main convective profile.

The temperature profile, with deposited heat added, can be divided into two cases, depending on and . If , or equivalently , The depiction in Figure 2 is accurate. As seen in Figure 2, in this case is restricted by a maximal value , limiting the reduction in internal luminosity, using Equation (11), to an insignificant effect .

In the case where , a secondary convective region emerges between and . This region is displayed in Figure 3. The modified luminosity can be calculated by finding . Using Equation (20), approximating , and noting that in this regime , we find


where is adiabatically related to . Combining this result with Equation (11), we get , which we rewrite as


in order to interpolate with the regime.

Equation (22) shows that the figure of merit of the heat source, in terms of its influence on the internal luminosity, is , which must satisfy for a significant effect. For example, if the heat is deposited at , then is needed. This result is intuitive, since the equilibrium luminosity itself is deposited at . If, on the other hand, the heat is deposited at the top of the convective region, , then, using Equation (11), is sufficient. Although not evident from Equation (22), is required for a significant effect (see also Guillot & Showman, 2002), even for deposition in the convective region . Otherwise, the deposited power simply replaces part of the internal luminosity as the convectively delivered heat which maintains the radiative profile at the envelope, and the temperature profile remains unchanged (with respect to ). In this case, the luminosity decreases by .

The secondary convective region, described in this section, was also noted in numerical works (Guillot & Showman, 2002; Batygin et al., 2011; Wu & Lithwick, 2013), for high dissipated power.

3.2. Cooling and Equilibrium

Substantial additional heating raises the effective equilibrium temperature, and therefore the inner radiative-convective boundary temperature above , as seen in Equation (21). Under the influence of such a heat source, using Equation (22), the cooling Equation (18) in the insolated regime changes to


slowing down the planet’s cooling, with the power calculated using the opacity values and , following the discussion in Section 2.

Equilibrium is reached when (see Section 2). According to Equations (9) and (21), the required heating intensity for this scenario is


with denoting, as before, the original radiative-convective transition, without additional heating. We may rewrite Equation (24), using Equation (11), as , showing that, as intuitively expected, a deposition of order forces equilibrium, when the deposition is at the center. For a given heat source, Equation (24) defines the equilibrium central temperature , for which evolutionary cooling not only slows down, but stops entirely


where we have again interpolated with the weak heating regime.

If, due to very intensive (or deep) heating, , the opacity behaves according to the Kramers law, with and , prohibiting convection (see Section 2). In this case, the planet becomes fully radiative from inward, even before equilibrium is reached. The internal luminosity in this case is given by Equation (6), . Combining this result with Equation (15), we find that the cooling time is , with hotter (more inflated) planets reaching equilibrium first (planets with radii close to cool even faster, due to the dependence of the opacity on the density). This regime corresponds to the equilibrium-inflation paradigm of hot Jupiters, which typically involves very deep heating (Guillot & Showman, 2002; Chabrier et al., 2004; Burrows et al., 2007; Liu et al., 2008; Spiegel & Burrows, 2013).

4. Effect of Heating on Planet Radius

In Sections 2 and 3 we found the effects of stellar irradiation and additional heat deposition on the planet luminosity. The diminished luminosity results in slower cooling, and therefore a higher bulk temperature and larger radius at a specific age, in addition to a larger equilibrium radius. For a rough estimate of this inflation, we introduce a simple model which relates the plant’s radius to its central temperature .

We start by writing a simple equation of state (EOS)


where is the electron degeneracy term, and is the electrostatic correction ( denotes the Planck constant, and denote the electron mass and charge, respectively, and is the Bohr radius; see, e.g., Padmanabhan, 2001, where we assume a mixture of hydrogen and helium, so the atomic weight and charge are approximately 1). The rightmost term is the ideal gas pressure of the non-degenerate ions. We are interested in the regime where the degeneracy parameter , which is adequate for HD 209458b (see Guillot, 2005). In this regime, the electrons are almost completely degenerate, and the ion thermal term is a small correction. Since the degeneracy pressure is inversely related to the particle mass, the ions are non-degenerate. The thermal contribution of the degenerate electrons is quadratic in (see, e.g., Chandrasekhar, 1939), and is therefore negligible. This EOS is also used by Arras & Bildsten (2006).

We proceed by replacing the profiles , with the characteristic values . This simplification is justified by noting that the entire planet (except for the radiative isothermal envelope) can be described approximately with a single almost degenerate polytrope. The distinction between a degenerate core and a non-degenerate envelope is insignificant, since the temperature profile is almost parallel to lines of constant (Guillot, 2005). A similar argument renders the distinction between a core without electrostatic corrections, and an envelope with significant corrections, also insignificant in comparison with the bulk electrostatic correction. For a given central temperature , the contribution of the isothermal () radiative envelope to the radius is negligible, because it scales with the temperature, and . We emphasize that the quasi-degenerate treatment in this section only holds for mild inflations, and should be taken only as an order of magnitude estimate for extremely inflated planets, with radii . In these cases, numerical radius-temperature relations are more adequate, as detailed below.

In equilibrium, the pressure of Equation (26) balances the gravitational pressure, of order ( is the gravitational constant),


For small variations of the temperature, differentiating Equation (27), while holding the planet mass constant, yields a linear relation


In the degenerate limit (), the radius increases with temperature, while in the non-degenerate limit (), the planet shrinks as the temperature increases. As mentioned above, the non-degenerate limit corresponds to the initial contraction phase, accompanied by an increase in temperature, while the degenerate limit describes the later cooling of the planet, which is accompanied by a decrease in radius. The inflated hot Jupiters are on the late-time cooling branch of this scenario (the transition is at a radius around , see Guillot et al., 1996).

We make a linear approximation for low temperatures, using Equations (27) and (28)


with denoting the inflation relative to the zero-temperature radius , is the surface gravity, and noting that Jupiter’s mass is close the well-known inversion of the zero-temperature mass-radius relation (Padmanabhan, 2001). Since we focus on hot Jupiters with , the electrostatic correction term does not influence our order of magnitude estimate. This relation fits well numerical radius-temperature curves (for example, Burrows et al., 1997, Table 1), up to a radius of , where the linear approximation breaks down (since it assumes ). See Arras & Bildsten (2006) for a similar linear relation between the radius change and the central temperature.

Although Equation (29) was derived for a constant planet mass, it is also valid (up to numerical factors) when the planet surface gravity , instead of the planet mass is kept constant. In this case, the gravitational pressure in Equation (27) may be written as , and the radius increase is found by differentiating . This scenario is interesting for comparison with the results of Spiegel & Burrows (2013), who studied models with constant surface gravity.

5. Comparison to Numerical Calculations

Guillot & Showman (2002) and Baraffe et al. (2003) find that less heat dissipation is needed to inflate hot Jupiters if the heat is deposited in deeper mass fractions of the planet (see also Wu & Lithwick, 2013). In this section we compare our estimates to the systematic work of Spiegel & Burrows (2013), henceforth SB2013, who studied the effects of different amounts of additional heat, deposited at different pressure levels.

In their Section 3.3, specifically Figure 4, SB2013 deposit different fractions of the incident flux at different pressures inside a planet similar to HD 209458b. For each value of additional power and deposition pressure, SB2013 present the inflated planet radius. We emphasize that, in their Section 3.3, SB2013 concentrate on planets which have not yet reached equilibrium, but are still contracting. Moreover, they compare models with the same effective temperature, defined as , and the same surface gravity, resulting in models with different ages and masses.

By combining Equations (14) and (22), we relate the central temperature to the deposited power,


The inflated planet radius can be related to the deposited power by combining Equation (30) with Equation (29)


with denoting the expansion without extra heating. This numeric value is estimated by substituting the parameters and , which are kept constant in SB2013, the observed temperature (Crossfield et al., 2012), and taking the adequate and . The model of SB2013 for HD209458b is indeed inflated by , due to the irradiation of the parent star, even without extra power sources (see also Baraffe et al., 2003; Burrows et al., 2007), in agreement with our estimate. The zero-temperature radius is estimated using Guillot (2005) and Baraffe et al. (2010), for the given surface gravity. An alternative estimate of the inflation may be given by relating the central temperature, found in Equation (30), to , using the numerically calculated curve given in Burrows et al. (1997). This analysis is more accurate at large expansions.

In Figure 4 we compare the results of SB2013 (their Figure 4) to the equation


with denoting the expansion relative to the radius without extra heat deposition ( in our case), , an order unity fitting parameter, and from Equation (31). The fraction of incident flux is given by SB2013, and the translation of deposition pressure to optical depth is done approximately by , with taken from Arras & Bildsten (2006), and uncertainties going into the fitting parameter . The parameters in Figure 4, , are not, strictly speaking, a best fit, but only an illustration. A second fit, with and the predicted , is also provided. This fit, however, incorporates the more accurate nonlinear radius correction (Burrows et al., 1997) to Equation (32).

While our model is not exact, as seen in Figure 4, it does manage to describe the main trends of atmospheric heating, and, given the general uncertainties and approximations of the model, it fits the numerical calculations reasonably well over several orders of magnitude in pressure level and flux. When the realistic radius-temperature relation is used, the analytically predicted power fits the results adequately. With the approximated linear radius-temperature relation, however, we fit to the numerical results.

Figure 4.— Comparison of our prediction for planet inflation (relative to ), due to heating (expressed as a fraction of the incident flux) at different pressure levels, with Spiegel & Burrows (2013), for constant surface gravity , effective temperature , and incident irradiation of HD 209458b. Our analytical model (dashed blue lines) is given by Equation (32), with fitted parameters . A second model (dot-dashed red lines), which takes into account the nonlinear radius-temperature relation, is also provided, with a fitted , and the analytically derived . The results of SB2013 (solid black lines, marked with circles) are taken from their Figure 4.

In the “weak heating” regime , it is possible to linearize Equation (32)


which has the advantage of a single fitting parameter , allowing for a more intuitive interpretation of the results.

6. Conclusions

In this work we presented a simplified analytical model for irradiated giant gas planets, which includes an additional heat source, deposited at an arbitrary depth inside the planet’s atmosphere.

For an irradiated planet with no extra heat sources, we found useful approximate scaling relations for the internal luminosity of the planet, and for the radiative-convective boundary, marked by an optical depth . Explicitly, the radiative-convective boundary and the internal luminosity are given by


with , denoting the equilibrium luminosity and radiation energy density imposed by the stellar irradiation, respectively ( is the equilibrium temperature). is the central radiation energy density of the planet, and represents its optical depth (with opacities extrapolated from the outer layers of the planet inward). The power is related to the adiabatic index of the planet and to the opacity near the planet’s edge. In addition, we identified the condition , as the boundary separating insolated (irradiated) planets () from isolated ones (). These relations serve as an intuitive interpretation of well-known results (Guillot et al., 1996; Burrows et al., 2000; Arras & Bildsten, 2006; Fortney & Nettelmann, 2010; Youdin & Mitchell, 2010; Spiegel & Burrows, 2013).

We addressed the effects of possible additional heat sources, parametrized by additional luminosity , deposited at an optical depth . Such heat sources lower the internal luminosity of the planet, slow its evolutionary cooling, and are therefore candidates for solving the puzzle of over-inflated hot Jupiters (see Baraffe et al., 2010; Fortney & Nettelmann, 2010; Spiegel & Burrows, 2013; Baraffe et al., 2014, and references within). According to our model, the figure of merit in this context is , with necessary for a significant effect. Concretely, an outer convective layer is formed in this case, between optical depths and , which is equivalent to an enhanced stellar irradiation


The internal luminosity is therefore reduced from a value of without extra heating to, roughly,


with a further threshold of required for a significant effect.

Since the inflation of a planet, relative to its zero-temperature radius, is proportional to its central temperature, these relations provide an intuitive explanation to the dependence of inflation on the deposition depth (Guillot & Showman, 2002; Baraffe et al., 2003; Wu & Lithwick, 2013), and to the large equilibrium radii, in the case of very deep heating (Guillot & Showman, 2002; Burrows et al., 2007; Liu et al., 2008; Spiegel & Burrows, 2013).

Quantitatively, we compared our model to the recent parameter survey of Spiegel & Burrows (2013). The model, though not exact, fits the numerical results adequately, considering its simplifying assumptions and approximations.

The scaling laws presented in this work, combined with potential observational correlations (see, e.g., Laughlin et al., 2011; Schneider et al., 2011), may be used to rule out or to favor specific energy sources, or even the entire additional heat source paradigm.

This research was partially supported by ISF, ISA and iCore grants, and a Packard Fellowship. RS would like to thank the Max Planck Institute for Extraterrestrial Physics and the Humbolt Foundation for support and worm hospitality while this research was completed. We thank Peter Goldreich for insightful discussions, and the anonymous referee for helpful comments.


  1. Anderson, D. R., Smith, A. M. S., Lanotte, A. A., et al. 2011, MNRAS, 416, 2108
  2. Arras, P., & Bildsten, L. 2006, ApJ, 650, 394
  3. Arras, P., & Socrates, A. 2009a, arXiv:0901.0735
  4. Arras, P., & Socrates, A. 2009b, arXiv:0912.2318
  5. Arras, P., & Socrates, A. 2010, ApJ, 714, 1
  6. Baraffe, I., Chabrier, G., & Barman, T. 2010, RPPh, 73, 016901
  7. Baraffe, I., Chabrier, G., Barman, T. S., Allard, F. & Hauschildt, P. H. 2003, A&A, 402, 701
  8. Baraffe, I., Chabrier, G., Fortney, J., & Sotin, C. 2014, Protostars and Planets VI (Tucson: Univ. Arizona Press), 763
  9. Barman, T. S., Hauschildt, P. H., & Allard, F. 2001, ApJ, 556, 885
  10. Batygin, K., & Stevenson, D. J. 2010, ApJ, 714, 238
  11. Batygin, K., Stevenson, D. J., & Bodenheimer, P. H. 2011, ApJ, 738, 1
  12. Bodenheimer, P., Laughlin, G., & Lin, D. N. C. 2003, ApJ, 592, 555
  13. Bodenheimer, P., Lin, D. N. C., & Mardling, R. A. 2001, ApJ, 548, 466
  14. Budaj, J., Hubeny, I., & Burrows, A. 2012, A&A, 537, A115
  15. Burrows, A., Guillot, T., Hubbard, W. B., et al. 2000, ApJ, 534, L97
  16. Burrows, A., Hubeny, I., Budaj, J., & Hubbard, W. B. 2007, ApJ, 661, 502
  17. Burrows, A., Marley, M., Hubbard, W. B., et al. 1997, ApJ, 491, 756
  18. Chabrier, G., & Baraffe, I. 2007, ApJ, 661, L81
  19. Chabrier, G., Barman, T., Baraffe, I., Allard, F., & Hauschildt, P. H. 2004, ApJ, 603, L53
  20. Chan, T., Ingemyr, M., Winn, J. N., et al. 2011, ApJ, 141, 179
  21. Chandrasekhar, S.  1939, An Introduction to the Study of Stellar Structure (Chicago, IL: Univ. Chicago Press)
  22. Crossfield, I. J. M., Knutson, H., Fortney, J., et al. 2012, ApJ, 752, 81
  23. Dobbs-Dixon, I., Agol, E., & Burrows A. 2012, ApJ, 751, 87
  24. Fortney, J. J., & Nettelmann, N. 2010, Space Sci. Rev., 152, 423
  25. Gu, P., Lin, D. N. C., & Bodenheimer, P. H. 2003, ApJ, 588, 509
  26. Guillot, T. 2005, AREPS, 33, 493
  27. Guillot, T. 2010, A&A, 520, A27
  28. Guillot, T., Burrows, A., Hubbard, W. B., Lunine, J. I., & Saumon, D. 1996, ApJ, 459, L35
  29. Guillot, T., Gautier, D., Chabrier, G., & Mosser, B. 1994, Icarus, 112, 337
  30. Guillot, T., & Havel, M. 2011, A&A, 527, A20
  31. Guillot, T., & Showman, A. P. 2002, A&A, 385, 156
  32. Hanel, R. A., Conrath, B. J., Herath, L. W., Kunde, V. G., & Pirraglia, J. A. 1981, J. Geophys. Res., 86, 8705
  33. Hansen, B. M. S. 2008, ApJS, 179, 484
  34. Hartman, J. D., Bakos, G. À., Torres, G., et al. 2011, ApJ, 742, 59
  35. Huang, X., & Cumming, A. 2012, ApJ, 757, 46
  36. Hubbard, W. B. 1968, ApJ, 152, 745
  37. Hubbard, W. B. 1977, Icarus, 30, 305
  38. Hubeny, I., Burrows, A., & Sudarsky, D. 2003, ApJ, 594, 1011
  39. Ibgui, L., & Burrows, A. 2009, ApJ, 700, 1921
  40. Ibgui, L., Burrows, A., & Spiegel, D. M. 2010, ApJ, 713, 751
  41. Ibgui, L., Spiegel, D. M., & Burrows, A. 2011, ApJ, 727, 75
  42. Jackson, B., Greenberg, R., & Barnes, R. 2008, ApJ, 681, 1631
  43. Kippenhahn, A., & Weigert, A. 1994, Stellar Structure and Evolution (Berlin: Springer)
  44. Laughlin, G., Crismani, M., & Adams, F. C. 2011, ApJ, 729, L7
  45. Laughlin, G., Wolf, A., Vanmunster, T., et al. 2005, ApJ, 621, 1072
  46. Leconte, J., & Chabrier, G. 2012, A&A, 540, A20
  47. Leconte, J., Chabrier, G., Baraffe, I., & Levard, B. 2010, A&A, 516, A64
  48. Liu, X., Burrows, A., & Ibgui, L. 2008, ApJ, 687, 1191
  49. Marleau, G. D., & Cumming, A. 2014, MNRAS, 437, 1378
  50. Matzner, C. D., & McKee, C. F. 1999, ApJ, 510, 379
  51. Mayne, N. J., Baraffe, I., Acreman, D. M., et al. 2014, A&A, 561, A1
  52. Mayor, M., & Queloz, D. 1995, Nature, 378, 355
  53. Miller, N., Fortney, J. J., & Jackson, B. 2009, ApJ, 702, 1413
  54. Padmanabhan, T. 2000, Theoretical Astrophysics, Volume 1: Astrophysical Processes (Cambridge: Cambridge Univ. Press)
  55. Padmanabhan, T. 2001, Theoretical Astrophysics, Volume 2: Stars and Stellar Systems (Cambridge: Cambridge Univ. Press)
  56. Peebles, P. J. E. 1964, ApJ, 140, 328
  57. Perez-Becker, D., & Showman, A. P. 2013, ApJ, 776, 134
  58. Perna, R., Heng, K., & Pont, F. 2012, ApJ, 751, 59
  59. Perna, R., Menou, K., & Rauscher, E. 2010, ApJ, 719, 1421
  60. Rauscher, E., & Menou K. 2013, ApJ, 764, 103
  61. Rauscher, E., & Showman A. P. 2014, ApJ, 784, 160
  62. Rogers, T. M., & Showman A. P. 2014, ApJ, 782, L4
  63. Saumon, D., Hubbard, W. B., Burrows A., et al. 1996, ApJ, 460, 993
  64. Schneider, J., Dedieu, C., Le Sidaner, P., Savalle, R., & Zolotukhin, I. 2011, A&A, 532, A79
  65. Showman, A. P., & Guillot, T. 2002, A&A, 385, 166
  66. Socrates, A. 2013, arXiv:1304.4121
  67. Spiegel, D. S., & Burrows, A. 2010, ApJ, 722, 871
  68. Spiegel, D. S., & Burrows, A. 2012, ApJ, 745, 174
  69. Spiegel, D. S., & Burrows, A. 2013, ApJ, 772, 76
  70. Stevenson, D. J., & Salpeter E. E. 1977, ApJS, 35, 221
  71. Trafton, L. M. 1964, ApJ, 140, 1340
  72. Trafton, L. M. 1967, ApJ, 147, 765
  73. Winn, J. N., & Holman, M. J. 2005, ApJ, 159, L159
  74. Wu, Y., & Lithwick, Y. 2013, ApJ, 763, 13
  75. Youdin, A. N., & Mitchell, J. L. 2010, ApJ, 721, 1113
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description