Search

Measured rotational speeds of giant planets and brown dwarfs frequently constitute appreciable fractions of the breakup limit, resulting in centrifugal expansion of these objects at the equator. According to models of internal energy transport, this expansion ought to make the poles of a rotator significantly hotter than the equator, so that inclination of the rotational axis greatly affects both spectral shape and total flux. In this article, we explore the dependence of a substellar object's observables on its rotational speed and axis inclination. To do so, we combine PICASO (a Planetary Intensity Code for Atmospheric Spectroscopy Observations) with software PARS (Paint the Atmospheres of Rotating Stars). The former computer program models radiative transfer within plane-parallel planetary atmospheres, while the latter computes disk-integrated spectra of centrifugally deformed gaseous masses. We find that the specific flux of a typical fast-rotating brown dwarf can increase by as much as a factor of 1.5 with movement from an equator-on to a pole-on view. On the other hand, the distinctive effect of rotation on spectral shape increases toward the equator-on view. The latter effect also increases with lower effective temperature. The bolometric luminosity estimate for a typical fast rotator at extreme inclinations has to be adjusted by as much as about 20% due to the anisotropy of the object's observed flux. We provide a general formula for the calculation of the corresponding adjustment factor in terms of rotational speed and inclination., 11 pages, 4 figures, MNRAS accepted

Planets orbiting M-dwarf stars are prime targets in the search for rocky exoplanet atmospheres. The small size of M dwarfs renders their planets exceptional targets for transmission spectroscopy, facilitating atmospheric characterization. However, it remains unknown whether their host stars' highly variable extreme-UV radiation environments allow atmospheres to persist. With JWST, we have begun to determine whether or not the most favorable rocky worlds orbiting M dwarfs have detectable atmospheres. Here, we present a 2.8-5.2 micron JWST NIRSpec/G395H transmission spectrum of the warm (700 K, 40.3x Earth's insolation) super-Earth GJ 486b (1.3 R$_{\oplus}$ and 3.0 M$_{\oplus}$). The measured spectrum from our two transits of GJ 486b deviates from a flat line at 2.2 - 3.3 $\sigma$, based on three independent reductions. Through a combination of forward and retrieval models, we determine that GJ 486b either has a water-rich atmosphere (with the most stringent constraint on the retrieved water abundance of H2O > 10% to 2$\sigma$) or the transmission spectrum is contaminated by water present in cool unocculted starspots. We also find that the measured stellar spectrum is best fit by a stellar model with cool starspots and hot faculae. While both retrieval scenarios provide equal quality fits ($\chi^2_\nu$ = 1.0) to our NIRSpec/G395H observations, shorter wavelength observations can break this degeneracy and reveal if GJ 486b sustains a water-rich atmosphere., Comment: 18 pages, 7 figures, 5 tables. Accepted in ApJ Letters. Co-First Authors

Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength coverage, moderate spectral resolution, and high precision that, together, are not achievable with previous observatories. Now that JWST has commenced science operations, we are able to observe exoplanets at previously uncharted wavelengths and spectral resolutions. Here we report time-series observations of the transiting exoplanet WASP-39b using JWST's Near InfraRed Camera (NIRCam). The long-wavelength spectroscopic and short-wavelength photometric light curves span 2.0 - 4.0 $\mu$m, exhibit minimal systematics, and reveal well-defined molecular absorption features in the planet's spectrum. Specifically, we detect gaseous H$_2$O in the atmosphere and place an upper limit on the abundance of CH$_4$. The otherwise prominent CO$_2$ feature at 2.8 $\mu$m is largely masked by H$_2$O. The best-fit chemical equilibrium models favour an atmospheric metallicity of 1-100$\times$ solar (i.e., an enrichment of elements heavier than helium relative to the Sun) and a sub-stellar carbon-to-oxygen (C/O) ratio. The inferred high metallicity and low C/O ratio may indicate significant accretion of solid materials during planet formation or disequilibrium processes in the upper atmosphere., Comment: 35 pages, 13 figures, 3 tables, Nature, accepted

