Search

The spectra of brown dwarfs are key to exploring the chemistry and physics that take place in their atmospheres. Late-T dwarf spectra are particularly diagnostic due to their relatively cloud-free atmospheres and deep molecular bands. With the use of powerful atmospheric retrieval tools applied to the spectra of these objects, direct constraints on molecular/atomic abundances, gravity, and vertical thermal profiles can be obtained enabling a broad exploration of the chemical/physical mechanisms operating in their atmospheres. We present a uniform retrieval analysis on low-resolution IRTF SpeX near-IR spectra of a sample of 50 T dwarfs, including new observations as part of a recent volume-limited survey. This analysis more than quadruples the sample of T dwarfs with retrieved temperature profiles and abundances (H$_2$O, CH$_4$, NH$_3$, K and subsequent C/O and metallicities). We are generally able to constrain effective temperatures to within 50K, volume mixing ratios for major species to within 0.25dex, atmospheric metallicities [M/H] to within 0.2, and C/O ratios to within 0.2. We compare our retrieved constraints on the thermal structure, chemistry, and gravities of these objects with predictions from self-consistent radiative-convective equilibrium models and find, in general though with substantial scatter, consistency with solar composition chemistry and thermal profiles of the neighboring stellar FGK population. Objects with notable discrepancies between the two modeling techniques and potential mechanisms for their differences, be they related to modeling approach or physically motivated, are discussed more thoroughly in the text., 21 pages, 11 figures

The detection of disk-integrated polarization from Luhman 16A and B in H-band, and subsequent modeling, has been interpreted in the framework of zonal cloud bands on these bodies. Recently, Tan and Showman (2021) investigated three-dimensional atmospheric circulation and cloud structures of brown dwarfs with general circulation models (GCMs), and their simulations yield complex cloud distributions showing some aspects of zonal jets, but also complex vortices that cannot be captured by a simple model. Here we use these 3D GCMs specific to Luhman 16A and B, along with the three-dimensional Monte Carlo radiative transfer code ARTES, to calculate their polarization signals. We adopt the 3D temperature-pressure and cloud profiles from the GCMs as our input atmospheric structures. Our polarization calculations at 1.6 $\mu$m agree well with the measured degree of linear polarization from both Luhman 16 A and B. Our calculations reproduce the measured polarization for both the objects with cloud particle sizes between 0.5-1 \,$\mu$m for Luhman 16 A and 5 \,$\mu$m for Luhman 16 B. We find that the degree of linear polarization can vary on hour-long timescales over the course of a rotation period. We also show that models with azimuthally symmetric band-like cloud geometries, typically used for interpreting polarimetry observations of brown dwarfs, over-predict the polarization signal if the cloud patterns do not include complex vortices within these bands. This exploratory work shows that GCMs are promising for modeling and interpreting polarization signals of brown dwarfs., Comment: 13 pages, 6 figures, Accepted for publication in The Astrophysical journal

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

We present $Spitzer$ IRS 5--14 $\mu$m spectra and 16 $\mu$m and 22 $\mu$m photometry of the T2.5 companion to the $\sim$300 Myr-old G0V star HN Peg. We incorporate previous 0.8--5 $\mu$m observations to obtain the most comprehensive spectral energy distribution of an intermediate-gravity L/T-transition dwarf which, together with an accurate Gaia EDR3 parallax of the primary, enable us to derive precise fundamental parameters. We find that young ($\approx$0.1--0.3 Gyr) early-T dwarfs on average have $\approx$140 K lower effective temperatures, $\approx$20% larger radii, and similar bolometric luminosities compared to $\gtrsim$1 Gyr-old field dwarfs with similar spectral types. Our accurate infrared spectrophotometry offers new detail at wavelengths where the dominant carbon-bearing molecules have their strongest transitions: at 3.4 $\mu$m for methane and at 4.6 $\mu$m for carbon monoxide. We assess the performance of various widely available photospheric models and find that models with condensates and/or clouds better reproduce the full SED of this moderately young early-T dwarf. However, cloud-free models incorporating a more general convective instability treatment reproduce at least the low-resolution near-IR spectrum similarly well. Our analysis of $R\approx2300$ $J$-band spectra shows that the near-infrared potassium absorption lines in HN Peg B have similar strengths to those seen in both younger and older T2-T3 dwarfs. We conclude that while alkali lines are well-established as surface gravity indicators for L-type or warmer stars, they are insensitive to surface gravity in early-T dwarfs, Comment: Accepted for publication in ApJ

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

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.

