WASP-39b is a hot Saturn-mass exoplanet with a predicted clear atmosphere based on observations in the optical and infrared. Here, we complete the transmission spectrum of the atmosphere with observations in the near-infrared (NIR) over three water absorption features with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) G102 (0.8-1.1 mu m) and G141 (1.1-1.7 mu m) spectroscopic grisms. We measure the predicted high-amplitude H2O feature centered at 1.4 mu m and the smaller amplitude features at 0.95 and 1.2 mu m, with a maximum water absorption amplitude of 2.4 planetary scale heights. We incorporate these new NIR measurements into previously published observational measurements to complete the transmission spectrum from 0.3 to 5 mu m. From these observed water features, combined with features in the optical and IR, we retrieve a well constrained temperature T-eq = 1030(20)(+30) K, and atmospheric metallicity 151(46) (+48) solar, which is relatively high with respect to the currently established mass-metallicity trends. This new measurement in the Saturn-mass range hints at further diversity in the planet formation process relative to our solar system giants.

We present simultaneous Hubble Space Telescope (HST) WFC3+Spitzer IRAC variability monitoring for the highly variable young (similar to 20 Myr) planetary-mass object PSO J318.5-22. Our simultaneous HST + Spitzer observations covered approximately two rotation periods with Spitzer and most of a rotation period with the HST. We derive a period of 8.6. +/-. 0.1 hr from the Spitzer light curve. Combining this period with the measuredvsinifor this object, we find an inclination of 56 degrees.2. +/-. 8 degrees.1. We measure peak-to-trough variability amplitudes of 3.4%. +/-. 0.1% for Spitzer Channel 2 and 4.4%-5.8% (typical 68% confidence errors of similar to 0.3%) in the near-IR bands (1.07-1.67 mu m) covered by the WFC3 G141 prism-the mid-IR variability amplitude for PSO J318.5-22 is one of the highest variability amplitudes measured in the mid-IR for any brown dwarf or planetary-mass object. Additionally, we detect phase offsets ranging from 200 degrees to 210 degrees (typical error of similar to 4 degrees) between synthesized near-IR light curves and the Spitzer mid-IR light curve, likely indicating depth-dependent longitudinal atmospheric structure in this atmosphere. The detection of similar variability amplitudes in wide spectral bands relative to absorption features suggests that the driver of the variability may be inhomogeneous clouds (perhaps a patchy haze layer over thick clouds), as opposed to hot spots or compositional inhomogeneities at the top-of-atmosphere level.

The determination of atmospheric structure and molecular abundances of planetary atmospheres via spectroscopy involves direct comparisons between models and data. While varying in sophistication, most model spectra comparisons fundamentally assume one-dimensional (1D) model physics. However, knowledge from general circulation models and of solar system planets suggests that planetary atmospheres are inherently three-dimensional in their structure and composition. We explore the potential biases resulting from standard "1D" assumptions within a Bayesian atmospheric retrieval framework. Specifically, we show how the assumption of a single 1D thermal profile can bias our interpretation of the thermal emission spectrum of a hot Jupiter atmosphere that is composed of two thermal profiles. We retrieve spectra of unresolved model planets as observed with a combination of the Hubble Space Telescope Wide Field Camera 3 (WFC3)+Spitzer Infrared Array Camera (IRAC) as well as the James Webb Space Telescope (JWST) under varying differences in the two thermal profiles. For WFC3+IRAC, there is a significantly biased estimate of CH4 abundance using a 1D model when the contrast is 80%. For JWST, two thermal profiles are required to adequately interpret the data and estimate the abundances when contrast is greater than 40%. We also apply this preliminary concept to the recent WFC3+IRAC phase curve data of the hot Jupiter WASP-43b. We see similar behavior as present in our simulated data: while the H2O abundance determination is robust, CH4 is artificially well-constrained to incorrect values under the 1D assumption. Our work demonstrates the need to evaluate model assumptions in order to extract meaningful constraints from atmospheric spectra and motivates exploration of optimal observational setups.

