Researcher profile

Shigeru Ida

Shigeru Ida contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 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

ALMA High-resolution Multiband Analysis for the Protoplanetary Disk around TW Hya

We present a high-resolution (2.5 au) multiband analysis of the protoplanetary disk around TW Hya using ALMA long baseline data at Bands 3, 4, 6, and 7. We aim to reconstruct a high-sensitivity millimeter continuum image and revisit the spectral index distribution. The imaging is performed by combining new ALMA data at Bands 4 and 6 with available archive data. Two methods are employed to reconstruct the images; multi-frequency synthesis (MFS) and the fiducial image-oriented method, where each band is imaged separately and the frequency dependence is fitted pixel by pixel. We find that the MFS imaging with the second order of Taylor expansion can reproduce the frequency dependence of the continuum emission between Bands 3 and 7 in a manner consistent with previous studies and is a reasonable method to reconstruct the spectral index map. The image-oriented method provides a spectral index map consistent with the MFS imaging, but with a two times lower resolution. Mock observations of an intensity model were conducted to validate the images from the two methods. We find that the MFS imaging provides a high-resolution spectral index distribution with an uncertainty of $<10$~\%. Using the submillimeter spectrum reproduced from our MFS images, we directly calculated the optical depth, power-law index of the dust opacity coefficient ($β$), and dust temperature. The derived parameters are consistent with previous works, and the enhancement of $β$ within the intensity gaps is also confirmed, supporting a deficit of millimeter-sized grains within the gaps.

preprint2022arXiv

Monte Carlo Simulation of Dust Particles in a Protoplanetary Disk: Crystalline to Amorphous Silicate Ratio in Comets

Observationally inferred crystalline abundance in silicates in comets, which should have been formed in the outer region of a protoplanetary disk, is relatively high (~ 10-60%), although crystalline silicates would be formed by annealing of amorphous precursors in the disk inner region. In order to quantitatively address this puzzle, we have performed Monte Carlo simulation of advection/diffusion of silicate particles in a turbulent disk, in the setting based on pebble accretion model: pebbles consisting of many small amorphous silicates embedded in icy mantle are formed in the disk outer region, silicate particles are released at the snow line, crystalline silicate particles are produced at the annealing line, the silicate particles diffused beyond the snow line, and they eventually stick to drifting pebbles to come back to the snow line. In a simple case without the sticking and with a steady pebble flux, we show through the simulations and analytical arguments that crystalline components in silicate materials beyond the snow line is robustly and uniformly ~ 5%. On the other hand, in a more realistic case with the sticking and with a decaying pebble flux, the crystalline abundance is raised up to ~ 20-25%, depending on the ratio of decay and diffusion timescales. This abundance is consistent with the observations. In this investigation, we assume a simple steady accretion disk. The simulations coupled with the disk evolution is needed for more detailed comparison with observed data.

preprint2022arXiv

Monte Carlo Simulation of Sugar Synthesis on Icy Dust Particles Intermittently Irradiated by UV in a Protoplanetary Disk

Context. While synthesis of organic molecules in molecular clouds or protoplanetary disks is complex, observations of interstellar grains, analyses of carbonaceous chondrites, and UV photochemistry experiments are rapidly developing and providing constraints on and clues to the complex organic molecule synthesis in space. It motivates us to construct a theoretical synthesis model. Aims. We develop a new code to simulate global reaction sequences of organic molecules to apply it for sugar synthesis by intermittent UV irradiation on the surface of icy particles in a protoplanetary disk. Here we show the first results of our new simulation. Methods. We apply a Monte Carlo method to select reaction sequences from all possible reactions, using the graph-theoretic matrix model for chemical reactions and modeling reactions on the icy particles during UV irradiation. Results. We here obtain the results consistent with the organic molecules in carbonaceous chondrites and obtained by the experiments, however, through a different pathway from the conventional formose reactions previously suggested. During UV irradiation, loosely-bonded O-rich large molecules are continuously created and destroyed. After UV irradiation is turned off, the ribose abundance rapidly increases, through the decomposition of the large molecules with break-ups of O-O bonds and replacements of C-OH by C-H to reach O/C = 1 for sugars. The sugar abundance is regulated mostly by the total atomic ratio H/O of starting materials, but not by their specific molecule forms. Deoxyribose is simultaneously synthesized, and most of the molecules end up with complex C-rich molecules.

preprint2021arXiv

