Showing 1–50 of 80 results for author: Lupu, R

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

Modeling the day-night temperature variations of ultra-hot Jupiters: confronting non-grey general circulation models and observations

Authors: Xianyu Tan, Thaddeus D. Komacek, Natasha E. Batalha, Drake Deming, Roxana Lupu, Vivien Parmentier, Raymond T. Pierrehumbert

Abstract: Ultra-hot Jupiters (UHJs) are natural laboratories to study extreme physics in planetary atmospheres and their rich observational data sets are yet to be confronted with models with varying complexities at a population level. In this work, we update the general circulation model of Tan & Komacek (2019) to include a non-grey radiative transfer scheme and apply it to simulate the realistic thermal s… ▽ More Ultra-hot Jupiters (UHJs) are natural laboratories to study extreme physics in planetary atmospheres and their rich observational data sets are yet to be confronted with models with varying complexities at a population level. In this work, we update the general circulation model of Tan & Komacek (2019) to include a non-grey radiative transfer scheme and apply it to simulate the realistic thermal structures, phase-dependent spectra, and wavelength-dependent phase curves of UHJs. We performed grids of models over a large range of equilibrium temperatures and rotation periods for varying assumptions, showing that the fractional day-night brightness temperature differences remain almost constant or slightly increase with increasing equilibrium temperature from the visible to mid-infrared wavelengths. This differs from previous work primarily due to the increasing planetary rotation rate with increasing equilibrium temperature for fixed host star type. Radiative effects of varying atmospheric compositions become more significant in dayside brightness temperature in longer wavelengths. Data-model comparisons of dayside brightness temperatures and phase curve amplitudes as a function of equilibrium temperature are in broad agreement. Observations show a large scatter compared to models even with a range of different assumptions, indicating significantly varying intrinsic properties in the hot Jupiter population. Our cloud-free models generally struggle to match all observations for individual targets with a single set of parameter choices, indicating the need for extra processes for understanding the heat transport of UHJs. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: Accepted to MNRAS; data underlying this article is available in Zenodo https://doi.org/10.5281/zenodo.10121933

Origin of Life Molecules in the Atmosphere After Big Impacts on the Early Earth

Authors: Nicholas F. Wogan, David C. Catling, Kevin J. Zahnle, Roxana Lupu

Abstract: The origin of life on Earth would benefit from a prebiotic atmosphere that produced nitriles, like HCN, which enable ribonucleotide synthesis. However, geochemical evidence suggests that Hadean air was relatively oxidizing with negligible photochemical production of prebiotic molecules. These paradoxes are resolved by iron-rich asteroid impacts that transiently reduced the entire atmosphere, allow… ▽ More The origin of life on Earth would benefit from a prebiotic atmosphere that produced nitriles, like HCN, which enable ribonucleotide synthesis. However, geochemical evidence suggests that Hadean air was relatively oxidizing with negligible photochemical production of prebiotic molecules. These paradoxes are resolved by iron-rich asteroid impacts that transiently reduced the entire atmosphere, allowing nitriles to form in subsequent photochemistry. Here, we investigate impact-generated reducing atmospheres using new time-dependent, coupled atmospheric chemistry and climate models, which account for gas-phase reactions and surface-catalysis. The resulting H$_2$-, CH$_4$- and NH$_3$-rich atmospheres persist for millions of years, until hydrogen escapes to space. HCN and HCCCN production and rainout to the surface can reach $10^9$ molecules cm$^{-2}$ s$^{-1}$ in hazy atmospheres with a mole ratio of $\mathrm{CH_4} / \mathrm{CO_2} > 0.1$. Smaller $\mathrm{CH_4} / \mathrm{CO_2}$ ratios produce HCN rainout rates $< 10^5$ molecules cm$^{-2}$ s$^{-1}$, and negligible HCCCN. The minimum impactor mass that creates atmospheric $\mathrm{CH_4} / \mathrm{CO_2} > 0.1$ is $4 \times 10^{20}$ to $5 \times 10^{21}$ kg (570 to 1330 km diameter), depending on how efficiently iron reacts with a steam atmosphere, the extent of atmospheric equilibration with an impact-induced melt pond, and the surface area of nickel that catalyzes CH$_4$ production. Alternatively, if steam permeates and deeply oxidizes crust, impactors $\sim 10^{20}$ kg could be effective. Atmospheres with copious nitriles have $> 360$ K surface temperatures, perhaps posing a challenge for RNA longevity, although cloud albedo can produce cooler climates. Regardless, post-impact cyanide can be stockpiled and used in prebiotic schemes after hydrogen has escaped to space. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: Accepted by Planetary Science Journal on 7/4/23

A reflective, metal-rich atmosphere for GJ 1214b from its JWST phase curve

Authors: Eliza M. -R. Kempton, Michael Zhang, Jacob L. Bean, Maria E. Steinrueck, Anjali A. A. Piette, Vivien Parmentier, Isaac Malsky, Michael T. Roman, Emily Rauscher, Peter Gao, Taylor J. Bell, Qiao Xue, Jake Taylor, Arjun B. Savel, Kenneth E. Arnold, Matthew C. Nixon, Kevin B. Stevenson, Megan Mansfield, Sarah Kendrew, Sebastian Zieba, Elsa Ducrot, Achrène Dyrek, Pierre-Olivier Lagage, Keivan G. Stassun, Gregory W. Henry , et al. (8 additional authors not shown)

Abstract: There are no planets intermediate in size between Earth and Neptune in our Solar System, yet these objects are found around a substantial fraction of other stars. Population statistics show that close-in planets in this size range bifurcate into two classes based on their radii. It is hypothesized that the group with larger radii (referred to as "sub-Neptunes") is distinguished by having hydrogen-… ▽ More There are no planets intermediate in size between Earth and Neptune in our Solar System, yet these objects are found around a substantial fraction of other stars. Population statistics show that close-in planets in this size range bifurcate into two classes based on their radii. It is hypothesized that the group with larger radii (referred to as "sub-Neptunes") is distinguished by having hydrogen-dominated atmospheres that are a few percent of the total mass of the planets. GJ 1214b is an archetype sub-Neptune that has been observed extensively using transmission spectroscopy to test this hypothesis. However, the measured spectra are featureless, and thus inconclusive, due to the presence of high-altitude aerosols in the planet's atmosphere. Here we report a spectroscopic thermal phase curve of GJ 1214b obtained with JWST in the mid-infrared. The dayside and nightside spectra (average brightness temperatures of 553 $\pm$ 9 and 437 $\pm$ 19 K, respectively) each show >3$σ$ evidence of absorption features, with H$_2$O as the most likely cause in both. The measured global thermal emission implies that GJ 1214b's Bond albedo is 0.51 $\pm$ 0.06. Comparison between the spectroscopic phase curve data and three-dimensional models of GJ 1214b reveal a planet with a high metallicity atmosphere blanketed by a thick and highly reflective layer of clouds or haze. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: Published online in Nature on May 10, 2023

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

A 1.46-2.48 $μ$m Spectroscopic Atlas of a T6 Dwarf (1060 K) Atmosphere with IGRINS: First Detections of H$_2$S and H$_2$, and Verification of H$_2$O, CH$_4$, and NH$_3$ Line Lists

Authors: Megan E. Tannock, Stanimir Metchev, Callie E. Hood, Gregory N. Mace, Jonathan J. Fortney, Caroline V. Morley, Daniel T. Jaffe, Roxana Lupu

