Source author record

Ryuki Hyodo

Ryuki Hyodo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

A "no-drift" runaway pile-up of pebbles in protoplanetary disks II. Characteristics of the resulting planetesimal belt

Forming planetesimals from pebbles is a major challenge in our current understanding of planet formation. In a protoplanetary disk, pebbles drift inward near the disk midplane via gas drag and they may enter a dead zone. In this context, we identified that the backreaction of the drag of pebbles onto the gas could lead to a runaway pile-up of pebbles, the so-called no-drift mechanism. We improve upon the previous study of the no-drift mechanism by investigating the nature and characteristics of the resultant planetesimal belt. We performed 1D diffusion-advection simulations of drifting pebbles in the outer region of a dead zone by including the backreaction to the radial drift of pebbles and including planetesimal formation via the streaming instability. We considered the parameters that regulate gas accretion and vertical stirring of pebbles in the disk midplane. In this study, the pebble-to-gas mass flux ($F_{\rm p/g}$) was fixed as a parameter. We find that planetesimals initially form within a narrow ring whose width expands as accumulating pebbles radially diffuse over time. The system finally reaches a steady-state where the width of the planetesimal belt no longer changes. A non-negligible total mass of planetesimals (more than one Earth mass) is formed for a disk having $F_{\rm p/g} \gtrsim 0.1$ for more than $\sim 10-100$ kyr with nominal parameters: a gas mass flux of $\gtrsim10^{-8} {\rm M}_\oplus$/yr, $τ_{\rm s} \simeq 0.01-0.1$, $α_{\rm mid} \lesssim 10^{-4}$, and $α_{\rm acc} \simeq 10^{-3}-10^{-2}$ at $r \lesssim 10$ au, where $r$, $τ_{\rm s}$, $α_{\rm mid}$, and $α_{\rm acc}$ are the heliocentric distance, the Stokes number, and the parameters in a dead zone controlling the efficiencies of vertical turbulent diffusion of pebbles (i.e., scale height of pebbles) and gas accretion of the $α$-disk (i.e., gas surface density), respectively.

preprint2022arXiv

Asteroid Flyby Cycler Trajectory Design Using Deep Neural Networks

Asteroid exploration has been attracting more attention in recent years. Nevertheless, we have just visited tens of asteroids while we have discovered more than one million bodies. As our current observation and knowledge should be biased, it is essential to explore multiple asteroids directly to better understand the remains of planetary building materials. One of the mission design solutions is utilizing asteroid flyby cycler trajectories with multiple Earth gravity assists. An asteroid flyby cycler trajectory design problem is a subclass of global trajectory optimization problems with multiple flybys, involving a trajectory optimization problem for a given flyby sequence and a combinatorial optimization problem to decide the sequence of the flybys. As the number of flyby bodies grows, the computation time of this optimization problem expands maliciously. This paper presents a new method to design asteroid flyby cycler trajectories utilizing a surrogate model constructed by deep neural networks approximating trajectory optimization results. Since one of the bottlenecks of machine learning approaches is the computation time to generate massive trajectory databases, we propose an efficient database generation strategy by introducing pseudo-asteroids satisfying the Karush-Kuhn-Tucker conditions. The numerical result applied to JAXA's DESTINY+ mission shows that the proposed method is practically applicable to space mission design and can significantly reduce the computational time for searching asteroid flyby sequences.

preprint2022arXiv

Challenges in forming Phobos and Deimos directly from a splitting of an ancestral single moon

