Showing 1–50 of 67 results for author: Freedman, R

High-precision atmospheric characterization of a Y dwarf with JWST NIRSpec G395H spectroscopy: isotopologue, C/O ratio, metallicity, and the abundances of six molecular species

Authors: Ben W. P. Lew, Thomas Roellig, Natasha E. Batalha, Michael Line, Thomas Greene, Sagnick Murkherjee, Richard Freedman, Michael Meyer, Charles Beichman, Catarina Alves De Oliveira, Matthew De Furio, Doug Johnstone, Alexandra Z. Greenbaum, Mark Marley, Jonathan J. Fortney, Erick T. Young, Jarron Leisenring, Martha Boyer, Klaus Hodapp, Karl Misselt, John Stansberry, Marcia Rieke

Abstract: The launch of the James Webb Space Telescope (JWST) marks a pivotal moment for precise atmospheric characterization of Y dwarfs, the coldest brown dwarf spectral type. In this study, we leverage moderate spectral resolution observations (R $\sim$ 2700) with the G395H grating of the Near-Infrared Spectrograph (NIRSpec) onboard of JWST to characterize the nearby (9.9 pc) Y dwarf WISEPA J182831.08+26… ▽ More The launch of the James Webb Space Telescope (JWST) marks a pivotal moment for precise atmospheric characterization of Y dwarfs, the coldest brown dwarf spectral type. In this study, we leverage moderate spectral resolution observations (R $\sim$ 2700) with the G395H grating of the Near-Infrared Spectrograph (NIRSpec) onboard of JWST to characterize the nearby (9.9 pc) Y dwarf WISEPA J182831.08+265037.8 (WISE 1828). With the NIRSpec G395H 2.88-5.12 $\mathrmμ$m spectrum, we measure the abundances of CO, CO$_2$, CH$_4$, H$_2$S, NH$_3$, and H$_2$O, which are the major carbon, nitrogen, oxygen, and sulfur bearing species in the atmosphere. Based on the retrieved volume mixing ratios with the atmospheric retrieval framework CHIMERA, we report that the C/O ratio is $0.45 \pm 0.01$, close to the solar C/O value of 0.55, and the metallicity to be +0.30 $\pm$ 0.02 dex. Comparison between the retrieval results with the forward modeling results suggests that the model bias for C/O and metallicity could be as high as 0.03 and 0.97 dex respectively. We also report a lower limit of the $^{12}$CO/$^{13}$CO ratio of $>40 $, being consistent with the nominal solar value of 90. Our results highlight the potential of JWST in measuring the C/O ratios down to percent-level precision and characterizing isotopologues of cold planetary atmospheres similar to WISE 1828. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 18 pages + references, including 11 figures, accepted for publication in The Astronomical Journal

The Sonora Substellar Atmosphere Models. III. Diamondback: Atmospheric Properties, Spectra, and Evolution for Warm Cloudy Substellar Objects

Authors: Caroline V. Morley, Sagnick Mukherjee, Mark S. Marley, Jonathan J. Fortney, Channon Visscher, Roxana Lupu, Ehsan Gharib-Nezhad, Daniel Thorngren, Richard Freedman, Natasha Batalha 7

Abstract: We present a new grid of cloudy atmosphere and evolution models for substellar objects. These models include the effect of refractory cloud species, including silicate clouds, on the spectra and evolution. We include effective temperatures from 900 to 2400 K and surface gravities from log g=3.5-5.5, appropriate for a broad range of objects with masses between 1 and 84 Jupiter masses. Model pressur… ▽ More We present a new grid of cloudy atmosphere and evolution models for substellar objects. These models include the effect of refractory cloud species, including silicate clouds, on the spectra and evolution. We include effective temperatures from 900 to 2400 K and surface gravities from log g=3.5-5.5, appropriate for a broad range of objects with masses between 1 and 84 Jupiter masses. Model pressure-temperature structures are calculated assuming radiative-convective and chemical equilibrium. We consider the effect of both clouds and metallicity on the atmospheric structure, resulting spectra, and thermal evolution of substellar worlds. We parameterize clouds using the Ackerman & Marley (2001) cloud model, including cloud parameter fsed values from 1-8; we include three metallicities (-0.5, 0.0, and +0.5). Refractory clouds and metallicity both alter the evolution of substellar objects, changing the inferred temperature at a given age by up to 100-200 K. We compare to the observed photometry of brown dwarfs, finding broad agreement with the measured photometry. We publish the spectra, evolution, and other data products online with open access. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 26 pages, 18 figures, submitted for publication in ApJ

The Sonora Substellar Atmosphere Models. IV. Elf Owl: Atmospheric Mixing and Chemical Disequilibrium with Varying Metallicity and C/O Ratios

Authors: Sagnick Mukherjee, Jonathan J. Fortney, Caroline V. Morley, Natasha E. Batalha, Mark S. Marley, Theodora Karalidi, Channon Visscher, Roxana Lupu, Richard Freedman, Ehsan Gharib-Nezhad

Abstract: Disequilibrium chemistry due to vertical mixing in the atmospheres of many brown dwarfs and giant exoplanets is well-established. Atmosphere models for these objects typically parameterize mixing with the highly uncertain $K_{\rm zz}$ diffusion parameter. The role of mixing in altering the abundances of C-N-O-bearing molecules has mostly been explored for solar composition atmospheres. However, at… ▽ More Disequilibrium chemistry due to vertical mixing in the atmospheres of many brown dwarfs and giant exoplanets is well-established. Atmosphere models for these objects typically parameterize mixing with the highly uncertain $K_{\rm zz}$ diffusion parameter. The role of mixing in altering the abundances of C-N-O-bearing molecules has mostly been explored for solar composition atmospheres. However, atmospheric metallicity and the C/O ratio also impact atmospheric chemistry. Therefore, we present the \texttt{Sonora Elf Owl} grid of self-consistent cloud-free 1D radiative-convective equilibrium model atmospheres for JWST observations, which includes a variation of $K_{\rm zz}$ across several orders of magnitude and also encompasses sub-solar to super-solar metallicities and C/O ratios. We find that the impact of $K_{\rm zz}$ on the $T(P)$ profile and spectra is a strong function of both $T_{\rm eff}$ and metallicity. For metal-poor objects $K_{\rm zz}$ has large impacts on the atmosphere at significantly higher $T_{\rm eff}$ compared to metal-rich atmospheres where the impact of $K_{\rm zz}$ is seen to occur at lower $T_{\rm eff}$. We identify significant spectral degeneracies between varying $K_{\rm zz}$ and metallicity in multiple wavelength windows, in particular at 3-5 $μ$m. We use the \texttt{Sonora Elf Owl} atmospheric grid to fit the observed spectra of a sample of 9 early to late T- type objects from $T_{\rm eff}=550-1150$ K. We find evidence for very inefficient vertical mixing in these objects with inferred $K_{\rm zz}$ values lying in the range between $\sim$ 10$^1$-10$^4$ cm$^2$s$^{-1}$. Using self-consistent models, we find that this slow vertical mixing is due to the observations probing mixing in the deep detached radiative zone in these atmospheres. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: Accepted for Publication in The Astrophysical Journal, 16 Figures, 3 Tables, 28 Pages

The impact of spectral line wing cut-off: Recommended standard method with application to MAESTRO opacity database

Authors: Ehsan, Gharib-Nezhad, Natasha E. Batalha, Katy Chubb, Richard Freedman, Iouli E. Gordon, Robert R. Gamache, Robert J. Hargreaves, Nikole K. Lewis, Jonathan Tennyson, Sergei N. Yurchenko

Abstract: When computing cross-sections from a line list, the result depends not only on the line strength, but also the line shape, pressure-broadening parameters, and line wing cut-off (i.e., the maximum distance calculated from each line centre). Pressure-broadening can be described using the Lorentz lineshape, but it is known to not represent the true absorption in the far wings. Both theory and experim… ▽ More When computing cross-sections from a line list, the result depends not only on the line strength, but also the line shape, pressure-broadening parameters, and line wing cut-off (i.e., the maximum distance calculated from each line centre). Pressure-broadening can be described using the Lorentz lineshape, but it is known to not represent the true absorption in the far wings. Both theory and experiment have shown that far from the line centre, non-Lorentzian behaviour controls the shape of the wings and the Lorentz lineshape fails to accurately characterize the absorption, leading to an underestimation or overestimation of the opacity continuum depending on the molecular species involved. The line wing cut-off is an often overlooked parameter when calculating absorption cross sections, but can have a significant effect on the appearance of the spectrum since it dictates the extent of the line wing that contributes to the calculation either side of every line centre. Therefore, when used to analyse exoplanet and brown dwarf spectra, an inaccurate choice for the line wing cut-off can result in errors in the opacity continuum, which propagate into the modeled transit spectra, and ultimately impact/bias the interpretation of observational spectra, and the derived composition and thermal structure. Here, we examine the different methods commonly utilized to calculate the wing cut-off and propose a standard practice procedure (i.e., absolute value of 25~cm$^{-1}$ for $P\leqslant$~200~bar and 100~cm$^{-1}$ for $P >$ ~200~bar) to generate molecular opacities which will be used by the open-access {\tt MAESTRO} (Molecules and Atoms in Exoplanet Science: Tools and Resources for Opacities) database. The pressing need for new measurements and theoretical studies of the far-wings is highlighted. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: 13 pages, 6 figures

Journal ref: RAS Techniques and Instruments, 2023, rzad058

The Missing Link: Testing Galactic Chemical Evolution Models with the First Multi-Isotopic Abundances in Solar Twin Stars

Authors: David R. Coria, Ian J. M. Crossfield, Joshua Lothringer, Becky Flores, Nikos Prantzos, Richard Freedman

Abstract: We present the first isotopic abundances of both $^{13}$CO and C$^{18}$O in solar twin stars and test the results against several galactic chemical evolution (GCE) models with different nucleosynthesis prescriptions. First, we compare M-band spectra from IRTF/iSHELL to synthetic spectra generated from custom solar atmosphere models using the PHOENIX atmosphere code. Next, we compare our calculated… ▽ More We present the first isotopic abundances of both $^{13}$CO and C$^{18}$O in solar twin stars and test the results against several galactic chemical evolution (GCE) models with different nucleosynthesis prescriptions. First, we compare M-band spectra from IRTF/iSHELL to synthetic spectra generated from custom solar atmosphere models using the PHOENIX atmosphere code. Next, we compare our calculated abundances to GCE models that consider isotopic yields from massive stars, asymptotic giant branch (AGB) stars and fast-rotating stars. The $^{12}$C/$^{13}$C ratios determined for this sample of solar twins are consistent with predictions from the selected GCE models; however, the $^{16}$O/$^{18}$O ratios tentatively contradict these predictions. This project constitutes the first in a stellar chemical abundance series seeking to: (1) support the James Webb Space Telescope (JWST) as it characterizes exoplanet atmospheres, interiors, and biosignatures by providing host star abundances (2) identify how unexplored stellar abundances reveal the process of galactic chemical evolution and correlate with star formation, interior, age, metallicity, and activity; and (3) provide improved stellar ages using stellar abundance measurements. By measuring elemental and isotopic abundances in a variety of stars, we not only supply refined host star parameters, but also provide the necessary foundations for complementary exoplanet characterization studies and ultimately contribute to the exploration of galactic, stellar, and planetary origins and evolution. △ Less

Submitted 30 August, 2023; originally announced August 2023.

Comments: 21 pages, 6 figures. Published in ApJ 8/29/2023

Journal ref: The Astrophysical Journal (ApJ 2023) Volume 954 Number 2 Page 121

The puzzle of the formation of T8 dwarf Ross 458c

Authors: Josefine Gaarn, Ben Burningham, Jacqueline K. Faherty, Channon Visscher, Mark S. Marley, Eileen C. Gonzales, Emily Calamari, Daniella Bardalez Gagliuffi, Roxana Lupu, Richard Freedman