Exoplanet atmospheric retrieval is a computational technique widely used to infer properties of planetary atmospheres from remote spectroscopic observations. Retrieval codes typically employ Bayesian sampling algorithms or machine learning approaches to explore the range of atmospheric properties (e.g., chemical composition, temperature structure, aerosols) compatible with an observed spectrum. However, despite the wide adoption of exoplanet retrieval techniques, there is currently no systematic summary of exoplanet retrieval codes in the literature. Here, we provide a catalogue of the atmospheric retrieval codes published to date, alongside links to their respective code repositories where available. Our catalogue will be continuously updated via a Zenodo archive., Comment: 5 pages, 1 giant Table. Published in RNAAS. Live catalogue will be updated at https://doi.org/10.5281/zenodo.7675743

Evidence of disequilibrium chemistry due to vertical mixing in the atmospheres of many T and Y-dwarfs has been inferred due to enhanced mixing ratios of CO and reduced NH$_3$. Atmospheric models of planets and brown dwarfs typically parameterize this vertical mixing phenomenon with the vertical eddy diffusion coefficient, $K_{\rm zz}$. While $K_{\rm zz}$ can perhaps be approximated in the convective regions in the atmosphere with mixing length theory, in radiative regions the strength of vertical mixing is uncertain by many orders of magnitude. With a new grid of self-consistent 1D model atmospheres from $T_{\rm eff}$ of 400 - 1000 K, computed with a new radiative-convective equilibrium python code PICASO 3.0, we aim to assess how molecular abundances and corresponding spectra can be used as a probe of depth-dependent $K_{\rm zz}$. At a given surface gravity, we find non-monotonic behavior in the CO abundance as a function of $T_{\rm eff}$, as chemical abundances are sometimes quenched in either of two potential atmospheric convective zones, or quenched in either of two possible radiative zones. The temperature structure and chemical quenching behavior also changes with gravity. We compare our models with available near-infrared and M-band spectroscopy of several T and Y-dwarfs and assess their atmospheric vertical mixing profiles. We also compare to color-magnitude diagrams and make predictions for JWST spectra. This work yields new constraints, and points the way to significant future gains, in determining $K_{\rm zz}$, a fundamental atmospheric parameter in substellar atmospheres, with significant implications for chemistry and cloud modeling., Comment: Accepted for Publication in The Astrophysical Journal, 31 Pages, 22 Figures. The model grid will be released via Zenodo

We present new functionality within PICASO, a state-of-the-art radiative transfer model for exoplanet and brown dwarf atmospheres, by developing a new pipeline that computes phase-resolved thermal emission (thermal phase curves) from three-dimensional (3D) models. Because PICASO is coupled to Virga, an open-source cloud code, we are able to produce cloudy phase curves with different sedimentation efficiencies ($f_{sed}$) and cloud condensate species. We present the first application of this new algorithm to hot Jupiter WASP-43b. Previous studies of the thermal emission of WASP-43b from Kataria et al. found good agreement between cloud-free models and dayside thermal emission, but an overestimation of the nightside flux, for which clouds have been suggested as a possible explanation. We use the temperature and vertical wind structure from the cloud-free 3D general circulation models of Kataria et al. and post-process it using PICASO, assuming that clouds form and affect the spectra. We compare our models to results from Kataria et al., including Hubble Space Telescope Wide-Field Camera 3 (WFC3) observations of WASP-43b from Stevenson et al. In addition, we compute phase curves for Spitzer at 3.6 and 4.5 $\mu m$ and compare them to observations from Stevenson et al. We are able to closely recover the cloud-free results, even though PICASO utilizes a coarse spatial grid. We find that cloudy phase curves provide much better agreement with the WFC3 and Spitzer nightside data, while still closely matching the dayside emission. This work provides the community with a convenient, user-friendly tool to interpret phase-resolved observations of exoplanet atmospheres using 3D models., Comment: Published ApJ; GitHub: https://github.com/natashabatalha/picaso/blob/master/picaso/ Docs and Tutorials: https://natashabatalha.github.io/picaso/tutorials.html#moving-to-3-dimensions

