BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Mark S. Marley

Mark S. Marley appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.IMastro-phastro-ph.GAastro-ph.CO

Catalog footprint

What is connected

87works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

87 published item(s)

preprint2022arXiv

PICASO 3.0: A One-Dimensional Climate Model for Giant Planets and Brown Dwarfs

Upcoming James Webb Space Telescope (JWST) observations will allow us to study exoplanet and brown dwarf atmospheres in great detail. The physical interpretation of these upcoming high signal-to-noise observations requires precise atmospheric models of exoplanets and brown dwarfs. While several one-dimensional and three-dimensional atmospheric models have been developed in the past three decades, these models have often relied on simplified assumptions like chemical equilibrium and are also often not open-source, which limits their usage and development by the wider community. We present a python-based one-dimensional atmospheric radiative-convective equilibrium model. This model has heritage from the Fortran-based code (Marley et al.,1996} which has been widely used to model the atmospheres of Solar System objects, brown dwarfs, and exoplanets. In short, the basic capability of the original model is to compute the atmospheric state of the object under radiative-convective equilibrium given its effective or internal temperature, gravity, and host--star properties (if relevant). In the new model, which has been included within the well-utilized code-base PICASO, we have added these original features as well as the new capability of self-consistently treating disequilibrium chemistry. This code is widely applicable to Hydrogen-dominated atmospheres (e.g., brown dwarfs and giant planets).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Polarization of Rotationally Oblate Self-Luminous Exoplanets with Anisotropic Atmospheres

Young self-luminous giant exoplanets are expected to be oblate in shape owing to the high rotational speeds observed for some objects. Similar to the case of brown dwarfs, the thermal emission from these planets should be polarized by scatterings of molecules and condensate cloud particles, and the rotation-induced asymmetry of the planet's disk would yield to net non-zero detectable polarization. Considering an anisotropic atmosphere, we present here a three-dimensional approach to estimate the disk-averaged polarization that arises due to the oblateness of the planets. We solve the multiple-scattering vector radiative transfer equations at each location on the planet's disk and calculate the local Stokes vectors and then calculate the disk-integrated flux and linear polarization. For a cloud-free atmosphere, the polarization signal is observable only in the visible wavelength region. However, the presence of clouds in the planetary atmospheres leads to a detectable amount of polarization in the infrared wavelength region where the planetary thermal emission peaks. Considering different broad-band filters of the SPHERE-IRDIS instrument of the Very Large Telescope, we present generic models for the polarization at different wavelength bands as a function of their rotation period. We also present polarization models for the Exoplanets $β$ Pic b and ROXs 42B b as two representative cases which can guide future observations. Our insights on the polarization of young giant planets presented here would be useful for the upcoming polarimetric observations of the directly imaged planets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A New Sedimentation Model for Greater Cloud Diversity in Giant Exoplanets and Brown Dwarfs

The observed atmospheric spectrum of exoplanets and brown dwarfs depends critically on the presence and distribution of atmospheric condensates. The Ackerman & Marley (2001) methodology for predicting the vertical distribution of condensate particles is widely used to study cloudy atmospheres and has recently been implemented in an open-source python package virga. The model relies upon input parameter $f_{\text{sed}}$, the sedimentation efficiency, which until now has been held constant. The relative simplicity of this model renders it useful for retrieval studies due to its rapidly attainable solutions. However, comparisons with more complex microphysical models such as CARMA have highlighted inconsistencies between the two approaches, namely that the cloud parameters needed for radiative transfer produced by virga are dissimilar to those produced by CARMA. To address these discrepancies, we have extended the original Ackerman and Marley methodology in virga to allow for non-constant $f_{\text{sed}}$ values, in particular those that vary with altitude. We discuss one such parameterization and compare the cloud mass mixing ratio produced by virga with constant and variable $f_{\text{sed}}$ profiles to that produced by CARMA. We find that the variable $f_{\text{sed}}$ formulation better captures the profile produced by CARMA with heterogeneous nucleation, yet performs comparatively to constant $f_{\text{sed}}$ for homogeneous nucleation. In general, virga has the capacity to handle any $f_{\text{sed}}$ with an explicit anti-derivative, permitting a plethora of alternative cloud profiles that are otherwise unattainable by constant $f_{\text{sed}}$ values. The ensuing flexibility has the potential to better agree with increasingly complex models and observed data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

HST/WFC3 Complete Phase-resolved Spectroscopy of White Dwarf-Brown Dwarf Binaries WD 0137 and EPIC 2122

Brown dwarfs in close-in orbits around white dwarfs offer an excellent opportunity to investigate properties of fast-rotating, tidally-locked, and highly-irradiated atmospheres. We present Hubble Space Telescope Wide Field Camera 3 G141 phase-resolved observations of two brown dwarf-white dwarf binaries: WD 0137-349 and EPIC 212235321. Their 1.1 to 1.7 $μ$m phase curves demonstrate rotational modulations with semi-amplitudes of $5.27\pm0.02$% and $29.1\pm0.1$%; both can be well fit by multi-order Fourier series models. The high-order Fourier components have the same phase as the first order and are likely caused by hot spots located at the substellar points, suggesting inefficient day/night heat transfer. Both brown dwarfs' phase-resolved spectra can be accurately represented by linear combinations of their day- and night-side spectra. Fitting the irradiated brown dwarf model grids to the day-side spectra require a filling factor of ~50%, further supporting a hot spot dominating the emission of the day-sides. The night-side spectrum of WD 0137-349B is reasonably well fit by non-irradiated substellar models and the one of EPIC 212235321B can be approximated by a Planck function. We find strong spectral variations in the brown dwarfs' day/night flux and brightness temperature contrasts, which highlights the limitations of band-integrated measurements in probing heat transfer in irradiated objects. On the color-magnitude diagram, WD 0137-349B evolves along a cloudless model track connecting the early-L and mid-T spectral types, demonstrating that clouds and disequilibrium chemistry have a negligible effect on this object. A full interpretation of these high-quality phase-resolved spectra calls for new models that couple atmospheric circulation and radiative transfer under high-irradiation conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Chemistry of Temperate Super-Earth and Mini-Neptune Atmospheric Hazes from Laboratory Experiments

Very little experimental work has been done to explore the properties of photochemical hazes formed in atmospheres with very different compositions or temperatures than that of the outer solar system or of early Earth. With extrasolar planet discoveries now numbering thousands, this untapped phase space merits exploration. This study presents measured chemical properties of haze particles produced in laboratory analogues of exoplanet atmospheres. We used very high resolution mass spectrometry to measure the chemical components of solid particles produced in atmospheric chamber experiments. Many complex molecular species with general chemical formulas C$_w$H$_x$N$_y$O$_z$ were detected. We detect molecular formulas of prebiotic interest in the data, including those for the monosaccharide glyceraldehyde, a variety of amino acids and nucleotide bases, and several sugar derivatives. Additionally, the experimental exoplanetary haze analogues exhibit diverse solubility characteristics, which provide insight into the possibility of further chemical or physical alteration of photochemical hazes in super-Earth and mini-Neptune atmospheres. These exoplanet analogue particles can help us better understand chemical atmospheric processes and suggest a possible source of in situ atmospheric prebiotic chemistry on distant worlds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

COol Companions ON Ultrawide orbiTS (COCONUTS). I. A High-Gravity T4 Benchmark around an Old White Dwarf and A Re-Examination of the Surface-Gravity Dependence of the L/T Transition

We present the first discovery from the COol Companions ON Ultrawide orbiTS (COCONUTS) program, a large-scale survey for wide-orbit planetary and substellar companions. We have discovered a co-moving system COCONUTS-1, composed of a hydrogen-dominated white dwarf (PSO J058.9855+45.4184; $d=31.5$ pc) and a T4 companion (PSO J058.9869+45.4296) at a $40.6''$ (1280 au) projected separation. We derive physical properties for COCONUTS-1B from (1) its near-infrared spectrum using cloudless Sonora atmospheric models, and (2) its luminosity and the white dwarf's age ($7.3_{-1.6}^{+2.8}$ Gyr) using Sonora evolutionary models. The two methods give consistent temperatures and radii, but atmospheric models infer a lower surface gravity and therefore an unphysically young age. Assuming evolutionary model parameters ($T_{\rm eff}=1255^{+6}_{-8}$ K, $\log{g}=5.44^{+0.02}_{-0.03}$ dex, $R=0.789^{+0.011}_{-0.005}$ R$_{\rm Jup}$), we find cloudless model atmospheres have brighter Y- and J-band fluxes than the data, suggesting condensate clouds have not fully dispersed around 1300 K. The W2 flux (4.6 $μ$m) of COCONUTS-1B is fainter than models, suggesting non-equilibrium mixing of CO. To investigate the gravity dependence of the L/T transition, we compile all 60 known L6-T6 benchmarks and derive a homogeneous set of temperatures, surface gravities, and masses. As is well-known, young, low-gravity late-L dwarfs have significantly fainter, redder near-infrared photometry and $\approx200-300$ K cooler temperatures than old, high-gravity objects. Our sample now reveals such gravity dependence becomes weaker for T dwarfs, with young objects having comparable near-infrared photometry and $\approx100$ K cooler temperatures compared to old objects. Finally, we find that young objects have a larger amplitude J-band brightening than old objects, and also brighten at H band as they cross the L/T transition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Debris Disk Results from the Gemini Planet Imager Exoplanet Survey's Polarimetric Imaging Campaign

We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $μ$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 °$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Enabling Effective Exoplanet / Planetary Collaborative Science

The field of exoplanetary science has emerged over the past two decades, rising up alongside traditional solar system planetary science. Both fields focus on understanding the processes which form and sculpt planets through time, yet there has been less scientific exchange between the two communities than is ideal. This white paper explores some of the institutional and cultural barriers which impede cross-discipline collaborations and suggests solutions that would foster greater collaboration. Some solutions require structural or policy changes within NASA itself, while others are directed towards other institutions, including academic publishers, that can also facilitate greater interdisciplinarity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Haze Formation in Warm H2-rich Exoplanet Atmospheres