We target the thermal emission spectrum of the non-transiting gas giant HD 88133 b with high-resolution near-infrared spectroscopy, by treating the planet and its host star as a spectroscopic binary. For sufficiently deep summed flux observations of the star and planet across multiple epochs, it is possible to resolve the signal of the hot gas giant's atmosphere compared to the brighter stellar spectrum, at a level consistent with the aggregate shot noise of the full data set. To do this, we first perform a principal component analysis to remove the contribution of the Earth's atmosphere to the observed spectra. Then, we use a cross-correlation analysis to tease out the spectra of the host star and HD 88133 b to determine its orbit and identify key sources of atmospheric opacity. In total, six epochs of Keck NIRSPEC L-band observations and three epochs of Keck NIRSPEC K-band observations of the HD 88133 system were obtained. Based on an analysis of the maximum likelihood curves calculated from the multi-epoch cross-correlation of the full data set with two atmospheric models, we report the direct detection of the emission spectrum of the non-transiting exoplanet HD 88133 b and measure a radial projection of the Keplerian orbital velocity of 40 +/- 15 km s(-1), a true mass of 1.02(-0.28)(+0.61) M-J, a nearly face-on orbital inclination of 15(-5)(+60), and an atmosphere opacity structure at high dispersion dominated by water vapor. This, combined with 11 years of radial velocity measurements of the system, provides the most up-to-date ephemeris for HD 88133.

The nature of the thermal structure of hot Jupiter atmospheres is one of the key questions raised by the characterization of transiting exoplanets over the past decade. There have been claims that many hot Jupiters exhibit atmospheric thermal inversions. However, these claims have been based on broadband photometry rather than the unambiguous identification of emission features with spectroscopy, and the chemical species that could cause the thermal inversions by absorbing stellar irradiation at high altitudes have not been identified despite extensive theoretical and observational effort. Here we present high-precision Hubble Space Telescope WFC3 observations of the dayside thermal emission spectrum of the hot Jupiter HD 209458b, which was the first exoplanet suggested to have a thermal inversion. In contrast to previous results for this planet, our observations detect water in absorption at 6.2 sigma confidence. When combined with Spitzer photometry, the data are indicative of a monotonically decreasing temperature with pressure over the range of 1-0.001 bars at 7.7 sigma confidence. We test the robustness of our results by exploring a variety of model assumptions, including the temperature profile parameterization, presence of a cloud, and choice of Spitzer data reduction. We also introduce a new analysis method to determine the elemental abundances from the spectrally retrieved mixing ratios with thermochemical self-consistency and find plausible abundances consistent with solar metallicity (0.06-10 x solar) and carbon-to oxygen ratios less than unity. This work suggests that high-precision spectrophotometric results are required to robustly infer thermal structures and compositions of extrasolar planet atmospheres and to perform comparative exoplanetology.

Planetesimals may form from the gravitational collapse of dense particle clumps initiated by the streaming instability. We use simulations of aerodynamically coupled gas-particle mixtures to investigate whether the properties of planetesimals formed in this way depend upon the sizes of the particles that participate in the instability. Based on three high-resolution simulations that span a range of dimensionless stopping times 6 X 10(-3) <= tau <= 2, no statistically significant differences in the initial planetesimal mass function are found. The mass functions are fit by a power law, dN/dM(p) proportional to M-p(-p), with p = 1.5-1.7 and errors of Delta p approximate to 0.1. Comparing the particle density fields prior to collapse, we find that the high-wavenumber power spectra are similarly indistinguishable, though the large-scale geometry of structures induced via the streaming instability is significantly different between all three cases. We interpret the results as evidence for a near-universal slope to the mass function, arising from the small-scale structure of streaming-induced turbulence.

Context. Irregular satellites of Jupiter with their highly eccentric, inclined and distant orbits suggest that their capture took place after the giant planet migration. Aims. We aim to improve our understanding of the surface composition of irregular satellites of Jupiter to gain insight into a narrow time window when our solar system was forming. Methods. We observed three Jovian irregular satellites, Himalia (JVI), Elara (JVII), and Carme (JXI), using a medium-resolution 0.8-5.5 mu m spectrograph, SpeX on the NASA Infrared Telescope Facility (IRTF). Using a linear spectral unmixing model we have constrained the major mineral phases on the surface of these three bodies. Results. Our results confirm that the surface of Himalia (JVI), Elara (JVII), and Carme (JXI) are dominated by opaque materials such as those seen in carbonaceous chondrite meteorites. Our spectral modeling of NIR spectra of Himalia and Elara confirm that their surface composition is the same and magnetite is the dominant mineral. A comparison of the spectral shape of Himalia with the two large main C-type asteroids, Themis (D similar to 176 km) and Europa (D similar to 352 km), suggests surface composition similar to Europa. The NIR spectrum of Carme exhibits blue slope up to 1.5 mu m and is spectrally distinct from those of Himalia and Elara. Our model suggests that it is compositionally similar to amorphous carbon. Conclusions. Himalia and Elara are compositionally similar but di ff er significantly from Carme. These results support the hypotheses that the Jupiter's irregular satellites are captured bodies that were subject to further breakup events and clustered as families based on their similar physical and surface compositions.

