Search

Exoplanets are abundant in our galaxy and yet characterizing them remains a technical challenge. Solar System planets provide an opportunity to test the practical limitations of exoplanet observations with high signal-to-noise data that we cannot access for exoplanets. However, data on Solar System planets differ from exoplanets in that Solar System planets are spatially resolved while exoplanets are unresolved point-sources. We present a novel instrument designed to observe Solar System planets as though they are exoplanets, the Planet as Exoplanet Analog Spectrograph (PEAS). PEAS consists of a dedicated 0.5-m telescope and off-the-shelf optics, located at Lick Observatory. PEAS uses an integrating sphere to disk-integrate light from the Solar System planets, producing spatially mixed light more similar to the spectra we can obtain from exoplanets. This paper describes the general system design and early results of the PEAS instrument., Comment: SPIE Astronomical Telescopes & Instrumentation 2020 Proceedings (11447-153), 9 pages, 9 figures

We present a transmission spectrum for the warm (500-600K) sub-Neptune HD3167c obtained using the Hubble Space Telescope Wide Field Camera 3 infrared spectrograph. We combine these data, which span the 1.125-1.643 micron wavelength range, with broadband transit measurements made using Kepler/K2 (0.6-0.9 micron) and Spitzer/IRAC (4-5 micron). We find evidence for absorption by at least one of H2O, HCN, CO2, and CH4 (Bayes factor 7.4; 2.5-sigma significance), although the data precision does not allow us to unambiguously discriminate between these molecules. The transmission spectrum rules out cloud-free hydrogen-dominated atmospheres with metallicities <100x solar at >5.8-sigma confidence. In contrast, good agreement with the data is obtained for cloud-free models assuming metallicities >700x solar. However, for retrieval analyses that include the effect of clouds, a much broader range of metallicities (including subsolar) is consistent with the data, due to the degeneracy with cloud-top pressure. Self-consistent chemistry models that account for photochemistry and vertical mixing are presented for the atmosphere of HD3167c. The predictions of these models are broadly consistent with our abundance constraints, although this is primarily due to the large uncertainties on the latter. Interior structure models suggest the core mass fraction is >40%, independent of a rock or water core composition, and independent of atmospheric envelope metallicity up to 1000x solar. We also report abundance measurements for fifteen elements in the host star, showing that it has a very nearly solar composition., Comment: Accepted for publication in the Astronomical Journal

The orbital architecture of the Solar System is thought to have been sculpted by a dynamical instability among the giant planets. During the instability a primordial outer disk of planetesimals was destabilized and ended up on planet-crossing orbits. Most planetesimals were ejected into interstellar space but a fraction were trapped on stable orbits in the Kuiper belt and Oort cloud. We use a suite of N-body simulations to map out the diversity of planetesimals' dynamical pathways. We focus on two processes: tidal disruption from very close encounters with a giant planet, and loss of surface volatiles from repeated passages close to the Sun. We show that the rate of tidal disruption is more than a factor of two higher for ejected planetesimals than for surviving objects in the Kuiper belt or Oort cloud. Ejected planetesimals are preferentially disrupted by Jupiter and surviving ones by Neptune. Given that the gas giants contracted significantly as they cooled but the ice giants did not, taking into account the thermal evolution of the giant planets decreases the disruption rate of ejected planetesimals. The frequency of volatile loss and extinction is far higher for ejected planetesimals than for surviving ones and is not affected by the giant planets' contraction. Even if all interstellar objects were ejected from Solar System-like systems, our analysis suggests that their physical properties should be more diverse than those of Solar System small bodies as a result of their divergent dynamical histories. This is consistent with the characteristics of the two currently-known interstellar objects., Comment: ApJ Letters, in press. 6 pages, 4 figures

Hot Jupiters have been predicted to have a strong day/night temperature contrast and a hot spot shifted eastward of the substellar point. This was confirmed by numerous phase curve observations probing the longitudinal brightness variation of the atmosphere. Global circulation models, however, systematically underestimate the phase curve amplitude and overestimate the shift of its maximum. We use a global circulation model including non-grey radiative transfer and realistic gas and cloud opacities to systematically investigate how the atmospheric circulation of hot Jupiters varies with equilibrium temperature from 1000 to 2200K. We show that the heat transport is very efficient for cloudless planets cooler than 1600K and becomes less efficient at higher temperatures. When nightside clouds are present, the day-to-night heat transport becomes extremely inefficient, leading to a good match to the observed low nightside temperatures. The constancy of this low temperature is, however, due to the strong dependence of the radiative timescale with temperature. We further show that nightside clouds increase the phase curve amplitude and decreases the phase curve offset at the same time. This change is very sensitive to the cloud chemical composition and particle size, meaning that the diversity in observed phase curves can be explained by a diversity of nightside cloud properties. Finally, we show that phase curve parameters do not necessarily track the day/night contrast nor the shift of the hot spot on isobars, and propose solutions to to recover the true hot-spot shift and day/night contrast., Comment: 32 pages, 21 figures, 4 years in the making. Accepted in MNRAS. See also the complementary study by Roman et al. today on ArXiv