TOI-1266c is a recently discovered super-Venus in the radius valley orbiting an early M dwarf. However, its notional bulk density ($\sim$2.2 g cm$^{-3}$) is consistent with a large volatile fraction, suggesting that it might have volatile reservoirs that have survived billions of years at more than twice the Earth's insolation. On the other hand, the upper mass limit paints a picture of a cool super Mercury dominated by >50\% iron core ($\sim$9.2 g cm$^{-3}$) that has tiptoed up to the collisional stripping limit and into the radius gap. Here, we examine several hypothetical states for TOI-1266c using a combination of new and updated open-source atmospheric escape, radiative-convective, and photochemical models. We find that water-rich atmospheres with trace amounts of H$_{2}$ and CO$_{2}$ are potentially detectable (SNR $>\sim 5$) in less than 20 hours of JWST observing time. We also find that water vapor spectral features are not substantially impacted by the presence of high-altitude water or ice clouds due the presence of a significant amount of water above the cloud-deck, although further work with self-consistent cloud models is needed. Regardless of its mass, however, TOI-1266c represents a unique proving ground for several hypotheses related to the evolution of sub-Neptunes and Venus-like worlds, particularly those near the radius valley., 40 pages (18 pages in the main text), 15 figures (10 figures in the main text), 2 appendices. Revised based on first round of reviews at PSJ. Comments welcome

The majority of exoplanets found to date have been discovered via the transit method, and transmission spectroscopy represents the primary method of studying these distant worlds. Currently, in-depth atmospheric characterization of transiting exoplanets entails the use of spectrographs on large telescopes, requiring significant observing time to study each planet. Previous studies have demonstrated trends for solar system worlds using color–color photometry of reflectance spectra, as well as trends within transmission spectra for hot Jupiters. Building on these concepts, we have investigated the use of transmission color photometric analysis for efficient, coarse categorization of exoplanets and for assessing the nature of these worlds, with a focus on resolving the bulk composition degeneracy to aid in discriminating super-Earths and sub-Neptunes. We present our methodology and first results, including spectrum models, model comparison frameworks, and wave band selection criteria. We present our results for different transmission “color” metrics, filter selection methods, and numbers of filters. Assuming noise-free spectra of isothermal atmospheres in chemical equilibrium, with our pipeline, we are able to constrain atmospheric mean molecular weight in order to distinguish between super-Earth and sub-Neptune atmospheres with >90% overall accuracy using specific low-resolution filter combinations, . We also found that increasing the number of filters does not substantially impact this performance. This method could allow for broad characterization of large numbers of planets much more efficiently than current methods permit, enabling population- and system-level studies. Additionally, data collected via this method could inform follow-up observing time by large telescopes for more detailed studies of worlds of interest.

Lithium is an important element for the understanding of ultracool dwarfs because it is lost to fusion at masses above $\sim 68\, M_{\rm J}$. Hence, the presence or absence of atomic Li has served as an indicator of the nearby H-burning boundary at about $75\,M_{\rm J}$ between brown-dwarfs and very low-mass stars. Historically the "Lithium test", a search for the presence and strength of the Li line at 670.8 nm, has been a marker if an object has a substellar mass with stellar-like spectral energy distribution (e.g., a late-type M dwarf). While the Li test could in principle also be used to distinguish masses of later-type L-T dwarfs, Li is predominantly no longer found as an atomic gas, but rather a molecular species such as LiH, LiF, LiOH, and LiCl in their cooler atmospheres. L- and T-type brown dwarfs are also quite faint at 670 nm and thus challenging targets for high resolution spectroscopy. But only recently have experimental molecular line lists become available for the molecular Li species, allowing molecular Li mass discrimination. In this study, we generated the latest opacity of each of these Li-bearing molecules and performed thermochemical equilibrium atmospheric composition calculation of the abundance of these molecules. Finally, we computed thermal emission spectra for a series of radiative-convective equilibrium models of cloudy and cloudless brown dwarf atmospheres (with $T_{\rm eff}=$ 500--2400~K, and $\log g$=4.0, 4.5, 5.0) to understand where the presence or absence of atmospheric lithium-bearing species is most easily detected as a function of brown dwarf mass and age. After atomic Li, the best spectral signatures were found to be LiF at $10.5-12.5$~\micron and LiCl at $14.5-18.5$ $\micron$. LiH also shows a narrow feature at $\sim 9.38$ $\micron$., 13 pages, 7 figures, 1 table, In final revision at ApJ. Comments welcome

