FORWARD AND INVERSE MODELING OF THE EMISSION AND TRANSMISSION SPECTRUM OF GJ 436B: INVESTIGATING METAL ENRICHMENT, TIDAL HEATING, AND CLOUDS

The Neptune-mass GJ 436b is one of the most-studied transiting exoplanets with repeated measurements of both its thermal emission and transmission spectra. We build on previous studies to answer outstanding questions about this planet, including its potentially high metallicity and tidal heating of its interior. We present new observations of GJ 436b's thermal emission at 3.6 and 4.5 micron, which reduce uncertainties in estimates of GJ 436b's flux at those wavelengths and demonstrate consistency between Spitzer observations spanning more than 7 years. We analyze the Spitzer thermal emission photometry and Hubble WFC3 transmission spectrum in tandem. We use a powerful dual-pronged modeling approach, comparing these data to both self-consistent and retrieval models. We vary the metallicity, intrinsic luminosity from tidal heating, disequilibrium chemistry, and heat redistribution. We also study the effect of clouds and photochemical hazes on the spectra, but do not find strong evidence for either. The self-consistent and retrieval modeling combine to suggest that GJ 436b has a high atmospheric metallicity, with best fits at or above several hundred times solar metallicity, tidal heating warming its interior with best-fit intrinsic effective effective temperatures around 300--350 K, and disequilibrium chemistry. High metal-enrichments (>600x solar) can only occur from the accretion of rocky, rather than icy, material. Assuming Tint~300--350 K, we find that Q'~2x10^5--10^6, larger than Neptune's Q', and implying a long tidal circularization timescale for the planet's orbit. We suggest that Neptune-mass planets may be a more diverse class than previously imagined, with metal-enhancements potentially spanning several orders of magnitude, to perhaps over 1000x solar metallicity. High fidelity observations with instruments like JWST will be critical for characterizing this diversity.

[1]  Drake Deming,et al.  Possible thermochemical disequilibrium in the atmosphere of the exoplanet GJ 436b , 2010, Nature.

[2]  K. Cahoy,et al.  EXOPLANET ALBEDO SPECTRA AND COLORS AS A FUNCTION OF PLANET PHASE, SEPARATION, AND METALLICITY , 2010, 1009.3071.

[3]  J. Beaulieu,et al.  METHANE IN THE ATMOSPHERE OF THE TRANSITING HOT NEPTUNE GJ436B? , 2010, 1007.0324.

[4]  I. Pater,et al.  Constraining the origins of Neptune’s carbon monoxide abundance with CARMA millimeter-wave observations , 2012, 1301.1990.

[5]  Caltech,et al.  THE GJ 436 SYSTEM: DIRECTLY DETERMINED ASTROPHYSICAL PARAMETERS OF AN M DWARF AND IMPLICATIONS FOR THE TRANSITING HOT NEPTUNE , 2012, 1202.0083.

[6]  AZ,et al.  Characterization of the hot Neptune GJ 436 b with Spitzer and ground-based observations , 2007, 0707.3809.

[7]  M. Marley,et al.  THE ATMOSPHERIC CIRCULATION OF A NINE-HOT-JUPITER SAMPLE: PROBING CIRCULATION AND CHEMISTRY OVER A WIDE PHASE SPACE , 2016, 1602.06733.

[8]  A. Merloni,et al.  X-ray spectral modelling of the AGN obscuring region in the CDFS: Bayesian model selection and catalogue , 2014, 1402.0004.

[9]  Debra A. Fischer,et al.  A Neptune-Mass Planet Orbiting the Nearby M Dwarf GJ 436 , 2004 .

[10]  David Charbonneau,et al.  KEPLER-93b: A TERRESTRIAL WORLD MEASURED TO WITHIN 120 km, AND A TEST CASE FOR A NEW SPITZER OBSERVING MODE , 2014, 1405.3659.

[11]  Christoph Mordasini,et al.  A FRAMEWORK FOR CHARACTERIZING THE ATMOSPHERES OF LOW-MASS LOW-DENSITY TRANSITING PLANETS , 2013, 1306.4329.

[12]  C P McKay,et al.  Thermal structure of Uranus' atmosphere. , 1999, Icarus.

[13]  A. Warn-Varnas,et al.  Yellow Sea ocean-acoustic solitary wave modeling studies , 2005 .

[14]  Drake Deming,et al.  Spitzer Transit and Secondary Eclipse Photometry of GJ 436b , 2007, 0707.2778.

[15]  M. Marley,et al.  ATMOSPHERIC CIRCULATION OF ECCENTRIC HOT NEPTUNE GJ436b , 2010, 1007.2942.