We have developed an open-source pipeline for the analysis of \textit{Spitzer}/IRAC channel 1 and 2 time-series photometry, incorporating some of the most popular decorrelation methods. We applied this pipeline to new phase curve observations of ultra-hot Jupiters MASCARA-1b and KELT-16b, and we performed the first comprehensive reanalysis of 15 phase curves. We find that MASCARA-1b and KELT-16b have phase offsets of $6^{+11}_{-11}~^{\circ}$W and $38^{+16}_{-15}~^{\circ}$W, dayside temperatures of $2952^{+100}_{-97}$ K and $3070^{+160}_{-150}$ K, and nightside temperatures of $1300^{+340}_{-340}$ K and $1900^{+430}_{-440}$ K, respectively. We confirm a strong correlation between dayside and irradiation temperatures with a shallower dependency for nightside temperature. We also find evidence that the normalized phase curve amplitude (peak-to-trough divided by eclipse depth) is correlated with stellar effective temperature. In addition, while our different models often retrieve similar parameters, significant differences occasionally arise between them, as well as between our preferred model and the literature values. Nevertheless, our preferred models are consistent with published phase offsets to within $-8\pm21$ degrees ($-1.6\pm3.2$ sigma), and normalized phase curve amplitudes are on average reproduced to within $-0.01\pm0.24$ ($-0.1\pm1.6$ sigma). Finally, we find that BLISS performs best in most cases, but not all; we therefore recommend future analyses consider numerous detector models to ensure an optimal fit and to assess model dependencies., Comment: Resubmitted to MNRAS after revision on February 25, 2021, originally submitted to MNRAS on October 1, 2020. 20 pages, 12 figures, and 5 tables. 3 pages of appendix, 16 pages of supplementary information

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., Comment: Accepted to AJ. No additional significant changes

We present a new optical (400-950nm) transmission spectrum of the hot Jupiter WASP-31b (M=0.48 MJ; R= 1.54 RJ; P=3.41 days), obtained by combining four transits observations. These transits were observed with IMACS on the Magellan Baade Telescope at Las Campanas Observatory as part of the ACCESS project. We investigate the presence of clouds/hazes in the upper atmosphere of this planet as well as the contribution of stellar activity on the observed features. In addition, we search for absorption features of the alkali elements Na I and K I, with particular focus on K I, for which there have been two previously published disagreeing results. Observations with HST/STIS detected K I, whereas ground-based low- and high-resolution observations did not. We use equilibrium and non-equilibrium chemistry retrievals to explore the planetary and stellar parameter space of the system with our optical data combined with existing near-IR observations. Our best-fit model is that with a scattering slope consistent with a Rayleigh slope (alpha=5.3+2.9-3.1), high-altitude clouds at a log cloud top pressure of -3.6+2.7-2.1 bars, and possible muted H2O features. We find that our observations support other ground-based claims of no K I. Clouds are likely why signals like H2O are extremely muted and Na or K cannot be detected. We then juxtapose our Magellan/IMACS transmission spectrum with existing VLT/FORS2, HST/WFC3, HST/STIS, and Spitzer observations to further constrain the optical-to-infrared atmospheric features of the planet. We find that a steeper scattering slope (alpha = 8.3+/-1.5) is anchored by STIS wavelengths blueward of 400 nm and only the original STIS observations show significant potassium signal., Comment: Accepted 14 September 2020 by AJ

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 $\mu$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., Comment: Submitted to AAS Journals, revised to reflect referee comments. Posting of this manuscript on the arXiv was coordinated with S. Ghandi et al

Uranus and Neptune form a distinct class of planets in our solar system. Given this fact, and ubiquity of similar-mass planets in other planetary systems, it is essential to understand their interior structure and composition. However, there are more open questions regarding these planets than answers. In this review we concentrate on the things we do not know about the interiors of Uranus and Neptune with a focus on why the planets may be different, rather than the same. We next summarize the knowledge about the planets' internal structure and evolution. Finally, we identify the topics that should be investigated further on the theoretical front as well as required observations from space missions., Comment: Accepted for publication in Philosophical Transactions of the Royal Society A

The nearby super-Earth 55 Cnc e orbits a bright (V = 5.95 mag) star with a period of ~ 18 hours and a mass of ~ 8 Earth masses. Its atmosphere may be water-rich and have a large scale-height, though attempts to characterize it have yielded ambiguous results. Here we present a sensitive search for water and TiO in its atmosphere at high spectral resolution using the Gemini North telescope and the GRACES spectrograph. We combine observations with previous observations from Subaru and CFHT, improving the constraints on the presence of water vapor. We adopt parametric models with an updated planet radius based on recent measurements, and use a cross-correlation technique to maximize sensitivity. Our results are consistent with atmospheres that are cloudy or contain minimal amounts of water and TiO. Using these parametric models, we rule out a water-rich atmosphere (VMR >= 0.1%) with a mean molecular weight of <= 15 g/mol at a 3 sigma confidence level, improving on the previous limit by a significant margin. For TiO, we rule out a mean molecular weight of <= 5 g/mol with a 3 sigma confidence level for a VMR greater than 10^-8; for a VMR of greater than 10^-7, the limit rises to a mean molecular weight of <= 10 g/mol. We can rule out low mean-molecular-weight chemical equilibrium models both including and excluding TiO/VO at very high confidence levels (> 10 sigma). Overall, our results are consistent with an atmosphere with a high mean molecular weight and/or clouds, or no atmosphere., Comment: 15 pages, 12 figures, 1 table

