-
To understand and compare the 3D atmospheric structure of HD 209458 b and HD
189733 b, focusing on the formation and distribution of cloud particles, as
well as their feedback on the dynamics and thermal profile. We couple the 3D
Met Office Unified Model (UM), including detailed treatments of atmospheric
radiative transfer and dynamics, to a kinetic cloud formation scheme. The
resulting model self--consistently solves for the formation of condensation
seeds, surface growth and evaporation, gravitational settling and advection,
cloud radiative feedback via absorption and, crucially, scattering. Fluxes
directly obtained from the UM are used to produce synthetic SEDs and phase
curves. Our simulations show extensive cloud formation in both planets,
however, cooler temperatures in the HD 189733 b result in higher cloud particle
number densities. Sub-micron particles are suspended by vertical flows leading
to extensive upper-atmosphere cloud cover. A combination of meridional
advection and efficient cloud formation in cooler high latitude regions, result
in enhanced cloud coverage for latitudes > 30 degrees and leads to a zonally
banded structure for all our simulations. The cloud bands extend around the
entire planet(s), as the temperatures, even on the day side, remain below the
condensation temperature of silicates and oxides. Therefore, our simulated
optical phase curve for HD 209458 b shows no `offset', in contrast to
observations. Efficient scattering by cloud results an atmospheric cooling of
up to 250K, and an advection-driven fluctuating cloud opacity causes temporal
variability in the thermal emission. The inclusion of this fundamental
cloud-atmosphere radiative feedback leads to significant differences with
approaches neglecting these physical elements and suggests both a note of
caution of interpretations neglecting such cloud feedback and scattering, and
merits further study.
-
Aims. We investigate the feasibility of planetesimal growth in circumbinary
protoplanetary disks around the observed systems Kepler- 16 and Kepler-34 under
the gravitational influence of a precessing eccentric gas disk. Methods. We
embed the results of our previous hydrodynamical simulations of protoplanetary
disks around binaries into an N-body code to perform 3D, high-resolution,
inter-particle gravity-enabled simulations of planetesimal growth and dynamics
that include the gravitational force imparted by the gas. Results. Including
the full, precessing asymmetric gas disk generates high eccentricity orbits for
planetesimals orbiting at the edge of the circumbinary cavity, where the gas
surface density and eccentricity have their largest values. The gas disk is
able to efficiently align planetesimal pericenters in some regions leading to
phased, non-interacting orbits. Outside of these areas eccentric planetesimal
orbits become misaligned and overlap leading to crossing orbits and high
relative velocities during planetesimal collisions. This can lead to an
increase in the number of erosive collisions that far outweighs the number of
collisions that result in growth. Gravitational focusing from the static
axisymmetric gas disk is weak and does not significantly alter collision
outcomes from the gas free case. Conclusions. Due to asymmetries in the gas
disk, planetesimals are strongly perturbed onto highly eccentric orbits. Where
planetesimals orbits are not well aligned, orbit crossings lead to an increase
in the number of erosive collisions. This makes it difficult for sustained
planetesimal accretion to occur at the location of Kepler-16b and Kepler-34b
and we therefore rule out in-situ growth. This adds further support to our
initial suggestions that most circumbinary planets should form further out in
the disk and migrate inwards.
-
The Kepler mission's discovery of a number of circumbinary planets orbiting
close (a_p < 1.1 au) to the stellar binary raises questions as to how these
planets could have formed given the intense gravitational perturbations the
dual stars impart on the disk. The gas component of circumbinary protoplanetary
disks is perturbed in a similar manner to the solid, planetesimal dominated
counterpart, although the mechanism by which disk eccentricity originates
differs. This is the first work of a series that aims to investigate the
conditions for planet formation in circumbinary protoplanetary disks. We
present a number of hydrodynamical simulations that explore the response of gas
disks around two observed binary systems: Kepler-16 and Kepler-34. We probe the
importance of disk viscosity, aspect-ratio, inner boundary condition, initial
surface density gradient, and self-gravity on the dynamical evolution of the
disk, as well as its quasi steady-state profile. We find there is a strong
influence of binary type on the mean disk eccentricity, e_d, leading to e_d =
0.02 - 0.08 for Kepler-16 and e_d = 0.10 - 0.15 in Kepler-34. The value of
alpha-viscosity has little influence on the disk, but we find a strong increase
in mean disk eccentricity with increasing aspect-ratio due to wave propagation
effects. The choice of inner boundary condition only has a small effect on the
surface density and eccentricity of the disk. Our primary finding is that
including disk self-gravity has little impact on the evolution or final state
of the disk for disks with masses less than 12.5 times that of the minimum-mass
solar nebula. This finding contrasts with the results of self-gravity relevance
in circumprimary disks, where its inclusion is found to be an important factor
in describing the disk evolution.