Optical secondary eclipse measurements made by \emph{Kepler} reveal a diverse set of geometric albedos for hot Jupiters with equilibrium temperatures between $1550-1700$ K. The presence or absence of high altitude condensates, such as Mg$_2$SiO$_4$, Fe, Al$_2$O$_3$, and TiO$_2$, can significantly alter optical albedos, but these clouds are expected to be confined to localized regions in the atmospheres of these tidally locked planets. Here, we present 3D general circulation models and corresponding cloud and albedo maps for six hot Jupiters with measured optical albedos in this temperature range. We find that the observed optical albedos of K2-31b and K2-107b are best matched by either cloud free models or models with relatively compact cloud layers, while Kepler-8b and Kepler-17b's optical albedos can be matched by moderately extended ($f_{sed}$ = 0.1) parametric cloud models. HATS-11b has a high optical albedo, corresponding to models with bright Mg$_2$SiO$_4$ clouds extending to very low pressures ($f_{sed}$ = 0.03). We are unable to reproduce Kepler-7b's high albedo, as our models predict that the dayside will be dominated by dark Al$_2$O$_3$ clouds at most longitudes. We compare our parametric cloud model with a two-zone microphysical cloud model (\texttt{CARMA}). We find that even after accounting for the 3D thermal structure, no single cloud model can explain the full range of observed albedos within the sample. We conclude that a better knowledge of the vertical mixing profiles, cloud radiative feedback, cloud condensate properties, and atmospheric metallicities is needed in order to explain the unexpected diversity of albedos in this temperature range., Comment: accepted, ApJ

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., Accepted ApJ; GitHub: http://github.com/natashabatalha/virga Docs: http://natashabatalha.github.io/virga/# Virga Tutorials: http://natashabatalha.github.io/virga/tutorials.html

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

The advantage of having a high-fidelity instrument simulation tool developed in tandem with novel instrumentation is having the ability to investigate, in isolation and in combination, the wide parameter space set by the instrument design. SCALES, the third generation thermal-infrared diffraction limited imager and low/med-resolution integral field spectrograph being designed for Keck, is an instrument unique in design in order to optimize for its driving science case of direct detection and characterization of thermal emission from cold exoplanets. This warranted an end-to-end simulation tool that systematically produces realistic mock data from SCALES to probe the recovery of injected signals under changes in instrument design parameters. In this paper, we quantify optomechanical tolerance and detector electronic requirements set by the fiducial science cases using information content analysis, and test the consequences of updates to the design of the instrument on meeting these requirements., Comment: 7 pages, 4 figures, SPIE Optical Engineering and Applications 2021

The exoplanet atmosphere characterization goals of future UV/Optical/Infrared flagship space missions will drive challenging design requirements for instrument wavefront controls, spatial and spectral sampling, spectral bandwidth, and detector performance. The new ExoSpec Project links four previously distinct research efforts at Goddard and Ames for enabling and enhancing the characterization of directly-imaged exoplanets. ExoSpec is comprised of three laboratory subsystem demonstrations: high-contrast integral field spectrographs, p-channel CCDs, and parabolic deformable mirrors. A fourth component, exoplanet spectral retrieval, is an iterative data simulation activity driven by the goal of optimizing the system-level instrument design for atmosphere retrieval metrics. The ExoSpec Project's emphasis on system-level spectroscopy performance complements the objectives of other technology demonstrations supported by NASA.

