Understanding Of Planet Formation Research Articles

FARGOCPT: 2D Multiphysics code for simulating disk interactions with stars, planets, and particles

Context. Planet-disk interactions play a crucial role in the understanding of planet formation and disk evolution. There are multiple numerical tools available to simulate these interactions, including the commonly used FARGO code and its variants. Many of the codes have been extended over time to include additional physical processes, with a focus on their accurate modeling. Aims. We introduce FARGOCPT, an updated version of FARGO that incorporates other previous enhancements to the code, to provide a simulation environment tailored to studies of the interactions between stars, planets, and disks. It is meant to ensure an accurate representation of planet systems, hydrodynamics, and dust dynamics, with a focus on usability. Methods. The radiation-hydrodynamics part of FARGOCPT uses a second-order upwind scheme in 2D polar coordinates, supporting multiple equations of state, radiation transport, heating and cooling, and self-gravity. Shocks are considered using artificial viscosity. The integration of the N-body system is achieved by leveraging the REBOUND code. The dust module utilizes massless tracer particles, adapted to drag laws for the Stokes and Epstein regimes. Moreover, FARGOCPT provides mechanisms to simulate accretion onto stars and planets. Results. The code has been tested in practice in the context of multiple studies. Additionally, it comes with an automated test suite for checking the physics modules. It is available online. Conclusions. FARGOCPT offers a unique set of simulation capabilities within the current landscape of publicly available planet-disk interaction simulation tools. Its structured interface and underlying technical updates are intended to assist researchers in ongoing explorations of planet formation.

Dust growth and pebble formation in the initial stages of protoplanetary disk evolution

Aims. The initial stages of planet formation may start concurrently with the formation of a gas-dust protoplanetary disk. This makes the study of the earliest stages of protoplanetary disk formation crucially important. Here we focus on dust growth and pebble formation in a protoplanetary disk that is still accreting from a parental cloud core. Methods. We have developed an original three-dimensional numerical hydrodynamics code, which computes the collapse of rotating clouds and disk formation on nested meshes using a novel hybrid Coarray Fortran-OpenMP approach for distributed and shared memory parallelization. Dust dynamics and growth are also included in the simulations. Results. We found that the dust growth from ~1 µm to 1–10 mm already occurs in the initial few thousand years of disk evolution but the Stokes number hardly exceeds 0.1 because of higher disk densities and temperatures compared to the minimum mass Solar nebular. The ratio of the dust-to-gas vertical scale heights remains rather modest, 0.2–0.5, which may be explained by the perturbing action of spiral arms that develop in the disk soon after its formation. The dust-to-gas mass ratio in the disk midplane is highly nonhomogeneous throughout the disk extent and is in general enhanced by a factor of several compared to the fiducial 1:100 value. Low St hinders strong dust accumulation in the spiral arms compared to the rest of the disk and the nonsteady nature of the spirals is also an obstacle. The spatial distribution of pebbles in the disk midplane exhibits a highly nonhomogeneous and patchy character. The total mass of pebbles in the disk increases with time and reaches a few tens of Earth masses after a few tens of thousand years of disk evolution. Conclusions. We found that protoplanetary disks with an age ≤20 kyr can possess notable amounts of pebbles and feature dust-togas density enhancements in the disk midplane. Hence, these young disks can already be ripe for the planet formation process to start. Multidimensional numerical models of disk formation that consider the coevolution of gas and dust including dust growth are important to improve our understanding of planet formation.

Synthetic populations of protoplanetary disks: Impact of magnetic fields and radiative transfer

