Researcher profile

Kevin B. Stevenson

Kevin B. Stevenson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

preprint2023arXiv

Eureka!: An End-to-End Pipeline for JWST Time-Series Observations

$\texttt{Eureka!}$ is a data reduction and analysis pipeline for exoplanet time-series observations, with a particular focus on JWST data. Over the next 1-2 decades, JWST will pursue four main science themes: Early Universe, Galaxies Over Time, Star Lifecycle, and Other Worlds. Our focus is on providing the astronomy community with an open source tool for the reduction and analysis of time-series observations of exoplanets in pursuit of the fourth of these themes, Other Worlds. The goal of $\texttt{Eureka!}$ is to provide an end-to-end pipeline that starts with uncalibrated FITS files and ultimately yields precise exoplanet spectra. The pipeline has a modular structure with six stages, and each stage uses a "Eureka! Control File" (ECF) to allow for easy control of the pipeline's behavior. Stage 5 also uses a "Eureka! Parameter File" (EPF) to control the fitted parameters. We provide template ECFs for the MIRI, NIRCam, NIRISS, and NIRSpec instruments on JWST and the WFC3 instrument on the Hubble Space Telescope (HST). These templates give users a good starting point for their analyses, but $\texttt{Eureka!}$ is not intended to be used as a black box tool, and users should expect to fine-tune some settings for each observation in order to achieve optimal results. At each stage, the pipeline creates intermediate figures and outputs that allow users to compare $\texttt{Eureka!}$'s performance using different parameter settings or to compare $\texttt{Eureka!}$ with an independent pipeline. The ECF used to run each stage is also copied into the output folder from each stage to enhance reproducibility. Finally, while $\texttt{Eureka!}$ has been optimized for exoplanet observations (especially the latter stages of the code), much of the core functionality could also be repurposed for JWST time-series observations in other research domains thanks to $\texttt{Eureka!}$'s modularity.

preprint2022arXiv

A new method to correct for host star variability in multi-epoch observations of exoplanet transmission spectra

Transmission spectra of exoplanets orbiting active stars suffer from wavelength-dependent effects due to stellar photospheric heterogeneity. WASP-19b, an ultra-hot Jupiter (T$_{eq}$ $\sim$ 2100 K), is one such strongly irradiated gas-giant orbiting an active solar-type star. We present optical (520-900 nm) transmission spectra of WASP-19b obtained across eight epochs using the Gemini Multi-Object Spectrograph (GMOS) on the Gemini-South telescope. We apply our recently developed Gaussian Processes regression based method to model the transit light curve systematics and extract the transmission spectrum at each epoch. We find that WASP-19b's transmission spectrum is affected by stellar variability at individual epochs. We report an observed anticorrelation between the relative slopes and offsets of the spectra across all epochs. This anticorrelation is consistent with the predictions from the forward transmission models, which account for the effect of unocculted stellar spots and faculae measured previously for WASP-19. We introduce a new method to correct for this stellar variability effect at each epoch by using the observed correlation between the transmission spectral slopes and offsets. We compare our stellar variability corrected GMOS transmission spectrum with previous contradicting MOS measurements for WASP-19b and attempt to reconcile them. We also measure the amplitude and timescale of broadband stellar variability of WASP-19 from TESS photometry, which we find to be consistent with the effect observed in GMOS spectroscopy and ground-based broadband photometric long-term monitoring. Our results ultimately caution against combining multi-epoch optical transmission spectra of exoplanets orbiting active stars before correcting each epoch for stellar variability.

preprint2022arXiv

Confirmation of Water Absorption in the Thermal Emission Spectrum of the Hot Jupiter WASP-77Ab with HST/WFC3

Secondary eclipse observations of hot Jupiters can reveal both their compositions and thermal structures. Previous observations have shown a diversity of hot Jupiter eclipse spectra, including absorption features, emission features, and featureless blackbody-like spectra. We present a secondary eclipse spectrum of the hot Jupiter WASP-77Ab observed between $1-5$ $μ$m with the Hubble Space Telescope (HST) and the Spitzer Space Telescope. The HST observations show signs of water absorption indicative of a non-inverted thermal structure. We fit the data with both a one-dimensional free retrieval and a grid of one-dimensional self-consistent forward models to confirm this non-inverted structure. The free retrieval places a $3σ$ lower limit on the atmospheric water abundance of $\log(n_\mathrm{H_2O})>-4.78$ and can not constrain the CO abundance. The grid fit produces a slightly super-stellar metallicity and constrains the carbon-to-oxygen ratio to less than or equal to the solar value. We also compare our data to recent high-resolution observations of WASP-77Ab taken with the Gemini-South/IGRINS spectrograph and find that our observations are consistent with the best-fit model to the high-resolution data. However, the metallicity derived from the IGRINS data is significantly lower than that derived from our self-consistent model fit. We find that this difference may be due to disequilibrium chemistry, and the varying results between the models applied here demonstrate the difficulty of constraining disequilibrium chemistry with low-resolution, low wavelength coverage data alone. Future work to combine observations from IGRINS, HST, and JWST will improve our estimate of the atmospheric composition of WASP-77Ab.