For solar-system objects, ultraviolet spectroscopy has been critical in identifying sources for stratospheric heating and measuring the abundances of a variety of hydrocarbon and sulfur-bearing species, produced via photochemical mechanisms, as well as oxygen and ozone. To date, less than 20 exoplanets have been probed in this critical wavelength range (0.2-0.4 um). Here we use data from Hubble's newly implemented WFC3 UVIS G280 grism to probe the atmosphere of the hot Jupiter HAT-P-41b in the ultraviolet through optical in combination with observations at infrared wavelengths. We analyze and interpret HAT-P-41b's 0.2-5.0 um transmission spectrum using a broad range of methodologies including multiple treatments of data systematics as well as comparisons with atmospheric forward, cloud microphysical, and multiple atmospheric retrieval models. Although some analysis and interpretation methods favor the presence of clouds or potentially a combination of Na, VO, AlO, and CrH to explain the ultraviolet through optical portions of HAT-P-41b's transmission spectrum, we find that the presence of a significant H- opacity provides the most robust explanation. We obtain a constraint for the abundance of H-, log(H-) = -8.65 +/- 0.62 in HAT-P-41b's atmosphere, which is several orders of magnitude larger than predictions from equilibrium chemistry for a 1700 - 1950 K hot Jupiter. We show that a combination of photochemical and collisional processes on hot hydrogen-dominated exoplanets can readily supply the necessary amount of H- and suggest that such processes are at work in HAT-P-41b and many other hot Jupiter atmospheres., Comment: 17 pages, 7 figures. Published in ApJL (October 8th 2020)

HD 21749 is a bright ($V=8.1$ mag) K dwarf at 16 pc known to host an inner terrestrial planet HD 21749c as well as an outer sub-Neptune HD 21749b, both delivered by TESS. Follow-up spectroscopic observations measured the mass of HD 21749b to be $22.7\pm2.2\ M_{\oplus}$ with a density of $7.0^{+1.6}_{-1.3}$ g~cm$^{-3}$, making it one of the densest sub-Neptunes. However, the mass measurement was suspected to be influenced by stellar rotation. Here we present new high-cadence PFS RV data to disentangle the stellar activity signal from the planetary signal. We find that HD 21749 has a similar rotational timescale as the planet's orbital period, and the amplitude of the planetary orbital RV signal is estimated to be similar to that of the stellar activity signal. We perform Gaussian Process (GP) regression on the photometry and RVs from HARPS and PFS to model the stellar activity signal. Our new models reveal that HD 21749b has a radius of $2.86\pm0.20\ R_{\oplus}$, an orbital period of $35.6133\pm0.0005$ d with a mass of $M_{b}=20.0\pm2.7\ M_{\oplus}$ and a density of $4.8^{+2.0}_{-1.4}$ g~cm$^{-3}$ on an eccentric orbit with $e=0.16\pm0.06$, which is consistent with the most recent values published for this system. HD 21749c has an orbital period of $7.7902\pm0.0006$ d, a radius of $1.13\pm0.10\ R_{\oplus}$, and a 3$\sigma$ mass upper limit of $3.5\ M_{\oplus}$. Our Monte Carlo simulations confirm that without properly taking stellar activity signals into account, the mass measurement of HD 21749b is likely to arrive at a significantly underestimated error bar., Comment: 20 pages, 16 figures, accepted to MNRAS

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