[16]  M. Line,et al.  A SYSTEMATIC RETRIEVAL ANALYSIS OF SECONDARY ECLIPSE SPECTRA. III. DIAGNOSING CHEMICAL DISEQUILIBRIUM IN PLANETARY ATMOSPHERES , 2013, 1309.6679.

[17]  R. Paul Butler,et al.  Refined stellar, orbital and planetary parameters of the eccentric HAT‐P‐2 planetary system , 2009, 0908.1705.

[18]  M. Holman,et al.  A SEARCH FOR ADDITIONAL PLANETS IN THE NASA EPOXI OBSERVATIONS OF THE EXOPLANET SYSTEM GJ 436 , 2009, 0909.2875.

[19]  Yuk L. Yung,et al.  Vertical transport and photochemistry in the terrestrial mesosphere and lower thermosphere (50–120 km) , 1981 .

[20]  Robert T. Zellem,et al.  THE 4.5 μm FULL-ORBIT PHASE CURVE OF THE HOT JUPITER HD 209458b , 2014, 1405.5923.

[21]  M. Marley,et al.  QUANTITATIVELY ASSESSING THE ROLE OF CLOUDS IN THE TRANSMISSION SPECTRUM OF GJ 1214b , 2013, 1305.4124.

[22]  Kyle L. Luther,et al.  CHARACTERIZING TRANSITING EXOPLANET ATMOSPHERES WITH JWST , 2015, 1511.05528.

[23]  B. Fegley,et al.  ATMOSPHERIC CHEMISTRY IN GIANT PLANETS, BROWN DWARFS, AND LOW-MASS DWARF STARS. III. IRON, MAGNESIUM, AND SILICON , 2010, 1001.3639.

[24]  A. Burrows,et al.  Chemical Equilibrium Abundances in Brown Dwarf and Extrasolar Giant Planet Atmospheres , 1999 .

[25]  S. Seager,et al.  Clouds and chemistry: Ultracool dwarf atmospheric properties from optical and infrared colors , 2002 .

[26]  M. Marley,et al.  GASEOUS MEAN OPACITIES FOR GIANT PLANET AND ULTRACOOL DWARF ATMOSPHERES OVER A RANGE OF METALLICITIES AND TEMPERATURES , 2014, 1409.0026.

[27]  F. Selsis,et al.  THE PUZZLING CHEMICAL COMPOSITION OF GJ 436B'S ATMOSPHERE: INFLUENCE OF TIDAL HEATING ON THE CHEMISTRY , 2013, 1312.3007.

[28]  Drake Deming,et al.  A Search for a Sub-Earth-Sized Companion to GJ 436 and a Novel Method to Calibrate Warm Spitzer IRAC Observations , 2010, 1009.0755.

[29]  M. R. Line,et al.  INFORMATION CONTENT OF EXOPLANETARY TRANSIT SPECTRA: AN INITIAL LOOK , 2011, 1111.2612.

[30]  I. P. Waldmann,et al.  REVISITING SPITZER TRANSIT OBSERVATIONS WITH INDEPENDENT COMPONENT ANALYSIS: NEW RESULTS FOR THE GJ 436 SYSTEM , 2015, 1501.05866.

[31]  Richard Greenberg,et al.  Tidal Heating of Extrasolar Planets , 2008, 0803.0026.

[32]  Sara Seager,et al.  HELIUM ATMOSPHERES ON WARM NEPTUNE- AND SUB-NEPTUNE-SIZED EXOPLANETS AND APPLICATIONS TO GJ 436b , 2015, 1505.02221.

[33]  Claudia Emde,et al.  New secondary-scattering correction in DISORT with increased efficiency for forward scattering , 2011 .

[34]  Steven Soter,et al.  Q in the solar system , 1966 .

[35]  R. P. Butler,et al.  The M Dwarf GJ 436 and its Neptune‐Mass Planet , 2006, astro-ph/0608260.

[36]  C. McKay,et al.  The thermal structure of Titan's atmosphere. , 1989, Icarus.

[37]  R. Redmer,et al.  Interior structure models of GJ 436b , 2010, 1002.4447.

[38]  G. Laughlin,et al.  A QUASI-STATIONARY SOLUTION TO GLIESE 436b's ECCENTRICITY , 2009, 0904.3146.

[39]  Gilles Chabrier,et al.  An Equation of State for Low-Mass Stars and Giant Planets , 1995 .

[40]  E. Agol,et al.  3.6 AND 4.5 μm PHASE CURVES OF THE HIGHLY IRRADIATED ECCENTRIC HOT JUPITER WASP-14b , 2015, 1505.03158.

