Showing 1–50 of 104 results for author: Bitsch, B

There is no disk mass budget problem of planet formation

Authors: Sofia Savvidou, Bertram Bitsch

Abstract: The inferred dust masses from Class II protoplanetary disk observations are lower than or equal to the masses of the observed exoplanet systems. This poses the question of how planets form if their natal environments do not contain enough mass. This hypothesis has entered the literature as the "mass budget problem" of planet formation. We utilize numerical simulations of planet formation via pebbl… ▽ More The inferred dust masses from Class II protoplanetary disk observations are lower than or equal to the masses of the observed exoplanet systems. This poses the question of how planets form if their natal environments do not contain enough mass. This hypothesis has entered the literature as the "mass budget problem" of planet formation. We utilize numerical simulations of planet formation via pebble and gas accretion, including migration, in a viscously evolving protoplanetary disk, while tracing the time evolution of the dust mass. As expected, we find that the presence of a giant planet in the disk can influence the evolution of the disk itself and prevent rapid dust mass loss by trapping the dust outside its orbit. Early formation is crucial for giant planet formation, as we found in our previous work; therefore, our findings strengthen the hypothesis that planet formation has already occurred or is ongoing in Class II disks. Most importantly, we find that the optically thin dust mass significantly underestimates the total dust mass in the presence of a dust-trapping deep gap. We also show that the beam convolution would smear out the feature from a deep gap, especially if the planet forms in the inner disk. Such hidden dust mass, along with early planet formation, could be the answer to the hypothetical mass budget problem. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 8 pages, 4 figures, submitted to A&A

A formation pathway for terrestrial planets with moderate water content involving atmospheric-volatile recycling

Authors: Jonas Müller, Bertram Bitsch, Aaron David Schneider

Abstract: Of the many recently discovered terrestrial exoplanets, some are expected to harbor moderate water mass fractions of a few percent. The formation pathways that can produce planets with these water mass fractions are not fully understood. Here, we use the code chemcomp, which consists of a semi-analytical 1D protoplanetary disk model harboring a migrating and accreting planet, to model the growth a… ▽ More Of the many recently discovered terrestrial exoplanets, some are expected to harbor moderate water mass fractions of a few percent. The formation pathways that can produce planets with these water mass fractions are not fully understood. Here, we use the code chemcomp, which consists of a semi-analytical 1D protoplanetary disk model harboring a migrating and accreting planet, to model the growth and composition of planets with moderate water mass fractions by pebble accretion in a protoplanetary disk around a TRAPPIST-1 analog star. This star is accompanied by seven terrestrial planets, of which the outer four planets likely contain water mass fractions of between 1\% and 10\%. We adopt a published model that considers the evaporation of pebbles in the planetary envelope, from where recycling flows can transport the volatile vapor back into the disk. We find that with this model, the planetary water content depends on the influx rate of pebbles onto the planet. A decreasing pebble influx with time reduces the envelope temperature and consequently allows the formation of planets with moderate water mass fractions as inferred for the outer TRAPPIST-1 planets for a number of different simulation configurations. This is further evidence that the recycling of vapor is an important component of planet formation needed to explain the vast and diverse population of exoplanets. △ Less

Submitted 30 June, 2024; v1 submitted 13 June, 2024; originally announced June 2024.

Comments: Accepted by A&A, 19 pages, 8 figures

Mind the gap: Distinguishing disc substructures and their impact on the inner disc composition

Authors: Jingyi Mah, Sofia Savvidou, Bertram Bitsch

Abstract: Improved observational technologies have enabled the resolution of substructures and the measurement of chemical abundances in discs. Understanding the chemical composition of the inner disc allows us to infer the building blocks available for planet formation. Recently, the depletion of water in the inner disc has been suggested to be linked to the presence of substructures, such as gaps and ring… ▽ More Improved observational technologies have enabled the resolution of substructures and the measurement of chemical abundances in discs. Understanding the chemical composition of the inner disc allows us to infer the building blocks available for planet formation. Recently, the depletion of water in the inner disc has been suggested to be linked to the presence of substructures, such as gaps and rings, further out. We investigate this hypothesis further by running 1D semi-analytical models of a disc with a gap to understand the combined effects of disc viscosity, gap depth, gap location, and gap formation timescales on the inner disc composition. Our results show that for a specific value of disc viscosity, the simulation outcome can be classified into three regimes: shallow gap, "traffic jam", and deep gap. While deep gaps may already be distinguishable with moderate-resolution techniques, it is still challenging to resolve shallow gaps with the current capabilities. On the other hand, discs with traffic jams have a higher chance of being resolved when observed with a high resolution, but they may appear as an intensity enhancement or even featureless when observed with moderate to low angular resolution. In this regard, information on the inner disc composition is useful because it can help to infer the existence of traffic jams or distinguish them from deep gaps: discs with deep gaps are expected to have a low water content and thus high C/O ratio in the inner disc due to the effective blocking of pebbles, while discs with shallow gaps would demonstrate the opposite trend. Furthermore, discs with a traffic jam would have a constant inward flux of water-rich pebbles resulting in a moderate water content and sub-stellar C/O ratios. Finally, we find that the effectiveness of gaps as pebble barriers diminishes quickly when they form late, as most of the pebbles have already drifted inwards. △ Less

Submitted 15 June, 2024; v1 submitted 10 June, 2024; originally announced June 2024.

Comments: Replace with language-edited version

Journal ref: A&A 686, L17 (2024)

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

A recipe for eccentricity and inclination damping for partial gap opening planets in 3D disks

Authors: Gabriele Pichierri, Bertram Bitsch, Elena Lega

Abstract: In a previous paper we showed that, like the migration speed, the eccentricity damping efficiency is modulated linearly by the depth of the partial gap a planet carves in the disk surface density profile, resulting in less efficient $e$-damping compared to the prescription commonly used in population synthesis works. Here, we extend our analysis to 3D, refining our $e$-damping formula and studying… ▽ More In a previous paper we showed that, like the migration speed, the eccentricity damping efficiency is modulated linearly by the depth of the partial gap a planet carves in the disk surface density profile, resulting in less efficient $e$-damping compared to the prescription commonly used in population synthesis works. Here, we extend our analysis to 3D, refining our $e$-damping formula and studying how the inclination damping efficiency is also affected. We perform high resolution 3D locally isothermal hydrodynamical simulations of planets with varying masses embedded in disks with varying aspect ratios and viscosities. We extract the gap profile and orbital damping timescales for fixed eccentricities and inclinations up to the disk scale height. The limit in gap depths below which vortices appear, in the low-viscosity case, happens roughly at the transition between classical type-I and type-II migration regimes. The orbital damping timescales can be described by two linear trends with a break around gap depths $\sim80\%$ and with slopes and intercepts depending on the eccentricity and inclination. These trends are understood on physical grounds and are reproduced by simple fitting formulas whose error is within the typically uncertainty of type-I torque formulas. Thus, our recipes for the gap depth and orbital damping efficiencies yield a simple description for planet-disk interactions to use in N-body codes in the case of partial gap opening planets that is consistent with high-resolution 3D hydro-simulations. Finally, we show examples of how our novel orbital damping prescription can affect the outcome of population synthesis experiments. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ

At least one in a dozen stars exhibits evidence of planetary ingestion

Authors: Fan Liu, Yuan-Sen Ting, David Yong, Bertram Bitsch, Amanda Karakas, Michael T. Murphy, Meridith Joyce, Aaron Dotter, Fei Dai

Abstract: Stellar chemical compositions can be altered by ingestion of planetary material and/or planet formation which removes refractory material from the proto-stellar disc. These "planet signatures" appear as correlations between elemental abundance differences and the dust condensation temperature. Detecting these planet signatures, however, is challenging due to unknown occurrence rates, small amplitu… ▽ More Stellar chemical compositions can be altered by ingestion of planetary material and/or planet formation which removes refractory material from the proto-stellar disc. These "planet signatures" appear as correlations between elemental abundance differences and the dust condensation temperature. Detecting these planet signatures, however, is challenging due to unknown occurrence rates, small amplitudes, and heterogeneous star samples with large differences in stellar ages, and therefore stars born together (i.e., co-natal) with identical compositions can facilitate such detections. While previous spectroscopic studies were limited to small number of binary stars, the Gaia satellite provides new opportunities for detecting stellar chemical signatures of planets among co-moving pairs of stars confirmed to be co-natal. Here we report high-precision chemical abundances for a homogeneous sample of 91 co-natal pairs of stars with a well-defined selection function and identify at least seven new instances of planetary ingestion, corresponding to an occurrence rate of 8%. An independent Bayesian indicator is deployed, which can effectively disentangle the planet signatures from other factors, such as random abundance variation and atomic diffusion. Our study provides new evidence of planet signatures and facilitates a deeper understanding of the star-planet-chemistry connection by providing new observational constraints on the mechanisms of planet engulfment, formation and evolution. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 29 pages, 11 figures. Author's submitted version before final edits. Published in Nature on March 21, 2024: https://www.nature.com/articles/s41586-024-07091-y

Changing disc compositions via internal photoevaporation

Authors: Julia Lena Lienert, Bertram Bitsch, Thomas Henning

Abstract: The chemical evolution of protoplanetary discs is not fully understood, several factors influence the final distribution of disc material. One such factor are inward drifting and evaporating pebbles that enrich the inner disc with vapour. In particular, it is first enriched with water vapour, resulting in a low C/O ratio, before carbon-rich gas from the outer disc is transported inwards elevating… ▽ More The chemical evolution of protoplanetary discs is not fully understood, several factors influence the final distribution of disc material. One such factor are inward drifting and evaporating pebbles that enrich the inner disc with vapour. In particular, it is first enriched with water vapour, resulting in a low C/O ratio, before carbon-rich gas from the outer disc is transported inwards elevating the C/O ratio again. However, it is unclear how internal photoevaporation, which carries away gas and opens gaps that block inward drifting pebbles, affects the chemical composition of the disc. We aim to study these effects in discs around solar-like stars, where we especially focus on the C/O ratio and the water content. The simulations are carried out using a semi-analytical 1D disc model. Our code chemcomp includes viscous evolution and heating, pebble growth and drift, pebble evaporation and condensation, and a simple chemical partitioning model. We show that internal photoevaporation plays a major role in the (chemical) evolution of protoplanetary discs: As it opens a gap, inward drifting pebbles are stopped and cannot contribute to the volatile content any more. In addition, gas from the outer disc is carried away by photoevaporative winds. Consequently, the C/O ratio in the inner disc is low. In contrast, gaps opened by giant planets allow the gas to pass, resulting in an elevated C/O ratio, similar to viscous discs without internal photoevaporation. This will enable us to distinguish observationally between these two scenarios when measuring the C/O ratio, implying that we can infer the cause of gap structures in disc observations. In the case of a photoevaporative disc, we additionally find an elevated water content in the inner disc as the water vapour and ice undergo a cycle of evaporation/re-condensation, preventing its inward accretion onto the star. △ Less

Submitted 15 February, 2024; v1 submitted 14 February, 2024; originally announced February 2024.

Which stars can form planets: Planetesimal formation at low metallicities

Authors: Geoffrey Andama, Jingyi Mah, Bertram Bitsch

Abstract: The exoplanet diversity has been linked to the disc environment in which they form, where the host star metallicity and the formation pathways play a crucial role. In the context of the core accretion paradigm, the initial stages of planet formation require the growth of dust material from micrometre size to planetesimal size bodies before core accretion can kick in. Although numerous studies have… ▽ More The exoplanet diversity has been linked to the disc environment in which they form, where the host star metallicity and the formation pathways play a crucial role. In the context of the core accretion paradigm, the initial stages of planet formation require the growth of dust material from micrometre size to planetesimal size bodies before core accretion can kick in. Although numerous studies have been conducted on planetesimal formation, it is still poorly understood how this process takes place in low metallicity stellar environments. We explore planetesimals formation in stellar environments primarily with low metallicity. We performed global 1D viscous disc evolution simulations including grain growth, evaporation and condensation of chemical species at ice lines. We followed the formation of planetesimals during disc evolution and tested different metallicities, disc sizes and turbulent viscosity strengths. We find that at solar and sub-solar metallicities, there is a significant enhancement in the midplane dust-to-gas mass ratios at the ice lines but this leads to planetesimal formation only at the water ice line. In our simulations, [Fe/H] = -0.6 is the lowest limit of metallicity for planetesimal formation where a few Earth masses of planetesimals could form. For such extreme disc environments, large discs are more conducive than small discs for forming large amounts of planetesimals at a fixed metallicity, because the pebble flux can be maintained for a longer time resulting in a longer and more efficient planetesimal formation phase. At lower metallicities, planetesimal formation is less supported in quiescent discs compared to turbulent discs, because the pebble flux can be maintained for a longer time. The amount of planetesimals formed at sub-solar metallicities in our simulations places a limit on core sizes that could possibly only result in the formation of super-Earths. △ Less

