Search

Here we present the first open-source radiative transfer model for computing the reflected light of exoplanets at any phase geometry, called PICASO: the 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, signal-to-noise ratio, and resolving power. [ABSTRACT FROM AUTHOR]

Seismology of the gas giants holds the potential to resolve long-standing questions about their internal structure and rotation state. We construct a family of Saturn interior models constrained by the gravity field and compute their adiabatic mode eigenfrequencies and corresponding Lindblad and vertical resonances in Saturn’s C ring, where more than 20 waves with pattern speeds faster than the ring mean motion have been detected and characterized using high-resolution Cassini Visual and Infrared Mapping Spectrometer stellar occultation data. We present identifications of the fundamental modes of Saturn that appear to be the origin of these observed ring waves, and we use their observed pattern speeds and azimuthal wavenumbers to estimate the bulk rotation period of Saturn’s interior to be (median and 5%/95% quantiles), significantly faster than Voyager and Cassini measurements of periods in Saturn’s kilometric radiation, the traditional proxy for Saturn’s bulk rotation period. The global fit does not exhibit any clear systematics indicating strong differential rotation in Saturn’s outer envelope. [ABSTRACT FROM AUTHOR]

The H2O abundance of a planetary atmosphere is a powerful indicator of formation conditions. Inferring H2O in the solar system giant planets is challenging, due to condensation depleting the upper atmosphere of water vapor. Substantially warmer hot Jupiter exoplanets readily allow detections of H2O via transmission spectroscopy, but such signatures are often diminished by the presence of clouds composed of other species. In contrast, highly scattering water clouds can brighten planets in reflected light, enhancing molecular signatures. Here, we present an extensive parameter space survey of the prominence of H2O absorption features in reflection spectra of cool (Teff < 400 K) giant exoplanetary atmospheres. The impact of effective temperature, gravity, metallicity, and sedimentation efficiency is explored. We find prominent H2O features around 0.94 μm, 0.83 μm, and across a wide spectral region from 0.4 to 0.73 μm. The 0.94 μm feature is only detectable where high-altitude water clouds brighten the planet: Teff ∼ 150 K, g ≳ 20 ms−2, fsed ≳ 3, m ≲ 10× solar. In contrast, planets with g ≲ 20 ms−2 and Teff ≳ 180 K display substantially prominent H2O features embedded in the Rayleigh scattering slope from 0.4 to 0.73 μm over a wide parameter space. High fsed enhances H2O features around 0.94 μm, and enables these features to be detected at lower temperatures. High m results in dampened H2O absorption features, due to water vapor condensing to form bright, optically thick clouds that dominate the continuum. We verify these trends via self-consistent modeling of the low-gravity exoplanet HD 192310c, revealing that its reflection spectrum is expected to be dominated by H2O absorption from 0.4 to 0.73 μm for m ≲ 10× solar. Our results demonstrate that H2O is manifestly detectable in reflected light spectra of cool giant planets only marginally warmer than Jupiter, providing an avenue to directly constrain the C/O and O/H ratios of a hitherto unexplored population of exoplanetary atmospheres. [ABSTRACT FROM AUTHOR]

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. presented a spectrum of WISE 0855 from 4.5–5.1 μm (M band), revealing water vapor features. Here, we present a spectrum of WISE 0855 in the L band, from 3.4–4.14 μm. 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. The James Webb Space Telescope 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. [ABSTRACT FROM AUTHOR]

