Planet Formation Research Articles

The TESS-Keck Survey. VII. A Superdense Sub-Neptune Orbiting TOI-1824* *Some of the data presented herein were obtained at Keck Observatory, which is a private 501(c)3 nonprofit organization operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of

We confirm a massive sub-Neptune-sized planet on a P = 22.8 days orbit around the star TOI-1824 (T eff = 5200 K, V = 9.7 mag). TESS first identified TOI-1824 b (formerly TOI-1824.01) as an object of interest in 2020 April after two transits in Sector 22 were matched with a single transit in Sector 21. TOI-1824 was subsequently targeted for ground-based Doppler monitoring with Keck-HIRES and APF-Levy. Using a joint model of the TESS photometry, radial velocities, and Ca ii H and K emission measurements as an activity indicator, we find that TOI-1824 b is an unusually dense sub-Neptune. The planet has a radius R p = 2.63 ± 0.15 R ⊕ and mass M p = 18.5 ± 3.2 M ⊕, implying a bulk density of 5.6 ± 1.4 g cm−3. TOI-1824 b's mass and radius situate it near a small group of “superdense sub-Neptunes” (R p ≲ 3 R ⊕ and M p ≳ 20 M ⊕). While the formation mechanism of superdense sub-Neptunes is a mystery, one possible explanation is the constructive collision of primordial icy cores; such giant impacts would drive atmospheric escape and could help explain these planets' apparent lack of massive envelopes. We discuss TOI-1824 b in the context of these overdense planets, whose unique location in the exoplanet mass–radius plane make them a potentially valuable tracer of planet formation.

Constraining Free–Free Emission and Photoevaporative Mass-loss Rates for Known Proplyds and New VLA–identified Candidate Proplyds in NGC 1977

We present NSF's Karl G. Jansky Very Large Array observations covering the NGC 1977 region at 3.0, 6.4, and 15.0 GHz. We search for compact radio sources and detect continuum emission from 34 NGC 1977 cluster members and 37 background objects. Of the 34 radio-detected cluster members, 3 are associated with known proplyds in NGC 1977, 22 are associated with additional young stellar objects in NGC 1977, and 9 are newly identified cluster members. We examine the radio spectral energy distributions, circular polarization, and variability of the detected NGC 1977 sources and identify 10 new candidate proplyds whose radio fluxes are dominated by optically thin free–free emission. We use measurements of free–free emission to calculate the mass-loss rates of known proplyds and new candidate proplyds in NGC 1977, and find values ∼10−9 to 10−8 M ⊙ yr−1, which are lower than the mass-loss rates measured toward proplyds in the Orion Nebula Cluster but consistent with the mass-loss rates predicted by external photoevaporation models for spatially extended disks that are irradiated by the typical external ultraviolet (UV) fields encountered in NGC 1977. Finally, we show that photoevaporative disk winds in NGC 1977 may be illuminated by internal or external sources of ionization, depending on their positions within the cluster. This study provides new constraints on disk properties in a clustered star-forming region with a weaker UV environment than the Orion Nebula Cluster but a stronger UV environment than low-mass star-forming regions like Taurus. Such intermediate UV environments represent the typical conditions of Galactic star and planet formation.

Probing the Kinematic Signatures of Planet Formation in the Gas Disk of MWC 480

The disk around MWC 480 has shown multiple substructures in both dust and gas observations, possibly suggesting ongoing planet formation in situ. In this paper, we explore the gas kinematics of the MWC 480 disk by analyzing the archival Atacama Large Millimeter/submillimeter Array observations of 12CO (J = 2-1), 13CO (J = 2-1), and C18O (J = 2-1). By modeling the line-of-sight velocities, inferred from the Doppler shifts of the emission lines, we are able to decompose the three-dimensional (3D) velocity field of the disk into rotational, radial, and vertical components. Further analysis reveals the presence of large-scale gas flows in the (r, z) plane. Notably, we identify potential meridional flows across various heights as traced by all three CO isotopologues in the 80–120 au region, possibly associated with ongoing planet formation activities in this region. Moreover, we find upward flows near 200 au for all three CO isotopologues, which may point to the presence of disk winds.