We use the Wide Field Camera 3 on the {\sl Hubble Space Telescope} to spectrophotometrically monitor the young L7.5 companion HD~203030B. Our time series reveal photometric variability at 1.27\,$\mu$m and 1.39\,$\mu$m on time scales compatible with rotation. We find a rotation period of $7.5^{+0.6}_{-0.5}$ h: comparable to those observed in other brown dwarfs and planetary-mass companions younger than 300 Myr. We measure variability amplitudes of $1.1\pm0.3\%$ (1.27\,$\mu$m) and $1.7\pm0.4\%$ (1.39\,$\mu$m), and a phase lag of 56$^\circ\pm$28$^\circ$ between the two light curves. We attribute the difference in photometric amplitudes and phases to a patchy cloud layer that is sinking below the level where water vapor becomes opaque. HD 203030B and the few other known variable young late-L dwarfs are unlike warmer (earlier-type and/or older) L dwarfs, for which variability is much less wavelength-dependent across the 1.1--1.7$\mu$m region. We further suggest that a sinking of the top-most cloud deck below the level where water or carbon monoxide gas become opaque may also explain the often enhanced variability amplitudes of even earlier-type low-gravity L dwarfs. Because these condensate and gas opacity levels are already well-differentiated in T dwarfs, we do not expect the same variability amplitude enhancement in young vs.\ old T dwarfs., Comment: Accepted for publication in The Astrophysical Journal

Measurements of photometric variability at different wavelengths provide insights into the vertical cloud structure of brown dwarfs and planetary-mass objects. In seven Hubble Space Telescope consecutive orbits, spanning $\sim$10 h of observing time}, we obtained time-resolved spectroscopy of the planetary-mass T8-dwarf Ross 458C using the near-infrared Wide Field Camera 3. We found spectrophotometric variability with a peak-to-peak signal of 2.62$\pm$0.02 % (in the 1.10-1.60~$\mu$m white light curve). Using three different methods, we estimated a rotational period of 6.75$\pm$1.58~h for the white light curve, and similar periods for narrow $J$- and $H$- band light curves. Sine wave fits to the narrow $J$- and $H$-band light curves suggest a tentative phase shift between the light curves with wavelength when we allow different periods between both light curves. If confirmed, this phase shift may be similar to the phase shift detected earlier for the T6.5 spectral type 2MASS J22282889-310262. We find that, in contrast with 2M2228, the variability of Ross~458C shows evidence for a {color trend} within the narrow $J$-band, but gray variations in the narrow $H$-band. The spectral time-resolved variability of Ross 458C might be potentially due to heterogeneous sulfide clouds in the atmosphere of the object. Our discovery extends the study of spectral modulations of condensate clouds to the coolest T dwarfs, planetary-mass companions., Comment: Accepted in ApJL

First generation optical coronagraphic telescopes will obtain images of cool gas and ice giant exoplanets around nearby stars. The albedo spectra of exoplanets at planet-star separations larger than about 1 AU are dominated by reflected light to beyond 1 {\mu}m and are punctuated by molecular absorption features. We consider how exoplanet albedo spectra and colors vary as a function of planet-star separation, metallicity, mass, and observed phase for Jupiter and Neptune analogs from 0.35 to 1 {\mu}m. We model Jupiter analogs with 1x and 3x the solar abundance of heavy elements, and Neptune analogs with 10x and 30x. Our model planets orbit a solar analog parent star at separations of 0.8 AU, 2 AU, 5 AU, and 10 AU. We use a radiative-convective model to compute temperature-pressure profiles. The giant exoplanets are cloud-free at 0.8 AU, have H2O clouds at 2 AU, and have both NH3 and H2O clouds at 5 AU and 10 AU. For each model planet we compute moderate resolution spectra as a function of phase. The presence and structure of clouds strongly influence the spectra. Since the planet images will be unresolved, their phase may not be obvious, and multiple observations will be needed to discriminate between the effects of planet-star separation, metallicity, and phase. We consider the range of these combined effects on spectra and colors. For example, we find that the spectral influence of clouds depends more on planet-star separation and hence temperature than metallicity, and it is easier to discriminate between cloudy 1x and 3x Jupiters than between 10x and 30x Neptunes. In addition to alkalis and methane, our Jupiter models show H2O absorption features near 0.94 {\mu}m. We also predict that giant exoplanets receiving greater insolation than Jupiter will exhibit higher equator to pole temperature gradients than are found on Jupiter and thus may have differing atmospheric dynamics., Comment: 62 pages, 19 figures, 6 tables Accepted for publication in ApJ

