ABSTRACT The composition of a planet's atmosphere is determined by its formation, evolution, and present-day insolation. A planet's spectrum therefore may hold clues on its origins. We present a ..."chain" of models, linking the formation of a planet to its observable present-day spectrum. The chain links include (1) the planet's formation and migration, (2) its long-term thermodynamic evolution, (3) a variety of disk chemistry models, (4) a non-gray atmospheric model, and (5) a radiometric model to obtain simulated spectroscopic observations with James Webb Space Telescope and ARIEL. In our standard chemistry model the inner disk is depleted in refractory carbon as in the Solar System and in white dwarfs polluted by extrasolar planetesimals. Our main findings are: (1) envelope enrichment by planetesimal impacts during formation dominates the final planetary atmospheric composition of hot Jupiters. We investigate two, under this finding, prototypical formation pathways: a formation inside or outside the water iceline, called "dry" and "wet" planets, respectively. (2) Both the "dry" and "wet" planets are oxygen-rich (C/O < 1) due to the oxygen-rich nature of the solid building blocks. The "dry" planet's C/O ratio is <0.2 for standard carbon depletion, while the "wet" planet has typical C/O values between 0.1 and 0.5 depending mainly on the clathrate formation efficiency. Only non-standard disk chemistries without carbon depletion lead to carbon-rich C/O ratios >1 for the "dry" planet. (3) While we consistently find C/O ratios <1, they still vary significantly. To link a formation history to a specific C/O, a better understanding of the disk chemistry is thus needed.
Context.
The anomalously large radii of hot Jupiters are still not fully understood, and all of the proposed explanations are based on the idea that these close-in giant planets possess hot ...interiors. Most of the mechanisms proposed have been tested on a handful of exoplanets.
Aims.
We approach the radius anomaly problem by adopting a statistical approach. We want to infer the internal luminosity for the sample of hot Jupiters, study its effect on the interior structure, and put constraints on which mechanism is the dominant one.
Methods.
We developed a flexible and robust hierarchical Bayesian model that couples the interior structure of exoplanets to the observed properties of close-in giant planets. We applied the model to 314 hot Jupiters and inferred the internal luminosity distribution for each planet and studied at the population level (i) the mass–luminosity–radius distribution and as a function of equilibrium temperature the distributions of the (ii) heating efficiency, (iii) internal temperature, and the (iv) pressure of the radiative–convective–boundary (RCB).
Results.
We find that hot Jupiters tend to have high internal luminosity with 10
4
L
J
for the largest planets. As a result, we show that all the inflated planets have hot interiors with an internal temperature ranging from 200 up to 800 K for the most irradiated ones. This has important consequences on the cooling rate and we find that the RCB is located at low pressures between 3 and 100 bar. Assuming that the ultimate source of the extra heating is the irradiation from the host star, we also illustrate that the heating efficiency increases with increasing equilibrium temperature and reaches a maximum of 2.5% at ~1860 K, beyond which the efficiency decreases, which is in agreement with previous results. We discuss our findings in the context of the proposed heating mechanisms and illustrate that ohmic dissipation, the advection of potential temperature, and thermal tides are in agreement with certain trends inferred from our analysis and thus all three models can explain various aspects of the observations.
Conclusions.
We provide new insights on the interior structure of hot Jupiters and show that with our current knowledge, it is still challenging to firmly identify the universal mechanism driving the inflated radii.
ABSTRACT Many parameters constraining the spectral appearance of exoplanets are still poorly understood. We therefore study the properties of irradiated exoplanet atmospheres over a wide parameter ...range including metallicity, C/O ratio, and host spectral type. We calculate a grid of 1D radiative-convective atmospheres and emission spectra. We perform the calculations with our new Pressure-Temperature Iterator and Spectral Emission Calculator for Planetary Atmospheres (PETIT) code, assuming chemical equilibrium. The atmospheric structures and spectra are made available online. We find that atmospheres of planets with C/O ratios ∼1 and 1500 K can exhibit inversions due to heating by the alkalis because the main coolants CH4, H2O, and HCN are depleted. Therefore, temperature inversions possibly occur without the presence of additional absorbers like TiO and VO. At low temperatures we find that the pressure level of the photosphere strongly influences whether the atmospheric opacity is dominated by either water (for low C/O) or methane (for high C/O), or both (regardless of the C/O). For hot, carbon-rich objects this pressure level governs whether the atmosphere is dominated by methane or HCN. Further we find that host stars of late spectral type lead to planetary atmospheres which have shallower, more isothermal temperature profiles. In agreement with prior work we find that for planets with K the transition between water or methane dominated spectra occurs at C/O ∼ 0.7, instead of ∼1, because condensation preferentially removes oxygen.
The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has recently been detected by direct imaging in the near-infrared. To study the interaction between ...the protoplanet and the disk itself, observations of the gas and dust surface density structure are needed. This paper reports new ALMA observations of the dust continuum at 1.3 mm, 12CO, 13CO, and C18O J = 2−1 emission from the system HD 169142 (which is observed almost face-on) at an angular resolution of ~\hbox{$0\farcs3 \,\times\, 0\farcs2$}0 .̋ 3 × 0 .̋ 2 (~35 × 20 au). The dust continuum emission reveals a double-ring structure with an inner ring between \hbox{$0\farcs17{-}0\farcs28$}0 .̋ 17−0 .̋ 28 (~20−35 au) and an outer ring between \hbox{$0\farcs48{-}0\farcs64$}0 .̋ 48−0 .̋ 64 (~56−83 au). The size and position of the inner ring is in good agreement with previous polarimetric observations in the near-infrared and is consistent with dust trapping by a massive planet. No dust emission is detected inside the inner dust cavity (R ≲ 20 au) or within the dust gap (~35−56 au) down to the noise level. In contrast, the channel maps of the J = 2−1 line of the three CO isotopologs reveal gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission; it extends to ~1\hbox{$\farcs5$}.̋ 5 (~180 au) in radius. This difference and the sharp drop of the continuum emission at large radii point to radial drift of large dust grains (>μm size). Using the thermo-chemical disk code dali, we modeled the continuum and the CO isotopolog emission to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of ~30−40 inward of the dust gap. The gas and dust distribution indicate that two giant planets shape the disk structure through dynamical clearing (dust cavity and gap) and dust trapping (double-ring dust distribution).
ABSTRACT The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887 m survey of ...the disk population around objects from ∼2 to 0.03 M in the nearby ∼2 Myr old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45 au in radius. Assuming isothermal and optically thin emission, we convert the 887 m flux densities into dust disk masses, hereafter Mdust. We find that the relation is steeper than linear and of the form Mdust ∝ (M*)1.3-1.9, where the range in the power-law index reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By reanalyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr old regions of Taurus, Lupus, and Chamaeleon I share the same relation, while the 10 Myr old Upper Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations.
Formation and evolution of water in the solar system and the origin of water on Earth constitute one of the most interesting questions in astronomy. In this paper, the isotopic and chemical evolution ...of water during the early history of the solar nebula, before the onset of planetesimal formation, is studied. A gas-grain chemical model that includes multiply deuterated species and nuclear spin-states is combined with a steady-state solar nebula model. To calculate initial abundances, we simulated 1 Myr of evolution of a cold and dark TMC-1-like prestellar core. Two time-dependent chemical models of the solar nebula are calculated over 1 Myr: (1) a laminar model and (2) a model with two-dimensional (2D) turbulent mixing. We find that with regards to the water isotopic composition and the origin of the comets, the mixing model seems to be favored over the laminar model.
Context. Through observations numerous giant molecular filaments (GMFs) have been discovered in the Milky Way. Their role in the Galactic star formation and Galaxy-scale evolution of dense gas is ...unknown. Aims. We investigate systematically the star-forming content of all currently known GMFs. This allows us to estimate the star formation rates (SFRs) of the GMFs and to establish relationships between the SFRs and the GMF properties. Methods. We identified and classified the young stellar object (YSO) population of each GMF using multiwavelength photometry from near- to far-infrared. We estimated the total SFRs assuming a universal and fully sampled initial mass function and luminosity function. Results. We uniformly estimate the physical properties of 57 GMFs. The GMFs show correlations between the 13CO line width, mass, and size, similar to Larson’s relations. We identify 36 394 infrared excess sources in 57 GMFs and obtain SFRs for 46 GMFs. The median SFR surface density (ΣSFR) and star formation efficiency (SFE) of GMFs are 0.62 M⊙ Myr−1 pc−2 and 1%, similar to the nearby star-forming clouds. The star formation rate per free-fall time of GMFs is between 0.002−0.05 with the median value of 0.02. We also find a strong correlation between SFR and dense gas mass that is defined as gas mass above a visual extinction of 7 mag, which suggests that the SFRs of the GMFs scale similarly with dense gas as those of nearby molecular clouds. We also find a strong correlation between the mean SFR per unit length and dense gas mass per unit length. The origin of this scaling remains unknown, calling for further studies that can link the structure of GMFs to their SF activity and explore the differences between GMFs and other molecular clouds.
The growth of solid particles towards meter sizes in protoplanetary disks has to circumvent at least two hurdles, namely the rapid loss of material due to radial drift and particle fragmentation due ...to destructive collisions. In this paper, we present the results of numerical simulations with more and more realistic physics involved. Step by step, we include various effects, such as particle growth, radial/vertical particle motion and dust particle fragmentation in our simulations. We demonstrate that the initial dust-to-gas ratio is essential for the particles to overcome the radial drift barrier. If this value is increased by a factor of 2 compared with the canonical value for the interstellar medium, km-sized bodies can form in the inner disk (<2 AU) within 104 yrs. However, we find that solid particles get destroyed through collisional fragmentation. Only with the unrealistically high-threshold velocities needed for fragmentation to occur (>30 m/s), particles are able to grow to larger sizes in disks with low α values. We also find that less than 5% of the small dust grains remain in the disk after 1 Myr due to radial drift, no matter whether fragmentation is included in the simulations or not. In this paper, we also present considerable improvements to existing algorithms for dust-particle coagulation, which speed up the coagulation scheme by a factor of ~ 104.
ABSTRACT The total gas mass of a protoplanetary disk is a fundamental, but poorly determined, quantity. A new technique has been demonstrated to assess directly the bulk molecular gas reservoir of ...molecular hydrogen using the HD J = 1-0 line at 112 m. In this work we present a Herschel Space Observatory10 survey of six additional T Tauri disks in the HD line. Line emission is detected at >3 significance in two cases: DM Tau and GM Aur. For the other four disks, we establish upper limits to the line flux. Using detailed disk structure and ray-tracing models, we calculate the temperature structure and dust mass from modeling the observed spectral energy distributions, and we include the effect of UV gas heating to determine the amount of gas required to fit the HD line. The ranges of gas masses are 1.0-4.7 × 10−2 for DM Tau and 2.5-20.4 × 10−2 for GM Aur. These values are larger than those found using CO for GM Aur, while the CO-derived gas mass for DM Tau is consistent with the lower end of our mass range. This suggests a CO chemical depletion from the gas phase of up to a factor of five for DM Tau and up to two orders of magnitude for GM Aur. We discuss how future analysis can narrow the mass ranges further.
ABSTRACT
The presence of CO gas around 10–50 Myr old A stars with debris discs has sparked debate on whether the gas is primordial or secondary. Since secondary gas released from planetesimals is ...poor in H2, it was thought that CO would quickly photodissociate never reaching the high levels observed around the majority of A stars with bright debris discs. Kral et al. showed that neutral carbon produced by CO photodissociation can effectively shield CO and potentially explain the high CO masses around 9 A stars with bright debris discs. Here, we present a new model that simulates the gas viscous evolution, accounting for carbon shielding and how the gas release rate decreases with time as the planetesimal disc loses mass. We find that the present gas mass in a system is highly dependant on its evolutionary path. Since gas is lost on long time-scales, it can retain a memory of the initial disc mass. Moreover, we find that gas levels can be out of equilibrium and quickly evolving from a shielded on to an unshielded state. With this model, we build the first population synthesis of gas around A stars, which we use to constrain the disc viscosity. We find a good match with a high viscosity (α ∼ 0.1), indicating that gas is lost on time-scales ∼1–10 Myr. Moreover, our model also shows that high CO masses are not expected around FGK stars since their planetesimal discs are born with lower masses, explaining why shielded discs are only found around A stars. Finally, we hypothesize that the observed carbon cavities could be due to radiation pressure or accreting planets.