Abstract: At the lowest masses, the distinction between brown dwarfs and giant exoplanets is often blurred and literature classifications rarely reflect the deuterium burning boundary. Atmospheric characterisation may reveal the extent to which planetary formation pathways contribute to the population of very-low mass brown dwarfs, by revealing if their abundance distributions differ from those of the local… ▽ More At the lowest masses, the distinction between brown dwarfs and giant exoplanets is often blurred and literature classifications rarely reflect the deuterium burning boundary. Atmospheric characterisation may reveal the extent to which planetary formation pathways contribute to the population of very-low mass brown dwarfs, by revealing if their abundance distributions differ from those of the local field population or, in the case of companions, their primary stars. The T8 dwarf Ross 458c is a possible planetary mass companion to a pair of M dwarfs, and previous work suggests that it is cloudy. We here present the results of the retrieval analysis of Ross 458c, using archival spectroscopic data in the 1.0 to 2.4 micron range. We test a cloud free model as well as a variety of cloudy models and find that the atmosphere of Ross 458c is best described by a cloudy model (strongly preferred). The CH4/H2O is higher than expected at 1.97 +0.13 -0.14. This value is challenging to understand in terms of equilibrium chemistry and plausible C/O ratios. Comparisons to thermochemical grid models suggest a C/O of ~ 1.35, if CH4 and H2O are quenched at 2000 K, requiring vigorous mixing. We find a [C/H] ratio of +0.18, which matches the metallicity of the primary system, suggesting that oxygen is missing from the atmosphere. Even with extreme mixing, the implied C/O is well beyond the typical stellar regime, suggesting a either non-stellar formation pathway, or the sequestration of substantial quantities of oxygen via hitherto unmodeled chemistry or condensation processes. △ Less

Submitted 3 April, 2023; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: 16 pages, 9 figures

The Sonora Substellar Atmosphere Models. II. Cholla: A Grid of Cloud-free, Solar Metallicity Models in Chemical Disequilibrium for the JWST Era

Authors: Theodora Karalidi, Mark Marley, Jonathan J. Fortney, Caroline Morley, Didier Saumon, Roxana Lupu, Channon Visscher, Richard Freedman

Abstract: Exoplanet and brown dwarf atmospheres commonly show signs of disequilibrium chemistry. In the James Webb Space Telescope era high resolution spectra of directly imaged exoplanets will allow the characterization of their atmospheres in more detail, and allow systematic tests for the presence of chemical species that deviate from thermochemical equilibrium in these atmospheres. Constraining the pres… ▽ More Exoplanet and brown dwarf atmospheres commonly show signs of disequilibrium chemistry. In the James Webb Space Telescope era high resolution spectra of directly imaged exoplanets will allow the characterization of their atmospheres in more detail, and allow systematic tests for the presence of chemical species that deviate from thermochemical equilibrium in these atmospheres. Constraining the presence of disequilibrium chemistry in these atmospheres as a function of parameters such as their effective temperature and surface gravity will allow us to place better constrains in the physics governing these atmospheres. This paper is part of a series of works presenting the Sonora grid of atmosphere models (Marley et al 2021, Morley et al in prep.). In this paper we present a grid of cloud-free, solar metallicity atmospheres for brown dwarfs and wide separation giant planets with key molecular species such as CH4, H2O, CO and NH3 in disequilibrium. Our grid covers atmospheres with Teff~[500 K,1300 K], logg~[3.0,5.5] (cgs) and an eddy diffusion parameter of logKzz=2, 4 and 7 (cgs). We study the effect of different parameters within the grid on the temperature and composition profiles of our atmospheres. We discuss their effect on the near-infrared colors of our model atmospheres and the detectability of CH4, H2O, CO and NH3 using the JWST. We compare our models against existing MKO and Spitzer observations of brown dwarfs and verify the importance of disequilibrium chemistry for T dwarf atmospheres. Finally, we discuss how our models can help constrain the vertical structure and chemical composition of these atmospheres. △ Less

Submitted 22 October, 2021; originally announced October 2021.

Comments: accepted to ApJ

The first retrieval of a substellar subdwarf: A cloud-free SDSS J125637.13-022452.4

Authors: Eileen C. Gonzales, Ben Burningham, Jacqueline K. Faherty, Channon Visscher, Mark Marley, Roxana Lupu, Richard Freedman, Nikole K. Lewis

Abstract: We present the first retrieval analysis of a substellar subdwarf, SDSS J125637.13-022452.4 (SDSS J1256-0224), using the Brewster retrieval code base. We find SDSS J1256-0224 is best fit by a cloud-free model with an ion (neutral H, H-, and electron) abundance corresponding to ion [Fe/H]=-1.5. However, this model is indistinguishable from a cloud-free model with ion [Fe/H]=-2.0 and a cloud-free mod… ▽ More We present the first retrieval analysis of a substellar subdwarf, SDSS J125637.13-022452.4 (SDSS J1256-0224), using the Brewster retrieval code base. We find SDSS J1256-0224 is best fit by a cloud-free model with an ion (neutral H, H-, and electron) abundance corresponding to ion [Fe/H]=-1.5. However, this model is indistinguishable from a cloud-free model with ion [Fe/H]=-2.0 and a cloud-free model with ion Fe/H]=-1.5 assuming a subsolar carbon-to-oxygen ratio. We are able to constrain abundances for water, FeH, and CrH, with an inability to constrain any carbon-bearing species likely due to the low-metallicity of SDSS J1256-0224. We also present an updated spectral energy distribution (SED) and semi-empirical fundamental parameters. Our retrieval- and SED-based fundamental parameters agree with the Baraffe low-metallicity evolutionary models. From examining our "rejected" models (those with $Δ$BIC>45), we find that we are able to retrieve gas abundances consistent with those of our best-fitting model. We find the cloud in these poorer fitting "cloudy" models is either pushed to the bottom of the atmosphere or made optically thin. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: 19 pages, 5 figures. arXiv admin note: text overlap with arXiv:2010.01224

The Sonora Brown Dwarf Atmosphere and Evolution Models I. Model Description and Application to Cloudless Atmospheres in Rainout Chemical Equilibrium

Authors: Mark S. Marley, Didier Saumon, Channon Visscher, Roxana Lupu, Richard Freedman, Caroline Morley, Jonathan J. Fortney, Christopher Seay, Adam J. R. W. Smith, D. J. Teal, Ruoyan Wang

Abstract: We present a new generation of substellar atmosphere and evolution models, appropriate for application to studies of L, T, and Y-type brown dwarfs and self-luminous extrasolar planets. The atmosphere models describe the expected temperature-pressure profiles and emergent spectra of atmospheres in radiative-convective equilibrium with effective temperatures and gravities within the ranges… ▽ More We present a new generation of substellar atmosphere and evolution models, appropriate for application to studies of L, T, and Y-type brown dwarfs and self-luminous extrasolar planets. The atmosphere models describe the expected temperature-pressure profiles and emergent spectra of atmospheres in radiative-convective equilibrium with effective temperatures and gravities within the ranges $200\le T_{\rm eff}\le2400\,\rm K$ and $2.5\le \log g \le 5.5$. These ranges encompass masses from about 0.5 to 85 Jupiter masses for a set of metallicities ($[{\rm M/H}] = -0.5$ to $+0.5$), C/O ratios (from 0.5 to 1.5 times that of solar), and ages. The evolution tables describe the cooling of these substellar objects through time. These models expand the diversity of model atmospheres currently available, notably to cooler effective temperatures and greater ranges in C/O. Notable improvements from past such models include updated opacities and atmospheric chemistry. Here we describe our modeling approach and present our initial tranche of models for cloudless, chemical equilibrium atmospheres. We compare the modeled spectra, photometry, and evolution to various datasets. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 27 pages, 16 figures, accepted for Astrophysical Journal. Models available at https://doi.org/10.5281/zenodo.5063476

Detection and Bulk Properties of the HR 8799 Planets with High Resolution Spectroscopy

Authors: Jason J. Wang, Jean-Baptiste Ruffio, Evan Morris, Jacques-Robert Delorme, Nemanja Jovanovic, Jacklyn Pezzato, Daniel Echeverri, Luke Finnerty, Callie Hood, J. J. Zanazzi, Marta L. Bryan, Charlotte Z. Bond, Sylvain Cetre, Emily C. Martin, Dimitri Mawet, Andy Skemer, Ashley Baker, Jerry W. Xuan, J. Kent Wallace, Ji Wang, Randall Bartos, Geoffrey A. Blake, Andy Boden, Cam Buzard, Benjamin Calvin , et al. (27 additional authors not shown)

Abstract: Using the Keck Planet Imager and Characterizer (KPIC), we obtained high-resolution (R$\sim$35,000) $K$-band spectra of the four planets orbiting HR 8799. We clearly detected \water{} and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and \water{} in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral reso… ▽ More Using the Keck Planet Imager and Characterizer (KPIC), we obtained high-resolution (R$\sim$35,000) $K$-band spectra of the four planets orbiting HR 8799. We clearly detected \water{} and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and \water{} in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral resolution to date when considering both their angular separations and flux ratios. We developed a forward modeling framework that allows us to jointly fit the spectra of the planets and the diffracted starlight simultaneously in a likelihood-based approach and obtained posterior probabilities on their effective temperatures, surface gravities, radial velocities, and spins. We measured $v\sin(i)$ values of $10.1^{+2.8}_{-2.7}$~km/s for HR 8799 d and $15.0^{+2.3}_{-2.6}$~km/s for HR 8799 e, and placed an upper limit of $< 14$~km/s of HR 8799 c. Under two different assumptions of their obliquities, we found tentative evidence that rotation velocity is anti-correlated with companion mass, which could indicate that magnetic braking with a circumplanetary disk at early times is less efficient at spinning down lower mass planets. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 31 pages, 12 figures, Accepted to AJ

Following the Lithium: Tracing Li-bearing Molecules Across Age, Mass, and Gravity in Brown Dwarfs

Authors: Ehsan Gharib-Nezhad, Mark S. Marley, Natasha E. Batalha, Channon Visscher, Richard S. Freedman, Roxana E. Lupu

Abstract: Lithium is an important element for the understanding of ultracool dwarfs because it is lost to fusion at masses above $\sim 68\, M_{\rm J}$. Hence, the presence or absence of atomic Li has served as an indicator of the nearby H-burning boundary at about $75\,M_{\rm J}$ between brown-dwarfs and very low-mass stars. Historically the "Lithium test", a search for the presence and strength of the Li l… ▽ More Lithium is an important element for the understanding of ultracool dwarfs because it is lost to fusion at masses above $\sim 68\, M_{\rm J}$. Hence, the presence or absence of atomic Li has served as an indicator of the nearby H-burning boundary at about $75\,M_{\rm J}$ between brown-dwarfs and very low-mass stars. Historically the "Lithium test", a search for the presence and strength of the Li line at 670.8 nm, has been a marker if an object has a substellar mass with stellar-like spectral energy distribution (e.g., a late-type M dwarf). While the Li test could in principle also be used to distinguish masses of later-type L-T dwarfs, Li is predominantly no longer found as an atomic gas, but rather a molecular species such as LiH, LiF, LiOH, and LiCl in their cooler atmospheres. L- and T-type brown dwarfs are also quite faint at 670 nm and thus challenging targets for high resolution spectroscopy. But only recently have experimental molecular line lists become available for the molecular Li species, allowing molecular Li mass discrimination. In this study, we generated the latest opacity of each of these Li-bearing molecules and performed thermochemical equilibrium atmospheric composition calculation of the abundance of these molecules. Finally, we computed thermal emission spectra for a series of radiative-convective equilibrium models of cloudy and cloudless brown dwarf atmospheres (with $T_{\rm eff}=$ 500--2400~K, and $\log g$=4.0, 4.5, 5.0) to understand where the presence or absence of atmospheric lithium-bearing species is most easily detected as a function of brown dwarf mass and age. After atomic Li, the best spectral signatures were found to be LiF at $10.5-12.5$~\micron and LiCl at $14.5-18.5$ $\micron$. LiH also shows a narrow feature at $\sim 9.38$ $\micron$. △ Less

Submitted 1 June, 2021; originally announced June 2021.

Comments: 13 pages, 7 figures, 1 table, In final revision at ApJ. Comments welcome

Cloud busting: enstatite and quartz clouds in the atmosphere of 2M2224-0158

Authors: Ben Burningham, Jacqueline K. Faherty, Eileen C. Gonzales, Mark S. Marley, Channon Visscher, Roxana Lupu, Josefine Gaarn, Michelle Fabienne Bieger, Richard Freedman, Didier Saumon

