Critical core mass for enriched envelopes: the role of H2O condensation

Context. Within the core accretion scenario of planetary formation, most simulations performed so far always assume the accreting envelope to have a solar composition. From the study of meteorite showers on Earth and numerical simulations, we know that planetesimals must undergo thermal ablation and disruption when crossing a protoplanetary envelope. Once the protoplanet has acquired an atmosphere, the primordial envelope gets enriched in volatiles and silicates from the planetesimals. This change of envelope composition during the formation can have a significant effect in the final atmospheric composition and on the formation timescale of giant planets. Aims. To investigate the physical implications of considering the envelope enrichment of protoplanets due to the disruption of icy planetesimals during their way to the core. Particular focus is placed on the effect on the critical core mass for envelopes where condensation of water can occur. Methods. Internal structure models are numerically solved with the implementation of updated opacities for all ranges of metallicities and the software CEA to compute the equation of state. CEA computes the chemical equilibrium for an arbitrary mixture of gases and allows the condensation of some species, including water. This means that the latent heat of phase transitions is consistently incorporated in the total energy budget. Results. The critical core mass is found to decrease significantly when an enriched envelope composition is considered in the internal structure equations. A particular strong reduction of the critical core mass is obtained for planets whose envelope metallicity is larger than Z=0.45 when the outer boundary conditions are suitable for condensation of water to occur in the top layers of the atmosphere. We show that this effect is qualitatively preserved when the atmosphere is out of chemical equilibrium.Comment: Accepted for publication in A&

Elemental abundances and minimum mass of heavy elements in the envelope of HD 189733b

Oxygen (O) and carbon (C) have been inferred recently to be subsolar in abundance from spectra of the atmosphere of the transiting hot Jupiter HD 189733b. Yet, the mass and radius of the planet coupled with structure models indicate a strongly supersolar abundance of heavy elements in the interior of this object. Here we explore the discrepancy between the large amount of heavy elements suspected in the planet's interior and the paucity of volatiles measured in its atmosphere. We describe the formation sequence of the icy planetesimals formed beyond the snow line of the protoplanetary disk and calculate the composition of ices ultimately accreted in the envelope of HD 189733b on its migration pathway. This allows us to reproduce the observed volatile abundances by adjusting the mass of ices vaporized in the envelope. The predicted elemental mixing ratios should be 0.15--0.3 times solar in the envelope of HD 189733b if they are fitted to the recent O and C determinations. However, our fit to the minimum mass of heavy elements predicted by internal structure models gives elemental abundances that are 1.2--2.4 times oversolar in the envelope of HD189733b. We propose that the most likely cause of this discrepancy is irradiation from the central star leading to development of a radiative zone in the planet's outer envelope which would induce gravitational settling of elements. Hence, all strongly irradiated extrasolar planets should present subsolar abundances of volatiles. We finally predict that the abundances of nitrogen (N), sulfur (S) and phosphorus (P) are of $\sim$ $2.8 \times 10^{-5}$, $5.3 \times 10^{-6}$ and $1.8 \times 10^{-7}$ relative to H$_2$, respectively in the atmosphere of HD 189733b.Comment: Accepted for publication in Astronomy & Astrophysic

A generalized bayesian inference method for constraining the interiors of super Earths and sub-Neptunes

We aim to present a generalized Bayesian inference method for constraining interiors of super Earths and sub-Neptunes. Our methodology succeeds in quantifying the degeneracy and correlation of structural parameters for high dimensional parameter spaces. Specifically, we identify what constraints can be placed on composition and thickness of core, mantle, ice, ocean, and atmospheric layers given observations of mass, radius, and bulk refractory abundance constraints (Fe, Mg, Si) from observations of the host star's photospheric composition. We employed a full probabilistic Bayesian inference analysis that formally accounts for observational and model uncertainties. Using a Markov chain Monte Carlo technique, we computed joint and marginal posterior probability distributions for all structural parameters of interest. We included state-of-the-art structural models based on self-consistent thermodynamics of core, mantle, high-pressure ice, and liquid water. Furthermore, we tested and compared two different atmospheric models that are tailored for modeling thick and thin atmospheres, respectively. First, we validate our method against Neptune. Second, we apply it to synthetic exoplanets of fixed mass and determine the effect on interior structure and composition when (1) radius, (2) atmospheric model, (3) data uncertainties, (4) semi-major axes, (5) atmospheric composition (i.e., a priori assumption of enriched envelopes versus pure H/He envelopes), and (6) prior distributions are varied. Our main conclusions are: [...]Comment: Astronomy & Astrophysics, 597, A37, 17 pages, 11 figure