Submitted 5 February, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A

Journal ref: A&A 683, A118 (2024)

chemcomp: Modeling the chemical composition of planets formed in protoplanetary disks

Authors: Aaron David Schneider, Bertram Bitsch

Abstract: Future observations of exoplanets will hopefully reveal detailed constraints on planetary compositions. Recently, we have developed and introduced chemcomp (Schneider & Bitsch 2021a), which simulates the formation of planets in viscously evolving protoplanetary disks by the accretion of pebbles and gas. The chemical composition of planetary building blocks (pebbles and gas) is traced by including… ▽ More Future observations of exoplanets will hopefully reveal detailed constraints on planetary compositions. Recently, we have developed and introduced chemcomp (Schneider & Bitsch 2021a), which simulates the formation of planets in viscously evolving protoplanetary disks by the accretion of pebbles and gas. The chemical composition of planetary building blocks (pebbles and gas) is traced by including a physical approach of the evaporation and condensation of volatiles at evaporation lines. We have now open-sourced the chemcomp code to enable comparisons between planet formation models and observational constraints by the community. The code can be found at https://github.com/AaronDavidSchneider/chemcomp, is easy to use (using configuration files) and comes with a detailed documentation and examples. △ Less

Submitted 22 November, 2023; originally announced January 2024.

Comments: url: https://github.com/AaronDavidSchneider/chemcomp

Constraining the formation history of the HAT-P-11 system by atmospheric abundances

Authors: Lena Chatziastros, Bertram Bitsch, Aaron David Schneider

Abstract: The chemical fingerprint of a planet reveals information about its formation history regarding when and where it formed. The water content of a planet can help to constrain its formation pathway: If the planet formed in the outer regions of the disk and migrated inward, it will be water-rich due to the accretion of water-ice-rich solids. Conversely, formation in the inner disk, where water-ice is… ▽ More The chemical fingerprint of a planet reveals information about its formation history regarding when and where it formed. The water content of a planet can help to constrain its formation pathway: If the planet formed in the outer regions of the disk and migrated inward, it will be water-rich due to the accretion of water-ice-rich solids. Conversely, formation in the inner disk, where water-ice is not available, will result in a smaller atmospheric water content due to the limited accretion of water vapor. This process complicates with the presence of gap-opening giant planets. A gas giant exerts a pressure bump exterior to its orbit, preventing the further influx of pebbles into the inner system, resulting in a water-poor environment. The different formation scenarios can help to constrain the formation of the HAT-P-11 system, which contains an inner sub-Neptune with a mass of 23.4 $\mathrm{M_{\oplus}}$ and substellar water abundances ($X_\mathrm{H_2O} \approx 0.11$), as well as an outer giant planet orbiting exterior to the water-ice line. Our planet formation model encompasses planetary growth through pebble and gas accretion, along with a pebble drift and evaporation module that enables us to track the chemical composition of the disk and the planets. We find that the presence of the gas giant is necessary to block water-ice-rich material, resulting in a substellar water content for the inner sub-Neptune. If the giant planet forms too early, not enough solid material can enter the inner disk regions, preventing the efficient growth of the inner planet. This highlights the importance of the timing of giant planet formation to explain the inner system structure. Our simulations predict a roughly stellar C/O ratio with superstellar C/H and O/H for HAT-P-11b, providing constraints for future observations that are essential for gaining a more detailed understanding of its formation. △ Less

Submitted 25 October, 2023; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 14 pages, accepted by A&A

Journal ref: A&A 681, A52 (2024)

Planetary system architectures with low-mass inner planets: Direct imaging exploration of mature systems beyond 1 au

Authors: Celia Desgrange, Julien Milli, Gael Chauvin, Thomas Henning, Anna Luashvili, Matthew Read, Mark Wyatt, Grant Kennedy, Remo Burn, Martin Schlecker, Flavien Kiefer, Valentina D'Orazi, Sergio Messina, Pascal Rubini, Anne-Marie Lagrange, Carine Babusiaux, Luca Matra, Bertram Bitsch, Mariangela Bonavita, Philippe Delorme, Elisabeth Matthews, Paulina Palma-Bifani, Arthur Vigan

Abstract: The discovery of planets orbiting at less than 1 au from their host star and less massive than Saturn in various exoplanetary systems revolutionized our theories of planetary formation. The fundamental question is whether these close-in low-mass planets could have formed in the inner disk interior to 1 au, or whether they formed further out in the planet-forming disk and migrated inward. Exploring… ▽ More The discovery of planets orbiting at less than 1 au from their host star and less massive than Saturn in various exoplanetary systems revolutionized our theories of planetary formation. The fundamental question is whether these close-in low-mass planets could have formed in the inner disk interior to 1 au, or whether they formed further out in the planet-forming disk and migrated inward. Exploring the role of additional giant planets in these systems may help us to pinpoint their global formation and evolution. We searched for additional substellar companions by using direct imaging in systems known to host close-in small planets. The use of direct imaging complemented by radial velocity and astrometric detection limits enabled us to explore the giant planet and brown dwarf demographics around these hosts to investigate the potential connection between both populations. We carried out a direct imaging survey with VLT/SPHERE to look for outer giant planets and brown dwarf companions in 27 systems hosting close-in low-mass planets discovered by radial velocity. Our sample is composed of very nearby (<20pc) planetary systems, orbiting G-, K-, and M-type mature (0.5-10Gyr) stellar hosts. We performed homogeneous direct imaging data reduction and analysis to search for and characterize point sources, and derived robust statistical detection limits. Of 337 point-source detections, we do not find any new bound companions. We recovered the emblematic very cool T-type brown dwarf GJ229B. Our typical sensitivities in direct imaging range from 5 to 30 MJup beyond 2 au. The non-detection of massive companions is consistent with predictions based on models of planet formation by core accretion. Our pilot study opens the way to a multi-technique approach for the exploration of very nearby exoplanetary systems with future ground-based and space observatories. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 49 pages including 31 pages of appendices and references, 31 figures, A&A, accepted

Journal ref: A&A 680, A64 (2023)

Composition of giant planets: the roles of pebbles and planetesimals

Authors: Claudia Danti, Bertram Bitsch, Jingyi Mah

Abstract: One of the current challenges of planet formation theory is to explain the enrichment of observed exoplanetary atmospheres. Past studies have focused on scenarios where either pebbles or planetesimals were the heavy element enrichment's drivers, we combine here both approaches to understand whether the composition of a planet can constrain its formation pathway. We study three different formation… ▽ More One of the current challenges of planet formation theory is to explain the enrichment of observed exoplanetary atmospheres. Past studies have focused on scenarios where either pebbles or planetesimals were the heavy element enrichment's drivers, we combine here both approaches to understand whether the composition of a planet can constrain its formation pathway. We study three different formation scenarios: pebble accretion, pebble accretion with planetesimal formation, combined pebble and planetesimal accretion. We use the chemcomp code to perform semi-analytical 1D simulations of protoplanetary discs, including viscous evolution, pebble drift, and simple chemistry to simulate the growth of planets from planetary embryos to gas giants as they migrate through the disc, while tracking their composition. Our simulations confirm that the composition of the planetary atmosphere is dominated by the accretion of gas enriched by inward drifting and evaporating pebbles. Including planetesimal formation hinders the enrichment, because the pebbles locked into planetesimals cannot evaporate and enrich the disc. This results in a big drop of the accreted heavy elements both in the planetesimal formation and accretion case, proving that planetesimal formation needs to be inefficient in order to explain planets with high heavy element content. Accretion of planetesimals enhances the refractory component of the atmosphere, leading to low volatile to refractory ratios, contrary to the pure pebble scenario. Such low volatile to refractory ratios can also be achieved by planets migrating in the inner disc in pure pebble scenario. Distinguishing these two scenarios requires knowledge about the planet's atmospheric C/H and O/H ratios, which are higher for pure pebble accretion. Therefore, a detailed knowledge of the composition of planetary atmospheres could help to constrain the planet's formation pathway. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 12 pages, 10 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 679, L7 (2023)

How to make giant planets via pebble accretion

Authors: Sofia Savvidou, Bertram Bitsch

Abstract: Planet formation is directly linked to the birthing environment that protoplanetary disks provide. The disk properties determine whether a giant planet will form and how it evolves. The number of exoplanet and disk observations is consistently rising, however, it is not yet possible to directly link these two populations. Therefore, a deep theoretical understanding of how planets form is crucial.… ▽ More Planet formation is directly linked to the birthing environment that protoplanetary disks provide. The disk properties determine whether a giant planet will form and how it evolves. The number of exoplanet and disk observations is consistently rising, however, it is not yet possible to directly link these two populations. Therefore, a deep theoretical understanding of how planets form is crucial. We performed numerical simulations of planet formation via pebble and gas accretion, while including migration, in a viscously evolving protoplanetary disk, with dust growing, drifting, and evaporating at the ice lines. In our investigation of the most favorable conditions for giant planet formation, we find that these are high disk masses, early formation, and a large enough disk to host a long-lasting pebble flux. However, small disks with the same mass allow more efficient gas accretion onto planetary cores, leading to more massive gas giants. Given the right conditions, high viscosity leads to more massive cores and it enhances gas accretion. It also causes faster type II migration rates, so the giants have a decreasing final position for increasing viscosity. Intermediate dust fragmentation velocities provide the necessary pebble sizes and radial drift velocities for maximized pebble accretion and pebble flux. An enhanced dust-to-gas ratio can compensate for lower disk masses, but early formation is still crucial. We conclude that there is no specific initial parameter that leads to giant planet formation; rather, it is the outcome of a combination of complementary factors. This also implies that the diversity of the exoplanet systems is the product of the intrinsic diversity of the protoplanetary disks and it is crucial to take advantage of the increasing number and quality of observations to constrain the disk population properties and ultimately devise planet formation theories. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: Accepted by A&A, 15 pages, 6 figures

Journal ref: A&A 679, A42 (2023)

C3PO: Towards a complete census of co-moving pairs of stars. I. High precision stellar parameters for 250 stars

Authors: David Yong, Fan Liu, Yuan-Sen Ting, Meridith Joyce, Bertram Bitsch, Fei Dai, Aaron Dotter, Amanda I. Karakas, Michael T. Murphy