Abstract: We present the most detailed data-driven exploration of cloud opacity in a substellar object to-date. We have tested over 60 combinations of cloud composition and structure, particle size distribution, scattering model, and gas phase composition assumptions against archival $1-15 {\rm μm}$ spectroscopy for the unusually red L4.5~dwarf 2MASSW~J2224438-015852 using the Brewster retrieval framework.… ▽ More We present the most detailed data-driven exploration of cloud opacity in a substellar object to-date. We have tested over 60 combinations of cloud composition and structure, particle size distribution, scattering model, and gas phase composition assumptions against archival $1-15 {\rm μm}$ spectroscopy for the unusually red L4.5~dwarf 2MASSW~J2224438-015852 using the Brewster retrieval framework. We find that, within our framework, a model that includes enstatite and quartz cloud layers at shallow pressures, combined with a deep iron cloud deck fits the data best. This models assumes a Hansen distribution for particle sizes for each cloud, and Mie scattering. We retrieved particle effective radii of $\log_{10} a {\rm (μm)} = -1.41^{+0.18}_{-0.17}$ for enstatite, $-0.44^{+0.04}_{-0.20}$ for quartz, and $-0.77^{+0.05}_{-0.06}$ for iron. Our inferred cloud column densities suggest ${\rm (Mg/Si)} = 0.69^{+0.06}_{-0.08}$ if there are no other sinks for magnesium or silicon. Models that include forsterite alongside, or in place of, these cloud species are strongly rejected in favour of the above combination. We estimate a radius of $0.75 \pm 0.02$ Rjup, which is considerably smaller than predicted by evolutionary models for a field age object with the luminosity of 2M2224-0158. Models which assume vertically constant gas fractions are consistently preferred over models that assume thermochemical equilibrium. From our retrieved gas fractions we infer ${\rm [M/H]} = +0.38^{+0.07}_{-0.06}$ and ${\rm C/O} = 0.83^{+0.06}_{-0.07}$. Both these values are towards the upper end of the stellar distribution in the Solar neighbourhood, and are mutually consistent in this context. A composition toward the extremes of the local distribution is consistent with this target being an outlier in the ultracool dwarf population. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 20 pages, 10 figures, MNRAS accepted

EXOPLINES: Molecular Absorption Cross-Section Database for Brown Dwarf and Giant Exoplanet Atmospheres

Authors: Ehsan Gharib-Nezhad, Aishwarya R. Iyer, Michael R. Line, Richard S. Freedman, Mark S. Marley, Natasha E. Batalha

Abstract: Stellar, substellar, and planetary atmosphere models are all highly sensitive to the input opacities. Generational differences between various state-of-the-art stellar/planetary models are primarily because of incomplete and outdated atomic/molecular line-lists. Here we present a database of pre-computed absorption cross-sections for all isotopologues of key atmospheric molecules relevant to late-… ▽ More Stellar, substellar, and planetary atmosphere models are all highly sensitive to the input opacities. Generational differences between various state-of-the-art stellar/planetary models are primarily because of incomplete and outdated atomic/molecular line-lists. Here we present a database of pre-computed absorption cross-sections for all isotopologues of key atmospheric molecules relevant to late-type stellar, brown dwarf, and planetary atmospheres: MgH, AlH, CaH, TiH, CrH, FeH, SiO, TiO, VO, and H2O. The pressure and temperature ranges of the computed opacities are between 10$^{-6}$--3000~bar and 75--4000~K, and their spectral ranges are 0.25--330~$μ$m for many cases where possible. For cases with no pressure-broadening data, we use collision theory to bridge the gap. We also probe the effect of absorption cross-sections calculated from different line lists in the context of Ultra-Hot Jupiter and M-dwarf atmospheres. Using 1-D self-consistent radiative-convective thermochemical equilibrium models, we report significant variations in the theoretical spectra and thermal profiles of substellar atmospheres. With a 2000 K representative Ultra-Hot Jupiter, we report variations of up to 320 and 80 ppm in transmission and thermal emission spectra, respectively. For a 3000 K M-dwarf, we find differences of up to 125$\%$ in the spectra. We find that the most significant differences arise due to the choice of TiO line-lists, primarily below 1$μ$m. In sum, we present (1) a database of pre-computed molecular absorption cross-sections, and (2) quantify biases that arise when characterizing substellar/exoplanet atmospheres due to line list differences, therefore highlighting the importance of correct and complete opacities for eventual applications to high precision spectroscopy and photometry. △ Less

Submitted 1 April, 2021; originally announced April 2021.

Comments: accepted in ApJS. 44 pages, 17 figures, 7 tables. Comments are welcome

Retrieval of SDSS J1416+1348AB

Authors: Eileen Gonzales, Ben Burningham, Jackie Faherty, Colleen Cleary, Channon Visscher, Mark Marley, Roxana Lupu, Richard Freedman

Abstract: We present the distance-calibrated spectral energy distribution (SED) of the d/sdL7 SDSS J14162408+1348263A (J1416A) and an updated SED for SDSS J14162408+1348263B (J1416B). We also present the first retrieval analysis of J1416A using the Brewster retrieval code base and the second retrieval of J1416B. We find that the primary is best fit by a non-grey cloud opacity with a power-law wavelength dep… ▽ More We present the distance-calibrated spectral energy distribution (SED) of the d/sdL7 SDSS J14162408+1348263A (J1416A) and an updated SED for SDSS J14162408+1348263B (J1416B). We also present the first retrieval analysis of J1416A using the Brewster retrieval code base and the second retrieval of J1416B. We find that the primary is best fit by a non-grey cloud opacity with a power-law wavelength dependence, but is indistinguishable between the type of cloud parameterization. J1416B is best fit by a cloud-free model, consistent with the results from Line et al. (2017). Most fundamental parameters derived via SEDs and retrievals are consistent within 1 sigma for both J1416A and J1416B. The exceptions include the radius of J1416A, where the retrieved radius is smaller than the evolutionary model-based radius from the SED for the deck cloud model, and the bolometric luminosity which is consistent within 2.5 sigma for both cloud models. The pair's metallicity and Carbon-to-Oxygen (C/O) ratio point towards formation and evolution as a system. By comparing the retrieved alkali abundances while using two opacity models, we are able to evaluate how the opacities behave for the L and T dwarf. Lastly, we find that relatively small changes in composition can drive major observable differences for lower temperature objects. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: 40 pages, 25 figures

Beyond Equilibrium Temperature: How the Atmosphere/Interior Connection Affects the Onset of Methane, Ammonia, and Clouds in Warm Transiting Giant Planets

Authors: Jonathan J. Fortney, Channon Visscher, Mark S. Marley, Callie E. Hood, Michael R. Line, Daniel P. Thorngren, Richard S. Freedman, Roxana Lupu

Abstract: The atmospheric pressure-temperature profiles for transiting giant planets cross a range of chemical transitions. Here we show that the particular shape of these irradiated profiles for warm giant planets below 1300 K lead to striking differences in the behavior of non-equilibrium chemistry compared to brown dwarfs of similar temperatures. Our particular focus is H$_2$O, CO, CH$_4$, CO$_2$, and NH… ▽ More The atmospheric pressure-temperature profiles for transiting giant planets cross a range of chemical transitions. Here we show that the particular shape of these irradiated profiles for warm giant planets below 1300 K lead to striking differences in the behavior of non-equilibrium chemistry compared to brown dwarfs of similar temperatures. Our particular focus is H$_2$O, CO, CH$_4$, CO$_2$, and NH$_3$ in Jupiter- and Neptune-class planets. We show the cooling history of a planet, which depends most significantly on planetary mass and age, can have a dominant effect on abundances in the visible atmosphere, often swamping trends one might expect based on Teq alone. The onset of detectable CH$_4$ in spectra can be delayed to lower Teq for some planets compared to equilibrium, or pushed to higher Teq. The detectability of NH$_3$ is typically enhanced compared to equilibrium expectations, which is opposite to the brown dwarf case. We find that both CH$_4$ and NH$_3$ can become detectable at around the same Teq (at Teq values that vary with mass and metallicity) whereas these "onset" temperatures are widely spaced for brown dwarfs. We suggest observational strategies to search for atmospheric trends and stress that non-equilibrium chemistry and clouds can serve as probes of atmospheric physics. As examples of atmospheric complexity, we assess three Neptune-class planets GJ 436b, GJ 3470b, and WASP-107, all around Teq=700 K. Tidal heating due to eccentricity damping in all three planets heats the deep atmosphere by thousands of degrees, and may explain the absence of CH$_4$ in these cool atmospheres. Atmospheric abundances must be interpreted in the context of physical characteristics of the planet. △ Less

Submitted 30 October, 2020; v1 submitted 30 September, 2020; originally announced October 2020.

Comments: Accepted to AJ. No additional significant changes

Prospects for Characterizing the Haziest Sub-Neptune Exoplanets with High Resolution Spectroscopy

Authors: Callie E. Hood, Jonathan J. Fortney, Michael R. Line, Emily C. Martin, Caroline V. Morley, Jayne L. Birkby, Zafar Rustamkulov, Roxana E. Lupu, Richard S. Freedman

Abstract: Observations to characterize planets larger than Earth but smaller than Neptune have led to largely inconclusive interpretations at low spectral resolution due to hazes or clouds that obscure molecular features in their spectra. However, here we show that high-resolution spectroscopy (R $\sim$ 25,000 to 100,000) enables one to probe the regions in these atmospheres above the clouds where the cores… ▽ More Observations to characterize planets larger than Earth but smaller than Neptune have led to largely inconclusive interpretations at low spectral resolution due to hazes or clouds that obscure molecular features in their spectra. However, here we show that high-resolution spectroscopy (R $\sim$ 25,000 to 100,000) enables one to probe the regions in these atmospheres above the clouds where the cores of the strongest spectral lines are formed. We present models of transmission spectra for a suite of GJ1214b-like planets with thick photochemical hazes covering 1 - 5 $μ$m at a range of resolutions relevant to current and future ground-based spectrographs. Furthermore, we compare the utility of the cross-correlation function that is typically used with a more formal likelihood-based approach, finding that only the likelihood based method is sensitive to the presence of haze opacity. We calculate the signal-to-noise of these spectra, including telluric contamination, required to robustly detect a host of molecules such as CO, CO$_{2}$, H$_{2}$O, and CH$_{4}$, and photochemical products like HCN, as a function of wavelength range and spectral resolution. Spectra in M band require the lowest S/N$_{res}$ to detect multiple molecules simultaneously. CH$_{4}$ is only observable for the coolest models ($T_{\rm{eff}} =$ 412 K) and only in the L band. We quantitatively assess how these requirements compare to what is achievable with current and future instruments, demonstrating that characterization of small cool worlds with ground-based high resolution spectroscopy is well within reach. △ Less

Submitted 25 August, 2020; originally announced August 2020.

Comments: Submitted to AAS Journals, revised to reflect referee comments. Posting of this manuscript on the arXiv was coordinated with S. Ghandi et al

ExoMol line lists -- XXXIX. Ro-vibrational molecular line list for CO$_2$

Authors: S. N. Yurchenko, Thomas M. Mellor, Richard S. Freedman, J. Tennyson

Abstract: A new hot line list for the main isotopologue of CO$_2$ is presented. The line list consists of almost 2.5 billion transitions between 3.5 million rotation-vibration states of CO$_2$ in its ground electronic state, covering the wavenumber range 0-20000 cm$^{-1}$ ($λ>0.5$ $μ$m) with the upper and lower energy thresholds of 36 000 cm$^{-1}$ and 16 000 cm$^{-1}$, respectively. The ro-vibrational ener… ▽ More A new hot line list for the main isotopologue of CO$_2$ is presented. The line list consists of almost 2.5 billion transitions between 3.5 million rotation-vibration states of CO$_2$ in its ground electronic state, covering the wavenumber range 0-20000 cm$^{-1}$ ($λ>0.5$ $μ$m) with the upper and lower energy thresholds of 36 000 cm$^{-1}$ and 16 000 cm$^{-1}$, respectively. The ro-vibrational energies and wavefunctions are computed variationally using the Ames-2 accurate empirical potential energy surface. The ro-vibrational transition probabilities in the form of Einstein coefficients are computed using an accurate ab initio dipole moment surface using variational program TROVE. A new implementation of TROVE which uses an exact nuclear-motion kinetic energy operator is employed. Comparisons with the existing hot line lists are presented. The line list should be useful for atmospheric retrievals of exoplanets and cool stars. The UCL-4000 line list is available from the CDS and ExoMol databases. △ Less

Submitted 4 July, 2020; originally announced July 2020.

Observations of Disequilibrium CO Chemistry in the Coldest Brown Dwarfs

Authors: Brittany E. Miles, Andrew J. I. Skemer, Caroline V. Morley, Mark S. Marley, Jonathan J. Fortney, Katelyn N. Allers, Jacqueline K. Faherty, Thomas R. Geballe, Channon Visscher, Adam C. Schneider, Roxana Lupu, Richard S. Freedman, Gordon L. Bjoraker