Context. Protostellar disks are the product of angular momentum conservation during protostellar collapse. Understanding their formation is crucial because they are the birthplace of planets and their formation is also tightly related to star formation. Unfortunately, the initial properties of Class 0 disks and their evolution are still poorly constrained both theoretically and observationally. Aims. We aim to better understand the mechanisms that set the statistics of disk properties as well as to study their formation in massive protostellar clumps. We also want to provide the community with synthetic disk populations to better interpret young disk observations. Methods. We used the ramses code to model star and disk formation in massive protostellar clumps with magnetohydrodynamics, including the effect of ambipolar diffusion and radiative transfer as well as stellar radiative feedback. Those simulations, resolved up to the astronomical unit scale, have allowed us to investigate the formation of disk populations. Results. Magnetic fields play a crucial role in disk formation. A weaker initial field leads to larger and massive disks and weakens the stellar radiative feedback by increasing fragmentation. We find that ambipolar diffusion impacts disk and star formation and leads to very different disk magnetic properties. The stellar radiative feedback also have a strong influence, increasing the temperature and reducing fragmentation. Comparing our disk populations with observations reveals that our models with a mass-to-flux ratio of 10 seems to better reproduce observed disk sizes. This also sheds light on a tension between models and observations for the disk masses. Conclusions. The clump properties and physical modeling significantly impact disk populations. It is critical to for the tension, with respect to disk mass estimates, between observations and models to be solved with synthetic observations. This is particularly important in the context of understanding planet formation.

Constraints on the Dust Size Distributions in the HD 163296 Disk from the Difference of the Apparent Dust Ring Widths between Two ALMA Bands

The dust size in protoplanetary disks is a crucial parameter for understanding planet formation, while the observational constraints on dust size distribution have large uncertainties. In this study, we present a new method to constrain the dust size distribution from the dust spatial distribution, utilizing the fact that larger dust grains are more spatially localized. We analyze Atacama Large Millimeter/submillimeter Array Band 6 (1.25 mm) and Band 4 (2.14 mm) high-resolution images and constrain the dust size distribution in the two rings of the HD 163296 disk. We find that the outer ring at 100 au appears narrower at longer wavelengths, while the inner ring at 67 au appears to have similar widths across the two wavelengths. We model dust rings trapped at gas pressure maxima, where the dust grains follow a power-law size distribution, and the dust grains of a specific size follow a Gaussian spatial distribution, with the width depending on the grain size. By comparing the observations with the models, we constrain the maximum dust size and the exponent of the dust size distribution p. We constrain and p < 3.3 in the inner ring, and and 3.4 < p < 3.7 in the outer ring. The larger maximum dust size in the outer ring implies a spatial dependency in dust growth, potentially influencing the formation location of the planetesimals. We further discuss the turbulence strength α derived from the constrained dust spatial distribution, assuming equilibrium between turbulent diffusion and the accumulation of dust grains.

Retrieval of the dayside atmosphere of WASP-43b with CRIRES+

Accurately estimating the C/O ratio of hot Jupiter atmospheres is a promising pathway towards understanding planet formation and migration, as well as the formation of clouds and the overall atmospheric composition. The atmosphere of the hot Jupiter WASP-43b has been extensively analysed using low-resolution observations with HST and Spitzer, but these previous observations did not cover the K band, which hosts prominent spectral features of major carbon-bearing species such as CO and CH4. As a result, the ability to establish precise constraints on the C/O ratio was limited. Moreover, the planet has not been studied at high spectral resolution, which can provide insights into the atmospheric dynamics. In this study, we present the first high-resolution dayside spectra of WASP-43b with the new CRIRES+ spectrograph. By observing the planet in the K band, we successfully detected the presence of CO and provide evidence for the existence of H2O using the cross-correlation method. This discovery represents the first direct detection of CO in the atmosphere of WASP-43b. Furthermore, we retrieved the temperature-pressure profile, abundances of CO and H2O, and a super-solar C/O ratio of 0.78 by applying a Bayesian retrieval framework to the data. Our findings also shed light on the atmospheric characteristics of WASP-43b. We found no evidence for a cloud deck on the dayside, and recovered a line broadening indicative of an equatorial super-rotation corresponding to a jet with a wind speed of ~5kms−1, matching the results of previous forward models and low-resolution atmospheric retrievals for this planet.