Spectroscopic phase curves provide unique access to the three-dimensional properties of transiting exoplanet atmospheres. However, a modeling framework must be developed to deliver accurate inferences of atmospheric properties for these complex data sets. Here, we develop an approach to retrieve temperature structures and molecular abundances from phase curve spectra at any orbital phase. In the context of a representative hot Jupiter with a large day-night temperature contrast, we examine the biases in typical one-dimensional (1D) retrievals as a function of orbital phase/geometry, compared to two-dimensional (2D) models that appropriately capture the disk-integrated phase geometry. We guide our intuition by applying our new framework on a simulated HST+Spitzer phase curve data set in which the "truth" is known, followed by an application to the spectroscopic phase curve of the canonical hot Jupiter, WASP-43b. We also demonstrate the retrieval framework on simulated JWST phase curve observations. We apply our new geometric framework to a joint-fit of all spectroscopic phases, assuming longitudinal molecular abundance homogeneity, resulting in an a factor of 2 improvement in abundances precision when compared to individual phase constraints. With a 1D retrieval model on simulated HST+Spitzer data, we find strongly biased molecular abundances for CH$_4$ and CO$_2$ at most orbital phases. With 2D, the day and night profiles retrieved from WASP-43b remain consistent throughout the orbit. JWST retrievals show that a 2D model is strongly favored at all orbital phases. Based on our new 2D retrieval implementation, we provide recommendations on when 1D models are appropriate and when more complex phase geometries involving multiple TP profiles are required to obtain an unbiased view of tidally locked planetary atmospheres., Comment: Accepted at AJ. 23 pages, 30 figures - see figure 2 for quick takeaway

We present a transmission spectrum for the Neptune-size exoplanet HD 106315 c from optical to infrared wavelengths based on transit observations from the Hubble Space Telescope/Wide Field Camera 3, K2, and Spitzer. The spectrum shows tentative evidence for a water absorption feature in the $1.1 - 1.7\mu$m wavelength range with a small amplitude of 30 ppm (corresponding to just $0.8 \pm 0.04$ atmospheric scale heights). Based on an atmospheric retrieval analysis, the presence of water vapor is tentatively favored with a Bayes factor of 1.7 - 2.6 (depending on prior assumptions). The spectrum is most consistent with either enhanced metallicity, high altitude condensates, or both. Cloud-free solar composition atmospheres are ruled out at $>5\sigma$ confidence. We compare the spectrum to grids of cloudy and hazy forward models and find that the spectrum is fit well by models with moderate cloud lofting or haze formation efficiency, over a wide range of metallicities ($1 - 100\times$ solar). We combine the constraints on the envelope composition with an interior structure model and estimate that the core mass fraction is $\gtrsim0.3$. With a bulk composition reminiscent of that of Neptune and an orbital distance of 0.15 AU, HD 106315 c hints that planets may form out of broadly similar material and arrive at vastly different orbits later in their evolution., Comment: Submitted to AAS journals; 19 pages, 12 figures

Aerosols are common in the atmospheres of exoplanets across a wide swath of temperatures, masses, and ages. These aerosols strongly impact observations of transmitted, reflected, and emitted light from exoplanets, obfuscating our understanding of exoplanet thermal structure and composition. Knowing the dominant aerosol composition would facilitate interpretations of exoplanet observations and theoretical understanding of their atmospheres. A variety of compositions have been proposed, including metal oxides and sulphides, iron, chromium, sulphur, and hydrocarbons. However, the relative contributions of these species to exoplanet aerosol opacity is unknown. Here we show that the aerosol composition of giant exoplanets observed in transmission is dominated by silicates and hydrocarbons. By constraining an aerosol microphysics model with trends in giant exoplanet transmission spectra, we find that silicates dominate aerosol opacity above planetary equilibrium temperatures of 950 K due to low nucleation energy barriers and high elemental abundances, while hydrocarbon aerosols dominate below 950 K due to an increase in methane abundance. Our results are robust to variations in planet gravity and atmospheric metallicity within the range of most giant transiting exoplanets. We predict that spectral signatures of condensed silicates in the mid-infrared are most prominent for hot (>1600 K), low-gravity (<10 m s$^{-2}$) objects., Comment: 32 pages, 4 main text figures, 7 supplementary figures, 4 supplementary tables, published in Nature Astronomy

We present an optical-to-infrared transmission spectrum of the inflated sub-Saturn KELT-11b measured with the Transiting Exoplanet Survey Satellite (TESS), the Hubble Space Telescope (HST) Wide Field Camera 3 G141 spectroscopic grism, and the Spitzer Space Telescope (Spitzer) at 3.6 $\mu$m, in addition to a Spitzer 4.5 $\mu$m secondary eclipse. The precise HST transmission spectrum notably reveals a low-amplitude water feature with an unusual shape. Based on free retrieval analyses with varying molecular abundances, we find strong evidence for water absorption. Depending on model assumptions, we also find tentative evidence for other absorbers (HCN, TiO, and AlO). The retrieved water abundance is generally $\lesssim 0.1\times$ solar (0.001--0.7$\times$ solar over a range of model assumptions), several orders of magnitude lower than expected from planet formation models based on the solar system metallicity trend. We also consider chemical equilibrium and self-consistent 1D radiative-convective equilibrium model fits and find they too prefer low metallicities ($[M/H] \lesssim -2$, consistent with the free retrieval results). However, all the retrievals should be interpreted with some caution since they either require additional absorbers that are far out of chemical equilibrium to explain the shape of the spectrum or are simply poor fits to the data. Finally, we find the Spitzer secondary eclipse is indicative of full heat redistribution from KELT-11b's dayside to nightside, assuming a clear dayside. These potentially unusual results for KELT-11b's composition are suggestive of new challenges on the horizon for atmosphere and formation models in the face of increasingly precise measurements of exoplanet spectra., Comment: Accepted to The Astronomical Journal. 31 pages, 20 figures, 7 tables

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., Comment: Accepted to the Astronomical Journal. 24 Pages, 16 Figures, 4 Tables. take care

We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~\mu$m and $0.30 \pm 0.02 \%$ at $4.5~\mu$m, corresponding to dayside brightness temperatures of $1542^{+32}_{-31}$~K and $1557^{+35}_{-36}$~K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are $1117^{+76}_{-71}$~K and $1167^{+69}_{-74}$~K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of $0.12_{-0.16}^{+0.14}$ and a moderate day-night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright-spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with \textit{Spitzer} full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: $11^{\circ}\pm 7^{\circ}$ at $3.6~\mu$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~\mu$m. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters., Comment: 14 pages, 8 figures. Accepted for publication in AJ

