Normalized to: Grimm, S.
[2]
oai:arXiv.org:1911.13158 [pdf] - 2115136
THOR 2.0: Major Improvements to the Open-Source General Circulation
Model
Submitted: 2019-11-29, last modified: 2020-05-25
THOR is the first open-source general circulation model (GCM) developed from
scratch to study the atmospheres and climates of exoplanets, free from Earth-
or Solar System-centric tunings. It solves the general non-hydrostatic Euler
equations (instead of the primitive equations) on a sphere using the
icosahedral grid. In the current study, we report major upgrades to THOR,
building upon the work of Mendon\c{c}a et al. (2016). First, while the
Horizontally Explicit Vertically Implicit (HEVI) integration scheme is the same
as that described in Mendon\c{c}a et al. (2016), we provide a clearer
description of the scheme and improved its implementation in the code. The
differences in implementation between the hydrostatic shallow (HSS),
quasi-hydrostatic deep (QHD) and non-hydrostatic deep (NHD) treatments are
fully detailed. Second, standard physics modules are added: two-stream,
double-gray radiative transfer and dry convective adjustment. Third, THOR is
tested on additional benchmarks: tidally-locked Earth, deep hot Jupiter,
acoustic wave, and gravity wave. Fourth, we report that differences between the
hydrostatic and non-hydrostatic simulations are negligible in the Earth case,
but pronounced in the hot Jupiter case. Finally, the effects of the so-called
"sponge layer", a form of drag implemented in most GCMs to provide numerical
stability, are examined. Overall, these upgrades have improved the flexibility,
user-friendliness, and stability of THOR.
[3]
oai:arXiv.org:2004.10106 [pdf] - 2115276
Information content of JWST-NIRSPEC transmission spectra of warm
Neptunes
Guzmán-Mesa, Andrea;
Kitzmann, Daniel;
Fisher, Chloe;
Burgasser, Adam J.;
Hoeijmakers, H. Jens;
Márquez-Neila, Pablo;
Grimm, Simon L.;
Mandell, Avi M.;
Sznitman, Raphael;
Heng, Kevin
Submitted: 2020-04-21, last modified: 2020-05-07
Warm Neptunes offer a rich opportunity for understanding exo-atmospheric
chemistry. With the upcoming James Webb Space Telescope (JWST), there is a need
to elucidate the balance between investments in telescope time versus
scientific yield. We use the supervised machine learning method of the random
forest to perform an information content analysis on a 11-parameter model of
transmission spectra from the various NIRSpec modes. The three bluest
medium-resolution NIRSpec modes (0.7 - 1.27 microns, 0.97 - 1.84 microns, 1.66
- 3.07 microns) are insensitive to the presence of CO. The reddest
medium-resolution mode (2.87 - 5.10 microns) is sensitive to all of the
molecules assumed in our model: CO, CO2, CH4, C2H2, H2O, HCN and NH3. It
competes effectively with the three bluest modes on the information encoded on
cloud abundance and particle size. It is also competitive with the
low-resolution prism mode (0.6 - 5.3 microns) on the inference of every
parameter except for the temperature and ammonia abundance. We recommend
astronomers to use the reddest medium-resolution NIRSpec mode for studying the
atmospheric chemistry of 800-1200 K warm Neptunes; its corresponding
high-resolution counterpart offers diminishing returns. We compare our findings
to previous JWST information content analyses that favor the blue orders, and
suggest that the reliance on chemical equilibrium could lead to biased outcomes
if this assumption does not apply. A simple, pressure-independent diagnostic
for identifying chemical disequilibrium is proposed based on measuring the
abundances of H2O, CO and CO2.
[5]
oai:arXiv.org:2004.08415 [pdf] - 2080208
High-resolution Transmission Spectroscopy of MASCARA-2 b with EXPRES
Hoeijmakers, H. Jens;
Cabot, Samuel H. C.;
Zhao, Lily;
Buchhave, Lars A.;
Tronsgaard, René;
Kitzmann, Daniel;
Grimm, Simon L.;
Cegla, Heather M.;
Bourrier, Vincent;
Ehrenreich, David;
Heng, Kevin;
Lovis, Christophe;
Fischer, Debra A.
Submitted: 2020-04-17
We report detections of atomic species in the atmosphere of MASCARA-2 b,
using the first transit observations obtained with the newly commissioned
EXPRES spectrograph. EXPRES is a highly stabilised optical echelle
spectrograph, designed to detect stellar reflex motions with amplitudes down to
30 cm/s, and was recently deployed at the Lowell Discovery Telescope. By
analysing the transmission spectrum of the ultra-hot Jupiter MASCARA-2 b using
the cross-correlation method, we confirm previous detections of Fe I, Fe II and
Na I, which likely originate in the upper regions of the inflated atmosphere.
In addition, we report significant detections of Mg I and Cr II. The absorption
strengths change slightly with time, possibly indicating different temperatures
and chemistry in the day-side and night-side terminators. Using the effective
stellar line-shape variation induced by the transiting planet, we constrain the
projected spin-orbit misalignment of the system to $1.6\pm3.1$ degrees,
consistent with an aligned orbit. We demonstrate that EXPRES joins a suite of
instruments capable of phase-resolved spectroscopy of exoplanet atmospheres.
[7]
oai:arXiv.org:1912.02622 [pdf] - 2072606
Science Case for the Einstein Telescope
Maggiore, Michele;
Broeck, Chris van den;
Bartolo, Nicola;
Belgacem, Enis;
Bertacca, Daniele;
Bizouard, Marie Anne;
Branchesi, Marica;
Clesse, Sebastien;
Foffa, Stefano;
García-Bellido, Juan;
Grimm, Stefan;
Harms, Jan;
Hinderer, Tanja;
Matarrese, Sabino;
Palomba, Cristiano;
Peloso, Marco;
Ricciardone, Angelo;
Sakellariadou, Mairi
Submitted: 2019-12-05, last modified: 2020-03-24
The Einstein Telescope (ET), a proposed European ground-based
gravitational-wave detector of third-generation, is an evolution of
second-generation detectors such as Advanced LIGO, Advanced Virgo, and KAGRA
which could be operating in the mid 2030s. ET will explore the universe with
gravitational waves up to cosmological distances. We discuss its main
scientific objectives and its potential for discoveries in astrophysics,
cosmology and fundamental physics.
[9]
oai:arXiv.org:1910.11627 [pdf] - 2076699
Interpreting High-Resolution Spectroscopy of Exoplanets Using
Cross-Correlations and Supervised Machine Learning
Submitted: 2019-10-25, last modified: 2020-02-29
We present a new method for performing atmospheric retrieval on ground-based,
high-resolution data of exoplanets. Our method combines cross-correlation
functions with a random forest, a supervised machine learning technique, to
overcome challenges associated with high-resolution data. A series of
cross-correlation functions are concatenated to give a "CCF-sequence" for each
model atmosphere, which reduces the dimensionality by a factor of ~100. The
random forest, trained on our grid of ~65,000 models, provides a
likelihood-free method of retrieval. The pre-computed grid spans 31 values of
both temperature and metallicity, and incorporates a realistic noise model. We
apply our method to HARPS-N observations of the ultra-hot Jupiter KELT-9b, and
obtain a metallicity consistent with solar (logM = $-0.2\pm0.2$). Our retrieved
transit chord temperature (T = $6000^{+0}_{-200}$K) is unreliable as the ion
cross-correlations lie outside of the training set, which we interpret as being
indicative of missing physics in our atmospheric model. We compare our method
to traditional nested-sampling, as well as other machine learning techniques,
such as Bayesian neural networks. We demonstrate that the likelihood-free
aspect of the random forest makes it more robust than nested-sampling to
different error distributions, and that the Bayesian neural network we tested
is unable to reproduce complex posteriors. We also address the claim in Cobb et
al. (2019) that our random forest retrieval technique can be over-confident but
incorrect. We show that this is an artefact of the training set, rather than
the machine learning method, and that the posteriors agree with those obtained
using nested-sampling.
[11]
oai:arXiv.org:1910.01070 [pdf] - 2057751
Helios-r2 -- A new Bayesian, open-source retrieval model for brown
dwarfs and exoplanet atmospheres
Submitted: 2019-10-02, last modified: 2020-02-10
We present an improved, hybrid CPU-GPU atmospheric retrieval code, Helios-r2,
which is applicable to medium-resolution emission spectra of brown dwarfs, in
preparation for precision atmospheric spectroscopy in the era of the James Webb
Space Telescope. The model is available as open-source code on the Exoclimes
Simulation Platform. We subject Helios-r2 to a battery of tests of varying
difficulty. The simplest test involves a mock retrieval on a forward model
generated using the same radiative transfer technique, the same implementation
of opacities, and the same chemistry model. The least trivial test involves a
mock retrieval on synthetic spectra from the Sonora model grid, which uses a
different radiative transfer technique, a different implementation of
opacities, and a different chemistry model. A calibration factor, which is
included to capture uncertainties in the brown dwarf radius, distance to the
brown dwarf and flux calibration of the spectrum, may compensate, sometimes
erroneously, for discrepancies in modeling choices and implementation. We
analyze spectra of the benchmark brown dwarf GJ 570 D and the binary brown
dwarf companions in the Epsilon Indi system. The retrieved surface gravities
are consistent with previous studies and/or values inferred from dynamical
masses (for Epsilon Indi Ba and Bb only). There remains no clear criterion on
how to reject unphysical values of the retrieved brown dwarf radii. The
inferred radii and corresponding masses should be taken with great caution. The
retrieved carbon-to-oxygen ratios and metallicity depend on whether chemical
equilibrium is assumed.
[12]
oai:arXiv.org:2002.02015 [pdf] - 2069033
On the impact of tides on the transit-timing fits to the TRAPPIST-1
system
Submitted: 2020-02-05
Transit Timing Variations, or TTVs, can be a very efficient way of
constraining masses and eccentricities of multi-planet systems. Recent
measurements of the TTVs of TRAPPIST-1 led to an estimate of the masses of the
planets, enabling an estimate of their densities. A recent TTV analysis using
data obtained in the past two years yields a 34% and 13% increase in mass for
TRAPPIST-1b and c, respectively. In most studies to date, a Newtonian N-body
model is used to fit the masses of the planets, while sometimes general
relativity is accounted for. Using the Posidonius N-body code, we show that in
the case of the TRAPPIST-1 system, non-Newtonian effects might be also relevant
to correctly model the dynamics of the system and the resulting TTVs. In
particular, using standard values of the tidal Love number $k_2$ (accounting
for the tidal deformation) and the fluid Love number $k_{2f}$ (accounting for
the rotational flattening) leads to differences in the TTVs of TRAPPIST-1b and
c similar to the differences caused by general relativity. We also show that
relaxing the values of tidal Love number $k_2$ and the fluid Love number
$k_{2f}$ can lead to TTVs which differ by as much as a few 10~s on a $3-4$-year
timescale, which is a potentially observable level. The high values of the Love
numbers needed to reach observable levels for the TTVs could be achieved for
planets with a liquid ocean, which, if detected, might then be interpreted as a
sign that TRAPPIST-1b and TRAPPIST-1c could have a liquid magma ocean. For
TRAPPIST-1 and similar systems, the models to fit the TTVs should potentially
account for general relativity, for the tidal deformation of the planets, for
the rotational deformation of the planets and, to a lesser extent, for the
rotational deformation of the star, which would add up to 7x2+1 = 15 additional
free parameters in the case of TRAPPIST-1.
[19]
oai:arXiv.org:1905.04587 [pdf] - 1925039
ExoMol Molecular linelists -- XXXIII. The spectrum of Titanium Oxide
Submitted: 2019-05-11, last modified: 2019-07-12
Accurate line lists are crucial for correctly modelling a variety of
astrophysical phenomena, including stellar photospheres and the atmospheres of
extra-solar planets. This paper presents a new line database Toto for the main
isotopologues of titanium oxide (TiO): $^{46}$Ti$^{16}$O, $^{47}$Ti$^{16}$O,
$^{48}$Ti$^{16}$O, $^{49}$Ti$^{16}$O and $^{50}$Ti$^{16}$O. The TiO line list
contains transitions with wave-numbers up to 30,000 cm$^{-1}$ ie long-wards of
0.33 $\mu$m. The Toto line list includes all dipole-allowed transitions between
13 low-lying electronic states (X $^3\Delta$, a $^1\Delta$, d $^1\Sigma^+$, E
$^3\Pi$, A $^3\Phi$ B $^3\Pi$, C $^3\Delta$, b $^1\Pi$, c $^1\Phi$, f
$^1\Delta$, e $^1\Sigma^+$). Ab initio potential energy curves (PECs) are
computed at the icMRCI level and combined with spin-orbit and other coupling
curves. These PECs and couplings are iteratively refined to match known
empirical energy levels. Accurate line intensities are generated using ab
initio dipole moment curves. The Toto line lists are appropriate for
temperatures below 5000 K and contain 30 million transitions for TiO; it is
made available in electronic form via the CDS data centre and via
www.exomol.com. Tests of the line lists show greatly improved agreement with
observed spectra for objects such as M-dwarfs GJ876 and GL581.
[21]
oai:arXiv.org:1905.08132 [pdf] - 1890541
Exoplanetary Monte Carlo Radiative Transfer with Correlated-k I.
Benchmarking Transit and Emission Observables
Submitted: 2019-05-20
Current observational data of exoplanets are providing increasing detail of
their 3D atmospheric structures. As characterisation efforts expand in scope,
the need to develop consistent 3D radiative-transfer methods becomes more
pertinent as the complex atmospheric properties of exoplanets are required to
be modelled together consistently. We aim to compare the transmission and
emission spectra results of a 3D Monte Carlo Radiative Transfer (MCRT) model to
contemporary radiative-transfer suites. We perform several benchmarking tests
of a MCRT code, Cloudy Monte Carlo Radiative Transfer (CMCRT), to transmission
and emission spectra model output. We add flexibility to the model through the
use of k-distribution tables as input opacities. We present a hybrid MCRT and
ray tracing methodology for the calculation of transmission spectra with a
multiple scattering component. CMCRT compares well to the transmission spectra
benchmarks at the 10s of ppm level. Emission spectra benchmarks are consistent
to within 10% of the 1D models. We suggest that differences in the benchmark
results are likely caused by geometric effects between plane-parallel and
spherical models. In a practical application, we post-process a cloudy 3DHD
189733b GCM model and compare to available observational data. Our results
suggest the core methodology and algorithms of CMCRT produce consistent results
to contemporary radiative transfer suites. 3D MCRT methods are highly suitable
for detailed post-processing of cloudy and non-cloudy 1D and 3D exoplanet
atmosphere simulations in instances where atmospheric inhomogeneities,
significant limb effects/geometry or multiple scattering components are
important considerations.
[23]
oai:arXiv.org:1905.02096 [pdf] - 1920916
A spectral survey of an ultra-hot Jupiter: Detection of metals in the
transmission spectrum of KELT-9 b
Hoeijmakers, H. J.;
Ehrenreich, D.;
Kitzmann, D.;
Allart, R.;
Grimm, S. L.;
Seidel, J. V.;
Wyttenbach, A.;
Pino, L.;
Nielsen, L. D.;
Fisher, C.;
Rimmer, P. B.;
Bourrier, V.;
Cegla, H. M.;
Lavie, B.;
Lovis, C.;
Patzer, A. B. C.;
Stock, J. W.;
Pepe, F. A.;
Heng, Kevin
Submitted: 2019-05-06
Context: KELT-9 b exemplifies a newly emerging class of short-period gaseous
exoplanets that tend to orbit hot, early type stars - termed ultra-hot
Jupiters. The severe stellar irradiation heats their atmospheres to
temperatures of $\sim 4,000$ K, similar to the photospheres of dwarf stars. Due
to the absence of aerosols and complex molecular chemistry at such
temperatures, these planets offer the potential of detailed chemical
characterisation through transit and day-side spectroscopy. Studies of their
chemical inventories may provide crucial constraints on their formation process
and evolution history.
Aims: To search the optical transmission spectrum of KELT-9 b for absorption
lines by metals using the cross-correlation technique.
Methods: We analyse 2 transits observed with the HARPS-N spectrograph. We use
an isothermal equilibrium chemistry model to predict the transmission spectrum
for each of the neutral and singly-ionized atoms with atomic numbers between 3
and 78. Of these, we identify the elements that are expected to have spectral
lines in the visible wavelength range and use those as cross-correlation
templates.
Results: We detect absorption of Na I, Cr II, Sc II and Y II, and confirm
previous detections of Mg I, Fe I, Fe II and Ti II. In addition, we find
evidence of Ca I, Cr I, Co I, and Sr II that will require further observations
to verify. The detected absorption lines are significantly deeper than model
predictions, suggesting that material is transported to higher altitudes where
the density is enhanced compared to a hydrostatic profile. There appears to be
no significant blue-shift of the absorption spectrum due to a net day-to-night
side wind. In particular, the strong Fe II feature is shifted by $0.18 \pm
0.27$ km~s$^{-1}$, consistent with zero. Using the orbital velocity of the
planet we revise the steller and planetary masses and radii.
[24]
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.
[25]
oai:arXiv.org:1901.07250 [pdf] - 1863958
Co-orbital exoplanets from close period candidates: The TOI-178 case
Submitted: 2019-01-22
Despite the existence of co-orbital bodies in the solar system, and the
prediction of the formation of co-orbital planets by planetary system formation
models, no co-orbital exoplanets (also called trojans) have been detected thus
far. Here we study the signature of co-orbital exoplanets in transit surveys
when two planet candidates in the system orbit the star with similar periods.
Such pair of candidates could be discarded as false positives because they are
not Hill-stable. However, horseshoe or long libration period tadpole co-orbital
configurations can explain such period similarity. This degeneracy can be
solved by considering the Transit Timing Variations (TTVs) of each planet. We
then focus on the three planet candidates system TOI-178: the two outer
candidates of that system have similar orbital period and had an angular
separation near $\pi/3$ during the TESS observation of sector 2. Based on the
announced orbits, the long-term stability of the system requires the two
close-period planets to be co-orbitals. Our independent detrending and transit
search recover and slightly favour the three orbits close to a 3:2:2 resonant
chain found by the TESS pipeline, although we cannot exclude an alias that
would put the system close to a 4:3:2 configuration. We then analyse in more
detail the co-orbital scenario. We show that despite the influence of an inner
planet just outside the 2:3 mean-motion resonance, this potential co-orbital
system can be stable on the Giga-year time-scale for a variety of planetary
masses, either on a trojan or a horseshoe orbit. We predict that large TTVs
should arise in such configuration with a period of several hundred days. We
then show how the mass of each planet can be retrieved from these TTVs.
[26]
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.
[27]
oai:arXiv.org:1807.10652 [pdf] - 1767545
ACCESS: A featureless optical transmission spectrum for WASP-19b from
Magellan/IMACS
Espinoza, Néstor;
Rackham, Benjamin V. \;
Jordán, Andrés;
Apai, Dániel;
López-Morales, Mercedes;
Osip, David J.;
Grimm, Simon L.;
Hoeijmakers, Jens;
Wilson, Paul A.;
Bixel, Alex;
McGruder, Chima;
Rodler, Florian;
Weaver, Ian;
Lewis, Nikole K.;
Fortney, Jonathan J.;
Fraine, Jonathan
Submitted: 2018-07-27, last modified: 2018-10-02
The short period ($0.94$-day) transiting exoplanet WASP-19b is an exceptional
target for transmission spectroscopy studies, due to its relatively large
atmospheric scale-height ($\sim 500$ km) and equilibrium temperature ($\sim
2100$ K). Here we report on six precise spectroscopic Magellan/IMACS
observations, five of which target the full optical window from $0.45-0.9\mu$m
and one targeting the $0.4-0.55\mu$m blue-optical range. Five of these datasets
are consistent with a transmission spectrum without any significant spectral
features, while one shows a significant slope as a function of wavelength,
which we interpret as arising from photospheric heterogeneities in the star.
Coupled with HST/WFC3 infrared observations, our optical/near-infrared
measurements point to the presence of high altitude clouds in WASP-19b's
atmosphere in agreement with previous studies. Using a semi-analytical
retrieval approach, considering both planetary and stellar spectral features,
we find a water abundance consistent with solar for WASP-19b and strong
evidence for sub-solar abundances for optical absorbers such as TiO and Na; no
strong optical slope is detected, which suggests that if hazes are present,
they are much weaker than previously suggested. In addition, two spot-crossing
events are observed in our datasets and analyzed, including one of the first
unambiguously detected bright spot-crossing events on an exoplanet host star.
[28]
oai:arXiv.org:1808.05653 [pdf] - 1743870
Atomic iron and titanium in the atmosphere of the exoplanet KELT-9b
Hoeijmakers, H. Jens;
Ehrenreich, David;
Heng, Kevin;
Kitzmann, Daniel;
Grimm, Simon L.;
Allart, Romain;
Deitrick, Russell;
Wyttenbach, Aurelien;
Oreshenko, Maria;
Pino, Lorenzo;
Rimmer, Paul B.;
Molinari, Emilio;
Di Fabrizio, Luca
Submitted: 2018-08-16
The chemical composition of an exoplanet is a key ingredient in constraining
its formation history. Iron is the most abundant transition metal, but has
never been directly detected in an exoplanet due to its highly refractory
nature. KELT-9b (HD 195689b) is the archetype of the class of ultra-hot
Jupiters that straddle the transition between stars and gas-giant exoplanets
and serve as distinctive laboratories for studying atmospheric chemistry,
because of its high equilibrium temperature of 4050 +/- 180 K. These properties
imply that its atmosphere is a tightly constrained chemical system that is
expected to be nearly in chemical equilibrium and cloud-free. It was previously
predicted that the spectral lines of iron will be detectable in the visible
range of wavelengths. At these high temperatures, iron and several other
transition metals are not sequestered in molecules or cloud particles and exist
solely in their atomic forms. Here, we report the direct detection of atomic
neutral and singly-ionized iron (Fe and Fe+), and singly-ionized titanium (Ti+)
in KELT-9b via the cross-correlation technique applied to high-resolution
spectra obtained during the primary transit of the exoplanet.
[29]
oai:arXiv.org:1808.01803 [pdf] - 1763035
Interior characterization in multiplanetary systems: TRAPPIST-1
Submitted: 2018-08-06
Interior characterization traditionally relies on individual planetary
properties, ignoring correlations between different planets of the same system.
For multi-planetary systems, planetary data are generally correlated. This is
because, the differential masses and radii are better constrained than absolute
planetary masses and radii. We explore such correlations and data specific to
the multiplanetary-system of TRAPPIST-1 and study their value for our
understanding of planet interiors. Furthermore, we demonstrate that the rocky
interior of planets in a multi-planetary system can be preferentially probed by
studying the most dense planet representing a rocky interior analogue. Our
methodology includes a Bayesian inference analysis that uses a Markov chain
Monte Carlo scheme. Our interior estimates account for the anticipated
variability in the compositions and layer thicknesses of core, mantle, water
oceans and ice layers, and a gas envelope. Our results show that (1) interior
estimates significantly depend on available abundance proxies and (2) that the
importance of inter-dependent planetary data for interior characterization is
comparable to changes in data precision by 30 %. For the interiors of
TRAPPIST-1 planets, we find that possible water mass fractions generally range
from 0-25 %. The lack of a clear trend of water budgets with orbital period or
planet mass challenges possible formation scenarios. While our estimates change
relatively little with data precision, they critically depend on data accuracy.
If planetary masses varied within ~24 %, interiors would be consistent with
uniform (~7 %) or an increasing water mass fractions with orbital period (~2-12
%).
[30]
oai:arXiv.org:1711.08492 [pdf] - 1724847
Towards Consistent Modeling of Atmospheric Chemistry and Dynamics in
Exoplanets: Validation and Generalization of Chemical Relaxation Method
Submitted: 2017-11-22, last modified: 2018-05-29
Motivated by the work of Cooper & Showman, we revisit the chemical relaxation
method, which seeks to enhance the computational efficiency of
chemical-kinetics calculations by replacing the chemical network with a handful
of independent source/sink terms. Chemical relaxation solves the evolution of
the system and can treat disequilibrium chemistry, as the source/sink terms are
driven towards chemical equilibrium on a prescribed chemical timescale, but it
has surprisingly never been validated. First, we generalize the treatment by
forgoing the use of a single chemical timescale, instead developing a pathway
analysis tool that allows us to identify the rate-limiting reaction as a
function of temperature and pressure. For the interconversion between methane
and carbon monoxide and between ammonia, and molecular nitrogen, we identify
the key rate-limiting reactions for conditions relevant to currently
characterizable exo-atmospheres (500-3000 K, 0.1 mbar to 1 kbar). Second, we
extend chemical relaxation to include carbon dioxide and water. Third, we
examine the role of metallicity and carbon-to-oxygen ratio in chemical
relaxation. Fourth, we apply our pathway analysis tool to diagnose the
differences between our chemical network and that of Moses and Venot. Finally,
we validate the chemical relaxation method against full chemical kinetics
calculations in one dimension. For WASP-18b-, HD 189733b- and GJ 1214-b-like
atmospheres, we show that chemical relaxation is mostly accurate to within an
order of magnitude, a factor of 2 and $\sim 10\%$, respectively. The level of
accuracy attained allows for the chemical relaxation method to be included in
three-dimensional general circulation models.
[31]
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.
[32]
oai:arXiv.org:1802.02250 [pdf] - 1640520
Atmospheric reconnaissance of the habitable-zone Earth-sized planets
orbiting TRAPPIST-1
de Wit, J.;
Wakeford, H. R.;
Lewis, N.;
Delrez, L.;
Gillon, M.;
Selsis, F.;
Leconte, J.;
Demory, B. -O.;
Bolmont, E.;
Bourrier, V.;
Burgasser, A. J.;
Grimm, S.;
Jehin, E.;
Lederer, S. M.;
Owen, J.;
Stamenkovic, V.;
Triaud, A. H. M. J.
Submitted: 2018-02-06
Seven temperate Earth-sized exoplanets readily amenable for atmospheric
studies transit the nearby ultracool dwarf star TRAPPIST-1 (refs 1,2). Their
atmospheric regime is unknown and could range from extended primordial
hydrogen-dominated to depleted atmospheres (refs 3-6). Hydrogen in particular
is a powerful greenhouse gas that may prevent the habitability of inner planets
while enabling the habitability of outer ones (refs 6-8). An atmosphere largely
dominated by hydrogen, if cloud-free, should yield prominent spectroscopic
signatures in the near-infrared detectable during transits. Observations of the
innermost planets have ruled out such signatures (ref 9). However, the
outermost planets are more likely to have sustained such a Neptune-like
atmosphere (refs 10,11). Here, we report observations for the four planets
within or near the system's habitable zone, the circumstellar region where
liquid water could exist on a planetary surface (refs 12-14). These planets do
not exhibit prominent spectroscopic signatures at near-infrared wavelengths
either, which rules out cloud-free hydrogen-dominated atmospheres for
TRAPPIST-1 d, e and f, with significance of 8, 6 and 4 sigma, respectively.
Such an atmosphere is instead not excluded for planet g. As high-altitude
clouds and hazes are not expected in hydrogen-dominated atmospheres around
planets with such insolation (refs 15,16), these observations further support
their terrestrial and potentially habitable nature.
[33]
oai:arXiv.org:1802.01377 [pdf] - 1629579
The nature of the TRAPPIST-1 exoplanets
Grimm, Simon L.;
Demory, Brice-Olivier;
Gillon, Michaël;
Dorn, Caroline;
Agol, Eric;
Burdanov, Artem;
Delrez, Laetitia;
Sestovic, Marko;
Triaud, Amaury H. M. J.;
Turbet, Martin;
Bolmont, Émeline;
Caldas, Anthony;
de Wit, Julien;
Jehin, Emmanuël;
Leconte, Jérémy;
Raymond, Sean N.;
Van Grootel, Valérie;
Burgasser, Adam J.;
Carey, Sean;
Fabrycky, Daniel;
Heng, Kevin;
Hernandez, David M.;
Ingalls, James G.;
Lederer, Susan;
Selsis, Franck;
Queloz, Didier
Submitted: 2018-02-05
Context. The TRAPPIST-1 system hosts seven Earth-sized, temperate exoplanets
orbiting an ultra-cool dwarf star. As such, it represents a remarkable setting
to study the formation and evolution of terrestrial planets that formed in the
same protoplanetary disk. While the sizes of the TRAPPIST-1 planets are all
known to better than 5% precision, their densities have significant
uncertainties (between 28% and 95%) because of poor constraints on the planet's
masses. Aims.The goal of this paper is to improve our knowledge of the
TRAPPIST-1 planetary masses and densities using transit-timing variations
(TTV). The complexity of the TTV inversion problem is known to be particularly
acute in multi-planetary systems (convergence issues, degeneracies and size of
the parameter space), especially for resonant chain systems such as TRAPPIST-1.
Methods. To overcome these challenges, we have used a novel method that employs
a genetic algorithm coupled to a full N-body integrator that we applied to a
set of 284 individual transit timings. This approach enables us to efficiently
explore the parameter space and to derive reliable masses and densities from
TTVs for all seven planets. Results. Our new masses result in a five- to
eight-fold improvement on the planetary density uncertainties, with precisions
ranging from 5% to 12%. These updated values provide new insights into the bulk
structure of the TRAPPIST-1 planets. We find that TRAPPIST-1\,c and e likely
have largely rocky interiors, while planets b, d, f, g, and h require envelopes
of volatiles in the form of thick atmospheres, oceans, or ice, in most cases
with water mass fractions less than 5%.
[34]
oai:arXiv.org:1508.00917 [pdf] - 1579678
Stochasticity & Predictability in Terrestrial Planet Formation
Submitted: 2015-08-04, last modified: 2017-09-04
Terrestrial planets are thought to be the result of a vast number of
gravitational interactions and collisions between smaller bodies. We use
numerical simulations to show that practically identical initial conditions
result in a wide array of final planetary configurations. This is a result of
the chaotic evolution of trajectories which are highly sensitive to minuscule
displacements. We determine that differences between systems evolved from
virtually identical initial conditions can be larger than the differences
between systems evolved from very different initial conditions. This implies
that individual simulations lack predictive power. For example, there is not a
reproducible mapping between the initial and final surface density profiles.
However, some key global properties can still be extracted if the statistical
spread across many simulations is considered. Based on these spreads, we
explore the collisional growth and orbital properties of terrestrial planets
which assemble from different initial conditions (we vary the initial
planetesimal distribution, planetesimal masses, and giant planet orbits).
Confirming past work, we find that the resulting planetary systems are sculpted
by sweeping secular resonances. Configurations with giant planets on eccentric
orbits produce fewer and more massive terrestrial planets on tighter orbits
than those with giants on circular orbits. This is further enhanced if the
initial mass distribution is biased to the inner regions. In all cases, the
outer edge of the system is set by the final location of the $\nu_6$ resonance
and we find that the mass distribution peaks at the $\nu_5$ resonance. Using
existing observations, we find that extrasolar systems follow similar trends.
Although differences between our numerical modelling and exoplanetary systems
remain, we suggest that CoRoT-7, HD 20003, and HD 20781 may host undetected
giant planets.
[35]
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.
[36]
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).
[37]
oai:arXiv.org:1705.10362 [pdf] - 1584013
Habitable Moist Atmospheres On Terrestrial Planets Near the Inner Edge
Of the Habitable Zone Around M-dwarfs
Submitted: 2017-05-29, last modified: 2017-07-04
Terrestrial planets in the habitable zones (HZs) of low-mass stars and cool
dwarfs have received significant scrutiny recently because their shorter
orbital periods increase their chances of detection and characterization
compared to planets around G-dwarfs. As these planets are likely tidal-locked,
improved 3D numerical simulations of such planetary atmospheres are needed to
guide target selection. Here we use a 3-D climate system model, updated with
new water-vapor absorption coefficients derived from the HITRAN 2012 database,
to study ocean covered planets at the inner edge of the HZ around late-M to
mid-K stars ($2600$ K <= Teff <= 4500K). Our results indicate that these
updated water-vapor coefficients result in significant warming compared to
previous studies, so the inner HZ around M-dwarfs is not as close as suggested
by earlier work. Assuming synchronously rotating planets with background 1 bar
N2 atmospheres, we find that planets at the inner HZ of stars with Teff > 3000K
undergo the classical "moist-greenhouse" (H2O mixing ratio > 10-3 in the
stratosphere) at significantly lower surface temperature (~ 280K) in our 3-D
model compared with 1-D climate models (~ 340K). This implies that some planets
around low mass stars can simultaneously undergo water-loss and remain
habitable. However, for star with Teff <= 3000K, planets at the inner HZ may
directly transition to a runaway state, while bypassing the moist greenhouse
water-loss entirely. We analyze transmission spectra of planets in a moist
green-house regime, and find that there are several prominent H2O features
within JWST instruments range. Thus, relying only upon standard Earth-analog
spectra with 24-hour rotation period around M-dwarfs for habitability studies
will miss the strong H2O features that one would expect to see on synchronously
rotating planets around M-dwarf stars, with JWST.
[38]
oai:arXiv.org:1703.04166 [pdf] - 1581897
A seven-planet resonant chain in TRAPPIST-1
Luger, Rodrigo;
Sestovic, Marko;
Kruse, Ethan;
Grimm, Simon L.;
Demory, Brice-Olivier;
Agol, Eric;
Bolmont, Emeline;
Fabrycky, Daniel;
Fernandes, Catarina S.;
Van Grootel, Valérie;
Burgasser, Adam;
Gillon, Michaël;
Ingalls, James G.;
Jehin, Emmanuël;
Raymond, Sean N.;
Selsis, Franck;
Triaud, Amaury H. M. J.;
Barclay, Thomas;
Barentsen, Geert;
Howell, Steve B.;
Delrez, Laetitia;
de Wit, Julien;
Foreman-Mackey, Daniel;
Holdsworth, Daniel L.;
Leconte, Jérémy;
Lederer, Susan;
Turbet, Martin;
Almleaky, Yaseen;
Benkhaldoun, Zouhair;
Magain, Pierre;
Morris, Brett;
Heng, Kevin;
Queloz, Didier
Submitted: 2017-03-12, last modified: 2017-05-22
The TRAPPIST-1 system is the first transiting planet system found orbiting an
ultra-cool dwarf star. At least seven planets similar to Earth in radius and in
mass were previously found to transit this host star. Subsequently, TRAPPIST-1
was observed as part of the K2 mission and, with these new data, we report the
measurement of an 18.77 d orbital period for the outermost planet, TRAPPIST-1h,
which was unconstrained until now. This value matches our theoretical
expectations based on Laplace relations and places TRAPPIST-1h as the seventh
member of a complex chain, with three-body resonances linking every member. We
find that TRAPPIST-1h has a radius of 0.727 Earth radii and an equilibrium
temperature of 173 K. We have also measured the rotational period of the star
at 3.3 d and detected a number of flares consistent with a low-activity,
middle-aged, late M dwarf.
[39]
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).
[40]
oai:arXiv.org:1607.05535 [pdf] - 1490798
THOR: A New and Flexible Global Circulation Model to Explore Planetary
Atmospheres
Submitted: 2016-07-19
We have designed and developed, from scratch, a global circulation model
named THOR that solves the three-dimensional non-hydrostatic Euler equations.
Our general approach lifts the commonly used assumptions of a shallow
atmosphere and hydrostatic equilibrium. We solve the "pole problem" (where
converging meridians on a sphere lead to increasingly smaller time steps near
the poles) by implementing an icosahedral grid. Irregularities in the grid,
which lead to grid imprinting, are smoothed using the "spring dynamics"
technique. We validate our implementation of spring dynamics by examining
calculations of the divergence and gradient of test functions. To prevent the
computational time step from being bottlenecked by having to resolve sound
waves, we implement a split-explicit method together with a horizontally
explicit and vertically implicit integration. We validate our global
circulation model by reproducing the Earth and also the hot Jupiter-like
benchmark tests. THOR was designed to run on Graphics Processing Units (GPUs),
which allows for physics modules (radiative transfer, clouds, chemistry) to be
added in the future, and is part of the open-source Exoclimes Simulation
Platform (ESP; www.exoclime.org).
[41]
oai:arXiv.org:1607.02494 [pdf] - 1528105
Optimal non-circular fiber geometries for image scrambling in
high-resolution spectrographs
Submitted: 2016-07-08
Optical fibers are a key component for high-resolution spectrographs to
attain high precision in radial velocity measurements. We present a custom
fiber with a novel core geometry - a 'D'-shape. From a theoretical standpoint,
such a fiber should provide superior scrambling and modal noise mitigation,
since unlike the commonly used circular and polygonal fiber cross sections, it
shows chaotic scrambling. We report on the fabrication process of a test fiber
and compare the optical properties, scrambling performance, and modal noise
behavior of the D-fiber with those of common polygonal fibers.
[42]
oai:arXiv.org:1503.03806 [pdf] - 1263925
HELIOS-K: An Ultrafast, Open-source Opacity Calculator for Radiative
Transfer
Submitted: 2015-03-12, last modified: 2015-06-18
We present an ultrafast opacity calculator that we name HELIOS-K. It takes a
line list as an input, computes the shape of each spectral line and provides an
option for grouping an enormous number of lines into a manageable number of
bins. We implement a combination of Algorithm 916 and Gauss-Hermite quadrature
to compute the Voigt profile, write the code in CUDA and optimise the
computation for graphics processing units (GPUs). We restate the theory of the
k-distribution method and use it to reduce $\sim 10^5$ to $10^8$ lines to $\sim
10$ to $10^4$ wavenumber bins, which may then be used for radiative transfer,
atmospheric retrieval and general circulation models. The choice of line-wing
cutoff for the Voigt profile is a significant source of error and affects the
value of the computed flux by $\sim 10\%$. This is an outstanding physical
(rather than computational) problem, due to our incomplete knowledge of
pressure broadening of spectral lines in the far line wings. We emphasize that
this problem remains regardless of whether one performs line-by-line
calculations or uses the k-distribution method and affects all calculations of
exoplanetary atmospheres requiring the use of wavelength-dependent opacities.
We elucidate the correlated-k approximation and demonstrate that it applies
equally to inhomogeneous atmospheres with a single atomic/molecular species or
homogeneous atmospheres with multiple species. Using a NVIDIA K20 GPU, HELIOS-K
is capable of computing an opacity function with $\sim 10^5$ spectral lines in
$\sim 1$ second and is publicly available as part of the Exoclimes Simulation
Platform (ESP; www.exoclime.org).
[43]
oai:arXiv.org:1404.2324 [pdf] - 893230
The GENGA Code: Gravitational Encounters in N-body simulations with GPU
Acceleration
Submitted: 2014-04-08, last modified: 2014-09-20
We describe an open source GPU implementation of a hybrid symplectic N-body
integrator, GENGA (Gravitational ENcounters with Gpu Acceleration), designed to
integrate planet and planetesimal dynamics in the late stage of planet
formation and stability analyses of planetary systems. GENGA uses a hybrid
symplectic integrator to handle close encounters with very good energy
conservation, which is essential in long-term planetary system integration. We
extended the second order hybrid integration scheme to higher orders. The GENGA
code supports three simulation modes: Integration of up to 2048 massive bodies,
integration with up to a million test particles, or parallel integration of a
large number of individual planetary systems. We compare the results of GENGA
to Mercury and pkdgrav2 in respect of energy conservation and performance, and
find that the energy conservation of GENGA is comparable to Mercury and around
two orders of magnitude better than pkdgrav2. GENGA runs up to 30 times faster
than Mercury and up to eight times faster than pkdgrav2. GENGA is written in
CUDA C and runs on all NVIDIA GPUs with compute capability of at least 2.0.
[44]
oai:arXiv.org:1305.4070 [pdf] - 1171441
Super Earths and Dynamical Stability of Planetary Systems: First
Parallel GPU Simulations Using GENGA
Submitted: 2013-05-17
We report on the stability of hypothetical Super-Earths in the habitable zone
of known multi-planetary systems. Most of them have not yet been studied in
detail concerning the existence of additional low-mass planets. The new N-body
code GENGA developed at the UZH allows us to perform numerous N-body
simulations in parallel on GPUs. With this numerical tool, we can study the
stability of orbits of hypothetical planets in the semi-major axis and
eccentricity parameter space in high resolution. Massless test particle
simulations give good predictions on the extension of the stable region and
show that HIP 14180 and HD 37124 do not provide stable orbits in the habitable
zone. Based on these simulations, we carry out simulations of 10 Earth mass
planets in several systems (HD 11964, HD 47186, HD 147018, HD 163607, HD
168443, HD 187123, HD 190360, HD 217107 and HIP 57274). They provide more exact
information about orbits at the location of mean motion resonances and at the
edges of the stability zones. Beside the stability of orbits, we study the
secular evolution of the planets to constrain probable locations of
hypothetical planets. Assuming that planetary systems are in general closely
packed, we find that apart from HD 168443, all of the systems can harbor 10
Earth mass planets in the habitable zone.
