99 research outputs found
Modelling the local and global cloud formation on HD 189733b
Context. Observations suggest that exoplanets such as HD 189733b form clouds
in their atmospheres which have a strong feedback onto their thermodynamical
and chemical structure, and overall appearance. Aims. Inspired by mineral cloud
modelling efforts for Brown Dwarf atmospheres, we present the first spatially
varying kinetic cloud model structures for HD 189733b. Methods. We apply a
2-model approach using results from a 3D global radiation-hydrodynamic
simulation of the atmosphere as input for a detailed, kinetic cloud formation
model. Sampling the 3D global atmosphere structure with 1D trajectories allows
us to model the spatially varying cloud structure on HD 189733b. The resulting
cloud properties enable the calculation of the scattering and absorption
properties of the clouds. Results. We present local and global cloud structure
and property maps for HD 189733b. The calculated cloud properties show
variations in composition, size and number density of cloud particles which are
strongest between the dayside and nightside. Cloud particles are mainly
composed of a mix of materials with silicates being the main component. Cloud
properties, and hence the local gas composition, change dramatically where
temperature inversions occur locally. The cloud opacity is dominated by
absorption in the upper atmosphere and scattering at higher pressures in the
model. The calculated 8{\mu}m single scattering Albedo of the cloud particles
are consistent with Spitzer bright regions. The cloud particles scattering
properties suggest that they would sparkle/reflect a midnight blue colour at
optical wavelengths.Comment: Accepted for publication (A&A) - 21/05/2015 (Low Resolution Maps
Critical velocity ionisation in substellar atmospheres
The observation of radio, X-ray and Hα emission from substellar objects indicates the presence of plasma regions and associated high-energy processes in their surrounding envelopes. This paper numerically simulates and characterises Critical Velocity Ionisation, a potential ionisation process, that can efficiently generate plasma as a result of neutral gas flows interacting with seed magnetized plasmas. By coupling a Gas-MHD interactions code (to simulate the ionisation mechanism) with a substellar global circulation model (to provide the required gas flows) we quantify the spatial extent of the resulting plasma regions, their degree of ionisation and their lifetime for a typical substellar atmosphere. It is found that the typical average ionisation fraction reached at equilibrium (where the ionisation and recombination rates are equal and opposite) ranges from 10-5 to 10-8, at pressures between 10-1 and 10-3 bar, with a trend of increasing ionisation fraction with decreasing atmospheric pressure. The ionisation fractions reached as a result of Critical Velocity Ionisation are sufficient to allow magnetic fields to couple to gas flows in the atmosphere
Dynamic mineral clouds on HD 189733b : I. 3D RHD with kinetic, non-equilibrium cloud formation
Funding: European community under the FP7 ERC starting grant 257431 (G.L., Ch. H.).Context. Observations of exoplanet atmospheres have revealed the presence of cloud particles in their atmospheres. 3D modelling of cloud formation in atmospheres of extrasolar planets coupled to the atmospheric dynamics has long been a challenge. Aims. We investigate the thermo-hydrodynamic properties of cloud formation processes in the atmospheres of hot Jupiter exoplanets. Methods. We simulate the dynamic atmosphere of HD 189733b with a 3D model that couples 3D radiative-hydrodynamics with a kinetic, microphysical mineral cloud formation module designed for RHD/GCM exoplanet atmosphere simulations. Our simulation includes the feedback effects of cloud advection and settling, gas phase element advection and depletion/replenishment and the radiative effects of cloud opacity. We model the cloud particles as a mix of mineral materials which change in size and composition as they travel through atmospheric thermo-chemical environments. All local cloud properties such as number density, grain size and material composition are time-dependently calculated. Gas phase element depletion as a result of cloud formation is included in the model. In situ effective medium theory and Mie theory is applied to calculate the wavelength dependent opacity of the cloud component. Results. We present a 3D cloud structure of a chemically complex, gaseous atmosphere of the hot Jupiter HD 189733b. Mean cloud particle sizes are typically sub-micron (0.01−0.5 μm) at pressures less than 1 bar with hotter equatorial regions containing the smallest grains. Denser cloud structures occur near terminator regions and deeper (~1 bar) atmospheric layers. Silicate materials such as MgSiO3[s] are found to be abundant at mid-high latitudes, while TiO2[s] and SiO2[s] dominate the equatorial regions. Elements involved in the cloud formation can be depleted by several orders of magnitude. Conclusions. The interplay between radiative-hydrodynamics and cloud kinetics leads to an inhomogeneous, wavelength dependent opacity cloud structure with properties differing in longitude, latitude and depth. This suggests that transit spectroscopy would sample a variety of cloud particles properties (sizes, composition, densities).Publisher PDFPeer reviewe
Dynamic mineral clouds on HD 189733b. II. Monte Carlo radiative transfer for 3D cloudy exoplanet atmospheres : combining scattering and emission spectra
G.L. and Ch.H. highlight the financial support of the European community under the FP7 ERC starting grant 257431.Context. As the 3D spatial properties of exoplanet atmospheres are being observed in increasing detail by current and new generations of telescopes, the modelling of the 3D scattering effects of cloud forming atmospheres with inhomogeneous opacity structures becomes increasingly important to interpret observational data. Aims. We model the scattering and emission properties of a simulated cloud forming, inhomogeneous opacity, hot Jupiter atmosphere of HD 189733b. We compare our results to available Hubble Space Telescope (HST) and Spitzer data and quantify the effects of 3D multiple scattering on observable properties of the atmosphere. We discuss potential observational properties of HD 189733b for the upcoming Transiting Exoplanet Survey Satellite (TESS) and CHaracterising ExOPlanet Satellite (CHEOPS) missions. Methods. We developed a Monte Carlo radiative transfer code and applied it to post-process output of our 3D radiative-hydrodynamic, cloud formation simulation of HD 189733b. We employed three variance reduction techniques, i.e. next event estimation, survival biasing, and composite emission biasing, to improve signal to noise of the output. For cloud particle scattering events, we constructed a log-normal area distribution from the 3D cloud formation radiative-hydrodynamic results, which is stochastically sampled in order to model the Rayleigh and Mie scattering behaviour of a mixture of grain sizes. Results. Stellar photon packets incident on the eastern dayside hemisphere show predominantly Rayleigh, single-scattering behaviour, while multiple scattering occurs on the western hemisphere. Combined scattered and thermal emitted light predictions are consistent with published HST and Spitzer secondary transit observations. Our model predictions are also consistent with geometric albedo constraints from optical wavelength ground-based polarimetry and HST B band measurements. We predict an apparent geometric albedo for HD 189733b of 0.205 and 0.229, in the TESS and CHEOPS photometric bands respectively. Conclusions. Modelling the 3D geometric scattering effects of clouds on observables of exoplanet atmospheres provides an important contribution to the attempt to determine the cloud properties of these objects. Comparisons between TESS and CHEOPS photometry may provide qualitative information on the cloud properties of nearby hot Jupiter exoplanets.Publisher PDFPeer reviewe
The mineral clouds on HD209 458b and HD189 733b
This is the final version of the article. Available from the publisher via the DOI in this record.3D atmosphere model results are used to comparatively study the kinetic, nonequilibrium
cloud formation in the atmospheres of two example planets guided by
the giant gas planets HD 209 458b and HD 189 733b. Rather independently of hydrodynamic
model differences, our cloud modelling suggests that both planets are covered
in mineral clouds throughout the entire modelling domain. Both planets harbour chemically
complex clouds that are made of mineral particles that have a height-dependent
material composition and size. The remaining gas-phase element abundances strongly
effects the molecular abundances of the atmosphere in the cloud forming regions. Hydrocarbon
and cyanopolyyne molecules can be rather abundant in the inner, dense
part of the atmospheres of HD 189 733b and HD 209 458b. No one value for metallicity
and the C/O ratio can be used to describe an extrasolar planet. Our results concerning
the presence and location of water in relation to the clouds explain some of the
observed differences between the two planets. In HD 189 733b, strong water features
have been reported while such features are less strong for HD 209 458b. By considering
the location of the clouds in the two atmospheres, we see that obscuring clouds exist
high in the atmosphere of HD 209 458b, but much deeper in HD 189 733b. We further
conclude that the (self-imposed) degeneracy of cloud parameters in retrieval methods
can only be lifted if the cloud formation processes are accurately modelled in contrast
to prescribing them by independent parametersWe highlight financial support of the European Community
under the FP7 by the ERC starting grant 257431 and by an
ERC advanced grant 247060. JK acknowledges the Rosen
fellowship from the Brooklyn College New York, US. Some
of the calculations for this paper were performed on the
DIRAC Facility jointly funded by STFC, the Large Facilities
Capital Fund of BIS, and the University of Exeter
Mapping Exoplanets
The varied surfaces and atmospheres of planets make them interesting places
to live, explore, and study from afar. Unfortunately, the great distance to
exoplanets makes it impossible to resolve their disk with current or near-term
technology. It is still possible, however, to deduce spatial inhomogeneities in
exoplanets provided that different regions are visible at different
times---this can be due to rotation, orbital motion, and occultations by a
star, planet, or moon. Astronomers have so far constructed maps of thermal
emission and albedo for short period giant planets. These maps constrain
atmospheric dynamics and cloud patterns in exotic atmospheres. In the future,
exo-cartography could yield surface maps of terrestrial planets, hinting at the
geophysical and geochemical processes that shape them.Comment: Updated chapter for Handbook of Exoplanets, eds. Deeg & Belmonte. 17
pages, including 6 figures and 4 pages of reference
Is tidal heating sufficient to explain bloated exoplanets? Consistent calculations accounting for finite initial eccentricity
In this paper, we present the consistent evolution of short-period exoplanets
coupling the tidal and gravothermal evolution of the planet. Contrarily to
previous similar studies, our calculations are based on the complete tidal
evolution equations of the Hut model, valid at any order in eccentricity,
obliquity and spin. We demonstrate, both analytically and numerically, that,
except if the system was formed with a nearly circular orbit (e<0.2), solving
consistently the complete tidal equations is mandatory to derive correct tidal
evolution histories. We show that calculations based on tidal models truncated
at second order in eccentricity, as done in all previous studies, lead to
erroneous tidal evolutions. As a consequence, tidal energy dissipation rates
are severely underestimated in all these calculations and the characteristic
timescales for the various orbital parameters evolutions can be wrong by up to
three orders in magnitude.
Based on these complete, consistent calculations, we revisit the viability of
the tidal heating hypothesis to explain the anomalously large radius of
transiting giant planets. We show that, even though tidal dissipation does
provide a substantial contribution to the planet's heat budget and can explain
some of the moderately bloated hot-Jupiters, this mechanism can not explain
alone the properties of the most inflated objects, including HD 209458b.
Indeed, solving the complete tidal equations shows that enhanced tidal
dissipation and thus orbit circularization occur too early during the planet's
evolution to provide enough extra energy at the present epoch. In that case
another mechanisms, such as stellar irradiation induced surface winds
dissipating in the planet's tidal bulges, or inefficient convection in the
planet's interior must be invoked, together with tidal dissipation, to provide
all the pieces of the abnormally large exoplanet puzzle.Comment: 14 pages, 10 figures, Accepted for publication in Astronomy and
Astrophysics
Tidal friction in close-in satellites and exoplanets. The Darwin theory re-visited
This report is a review of Darwin's classical theory of bodily tides in which
we present the analytical expressions for the orbital and rotational evolution
of the bodies and for the energy dissipation rates due to their tidal
interaction. General formulas are given which do not depend on any assumption
linking the tidal lags to the frequencies of the corresponding tidal waves
(except that equal frequency harmonics are assumed to span equal lags).
Emphasis is given to the cases of companions having reached one of the two
possible final states: (1) the super-synchronous stationary rotation resulting
from the vanishing of the average tidal torque; (2) the capture into a 1:1
spin-orbit resonance (true synchronization). In these cases, the energy
dissipation is controlled by the tidal harmonic with period equal to the
orbital period (instead of the semi-diurnal tide) and the singularity due to
the vanishing of the geometric phase lag does not exist. It is also shown that
the true synchronization with non-zero eccentricity is only possible if an
extra torque exists opposite to the tidal torque. The theory is developed
assuming that this additional torque is produced by an equatorial permanent
asymmetry in the companion. The results are model-dependent and the theory is
developed only to the second degree in eccentricity and inclination
(obliquity). It can easily be extended to higher orders, but formal accuracy
will not be a real improvement as long as the physics of the processes leading
to tidal lags is not better known.Comment: 30 pages, 7 figures, corrected typo
- …