The James Webb Space Telescope (JWST) is expected to revolutionize our understanding of Jovian worlds over the coming decade. However, as we push towards characterizing cooler, smaller, "terrestrial-like" planets, dedicated next-generation facilities will be required to tease out the small spectral signatures indicative of biological activity. Here, we evaluate the feasibility of determining atmospheric properties, from near-to-mid-infrared transmission spectra, of transiting temperate terrestrial M-dwarf companions. Specifically, we utilize atmospheric retrievals to explore the trade space between spectral resolution, wavelength coverage, and signal-to-noise on our ability to both detect molecular species and constrain their abundances. We find that increasing spectral resolution beyond R=100 for near-infrared wavelengths, shorter than 5$\mu$m, proves to reduce the degeneracy between spectral features of different molecules and thus greatly benefits the abundance constraints. However, this benefit is greatly diminished beyond 5$\mu$m as any overlap between broad features in the mid-infrared does not deconvolve with higher resolutions. Additionally, our findings revealed that the inclusion of features beyond 11$\mu$m did not meaningfully improve the detection significance nor abundance constraints results. We conclude that an instrument with continuous wavelength coverage from $\sim$2-11$\mu$m, spectral resolution of R$\simeq$50-300, and a 25m$^2$ collecting area, would be capable of detecting H$_2$O, CO$_2$, CH$_4$, O$_3$, and N$_2$O in the atmosphere of an Earth-analog transiting a M-dwarf (mag$_{K}=8.0$) within 50 transits, and obtain better than an order-of-magnitude constraint on each of their abundances., Comment: Accepted to AJ 12/19/2019

The gravity field of a giant planet is typically our best window into its interior structure and composition. Through comparison of a model planet's calculated gravitational potential with the observed potential, inferences can be made about interior quantities, including possible composition and the existence of a core. Necessarily, a host of assumptions go into such calculations, making every inference about a giant planet's structure strongly model dependent. In this work we present a more general picture by setting Saturn's gravity field, as measured during the \emph{Cassini} Grand Finale, as a likelihood function driving a Markov-Chain Monte Carlo exploration of the possible interior density profiles. The result is a posterior distribution of the interior structure that is not tied to assumed composition, thermal state, or material equations of state. Constraints on interior structure derived in this Bayesian framework are necessarily less informative, but are also less biased and more general. These empirical and probabilistic constraints on the density structure are our main data product which we archive for continued analysis. We find that the outer half of Saturn's radius is relatively well constrained, and we interpret our findings as suggesting a significant metal enrichment, in line with atmospheric abundances from remote sensing. As expected, the inner half of Saturn's radius is less well-constrained by gravity, but we generally find solutions that include a significant density enhancement, which can be interpreted as a core, although this core is often lower in density and larger in radial extent than typically found by standard models. This is consistent with a dilute core and/or composition gradients.

We examine the comparative thermal evolution of Jupiter and Saturn applying recent theoretical results for helium's immiscibility in fluid metallic hydrogen. The redistribution of helium in their interiors proceeds very differently for the two planets. We confirm that based on Jupiter's atmospheric helium depletion as observed in situ by the Galileo entry probe, Jupiter's interior helium has differentiated modestly, and we present models reconciling Jupiter's helium depletion, radius, and heat flow at the solar age. Jupiter's recently revised Bond albedo implies a lower intrinsic flux for the planet, accommodating more luminosity from helium differentiation such that mildly superadiabatic interiors can satisfy all constraints. The same phase diagram applied to the less massive Saturn predicts dramatic helium differentiation to the degree that Saturn inevitably forms a helium-rich shell or core, an outcome previously proposed by Stevenson & Salpeter and others. The luminosity from Saturn's helium differentiation is sufficient to extend its cooling time to the solar age, even for adiabatic interiors. This model predicts Saturn's atmospheric helium to be depleted to Y=0.07+/-0.01, corresponding to a He/H_2 mixing ratio 0.036+/-0.006. We also show that neon differentiation may have contributed to both planets' luminosity in the past. These results demonstrate that Jupiter and Saturn's thermal evolution can be explained self-consistently with a single physical model, and emphasize that nontrivial helium distributions should be considered in future models for Saturn's internal structure and dynamo., Comment: Accepted for publication in ApJ

The relationship between the compositions of giant planets and their host stars is of fundamental interest in understanding planet formation. The solar system giant planets are enhanced above solar composition in metals, both in their visible atmospheres and bulk compositions. A key question is whether the metal enrichment of giant exoplanets is correlated with that of their host stars. Thorngren et al. (2016) showed that in cool (Teq < 1000 K) giant exoplanets, the total heavy-element mass increases with total Mp and the heavy element enrichment relative to the parent star decreases with total Mp. In their work, the host star metallicity was derived from literature [Fe/H] measurements. Here we conduct a more detailed and uniform study to determine whether different host star metals (C, O, Mg, Si, Fe, and Ni) correlate with the bulk metallicity of their planets, using correlation tests and Bayesian linear fits. We present new host star abundances of 19 cool giant planet systems, and combine these with existing host star data for a total of 22 cool giant planet systems (24 planets). Surprisingly, we find no clear correlation between stellar metallicity and planetary residual metallicity (the relative amount of metal versus that expected from the planet mass alone), which is in conflict with common predictions from formation models. We also find a potential correlation between residual planet metals and stellar volatile-to-refractory element ratios. These results provide intriguing new relationships between giant planet and host star compositions for future modeling studies of planet formation., Comment: Main results in Figure 3-5. Published in AJ on 2019-11-20. 24 pages (single column style)

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$\,\mu$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\,\mu$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\sigma$ 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$\,\mu$m photometric point, instead of spectroscopy, combined with the green-optical point, fails to produce meaningful information about atmospheric water content., Comment: Accepted to AJ