Estimating the number of planets that PLATO can detect

Context. The PLATO mission is scheduled for launch in 2026. It will monitor more than 245 000 FGK stars of magnitude 13 or brighter for planet transit events. Among the key scientific goals are the detection of Earth-Sun analogs; the detailed characterization of stars and planets in terms of mass, radius, and ages; the detection of planetary systems with longer orbital periods than are detected in current surveys; and to advance our understanding of planet formation and evolution processes. Aims. This study aims to estimate the number of exoplanets that PLATO can detect as a function of planetary size and period, stellar brightness, and observing strategy options. Deviations from these estimates will be informative of the true occurrence rates of planets, which helps constraining planet formation models. Methods. For this purpose, we developed the Planet Yield for PLATO estimator (PYPE), which adopts a statistical approach. We apply given occurrence rates from planet formation models and from different search and vetting pipelines for the Kepler data. We estimate the stellar sample to be observed by PLATO using a fraction of the all-sky PLATO stellar input catalog (PIC). PLATO detection efficiencies are calculated under different assumptions that are presented in detail in the text. Results. The results presented here primarily consider the current baseline observing duration of 4 yr. We find that the expected PLATO planet yield increases rapidly over the first year and begins to saturate after 2 yr. A nominal (2+2) 2-yr mission could yield about several thousand to several tens of thousands of planets, depending on the assumed planet occurrence rates. We estimate a minimum of 500 Earth-size (0.8−1.25 RE) planets, about a dozen of which would reside in a 250–500 days period bin around G stars. We find that one-third of the detected planets are around stars bright enough (V ≤11) for RV-follow-up observations. We find that a 3-yr-long observation followed by 6 two-month short observations (3+1 yr) yield roughly twice as many planets as two long observations of 2 yr (2+2 yr). The former strategy is dominated by short-period planets, while the latter is more beneficial for detecting earths in the habitable zone. Conclusions. Of the many sources of uncertainties for the PLATO planet yield, the real occurrence rates matters most. Knowing the latter is crucial for using PLATO observations to constrain planet formation models by comparing their statistical yields.

An inflationary disk phase to explain extended protoplanetary dust disks

Context. Understanding planetesimal formation is an essential first step towards understanding planet formation. The distribution of these first solid bodies drives the locations where planetary embryos, which eventually form fully-fledged planets, grow. Aims. We seek to understand the parameter space of possible protoplanetary disk formation and evolution models of our Solar System. A good protoplanetary disk scenario for the Solar System must meet at least the following three criteria: (1) It must produce an extended gas and dust disk (e.g. 45 au for the dust); (2) within the disk, the local dust-to-gas ratio in at least two distinct locations must sufficiently increase to explain the early formation of the parent bodies of non-carbonaceous and carbonaceous iron meteorites; and (3) dust particles, which have condensed at high temperatures (i.e. calcium–aluminium-rich inclusions), must be transported to the outer disk. Though current protoplanetary disk models are able to satisfy one or two of these criteria, none have been successful in recreating all three. We aim to find scenarios that satisfy all three. Methods. In this study we used a 1D disk model that tracks the evolution of the gas and dust disks. Planetesimals are formed within the disk at locations where the streaming instability can be triggered. We explored a large parameter space to study the effect of the disk viscosity, the timescale of infall of material into the disk, the distance within which material is deposited into the disk, and the fragmentation threshold of dust particles. Results. We find that scenarios with a large initial disk viscosity (α &gt; 0.05), a relatively short infall timescale (Tinfall &lt; 100–200kyr), and a small centrifugal radius (RC ~ 0.4 au; i.e. the distance within which material falls into the disk) result in disks that satisfy all three criteria needed to represent the protoplanetary disk of the Solar System. The large initial viscosity and short infall timescale result in a rapid initial expansion of the disk, which we dub the ‘inflationary phase’ of the disk. Furthermore, a temperature-dependent fragmentation threshold, which accounts for cold icy particles breaking more easily, results in larger and more massive disks. This, in turn, results in more ‘icy’ than ‘rocky’ planetesimals. Such scenarios are also better in line with our Solar System, which has small terrestrial planets and massive giant planet cores. Finally, we find that scenarios with large RC cannot transport calcium–aluminium-rich inclusions to the outer disk and do not produce planetesimals at two locations within the disk.

