Normalized to: Malik, M.
[1]
oai:arXiv.org:1910.01567 [pdf] - 2027894
Evidence for H$_{2}$ Dissociation and Recombination Heat Transport in
the Atmosphere of KELT-9b
Mansfield, Megan;
Bean, Jacob L.;
Stevenson, Kevin B.;
Komacek, Thaddeus D.;
Bell, Taylor J.;
Tan, Xianyu;
Malik, Matej;
Beatty, Thomas G.;
Wong, Ian;
Cowan, Nicolas B.;
Dang, Lisa;
Désert, Jean-Michel;
Fortney, Jonathan J.;
Gaudi, B. Scott;
Keating, Dylan;
Kempton, Eliza M. -R.;
Kreidberg, Laura;
Line, Michael R.;
Parmentier, Vivien;
Stassun, Keivan G.;
Swain, Mark R.;
Zellem, Robert T.
Submitted: 2019-10-03, last modified: 2020-01-07
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 $\mu$m. KELT-9b is the hottest
known transiting planet, with a 4.5-$\mu$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\sigma$ confidence. This discrepancy may
be due to magnetic effects in the planet's highly ionized atmosphere.
[2]
oai:arXiv.org:1910.11795 [pdf] - 2026138
Supervised Machine Learning for Intercomparison of Model Grids of Brown
Dwarfs: Application to GJ 570D and the Epsilon Indi B Binary System
Submitted: 2019-10-25, last modified: 2019-12-18
Self-consistent model grids of brown dwarfs involve complex physics and
chemistry, and are often computed using proprietary computer codes, making it
challenging to identify the reasons for discrepancies between model and data as
well as between the models produced by different research groups. In the
current study, we demonstrate a novel method for analyzing brown dwarf spectra,
which combines the use of the Sonora, AMES-Cond and HELIOS model grids with the
supervised machine learning method of the random forest. Besides performing
atmospheric retrieval, the random forest enables information content analysis
of the three model grids as a natural outcome of the method, both individually
on each grid and by comparing the grids against one another, via computing
large suites of mock retrievals. Our analysis reveals that the different
choices made in modelling the alkali line shapes hinder the use of the alkali
lines as gravity indicators. Nevertheless, the spectrum longward of 1.2 micron
encodes enough information on the surface gravity to allow its inference from
retrieval. Temperature may be accurately and precisely inferred independent of
the choice of model grid, but not the surface gravity. We apply random forest
retrieval to three objects: the benchmark T7.5 brown dwarf GJ 570D; and Epsilon
Indi Ba (T1.5 brown dwarf) and Bb (T6 brown dwarf), which are part of a binary
system and have measured dynamical masses. For GJ 570D, the inferred effective
temperature and surface gravity are consistent with previous studies. For
Epsilon Indi Ba and Bb, the inferred surface gravities are broadly consistent
with the values informed by the dynamical masses.
[3]
oai:arXiv.org:1907.13135 [pdf] - 2010011
Analyzing Atmospheric Temperature Profiles and Spectra of M dwarf Rocky
Planets
Submitted: 2019-07-30, last modified: 2019-12-05
The James Webb Space Telescope (JWST) will open up the possibility of
comprehensively measuring the emission spectra of rocky exoplanets orbiting M
dwarfs to detect and characterize their atmospheres. In preparation for this
opportunity, we present model atmospheres for three M-dwarf planets
particularly amenable to secondary eclipse spectroscopy -- TRAPPIST-1b, GJ
1132b, and LHS 3844b. Using three limiting cases of candidate atmospheric
compositions (pure H2O, pure CO2 and solar abundances) we calculate
temperature-pressure profiles and emission and reflection spectra in
radiative-convective equilibrium, including the effects of a solid surface. We
find that the atmospheric radiative transfer is significantly influenced by the
cool M-star irradiation; H2O and CO2 absorption bands in the near-infrared are
strong enough to absorb a sizeable fraction of the incoming stellar light at
low pressures, which leads to temperature inversions in the upper atmosphere.
The non-gray band structure of gaseous opacities in the infrared is hereby an
important factor. Opacity windows are muted at higher atmospheric temperatures,
so we expect temperature inversions to be common only for sufficiently cool
planets. We also find that pure CO2 atmospheres exhibit lower overall
temperatures and stronger reflection spectra compared to models of the other
two compositions. We estimate that for GJ 1132b and LHS 3844b we should be able
to distinguish between different atmospheric compositions with JWST. The
emission lines from the predicted temperature inversions are currently hard to
measure, but high resolution spectroscopy with future Extremely Large
Telescopes may be able to detect them.
[4]
oai:arXiv.org:1907.13150 [pdf] - 2025658
Identifying Atmospheres on Rocky Exoplanets Through Inferred High Albedo
Submitted: 2019-07-30, last modified: 2019-11-22
The upcoming launch of the James Webb Space Telescope (JWST) means that we
will soon have the capability to characterize the atmospheres of rocky
exoplanets. However, it is still unknown whether such planets orbiting close to
M dwarf stars can retain their atmospheres, or whether high-energy irradiation
from the star will strip the gaseous envelopes from these objects. We present a
new method to detect an atmosphere on a synchronously rotating rocky exoplanet
around a K/M dwarf, by using thermal emission during secondary eclipse to infer
a high dayside albedo that could only be explained by bright clouds. Based on
calculations for plausible surface conditions, we conclude that a high albedo
could be unambiguously interpreted as a signal of an atmosphere for planets
with substellar temperatures of $T_{sub}=$ 410-1250 K. This range corresponds
to equilibrium temperatures of $T_{eq}=$ 300-880 K. We compare the inferred
albedos of eight possible planet surface compositions to cloud albedo
calculations. We determine that a layer of clouds with optical depths greater
than $\tau=0.5$-$7$, would have high enough albedos to be distinguishable from
a bare rock surface. This method of detecting an atmosphere on a rocky planet
is complementary to existing methods for detecting atmospheres, because it
provides a way to detect atmospheres with pressures below 1 bar (e.g. Mars),
which are too tenuous to transport significant heat but thick enough to host
high-albedo clouds.
[5]
oai:arXiv.org:1907.13138 [pdf] - 2025657
Identifying Candidate Atmospheres on Rocky M dwarf Planets via Eclipse
Photometry
Submitted: 2019-07-30
Most rocky planets in the galaxy orbit a cool host star, and there is large
uncertainty among theoretical models whether these planets can retain an
atmosphere. The James Webb Space Telescope (JWST) might be able to settle this
question empirically, but most proposals for doing so require large
observational effort because they are based on spectroscopy. Here we show that
infrared photometry of secondary eclipses could quickly identify "candidate"
atmospheres, by searching for rocky planets with atmospheres thick enough that
atmospheric heat transport noticeably reduces their dayside thermal emission
compared to that of a bare rock. For a planet amenable to atmospheric
follow-up, we find that JWST should be able to confidently detect the heat
redistribution signal of an O(1) bar atmosphere with one to two eclipses. One
to two eclipses is generally much less than the effort needed to infer an
atmosphere via transmission or emission spectroscopy. Candidate atmospheres can
be further validated via follow-up spectroscopy or phase curves. In addition,
because this technique is fast it could enable a first atmospheric survey of
rocky exoplanets with JWST. We estimate that the TESS mission will find ~100
planets that are too hot to be habitable but that can be quickly probed via
eclipse photometry. Knowing whether hot, rocky planets around M dwarfs have
atmospheres is important not only for understanding the evolution of
uninhabitable worlds: if atmospheres are common on hot planets, then cooler,
potentially habitable planets around M dwarfs are also likely to have
atmospheres.
[6]
oai:arXiv.org:1812.02027 [pdf] - 1859942
Evolutionary models of cold and low-mass planets: Cooling curves,
magnitudes, and detectability
Submitted: 2018-12-05, last modified: 2019-03-27
Future instruments like NIRCam and MIRI on JWST or METIS at the ELT will be
able to image exoplanets that are too faint for current direct imaging
instruments. Evolutionary models predicting the planetary intrinsic luminosity
as a function of time have traditionally concentrated on gas-dominated giant
planets. We extend these cooling curves to Saturnian and Neptunian planets. We
simulate the cooling of isolated core-dominated and gas giant planets with
masses of 5 Earthmasses to 2 Jupitermasses. The luminosity includes the
contribution from the cooling and contraction of the core and of the H/He
envelope, as well as radiogenic decay. For the atmosphere we use grey,
AMES-Cond, petitCODE, and HELIOS models. We consider solar and non-solar
metallicities as well as cloud-free and cloudy atmospheres. The most important
initial conditions, namely the core-to-envelope ratio and the initial
luminosity are taken from planet formation simulations based on the core
accretion paradigm. We first compare our cooling curves for Uranus, Neptune,
Jupiter, Saturn, GJ 436b, and a 5 Earthmass-planet with a 1% H/He envelope with
other evolutionary models. We then present the temporal evolution of planets
with masses between 5 Earthmasses and 2 Jupitermasses in terms of their
luminosity, effective temperature, radius, and entropy. We discuss the impact
of different post formation entropies. For the different atmosphere types and
initial conditions magnitudes in various filter bands between 0.9 and 30
micrometer wavelength are provided. Using black body fluxes and non-grey
spectra, we estimate the detectability of such planets with JWST. It is found
that a 20 (100) Earthmass-planet can be detected with JWST in the background
limit up to an age of about 10 (100) Myr with NIRCam and MIRI, respectively.
[7]
oai:arXiv.org:1903.06794 [pdf] - 1868170
Self-luminous and irradiated exoplanetary atmospheres explored with
HELIOS
Submitted: 2019-03-15
We present new methodological features and physical ingredients included in
the 1D radiative transfer code HELIOS, improving the hemispheric two-stream
formalism. We conduct a thorough intercomparison survey with several
established forward models, including COOLTLUSTY, PHOENIX, and find
satisfactory consistency with their results. Then, we explore the impact of (i)
different groups of opacity sources, (ii) a stellar path length adjustment, and
(iii) a scattering correction on self-consistently calculated atmospheric
temperatures and planetary emission spectra. First, we observe that
temperature-pressure (T-P) profiles are very sensitive to the opacities
included, with metal oxides, hydrides, the alkali atoms (and ionized hydrogen)
playing an important role for the absorption of shortwave radiation (in very
hot surroundings). Moreover, if these species are sufficiently abundant, they
are likely to induce non-monotonic T-P profiles. Second, without the stellar
path length adjustment, the incoming stellar flux is significantly
underestimated for zenith angles above 80{\deg}, which somewhat affects the
upper atmospheric temperatures and the planetary emission. Third, the
scattering correction improves the accuracy of the computation of the reflected
stellar light by ~10%. We use HELIOS to calculate a grid of cloud-free
atmospheres in radiative-convective equilibrium for self-luminous planets for a
range of effective temperatures, surface gravities, metallicities, and C/O
ratios, to be used by planetary evolution studies. Furthermore, we calculate
dayside temperatures and secondary eclipse spectra for a sample of exoplanets
for varying chemistry and heat redistribution. These results may be used to
make predictions on the feasibility of atmospheric characterizations with
future observations.
[8]
oai:arXiv.org:1808.00501 [pdf] - 1806051
Three-Dimensional Circulation Driving Chemical Disequilibrium in
WASP-43b
Submitted: 2018-08-01, last modified: 2018-10-30
Spectral features in the observed spectra of exoplanets depend on the
composition of their atmospheres. A good knowledge of the main atmospheric
processes that drive the chemical distribution is therefore essential to
interpret exoplanetary spectra. An atmosphere reaches chemical equilibrium if
the rates of the forward and backward chemical reactions converge to the same
value. However, there are atmospheric processes, such as atmospheric transport,
that destabilize this equilibrium. In this work we study the changes in
composition driven by a 3D wind field in WASP-43b using our Global Circulation
Model, THOR. Our model uses validated temperature- and pressure-dependent
chemical timescales that allow us to explore the disequilibrium chemistry of
CO, CO$_2$, H$_2$O and CH$_4$. In WASP-43b the formation of the equatorial jet
has an important impact in the chemical distribution of the different species
across the atmosphere. At low latitudes the chemistry is longitudinally
quenched, except for CO$_2$ at solar abundances. The polar vortexes have a
distinct chemical distribution since these are regions with lower temperature
and atmospheric mixing. Vertical and latitudinal mixing have a secondary impact
in the chemical transport. We determine graphically the effect of
disequilibrium on observed emission spectra. Our results do not show any
significant differences in the emission spectra between the equilibrium and
disequilibrium solutions for C/O = 0.5. However, if C/O is increased to 2.0,
differences in the spectra due to the disequilibrium chemistry of CH$_4$ become
non-negligible. In some spectral ranges the emission spectra can have more than
15$\%$ departures from the equilibrium solution.
[9]
oai:arXiv.org:1803.07087 [pdf] - 1761472
Confronting phantom inflation with Planck data
Submitted: 2018-03-19, last modified: 2018-09-25
The latest Planck results are in excellent agreement with the theoretical
expectations predicted from standard normal inflation based on slow-roll
approximation which assumes equation-of-state $\omega \geq-1$. In this work, we
study the phantom inflation ($\omega<-1$) as an alternative cosmological model
within the slow-climb approximation using two hybrid inflationary fields. We
perform Chain Monte Carlo analysis to determine the posterior distribution and
best fit values for the cosmological parameters using Planck data and show that
current CMB data does not discriminate between normal and phantom inflation.
Interestingly, unlike in normal inflation, $\omega$ in phantom induced
inflation evolves very slowly away from $-1$ during the inflation. Furthermore,
in contrast to the standard normal inflation for which only upper bound on
tensor-to-scalar ratio $r$ are possible, we obtain both upper and lower bounds
for the two hybrid fields in the phantom scenario. Finally, we discuss
prospects of future high precision polarization measurements and show that it
may be possible to establish the dominance of one model over the other.
[10]
oai:arXiv.org:1807.10769 [pdf] - 1775614
Diagnosing aerosols in extrasolar giant planets with cross-correlation
function of water bands
Submitted: 2018-07-27
Transmission spectroscopy with ground-based, high-resolution instruments
provides key insight into the composition of exoplanetary atmospheres.
Molecules such as water and carbon monoxide have been unambiguously identified
in hot gas giants through cross- correlation techniques. A maximum in the
cross-correlation function (CCF) is found when the molecular absorption lines
in a binary mask or model template match those contained in the planet. Here,
we demonstrate how the CCF method can be used to diagnose broad-band
spectroscopic features such as scattering by aerosols in high-resolution
transit spectra. The idea consists in exploiting the presence of multiple water
bands from the optical to the near-infrared. We produced a set of models of a
typical hot Jupiter spanning various conditions of temperature and aerosol
coverage. We demonstrate that comparing the CCFs of individual water bands for
the models constrains the presence and the properties of the aerosol layers.
The contrast difference between the CCFs of two bands can reach about 100 ppm,
which could be readily detectable with current or upcoming high-resolution
stabilized spectrographs spanning a wide spectral range, such as ESPRESSO,
CARMENES, HARPS-N+GIANO, HARPS+NIRPS, SPIRou or CRIRES+.
[11]
oai:arXiv.org:1804.07137 [pdf] - 1739927
The Peculiar Atmospheric Chemistry of KELT-9b
Submitted: 2018-04-19, last modified: 2018-07-12
The atmospheric temperatures of the ultra-hot Jupiter KELT-9b straddle the
transition between gas giants and stars, and therefore between two
traditionally distinct regimes of atmospheric chemistry. Previous theoretical
studies assume the atmosphere of KELT-9b to be in chemical equilibrium. Despite
the high ultraviolet flux from KELT-9, we show using photochemical kinetics
calculations that the observable atmosphere of KELT-9b is predicted to be close
to chemical equilibrium, which greatly simplifies any theoretical
interpretation of its spectra. It also makes the atmosphere of KELT-9b, which
is expected to be cloudfree, a tightly constrained chemical system that lends
itself to a clean set of theoretical predictions. Due to the lower pressures
probed in transmission (compared to emission) spectroscopy, we predict the
abundance of water to vary by several orders of magnitude across the
atmospheric limb depending on temperature, which makes water a sensitive
thermometer. Carbon monoxide is predicted to be the dominant molecule under a
wide range of scenarios, rendering it a robust diagnostic of the metallicity
when analyzed in tandem with water. All of the other usual suspects (acetylene,
ammonia, carbon dioxide, hydrogen cyanide, methane) are predicted to be
subdominant at solar metallicity, while atomic oxygen, iron and magnesium are
predicted to have relative abundances as high as 1 part in 10,000. Neutral
atomic iron is predicted to be seen through a forest of optical and
near-infrared lines, which makes KELT-9b suitable for high-resolution
ground-based spectroscopy with HARPS-N or CARMENES. We summarize future
observational prospects of characterizing the atmosphere of KELT-9b.
[12]
oai:arXiv.org:1804.04961 [pdf] - 1728934
Analytical Models of Exoplanetary Atmospheres. VI. Full Solutions for
Improved Two-stream Radiative Transfer Including Direct Stellar Beam
Submitted: 2018-04-13, last modified: 2018-07-03
Two-stream radiative transfer is a workhorse in the Earth, planetary and
exoplanetary sciences due to its simplicity and ease of implementation.
However, a longstanding limitation of the two-stream approximation is its
inaccuracy in the presence of medium-sized or large aerosols. This limitation
was lifted by the discovery of the improved two-stream technique, where the
accuracy of the scattering greenhouse effect is matched to that of multi-stream
calculations by construction. In the present study, we derive the full
solutions for improved two-stream radiative transfer, following its
introduction by Heng & Kitzmann (2017), and include contributions from the
direct stellar beam. The generalization of the original two-stream flux
solutions comes in the form of a correction factor, traditionally set to unity,
which is the ratio of a pair of first Eddington coefficients. We derive an
analytical expression for this correction factor and also provide a simple
fitting function for its ease of use by other workers. We prove that the direct
stellar beam is associated with a second Eddington coefficient that is on the
order of unity. Setting this second Eddington coefficient to 2/3 and
$1/\sqrt{3}$ reproduces the Eddington and quadrature closures, respectively,
associated with the direct beam. We use our improved two-stream solutions for
the fluxes to derive two-stream source function solutions for both the
intensity and fluxes.
[13]
oai:arXiv.org:1709.09678 [pdf] - 1678459
Combining low- to high-resolution transit spectroscopy of HD 189733b.
Linking the troposphere and the thermosphere of a hot gas giant
Pino, Lorenzo;
Ehrenreich, David;
Wyttenbach, Aurelien;
Bourrier, Vincent;
Nascimbeni, Valerio;
Heng, Kevin;
Grimm, Simon;
Lovis, Christophe;
Malik, Matej;
Pepe, Francesco;
Piotto, Giampaolo
Submitted: 2017-09-27, last modified: 2018-05-08
Space-borne low-to medium-resolution (R~10^2-10^3) transmission spectroscopy
of atmospheres detect the broadest spectral features (alkali doublets,
molecular bands, scattering), while high-resolution (R~10^5), ground-based
observations probe the sharpest features (cores of the alkali lines, molecular
lines).The two techniques differ by:(1) The LSF of ground-based observations is
10^3 times narrower than for space-borne observations;(2)Space-borne
transmission spectra probe up to the base of thermosphere, while ground-based
observations can reach pressures down to 10^(-11);(3)Space-borne observations
directly yield the transit depth of the planet, while ground-based observations
measure differences in the radius of the planet at different wavelengths.It is
challenging to combine both techniques.We develop a method to compare
theoretical models with observations at different resolutions.We introduce
PyETA, a line-by-line 1D radiative transfer code to compute transmission
spectra at R~10^6 (0.01 A) over a broad wavelength range.An hybrid forward
modeling/retrieval optimization scheme is devised to deal with the large
computational resources required by modeling a broad wavelength range (0.3-2
$\mu$m) at high resolution.We apply our technique to HD189733b.Here, HST
observations reveal a flattened spectrum due to scattering by aerosols, while
high-resolution ground-based HARPS observations reveal the sharp cores of
sodium lines.We reconcile these results by building models that reproduce
simultaneously both data sets, from the troposphere to the thermosphere. We
confirm:(1)the presence of scattering by tropospheric aerosols;(2)that the
sodium core feature is of thermospheric origin.Accounting for aerosols, the
sodium cores indicate T up to 10000K in the thermosphere.The precise value of
the thermospheric temperature is degenerate with the abundance of sodium and
altitude of the aerosol deck.
[14]
oai:arXiv.org:1802.03047 [pdf] - 1656310
Revisiting the Phase Curves of WASP-43b: Confronting Reanalyzed Spitzer
Data with Cloudy Atmospheres
Submitted: 2018-02-08
Recently acquired Hubble and Spitzer phase curves of the short-period hot
Jupiter WASP-43b make it an ideal target for confronting theory with data. On
the observational front, we re-analyze the 3.6 and 4.5 $\mu$m Spitzer phase
curves and demonstrate that our improved analysis better removes residual red
noise due to intra-pixel sensitivity, which leads to greater fluxes emanating
from the nightside of WASP-43b, thus reducing the tension between theory and
data. On the theoretical front, we construct cloudfree and cloudy atmospheres
of WASP-43b using our Global Circulation Model (GCM), THOR, which solves the
non-hydrostatic Euler equations (compared to GCMs that typically solve the
hydrostatic primitive equations). The cloudfree atmosphere produces a
reasonable fit to the dayside emission spectrum. The multi-phase emission
spectra constrain the cloud deck to be confined to the nightside and have a
finite cloud-top pressure. The multi-wavelength phase curves are naturally
consistent with our cloudy atmospheres, except for the 4.5 $\mu$m phase curve,
which requires the presence of enhanced carbon dioxide in the atmosphere of
WASP-43b. Multi-phase emission spectra at higher spectral resolution, as may be
obtained using the James Webb Space Telescope, and a reflected-light phase
curve at visible wavelengths would further constrain the properties of clouds
in WASP-43b.
[15]
oai:arXiv.org:1705.04444 [pdf] - 1604213
AGN feedback with the Square Kilometer Array (SKA) and implications for
cluster physics and cosmology
Submitted: 2017-05-12, last modified: 2017-10-17
AGN feedback is regarded as an important non-gravitational process in galaxy
clusters, providing useful constraints on large-scale structure formation. It
modifies the structure and energetics of the intra-cluster medium (ICM) and
hence its understanding is crucially needed in order to use clusters as high
precision cosmological probes. In this context, particularly keeping in mind
the upcoming high quality radio data expected from radio surveys like SKA with
its higher sensitivity, high spatial and spectral resolutions, we review our
current understanding of AGN feedback, its cosmological implications and the
impact that SKA can have in revolutionizing our understanding of AGN feedback
in large-scale structures. Recent developments regarding the AGN outbursts and
its possible contribution to excess entropy in the hot atmospheres of groups
and clusters, its correlation with the feedback energy in ICM, quenching of
cooling flows and the possible connection between cool core clusters and radio
mini-halos, are discussed. We describe current major issues regarding modeling
of AGN feedback and its impact on the surrounding medium. With regard to the
future of AGN feedback studies, we examine the possible breakthroughs that can
be expected from SKA observations. In the context of cluster cosmology, for
example, we point out the importance of SKA observations for cluster mass
calibration by noting that most of $z>1$ clusters discovered by eROSITA X-ray
mission can be expected to be followed up through a 1000 hour SKA-1 mid
programme. Moreover, approximately $1000$ radio mini halos and $\sim 2500$
radio halos at $z<0.6$ can be potentially detected by SKA1 and SKA2 and used as
tracers of galaxy clusters and determination of cluster selection function.
[16]
oai:arXiv.org:1710.02427 [pdf] - 1589389
Aerosol Constraints on the Atmosphere of the Hot Saturn-mass planet
WASP-49b
Submitted: 2017-10-06
The strong, nearly wavelength-independent absorption cross section of
aerosols produces featureless exoplanet transmission spectra, limiting our
ability to characterize their atmospheres. Here we show that even in the
presence of featureless spectra, we can still characterize certain atmospheric
properties. Specifically, we constrain the upper and lower pressure boundaries
of aerosol layers, and present plausible composition candidates. We study the
case of the bloated Saturn-mass planet WASP-49b, where near-infrared
observations reveal a flat transmission spectrum between 0.7 and 1.0
{\microns}. First, we use a hydrodynamic upper-atmosphere code to estimate the
pressure reached by the ionizing stellar high-energy photons at $10^{-8}$ bar,
setting the upper pressure boundary where aerosols could exist. Then, we
combine HELIOS and Pyrat Bay radiative-transfer models to constrain the
temperature and photospheric pressure of atmospheric aerosols, in a Bayesian
framework. For WASP-49b, we constrain the transmission photosphere (hence, the
aerosol deck boundaries) to pressures above $10^{-5}$ bar (100$\times$ solar
metallicity), $10^{-4}$ bar (solar), and $10^{-3}$ bar (0.1$\times$ solar) as
lower boundary, and below $10^{-7}$ bar as upper boundary. Lastly, we compare
condensation curves of aerosol compounds with the planet's pressure-temperature
profile to identify plausible condensates responsible for the absorption. Under
these circumstances, we find as candidates: Na$_{2}$S (at 100$\times$ solar
metallicity); Cr and MnS (at solar and 0.1$\times$ solar); and forsterite,
enstatite, and alabandite (at 0.1$\times$ solar).
[17]
oai:arXiv.org:1709.00338 [pdf] - 1587807
Retrieval Analysis of the Emission Spectrum of WASP-12b: Sensitivity of
Outcomes to Prior Assumptions and Implications for Formation History
Oreshenko, Maria;
Lavie, Baptiste;
Grimm, Simon L.;
Tsai, Shang-Min;
Malik, Matej;
Demory, Brice-Olivier;
Mordasini, Christoph;
Alibert, Yann;
Benz, Willy;
Quanz, Sascha P.;
Trotta, Roberto;
Heng, Kevin
Submitted: 2017-09-01
We analyze the emission spectrum of the hot Jupiter WASP-12b using our
HELIOS-R retrieval code and HELIOS-K opacity calculator. When interpreting
Hubble and Spitzer data, the retrieval outcomes are found to be
prior-dominated. When the prior distributions of the molecular abundances are
assumed to be log-uniform, the volume mixing ratio of HCN is found to be
implausibly high. A VULCAN chemical kinetics model of WASP-12b suggests that
chemical equilibrium is a reasonable assumption even when atmospheric mixing is
implausibly rigorous. Guided by (exo)planet formation theory, we set Gaussian
priors on the elemental abundances of carbon, oxygen and nitrogen with the
Gaussian peaks being centered on the measured C/H, O/H and N/H values of the
star. By enforcing chemical equilibrium, we find substellar O/H and stellar to
slightly superstellar C/H for the dayside atmosphere of WASP-12b. The
superstellar carbon-to-oxygen ratio is just above unity, regardless of whether
clouds are included in the retrieval analysis, consistent with Madhusudhan et
al. (2011). Furthermore, whether a temperature inversion exists in the
atmosphere depends on one's assumption for the Gaussian width of the priors.
Our retrieved posterior distributions are consistent with the formation of
WASP-12b in a solar-composition protoplanetary disk, beyond the water iceline,
via gravitational instability or pebble accretion (without core erosion) and
migration inwards to its present orbital location via a disk-free mechanism,
and are inconsistent with both in-situ formation and core accretion with disk
migration, as predicted by Madhusudhan et al. (2017). We predict that the
interpretation of James Webb Space Telescope WASP-12b data will not be
prior-dominated.
[18]
oai:arXiv.org:1610.03216 [pdf] - 1580376
HELIOS-Retrieval: An Open-source, Nested Sampling Atmospheric Retrieval
Code, Application to the HR 8799 Exoplanets and Inferred Constraints for
Planet Formation
Lavie, Baptiste;
Mendonça, João M.;
Mordasini, Christoph;
Malik, Matej;
Bonnefoy, Mickaël;
Demory, Brice-Olivier;
Oreshenko, Maria;
Grimm, Simon L.;
Ehrenreich, David;
Heng, Kevin
Submitted: 2016-10-11, last modified: 2017-08-21
We present an open-source retrieval code named HELIOS-Retrieval (hereafter
HELIOS-R), designed to obtain chemical abundances and temperature-pressure
profiles from inverting the measured spectra of exoplanetary atmospheres. In
the current implementation, we use an exact solution of the radiative transfer
equation, in the pure absorption limit, in our forward model, which allows us
to analytically integrate over all of the outgoing rays (instead of performing
Gaussian quadrature). Two chemistry models are considered: unconstrained
chemistry (where the mixing ratios are treated as free parameters) and
equilibrium chemistry (enforced via analytical formulae, where only the
elemental abundances are free parameters). The nested sampling algorithm allows
us to formally implement Occam's Razor based on a comparison of the Bayesian
evidence between models. We perform a retrieval analysis on the measured
spectra of the HR 8799b, c, d and e directly imaged exoplanets. Chemical
equilibrium is disfavored by the Bayesian evidence for HR 8799b, c and d. We
find supersolar C/O, C/H and O/H values for the outer HR 8799b and c
exoplanets, while the inner HR 8799d and e exoplanets have substellar C/O,
substellar C/H and superstellar O/H values. If these retrieved properties are
representative of the bulk compositions of the exoplanets, then they are
inconsistent with formation via gravitational instability (without late-time
accretion) and consistent with a core accretion scenario in which late-time
accretion of ices occurred differently for the inner and outer exoplanets. For
HR 8799e, we find that spectroscopy in the K band is crucial for constraining
C/O and C/H. HELIOS-R is publicly available as part of the Exoclimes Simulation
Platform (ESP; www.exoclime.org).
[19]
oai:arXiv.org:1703.00028 [pdf] - 1581687
Excess entropy and energy feedback from within cluster cores up to
r$_{200}$
Submitted: 2017-02-28, last modified: 2017-08-03
We estimate the "non-gravitational" entropy-injection profiles, $\Delta K$,
and the resultant energy feedback profiles, $\Delta E$, of the intracluster
medium for 17 clusters using their Planck SZ and ROSAT X-Ray observations,
spanning a large radial range from $0.2r_{500}$ up to $r_{200}$. The feedback
profiles are estimated by comparing the observed entropy, at fixed gas mass
shells, with theoretical entropy profiles predicted from non-radiative
hydrodynamic simulations. We include non-thermal pressure and gas clumping in
our analysis. The inclusion of non-thermal pressure and clumping results in
changing the estimates for $r_{500}$ and $r_{200}$ by 10\%-20\%. When
clumpiness is not considered it leads to an under-estimation of $\Delta
K\approx300$ keV cm$^2$ at $r_{500}$ and $\Delta K\approx1100$ keV cm$^2$ at
$r_{200}$. On the other hand, neglecting non-thermal pressure results in an
over-estimation of $\Delta K\approx 100$ keV cm$^2$ at $r_{500}$ and
under-estimation of $\Delta K\approx450$ keV cm$^2$ at $r_{200}$. For the
estimated feedback energy, we find that ignoring clumping leads to an
under-estimation of energy per particle $\Delta E\approx1$ keV at $r_{500}$ and
$\Delta E\approx1.5$ keV at $r_{200}$. Similarly, neglect of the non-thermal
pressure results in an over-estimation of $\Delta E\approx0.5$ keV at $r_{500}$
and under-estimation of $\Delta E\approx0.25$ keV at $r_{200}$. We find entropy
floor of $\Delta K\approx300$ keV cm$^2$ is ruled out at $\approx3\sigma$
throughout the entire radial range and $\Delta E\approx1$ keV at more than
3$\sigma$ beyond $r_{500}$, strongly constraining ICM pre-heating scenarios. We
also demonstrate robustness of results w.r.t sample selection, X-Ray analysis
procedures, entropy modeling etc.
[20]
oai:arXiv.org:1610.05776 [pdf] - 1560179
Low-$\ell$ power suppression in punctuated inflation
Submitted: 2016-10-18, last modified: 2017-04-03
Motivated by Planck confirmation of an anomalously low value of the CMB
temperature fluctuations up to multipole $\ell<40$, we in this paper try to
explain such feature by investigating case of punctuated inflation scenario.
This form of inflation potential is inspired by Minimal Super-symmetric
Standard Model (MSSM) wherein suppression of curvature perturbation power at
large scales is produced by introducing period of fast-roll phase of the
inflation sandwiched between two stages of slow-roll phase.
We apply Markov Chain Monte Carlo analysis to determine posterior
distribution and the best fit values of the model parameters using recent WMAP9
and Planck data. We show that WMAP9 and Planck results are consistent with each
other and that with Planck data we obtain tighter constraints for punctuated
inflation parameters. We find that punctuated inflation leads to better fit in
CMB data compared to simple power law model. The improvement in the fit to the
WMAP9 data is $\Delta \chi^2\sim 3.6$ and for Planck the improvement is $\Delta
\chi^2\sim 5.4$. We find that $AIC$ does not discriminate between punctuated
inflation and simple power law model for WMAP9 data. However, for Planck data
we find that punctuated inflation is moderately preferred over a simple power
law model.
[21]
oai:arXiv.org:1606.05474 [pdf] - 1530855
HELIOS: An Open-source, GPU-accelerated Radiative Transfer Code For
Self-consistent Exoplanetary Atmospheres
Malik, Matej;
Grosheintz, Luc;
Mendonça, João M.;
Grimm, Simon L.;
Lavie, Baptiste;
Kitzmann, Daniel;
Tsai, Shang-Min;
Burrows, Adam;
Kreidberg, Laura;
Bedell, Megan;
Bean, Jacob L.;
Stevenson, Kevin B.;
Heng, Kevin
Submitted: 2016-06-17, last modified: 2016-11-19
We present the open-source radiative transfer code named HELIOS, which is
constructed for studying exoplanetary atmospheres. In its initial version, the
model atmospheres of HELIOS are one-dimensional and plane-parallel, and the
equation of radiative transfer is solved in the two-stream approximation with
non-isotropic scattering. A small set of the main infrared absorbers is
employed, computed with the opacity calculator HELIOS-K and combined using a
correlated-$k$ approximation. The molecular abundances originate from validated
analytical formulae for equilibrium chemistry. We compare HELIOS with the work
of Miller-Ricci & Fortney using a model of GJ 1214b, and perform several tests,
where we find: model atmospheres with single-temperature layers struggle to
converge to radiative equilibrium; $k$-distribution tables constructed with
$\gtrsim 0.01$ cm$^{-1}$ resolution in the opacity function ($ \lesssim 10^3$
points per wavenumber bin) may result in errors $\gtrsim 1$-10 % in the
synthetic spectra; and a diffusivity factor of 2 approximates well the exact
radiative transfer solution in the limit of pure absorption. We construct
"null-hypothesis" models (chemical equilibrium, radiative equilibrium and solar
element abundances) for 6 hot Jupiters. We find that the dayside emission
spectra of HD 189733b and WASP-43b are consistent with the null hypothesis,
while it consistently under-predicts the observed fluxes of WASP-8b, WASP-12b,
WASP-14b and WASP-33b. We demonstrate that our results are somewhat insensitive
to the choice of stellar models (blackbody, Kurucz or PHOENIX) and metallicity,
but are strongly affected by higher carbon-to-oxygen ratios. The code is
publicly available as part of the Exoclimes Simulation Platform (ESP;
exoclime.net).
[22]
oai:arXiv.org:1606.00014 [pdf] - 1522507
Little evidence for entropy and energy excess beyond $r_{500}$ - An end
to ICM preheating?
Submitted: 2016-05-31, last modified: 2016-10-27
Non-gravitational feedback affects the nature of the intra-cluster medium
(ICM). X-ray cooling of the ICM and in situ energy feedback from AGN's and SNe
as well as {\it preheating} of the gas at epochs preceding the formation of
clusters are proposed mechanisms for such feedback. While cooling and AGN
feedbacks are dominant in cluster cores, the signatures of a preheated ICM are
expected to be present even at large radii. To estimate the degree of
preheating, with minimum confusion from AGN feedback/cooling, we study the
excess entropy and non-gravitational energy profiles upto $r_{200}$ for a
sample of 17 galaxy clusters using joint data sets of {\it Planck} SZ pressure
and {\it ROSAT/PSPC} gas density profiles. The canonical value of preheating
entropy floor of $\gtrsim 300$ keV cm$^2$, needed in order to match cluster
scalings, is ruled out at $\approx 3\sigma$. We also show that the feedback
energy of 1 keV/particle is ruled out at 5.2$\sigma$ beyond $r_{500}$. Our
analysis takes both non-thermal pressure and clumping into account which can be
important in outer regions. Our results based on the direct probe of the ICM in
the outermost regions do not support any significant preheating.
[23]
oai:arXiv.org:1501.02647 [pdf] - 1126203
Joint Planck and WMAP Assessment of Low CMB Multipoles
Submitted: 2015-01-12, last modified: 2015-06-08
The remarkable progress in cosmic microwave background (CMB) studies over
past decade has led to the era of precision cosmology in striking agreement
with the $\Lambda$CDM model. However, the lack of power in the CMB temperature
anisotropies at large angular scales (low-$\ell$), as has been confirmed by the
recent Planck data also (up to $\ell=40$), although statistically not very
strong (less than $3\sigma$), is still an open problem. One can avoid to seek
an explanation for this problem by attributing the lack of power to cosmic
variance orcan look for explanations i.e., different inflationary potentials or
initial conditions for infl ation to begin with, non-trivial topology, ISW
effect etc. Features in the primordial power spectrum (PPS) motivated by the
early universe physics has been the most common solution to address this
problem. In the present work we also follow this approach and consider a set of
PPS which have features and constrain the parameters of those using WMAP 9 year
and Planck data employing Markov-Chain Monte Carlo (MCMC) analysis. The
prominent feature of all the models of PPS that we consider is an infra-red cut
off which leads to suppression of power at large angular scales. We consider
models of PPS with maximum three extra parameters and use Akaike information
criterion ($AIC$) and Bayesian information criterion ($BIC$) of model and
Bayesian information criterion ($BIC$) of model selection to compare the
models. For most models, we find good constraints for the cut off scale $k_c$,
however, for other parameters our constraints are not that good. We find that
sharp cut off model gives best likelihood value for the WMAP 9 year data, but
is as good as power law model according to $AIC$. For the joint WMAP 9+Planck
data set, Starobinsky model is slightly preferred by $AIC$ which is also able
to produce CMB power suppression up to $\ell\leq30$ to some extent.
[24]
oai:arXiv.org:1502.06597 [pdf] - 1231836
On the gap-opening criterion of migrating planets in protoplanetary
disks
Submitted: 2015-02-23
We perform two-dimensional hydrodynamical simulations to quantitatively
explore the torque balance criterion for gap-opening (as formulated by Crida et
al. 2006) in a variety of disks when considering a migrating planet. We find
that even when the criterion is satisfied, there are instances when planets
still do not open gaps. We stress that gap-opening is not only dependent on
whether a planet has the ability to open a gap, but whether it can do so
quickly enough. This can be expressed as an additional condition on the
gap-opening timescale versus the crossing time, i.e. the time it takes the
planet to cross the region which it is carving out. While this point has been
briefly made in the previous literature, our results quantify it for a range of
protoplanetary disk properties and planetary masses, demonstrating how crucial
it is for gap-opening. This additional condition has important implications for
the survival of planets formed by core accretion in low mass disks as well as
giant planets or brown dwarfs formed by gravitational instability in massive
disks. It is particularly important for planets with intermediate masses
susceptible to Type III-like migration. For some observed transition disks or
disks with gaps, we expect that estimates on the potential planet masses based
on the torque balance gap-opening criterion alone may not be sufficient. With
consideration of this additional timescale criterion theoretical studies may
find a reduced planet survivability or that planets may migrate further inwards
before opening a gap.
[25]
oai:arXiv.org:1304.3125 [pdf] - 1165885
Debris discs around M stars: non-existence versus non-detection
Submitted: 2013-04-10
Motivated by the reported dearth of debris discs around M stars, we use
survival models to study the occurrence of planetesimal discs around them.
These survival models describe a planetesimal disc with a small number of
parameters, determine if it may survive a series of dynamical processes and
compute the associated infrared excess. For the WISE satellite, we demonstrate
that the dearth of debris discs around M stars may be attributed to the small
semi-major axes generally probed if either: 1. the dust grains behave like
blackbodies emitting at a peak wavelength coincident with the observed one; 2.
or the grains are hotter than predicted by their blackbody temperatures and
emit at peak wavelengths that are shorter than the observed one. At these small
distances from the M star, planetesimals are unlikely to survive or persist for
time scales of 300 Myr or longer if the disc is too massive. Conversely, our
survival models allow for the existence of a large population of low-mass
debris discs that are too faint to be detected with current instruments.
However, our interpretation becomes less clear and large infrared excesses are
allowed if only one of these scenarios holds: 3. the dust grains are hotter
than blackbody and predominantly emit at the observed wavelength; 4. or are
blackbody in nature and emit at peak wavelengths longer than the observed one.
Both scenarios imply that the parent planetesimals reside at larger distances
from the star than inferred if the dust grains behaved like blackbodies. In all
scenarios, we show that the infrared excesses detected at 22 and 70 microns
from AU Mic are easily reconciled with its young age. We elucidate the
conditions under which stellar wind drag may be neglected when considering dust
populations around M stars. The WISE satellite should be capable of detecting
debris discs around young M stars with ages on the order of 10 Myr.
[26]
oai:arXiv.org:1011.5170 [pdf] - 1042156
MARVELS-1b: A Short-Period, Brown Dwarf Desert Candidate from the
SDSS-III MARVELS Planet Search
Lee, Brian L.;
Ge, Jian;
Fleming, Scott W.;
Stassun, Keivan G.;
Gaudi, B. Scott;
Barnes, Rory;
Mahadevan, Suvrath;
Eastman, Jason D.;
Wright, Jason;
Siverd, Robert J.;
Gary, Bruce;
Ghezzi, Luan;
Laws, Chris;
Wisniewski, John P.;
de Mello, G. F. Porto;
Ogando, Ricardo L. C.;
Maia, Marcio A. G.;
da Costa, Luiz Nicolaci;
Sivarani, Thirupathi;
Pepper, Joshua;
Nguyen, Duy Cuong;
Hebb, Leslie;
De Lee, Nathan;
Wang, Ji;
Wan, Xiaoke;
Zhao, Bo;
Chang, Liang;
Groot, John;
Varosi, Frank;
Hearty, Fred;
Hanna, Kevin;
van Eyken, J. C.;
Kane, Stephen R.;
Agol, Eric;
Bizyaev, Dmitry;
Bochanski, John J.;
Brewington, Howard;
Chen, Zhiping;
Costello, Erin;
Dou, Liming;
Eisenstein, Daniel J.;
Fletcher, Adam;
Ford, Eric B.;
Guo, Pengcheng;
Holtzman, Jon A.;
Jiang, Peng;
Leger, R. French;
Liu, Jian;
Long, Daniel C.;
Malanushenko, Elena;
Malanushenko, Viktor;
Malik, Mohit;
Oravetz, Daniel;
Pan, Kaike;
Rohan, Pais;
Schneider, Donald P.;
Shelden, Alaina;
Snedden, Stephanie A.;
Simmons, Audrey;
Weaver, B. A.;
Weinberg, David H.;
Xie, Ji-Wei
Submitted: 2010-11-23
We present a new short-period brown dwarf candidate around the star TYC
1240-00945-1. This candidate was discovered in the first year of the
Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which
is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we
designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of
dispersed fixed-delay interferometery to simultaneously obtain radial velocity
measurements for 60 objects per field using a single, custom-built instrument
that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity
measurements spread over a ~370 d time baseline, we derive a Keplerian orbital
fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and
eccentricity consistent with circular. Independent follow-up radial velocity
data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly
evolved F9 host star, we infer that the companion has a minimum mass of
28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit,
and is probably tidally synchronized. We find no evidence for coherent
instrinsic variability of the host star at the period of the companion at
levels greater than a few millimagnitudes. The companion has an a priori
transit probability of ~14%. Although we find no evidence for transits, we
cannot definitively rule them out for companion radii ~<1 R_Jup.