Light scattering by atmospheric dust particles is responsible for the polarization observed in some L dwarfs. Whether this polarization arises from an inhomogeneous distribution of dust across the disk or an oblate shape induced by rotation remains unclear. Here we argue that the latter case is plausible and, for many L dwarfs, the more likely one. Furthermore evolutionary models of mature field L dwarfs predict surface gravities ranging from about 200 to 2500 m/s^2 (corresponding to masses of about 15 to 70 times of Jupiter mass). Yet comparison of observed spectra to available synthetic spectra often does not permit more precise determination of the surface gravity of individual field L dwarfs, leading to important uncertainties in their properties. Since rotationally-induced non-sphericity, which gives rise to non-zero disk-integrated polarization, is more pronounced at lower gravities, polarization is a promising low gravity indicator. Here we combine a rigorous multiple scattering analysis with a self-consistent cloudy atmospheric model and observationally inferred rotational velocities and find that the observed optical polarization can be explained if the surface gravity of the polarized objects is about 300 m/s^2 or less, potentially providing a new method for constraining L dwarf masses., 14 pages (AASTEX) including 3 eps figures. Accepted for publication in The Astrophysical Journal Letters

We present 7.5-14.2um low-resolution spectroscopy, obtained with the Spitzer Infrared Spectrograph, of the T8.5 dwarf Wolf 940 B, which is a companion to an M4 dwarf with a projected separation of 400 AU. We combine these data with previously published near-infrared spectroscopy and mid-infrared photometry, to produce the spectral energy distribution for the very low-temperature T dwarf. We use atmospheric models to derive the bolometric correction and obtain a luminosity of log L/Lsun = -6.01 +/- 0.05. Evolutionary models are used with the luminosity to constrain the values of effective temperature (T_eff) and surface gravity, and hence mass and age for the T dwarf. We further restrict the allowed range of T_eff and gravity using age constraints implied by the M dwarf primary, and refine the physical properties of the T dwarf by comparison of the observed and modelled spectroscopy and photometry. This comparison indicates that Wolf 940 B has a metallicity within 0.2 dex of solar, as more extreme values give poor fits to the data - lower metallicity produces a poor fit at lambda > 2um while higher metallicity produces a poor fit at lambda < 2um. This is consistent with the independently derived value of [m/H] = +0.24 +/- 0.09 for the primary star, using the Johnson & Apps (2008) M_K:V-K relationship. We find that the T dwarf atmosphere is undergoing vigorous mixing, with an eddy diffusion coefficient K_zz of 10^4 to 10^6 cm^2 s^-1. We derive an effective temperature of 585 K to 625 K, and surface gravity log g = 4.83 to 5.22 (cm s^-2), for an age range of 3 Gyr to 10 Gyr, as implied by the kinematic and H alpha properties of the M dwarf primary. The lower gravity corresponds to the lower temperature and younger age for the system, and the higher value to the higher temperature and older age. The mass of the T dwarf is 24 M_Jupiter to 45 M_Jupiter for the younger to older age limit., 24 pages which include 5 Figures and 3 Tables. Accepted for publication in the Astrophysical Journal July 2 2010