New observing capabilities coming online over the next few years will provide opportunities for characterization of exoplanet atmospheres. However, clouds/hazes could be present in the atmospheres of many exoplanets, muting the amplitude of spectral features. We use laboratory simulations to explore photochemical haze formation in H2-rich exoplanet atmospheres at 800 K with metallicity either 100 and 1000 times solar. We find that haze particles are produced in both simulated atmospheres with small particle size (20 to 140 nm) and relative low production rate (2.4 x 10-5 to 9.7 x 10-5 mg cm-3 h-1), but the particle size and production rate is dependent on the initial gas mixtures and the energy sources used in the simulation experiments. The gas phase mass spectra show that complex chemical processes happen in these atmospheres and generate new gas products that can further react to form larger molecules and solid haze particles. Two H2-rich atmospheres with similar C/O ratios (~0.5) yield different haze particles size, haze production rate, and gas products, suggesting both the elemental abundances and their bonding environments in an atmosphere can significantly affect the photochemistry. There is no methane (CH4) in our initial gas mixtures, although CH4 is often believed to be required to generate organic hazes. However, haze production rates from our experiments with different initial gas mixtures indicate that CH4 is neither required to generate organic hazes nor necessary to promote the organic haze formation. The variety and relative yield of the gas products indicate that CO and N2 enrich chemical reactions in H2-rich atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

HD 165054: an astrometric calibration field for high-contrast imagers in Baade's Window

We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square arcsecond with $H<22$), necessary because of the small field-of-view of the Gemini Planet Imager. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (angular separation $< 2$ arcsec). We achieved a proper motion precision of $\sim 0.3$ mas/yr, and constrained the parallax of each star to be $\lesssim 1$ mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with an excellent astrometric calibration. We find that predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Helios-r2 -- A new Bayesian, open-source retrieval model for brown dwarfs and exoplanet atmospheres

We present an improved, hybrid CPU-GPU atmospheric retrieval code, Helios-r2, which is applicable to medium-resolution emission spectra of brown dwarfs, in preparation for precision atmospheric spectroscopy in the era of the James Webb Space Telescope. The model is available as open-source code on the Exoclimes Simulation Platform. We subject Helios-r2 to a battery of tests of varying difficulty. The simplest test involves a mock retrieval on a forward model generated using the same radiative transfer technique, the same implementation of opacities, and the same chemistry model. The least trivial test involves a mock retrieval on synthetic spectra from the Sonora model grid, which uses a different radiative transfer technique, a different implementation of opacities, and a different chemistry model. A calibration factor, which is included to capture uncertainties in the brown dwarf radius, distance to the brown dwarf and flux calibration of the spectrum, may compensate, sometimes erroneously, for discrepancies in modeling choices and implementation. We analyze spectra of the benchmark brown dwarf GJ 570 D and the binary brown dwarf companions in the Epsilon Indi system. The retrieved surface gravities are consistent with previous studies and/or values inferred from dynamical masses (for Epsilon Indi Ba and Bb only). There remains no clear criterion on how to reject unphysical values of the retrieved brown dwarf radii. The inferred radii and corresponding masses should be taken with great caution. The retrieved carbon-to-oxygen ratios and metallicity depend on whether chemical equilibrium is assumed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multiband Polarimetric Imaging of HR 4796A with the Gemini Planet Imager

HR4796A hosts a well-studied debris disk with a long history due to its high fractional luminosity and favorable inclination lending itself well to both unresolved and resolved observations. We present new J- and K1-band images of the resolved debris disk HR4796A taken in the polarimetric mode of the Gemini Planet Imager (GPI). The polarized intensity features a strongly forward scattered brightness distribution and is undetected at the far side of the disk. The total intensity is detected at all scattering angles and also exhibits a strong forward scattering peak. We use a forward modelled geometric disk in order to extract geometric parameters, polarized fraction and total intensity scattering phase functions for these data as well as H-band data previously taken by GPI. We find the polarized phase function becomes increasingly more forward scattering as wavelength increases. We fit Mie and distribution of hollow spheres grain (DHS) models to the extracted functions. We find that while it is possible to describe generate a satisfactory model for the total intensity using a DHS model, but not with a Mie model. We find that no single grain population of DHS or Mie grains of arbitrary composition can simultaneously reproduce the polarized fraction and total intensity scattering phase functions, indicating the need for more sophisticated grain models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Observations of Disequilibrium CO Chemistry in the Coldest Brown Dwarfs

Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects discovered by future direct imaging surveys. We present a low resolution R $\sim$ 300 $M$-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 K and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K$_{zz}$) to estimate the strength of vertical mixing in each object. The K$_{zz}$ values of cooler gaseous objects are close to their theoretical maximum and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K$_{zz}$ values imply that disequilibrium phosphine (PH$_{3}$) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH$_{3}$ absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855's $M$-band spectrum by including both CO and water clouds in the atmospheric model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Revised Astrometric Calibration of the Gemini Planet Imager

We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in reduced GPI images. As well as introducing a small error in position angle measurements for targets observed at small zenith distances, this error led to a significant error in the previous astrometric calibration that has affected all subsequent astrometric measurements. We present a detailed description of these issues, and how they were corrected. We reduced GPI observations of calibration binaries taken periodically since the instrument was commissioned in 2014 using an updated version of the DRP. These measurements were compared to observations obtained with the NIRC2 instrument on Keck II, an instrument with an excellent astrometric calibration, allowing us to derive an updated plate scale and north offset angle for GPI. This revised astrometric calibration should be used to calibrate all measurements obtained with GPI for the purposes of precision astrometry.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sulfur-driven Haze Formation in Warm CO2-rich Exoplanet Atmospheres

Sulfur gases significantly affect the photochemistry of planetary atmospheres in our Solar System, and are expected to be important components in exoplanet atmospheres. However, sulfur photochemistry in the context of exoplanets is poorly understood due to a lack of chemical-kinetics information for sulfur species under relevant conditions. Here, we study the photochemical role of hydrogen sulfide (H2S) in warm CO2-rich exoplanet atmospheres (800 K) by carrying out laboratory simulations. We find that H2S plays a significant role in photochemistry, even when present in the atmosphere at relatively low concentrations (1.6%). It participates in both gas and solid phase chemistry, leading to the formation of other sulfur gas products (CH3SH/SO, C2H4S/OCS, SO2/S2, and CS2) and to an increase in solid haze particle production and compositional complexity. Our study shows that we may expect thicker haze with small particle sizes (20 to 140 nm) for warm CO2-rich exoplanet atmospheres that possess H2S.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Gemini Planet Imager view of the HD 32297 debris disk

We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Toward Complete Characterization: Prospects for Directly Imaging Transiting Exoplanets

High contrast direct imaging of exoplanets can provide many important observables, including measurements of the orbit, spectra that probe the lower layers of the atmosphere, and phase variations of the planet, but cannot directly measure planet radius or mass. Our future understanding of directly imaged exoplanets will therefore rely on extrapolated models of planetary atmospheres and bulk composition, which need robust calibration. We estimate the population of extrasolar planets that could serve as calibrators for these models. Critically, this population of "standard planets" must be accessible to both direct imaging and the transit method, allowing for radius measurement. We show that the search volume of a direct imaging mission eventually overcomes the transit probability falloff with semi-major axis, so that as long as cold planets are not exceedingly rare, the population of transiting planets and directly imageable planets overlaps. Using current extrapolations of Kepler occurrence rates, we estimate that ~8 standard planets could be characterized shortward of 800 nm with an ambitious future direct imaging mission like LUVOIR-A and several dozen could be detected at V band. We show the design space that would expand the sample size and discuss the extent to which ground- and space-based surveys could detect this small but crucial population of planets.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An updated visual orbit of the directly-imaged exoplanet 51 Eridani b and prospects for a dynamical mass measurement with Gaia