Abstract: We conduct a line-by-line differential analysis of a sample of 125 co-moving pairs of stars (dwarfs and subgiants near solar metallicity). We obtain high precision stellar parameters with average uncertainties in effective temperature, surface gravity and metallicity of 16.5 K, 0.033 dex and 0.014 dex, respectively. We classify the co-moving pairs of stars into two groups, chemically homogeneous (… ▽ More We conduct a line-by-line differential analysis of a sample of 125 co-moving pairs of stars (dwarfs and subgiants near solar metallicity). We obtain high precision stellar parameters with average uncertainties in effective temperature, surface gravity and metallicity of 16.5 K, 0.033 dex and 0.014 dex, respectively. We classify the co-moving pairs of stars into two groups, chemically homogeneous (conatal; |Delta[Fe/H]| $\le$ 0.04 dex) and inhomogeneous (non-conatal), and examine the fraction of chemically homogeneous pairs as a function of separation and effective temperature. The four main conclusions from this study are: (1) A spatial separation of \ds = 10$^6$ AUえーゆー is an approximate boundary between homogeneous and inhomogeneous pairs of stars, and we restrict our conclusions to only consider the 91 pairs with \ds $\le$ 10$^6$ AUえーゆー; (2) There is no trend between velocity separation and the fraction of chemically homogeneous pairs in the range \dv $\le$ 4 \kms; (3) We confirm that the fraction of chemically inhomogeneous pairs increases with increasing \teff\ and the trend matches a toy model of that expected from planet ingestion; (4) Atomic diffusion is not the main cause of the chemical inhomogeneity. A major outcome from this study is a sample of 56 bright co-moving pairs of stars with chemical abundance differences $\leq$ 0.02 dex (5\%) which is a level of chemical homogeneity comparable to that of the Hyades open cluster. These important objects can be used, in conjunction with star clusters and the \gaia\ ``benchmark'' stars, to calibrate stellar abundances from large-scale spectroscopic surveys. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: MNRAS in press (see source file for full versions of long tables)

Enriching inner discs and giant planets with heavy elements

Authors: Bertram Bitsch, Jingyi Mah

Abstract: Giant exoplanets seem to have on average a much larger heavy element content than the solar system giants. Past attempts to explain these heavy element contents include collisions between planets, accretion of volatile rich gas and accretion of gas enriched in micro-metre sized solids. However, these different theories individually could not explain the heavy element content of giants and the vola… ▽ More Giant exoplanets seem to have on average a much larger heavy element content than the solar system giants. Past attempts to explain these heavy element contents include collisions between planets, accretion of volatile rich gas and accretion of gas enriched in micro-metre sized solids. However, these different theories individually could not explain the heavy element content of giants and the volatile to refractory ratios in atmospheres of giant planets at the same time. Here we combine the approaches of gas accretion enhanced with vapor and small micro-meter sized dust grains. As pebbles drift inwards, the volatile component evaporates and enriches the disc, while the smaller silicate core of the pebble continues to move inwards. The smaller silicate pebbles drift slower, leading to a pile-up of material interior to the water ice line, increasing the dust-to-gas ratio interior to the ice line. Under the assumption that these small dust grains follow the motion of the gas, gas accreting giants accrete large fractions of small solids in addition to the volatile vapor. The effectiveness of the solid enrichment requires a large disc radius to maintain the pebble flux for a long time and a large viscosity that reduces the size and inward drift of the small dust grains. However, this process depends crucially on the debated size difference of the pebbles interior and exterior of the water ice line. On the other hand, the volatile component released by the inward drifting pebbles can lead to a large enrichment with heavy element vapor, independently of a size difference of pebbles interior and exterior to the water ice line. Our model stresses the importance of the disc's radius and viscosity on the enrichment of dust and vapor. Consequently we show how our model could explain the heavy element content of the majority of giant planets by using combined estimates of dust and vapor enrichment. △ Less

Submitted 2 October, 2023; v1 submitted 1 September, 2023; originally announced September 2023.

Comments: Accepted by A&A, 10 pages, 9 figures, now with language corrections

Journal ref: A&A 679, A11 (2023)

Close-in ice lines and the super-stellar C/O ratio in discs around very low-mass stars

Authors: J. Mah, B. Bitsch, I. Pascucci, T. Henning

Abstract: The origin of the elevated C/O ratios in discs around late M dwarfs compared to discs around solar-type stars is not well understood. Here we endeavour to reproduce the observed differences in the disc C/O ratios as a function of stellar mass using a viscosity-driven disc evolution model and study the corresponding atmospheric composition of planets that grow inside the water-ice line in these dis… ▽ More The origin of the elevated C/O ratios in discs around late M dwarfs compared to discs around solar-type stars is not well understood. Here we endeavour to reproduce the observed differences in the disc C/O ratios as a function of stellar mass using a viscosity-driven disc evolution model and study the corresponding atmospheric composition of planets that grow inside the water-ice line in these discs. We carried out simulations using a coupled disc evolution and planet formation code that includes pebble drift and evaporation. We used a chemical partitioning model for the dust composition in the disc midplane. Inside the water-ice line, the disc's C/O ratio initially decreases to sub-stellar due to the inward drift and evaporation of water-ice-rich pebbles before increasing again to super-stellar values due to the inward diffusion of carbon-rich vapour. We show that this process is more efficient for very low-mass stars compared to solar-type stars due to the closer-in ice lines and shorter disc viscous timescales. In high-viscosity discs, the transition from sub-stellar to super-stellar takes place faster due to the fast inward advection of carbon-rich gas. Our results suggest that planets accreting their atmospheres early (when the disc C/O is still sub-stellar) will have low atmospheric C/O ratios, while planets that accrete their atmospheres late (when the disc C/O has become super-stellar) can obtain high C/O ratios. Our model predictions are consistent with observations, under the assumption that all stars have the same metallicity and chemical composition, and that the vertical mixing timescales in the inner disc are much shorter than the radial advection timescales. This further strengthens the case for considering stellar abundances alongside disc evolution in future studies that aim to link planet (atmospheric) composition to disc composition. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 677, L7 (2023)

Do all gaps in protoplanetary discs host planets?

Authors: Anastasia Tzouvanou, Bertram Bitsch, Gabriele Pichierri

Abstract: Following the assumption that the disc substructures observed in protoplanetary discs originate from the interaction between the disc and the forming planets embedded therein, we aim to test if these putative planets could represent the progenitors of the currently observed giant exoplanets. We performed N-body simulations assuming initially three, four, five or seven planets. Our model includes… ▽ More Following the assumption that the disc substructures observed in protoplanetary discs originate from the interaction between the disc and the forming planets embedded therein, we aim to test if these putative planets could represent the progenitors of the currently observed giant exoplanets. We performed N-body simulations assuming initially three, four, five or seven planets. Our model includes pebble and gas accretion, migration, damping of eccentricities and inclinations, disc-planet interaction and disc evolution. We locate the planets in the positions where the gaps in protoplanetary discs have been observed and we evolve the systems for 100Myr including a few Myr of gas disc evolution, while also testing three values of $αあるふぁ$ viscosity. For planetary systems with initially three and four planets we find that most of the growing planets lie beyond the RV detection limit of 5AUえーゆー and only a small fraction of them migrate into the inner region. We also find that these systems have too low final eccentricities to be in agreement with the observed giant planet population. Systems initially consisting of five or seven planets become unstable after $\approx$40Kyr of integration time. This clearly shows that not every gap can host a planet. The general outcome of our simulations - too low eccentricities - is independent of the disc's viscosity and surface density. Further observations could either confirm the existence of an undetected population of wide-orbit giants or exclude the presence of such undetected population to constrain how many planets hide in gaps even further. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 677, A82 (2023)

How does accretion of planet-forming disks influence stellar abundances?

Authors: León-Alexander Hühn, Bertram Bitsch

Abstract: Millimeter sized dust grains experience radial velocities exceeding the gas velocities by orders of magnitude. The viscous evolution of the accretion disk adds disk material onto the central star's convective envelope, influencing its elemental abundances, [X/H]. At the same time, the envelope mass shrinks over time, amplifying the rate of abundance change. Therefore, the elemental abundances of t… ▽ More Millimeter sized dust grains experience radial velocities exceeding the gas velocities by orders of magnitude. The viscous evolution of the accretion disk adds disk material onto the central star's convective envelope, influencing its elemental abundances, [X/H]. At the same time, the envelope mass shrinks over time, amplifying the rate of abundance change. Therefore, the elemental abundances of the star are sensitive to disk processes that alter the composition and timing of disk accretion. We perform numerical 1D log-radial simulations integrating the disk advection-diffusion equation, accounting for phase transitions of chemical species at the evaporation fronts. They reveal a peak of refractory abundance within the first 2 Myr of $Δでるた\mathrm{[X/H]}\sim 5\times 10^{-2}$ if grain growth is significant, but subsequent accretion diminishes previous refractory abundance increases for long-lived disks. Planet formation can reduce the abundance of dust species whose evaporation fronts lie within the planet's orbit by opening a gap and consequently blocking inward drifting pebbles. We expect the accretion of the Solar protoplanetary disk with Jupiter present to have changed the Sun's elemental abundances by ${\sim}10^{-2}$ throughout its lifetime. These considerations are also applied to the HD106515 wide binary system. We find that measurements of $Δでるた\mathrm{[X/H]}$ are in reasonable agreement with results from simulations where the observed giant planet around HD106515 A is included and if HD106515B's disk formed planetesimals more efficiently. Simulations where the planet formed inside the water ice line are more favorable to agree with observations. Even though the general changes in the stellar abundances due to disk accretion are small, they are detectable at current sensitivities, indicating that the here presented methods can be used to constrain the planet formation pathway. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 23 pages, 17 figures, 2 tables. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 676, A87 (2023)

Giants are bullies: how their growth influences systems of inner sub-Neptunes and super-Earths

Authors: Bertram Bitsch, Andre Izidoro

Abstract: Observations point to a correlation between outer giants and inner sub-Neptunes, unexplained by simulations so far. We utilize N-body simulations including pebble and gas accretion as well as planetary migration to investigate how the gas accretion rates influence the formation of systems of inner sub-Neptunes and outer gas giants as well as the eccentricity distribution of the outer giant planets… ▽ More Observations point to a correlation between outer giants and inner sub-Neptunes, unexplained by simulations so far. We utilize N-body simulations including pebble and gas accretion as well as planetary migration to investigate how the gas accretion rates influence the formation of systems of inner sub-Neptunes and outer gas giants as well as the eccentricity distribution of the outer giant planets. Less efficient envelope contraction rates allow a more efficient formation of systems with inner sub-Neptunes and outer giants. This is caused by the fact that the cores formed in the inner disc are too small to accrete large envelopes and only cores growing in the outer disc can become giants. As a result, instabilities between the outer giant planets do not necessarily destroy the inner systems of sub-Neptunes unlike simulations where giant planets can form closer in. Our simulations show that up to 50% of the systems of cold Jupiters could have inner sub-Neptunes, in agreement with observations. Our simulations show a good agreement with the eccentricity distribution of giants, even though we find a slight mismatch to the mass and semi-major axes distributions. Synthetic transit observations of the inner systems (r<0.7 AUえーゆー) reveal an excellent match to the Kepler observations, where our simulations match the period ratios of adjacent planet pairs. Thus, the breaking the chains model for super-Earth and sub-Neptune formation remains consistent with observations even when outer giant planets are present. However, simulations with outer giant planets produce more systems with mostly only one inner planet and with larger eccentricities, in contrast to simulations without outer giants. We thus predict that systems with truly single close-in planets are more likely to host outer gas giants and we consequently suggest RV follow-up observations of these systems to constrain the formation pathway. △ Less

Submitted 11 May, 2023; v1 submitted 25 April, 2023; originally announced April 2023.

Comments: 21 pages, 17 figures, accepted for publication by A&A, now with language correction

Journal ref: A&A 674, A178 (2023)

Forming super-Mercuries: The role of stellar abundances

Authors: J. Mah, B. Bitsch

Abstract: Super-Mercuries, rocky exoplanets with bulk iron mass fraction of more than 60 per cent, appear to be preferentially hosted by stars with higher iron mass fraction than the Earth. It is unclear whether these iron-rich planets can form in the disc, or if giant impacts are necessary. Here we investigate the formation of super-Mercuries in their natal protoplanetary discs by taking into account their… ▽ More Super-Mercuries, rocky exoplanets with bulk iron mass fraction of more than 60 per cent, appear to be preferentially hosted by stars with higher iron mass fraction than the Earth. It is unclear whether these iron-rich planets can form in the disc, or if giant impacts are necessary. Here we investigate the formation of super-Mercuries in their natal protoplanetary discs by taking into account their host stars' abundances (Fe, Mg, Si, S). We employ a disc evolution model which includes the growth, drift, evaporation and recondensation of pebbles to compute the pebble iron mass fraction. The recondensation of outward-drifting iron vapour near the iron evaporation front is the key mechanism that facilitates an increase in the pebble iron mass fraction. We also simulate the growth of planetary seeds around the iron evaporation front using a planet formation model which includes pebble accretion and planet migration, and compute the final composition of the planets. Our simulations are able to reproduce the observed iron compositions of the super-Mercuries provided that all the iron in the disc are locked in pure Fe grains and that the disc viscosity is low. The combined effects of slow orbital migration of planets and long retention time of iron vapour in low-viscosity discs makes it easier to form iron-rich planets. Furthermore, we find that decreasing the stellar Mg/Si ratio results in an increase in the iron mass fraction of the planet due to a reduction in the abundance of Mg2SiO4, which has a very similar condensation temperature as iron, in the disc. Our results thus imply that super-Mercuries are more likely to form around stars with low Mg/Si, in agreement with observational data. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 9 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 673, A17 (2023)

Assessing the spin-orbit obliquity of low-mass planets in the breaking the chain formation model: A story of misalignment

Authors: Leandro Esteves, André Izidoro, Othon C. Winter, Bertram Bitsch, Andrea Isella

Abstract: The spin-orbit obliquity of a planetary system constraints its formation history. A large obliquity may either indicate a primordial misalignment between the star and its gaseous disk or reflect the effect of different mechanisms tilting planetary systems after formation. Observations and statistical analysis suggest that system of planets with sizes between 1 and 4 R$_{\oplus}$ have a wide range… ▽ More The spin-orbit obliquity of a planetary system constraints its formation history. A large obliquity may either indicate a primordial misalignment between the star and its gaseous disk or reflect the effect of different mechanisms tilting planetary systems after formation. Observations and statistical analysis suggest that system of planets with sizes between 1 and 4 R$_{\oplus}$ have a wide range of obliquities ($\sim0-30^{\circ}$), and that single- and multi-planet transiting have statistically indistinguishable obliquity distributions. Here, we revisit the ``breaking the chains'' formation model with focus in understanding the origin of spin-orbit obliquities. This model suggests that super-Earths and mini-Neptunes migrate close to their host stars via planet-disk gravitational interactions, forming chain of planets locked in mean-motion resonances. After gas-disk dispersal, about 90-99\% of these planetary systems experience dynamical instabilities, which spread the systems out. Using synthetic transit observations, we show that if planets are born in disks where the disk angular momentum is virtually aligned with the star's rotation spin, their final obliquity distributions peak at about $\sim$5 degrees or less, and the obliquity distributions of single and multi-planet transiting systems are statistically distinct. By treating the star-disk alignment as a free-parameter, we show that the obliquity distributions of single and multi-planet transiting systems only become statistically indistinguishable if planets are assumed to form in primordially misaligned natal disks with a ``tilt'' distribution peaking at $\gtrsim$10-20 deg. We discuss the origin of these misalignments in the context of star formation and potential implications of this scenario for formation models. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: 11 pages, 9 figures, 2 tables. Accepted for publication in MNRAS

MSC Class: 85-10 ACM Class: J.2

How the origin of stars in the Galaxy impacts the composition of planetary building blocks

Authors: N. Cabral, A. Guilbert-Lepoutre, B. Bitsch, N. Lagarde, S. Diakite

Abstract: Context. Our Galaxy is composed of different stellar populations with varying chemical abundances, which are thought to imprint the composition of planet building blocks (PBBs). As such, the properties of stars should affect the properties of planets and small bodies formed in their systems. In this context, high-resolution spectroscopic surveys open a window into the chemical links between and th… ▽ More Context. Our Galaxy is composed of different stellar populations with varying chemical abundances, which are thought to imprint the composition of planet building blocks (PBBs). As such, the properties of stars should affect the properties of planets and small bodies formed in their systems. In this context, high-resolution spectroscopic surveys open a window into the chemical links between and their host stars. Aims. We aim to determine the PBB composition trends for various stellar populations across the Galaxy by comparing the two large spectroscopic surveys APOGEE and GALAH. We assess the reliability of the PBB composition as determined with these surveys with a propagation error study. Methods. Stellar spectroscopic abundances from the large surveys GALAH-DR3 and APOGEE-DR17 were used as input with a stoichiometric condensation model. We classified stars into different Galactic components and we quantified the PBB composition trends as a function of [Fe/H]. We also analysed the distribution composition patterns in the [$αあるふぁ$/Fe]-[Fe/H] diagram. Results. Our propagation error study suggests that the overall trends with [Fe/H] and [$αあるふぁ$/Fe] are robust, which is supported by the double study of both APOGEE and GALAH. We therefore confirm the existence of a bimodal PBB composition separating the thin disc stars from the thick disc stars. Furthermore, we confirm that the stoichiometric water PBB content is anti-correlated with [Fe/H]. Conclusions. Our results imply that metal-poor stars both in the thin and thick disks are suitable hosts for water-rich PBBs and for ice-rich small bodies. However, for metal-poor stars ([Fe/H]<0), the PBBs around thick disc stars should have a higher water content than that around thin disc stars because of the $αあるふぁ$-contentdependence of the water mass fraction. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Journal ref: A&A 673, A117 (2023)

A recipe for orbital eccentricity damping in the type-I regime for low viscosity 2D-discs

Authors: Gabriele Pichierri, Bertram Bitsch, Elena Lega

Abstract: It is known that gap opening depends on the disc's viscosity; however, eccentricity damping formulas have only been derived at high viscosities, ignoring partial gap opening. We aim at obtaining a simple formula to model $e$-damping of the type-I regime in low viscosity discs, where even small planets may start opening partial. We perform high resolution 2D locally isothermal hydrodynamical simula… ▽ More It is known that gap opening depends on the disc's viscosity; however, eccentricity damping formulas have only been derived at high viscosities, ignoring partial gap opening. We aim at obtaining a simple formula to model $e$-damping of the type-I regime in low viscosity discs, where even small planets may start opening partial. We perform high resolution 2D locally isothermal hydrodynamical simulations of planets with varying masses on fixed orbits in discs with varying aspect ratios and viscosities. We determine the torque and power felt by the planet to derive migration and eccentricity damping timescales. We first find a lower limit to the gap depths below which vortices appear; this happens roughly at the transition between type-I and type-II regimes. For the simulations that remain stable, we obtain a fit to the observed gap depth in the limit of vanishing eccentricities that is similar to the one currently used in the literature but is accurate down to $αあるふぁ=3.16\times 10^{-5}$. We record the $e$-damping efficiency as a function of the observed gap depth and $e$: when the planet has opened a deep enough gap, a linear trend is observed independently of $e$; at shallower gaps this linear trend is preserved at low $e$, while it deviates to more efficient damping when $e$ is comparable to the disc's scale height. Both trends can be understood on theoretical grounds and are reproduced by a simple fitting formula. Our combined fits yield a simple recipe to implement type-I $e$-damping in $N$-body for partial gap opening planets that is consistent with high-resolution 2D hydro-simulations. The typical error of the fit is of the order of a few percent, and lower than the error of type-I torque formulas widely used in the literature. This will allow a more self-consistent treatment of planet-disc interactions of the type-I regime for population synthesis models at low viscosities. △ Less

Submitted 7 December, 2022; originally announced December 2022.

Comments: Submitted to A&A

Journal ref: A&A 670, A148 (2023)

Simultaneous gas accretion onto a pair of giant planets: Impact on their final mass and on the protoplanetary disk structure

Authors: Camille Bergez-Casalou, Bertram Bitsch, Sean N. Raymond

Abstract: Several planetary systems are known to host multiple giant planets. However, when two giant planets are accreting from the same disk, it is unclear what effect the presence of the second planet has on the gas accretion process of both planets. In this paper we perform long-term 2D isothermal hydrodynamical simulations (over more than 0.5 Myrs) with the FARGO-2D1D code, considering two non-migratin… ▽ More Several planetary systems are known to host multiple giant planets. However, when two giant planets are accreting from the same disk, it is unclear what effect the presence of the second planet has on the gas accretion process of both planets. In this paper we perform long-term 2D isothermal hydrodynamical simulations (over more than 0.5 Myrs) with the FARGO-2D1D code, considering two non-migrating planets accreting from the same gaseous disk. We find that the evolution of the planets' mass ratio depends on gap formation. However, in all cases, when the planets start accreting at the same time, they end up with very similar masses (0.9 $<m_{p,out}/m_{p,in}<$ 1.1 after 0.5 Myrs). Delaying the onset of accretion of one planet allows the planets' mass ratio to reach larger values initially, but they quickly converge to similar masses afterward (0.8 $<m_{p,out}/m_{p,in}<$ 2 in $10^5$ yrs). In order to reproduce the more diverse observed mass ratios of exoplanets, the planets must start accreting gas at different times, and their accretion must be stopped quickly after the beginning of runaway gas accretion (less than 0.5 Myrs), for example via disk dispersal. The evolution of the planets' mass ratio can have an important impact on the dynamics of the system and may constrain the formation history of Jupiter and Saturn. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted to A&A

Journal ref: A&A 669, A129 (2023)

The Exoplanet Radius Valley from Gas-driven Planet Migration and Breaking of Resonant Chains

Authors: Andre Izidoro, Hilke E. Schlichting, Andrea Isella, Rajdeep Dasgupta, Christian Zimmermann, Bertram Bitsch

Abstract: The size frequency distribution of exoplanet radii between 1 and 4$R_{\oplus}$ is bimodal with peaks at $\sim$1.4 $R_{\oplus}$ and $\sim$2.4 $R_{\oplus}$, and a valley at $\sim$1.8$R_{\oplus}$. This radius valley separates two classes of planets -- usually referred to as "super-Earths" and "mini-Neptunes" -- and its origin remains debated. One model proposes that super-Earths are the outcome of ph… ▽ More The size frequency distribution of exoplanet radii between 1 and 4$R_{\oplus}$ is bimodal with peaks at $\sim$1.4 $R_{\oplus}$ and $\sim$2.4 $R_{\oplus}$, and a valley at $\sim$1.8$R_{\oplus}$. This radius valley separates two classes of planets -- usually referred to as "super-Earths" and "mini-Neptunes" -- and its origin remains debated. One model proposes that super-Earths are the outcome of photo-evaporation or core-powered mass-loss stripping the primordial atmospheres of the mini-Neptunes. A contrasting model interprets the radius valley as a dichotomy in the bulk compositions, where super-Earths are rocky planets and mini-Neptunes are water-ice rich worlds. In this work, we test whether the migration model is consistent with the radius valley and how it distinguishes these views. In the migration model, planets migrate towards the disk inner edge forming a chain of planets locked in resonant configurations. After the gas disk dispersal, orbital instabilities "break the chains" and promote late collisions. This model broadly matches the period-ratio and planet-multiplicity distributions of Kepler planets, and accounts for resonant chains such as TRAPPIST-1, Kepler-223, and TOI-178. Here, by combining the outcome of planet formation simulations with compositional mass-radius relationships, and assuming complete loss of primordial H-rich atmospheres in late giant-impacts, we show that the migration model accounts for the exoplanet radius valley and the intra-system uniformity ("peas-in-a-pod") of Kepler planets. Our results suggest that planets with sizes of $\sim$1.4 $R_{\oplus}$ are mostly rocky, whereas those with sizes of $\sim$2.4 $R_{\oplus}$ are mostly water-ice rich worlds. Our results do not support an exclusively rocky composition for the cores of mini-Neptunes. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: Accepted for publicationn in ApJL

How drifting and evaporating pebbles shape giant planets III: The formation of WASP-77A b and $τたう$ Boötis b

Authors: Bertram Bitsch, Aaron David Schneider, Laura Kreidberg

Abstract: Atmospheric abundances are thought to constrain the planet formation pathway, because different species evaporate at different temperatures leaving distinct signatures in the accreted atmosphere. The planetary C/O ratio is thought to constrain the planet formation pathway, because of the condensation sequence of H$_2$O, CO$_2$, CH$_4$, and CO, resulting in an increase of the gas phase C/O ratio wi… ▽ More Atmospheric abundances are thought to constrain the planet formation pathway, because different species evaporate at different temperatures leaving distinct signatures in the accreted atmosphere. The planetary C/O ratio is thought to constrain the planet formation pathway, because of the condensation sequence of H$_2$O, CO$_2$, CH$_4$, and CO, resulting in an increase of the gas phase C/O ratio with increasing distance. Here we use a disc evolution model including pebble growth, drift and evaporation coupled with a planet formation model that includes pebble and gas accretion as well as planet migration to compute the atmospheric compositions of giant planets. We compare our results to the recent observations of the hot Jupiters WASP-77A b and $τたう$ Boötis b, which feature sub-solar and super-solar C/H and O/H values, respectively. Our simulations reproduce these measurements and show that giants like WASP-77A b should start to form beyond the CO$_2$ evaporation front, while giants like $τたう$ Boötis b should originate from beyond the H$_2$O line. Our model allows the formation of sub- and super-solar atmospheric compositions. However simulations without pebble evaporation can not reproduce the super-solar C/H and O/H ratios of $τたう$ Boötis b's atmosphere without additional accretion of solids. We identify the $αあるふぁ$ viscosity parameter of the disc as a key ingredient, because the viscosity drives the inward motion of volatile enriched vapor, responsible for the accretion of gaseous C and O. Depending on the planet's migration history order-of-magnitude differences in atmospheric C/H and O/H are expected. Our simulations also predict super-solar N/H for $τたう$ Boötis b and solar N/H for WASP-77A b. We conclude that pebble evaporation is a key ingredient to explain the variety of exoplanet atmospheres, because it can explain both, sub- and super-solar atmospheric abundances. △ Less

Submitted 14 July, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: Accepted by A&A, 7 pages, 4 figures

Journal ref: A&A 665, A138 (2022)

A new third planet and the dynamical architecture of the HD33142 planetary system

Authors: Trifon Trifonov, Anna Wollbold, Martin Kürster, Jan Eberhardt, Stephan Stock, Thomas Henning, Sabine Reffert, R. Paul Butler, Steven S. Vogt, Ansgar Reiners, Man Hoi Lee, Bertram Bitsch, Mathias Zechmeister, Florian Rodler, Volker Perdelwitz, Lev Tal-Or, Jan Rybizki, Paul Heeren, Davide Gandolfi, Oscar Barragán, Olga Zakhozhay, Paula Sarkis, Marcelo Tala Pinto, Diana Kossakowski, Vera Wolthoff , et al. (2 additional authors not shown)

Abstract: Based on recently-taken and archival HARPS, FEROS and HIRES radial velocities (RVs), we present evidence for a new planet orbiting the first ascent red giant star HD33142 (with an improved mass estimate of 1.52$\pm$0.03 M$_\odot$), already known to host two planets. We confirm the Jovian mass planets HD33142 b and c with periods of $P_{\rm b}$ = 330.0$_{-0.4}^{+0.4}$ d and $P_{\rm c}$ = 810.2… ▽ More Based on recently-taken and archival HARPS, FEROS and HIRES radial velocities (RVs), we present evidence for a new planet orbiting the first ascent red giant star HD33142 (with an improved mass estimate of 1.52$\pm$0.03 M$_\odot$), already known to host two planets. We confirm the Jovian mass planets HD33142 b and c with periods of $P_{\rm b}$ = 330.0$_{-0.4}^{+0.4}$ d and $P_{\rm c}$ = 810.2$_{-4.2}^{+3.8}$ d and minimum dynamical masses of $m_{\rm b}\sin{i}$ = 1.26$_{-0.05}^{+0.05}$ M$_{\rm Jup}$ and $m_{\rm c}\sin{i}$ = 0.89$_{-0.05}^{+0.06}$ M$_{\rm Jup}$. Furthermore, our periodogram analysis of the precise RVs shows strong evidence for a short-period Doppler signal in the residuals of a two-planet Keplerian fit, which we interpret as a third, Saturn-mass planet with $m_\mathrm{d}\sin{i}$ = 0.20$_{-0.03}^{+0.02}$ M$_{\rm Jup}$ on a close-in orbit with an orbital period of $P_{\rm d}$ =89.9$_{-0.1}^{+0.1}$ d. We study the dynamical behavior of the three-planet system configurations with an N-body integration scheme, finding it long-term stable with the planets alternating between low and moderate eccentricities episodes. We also performed N-body simulations, including stellar evolution and second-order dynamical effects such as planet-stellar tides and stellar mass-loss on the way to the white dwarf phase. We find that planets HD33142 b, c and d are likely to be engulfed near the tip of the red giant branch phase due to tidal migration. These results make the HD33142 system an essential benchmark for the planet population statistics of the multiple-planet systems found around evolved stars. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: Accepted for publication in AJ

Interpreting the atmospheric composition of exoplanets: sensitivity to planet formation assumptions

Authors: Paul Mollière, Tamara Molyarova, Bertram Bitsch, Thomas Henning, Aaron Schneider, Laura Kreidberg, Christian Eistrup, Remo Burn, Evert Nasedkin, Dmitry Semenov, Christoph Mordasini, Martin Schlecker, Kamber R. Schwarz, Sylvestre Lacour, Mathias Nowak, Matthäus Schulik

Abstract: Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric comp… ▽ More Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric compositions into planetary formation inferences. In this study we summarize the complexities and uncertainties of state-of-the-art planet formation models and how they influence planetary atmospheric compositions. We introduce a methodology that explores the effect of different formation model assumptions when interpreting atmospheric compositions. We apply this framework to the directly imaged planet HR 8799e. Based on its atmospheric composition, this planet may have migrated significantly during its formation. We show that including the chemical evolution of the protoplanetary disk leads to a reduced need for migration. Moreover, we find that pebble accretion can reproduce the planet's composition, but some of our tested setups lead to too low atmospheric metallicities, even when considering that evaporating pebbles may enrich the disk gas. We conclude that the definitive inversion from atmospheric abundances to planet formation for a given planet may be challenging, but a qualitative understanding of the effects of different formation models is possible, opening up pathways for new investigations. △ Less

Submitted 28 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ, base formation inversion code is available at: https://gitlab.com/mauricemolli/formation-inversion

Emerging population of gap-opening planets around type-A stars -- Long-term evolution of the forming planets around HD 163296

Authors: Johanna Müller-Horn, Gabriele Pichierri, Bertram Bitsch

Abstract: The presence of forming planets embedded in their protoplanetary disks has been inferred from the detection of multiring structures in such disks. Most of these suspected planets are undetectable by direct imaging observations at current measurement sensitivities. Inward migration and accretion might make these putative planets accessible to the Doppler method, but the actual extent of growth and… ▽ More The presence of forming planets embedded in their protoplanetary disks has been inferred from the detection of multiring structures in such disks. Most of these suspected planets are undetectable by direct imaging observations at current measurement sensitivities. Inward migration and accretion might make these putative planets accessible to the Doppler method, but the actual extent of growth and orbital evolution remains unconstrained. Under the premise that the gaps in the disk around HD 163296 originate from new-born planets, we investigate if and under which circumstances the gap-opening planets could represent progenitors of the exoplanet population detected around A-type stars. In particular, we study the dependence of final planetary masses and orbital parameters on the viscosity of the disk. The evolution of the embedded planets was simulated throughout the disk lifetime and up to 100 Myr after the dispersal of the disk, taking the evolving disk structure and a likely range of disk lifetimes into account. We find that the final configuration of the planets is largely determined by the $αあるふぁ$ viscosity parameter of the disk and less dependent on the choice for the disk lifetime and the initial planetary parameters. If we assume that planets such as those in HD 163296 evolve to form the observed exoplanet population of A-type stars, a $αあるふぁ$ parameter on the order of $3.16 \times 10^{-4} \lesssim αあるふぁ\lesssim 10^{-3}$ is required for the disks to induce sufficiently high migration rates. Depending on whether or not future direct imaging surveys will uncover a larger number of planets with $m_\mathrm{pl} \lesssim 3 M_\mathrm{Jup}$ and $a_\mathrm{pl} \gtrsim 10 \mathrm{AUえーゆー}$ we expect the $αあるふぁ$ parameter to be at the lower or upper end of this range, always under the assumption that such disks indeed harbor wide orbit planets. △ Less

Submitted 23 April, 2022; originally announced April 2022.

Comments: Accepted for publication in Astronomy and Astrophysics, 13 pages, 7 figures

Journal ref: A&A 663, A163 (2022)

Planet Formation Theory in the Era of ALMA and Kepler: from Pebbles to Exoplanets

Authors: Joanna Drazkowska, Bertram Bitsch, Michiel Lambrechts, Gijs D. Mulders, Daniel Harsono, Allona Vazan, Beibei Liu, Chris W. Ormel, Katherine Kretke, Alessandro Morbidelli

Abstract: Our understanding of planet formation has been rapidly evolving in recent years. The classical planet formation theory, developed when the only known planetary system was our own Solar System, has been revised to account for the observed diversity of the exoplanetary systems. At the same time, the increasing observational capabilities of the young stars and their surrounding disks bring new constr… ▽ More Our understanding of planet formation has been rapidly evolving in recent years. The classical planet formation theory, developed when the only known planetary system was our own Solar System, has been revised to account for the observed diversity of the exoplanetary systems. At the same time, the increasing observational capabilities of the young stars and their surrounding disks bring new constraints on the planet formation process. In this chapter, we summarize the new information derived from the exoplanets population and the circumstellar disks observations. We present the new developments in planet formation theory, from dust evolution to the growth of planetary cores by accretion of planetesimals, pebbles, and gas. We review the state-of-the-art models for the formation of diverse planetary systems, including the population synthesis approach which is necessary to compare theoretical model outcomes to the exoplanet population. We emphasize that the planet formation process may not be spatially uniform in the disk and there are preferential locations for the formation of planetesimals and planets. Outside of these locations, a significant fraction of solids is not growing past the pebble-sizes. The reservoir of pebbles plays an important role in the growth of planetary cores in the pebble accretion process. The timescale of the emergence of massive planetary cores is an important aspect of the present models and it is likely that the cores within one disk form at different times. In addition, there is growing evidence that the first planetary cores start forming early, during the circumstellar disk buildup process. △ Less

Submitted 24 May, 2023; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 38 pages, 13 figures, review chapter accepted to Protostars and Planets VII, editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura, corrected label in Fig. 1

Shallowness of circulation in hot Jupiters -- Advancing the Ohmic dissipation model

Authors: Henrik Knierim, Konstantin Batygin, Bertram Bitsch

Abstract: The inflated radii of giant short-period extrasolar planets collectively indicate that the interiors of hot Jupiters are heated by some anomalous energy dissipation mechanism. Although a variety of physical processes have been proposed to explain this heating, recent statistical evidence points to the confirmation of explicit predictions of the Ohmic dissipation theory, elevating this mechanism as… ▽ More The inflated radii of giant short-period extrasolar planets collectively indicate that the interiors of hot Jupiters are heated by some anomalous energy dissipation mechanism. Although a variety of physical processes have been proposed to explain this heating, recent statistical evidence points to the confirmation of explicit predictions of the Ohmic dissipation theory, elevating this mechanism as the most promising candidate for resolving the radius inflation problem. In this work, we present an analytic model for the dissipation rate and derive a simple scaling law that links the magnitude of energy dissipation to the thickness of the atmospheric weather layer. From this relation, we find that the penetration depth influences the Ohmic dissipation rate by an order of magnitude. We further investigate the weather layer depth of hot Jupiters from the extent of their inflation and show that, depending on the magnetic field strength, hot Jupiter radii can be maintained even if the circulation layer is relatively shallow. Additionally, we explore the evolution of zonal wind velocities with equilibrium temperature by matching our analytic model to statistically expected dissipation rates. From this analysis, we deduce that the wind speed scales approximately as $1/\sqrt{T_\mathrm{eq}-T_0}$, where $T_0$ is a constant that equals $T_0 \sim 1000~\mathrm{K}-1800~\mathrm{K}$ depending on planet-specific parameters (radius, mass, etc.). This work outlines inter-related constraints on the atmospheric flow and the magnetic field of hot Jupiters and provides a foundation for future work on the Ohmic heating mechanism. △ Less

Submitted 10 February, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: published in A&A, 5 pages, including 5 figures, plus 2 pages appendix, including 2 figures; v2 follows the A&A language guidelines

Journal ref: A&A 658, L7 (2022)

Constraining giant planet formation with synthetic ALMA images of the Solar System's natal protoplanetary disk

Authors: C. Bergez-Casalou, B. Bitsch, N. T. Kurtovic, P. Pinilla

Abstract: New ALMA observations of protoplanetary disks allow us to probe planet formation in other systems, giving us new constraints on planet formation processes. Meanwhile, studies of our own Solar System rely on constraints derived in a completely different way. However, it is still unclear what features the Solar System's disk could have produced during its gas phase. By running 2D isothermal hydro-si… ▽ More New ALMA observations of protoplanetary disks allow us to probe planet formation in other systems, giving us new constraints on planet formation processes. Meanwhile, studies of our own Solar System rely on constraints derived in a completely different way. However, it is still unclear what features the Solar System's disk could have produced during its gas phase. By running 2D isothermal hydro-simulations and a dust evolution model, we derive synthetic images at 1.3 mm wavelength using the radiative transfer code RADMC3D. We find that the embedded multiple giant planets strongly perturb the radial gas velocities of the disk, creating traffic jams in the dust. They produce over-densities different from the ones created by pressure traps and located away from the planets' positions in the disk. By deriving the images at 1.3mm from these dust distributions, we show that the traffic jams, observable with a high resolution, further blur the link between the number of gaps and rings in disks and the number of embedded planets. We additionally show that a system of 3 compact giant planets does not automatically produce bright outer rings at large radii in the disk. This means that high resolution observations of disks of various sizes are needed to distinguish between different giant planet formation scenarios during the disk phase, where the giants form either in the outer regions of the disks or in the inner regions. Finally, we find that, even when the dust temperature is determined self-consistently, the dust masses derived observationally might be off by up to a factor of ten compared to the dust contained in our simulations due to the creation of optically thick regions. Our study clearly shows that in addition to the constraints from exoplanets and the Solar System, ALMA has the power to constrain different stages of planet formation already during the first few million years. △ Less

Submitted 20 January, 2022; v1 submitted 10 January, 2022; originally announced January 2022.

Comments: Accepted to A&A

Journal ref: A&A 659, A6 (2022)

Planetesimal rings as the cause of the Solar System's planetary architecture

Authors: Andre Izidoro, Rajdeep Dasgupta, Sean N. Raymond, Rogerio Deienno, Bertram Bitsch, Andrea Isella

Abstract: Astronomical observations reveal that protoplanetary disks around young stars commonly have ring- and gap-like structures in their dust distributions. These features are associated with pressure bumps trapping dust particles at specific locations, which simulations show are ideal sites for planetesimal formation. Here we show that our Solar System may have formed from rings of planetesimals -- cre… ▽ More Astronomical observations reveal that protoplanetary disks around young stars commonly have ring- and gap-like structures in their dust distributions. These features are associated with pressure bumps trapping dust particles at specific locations, which simulations show are ideal sites for planetesimal formation. Here we show that our Solar System may have formed from rings of planetesimals -- created by pressure bumps -- rather than a continuous disk. We model the gaseous disk phase assuming the existence of pressure bumps near the silicate sublimation line (at $T \sim$1400~K), water snowline (at $T \sim$170~K), and CO-snowline (at $T \sim$30~K). Our simulations show that dust piles up at the bumps and forms up to three rings of planetesimals: a narrow ring near 1~au, a wide ring between $\sim$3-4~au and $\sim$10-20~au, and a distant ring between $\sim$20~au and $\sim$45~au. We use a series of simulations to follow the evolution of the innermost ring and show how it can explain the orbital structure of the inner Solar System and provides a framework to explain the origins of isotopic signatures of Earth, Mars and different classes of meteorites. The central ring contains enough mass to explain the rapid growth of the giant planets' cores. The outermost ring is consistent with dynamical models of Solar System evolution proposing that the early Solar System had a primordial planetesimal disk beyond the current orbit of Uranus. △ Less

Submitted 31 December, 2021; originally announced December 2021.

Comments: Published in Nature Astronomy on Dec 30, 2021. Authors' version including Methods and Supplementary Information

Composition of super-Earths, super-Mercuries, and their host stars

Authors: V. Adibekyan, N. C. Santos, C. Dorn, S. G. Sousa, A. A. Hakobyan, B. Bitsch, Ch. Mordasini, S. C. C. Barros, E. Delgado Mena, O. D. S. Demangeon, J. P. Faria, P. Figueira, B. M. T. B. Soares, G. Israelian

Abstract: Because of their common origin, it was assumed that the composition of planet building blocks should, to a first order, correlate with stellar atmospheric composition, especially for refractory elements. In fact, information on the relative abundance of refractory and major rock-forming elements such as Fe, Mg, Si has been commonly used to improve interior estimates for terrestrial planets. Recent… ▽ More Because of their common origin, it was assumed that the composition of planet building blocks should, to a first order, correlate with stellar atmospheric composition, especially for refractory elements. In fact, information on the relative abundance of refractory and major rock-forming elements such as Fe, Mg, Si has been commonly used to improve interior estimates for terrestrial planets. Recently Adibekyan et al. (2021) presented evidence of a tight chemical link between rocky planets and their host stars. In this study we add six recently discovered exoplanets to the sample of Adibekyan et al and re-evaluate their findings in light of these new data. We confirm that i) iron-mass fraction of rocky exoplanets correlates (but not a 1:1 relationship) with the composition of their host stars, ii) on average the iron-mass fraction of planets is higher than that of the primordial iron-mass fraction of the protoplanetary disk, iii) super-Mercuries are formed in disks with high iron content. Based on these results we conclude that disk-chemistry and planet formation processes play an important role in the composition, formation, and evolution of super-Earths and super-Mercuries. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Comments: Peer-reviewed conference (Astronomy in the Crossroads of Interdisciplinary and Multidisciplinary Sciences) proceeding