Planetary mass-radius relations across the galaxy

Planet formation theory suggests that planet bulk compositions are likely to reflect the chemical abundance ratios of their host star's photosphere. Variations in the abundance of particular chemical species in stellar photospheres between different galactic stellar populations demonstrate that there are differences among the expected solid planet bulk compositions. We aim to present planetary mass-radius relations of solid planets for kinematically differentiated stellar populations, namely, the thin disc, thick disc, and halo. Using two separate internal structure models, we generated synthetic planets using bulk composition inputs derived from stellar abundances. We explored two scenarios, specifically iron-silicate planets at 0.1 AU and silicate-iron-water planets at 4 AU. We show that there is a persistent statistical difference in the expected mass-radius relations of solid planets among the different galactic stellar populations. At 0.1 AU for silicate-iron planets, there is a 1.51 to 2.04\% mean planetary radius difference between the thick and thin disc stellar populations, whilst for silicate-iron-water planets past the ice line at 4 AU, we calculate a 2.93 to 3.26\% difference depending on the models. Between the halo and thick disc, we retrieve at 0.1 AU a 0.53 to 0.69\% mean planetary radius difference, and at 4 AU we find a 1.24 to 1.49\% difference depending on the model. Future telescopes (such as PLATO) will be able to precisely characterize solid exoplanets and demonstrate the possible existence of planetary mass-radius relationship variability between galactic stellar populations.Comment: 11 pages, 9 figures, accepted for publication in Astronomy & Astrophysic

Non-uniform integrability and generalized Young measures,

Given a bounded sequence (u n ) in L 1 (Ω, µ; IR d ), we describe the weak limits in the sense of measures of f (x, u n ) µ for a class of continuous integrands with linear growth at infinity. The defect of uniform integrability of the sequence f (x, u n ) is described by a measure m and a family of probability measures on S d−1 whereas the classical Young measure is associated with the biting limits in the sense of Chacon's lemma. Some consequences of this new approach are given in Calculus of Variations

Beat Cepheids as Probes of Stellar and Galactic Metallicity