We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$\,au and a period of $28.1_{-4.9}^{+17.2}$\,yr, assuming a mass of 1.75\,M$_{\odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {\it Hipparcos} and {\it Gaia} catalogues was found to be inconsistent at the 2--3$σ$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {\it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {\it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$\% probability) given the most optimistic predictions for the {\it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($\gtrsim3.6$\,M$_{\rm Jup}$).

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Detection of a low-mass stellar companion to the accelerating A2IV star HR 1645

The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to be resolved with ground-based high contrast imaging instruments; long-period equal mass ratio stellar companions that are also consistent with the measured acceleration are excluded with previous imaging observations. The small but significant amplitude of the acceleration made HR 1645 a promising candidate for targeted searches for brown dwarf and planetary-mass companions around nearby, young stars. In this paper we explore the origin of the astrometric acceleration by modelling the signal induced by a wide-orbit M8 companion discovered with the Gemini Planet Imager, as well as the effects of an inner short-period spectroscopic companion discovered a century ago but not since followed-up. We present the first constraints on the orbit of the inner companion, and demonstrate that it is a plausible cause of the astrometric acceleration. This result demonstrates the importance of vetting of targets with measured astrometric acceleration for short-period stellar companions prior to conducting targeted direct imaging surveys for wide-orbit substellar companions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

First Resolved Scattered-Light Images of Four Debris Disks in Scorpius-Centaurus with the Gemini Planet Imager

We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet-disk system. Uniquely, HD 98363 has a wide co-moving stellar companion Wray 15-788 with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one non-detection to one of the brightest disks in the region. With the new results, a total of 15 A/F-stars in Sco-Cen have resolved scattered light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F-stars at this important stage in the development of planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet beta Pictoris b from Combined Direct Imaging and Astrometry

We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: 12.8 [+5.3, -3.2] M_Jup. Our eccentricity determination of 0.12 [+0.04, -0.03] disfavors circular orbits. We consider orbit fits to several different imaging datasets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well constrained within 2.8 days of 2017 September 13, with the star behind the planet's Hill sphere between 2017 April 11 and 2018 February 16 (+/- 18 days). Following the recent radial velocity detection of a second planet in the system, beta Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass for the imaged planet beta Pic b, of 8.0 +/- 2.6 M_Jup, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Possible Astrometric Signature of a Planetary-mass Companion to the Nearby Young Star TW Piscis Austrini (Fomalhaut B): Constraints from Astrometry, Radial Velocities, and Direct Imaging

We present constraints on the presence of substellar companions to the nearby ($d\sim7.6$ pc), young ($440\pm40$ Myr) K4Ve star TW Piscis Austrini, the wide ($\sim$0.3 pc) companion to the A4V star Fomalhaut. We combined absolute astrometry from Hipparcos and Gaia with literature radial velocity measurements and dedicated high-contrast imaging observations obtained with Keck/NIRC2 to achieve sensitivity to brown dwarf and planetary-mass companions ($\gtrsim2$ M$_{\rm Jup}$) over many decades of orbital period ($\lesssim 10^3$ yr). The significant astrometric acceleration measured between the Hipparcos and Gaia catalogues, reported previously in the literature, cannot be explained by the orbital motion of TW PsA around the barycenter of the Fomalhaut triple system. Instead, we find that it is consistent with the reflex motion induced by an orbiting substellar companion. The combination of astrometry, radial velocities, and a deep $L^{\prime}$ imaging dataset leads to a constraint on the companion mass of $1.2_{-0.6}^{+0.7}$ M$_{\rm Jup}$. However, the period of the companion is poorly constrained, with a highly multi-modal period posterior distribution due to aliasing with the 24.25 year baseline between Hipparcos and Gaia. If confirmed through continued astrometric or spectroscopic monitoring, or via direct detection, the companion to TW PsA would represent a choice target for detailed atmospheric characterization with high-contrast instruments on the upcoming James Webb Space Telescope and Wide Field Infrared Survey Telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Transmission Spectroscopy of WASP-79b from 0.6 to 5.0 $μ$m

As part of the PanCET program, we have conducted a spectroscopic study of WASP-79b, an inflated hot Jupiter orbiting an F-type star in Eridanus with a period of 3.66 days. Building on the original WASP and TRAPPIST photometry of Smalley et al (2012), we examine HST/WFC3 (1.125 - 1.650 $μ$m), Magellan/LDSS-3C (0.6 - 1 $μ$m) data, and Spitzer data (3.6 and 4.5 $μ$m). Using data from all three instruments, we constrain the water abundance to be --2.20 $\leq$ log(H$_2$O) $\leq$ --1.55. We present these results along with the results of an atmospheric retrieval analysis, which favor inclusion of FeH and H$^-$ in the atmospheric model. We also provide an updated ephemeris based on the Smalley, HST/WFC3, LDSS-3C, Spitzer, and TESS transit times. With the detectable water feature and its occupation of the clear/cloudy transition region of the temperature/gravity phase space, WASP-79b is a target of interest for the approved JWST Director's Discretionary Early Release Science (DD ERS) program, with ERS observations planned to be the first to execute in Cycle 1. Transiting exoplanets have been approved for 78.1 hours of data collection, and with the delay in the JWST launch, WASP-79b is now a target for the Panchromatic Transmission program. This program will observe WASP-79b for 42 hours in 4 different instrument modes, providing substantially more data by which to investigate this hot Jupiter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

$Extrasolar~Storms$: Pressure-dependent Changes In Light Curve Phase In Brown Dwarfs From Simultaneous $Hubble$ and $Spitzer$ Observations

We present $Spitzer$/IRAC Ch1 and Ch2 monitoring of six brown dwarfs during 8 different epochs over the course of 20 months. For four brown dwarfs, we also obtained simulataneous $HST$/WFC3 G141 Grism spectra during two epochs and derived light curves in five narrow-band filters. Probing different pressure levels in the atmospheres, the multi-wavelength 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 h to 13 h. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cloud Atlas: Discovery of Patchy Clouds and High-amplitude Rotational Modulations In a Young, Extremely Red L-type Brown Dwarf

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 μ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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detecting Exomoons Around Self-luminous Giant Exoplanets Through Polarization

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 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 which are appropriate for self-luminous exoplanets, we present the polarization profiles of these objects in the infrared during 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Discovery of a Substellar Companion to the Nearby Debris Disk Host HR 2562

We present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observations to more than $5σ$. Spectral data in $J$, $H$, and $K$ bands show morphological similarity to L/T transition objects. We assign a spectral type of L7$\pm$3 to HR 2562B, and derive a luminosity of $\log$(L$_{\rm bol}$/L$_{\odot}$)=-4.62$\pm$0.12, corresponding to a mass of 30$\pm$15 M$_{\rm Jup}$ from evolutionary models at an estimated age of the system of 300-900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Exploring the Role of Sub-micron Sized Dust Grains in the Atmospheres of Red L0 - L6 Dwarfs

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 which 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 submicron 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High temperature condensate clouds in super-hot Jupiter atmospheres

Deciphering the role of clouds is central to our understanding of exoplanet atmospheres, as they have a direct impact on the temperature and pressure structure, and observational properties of the planet. Super-hot Jupiters occupy a temperature regime similar to low mass M-dwarfs, where minimal cloud condensation is expected. However, observations of exoplanets such as WASP-12b (Teq ~ 2500 K) result in a transmission spectrum indicative of a cloudy atmosphere. We re-examine the temperature and pressure space occupied by these super-hot Jupiter atmospheres, to explore the role of the initial Al- and Ti-bearing condensates as the main source of cloud material. Due to the high temperatures a majority of the more common refractory material is not depleted into deeper layers and would remain in the vapor phase. The lack of depletion into deeper layers means that these materials with relatively low cloud masses can become significant absorbers in the upper atmosphere. We provide condensation curves for the initial Al- and Ti-bearing condensates that may be used to provide quantitative estimates of the effect of metallicity on cloud masses, as planets with metal-rich hosts potentially form more opaque clouds because more mass is available for condensation. Increased metallicity also pushes the point of condensation to hotter, deeper layers in the planetary atmosphere further increasing the density of the cloud. We suggest that planets around metal-rich hosts are more likely to have thick refractory clouds, and discuss the implication on the observed spectra of WASP-12b.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Imaging an 80 AU Radius Dust Ring Around the F5V Star HD 157587

We present $H$-band near-infrared polarimetric imaging observations of the F5V star HD~157587 obtained with the Gemini Planet Imager (GPI) that reveal the debris disk as a bright ring structure at a separation of $\sim$80$-$100~AU. The new GPI data complement recent HST/STIS observations that show the disk extending out to over 500~AU. The GPI image displays a strong asymmetry along the projected minor axis as well as a fainter asymmetry along the projected major axis. We associate the minor and major axis asymmetries with polarized forward scattering and a possible stellocentric offset, respectively. To constrain the disk geometry we fit two separate disk models to the polarized image, each using a different scattering phase function. Both models favor a disk inclination of $\sim 70\degr$ and a $1.5\pm0.6$ AU stellar offset in the plane of the sky along the projected major axis of the disk. We find that the stellar offset in the disk plane, perpendicular to the projected major axis is degenerate with the form of the scattering phase function and remains poorly constrained. The disk is not recovered in total intensity due in part to strong adaptive optics residuals, but we recover three point sources. Considering the system's proximity to the galactic plane and the point sources' positions relative to the disk, we consider it likely that they are background objects and unrelated to the disk's offset from the star.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Maps of Evolving Cloud Structures in Luhman 16AB from HST Time-Resolved Spectroscopy

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 & 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 structure of brown dwarf and other ultracool atmospheres, to light curves of Luhman 16A & B using the Hubble Space Telescope's G141 and G102 grisms. Aeolus retrieved three or four spots in the top-of-the-atmosphere of Luhman 16A & 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 top-of-the-atmosphere. We compared our Luhman 16B map with the only previously published map. Interestingly, our map contained a large, cooler (DT~51 K) than the background top-of-the-atmosphere spot that 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 this feature is related to top-of-the-atmosphere structures of Luhman 16B.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Photolytic Hazes in the Atmosphere of 51 Eri b

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Point Source Polarimetry with the Gemini Planet Imager: Sensitivity Characterization with T5.5 Dwarf Companion HD 19467 B

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 $p_{\text{CL}99.73\%} \leq 2.4\%$. We discuss our results in the context of T dwarf cloud models and photometric variability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The atmospheric circulation of a nine-hot Jupiter sample: Probing circulation and chemistry over a wide phase space

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The First Detection of Photometric Variability in a Y Dwarf: WISE J140518.39+553421.3

We present the first detection of photometric variability of a spectroscopically-confirmed Y dwarf. The Infrared Array Camera on board the Spitzer Space Telescope was used to obtain times series photometry at 3.6 and 4.5 microns over a twenty four hour period at two different epochs separated by 149 days. Variability is evident at 4.5 um in the first epoch and at 3.6 and 4.5 um 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 hours, 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. More advanced two- or three-dimensional models that include time-dependent phenomena like vertical mixing, cloud formation, and thermal relaxation are therefore sorely needed in order to properly interpret our observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Interior of Saturn

We review our current understanding of the interior structure and thermal evolution of Saturn, with a focus on recent results in the Cassini era. There has been important progress in understanding physical inputs, including equations of state of planetary materials and their mixtures, physical parameters like the gravity field and rotation rate, and constraints on Saturnian free oscillations. At the same time, new methods of calculation, including work on the gravity field of rotating fluid bodies, and the role of interior composition gradients, should help to better constrain the state of Saturn's interior, now and earlier in its history. However, a better appreciation of modeling uncertainties and degeneracies, along with a greater exploration of modeling phase space, still leave great uncertainties in our understanding of Saturn's interior. Further analysis of Cassini data sets, as well as precise gravity field measurements from the Cassini Grand Finale orbits, will further revolutionize our understanding of Saturn's interior over the next few years.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Need for Laboratory Work to Aid in The Understanding of Exoplanetary Atmospheres

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Orbit and Transit Prospects for $β$ Pictoris b constrained with One Milliarcsecond Astrometry

A principal scientific goal of the Gemini Planet Imager (GPI) is obtaining milliarcsecond astrometry to constrain exoplanet orbits. However, astrometry of directly imaged exoplanets is subject to biases, systematic errors, and speckle noise. Here we describe an analytical procedure to forward model the signal of an exoplanet that accounts for both the observing strategy (angular and spectral differential imaging) and the data reduction method (Karhunen-Loève Image Projection algorithm). We use this forward model to measure the position of an exoplanet in a Bayesian framework employing Gaussian processes and Markov chain Monte Carlo (MCMC) to account for correlated noise. In the case of GPI data on $β$ Pic b, this technique, which we call Bayesian KLIP-FM Astrometry (BKA), outperforms previous techniques and yields 1$σ$-errors at or below the one milliarcsecond level. We validate BKA by fitting a Keplerian orbit to twelve GPI observations along with previous astrometry from other instruments. The statistical properties of the residuals confirm that BKA is accurate and correctly estimates astrometric errors. Our constraints on the orbit of $β$ Pic b firmly rule out the possibility of a transit of the planet at 10-$σ$ significance. However, we confirm that the Hill sphere of $β$ Pic b will transit, giving us a rare chance to probe the circumplanetary environment of a young, evolving exoplanet. We provide an ephemeris for photometric monitoring of the Hill sphere transit event, which will begin at the start of April in 2017 and finish at the end of January in 2018.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Transitions in the cloud composition of hot Jupiters

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 lightcurves of some hot Jupiters are asymmetric: for the hottest planets, the lightcurve peaks before secondary eclipse, whereas for planets cooler than $\sim1900\rm\,K$, it peaks after secondary eclipse. We use the thermal structure from 3D global circulation models to determine the expected cloud distribution and Kepler lightcurves of hot Jupiters. We demonstrate that the change from an optical lightcurve 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 lightcurve 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\rm\,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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

$β$ Pictoris' inner disk in polarized light and new orbital parameters for $β$ Pictoris b

We present $H$-band observations of $β$ Pic with the Gemini Planet Imager's (GPI's) polarimetry mode that reveal the debris disk between ~0.3" (~6 AU) and ~1.7" (~33 AU), while simultaneously detecting $β$ Pic $b$. The polarized disk image was fit with a dust density model combined with a Henyey-Greenstein scattering phase function. The best fit model indicates a disk inclined to the line of sight ($ϕ=85.27°^{+0.26}_{-0.19}$) with a position angle $θ_{PA}=30.35°^{+0.29}_{-0.28}$ (slightly offset from the main outer disk, $θ_{PA}\approx29°$), that extends from an inner disk radius of $23.6^{+0.9}_{-0.6}$ AU to well outside GPI's field of view. In addition, we present an updated orbit for $β$ Pic $b$ based on new astrometric measurements taken in GPI's spectroscopic mode spanning 14 months. The planet has a semi-major axis of $a=9.2^{+1.5}_{-0.4}$AU, with an eccentricity $e\leq 0.26$. The position angle of the ascending node is $Ω=31.75°\pm0.15$, offset from both the outer main disk and the inner disk seen in the GPI image. The orbital fit constrains the stellar mass of $β$ Pic to $1.60\pm0.05 M_{\odot}$. Dynamical sculpting by $β$ Pic $b$ cannot easily account for the following three aspects of the inferred disk properties: 1) the modeled inner radius of the disk is farther out than expected if caused by $β$ Pic b; 2) the mutual inclination of the inner disk and $β$ Pic $b$ is $4°$, when it is expected to be closer to zero; and 3) the aspect ratio of the disk ($h_0 = 0.137^{+0.005}_{-0.006}$) is larger than expected from interactions with $β$ Pic $b$ or self-stirring by the disk's parent bodies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Characterizing Rocky and Gaseous Exoplanets with 2-meter Class Space-based Coronagraphs