Hot Jupiters receive strong stellar irradiation, producing equilibrium temperatures of 1000-2500 K. Incoming irradiation directly heats just their thin outer layer, down to pressures of similar to 0.1 bars. In standard irradiated evolution models of hot Jupiters, predicted transit radii are too small. Previous studies have shown that deeper heating-at a small fraction of the heating rate from irradiation-can explain observed radii. Here we present a suite of evolution models for HD 209458b, where we systematically vary both the depth and intensity of internal heating, without specifying the uncertain heating mechanism(s). Our models start with a hot, high-entropy planet whose radius decreases as the convective interior cools. The applied heating suppresses this cooling. We find that very shallow heating-at pressures of 1-10 bars-does not significantly suppress cooling, unless the total heating rate is greater than or similar to 10% of the incident stellar power. Deeper heating, at 100 bars, requires heating at only 1% of the stellar irradiation to explain the observed transit radius of 1.4R(Jup) after 5 Gyr of cooling. In general, more intense and deeper heating results in larger hot-Jupiter radii. Surprisingly, we find that heat deposited at 10(4) bars-which is exterior to approximate to 99% of the planet's mass-suppresses planetary cooling as effectively as heating at the center. In summary, we find that relatively shallow heating is required to explain the radii of most hot Jupiters, provided that this heat is applied early and persists throughout their evolution.

The Upsilon Andromedae system was the first multi-planet system discovered orbiting a main-sequence star. We describe the detection of water vapor in the atmosphere of the innermost non-transiting gas giant ups. And. b by treating the star-planet system as a spectroscopic binary with high-resolution, ground-based spectroscopy. We resolve the signal of the planet's motion and break the mass-inclination degeneracy for this non-transiting planet via deep combined flux observations of the star and the planet. In total, seven epochs of Keck NIRSPEC L band observations, three epochs of Keck NIRSPEC short-wavelength K band observations, and three epochs of Keck NIRSPEC long wavelength K band observations of the ups. And. system were obtained. We perform a multi-epoch cross-correlation of the full data set with an atmospheric model. We measure the radial projection of the Keplerian velocity (K-P = 55 +/- 9 km s(-1)), true mass (M-b = 1.7(-0.24)(+0.33)M(J)), and orbital inclination (i(b) 24 degrees +/- 4 degrees), and determine that the planet's opacity structure is dominated by water vapor at the probed wavelengths. Dynamical simulations of the planets in the ups. And. system with these orbital elements for ups. And. b show that stable, long-term (100 Myr) orbital configurations exist. These measurements will inform future studies of the stability and evolution of the ups. And. system, as well as the atmospheric structure and composition of the hot Jupiter.

As followup to our recent detection of a pre-transit signal around HD 189733 b, we obtained full pre-transit phase coverage of a single planetary transit. The pre-transit signal is again detected in the Balmer lines but with variable strength and timing, suggesting that the bow shock geometry reported in our previous work does not describe the signal from the latest transit. We also demonstrate the use of the Ca II H and K residual core flux as a proxy for the stellar activity level throughout the transit. A moderate trend is found between the pre-transit absorption signal in the 2013 data and the Ca II H flux. This suggests that some of the 2013 pre-transit hydrogen absorption can be attributed to varying stellar activity levels. A very weak correlation is found between the Ca II H core flux and the Balmer line absorption in the 2015 transit, hinting at a smaller contribution from stellar activity compared to the 2013 transit. We simulate how varying stellar activity levels can produce changes in the Balmer line transmission spectra. These simulations show that the strength of the 2013 and 2015 pre-transit signals can be reproduced by stellar variability. If the pre-transit signature is attributed to circumplanetary material, its evolution in time can be described by accretion clumps spiraling toward the star, although this interpretation has serious limitations. Further high-cadence monitoring at H alpha is necessary to distinguish between true absorption by transiting material and short-term variations in the stellar activity level.