preprint2022arXiv

Hierarchical Bayesian Atmospheric Retrieval Modeling for Population Studies of Exoplanet Atmospheres: A Case Study on the Habitable Zone

With the growing number of spectroscopic observations and observational platforms capable of exoplanet atmospheric characterization, there is a growing need for analysis techniques that can distill information about a large population of exoplanets into a coherent picture of atmospheric trends expressed within the statistical sample. In this work, we develop a Hierarchical Bayesian Atmospheric Retrieval (HBAR) model to infer population-level trends in exoplanet atmospheric characteristics. We demonstrate HBAR on the case of inferring a trend in atmospheric CO2 with incident stellar flux, predicted by the presence of a functioning carbonate-silicate weathering negative feedback cycle, an assumption upon which all calculations of the habitable zone (HZ) rest. Using simulated transmission spectra and JWST-quality observations of rocky planets with H2O, CO2, and N2 bearing atmospheres, we find that the predicted trend in CO2 causes subtle differences in the spectra of order 10 ppm in the 1-5 um range, underscoring the challenge inherent to testing this hypothesis. In the limit of highly precise data (100 stacked transits per planet), we show that our HBAR model is capable of inferring the population-level parameters that characterize the trend in CO2, and we demonstrate that the null hypothesis and other simpler trends can be rejected at high confidence. Although we find that this specific empirical test of the HZ may be prohibitively challenging in the JWST era, the HBAR framework developed in this work may find a more immediate usage for the analysis of gas giant spectra observed with JWST, Ariel, and other upcoming missions.

preprint2022arXiv

UV Absorption by Silicate Cloud Precursors in Ultra-hot Jupiter WASP-178b

Aerosols have been found to be nearly ubiquitous in substellar atmospheres. The precise temperature at which these aerosols begin to form in exoplanets has yet to be observationally constrained. Theoretical models and observations of muted spectral features suggest that silicate clouds play an important role in exoplanets between at least 950 and 2,100 K. However, some giant planets are thought to be hot enough to avoid condensation altogether. Here, we present the near-UV transmission spectrum of an ultra-hot Jupiter, WASP-178b ($\sim$2,450~K), that exhibits significant NUV absorption. This short-wavelength absorption is among the largest spectral features ever observed in an exoplanet in terms of atmospheric scale heights. Bayesian retrievals indicate the presence of gaseous refractory species containing silicon and magnesium, which are the precursors to condensate clouds at lower temperatures. SiO in particular has not been detected in exoplanets before, but the presence of SiO in WASP-178b is consistent with theoretical expectation as the dominant Si-bearing species at high temperatures. These observations allow us to re-interpret previous observations of HAT-P-41b and WASP-121b that did not consider SiO to suggest that silicate cloud formation begins on exoplanets with equilibrium temperatures between 1,950 and 2,450~K.

preprint2021arXiv

No umbrella needed: Confronting the hypothesis of iron rain on WASP-76b with post-processed general circulation models

High-resolution spectra are unique indicators of three-dimensional processes in exoplanetary atmospheres. For instance, in 2020, Ehrenreich et al. reported transmission spectra from the ESPRESSO spectrograph yielding an anomalously large Doppler blueshift from the ultra-hot Jupiter WASP-76b. Interpretations of these observations invoke toy model depictions of gas-phase iron condensation in lower-temperature regions of the planet's atmosphere. In this work, we forward model the atmosphere of WASP-76b with double-gray general circulation models (GCMs) and ray-striking radiative transfer to diagnose the planet's high-resolution transmission spectrum. We confirm that a physical mechanism driving strong east-west asymmetries across the terminator must exist to reproduce large Doppler blueshifts in WASP-76b's transmission spectrum. We identify low atmospheric drag and a deep radiative-convective boundary as necessary components of our GCM to produce this asymmetry (the latter is consistent with existing Spitzer phase curves). However, we cannot reproduce either the magnitude or the time-dependence of the WASP-76b Doppler signature with gas-phase iron condensation alone. Instead, we find that high-altitude, optically thick clouds composed of $\rm Al_2O_3$, Fe, or $\rm Mg_2SiO_4$ provide reasonable fits to the Ehrenreich et al. observations -- with marginal contributions from condensation. This fit is further improved by allowing a small orbital eccentricity ($e \approx 0.017$), consistent with prior WASP-76b orbital constraints. We additionally validate our forward-modeled spectra by reproducing lines of nearly all species detected in WASP-76b by Tabernero et al. 2021. Our procedure's success in diagnosing phase-resolved Doppler shifts demonstrates the benefits of physical, self-consistent, three-dimensional simulations in modeling high-resolution spectra of exoplanet atmospheres.