Condensation clouds in substellar atmospheres have been widely inferred from spectra and photometric variability. Up until now, their horizontally averaged vertical distribution and mean particle size have been largely characterized using models, one of which is the eddy diffusion–sedimentation model from Ackerman and Marley that relies on a sedimentation efficiency parameter, fsed, to determine the vertical extent of clouds in the atmosphere. However, the physical processes controlling the vertical structure of clouds in substellar atmospheres are not well understood. In this work, we derive trends in fsed across a large range of eddy diffusivities (Kzz), gravities, material properties, and cloud formation pathways by fitting cloud distributions calculated by a more detailed cloud microphysics model. We find that fsed is dependent on Kzz, but not gravity, when Kzz is held constant. fsed is most sensitive to the nucleation rate of cloud particles, as determined by material properties like surface energy and molecular weight. High surface energy materials form fewer, larger cloud particles, leading to large fsed (>1), and vice versa for materials with low surface energy. For cloud formation via heterogeneous nucleation, fsed is sensitive to the condensation nuclei flux and radius, connecting cloud formation in substellar atmospheres to the objects’ formation environments and other atmospheric aerosols. These insights could lead to improved cloud models that help us better understand substellar atmospheres. For example, we demonstrate that fsed could increase with increasing cloud base depth in an atmosphere, shedding light on the nature of the brown dwarf L/T transition. [ABSTRACT FROM AUTHOR]

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–800 K) near-infrared (1–2.5 μm) spectra. With these spectra we are able to directly constrain the molecular abundances for the first time of H2O, CH4, CO, CO2, NH3, H2S, and Na+K, surface gravity, effective temperature, thermal structure, photometric radius, and cloud optical depths. We find that ammonia, water, methane, and the alkali metals are present and that their abundances are well constrained in all 11 objects. We find no significant trend in the water, methane, or ammonia abundances with temperature, but find a very strong (>25σ) decreasing trend in the alkali metal abundances with decreasing effective temperature, indicative of alkali rainout. As expected from previous work, 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 subsolar (−0.4 < [M/H] < 0.1 dex) and carbon-to-oxygen ratios are somewhat supersolar (0.4 < C/O < 1.2), 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 some of our radii are larger than expectations from evolutionary models, possibly indicative of unresolved binaries. This investigation and method represent a new and powerful paradigm for using spectra to determine the fundamental chemical and physical processes governing cool brown dwarf atmospheres. [ABSTRACT FROM AUTHOR]

We present a new matched-filter algorithm for direct detection of point sources in the immediate vicinity of bright stars. The stellar point-spread function (PSF) is first subtracted using a Karhunen-Loéve image processing (KLIP) algorithm with angular and spectral differential imaging (ADI and SDI). The KLIP-induced distortion of the astrophysical signal is included in the matched-filter template by computing a forward model of the PSF at every position in the image. To optimize the performance of the algorithm, we conduct extensive planet injection and recovery tests and tune the exoplanet spectra template and KLIP reduction aggressiveness to maximize the signal-to-noise ratio (S/N) of the recovered planets. We show that only two spectral templates are necessary to recover any young Jovian exoplanets with minimal S/N loss. We also developed a complete pipeline for the automated detection of point-source candidates, the calculation of receiver operating characteristics (ROC), contrast curves based on false positives, and completeness contours. We process in a uniform manner more than 330 data sets from the Gemini Planet Imager Exoplanet Survey and assess GPI typical sensitivity as a function of the star and the hypothetical companion spectral type. This work allows for the first time a comparison of different detection algorithms at a survey scale accounting for both planet completeness and false-positive rate. We show that the new forward model matched filter allows the detection of 50% fainter objects than a conventional cross-correlation technique with a Gaussian PSF template for the same false-positive rate. [ABSTRACT FROM AUTHOR]

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 the signal-to-noise ratio (S/N) of the observed spectra, particularly for low metallicity atmospheres, which lack deep absorption signatures. Meaningful results may only be possible above certain S/N thresholds; for cases across a metallicity range of 1–50 times solar, we find that only an S/N of 20 systematically reproduces a value 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. [ABSTRACT FROM AUTHOR]

We examine the hypothesis that the red near-infrared colors of some L dwarfs could be explained by a “dust haze” of small particles in their upper atmospheres. This dust haze would exist in conjunction with the clouds found in dwarfs with more typical colors. We developed a model that uses Mie theory and the Hansen particle size distributions to reproduce the extinction due to the proposed dust haze. We apply our method to 23 young L dwarfs and 23 red field L dwarfs. We constrain the properties of the dust haze including particle size distribution and column density using Markov Chain Monte Carlo methods. We find that sub-micron-range silicate grains reproduce the observed reddening. Current brown dwarf atmosphere models include large-grain (1–100 μm) dust clouds but not sub-micron dust grains. Our results provide a strong proof of concept and motivate a combination of large and small dust grains in brown dwarf atmosphere models. [ABSTRACT FROM AUTHOR]