The Kepler mission revealed a class of planets known as ''super-puffs,'' with masses only a few times larger than Earth's but radii larger than Neptune, giving them very low mean densities. All three of the known planets orbiting the young solar-type star Kepler 51 are super-puffs. The Kepler 51 system thereby provides an opportunity for a comparative study of the structures and atmospheres of this mysterious class of planets, which may provide clues about their formation and evolution. We observed two transits each of Kepler 51b and 51d with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope. Combining new WFC3 transit times with re-analyzed Kepler data and updated stellar parameters, we confirmed that all three planets have densities lower than 0.1 g/cm$^{3}$. We measured the WFC3 transmission spectra to be featureless between 1.15 and 1.63 $\mu$m, ruling out any variations greater than 0.6 scale heights (assuming a H/He dominated atmosphere), thus showing no significant water absorption features. We interpreted the flat spectra as the result of a high-altitude aerosol layer (pressure $<$3 mbar) on each planet. Adding this new result to the collection of flat spectra that have been observed for other sub-Neptune planets, we find support for one of the two hypotheses introduced by Crossfield and Kreidberg (2017), that planets with cooler equilibrium temperatures have more high-altitude aerosols. We strongly disfavor their other hypothesis that the H/He mass fraction drives the appearance of large amplitude transmission features., Comment: 29 pages, 16 figures, Re-submitted to AJ

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $\mu$m. KELT-9b is the hottest known transiting planet, with a 4.5-$\mu$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5\sigma$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.

Two of TESS's major science goals are to measure masses for 50 planets smaller than 4 Earth radii and to discover high-quality targets for atmospheric characterization efforts. It is important that these two goals are linked by quantifying what precision of mass constraint is required to yield robust atmospheric properties of planets. Here, we address this by conducting retrievals on simulated JWST transmission spectra under various assumptions for the degree of uncertainty in the planet's mass for a representative population of seven planets ranging from terrestrials to warm Neptunes to hot Jupiters. Only for the cloud-free, low metallicity gas giants are we able to infer exoplanet mass from transmission spectroscopy alone, to ~10% accuracy. For low metallicity cases (<4xSolar) we are able to accurately constrain atmospheric properties without prior knowledge of the planet's mass. For all other cases (including terrestrial-like planets), atmospheric properties can only be inferred with a mass precision of better than ~50%. At this level, though, the widths of the posterior distributions of the atmospheric properties are dominated by the uncertainties in mass. With a precision of ~20%, the widths of the posterior distributions are dominated by the spectroscopic data quality. Therefore, as a rule-of-thumb, we recommend: a ~50% mass precision for initial atmospheric characterization and a ~20% mass precision for more detailed atmospheric analyses., Comment: Accepted ApJL, 9 pages, 4 figures, 1 table

Probabilistic relationships between the mass and radius of planets are useful for a variety of purposes, including estimating the yields of planet discovery efforts and the radii of discovered planets given their masses. Previous work on giant planets often do not make the important distinction between cool and (inflated) hot giants. This can lead to needlessly large uncertainties and overestimated radii for the cooler objects. Using hierarchical bayesian modeling and the AIC model selection criterion, we identify an empirical mass-radius relationship for cool giant planets above 15 Earth masses. It takes the form of a quadratic function in log mass, with a power-law dispersion. Parametric uncertainties from MCMC are provided, and we compare our results with the outputs of interior structure models and the radii of the solar system giant planets, which fit comfortably within our relation., Comment: Two pages, one figure

Results from the Kepler mission indicate that the occurrence rate of small planets ($<3$ $R_\oplus$) in the habitable zone of nearby low-mass stars may be as high as 80%. Despite this abundance, probing the conditions and atmospheric properties on any habitable-zone planet is extremely difficult and has remained elusive to date. Here, we report the detection of water vapor and the likely presence of liquid and icy water clouds in the atmosphere of the $2.6$ $R_\oplus$ habitable-zone planet K2-18b. The simultaneous detection of water vapor and clouds in the mid-atmosphere of K2-18b is particularly intriguing because K2-18b receives virtually the same amount of total insolation from its host star ($1368_{-107}^{+114}$ W m$^{-2}$) as the Earth receives from the Sun (1361 W m$^{-2}$), resulting in the right conditions for water vapor to condense and explain the detected clouds. In this study, we observed nine transits of K2-18b using HST/WFC3 in order to achieve the necessary sensitivity to detect the water vapor, and we supplement this data set with Spitzer and K2 observations to obtain a broader wavelength coverage. While the thick hydrogen-dominated envelope we detect on K2-18b means that the planet is not a true Earth analog, our observations demonstrate that low-mass habitable-zone planets with the right conditions for liquid water are accessible with state-of-the-art telescopes., Comment: Published in ApJL, includes important updates to stellar and planet parameters