How the planetary eccentricity influences the pebble isolation mass

Authors: Raúl O. Chametla, Frédéric S. Masset, Clément Baruteau, Bertram Bitsch

Abstract: We investigate the pebble isolation mass for a planet on a fixed eccentric orbit in its protoplanetary disc by conducting a set of 2D hydrodynamical simulations including dust turbulent diffusion. A range of planet eccentricities up to $e=0.2$ is adopted. Our simulations also cover a range of $αあるふぁ-$turbulent viscosities, and for each pair $\{αあるふぁ,e\}$ the pebble isolation mass is estimated as the minim… ▽ More We investigate the pebble isolation mass for a planet on a fixed eccentric orbit in its protoplanetary disc by conducting a set of 2D hydrodynamical simulations including dust turbulent diffusion. A range of planet eccentricities up to $e=0.2$ is adopted. Our simulations also cover a range of $αあるふぁ-$turbulent viscosities, and for each pair $\{αあるふぁ,e\}$ the pebble isolation mass is estimated as the minimum planet mass in our simulations such that solids with a Stokes number $\gtrsim 0.05$ do not flow across the planet orbit and remain trapped around a pressure bump outside the planet gap. For $αあるふぁ<10^{-3}$, we find that eccentric planets reach a well-defined pebble isolation mass, which can be smaller than for planets on circular orbits when the eccentricity remains smaller than the disc's aspect ratio. We provide a fitting formula for how the pebble isolation mass depends on planet eccentricity. However, for $αあるふぁ> 10^{-3}$, eccentric planets cannot fully stall the pebbles flow, and thus do not reach a well-defined pebble isolation mass. Our results suggest that the maximum mass reached by rocky cores should exhibit a dichotomy depending on the disc turbulent viscosity. While being limited to ${\cal O}(10\,M_\oplus)$ in low-viscosity discs, this maximum mass could reach much larger values in discs with a high turbulent viscosity in the planet vicinity. Our results further highlight that pebble filtering by growing planets might not be as effective as previously thought, especially in high-viscosity discs, with important implications to protoplanetary discs observations. △ Less