While there have been multiple observational programs aimed at detecting linear polarization of optical radiation emitted by ultracool dwarfs, there has been comparatively less rigorous theoretical analysis of the problem. The general expectation has been that the atmospheres of those substellar-mass objects with condensate clouds would give rise to linear polarization due to scattering. Because of rotation-induced non-sphericity, there is expected to be incomplete cancellation of disk-integrated net polarization and thus a finite polarization. For cloudless objects, however, only molecular Rayleigh scattering will contribute to any net polarization and this limit has not been well studied. Hence in this paper we present a detailed multiple scattering analysis of the polarization expected from those T-dwarfs whose spectra show absence of condensates. For this, we develop and solve the full radiative transfer equations for linearly polarized radiation. Only atomic and molecular Rayleigh scattering are considered to be the source of polarization. We compute the local polarization at different angular directions in a plane-parallel atmospheres calculated for the range of effective temperatures of T dwarfs and then average over the whole surface of the object. The effects of gravity and limb darkening as well as rotation induced non-sphericity are included. It is found that the amount of polarization decreases with the increase in effective temperature. It is also found that significant polarization at any local point in the atmosphere arises only in the optical (B-band). However, the disk integrated polarization--even in the B-band--is negligible. Hence we conclude that, unlike the case for cloudy L dwarfs, polarization of cloudless T-dwarfs by atomic and molecular scattering may not be detectable., Comment: 31 pages (aastex) including 11 figures (eps), accepted for publication in the Astrophysical Journal

Determinations of the luminosity and atmospheric properties of the T8 brown dwarf 2MASS J09393548 - 2448279 are presented, based on Spitzer IRAC and IRS observations and ground-based astrometry. We find log_(10)(L_(bol)/L⊙) = -5.69 ± 0.03 for this source, comparable to the current low-luminosity record holder 2MASS J04151954 - 0935066. However, modeling of near- and mid-infrared spectral data indicates an effective temperature of 600 ± 35 K, roughly 100 K cooler than 2M0415. These parameters require a highly inflated radius for 2M0939 (R ≈ 0.13 R⊙) which cannot be reconciled with brown dwarf structure models. However, if this source is an unresolved, equal-mass binary, then the reduced luminosity of each component (L_(bol) ≈ 10^-6 L⊙) can be brought into agreement with the inferred atmospheric parameters for an age of 0.4-12 Gyr and component masses of 0.01-0.05 M⊙. This hypothesis can be tested through future high-resolution imaging and/or spectroscopic observations.

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

We present 10 - 19 um moderate resolution spectra of ten M dwarfs, one L dwarf, and two T dwarf systems obtained with the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope. The IRS allows us to examine molecular spectroscopic features/lines at moderate spectral resolution in a heretofore untapped wavelength regime. These R~600 spectra allow for a more detailed examination of clouds, non-equilibrium chemistry, as well as the molecular features of H2O, NH3, and other trace molecular species that are the hallmarks of these objects. A cloud-free model best fits our mid-infrared spectrum of the T1 dwarf epsilon Indi Ba, and we find that the NH3 feature in epsilon Indi Bb is best explained by a non-equilibrium abundance due to vertical transport in its atmosphere. We examined a set of objects (mostly M dwarfs) in multiple systems to look for evidence of emission features, which might indicate an atmospheric temperature inversion, as well as trace molecular species; however, we found no evidence of either., Comment: 19 pages, 7 figures, accepted ApJ 1/12/07

We present new L' (3.75 um) photometry of six L and T dwarfs, and M' (4.70 um) photometry of ten L and T dwarfs, observed at Gemini Observatory, and new 3.55, 4.49, 5.73 and 7.87 um photometry of nine L and T dwarfs, obtained with the Spitzer Space Telescope. The sample includes unusually blue and red dwarfs from our near-infrared studies. The data are combined with published L', M' and Spitzer photometry of L and T dwarfs, and trends of colors with spectral type and other colors are examined. Model atmospheres by Marley and Saumon are used to generate synthetic colors for ranges of effective temperature, gravity, grain sedimentation efficiency, metallicity and vertical mixing efficiency. We explore how these parameters affect the mid-infrared colors of L and T dwarfs and find that the data are modelled satisfactorily only if substantial vertical mixing occurs in both L- and T-dwarf atmospheres. The location and range of the L and T dwarf sequences in IRAC color-color and color - magnitude diagrams is also only reproduced if this mixing occurs, with a range of efficiency described by K_zz ~ 10^2-10^6 cm^2/s. The colors of the unusually red dwarfs are best reproduced by non-equilibrium models with low sedimentation efficiency, i.e. thick cloud decks, and those of the unusually blue dwarfs by non-equilibrium models with high sedimentation efficiency, i.e. thin cloud decks. The K-L' and Spitzer [3.55]-[4.49] colors can be used as indicators of effective temperature for L and T dwarfs, but care must be taken to include gravity and metallicity effects for late-T dwarfs and vertical mixing for both late-L and T dwarfs., Accepted for publication in the Astrophysical Journal; 21 pages, 8 Figures, 3 Tables