A &#34;no-drift&#34; 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 &#34;no-drift&#34; 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

A new and simple prescription for planet orbital migration and eccentricity damping by planet-disc interactions based on dynamical friction

During planet formation gravitational interaction between a planetary embryo and the protoplanetary gas disc causes orbital migration of the planetary embryo, which plays an important role in shaping the final planetary system. While migration sometimes occurs in the supersonic regime, wherein the relative velocity between the planetary embryo and the gas is higher than the sound speed, migration prescriptions proposed thus far describing the planet-disc interaction force and the timescales of orbital change in the supersonic regime are inconsistent with one another. Here we discuss the details of existing prescriptions in the literature and derive a new simple and intuitive formulation for planet-disc interactions based on dynamical friction that can be applied in both supersonic and subsonic cases. While the existing prescriptions assume particular disc models, ours include the explicit dependence on the disc parameters; hence it can be applied to discs with any radial surface density and temperature dependence (except for the local variations with radial scales less than the disc scale height). Our prescription will reduce the uncertainty originating from different literature formulations of planet migration and will be an important tool to study planet accretion processes, especially when studying the formation of close-in low-mass planets that are commonly found in exoplanetary systems.

preprint2020arXiv

Orbital evolution of Saturn&#39;s satellites due to the interaction between the moons and massive Saturn&#39;s rings

Saturn&#39;s mid-sized moons (satellites) have a puzzling orbital configuration with trapping in mean-motion resonances with every other pairs (Mimas-Tethys 4:2 and Enceladus-Dione 2:1). To reproduce their current orbital configuration on the basis of Crida & Charnoz&#39;s model of satellite formation from a hypothetical ancient massive rings, adjacent pairs must pass 1st-order mean-motion resonances without being trapped. The trapping could be avoided by fast orbital migration and/or excitation of the satellite&#39;s eccentricity caused by gravitational interactions between the satellites and the rings (the disk), which are still unknown. In our research, we investigate the satellite orbital evolution due to interactions with the disk through full N-body simulations. We performed global high-resolution N-body simulations of a self-gravitating particle disk interacting with a single satellite. We used $N \sim 10^5$ particles for the disk. Gravitational forces of all the particles and their inelastic collisions are taken into account. As a result, dense short-wavelength wake structure is created by the disk self-gravity and global spiral arms with $m \sim$ a few is induced by the satellite. The self-gravity wakes regulate the orbital evolution of the satellite, which has been considered as a disk spreading mechanism but not as a driver for the orbital evolution. The self-gravity wake torque to the satellite is so effective that the satellite migration is much faster than that was predicted with the spiral arms torque. It provides a possible model to avoid the resonance capture of adjacent satellite pairs and establish the current orbital configuration of Saturn&#39;s mid-sized satellites.

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.

preprint2020arXiv

Uranian Satellite Formation by Evolution of a Water Vapor Disk Generated by a Giant Impact

The ice-giant planet Uranus likely underwent a giant impact, given that its spin axis is tilted by 98 degrees. That its satellite system is equally inclined and prograde suggests that it was formed as a consequence of the impact. However, the disks predicted by the impact simulations generally have sizes one order smaller and masses two orders larger than those of the observed system at present. Here we show, by means of a theoretical model, that the Uranian satellite formation is regulated by the evolution of the impact-generated disk. Because the vaporization temperature of water ice is low and both Uranus and the impactor are assumed to be ice-dominated, we can conclude that the impact-generated disk has mostly vaporized. We predict that the disk lost a significant amount of water vapour mass and spread to the levels of the current system until the disk cooled down enough for ice condensation and accretion of icy particles to begin. From the predicted distribution of condensed ices, our N-body simulation is able to reproduce the observed mass-orbit configuration of Uranian satellites. This scenario contrasts with the giant-impact model for the Earth&#39;s Moon, in which about half of the compact, impact-generated, solid or liquid disk is immediately incorporated into the Moon on impact.

preprint2019arXiv

The Galilean Satellites Formed Slowly from Pebbles

It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter&#39;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&#39;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.

preprint2012arXiv

Dynamics and Accretion of Planetesimals