[41]  T. Barman,et al.  COMPOSITIONAL DIVERSITY IN THE ATMOSPHERES OF HOT NEPTUNES, WITH APPLICATION TO GJ 436b , 2013, The Astrophysical journal.

[42]  Carl Sagan,et al.  Physical properties of the particles composing the Martian dust storm of 1971–1972 , 1977 .

[43]  Andrew S. Ackerman,et al.  Precipitating Condensation Clouds in Substellar Atmospheres , 2001, astro-ph/0103423.

[44]  Yuk L. Yung,et al.  High-temperature Photochemistry in the Atmosphere of HD 189733b , 2010 .

[45]  M. Marley,et al.  Line and Mean Opacities for Ultracool Dwarfs and Extrasolar Planets , 2007, 0706.2374.

[46]  M. Holman,et al.  Transit infrared spectroscopy of the hot Neptune around GJ 436 with the Hubble Space Telescope , 2008, 0810.5731.

[47]  Portugal,et al.  Accurate Spitzer infrared radius measurement for the hot Neptune GJ 436b , 2007, 0707.2261.

[48]  G. Vasisht,et al.  THERMOCHEMICAL AND PHOTOCHEMICAL KINETICS IN COOLER HYDROGEN-DOMINATED EXTRASOLAR PLANETS: A METHANE-POOR GJ436b? , 2011, 1104.3183.

[49]  D. Saumon,et al.  The Evolution of L and T Dwarfs in Color-Magnitude Diagrams , 2008, 0808.2611.

[50]  Ignasi Ribas,et al.  A ~5 M⊕ Super-Earth Orbiting GJ 436? The Power of Near-Grazing Transits , 2008, 0801.3230.

[51]  F. Feroz,et al.  MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics , 2008, 0809.3437.

[52]  M. Marley,et al.  METHANE, CARBON MONOXIDE, AND AMMONIA IN BROWN DWARFS AND SELF-LUMINOUS GIANT PLANETS , 2014, 1408.6283.

[53]  Kerri Cahoy,et al.  THERMAL EMISSION AND REFLECTED LIGHT SPECTRA OF SUPER EARTHS WITH FLAT TRANSMISSION SPECTRA , 2015, 1511.01492.

[54]  S. L. Thompson.  ANEOS analytic equations of state for shock physics codes input manual , 1990 .

[55]  T. Guillot,et al.  A Nongray Theory of Extrasolar Giant Planets and Brown Dwarfs , 1997, astro-ph/9705201.

[56]  Adam Burrows,et al.  ALBEDO AND REFLECTION SPECTRA OF EXTRASOLAR GIANT PLANETS , 1999 .

[57]  Thomas P. Greene,et al.  TRANSMISSION SPECTRA OF TRANSITING PLANET ATMOSPHERES: MODEL VALIDATION AND SIMULATIONS OF THE HOT NEPTUNE GJ 436b FOR THE JAMES WEBB SPACE TELESCOPE , 2010, 1010.2451.

[58]  J. Fortney,et al.  THE MASS–METALLICITY RELATION FOR GIANT PLANETS , 2015, 1511.07854.

[59]  Sanford Gordon,et al.  Computer program for calculation of complex chemical equilibrium compositions , 1972 .

[60]  E. Agol,et al.  Analytic Light Curves for Planetary Transit Searches , 2002, astro-ph/0210099.

[61]  Comparative Planetary Atmospheres: Models of TrES-1 and HD 209458b , 2005, astro-ph/0505359.

[62]  France.,et al.  Limits to the planet candidate GJ 436c , 2008, 0804.3030.

[63]  M. Marley,et al.  Multiwaveband photometry of the irradiated brown dwarf WD0137−349B , 2014, 1412.6363.

[64]  Nikole K. Lewis,et al.  ORBITAL PHASE VARIATIONS OF THE ECCENTRIC GIANT PLANET HAT-P-2b , 2013, 1302.5084.

[65]  S. Hinkley,et al.  FRIENDS OF HOT JUPITERS. I. A RADIAL VELOCITY SEARCH FOR MASSIVE, LONG-PERIOD COMPANIONS TO CLOSE-IN GAS GIANT PLANETS , 2013, 1312.2954.

[66]  Nikole K. Lewis,et al.  DISEQUILIBRIUM CARBON, OXYGEN, AND NITROGEN CHEMISTRY IN THE ATMOSPHERES OF HD 189733b AND HD 209458b , 2011, 1102.0063.

[67]  D. Crisp,et al.  A SYSTEMATIC RETRIEVAL ANALYSIS OF SECONDARY ECLIPSE SPECTRA. I. A COMPARISON OF ATMOSPHERIC RETRIEVAL TECHNIQUES , 2013, 1304.5561.