Abstract: Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects discovered by future direct imaging surveys. We present a low resolution R $\sim$ 300 $M$-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 K and 250 K along with Jupiter. Thes… ▽ More Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects discovered by future direct imaging surveys. We present a low resolution R $\sim$ 300 $M$-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 K and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K$_{zz}$) to estimate the strength of vertical mixing in each object. The K$_{zz}$ values of cooler gaseous objects are close to their theoretical maximum and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K$_{zz}$ values imply that disequilibrium phosphine (PH$_{3}$) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH$_{3}$ absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855's $M$-band spectrum by including both CO and water clouds in the atmospheric model. △ Less

Submitted 27 April, 2020; v1 submitted 22 April, 2020; originally announced April 2020.

Comments: Accepted to the Astronomical Journal. 24 Pages, 16 Figures, 4 Tables. take care

The Need for Laboratory Measurements and Ab Initio Studies to Aid Understanding of Exoplanetary Atmospheres

Authors: Jonathan J. Fortney, Tyler D. Robinson, Shawn Domagal-Goldman, Anthony D. Del Genio, Iouli E. Gordon, Ehsan Gharib-Nezhad, Nikole Lewis, Clara Sousa-Silva, Vladimir Airapetian, Brian Drouin, Robert J. Hargreaves, Xinchuan Huang, Tijs Karman, Ramses M. Ramirez, Gregory B. Rieker, Jonathan Tennyson, Robin Wordsworth, Sergei N Yurchenko, Alexandria V Johnson, Timothy J. Lee, Chuanfei Dong, Stephen Kane, Mercedes Lopez-Morales, Thomas Fauchez, Timothy Lee , et al. (63 additional authors not shown)

Abstract: We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as c… ▽ More We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Comments: Submitted as an Astro2020 Science White Paper

Exploring A Photospheric Radius Correction to Model Secondary Eclipse Spectra for Transiting Exoplanets

Authors: Jonathan J. Fortney, Roxana E. Lupu, Caroline V. Morley, Richard S. Freedman, Callie Hood

Abstract: We highlight a physical effect that is often not considered that impacts the calculation of model spectra of planets at secondary eclipse, affecting both emission and reflection spectra. The radius of the emitting surface of the planet is not merely one value measured from a transit light curve, but is itself a function of wavelength, yet it is not directly measurable. At high precision, a similar… ▽ More We highlight a physical effect that is often not considered that impacts the calculation of model spectra of planets at secondary eclipse, affecting both emission and reflection spectra. The radius of the emitting surface of the planet is not merely one value measured from a transit light curve, but is itself a function of wavelength, yet it is not directly measurable. At high precision, a similar effect is well-known in transit "transmission spectroscopy" but this related effect also impacts emission and reflection. As is well-appreciated, the photospheric radius can vary across $\sim$4-8 atmospheric scale heights, depending on atmospheric opacity and spectral resolution. This effect leads to a decreased weighting in model calculations at wavelengths where atmospheric opacity is low, and one sees more deeply into the atmosphere, to a smaller radius. The overall effect serves to mute emission spectra features for atmospheres with no thermal inversion but to enhance features for atmospheres with a thermal inversion. While this effect can be ignored for current \emph{Hubble} observations, it can lead to wavelength-dependent 10-20\% changes in planet-to-star flux ratios in the infrared at $R\sim~200-1000$ (readily achievable for JWST) for low-gravity hot Jupiters, although values of 5\% are more typical for the population. The effect is mostly controlled by the ratio of the atmospheric scale height to the planet radius, and can be important at any planetary temperature. Of known planets, the effect is largest for the cool "super-puffs" at very low surface gravity, where it can alter calculated flux ratios by over 100\%. We discuss complexities of including this photospheric radius effect in 1D and 3D atmosphere models. △ Less

Submitted 29 March, 2019; originally announced April 2019.

Comments: Submitted to ApJ Letters, revised to reflect referee comments

Imaging Cool Giant Planets in Reflected Light: Science Investigations and Synergy with Habitable Planets

Authors: Mark Marley, Nikole Lewis, Giada Arney, Vanessa Bailey, Natasha Batalha, Charles Beichman, Björn Benneke, Jasmina Blecic, Kerri Cahoy, Jeffrey Chilcote, Shawn Domagal-Goldman, Courtney Dressing, Michael Fitzgerald, Jonathan Fortney, Richard Freedman, Dawn Gelino, John Gizis, Olivier Guyon, Thomas Greene, Heidi Hammel, Yasuhiro Hasegawa, Nemanja Jovanovic, Quinn Konopacky, Ravi Kopparapu, Michael Liu , et al. (16 additional authors not shown)

Abstract: Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the… ▽ More Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the decades to follow. △ Less

Submitted 26 March, 2019; v1 submitted 21 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro 2020 Decadal Survey on Astronomy and Astrophysics. Replace version to fix typo in co-signer name and add figure credits

Unusual Isotopic Abundances in a Fully-Convective Stellar Binary

Authors: I. J. M. Crossfield, J. D. Lothringer, B. Flores, E. A. C. Mills, R. Freedman, J. Valverde, B. Miles, X. Guo, A. Skemer

Abstract: Low-mass M dwarfs represent the most common outcome of star formation, but their complex emergent spectra hinder detailed studies of their composition and initial formation. The measurement of isotopic ratios is a key tool that has been used to unlock the formation of our Solar System, the Sun, and the nuclear processes within more massive stars. We observed GJ 745AB, two M dwarfs orbiting in a wi… ▽ More Low-mass M dwarfs represent the most common outcome of star formation, but their complex emergent spectra hinder detailed studies of their composition and initial formation. The measurement of isotopic ratios is a key tool that has been used to unlock the formation of our Solar System, the Sun, and the nuclear processes within more massive stars. We observed GJ 745AB, two M dwarfs orbiting in a wide binary, with the IRTF/iSHELL spectrograph. Our spectroscopy of CO in these stars at the 4.7 micron fundamental and 2.3 micron first-overtone rovibrational bandheads reveals 12C16O, 13C16O, and 12C18O in their photospheres. Since the stars are fully convective, the atomic constituents of these isotopologues should be uniformly mixed throughout the stars' interiors. We find that in these M dwarfs, both 12C/13C and 16O/18O greatly exceed the Solar values. These measurements cannot be explained solely by models of Galactic chemical evolution, but require that the stars formed from an ISM significantly enriched by material ejected from an exploding core-collape supernova. These isotopic measurements complement the elemental abundances provided by large-scale spectroscopic surveys, and open a new window onto studies of Galactic evolution, stellar populations, and individual systems. △ Less

Submitted 9 January, 2019; originally announced January 2019.

Comments: 11 pages, 4 data files, 3 figures, 2 tables. ApJ in press

Journal ref: ApJL 2019, 817, 1, L3

Perspectives on Astrophysics Based on Atomic, Molecular, and Optical (AMO) Techniques

Authors: Daniel Wolf Savin, James F. Babb, Paul M. Bellan, Crystal Brogan, Jan Cami, Paola Caselli, Lia Corrales, Gerardo Dominguez, Steven R. Federman, Chris J. Fontes, Richard Freedman, Brad Gibson, Leon Golub, Thomas W. Gorczyca, Michael Hahn, Sarah M. Hörst, Reggie L. Hudson, Jeffrey Kuhn, James E. Lawler, Maurice A. Leutenegger, Joan P. Marler, Michael C. McCarthy, Brett A. McGuire, Stefanie N. Milam, Nicholas A. Murphy , et al. (13 additional authors not shown)

Abstract: About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides i… ▽ More About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides in this intermediate region. We put forward that the next frontier for AMO science is to explore the AMO complexity that describes most of the Cosmos. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Comments: White paper submission to the Decadal Assessment and Outlook Report on Atomic, Molecular, and Optical (AMO) Science (AMO 2020)

Measuring the D/H Ratios of Exoplanets and Brown Dwarfs

Authors: Caroline V. Morley, Andrew J. Skemer, Brittany E. Miles, Michael R. Line, Eric D. Lopez, Matteo Brogi, Richard S. Freedman, Mark S. Marley

Abstract: The relative abundance of deuterium and hydrogen is a potent tracer of planet formation and evolution. Jupiter and Saturn have protosolar atmospheric D/H ratios, a relic of substantial gas accretion from the nebula, while the atmospheres of Neptune and Uranus are enhanced in D by accretion of ices into their envelopes. For terrestrial planets, D/H ratios are used to determine the mechanisms of vol… ▽ More The relative abundance of deuterium and hydrogen is a potent tracer of planet formation and evolution. Jupiter and Saturn have protosolar atmospheric D/H ratios, a relic of substantial gas accretion from the nebula, while the atmospheres of Neptune and Uranus are enhanced in D by accretion of ices into their envelopes. For terrestrial planets, D/H ratios are used to determine the mechanisms of volatile delivery and subsequent atmosphere loss over the lifetime of the planet. Planets and brown dwarfs more massive than 13 M_J quickly fuse their initial D reservoir. Here, we simulate spectra for giant exoplanets and brown dwarfs (2 M_Neptune to ~10 M_Jupiter) from Teff=200--1800 K including both CH3D and HDO to determine the observability of these dominant deuterium isotopologues in mid-infrared thermal emission spectra. Colder objects have stronger molecular features in their spectra, due to the temperature-dependence of molecular cross sections. CH3D is easier to observe than HDO at all temperatures considered, due to the strength of its absorption bands and locations of features at wavelengths with few other strong absorption features. We predict that for nearby cool brown dwarfs, the CH3D abundance can be measured at high signal-to-noise with the James Webb Space Telescope; for objects from 200--800 K closer than 10 pc, a protosolar D/H ratio would be readily observable in 2.5 hours. Moderately young Jupiter-mass planets (100--300 Myr) and young Neptunes (10 Myr) may be discovered with JWST and provide the best targets for detecting deuterium on an exoplanet in the coming decade. Future telescope designs should consider the importance of isotopes for understanding the formation and evolution of planetary atmospheres. △ Less

Submitted 11 October, 2019; v1 submitted 9 October, 2018; originally announced October 2018.

Comments: 10 pages, 5 figures, published in ApJ Letters

Journal ref: The Astrophysical Journal Letters, Volume 882, Issue 2, article id. L29, 8 pp. (2019)

An L Band Spectrum of the Coldest Brown Dwarf

Authors: Caroline V. Morley, Andrew J. Skemer, Katelyn N. Allers, Mark. S. Marley, Jacqueline K. Faherty, Channon Visscher, Samuel A. Beiler, Brittany E. Miles, Roxana Lupu, Richard S. Freedman, Jonathan J. Fortney, Thomas R. Geballe, Gordon L. Bjoraker

Abstract: The coldest brown dwarf, WISE 0855, is the closest known planetary-mass, free-floating object and has a temperature nearly as cold as the solar system gas giants. Like Jupiter, it is predicted to have an atmosphere rich in methane, water, and ammonia, with clouds of volatile ices. WISE 0855 is faint at near-infrared wavelengths and emits almost all its energy in the mid-infrared. Skemer et al. 201… ▽ More The coldest brown dwarf, WISE 0855, is the closest known planetary-mass, free-floating object and has a temperature nearly as cold as the solar system gas giants. Like Jupiter, it is predicted to have an atmosphere rich in methane, water, and ammonia, with clouds of volatile ices. WISE 0855 is faint at near-infrared wavelengths and emits almost all its energy in the mid-infrared. Skemer et al. 2016 presented a spectrum of WISE 0855 from 4.5-5.1 micron (M band), revealing water vapor features. Here, we present a spectrum of WISE 0855 in L band, from 3.4-4.14 micron. We present a set of atmosphere models that include a range of compositions (metallicities and C/O ratios) and water ice clouds. Methane absorption is clearly present in the spectrum. The mid-infrared color can be better matched with a methane abundance that is depleted relative to solar abundance. We find that there is evidence for water ice clouds in the M band spectrum, and we find a lack of phosphine spectral features in both the L and M band spectra. We suggest that a deep continuum opacity source may be obscuring the near-infrared flux, possibly a deep phosphorous-bearing cloud, ammonium dihyrogen phosphate. Observations of WISE 0855 provide critical constraints for cold planetary atmospheres, bridging the temperature range between the long-studied solar system planets and accessible exoplanets. JWST will soon revolutionize our understanding of cold brown dwarfs with high-precision spectroscopy across the infrared, allowing us to study their compositions and cloud properties, and to infer their atmospheric dynamics and formation processes. △ Less

Submitted 20 April, 2018; originally announced April 2018.

Comments: 19 pages, 21 figures. Accepted for publication in ApJ

Retrieval of atmospheric properties of cloudy L dwarfs

Authors: Ben Burningham, Mark S. Marley, Michael R. Line, Roxana Lupu, Channon Visscher, Caroline V. Morley, Didier Saumon, Richard Freedman