The stability of unevenly spaced planetary systems

Studying the orbital stability of multi-planet systems is essential to understand planet formation, estimate the stable time of an observed planetary system, and advance population synthesis models. Although previous studies have primarily focused on ideal systems characterized by uniform orbital separations, in reality a diverse range of orbital separations exists among planets within the same system. This study focuses on investigating the dynamical stability of systems with non-uniform separation. We considered a system with 10 planets with masses of 10−7 solar masses around a central star with a mass of 1 solar mass. We performed more than 100,000 runs of N-body simulations with different parameters. Results demonstrate that reducing merely one pair of planetary spacing leads to an order of magnitude shorter orbital crossing times that could be formulated based on the Keplerian periods of the closest separation pair. Furthermore, the first collisions are found to be closely associated with the first encounter pair that is likely to be the closest separation pair initially. We conclude that when estimating the orbital crossing time and colliding pairs in a realistic situation, updating the formula derived for evenly spaced systems would be necessary.

Inner Planetary System Gap Complexity is a Predictor of Outer Giant Planets

The connection between inner small planets and outer giant planets is crucial to our understanding of planet formation across a wide range of orbital separations. While Kepler provided a plethora of compact multiplanet systems at short separations (≲1 au), relatively little is known about the occurrence of giant companions at larger separations and how they impact the architectures of the inner systems. Here, we use the catalog of systems from the Kepler Giant Planet Search to study how the architectures of the inner transiting planets correlate with the presence of outer giant planets. We find that for systems with at least three small transiting planets, the distribution of inner-system gap complexity (), a measure of the deviation from uniform spacings, appears to differ (p ≲ 0.02) between those with an outer giant planet () and those without any outer giants. All four inner systems (with three or more transiting planets) with outer giant(s) have a higher gap complexity () than 79% (19/24) of the inner systems without any outer giants (median ). This suggests that one can predict the occurrence of outer giant companions by selecting multitransiting systems with highly irregular spacings. We do not find any correlation between the outer giant occurrence and the size (similarity or ordering) patterns of the inner planets. The higher gap complexities of inner systems with an outer giant hints that massive external planets play an important role in the formation and/or disruption of the inner systems.

Threshold velocity for the collisional growth of porous dust aggregates consisting of cohesive frictionless spheres

Understanding the collisional outcomes of dust aggregates and their dependence on the material properties of the constituting particles is of great importance for understanding planet formation. Recent numerical simulations have revealed that interparticle tangential friction plays a crucial role in energy dissipation during collisions between porous dust aggregates, but the importance of friction for the collisional growth of dust aggregates remains poorly understood. Here we demonstrate the effects of interparticle tangential friction on the collisional growth of dust aggregates. We performed numerical simulations of collisions between equal-mass porous dust aggregates consisting of cohesive and frictionless spheres. We changed the collision velocity and impact angle systematically and calculated the collisional growth efficiency as a function of the collision velocity. We found that the threshold velocity for collisional growth decreases when dust aggregates are made of frictionless spheres compared to frictional spheres. Our results highlight the importance of tangential interactions for the collisional behavior of dust aggregates and indicate that the predictive equation for the threshold velocity should be reconstructed.

Stable Sn isotope signatures of Mid-ocean ridge basalts