Abstract: We present Gemini South/IGRINS observations of the 1060 K T6 dwarf 2MASS J08173001$-$6155158 with unprecedented resolution ($R\equivλ/Δλ=45\,000$) and signal-to-noise ratio (SNR > 200) for a late-type T dwarf. We use this benchmark observation to test the reliability of molecular line lists used up-to-date atmospheric models. We determine which spectroscopic regions should be used to estimate the… ▽ More We present Gemini South/IGRINS observations of the 1060 K T6 dwarf 2MASS J08173001$-$6155158 with unprecedented resolution ($R\equivλ/Δλ=45\,000$) and signal-to-noise ratio (SNR > 200) for a late-type T dwarf. We use this benchmark observation to test the reliability of molecular line lists used up-to-date atmospheric models. We determine which spectroscopic regions should be used to estimate the parameters of cold brown dwarfs and, by extension, exoplanets. We present a detailed spectroscopic atlas with molecular identifications across the $H$ and $K$ bands of the near-infrared. We find that water (H$_2$O) line lists are overall reliable. We find the most discrepancies amongst older methane (CH$_4$) line lists, and that the most up-to-date CH$_4$ line lists correct many of these issues. We identify individual ammonia (NH$_3$) lines, a hydrogen sulfide (H$_2$S) feature at 1.5900 $μ$m, and a molecular hydrogen (H$_2$) feature at 2.1218 $μ$m. These are the first unambiguous detections of H$_2$S and H$_2$ absorption features in an extra-solar atmosphere. With the H$_2$ detection, we place an upper limit on the atmospheric dust concentration of this T6 dwarf: at least 500 times less than the interstellar value, implying that the atmosphere is effectively dust-free. We additionally identify several features that do not appear in the model spectra. Our assessment of the line lists is valuable for atmospheric model applications to high-dispersion, low-SNR, high-background spectra, such as an exoplanet around a star. We demonstrate a significant enhancement in the detection of the CH$_4$ absorption signal in this T6 dwarf with the most up-to-date line lists. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 45 pages, 63 figures, 5 tables, Accepted to MNRAS

The impact of ultraviolet heating and cooling on the dynamics and observability of lava planet atmospheres

Authors: T. Giang Nguyen, Nicolas B. Cowan, Raymond T. Pierrehumbert, Roxana E. Lupu, John E. Moores

Abstract: Lava planets have non-global, condensible atmospheres similar to icy bodies within the solar system. Because they depend on interior dynamics, studying the atmospheres of lava planets can lead to understanding unique geological processes driven by their extreme environment. Models of lava planet atmospheres have thus far focused on either radiative transfer or hydrodynamics. In this study, we coup… ▽ More Lava planets have non-global, condensible atmospheres similar to icy bodies within the solar system. Because they depend on interior dynamics, studying the atmospheres of lava planets can lead to understanding unique geological processes driven by their extreme environment. Models of lava planet atmospheres have thus far focused on either radiative transfer or hydrodynamics. In this study, we couple the two processes by introducing ultraviolet and infrared radiation to a turbulent boundary layer model. We also test the effect of different vertical temperature profiles on atmospheric dynamics. Results from the model show that UV radiation affects the atmosphere much more than IR. UV heating and cooling work together to produce a horizontally isothermal atmosphere away from the sub-stellar point regardless of the vertical temperature profile. We also find that stronger temperature inversions induce stronger winds and hence cool the atmosphere. Our simulated transmission spectra of the bound atmosphere show a strong SiO feature in the UV that would be challenging to observe in the planet's transit spectrum due to the precision required. Our simulated emission spectra are more promising, with significant SiO spectral features at 4.5 and 9 $μ$m that can be observed with the James Webb Space Telescope. Different vertical temperature profiles produce discernible dayside emission spectra, but not in the way one would expect. △ Less

Submitted 9 May, 2022; v1 submitted 28 April, 2022; originally announced April 2022.

Comments: Accepted to MNRAS

K2 and Spitzer phase curves of the rocky ultra-short-period planet K2-141 b hint at a tenuous rock vapor atmosphere

Authors: S. Zieba, M. Zilinskas, L. Kreidberg, T. G. Nguyen, Y. Miguel, N. B. Cowan, R. Pierrehumbert, L. Carone, L. Dang, M. Hammond, T. Louden, R. Lupu, L. Malavolta, K. B. Stevenson

Abstract: K2-141 b is a transiting, small (1.5 Re) ultra-short-period (USP) planet discovered by Kepler orbiting a K-dwarf host star every 6.7 hours. The planet's high surface temperature makes it an excellent target for thermal emission observations. Here we present 65 hours of continuous photometric observations of K2-141 b collected with Spitzer's IRAC Channel 2 at 4.5 micron spanning 10 full orbits of t… ▽ More K2-141 b is a transiting, small (1.5 Re) ultra-short-period (USP) planet discovered by Kepler orbiting a K-dwarf host star every 6.7 hours. The planet's high surface temperature makes it an excellent target for thermal emission observations. Here we present 65 hours of continuous photometric observations of K2-141 b collected with Spitzer's IRAC Channel 2 at 4.5 micron spanning 10 full orbits of the planet. We measure an infrared eclipse depth of 143 +/- 39 ppm and a peak to trough amplitude variation of 121 +/- 43 ppm. The best fit model to the Spitzer data shows no significant thermal hotspot offset, in contrast to the previously observed offset for the well-studied USP planet 55 Cnc e. We also jointly analyze the new Spitzer observations with the photometry collected by Kepler during two separate K2 campaigns. We model the planetary emission with a range of toy models that include a reflective and a thermal contribution. With a two-temperature model, we measure a dayside temperature of 2049 +/- 361 K and a night-side temperature that is consistent with zero (Tp,n < 1712 K at 2 sigma). Models with a steep dayside temperature gradient provide a better fit to the data than a uniform dayside temperature (DeltaBIC = 22.2). We also find evidence for a non-zero geometric albedo Ag = 0.28 +/- 0.07. We also compare the data to a physically motivated, pseudo-2D rock vapor model and a 1D turbulent boundary layer model. Both models fit the data well. Notably, we find that the optical eclipse depth can be explained by thermal emission from a hot inversion layer, rather than reflected light. A thermal inversion may also be responsible for the deep optical eclipse observed for another USP, Kepler-10 b. Finally, we significantly improve the ephemerides for K2-141 b and c, which will facilitate further follow-up observations of this interesting system with state-of-the-art observatories like JWST. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: 31 pages, 27 figures, 12 tables, accepted for publication in A&A

Journal ref: A&A 664, A79 (2022)

Observability of Evaporating Lava Worlds

Authors: Mantas Zilinskas, Christiaan van Buchem, Yamila Miguel, Amy Louca, Roxana Lupu, Sebastian Zieba, Wim van Westrenen

Abstract: Lava worlds belong to a class of short orbital period planets reaching dayside temperatures high enough to melt their silicate crust. Theory predicts that the resulting lava oceans outgas their volatile components, attaining equilibrium with the overlying vapour. This creates a tenuous, silicate-rich atmosphere that may be confined to the permanent dayside of the planet. The launch of JWST will pr… ▽ More Lava worlds belong to a class of short orbital period planets reaching dayside temperatures high enough to melt their silicate crust. Theory predicts that the resulting lava oceans outgas their volatile components, attaining equilibrium with the overlying vapour. This creates a tenuous, silicate-rich atmosphere that may be confined to the permanent dayside of the planet. The launch of JWST will provide the much needed sensitivity and spectral coverage to characterise these worlds. In this paper, we assess the observability of characterisable spectral features by self-consistently modelling silicate atmospheres for all the currently confirmed targets having sufficient substellar temperatures (> 1500 K). To achieve this we use outgassed equilibrium chemistry and radiative transfer methods to compute temperature-pressure profiles, atmospheric chemical compositions and emission spectra. We explore varying melt compositions, free of highly-volatile elements, accounting for possible atmospheric evolution. Our models include a large number of neutral and ionic species, as well as all up-to-date opacities. The results indicate that SiO and SiO2 infrared features are the best, unique identifiers of silicate atmospheres, detectable using MIRI instrument of JWST. Detection of these two species in emission would allow for strong constraints on atmospheric thermal structure and possibly the composition of the melt. We also propose that certain species, e.g., TiO, may be directly tied to different classes of melts, possibly revealing surface and interior dynamics. Currently, there are nearly a dozen confirmed lava planets ideal for characterisation of silicate atmospheres using JWST, with two of these already accepted for the initial General Observers program. △ Less