Submitted 28 December, 2021; originally announced December 2021.

Comments: 9 pages, 7 figures. Accepted for publication in MNRAS

The "Breaking The Chains" migration model for super-Earths formation: the effect of collisional fragmentation

Authors: Leandro Esteves, André Izidoro, Bertram Bitsch, Seth A. Jacobson, Sean N. Raymond, Rogerio Deienno, Othon C. Winter

Abstract: Planets between 1-4 Earth radii with orbital periods <100 days are strikingly common. The migration model proposes that super-Earths migrate inwards and pile up at the disk inner edge in chains of mean motion resonances. After gas disk dispersal, simulations show that super-Earth's gravitational interactions can naturally break their resonant configuration leading to a late phase of giant impacts.… ▽ More Planets between 1-4 Earth radii with orbital periods <100 days are strikingly common. The migration model proposes that super-Earths migrate inwards and pile up at the disk inner edge in chains of mean motion resonances. After gas disk dispersal, simulations show that super-Earth's gravitational interactions can naturally break their resonant configuration leading to a late phase of giant impacts. The instability phase is key to matching the orbital spacing of observed systems. Yet, most previous simulations have modelled collisions as perfect accretion events, ignoring fragmentation. In this work, we investigate the impact of imperfect accretion on the breaking the chains scenario. We performed N-body simulations starting from distributions of planetary embryos and modelling the effects of pebble accretion and migration in the gas disk. Our simulations also follow the long-term dynamical evolution of super-Earths after the gas disk dissipation. We compared the results of simulations where collisions are treated as perfect merging events with those where imperfect accretion and fragmentation are allowed. We concluded that the perfect accretion is a suitable approximation in this regime, from a dynamical point of view. Although fragmentation events are common, only ~10% of the system mass is fragmented during a typical "late instability phase", with fragments being mostly reacreted by surviving planets. This limited total mass in fragments proved to be insufficient to alter qualitatively the final system dynamical configuration -- e.g. promote strong dynamical friction or residual migration -- compared to simulations where fragmentation is neglected. △ Less