preprint2020arXiv

A New Method For Studying Exoplanet Atmospheres Using Planetary Infrared Excess

To date, the ability for observers to reveal the composition or thermal structure of an exoplanet's atmosphere has rested on two techniques: high-contrast direct imaging and time-series observations of transiting exoplanets. The former is currently limited to characterizing young, massive objects while the latter requires near 90 degree orbital inclinations, thus limiting atmospheric studies to a small fraction of the total exoplanet population. Here we present an observational and analysis technique for studying the atmospheres of non-transiting exoplanets that relies on acquiring simultaneous, broad-wavelength spectra and resolving planetary infrared emission from the stellar spectrum. This method could provide an efficient means to study exoplanet atmospheric dynamics using sparsely-sampled phase curve observations or a mechanism to search for signs of life on non-transiting exoplanets orbiting the nearest M-dwarf stars (such as Proxima Centauri). If shown to be effective with James Webb Space Telescope (JWST) observations, the method of measuring planetary infrared excess (PIE) would open up the large population of nearby, non-transiting exoplanets for atmospheric characterization.

preprint2020arXiv

Clouds in Three-Dimensional Models of Hot Jupiters Over a Wide Range of Temperatures I: Thermal Structures and Broadband Phase Curve Predictions

Using a general circulation model (GCM), we investigate trends in simulated hot Jupiter atmospheres for a range of irradiation temperatures (1,500 - 4,000 K), surface gravities (10 and 40 m s-2), and cloud conditions. Our models include simplified temperature-dependent clouds with radiative feedback and show how different cloud compositions, vertical thicknesses, and opacities shape hot Jupiters atmospheres by potentially increasing planetary albedos, decreasing photospheric pressures and nightside temperatures, and in some cases producing strong dayside thermal inversions. With decreasing irradiation, clouds progressively form on the nightside and cooler western limb followed by the eastern limb and central dayside. We find that clouds significantly modify the radiative transport and affect the observable properties of planets colder than T_irr ~ 3,000~K (T_eq~2,100 K) depending on the clouds&#39; vertical extent. The precise strength of expected effects depends on the assumed parameters, but trends in predicted phase curves emerge from an ensemble of simulations. Clouds lead to larger phase curve amplitudes and smaller phase curve offsets at IR wavelengths, compared to cloud-free models. At optical wavelengths, we predict mostly westward phase curve offsets at intermediate temperatures (T_irr ~ 2,000 - 3,500 K) with clouds confined to the nightside and western limb. If clouds are vertically compact (i.e. on order of a pressure scale height in thickness), their distributions and effects become more complicated as different condensates form at different heights -- some too deep to significantly affect the observable atmosphere. Our results have implications for interpreting the diversity of phase curve observations of planets with T_irr <~3,000~K.

preprint2020arXiv

Evidence for H$_{2}$ Dissociation and Recombination Heat Transport in the Atmosphere of KELT-9b

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $μ$m. KELT-9b is the hottest known transiting planet, with a 4.5-$μ$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5σ$ confidence. This discrepancy may be due to magnetic effects in the planet&#39;s highly ionized atmosphere.

preprint2020arXiv

Global Chemistry and Thermal Structure Models for the Hot Jupiter WASP-43b and Predictions for JWST

The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5--12 $μ$m range with JWST/MIRI. To prepare for these observations, we performed an extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 $μ$m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the night-side emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, providing eventually the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition.

preprint2020arXiv

Introducing a New Spitzer Master BLISS Map to Remove the Instrument-Systematic -- Phase-Curve-Parameter Degeneracy, as Demonstrated by a Reanalysis of the 4.5$μm$ WASP-43b Phase Curve