Submitted 19 April, 2022; v1 submitted 9 February, 2022; originally announced February 2022.

Comments: Accepted for publication in A&A Preprint: https://doi.org/10.1051/0004-6361/202142984

Journal ref: A&A 661, A126 (2022)

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

Impacts of Water Latent Heat on the Thermal Structure of Ultra-Cool Objects: Brown Dwarfs and Free-Floating Planets

Authors: Shih-Yun Tang, Tyler D. Robinson, Mark S. Marley, Natasha E. Batalha, Roxana Lupu, L. Prato

Abstract: Brown dwarfs are essential targets for understanding planetary and sub-stellar atmospheres across a wide range of thermal and chemical conditions. As surveys continue to probe ever deeper, and as observing capabilities continue to improve, the number of known Y dwarfs -- the coldest class of sub-stellar objects, with effective temperatures below about 600 K -- is rapidly growing. Critically, this… ▽ More Brown dwarfs are essential targets for understanding planetary and sub-stellar atmospheres across a wide range of thermal and chemical conditions. As surveys continue to probe ever deeper, and as observing capabilities continue to improve, the number of known Y dwarfs -- the coldest class of sub-stellar objects, with effective temperatures below about 600 K -- is rapidly growing. Critically, this class of ultra-cool objects has atmospheric conditions that overlap with Solar System worlds and, as a result, tools and ideas developed from studying Earth, Jupiter, Saturn and other nearby worlds are well-suited for application to sub-stellar atmospheres. To that end, we developed a one-dimensional (vertical) atmospheric structure model for ultra-cool objects that includes moist adiabatic convection, as this is an important process for many Solar System planets. Application of this model across a range of effective temperatures (350, 300, 250, 200 K), metallicities ([M/H] of 0.0, 0.5, 0.7, 1.5), and gravities (log $g$ of 4.0, 4.5, 4.7, 5.0) demonstrates strong impacts of water latent heat release on simulated temperature-pressure profiles. At the highest metallicities, water vapor mixing ratios reach an Earth-like 3%, with associated major alterations to the thermal structure in the atmospheric regions where water condenses. Spectroscopic and photometric signatures of metallicity and moist convection should be readily detectable at near- and mid-infrared wavelengths, especially with James Webb Space Telescope observations, and can help indicate the formation history of an object. △ Less

Submitted 16 August, 2021; v1 submitted 14 May, 2021; originally announced May 2021.

Comments: 16 pages, 10 figures, ApJ accepted version, data public available on zenodo

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

The Nancy Grace Roman Space Telescope Coronagraph Instrument (CGI) Technology Demonstration

Authors: N. Jeremy Kasdin, Vanessa P. Bailey, Bertrand Mennesson, Robert T. Zellem, Marie Ygouf, Jason Rhodes, Thomas Luchik, Feng Zhao, A J Eldorado Riggs, Young-Joon Seo, John Krist, Brian Kern, Hong Tang, Bijan Nemati, Tyler D. Groff, Neil Zimmerman, Bruce Macintosh, Margaret Turnbull, John Debes, Ewan S. Douglas, Roxana E. Lupu

Abstract: The Coronagraph Instrument (CGI) on the Nancy Grace Roman Space Telescope will demonstrate the high-contrast technology necessary for visible-light exoplanet imaging and spectroscopy from space via direct imaging of Jupiter-size planets and debris disks. This in-space experience is a critical step toward future, larger missions targeted at direct imaging of Earth-like planets in the habitable zone… ▽ More The Coronagraph Instrument (CGI) on the Nancy Grace Roman Space Telescope will demonstrate the high-contrast technology necessary for visible-light exoplanet imaging and spectroscopy from space via direct imaging of Jupiter-size planets and debris disks. This in-space experience is a critical step toward future, larger missions targeted at direct imaging of Earth-like planets in the habitable zones of nearby stars. This paper presents an overview of the current instrument design and requirements, highlighting the critical hardware, algorithms, and operations being demonstrated. We also describe several exoplanet and circumstellar disk science cases enabled by these capabilities. A competitively selected Community Participation Program team will be an integral part of the technology demonstration and could perform additional CGI observations beyond the initial tech demo if the instrument performance warrants it. △ Less

Submitted 2 March, 2021; originally announced March 2021.

Comments: SPIE Space Telescopes and Instrumentation 2020: Optical, Infrared, and Millimeter Wave, 114431U (15 December 2020)

Report number: Proc. SPIE 11443

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

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

Creation and Evolution of Impact-generated Reduced Atmospheres of Early Earth

Authors: Kevin Zahnle, Roxana Lupu, David Catling, Nick Wogan

Abstract: The origin of life on Earth seems to demand a highly reduced early atmosphere, rich in CH4, H2, and NH3, but geological evidence suggests that Earth's mantle has always been relatively oxidized and its emissions dominated by CO2 H2O, and N2. The paradox can be resolved by exploiting the reducing power inherent in the "late veneer," i.e., material accreted by Earth after the Moon-forming impact. Is… ▽ More The origin of life on Earth seems to demand a highly reduced early atmosphere, rich in CH4, H2, and NH3, but geological evidence suggests that Earth's mantle has always been relatively oxidized and its emissions dominated by CO2 H2O, and N2. The paradox can be resolved by exploiting the reducing power inherent in the "late veneer," i.e., material accreted by Earth after the Moon-forming impact. Isotopic evidence indicates that the late veneer consisted of extremely dry, highly reduced inner solar system materials, suggesting that Earth's oceans were already present when the late veneer came. The major primary product of reaction between the late veneer's iron and Earth's water was H2. Ocean vaporizing impacts generate high pressures and long cooling times that favor CH4 and NH3. Impacts too small to vaporize the oceans are much less productive of CH4 and NH3, unless (i) catalysts were available to speed their formation, or (ii) additional reducing power was extracted from pre-existing crustal or mantle materials. The transient H2-CH4 atmospheres evolve photochemically to generate nitrogenated hydrocarbons at rates determined by solar radiation and hydrogen escape, on timescales ranging up to tens of millions of years and with cumulative organic production ranging up to half a kilometer. Roughly one ocean of hydrogen escapes. The atmosphere after the methane's gone is typically H2 and CO rich, with eventual oxidation to CO2 rate-limited by water photolysis and hydrogen escape. △ Less

Submitted 25 February, 2020; v1 submitted 31 December, 2019; originally announced January 2020.

Comments: Submitted 31 Dec 2019 to "The Planetary Science Journal." Accepted 10 Feb 2020

Detecting and Characterizing Water Vapor in the Atmospheres of Earth Analogs through Observation of the 0.94 Micron Feature in Reflected Light

Authors: Adam J. R. W. Smith, Y. Katherina Feng, Jonathan J. Fortney, Tyler D. Robinson, Mark S. Marley, Roxana E. Lupu, Nikole K. Lewis