Over a large range of equilibrium temperatures, clouds shape the transmission spectrum of hot Jupiter atmospheres, yet their composition remains unknown. Recent observations show that the Kepler light curves of some hot Jupiters are asymmetric: for the hottest planets, the light curve peaks before secondary eclipse, whereas for planets cooler than ∼1900 K, it peaks after secondary eclipse. We use the thermal structure from 3D global circulation models to determine the expected cloud distribution and Kepler light curves of hot Jupiters. We demonstrate that the change from an optical light curve dominated by thermal emission to one dominated by scattering (reflection) naturally explains the observed trend from negative to positive offset. For the cool planets the presence of an asymmetry in the Kepler light curve is a telltale sign of the cloud composition, because each cloud species can produce an offset only over a narrow range of effective temperatures. By comparing our models and the observations, we show that the cloud composition of hot Jupiters likely varies with equilibrium temperature. We suggest that a transition occurs between silicate and manganese sulfide clouds at a temperature near 1600 K, analogous to the L/T transition on brown dwarfs. The cold trapping of cloud species below the photosphere naturally produces such a transition and predicts similar transitions for other condensates, including TiO. We predict that most hot Jupiters should have cloudy nightsides, that partial cloudiness should be common at the limb, and that the dayside hot spot should often be cloud-free. [ABSTRACT FROM AUTHOR]

We present Spitzer/Infrared Array Camera Ch1 and Ch2 monitoring of six brown dwarfs during eight different epochs over the course of 20 months. For four brown dwarfs, we also obtained simulataneous Hubble Space Telescope (HST)/WFC3 G141 grism spectra during two epochs and derived light curves in five narrowband filters. Probing different pressure levels in the atmospheres, the multiwavelength light curves of our six targets all exhibit variations, and the shape of the light curves evolves over the timescale of a rotation period, ranging from 1.4 to 13 hr. We compare the shapes of the light curves and estimate the phase shifts between the light curves observed at different wavelengths by comparing the phase of the primary Fourier components. We use state-of-the-art atmosphere models to determine the flux contribution of different pressure layers to the observed flux in each filter. We find that the light curves that probe higher pressures are similar and in phase, but are offset and often different from the light curves that probe lower pressures. The phase differences between the two groups of light curves suggest that the modulations seen at lower and higher pressures may be introduced by different cloud layers. [ABSTRACT FROM AUTHOR]

WISE J104915.57-531906.1 is the nearest brown dwarf binary to our solar system, consisting of two brown dwarfs in the L/T transition: Luhman 16A and B. In this paper, we present the first map of Luhman 16A, and maps of Luhman 16B for two epochs. Our maps were created by applying Aeolus, a Markov-Chain Monte Carlo code that maps the top-of-the-atmosphere (TOA) structure of brown dwarf and other ultracool atmospheres, to light curves of Luhman 16A and B using the Hubble Space Telescope’s G141 and G102 grisms. Aeolus retrieved three or four spots in the TOA of Luhman 16A and B, with a surface coverage of 19%–32% (depending on an assumed rotational period of 5 hr or 8 hr) or 21%–38.5% (depending on the observational epoch), respectively. The brightness temperature of the spots of the best-fit models was ∼200 K hotter than the background TOA. We compared our Luhman 16B map with the only previously published map. Interestingly, our map contained a large TOA spot that was cooler (ΔT ∼ 51 K) than the background, which lay at low latitudes, in agreement with the previous Luhman 16B map. Finally, we report the detection of a feature reappearing in Luhman 16B light curves that are separated by tens of hundreds of rotations from each other. We speculate that this feature is related to TOA structures of Luhman 16B. [ABSTRACT FROM AUTHOR]