Abstract: We present the first results from applying the spectral inversion technique in the cloudy L dwarf regime. Our new framework provides a flexible approach to modelling cloud opacity which can be built incrementally as the data requires, and improves upon previous retrieval experiments in the brown dwarf regime by allowing for scattering in two stream radiative transfer. Our first application of the… ▽ More We present the first results from applying the spectral inversion technique in the cloudy L dwarf regime. Our new framework provides a flexible approach to modelling cloud opacity which can be built incrementally as the data requires, and improves upon previous retrieval experiments in the brown dwarf regime by allowing for scattering in two stream radiative transfer. Our first application of the tool to two mid-L dwarfs is able to reproduce their near-infrared spectra far more closely than grid models. Our retrieved thermal, chemical, and cloud profiles allow us to estimate $T_{\rm eff} = 1796^{+23}_{-25}$ K and $\log g = 5.21^{+0.05}_{-0.08}$ for 2MASS J05002100+0330501 and for 2MASSW J2224438-015852 we find $T_{\rm eff} = 1723^{+18}_{-19}$ K and $\log g = 5.31^{+0.04}_{-0.08}$, in close agreement with previous empirical estimates. Our best model for both objects includes an optically thick cloud deck which passes $τ_{cloud} \geq 1$ (looking down) at a pressure of around 5 bar. The temperature at this pressure is too high for silicate species to condense, and we argue that corundum and/or iron clouds are responsible for this cloud opacity. Our retrieved profiles are cooler at depth, and warmer at altitude than the forward grid models that we compare, and we argue that some form of heating mechanism may be at work in the upper atmospheres of these L dwarfs. We also identify anomalously high CO abundance in both targets, which does not correlate with the warmth of our upper atmospheres or our choice of cloud model, and find similarly anomalous alkali abundance for one of our targets. These anomalies may reflect unrecognised shortcomings in our retrieval model, or inaccuracies in our gas phase opacities. △ Less

Submitted 18 May, 2017; v1 submitted 5 January, 2017; originally announced January 2017.

Comments: Accepted for publication in MNRAS. Changes in review: additional validation against T dwarf test case, posteriors for cloud free retrievals and contribution functions added to Appendix. Also, entire investigation re-run with new UCL H2O opacities, with minimal impact on results

Uniform Atmospheric Retrieval Analysis of Ultracool Dwarfs II: Properties of 11 T-dwarfs

Authors: Michael R. Line, Mark S. Marley, Michael C. Liu, Ben Burningham, Caroline V. Morley, Natalie R. Hinkel, Johanna Teske, Jonathan J. Fortney, Roxana Lupu, Richard Freedman

Abstract: Brown dwarf spectra are rich in information revealing of the chemical and physical processes operating in their atmospheres. We apply a recently developed atmospheric retrieval tool to an ensemble of late T-dwarf (600-800K) near infrared spectra. With these spectra we are able to place direct constraints the molecular abundances of H$_2$O, CH$_4$, CO, CO$_2$, NH$_3$, H$_2$S, and Na+K, gravity, the… ▽ More Brown dwarf spectra are rich in information revealing of the chemical and physical processes operating in their atmospheres. We apply a recently developed atmospheric retrieval tool to an ensemble of late T-dwarf (600-800K) near infrared spectra. With these spectra we are able to place direct constraints the molecular abundances of H$_2$O, CH$_4$, CO, CO$_2$, NH$_3$, H$_2$S, and Na+K, gravity, thermal structure (and effective temperature), photometric radius, and cloud optical depths. We find that ammonia, water, methane, and the alkali metals are present and well constrained in all 11 objects. From the abundance constraints we find no significant trend in the water, methane, or ammonia abundances with temperature, but find a very strong ($>$25$σ$) increasing trend in the alkali metal abundances with effective temperature, indicative of alkali rainout. We also find little evidence for optically thick clouds. With the methane and water abundances, we derive the intrinsic atmospheric metallicity and carbon-to-oxygen ratios. We find in our sample, that metallicities are typically sub solar and carbon-to-oxygen ratios are somewhat super solar, different than expectations from the local stellar population. We also find that the retrieved vertical thermal profiles are consistent with radiative equilibrium over the photospheric regions. Finally, we find that our retrieved effective temperatures are lower than previous inferences for some objects and that our radii are larger than expectations from evolutionary models, possibly indicative of un-resolved binaries. This investigation and methodology represents a paradigm in linking spectra to the determination of the fundamental chemical and physical processes governing cool brown dwarf atmospheres. △ Less

Submitted 24 July, 2017; v1 submitted 8 December, 2016; originally announced December 2016.

Comments: Accepted ApJ. Supplemental material including full posteriors will be included through the link in the published ApJ article

The Need for Laboratory Work to Aid in The Understanding of Exoplanetary Atmospheres

Authors: Jonathan J. Fortney, Tyler D. Robinson, Shawn Domagal-Goldman, David Skålid Amundsen, Matteo Brogi, Mark Claire, David Crisp, Eric Hebrard, Hiroshi Imanaka, Remco de Kok, Mark S. Marley, Dillon Teal, Travis Barman, Peter Bernath, Adam Burrows, David Charbonneau, Richard S. Freedman, Dawn Gelino, Christiane Helling, Kevin Heng, Adam G. Jensen, Stephen Kane, Eliza M. -R. Kempton, Ravi Kumar Kopparapu, Nikole K. Lewis , et al. (9 additional authors not shown)

Abstract: Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements… ▽ More Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr) △ Less

Submitted 23 February, 2016; v1 submitted 19 February, 2016; originally announced February 2016.

Comments: This is a white paper originally initiated within NASA's Nexus for Exoplanet System Science (NExSS), with additional community input. Community comment is welcome. Please submit comments and suggestions by going to the Google Doc version of the paper: https://docs.google.com/document/d/1109uy_5n9VWuOOn-huTVv-gJgBFRO9QH0Ojkao1a6Lk/edit?usp=sharing

The LEECH Exoplanet Imaging Survey: Characterization of the Coldest Directly Imaged Exoplanet, GJ 504 b, and Evidence for Super-Stellar Metallicity

Authors: Andrew J. Skemer, Caroline V. Morley, Neil T. Zimmerman, Michael F. Skrutskie, Jarron Leisenring, Esther Buenzli, Mickael Bonnefoy, Vanessa Bailey, Philip Hinz, Denis Defrére, Simone Esposito, Dániel Apai, Beth Biller, Wolfgang Brandner, Laird Close, Justin R. Crepp, Robert J. De Rosa, Silvano Desidera, Josh Eisner, Jonathan Fortney, Richard Freedman, Thomas Henning, Karl-Heinz Hofmann, Taisiya Kopytova, Roxana Lupu , et al. (17 additional authors not shown)