The constraint of Sn isotope composition of the Earth's mantle is the prerequisite for understanding planet formation with Sn isotopes and Sn isotope fractionation in terrestrial reservoirs. However, previous estimates of the Sn isotope composition of the Earth's mantle are hampered by the small number of mantle-related samples. In this study, we analyzed Sn stable isotope compositions of 27 basalts from geochemically and geologically diverse mid-ocean ridge segments, to constrain the Sn isotope composition of the depleted mantle source. The spread of Sn isotope compositions for the 27 different mid-ocean ridge basalts (MORBs) is limited, with δ122/118SnBHVO-2 (delta notation of 122Sn/118Sn ratios relative to USGS BHVO-2 basalt standard) values ranging from −0.051 ± 0.015‰ to 0.209 ± 0.016‰ (δ122/118Sn3161a from 0.280 ± 0.015 to 0.540 ± 0.016‰), despite great geochemical and geographical diversities of the samples. The slightly higher δ122/118SnBHVO-2 value (up to 0.209 ± 0.016‰ (δ122/118Sn3161a up to 0.540 ± 0.016‰)) of one MORB sample appears to be caused by post-eruptive alteration of the basalts. Other than that, the δ122/118SnBHVO-2 values of the MORB samples do not show correlations with sample latitude, spreading rate of mid-ocean ridge, MgO content of the bulk rock samples, and partial melting index Na8.0 (defined as [Na2O] + 0.373 × [MgO] − 2.98), implying that Sn isotopes do not fractionate significantly during MORB melt generation and evolution processes. This study confirms the limited Sn isotopic variability between fresh MORBs globally, pointing to the Sn isotopic homogeneity of the depleted mantle source. Using these new MORB data, we proposed an estimate of the δ122/118SnBHVO-2 value for the Earth's depleted mantle to be 0.036 ± 0.087‰ (δ122/118Sn3161a of 0.367 ± 0.087‰) (2SD, N = 12). This value provides a reference point for understanding the planetary and magmatic processes of Earth from a Sn isotope perspective.

The chemical inventory of the inner regions of planet-forming disks - the JWST/MINDS program.

The understanding of planet formation has changed recently, embracing the new idea of pebble accretion. This means that the influx of pebbles from the outer regions of planet-forming disks to their inner zones could determine the composition of planets and their atmospheres. The solid and molecular components delivered to the planet-forming region can be best characterized by mid-infrared spectroscopy. With Spitzer low-resolution (R = 100, 600) spectroscopy, this approach was limited to the detection of abundant molecules, such as H2O, C2H2, HCN and CO2. This contribution will present the first results of the MINDS (MIRI mid-INfrared Disk Survey, PI:Th Henning) project. Due do the sensitivity and spectral resolution provided by the James Webb Space Telescope (JWST), we now have a unique tool to obtain the full inventory of chemistry in the inner disks of solar-type stars and brown dwarfs, including also less-abundant hydrocarbons and isotopologues. The Integral Field Unit (IFU) capabilities will enable at the same time spatial studies of the continuum and line emission in extended sources such as debris disks, the flying saucer and also the search for mid-IR signatures of forming planets in systems such as PDS 70. These JWST observations are complementary to ALMA and NOEMA observations of outer-disk chemistry; together these datasets will provide an integral view of the processes occurring during the planet-formation phase.

Importance of Sample Selection in Exoplanet-atmosphere Population Studies

Understanding planet formation requires robust population studies, which are designed to reveal trends in planet properties. In this work we aim to determine if and how different methods for selecting populations of exoplanets for atmospheric characterization with JWST could influence population-level inferences. We generate three hypothetical surveys of super-Earths/sub-Neptunes, with each survey designed to span a similar radius-insolation flux space. The survey samples are constructed based on three different selection criteria (evenly spaced by eye, binned, and a quantitative selection function). Using an injection-recovery technique, we test how robustly individual-planet atmospheric parameters and population-level parameters can be retrieved. We find that all three survey designs result in equally suitable targets for individual atmospheric characterization, but not equally suitable targets for constraining population parameters. Only samples constructed with a quantitative method or that are sufficiently evenly spaced-by-eye result in robust population parameter constraints. Furthermore, we find that the sample with the best targets for individual atmospheric study does not necessarily result in the best-constrained population parameters. The method of sample selection must be considered. We also find that there may be large variability in population-level results with a sample that is small enough to fit in a single JWST cycle (∼12 planets), suggesting that the most successful population-level analyses will be multicycle. Lastly, we infer that our exploration of sample selection is limited by the small number of transiting planets with measured masses around bright stars. Our results can guide future development of programs that aim to determine underlying trends in exoplanet-atmospheric properties, and, by extension, formation and evolution processes.