Many of the directly imaged self-luminous gas-giant exoplanets have been found to have cloudy atmospheres. Scattering of the emergent thermal radiation from these planets by the dust grains in their atmospheres should locally give rise to significant linear polarization of the emitted radiation. However, the observable disk-averaged polarization should be zero if the planet is spherically symmetric. Rotation-induced oblateness may yield a net non-zero disk-averaged polarization if the planets have sufficiently high spin rotation velocity. On the other hand, when a large natural satellite or exomoon transits a planet with a cloudy atmosphere along the line of sight, the asymmetry induced during the transit should give rise to a net non-zero, time-resolved linear polarization signal. The peak amplitude of such time-dependent polarization may be detectable even for slowly rotating exoplanets. Therefore, we suggest that large exomoons around directly imaged self-luminous exoplanets may be detectable through time-resolved imaging polarimetry. Adopting detailed atmospheric models for several values of effective temperature and surface gravity that are appropriate for self-luminous exoplanets, we present the polarization profiles of these objects in the infrared during the transit phase and estimate the peak amplitude of polarization that occurs during the inner contacts of the transit ingress/egress phase. The peak polarization is predicted to range between 0.1% and 0.3% in the infrared. [ABSTRACT FROM AUTHOR]

We present the first detection of the photometric variability in a spectroscopically confirmed Y dwarf. The Infrared Array Camera on board the Spitzer Space Telescope was used to obtain time series photometry of WISE J140518.39+553421.3 at 3.6 and 4.5 μm over a 24-hr period at two different epochs separated by 149 days. Variability is evident at 4.5 μm in the first epoch and at 3.6 and 4.5 μm in the second epoch, which suggests that the underlying cause or causes of this variability change on the timescales of months. The second-epoch [3.6] and [4.5] light curves are nearly sinusoidal in form, in phase, have periods of roughly 8.5 hr, and have semi-amplitudes of 3.5%. We find that a simple geometric spot model with a single bright spot reproduces these observations well. We also compare our measured semi-amplitudes of the second-epoch light curves to predictions of the static, one-dimensional, partly cloudy, and hot spot models of Morley and collaborators, and find that neither set of models can reproduce the observed [3.6] and [4.5] semi-amplitudes simultaneously. Therefore, more advanced two-dimensional or three-dimensional models that include time-dependent phenomena like vertical mixing, cloud formation, and thermal relaxation are sorely needed in order to properly interpret our observations. [ABSTRACT FROM AUTHOR]

We present results from an atmospheric circulation study of nine hot Jupiters that compose a large transmission spectral survey using the Hubble and Spitzer Space Telescopes. These observations exhibit a range of spectral behavior over optical and infrared wavelengths, suggesting diverse cloud and haze properties in their atmospheres. By utilizing the specific system parameters for each planet, we naturally probe a wide phase space in planet radius, gravity, orbital period, and equilibrium temperature. First, we show that our model “grid” recovers trends shown in traditional parametric studies of hot Jupiters, particularly equatorial superrotation and increased day–night temperature contrast with increasing equilibrium temperature. We show how spatial temperature variations, particularly between the dayside and nightside and west and east terminators, can vary by hundreds of kelvin, which could imply large variations in Na, K, CO and abundances in those regions. These chemical variations can be large enough to be observed in transmission with high-resolution spectrographs, such as ESPRESSO on VLT, METIS on the E-ELT, or MIRI and NIRSpec aboard JWST. We also compare theoretical emission spectra generated from our models to available Spitzer eclipse depths for each planet and find that the outputs from our solar-metallicity, cloud-free models generally provide a good match to many of the data sets, even without additional model tuning. Although these models are cloud-free, we can use their results to understand the chemistry and dynamics that drive cloud formation in their atmospheres. [ABSTRACT FROM AUTHOR]

Detecting polarized light from self-luminous exoplanets has the potential to provide key information about rotation, surface gravity, cloud grain size, and cloud coverage. While field brown dwarfs with detected polarized emission are common, no exoplanet or substellar companion has yet been detected in polarized light. With the advent of high contrast imaging spectro-polarimeters such as GPI and SPHERE, such a detection may now be possible with careful treatment of instrumental polarization. In this paper, we present 28 minutes of H-band GPI polarimetric observations of the benchmark T5.5 companion HD 19467 B. We detect no polarization signal from the target, and place an upper limit on the degree of linear polarization of . We discuss our results in the context of T dwarf cloud models and photometric variability. [ABSTRACT FROM AUTHOR]

