We present new equations of state (EOS) for hydrogen and helium covering a wide range of temperatures from 60 K to 10$^7$ K and densities from $10^{-10}$ g/cm$^3$ to $10^3$ g/cm$^3$. They include an extended set of ab initio EOS data for the strongly
correlated quantum regime with an accurate connection to data derived from other approaches for the neighboring regions. We compare linear-mixing isotherms based on our EOS tables with available real-mixture data. A first important astrophysical application of this new EOS data is the calculation of interior models for Jupiter and the comparison with recent results. Secondly, mass-radius relations are calculated for Brown Dwarfs which we compare with predictions derived from the widely used EOS of Saumon, Chabrier and van Horn. Furthermore, we calculate interior models for typical Brown Dwarfs with different masses, namely Corot-3b, Gliese-229b and Corot-15b, and the Giant Planet KOI-889b. The predictions for the central pressures and densities differ by up to 10$%$ dependent on the EOS used. Our EOS tables are made available in the supplemental material of this paper.
The hot-Jupiter HAT-P-2b has become a prime target for Spitzer Space Telescope observations aimed at understanding the atmospheric response of exoplanets on highly eccentric orbits. Here we present a suite of three-dimensional atmospheric circulation
models for HAT-P-2b that investigate the effects of assumed atmospheric composition and rotation rate on global scale winds and thermal patterns. We compare and contrast atmospheric models for HAT-P-2b, which assume one and five times solar metallicity, both with and without TiO/VO as atmospheric constituents. Additionally we compare models that assume a rotation period of half, one, and two times the nominal pseudo-synchronous rotation period. We find that changes in assumed atmospheric metallicity and rotation rate do not significantly affect model predictions of the planetary flux as a function of orbital phase. However, models in which TiO/VO are present in the atmosphere develop a transient temperature inversion between the transit and secondary eclipse events that results in significant variations in the timing and magnitude of the peak of the planetary flux compared with models in which TiO/VO are omitted from the opacity tables. We find that no one single atmospheric model can reproduce the recently observed full orbit phase curves at 3.6, 4.5 and 8.0 microns, which is likely due to a chemical process not captured by our current atmospheric models for HAT-P-2b. Further modeling and observational efforts focused on understanding the chemistry of HAT-P-2bs atmosphere are needed and could provide key insights into the interplay between radiative, dynamical, and chemical processes in a wide range of exoplanet atmospheres.
We present new calculations of Rosseland and Planck gaseous mean opacities relevant to the atmospheres of giant planets and ultracool dwarfs. Such calculations are used in modeling the atmospheres, interiors, formation, and evolution of these objects
. Our calculations are an expansion of those presented in Freedman et al. (2008) to include lower pressures, finer temperature resolution, and also the higher metallicities most relevant for giant planet atmospheres. Calculations span 1 microbar to 300 bar, and 75 K to 4000 K, in a nearly square grid. Opacities at metallicities from solar to 50 times solar abundances are calculated. We also provide an analytic fit to the Rosseland mean opacities over the grid in pressure, temperature, and metallicity. In addition to computing mean opacities at these local temperatures, we also calculate them with weighting functions up to 7000 K, to simulate the mean opacities for incident stellar intensities, rather than locally thermally emitted intensities. The chemical equilibrium calculations account for the settling of condensates in a gravitational field and are applicable to cloud-free giant planet and ultracool dwarf atmospheres, but not circumstellar disks. We provide our extensive opacity tables for public use.
The formation of clouds affects brown dwarf and planetary atmospheres of nearly all effective temperatures. Iron and silicate condense in L dwarf atmospheres and dissipate at the L/T transition. Minor species such as sulfides and salts condense in mi
d-late T dwarfs. For brown dwarfs below Teff=450 K, water condenses in the upper atmosphere to form ice clouds. Currently over a dozen objects in this temperature range have been discovered, and few previous theoretical studies have addressed the effect of water clouds on brown dwarf or exoplanetary spectra. Here we present a new grid of models that include the effect of water cloud opacity. We find that they become optically thick in objects below Teff=350-375 K. Unlike refractory cloud materials, water ice particles are significantly non-gray absorbers; they predominantly scatter at optical wavelengths through J band and absorb in the infrared with prominent features, the strongest of which is at 2.8 microns. H2O, NH3, CH4, and H2 CIA are dominant opacity sources; less abundant species such as may also be detectable, including the alkalis, H2S, and PH3. PH3, which has been detected in Jupiter, is expected to have a strong signature in the mid-infrared at 4.3 microns in Y dwarfs around Teff=450 K; if disequilibrium chemistry increases the abundance of PH3, it may be detectable over a wider effective temperature range than models predict. We show results incorporating disequilibrium nitrogen and carbon chemistry and predict signatures of low gravity in planetary- mass objects. Lastly, we make predictions for the observability of Y dwarfs and planets with existing and future instruments including the James Webb Space Telescope and Gemini Planet Imager.
The hundreds of exoplanets that have been discovered in the past two decades offer a new perspective on planetary structure. Instead of being the archetypal examples of planets, those of our Solar System are merely possible outcomes of planetary syst
em formation and evolution, and conceivably not even terribly common outcomes (although this remains an open question). Here, we review the diverse range of interior structures that are known to, and speculated to, exist in exoplanetary systems -- from mostly degenerate objects that are more than 10 times as massive as Jupiter, to intermediate-mass Neptune-like objects with large cores and moderate hydrogen/helium envelopes, to rocky objects with roughly the mass of the Earth.
Recent observations of the super-Earth GJ 1214b show that it has a relatively featureless transmission spectrum. One suggestion is that these observations indicate that the planets atmosphere is vertically compact, perhaps due to a water-rich composi
tion that yields a large mean molecular weight. Another suggestion is that the atmosphere is hydrogen/helium-rich with clouds that obscure predicted absorption features. Previous models that incorporate clouds have included their effect without a strong physical motivation for their existence. Here, we present model atmospheres of GJ 1214b that include physically-motivated clouds of two types. We model the clouds that form as a result of condensation in chemical equilibrium, as they likely do on brown dwarfs, which include KCl and ZnS for this planet. We also include clouds that form as a result of photochemistry, forming a hydrocarbon haze layer. We use a photochemical kinetics model to understand the vertical distribution and available mass of haze-forming molecules. We model both solar and enhanced-metallicity cloudy models and determine the cloud properties necessary to match observations. In enhanced-metallicity atmospheres, we find that the equilibrium clouds can match the observations of GJ 1214b if they are lofted high into the atmosphere and have a low sedimentation efficiency (fsed=0.1). We find that models with a variety of hydrocarbon haze properties can match the observations. Particle sizes from 0.01 to 0.25 micron can match the transmission spectrum with haze-forming efficiencies as low as 1-5%.
We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 micron bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 a
nd 24 micron, these data allow us to characterize the exoplanets emission spectrum as a function of planetary longitude. We utilize improved methods for removing the effects of intrapixel sensitivity variations and accounting for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% +/- 0.0061% in the 3.6 micron band and 0.0982% +/- 0.0089% in the 4.5 micron band. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 micron, and we present new evidence indicating that the flux minimum observed in the 8 micron is likely caused by an over-shooting effect in the 8 micron array. We obtain improved estimates for HD 189733bs dayside planet-star flux ratio of 0.1466% +/- 0.0040% at 3.6 micron and 0.1787% +/- 0.0038% at 4.5 micron; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of models from Burrows et al. (2008) and Showman et al. (2009). We find that HD 189733bs 4.5 micron nightside flux is 3.3 sigma smaller than predicted by the Showman et al. models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 micron absorption in this region. [abridged]
GJ436b is a unique member of the transiting extrasolar planet population being one of the smallest and least irradiated and possessing an eccentric orbit. Because of its size, mass and density, GJ436b could plausibly have an atmospheric metallicity s
imilar to Neptune (20-60 times solar abundances), which makes it an ideal target to study the effects of atmospheric metallicity on dynamics and radiative transfer in an extrasolar planetary atmosphere. We present three-dimensional atmospheric circulation models that include realistic non-gray radiative transfer for 1, 3, 10, 30, and 50 times solar atmospheric metallicity cases of GJ436b. Low metallicity models (1 and 3 times solar) show little day/night temperature variation and strong high-latitude jets. In contrast, higher metallicity models (30 and 50 times solar) exhibit day/night temperature variations and a strong equatorial jet. Spectra and light curves produced from these simulations show strong orbital phase dependencies in the 50 times solar case and negligible variations with orbital phase in the 1 times solar case. Comparisons between the predicted planet/star flux ratio from these models and current secondary eclipse measurements support a high metallicity atmosphere (30-50 times solar abundances) with disequilibrium carbon chemistry at play for GJ436b. Regardless of the actual atmospheric composition of GJ436b, our models serve to illuminate how metallicity influences the atmospheric circulation for a broad range of warm extrasolar planets.
The determination of Saturns atmospheric noble gas abundances are critical to understanding the formation and evolution of Saturn, and giant planets in general. These measurements can only be performed with an entry probe. A Saturn probe will address
whether enhancement in heavy noble gases, as was found in Jupiter, are a general feature of giant planets, and their ratios will be a powerful constraint on how they form. The helium abundance will show the extent to which helium has phase separated from hydrogen in the planets deep interior. Jupiters striking neon depletion may also be tied to its helium depletion, and must be confirmed or refuted in Saturn. Together with Jupiters measured atmospheric helium abundance, a consistent evolutionary theory for both planets, including helium rain will be possible. We will then be able to calibrate the theory of the evolution of all giant planets, including exoplanets. In addition, high pressure H/He mixtures under giant planet conditions are an important area of condensed matter physics that are beyond the realm of experiment.
We present global, three-dimensional numerical simulations of HD 189733b and HD 209458b that couple the atmospheric dynamics to a realistic representation of non-gray cloud-free radiative transfer. The model, which we call the Substellar and Planetar
y Atmospheric Radiation and Circulation (SPARC) model, adopts the MITgcm for the dynamics and uses the radiative model of McKay, Marley, Fortney, and collaborators for the radiation. Like earlier work with simplified forcing, our simulations develop a broad eastward equatorial jet, mean westward flow at higher latitudes, and substantial flow over the poles at low pressure. For HD 189733b, our simulations without TiO and VO opacity can explain the broad features of the observed 8 and 24-micron light curves, including the modest day-night flux variation and the fact that the planet/star flux ratio peaks before the secondary eclipse. Our simulations also provide reasonable matches to the Spitzer secondary-eclipse depths at 4.5, 5.8, 8, 16, and 24 microns and the groundbased upper limit at 2.2 microns. However, we substantially underpredict the 3.6-micron secondary-eclipse depth, suggesting that our simulations are too cold in the 0.1-1 bar region. Predicted temporal variability in secondary-eclipse depths is ~1% at Spitzer bandpasses, consistent with recent observational upper limits at 8 microns. We also show that nonsynchronous rotation can significantly alter the jet structure. For HD 209458b, we include TiO and VO opacity; these simulations develop a hot (>2000 K) dayside stratosphere. Despite this stratosphere, we do not reproduce current Spitzer photometry of this planet. Light curves in Spitzer bandpasses show modest phase variation and satisfy the observational upper limit on day-night phase variation at 8 microns. (abridged)
