-
Outflow Driven by a Protoplanet Embedded in the TW Hya Disk
Authors:
Tomohiro C. Yoshida,
Hideko Nomura,
Charles J. Law,
Richard Teague,
Yuhito Shibaike,
Kenji Furuya,
Takashi Tsukagoshi
Abstract:
Gas giant planets are formed by gas accretion onto planetary cores in protoplanetary disks. However, direct evidence of this process is still lacking, limiting our understanding of planetary formation processes. During mass accretion, planet-driven outflows may be launched, which could be observable by shock tracers such as sulfur monoxide (SO). We report the detection of SO gas in the protoplanet…
▽ More
Gas giant planets are formed by gas accretion onto planetary cores in protoplanetary disks. However, direct evidence of this process is still lacking, limiting our understanding of planetary formation processes. During mass accretion, planet-driven outflows may be launched, which could be observable by shock tracers such as sulfur monoxide (SO). We report the detection of SO gas in the protoplanetary disk around TW Hya in archival Atacama Large Millimeter/sub-millimeter Array (ALMA) observations. The $\rm SO\ J=8_7 - 7_6\ $ emission line is detected at a $6σ$ significance and localized to the southeast region of the disk with an arc-like morphology. The line center is red-shifted with respect to the systemic velocity by $\sim5\ \rm km\ s^{-1}$. The starting point of the SO emission is located at a planet-carved dust gap at $42$ au. We attribute this to an outflow driven by an embedded protoplanet. Indeed, the observed morphology is well reproduced by a ballistic outflow model. The outflow velocity suggests that the outflow launching source has a mass of $\sim 4 M_\oplus\ (0.012 M_{\rm Jup})$ and the mass-loss rate is $3\times10^{-8} - 1\times10^{-6}\ M_{\rm Jup}\ {\rm yr^{-1}}$. With the relation of mass-loss and mass-accretion rates established for protostars, we estimated the mass-accretion rate onto the protoplanet to be $3\times10^{-7} - 1\times10^{-5}\ M_{\rm Jup}\ {\rm yr^{-1}}$, which matches theoretical predictions for a $\sim 4 M_\oplus$ planet at this separation. The detection of planet-driven outflow provides us a unique opportunity to directly probe the earliest phase of gas giant planet formation.
△ Less
Submitted 19 July, 2024;
originally announced July 2024.
-
Constraints on PDS 70 b and c from the dust continuum emission of the circumplanetary discs considering in situ dust evolution
Authors:
Yuhito Shibaike,
Christoph Mordasini
Abstract:
The young T Tauri star PDS 70 has two gas accreting planets sharing one large gap in a pre-transitional disc. Dust continuum emission from PDS 70 c has been detected by Atacama Large Millimeter/submillimeter Array (ALMA) Band 7, considered as the evidence of a circumplanetary disc, but the emission from PDS 70 b has not. We constrain the planet mass and the gas accretion rate of the planets by int…
▽ More
The young T Tauri star PDS 70 has two gas accreting planets sharing one large gap in a pre-transitional disc. Dust continuum emission from PDS 70 c has been detected by Atacama Large Millimeter/submillimeter Array (ALMA) Band 7, considered as the evidence of a circumplanetary disc, but the emission from PDS 70 b has not. We constrain the planet mass and the gas accretion rate of the planets by introducing a model of dust evolution in the CPDs and reproducing the detection and non-detection. We first develop a 1D steady gas disc model of the CPDs reflecting the planet properties. We then calculate the radial distribution of the dust profiles considering the dust evolution in the gas disc and calculate the total flux density of dust thermal emission from the CPDs. We find positive correlations between the flux density and three planet properties, the planet mass, gas accretion rate, and their product called "MMdot". We then find that the MMdot of PDS 70 c is larger than 0.4 MJ^2/Myr, and the planet mass and the gas accretion rate are larger than 5MJ and 0.02 MJ/Myr, respectively. This is the first case to succeed in obtaining constraints on planet properties from the flux density of dust continuum from a CPD. We also find some loose constraints on the properties of PDS 70 b from the non-detection. We propose possible scenarios for the planets explaining the non-detection respectively detection of the dust emission. The first explanation is that planet c has larger planet mass and/or larger gas accretion rate than planet b. The other possibility is that the CPD of planet c has a larger amount of dust supply and/or weaker turbulence than that of planet b. If the dust supply to planet c is larger than b due to its closeness to the outer dust ring, it is also quantitatively consistent with that planet c has weaker H-alpha line emission than planet b considering the dust extinction effect.
△ Less
Submitted 4 May, 2024; v1 submitted 22 April, 2024;
originally announced April 2024.
-
Delivery of Dust Particles from Protoplanetary Disks onto Circumplanetary Disks of Giant Planets
Authors:
Natsuho Maeda,
Keiji Ohtsuki,
Ryo Suetsugu,
Yuhito Shibaike,
Takayuki Tanigawa,
Masahiro N. Machida
Abstract:
Principal regular satellites of gas giants are thought to be formed by the accumulation of solid materials in circumplanetary disks (CPDs). While there has been significant progress in the study of satellite formation in CPDs, details of the supply of satellite building blocks to CPDs remain unclear. We performed orbital integration of solid particles in the protoplanetary disk (PPD) approaching a…
▽ More
Principal regular satellites of gas giants are thought to be formed by the accumulation of solid materials in circumplanetary disks (CPDs). While there has been significant progress in the study of satellite formation in CPDs, details of the supply of satellite building blocks to CPDs remain unclear. We performed orbital integration of solid particles in the protoplanetary disk (PPD) approaching a planet, considering the gas drag force using the results of three-dimensional hydrodynamical simulations of a local region around the planet. We investigated planetary-mass dependence of the capture positions and capture rates of dust particles accreting onto the CPD. We also examined the degree of dust retention in accreting gas onto the CPD, which is important for determining the ratio of dust-to-gas inflow rates, a key parameter in satellite formation. We found that the degree of dust retention increases with increasing planetary mass for a given dust scale height in the PPD. In the case of a small planet ($M_{\rm p}=0.2M_{\rm Jup}$), most particles with insufficient initial altitudes in the PPD are isolated from the gas in the accreting region. On the other hand, in the case of a massive planet ($M_{\rm p}=1M_{\rm Jup}$), dust particles can be coupled to the vertically accreting gas, even when the dust scale height is about $10-30$\% of the gas scale height. The results of this study can be used for models of dust delivery and satellite formation in the CPDs of gas giants of various masses, including exoplanets.
△ Less
Submitted 17 April, 2024;
originally announced April 2024.
-
Planetesimal formation at the gas pressure bump following a migrating planet II. Effects of dust growth
Authors:
Yuhito Shibaike,
Yann Alibert
Abstract:
Planetesimal formation is still mysterious. One of the ways to form planetesimals is to invoke a gas pressure bump in a protoplanetary disc. In our previous paper, we propose a new scenario in which the piled-up dust at a gas pressure bump created by a migrating planet form planetesimals by streaming instability in a wide region of the disc as the planet migrates inward. In this work, we consider…
▽ More
Planetesimal formation is still mysterious. One of the ways to form planetesimals is to invoke a gas pressure bump in a protoplanetary disc. In our previous paper, we propose a new scenario in which the piled-up dust at a gas pressure bump created by a migrating planet form planetesimals by streaming instability in a wide region of the disc as the planet migrates inward. In this work, we consider the global time evolution of dust and investigate the detailed conditions and results of the planetesimal formation in our scenario. We use a 1D grid single-sized dust evolution model, which can follow the growth of the particles by their mutual collision and their radial drift and diffusion. We calculate the time-evolution of the radial distribution of the peak mass and surface density of the dust in a gas disc perturbed by an embedded migrating planet and investigate if the dust satisfies the condition for planetesimal formation. We find that planetesimals form in a belt-like region between the snowline and the position where the planet reaches its pebble-isolation mass when the strength of turbulence is $10^{-4}\leqα\leq10^{-3}$, which is broadly consistent with observed value. The mechanism of the formation, streaming instability or mutual collision, depends on the timescale of the streaming instability. The total mass of planetesimals also depends on $α$ and is about $30-100~M_{\rm E}$ if the planetary core has already existed at the beginning and grows by gas accretion, but it decreases as the timing of the formation of the planetary core is later. We also provide simple approximate expressions of the surface density and total mass of the planetesimals and find that the total mass strongly depends on the dust mass. We show that planetesimals form in a belt-like region by the combination of the dust pile-up at the gas pressure bump formed by a planet and its inward migration.
△ Less
Submitted 20 June, 2023;
originally announced June 2023.
-
Effective dust growth in laminar circumplanetary discs with magnetic wind-driven accretion
Authors:
Yuhito Shibaike,
Shoji Mori
Abstract:
It has been considered that large satellites around gas planets form in-situ circumplanetary discs (CPDs). However, dust particles supplied into CPDs drift toward the central planets before they grow into satellitesimals, building blocks of the satellites. We investigate the dust growth in laminar CPDs with magnetic wind-driven accretion. In such laminar discs, dust particles can settle onto the m…
▽ More
It has been considered that large satellites around gas planets form in-situ circumplanetary discs (CPDs). However, dust particles supplied into CPDs drift toward the central planets before they grow into satellitesimals, building blocks of the satellites. We investigate the dust growth in laminar CPDs with magnetic wind-driven accretion. In such laminar discs, dust particles can settle onto the mid-plane and grow large by mutual collision more efficient than in classical turbulent CPDs. First, we carry out 3D local MHD simulations of a CPD including all the nonideal MHD effects (Ohmic resistivity, Hall effect and ambipolar diffusion). We investigate if the disc accretion can be governed by magnetic wind-driven accretion and how laminar the disc can be, in a situation where the magnetic disc wind can be launched from the disc. Second, we model 1D steady CPDs consistent with the results of the MHD simulations and calculate the steady radial distributions of the dust profiles in the modeled discs, taking account of the collisional growth, radial drift, fragmentation, and vertical stirring by the Kelvin-Helmholtz instability. We show that satellitesimals can form in such CPDs if the dust-to-gas mass ratio of the inflow to the discs is larger than 0.02, which is 50 times smaller than the critical value in turbulent CPDs. This condition can be satisfied when enough amount of dust piles up at the gas pressure bump created by the planets. This result shows that satellitesimals would form in laminar CPDs with magnetic wind-driven accretion.
△ Less
Submitted 9 December, 2022; v1 submitted 16 November, 2022;
originally announced November 2022.
-
Planetesimal formation at the gas pressure bump following a migrating planet I. Basic characteristics of the new formation model
Authors:
Yuhito Shibaike,
Yann Alibert
Abstract:
To avoid known difficulties in planetesimal formation such as the drift or fragmentation barriers, many scenarios have been proposed. However, in these scenarios, planetesimals form in general only at some specific locations in protoplanetary discs. On the other hand, it is generally assumed in planet formation models and population synthesis models, that planetesimals are broadly distributed in t…
▽ More
To avoid known difficulties in planetesimal formation such as the drift or fragmentation barriers, many scenarios have been proposed. However, in these scenarios, planetesimals form in general only at some specific locations in protoplanetary discs. On the other hand, it is generally assumed in planet formation models and population synthesis models, that planetesimals are broadly distributed in the protoplanetary disc. Here we propose a new scenario in which planetesimals can form in broad areas of the discs. Planetesimals form at the gas pressure bump formed by a first-generation planet (e.g. formed by pebble accretion) and the formation region spreads inward in the disc as the planet migrates. We use a simple 1D Lagrangian particle model to calculate the radial distribution of pebbles in the gas disc perturbed by a migrating embedded planet. We consider that planetesimals form by streaming instability at the points where the pebble-to-gas density ratio on the mid-plane becomes larger than unity. We also study the effect of some key parameters like the ones of the gas disc model, the pebble mass flux, the migration speed of the planet, and the strength of turbulence. We find that planetesimals form in wide areas of the discs provided the flux of pebbles is typical and the turbulence is not too strong. The planetesimal surface density depends on the pebble mass flux and the migration speed of the planet. The total mass of the planetesimals and the orbital position of the formation area depend strongly on the pebble mass flux. We also find that the profile of the planetesimal surface density and its slope can be estimated by very simple equations. We show that our new scenario can explain the formation of planetesimals in broad areas. The simple estimates we provide for the planetesimal surface density profile can be used as initial conditions for population synthesis models.
△ Less
Submitted 20 October, 2020;
originally announced October 2020.
-
The Galilean Satellites Formed Slowly from Pebbles
Authors:
Yuhito Shibaike,
Chris W. Ormel,
Shigeru Ida,
Satoshi Okuzumi,
Takanori Sasaki
Abstract:
It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter's formation. However, understanding the specifics of the formation process is challenging as both small particles (pebbles) as well as the satellites are subject to fast migration processes. Here, we hypothesize a new scenario for the origin of the Galilean system, based on the ca…
▽ More
It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter's formation. However, understanding the specifics of the formation process is challenging as both small particles (pebbles) as well as the satellites are subject to fast migration processes. Here, we hypothesize a new scenario for the origin of the Galilean system, based on the capture of several planetesimal seeds and subsequent slow accretion of pebbles. To halt migration, we invoke an inner disk truncation radius, and other parameters are tuned for the model to match physical, dynamical, compositional, and structural constraints. In our scenario it is natural that Ganymede's mass is determined by pebble isolation. Our slow-pebble-accretion scenario then reproduces the following characteristics: (1) the mass of all the Galilean satellites; (2) the orbits of Io, Europa, and Ganymede captured in mutual 2:1 mean motion resonances; (3) the ice mass fractions of all the Galilean satellites; (4) the unique ice-rock partially differentiated Callisto and the complete differentiation of the other satellites. Our scenario is unique to simultaneously reproduce these disparate properties.
△ Less
Submitted 28 October, 2019; v1 submitted 31 August, 2019;
originally announced September 2019.
-
Photophoresis in the circumjovian disk and its impact on the orbital configuration of the Galilean satellites
Authors:
Sota Arakawa,
Yuhito Shibaike
Abstract:
Jupiter has four large regular satellites called the Galilean satellites: Io, Europa, Ganymede, and Callisto. The inner three of the Galilean satellites orbit in a 4:2:1 mean motion resonance; therefore their orbital configuration may originate from the stopping of the migration of Io near the bump in the surface density distribution and following resonant trapping of Europa and Ganymede. The form…
▽ More
Jupiter has four large regular satellites called the Galilean satellites: Io, Europa, Ganymede, and Callisto. The inner three of the Galilean satellites orbit in a 4:2:1 mean motion resonance; therefore their orbital configuration may originate from the stopping of the migration of Io near the bump in the surface density distribution and following resonant trapping of Europa and Ganymede. The formation mechanism of the bump near the orbit of the innermost satellite, Io, is not yet understood, however. Here, we show that photophoresis in the circumjovian disk could be the cause of the bump, using analytic calculations of steady-state accretion disks. We propose that photophoresis in the circumjovian disk could stop the inward migration of dust particles near the orbit of Io. The resulting dust depleted inner region would have a higher ionization fraction, and thus admit increased magnetorotational-instability-driven accretion stress than the outer region. The increase of the accretion stress at the photophoretic dust barrier would form a bump in the surface density distribution, halting the migration of Io.
△ Less
Submitted 8 August, 2019;
originally announced August 2019.
-
Possibility of Concentration of Non-volatile Species near the Surface of Comet 67P/Churyumov-Gerasimenko
Authors:
Taiki Suzuki,
Yoshiharu Shinnaka,
Liton Majumdar,
Takashi Shibata,
Yuhito Shibaike,
Hideko Nomura,
Harumi Minamoto
Abstract:
The cometary materials are thought to be the reservoir of primitive materials in the Solar System. The recent detection of glycine and CH$_3$NH$_2$ by the ROSINA mass spectrometer in the coma of 67P/Churyumov-Gerasimenko suggests that amino acids and their precursors may have been formed in an early evolutionary phase of the Solar System. We investigate the evolution of cometary interior consideri…
▽ More
The cometary materials are thought to be the reservoir of primitive materials in the Solar System. The recent detection of glycine and CH$_3$NH$_2$ by the ROSINA mass spectrometer in the coma of 67P/Churyumov-Gerasimenko suggests that amino acids and their precursors may have been formed in an early evolutionary phase of the Solar System. We investigate the evolution of cometary interior considering the evaporation process of water followed by the concentration of non-volatile species. We develop a Simplified Cometary Concentration Model (SCCM) to simulate the evaporation and concentration processes on the cometary surface.We use 67P/Churyumov-Gerasimenko as the benchmark of SCCM. We investigate the depth of the layer where non-volatile species concentrate after the numerous passages of perihelion after the formation of the Solar System. As a result, the SCCM explains the observed production rates of water and CH$_3$NH$_2$ at 100 comet years. SCCM results suggest that the non-volatile species would concentrate at depths between 0 and 100cm of comet surface within 10 comet years. Our results also suggest that the non-volatile species would concentrate several meters beneath the surface before it hit the early Earth. This specific mass of non-volatile species may provide unique chemical condition to the volcanic hot spring pools.
△ Less
Submitted 3 December, 2020; v1 submitted 16 January, 2019;
originally announced January 2019.
-
Impacts of dust feedback on a dust ring induced by a planet in a protoplanetary disk
Authors:
Kazuhiro D. Kanagawa,
Takayuki Muto,
Satoshi Okuzumi,
Takayuki Tanigawa,
Tetsuo Taki,
Yuhito Shibaike
Abstract:
When a planet forms a deep gap in a protoplanetary disk, dust grains cannot pass through the gap. As a consequence, the density of the dust grains can increase up to the same level of the density of the gas at the outer edge. The feedback on the gas from the drifting dust grains is not negligible, in such a dusty region. We carried out two-dimensional two-fluid (gas and dust) hydrodynamic simulati…
▽ More
When a planet forms a deep gap in a protoplanetary disk, dust grains cannot pass through the gap. As a consequence, the density of the dust grains can increase up to the same level of the density of the gas at the outer edge. The feedback on the gas from the drifting dust grains is not negligible, in such a dusty region. We carried out two-dimensional two-fluid (gas and dust) hydrodynamic simulations. We found that when the radial flow of the dust grains across the gap is halted, a broad ring of the dust grains can be formed because of the dust feedback and the diffusion of the dust grains. The minimum mass of the planet to form the broad dust ring is consistent with the pebble-isolation mass, in the parameter range of our simulations. The broad ring of the dust grains is good environment for the formation of the protoplanetary solid core. If the ring is formed in the disk around the sun-like star at $\sim 2\ \mbox{AU}$, a massive solid core ($\sim 50M_{\oplus}$) can be formed within the ring, which may be connected to the formation of Hot Jupiters holding a massive solid core such as HD 149026b. In the disk of the dwarf star, a number of Earth-sized planets can be formed within the dust ring around $\sim 0.5\ \mbox{AU}$, which potentially explain the planet system made of multiple Earth-sized planets around the dwarf star such as TRAPPIST-1.
△ Less
Submitted 19 November, 2018; v1 submitted 12 October, 2018;
originally announced October 2018.
-
Satellitesimal Formation via Collisional Dust Growth in Steady Circumplanetary Disks
Authors:
Yuhito Shibaike,
Satoshi Okuzumi,
Takanori Sasaki,
Shigeru Ida
Abstract:
The icy satellites around Jupiter are considered to have formed in a circumplanetary disk. While previous models focused on the formation of satellites starting from satellitesimals, the question of how satellitesimals form from smaller dust particles has not been addressed so far. In this work, we study the possibility that satellitesimals form in situ in a circumplanetary disk. We calculate the…
▽ More
The icy satellites around Jupiter are considered to have formed in a circumplanetary disk. While previous models focused on the formation of satellites starting from satellitesimals, the question of how satellitesimals form from smaller dust particles has not been addressed so far. In this work, we study the possibility that satellitesimals form in situ in a circumplanetary disk. We calculate the radial distribution of the surface density and representative size of icy dust particles that grow by colliding with each other and drift toward the central planet in a steady circumplanetary disk with a continuous supply of gas and dust from the parent protoplanetary disk. The radial drift barrier is overcome if the ratio of the dust to gas accretion rates onto the circumplanetary disk, $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$, is high and the strength of turbulence, $α$, is not too low. The collision velocity is lower than the critical velocity of fragmentation when $α$ is low. Taken together, we find that the conditions for satellitesimal formation via dust coagulation are given by $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}\ge1$ and $10^{-4}\leα<10^{-2}$. The former condition is generally difficult to achieve, suggesting that the in-situ satellitesimal formation via particle sticking is viable only under an extreme condition. We also show that neither satellitesimal formation via the collisional growth of porous aggregates nor via streaming instability is viable as long as $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$ is low.
△ Less
Submitted 3 August, 2017;
originally announced August 2017.
-
Excavation and Melting of the Hadean Continental Crust by Late Heavy Bombardment
Authors:
Yuhito Shibaike,
Takanori Sasaki,
Shigeru Ida
Abstract:
No Hadean rocks have ever been found on Earth's surface except for zircons---evidence of continental crust, suggesting that Hadean continental crust existed but later disappeared. One hypothesis for the disappearance of the continental crust is excavation/melting by the Late Heavy Bombardment (LHB), a concentration of impacts in the last phase of the Hadean eon. In this paper, we calculate the eff…
▽ More
No Hadean rocks have ever been found on Earth's surface except for zircons---evidence of continental crust, suggesting that Hadean continental crust existed but later disappeared. One hypothesis for the disappearance of the continental crust is excavation/melting by the Late Heavy Bombardment (LHB), a concentration of impacts in the last phase of the Hadean eon. In this paper, we calculate the effects of LHB on Hadean continental crust in order to investigate this hypothesis. Approximating the size-frequency distribution of the impacts by a power-law scaling with an exponent α as a parameter, we have derived semi-analytical expressions for the effects of LHB impacts. We calculated the total excavation/melting volume and area affected by the LHB from two constraints of LHB on the moon, the size of the largest basin during LHB, and the density of craters larger than 20 km. We also investigated the effects of the value of α. Our results show that LHB does not excavate/melt all of Hadean continental crust directly, but over 70% of the Earth's surface area can be covered by subsequent melts in a broad range of α. If there have been no overturns of the continental crust until today, LHB could be responsible for the absence of Hadean rocks because most of Hadean continental crust is not be exposed on the Earth's surface in this case.
△ Less
Submitted 2 November, 2015;
originally announced November 2015.