Abstract: The characterization of rocky, Earth-like planets is an important goal for future large ground- and space-based telescopes. In support of developing an efficient observational strategy, we have applied Bayesian statistical inference to interpret the albedo spectrum of cloudy true-Earth analogs that include a diverse spread in their atmospheric water vapor mixing ratios. We focus on detecting water… ▽ More The characterization of rocky, Earth-like planets is an important goal for future large ground- and space-based telescopes. In support of developing an efficient observational strategy, we have applied Bayesian statistical inference to interpret the albedo spectrum of cloudy true-Earth analogs that include a diverse spread in their atmospheric water vapor mixing ratios. We focus on detecting water-bearing worlds by characterizing their atmospheric water vapor content via the strong 0.94$\,μ$m H$_2$O absorption feature, with several observational configurations. Water vapor is an essential signpost when assessing planetary habitability, and determining its presence is important in vetting whether planets are suitable for hosting life. We find that R=140 spectroscopy of the absorption feature combined with a same-phase green optical photometric point at $0.525-0.575\,μ$m is capable of distinguishing worlds with less than $0.1\times$ Earth-like water vapor levels from worlds with $1\times$ Earth-like levels or greater at a signal-to-noise ratio of 5 or better with $2σ$ confidence. This configuration can differentiate between $0.01\times$ and $0.1\times$ Earth-like levels when the signal-to-noise ratio is 10 or better at the same confidence. However, strong constraints on the water vapor mixing ratio remained elusive with this configuration even at signal-to-noise of 15. We find that adding the same-phase optical photometric point does not significantly help characterize the H$_2$O mixing ratio, but does enable an upper limit on atmospheric ozone levels. Finally, we find that a 0.94$\,μ$m photometric point, instead of spectroscopy, combined with the green-optical point, fails to produce meaningful information about atmospheric water content. △ Less

Submitted 20 November, 2019; originally announced November 2019.

Comments: Accepted to AJ

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

Characterizing the Atmospheres of Irradiated Exoplanets at High Spectral Resolution

Authors: Diana Dragomir, Eliza Kempton, Jacob Bean, Ian Crossfield, Eric Gaidos, Nikole Lewis, Michael Line, Roxana Lupu, George Zhou

Abstract: The best-characterized exoplanets to date are planets on close-in transiting orbits around their host stars. The high level of irradiation and transiting geometry of these objects make them ideal targets for atmospheric investigations. However, the modest apertures of many current telescopes allow mostly low resolution spectra to be observed for transiting planets, failing to extract key physical… ▽ More The best-characterized exoplanets to date are planets on close-in transiting orbits around their host stars. The high level of irradiation and transiting geometry of these objects make them ideal targets for atmospheric investigations. However, the modest apertures of many current telescopes allow mostly low resolution spectra to be observed for transiting planets, failing to extract key physical and chemical properties of their atmospheres. Ground-based 30-meter class telescopes will set the stage for a substantial leap in our understanding of exoplanet atmospheres. We outline a two-pronged survey, recently submitted as a Key Science Project (KSP) for the US ELTs, which would yield unprecedented insight into the atmospheres of close-in exoplanets via combined observations with the GMT and TMT. (1) The first opportunity involves measuring the global-scale atmospheric circulation and planetary rotation for a sample of 40 hot Jupiters to glean insight into the unique radiative forcing regime governing highly-irradiated, tidally-locked giant planets. (2) The second opportunity involves measuring atmospheric mass-loss and extracting atmospheric composition and abundance ratio information for $\sim60$ sub-Neptunes and super-Earths (including candidate disintegrating planets) to constrain their formation and evolution histories. These efforts would be made possible by the unparalleled combination of high spectral resolution instrumentation and large aperture size of the ELTs. This survey would enable the first statistical study of atmospheric circulation in extrasolar giant planets, and would provide detections of trace gases and measurements of atmospheric escape in small-planet atmospheres, far exceeding the reach of \textit{JWST}. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: Science White Paper submitted to the Astro2020 Decadal Survey; 7 pages, 3 figures

New Frontiers for Terrestrial-sized to Neptune-sized Exoplanets In the Era of Extremely Large Telescopes

Authors: Ji Wang, Michael R. Meyer, Alan Boss, Laird Close, Thayne Currie, Diana Dragomir, Jonathan Fortney, Eric Gaidos, Yasuhiro Hasegawa, Irina Kitiashvili, Quinn Konopacky, Chien-Hsiu Lee, Nikole K. Lewis, Michael Liu, Roxana Lupu, Dimitri Mawet, Carl Melis, Mercedes Lopez-Morales, Caroline V. Morley, Chris Packham, Eliad Peretz, Andy Skemer, Mel Ulmer

Abstract: Surveys reveal that terrestrial- to Neptune-sized planets (1 $< R <$ 4 R$_{\rm{Earth}}$) are the most common type of planets in our galaxy. Detecting and characterizing such small planets around nearby stars holds the key to understanding the diversity of exoplanets and will ultimately address the ubiquitousness of life in the universe. The following fundamental questions will drive research in th… ▽ More Surveys reveal that terrestrial- to Neptune-sized planets (1 $< R <$ 4 R$_{\rm{Earth}}$) are the most common type of planets in our galaxy. Detecting and characterizing such small planets around nearby stars holds the key to understanding the diversity of exoplanets and will ultimately address the ubiquitousness of life in the universe. The following fundamental questions will drive research in the next decade and beyond: (1) how common are terrestrial to Neptune-sized planets within a few AU of their host star, as a function of stellar mass? (2) How does planet composition depend on planet mass, orbital radius, and host star properties? (3) What are the energy budgets, atmospheric dynamics, and climates of the nearest worlds? Addressing these questions requires: a) diffraction-limited spatial resolution; b) stability and achievable contrast delivered by adaptive optics; and c) the light-gathering power of extremely large telescopes (ELTs), as well as multi-wavelength observations and all-sky coverage enabled by a comprehensive US ELT Program. Here we provide an overview of the challenge, and promise of success, in detecting and comprehensively characterizing small worlds around the very nearest stars to the Sun with ELTs. This white paper extends and complements the material presented in the findings and recommendations published in the National Academy reports on Exoplanet Science Strategy and Astrobiology Strategy for the Search for Life in the Universe. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Astro2020 Science White Paper

The Effect of 3D Transport-Induced Disequilibrium Carbon Chemistry on the Atmospheric Structure and Phase Curves and Emission Spectra of Hot Jupiter HD 189733b

Authors: Maria E Steinrueck, Vivien Parmentier, Adam P Showman, Joshua D Lothringer, Roxana E Lupu

Abstract: On hot Jupiter exoplanets, strong horizontal and vertical winds should homogenize the abundances of the important absorbers CH$_4$ and CO much faster than chemical reactions restore chemical equilibrium. This effect, typically neglected in general circulation models (GCMs), has been suggested as explanation for discrepancies between observed infrared lightcurves and those predicted by GCMs: On the… ▽ More On hot Jupiter exoplanets, strong horizontal and vertical winds should homogenize the abundances of the important absorbers CH$_4$ and CO much faster than chemical reactions restore chemical equilibrium. This effect, typically neglected in general circulation models (GCMs), has been suggested as explanation for discrepancies between observed infrared lightcurves and those predicted by GCMs: On the nightsides of several hot Jupiters, GCMs predict outgoing fluxes that are too large, especially in the Spitzer 4.5 $μ$m band. We modified the SPARC/MITgcm to include disequilibrium abundances of CH$_4$, CO and H$_2$O by assuming that the CH$_4$/CO ratio is constant throughout the simulation domain. We ran simulations of hot Jupiter HD 189733b with 8 CH$_4$/CO ratios. In the more likely CO-dominated regime, we find temperature changes $\geq$50-100 K compared to the equilibrium chemistry case across large regions. This effect is large enough to affect predicted emission spectra and should thus be included in GCMs of hot Jupiters with equilibrium temperatures between 600K and 1300K. We find that spectra in regions with strong methane absorption, including the Spitzer 3.6 and 8 $μ$m bands, are strongly impacted by disequilibrium abundances. We expect chemical quenching to result in much larger nightside fluxes in the 3.6 $μ$m band, in stark contrast to observations. Meanwhile, we find almost no effect on predicted observations in the 4.5 $μ$m band, as the opacity changes due to CO and H$_2$O offset each other. We thus conclude that disequilibrium carbon chemistry cannot explain the observed low nightside fluxes in the 4.5 $μ$m band. △ Less