A giant planet shaping the disk around the very low-mass star CIDA 1

Context. Exoplanetary research has provided us with exciting discoveries of planets around very low-mass (VLM) stars (0.08 M⊙ ≲ M* ≲ 0.3 M⊙; e.g., TRAPPIST-1 and Proxima Centauri). However, current theoretical models still strive to explain planet formation in these conditions and do not predict the development of giant planets. Recent high-resolution observations from the Atacama Large Millimeter/submillimeter Array (ALMA) of the disk around CIDA 1, a VLM star in Taurus, show substructures that hint at the presence of a massive planet. Aims. We aim to reproduce the dust ring of CIDA 1, observed in the dust continuum emission in ALMA Band 7 (0.9 mm) and Band 4 (2.1 mm), along with its 12CO (J = 3−2) and 13CO (J = 3−2) channel maps, assuming the structures are shaped by the interaction of the disk with a massive planet. We seek to retrieve the mass and position of the putative planet, through a global simulation that assesses planet-disk interactions to quantitatively reproduce protoplanetary disk observations of both dust and gas emission in a self-consistent way. Methods. Using a set of hydrodynamical simulations, we model a protoplanetary disk that hosts an embedded planet with a starting mass of between 0.1 and 4.0 MJup and initially located at a distance of between 9 and 11 au from the central star. We compute the dust and gas emission using radiative transfer simulations, and, finally, we obtain the synthetic observations, treating the images as the actual ALMA observations. Results. Our models indicate that a planet with a minimum mass of ~1.4 MJup orbiting at a distance of ~9−10 au can explain the morphology and location of the observed dust ring in Band 7 and Band 4. We match the flux of the dust emission observation with a dust-to-gas mass ratio in the disk of ~10−2. We are able to reproduce the low spectral index (~2) observed where the dust ring is detected, with a ~40−50% fraction of optically thick emission. Assuming a 12CO abundance of 5 × 10−5 and a 13CO abundance 70 times lower, our synthetic images reproduce the morphology of the 12CO (J = 3−2) and 13CO (J = 3−2) observed channel maps where the cloud absorption allowed a detection. From our simulations, we estimate that a stellar mass M* = 0.2 M⊙ and a systemic velocity vsys = 6.25 km s−1 are needed to reproduce the gas rotation as retrieved from molecular line observations. Applying an empirical relation between planet mass and gap width in the dust, we predict a maximum planet mass of ~4−8 MJup. Conclusions. Our results suggest the presence of a massive planet orbiting CIDA 1, thus challenging our understanding of planet formation around VLM stars.

Discovery of Line Pressure Broadening and Direct Constraint on Gas Surface Density in a Protoplanetary Disk

The gas surface density profile of protoplanetary disks is one of the most fundamental physical properties to understanding planet formation. However, it is challenging to determine the surface density profile observationally, because the H2 emission cannot be observed in low-temperature regions. We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) archival data of the 12CO J = 3 − 2 line toward the protoplanetary disk around TW Hya and discovered extremely broad line wings due to the pressure broadening. In conjunction with a previously reported optically thin CO isotopologue line, the pressure broadened line wings enabled us to directly determine the midplane gas density for the first time. The gas surface density at ∼5 au from the central star reaches ∼103 g cm−2, which suggests that the inner region of the disk has enough mass to form a Jupiter-mass planet. Additionally, the gas surface density drops at the inner cavity by ∼2 orders of magnitude compared to outside the cavity. We also found a low CO abundance of ∼10−6 with respect to H2, even inside the CO snow line, which suggests conversion of CO to less volatile species. Combining our results with previous studies, the gas surface density jumps at r ∼ 20 au, suggesting that the inner region (3 < r < 20 au) might be the magnetorotational instability dead zone. This study sheds light on the direct gas surface density constraint without assuming the CO/H2 ratio using ALMA.