Traditional thermal evolution models of giant planets employ arbitrary initial conditions selected more for computational expediency than physical accuracy. Since the initial conditions are eventually forgotten by the evolving planet, this approach is valid for mature planets, if not young ones. To explore the evolution at young ages of jovian mass planets we have employed model planets created by one implementation of the core accretion mechanism as initial conditions for evolutionary calculations. We find that young jovian planets are smaller, cooler, and several to 100 times less luminous than predicted by earlier models. Furthermore the time interval during which the young jupiters are fainter than expected depends on the mass of planet. Jupiter mass planets (1 M_J) align with the conventional model luminosity in as little at 20 million years, but 10 M_J planets can take up to 1 billion years to match commonly cited luminosities, given our implementation of the core accretion mechanism. If our assumptions, especially including our treatment of the accretion shock, are correct, then young jovian planets are substantially fainter at young ages than currently believed. These results have important consequences both for detection strategies and for assigning masses to young jovian planets based on observed luminosities., 20 pages, 6 figures, Astrophys. J. in press

We present the first mid-infrared spectra of brown dwarfs, together with observations of a low-mass star. Our targets are the M3.5 dwarf GJ 1001A, the L8 dwarf DENIS-P J0255-4700, and the T1/T6 binary system epsilon Indi Ba/Bb. As expected, the mid-infrared spectral morphology of these objects changes rapidly with spectral class due to the changes in atmospheric chemistry resulting from their differing effective temperatures and atmospheric structures. By taking advantage of the unprecedented sensitivity of the Infrared Spectrograph on the Spitzer Space Telescope we have detected the 7.8 micron methane and 10 micron ammonia bands for the first time in brown dwarf spectra.

We present the first analysis of an optical to mid-infrared spectrum of the T7.5 dwarf Gliese 570D with model atmospheres, synthetic spectra, and brown dwarf evolution sequences. We obtain precise values for the basic parameters of Gl 570D: Teff=800 - 820K, log g (cm/s^2)=5.09 - 5.23, and log L/Lsun= -5.525 to -5.551. The Spitzer IRS spectrum shows prominent features of ammonia (NH3) that can only be fitted by reducing the abundance of NH3 by about one order of magnitude from the value obtained with chemical equilibrium models. We model departures from chemical equilibrium in the atmosphere of Gl 570D by considering the kinetics of nitrogen and carbon chemistry in the presence of vertical mixing. The resulting model spectrum reproduces the data very well., Comment: Accepted for publication in the ApJ. 10 pages, including 3 figures

We investigate the atmosphere and interior of the new transiting planet HD 149026b, which appears to be very rich in heavy elements. We first compute model atmospheres at metallicities ranging from solar to ten times solar, and show how for cases with high metallicity or inefficient redistribution of energy from the day side, the planet may develop a hot stratosphere due to absorption of stellar flux by TiO and VO. The spectra predicted by these models are very different than cooler atmosphere models without stratospheres. The spectral effects are potentially detectable with the Spitzer Space Telescope. In addition the models with hot stratospheres lead to a large limb brightening, rather than darkening. We compare the atmosphere of HD 149026b to other well-known transiting planets, including the recently discovered HD 189733b, which we show have planet-to-star flux ratios twice that of HD 209458 and TrES-1. The methane abundance in the atmosphere of HD 189733b is a sensitive indicator of atmospheric temperature and metallicity and can be constrained with Spitzer IRAC observations. We then turn to interior studies of HD 149026b and use a grid of self-consistent model atmospheres and high-pressure equations of state for all components to compute thermal evolution models of the planet. We estimate that the mass of heavy elements within the planet is in the range of 60 to 93 M_earth. Finally, we discuss trends in the radii of transiting planets with metallicity in light of this new member of the class., Comment: Accepted to the Astrophysical Journal. 18 pages, including 10 figures. New section on the atmosphere of planet HD 189733b. Enhanced discussion of atmospheric Ti chemistry and core mass for HD 149026b