Stellar, substellar, and planetary atmosphere models are all highly sensitive to the input opacities. Generational differences between various state-of-the-art stellar/planetary models are primarily because of incomplete and outdated atomic/molecular line-lists. Here we present a database of pre-computed absorption cross-sections for all isotopologues of key atmospheric molecules relevant to late-type stellar, brown dwarf, and planetary atmospheres: MgH, AlH, CaH, TiH, CrH, FeH, SiO, TiO, VO, and H2O. The pressure and temperature ranges of the computed opacities are between 10$^{-6}$--3000~bar and 75--4000~K, and their spectral ranges are 0.25--330~$\mu$m for many cases where possible. For cases with no pressure-broadening data, we use collision theory to bridge the gap. We also probe the effect of absorption cross-sections calculated from different line lists in the context of Ultra-Hot Jupiter and M-dwarf atmospheres. Using 1-D self-consistent radiative-convective thermochemical equilibrium models, we report significant variations in the theoretical spectra and thermal profiles of substellar atmospheres. With a 2000 K representative Ultra-Hot Jupiter, we report variations of up to 320 and 80 ppm in transmission and thermal emission spectra, respectively. For a 3000 K M-dwarf, we find differences of up to 125$\%$ in the spectra. We find that the most significant differences arise due to the choice of TiO line-lists, primarily below 1$\mu$m. In sum, we present (1) a database of pre-computed molecular absorption cross-sections, and (2) quantify biases that arise when characterizing substellar/exoplanet atmospheres due to line list differences, therefore highlighting the importance of correct and complete opacities for eventual applications to high precision spectroscopy and photometry., Comment: accepted in ApJS. 44 pages, 17 figures, 7 tables. Comments are welcome

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraphic Instrument (CGI) technology demonstration and potential science mission, as well as mission concepts with exoplanet imaging capabilities such as HabEx and LUVOIR, all require the identification of the best targets for exoplanet observations. To date, the focus has been primarily on two classes of targets: those with known exoplanets discovered indirectly that may be observable by these imaging missions, and those targets with no known planets that have high completeness values (probabilities of planet detections) under some assumptions of the instrument performance and the overall population of exoplanets. A third class of target, however, has received much less scrutiny: stars with known exoplanets that could not possibly be directly imaged due to the size of their orbits. These are planets that would be guaranteed to spend all of their time inside the inner working angles or outside the outer working angles of all currently proposed coronagraphs. However, these systems could potentially harbor additional planets, either exterior to or interior to the currently known planets, but not yet detectable by indirect means. Here, we discuss how to identify systems from all three categories that would be good targets for WFIRST and other, future, space-based imagers. We present a method for assessing the utility of these targets based on an exploration of the available dynamical phase space of the systems that would result in long-term stable orbits for both the currently known and the potentially discoverable companions, and show how it augments existing methods for assessing target completeness and utility.

Vortex fiber nulling (VFN) is a method that may enable the detection and characterization of exoplanets at small angular separations (0.5-2 $\lambda/D$) with ground- and space-based telescopes. Since the field of view is within the inner working angle of most coronagraphs, nulling accesses non-transiting planets that are otherwise too close to their star for spectral characterization by other means, thereby significantly increasing the number of known exoplanets available for direct spectroscopy in the near-infrared. Furthermore, VFN targets planets on closer-in orbits which tend to have more favorable planet-to-star flux ratios in reflected light. Here, we present the theory and applications of VFN, show that the optical performance is approximately equivalent for a variety of implementations and aperture shapes, and discuss the trade-offs between throughput and engineering requirements using numerical simulations. We compare vector and scalar approaches and, finally, show that beam shaping optics may be used to significantly improve the throughput for planet light. Based on theoretical performance, we estimate the number of known planets and theoretical exoEarths accessible with a VFN instrument linked to a high-resolution spectrograph on the future Thirty Meter Telescope., Comment: To appear in Proceedings of the SPIE, Techniques and Instrumentation for Detection of Exoplanets IX, Vol. 11117

The reflected light signal from a planet throughout its orbit is a powerful probe of a planet's atmospheric properties. There are a number of planets that are amenable to reflected light phase curve studies with present and future space-based instrumentation and here we assess our ability to characterize these worlds. Using simulated TESS populations we identify the Nine, a set of archetypal exoplanets with the potential to be bright in reflected light, because of their radii and proximity to their star, while still being cool enough to have minimal thermal contamination at optical wavelengths. For each planet we compute albedo spectra for several cloud and atmosphere assumptions (e.g. thermochemical equilibrium, solar composition). We find that in the TESS bandpass the estimated contrast at optical wavelengths is typically, 28 pages, 9 Figures, 3 Tables, Accepted to AJ

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