Depletion of gaseous CO in protoplanetary disks by surface-energy-regulated ice formation

Empirical constraints of fundamental properties of protoplanetary disks are essential for understanding planet formation and planetary properties1,2. Carbon monoxide (CO) gas is often used to constrain disk properties3. However, estimates show that the CO gas abundance in disks is depleted relative to expected values4,5,6,7, and models of various disk processes impacting the CO abundance could not explain this depletion on observed ~1 Myr timescales8,9,10,11,12,13,14. Here we demonstrate that surface energy effects on particles in disks, such as the Kelvin effect, that arise when ice heterogeneously nucleates onto an existing particle can efficiently trap CO in its ice phase. In previous ice formation models, CO gas was released when small ice-coated particles were lofted to warmed disk layers. Our model can reproduce the observed abundance, distribution and time evolution of gaseous CO in the four most studied protoplanetary disks7. We constrain the solid and gaseous CO inventory at the midplane and disk diffusivities and resolve inconsistencies in estimates of the disk mass—three crucial parameters that control planetary formation.

The Ariel Target List: The Impact of TESS and the Potential for Characterizing Multiple Planets within a System

The ESA Ariel mission has been adopted for launch in 2029 and will conduct a survey of around 1000 exoplanetary atmospheres during its primary mission life. By providing homogeneous data sets with a high signal-to-noise ratio and wide wavelength coverage, Ariel will unveil the atmospheric demographics of these faraway worlds, helping to constrain planet formation and evolution processes on a galactic scale. Ariel seeks to undertake a statistical survey of a diverse population of planets; therefore, the sample of planets from which this selection can be made is of the utmost importance. While many suitable targets have already been found, hundreds more will be discovered before the mission is operational. Previous studies have used predictions of exoplanet detections to forecast the available planet population by the launch date of Ariel, with the most recent noting that the Transiting Exoplanet Survey Satellite (TESS) alone should provide over 1000 potential targets. In this work, we consider the planet candidates found to date by TESS to show that, with the addition of already confirmed planets, Ariel will already have a more than sufficient sample to choose its target list from once these candidates are validated. We showcase the breadth of this population, as well as exploring, for the first time, the ability of Ariel to characterize multiple planets within a single system. Comparative planetology of worlds orbiting the same star, as well as across the wider population, will undoubtedly revolutionize our understanding of planet formation and evolution.

Asteroids accretion, differentiation, and break-up in the Vesta source region: Evidence from cosmochemistry of mesosiderites

The cosmochemistry of meteorites provides unique clues on asteroids accretion, differentiation, collisional break-up and reassembly - processes of critical importance for understanding planet formation in the early solar system. Mesosiderites are a complex group of achondrites whose nearly 50:50 metal-silicate composition is interpreted in the literature as resulting from the mixing of core and crustal materials derived from differentiated asteroids. Because of their complex nature and contrasting geochemical, isotopic, and spectroscopic data, the formation mechanism of mesosiderites is still poorly understood and is open to a large variety of planetary differentiation scenarios and collisional histories. In this study, based on new petrographic and geochemical data of 16 mesosiderites, we investigate in detail the proposal that mesosiderites are related to the howardite-eucrite-diogenite (HED) meteorite group, whose parent body is widely considered to be asteroid 4 Vesta (∼500 km diameter), the target of the recent NASA’s Dawn mission. We present the first high precision oxygen isotope analyses on the matrix of a set of mesosiderite samples, coupled with new chemical and petrographic analyses of mesosiderites Um Hadid, Estherville, and Mount Padbury. Concordant Δ17O values between mesosiderites (–0.241 ± 0.015 (2σ)) and howardite-eucrite-diogenites (−0.241 ± 0.017‰ (2σ)) indicate that they derived from the same oxygen isotope reservoir, but petrological evidence, in particular the distinctly lower Fe/Mn ratios and the larger lithological diversity in mesosiderites, indicates that they formed within different parent bodies. This suggests that mesosiderites and howardite-eucrite-diogenites originated in distinct parent bodies that accreted in the 4 Vesta source region but experienced different geologic evolution in terms of crustal differentiation and impact history, which were more complex and catastrophic in the mesosiderite parent body.