We present K- and L-band photometry obtained with the Keck I telescope for a representative sample of L and T dwarfs. These observations were motivated in part by the dominant role H2O and CH4 play in shaping the flux near 2 and 3 μm and by the potential use of these bands as indicators of spectral class in the infrared. In addition, these observations aid the determination of the bolometric luminosity of L and T dwarfs. Here we report the K, L', and Ls magnitudes of our objects and the trends observed in the (K-L') and (K-Ls) colors as a function of L and T dwarf spectral class. We compare these colors with theoretical models, derive a relationship between effective temperature and L spectral class, and compare our temperature estimates with others.

We have detected weak absorption features produced by the strong nu(3) methane band at 3.3 microns in two L dwarfs, 2MASSW J1507476-162738 and 2MASSI J0825196+211552, classified by Kirkpatrick et al. (2000) as spectral types L5 and L7.5 respectively. These absorptions occur in objects warmer than any in which methane previously has been detected, and mark the first appearance of methane in the ultracool star-to-brown dwarf spectral sequence., Comment: 12 pages, 2 figures; accepted by ApJ Letters

We report the detection of carbon monoxide in the atmosphere of the brown dwarf Gliese 229B in spectra obtained at wavelengths near 4.7 μm. The spectrum is consistent with mole fractions of CO of more than 50 parts per million (ppm), much higher than expected for gas in thermochemical equilibrium at the temperature of the atmosphere we observe. As in Jupiter, mixing in the upper atmosphere of Gl 229B may prevent equilibrium between CO and CH4 at low temperatures. The average flux from 4.55 to 5.10 μm in Gl 229B is measured to be 3.6±0.2 mJy, lower than previous broadband measurements.

We have monitored photometrically the Y0 brown dwarf WISEP J173835.52+273258.9 (W1738) at both near- and mid-infrared wavelengths. This ~1 Gyr-old 400K dwarf is at a distance of 8pc and has a mass around 5 M_Jupiter. We observed W1738 using two near-infrared filters at lambda~1um, Y and J, on Gemini observatory, and two mid-infrared filters at lambda~4um, [3.6] and [4.5], on the Spitzer observatory. Twenty-four hours were spent on the source by Spitzer on each of June 30 and October 30 2013 UT. Between these observations, around 5 hours were spent on the source by Gemini on each of July 17 and August 23 2013 UT. The mid-infrared light curves show significant evolution between the two observations separated by four months. We find that a double sinusoid can be fit to the [4.5] data, where one sinusoid has a period of 6.0 +/- 0.1 hours and the other a period of 3.0 +/- 0.1 hours. The near-infrared observations suggest variability with a ~3.0 hour period, although only at a, Comment: Accepted by ApJ July 26 2016. Twenty-six pages include 8 Figures and 5 Tables

Condensate clouds fundamentally impact the atmospheric structure and spectra of exoplanets and brown dwarfs but the connections between surface gravity, cloud structure, dust in the upper atmosphere, and the red colors of some brown dwarfs remain poorly understood. Rotational modulations enable the study of different clouds in the same atmosphere, thereby providing a method to isolate the effects of clouds. Here we present the discovery of high peak-to-peak amplitude (8%) rotational modulations in a low-gravity, extremely red (J-Ks=2.55) L6 dwarf WISEP J004701.06+680352.1 (W0047). Using the Hubble Space Telescope (HST) time-resolved grism spectroscopy we find a best-fit rotational period (13.20$\pm$0.14 hours) with a larger amplitude at 1.1 micron than at 1.7 micron. This is the third largest near-infrared variability amplitude measured in a brown dwarf, demonstrating that large-amplitude variations are not limited to the L/T transition but are present in some extremely red L-type dwarfs. We report a tentative trend between the wavelength dependence of relative amplitude, possibly proxy for small dust grains lofted in the upper atmosphere, and the likelihood of large-amplitude variability. By assuming forsterite as haze particle, we successfully explain the wavelength dependent amplitude with submicron-sized haze particles sizes of around 0.4 {\mu}m. W0047 links the earlier spectral and later spectral type brown dwarfs in which rotational modulations have been observed, the large amplitude variations in this object make this a benchmark brown dwarf for the study of cloud properties close to the L/T transition., Comment: 9 pages, 4 figures, Accepted for publication in ApJL on September 15th

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