Abstract: As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the first generation of directly-imaged exoplanets were all L-type. Recently, Kuzuhara et al. (2… ▽ More As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the first generation of directly-imaged exoplanets were all L-type. Recently, Kuzuhara et al. (2013) announced the discovery of GJ 504 b, the first T dwarf exoplanet. GJ 504 b provides a unique opportunity to study the atmosphere of a new type of exoplanet with a ~500 K temperature that bridges the gap between the first directly imaged planets (~1000 K) and our own Solar System's Jupiter (~130 K). We observed GJ 504 b in three narrow L-band filters (3.71, 3.88, and 4.00 microns), spanning the red end of the broad methane fundamental absorption feature (3.3 microns) as part of the LEECH exoplanet imaging survey. By comparing our new photometry and literature photometry to a grid of custom model atmospheres, we were able to fit GJ 504 b's unusual spectral energy distribution for the first time. We find that GJ 504 b is well-fit by models with the following parameters: T_eff=544+/-10 K, g<600 m/s^2, [M/H]=0.60+/-0.12, cloud opacity parameter of f_sed=2-5, R=0.96+/-0.07 R_Jup, and log(L)=-6.13+/-0.03 L_Sun, implying a hot start mass of 3-30 M_jup for a conservative age range of 0.1-6.5 Gyr. Of particular interest, our model fits suggest that GJ 504 b has a super-stellar metallicity. Since planet formation can create objects with non-stellar metallicities, while binary star formation cannot, this result suggests that GJ 504 b formed like a planet, not like a binary companion. △ Less

Submitted 30 November, 2015; originally announced November 2015.

Comments: resubmitted to ApJ with referee's comments

Gaseous Mean Opacities for Giant Planet and Ultracool Dwarf Atmospheres over a Range of Metallicities and Temperatures

Authors: Richard S. Freedman, Jacob Lustig-Yaeger, Jonathan J. Fortney, Roxana E. Lupu, Mark S. Marley, Katharina Lodders

Abstract: 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 th… ▽ More 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. △ Less

Submitted 24 June, 2021; v1 submitted 29 August, 2014; originally announced September 2014.

Comments: This version fixes a units description error that was only in the arxiv version, but is correct in the ApJS version. In equations 4 and 5 pressure is in dyne/cm2. Online tables available at the ApJS website: https://iopscience.iop.org/article/10.1088/0067-0049/214/2/25

Warming early Mars with CO2 and H2

Authors: Ramses M. Ramirez, Ravi Kopparapu, Michael E. Zugger, Tyler D. Robinson, Richard Freedman, James F. Kasting

Abstract: The presence of valleys on ancient terrains of Mars suggest that liquid water flowed on the martian surface 3.8 billion years ago or before. The above-freezing temperatures required to explain valley formation could have been transient, in response to frequent large meteorite impacts on early Mars, or they could have been caused by long-lived greenhouse warming. Climate models that consider only t… ▽ More The presence of valleys on ancient terrains of Mars suggest that liquid water flowed on the martian surface 3.8 billion years ago or before. The above-freezing temperatures required to explain valley formation could have been transient, in response to frequent large meteorite impacts on early Mars, or they could have been caused by long-lived greenhouse warming. Climate models that consider only the greenhouse gases carbon dioxide and water vapor have been unable to recreate warm surface conditions, given the lower solar luminosity at that time. Here we use a one-dimensional climate model to demonstrate that an atmosphere containing 1.3-4 bar of CO2 and water vapor, along with 5 to 20 percent H2, could have raised the mean surface temperature of early Mars above the freezing point of water. Vigorous volcanic outgassing from a highly reduced early martian mantle is expected to provide sufficient atmospheric H2 and CO2, the latter from the photochemical oxidation of outgassed CH4 and CO, to form a CO2-H2 greenhouse. Such a dense early martian atmosphere is consistent with independent estimates of surface pressure based on cratering data. △ Less

Submitted 26 May, 2014; originally announced May 2014.

Comments: 48 pages, 12 figures, and 4 tables. The official Nature Geoscience version is found here: http://www.nature.com/ngeo/journal/v7/n1/full/ngeo2000.html

Journal ref: Nature Geoscience 7, 59 to 63 (2014)

The atmospheric circulation of the super Earth GJ 1214b: Dependence on composition and metallicity

Authors: Tiffany Kataria, Adam P. Showman, Jonathan J. Fortney, Mark S. Marley, Richard S. Freedman

Abstract: We present three-dimensional atmospheric circulation models of GJ 1214b, a 2.7 Earth-radius, 6.5 Earth-mass super Earth detected by the MEarth survey. Here we explore the planet's circulation as a function of atmospheric metallicity and atmospheric composition, modeling atmospheres with a low mean-molecular weight (i.e., H2-dominated) and a high mean-molecular weight (i.e. water- and CO2-dominated… ▽ More We present three-dimensional atmospheric circulation models of GJ 1214b, a 2.7 Earth-radius, 6.5 Earth-mass super Earth detected by the MEarth survey. Here we explore the planet's circulation as a function of atmospheric metallicity and atmospheric composition, modeling atmospheres with a low mean-molecular weight (i.e., H2-dominated) and a high mean-molecular weight (i.e. water- and CO2-dominated). We find that atmospheres with a low mean-molecular weight have strong day-night temperature variations at pressures above the infrared photosphere that lead to equatorial superrotation. For these atmospheres, the enhancement of atmospheric opacities with increasing metallicity lead to shallower atmospheric heating, larger day-night temperature variations and hence stronger superrotation. In comparison, atmospheres with a high mean-molecular weight have larger day-night and equator-to-pole temperature variations than low mean-molecular weight atmospheres, but differences in opacity structure and energy budget lead to differences in jet structure. The circulation of a water-dominated atmosphere is dominated by equatorial superrotation, while the circulation of a CO2-dominated atmosphere is instead dominated by high-latitude jets. By comparing emergent flux spectra and lightcurves for 50x solar and water-dominated compositions, we show that observations in emission can break the degeneracy in determining the atmospheric composition of GJ 1214b. The variation in opacity with wavelength for the water-dominated atmosphere leads to large phase variations within water bands and small phase variations outside of water bands. The 50x solar atmosphere, however, yields small variations within water bands and large phase variations at other characteristic wavelengths. These observations would be much less sensitive to clouds, condensates, and hazes than transit observations. △ Less

Submitted 4 March, 2014; v1 submitted 9 January, 2014; originally announced January 2014.

Comments: 12 pages, 11 figures, 2 tables, accepted to ApJ

The Atmospheres of Earth-like Planets after Giant Impact Events

Authors: R. E. Lupu, Kevin Zahnle, Mark S. Marley, Laura Schaefer, Bruce Fegley, Caroline Morley, Kerri Cahoy, Richard Freedman, Jonathan J. Fortney

Abstract: It is now understood that the accretion of terrestrial planets naturally involves giant collisions, the moon-forming impact being a well known example. In the aftermath of such collisions the surface of the surviving planet is very hot and potentially detectable. Here we explore the atmospheric chemistry, photochemistry, and spectral signatures of post-giant-impact terrestrial planets enveloped by… ▽ More It is now understood that the accretion of terrestrial planets naturally involves giant collisions, the moon-forming impact being a well known example. In the aftermath of such collisions the surface of the surviving planet is very hot and potentially detectable. Here we explore the atmospheric chemistry, photochemistry, and spectral signatures of post-giant-impact terrestrial planets enveloped by thick atmospheres consisting predominantly of CO2, and H2O. The atmospheric chemistry and structure are computed self-consistently for atmospheres in equilibrium with hot surfaces with composition reflecting either the bulk silicate Earth (which includes the crust, mantle, atmosphere and oceans) or Earth's continental crust. We account for all major molecular and atomic opacity sources including collision-induced absorption. We find that these atmospheres are dominated by H2O and CO2, while the formation of CH4, and NH3 is quenched due to short dynamical timescales. Other important constituents are HF, HCl, NaCl, and SO2. These are apparent in the emerging spectra, and can be indicative that an impact has occurred. The use of comprehensive opacities results in spectra that are a factor of 2 lower in surface brightness in the spectral windows than predicted by previous models. The estimated luminosities show that the hottest post-giant-impact planets will be detectable with near-infrared coronagraphs on the planned 30m-class telescopes. The 1-4um region will be most favorable for such detections, offering bright features and better contrast between the planet and a potential debris disk. We derive cooling timescales on the order of 10^5-10^6 Myrs, based on the modeled effective temperatures. This leads to the possibility of discovering tens of such planets in future surveys. △ Less

Submitted 6 January, 2014; originally announced January 2014.

Comments: 51 pages, 16 figures, preprint format; ApJ submitted

Bulk Composition of GJ 1214b and other sub-Neptune exoplanets

Authors: Diana Valencia, Tristan Guillot, Vivien Parmentier, Richard S. Freedman

Abstract: GJ1214b stands out among the detected low-mass exoplanets, because it is, so far, the only one amenable to transmission spectroscopy. Up to date there is no consensus about the composition of its envelope although most studies suggest a high molecular weight atmosphere. In particular, it is unclear if hydrogen and helium are present or if the atmosphere is water dominated. Here, we present results… ▽ More GJ1214b stands out among the detected low-mass exoplanets, because it is, so far, the only one amenable to transmission spectroscopy. Up to date there is no consensus about the composition of its envelope although most studies suggest a high molecular weight atmosphere. In particular, it is unclear if hydrogen and helium are present or if the atmosphere is water dominated. Here, we present results on the composition of the envelope obtained by using an internal structure and evolutionary model to fit the mass and radius data. By examining all possible mixtures of water and H/He, with the corresponding opacities, we find that the bulk amount of H/He of GJ1214b is at most 7% by mass. In general, we find the radius of warm sub-Neptunes to be most sensitive to the amount of H/He. We note that all (Kepler-11b,c,d,f, Kepler-18b, Kepler-20b, 55Cnc-e, Kepler-36c and Kepler-68b) but two (Kepler-11e and Kepler-30b) of the discovered low-mass planets so far have less than 10% H/He. In fact, Kepler-11e and Kepler-30b have 10-18% and 5-15% bulk H/He. Conversely, little can be determined about the H2O or rocky content of sub-Neptune planets. We find that although a 100% water composition fits the data for GJ1214b, based on formation constraints the presence of heavier refractory material on this planet is expected, and hence, so is a component lighter than water required. A robust determination by transmission spectroscopy of the composition of the upper atmosphere of GJ1214b will help determine the extent of compositional segregation between the atmosphere and envelope. △ Less

Submitted 14 June, 2013; v1 submitted 12 May, 2013; originally announced May 2013.

Comments: Updated the masses and radii of the Kepler-11 system, added Kepler-30b as well in the analysis. Accepted in ApJ, 39 pages, 9 figures

Chemical consequences of the C/O ratio on hot Jupiters: Examples from WASP-12b,CoRoT-2b, XO-1b, and HD 189733b

Authors: J. I. Moses, N. Madhusudhan, C. Visscher, R. S. Freedman

Abstract: Motivated by recent spectroscopic evidence for carbon-rich atmospheres on some transiting exoplanets, we investigate the influence of the C/O ratio on the chemistry, composition, and spectra of extrasolar giant planets both from a thermochemical-equilibrium perspective and from consideration of disequilibrium processes like photochemistry and transport-induced quenching. We find that although CO i… ▽ More Motivated by recent spectroscopic evidence for carbon-rich atmospheres on some transiting exoplanets, we investigate the influence of the C/O ratio on the chemistry, composition, and spectra of extrasolar giant planets both from a thermochemical-equilibrium perspective and from consideration of disequilibrium processes like photochemistry and transport-induced quenching. We find that although CO is predicted to be a major atmospheric constituent on hot Jupiters for all C/O ratios, other oxygen-bearing molecules like H2O and CO2 are much more abundant when C/O < 1, whereas CH4, HCN, and C2H2 gain significantly in abundance when C/O > 1. Disequilibrium processes tend to enhance the abundance of CH4, NH3, HCN, and C2H2 over a wide range of C/O ratios. We compare the results of our models with secondary-eclipse photometric data from the Spitzer Space Telescope and conclude that (1) disequilibrium models with C/O ~ 1 are consistent with spectra of WASP-12b, XO-1b, and CoRoT-2b, confirming the possible carbon-rich nature of these planets, (2) spectra from HD 189733b are consistent with C/O ~< 1, but as the assumed metallicity is increased above solar, the required C/O ratio must increase toward 1 to prevent too much H2O absorption, (3) species like HCN can have a significant influence on spectral behavior in the 3.6 and 8.0 um Spitzer channels, potentially providing even more opacity than CH4 when C/O > 1, and (4) the very high CO2 abundance inferred for HD 189733b from near-infrared observations cannot be explained through equilibrium or disequilibrium chemistry in a H2-dominated atmosphere. We discuss possible formation mechanisms for carbon-rich hot Jupiters. The C/O ratio and bulk atmospheric metallicity provide important clues regarding the formation and evolution of the giant planets. △ Less

Submitted 13 November, 2012; originally announced November 2012.

Comments: accepted in the Astrophysical Journal

Journal ref: Astrophys. J., 763, 25 (2013)

Three-dimensional atmospheric circulation of hot Jupiters on highly eccentric orbits

Authors: Tiffany Kataria, Adam P. Showman, Nikole K. Lewis, Jonathan J. Fortney, Mark S. Marley, Richard S. Freedman

Abstract: Of the over 800 exoplanets detected to date, over half are on non-circular orbits, with eccentricities as high as 0.93. Such orbits lead to time-variable stellar heating, which has implications for the planet's atmospheric dynamical regime. However, little is known about this dynamical regime, and how it may influence observations. Therefore, we present a systematic study of hot Jupiters on highly… ▽ More Of the over 800 exoplanets detected to date, over half are on non-circular orbits, with eccentricities as high as 0.93. Such orbits lead to time-variable stellar heating, which has implications for the planet's atmospheric dynamical regime. However, little is known about this dynamical regime, and how it may influence observations. Therefore, we present a systematic study of hot Jupiters on highly eccentric orbits using the SPARC/MITgcm, a model which couples a three-dimensional general circulation model with a plane-parallel, two-stream, non-grey radiative transfer model. In our study, we vary the eccentricity and orbit-average stellar flux over a wide range. We demonstrate that the eccentric hot Jupiter regime is qualitatively similar to that of planets on circular orbits; the planets possess a superrotating equatorial jet and exhibit large day-night temperature variations. We show that these day-night heating variations induce momentum fluxes equatorward to maintain the superrotating jet throughout its orbit. As the eccentricity and/or stellar flux is increased, the superrotating jet strengthens and narrows, due to a smaller Rossby deformation radius. For a select number of model integrations, we generate full-orbit lightcurves and find that the timing of transit and secondary eclipse viewed from Earth with respect to periapse and apoapse can greatly affect what we see in infrared (IR) lightcurves; the peak in IR flux can lead or lag secondary eclipse depending on the geometry. For those planets that have large day-night temperature variations and rapid rotation rates, we find that the lightcurves exhibit "ringing" as the planet's hottest region rotates in and out of view from Earth. These results can be used to explain future observations of eccentric transiting exoplanets. △ Less

Submitted 6 March, 2013; v1 submitted 18 August, 2012; originally announced August 2012.

Comments: 20 pages, 18 figures, 2 tables; Accepted to ApJ

Masses, Radii, and Cloud Properties of the HR 8799 Planets

Authors: Mark S. Marley, Didier Saumon, Michael Cushing, Andrew S. Ackerman, Jonathan J. Fortney, Richard Freedman

Abstract: The near-infrared colors of the planets directly imaged around the A star HR 8799 are much redder than most field brown dwarfs of the same effective temperature. Previous theoretical studies of these objects have concluded that the atmospheres of planets b, c, and d are unusually cloudy or have unusual cloud properties. Some studies have also found that the inferred radii of some or all of the pla… ▽ More The near-infrared colors of the planets directly imaged around the A star HR 8799 are much redder than most field brown dwarfs of the same effective temperature. Previous theoretical studies of these objects have concluded that the atmospheres of planets b, c, and d are unusually cloudy or have unusual cloud properties. Some studies have also found that the inferred radii of some or all of the planets disagree with expectations of standard giant planet evolution models. Here we compare the available data to the predictions of our own set of atmospheric and evolution models that have been extensively tested against observations of field L and T dwarfs, including the reddest L dwarfs. Unlike some previous studies we require mutually consistent choices for effective temperature, gravity, cloud properties, and planetary radius. This procedure thus yields plausible values for the masses, effective temperatures, and cloud properties of all three planets. We find that the cloud properties of the HR 8799 planets are not unusual but rather follow previously recognized trends, including a gravity dependence on the temperature of the L to T spectral transition--some reasons for which we discuss. We find the inferred mass of planet b is highly sensitive to whether or not we include the H and K band spectrum in our analysis. Solutions for planets c and d are consistent with the generally accepted constraints on the age of the primary star and orbital dynamics. We also confirm that, like in L and T dwarfs and solar system giant planets, non-equilibrium chemistry driven by atmospheric mixing is also important for these objects. Given the preponderance of data suggesting that the L to T spectral type transition is gravity dependent, we present an exploratory evolution calculation that accounts for this effect. Finally we recompute the the bolometric luminosity of all three planets. △ Less

Submitted 18 June, 2012; v1 submitted 29 May, 2012; originally announced May 2012.

Comments: 52 pages, 12 figures, Astrophysical Journal, in press. v2 features minor editorial updates and corrections

New H2 collision-induced absorption and NH3 opacity and the spectra of the coolest brown dwarfs

Authors: D. Saumon, M. S. Marley, M. Abel, L. Frommhold, R. S. Freedman

Abstract: We present new cloudy and cloudless model atmospheres for brown dwarfs using recent ab initio calculations of the line list of ammonia (NH3) and of the collision-induced absorption of molecular hydrogen (H2). We compare the new synthetic spectra with models based on an earlier description of the H2 and NH3 opacities. We find a significant improvement in fitting the nearly complete spectral energy… ▽ More We present new cloudy and cloudless model atmospheres for brown dwarfs using recent ab initio calculations of the line list of ammonia (NH3) and of the collision-induced absorption of molecular hydrogen (H2). We compare the new synthetic spectra with models based on an earlier description of the H2 and NH3 opacities. We find a significant improvement in fitting the nearly complete spectral energy distribution of the T7p dwarf Gliese 570D and in near infrared color-magnitude diagrams of field brown dwarfs. We apply these new models to the identification of NH3 absorption in the H band peak of very late T dwarfs and the new Y dwarfs and discuss the observed trend in the NH3-H spectral index. The new NH3 line list also allows a detailed study of the medium resolution spectrum of the T9/T10 dwarf UGPS J072227.51-054031.2 where we identify several specific features caused by NH3. △ Less

Submitted 28 February, 2012; originally announced February 2012.

Comments: 37 pages, 13 figures. Accepted for publication in the Astrophysical J

Report number: LA-UR-11-07065

The Discovery of Y Dwarfs Using Data from the Wide-field Infrared Survey Explorer (WISE)

Authors: Michael C. Cushing, J. Davy Kirkpatrick, Christopher R. Gelino, Roger L. Griffith, Michael F. Skrutskie, Amanda K. Mainzer, Kenneth A. Marsh, Charles A. Beichman, Adam J. Burgasser, Lisa A. Prato, Robert A. Simcoe, Mark S. Marley, D. Saumon, Richard S. Freedman, Peter R. Eisenhardt, Edward L. Wright

Abstract: We present the discovery of seven ultracool brown dwarfs identified with the Wide-field Infrared Survey Explorer (WISE). Near-infrared spectroscopy reveals deep absorption bands of H_2O and CH_4 that indicate all seven of the brown dwarfs have spectral types later than UGPS J072227.51-054031.2, the latest type T dwarf currently known. The spectrum of WISEP J182831.08+265037.8 is distinct in that t… ▽ More We present the discovery of seven ultracool brown dwarfs identified with the Wide-field Infrared Survey Explorer (WISE). Near-infrared spectroscopy reveals deep absorption bands of H_2O and CH_4 that indicate all seven of the brown dwarfs have spectral types later than UGPS J072227.51-054031.2, the latest type T dwarf currently known. The spectrum of WISEP J182831.08+265037.8 is distinct in that the heights of the J- and H-band peaks are approximately equal in units of f_lambda, so we identify it as the archetypal member of the Y spectral class. The spectra of at least two of the other brown dwarfs exhibit absorption on the blue wing of the H-band peak that we tentatively ascribe to NH_3. These spectral morphological changes provide a clear transition between the T dwarfs and the Y dwarfs. In order to produce a smooth near-infrared spectral sequence across the T/Y dwarf transition, we have reclassified UGPS J0722-0540 as the T9 spectral standard and tentatively assign WISEP J173835.52+273258.9 as the Y0 spectral standard. In total, six of the seven new brown dwarfs are classified as Y dwarfs: four are classified as Y0, one is classified as Y0 (pec?), and WISEP J1828+2650 is classified as >Y0. We have also compared the spectra to the model atmospheres of Marley and Saumon and infer that the brown dwarfs have effective temperatures ranging from 300 K to 500 K, making them the coldest spectroscopically confirmed brown dwarfs known to date. △ Less

Submitted 23 August, 2011; originally announced August 2011.

Comments: Submitted June 11 and accepted August 2 for publication in the Astrophysical Journal

Disequilibrium Carbon, Oxygen, and Nitrogen Chemistry in the Atmospheres of HD 189733b and HD 209458b

Authors: Julianne I. Moses, Channon Visscher, Jonathan J. Fortney, Adam P. Showman, Nikole K. Lewis, Caitlin A. Griffith, Stephen J. Klippenstein, Megan Shabram, A. James Friedson, Mark S. Marley, Richard S. Freedman

Abstract: We have developed 1-D photochemical and thermochemical kinetics and diffusion models for the transiting exoplanets HD 189733b and HD 209458b to study the effects of disequilibrium chemistry on the atmospheric composition of "hot Jupiters." Here we investigate the coupled chemistry of neutral carbon, hydrogen, oxygen, and nitrogen species, and we compare the model results with existing transit and… ▽ More We have developed 1-D photochemical and thermochemical kinetics and diffusion models for the transiting exoplanets HD 189733b and HD 209458b to study the effects of disequilibrium chemistry on the atmospheric composition of "hot Jupiters." Here we investigate the coupled chemistry of neutral carbon, hydrogen, oxygen, and nitrogen species, and we compare the model results with existing transit and eclipse observations. We find that the vertical profiles of molecular constituents are significantly affected by transport-induced quenching and photochemistry, particularly on cooler HD 189733b; however, the warmer stratospheric temperatures on HD 209458b can help maintain thermochemical equilibrium and reduce the effects of disequilibrium chemistry. For both planets, the methane and ammonia mole fractions are found to be enhanced over their equilibrium values at pressures of a few bar to less than a mbar due to transport-induced quenching, but CH4 and NH3 are photochemically removed at higher altitudes. Atomic species, unsaturated hydrocarbons (particularly C2H2), some nitriles (particularly HCN), and radicals like OH, CH3, and NH2 are enhanced overequilibrium predictions because of quenching and photochemistry. In contrast, CO, H2O, N2, and CO2 more closely follow their equilibrium profiles, except at pressures < 1 microbar, where CO, H2O, and N2 are photochemically destroyed and CO2 is produced before its eventual high-altitude destruction. The enhanced abundances of HCN, CH4, and NH3 in particular are expected to affect the spectral signatures and thermal profiles HD 189733b and other, relatively cool, close-in transiting exoplanets. We examine the sensitivity of our results to the assumed temperature structure and eddy diffusion coefficientss and discuss further observational consequences of these models. △ Less

Submitted 16 May, 2011; v1 submitted 31 January, 2011; originally announced February 2011.

Comments: 40 pages, 16 figures, accepted for publication in Astrophysical Journal

Journal ref: Astrophys. J., 737, 15 (2011)

Transmission Spectra of Transiting Planet Atmospheres: Model Validation and Simulations of the Hot Neptune GJ 436b for JWST

Authors: Megan Shabram, Jonathan J. Fortney, Thomas P. Greene, Richard S. Freedman

Abstract: We explore the transmission spectrum of the Neptune-class exoplanet GJ 436b, including the possibility that its atmospheric opacity is dominated by a variety of non- equilibrium chemical products. We also validate our transmission code by demonstrating close agreement with analytic models that use only Rayleigh scattering or water vapor opacity. We find broad disagreement with radius variations pr… ▽ More We explore the transmission spectrum of the Neptune-class exoplanet GJ 436b, including the possibility that its atmospheric opacity is dominated by a variety of non- equilibrium chemical products. We also validate our transmission code by demonstrating close agreement with analytic models that use only Rayleigh scattering or water vapor opacity. We find broad disagreement with radius variations predicted by another published model. For GJ 436b, the relative coolness of the planet's atmosphere, along with its implied high metallicity, may make it dissimilar in character compared to "hot Jupiters." Some recent observational and modeling efforts suggest low relative abundances of H2O and CH4 present in GJ 436b's atmosphere, compared to calculations from equilibrium chemistry. We include these characteristics in our models and examine the effects of absorption from methane-derived higher order hydrocarbons. Significant absorption from HCN and C2H2 are found throughout the infrared, while C2H4 and C2H6 are less easily seen. We perform detailed simulations of JWST observations, including all likely noise sources, and find that we will be able to constrain chemical abundance regimes from this planet's transmission spectrum. For instance, the width of the features at 1.5, 3.3, and 7 microns indicates the amount of HCN versus C2H2 present. The NIRSpec prism mode will be useful due to its large spectral range and the relatively large number of photo-electrons recorded per spectral resolution element. However, extremely bright host stars like GJ 436 may be better observed with a higher spectroscopic resolution mode in order to avoid detector saturation. We find that observations with the MIRI low resolution spectrograph should also have high signal-to-noise in the 5 - 10 micron range due to the brightness of the star and the relatively low spectral resolution (R ~ 100) of this mode. △ Less

Submitted 23 November, 2010; v1 submitted 12 October, 2010; originally announced October 2010.

Comments: 33 pages, 12 figures, Accepted to ApJ

Atmospheric Circulation of Eccentric Hot Neptune GJ436b

Authors: Nikole K. Lewis, Adam P. Showman, Jonathan J. Fortney, Mark S. Marley, Richard S. Freedman, Katharina Lodders

Abstract: 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 similar to Neptune (20-60 times solar abundances), which makes it an ideal target to study the effects of atmospheric metallicity on dynamics and radia… ▽ More 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 similar 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. △ Less

Submitted 5 October, 2010; v1 submitted 17 July, 2010; originally announced July 2010.

Comments: 25 pages, 13 figures

Journal ref: The Astrophysical Journal, Volume 720, Issue 1, pp. 344-356 (2010)

Transmission Spectra of Three-Dimensional Hot Jupiter Model Atmospheres

Authors: J. J. Fortney, M. Shabram, A. P. Showman, Y. Lian, R. S. Freedman, M. S. Marley, N. K. Lewis

Abstract: We compute models of the transmission spectra of planets HD 209458b, HD 189733b, and generic hot Jupiters. We examine the effects of temperature, surface gravity, and metallicity for the generic planets as a guide to understanding transmission spectra in general. We find that carbon dioxide absorption at 4.4 and 15 microns is prominent at high metallicity, and is a clear metallicity indicator. F… ▽ More We compute models of the transmission spectra of planets HD 209458b, HD 189733b, and generic hot Jupiters. We examine the effects of temperature, surface gravity, and metallicity for the generic planets as a guide to understanding transmission spectra in general. We find that carbon dioxide absorption at 4.4 and 15 microns is prominent at high metallicity, and is a clear metallicity indicator. For HD 209458b and HD 189733b, we compute spectra for both one-dimensional and three-dimensional model atmospheres and examine the differences between them. The differences are usually small, but can be large if atmospheric temperatures are near important chemical abundance boundaries. The calculations for the 3D atmospheres, and their comparison with data, serve as constraints on these dynamical models that complement the secondary eclipse and light curve data sets. For HD 209458b, even if TiO and VO gases are abundant on the day side, their abundances can be considerably reduced on the cooler planetary limb. However, given the predicted limb temperatures and TiO abundances, the model's optical opacity is too high. For HD 189733b we find a good match with some infrared data sets and constrain the altitude of a postulated haze layer. For this planet, substantial differences can exist between the transmission spectra of the leading and trailing hemispheres, which is an excellent probe of carbon chemistry. In thermochemical equilibrium, the cooler leading hemisphere is methane-dominated, and the hotter trailing hemisphere is CO-dominated, but these differences may be eliminated by non-equilibrium chemistry due to vertical mixing. It may be possible to constrain the carbon chemistry of this planet, and its spatial variation, with JWST. △ Less

Submitted 11 December, 2009; originally announced December 2009.

Comments: Revised for ApJ. Submitted 9/16/2009. 20 pages, including figures

Journal ref: Astrophys.J.709:1396-1406,2010

Atmospheric Sulfur Photochemistry on Hot Jupiters

Authors: K. Zahnle, M. S. Marley, R. S. Freedman, K. Lodders, J. J. Fortney

Abstract: We develop a new 1D photochemical kinetics code to address stratospheric chemistry and stratospheric heating in hot Jupiters. Here we address optically active S-containing species and CO2 at 1200 < T < 2000 K. HS (mercapto) and S2 are highly reactive species that are generated photochemically and thermochemically from H2S with peak abundances between 1-10 mbar. S2 absorbs UV between 240 and 340… ▽ More We develop a new 1D photochemical kinetics code to address stratospheric chemistry and stratospheric heating in hot Jupiters. Here we address optically active S-containing species and CO2 at 1200 < T < 2000 K. HS (mercapto) and S2 are highly reactive species that are generated photochemically and thermochemically from H2S with peak abundances between 1-10 mbar. S2 absorbs UV between 240 and 340 nm and is optically thick for metallicities [SH] > 0 at T > 1200 K. HS is probably more important than S2, as it is generally more abundant than S2 under hot Jupiter conditions and it absorbs at somewhat redder wavelengths. We use molecular theory to compute an HS absorption spectrum from sparse available data and find that HS should absorb strongly between 300 and 460 nm, with absorption at the longer wavelengths being temperature sensitive. When the two absorbers are combined, radiative heating (per kg of gas) peaks at 100 microbars, with a total stratospheric heating of about 8 x 10^4 W/m^2 for a jovian planet orbiting a solar-twin at 0.032 AU. Total heating is insensitive to metallicity. The CO2 mixing ratio is a well-behaved quadratic function of metallicity, ranging from 1.6 x 10^-8 to 1.6 x 10^-4 for -0.3 < [M/H] < 1.7. CO2 is insensitive to insolation, vertical mixing, temperature (1200 < T <2000 K), and gravity. The photochemical calculations confirm that CO2 should prove a useful probe of planetary metallicity. △ Less

Submitted 26 June, 2009; v1 submitted 9 March, 2009; originally announced March 2009.

Comments: Astrophysical Journal Lett. in press; important revision includes effect of updated thermodynamic data and a new opacity source

Journal ref: Astrophys.J.701:L20-L24,2009

Atmospheric circulation of hot Jupiters: Coupled radiative-dynamical general circulation model simulations of HD 189733b and HD 209458b

Authors: Adam P. Showman, Jonathan J. Fortney, Yuan Lian, Mark S. Marley, Richard S. Freedman, Heather A. Knutson, David Charbonneau

Abstract: 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 Planetary Atmospheric Radiation and Circulation (SPARC) model, adopts the MITgcm for the dynamics and uses the radiative model of McKay, Marley, Fortney, an… ▽ More 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 Planetary 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) △ Less