While Spitzer IRAC systematics are generally well understood, each data set can provide its own challenges that continue to teach us about the underlying functional form of these systematics. Multiple groups have analyzed the phase curves of WASP-43b with varying detrending techniques, each obtaining different results. In this work, we take another look at WASP-43b while exploring the degenerate relation between Bilinearly-Interpolated Subpixel Sensitivity (BLISS) mapping, point response function full width at half-maximum (PRF-FWHM) detrending, and phase curve parameters. We find that there is a strong correlation between the detrending parameters in the two models, and best fit phase curve amplitudes vary strongly when the data are temporally binned. To remove this degeneracy, we present a new Gaussian centroided intrapixel sensitivity map (hereafter fixed sensitivity map), generated using 3,712,830 exposures spanning 5 years, for a variety of aperture sizes at 4.5$μm$. We find evidence for time variability in the sensitivity at 3.6$μm$ and do not generate a visit-independent map for this channel. With the fixed 4.5 $μm$ intrapixel sensitivity map, the best fits for WASP-43b no longer vary strongly with bin size and PRF-FWHM detrending is no longer required to remove correlated noise. For data sets that do not fall completely within the sweet spot, temporal binning should not be used in the analysis of Spitzer phase curves. We confirm night side emission for WASP-43b with a disk integrated nightside temperature of 806 $\pm$ 48 K at 4.5$μm$. The 4.5$μm$ maps are available at \url{github.com/kevin218/POET}

preprint2020arXiv

Smaller than expected bright-spot offsets in Spitzer phase curves of the hot Jupiter Qatar-1b

We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~μ$m and $0.30 \pm 0.02 \%$ at $4.5~μ$m, corresponding to dayside brightness temperatures of $1542^{+32}_{-31}$~K and $1557^{+35}_{-36}$~K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are $1117^{+76}_{-71}$~K and $1167^{+69}_{-74}$~K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of $0.12_{-0.16}^{+0.14}$ and a moderate day-night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright-spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with \textit{Spitzer} full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: $11^{\circ}\pm 7^{\circ}$ at $3.6~μ$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~μ$m. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters.

preprint2020arXiv

Statistical Characterization of Hot Jupiter Atmospheres using Spitzer&#39;s Secondary Eclipses

We report 78 secondary eclipse depths for a sample of 36 transiting hot Jupiters observed at 3.6- and 4.5 microns using the Spitzer Space Telescope. Our eclipse results for 27 of these planets are new, and include highly irradiated worlds such as KELT-7b, WASP-87b, WASP-76b, and WASP-64b, and important targets for JWST such as WASP-62b. We find that WASP-62b has a slightly eccentric orbit e cos(omega) = 0.00614+/- 0.00064, and we confirm the eccentricity of HAT-P-13b and WASP-14b. The remainder are individually consistent with circular orbits, but we find statistical evidence for eccentricity increasing with orbital period in our range from 1 to 5 days. Our day-side brightness temperatures for the planets yield information on albedo and heat redistribution, following Cowan and Agol (2011). Planets having maximum day side temperatures exceeding ~ 2200K are consistent with zero albedo and distribution of stellar irradiance uniformly over the day-side hemisphere. Our most intriguing result is that we detect a systematic difference between the emergent spectra of these hot Jupiters as compared to blackbodies. The ratio of observed brightness temperatures, Tb(4.5)/Tb(3.6), increases with equilibrium temperature by 100 +/- 24 parts-per-million per Kelvin, over the entire temperature range in our sample (800K to 2500K). No existing model predicts this trend over such a large range of temperature. We suggest that this may be due to a structural difference in the atmospheric temperature profile between the real planetary atmospheres as compared to models.

preprint2019arXiv

Impact of Clouds and Hazes on the Simulated JWST Transmission Spectra of Habitable Zone Planets in the TRAPPIST-1 System

The TRAPPIST-1 system, consisting of an ultra-cool host star having seven known Earth-size planets will be a prime target for atmospheric characterization with JWST. However, the detectability of atmospheric molecular species may be severely impacted by the presence of clouds and/or hazes. In this work, we perform 3-D General Circulation Model (GCM) simulations with the LMD Generic model supplemented by 1-D photochemistry simulations at the terminator with the Atmos model to simulate several possible atmospheres for TRAPPIST-1e, 1f and 1g: 1) modern Earth, 2) Archean Earth, and 3) CO2-rich atmospheres. JWST synthetic transit spectra were computed using the GSFC Planetary Spectrum Generator (PSG). We find that TRAPPIST-1e, 1f and 1g atmospheres, with clouds and/or hazes, could be detected using JWST&#39;s NIRSpec prism from the CO2 absorption line at 4.3 um in less than 15 transits at 3 sigma or less than 35 transits at 5 sigma. However, our analysis suggests that other gases would require hundreds (or thousands) of transits to be detectable. We also find that H2O, mostly confined in the lower atmosphere, is very challenging to detect for these planets or similar systems if the planets&#39; atmospheres are not in a moist greenhouse state. This result demonstrates that the use of GCMs, self-consistently taking into account the effect of clouds and sub-saturation, is crucial to evaluate the detectability of atmospheric molecules of interest as well as for interpreting future detections in a more global (and thus robust and relevant) approach.

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&#39;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.