Submitted 29 October, 2021; originally announced November 2021.

Comments: 13 pages, 8 figures and 1 table. Accepted for publication in MNRAS

MSC Class: 85-10

The Kepler-223 resonance holds information on turbulence during the gas disk phase

Authors: L. -A. Hühn, G. Pichierri, B. Bitsch, K. Batygin

Abstract: Many fundamental physical processes regarding planetary formation in protoplanetary disks are still imperfectly understood, with an elusive phenomenon being turbulence in such disks. Observations are available of planetary systems and of some protoplanetary disks, which can serve as a starting point for these investigations. Detected systems reveal different architectures of planets. One particula… ▽ More Many fundamental physical processes regarding planetary formation in protoplanetary disks are still imperfectly understood, with an elusive phenomenon being turbulence in such disks. Observations are available of planetary systems and of some protoplanetary disks, which can serve as a starting point for these investigations. Detected systems reveal different architectures of planets. One particularly interesting case to consider is the Kepler-223 system, which contains a rare configuration of four planets in a resonance chain, implying a certain migration history. We aim to use the orbital configuration of Kepler-223's planets to constrain the parameters of the protoplanetary disk that allow for the formation of a chain of mean-motion resonances that resembles the one of Kepler-223. The parameters we aim to investigate are primarily the disk viscosity and surface density. We use the swift_symba N-body integrator with additional dissipative forces to mimic planet-disk interactions. We constrain the surface densities and viscosities that allow the formation of a resonant chain like that of Kepler-223. We find that surface densities of up to a few minimum mass solar nebula (MMSN) surface densities and disk viscosity parameters $αあるふぁ$ of a few $10^{-3}$ up to $10^{-2}$ are most successful at reproducing the architecture of this particular planetary system. We describe how these two quantities are linked to each other considering the success of reproducing the chain and find that higher disk surface densities in turn require lower viscosities to build the chain. Our results show that well characterized observed planetary systems hold information about their formation conditions in the protoplanetary disks and that it is possible to extract this information, namely the initial disk surface density and viscosity, helping to constrain planet formation. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 12 pages, 7 figues, accepted for publication in A&A

Journal ref: A&A 656, A115 (2021)

How drifting and evaporating pebbles shape giant planets II: Volatiles and refractories in atmospheres

Authors: Aaron David Schneider, Bertram Bitsch

Abstract: Upcoming studies of extrasolar gas giants will give precise insights into the composition of planetary atmospheres with the ultimate goal to link it to the formation history of the planet. Here, we investigate how drifting and evaporating pebbles that enrich the gas phase of the disk influence the chemical composition of growing and migrating gas giants. To achieve this goal, we perform semi analy… ▽ More Upcoming studies of extrasolar gas giants will give precise insights into the composition of planetary atmospheres with the ultimate goal to link it to the formation history of the planet. Here, we investigate how drifting and evaporating pebbles that enrich the gas phase of the disk influence the chemical composition of growing and migrating gas giants. To achieve this goal, we perform semi analytical 1D models of protoplanetary disks including viscous evolution, pebble drift and evaporation to simulate the growth of planets from planetary embryos to Jupiter mass objects by the accretion of pebbles and gas while they migrate through the disk. The gas phase of the protoplanetary disk is enriched due to the evaporation of inward drifting pebbles crossing evaporation lines, leading to the accretion of large amounts of volatiles into the planetary atmosphere. As a consequence, gas accreting planets are enriched in volatiles (C, O, N) compared to refractories (e.g., Mg, Si, Fe) by up to a factor of 100, depending on the chemical species, its exact abundance and volatility as well as the disk's viscosity. A simplified model for the formation of Jupiter reveals that its nitrogen content can be explained by inward diffusing nitrogen rich vapor, implying that Jupiter does not need to form close to the N2 evaporation front as indicated by previous simulations. However, our model predicts a too low oxygen abundance for Jupiter, implying either Jupiter's migration across the water ice line or an additional accretion of solids into the atmosphere, which can also increase Jupiter's carbon abundance. The accretion of solids will increase the refractory to volatile ratio in planetary atmospheres substantially. We thus conclude that the volatile to refractory ratio can place a strong constraint on planet formation theories making it an important target for future observations. △ Less

Submitted 15 September, 2021; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: accepted for publication in A&A

Journal ref: A&A 654, A72 (2021)

Detailed elemental abundances of binary stars: Searching for signatures of planet formation and atomic diffusion