We use a 1D model to address photochemistry and possible haze formation in the irradiated warm Jupiter, 51 Eridani b. The intended focus was to be carbon, but sulfur photochemistry turns out to be important. The case for organic photochemical hazes is intriguing but falls short of being compelling. If organic hazes form, they are likeliest to do so if vertical mixing in 51 Eri b is weaker than in Jupiter, and they would be found below the altitudes where methane and water are photolyzed. The more novel result is that photochemistry turns H$_2$S into elemental sulfur, here treated as S$_8$. In the cooler models, S$_8$ is predicted to condense in optically thick clouds of solid sulfur particles, whilst in the warmer models S$_8$ remains a vapor along with several other sulfur allotropes that are both visually striking and potentially observable. For 51 Eri b, the division between models with and without condensed sulfur is at an effective temperature of 700 K, which is within error its actual effective temperature; the local temperature where sulfur condenses is between 280 and 320 K. The sulfur photochemistry we have discussed is quite general and ought to be found in a wide variety of worlds over a broad temperature range, both colder and hotter than the 650-750 K range studied here, and we show that products of sulfur photochemistry will be nearly as abundant on planets where the UV irradiation is orders of magnitude weaker than it is on 51 Eri b., Comment: 24 pages including 11 figures and a table

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

We present results from an atmospheric circulation study of nine hot Jupiters that comprise 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 which suggest 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 K, which could imply large variations in Na, K, CO and CH4 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 with 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 datasets, 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., Comment: 16 pages, 10 figures; accepted to the Astrophysical Journal

As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the first generation of directly-imaged exoplanets were all L-type. Recently, Kuzuhara et al. (2013) announced the discovery of GJ 504 b, the first T dwarf exoplanet. GJ 504 b provides a unique opportunity to study the atmosphere of a new type of exoplanet with a ~500 K temperature that bridges the gap between the first directly imaged planets (~1000 K) and our own Solar System's Jupiter (~130 K). We observed GJ 504 b in three narrow L-band filters (3.71, 3.88, and 4.00 microns), spanning the red end of the broad methane fundamental absorption feature (3.3 microns) as part of the LEECH exoplanet imaging survey. By comparing our new photometry and literature photometry to a grid of custom model atmospheres, we were able to fit GJ 504 b's unusual spectral energy distribution for the first time. We find that GJ 504 b is well-fit by models with the following parameters: T_eff=544+/-10 K, g, resubmitted to ApJ with referee's comments