The origin and evolution of Martian moons have been intensively debated in recent years. It is proposed that Phobos and Deimos may originate directly from a splitting of an ancestral moon orbiting at around the Martian synchronous orbit. At this hypothetical splitting, the apocenter of the inner moon (presumed as Phobos) and the pericenter of the outer moon (presumed as Deimos) are reported to coincide, in that, their semi-major axes reside inside and outside the Martian synchronous orbit with non-zero eccentricities, respectively. However, the successive orbital evolution of the two moons is not studied. Here, we perform direct $N$-body orbital integrations of the moons, including the Martian oblateness of the $J_2$ and $J_4$ terms. We show that the two moons, while they precess, likely collide within $\sim 10^4$ years with an impact velocity of $v_{\rm imp} \sim 100-300$ m s$^{-1}$ ($\sim 10-30$ times moons' escape velocity) and with an isotropic impact direction. The impact occurs around the apocenter and the pericenter of the inner and outer moons, respectively, where the timescale of this periodic orbital alignment is regulated by the precession. By performing additional impact simulations, we show that such a high-velocity impact likely results in a disruptive outcome, forming a debris ring at around the Martian synchronous orbit, from which several small moons would accrete. Such an evolutionary path would eventually form a different Martian moons system from the one we see today. Therefore, it seems unlikely that Phobos and Deimos are split directly from a single ancestral moon.

preprint2021arXiv

A "no-drift" runaway pile-up of pebbles in protoplanetary disks in which midplane turbulence increases with radius

A notable challenge of planet formation is to find a path to directly form planetesimals from small particles. We aim to understand how drifting pebbles pile up in a protoplanetary disk with a non-uniform turbulence structure. We consider a disk structure in which the midplane turbulence viscosity is increasing with radius in protoplanetary disks as in the outer region of a dead zone. We perform 1D diffusion-advection simulations of pebbles that include back-reaction (the inertia) to radial drift and vertical/radial diffusion of pebbles for a given pebble-to-gas mass flux. We report a new mechanism, the "no-drift" runaway pile-up, leading to a runaway accumulation of pebbles in disks, thus favoring the formation of planetesimals by streaming and/or gravitational instabilities. This occurs when pebbles drifting in from the outer disk and entering a dead zone experience a decrease in vertical turbulence. The scale height of the pebble subdisk then decreases, and for small enough values of the turbulence in the dead zone and high values of the pebble to gas flux ratio, the back-reaction of pebbles on gas leads to a significant decrease in their drift velocity and thus their progressive accumulation. This occurs when the ratio of the flux of pebbles to that of the gas is large enough so that the effect dominates over any Kelvin-Helmholtz shear instability. This process is independent of the existence of a pressure bump.

preprint2021arXiv

Planetesimal formation around the snow line. II. Dust or pebbles?

Around the snow line, icy pebbles and silicate dust may locally pile-up and form icy and rocky planetesimals via streaming instability and/or gravitational instability. We perform 1D diffusion-advection simulations that include the back-reaction to radial drift and diffusion of icy pebbles and silicate dust, ice sublimation, release of silicate dust, and their recycling through recondensation and sticking onto pebbles outside the snow line. We use a realistic description of the scale height of silicate dust obtained from Ida et al. and that of pebbles including the effects of a Kelvin-Helmholtz instability. We study the dependence of solid pile-up on distinct effective viscous parameters for turbulent diffusions in the radial and vertical directions ($α_{\rm Dr}$ and $α_{\rm Dz}$) and for the gas accretion to the star ($α_{\rm acc}$) as well as that on the pebble-to-gas mass flux ($F_{\rm p/g}$). We derive the sublimation width of drifting icy pebbles which is a critical parameter to characterize the pile-up of silicate dust and pebbles around the snow line. We identify a parameter space (in the $F_{\rm p/g}-α_{\rm acc}-α_{\rm Dz}(=α_{\rm Dr})$ space) where pebbles no longer drift inward to reach the snow line due to the back-reaction that slows down radial velocity of pebbles. We show that the pile-up of solids around the snow line occurs in a broader range of parameters for $α_{\rm acc}=10^{-3}$ than for $α_{\rm acc}=10^{-2}$. Above a critical $F_{\rm p/g}$ value, the runaway pile-up of silicate dust inside the snow line is favored for $α_{\rm Dr}/α_{\rm acc} \ll 1$, while that of pebbles outside the snow line is favored for $α_{\rm Dr}/α_{\rm acc} \sim 1$. Our results imply that a distinct evolutionary path could produce a diversity of outcomes in terms of planetesimal formation around the snow line.

preprint2020arXiv

Escape and accretion by cratering impacts: Formulation of scaling relations for high-speed ejecta

Numerous small bodies inevitably lead to cratering impacts on large planetary bodies during planet formation and evolution. As a consequence of these small impacts, a fraction of the target material escapes from the gravity of the large body, and a fraction of the impactor material accretes onto the target surface, depending on the impact velocities and angles. Here, we study the mass of the high-speed ejecta that escapes from the target gravity by cratering impacts when material strength is neglected. We perform a large number of cratering impact simulations onto a planar rocky target using the smoothed particle hydrodynamics method. We show that the escape mass of the target material obtained from our numerical simulations agrees with the prediction of a scaling law under a point-source assumption when $v_{\rm imp} \gtrsim 12 v_{\rm esc}$, where $v_{\rm imp}$ is the impact velocity and $v_{\rm esc}$ is the escape velocity of the target. However, we find that the point-source scaling law overestimates the escape mass up to a factor of $\sim 70$, depending on the impact angle, when $v_{\rm imp} \lesssim 12 v_{\rm esc}$. Using data obtained from numerical simulations, we derive a new scaling law for the escape mass of the target material for $v_{\rm imp} \lesssim 12 v_{\rm esc}$. We also derive a scaling law that predicts the accretion mass of the impactor material onto the target surface upon cratering impacts by numerically evaluating the escape mass of the impactor material. Our newly derived scaling laws are useful for predicting the escape mass of the target material and the accretion mass of the impactor material for a variety of cratering impacts that would occur on large planetary bodies during planet formation.

preprint2020arXiv

Planetesimal formation around the snow line: I. Monte Carlo simulations of silicate dust pile-up in a turbulent disk

Context: The formation of rocky planetesimals is a long-standing problem in planet formation theory. One of the possibilities is that it results from gravitational instability as a result of pile-up of small silicate dust particles released from sublimating icy pebbles that pass the snow line. Aims: We want to understand and quantify the role of the water snow line for the formation of rock-rich and ice-rich planetesimals. In this paper, we focus on the formation of rock-rich planetesimals. A companion paper examines the combined formation of both rock-rich and ice-rich planetesimals. Methods: We develop a new Monte Carlo code to calculate the radial evolution of silicate particles in a turbulent accretion disk, accounting for the back-reaction (i.e., inertia) of the particles on their radial drift velocity and diffusion. Results depend in particular on the particle injection width (determined from the radial sublimation width of icy pebbles), the pebble scale height and the pebble mass flux through the disk. The scale height evolution of the silicate particles, which is the most important factor for the runaway pile-up, is automatically calculated in this Lagrange method. Results: From the numerical results, we derive semi-analytical relations for the scale height of the silicate dust particles and the particles-to-gas density ratio at the midplane, as functions of a pebble-to-gas mass flux ratio and the $α$ parameters for disk gas accretion and vertical/radial diffusion. We find that the runaway pile-up of the silicate particles (formation of rocky planetesimals) occurs if the pebble-to-gas mass flux ratio is $> [(α_{Dz}/α_{acc})/3 \times 10^{-2}]^{1/2}$ where $α_{Dz}$ and $α_{acc}$ are the $α$ parameters for vertical turbulent diffusion and disk gas accretion.

preprint2016arXiv

Formation of Centaurs' rings through their partial tidal disruption during planetary encounters

Centaurs are minor planets orbiting between Jupiter and Neptune that have or had crossing orbits with one or more giant planets. Recent observations and reinterpretation of previous observations have revealed the existence of ring systems around 10199 Chariklo and 2060 Chiron. However, the origin of the ring systems around such a minor planet is still an open question. Here, we propose that the tidal disruption of a differentiated object that experiences a close encounter with a giant planet could naturally form diverse ring-satellite systems around the Centaurs. During the close encounter, the icy mantle of the passing object is preferentially ripped off by the planet's tidal force and the debris is distributed mostly within the Roche limit of the largest remnant body. Assuming the existence of $20-50$wt% silicate core below the icy mantle, a disk of particles is formed when the objects pass within $0.4-0.8$ of the planet's Roche limit with the relative velocity at infinity $3-6$km s$^{-1}$ and 8h initial spin period of the body. The resultant ring mass is $0.1-10$% of the central object's mass. Such particle disks are expected to spread radially, and materials spreading beyond the Roche limit would accrete into satellite(s). Our numerical results suggest that ring formation would be a natural outcome of such extreme close encounters and Centaurs can naturally have such ring systems because they cross the orbits of the giant planets.

preprint2016arXiv

Ring Formation around Giant Planets by Tidal Disruption of a Single Passing Large Kuiper Belt Object

The origin of rings around giant planets remains elusive. Saturn's rings are massive and made of 90-95% of water ice. In contrast, the much less massive rings of Uranus and Neptune are dark and likely to have higher rock fraction. Here we investigate, for the first time, the tidal disruption of a passing object, including the subsequent formation of planetary rings. First, we perform SPH simulations of the tidal destruction of big differentiated objects ($M_{\rm body}=10^{21-23}$) that experience close encounters with Saturn or Uranus. We find that about $0.1-10$% of the mass of the passing body is gravitationally captured around the planet. However, these fragments are initially big chunks and have highly eccentric orbits around the planet. Then, we perform N-body simulations including the planet's oblateness, starting with data obtained from the SPH simulations. Our N-body simulations show that the chunks are tidally destroyed during their next several orbits. Their individual orbits then start to precess incoherently around the planet's equator, which enhances their encounter velocities on longer-term evolution, resulting in more destructive impacts. These collisions would damp their eccentricities resulting in a progressive collapse of the debris cloud into a thin equatorial and low-eccentricity ring. These high energy impacts are expected to be catastrophic enough to produce small particles. Our numerical results also show that the mass of formed rings is large enough to explain current rings including inner regular satellites around Saturn and Uranus. In the case of Uranus, a body can go deeper inside the planet's Roche limit resulting in a more efficient capture of rocky material compared to Saturn's case in which mostly ice is captured. Thus, our results can naturally explain the compositional difference between the rings of Saturn, Uranus and Neptune.

preprint2014arXiv

Formation of Multiple-Satellite Systems From Low-Mass Circumplanetary Particle Disks

Circumplanetary particle disks would be created in the late stage of planetary formation either by impacts of planetary bodies or disruption of satellites or passing bodies, and satellites can be formed by accretion of disk particles spreading across the Roche limit. Previous N-body simulation of lunar accretion focused on the formation of single-satellite systems from disks with large disk-to-planet mass ratios, while recent models of the formation of multiple-satellite systems from disks with smaller mass ratios do not take account of gravitational interaction between formed satellites. In the present work, we investigate satellite accretion from particle disks with various masses, using N-body simulation. In the case of accretion from somewhat less massive disks than the case of lunar accretion, formed satellites are not massive enough to clear out the disk, but can become massive enough to gravitationally shepherd the disk outer edge and start outward migration due to gravitational interaction with the disk. When the radial location of the 2:1 mean motion resonance of the satellite reaches outside the Roche limit, the second satellite can be formed near the disk outer edge, and then the two satellites continue outward migration while being locked in the resonance. Co-orbital satellites are found to be occasionally formed on the orbit of the first satellite. Our simulations also show that stochastic nature involved in gravitational interaction and collision between aggregates in the tidal environment can lead to diversity in the final mass and orbital architecture, which would be expected in satellite systems of exoplanets.