The mere location of a Beat Cepheid model in a Period Ratio vs. Period diagram (Petersen diagram) puts very tight constraints on its metallicity Z. The Beat Cepheid Peterson diagrams are revisited with linear nonadiabatic turbulent convective models, and their accuracy as a probe for stellar metallicity is evaluated. They are shown to be largely independent of the helium content Y, and they are also only weakly dependent on the mass-luminosity relation that is used in their construction. However, they are found to show sensitivity to the relative abundances of the elements that are lumped into the metallicity parameter Z. Rotation is estimated to have but a small effect on the 'pulsation metallicities'. A composite Petersen diagram is presented that allows one to read off upper and lower limits on the metallicity Z from the measured period P0 and period ratio P1/P0.Comment: 9 pages, 12 color figures (black and white version available from 1st author's website). With minor revisions. to appear in Ap

The photophoretic sweeping of dust in transient protoplanetary disks

Context: Protoplanetary disks start their lives with a dust free inner region where the temperatures are higher than the sublimation temperature of solids. As the star illuminates the innermost particles, which are immersed in gas at the sublimation edge, these particles are subject to a photophoretic force. Aims: We examine the motion of dust particles at the inner edge of protoplanetary disks due to photophoretic drag. Methods: We give a detailed treatment of the photophoretic force for particles in protoplanetary disks. The force is applied to particles at the inner edge of a protoplanetary disk and the dynamical behavior of the particles is analyzed. Results: We find that, in a laminar disk, photophoretic drag increases the size of the inner hole after accretion onto the central body has become subdued. This region within the hole becomes an optically transparent zone containing gas and large dusty particles (>>10 cm), but devoid of, or strongly depleted in, smaller dust aggregates. Photophoresis can clear the inner disk of dust out to 10 AU in less than 1 Myr. The details of this clearance depend on the size distribution of the dust. Any replenishment of the dust within the cleared region will be continuously and rapidly swept out to the edge. At late times, the edge reaches a stable equilibrium between inward drift and photophoretic outward drift, at a distance of some tens of AU. Eventually, the edge will move inwards again as the disk disperses, shifting the equilibrium position back from about 40 AU to below 30 AU in 1-2 Myr in the disk model. In a turbulent disk, diffusion can delay the clearing of a disk by photophoresis. Smaller and/or age-independent holes of radii of a few AU are also possible outcomes of turbulent diffusion counteracting photophoresis. Conclusions: This outward and then inward moving edge marks a region of high dust concentration. This density enhancement, and the efficient transport of particles from close to the star to large distances away, can explain features of comets such as high measured ratios of crystalline to amorphous silicates, and has a large number of other applications

Impact of the measured parameters of exoplanets on the inferred internal structure

Exoplanet characterization is one of the main foci of current exoplanetary science. For super-Earths and sub-Neptunes, we mostly rely on mass and radius measurements, which allow to derive the body's mean density and give a rough estimate of the planet's bulk composition. However, the determination of planetary interiors is a very challenging task. In addition to the uncertainty in the observed fundamental parameters, theoretical models are limited due to the degeneracy in determining the planetary composition. We aim to study several aspects that affect internal characterization of super-Earths and sub-Neptunes: observational uncertainties, location on the M-R diagram, impact of additional constraints as bulk abundances or irradiation, and model assumptions. We use a full probabilistic Bayesian inference analysis that accounts for observational and model uncertainties. We employ a Nested Sampling scheme to efficiently produce the posterior probability distributions for all the planetary structural parameter of interest. We include a structural model based on self-consistent thermodynamics of core, mantle, high-pressure ice, liquid water, and H-He envelope. Regarding the effect of mass and radius uncertainties on the determination of the internal structure, we find three different regimes: below the Earth-like composition line and above the pure-water composition line smaller observational uncertainties lead to better determination of the core and atmosphere mass respectively, and between them structure characterization only weakly depends on the observational uncertainties. We show that small variations in the temperature or entropy profiles lead to radius variations that are comparable to the observational uncertainty, suggesting that uncertainties linked to model assumptions can become more relevant to determine the internal structure than observational uncertainties.Comment: 12 pages, 12 figure

A delta Scuti distance to the Large Magellanic Cloud

We present results from a well studied delta Scuti star discovered in the LMC. The absolute magnitude of the variable was determined from the PL relation for Galactic delta Scuti stars and from the theoretical modeling of the observed B,V,I light curves. The two methods give distance moduli for the LMC of 18.46+-0.19 and 18.48+-0.15, respectively, for a consistent value of the stellar reddening of E(B-V)=0.08+-0.02. We have also analyzed 24 delta Scuti candidates discovered in the OGLE II survey of the LMC, and 7 variables identified in the open cluster LW 55 and in the galaxy disk by Kaluzny et al. (2003, 2006). We find that the LMC delta Scuti stars define a PL relation whose slope is very similar to that defined by the Galactic delta Scuti variables, and yield a distance modulus for the LMC of 18.50+-0.22 mag. We compare the results obtained from the delta Scuti variables with those derived from the LMC RR Lyrae stars and Cepheids. Within the observational uncertainties, the three groups of pulsating stars yield very similar distance moduli. These moduli are all consistent with the "long" astronomical distance scale for the Large Magellanic Cloud.Comment: Accepted for publication on A