We present new 3.6 and 4.5 micron secondary eclipse measurements for five cool (less than approximately 1000 K) transiting gas giant planets: HAT-P-15b, HAT-P-17b, HAT-P-18b, HAT-P-26b, and WASP-69b. We detect eclipses in at least one bandpass for all planets except HAT-P-15b. We confirm and refine the orbital eccentricity of HAT-P-17b, which is also the only planet in our sample with a known outer companion. We compare our measured eclipse depths in these two bands, which are sensitive to the relative abundances of methane versus carbon monoxide and carbon dioxide, respectively, to predictions from 1D atmosphere models for each planet. For planets with hydrogen-dominated atmospheres and equilibrium temperatures cooler than approximately 1000 K, this ratio should vary as a function of both atmospheric metallicity and the carbon-to-oxygen ratio. For HAT-P-26b, our observations are in good agreement with the low atmospheric metallicity inferred from transmission spectroscopy. We find that all four of the planets with detected eclipses are best matched by models with relatively efficient circulation of energy to the nightside. We see no evidence for a solar-system-like correlation between planet mass and atmospheric metallicity, but instead identify a potential (1.9 sigma) correlation between the inferred methane/(carbon monoxide + carbon dioxide) ratio and stellar metallicity. Our ability to characterize this potential trend is limited by the relatively large uncertainties in the stellar metallicity values. Our observations provide a first look at the brightness of these planets at wavelengths accessible to the James Webb Space Telescope, which will be able to resolve individual methane, carbon monoxide, and carbon dioxide bands and provide much stronger constraints on their atmospheric compositions., Comment: 24 pages, 25 figures, accepted for publication in AJ

In giant planet atmosphere modelling, the intrinsic temperature $T_\mathrm{int}$ and radiative-convective boundary (RCB) are important lower boundary conditions. Often in one-dimensional radiative-convective models and in three-dimensional general circulation models it is assumed that $T_\mathrm{int}$ is similar to that of Jupiter itself, around 100 K, which yields a RCB around 1 kbar for hot Jupiters. In this work, we show that the inflated radii, and hence high specific entropy interiors, of hot Jupiters suggest much higher $T_\mathrm{int}$ values. Assuming the effect is primarily due to current heating (rather than delayed cooling), we derive an equilibrium relation between $T_\mathrm{eq}$ and $T_\mathrm{int}$, showing that the latter can take values as high as 700 K. In response, the RCB moves upward in the atmosphere. Using one-dimensional radiative-convective atmosphere models, we find RCBs of only a few bars, rather than the kilobar typically supposed. This much shallower RCB has important implications for the atmospheric structure, vertical and horizontal circulation, interpretations of phase curves, and the effect of deep cold traps on cloud formation., Comment: 6 pages, 4 figures, accepted to ApJ Letters. This version corrects the constants in equation 3; our calculations used the form in equation 2, so none of our other results are changed. An corresponding erratum has been submitted to ApJ Letters

The exoplanet WASP-12b is the prototype for the emerging class of ultra-hot, Jupiter-mass exoplanets. Past models have predicted---and near ultra-violet observations have shown---that this planet is losing mass. We present an analysis of two sets of 3.6 $\mu$m and 4.5 $\mu$m $\textit{Spitzer}$ phase curve observations of the system which show clear evidence of infrared radiation from gas stripped from the planet, and the gas appears to be flowing directly toward or away from the host star. This accretion signature is only seen at 4.5 $\mu$m, not at 3.6 $\mu$m, which is indicative either of CO emission at the longer wavelength or blackbody emission from cool, $\lesssim$ 600 K gas. It is unclear why WASP-12b is the only ultra-hot Jupiter to exhibit this mass loss signature, but perhaps WASP-12b's orbit is decaying as some have claimed, while the orbits of other exoplanets may be more stable; alternatively, the high energy irradiation from WASP-12A may be stronger than the other host stars. We also find evidence for phase offset variability at the level of $6.4\sigma$ ($46.2^{\circ}$) at 3.6 $\mu$m., Comment: 9 pages plus appendices and supplementary information, published in MNRAS

Launching in 2028, ESA's Atmospheric Remote-sensing Exoplanet Large-survey (ARIEL) survey of $\sim$1000 transiting exoplanets will build on the legacies of Kepler and TESS and complement JWST by placing its high precision exoplanet observations into a large, statistically-significant planetary population context. With continuous 0.5--7.8~$\mu$m coverage from both FGS (0.50--0.55, 0.8--1.0, and 1.0--1.2~$\mu$m photometry; 1.25--1.95~$\mu$m spectroscopy) and AIRS (1.95--7.80~$\mu$m spectroscopy), ARIEL will determine atmospheric compositions and probe planetary formation histories during its 3.5-year mission. NASA's proposed Contribution to ARIEL Spectroscopy of Exoplanets (CASE) would be a subsystem of ARIEL's FGS instrument consisting of two visible-to-infrared detectors, associated readout electronics, and thermal control hardware. FGS, to be built by the Polish Academy of Sciences' Space Research Centre, will provide both fine guiding and visible to near-infrared photometry and spectroscopy, providing powerful diagnostics of atmospheric aerosol contribution and planetary albedo, which play a crucial role in establishing planetary energy balance. The CASE team presents here an independent study of the capabilities of ARIEL to measure exoplanetary metallicities, which probe the conditions of planet formation, and FGS to measure scattering spectral slopes, which indicate if an exoplanet has atmospheric aerosols (clouds and hazes), and geometric albedos, which help establish planetary climate. Our design reference mission simulations show that ARIEL could measure the mass-metallicity relationship of its 1000-planet single-visit sample to $>7.5\sigma$ and that FGS could distinguish between clear, cloudy, and hazy skies and constrain an exoplanet's atmospheric aerosol composition to $>5\sigma$ for hundreds of targets, providing statistically-transformative science for exoplanet atmospheres., Comment: accepted to PASP; 23 pages, 6 figures

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., Comment: Submitted as an Astro2020 Science White Paper