Submitted 21 May, 2019; v1 submitted 6 August, 2018; originally announced August 2018.

Comments: 28 pages, 12 figures; accepted for publication in ApJ. Key changes since v1: updated introduction and discussion to include recent literature more thoroughly; added paragraph and figure to clarify change in water abundance; ran simulations for additional 600 days simulation time (leading to small quantitative changes in some figures but no qualitative changes)

From thermal dissociation to condensation in the atmospheres of ultra hot Jupiters: WASP-121b in context

Authors: Vivien Parmentier, Mike R. Line, Jacob L. Bean, Megan Mansfield, Laura Kreidberg, Roxana Lupu, Channon Visscher, Jean-Michel Desert, Jonathan J. Fortney, Magalie Deleuil, Jacob Arcangeli, Adam P. Showman, Mark S. Marley

Abstract: A new class of exoplanets has emerged: the ultra hot Jupiters, the hottest close-in gas giants. Most of them have weaker than expected spectral features in the $1.1-1.7μm$ bandpass probed by HST/WFC3 but stronger spectral features at longer wavelengths probed by Spitzer. This led previous authors to puzzling conclusions about the thermal structures and chemical abundances of these planets. Using t… ▽ More A new class of exoplanets has emerged: the ultra hot Jupiters, the hottest close-in gas giants. Most of them have weaker than expected spectral features in the $1.1-1.7μm$ bandpass probed by HST/WFC3 but stronger spectral features at longer wavelengths probed by Spitzer. This led previous authors to puzzling conclusions about the thermal structures and chemical abundances of these planets. Using the SPARC/MITgcm, we investigate how thermal dissociation, ionization, H$^-$ opacity and clouds shape the thermal structures and spectral properties of ultra hot Jupiters with a special focus on WASP-121b. We expand our findings to the whole population of ultra hot Jupiters through analytical quantification of the thermal dissociation and its influence on the strength of spectral features. We predict that most molecules are thermally dissociated and alkalies are ionized in the dayside photospheres of ultra hot Jupiters. This includes H$_{\rm 2}$O, TiO, VO, and H$_{\rm 2}$ but not CO, which has a stronger molecular bond. The vertical molecular gradient created by the dissociation significantly weakens the spectral features from water while the $4.5μm$ CO feature remain unchanged. The water band in the HST/WFC3 bandpass is further weakened by H$^-$ continuum opacity. Molecules are expected to recombine before reaching the limb, leading to order of magnitude variations of the chemical composition and cloud coverage between the limb and the dayside. Overall, molecular dissociation provides a qualitative understanding of the lack of strong spectral feature of water in the $1-2μm$ bandpass observed in most ultra hot Jupiters. Quantitatively, however, our model does not provide a satisfactory match to the WASP-121b emission spectrum. Together with WASP-33b and Kepler-33Ab, they seem the outliers among the population of ultra hot Jupiters in need of a more thorough understanding. △ Less

Submitted 6 August, 2018; v1 submitted 30 April, 2018; originally announced May 2018.

Comments: Accepted in A&A

Journal ref: A&A 617, A110 (2018)

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

Characterizing Earth Analogs in Reflected Light: Atmospheric Retrieval Studies for Future Space Telescopes

Authors: Y. Katherina Feng, Tyler D. Robinson, Jonathan J. Fortney, Roxana E. Lupu, Mark S. Marley, Nikole K. Lewis, Bruce Macintosh, Michael R. Line

Abstract: Space-based high contrast imaging mission concepts for studying rocky exoplanets in reflected light are currently under community study. We develop an inverse modeling framework to estimate the science return of such missions given different instrument design considerations. By combining an exoplanet albedo model, an instrument noise model, and an ensemble Markov chain Monte Carlo sampler, we expl… ▽ More Space-based high contrast imaging mission concepts for studying rocky exoplanets in reflected light are currently under community study. We develop an inverse modeling framework to estimate the science return of such missions given different instrument design considerations. By combining an exoplanet albedo model, an instrument noise model, and an ensemble Markov chain Monte Carlo sampler, we explore retrievals of atmospheric and planetary properties for Earth twins as a function of signal-to-noise ratio (SNR) and resolution ($R$). Our forward model includes Rayleigh scattering, single-layer water clouds with patchy coverage, and pressure-dependent absorption due to water vapor, oxygen, and ozone. We simulate data at $R = 70$ and $R = 140$ from 0.4-1.0 $μ$m with SNR $ = 5, 10, 15, 20$ at 550 nm (i.e., for HabEx/LUVOIR-type instruments). At these same SNR, we simulate data for WFIRST paired with a starshade, which includes two photometric points between 0.48-0.6 $μ$m and $R = 50$ spectroscopy from 0.6-0.97 $μ$m. Given our noise model for WFIRST-type detectors, we find that weak detections of water vapor, ozone, and oxygen can be achieved with observations with at least $R = 70$ / SNR$\ = 15$, or $R = 140$ / SNR$\ = 10$ for improved detections. Meaningful constraints are only achieved with $R = 140$ / SNR$\ = 20$ data. The WFIRST data offer limited diagnostic information, needing at least SNR = 20 to weakly detect gases. Most scenarios place limits on planetary radius, but cannot constrain surface gravity and, thus, planetary mass. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Comments: Resubmitted to AAS Journals after incorporating reviewer feedback. 26 pages, 18 figure, 9 tables

Direct detection and characterization of M-dwarf planets using light echoes

Authors: William B. Sparks, Richard L. White, Roxana E. Lupu, Holland C. Ford

Abstract: Exoplanets orbiting M dwarf stars are a prime target in the search for life in the Universe. M dwarf stars are active, with powerful flares that could adversely impact prospects for life, though there are counter-arguments. Here, we turn flaring to advantage and describe ways in which it can be used to enhance the detectability of planets, in the absence of transits or a coronagraph, significantly… ▽ More Exoplanets orbiting M dwarf stars are a prime target in the search for life in the Universe. M dwarf stars are active, with powerful flares that could adversely impact prospects for life, though there are counter-arguments. Here, we turn flaring to advantage and describe ways in which it can be used to enhance the detectability of planets, in the absence of transits or a coronagraph, significantly expanding the accessible discovery and characterization space. Flares produce brief bursts of intense luminosity, after which the star dims. Due to the light travel time between the star and planet, the planet receives the high intensity pulse, which it re-emits through scattering (a light echo) or intrinsic emission when the star is much fainter, thereby increasing the planet's detectability. The planet's light echo emission can potentially be discriminated from that of the host star by means of a time delay, Doppler shift, spatial shift, and polarization, each of which can improve the contrast of the planet to the star. Scattered light can reveal the albedo spectrum of the planet to within a size scale factor, and is likely to be polarized. Intrinsic emission mechanisms include fluorescent pumping of multiple molecular hydrogen and neutral oxygen lines by intense LyAlpha and LyBeta flare emission, recombination radiation of ionized and photodissociated species, and atmospheric processes such as terrestrial upper atmosphere airglow and near infrared hydroxyl emission. We discuss the feasibility of detecting light echoes and find that under favorable circumstances, echo detection is possible. △ Less