Authors: Fan Liu, Bertram Bitsch, Martin Asplund, Bei-Bei Liu, Michael T. Murphy, David Yong, Yuan-Sen Ting, Sofia Feltzing

Abstract: Binary star systems are assumed to be co-natal and coeval, thus to have identical chemical composition. In this work we aim to test the hypothesis that there is a connection between observed element abundance patterns and the formation of planets using binary stars. Moreover, we also want to test how atomic diffusion might influence the observed abundance patterns. We conduct a strictly line-by-li… ▽ More Binary star systems are assumed to be co-natal and coeval, thus to have identical chemical composition. In this work we aim to test the hypothesis that there is a connection between observed element abundance patterns and the formation of planets using binary stars. Moreover, we also want to test how atomic diffusion might influence the observed abundance patterns. We conduct a strictly line-by-line differential chemical abundance analysis of 7 binary systems. Stellar atmospheric parameters and elemental abundances are obtained with extremely high precision (< 3.5%) using the high quality spectra from VLT/UVES and Keck/HIRES. We find that 4 of 7 binary systems show subtle abundance differences (0.01 - 0.03 dex) without clear correlations with the condensation temperature, including two planet-hosting pairs. The other 3 binary systems exhibit similar degree of abundance differences correlating with the condensation temperature. We do not find any clear relation between the abundance differences and the occurrence of known planets in our systems. Instead, the overall abundance offsets observed in the binary systems (4 of 7) could be due to the effects of atomic diffusion. Although giant planet formation does not necessarily imprint chemical signatures onto the host star, the differences in the observed abundance trends with condensation temperature, on the other hand, are likely associated with diverse histories of planet formation (e.g., formation location). Furthermore, we find a weak correlation between abundance differences and binary separation, which may provide a new constraint on the formation of binary systems. △ Less

Submitted 24 August, 2021; originally announced August 2021.

Comments: 15 pages, 16 figures; accepted by MNRAS

A pair of warm giant planets near the 2:1 mean motion resonance around the K-dwarf star TOI-2202

Authors: Trifon Trifonov, Rafael Brahm, Nestor Espinoza, Thomas Henning, Andrés Jordán, David Nesvorny, Rebekah I. Dawson, Jack J. Lissauer, Man Hoi Lee, Diana Kossakowski, Felipe I. Rojas, Melissa J. Hobson, Paula Sarkis, Martin Schlecker, Bertram Bitsch, Gaspar Á. Bakos, Mauro Barbieri, Waqas Bhatti, R. Paul Butler, Jeffrey D. Crane, Sangeetha Nandakumar, Matías R. Díaz, Stephen Shectman, Johanna Teske, Pascal Torres , et al. (15 additional authors not shown)

Abstract: TOI-2202 b is a transiting warm Jovian-mass planet with an orbital period of P=11.91 days identified from the Full Frame Images data of five different sectors of the TESS mission. Ten TESS transits of TOI-2202 b combined with three follow-up light curves obtained with the CHAT robotic telescope show strong transit timing variations (TTVs) with an amplitude of about 1.2 hours. Radial velocity follo… ▽ More TOI-2202 b is a transiting warm Jovian-mass planet with an orbital period of P=11.91 days identified from the Full Frame Images data of five different sectors of the TESS mission. Ten TESS transits of TOI-2202 b combined with three follow-up light curves obtained with the CHAT robotic telescope show strong transit timing variations (TTVs) with an amplitude of about 1.2 hours. Radial velocity follow-up with FEROS, HARPS and PFS confirms the planetary nature of the transiting candidate (a$_{\rm b}$ = 0.096 $\pm$ 0.002 au, m$_{\rm b}$ = 0.98 $\pm$ 0.06 M$_{\rm Jup}$), and dynamical analysis of RVs, transit data, and TTVs points to an outer Saturn-mass companion (a$_{\rm c}$ = 0.155 $\pm$ 0.003 au, m$_{\rm c}$= $0.37 \pm 0.10$ M$_{\rm Jup}$) near the 2:1 mean motion resonance. Our stellar modeling indicates that TOI-2202 is an early K-type star with a mass of 0.82 M$_\odot$, a radius of 0.79 R$_\odot$, and solar-like metallicity. The TOI-2202 system is very interesting because of the two warm Jovian-mass planets near the 2:1 MMR, which is a rare configuration, and their formation and dynamical evolution are still not well understood. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: Accepted for publication in AJ

How drifting and evaporating pebbles shape giant planets I: Heavy element content and atmospheric C/O

Authors: Aaron David Schneider, Bertram Bitsch

Abstract: Recent observations of extrasolar gas giants suggest super-stellar C/O ratios in planetary atmospheres, while interior models of observed extrasolar giant planets additionally suggest high heavy element contents. Furthermore, recent observations of protoplanetary disks revealed super-solar C/H ratios, which are explained by inward drifting and evaporating pebbles, enhancing the volatile content of… ▽ More Recent observations of extrasolar gas giants suggest super-stellar C/O ratios in planetary atmospheres, while interior models of observed extrasolar giant planets additionally suggest high heavy element contents. Furthermore, recent observations of protoplanetary disks revealed super-solar C/H ratios, which are explained by inward drifting and evaporating pebbles, enhancing the volatile content of the disk. We investigate how the inward drift and evaporation of volatile rich pebbles influences the atmospheric C/O ratio and heavy element content of giant planets growing by pebble and gas accretion. To achieve this goal, we perform semi analytical 1D models of protoplanetary disks including the treatment of viscous evolution and heating, pebble drift and simple chemistry to simulate the growth of planets from planetary embryos to Jupiter mass objects by accretion of pebbles and gas while they migrate through the disk. Our simulations show that the composition of the planetary gas atmosphere is dominated by the accretion of vapour, originating from inward drifting evaporating pebbles. This process allows the giant planets to harbour large heavy element contents. In addition, our model reveals that giant planets originating further away from the central star have a higher C/O ratio on average due to the evaporation of methane rich pebbles in the outer disk. These planets can then also harbour super-solar C/O ratios, in line with exoplanet observations. However, planets formed in the outer disk harbour a smaller heavy element content, due to a smaller vapour enrichment of the outer disk. Our model predicts that giant planets with low/large atmospheric C/O should harbour a large/low total heavy element content. We further conclude that the inclusion of pebble evaporation at evaporation lines is a key ingredient to determine the heavy element content and composition of giant planets. △ Less

Submitted 9 September, 2021; v1 submitted 27 May, 2021; originally announced May 2021.

Comments: accepted for publication in A&A

Journal ref: A&A 654, A71 (2021)

The effect of a strong pressure bump in the Sun's natal disk: Terrestrial planet formation via planetesimal accretion rather than pebble accretion

Authors: André Izidoro, Bertram Bitsch, Rajdeep Dasgupta

Abstract: Mass-independent isotopic anomalies of carbonaceous and non-carbonaceous meteorites show a clear dichotomy suggesting an efficient separation of the inner and outer solar system. Observations show that ring-like structures in the distribution of mm-sized pebbles in protoplanetary disks are common. These structures are often associated with drifting pebbles being trapped by local pressure maxima in… ▽ More Mass-independent isotopic anomalies of carbonaceous and non-carbonaceous meteorites show a clear dichotomy suggesting an efficient separation of the inner and outer solar system. Observations show that ring-like structures in the distribution of mm-sized pebbles in protoplanetary disks are common. These structures are often associated with drifting pebbles being trapped by local pressure maxima in the gas disk. Similar structures may also have existed in the sun's natal disk, which could naturally explain the meteorite/planetary isotopic dichotomy. Here, we test the effects of a strong pressure bump in the outer disk (e.g. $\sim$5~au) on the formation of the inner solar system. We model dust coagulation and evolution, planetesimal formation, as well as embryo's growth via planetesimal and pebble accretion. Our results show that terrestrial embryos formed via planetesimal accretion rather than pebble accretion. In our model, the radial drift of pebbles foster planetesimal formation. However, once a pressure bump forms, pebbles in the inner disk are lost via drift before they can be efficiently accreted by embryos growing at $\gtrapprox$1~au. Embryos inside $\sim$0.5-1.0au grow relatively faster and can accrete pebbles more efficiently. However, these same embryos grow to larger masses so they should migrate inwards substantially, which is inconsistent with the current solar system. Therefore, terrestrial planets most likely accreted from giant impacts of Moon to roughly Mars-mass planetary embryos formed around $\gtrapprox$1.0~au. Finally, our simulations produce a steep radial mass distribution of planetesimals in the terrestrial region which is qualitatively aligned with formation models suggesting that the asteroid belt was born low-mass. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: Accepted in ApJ

Dry or water world? How the water contents of inner sub-Neptunes constrain giant planet formation and the location of the water ice line

Authors: Bertram Bitsch, Sean N. Raymond, Lars A. Buchhave, Aaron Bello-Arufe, Alexander D. Rathcke, Aaron David Schneider

Abstract: In the pebble accretion scenario, the pebbles that form planets drift inward from the outer disk regions, carrying water ice with them. At the water ice line, the water ice on the inward drifting pebbles evaporates and is released into the gas phase, resulting in water-rich gas and dry pebbles that move into the inner disk regions. Large planetary cores can block the inward drifting pebbles by for… ▽ More In the pebble accretion scenario, the pebbles that form planets drift inward from the outer disk regions, carrying water ice with them. At the water ice line, the water ice on the inward drifting pebbles evaporates and is released into the gas phase, resulting in water-rich gas and dry pebbles that move into the inner disk regions. Large planetary cores can block the inward drifting pebbles by forming a pressure bump outside their orbit in the protoplanetary disk. Depending on the relative position of a growing planetary core relative to the water ice line, water-rich pebbles might be blocked outside or inside the water ice line. Pebbles blocked outside the water ice line do not evaporate and thus do not release their water vapor into the gas phase, resulting in a dry inner disk, while pebbles blocked inside the water ice line release their water vapor into the gas phase, resulting in water vapor diffusing into the inner disk. As a consequence, close-in sub-Neptunes that accrete some gas from the disk should be dry or wet, respectively, if outer gas giants are outside or inside the water ice line, assuming that giant planets form fast, as has been suggested for Jupiter in our Solar System. Alternatively, a sub-Neptune could form outside the water ice line, accreting a large amount of icy pebbles and then migrating inward as a very wet sub-Neptune. We suggest that the water content of inner sub-Neptunes in systems with giant planets that can efficiently block the inward drifting pebbles could constrain the formation conditions of these systems, thus making these sub-Neptunes exciting targets for detailed characterization (e.g., with JWST, ELT, or ARIEL). In addition, the search for giant planets in systems with already characterized sub-Neptunes can be used to constrain the formation conditions of giant planets as well. △ Less

Submitted 3 May, 2021; v1 submitted 23 April, 2021; originally announced April 2021.

Comments: Accepted by A&A Letters

Journal ref: A&A 649, L5 (2021)

The water-ice line as a birthplace of planets: Implications of a species-dependent dust fragmentation threshold

Authors: Jonas Müller, Sofia Savvidou, Bertram Bitsch

Abstract: The thermodynamic structure of protoplanetary discs is determined by dust opacities, which depend on the size of the dust grains and their chemical composition. In the inner regions, the grain sizes are regulated by the level of turbulence (e.g. $αあるふぁ$ viscosity) and by the dust fragmentation velocity that represents the maximal velocity that grains can have at a collision before they fragment. Here,… ▽ More The thermodynamic structure of protoplanetary discs is determined by dust opacities, which depend on the size of the dust grains and their chemical composition. In the inner regions, the grain sizes are regulated by the level of turbulence (e.g. $αあるふぁ$ viscosity) and by the dust fragmentation velocity that represents the maximal velocity that grains can have at a collision before they fragment. Here, we perform self-consistently calculated 2D hydrodynamical simulations that consider a full grain size distribution of dust grains with a transition in the dust fragmentation velocity at the water-ice line. This approach accounts for the results of previous particle collision laboratory experiments, in which silicate particles typically have a lower dust fragmentation velocity than water-ice particles. Furthermore, we probe the effects of variations in the water abundance, the dust-to-gas ratio, and the turbulence parameter on the disc structure. For the discs with a transition in the dust fragmentation velocity at the water-ice line, we find a narrow but striking zone of planetary outward migration, including for low viscosities. In addition, we find a bump in the radial pressure gradient profile that tends to be located slightly inside the ice line. Both of these features are present for all tested disc parameters. Thus, we conclude that the ice line can function both as a migration trap, which can extend the growth times of planets before they migrate to the inner edge of the protoplanetary disc, and as a pressure trap, where planetesimal formation can be initiated or enhanced. △ Less