Five Key Exoplanet Questions Answered via the Analysis of 25 Hot-Jupiter Atmospheres in Eclipse

Population studies of exoplanets are key to unlocking their statistical properties. So far, the inferred properties have been mostly limited to planetary, orbital, and stellar parameters extracted from, e.g., Kepler, radial velocity, and Gaia data. More recently an increasing number of exoplanet atmospheres have been observed in detail from space and the ground. Generally, however, these atmospheric studies have focused on individual planets, with the exception of a couple of works that have detected the presence of water vapor and clouds in populations of gaseous planets via transmission spectroscopy. Here, using a suite of retrieval tools, we analyze spectroscopic and photometric data of 25 hot Jupiters, obtained with the Hubble and Spitzer Space Telescopes via the eclipse technique. By applying the tools uniformly across the entire set of 25 planets, we extract robust trends in the thermal structure and chemical properties of hot Jupiters not obtained in past studies. With the recent launch of the James Webb Space Telescope and the upcoming missions Twinkle and Ariel, population-based studies of exoplanet atmospheres, such as the one presented here, will be a key approach to understanding planet characteristics, formation, and evolution in our galaxy.

The mass and size of Herbig disks as seen by ALMA

Context. Many population studies have been performed over the past decade with the Atacama Large millimeter/submillimeter Array (ALMA) to understand the bulk properties of protoplanetary disks around young stars. The studied populations have mostly consisted of late spectral type (i.e., G, K &amp; M) stars, with relatively few more massive Herbig stars (spectral types B, A &amp; F). With Gaia-updated distances, now is a good time to use ALMA archival data for a Herbig disk population study and take an important step forward in our understanding of planet formation. Aims. The aim of this work is to determine the masses and sizes of all Herbig dust disks observed with ALMA to date in a volume-limited sample out to 450 pc. These masses and sizes are put in the context of the Lupus and Upper Sco T Tauri disk populations. Methods. ALMA Band 6 and Band 7 archival data of 36 Herbig stars are used, making this work 64% complete out to 225 pc, and 38% complete out to 450 pc also including Orion. Using stellar parameters and distances, the dust masses and sizes of the disks are determined via a curve-of-growth method. Survival analysis is used to obtain cumulative distributions of the dust masses and radii. Results. Herbig disks have a higher dust mass than the T Tauri disk populations of Lupus and Upper Sco by factors of ~3 and ~7 respectively. In addition, Herbig disks are often larger than the typical T Tauri disk. Although the masses and sizes of Herbig disks extend over a similar range to those of T Tauri disks, the distributions of masses and sizes of Herbig disks are significantly skewed toward higher values. Lastly, group I disks are more massive than group II disks. An insufficient number of group II disks are observed at sufficient angular resolution to determine whether or not they are also small in size compared to group I disks. Conclusions. Herbig disks are skewed towards more massive and larger dust disks compared to T Tauri disks. Based on this we speculate that these differences find their origin in an initial disk mass that scales with the stellar mass, and that subsequent disk evolution enlarges the observable differences, especially if (sub)millimeter continuum optical depth plays a role. Moreover, the larger disk masses and sizes of Herbig stars could be linked to the increasing prevalence of giant planets with host star mass.