Submitted 3 January, 2018; originally announced January 2018.

An ultrahot gas-giant exoplanet with a stratosphere

Authors: Thomas M. Evans, David K. Sing, Tiffany Kataria, Jayesh Goyal, Nikolay Nikolov, Hannah R. Wakeford, Drake Deming, Mark S. Marley, David S. Amundsen, Gilda E. Ballester, Joanna K. Barstow, Lotfi Ben-Jaffel, Vincent Bourrier, Lars A. Buchhave, Ofer Cohen, David Ehrenreich, Antonio Garcia Munoz, Gregory W. Henry, Heather Knutson, Panayotis Lavvas, Alain Lecavelier des Etangs, Nikole K. Lewis, Mercedes Lopez-Morales, Avi M. Mandell, Jorge Sanz-Forcada , et al. (2 additional authors not shown)

Abstract: Infrared radiation emitted from a planet contains information about the chemical composition and vertical temperature profile of its atmosphere. If upper layers are cooler than lower layers, molecular gases will produce absorption features in the planetary thermal spectrum. Conversely, if there is a stratosphere - where temperature increases with altitude - these molecular features will be observe… ▽ More Infrared radiation emitted from a planet contains information about the chemical composition and vertical temperature profile of its atmosphere. If upper layers are cooler than lower layers, molecular gases will produce absorption features in the planetary thermal spectrum. Conversely, if there is a stratosphere - where temperature increases with altitude - these molecular features will be observed in emission. It has been suggested that stratospheres could form in highly irradiated exoplanets, but the extent to which this occurs is unresolved both theoretically and observationally. A previous claim for the presence of a stratosphere remains open to question, owing to the challenges posed by the highly variable host star and the low spectral resolution of the measurements. Here we report a near-infrared thermal spectrum for the ultrahot gas giant WASP-121b, which has an equilibrium temperature of approximately 2,500 kelvin. Water is resolved in emission, providing a detection of an exoplanet stratosphere at 5-sigma confidence. These observations imply that a substantial fraction of incident stellar radiation is retained at high altitudes in the atmosphere, possibly by absorbing chemical species such as gaseous vanadium oxide and titanium oxide. △ Less

Submitted 3 August, 2017; originally announced August 2017.

Comments: This is the authors' version of the manuscript. 23 pages, 9 figures, 4 tables

Journal ref: Nature, 548, 58-61 (3 August 2017)

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

High Efficiency UV/Optical/NIR Detectors for Large Aperture Telescopes and UV Explorer Missions: Development of and Field Observations with Delta-doped Arrays

Authors: Shouleh Nikzad, April D. Jewell, Michael E. Hoenk, Todd Jones, John Hennessy, Tim Goodsall, Alexander Carver, Charles Shapiro, Samuel R. Cheng, Erika Hamden, Gillian Kyne, D. Christopher Martin, David Schiminovich, Paul Scowen, Kevin France, Stephan McCandliss, Roxana E. Lupu

Abstract: A number of exciting concepts are under development for Flagship, Probe class, Explorer class, and Suborbital class NASA missions in the ultraviolet/optical spectral ranges. These missions will depend on high performance silicon detector arrays being delivered affordably and in high numbers. In a focused effort we have advanced delta-doping technology to high throughput and high yield wafer-scale… ▽ More A number of exciting concepts are under development for Flagship, Probe class, Explorer class, and Suborbital class NASA missions in the ultraviolet/optical spectral ranges. These missions will depend on high performance silicon detector arrays being delivered affordably and in high numbers. In a focused effort we have advanced delta-doping technology to high throughput and high yield wafer-scale processing, encompassing a multitude of state-of-the-art silicon-based detector formats and designs. As part of this technology advancement and in preparation for upcoming missions, we have embarked on a number of field observations, instrument integrations, and independent evaluations of delta-doped arrays. In this paper, we present recent data and innovations from the Advanced Detectors and Systems program at JPL, including two-dimensional doping technology; our end-to-end post-fabrication processing of high performance UV/Optical/NIR arrays; and advanced coatings for detectors and optical elements. Additionally, we present examples of past, in-progress, and planned observations and deployments of delta-doped arrays. △ Less

Submitted 12 January, 2017; v1 submitted 14 December, 2016; originally announced December 2016.

Comments: submitted to JATIS

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

Atmospheric Retrieval for Direct Imaging Spectroscopy of Gas Giants In Reflected Light II: Orbital Phase and Planetary Radius

Authors: Michael Nayak, Roxana Lupu, Mark Marley, Jonathan Fortney, Tyler Robinson, Nikole Lewis

Abstract: Future space-based telescopes, such as the Wide-Field Infrared Survey Telescope (WFIRST), will observe the reflected-light spectra of directly imaged extrasolar planets. Interpretation of such data presents a number of novel challenges, including accounting for unknown planet radius and uncertain stellar illumination phase angle. Here we report on our continued development of Markov Chain Monte Ca… ▽ More Future space-based telescopes, such as the Wide-Field Infrared Survey Telescope (WFIRST), will observe the reflected-light spectra of directly imaged extrasolar planets. Interpretation of such data presents a number of novel challenges, including accounting for unknown planet radius and uncertain stellar illumination phase angle. Here we report on our continued development of Markov Chain Monte Carlo retrieval methods for addressing these issues in the interpretation of such data. Specifically we explore how the unknown planet radius and potentially poorly known observer-planet-star phase angle impacts retrievals of parameters of interest such as atmospheric methane abundance, cloud properties and surface gravity. As expected, the uncertainty in retrieved values is a strong function of signal-to-noise ratio (SNR) of the observed spectra, particularly for low metallicity atmospheres, which lack deep absorption signatures. Meaningful results may only be possible above certain SNR thresholds; for cases across a metallicity range of 1-50 times solar, we find that only an SNR of 20 systematically reproduces close to the correct methane abundance at all phase angles. However, even in cases where the phase angle is poorly known we find that the planet radius can be constrained to within a factor of two. We find that uncertainty in planet radius decreases at phase angles past quadrature, as the highly forward scattering nature of the atmosphere at these geometries limits the possible volume of phase space that relevant parameters can occupy. Finally, we present an estimation of possible improvement that can result from combining retrievals against observations at multiple phase angles. △ Less

Submitted 1 December, 2016; originally announced December 2016.

Comments: Accepted to PASP

Exo-Transmit: An Open-Source Code for Calculating Transmission Spectra for Exoplanet Atmospheres of Varied Composition

Authors: Eliza M. -R. Kempton, Roxana E. Lupu, Albert Owusu-Asare, Patrick Slough, Bryson Cale

Abstract: We present Exo-Transmit, a software package to calculate exoplanet transmission spectra for planets of varied composition. The code is designed to generate spectra of planets with a wide range of atmospheric composition, temperature, surface gravity, and size, and is therefore applicable to exoplanets ranging in mass and size from hot Jupiters down to rocky super-Earths. Spectra can be generated w… ▽ More We present Exo-Transmit, a software package to calculate exoplanet transmission spectra for planets of varied composition. The code is designed to generate spectra of planets with a wide range of atmospheric composition, temperature, surface gravity, and size, and is therefore applicable to exoplanets ranging in mass and size from hot Jupiters down to rocky super-Earths. Spectra can be generated with or without clouds or hazes with options to (1) include an optically thick cloud deck at a user-specified atmospheric pressure or (2) to augment the nominal Rayleigh scattering by a user-specified factor. The Exo-Transmit code is written in C and is extremely easy to use. Typically the user will only need to edit parameters in a single user input file in order to run the code for a planet of their choosing. Exo-Transmit is available publicly on Github with open-source licensing at https://github.com/elizakempton/Exo_Transmit . △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: 8 pages, 7 figures, 1 table, submitted to PASP