Explaining the evidence for surface liquid water on early Mars has been a challenge for climate modelers, as the sun was ∼ 30% less luminous during the late-Noachian. We propose that the additional greenhouse forcing of C O 2 - H 2 collision-induced absorption is capable of bringing the surface temperature above freezing and can put early Mars into a limit-cycling regime. Limit cycles occur when insolation is low and C O 2 outgassing rates are unable to balance with the rapid drawdown of C O 2 during warm weathering periods. Planets in this regime will alternate between global glaciation and transient warm climate phases. This mechanism is capable of explaining the geomorphological evidence for transient warm periods in the martian record. Previous work has shown that collision-induced absorption of C O 2 - H 2 was capable of deglaciating early Mars, but only with high H 2 outgassing rates (greater than ∼ 600 Tmol/yr) and at high surface pressures (between 3 to 4 bars). We used new theoretically derived collision-induced absorption coefficients for C O 2 - H 2 to reevaluate the climate limit cycling hypothesis for early Mars. Using the new and stronger absorption coefficients in our 1-dimensional radiative convective model as well as our energy balance model, we find that limit cycling can occur with an H 2 outgassing rate as low as ∼ 300 Tmol/yr at surface pressures below 3 bars. Our results agree more closely with paleoparameters for early martian surface pressure and hydrogen abundance.

A key legacy of the recently launched TESS mission will be to provide the astronomical community with many of the best transiting exoplanet targets for atmospheric characterization. However, time is of the essence to take full advantage of this opportunity. JWST, although delayed, will still complete its nominal five year mission on a timeline that motivates rapid identification, confirmation, and mass measurement of the top atmospheric characterization targets from TESS. Beyond JWST, future dedicated missions for atmospheric studies such as ARIEL require the discovery and confirmation of several hundred additional sub-Jovian size planets (R_p < 10 R_Earth) orbiting bright stars, beyond those known today, to ensure a successful statistical census of exoplanet atmospheres. Ground-based ELTs will also contribute to surveying the atmospheres of the transiting planets discovered by TESS. Here we present a set of two straightforward analytic metrics, quantifying the expected signal-to-noise in transmission and thermal emission spectroscopy for a given planet, that will allow the top atmospheric characterization targets to be readily identified among the TESS planet candidates. Targets that meet our proposed threshold values for these metrics would be encouraged for rapid follow-up and confirmation via radial velocity mass measurements. Based on the catalog of simulated TESS detections by Sullivan et al. (2015), we determine appropriate cutoff values of the metrics, such that the TESS mission will ultimately yield a sample of $\sim300$ high-quality atmospheric characterization targets across a range of planet size bins, extending down to Earth-size, potentially habitable worlds., accepted to PASP

A recent study by Ramirez et al. (2014) demonstrated that an atmosphere with 1.3-4 bar of CO2 and H2O, in addition to 5-20% H2, could have raised the mean annual and global surface temperature of early Mars above the freezing point of water. Such warm temperatures appear necessary to generate the rainfall (or snowfall) amounts required to carve the ancient martian valleys. Here, we use our best estimates for early martian outgassing rates, along with a 1-D photochemical model, to assess the conversion efficiency of CO, CH4, and H2S to CO2, SO2, and H2. Our outgassing estimates assume that Mars was actively recycling volatiles between its crust and interior, as Earth does today. H2 production from serpentinization and deposition of banded iron-formations is also considered. Under these assumptions, maintaining an H2 concentration of ~1-2% by volume is achievable, but reaching 5% H2 requires additional H2 sources or a slowing of the hydrogen escape rate below the diffusion limit. If the early martian atmosphere was indeed H2-rich, we might be able to see evidence of this in the rock record. The hypothesis proposed here is consistent with new data from the Curiosity Rover, which show evidence for a long-lived lake in Gale Crater near Mt. Sharp. It is also consistent with measured oxygen fugacities of martian meteorites, which show evidence for progressive mantle oxidation over time., Comment: 23 pages, 9 figures, published Icarus