Age Analysis of Extrasolar Planets: Insight from Stellar Isochrone Models

There is growing evidence from stellar kinematics and galactic chemical evolution suggesting that giant planets (M P ≥ 0.3M J ) are relatively young compared to the most commonly occurring population of small planets (M P < 0.3M J ). To further test the validity of these results, we analyzed the ages for a large number of 2336 exoplanet hosting stars determined using three different but well-established isochrone fitting models, namely, PARSEC, MIST, and Yonsei Yale. As input parameters, we used Gaia DR3 parallaxes, magnitudes, and photometric temperature, as well as spectroscopically determined more accurate temperatures and metallicities from the Sweet Catalog. Our analysis suggests that ∼50%–70% of stars with planets are younger than the Sun. We also find that, among the confirmed exoplanetary systems, stars hosting giant planets are even younger compared to small planet hosts. The median age of ∼2.61–3.48 Gyr estimated for the giant planet-hosting stars (depending on the model input parameters) suggests that the later chemical enrichment of the galaxy by the iron-peak elements, largely produced from Type Ia supernovae, may have paved the way for the formation of gas giants. Furthermore, within the giant planet population itself, stars hosting hot Jupiters (orbital period ≤10 days) are found to be younger compared to the stellar hosts of cool and warm Jupiters (orbital period >10 days), implying that hot Jupiters could be the youngest systems to emerge in the progression of planet formation.

HED zircons as a window into the solar system’s first crust: Decoupling primordial differentiation, metamorphism and impact events through textural and chemical studies

The study of Howardite-Eucrite-Diogenite (HED) meteorites provides unique insights into early planet formation and the impact events that shaped the early Solar System. However, unraveling the complex history of the HED parent body (hypothesized to be the asteroid 4 Vesta) from whole-rock samples is challenging since most HEDs are impact-related breccias comprising mixed lithic and mineral fragments that experienced variable deformation and alteration. Combining U-Pb geochronology, trace element geochemistry, and microstructural analysis of zircon can unravel magmatic, metamorphic and impact processes through time to decipher the HED parent body evolution. Here we present textural (EBSD), geochronological (207Pb/206Pb SIMS dating) and geochemical data (Th/U, REE, Ti-in-zircon thermometry) on 61 zircon grains from melt breccia eucrites, unbrecciated/monomict/polymict eucrites, howardites and diogenites. Diverse textures indicate variable histories of impact deformation and high-temperature recrystallization. Undeformed, fractured zircons preserve primary zoning (CL, Th/U, REE) indicating magmatic and metamorphic origins. At least three magmatic zircon grains (Th/U > 0.3) give 207Pb/206Pb ages of 4558–4565 Ma, suggesting primary differentiation in the parent body first million years. Twenty metamorphic zircon grains (Th/U < 0.3) date to 4420–4568 Ma, indicating prolonged thermal metamorphism from impact heating and/or crustal cooling. Impact-recrystallized granular zircon grains reveal major impacts during and just after the parent body differentiation (4500–4560 Ma), plus later events potentially linked to synchronous impacts in the Solar System (e.g. the Moon). Similarity of metamorphic and shocked zircon ages (circa 4550–4450 Ma) suggests impacts occurred for ≥100 million years after the parent body formed.

Introducing two improved methods for approximating radiative cooling in hydrodynamical simulations of accretion discs

ABSTRACT The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate at which the system loses thermal energy. To study the evolution of these systems, we require radiative cooling approximations because full radiative transfer is generally too expensive to be coupled to hydrodynamical models. Here, we present two new approximate methods for computing radiative cooling that make use of the polytropic cooling approximation. This approach invokes the assumption that each parcel of gas is located within a spherical pseudo-cloud, which can then be used to approximate the optical depth. The first method combines the methods introduced by Stamatellos et al. and Lombardi et al. to overcome the limitations of each method at low and high optical depths, respectively. The second method, the ‘modified Lombardi’ method, is specifically tailored for self-gravitating discs. This modifies the scale height estimate from the method of Lombardi et al. using the analytical scale height for a self-gravitating disc. We show that the modified Lombardi method provides an excellent approximation for the column density in a fragmenting disc, a regime in which the existing methods fail to recover the clumps and spiral structures. We therefore recommend this improved radiative cooling method for more realistic simulations of self-gravitating discs.