Forward and Inverse Modeling of the Emission and Transmission Spectrum of GJ 436b: Investigating Metal Enrichment, Tidal Heating, and Clouds

Authors: Caroline V. Morley, Heather Knutson, Michael Line, Jonathan J. Fortney, Daniel Thorngren, Mark S. Marley, Dillon Teal, Roxana Lupu

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

Submitted 24 October, 2016; originally announced October 2016.

Comments: 15 pages, 18 figures. Revised for publication in ApJ

The First Spectrum of the Coldest Brown Dwarf

Authors: Andrew Skemer, Caroline Morley, Katelyn Allers, Thomas Geballe, Mark Marley, Jonathan Fortney, Jacqueline Faherty, Gordon Bjoraker, Roxana Lupu

Abstract: The recently discovered brown dwarf WISE 0855 presents our first opportunity to directly study an object outside the Solar System that is nearly as cold as our own gas giant planets. However the traditional methodology for characterizing brown dwarfs---near infrared spectroscopy---is not currently feasible as WISE 0855 is too cold and faint. To characterize this frozen extrasolar world we obtained… ▽ More The recently discovered brown dwarf WISE 0855 presents our first opportunity to directly study an object outside the Solar System that is nearly as cold as our own gas giant planets. However the traditional methodology for characterizing brown dwarfs---near infrared spectroscopy---is not currently feasible as WISE 0855 is too cold and faint. To characterize this frozen extrasolar world we obtained a 4.5-5.2 $μ$m spectrum, the same bandpass long used to study Jupiter's deep thermal emission. Our spectrum reveals the presence of atmospheric water vapor and clouds, with an absorption profile that is strikingly similar to Jupiter. The spectrum is high enough quality to allow the investigation of dynamical and chemical processes that have long been studied in Jupiter's atmosphere, but now on an extrasolar world. △ Less

Submitted 16 May, 2016; originally announced May 2016.

Comments: submitted to ApJL

The Hunt for Planet Nine: Atmosphere, Spectra, Evolution, and Detectability

Authors: Jonathan J. Fortney, Mark S. Marley, Gregory Laughlin, Nadine Nettelmann, Caroline V. Morley, Roxana E. Lupu, Channon Visscher, Pavle Jeremic, Wade G. Khadder, Mason Hargrave

Abstract: We investigate the physical characteristics of the Solar System's proposed Planet Nine using modeling tools with a heritage in studying Uranus and Neptune. For a range of plausible masses and interior structures, we find upper limits on the intrinsic Teff, from ~35-50 K for masses of 5-20 M_Earth, and we also explore lower Teff values. Possible planetary radii could readily span from 3 to 6 R_Eart… ▽ More We investigate the physical characteristics of the Solar System's proposed Planet Nine using modeling tools with a heritage in studying Uranus and Neptune. For a range of plausible masses and interior structures, we find upper limits on the intrinsic Teff, from ~35-50 K for masses of 5-20 M_Earth, and we also explore lower Teff values. Possible planetary radii could readily span from 3 to 6 R_Earth depending on the mass fraction of any H/He envelope. Given its cold temperature, the planet encounters significant methane condensation, which dramatically alters the atmosphere away from simple Neptune-like expectations. We find the atmosphere is strongly depleted in molecular absorption at visible wavelengths, suggesting a Rayleigh scattering atmosphere with a high geometric albedo approaching 0.75. We highlight two diagnostics for the atmosphere's temperature structure, the first being the value of the methane mixing ratio above the methane cloud. The second is the wavelength at which cloud scattering can be seen, which yields the cloud-top pressure. Surface reflection may be seen if the atmosphere is thin. Due to collision-induced opacity of H2 in the infrared, the planet would be extremely blue (instead of red) in the shortest wavelength WISE colors if methane is depleted, and would, in some cases, exist on the verge of detectability by WISE. For a range of models, thermal fluxes from ~3-5 microns are ~20 orders of magnitude larger than blackbody expectations. We report a search of the AllWISE Source Catalog for Planet Nine, but find no detection. △ Less

Submitted 6 June, 2016; v1 submitted 25 April, 2016; originally announced April 2016.

Comments: Accepted to ApJ Letters, with minor changes from the submitted version

Developing Atmospheric Retrieval Methods for Direct Imaging Spectroscopy of Gas Giants in Reflected Light I: Methane Abundances and Basic Cloud Properties

Authors: Roxana E. Lupu, Mark S. Marley, Nikole Lewis, Michael Line, Wesley A. Traub, Kevin Zahnle

Abstract: Upcoming space-based coronagraphic instruments in the next decade will perform reflected light spectroscopy and photometry of cool, directly imaged extrasolar giant planets. We are developing a new atmospheric retrieval methodology to help assess the science return and inform the instrument design for such future missions, and ultimately interpret the resulting observations. Our retrieval techniqu… ▽ More Upcoming space-based coronagraphic instruments in the next decade will perform reflected light spectroscopy and photometry of cool, directly imaged extrasolar giant planets. We are developing a new atmospheric retrieval methodology to help assess the science return and inform the instrument design for such future missions, and ultimately interpret the resulting observations. Our retrieval technique employs a geometric albedo model coupled with both a Markov chain Monte Carlo Ensemble Sampler (emcee) and a multimodal nested sampling algorithm (MultiNest) to map the posterior distribution. This combination makes the global evidence calculation more robust for any given model, and highlights possible discrepancies in the likelihood maps. As a proof-of-concept, our current atmospheric model contains 1 or 2 cloud layers, methane as a major absorber, and a H$_2$-He background gas. This 6-to-9 parameter model is appropriate for Jupiter-like planets and can be easily expanded in the future. In addition to deriving the marginal likelihood distribution and confidence intervals for the model parameters, we perform model selection to determine the significance of methane and cloud detection as a function of expected signal-to-noise in the presence of spectral noise correlations. After internal validation, the method is applied to realistic spectra of Jupiter, Saturn, and HD 99492 c, a model observing target. We find that the presence or absence of clouds and methane can be determined with high confidence, while parameter uncertainties are model-dependent and correlated. Such general methods will also be applicable to the interpretation of direct imaging spectra of cloudy terrestrial planets. △ Less

Submitted 16 September, 2016; v1 submitted 18 April, 2016; originally announced April 2016.

Comments: 37 pages, 30 figures, accepted for publication in AJ

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

The Tethered Moon

Authors: Kevin J. Zahnle, Roxana Lupu, Anthony Dobrovolskis, Norman H. Sleep

Abstract: We address the thermal history of the Earth after the Moon-forming impact, taking tidal heating and thermal blanketing by the atmosphere into account. The atmosphere sets an upper bound of ~100 W/m^2 on how quickly the Earth can cool. The liquid magma ocean cools over 2-10 Myrs, with longer times corresponding to high angular-momentum events. Tidal heating is focused mostly in mantle materials tha… ▽ More We address the thermal history of the Earth after the Moon-forming impact, taking tidal heating and thermal blanketing by the atmosphere into account. The atmosphere sets an upper bound of ~100 W/m^2 on how quickly the Earth can cool. The liquid magma ocean cools over 2-10 Myrs, with longer times corresponding to high angular-momentum events. Tidal heating is focused mostly in mantle materials that are just beginning to freeze. The atmosphere's control over cooling sets up a negative feedback between viscosity-dependent tidal heating and temperature-dependent viscosity of the magma ocean. While the feedback holds, evolution of the Moon's orbit is limited by the modest radiative cooling rate of Earth's atmosphere. Orbital evolution is orders of magnitude slower than in conventional constant Q models, which promotes capture by resonances. The evection resonance is encountered early, when the Earth is molten. Capture by the evection resonance appears certain but unlikely to generate much eccentricity because it is encountered early when the Earth is molten and Q_Earth >> Q_Moon. Tidal dissipation in the Earth becomes more efficient (Q_Earth << Q_Moon) later when the Moon is between ~20 R_Earth and ~40 R_Earth. If lunar eccentricity grew great, this was when it did so, perhaps setting the table for some other process to leave its mark on the inclination of the Moon. △ Less