Rotational modulations of brown dwarfs have recently provided powerful constraints on the properties of ultra-cool atmospheres, including longitudinal and vertical cloud structures and cloud evolution. Furthermore, periodic light curves directly probe the rotational periods of ultra-cool objects. We present here, for the first time, time-resolved high-precision photometric measurements of a planetary-mass companion, 2M1207b. We observed the binary system with Hubble Space Telescope/Wide Field Camera 3 in two bands and with two spacecraft roll angles. Using point-spread function-based photometry, we reach a nearly photon-noise limited accuracy for both the primary and the secondary. While the primary is consistent with a flat light curve, the secondary shows modulations that are clearly detected in the combined light curve as well as in different subsets of the data. The amplitudes are 1.36% in the F125W and 0.78% in the F160W filters, respectively. By fitting sine waves to the light curves, we find a consistent period of hr and similar phases in both bands. The J- and H-band amplitude ratio of 2M1207b is very similar to a field brown dwarf that has identical spectral type but different J–H color. Importantly, our study also measures, for the first time, the rotation period for a directly imaged extra-solar planetary-mass companion. [ABSTRACT FROM AUTHOR]

Observations of Neptune with the Kepler Space Telescope yield a 49 day light curve with 98% coverage at a 1 minute cadence. A significant signature in the light curve comes from discrete cloud features. We compare results extracted from the light curve data with contemporaneous disk-resolved imaging of Neptune from the Keck 10-m telescope at 1.65 microns and Hubble Space Telescope visible imaging acquired nine months later. This direct comparison validates the feature latitudes assigned to the K2 light curve periods based on Neptune's zonal wind profile, and confirms observed cloud feature variability. Although Neptune's clouds vary in location and intensity on short and long timescales, a single large discrete storm seen in Keck imaging dominates the K2 and Hubble light curves; smaller or fainter clouds likely contribute to short-term brightness variability. The K2 Neptune light curve, in conjunction with our imaging data, provides context for the interpretation of current and future brown dwarf and extrasolar planet variability measurements. In particular we suggest that the balance between large, relatively stable, atmospheric features and smaller, more transient, clouds controls the character of substellar atmospheric variability. Atmospheres dominated by a few large spots may show inherently greater light curve stability than those which exhibit a greater number of smaller features. [ABSTRACT FROM AUTHOR]

Several concepts now exist for small, space-based missions to directly characterize exoplanets in reflected light. While studies have been performed that investigate the potential detection yields of such missions, little work has been done to understand how instrumental and astrophysical parameters will affect the ability of these missions to obtain spectra that are useful for characterizing their planetary targets. Here, we develop an instrument noise model suitable for studying the spectral characterization potential of a coronagraph-equipped, space-based telescope. We adopt a baseline set of telescope and instrument parameters appropriate for near-future planned missions like WFIRST-AFTA, including a 2 m diameter primary aperture, an operational wavelength range of 0.4–1.0 μm, and an instrument spectral resolution of λ/Δλ = 70, and apply our baseline model to a variety of spectral models of different planet types, including Earth twins, Jupiter twins, and warm and cool Jupiters and Neptunes. With our exoplanet spectral models, we explore wavelength-dependent planet–star flux ratios for main-sequence stars of various effective temperatures and discuss how coronagraph inner and outer working angle constraints will influence the potential to study different types of planets. For planets most favorable to spectroscopic characterization—cool Jupiters and Neptunes as well as nearby super-Earths—we study the integration times required to achieve moderate signal-to-noise ratio spectra. We also explore the sensitivity of the integration times required to either detect the bottom or presence of key absorption bands (for methane, water vapor, and molecular oxygen) to coronagraph raw contrast performance, exozodiacal light levels, and the distance to the planetary system. Decreasing detector quantum efficiency at longer visible wavelengths makes the detection of water vapor in the atmospheres of Earth-like planets extremely challenging, and also hinders detections of the 0.89 μm methane band. Additionally, most modeled observations have noise dominated by dark currents, indicating that improving CCD performance could substantially drive down requisite integration times. Finally, we briefly discuss the extension of our models to a more distant future Large UV-Optical-InfraRed (LUVOIR) mission. [ABSTRACT FROM AUTHOR]