Time-dependent, long-term hydrodynamic simulations of the inner protoplanetary disk. III. The influence of photoevaporation

The final stages of a protoplanetary disk are essential for our understanding of the formation and evolution of planets. Photoevaporation is an important mechanism that contributes to the dispersal of an accretion disk and has significant consequences for the disk's lifetime. However, the combined effects of photoevaporation and star-disk interaction have not been investigated in previous studies. A photoevaporative disk evolution model including a detailed formulation of the inner star-disk interaction region will improve the understanding of the final stages of disk evolution. We combined an implicit disk evolution model with a photoevaporative mass-loss profile. By including the innermost disk regions down to 0.01 AU, we could calculate the star-disk interaction, the stellar spin evolution, and the transition from an accreting disk to the propeller regime self-consistently. Starting from an early Class II star-disk system (with an age of 1 Myr), we calculated the long-term evolution of the system until the disk becomes almost completely dissolved. Photoevaporation has a significant effect on disk structure and evolution. The radial extent of the dead zone decreases, and the number of episodic accretion events (outbursts) is reduced by high stellar X-ray luminosities. Reasonable accretion rates ($&gt; 10^ M_ / yr $) in combination with photoevaporative gaps are possible for a dead zone that is still massive enough to develop episodic accretion events. Furthermore, the stellar spin evolution during the Class II evolution is less affected by the star-disk interaction in the case of high X-ray luminosities. Our results suggest that the formation of planets, especially habitable planets, in the dead zone is strongly impaired in the case of strong X-ray luminosities. Additionally, the importance of the star-disk interaction during the Class II phase with respect to the stellar spin evolution is reduced.

Terrestrial planet formation from a ring: Long-term simulations accounting for the giant planet instability