[68]  Michael R. Line,et al.  THE INFLUENCE OF NONUNIFORM CLOUD COVER ON TRANSIT TRANSMISSION SPECTRA , 2015, 1511.09443.

[69]  D. Hamilton,et al.  Orbital resonances in the inner neptunian system. II. Resonant history of Proteus, Larissa, Galatea, and Despina , 2008 .

[70]  K. Stamnes,et al.  Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media. , 1988, Applied optics.

[71]  T. Guillot,et al.  A non-grey analytical model for irradiated atmospheres - I. Derivation , 2013, 1311.6597.

[72]  D. Saumon,et al.  NEGLECTED CLOUDS IN T AND Y DWARF ATMOSPHERES , 2012, 1206.4313.

[73]  S. Seager,et al.  HIGH METALLICITY AND NON-EQUILIBRIUM CHEMISTRY IN THE DAYSIDE ATMOSPHERE OF HOT-NEPTUNE GJ 436b , 2010, 1004.5121.

[74]  Heather Knutson,et al.  A SYSTEMATIC RETRIEVAL ANALYSIS OF SECONDARY ECLIPSE SPECTRA. II. A UNIFORM ANALYSIS OF NINE PLANETS AND THEIR C TO O RATIOS , 2013, 1309.6663.

[75]  Mark S. Marley,et al.  Planetary Radii across Five Orders of Magnitude in Mass and Stellar Insolation: Application to Transits , 2006 .

[76]  E. Observatory,et al.  High cadence near infrared timing observations of extrasolar planets , 2009, 0905.1728.

[77]  C. McKay,et al.  Rapid calculation of radiative heating rates and photodissociation rates in inhomogeneous multiple scattering atmospheres , 1989 .

[78]  M. Allen,et al.  Photochemistry of the atmosphere of Titan: comparison between model and observations. , 1984, The Astrophysical journal. Supplement series.

[79]  G. Montagnier,et al.  Dynamical evolution of the Gliese 436 planetary system - Kozai migration as a potential source for Gliese 436b's eccentricity , 2012, 1208.0237.

[80]  Drake Deming,et al.  A featureless transmission spectrum for the Neptune-mass exoplanet GJ 436b , 2014, Nature.

[81]  Drake Deming,et al.  A SPITZER TRANSMISSION SPECTRUM FOR THE EXOPLANET GJ 436b, EVIDENCE FOR STELLAR VARIABILITY, AND CONSTRAINTS ON DAYSIDE FLUX VARIATIONS , 2011, 1104.2901.

[82]  E. Agol,et al.  Spitzer Secondary Eclipse Observations of Five Cool Gas Giant Planets and Empirical Trends in Cool Planet Emission Spectra , 2015, 1508.00902.

[83]  Drake Deming,et al.  SPITZER SECONDARY ECLIPSES OF THE DENSE, MODESTLY-IRRADIATED, GIANT EXOPLANET HAT-P- 20 b ?> USING PIXEL-LEVEL DECORRELATION , 2014, 1411.7404.

[84]  T. Guillot,et al.  Atmospheric, Evolutionary, and Spectral Models of the Brown Dwarf Gliese 229 B , 1996, Science.

[85]  D. Queloz,et al.  A global analysis of Spitzer and new HARPS data confirms the loneliness and metal-richness of GJ 436 b , 2014, 1409.4038.

[86]  Aisey M Andel.  ANALYTIC LIGHTCURVES FOR PLANETARY TRANSIT SEARCHES , 2002 .

[87]  Mark S. Marley,et al.  Synthetic Spectra and Colors of Young Giant Planet Atmospheres: Effects of Initial Conditions and Atmospheric Metallicity , 2008, 0805.1066.

[88]  J. Lunine,et al.  Reflected Spectra and Albedos of Extrasolar Giant Planets. I. Clear and Cloudy Atmospheres , 1998, astro-ph/9810073.

[89]  Jacob L. Bean,et al.  A DETECTION OF WATER IN THE TRANSMISSION SPECTRUM OF THE HOT JUPITER WASP-12b AND IMPLICATIONS FOR ITS ATMOSPHERIC COMPOSITION , 2015, 1504.05586.

[90]  Christopher J. Campo,et al.  TWO NEARBY SUB-EARTH-SIZED EXOPLANET CANDIDATES IN THE GJ 436 SYSTEM , 2012, 1207.4245.

[91]  D. Queloz,et al.  Detection of transits of the nearby hot Neptune GJ 436 b , 2007, Astronomy & Astrophysics.