Planets larger than Earth and smaller than Neptune are some of the most numerous in the galaxy, but observational efforts to understand this population have proved challenging because optically thick clouds or hazes at high altitudes obscure molecular features. We present models of super Earths that include thick clouds and hazes and predict their transmission, thermal emission, and reflected light spectra. Very thick, lofted clouds of salts or sulfides in high metallicity (1000× solar) atmospheres create featureless transmission spectra in the near-infrared. Photochemical hazes with a range of particle sizes also create featureless transmission spectra at lower metallicities. Cloudy thermal emission spectra have muted features more like blackbodies, and hazy thermal emission spectra have emission features caused by an inversion layer at altitudes where the haze forms. Close analysis of reflected light from warm (∼400–800 K) planets can distinguish cloudy spectra, which have moderate albedos (0.05–0.20), from hazy models, which are very dark (0.0–0.03). Reflected light spectra of cold planets (∼200 K) accessible to a space-based visible light coronagraph will have high albedos and large molecular features that will allow them to be more easily characterized than the warmer transiting planets. We suggest a number of complementary observations to characterize this population of planets, including transmission spectra of hot ( K) targets, thermal emission spectra of warm targets using the James Webb Space Telescope, high spectral resolution (R ∼ 105) observations of cloudy targets, and reflected light spectral observations of directly imaged cold targets. Despite the dearth of features observed in super Earth transmission spectra to date, different observations will provide rich diagnostics of their atmospheres. [ABSTRACT FROM AUTHOR]

The re-emergence of the 0.99 μm FeH feature in brown dwarfs of early- to mid-T spectral type has been suggested as evidence for cloud disruption where flux from deep, hot regions below the Fe cloud deck can emerge. The same mechanism could account for color changes at the L/T transition and photometric variability. We present the first observations of spectroscopic variability of brown dwarfs covering the 0.99 μm FeH feature. We observed the spatially resolved very nearby brown dwarf binary WISE J104915.57–531906.1 (Luhman 16AB), a late-L and early-T dwarf, with Hubble Space Telescope/WFC3 in the G102 grism at 0.8–1.15 μm. We find significant variability at all wavelengths for both brown dwarfs, with peak-to-valley amplitudes of 9.3% for Luhman 16B and 4.5% for Luhman 16A. This represents the first unambiguous detection of variability in Luhman 16A. We estimate a rotational period between 4.5 and 5.5 hr, very similar to Luhman 16B. Variability in both components complicates the interpretation of spatially unresolved observations. The probability for finding large amplitude variability in any two brown dwarfs is less than 10%. Our finding may suggest that a common but yet unknown feature of the binary is important for the occurrence of variability. For both objects, the amplitude is nearly constant at all wavelengths except in the deep K i feature below 0.84 μm. No variations are seen across the 0.99 μm FeH feature. The observations lend strong further support to cloud height variations rather than holes in the silicate clouds, but cannot fully rule out holes in the iron clouds. We re-evaluate the diagnostic potential of the FeH feature as a tracer of cloud patchiness. [ABSTRACT FROM AUTHOR]

Interpreting the spectra of brown dwarfs is key to determining the fundamental physical and chemical processes occurring in their atmospheres. Powerful Bayesian atmospheric retrieval tools have recently been applied to both exoplanet and brown dwarf spectra to tease out the thermal structures and molecular abundances to understand those processes. In this manuscript we develop a significantly upgraded retrieval method and apply it to the SpeX spectral library data of two benchmark late T dwarfs, Gl 570D and HD 3651B, to establish the validity of our upgraded forward model parameterization and Bayesian estimator. Our retrieved metallicities, gravities, and effective temperatures are consistent with the metallicity and presumed ages of the systems. We add the carbon-to-oxygen ratio as a new dimension to benchmark systems and find good agreement between carbon-to-oxygen ratios derived in the brown dwarfs and the host stars. Furthermore, we have for the first time unambiguously determined the presence of ammonia in the low-resolution spectra of these two late T dwarfs. We also show that the retrieved results are not significantly impacted by the possible presence of clouds, though some quantities are significantly impacted by uncertainties in photometry. This investigation represents a watershed study in establishing the utility of atmospheric retrieval approaches on brown dwarf spectra. [ABSTRACT FROM AUTHOR]