Several concepts now exist for small, space-based missions to directly characterize exoplanets in reflected light. 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 um, 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, 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 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 um methane band. Additionally, most modeled observations have noise dominated by dark current, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Discovery of Rotational Modulations in the Planetary-Mass Companion 2M1207b: Intermediate Rotation Period and Heterogeneous Clouds in a Low Gravity Atmosphere

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 HST/WFC3 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 $10.7^{+1.2}_{-0.6}$ hours 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neptune's Dynamic Atmosphere from Kepler K2 Observations: Implications for Brown Dwarf Light Curve Analyses

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-meter telescope at 1.65 microns and Hubble Space Telescope visible imaging acquired 9 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 time scales, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A non-grey analytical model for irradiated atmospheres. II: Analytical vs. numerical solutions

The recent discovery and characterization of the diversity of the atmospheres of exoplanets and brown dwarfs calls for the development of fast and accurate analytical models. We quantify the accuracy of the analytical solution derived in paper I for an irradiated, non-grey atmosphere by comparing it to a state-of-the-art radiative transfer model. Then, using a grid of numerical models, we calibrate the different coefficients of our analytical model for irradiated solar-composition atmospheres of giant exoplanets and brown dwarfs. We show that the so-called Eddington approximation used to solve the angular dependency of the radiation field leads to relative errors of up to 5% on the temperature profile. We show that for realistic non-grey planetary atmospheres, the presence of a convective zone that extends to optical depths smaller than unity can lead to changes in the radiative temperature profile on the order of 20% or more. When the convective zone is located at deeper levels (such as for strongly irradiated hot Jupiters), its effect on the radiative atmosphere is smaller. We show that the temperature inversion induced by a strong absorber in the optical, such as TiO or VO is mainly due to non-grey thermal effects reducing the ability of the upper atmosphere to cool down rather than an enhanced absorption of the stellar light as previously thought. Finally, we provide a functional form for the coefficients of our analytical model for solar-composition giant exoplanets and brown dwarfs. This leads to fully analytical pressure-temperature profiles for irradiated atmospheres with a relative accuracy better than 10% for gravities between 2.5m/s^2 and 250 m/s^2 and effective temperatures between 100 K and 3000 K. This is a great improvement over the commonly used Eddington boundary condition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Approach for Retrieving Temperatures and Abundances in Brown Dwarf Atmospheres

Brown dwarf spectra contain a wealth of information about their molecular abundances, temperature structure, and gravity. We present a new data driven retrieval approach, previously used in planetary atmosphere studies, to extract the molecular abundances and temperature structure from brown dwarf spectra. The feasibility of the approach is first demonstrated on a synthetic brown dwarf spectrum. Given typical spectral resolutions, wavelength coverage, and noise properties precisions of tens of percent can be obtained for the molecular abundances and 10s-100s K on the temperature profile. The technique is then applied to the well studied brown dwarf, Gl 570D. From this spectral retrieval the spectroscopic radius is constrained to be 0.75 - 0.83 $R_{\mathrm{J}}$, $\log(g)$ to be 5.13 - 5.46 and $T_{\mathrm{eff}}$ to be between 804 and 849 K. Estimates for the range of abundances and allowed temperature profiles are also derived. The results from our retrieval approach are in agreement with the self-consistent grid modeling results of Saumon et al (2006). This new approach will allow us to address issues of compositional differences between brown dwarfs and possibly their formation environments, disequilibrium chemistry, missing physics in current grid modeling approaches as well as a many other issues.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Atmospheric Circulation of Eccentric Hot Jupiter HAT-P-2b

The hot-Jupiter HAT-P-2b has become a prime target for Spitzer Space Telescope observations aimed at understanding the atmospheric response of exoplanets on highly eccentric orbits. Here we present a suite of three-dimensional atmospheric circulation models for HAT-P-2b that investigate the effects of assumed atmospheric composition and rotation rate on global scale winds and thermal patterns. We compare and contrast atmospheric models for HAT-P-2b, which assume one and five times solar metallicity, both with and without TiO/VO as atmospheric constituents. Additionally we compare models that assume a rotation period of half, one, and two times the nominal pseudo-synchronous rotation period. We find that changes in assumed atmospheric metallicity and rotation rate do not significantly affect model predictions of the planetary flux as a function of orbital phase. However, models in which TiO/VO are present in the atmosphere develop a transient temperature inversion between the transit and secondary eclipse events that results in significant variations in the timing and magnitude of the peak of the planetary flux compared with models in which TiO/VO are omitted from the opacity tables. We find that no one single atmospheric model can reproduce the recently observed full orbit phase curves at 3.6, 4.5 and 8.0 microns, which is likely due to a chemical process not captured by our current atmospheric models for HAT-P-2b. Further modeling and observational efforts focused on understanding the chemistry of HAT-P-2b's atmosphere are needed and could provide key insights into the interplay between radiative, dynamical, and chemical processes in a wide range of exoplanet atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cloud structure of brown dwarfs from spectroscopic variability observations

Recent discoveries of variable brown dwarfs have provided us with a new window into their three-dimensional cloud structure. The highest variables are found at the L/T transition, where the cloud cover is thought to break up, but variability has been found to occur also for both cloudy L dwarfs and (mostly) cloud-free mid T dwarfs. We summarize results from recent HST programs measuring the spectral variability of brown dwarfs in the near-infrared and compare to results from ground-based programs. We discuss the patchy cloud structure of L/T transition objects, for which it is becoming increasingly certain that the variability does not arise from cloud holes into the deep hot regions but from varying cloud thickness. We present a new patchy cloud model to explain the spectral variability of 2MASSJ21392676+0220226. We also discuss the curious multi-wavelength variability behavior of the recently discovered very nearby early T dwarf WISE J104915.57-531906.1B (Luhman 16B) and the mid T dwarf 2MASS J22282889-431026.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Directly Imaged L-T Transition Exoplanets in the Mid-Infrared

Gas-giant planets emit a large fraction of their light in the mid-infrared ($\gtrsim$3$μ$m), where photometry and spectroscopy are critical to our understanding of the bulk properties of extrasolar planets. Of particular importance are the L and M-band atmospheric windows (3-5$μ$m), which are the longest wavelengths currently accessible to ground-based, high-contrast imagers. We present binocular LBT AO images of the HR 8799 planetary system in six narrow-band filters from 3-4$μ$m, and a Magellan AO image of the 2M1207 planetary system in a broader 3.3$μ$m band. These systems encompass the five known exoplanets with luminosities consistent with L$\rightarrow$T transition brown dwarfs. Our results show that the exoplanets are brighter and have shallower spectral slopes than equivalent temperature brown dwarfs in a wavelength range that contains the methane fundamental absorption feature (spanned by the narrowband filters and encompassed by the broader 3.3$μ$m filter). For 2M1207 b, we find that thick clouds and non-equilibrium chemistry caused by vertical mixing can explain the object's appearance. For the HR 8799 planets, we present new models that suggest the atmospheres must have patchy clouds, along with non-equilibrium chemistry. Together, the presence of a heterogeneous surface and vertical mixing presents a picture of dynamic planetary atmospheres in which both horizontal and vertical motions influence the chemical and condensate profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gemini Planet Imager Spectroscopy of the HR 8799 planets c and d