Submitted 27 April, 2021; v1 submitted 14 April, 2021; originally announced April 2021.

Comments: Accepted by A&A, 13 pages, 17 figures

Journal ref: A&A 650, A185 (2021)

The growth of super-Earths: the importance of a self-consistent treatment of disc structures and pebble accretion

Authors: Sofia Savvidou, Bertram Bitsch

Abstract: The conditions in the protoplanetary disc are determinant for the various planet formation mechanisms. We present a framework which combines self-consistent disc structures with the calculations of the growth rates of planetary embryos via pebble accretion, in order to study the formation of Super-Earths. We first perform 2D hydrodynamical simulations of the inner discs, considering a grain size d… ▽ More The conditions in the protoplanetary disc are determinant for the various planet formation mechanisms. We present a framework which combines self-consistent disc structures with the calculations of the growth rates of planetary embryos via pebble accretion, in order to study the formation of Super-Earths. We first perform 2D hydrodynamical simulations of the inner discs, considering a grain size distribution with multiple chemical species and their corresponding size and composition dependent opacities. The resulting aspect ratios are almost constant with orbital distance, resulting in radially constant pebble isolation masses, the mass where pebble accretion stops. This supports the "peas-in-a-pod" constraint from the Kepler observations. The derived pebble sizes are used to calculate the growth rates of planetary embryos via pebble accretion. Discs with low levels of turbulence (expressed through the $αあるふぁ$-viscosity) and/or high dust fragmentation velocities allow larger particles, hence lead to smaller pebble isolation masses, and the contrary. At the same time, small pebble sizes lead to low accretion rates. We find that there is a trade-off between the pebble isolation mass and the growth timescale with the best set of parameters being an $αあるふぁ$-viscosity of $10^{-3}$ and a dust fragmentation velocity of 10 m/s, mainly for an initial gas surface density (at 1 AUえーゆー) greater than 1000 $\rm{g/cm^2}$. A self-consistent treatment between the disc structures and the pebble sizes is thus of crucial importance for planet formation simulations. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: Accepted by A&A, 15 pages, 13 figures

Journal ref: A&A 650, A132 (2021)

A compositional link between rocky exoplanets and their host stars

Authors: Vardan Adibekyan, Caroline Dorn, Sérgio G. Sousa, Nuno C. Santos, Bertram Bitsch, Garik Israelian, Christoph Mordasini, Susana C. C. Barros, Elisa Delgado Mena, Olivier D. S. Demangeon, João P. Faria, Pedro Figueira, Artur A. Hakobyan, Mahmoudreza Oshagh, Barbara M. T. B. Soares, Masanobu Kunitomo, Yoichi Takeda, Emiliano Jofré, Romina Petrucci, Eder Martioli

Abstract: Stars and planets both form by accreting material from a surrounding disk. Because they grow from the same material, theory predicts that there should be a relationship between their compositions. In this study, we search for a compositional link between rocky exoplanets and their host stars. We estimate the iron-mass fraction of rocky exoplanets from their masses and radii and compare it with the… ▽ More Stars and planets both form by accreting material from a surrounding disk. Because they grow from the same material, theory predicts that there should be a relationship between their compositions. In this study, we search for a compositional link between rocky exoplanets and their host stars. We estimate the iron-mass fraction of rocky exoplanets from their masses and radii and compare it with the compositions of their host stars, which we assume reflect the compositions of the protoplanetary disks. We find a correlation (but not a 1:1 relationship) between these two quantities, with a slope of >4, which we interpret as being attributable to planet formation processes. Super-Earths and super-Mercuries appear to be distinct populations with differing compositions, implying differences in their formation processes. △ Less

Submitted 15 October, 2021; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: Authors' version of the manuscript. Published in Science

Journal ref: Science, vol 374, issue 6565, pp. 330-332 (2021)

Influence of grain sizes and composition on the contraction rates of planetary envelopes and on planetary migration

Authors: Bertram Bitsch, Sofia Savvidou

Abstract: A crucial phase during planetary growth is the migration when the planetary core has been assembled, but the planet did not open a deep gap yet. During this phase the planet is subject to fast type-I migration, which is mostly directed inwards, and the planet can lose a significant fraction of its semi-major axis. The duration of this phase is set by how long the planetary envelope needs to contra… ▽ More A crucial phase during planetary growth is the migration when the planetary core has been assembled, but the planet did not open a deep gap yet. During this phase the planet is subject to fast type-I migration, which is mostly directed inwards, and the planet can lose a significant fraction of its semi-major axis. The duration of this phase is set by how long the planetary envelope needs to contract until it reaches a mass similar to the mass of the planetary core, which is when runaway gas accretion can set in and the planet can open a deeper gap in the disc, transitioning into the slower type-II migration. This envelope contraction phase depends crucially on the planetary mass and on the opacity inside the planetary envelope. Here we study how different opacity prescriptions influence the envelope contraction time and how this in turn influences how far the planet migrates through the disc. We find within our simulations that the size distribution of the grains as well as the chemical composition of the grains crucially influences how far the planet migrates before it can reach the runaway gas accretion phase. Grain size distributions with larger grain sizes result in less inward migration of the growing planet, due to faster gas accretion enabled by more efficient cooling. In addition, we find that planets forming in water poor environments can contract their envelope faster and thus migrate less, implying that gas giants forming in water poor environments might be located further away from their central star compared to gas giants forming in water rich environments. Future studies of planet formation that aim to investigate the chemical composition of formed gas giants need to take these effects self consistently into account. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: Accepted by A&A, 11 pages

Journal ref: A&A 647, A96 (2021)

Probing the impact of varied migration and gas accretion rates for the formation of giant planets in the pebble accretion scenario

Authors: Nelson Ndugu, Bertram Bitsch, Alessandro Morbidelli, Aurélien Crida, Edward Jurua

Abstract: The final orbital position of growing planets is determined by their migration speed, which is essentially set by the planetary mass. Small mass planets migrate in type I migration, while more massive planets migrate in type II migration, which is thought to depend mostly on the viscous evolution rate of the disc. A planet is most vulnerable to inward migration before it reaches type II migration… ▽ More The final orbital position of growing planets is determined by their migration speed, which is essentially set by the planetary mass. Small mass planets migrate in type I migration, while more massive planets migrate in type II migration, which is thought to depend mostly on the viscous evolution rate of the disc. A planet is most vulnerable to inward migration before it reaches type II migration and can lose a significant fraction of its semi-major axis at this stage. We investigated the influence of different disc viscosities, the dynamical torque and gas accretion from within the horseshoe region as mechanisms for slowing down planet migration. Our study confirms that planets growing in low viscosity environments migrate less, due to the earlier gap opening and slower type II migration rate. We find that taking the gas accretion from the horseshoe region into account allows an earlier gap opening and this results in less inward migration of growing planets. Furthermore, this effect increases the planetary mass compared to simulations that do not take the effect of gas accretion from the horseshoe region. Moreover, combining the effect of the dynamical torque with the effect of gas accretion from the horseshoe region, significantly slows down inward migration. Taking these effects into account could allow the formation of cold Jupiters (a > 1 au) closer to the water ice line region compared to previous simulations that did not take these effects into account. We thus conclude that gas accretion from within the horseshoe region and the dynamical torque play crucial roles in shaping planetary systems. △ Less

Submitted 11 December, 2020; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: 13 pages, 8 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

The CARMENES search for exoplanets around M dwarfs -- LP 714-47b (TOI 442.01): Populating the Neptune desert

Authors: S. Dreizler, I., J., M. Crossfield, D. Kossakowski, P. Plavchan, S., V. Jeffers, J. Kemmer, R. Luque, N. Espinoza, E. Pallé, K. Stassun, E. Matthews, B. Cale, J., A. Caballero, M. Schlecker, J. Lillo-Box, M. Zechmeister, S. Lalitha, A. Reiners, A. Soubkiou, B. Bitsch, M. , et al. (130 additional authors not shown)

Abstract: We report the discovery of a Neptune-like planet (LP 714-47 b, P = 4.05204 d, m_b = 30.8 +/- 1.5 M_earth , R_b = 4.7 +/- 0.3 R_earth ) located in the 'hot Neptune desert'. Confirmation of the TESS Object of Interest (TOI 442.01) was achieved with radial-velocity follow-up using CARMENES, ESPRESSO, HIRES, iSHELL, and PFS, as well as from photometric data using TESS, Spitzer, and ground-based photom… ▽ More We report the discovery of a Neptune-like planet (LP 714-47 b, P = 4.05204 d, m_b = 30.8 +/- 1.5 M_earth , R_b = 4.7 +/- 0.3 R_earth ) located in the 'hot Neptune desert'. Confirmation of the TESS Object of Interest (TOI 442.01) was achieved with radial-velocity follow-up using CARMENES, ESPRESSO, HIRES, iSHELL, and PFS, as well as from photometric data using TESS, Spitzer, and ground-based photometry from MuSCAT2, TRAPPIST- South, MONET-South, the George Mason University telescope, the Las Cumbres Observatory Global Telescope network, the El Sauce telescope, the TUBITAK National Observatory, the University of Louisville Manner Telescope, and WASP-South. We also present high-spatial resolution adaptive optics imaging with the Gemini Near-Infrared Imager. The low uncertainties in the mass and radius determination place LP 714-47 b among physically well-characterised planets, allowing for a meaningful comparison with planet structure models. The host star LP 714-47 is a slowly rotating early M dwarf (T_eff = 3950 +/- 51 K) with a mass of 0.59 +/- 0.02 M_sun and a radius of 0.58 +/- 0.02 R_sun. From long-term photometric monitoring and spectroscopic activity indicators, we determine a stellar rotation period of about 33 d. The stellar activity is also manifested as correlated noise in the radial-velocity data. In the power spectrum of the radial-velocity data, we detect a second signal with a period of 16 days in addition to the four-day signal of the planet. This could be shown to be a harmonic of the stellar rotation period or the signal of a second planet. It may be possible to tell the difference once more TESS data and radial-velocity data are obtained. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: Accepted for publication in A&A. 24 pages, 21 figures, 5 tables

Journal ref: A&A 644, A127 (2020)

Influence of planetary gas accretion on the shape and depth of gaps in protoplanetary discs

Authors: C. Bergez-Casalou, B. Bitsch, A. Pierens, A. Crida, S. N. Raymond

Abstract: It is widely known that giant planets have the capacity to open deep gaps in their natal gaseous protoplanetary discs. It is unclear, however, how gas accretion onto growing planets influences the shape and depth of their growing gaps. We performed isothermal hydrodynamical simulations with the Fargo-2D1D code, which assumes planets accreting gas within full discs that range from 0.1 to 260 AUえーゆー. Th… ▽ More It is widely known that giant planets have the capacity to open deep gaps in their natal gaseous protoplanetary discs. It is unclear, however, how gas accretion onto growing planets influences the shape and depth of their growing gaps. We performed isothermal hydrodynamical simulations with the Fargo-2D1D code, which assumes planets accreting gas within full discs that range from 0.1 to 260 AUえーゆー. The gas accretion routine uses a sink cell approach, in which different accretion rates are used to cope with the broad range of gas accretion rates cited in the literature. We find that the planetary gas accretion rate increases for larger disc aspect ratios and greater viscosities. Our main results show that gas accretion has an important impact on the gap-opening mass: we find that when the disc responds slowly to a change in planetary mass (i.e., at low viscosity), the gap-opening mass scales with the planetary accretion rate, with a higher gas accretion rate resulting in a larger gap-opening mass. On the other hand, if the disc response time is short (i.e., at high viscosity), then gas accretion helps the planet carve a deep gap. As a consequence, higher planetary gas accretion rates result in smaller gap-opening masses. Our results have important implications for the derivation of planet masses from disc observations: depending on the planetary gas accretion rate, the derived masses from ALMA observations might be off by up to a factor of two. We discuss the consequences of the change in the gap-opening mass on the evolution of planetary systems based on the example of the grand tack scenario. Planetary gas accretion also impacts stellar gas accretion, where the influence is minimal due to the presence of a gas-accreting planet. △ Less

Submitted 29 November, 2020; v1 submitted 1 October, 2020; originally announced October 2020.

Journal ref: A&A 643, A133 (2020)