We review the basic dynamics and accretion of planetesimals by showing N-body simulations. The orbits of planetesimals evolve through two-body gravitational relaxation: viscous stirring increases the random velocity and dynamical friction realizes the equiparation of the random energy. In the early stage of planetesimal accretion the growth mode of planetesimals is runaway growth where larger planetesimals grow faster than smaller ones. When a protoplanet (runaway-growing planetesimal) exceeds a critical mass the growth mode shifts to oligarchic growth where similar-sized protoplanets grow keeping a certain orbital separation. The final stage of terrestrial planet formation is collision among protoplanets known as giant impacts. We also summarize the dynamical effects of disk gas on planets and the core accretion model for formation of gas giants and discuss the diversity of planetary systems.

preprint2010arXiv

Planet Migration

Planet migration is the process by which a planet&#39;s orbital radius changes in time. The main agent for causing gas giant planet migration is the gravitational interaction of the young planet with the gaseous disk from which it forms. We describe the migration rates resulting from these interactions based on a simple model for disk properties. These migration rates are higher than is reasonable for planet survival. We discuss some proposed models for which the migration rates are lower. There are major uncertainties in migration rates due to a lack of knowledge about the detailed physical properties of disks. We also describe some additional forms of migration.

preprint2010arXiv

The Occurrence and Mass Distribution of Close-in Super-Earths, Neptunes, and Jupiters

The questions of how planets form and how common Earth-like planets are can be addressed by measuring the distribution of exoplanet masses and orbital periods. We report the occurrence rate of close-in planets (with orbital periods less than 50 days) based on precise Doppler measurements of 166 Sun-like stars. We measured increasing planet occurrence with decreasing planet mass (M). Extrapolation of a power law mass distribution fitted to our measurements, df/dlogM = 0.39M^-0.48, predicts that 23% of stars harbor a close-in Earth-mass planet (ranging from 0.5 to 2.0 Earth masses). Theoretical models of planet formation predict a deficit of planets in the domain from 5 to 30 Earth masses and with orbital periods less than 50 days. This region of parameter space is in fact well populated, implying that such models need substantial revision.

preprint2009arXiv

N-body simulations of planetary accretion around M dwarf stars

We have investigated planetary accretion from planetesimals in terrestrial planet regions inside the ice line around M dwarf stars through N-body simulations including tidal interactions with disk gas. Because of low luminosity of M dwarfs, habitable zones (HZs) are located in inner regions. In the close-in HZ, type-I migration and the orbital decay induced by eccentricity damping are efficient according to the high disk gas density in the small orbital radii. In the case of full efficiency of type-I migration predicted by the linear theory, we found that protoplanets that migrate to the vicinity of the host star undergo close scatterings and collisions, and 4 to 6 planets eventually remain in mutual mean motion resonances and their orbits have small eccentricities and they are stable both before and after disk gas decays. In the case of slow migration, the resonant capture is so efficient that densely-packed ~ 40 small protoplanets remain in mutual mean motion resonances. In this case, they start orbit crossing, after the disk gas decays and eccentricity damping due to tidal interaction with gas is no more effective. Through merging of the protoplanets, several planets in widely-separated non-resonant orbits with relatively large eccentricities are formed. Thus, the final orbital configurations of the terrestrial planets around M dwarfs sensitively depend on strength of type-I migration. We also found that large amount of water-ice is delivered by type-I migration from outer regions and final planets near the inner disk edge around M dwarfs are generally abundant in water-ice except for the innermost one that is shielded by the outer planets, unless type-I migration speed is reduced by a factor of more than 100 from that predicted by the linear theory.

preprint2006arXiv

The Separation/Period Gap in the Distribution of Extrasolar Planets around Stars with Masses M >= 1.2 Msolar

The evidence for a shortage of exosolar planets with semimajor axes -1.1 <= log (a/AU) <= -0.2 is investigated. It is shown that this valley results from a gap in the radial distribution of planets, orbiting stars with masses M >= 1.2 Msolar (the high-mass sample, HMS). No underabundance is found for planets orbiting stars with smaller masses. The observational data also indicate that within the HMS population it is preferentially the more massive planets with M sin(i) >= 0.8 M_J that are missing. Monte-Carlo simulations of planet formation and migration are presented that reproduce the observed shortage of planets in the observed radius regime. A dependence on the disk depletion timescale tau_dep is found. The gap is more pronounced for tau_dep = 10^6 - 10^7 yrs than for tau_dep = 3*10^6 - 3*10^7 yrs. This might explain the observed trend with stellar mass if disks around stars with masses M* >= 1.2 Msolar have shorter depletion timescales than those around less massive stars. Possible reasons for such a dependence are a decrease of disk size and an increase of stellar EUV flux with stellar mass.