During the first-light run of the Gemini Planet Imager (GPI) we obtained K-band spectra of exoplanets HR 8799 c and d. Analysis of the spectra indicates that planet d may be warmer than planet c. Comparisons to recent patchy cloud models and previously obtained observations over multiple wavelengths confirm that thick clouds combined with horizontal variation in the cloud cover generally reproduce the planets' spectral energy distributions. When combined with the 3 to 4 um photometric data points, the observations provide strong constraints on the atmospheric methane content for both planets. The data also provide further evidence that future modeling efforts must include cloud opacity, possibly including cloud holes, disequilibrium chemistry, and super-solar metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

HST Rotational Spectral Mapping of Two L-Type Brown Dwarfs: Variability In and Out of Water Bands Indicates High-Altitude Haze Layers

We present time-resolved near-infrared spectroscopy of two L5 dwarfs, 2MASS J18212815+1414010 and 2MASS J15074759-1627386, observed with the Wide Field Camera 3 instrument on the Hubble Space Telescope (HST). We study the wavelength dependence of rotation-modulated flux variations between 1.1 $μ$m and 1.7 $μ$m. We find that the water absorption bands of the two L5 dwarfs at 1.15 $μ$m and 1.4 $μ$m vary at similar amplitudes as the adjacent continuum. This differs from the results of previous HST observations of L/T transition dwarfs, in which the water absorption at 1.4 $μ$m displays variations of about half of the amplitude at other wavelengths. We find that the relative amplitude of flux variability out of the water band with respect to that in the water band shows a increasing trend from the L5 dwarfs toward the early T dwarfs. We utilize the models of Saumon & Marley (2008) and find that the observed variability of the L5 dwarfs can be explained by the presence of spatially varying high-altitude haze layers above the condensate clouds. Therefore, our observations show that the heterogeneity of haze layers - the driver of the variability - must be located at very low pressures, where even the water opacity is negligible. In the near future, the rotational spectral mapping technique could be utilized for other atomic and molecular species to probe different pressure levels in the atmospheres of brown dwarfs and exoplanets and uncover both horizontal and vertical cloud structures.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Methane, Carbon Monoxide, and Ammonia in Brown Dwarfs and Self-Luminous Giant Planets

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).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Cool Side: Modeling the Atmospheres of Brown Dwarfs and Giant Planets

The atmosphere of a brown dwarf or extrasolar giant planet controls the spectrum of radiation emitted by the object and regulates its cooling over time. While the study of these atmospheres has been informed by decades of experience modeling stellar and planetary atmospheres, the distinctive characteristics of these objects present unique challenges to forward modeling. In particular, complex chemistry arising from molecule-rich atmospheres, molecular opacity line lists (sometimes running to 10 billion absorption lines or more) multiple cloud-forming condensates, and disequilibrium chemical processes all combine to create a challenging task for any modeling effort. This review describes the process of incorporating these complexities into one-dimensional radiative-convective equilibrium models of sub-stellar objects. We discuss the underlying mathematics as well as the techniques used to model the physics, chemistry, radiative transfer, and other processes relevant to understanding these atmospheres. The review focuses on the process of the creation of atmosphere models and briefly presents some comparisons of model predictions to data. Current challenges in the field and some comments on the future conclude the review.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Saturn Ring Seismology: Looking Beyond First Order Resonances

Some wave features found in the C-ring of Saturn appear to be excited by resonances with normal mode oscillations of the planet. The waves are found at locations in the rings where the ratio of orbital to oscillation frequencies is given by m : m+1 where m is a small integer. I suggest here that it is plausible that ring waves may also be launched at second order resonances, where the frequency ratio would be m : m+2. Indeed, otherwise unassociated wave features are found at such locations in the C-ring. If confirmed the association of planetary modes with additional C-ring wave features would measure additional oscillation frequencies of Saturn and improve the utility of the waves for constraining the internal structure of the planet. Second-order resonances in general do not lie near first order ring resonance locations and thus are not the explanation for the apparent frequency splitting of modes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spectral Variability from the Patchy Atmospheres of T and Y Dwarfs

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Temperature Fluctuations as a Source of Brown Dwarf Variability

A number of brown dwarfs are now known to be variable with observed amplitudes as large as 10-30% at some wavelengths. While spatial inhomogeneities in cloud coverage and thickness are likely responsible for much of the observed variability, it is possible that some of the variations arise from atmospheric temperature fluctuations instead of, or in addition to, clouds. To better understand the role that thermal variability might play we present a case study of brown dwarf variability using a newly-developed one-dimensional, time-stepping model of atmospheric thermal structure. We focus on the effects of thermal perturbations, intentionally simplifying the problem through omission of clouds and atmospheric circulation. Model results demonstrate that thermal perturbations occurring deep in the atmosphere (at pressures greater than 10 bar) of a model T-dwarf can be communicated to the upper atmosphere through radiative heating via the windows in near-infrared water opacity. The response time depends on where in the atmosphere a thermal perturbation is introduced. We show that, for certain periodic perturbations, the emission spectrum can have complex, time- and wavelength-dependent behaviors, including phase shifts in times of maximum flux observed at different wavelengths. Since different wavelengths probe different levels in the atmosphere, these variations track a wavelength-dependent set of radiative exchanges happening between different atmospheric levels as a perturbation evolves in time. We conclude that thermal--as well as cloud--fluctuations must be considered as possible contributors to the observed brown dwarf variability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Atmospheres of Earth-like Planets after Giant Impact Events

It is now understood that the accretion of terrestrial planets naturally involves giant collisions, the moon-forming impact being a well known example. In the aftermath of such collisions the surface of the surviving planet is very hot and potentially detectable. Here we explore the atmospheric chemistry, photochemistry, and spectral signatures of post-giant-impact terrestrial planets enveloped by thick atmospheres consisting predominantly of CO2, and H2O. The atmospheric chemistry and structure are computed self-consistently for atmospheres in equilibrium with hot surfaces with composition reflecting either the bulk silicate Earth (which includes the crust, mantle, atmosphere and oceans) or Earth's continental crust. We account for all major molecular and atomic opacity sources including collision-induced absorption. We find that these atmospheres are dominated by H2O and CO2, while the formation of CH4, and NH3 is quenched due to short dynamical timescales. Other important constituents are HF, HCl, NaCl, and SO2. These are apparent in the emerging spectra, and can be indicative that an impact has occurred. The use of comprehensive opacities results in spectra that are a factor of 2 lower in surface brightness in the spectral windows than predicted by previous models. The estimated luminosities show that the hottest post-giant-impact planets will be detectable with near-infrared coronagraphs on the planned 30m-class telescopes. The 1-4um region will be most favorable for such detections, offering bright features and better contrast between the planet and a potential debris disk. We derive cooling timescales on the order of 10^5-10^6 Myrs, based on the modeled effective temperatures. This leads to the possibility of discovering tens of such planets in future surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The atmospheric circulation of the super Earth GJ 1214b: Dependence on composition and metallicity

We present three-dimensional atmospheric circulation models of GJ 1214b, a 2.7 Earth-radius, 6.5 Earth-mass super Earth detected by the MEarth survey. Here we explore the planet's circulation as a function of atmospheric metallicity and atmospheric composition, modeling atmospheres with a low mean-molecular weight (i.e., H2-dominated) and a high mean-molecular weight (i.e. water- and CO2-dominated). We find that atmospheres with a low mean-molecular weight have strong day-night temperature variations at pressures above the infrared photosphere that lead to equatorial superrotation. For these atmospheres, the enhancement of atmospheric opacities with increasing metallicity lead to shallower atmospheric heating, larger day-night temperature variations and hence stronger superrotation. In comparison, atmospheres with a high mean-molecular weight have larger day-night and equator-to-pole temperature variations than low mean-molecular weight atmospheres, but differences in opacity structure and energy budget lead to differences in jet structure. The circulation of a water-dominated atmosphere is dominated by equatorial superrotation, while the circulation of a CO2-dominated atmosphere is instead dominated by high-latitude jets. By comparing emergent flux spectra and lightcurves for 50x solar and water-dominated compositions, we show that observations in emission can break the degeneracy in determining the atmospheric composition of GJ 1214b. The variation in opacity with wavelength for the water-dominated atmosphere leads to large phase variations within water bands and small phase variations outside of water bands. The 50x solar atmosphere, however, yields small variations within water bands and large phase variations at other characteristic wavelengths. These observations would be much less sensitive to clouds, condensates, and hazes than transit observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The First H-band Spectrum of the Massive Gas Giant Planet beta Pictoris b with the Gemini Planet Imager

Using the recently installed Gemini Planet Imager (GPI), we have taken the first H-band spectrum of the planetary companion to the nearby young star beta Pictoris. GPI is designed to image and provide low-resolution spectra of Jupiter sized, self-luminous planetary companions around young nearby stars. These observations were taken covering the H-band (1.65 microns). The spectrum has a resolving power of $\sim$ 45 and demonstrates the distinctive triangular shape of a cool substellar object with low surface gravity. Using atmospheric models, we find an effective temperature of $1650 \pm 50$ K and a surface gravity of $\log(g) = 4.0 \pm 0.25$ (cgs units). These values agree well with predictions from planetary evolution models for a gas giant with mass between 10 and 12 $M_{\rm Jup}$ and age between 10 and 20 Myrs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Titan solar occultation observations reveal transit spectra of a hazy world

High altitude clouds and hazes are integral to understanding exoplanet observations, and are proposed to explain observed featureless transit spectra. However, it is difficult to make inferences from these data because of the need to disentangle effects of gas absorption from haze extinction. Here, we turn to the quintessential hazy world -- Titan -- to clarify how high altitude hazes influence transit spectra. We use solar occultation observations of Titan's atmosphere from the Visual and Infrared Mapping Spectrometer (VIMS) aboard NASA's Cassini spacecraft to generate transit spectra. Data span 0.88-5 microns at a resolution of 12-18 nm, with uncertainties typically smaller than 1%. Our approach exploits symmetry between occultations and transits, producing transit radius spectra that inherently include the effects of haze multiple scattering, refraction, and gas absorption. We use a simple model of haze extinction to explore how Titan's haze affects its transit spectrum. Our spectra show strong methane absorption features, and weaker features due to other gases. Most importantly, the data demonstrate that high altitude hazes can severely limit the atmospheric depths probed by transit spectra, bounding observations to pressures smaller than 0.1-10 mbar, depending on wavelength. Unlike the usual assumption made when modeling and interpreting transit observations of potentially hazy worlds, the slope set by haze in our spectra is not flat, and creates a variation in transit height whose magnitude is comparable to those from the strongest gaseous absorption features. These findings have important consequences for interpreting future exoplanet observations, including those from NASA's James Webb Space Telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Water Clouds in Y Dwarfs and Exoplanets