Submitted 24 June, 2009; v1 submitted 11 September, 2008; originally announced September 2008.

Comments: 20 pages (emulate-apj format), 21 figures, final version now published in ApJ. Includes expanded discussion of radiative-transfer methods and two new figures

Journal ref: Astrophys.J.699:564-584,2009

Atmospheric Parameters of Field L and T Dwarfs

Authors: Michael C. Cushing, Mark S. Marley, D. Saumon, Brandon C. Kelly, William D. Vacca, John T. Rayner, Richard S. Freedman, Katharina Lodders, Thomas L. Roellig

Abstract: We present an analysis of the 0.95-14.5 micron spectral energy distributions of nine field ultracool dwarfs with spectral types ranging from L1 to T4.5. Effective temperatures, gravities, and condensate cloud sedimentation efficiencies are derived by comparing the data to synthetic spectra computed from atmospheric models that self-consistently include the formation of condensate clouds. Derived… ▽ More We present an analysis of the 0.95-14.5 micron spectral energy distributions of nine field ultracool dwarfs with spectral types ranging from L1 to T4.5. Effective temperatures, gravities, and condensate cloud sedimentation efficiencies are derived by comparing the data to synthetic spectra computed from atmospheric models that self-consistently include the formation of condensate clouds. Derived effective temperatures decrease steadily through the L1 to T4.5 spectral types and we confirm that the effective temperatures of ultracool dwarfs at the L/T transition are nearly constant, decreasing by only ~200 K from spectral types L7.5 to T4.5. The two objects in our sample with very red J-Ks colors are best fitted with synthetic spectra that have thick clouds which hints at a possible correlation between the near-infrared colors of L dwarfs and the condensate cloud properties. The fits to the two T dwarfs in our sample (T2 and T4.5) also suggest that the clouds become thinner in this spectral class, in agreement with previous studies. Restricting the fits to narrower wavelength ranges (i.e., individual photometric bands) almost always yields excellent agreement between the data and models. Limitations in our knowledge of the opacities of key absorbers such as FeH, VO, and CH4 at certain wavelengths remain obvious, however. The effective temperatures obtained by fitting the narrower wavelength ranges can show a large scatter compared to the values derived by fitting the full spectral energy distributions; deviations are typically ~200 K and in the worst cases, up to 700 K. △ Less