The process leading to the formation of the terrestrial planet remains elusive. In a previous publication, we have shown that, if the first generation of planetesimals forms in a ring at ∼1 AU and the gas disk's density peaks at the ring location, planetary embryos of a few martian masses can grow and remain in the ring. In this work, we extend our simulations beyond the gas-disk stage, covering ∼200 Myr and accounting for the phase of giant planet instability, assumed to happen at different times. About half of the simulations form a pair of Venus and Earth analogues and, independently, about 10% form a Mars analogue. We find that the timing of the giant planet instability affects statistically the terrestrial system's excitation state and the timing of the last giant impacts. Hence a late instability (∼60 to 100 Myr after the Solar system's birth) is more consistent with a late Moon-formation time, as suggested by radioactive chronometers. However, the late veneer mass (LVM: mass accreted after the last giant impact) of Earth-sized planets suffering a giant impact after 80 My is usually an order of magnitude lower than the value inferred from geochemistry. In addition, the final angular momentum deficit (AMD) of the terrestrial planets tends to be too high. We tested the effect on the final AMD of the generation of debris during collisions and found that it is too small to change these conclusions. We argue that the best-case scenario is that the Moon-forming event occurred between 50 and 80 My, possibly just following the giant planet instability.

Condition for metal fragmentation during Earth-forming collisions

The long-term evolution of the deep Earth depends on its initial temperature and composition. These were set by the large planetary collisions that formed the Earth. After each collision, the metallic core of the impactor fell into a molten silicate magma ocean. Previous investigations showed that, as it sank, the impactor core fragmented into drops. The overall fragmentation of the core controlled the efficiency of chemical transfers between the impactor metal and the magma ocean, and, as a consequence, the composition of the Earth's core and mantle. However, because previous studies lack an impact stage, it is unclear whether the projectile core fragmented during the impact at Earth's surface, or deeper in the magma ocean.To answer this question, we conduct laboratory experiments modeling the collision of single-phase and two-phase impactors. In a first series of experiments, we investigate the impact of a single-phase centimetric liquid volume, representing the impactor core, onto a lighter immiscible liquid, representing the magma ocean. Our experiments approach the dynamical regime of planetary collisions for which inertia is large compared to surface tension. Varying the velocity and size of the impactor, we determine the conditions under which the impactor fragments into drops. We find that fragmentation occurs when the Froude number, which measures the relative importance of inertia to gravity, is larger than 40, regardless of surface tension. This fragmentation results from the growth of a turbulent Rayleigh-Taylor instability at the interface between the impacting liquid and the target pool. In contrast, when Fr<10, the impactor remains coherent. In a second series of experiments, we use two-phase impactors to show that these results hold for impactors that are differentiated into a core and a mantle.Applied to planet formation, our results suggest that the core of impactors less than 330 km in radius impacting at the escape velocity onto an Earth-sized planet fully fragments into droplets during the impact process, whereas the core of a giant Mars-sized impactor remains coherent. We derive a model for the depth at which the impactor core fragments in the magma ocean as a function of the impactor size and velocity. This model predicts that impactors with a radius less than 800 km fully fragment before reaching the bottom of the magma ocean. For velocities higher than twice the escape speed, some degree of fragmentation is unavoidable for any impactor size.

Mixing is easy: New insights for cosmochemical evolution from pre-stellar core collapse

Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot, and high-pressure regions that formed during the gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas and dust interaction and resolve potential zones of dust concentration during star and disc formation stages. We explore whether the thermal and dynamical conditions that developed during protostellar disc formation can generate gas flows that efficiently mix and transport the well-coupled gas and dust components. We simulated the collapse of dusty molecular cloud cores with the hydrodynamics code PLUTO augmented with radiation transport and self-gravity. We used a two-dimensional axisymmetric geometry and followed the azimuthal component of the velocity. The dust was treated as Lagrangian particles that are subject to drag from the gas, whose motion is computed on a Eulerian grid. We considered 1, 10, and 100 micron-sized neutral, spherical dust grains. Importantly, the equation of state accurately includes molecular hydrogen dissociation. We focus on molecular cloud core masses of 1 and 3 $M_ odot $ and explore the effects of different initial rotation rates and cloud core sizes. Our study underlines mechanisms for the early transport of dust from the inner hot disc regions via the occurrence of two transient gas motions, namely meridional flow and outflow. The vortical flow fosters dynamical mixing and retention of dust, while the thermal pressure driven outflow replenishes dust in the outer disc. Notably, these phenomena occur regardless of the initial cloud core mass, size, and rotation rate. Young dynamical precursors to planet-forming discs exhibit regions with complex hydrodynamical gas features and high-temperature structures. These can play a crucial role in concentrating dust for subsequent growth into protoplanets. Dust transport, especially, from sub-au scales surrounding the protostar to the outer relatively cooler parts, offers an efficient pathway for thermal reprocessing during pre-stellar core collapse.

PDRs4All

Context. JWST has taken the sharpest and most sensitive infrared (IR) spectral imaging observations ever of the Orion Bar photodis-sociation region (PDR), which is part of the nearest massive star-forming region the Orion Nebula, and often considered to be the ‘prototypical’ strongly illuminated PDR. Aims. We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the H II region to the atomic PDR – crossing the ionisation front (IF) –, and the subsequent transition to the molecular PDR – crossing the dissociation front (DF). Given the prevalence of PDRs in the interstellar medium and their dominant contribution to IR radiation, understanding the response of the PDR gas to far-ultraviolet (FUV) photons and the associated physical and chemical processes is fundamental to our understanding of star and planet formation and for the interpretation of any unresolved PDR as seen by JWST. Methods. We used high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion Bar PDR as part of the PDRs4All JWST Early Release Science programme. We constructed a 3″ × 25″’ spatio-spectral mosaic covering 0.97– 5.27 μm at a spectral resolution R of ~2700 and an angular resolution of 0.075″–0.173″. To study the properties of key regions captured in this mosaic, we extracted five template spectra in apertures centred on the three H2 dissociation fronts, the atomic PDR, and the H II region. This wealth of detailed spatial-spectral information was analysed in terms of variations in the physical conditions-incident UV field, density, and temperature – of the PDR gas. Results. The NIRSpec data reveal a forest of lines including, but not limited to, He I, H I , and C I recombination lines; ionic lines (e.g. Fe III and Fe II); O I and N I fluorescence lines; aromatic infrared bands (AIBs, including aromatic CH, aliphatic CH, and their CD counterparts); pure rotational and ro-vibrational lines from H2; and ro-vibrational lines from HD, CO, and CH+, with most of them having been detected for the first time towards a PDR. Their spatial distribution resolves the H and He ionisation structure in the Huygens region, gives insight into the geometry of the Bar, and confirms the large-scale stratification of PDRs. In addition, we observed numerous smaller-scale structures whose typical size decreases with distance from θ1 Ori C and IR lines from C I , if solely arising from radiative recombination and cascade, reveal very high gas temperatures (a few 1000 K) consistent with the hot irradiated surface of small-scale dense clumps inside the PDR. The morphology of the Bar, in particular that of the H2 lines, reveals multiple prominent filaments that exhibit different characteristics. This leaves the impression of a ‘terraced’ transition from the predominantly atomic surface region to the CO-rich molecular zone deeper in. We attribute the different characteristics of the H2 filaments to their varying depth into the PDR and, in some cases, not reaching the C+/C/CO transition. These observations thus reveal what local conditions are required to drive the physical and chemical processes needed to explain the different characteristics of the DFs and the photochemical evolution of the AIB carriers. Conclusions. This study showcases the discovery space created by JWST to further our understanding of the impact radiation from young stars has on their natal molecular cloud and proto-planetary disk, which touches on star and planet formation as well as galaxy evolution.

MINDS: The JWST MIRI Mid-INfrared Disk Survey

The study of protoplanetary disks has become increasingly important with the Kepler satellite finding that exoplanets are ubiquitous around stars in our galaxy and the discovery of enormous diversity in planetary system architectures and planet properties. High-resolution near-IR and ALMA images show strong evidence for ongoing planet formation in young disks. The JWST MIRI mid-INfrared Disk Survey (MINDS) aims to (1) investigate the chemical inventory in the terrestrial planet-forming zone across stellar spectral type, (2) follow the gas evolution into the disk dispersal stage, and (3) study the structure of protoplanetary and debris disks in the thermal mid-IR. The MINDS survey will thus build a bridge between the chemical inventory of disks and the properties of exoplanets. The survey comprises 52 targets (Herbig Ae stars, T Tauri stars, very low-mass stars and young debris disks). We primarily obtain MIRI/MRS spectra with high signal-to-noise ratio (∼100–500) covering the complete wavelength range from 4.9 to 27.9 μm. For a handful of selected targets we also obtain NIRSpec IFU high resolution spectroscopy (2.87–5.27 μm). We will search for signposts of planet formation in thermal emission of micron-sized dust—information complementary to near-IR scattered light emission from small dust grains and emission from large dust in the submillimeter wavelength domain. We will also study the spatial structure of disks in three key systems that have shown signposts for planet formation, TW Hya and HD 169142 using the MIRI coronagraph at 15.5 μm and 10.65 μm respectively and PDS 70 using NIRCam imaging in the 1.87 μm narrow and the 4.8 μm medium band filter. We provide here an overview of the MINDS survey and showcase the power of the new JWST mid-IR molecular spectroscopy with the TW Hya disk spectrum where we report the detection of the molecular ion CH3+ and the robust confirmation of HCO+ earlier detected with Spitzer.

Evolution of the relation between the mass accretion rate and the stellar and disk mass from brown dwarfs to stars

The time evolution of the dependence of the mass accretion rate with the stellar mass and the disk mass represents a fundamental way to understand the evolution of protoplanetary disks and the formation of planets. In this work, we present observations with X-shooter of 26 Class II very low-mass stars (&lt; 0.2 M⊙) and brown dwarfs in the Ophiuchus, Chamaeleon-I, and Upper Scorpius star-forming regions. These new observations extend the measurement of the mass accretion rate down to spectral type (SpT) M9 (∼0.02 M⊙) in Ophiuchus and Chamaeleon-I and add 11 very-low-mass stars to the sample of objects previously studied with broadband spectroscopy in Upper Scorpius. We obtained the spectral type and extinction, as well as the physical parameters of the sources. We used the intensity of various emission lines in the spectra of these sources to derive the accretion luminosity and mass accretion rates for the entire sample. Combining these new observations with data from the literature, we compare relations between accretion and stellar and disk properties of four different star-forming regions with different ages: Ophiuchus (∼1 Myr), Lupus (∼2 Myr), Chamaeleon-I (∼3 Myr), and Upper Scorpius (5−12 Myr). We find the slopes of the accretion relationships (L* − Lacc, M∗ − Ṁacc) to steepen in the 1−3 Myr age range (i.e., between Ophiuchus, Lupus, and Chamaeleon-I) and that both relationships may be better described with a single power law. We find that previous claims for a double power-law behavior of the M∗ − Ṁacc relationship may have been triggered by the use of a different SpT–Teff scale. We also find the relationship between the protoplanetary disk mass and the mass accretion rate of the stellar population to steepen with time down to the age of Upper Scorpius. Overall, we observe hints of a faster evolution into low accretion rates of low-mass stars and brown dwarfs. At the same time, we also find that brown dwarfs present higher Mdisk/Ṁacc ratios (i.e., longer accretion depletion timescales) than stars in Ophiuchus, Lupus, and Cha-I. This apparently contradictory result may imply that the evolution of protoplanetary disks around brown dwarfs may be different than what is seen in the stellar regime.

Methods for Incorporating Model Uncertainty into Exoplanet Atmospheric Analysis

A key goal of exoplanet spectroscopy is to measure atmospheric properties, such as abundances of chemical species, in order to connect them to our understanding of atmospheric physics and planet formation. In this new era of high-quality JWST data, it is paramount that these measurement methods are robust. When comparing atmospheric models to observations, multiple candidate models may produce reasonable fits to the data. Typically, conclusions are reached by selecting the best-performing model according to some metric. This ignores model uncertainty in favor of specific model assumptions, potentially leading to measured atmospheric properties that are overconfident and/or incorrect. In this paper, we compare three ensemble methods for addressing model uncertainty by combining posterior distributions from multiple analyses: Bayesian model averaging, a variant of Bayesian model averaging using leave-one-out predictive densities, and stacking of predictive distributions. We demonstrate these methods by fitting the Hubble Space Telescope (HST) + Spitzer transmission spectrum of the hot Jupiter HD 209458b using models with different cloud and haze prescriptions. All of our ensemble methods lead to uncertainties on retrieved parameters that are larger but more realistic and consistent with physical and chemical expectations. Since they have not typically accounted for model uncertainty, uncertainties of retrieved parameters from HST spectra have likely been underreported. We recommend stacking as the most robust model combination method. Our methods can be used to combine results from independent retrieval codes and from different models within one code. They are also widely applicable to other exoplanet analysis processes, such as combining results from different data reductions.

Fractionation in young cores: Direct determinations of nitrogen and carbon fractionation in HCN

Context. Nitrogen fractionation is a powerful tracer of the chemical evolution during star and planet formation. It requires robust determinations of the nitrogen fractionation across different evolutionary stages. Aims. We aim to determine the 14N/15N and 12C/13C ratios for HCN in six starless and prestellar cores and to compare the results between the direct method using radiative transfer modeling and the indirect double isotope method, assuming a fixed 12C/13C ratio. Methods. We present IRAM observations of the HCN 1–0, HCN 3–2, HC15N 1–0 and H13CN 1–0 transitions toward six embedded cores. The 14N/15N ratio was derived using both the indirect double isotope method and directly through non-local thermodynamic equilibrium (NLTE) 1D radiative transfer modeling of the HCN emission. The latter also provides the 12C/13C ratio, which we compared to the local interstellar value. Results. The derived 14N/15N ratios using the indirect method are generally in the range of 300-550. This result could suggest an evolutionary trend in the nitrogen fractionation of HCN between starless cores and later stages of the star formation process. However, the direct method reveals lower fractionation ratios of around ~250, mainly resulting from a lower 12C/13C ratio in the range ~20–40, as compared to the local interstellar medium value of 68. Conclusions. This study reveals a significant difference between the nitrogen fractionation ratio in HCN derived using direct and indirect methods. This can influence the interpretation of the chemical evolution and reveal the pitfalls of the indirect double isotope method for fractionation studies. However, the direct method is challenging, as it requires well-constrained source models to produce accurate results. No trend in the nitrogen fractionation of HCN between earlier and later stages of the star formation process is evident when the results of the direct method are considered.

Where Are the Water Worlds? Identifying Exo-water-worlds Using Models of Planet Formation and Atmospheric Evolution

Planet formation models suggest that the small exoplanets that migrate from beyond the snowline of the protoplanetary disk likely contain water-ice-rich cores (∼50% by mass), also known as water worlds. While the observed radius valley of the Kepler planets is well explained by the atmospheric dichotomy of the rocky planets, precise measurements of the mass and radius of the transiting planets hint at the existence of these water worlds. However, observations cannot confirm the core compositions of those planets, owing to the degeneracy between the density of a bare water-ice-rich planet and the bulk density of a rocky planet with a thin atmosphere. We combine different formation models from the Genesis library with atmospheric escape models, such as photoevaporation and impact stripping, to simulate planetary systems consistent with the observed radius valley. We then explore the possibility of water worlds being present in the currently observed sample by comparing them with simulated planets in the mass–radius–orbital period space. We find that the migration models suggest ≳10% and ≳20% of the bare planets, i.e., planets without primordial H/He atmospheres, to be water-ice-rich around G- and M-type host stars, respectively, consistent with the mass–radius distributions of the observed planets. However, most of the water worlds are predicted to be outside a period of 10 days. A unique identification of water worlds through radial velocity and transmission spectroscopy is likely to be more successful when targeting such planets with longer orbital periods.

The Evolution of Accreting Binaries: From Brown Dwarfs to Supermassive Black Holes

Circumbinary accretion occurs throughout the universe, from the formation of stars and planets to the aftermath of major galactic mergers. We present an extensive investigation of circumbinary accretion disks, studying circular binaries with mass ratios (q ≡ M 2/M 1) from 0.01 to 1 and at each mass ratio probing the effects of disk thickness and viscosity. We study disks with aspect ratios H/r ∈ {0.1, 0.05, 0.03} and vary both the magnitude and spatial dependence of viscosity. Although thin accretion disks have previously been found to promote rapid inspirals of equal-mass binaries, we find that gravitational torques become weaker at lower mass ratios and most binaries with 0.01 ≤ q ≤ 0.04 outspiral, which may delay the coalescence of black hole binaries formed from minor mergers and cause high-mass exoplanets to migrate outward. However, in a number of cases, the disks accreting onto binaries with mass ratios ∼0.07 fail to develop eccentric modes, leading to extremely rapid inspirals. Variability in black hole accretion correlates with disk eccentricity, and we observe variability above the ∼10% level even for mass ratios of 0.01. We demonstrate that the spatial dependence of the viscosity (e.g., α vs. constant ν) significantly affects the degree of preferential accretion onto the secondary, resolving discrepancies between previous studies. Colder circumbinary disks remain eccentric even at q ∼ 0.01 and sustain deep, asymmetric cavities.

Constraining the formation of WASP-39b using JWST transit spectroscopy

Context. Understanding the formation history of planets is one of the goals of studying exoplanet atmospheres. The atmospheric composition of planets can provide insights into the formation pathways of planets. Even though the mapping of the atmospheric composition onto a formation pathway is not unambiguous, with the increasing sensitivity of modern instruments, we can derive promising constraints. Aims. In this work, we aim to understand the formation pathway of WASP-39b. We discuss whether the detection of SO2 in its atmosphere would impact our understanding of the formation of the planet and whether it enables us to determine the formation pathway of the planet with greater accuracy. Methods. We used the JWST transit observation of the planet together with the available HST and Spitzer observations. We used a formation model coupled with a radiative transfer retrieval model to derive the planet's atmospheric characteristics and formation history. Furthermore, we used a photochemical model to derive the impact of photochemistry on the atmosphere of the planet. Results. In this work, we show that the planet is most likely to have initiated beyond the CO2 ice line of its natal disk. Furthermore, the planet is likely to have have accreted some planetesimals during its formation. We show that the sulfur abundance in the atmosphere of the planet is probably lower than 2.27 × 10−4. This abundance indicates that the planet is likely to exhibit a lower metallicity than suggested by the retrievals. Furthermore, such an abundance for sulfur is more likely if WASP-39b had been formed beyond the CO ice line of its natal disk.

Study of type II migration under the framework of the disk instability model for giant planet formation

Context. Hydrodynamic simulations of the migration of planets formed by gravitational instability suggest that after an initial phase of fast migration, planets can open gaps and continue to migrate on a type II migration timescale. The simulation time length is typically on the order of 104 yr. Aims. We study the effects of the subsequent type II migration during the disk lifetime on the final orbital radii of planets. Methods. We used a numerical disk model that follows the disk formation and evolution. The disk acquires mass through the mass influx from the collapse of its parent molecular cloud core. The model reflects the influence of the properties of the parent core on the disk. Considering clumps forming at different times in a disk and also in different disks with different parent core properties, we used the type II migration rate to follow the clump migration from the formation location. We studied the dependence of the clump migration on the properties of the parent core. Results. The mass influx drag enhances the migration process. The duration and viscosity of gravitational instability, viscosity in the dead zone, and the collapse time of the parent core play important roles in planet migration. As the angular momentum and mass of the parent core increase, migration is enhanced. The final radius is sensitive to the initial radius. Clumps forming at large radii might migrate outward with the disk expansion. Conclusions. Even though type II migration is slow, clumps can migrate over significant distances. A considerable proportion of clumps migrate to the central protostar via type II migration. Our calculations support the idea that the observed pile-up of planets at &lt;0.3 AU is explained by a scenario where planets might form at large radii, then migrate to orbits of &lt;0.3 AU, and halt by a stopping mechanism at this location.

Planets similar in size are often dissimilar in interior

The number of discovered exoplanets now exceeds 5500 allowing statistical analyses of planetary systems. Multi-planet systems are mini-laboratories of planet formation and evolution, and analysing their system architectures can help us to constrain the physics of these processes. Recent works have found evidence of significant intrasystem uniformity in planet properties such as radius, mass, and orbital spacing, collectively termed ‘peas in a pod’ trends. In particular, correlations in radius and mass have been interpreted as implying uniformity in planet bulk density and composition within a system. However, the samples used to assess trends in mass tend to be small and biased. In this paper, we re-evaluate correlations in planet properties in a large sample of systems with at least two planets for which mass and radius have been directly measured, and therefore bulk density can be calculated. Our sample was assembled using the most up-to-date exoplanet catalogue data, and we compute the relevant statistics while using a procedure to ‘weight’ the data points according to measurement precision. We find a moderate correlation in radius and a weak correlation in the densities of adjacent planets. However, masses of neighbouring planets show no overall correlation in our main sample and a weak correlation among pairs of planets similar in size or pairs restricted to Mp&lt;100 M⊕, Rp&lt;10 R⊕. Similarly, we show that the intrasystem dispersion in radius is typically less than that in mass and density. We identify ranges in stellar host properties that correlate with stronger uniformity in pairs of adjacent planets: low Teff for planet masses, and low metallicity and old age for planet densities. Furthermore, we explore whether peas in a pod trends extend into planet compositions or interior structures. For small neighbouring planets with similar radii, we show that their masses and interior structures are often disparate, indicating that even within the same system, similarity in radii is not necessarily a good proxy for similarity in composition or the physical nature of the planets.