The formation of clouds affects brown dwarf and planetary atmospheres of nearly all effective temperatures. Iron and silicate condense in L dwarf atmospheres and dissipate at the L/T transition. Minor species such as sulfides and salts condense in mid-late T dwarfs. For brown dwarfs below Teff=450 K, water condenses in the upper atmosphere to form ice clouds. Currently over a dozen objects in this temperature range have been discovered, and few previous theoretical studies have addressed the effect of water clouds on brown dwarf or exoplanetary spectra. Here we present a new grid of models that include the effect of water cloud opacity. We find that they become optically thick in objects below Teff=350-375 K. Unlike refractory cloud materials, water ice particles are significantly non-gray absorbers; they predominantly scatter at optical wavelengths through J band and absorb in the infrared with prominent features, the strongest of which is at 2.8 microns. H2O, NH3, CH4, and H2 CIA are dominant opacity sources; less abundant species such as may also be detectable, including the alkalis, H2S, and PH3. PH3, which has been detected in Jupiter, is expected to have a strong signature in the mid-infrared at 4.3 microns in Y dwarfs around Teff=450 K; if disequilibrium chemistry increases the abundance of PH3, it may be detectable over a wider effective temperature range than models predict. We show results incorporating disequilibrium nitrogen and carbon chemistry and predict signatures of low gravity in planetary- mass objects. Lastly, we make predictions for the observability of Y dwarfs and planets with existing and future instruments including the James Webb Space Telescope and Gemini Planet Imager.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Weather on Other Worlds. II. Survey Results: Spots Are Ubiquitous on L and T Dwarfs

We present results from the "Weather on Other Worlds" Spitzer Exploration Science program to investigate photometric variability in L and T dwarfs, usually attributed to patchy clouds. We surveyed 44 L3-T8 dwarfs, spanning a range of $J-K_s$ colors and surface gravities. We find that 14/23 (61%; 95% confidence interval: 41%-78%) of our single L3-L9.5 dwarfs are variable with peak-to-peak amplitudes between 0.2% and 1.5%, and 5/16 (31%; 95% confidence interval: 14%-56%) of our single T0-T8 dwarfs are variable with amplitudes between 0.8% and 4.6%. After correcting for sensitivity, we find that 80% (95% confidence interval: 53%-100%) of L dwarfs vary by >0.2%, and 36% (95% confidence interval: 19%-52%) of T dwarfs vary by >0.4%. Given viewing geometry considerations, we conclude that photospheric heterogeneities causing >0.2% 3-5-micron flux variations are present on virtually all L dwarfs, and probably on most T dwarfs. A third of L dwarf variables show irregular light curves, indicating that L dwarfs may have multiple spots that evolve over a single rotation. Also, approximately a third of the periodicities are on time scales >10 h, suggesting that slowly-rotating brown dwarfs may be common. We observe an increase in the maximum amplitudes over the entire spectral type range, revealing a potential for greater temperature contrasts in T dwarfs than in L dwarfs. We find a tentative association (92% confidence) between low surface gravity and high-amplitude variability among L3-L5.5 dwarfs. Although we can not confirm whether lower gravity is also correlated with a higher incidence of variables, the result is promising for the characterization of directly imaged young extrasolar planets through variability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

WISE Y Dwarfs As Probes of the Brown Dwarf-Exoplanet Connection

We have determined astrometric positions for 15 WISE-discovered late-type brown dwarfs (6 T8-9 and 9 Y dwarfs) using the Keck II telescope, the Spitzer Space Telescope, and the Hubble Space Telescope. Combining data from 8 to 20 epochs we derive parallactic and proper motions for these objects which put the majority within 15 parsecs. For ages greater than a few Gyr, as suggested from kinematic considerations, we find masses of 10-30 M$_{Jup}$ based on standard models for the evolution of low mass objects with a range of mass estimates for individual objects depending on the model in question. Three of the coolest objects have effective temperatures $\sim$ 350 K and inferred masses of 10-15 M$_{Jup}$. Our parallactic distances confirm earlier photometric estimates (Kirkpatrick et al. 2012) and direct measurements (Marsh et al. 2013, Beichman et al. 2013, Dupuy and Krauss 2013) and suggest that the number of objects with masses below about 15 M$_{Jup}$ must be flat or declining relative to higher mass objects. The masses of the coldest Y dwarfs may be similar to those inferred for recently imaged planet-mass companions to nearby young stars. Objects in this mass range, which appear to be rare in both the interstellar and proto-planetary environments, may both have formed via gravitational fragmentation: the brown dwarfs in interstellar clouds and companion objects in a protoplanetary disk. In both cases, however, the fact that objects in this mass range are relatively infrequent, suggests that this mechanism must be inefficient in both environments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Study of the Diverse T Dwarf Population Revealed by WISE

We report the discovery of 87 new T dwarfs uncovered with the Wide-field Infrared Survey Explorer (WISE) and three brown dwarfs with extremely red near-infrared colors that exhibit characteristics of both L and T dwarfs. Two of the new T dwarfs are likely binaries with L7+/-1 primaries and mid-type T secondaries. In addition, our follow-up program has confirmed 10 previously identified T dwarfs and four photometrically-selected L and T dwarf candidates in the literature. This sample, along with the previous WISE discoveries, triples the number of known brown dwarfs with spectral types later than T5. Using the WISE All-Sky Source Catalog we present updated color-color and color-type diagrams for all the WISE-discovered T and Y dwarfs. Near-infrared spectra of the new discoveries are presented, along with spectral classifications. To accommodate later T dwarfs we have modified the integrated flux method of determining spectral indices to instead use the median flux. Furthermore, a newly defined J-narrow index differentiates the early-type Y dwarfs from late-type T dwarfs based on the J-band continuum slope. The K/J indices for this expanded sample show that 32% of late-type T dwarfs have suppressed K-band flux and are blue relative to the spectral standards, while only 11% are redder than the standards. Comparison of the Y/J and K/J index to models suggests diverse atmospheric conditions and supports the possible re-emergence of clouds after the L/T transition. We also discuss peculiar brown dwarfs and candidates that were found not to be substellar, including two Young Stellar Objects and two Active Galactic Nuclei. The coolest WISE-discovered brown dwarfs are the closest of their type and will remain the only sample of their kind for many years to come.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Clouds and Hazes in Exoplanet Atmospheres

Clouds and hazes are commonplace in the atmospheres of solar system planets and are likely ubiquitous in the atmospheres of extrasolar planets as well. Clouds affect every aspect of a planetary atmosphere, from the transport of radiation, to atmospheric chemistry, to dynamics and they influence - if not control - aspects such as surface temperature and habitability. In this review we aim to provide an introduction to the role and properties of clouds in exoplanetary atmospheres. We consider the role clouds play in influencing the spectra of planets as well as their habitability and detectability. We briefly summarize how clouds are treated in terrestrial climate models and consider the far simpler approaches that have been taken so far to model exoplanet clouds, the evidence for which we also review. Since clouds play a major role in the atmospheres of certain classes of brown dwarfs we briefly discuss brown dwarf cloud modeling as well. We also review how the scattering and extinction efficiencies of cloud particles may be approximated in certain limiting cases of small and large particles in order to facilitate physical understanding. Since clouds play such important roles in planetary atmospheres, cloud modeling may well prove to be the limiting factor in our ability to interpret future observations of extrasolar planets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantitatively Assessing the Role of Clouds in the Transmission Spectrum of GJ 1214b

Recent observations of the super-Earth GJ 1214b show that it has a relatively featureless transmission spectrum. One suggestion is that these observations indicate that the planet's atmosphere is vertically compact, perhaps due to a water-rich composition that yields a large mean molecular weight. Another suggestion is that the atmosphere is hydrogen/helium-rich with clouds that obscure predicted absorption features. Previous models that incorporate clouds have included their effect without a strong physical motivation for their existence. Here, we present model atmospheres of GJ 1214b that include physically-motivated clouds of two types. We model the clouds that form as a result of condensation in chemical equilibrium, as they likely do on brown dwarfs, which include KCl and ZnS for this planet. We also include clouds that form as a result of photochemistry, forming a hydrocarbon haze layer. We use a photochemical kinetics model to understand the vertical distribution and available mass of haze-forming molecules. We model both solar and enhanced-metallicity cloudy models and determine the cloud properties necessary to match observations. In enhanced-metallicity atmospheres, we find that the equilibrium clouds can match the observations of GJ 1214b if they are lofted high into the atmosphere and have a low sedimentation efficiency (fsed=0.1). We find that models with a variety of hydrocarbon haze properties can match the observations. Particle sizes from 0.01 to 0.25 micron can match the transmission spectrum with haze-forming efficiencies as low as 1-5%.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Three-dimensional atmospheric circulation of hot Jupiters on highly eccentric orbits

Of the over 800 exoplanets detected to date, over half are on non-circular orbits, with eccentricities as high as 0.93. Such orbits lead to time-variable stellar heating, which has implications for the planet's atmospheric dynamical regime. However, little is known about this dynamical regime, and how it may influence observations. Therefore, we present a systematic study of hot Jupiters on highly eccentric orbits using the SPARC/MITgcm, a model which couples a three-dimensional general circulation model with a plane-parallel, two-stream, non-grey radiative transfer model. In our study, we vary the eccentricity and orbit-average stellar flux over a wide range. We demonstrate that the eccentric hot Jupiter regime is qualitatively similar to that of planets on circular orbits; the planets possess a superrotating equatorial jet and exhibit large day-night temperature variations. We show that these day-night heating variations induce momentum fluxes equatorward to maintain the superrotating jet throughout its orbit. As the eccentricity and/or stellar flux is increased, the superrotating jet strengthens and narrows, due to a smaller Rossby deformation radius. For a select number of model integrations, we generate full-orbit lightcurves and find that the timing of transit and secondary eclipse viewed from Earth with respect to periapse and apoapse can greatly affect what we see in infrared (IR) lightcurves; the peak in IR flux can lead or lag secondary eclipse depending on the geometry. For those planets that have large day-night temperature variations and rapid rotation rates, we find that the lightcurves exhibit "ringing" as the planet's hottest region rotates in and out of view from Earth. These results can be used to explain future observations of eccentric transiting exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Masses, Radii, and Cloud Properties of the HR 8799 Planets