Submitted 6 November, 2007; originally announced November 2007.

Comments: 44 pages, 24 figures, accepted to ApJ on 2007 Nov 5

Journal ref: Astrophys.J.678:1372-1395,2008

A Unified Theory for the Atmospheres of the Hot and Very Hot Jupiters: Two Classes of Irradiated Atmospheres

Authors: Jonathan J. Fortney, Katharina Lodders, Mark S. Marley, Richard S. Freedman

Abstract: We highlight the importance of gaseous TiO and VO opacity on the highly irradiated close-in giant planets. The atmospheres of these planets naturally fall into two classes that are somewhat analogous to the M- and L-type dwarfs. Those that are warm enough to have appreciable opacity due to TiO and VO gases we term the ``pM Class'' planets, and those that are cooler we term ``pL Class'' planets.… ▽ More We highlight the importance of gaseous TiO and VO opacity on the highly irradiated close-in giant planets. The atmospheres of these planets naturally fall into two classes that are somewhat analogous to the M- and L-type dwarfs. Those that are warm enough to have appreciable opacity due to TiO and VO gases we term the ``pM Class'' planets, and those that are cooler we term ``pL Class'' planets. We calculate model atmospheres for these planets, including pressure-temperature profiles, spectra, and characteristic radiative time constants. We show that pM Class planets have hot stratospheres $\sim$2000 K and appear ``anomalously'' bright in the mid infrared secondary eclipse, as was recently found for planets HD 149026b and HD 209458b. This class of planets absorbs incident flux and emits thermal flux from high in their atmospheres. Consequently, they will have large day/night temperature contrasts and negligible phase shifts between orbital phase and thermal emission light curves, because radiative timescales are much shorter than possible dynamical timescales. The pL Class planets absorb incident flux deeper in the atmosphere where atmospheric dynamics will more readily redistribute absorbed energy. This will lead to cooler day sides, warmer night sides, and larger phase shifts in thermal emission light curves. Around a Sun-like primary this boundary occurs at $\sim$0.04-0.05 AU. The eccentric transiting planets HD 147506b and HD 17156b alternate between the classes. Thermal emission in the optical from pM Class planets is significant red-ward of 400 nm, making these planets attractive targets for optical detection. The difference in the observed day/night contrast between ups Andromeda b (pM Class) and HD 189733b (pL Class) is naturally explained in this scenario. (Abridged.) △ Less

Submitted 8 December, 2007; v1 submitted 12 October, 2007; originally announced October 2007.

Comments: Accepted to the Astrophysical Journal

Journal ref: Astrophys.J.678:1419-1435,2008

Line and Mean Opacities for Ultracool Dwarfs and Extrasolar Planets

Authors: Richard S. Freedman, Mark S. Marley, Katharina Lodders

Abstract: Opacities and chemical abundance data are crucial ingredients of ultracool dwarf and extrasolar giant planet atmosphere models. We report here on the detailed sources of molecular opacity data employed by our group for this application. We also present tables of Rosseland and Planck mean opacities which are of use in some studies of the atmospheres, interiors, and evolution of planets and brown… ▽ More Opacities and chemical abundance data are crucial ingredients of ultracool dwarf and extrasolar giant planet atmosphere models. We report here on the detailed sources of molecular opacity data employed by our group for this application. We also present tables of Rosseland and Planck mean opacities which are of use in some studies of the atmospheres, interiors, and evolution of planets and brown dwarfs. For the tables presented here we have included the opacities of important atomic and molecular species, including the alkali elements, pressure induced absorption by hydrogen, and other significant opacity sources but neglect opacity from condensates. We report for each species how we have assembled molecular line data from a combination of public databases, laboratory data that is not yet in the public databases, and our own numerical calculations. We combine these opacities with abundances computed from a chemical equilibrium model using recently revised solar abundances to compute mean opacities. The chemical equilibrium calculation accounts for the settling of condensates in a gravitational field, and is applicable to ultracool dwarf and extrasolar planetary atmospheres, but not circumstellar disks. We find that the inclusion of alkali atomic opacity substantially increases the mean opacities over those currently in the literature at densities relevant to the atmospheres and interiors of giant planets and brown dwarfs. We provide our opacity tables for public use and discuss their limitations. △ Less

Submitted 19 July, 2007; v1 submitted 15 June, 2007; originally announced June 2007.

Comments: 32 pages, 6 figures, accepted for publication in Astrphysical Journal Supplement, units revised, now includes on-line only data files (available in "source" package from "other formats" link)

Physical and Spectral Characteristics of the T8 and Later-Type Dwarfs

Authors: S. K. Leggett, M. S. Marley, R. Freedman, D. Saumon, Michael C. Liu, T. R. Geballe, D. A. Golimowski, D. C. Stephens

Abstract: We use new and published near-IR spectra, with synthetic spectra, to derive physical properties of three of the latest-type T dwarfs. A new R~1700 spectrum of the T7.5 dwarf HD 3651B, with existing data, allows a detailed comparison to the well-studied and very similar dwarf, Gl 570D. We find that HD 3651B has both higher gravity and metallicity than Gl 570D, with Teff=820-830K, log g= 5.4-5.5,… ▽ More We use new and published near-IR spectra, with synthetic spectra, to derive physical properties of three of the latest-type T dwarfs. A new R~1700 spectrum of the T7.5 dwarf HD 3651B, with existing data, allows a detailed comparison to the well-studied and very similar dwarf, Gl 570D. We find that HD 3651B has both higher gravity and metallicity than Gl 570D, with Teff=820-830K, log g= 5.4-5.5, [m/H]= +0.2 and Kzz=10^4cm^2/s. Its age is 8-12 Gyr and its implied mass is 60-70 M_Jup. We perform a similar analyis of the T8 and T7.5 dwarfs 2MASS J09393548-2448279 and 2MASS J11145133-2618235 using published data, comparing them to the well-studied T8, 2MASS J04151954-0935066. We find that the two dwarfs have the same Teff as the reference dwarf, and similar gravities, but lower metallicities. The parameters are Teff=725-775K and [m/H]= -0.3; log g=5.3-5.45 for 2MASS J09393548-2448279 and log g=5.0-5.3 for 2MASS J11145133- 261823. The age and mass are ~10Gyr and 60M_Jup for 2MASS J09393548-2448279, and ~5 Gyr and 40M_Jup for 2MASS J11145133-261823. A serious limitation is the incompleteness of the line lists of CH4 and NH3 at lambda <1.7um. Spectra of Saturn and Jupiter, and of laboratory CH4 and NH3 gas, suggest that NH3 features in the Y- and J-bands may be useful as indicators of the next cooler spectral type, and not features in the H- and K-bands as previously thought. However large uncertainties remain, as the abundance of NH3 is likely to be significantly below the chemical equilibrium value, and inclusion of laboratory NH3 opacities predicts band shapes that are discrepant with existing data. It is possible that the T spectral class will have to be extended to low temperatures around 400K, when water clouds condense in the atmosphere [abridged]. △ Less

Submitted 17 May, 2007; originally announced May 2007.

Comments: 34 pages including 10 figures and two tables; accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.667:537-548,2007

Physical parameters of two very cool T dwarfs

Authors: D. Saumon, M. S. Marley, S. K. Leggett, T. R. Geballe, D. Stephens, D. A. Golimowski, M. C. Cushing, X. Fan, J. T. Rayner, K. Lodders, R. S. Freedman

Abstract: We present new infrared spectra of the T8 brown dwarf 2MASS J04151954-0935066: 2.9-4.1 micron spectra obtained with the Infrared Camera and Spectrograph on the Subaru Telescope, and 5.2-14.5 micron spectra obtained with the Infrared Spectrograph on the Spitzer Space Telescope. We use these data and models to determine an accurate bolometric luminosity of log L_bol/L_sun=-5.67 and to constrain th… ▽ More We present new infrared spectra of the T8 brown dwarf 2MASS J04151954-0935066: 2.9-4.1 micron spectra obtained with the Infrared Camera and Spectrograph on the Subaru Telescope, and 5.2-14.5 micron spectra obtained with the Infrared Spectrograph on the Spitzer Space Telescope. We use these data and models to determine an accurate bolometric luminosity of log L_bol/L_sun=-5.67 and to constrain the effective temperature, gravity, mass and age to 725-775K, log g=5.00-5.37, M=33-58 M_Jupiter and age=3-10Gyr. We perform the same analysis using published 0.6-15 micron spectra for the T7.5 dwarf 2MASS J12171110-0311131, for which we find a metal-rich composition ([Fe/H]~0.3) and log L_bol/L_sun=-5.31, T_eff=850-950K, log g=4.80-5.42, M=25-66M_Jupiter and age=1-10Gyr. These luminosities and effective temperatures straddle those determined with the same method and models for Gl 570D by Saumon et al. (2006) and make 2MASS J04151954-0935066 the coolest and least luminous T dwarf with well-determined properties. We find that synthetic spectra generated by the models reproduce the observed red through mid-infrared spectra of 2MASS J04151954-0935066 and 2MASS J12171110-0311131 very well, except for known discrepancies which are most likely due to the incomplete CH4 opacities. Both objects show evidence of departures from strict chemical equilibrium and we discuss this result in the context of other late T dwarfs where disequilibrium phenomena have been observed. △ Less

Submitted 2 November, 2006; originally announced November 2006.

Comments: Accepted for publication in the Ap. J. 35 pages, including 12 figures (10 in color)

Report number: LA-UR-06-6933

Journal ref: Astrophys.J.656:1136-1149,2007