Submitted 6 August, 2015; originally announced August 2015.

Comments: 34 pages, 10 figures, preprint format

Journal ref: Earth and Planetary Science Letters 427, 74-82 (2015)

A Quick Study of the Characterization of Radial Velocity Giant Planets in Reflected Light by Forward and Inverse Modeling

Authors: Mark Marley, Roxana Lupu, Nikole Lewis, Michael Line, Caroline Morley, Jonathan Fortney

Abstract: We explored two aspects of the problem of characterizing cool extrasolar giant planets in scattered optical light with a space based coronagraph. First, for a number of the known radial velocity (RV) giants we computed traditional forward models of their atmospheric structure and clouds, given various input assumptions, and computed model albedo spectra. Such models have been computed before, but… ▽ More We explored two aspects of the problem of characterizing cool extrasolar giant planets in scattered optical light with a space based coronagraph. First, for a number of the known radial velocity (RV) giants we computed traditional forward models of their atmospheric structure and clouds, given various input assumptions, and computed model albedo spectra. Such models have been computed before, but mostly for generic planets. Our new models demonstrate that there is likely interesting spectral diversity among those planets that are most favorable for direct detection. Second, we applied a powerful Markov Chain Monte Carlo (MCMC) retrieval technique to synthetic noisy data of cool giants to better understand how well various atmospheric parameters--particularly molecular abundances and cloud properties--could be constrained. This is the first time such techniques have been applied to this problem. The process is time consuming, so only a dozen or so cases could be completed in the limited time available. Nevertheless the results clearly show that even at S/N ~ 5, scientifically interesting and valuable conclusions can be drawn about the properties of giant planet atmospheres from noisy spectra. We find that atmospheric abundances are best constrained when the planet gravity is bounded. Thus direct imaging observations of known radial velocity planets are extremely valuable as limits on a target planet's gravity can be obtained from astrometric imaging and reflectivity and mass-radius relationship arguments (for mass and radius, respectively). Further retrieval studies are clearly warranted and would be valuable to help guide instrument design decisions. △ Less

Submitted 29 December, 2014; originally announced December 2014.

Comments: Report of a quick study of science return from direct-imaging exoplanet missions, commissioned by the NASA Exoplanet Exploration Program on behalf of the WFIRST/AFTA Science Definition Team and the Exo-S and Exo-C Science and Technology Definition Teams

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

Herschel-ATLAS and ALMA: HATLAS J142935.3-002836, a lensed major merger at redshift 1.027

Authors: Hugo Messias, Simon Dye, Neil Nagar, Gustavo Orellana, R. Shane Bussmann, Jae Calanog, Helmut Dannerbauer, Hai Fu, Edo Ibar, Andrew Inohara, R. J. Ivison, Mattia Negrello, Dominik A. Riechers, Yun-Kyeong Sheen, Simon Amber, Mark Birkinshaw, Nathan Bourne, Dave L. Clements, Asantha Cooray, Gianfranco De Zotti, Ricardo Demarco, Loretta Dunne, Stephen Eales, Simone Fleuren, Roxana E. Lupu , et al. (7 additional authors not shown)

Abstract: [Abridged] Aims: This work focuses on one lensed system, HATLAS J142935.3-002836 (H1429-0028), selected in the Herschel-ATLAS field. Gathering a rich, multi-wavelength dataset, we aim to confirm the lensing hypothesis and model the background source's morphology and dynamics, as well as to provide a full physical characterisation. Methods: Multi-wavelength high-resolution data is utilised to asses… ▽ More [Abridged] Aims: This work focuses on one lensed system, HATLAS J142935.3-002836 (H1429-0028), selected in the Herschel-ATLAS field. Gathering a rich, multi-wavelength dataset, we aim to confirm the lensing hypothesis and model the background source's morphology and dynamics, as well as to provide a full physical characterisation. Methods: Multi-wavelength high-resolution data is utilised to assess the nature of the system. A lensing-analysis algorithm which simultaneously fits different wavebands is adopted to characterise the lens. The background galaxy dynamical information is studied by reconstructing the 3-D source-plane of the ALMA CO(J:4-3) transition. Near-IR imaging from HST and Keck-AO allows to constrain rest-frame optical photometry independently for the foreground and background systems. Physical parameters (such as stellar and dust masses) are estimated via modelling of the spectral energy distribution taking into account source blending, foreground obscuration, and differential magnification. Results: The system comprises a foreground edge-on disk galaxy (at z_sp=0.218) with an almost complete Einstein ring around it. The background source (at z_sp=1.027) is magnified by a factor of ~8-10 depending on wavelength. It is comprised of two components and a tens of kpc long tidal tail resembling the Antennae merger. As a whole, the system is a massive stellar system (1.32[-0.41,+0.63] x1E11 Mo) forming stars at a rate of 394+-90 Mo/yr, and has a significant gas reservoir M_ISM = 4.6+-1.7 x1E10 Mo. Its depletion time due to star formation alone is thus expected to be tau_SF=M_ISM/SFR=117+-51 Myr. The dynamical mass of one of the components is estimated to be 5.8+-1.7 x1E10 Mo, and, together with the photometric total mass estimate, it implies that H1429-0028 is a major merger system (1:2.8[-1.5,+1.8]). △ Less

Submitted 18 June, 2014; originally announced June 2014.

Comments: submitted to A&A, 20 pages, 12 figures

Journal ref: A&A 568, A92 (2014)

Spectral Variability from the Patchy Atmospheres of T and Y Dwarfs

Authors: Caroline V. Morley, Mark S. Marley, Jonathan J. Fortney, Roxana Lupu

Abstract: Brown dwarfs of a variety of spectral types have been observed to be photometrically variable. Previous studies have focused on objects at the L/T transition, where the iron and silicate clouds in L dwarfs break up or dissipate. However, objects outside of this transitional effective temperature regime also exhibit variability. Here, we present models for mid-late T dwarfs and Y dwarfs. We present… ▽ More Brown dwarfs of a variety of spectral types have been observed to be photometrically variable. Previous studies have focused on objects at the L/T transition, where the iron and silicate clouds in L dwarfs break up or dissipate. However, objects outside of this transitional effective temperature regime also exhibit variability. Here, we present models for mid-late T dwarfs and Y dwarfs. We present models that include patchy salt and sulfide clouds as well as water clouds for the Y dwarfs. We find that for objects over 375 K, patchy cloud opacity would generate the largest amplitude variability within near-infrared spectral windows. For objects under 375 K, water clouds also become important and generate larger amplitude variability in the mid-infrared. We also present models in which we perturb the temperature structure at different pressure levels of the atmosphere to simulate hot spots. These models show the most variability in the absorption features between spectral windows. The variability is strongest at wavelengths that probe pressure levels at which the heating is the strongest. The most illustrative types of observations for understanding the physical processes underlying brown dwarf variability are simultaneous, multi-wavelength observations that probe both inside and outside of molecular absorption features. △ Less

Submitted 3 June, 2014; originally announced June 2014.

Comments: 6 pages, 5 figures, Accepted for publication in ApJ Letters