The near-infrared colors of the planets directly imaged around the A star HR 8799 are much redder than most field brown dwarfs of the same effective temperature. Previous theoretical studies of these objects have concluded that the atmospheres of planets b, c, and d are unusually cloudy or have unusual cloud properties. Some studies have also found that the inferred radii of some or all of the planets disagree with expectations of standard giant planet evolution models. Here we compare the available data to the predictions of our own set of atmospheric and evolution models that have been extensively tested against observations of field L and T dwarfs, including the reddest L dwarfs. Unlike some previous studies we require mutually consistent choices for effective temperature, gravity, cloud properties, and planetary radius. This procedure thus yields plausible values for the masses, effective temperatures, and cloud properties of all three planets. We find that the cloud properties of the HR 8799 planets are not unusual but rather follow previously recognized trends, including a gravity dependence on the temperature of the L to T spectral transition--some reasons for which we discuss. We find the inferred mass of planet b is highly sensitive to whether or not we include the H and K band spectrum in our analysis. Solutions for planets c and d are consistent with the generally accepted constraints on the age of the primary star and orbital dynamics. We also confirm that, like in L and T dwarfs and solar system giant planets, non-equilibrium chemistry driven by atmospheric mixing is also important for these objects. Given the preponderance of data suggesting that the L to T spectral type transition is gravity dependent, we present an exploratory evolution calculation that accounts for this effect. Finally we recompute the the bolometric luminosity of all three planets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Neglected Clouds in T and Y Dwarf Atmospheres

As brown dwarfs cool, a variety of species condense in their atmospheres, forming clouds. Iron and silicate clouds shape the emergent spectra of L dwarfs, but these clouds dissipate at the L/T transition. A variety of other condensates are expected to form in cooler T dwarf atmospheres. These include Cr, MnS, Na2S, ZnS, and KCl, but the opacity of these optically thinner clouds has not been included in previous atmosphere models. Here, we examine their effect on model T and Y dwarf atmospheres. The cloud structures and opacities are calculated using the Ackerman & Marley (2001) cloud model, which is coupled to an atmosphere model to produce atmospheric pressure-temperature profiles in radiative-convective equilibrium. We generate a suite of models between Teff = 400 and 1300 K, log g=4.0 and 5.5, and condensate sedimentation efficiencies from fsed=2 to 5. Model spectra are compared to two red T dwarfs, Ross 458C and UGPS 0722-05; models that include clouds are found to match observed spectra significantly better than cloudless models. The emergence of sulfide clouds in cool atmospheres, particularly Na2S, may be a more natural explanation for the "cloudy" spectra of these objects, rather than the re-emergence of silicate clouds that wane at the L-to-T transition. We find that sulfide clouds provide a mechanism to match the near- and mid-infrared colors of observed T dwarfs. Our results indicate that including the opacity of condensates in T dwarf atmospheres is necessary to accurately determine the physical characteristics of many of the observed objects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Vertical Atmospheric Structure in a Variable Brown Dwarf: Pressure-dependent Phase Shifts in Simultaneous Hubble Space Telescope-Spitzer Light Curves

Heterogeneous clouds or temperature perturbations in rotating brown dwarfs produce variability in the observed flux. We report time-resolved simultaneous observations of the variable T6.5 brown dwarf 2MASSJ22282889-431026 over the wavelength ranges 1.1-1.7 microns and broadband 4.5 microns. Spectroscopic observations were taken with Wide Field Camera 3 on board the Hubble Space Telescope and photometry with the Spitzer Space Telescope. The object shows sinusoidal infrared variability with a period of 1.4 hr at most wavelengths with peak-to-peak amplitudes between 1.45% and 5.3% of the mean flux. While the light curve shapes are similar at all wavelengths, their phases differ from wavelength to wavelength with a maximum difference of more than half of a rotational period. We compare the spectra with atmospheric models of different cloud prescriptions, from which we determine the pressure levels probed at different wavelengths. We find that the phase lag increases with decreasing pressure level, or higher altitude. We discuss a number of plausible scenarios that could cause this trend of light curve phase with probed pressure level. These observations are the first to probe heterogeneity in an ultracool atmosphere in both horizontal and vertical directions, and thus are an ideal test case for realistic three dimensional simulations of the atmospheric structure with clouds in brown dwarfs and extrasolar planets.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Disequilibrium Carbon, Oxygen, and Nitrogen Chemistry in the Atmospheres of HD 189733b and HD 209458b

We have developed 1-D photochemical and thermochemical kinetics and diffusion models for the transiting exoplanets HD 189733b and HD 209458b to study the effects of disequilibrium chemistry on the atmospheric composition of "hot Jupiters." Here we investigate the coupled chemistry of neutral carbon, hydrogen, oxygen, and nitrogen species, and we compare the model results with existing transit and eclipse observations. We find that the vertical profiles of molecular constituents are significantly affected by transport-induced quenching and photochemistry, particularly on cooler HD 189733b; however, the warmer stratospheric temperatures on HD 209458b can help maintain thermochemical equilibrium and reduce the effects of disequilibrium chemistry. For both planets, the methane and ammonia mole fractions are found to be enhanced over their equilibrium values at pressures of a few bar to less than a mbar due to transport-induced quenching, but CH4 and NH3 are photochemically removed at higher altitudes. Atomic species, unsaturated hydrocarbons (particularly C2H2), some nitriles (particularly HCN), and radicals like OH, CH3, and NH2 are enhanced overequilibrium predictions because of quenching and photochemistry. In contrast, CO, H2O, N2, and CO2 more closely follow their equilibrium profiles, except at pressures < 1 microbar, where CO, H2O, and N2 are photochemically destroyed and CO2 is produced before its eventual high-altitude destruction. The enhanced abundances of HCN, CH4, and NH3 in particular are expected to affect the spectral signatures and thermal profiles HD 189733b and other, relatively cool, close-in transiting exoplanets. We examine the sensitivity of our results to the assumed temperature structure and eddy diffusion coefficientss and discuss further observational consequences of these models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Probing the Physical Properties of Directly Imaged Gas Giant Exoplanets Through Polarization

It is becoming clear that the atmospheres of the young, self-luminous extrasolar giant planets imaged to date are dusty. Planets with dusty atmospheres may exhibit detectable amounts of linear polarization in the near-infrared, as has been observed from some field L dwarfs. The asymmetry required in the thermal radiation field to produce polarization may arise either from the rotation-induced oblateness or from surface inhomogeneities, such as partial cloudiness. While it is not possible at present to predict the extent to which atmospheric dynamics on a given planet may produce surface inhomogeneities substantial enough to produce net non-zero disk integrated polarization, the contribution of rotation-induced oblateness can be estimated. Using a self-consistent, spatially homogeneous atmospheric model and a multiple scattering polarization formalism for this class of exoplanets, we show that polarization on the order of 1% may arise due to the rotation-induced oblateness of the planets. The degree of polarization for cloudy planets should peak at the same wavelengths at which the planets are brightest in the near-infrared. The observed polarization may be even higher if surface inhomogeneities exist and play a significant role. Polarized radiation from self-luminous gas giant exoplanets, if detected, provides an additional tool to characterize these young planets and a new method to constrain their surface gravity and masses.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Discovery of Y Dwarfs Using Data from the Wide-field Infrared Survey Explorer (WISE)

We present the discovery of seven ultracool brown dwarfs identified with the Wide-field Infrared Survey Explorer (WISE). Near-infrared spectroscopy reveals deep absorption bands of H_2O and CH_4 that indicate all seven of the brown dwarfs have spectral types later than UGPS J072227.51-054031.2, the latest type T dwarf currently known. The spectrum of WISEP J182831.08+265037.8 is distinct in that the heights of the J- and H-band peaks are approximately equal in units of f_lambda, so we identify it as the archetypal member of the Y spectral class. The spectra of at least two of the other brown dwarfs exhibit absorption on the blue wing of the H-band peak that we tentatively ascribe to NH_3. These spectral morphological changes provide a clear transition between the T dwarfs and the Y dwarfs. In order to produce a smooth near-infrared spectral sequence across the T/Y dwarf transition, we have reclassified UGPS J0722-0540 as the T9 spectral standard and tentatively assign WISEP J173835.52+273258.9 as the Y0 spectral standard. In total, six of the seven new brown dwarfs are classified as Y dwarfs: four are classified as Y0, one is classified as Y0 (pec?), and WISEP J1828+2650 is classified as >Y0. We have also compared the spectra to the model atmospheres of Marley and Saumon and infer that the brown dwarfs have effective temperatures ranging from 300 K to 500 K, making them the coldest spectroscopically confirmed brown dwarfs known to date.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Patchy Cloud Model for the L to T Dwarf Transition

One of the mechanisms suggested for the L to T dwarf spectral type transition is the appearance of relatively cloud-free regions across the disk of brown dwarfs as they cool. The existence of partly cloudy regions has been supported by evidence for variability in dwarfs in the late L to early T spectral range, but no self-consistent atmosphere models of such partly cloudy objects have yet been constructed. Here we present a new approach for consistently modeling partly cloudy brown dwarfs and giant planets. We find that even a small fraction of cloud holes dramatically alter the atmospheric thermal profile, spectra, and photometric colors of a given object. With decreasing cloudiness objects briskly become bluer in J - K and brighten in J band, as is observed at the L/T transition. Model spectra of partly cloudy objects are similar to our models with globally homogenous, but thinner, clouds. Hence spectra alone may not be sufficient to distinguish partial cloudiness although variability and polarization measurements are potential observational signatures. Finally we note that partial cloud cover may be an alternative explanation for the blue L dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Atmospheric Circulation of Eccentric Hot Neptune GJ436b

