Normalized to: Blecic, J.
[1]
oai:arXiv.org:2001.04759 [pdf] - 2057869
Global Chemistry and Thermal Structure Models for the Hot Jupiter
WASP-43b and Predictions for JWST
Venot, Olivia;
Parmentier, Vivien;
Blecic, Jasmina;
Cubillos, Patricio E.;
Waldmann, Ingo P.;
Changeat, Quentin;
Moses, Julianne I.;
Tremblin, Pascal;
Crouzet, Nicolas;
Gao, Peter;
Powell, Diana;
Lagage, Pierre-Olivier;
Dobbs-Dixon, Ian;
Steinrueck, Maria E.;
Kreidberg, Laura;
Batalha, Natalie;
Bean, Jacob L.;
Stevenson, Kevin B.;
Casewell, Sarah;
Carone, Ludmila
Submitted: 2020-01-14
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 $\mu$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 $\mu$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.
[2]
oai:arXiv.org:1905.07064 [pdf] - 1884966
The Need for Laboratory Measurements and Ab Initio Studies to Aid
Understanding of Exoplanetary Atmospheres
Fortney, Jonathan J.;
Robinson, Tyler D.;
Domagal-Goldman, Shawn;
Del Genio, Anthony D.;
Gordon, Iouli E.;
Gharib-Nezhad, Ehsan;
Lewis, Nikole;
Sousa-Silva, Clara;
Airapetian, Vladimir;
Drouin, Brian;
Hargreaves, Robert J.;
Huang, Xinchuan;
Karman, Tijs;
Ramirez, Ramses M.;
Rieker, Gregory B.;
Tennyson, Jonathan;
Wordsworth, Robin;
Yurchenko, Sergei N;
Johnson, Alexandria V;
Lee, Timothy J.;
Dong, Chuanfei;
Kane, Stephen;
Lopez-Morales, Mercedes;
Fauchez, Thomas;
Lee, Timothy;
Marley, Mark S.;
Sung, Keeyoon;
Haghighipour, Nader;
Robinson, Tyler;
Horst, Sarah;
Gao, Peter;
Kao, Der-you;
Dressing, Courtney;
Lupu, Roxana;
Savin, Daniel Wolf;
Fleury, Benjamin;
Venot, Olivia;
Ascenzi, Daniela;
Milam, Stefanie;
Linnartz, Harold;
Gudipati, Murthy;
Gronoff, Guillaume;
Salama, Farid;
Gavilan, Lisseth;
Bouwman, Jordy;
Turbet, Martin;
Benilan, Yves;
Henderson, Bryana;
Batalha, Natalie;
Jensen-Clem, Rebecca;
Lyons, Timothy;
Freedman, Richard;
Schwieterman, Edward;
Goyal, Jayesh;
Mancini, Luigi;
Irwin, Patrick;
Desert, Jean-Michel;
Molaverdikhani, Karan;
Gizis, John;
Taylor, Jake;
Lothringer, Joshua;
Pierrehumbert, Raymond;
Zellem, Robert;
Batalha, Natasha;
Rugheimer, Sarah;
Lustig-Yaeger, Jacob;
Hu, Renyu;
Kempton, Eliza;
Arney, Giada;
Line, Mike;
Alam, Munazza;
Moses, Julianne;
Iro, Nicolas;
Kreidberg, Laura;
Blecic, Jasmina;
Louden, Tom;
Molliere, Paul;
Stevenson, Kevin;
Swain, Mark;
Bott, Kimberly;
Madhusudhan, Nikku;
Krissansen-Totton, Joshua;
Deming, Drake;
Kitiashvili, Irina;
Shkolnik, Evgenya;
Rustamkulov, Zafar;
Rogers, Leslie;
Close, Laird
Submitted: 2019-05-16
We are now on a clear trajectory for improvements in exoplanet observations
that will revolutionize our ability to characterize their atmospheric
structure, composition, and circulation, from gas giants to rocky planets.
However, exoplanet atmospheric models capable of interpreting the upcoming
observations are often limited by insufficiencies in the laboratory and
theoretical data that serve as critical inputs to atmospheric physical and
chemical tools. Here we provide an up-to-date and condensed description of
areas where laboratory and/or ab initio investigations could fill critical gaps
in our ability to model exoplanet atmospheric opacities, clouds, and chemistry,
building off a larger 2016 white paper, and endorsed by the NAS Exoplanet
Science Strategy report. Now is the ideal time for progress in these areas, but
this progress requires better access to, understanding of, and training in the
production of spectroscopic data as well as a better insight into chemical
reaction kinetics both thermal and radiation-induced at a broad range of
temperatures. Given that most published efforts have emphasized relatively
Earth-like conditions, we can expect significant and enlightening discoveries
as emphasis moves to the exotic atmospheres of exoplanets.
[3]
oai:arXiv.org:1903.09322 [pdf] - 1857235
Imaging Cool Giant Planets in Reflected Light: Science Investigations
and Synergy with Habitable Planets
Marley, Mark;
Lewis, Nikole;
Arney, Giada;
Bailey, Vanessa;
Batalha, Natasha;
Beichman, Charles;
Benneke, Björn;
Blecic, Jasmina;
Cahoy, Kerri;
Chilcote, Jeffrey;
Domagal-Goldman, Shawn;
Dressing, Courtney;
Fitzgerald, Michael;
Fortney, Jonathan;
Freedman, Richard;
Gelino, Dawn;
Gizis, John;
Guyon, Olivier;
Greene, Thomas;
Hammel, Heidi;
Hasegawa, Yasuhiro;
Jovanovic, Nemanja;
Konopacky, Quinn;
Kopparapu, Ravi;
Liu, Michael;
Lopez, Eric;
Lunine, Jonathan;
Lupu, Roxana;
Macintosh, Bruce;
Mandt, Kathleen;
Marois, Christian;
Mawet, Dimitri;
Mayorga, Laura;
Morley, Caroline;
Nielsen, Eric;
Roberge, Aki;
Serabyn, Eugene;
Skemer, Andrew;
Stapelfeldt, Karl;
Visscher, Channon;
Wang, Jason
Submitted: 2019-03-21, last modified: 2019-03-26
Planned astronomical observatories of the 2020s will be capable of obtaining
reflected light photometry and spectroscopy of cool extrasolar giant planets.
Here we explain that such data are valuable both for understanding the origin
and evolution of giant planets as a whole and for preparing for the
interpretation of similar datasets from potentially habitable extrasolar
terrestrial planets in the decades to follow.
[4]
oai:arXiv.org:1901.03764 [pdf] - 1838327
Toward More Reliable Analytic Thermochemical-equilibrium Abundances
Submitted: 2019-01-11, last modified: 2019-01-31
Heng & Tsai (2016) developed an analytic framework to obtain
thermochemical-equilibrium abundances for H$_{2}$O, CO, CO$_2$, CH$_4$,
C$_{2}$H$_{2}$, C$_{2}$H$_{4}$, HCN, NH$_3$, and N$_2$ for a system with known
temperature, pressure, and elemental abundances (hydrogen, carbon, nitrogen,
and oxygen). However, the implementation of their approach can become
numerically unstable under certain circumstances, leading to inaccurate
solutions (e.g., ${\rm C/O} \ge 1$ atmospheres at low pressures). Building up
on their approach, we identified the conditions that prompt inaccurate
solutions, and developed a new framework to avoid them, providing a reliable
implementation for arbitrary values of temperature (200 to $\sim$2000 K),
pressure ($10^{-8}$ to $10^{3}$ bar), and CNO abundances ($10^{-3}$ to $\sim
10^{2}\times$ solar elemental abundances), for hydrogen-dominated atmospheres.
The accuracy our analytic framework is better than 10% for the more abundant
species that have mixing fractions larger than $10^{-10}$, whereas the accuracy
is better than 50% for the less abundant species. Additionally, we added the
equilibrium-abundance calculation of atomic and molecular hydrogen into the
system, and explored the physical limitations of this approach. Efficient and
reliable tools, such as this one, are highly valuable for atmospheric Bayesian
studies, which need to evaluate a large number of models. We implemented our
analytic framework into the \textsc{rate} Python open-source package, available
at https://github.com/pcubillos/rate .
[5]
oai:arXiv.org:1803.04985 [pdf] - 1767381
The Transiting Exoplanet Community Early Release Science Program for
JWST
Bean, Jacob L.;
Stevenson, Kevin B.;
Batalha, Natalie M.;
Berta-Thompson, Zachory;
Kreidberg, Laura;
Crouzet, Nicolas;
Benneke, Björn;
Line, Michael R.;
Sing, David K.;
Wakeford, Hannah R.;
Knutson, Heather A.;
Kempton, Eliza M. -R.;
Désert, Jean-Michel;
Crossfield, Ian;
Batalha, Natasha E.;
de Wit, Julien;
Parmentier, Vivien;
Harrington, Joseph;
Moses, Julianne I.;
Lopez-Morales, Mercedes;
Alam, Munazza K.;
Blecic, Jasmina;
Bruno, Giovanni;
Carter, Aarynn L.;
Chapman, John W.;
Decin, Leen;
Dragomir, Diana;
Evans, Thomas M.;
Fortney, Jonathan J.;
Fraine, Jonathan D.;
Gao, Peter;
Muñoz, Antonio García;
Gibson, Neale P.;
Goyal, Jayesh M.;
Heng, Kevin;
Hu, Renyu;
Kendrew, Sarah;
Kilpatrick, Brian M.;
Krick, Jessica;
Lagage, Pierre-Olivier;
Lendl, Monika;
Louden, Tom;
Madhusudhan, Nikku;
Mandell, Avi M.;
Mansfield, Megan;
May, Erin M.;
Morello, Giuseppe;
Morley, Caroline V.;
Nikolov, Nikolay;
Redfield, Seth;
Roberts, Jessica E.;
Schlawin, Everett;
Spake, Jessica J.;
Todorov, Kamen O.;
Tsiaras, Angelos;
Venot, Olivia;
Waalkes, William C.;
Wheatley, Peter J.;
Zellem, Robert T.;
Angerhausen, Daniel;
Barrado, David;
Carone, Ludmila;
Casewell, Sarah L.;
Cubillos, Patricio E.;
Damiano, Mario;
de Val-Borro, Miguel;
Drummond, Benjamin;
Edwards, Billy;
Endl, Michael;
Espinoza, Nestor;
France, Kevin;
Gizis, John E.;
Greene, Thomas P.;
Henning, Thomas K.;
Hong, Yucian;
Ingalls, James G.;
Iro, Nicolas;
Irwin, Patrick G. J.;
Kataria, Tiffany;
Lahuis, Fred;
Leconte, Jérémy;
Lillo-Box, Jorge;
Lines, Stefan;
Lothringer, Joshua D.;
Mancini, Luigi;
Marchis, Franck;
Mayne, Nathan;
Palle, Enric;
Rauscher, Emily;
Roudier, Gaël;
Shkolnik, Evgenya L.;
Southworth, John;
Swain, Mark R.;
Taylor, Jake;
Teske, Johanna;
Tinetti, Giovanna;
Tremblin, Pascal;
Tucker, Gregory S.;
van Boekel, Roy;
Waldmann, Ingo P.;
Weaver, Ian C.;
Zingales, Tiziano
Submitted: 2018-03-13, last modified: 2018-09-03
The James Webb Space Telescope (JWST) presents the opportunity to transform
our understanding of planets and the origins of life by revealing the
atmospheric compositions, structures, and dynamics of transiting exoplanets in
unprecedented detail. However, the high-precision, time-series observations
required for such investigations have unique technical challenges, and prior
experience with other facilities indicates that there will be a steep learning
curve when JWST becomes operational. In this paper we describe the science
objectives and detailed plans of the Transiting Exoplanet Community Early
Release Science (ERS) Program, which is a recently approved program for JWST
observations early in Cycle 1. The goal of this project, for which the obtained
data will have no exclusive access period, is to accelerate the acquisition and
diffusion of technical expertise for transiting exoplanet observations with
JWST, while also providing a compelling set of representative datasets that
will enable immediate scientific breakthroughs. The Transiting Exoplanet
Community ERS Program will exercise the time-series modes of all four JWST
instruments that have been identified as the consensus highest priorities,
observe the full suite of transiting planet characterization geometries
(transits, eclipses, and phase curves), and target planets with host stars that
span an illustrative range of brightnesses. The observations in this program
were defined through an inclusive and transparent process that had
participation from JWST instrument experts and international leaders in
transiting exoplanet studies. Community engagement in the project will be
centered on a two-phase Data Challenge that culminates with the delivery of
planetary spectra, time-series instrument performance reports, and open-source
data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST
mission.
[6]
oai:arXiv.org:1704.07421 [pdf] - 1741983
Community targets for JWST's early release science program: evaluation
of WASP-63b
Kilpatrick, Brian M.;
Cubillos, Patricio E.;
Stevenson, Kevin B.;
Lewis, Nikole K.;
Wakeford, Hannah;
Macdonald, Ryan J.;
Madhusudhan, Nikku;
Blecic, Jasmina;
Bruno, Giovanni;
Burrows, Adam;
Deming, Drake;
Heng, Kevin;
Line, Michael R.;
Morley, Caroline V.;
Parmentier, Vivien;
Tucker, Gregory S.;
Valenti, Jeff A.;
Waldmann, Ingo P.;
Bean, Jacob L.;
Beichman, Charles;
Fraine, Jonathan;
Krick, J. E.;
Lothringer, Joshua D.;
Mandell, Avi M.
Submitted: 2017-04-24, last modified: 2018-08-30
We present observations of WASP-63b by the Hubble Space Telescope (HST) as
part of "A Preparatory Program to Identify the Single Best Transiting Exoplanet
for JWST Early Release Science". WASP-63b is one of the community targets under
consideration for the James Webb Space Telescope (JWST) Early Release Science
(ERS) program. We present a spectrum derived from a single observation by HST
Wide Field Camera 3 in the near infrared. We engaged groups across the
transiting exoplanet community to participate in the analysis of the data and
present results from each. There is general agreement amongst all results that
we find an H2O absorption feature with 3.5-4.0 sigma significance. However, the
feature is muted in comparison to a clear atmosphere at solar composition.
Although the detection of the water feature is robust, the reasons for the
muting of this feature are ambiguous due to a degeneracy between clouds and
composition. The data does not yield robust detections of any molecular species
other than H2O. The group was motivated to perform an additional set of
retrieval exercises to investigate an apparent bump in the spectrum at ~ 1.55
um. We explore possible disequilibrium chemistry and find this feature is
consistent with super-solar HCN abundance but it is questionable if the
required mixing ratio of HCN is chemically and physically plausible. The
ultimate goal of this study is to vet WASP-63b as a potential community target
to best demonstrate the capabilities and systematics of JWST instruments for
transiting exoplanet science. In the case of WASP-63b, the presence of a
detectable water feature indicates that WASP-63b remains a plausible target for
ERS observations.
[7]
oai:arXiv.org:1803.06678 [pdf] - 1651725
The Implications of 3D Thermal Structure on 1D Atmospheric Retrieval
Submitted: 2018-03-18
Using the atmospheric structure from a 3D global radiation-hydrodynamic
simulation of HD 189733b and the open-source BART code, we investigate the
difference between the secondary-eclipse temperature structure produced with a
3D simulation and the best-fit 1D retrieved model. Synthetic data are generated
by integrating the 3D models over the Spitzer, HST, and JWST bandpasses,
covering the wavelength range between 1 and 11 um. Using the data from
different observing instruments, we present detailed comparisons between the
temperature-pressure profiles recovered by BART and those from the 3D
simulations. We calculate several averages of the 3D thermal structure and
implement two temperature parameterizations to investigate different thermal
profile shapes. To assess which part of the thermal structure is best
constrained by the data, we generate contribution functions for both our
theoretical model and each of our retrieved models. Our conclusions are
strongly affected by the spectral resolution of the instruments included, their
wavelength coverage, and the number of data points combined. We also see some
limitations in each of the temperature parametrizations. The results show that
our 1D retrieval is recovering a temperature and pressure profile that most
closely matches the arithmetic average of the 3D thermal structure. When we use
a higher resolution, more data points, and a parametrized temperature profile
that allows more flexibility in the middle part of the atmosphere, we find a
better match between the retrieved temperature and pressure profile and the
arithmetic average.
[8]
oai:arXiv.org:1801.08482 [pdf] - 1709326
Dust in brown dwarfs and extra-solar planets VI. Assessing seed
formation across the brown dwarf and exoplanet regimes
Submitted: 2018-01-25
The cloud formation process starts with the formation of seed particles,
after which, surface chemical reactions grow or erode the cloud particles. We
investigate which materials may form cloud condensation seeds in the gas
temperature and pressure regimes (T$_{\rm gas}$ = 100-2000 K, p$_{\rm gas}$ =
10$^{-8}$-100 bar) expected to occur in planetary and brown dwarf atmospheres.
We apply modified classical nucleation theory which requires surface tensions
and vapour pressure data for each solid species, which are taken from the
literature. We calculate the seed formation rates of TiO$_{2}$[s] and SiO[s]
and find that they efficiently nucleate at high temperatures of T$_{\rm gas}$ =
1000-1750 K. Cr[s], KCl[s] and NaCl[s] are found to efficiently nucleate across
an intermediate temperature range of T$_{\rm gas}$ = 500-1000 K. We find
CsCl[s] may serve as the seed particle for the water cloud layers in cool
sub-stellar atmospheres. Four low temperature ice species, H$_{2}$O[s/l],
NH$_{3}$[s], H$_{2}$S[s/l] and CH$_{4}$[s], nucleation rates (T$_{\rm gas}$ =
100-250 K) are also investigated for the coolest sub-stellar/planetary
atmospheres. Our results suggest a possibly, (T$_{\rm gas}$, p$_{\rm gas}$)
distributed hierarchy of seed particle formation regimes throughout the
sub-stellar and planetary atmospheric temperature-pressure space. In order to
improve the accuracy of the nucleation rate calculation, further research into
the small cluster thermochemical data for each cloud species is warranted. The
validity of these seed particle scenarios will be tested by applying it to more
complete cloud models in the future.
[9]
oai:arXiv.org:1712.01010 [pdf] - 1693956
Equilibrium chemistry down to 100 K - Impact of silicates and
phyllosilicates on carbon/oxygen ratio
Submitted: 2017-12-04
We introduce a fast and versatile computer code, GGchem, to determine the
chemical composition of gases in thermo-chemical equilibrium down to 100 K,
with or without equilibrium condensation. We review the data for molecular
equilibrium constants, kp(T), from several sources and discuss which functional
fits are most suitable for low temperatures. We benchmark our results against
another chemical equilibrium code. We collect Gibbs free energies, dG(T), for
about 200 solid and liquid species from the NIST-JANAF database and the
geophysical database SUPCRTBL. We discuss the condensation sequence of the
elements with solar abundances in phase equilibrium down to 100 K. Once the
major magnesium silicates Mg2SiO4[s] and MgSiO3[s] have formed, the dust/gas
mass ratio jumps to a value of about 0.0045 which is significantly lower than
the often assumed value of 0.01. Silicate condensation is found to increase the
carbon/oxygen ratio (C/O) in the gas from its solar value of ~0.55 up to ~0.71,
and, by the additional intake of water and hydroxyl into the solid matrix, the
formation of phyllosilicates at temperatures below ~400 K increases the gaseous
C/O further to about 0.83. Metallic tungsten (W) is the first condensate found
to become thermodynamically stable around 1600 - 2200 K (depending on
pressure), several hundreds of Kelvin before subsequent materials like
zirconium dioxide (ZrO2) or corundum (Al2O3) can condense. We briefly discuss
whether tungsten, despite its low abundance of ~2.E-7 times the silicon
abundance, could provide the first seed particles for astrophysical dust
formation. The GGchem code is publicly available at
https://github.com/pw31/GGchem.
[10]
oai:arXiv.org:1701.00828 [pdf] - 1534094
The Atmosphere and Interior Structure of HAT-P-13b from Spitzer
Secondary Eclipses
Submitted: 2017-01-03
We present {\em Spitzer} secondary-eclipse observations of the hot Jupiter
HAT-P-13 b in the 3.6 {\micron} and 4.5 {\micron} bands. HAT-P-13 b inhabits a
two-planet system with a configuration that enables constraints on the planet's
second Love number, \math{k\sb{2}}, from precise eccentricity measurements,
which in turn constrains models of the planet's interior structure. We exploit
the direct measurements of \math{e \cos \omega} from our secondary-eclipse data
and combine them with previously published radial velocity data to generate a
refined model of the planet's orbit and thus an improved estimate on the
possible interval for \math{k\sb{2}}. We report eclipse phases of \math{0.49154
\pm 0.00080} and \math{0.49711 \pm 0.00083} and corresponding \math{e \cos
\omega} estimates of \math{-0.0136 \pm 0.0013} and \math{-0.0048 \pm 0.0013}.
Under the assumptions of previous work, our estimate of \math{k\sb{2}} of 0.81
{\pm} 0.10 is consistent with the lower extremes of possible core masses found
by previous models, including models with no solid core. This anomalous result
challenges both interior models and the dynamical assumptions that enable them,
including the essential assumption of apsidal alignment. We also report eclipse
depths of 0.081\% {\pm} 0.008\% in the 3.6 {\micron} channel and 0.088 \% {\pm}
0.028 \% in the 4.5 {\micron} channel. These photometric results are
non-uniquely consistent with solar-abundance composition without any thermal
inversion.
[11]
oai:arXiv.org:1610.01336 [pdf] - 1531883
On Correlated-noise Analyses Applied To Exoplanet Light Curves
Submitted: 2016-10-05
Time-correlated noise is a significant source of uncertainty when modeling
exoplanet light-curve data. A correct assessment of correlated noise is
fundamental to determine the true statistical significance of our findings.
Here we review three of the most widely used correlated-noise estimators in the
exoplanet field, the time-averaging, residual-permutation, and
wavelet-likelihood methods. We argue that the residual-permutation method is
unsound in estimating the uncertainty of parameter estimates. We thus recommend
to refrain from this method altogether. We characterize the behavior of the
time averaging's rms-vs.-bin-size curves at bin sizes similar to the total
observation duration, which may lead to underestimated uncertainties. For the
wavelet-likelihood method, we note errors in the published equations and
provide a list of corrections. We further assess the performance of these
techniques by injecting and retrieving eclipse signals into synthetic and real
Spitzer light curves, analyzing the results in terms of the relative-accuracy
and coverage-fraction statistics. Both the time-averaging and
wavelet-likelihood methods significantly improve the estimate of the eclipse
depth over a white-noise analysis (a Markov-chain Monte Carlo exploration
assuming uncorrelated noise). However, the corrections are not perfect, when
retrieving the eclipse depth from Spitzer datasets, these methods covered the
true (injected) depth within the 68\% credible region in only $\sim$45--65\% of
the trials. Lastly, we present our open-source model-fitting tool, Multi-Core
Markov-Chain Monte Carlo ({MC$^3$}). This package uses Bayesian statistics to
estimate the best-fitting values and the credible regions for the parameters
for a (user-provided) model. {MC$^3$} is a Python/C code, available at
https://github.com/pcubillos/MCcubed.
[12]
oai:arXiv.org:1604.02692 [pdf] - 1388317
Observations, Thermochemical Calculations, and Modeling of Exoplanetary
Atmospheres
Submitted: 2016-04-10
This dissertation as a whole aims to provide means to better understand
hot-Jupiter planets through observing, performing thermochemical calculations,
and modeling their atmospheres. We used Spitzer multi-wavelength
secondary-eclipse observations and targets with high signal-to-noise ratios, as
their deep eclipses allow us to detect signatures of spectral features and
assess planetary atmospheric structure and composition with greater certainty.
Chapter 1 gives a short introduction. Chapter 2 presents the Spitzer
secondary-eclipse analysis and atmospheric characterization of WASP-14b.
WASP-14b is a highly irradiated, transiting hot Jupiter. By applying a Bayesian
approach in the atmospheric analysis, we found an absence of thermal inversion
contrary to theoretical predictions. Chapter 3 describes the infrared
observations of WASP-43b Spitzer secondary eclipses, data analysis, and
atmospheric characterization. WASP-43b is one of the closest-orbiting hot
Jupiters, orbiting one of the coolest stars with a hot Jupiter. The atmospheric
analysis ruled out a strong thermal inversion in its dayside atmosphere.
Chapter 4 presents an open-source Thermochemical Equilibrium Abundances (TEA)
code and its application to several hot Jupiters. TEA calculates the abundances
of gaseous species using the Gibbs free-energy minimization method within an
iterative Lagrangian optimization scheme. The code is written in Python and
available to the community via http://github.com/dzesmin/TEA. Chapter 5
presents my contributions to an open-source Bayesian Atmospheric Radiative
Transfer (BART) code, and its application to WASP-43b. BART characterizes
planetary atmospheres based on the observed spectroscopic information. It
initializes a planetary atmospheric model, performs radiative-transfer
calculations to produce models of planetary spectra, and using a statistical
module compares models with observations.
[13]
oai:arXiv.org:1505.06392 [pdf] - 1439690
TEA: A Code for Calculating Thermochemical Equilibrium Abundances
Submitted: 2015-05-23
We present an open-source Thermochemical Equilibrium Abundances (TEA) code
that calculates the abundances of gaseous molecular species. The code is based
on the methodology of White et al. (1958) and Eriksson (1971). It applies Gibbs
free-energy minimization using an iterative, Lagrangian optimization scheme.
Given elemental abundances, TEA calculates molecular abundances for a
particular temperature and pressure or a list of temperature-pressure pairs. We
tested the code against the method of Burrows & Sharp (1999), the free
thermochemical equilibrium code CEA (Chemical Equilibrium with Applications),
and the example given by White et al. (1958). Using their thermodynamic data,
TEA reproduces their final abundances, but with higher precision. We also
applied the TEA abundance calculations to models of several hot-Jupiter
exoplanets, producing expected results. TEA is written in Python in a modular
format. There is a start guide, a user manual, and a code document in addition
to this theory paper. TEA is available under a reproducible-research,
open-source license via https://github.com/dzesmin/TEA.
[14]
oai:arXiv.org:1302.7003 [pdf] - 785169
Spitzer observations of the thermal emission from WASP-43b
Blecic, Jasmina;
Harrington, Joseph;
Madhusudhan, Nikku;
Stevenson, Kevin B.;
Hardy, Ryan A.;
Cubillos, Patricio E.;
Hardin, Matthew;
Nymeyer, Sarah;
Anderson, David R.;
Hellier, Coel;
Smith, Alexis M. S.;
Cameron, Andrew Collier
Submitted: 2013-02-27, last modified: 2014-02-18
WASP-43b is one of the closest-orbiting hot Jupiters, with a semimajor axis
of a = 0.01526 +/- 0.00018 AU and a period of only 0.81 days. However, it
orbits one of the coolest stars with a hot Jupiter (Tstar = 4520 +/- 120 K),
giving the planet a modest equilibrium temperature of Teq = 1440 +/- 40 K,
assuming zero Bond albedo and uniform planetary energy redistribution. The
eclipse depths and brightness temperatures from our jointly fit model are
0.347% +/- 0.013% and 1670 +/- 23 K at 3.6 {\mu}m and 0.382% +/- 0.015% and
1514 +/- 25 K at 4.5 {\mu}m. The eclipse timings improved the estimate of the
orbital period, P, by a factor of three (P = 0.81347436 +/- 1.4*10-7 days) and
put an upper limit on the eccentricity (e = 0.010+0.010 -0.007). We use our
Spitzer eclipse depths along with four previously reported ground-based
photometric observations in the near-infrared to constrain the atmospheric
properties of WASP-43b. The data rule out a strong thermal inversion in the
dayside atmosphere of WASP-43b. Model atmospheres with no thermal inversions
and fiducial oxygen-rich compositions are able to explain all the available
data. However, a wide range of metallicities and C/O ratios can explain the
data. The data suggest low day-night energy redistribution in the planet,
consistent with previous studies, with a nominal upper limit of about 35% for
the fraction of energy incident on the dayside that is redistributed to the
nightside.
[15]
oai:arXiv.org:1111.2363 [pdf] - 748009
Thermal Emission of WASP-14b Revealed with Three Spitzer Eclipses
Blecic, Jasmina;
Harrington, Joseph;
Madhusudhan, Nikku;
Stevenson, Kevin B.;
Hardy, Ryan A.;
Cubillos, Patricio;
Hardin, Matthew;
Campo, Christopher J.;
Bowman, William C.;
Nymeyer, Sarah;
Loredo, Thomas J.;
Anderson, David R.;
Maxted, Pierre F. L.
Submitted: 2011-11-09, last modified: 2013-11-19
Exoplanet WASP-14b is a highly irradiated, transiting hot Jupiter. Joshi et
al. calculate an equilibrium temperature Teq of 1866 K for zero albedo and
reemission from the entire planet, a mass of 7.3 +/- 0.5 Jupiter masses and a
radius of 1.28 +/- 0.08 Jupiter radii. Its mean density of 4.6 g/cm3 is one of
the highest known for planets with periods less than 3 days. We obtained three
secondary eclipse light curves with the Spitzer Space Telescope. The eclipse
depths from the best jointly fit model are $0.224\%$ +/- $0.018\%$ at 4.5
{\mu}m and $0.181\%$ +/- $0.022\%$ at 8.0 {\mu}m. The corresponding brightness
temperatures are 2212 +/- 94 K and 1590 +/- 116 K. A slight ambiguity between
systematic models suggests a conservative 3.6 {\mu}m eclipse depth of $0.19\%$
+/- $0.01\%$ and brightness temperature of 2242 +/- 55 K. Although extremely
irradiated, WASP-14b does not show any distinct evidence of a thermal
inversion. In addition, the present data nominally favor models with day night
energy redistribution less than $~30\%$. The current data are generally
consistent with oxygen-rich as well as carbon-rich compositions, although an
oxygen-rich composition provides a marginally better fit. We confirm a
significant eccentricity of e = 0.087 +/- 0.002 and refine other orbital
parameters.
[16]
oai:arXiv.org:1303.5468 [pdf] - 1165462
WASP-8b: Characterization of a Cool and Eccentric Exoplanet with Spitzer
Submitted: 2013-03-21
WASP-8b has 2.18 times Jupiter's mass and is on an eccentric ($e=0.31$)
8.16-day orbit. With a time-averaged equilibrium temperature of 948 K, it is
one of the least-irradiated hot Jupiters observed with the Spitzer Space
Telescope. We have analyzed six photometric light curves of WASP-8b during
secondary eclipse observed in the 3.6, 4.5, and 8.0 {\microns} Infrared Array
Camera bands. The eclipse depths are $0.113\pm 0.018$%, $0.069\pm 0.007$%, and
$0.093\pm 0.023$%, respectively, giving respective brightness temperatures of
1552, 1131, and 938 K. We characterized the atmospheric thermal profile and
composition of the planet using a line-by-line radiative transfer code and a
Markov-chain Monte Carlo sampler. The data indicated no thermal inversion,
independently of any assumption about chemical composition. We noted an
anomalously high 3.6-{\microns} brightness temperature (1552 K); by modeling
the eccentricity-caused thermal variation, we found that this temperature is
plausible for radiative time scales less than $\sim 10^{2}$ hours. However, as
no model spectra fit all three data points well, the temperature discrepancy
remains as an open question.
[17]
oai:arXiv.org:1203.6017 [pdf] - 1511785
Thermal emission from WASP-24b at 3.6 and 4.5 {\mu}m
Smith, A. M. S.;
Anderson, D. R.;
Madhusudhan, N.;
Southworth, J.;
Cameron, A. Collier;
Blecic, J.;
Harrington, J.;
Hellier, C.;
Maxted, P. F. L.;
Pollacco, D.;
Queloz, D.;
Smalley, B.;
Triaud, A. H. M . J.;
Wheatley, P. J.
Submitted: 2012-03-27, last modified: 2012-08-20
Aims. We observe occultations of WASP-24b to measure brightness temperatures
and to determine whether or not its atmosphere exhibits a thermal inversion
(stratosphere). Methods. We observed occultations of WASP-24b at 3.6 and 4.5
{\mu}m using the Spitzer Space Telescope. It has been suggested that there is a
correlation between stellar activity and the presence of inversions, so we
analysed existing HARPS spectra in order to calculate log R'HK for WASP-24 and
thus determine whether or not the star is chromospherically active. We also
observed a transit of WASP-24b in the Str\"{o}mgren u and y bands, with the
CAHA 2.2-m telescope. Results. We measure occultation depths of 0.159 \pm 0.013
per cent at 3.6 {\mu}m and 0.202 \pm 0.018 per cent at 4.5 {\mu}m. The
corresponding planetary brightness temperatures are 1974 \pm 71 K and 1944 \pm
85 K respectively. Atmosphere models with and without a thermal inversion fit
the data equally well; we are unable to constrain the presence of an inversion
without additional occultation measurements in the near-IR. We find log R'HK =
-4.98 \pm 0.12, indicating that WASP-24 is not a chromospherically active star.
Our global analysis of new and previously-published data has refined the system
parameters, and we find no evidence that the orbit of WASP-24b is non-circular.
Conclusions. These results emphasise the importance of complementing Spitzer
measurements with observations at shorter wavelengths to gain a full
understanding of hot Jupiter atmospheres.
[18]
oai:arXiv.org:1207.4245 [pdf] - 859324
Two nearby sub-Earth-sized exoplanet candidates in the GJ 436 system
Stevenson, Kevin B.;
Harrington, Joseph;
Lust, Nate B.;
Lewis, Nikole K.;
Montagnier, Guillaume;
Moses, Julianne I.;
Visscher, Channon;
Blecic, Jasmina;
Hardy, Ryan A.;
Cubillos, Patricio;
Campo, Christopher J.
Submitted: 2012-07-17
We report the detection of UCF-1.01, a strong exoplanet candidate with a
radius 0.66 +/- 0.04 times that of Earth (R_{\oplus}). This sub-Earth-sized
planet transits the nearby M-dwarf star GJ 436 with a period of 1.365862 +/-
8x10^{-6} days. We also report evidence of a 0.65 +/- 0.06 R_{\oplus} exoplanet
candidate (labeled UCF-1.02) orbiting the same star with an undetermined
period. Using the Spitzer Space Telescope, we measure the dimming of light as
the planets pass in front of their parent star to assess their sizes and
orbital parameters. If confirmed, UCF-1.01 and UCF-1.02 would be called GJ 436c
and GJ 436d, respectively, and would be part of the first
multiple-transiting-planet system outside of the Kepler field. Assuming
Earth-like densities of 5.515 g/cm^3, we predict both candidates to have
similar masses (~0.28 Earth-masses, M_{\oplus}, 2.6 Mars-masses) and surface
gravities of ~0.65 g (where g is the gravity on Earth). UCF-1.01's equilibrium
temperature (T_{eq}, where emitted and absorbed radiation balance for an
equivalent blackbody) is 860 K, making the planet unlikely to harbor life as on
Earth. Its weak gravitational field and close proximity to its host star imply
that UCF-1.01 is unlikely to have retained its original atmosphere; however, a
transient atmosphere is possible if recent impacts or tidal heating were to
supply volatiles to the surface. We also present additional observations of GJ
436b during secondary eclipse. The 3.6-micron light curve shows indications of
stellar activity, making a reliable secondary eclipse measurement impossible. A
second non-detection at 4.5 microns supports our previous work in which we find
a methane-deficient and carbon monoxide-rich dayside atmosphere.
[19]
oai:arXiv.org:1206.0774 [pdf] - 1123887
Infrared Eclipses of the Strongly Irradiated Planet WASP-33b, and
Oscillations of its Host Star
Deming, Drake;
Fraine, Jonathan D.;
Sada, Pedro V.;
Madhusudhan, Nikku;
Knutson, Heather A.;
Harrington, Joseph;
Blecic, Jasmina;
Nymeyer, Sarah;
Smith, Alexis M. S.;
Jackson, Brian
Submitted: 2012-06-04
We observe two secondary eclipses of the strongly irradiated transiting
planet WASP-33b in the Ks band, and one secondary eclipse each at 3.6- and 4.5
microns using Warm Spitzer. This planet orbits an A5V delta-Scuti star that is
known to exhibit low amplitude non-radial p-mode oscillations at about
0.1-percent semi-amplitude. We detect stellar oscillations in all of our
infrared eclipse data, and also in one night of observations at J-band out of
eclipse. The oscillation amplitude, in all infrared bands except Ks, is about
the same as in the optical. However, the stellar oscillations in Ks band have
about twice the amplitude as seen in the optical, possibly because the
Brackett-gamma line falls in this bandpass. We use our best-fit values for the
eclipse depth, as well as the 0.9 micron eclipse observed by Smith et al., to
explore possible states of the exoplanetary atmosphere, based on the method of
Madhusudhan and Seager. On this basis we find two possible states for the
atmospheric structure of WASP-33b. One possibility is a non-inverted
temperature structure in spite of the strong irradiance, but this model
requires an enhanced carbon abundance (C/O>1). The alternative model has solar
composition, but an inverted temperature structure. Spectroscopy of the planet
at secondary eclipse, using a spectral resolution that can resolve the water
vapor band structure, should be able to break the degeneracy between these very
different possible states of the exoplanetary atmosphere. However, both of
those model atmospheres absorb nearly all of the stellar irradiance with
minimal longitudinal re-distribution of energy, strengthening the hypothesis of
Cowan et al. that the most strongly irradiated planets circulate energy poorly.
Our measurement of the central phase of the eclipse yields e*cos(omega)=0.0003
+/-0.00013, which we regard as being consistent with a circular orbit.
[20]
oai:arXiv.org:1005.1017 [pdf] - 955186
Spitzer Secondary Eclipses of WASP-18b
Nymeyer, Sarah;
Harrington, Joseph;
Hardy, Ryan A.;
Stevenson, Kevin B.;
Campo, Christopher J.;
Madhusudhan, Nikku;
Collier-Cameron, Andrew;
Loredo, Thomas J.;
Blecic, Jasmina;
Bowman, William C.;
Britt, Christopher B. T.;
Cubillos, Patricio;
Hellier, Coel;
Gillon, Michael;
Maxted, Pierre F. L.;
Hebb, Leslie;
Wheatley, Peter J.;
Pollacco, Don;
Anderson, David R.
Submitted: 2010-05-06, last modified: 2011-08-16
The transiting exoplanet WASP-18b was discovered in 2008 by the Wide Angle
Search for Planets (WASP) project. The Spitzer Exoplanet Target of Opportunity
Program observed secondary eclipses of WASP-18b using Spitzer's Infrared Array
Camera (IRAC) in the 3.6 micron and 5.8 micron bands on 2008 December 20, and
in the 4.5 micron and 8.0 micron bands on 2008 December 24. We report eclipse
depths of 0.30 +/- 0.02%, 0.39 +/- 0.02%, 0.37 +/- 0.03%, 0.41 +/- 0.02%, and
brightness temperatures of 3100 +/- 90, 3310 +/- 130, 3080 +/- 140 and 3120 +/-
110 K in order of increasing wavelength. WASP-18b is one of the hottest planets
yet discovered - as hot as an M-class star. The planet's pressure-temperature
profile most likely features a thermal inversion. The observations also require
WASP-18b to have near-zero albedo and almost no redistribution of energy from
the day-side to the night side of the planet.
[21]
oai:arXiv.org:1003.2763 [pdf] - 1025697
On the Orbit of Exoplanet WASP-12b
Campo, Christopher J.;
Harrington, Joseph;
Hardy, Ryan A.;
Stevenson, Kevin B.;
Nymeyer, Sarah;
Ragozzine, Darin;
Lust, Nate B.;
Anderson, David R.;
Collier-Cameron, Andrew;
Blecic, Jasmina;
Britt, Christopher B. T.;
Bowman, William C.;
Wheatley, Peter J.;
Loredo, Thomas J.;
Deming, Drake;
Hebb, Leslie;
Hellier, Coel;
Maxted, Pierre F. L.;
Pollaco, Don;
West, Richard G.
Submitted: 2010-03-14, last modified: 2010-12-09
We observed two secondary eclipses of the exoplanet WASP-12b using the
Infrared Array Camera on the Spitzer Space Telescope. The close proximity of
WASP-12b to its G-type star results in extreme tidal forces capable of inducing
apsidal precession with a period as short as a few decades. This precession
would be measurable if the orbit had a significant eccentricity, leading to an
estimate of the tidal Love number and an assessment of the degree of central
concentration in the planetary interior. An initial ground-based secondary
eclipse phase reported by Lopez-Morales et al. (0.510 +/- 0.002) implied
eccentricity at the 4.5 sigma level. The spectroscopic orbit of Hebb et al. has
eccentricity 0.049 +/- 0.015, a 3 sigma result, implying an eclipse phase of
0.509 +/- 0.007. However, there is a well documented tendency of spectroscopic
data to overestimate small eccentricities. Our eclipse phases are 0.5010 +/-
0.0006 (3.6 and 5.8 microns) and 0.5006 +/- 0.0007 (4.5 and 8.0 microns). An
unlikely orbital precession scenario invoking an alignment of the orbit during
the Spitzer observations could have explained this apparent discrepancy, but
the final eclipse phase of Lopez-Morales et al. (0.510 -0.006 / +0.007) is
consistent with a circular orbit at better than 2 sigma. An orbit fit to all
the available transit, eclipse, and radial-velocity data indicates precession
at <1 sigma; a non-precessing solution fits better. We also comment on analysis
and reporting for Spitzer exoplanet data in light of recent re-analyses.
[22]
oai:arXiv.org:1012.1603 [pdf] - 1042430
High C/O Ratio and Weak Thermal Inversion in the Very Hot Atmosphere of
Exoplanet WASP-12b
Madhusudhan, Nikku;
Harrington, Joseph;
Stevenson, Kevin B.;
Nymeyer, Sarah;
Campo, Christopher J.;
Wheatley, Peter J.;
Deming, Drake;
Blecic, Jasmina;
Hardy, Ryan A.;
Lust, Nate B.;
Anderson, David R.;
Collier-Cameron, Andrew;
Britt, Christopher B. T.;
Bowman, William C.;
Hebb, Leslie;
Hellier, Coel;
Maxted, Pierre F. L.;
Pollacco, Don;
West, Richard G.
Submitted: 2010-12-07
The carbon-to-oxygen ratio (C/O) in a planet provides critical information
about its primordial origins and subsequent evolution. A primordial C/O greater
than 0.8 causes a carbide-dominated interior as opposed to a silicate-dominated
composition as found on Earth, and the atmospheres can also differ from those
in the Solar System. The solar C/O is 0.54. Here we report an analysis of
dayside multi-wavelength photometry of the transiting hot-Jupiter WASP-12b that
reveals C/O >= 1 in its atmosphere. The atmosphere is abundant in CO. It is
depleted in water vapor and enhanced in methane by over two orders of magnitude
each compared to a solar-abundance chemical-equilibrium model at the expected
temperatures. We also find that the extremely irradiated atmosphere (T > 2,500
K) of WASP-12b lacks a prominent thermal inversion, or a stratosphere, and has
very efficient day- night energy circulation. The absence of a strong thermal
inversion is in stark contrast to theoretical predictions for the most highly
irradiated hot-Jupiter atmospheres.