Here we present the first open-source radiative transfer model for computing the reflected light of exoplanets at any phase geometry, called PICASO: Planetary Intensity Code for Atmospheric Scattering Observations. This code, written in Python, has heritage from a decades old, well-known Fortran model used for several studies of planetary objects within the Solar System and beyond. We have adopted it to include several methodologies for computing both direct and diffuse scattering phase functions, and have added several updates including the ability to compute Raman scattering spectral features. Here we benchmark PICASO against two independent codes and discuss the degree to which the model is sensitive to a user's specification for various phase functions. Then, we conduct a full information content study of the model across a wide parameter space in temperature, cloud profile, SNR and resolving power., Comment: Accepted ApJ; GitHub: https://github.com/natashabatalha/picaso Docs: https://natashabatalha.github.io/picaso/ Radiative Transfer Tutorial: https://natashabatalha.github.io/picaso_dev

Several of NASA missions (TESS, JWST, WFIRST, etc.) and mission concepts (LUVOIR, HabEx, and OST) emphasize the exploration and characterization of exoplanets, and the study of the interstellar medium. We anticipate that a much broader set of chemical environments exists on exoplanets, necessitating data from a correspondingly broader set of chemical reactions. Similarly, the conditions that exist in astrophysical environments are very different from those traditionally probed in laboratory chemical kinetics studies. These are areas where quantum mechanical theory, applied to important reactions via well-validated chemical kinetics models, can fill a critical knowledge gap. Quantum chemical calculations are also introduced to study interior of planets, photochemical escape, and many critical chemical pathways (e.g. prebiotic environments, contaminations, etc.) After years of development of the relevant quantum chemical theories and significant advances in computational power, quantum chemical simulations have currently matured enough to describe real systems with an accuracy that competes with experiments. These approaches, therefore, have become the best possible alternative in many circumstances where performing experiments is too difficult, too expensive, or too dangerous, or simply not possible. In this white paper, several existing quantum chemical studies supporting exoplanetary science, planetary astronomy, and astrophysics are described, and the potential positive impacts of improved models associated with scientific goals of missions are addressed. In the end, a few recommendations from the scientific community to strengthen related research efforts at NASA are provided.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^\alpha a^\beta$, finding $\alpha$ = -2.4 $\pm$ 0.8 and $\beta$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability., Comment: 52 pages, 18 figures. AJ in press

Here we present the analysis of multi-epoch secondary eclipse observations of HD 189733b and HD 209458b as a probe of temporal variability in the planetary climate using both Spitzer channels 1 and 2 (3.6 and 4.5 um). Constraining temporal variability will inform models and identify physical processes occurring at either length scales too small to directly observe or at pressure levels that are inaccessible to transit observations. We do not detect statistically significant variability and are able to place useful upper limits on the IR variability amplitudes in these atmospheres. There are very few planets with multi-epoch observations at the required precision to probe variability in dayside emission. The observations considered in this study span several years, providing insight into temporal variability at multiple timescales. In the case of HD 189733b, the best fit eclipse depths for the channel 2 observations exhibit a scatter of 102 ppm about a median depth of 1827 ppm and in channel 1 exhibit a scatter of 88 ppm about a median depth of 1481 ppm. For HD 209458b, the best fit eclipse depths for the channel 2 observations exhibit a scatter of 22 ppm about a median depth of 1406 ppm and in channel 1 exhibit a scatter of 131 ppm about a median depth of 1092 ppm. The precision and scatter in these observations allow us to constrain variability to less than (5.6% and 6.0%) and (12% and 1.6%) for channels (1,2) of HD 189733b and HD 209458b respectively. There is a difference in the best fit eclipse timing compared to the predicted time consistent with an offset hotspot as predicted by GCMs and confirmed in previous phase curve observations., Comment: Accepted for Publication in AJ 12/17/2019, 9 pages, 8 figures

We present the analysis of a full-orbit, spectroscopic phase curve of the ultra hot Jupiter WASP-18b, obtained with the Wide Field Camera 3 aboard the Hubble Space Telescope. We measure the planet's normalized day-night contrast as >0.96 in luminosity: the disk-integrated dayside emission from the planet is at 964+-25 ppm, corresponding to 2894+-30 K, and we place an upper limit on the nightside emission of <32ppm or 1430K at the 3-sigma level. We also find that the peak of the phase curve exhibits a small, but significant offset in brightness of 4.5+-0.5 degrees eastward. We compare the extracted phase curve and phase resolved spectra to 3D Global Circulation Models and find that broadly the data can be well reproduced by some of these models. We find from this comparison several constraints on the atmospheric properties of the planet. Firstly we find that we need efficient drag to explain the very inefficient day-night re-circulation observed. We demonstrate that this drag could be due to Lorentz-force drag by a magnetic field as weak as 10 Gauss. Secondly, we show that a high metallicity is not required to match the large day-night temperature contrast. In fact, the effect of metallicity on the phase curve is different from cooler gas-giant counterparts, due to the high-temperature chemistry in WASP-18b's atmosphere. Additionally, we compare the current UHJ spectroscopic phase curves, WASP-18b and WASP-103b, and show that these two planets provide a consistent picture with remarkable similarities in their measured and inferred properties. However, key differences in these properties, such as their brightness offsets and radius anomalies, suggest that UHJ could be used to separate between competing theories for the inflation of gas-giant planets., Comment: Accepted for publication in A&A

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., Comment: Submitted to ApJ Letters, revised to reflect referee comments