GJ436b is a unique member of the transiting extrasolar planet population being one of the smallest and least irradiated and possessing an eccentric orbit. Because of its size, mass and density, GJ436b could plausibly have an atmospheric metallicity similar to Neptune (20-60 times solar abundances), which makes it an ideal target to study the effects of atmospheric metallicity on dynamics and radiative transfer in an extrasolar planetary atmosphere. We present three-dimensional atmospheric circulation models that include realistic non-gray radiative transfer for 1, 3, 10, 30, and 50 times solar atmospheric metallicity cases of GJ436b. Low metallicity models (1 and 3 times solar) show little day/night temperature variation and strong high-latitude jets. In contrast, higher metallicity models (30 and 50 times solar) exhibit day/night temperature variations and a strong equatorial jet. Spectra and light curves produced from these simulations show strong orbital phase dependencies in the 50 times solar case and negligible variations with orbital phase in the 1 times solar case. Comparisons between the predicted planet/star flux ratio from these models and current secondary eclipse measurements support a high metallicity atmosphere (30-50 times solar abundances) with disequilibrium carbon chemistry at play for GJ436b. Regardless of the actual atmospheric composition of GJ436b, our models serve to illuminate how metallicity influences the atmospheric circulation for a broad range of warm extrasolar planets.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Clouds in the Coldest Brown Dwarfs: FIRE Spectroscopy of Ross 458C

Condensate clouds are a salient feature of L dwarf atmospheres, but have been assumed to play little role in shaping the spectra of the coldest T-type brown dwarfs. Here we report evidence of condensate opacity in the near-infrared spectrum of the brown dwarf candidate Ross 458C, obtained with the Folded-Port Infrared Echellette (FIRE) spectrograph at the Magellan Telescopes. These data verify the low-temperature nature of this source, indicating a T8 spectral classification, log Lbol/Lsun = -5.62+/-0.03, Teff = 650+/-25 K, and a mass at or below the deuterium burning limit. The data also reveal enhanced emission at K-band associated with youth (low surface gravity) and supersolar metallicity, reflecting the properties of the Ross 458 system (age = 150-800 Myr, [Fe/H] = +0.2 to +0.3). We present fits of FIRE data for Ross 458C, the T9 dwarf ULAS J133553.45+113005.2, and the blue T7.5 dwarf SDSS J141624.08+134826.7B, to cloudless and cloudy spectral models from Saumon & Marley. For Ross 458C we confirm a low surface gravity and supersolar metallicity, while the temperature differs depending on the presence (635 [+25,-35] K) or absence (760 [+70,-45] K) of cloud extinction. ULAS J1335+1130 and SDSS J1416+1348B have similar temperatures (595 [+25,-45] K), but distinct surface gravities (log g = 4.0-4.5 cgs versus 5.0-5.5 cgs) and metallicities ([M/H] ~ +0.2 versus -0.2). In all three cases, cloudy models provide better fits to the spectral data, significantly so for Ross 458C. These results indicate that clouds are an important opacity source in the spectra of young cold T dwarfs, and should be considered when characterizing the spectra of planetary-mass objects in young clusters and directly-imaged exoplanets. The characteristics of Ross 458C suggest it could itself be regarded as a planet, albeit one whose cosmogony does not conform with current planet formation theories.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exoplanet albedo spectra and colors as a function of planet phase, separation, and metallicity

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 μ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 μ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 μ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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Observed Polarization of Brown Dwarfs Suggests Low Surface Gravity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

SDSS J141624.08+134826.7: Blue L Dwarfs and Non-Equilibrium Chemistry

We present an analysis of the recently discovered blue L dwarf SDSS J141624.08+134826.7. We extend the spectral coverage of its published spectrum to ~4 microns by obtaining a low-resolution L band spectrum with SpeX on the NASA IRTF. The spectrum exhibits a tentative weak CH4 absorption feature at 3.3 microns but is otherwise featureless. We derive the atmospheric parameters of SDSS J141624.08+134826.7 by comparing its 0.7-4.0 micron spectrum to the atmospheric models of Marley and Saumon which include the effects of both condensate cloud formation and non-equilibrium chemistry due to vertical mixing and find the best fitting model has Teff=1700 K, log g=5.5 [cm s-2], fsed=4, and Kzz=10^4 cm2 s-1. The derived effective temperature is significantly cooler than previously estimated but we confirm the suggestion by Bowler et al. that the peculiar spectrum of SDSS J141624.08+134826.7 is primarily a result of thin condensate clouds. In addition, we find strong evidence of vertical mixing in the atmosphere of SDSS J141624.08+134826.7 based on the absence of the deep 3.3 micron CH4 absorption band predicted by models computed in chemical equilibrium. This result suggests that observations of blue L dwarfs are an appealing way to quantitatively estimate the vigor of mixing in the atmospheres of L dwarfs because of the dramatic impact such mixing has on the strength of the 3.3 micron CH4 band in the emergent spectra of L dwarfs with thin condensate clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Atmospheric circulation of hot Jupiters: Coupled radiative-dynamical general circulation model simulations of HD 189733b and HD 209458b

We present global, three-dimensional numerical simulations of HD 189733b and HD 209458b that couple the atmospheric dynamics to a realistic representation of non-gray cloud-free radiative transfer. The model, which we call the Substellar and Planetary Atmospheric Radiation and Circulation (SPARC) model, adopts the MITgcm for the dynamics and uses the radiative model of McKay, Marley, Fortney, and collaborators for the radiation. Like earlier work with simplified forcing, our simulations develop a broad eastward equatorial jet, mean westward flow at higher latitudes, and substantial flow over the poles at low pressure. For HD 189733b, our simulations without TiO and VO opacity can explain the broad features of the observed 8 and 24-micron light curves, including the modest day-night flux variation and the fact that the planet/star flux ratio peaks before the secondary eclipse. Our simulations also provide reasonable matches to the Spitzer secondary-eclipse depths at 4.5, 5.8, 8, 16, and 24 microns and the groundbased upper limit at 2.2 microns. However, we substantially underpredict the 3.6-micron secondary-eclipse depth, suggesting that our simulations are too cold in the 0.1-1 bar region. Predicted temporal variability in secondary-eclipse depths is ~1% at Spitzer bandpasses, consistent with recent observational upper limits at 8 microns. We also show that nonsynchronous rotation can significantly alter the jet structure. For HD 209458b, we include TiO and VO opacity; these simulations develop a hot (>2000 K) dayside stratosphere. Despite this stratosphere, we do not reproduce current Spitzer photometry of this planet. Light curves in Spitzer bandpasses show modest phase variation and satisfy the observational upper limit on day-night phase variation at 8 microns. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Multiple Scattering Polarization of Substellar-mass Objects : T-dwarfs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Future of Ultracool Dwarf Science with JWST

Ultracool dwarfs exhibit a remarkably varied set of characteristics which hint at the complex physical processes acting in their atmospheres and interiors. Spectra of these objects not only depend upon their mass and effective temperature, but also their atmospheric chemistry, weather, and dynamics. As a consequence divining their mass, metallicity and age solely from their spectra has been a challenge. JWST, by illuminating spectral blind spots and observing objects with constrained masses and ages should finally unearth a sufficient number of ultracool dwarf Rosetta Stones to allow us to decipher the processes underlying the complex brown dwarf cooling sequence. In addition the spectra of objects invisible from the ground, including very low mass objects in clusters and nearby cold dwarfs from the disk population, will be seen for the first time. In combination with other ground- and space-based assets and programs, JWST will usher in a new golden era of brown dwarf science and discovery.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A Unified Theory for the Atmospheres of the Hot and Very Hot Jupiters: Two Classes of Irradiated Atmospheres

We highlight the importance of gaseous TiO and VO opacity on the highly irradiated close-in giant planets. The atmospheres of these planets naturally fall into two classes that are somewhat analogous to the M- and L-type dwarfs. Those that are warm enough to have appreciable opacity due to TiO and VO gases we term the ``pM Class'' planets, and those that are cooler we term ``pL Class'' planets. We calculate model atmospheres for these planets, including pressure-temperature profiles, spectra, and characteristic radiative time constants. We show that pM Class planets have hot stratospheres $\sim$2000 K and appear ``anomalously'' bright in the mid infrared secondary eclipse, as was recently found for planets HD 149026b and HD 209458b. This class of planets absorbs incident flux and emits thermal flux from high in their atmospheres. Consequently, they will have large day/night temperature contrasts and negligible phase shifts between orbital phase and thermal emission light curves, because radiative timescales are much shorter than possible dynamical timescales. The pL Class planets absorb incident flux deeper in the atmosphere where atmospheric dynamics will more readily redistribute absorbed energy. This will lead to cooler day sides, warmer night sides, and larger phase shifts in thermal emission light curves. Around a Sun-like primary this boundary occurs at $\sim$0.04-0.05 AU. The eccentric transiting planets HD 147506b and HD 17156b alternate between the classes. Thermal emission in the optical from pM Class planets is significant red-ward of 400 nm, making these planets attractive targets for optical detection. The difference in the observed day/night contrast between ups Andromeda b (pM Class) and HD 189733b (pL Class) is naturally explained in this scenario. (Abridged.)

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Atmospheric Parameters of Field L and T Dwarfs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Moderate Resolution Spitzer Infrared Spectrograph (IRS) Observations of M, L, and T Dwarfs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Planetary Radii across Five Orders of Magnitude in Mass and Stellar Insolation: Application to Transits

To aid in the physical interpretation of planetary radii constrained through observations of transiting planets, or eventually direct detections, we compute model radii of pure hydrogen-helium, water, rock, and iron planets, along with various mixtures. Masses ranging from 0.01 Earth masses to 10 Jupiter masses at orbital distances of 0.02 to 10 AU are considered. For hydrogen-helium rich planets, our models are the first to couple planetary evolution to stellar irradiation over a wide range of orbital separations (0.02 to 10 AU) through a non-gray radiative-convective equilibrium atmosphere model. Stellar irradiation retards the contraction of giant planets, but its effect is not a simple function of the irradiation level: a planet at 1 AU contracts as slowly as a planet at 0.1 AU. For hydrogen-helium planets, we consider cores up to 90% of the total planet mass, comparable to those of Uranus and Neptune. If "hot Neptunes" have maintained their original masses and are not remnants of more massive planets, radii of 0.30-0.45 times Jupiter's radius are expected. Water planets are ~40-50% larger than rocky planets, independent of mass. Finally, we provide tables of planetary radii at various ages and compositions, and for ice-rock-iron planets we fit our results to analytic functions, which will allow for quick composition estimates, given masses and radii, or mass estimates, given only planetary radii. These results will assist in the interpretation of observations for both the current transiting planet surveys as well as upcoming space missions, including CoRoT and Kepler.

0 citations0 reviewsNo code linkedOpen access