Author(s): Morley, CV; Fortney, JJ; Marley, MS; Zahnle, K; Line, M; Kempton, E; Lewis, N; Cahoy, K | Abstract: 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 (≳1000 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.

The re-emergence of the 0.99 $\mu$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 $\mu$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 HST/WFC3 in the G102 grism at 0.8-1.15 $\mu$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 h, 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 $\mu$m. No variations are seen across the 0.99 $\mu$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., Comment: 11 pages, 7 figures. Accepted for publication in the Astrophysical Journal

We use a planetary albedo model to investigate variations in visible wavelength phase curves of exoplanets. Thermal and cloud properties for these exoplanets are derived using one-dimensional radiative-convective and cloud simulations. The presence of clouds on these exoplanets significantly alters their planetary albedo spectra. We confirm that non-uniform cloud coverage on the dayside of tidally locked exoplanets will manifest as changes to the magnitude and shift of the phase curve. In this work, we first investigate a test case of our model using a Jupiter-like planet, at temperatures consistent to 2.0 AU insolation from a solar type star, to consider the effect of H2O clouds. We then extend our application of the model to the exoplanet Kepler-7b and consider the effect of varying cloud species, sedimentation efficiency, particle size, and cloud altitude. We show that, depending on the observational filter, the largest possible shift of the phase curve maximum will be ~2°–10° for a Jupiter-like planet, and up to ~30° (~0.08 in fractional orbital phase) for hot-Jupiter exoplanets at visible wavelengths as a function of dayside cloud distribution with a uniformly averaged thermal profile. The models presented in this work can be adapted for a variety of planetary cases at visible wavelengths to include variations in planet–star separation, gravity, metallicity, and source-observer geometry. Finally, we tailor our model for comparison with, and confirmation of, the recent optical phase-curve observations of Kepler-7b with the Kepler space telescope. The average planetary albedo can vary between 0.1 and 0.6 for the 1300 cloud scenarios that were compared to the observations. Many of these cases cannot produce a high enough albedo to match the observations. We observe that smaller particle size and increasing cloud altitude have a strong effect on increasing albedo. In particular, we show that a set of models where Kepler-7b has roughly half of its dayside covered in small-particle clouds high in the atmosphere, made of bright minerals like MgSiO3 and Mg2SiO4, provide the best fits to the observed offset and magnitude of the phase-curve, whereas Fe clouds are found to be too dark to fit the observations.

We address disequilibrum abundances of some simple molecules in the atmospheres of solar composition brown dwarfs and self-luminous extrasolar giant planets using a kinetics-based 1D atmospheric chemistry model. Our approach is to use the full kinetics model to survey the parameter space with effective temperatures between 500 K and 1100 K. In all of these worlds equilibrium chemistry favors CH4 over CO in the parts of the atmosphere that can be seen from Earth, but in most disequilibrium favors CO. The small surface gravity of a planet strongly discriminates against CH4 when compared to an otherwise comparable brown dwarf. If vertical mixing is like Jupiter's, the transition from methane to CO occurs at 500 K in a planet. Sluggish vertical mixing can raise this to 600 K; but clouds or more vigorous vertical mixing could lower this to 400 K. The comparable thresholds in brown dwarfs are $1100\pm100$ K. Ammonia is also sensitive to gravity, but unlike CH4/CO, the NH3/N2 ratio is insensitive to mixing, which makes NH3 a potential proxy for gravity. HCN may become interesting in high gravity brown dwarfs with very strong vertical mixing. Detailed analysis of the CO-CH4 reaction network reveals that the bottleneck to CO hydrogenation goes through methanol, in partial agreement with previous work. Simple, easy to use quenching relations are derived by fitting to the complete chemistry of the full ensemble of models. These relations are valid for determining CO, CH4, NH3, HCN, and CO2 abundances in the range of self-luminous worlds we have studied but may not apply if atmospheres are strongly heated at high altitudes by processes not considered here (e.g., wave breaking)., Comment: Astrophysical Journal, in press. Clarity improvements throughout and one new figure. 17 figures, 20 pages

We present 0.9 - 2.5 um resolved spectra for the ultracool binary WISEPC J121756.91+162640.2AB. The system consists of a pair of brown dwarfs that straddles the currently defined T/Y spectral type boundary. We use synthetic spectra generated by model atmospheres that include chloride and sulfide clouds (Morley et al.), the distance to the system (Dupuy & Kraus), and the radius of each component based on evolutionary models (Saumon & Marley) to determine a probable range of physical properties for the binary. The effective temperature of the T8.5 primary is 550 - 600 K, and that of the Y0 - Y0.5 secondary is 450 K. The atmospheres of both components are either free of clouds or have extremely thin cloud layers. We find that the masses of the primary and secondary are 30 and 22 M_Jup, respectively, and that the age of the system is 4 - 8 Gyr. This age is consistent with astrometric measurements (Dupuy & Kraus) that show that the system has kinematics intermediate between those of the thin and thick disks of the Galaxy. An older age is also consistent with an indication by the H - K colors that the system is slightly metal-poor., 21 pages which include 6 Figures and 3 Tables. Accepted on November 8 2013 for publication in ApJ

The recent discovery that the close-in extrasolar giant planet, HD209458b, transits its star has provided a first-of-its-kind measurement of the planet's radius and mass. In addition, there is a provocative detection of the light reflected off of the giant planet, $\tau$ Boo b. Including the effects of stellar irradiation, we estimate the general behavior of radius/age trajectories for such planets and interpret the large measured radii of HD209458b and $\tau$ Boo b in that context. We find that HD209458b must be a hydrogen-rich gas giant. Furthermore, the large radius of close-in gas giant is not due to the thermal expansion of its atmosphere, but to the high residual entropy that remains throughout its bulk by dint of its early proximity to a luminous primary. The large stellar flux does not inflate the planet, but retards its otherwise inexorable contraction from a more extended configuration at birth. This implies either that such a planet was formed near its current orbital distance or that it migrated in from larger distances ($\geq$0.5 A.U.), no later than a few times $10^7$ years of birth., Comment: aasms4 LaTeX, 1 figure, accepted to Ap.J. Letters