Time-series spectrophotometric studies of exoplanets during transit using ground-based facilities are a promising approach to characterize their atmospheric compositions. We aim to investigate the transit spectrum of the hot Jupiter HAT-P-1b. We compare our results to those obtained at similar wavelengths by previous space-based observations. We observed two transits of HAT-P-1b with the Gemini Multi-Object Spectrograph (GMOS) instrument on the Gemini North telescope using two instrument modes covering the 320 - 800 nm and 520 - 950 nm wavelength ranges. We used time-series spectrophotometry to construct transit light curves in individual wavelength bins and measure the transit depths in each bin. We accounted for systematic effects. We addressed potential photometric variability due to magnetic spots in the planet's host star with long-term photometric monitoring. We find that the resulting transit spectrum is consistent with previous Hubble Space Telescope (HST) observations. We compare our observations to transit spectroscopy models that marginally favor a clear atmosphere. However, the observations are also consistent with a flat spectrum, indicating high-altitude clouds. We do not detect the Na resonance absorption line (589 nm), and our observations do not have sufficient precision to study the resonance line of K at 770 nm. We show that even a single Gemini/GMOS transit can provide constraining power on the properties of the atmosphere of HAT-P-1b to a level comparable to that of HST transit studies in the optical when the observing conditions and target and reference star combination are suitable. Our 520 - 950 nm observations reach a precision comparable to that of HST transit spectra in a similar wavelength range of the same hot Jupiter, HAT-P-1b. However, our GMOS transit between 320 - 800 nm suffers from strong systematic effects and yields larger uncertainties., Comment: A&A, accepted, 16 pages, 8 figures, 5 tables

We report 78 secondary eclipse depths for a sample of 36 transiting hot Jupiters observed at 3.6- and 4.5 microns using the Spitzer Space Telescope. Our eclipse results for 27 of these planets are new, and include highly irradiated worlds such as KELT-7b, WASP-87b, WASP-76b, and WASP-64b, and important targets for JWST such as WASP-62b. We find that WASP-62b has a slightly eccentric orbit e cos(omega) = 0.00614+/- 0.00064, and we confirm the eccentricity of HAT-P-13b and WASP-14b. The remainder are individually consistent with circular orbits, but we find statistical evidence for eccentricity increasing with orbital period in our range from 1 to 5 days. Our day-side brightness temperatures for the planets yield information on albedo and heat redistribution, following Cowan and Agol (2011). Planets having maximum day side temperatures exceeding ~ 2200K are consistent with zero albedo and distribution of stellar irradiance uniformly over the day-side hemisphere. Our most intriguing result is that we detect a systematic difference between the emergent spectra of these hot Jupiters as compared to blackbodies. The ratio of observed brightness temperatures, Tb(4.5)/Tb(3.6), increases with equilibrium temperature by 100 +/- 24 parts-per-million per Kelvin, over the entire temperature range in our sample (800K to 2500K). No existing model predicts this trend over such a large range of temperature. We suggest that this may be due to a structural difference in the atmospheric temperature profile between the real planetary atmospheres as compared to models., Comment: accepted for publication in AJ

We present high-resolution optical transmission spectroscopy of two sub-Saturn mass transiting exoplanets, HAT-P-12b and WASP-69b. With relatively low densities and high atmospheric scale heights, these planets are particularly well-suited to characterization through transit spectroscopy, and serve as ideal candidates for extending previously-tested methods to lower planetary masses. Using a single transit for each planet, we take advantage of the Doppler cross-correlation technique to search for sodium, potassium, and water absorption features. Our analysis reveals a likely ($3.2\sigma$) detection of sodium absorption features in the atmosphere of HAT-P-12b, and enables us to place constraints on the presence of alkaline and molecular species in the atmospheres of both planets. With our results, we highlight the efficacy of ground-based campaigns for characterizing exoplanetary atmospheres, and pave the way for future analyses of low-mass planets., Comment: 15 pages, 11 figures

The helium absorption triplet at a wavelength of 10,833 \AA\ has been proposed as a way to probe the escaping atmospheres of exoplanets. Recently this feature was detected for the first time using Hubble Space Telescope (HST) WFC3 observations of the hot Jupiter WASP-107b. We use similar HST/WFC3 observations to detect helium in the atmosphere of the hot Neptune HAT-P-11b at the $4\sigma$ confidence level. We compare our observations to a grid of 1D models of hydrodynamic escape to constrain the thermospheric temperatures and mass loss rate. We find that our data are best fit by models with high mass loss rates of $\dot{M} \approx 10^{9}$ - $10^{11}$ g s$^{-1}$. Although we do not detect the planetary wind directly, our data are consistent with the prediction that HAT-P-11b is experiencing hydrodynamic atmospheric escape. Nevertheless, the mass loss rate is low enough that the planet has only lost up to a few percent of its mass over its history, leaving its bulk composition largely unaffected. This matches the expectation from population statistics, which indicate that close-in planets with radii greater than 2 R$_{\oplus}$ form and retain H/He-dominated atmospheres. We also confirm the independent detection of helium in HAT-P-11b obtained with the CARMENES instrument, making this the first exoplanet with the detection of the same signature of photoevaporation from both ground- and space-based facilities., Comment: 12 pages, 9 figures, accepted for publication in ApJL

Atmospheric characterization through spectroscopic analysis, an essential tool of modern exoplanet science, can benefit significantly from the context provided by interior structure models. In particular, the planet's bulk metallicity $Z_p$ places an upper limit on potential atmospheric metallicity. Here we construct interior structure models to derive $Z_p$ and atmospheric metallicity upper limits for 403 known transiting giant exoplanets. These limits are low enough that they can usefully inform atmosphere models. Additionally, we argue that comparing $Z_p$ to the observed atmospheric metallicity gives a useful measure of how well-mixed metals are within the planet. This represents a new avenue for learning about planetary interiors. To aid in the future characterization of new planet discoveries we derive analytic prior predictions of atmosphere metallicity as a function of planet mass, and evaluate the effectiveness of our approach on Jupiter and Saturn. We we include log-linear fits for approximating the metallicities of planets not in our catalog., Comment: Submitted to ApJ. Revision includes more planets and clarity improvements, scientific results unchanged