Source author record

Ryuji Morishima

Ryuji Morishima 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

11works
2topics
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

11 published item(s)

preprint2022arXiv

Disruption of Saturn's ring particles by thermal stress

Spacecraft and ground-based observations show that the main rings of Saturn lack particles larger than 10 m. Tidal or collisional destruction of satellites/comets have been proposed as the origin of the main rings; however, Saturn's tide alone cannot grind km-sized fragments into submeter-sized particles because of the high mechanical strength of water ice and rock. The question arises as to why such large particles are not left in the current ring. It is known that thermal stress induced by diurnal and seasonal temperature variations can cause weathering and fragmentation of boulders and contribute to dust and regolith production on the Moon and terrestrial planets, and then such thermal stress can break particles larger than a critical radius while cannot smaller than the critical radius. In this study, we examined the role of thermal stress acting on Saturn's ring particles. We found that thermal stress can grind porous ring particles larger than 10-20 m, which explains the lack of particles larger than 10 m in Saturn's ring. Also, fragmentation by thermal stress can be adoptable for the Epsilon rings of Uranus. Furthermore, thermal stress caused by diurnal or seasonal temperature variation acting on boulders on surfaces of icy satellites and asteroids may play an important role in the evolution of their sizes. Our calculations explain the lack of boulders on icy satellites, except in the geologically active provinces such as the tiger stripes of Enceladus, where boulders are supplied by recent geological activity. We predict that future observations can find numerous boulders around Europa's geologically active cracks.

preprint2016arXiv

Onset of oligarchic growth and implication for accretion histories of dwarf planets

We investigate planetary accretion that starts from equal-mass planetesimals using an analytic theory and numerical simulations. We particularly focus on how the planetary mass $M_{\rm oli}$ at the onset of oligarchic growth depends on the initial mass $m_0$ of a planetesimal. Oligarchic growth commences when the velocity dispersion relative to the Hill velocity of the protoplanet takes its minimum. We find that if $m_0$ is small enough, this normalized velocity dispersion becomes as low as unity during the intermediate stage between the runaway and oligarchic growth stages. In this case, $M_{\rm oli}$ is independent of $m_0$. If $m_0$ is large, on the other hand, oligarchic growth commences directly after runaway growth, and $M_{\rm oli} \propto m_0^{3/7}$. The planetary mass $M_{\rm oli}$ for the solid surface density of the Minimum Mass Solar Nebula is close to the masses of the dwarf planets in a reasonable range of $m_0$. This indicates that they are likely to be the largest remnant planetesimals that failed to become planets. The power-law exponent $q$ of the differential mass distribution of remnant planetesimals is typically $-2.0$ and $-2.7$ to $-2.5$ for small and large $m_0$. The slope, $q \simeq -2.7$, and the bump at $10^{21}$ g (or 50 km in radius) for the mass distribution of hot Kuiper belt objects are reproduced if $m_0$ is the bump mass. On the other hand, small initial planetesimals with $m_0 \sim 10^{13}$ g or less are favored to explain the slope of large asteroids, $q \simeq -2.0$, while the bump at $10^{21}$ g can be reproduced by introducing a small number of asteroid seeds each with mass of $10^{19} $ g.

preprint2015arXiv

A particle-based hybrid code for planet formation

We introduce a new particle-based hybrid code for planetary accretion. The code uses an $N$-body routine for interactions with planetary embryos while it can handle a large number of planetesimals using a super-particle approximation, in which a large number of small planetesimals are represented by a small number of tracers. Tracer-tracer interactions are handled by a statistical routine which uses the phase-averaged stirring and collision rates. We compare hybrid simulations with analytic predictions and pure $N$-body simulations for various problems in detail and find good agreements for all cases. The computational load on the portion of the statistical routine is comparable to or less than that for the $N$-body routine. The present code includes an option of hit-and-run bouncing but not fragmentation, which remains for future work.

preprint2015arXiv

Incomplete cooling down of Saturn's A ring at solar equinox: Implication for seasonal thermal inertia and internal structure of ring particles

At the solar equinox in August 2009, the Composite Infrared Spectrometer (CIRS) onboard Cassini showed the lowest Saturn's ring temperatures ever observed. Detailed radiative transfer models show that the observed equinox temperatures of Saturn's A ring are much higher than model predictions as long as only the flux from Saturn is taken into account. This indicates that the A ring was not completely cooled down at the equinox. We develop a simple seasonal model for ring temperatures and first assume that the internal density and the thermal inertia of a ring particle are uniform with depth. The particle size is estimated to be 1-2 m. The seasonal thermal inertia is found to be 30-50 Jm$^{-2}$K$^{-1}$s$^{-1/2}$ in the middle A ring whereas it is $\sim$ 10 Jm$^{-2}$K$^{-1}$s$^{-1/2}$ or as low as the diurnal thermal inertia in the inner and outermost regions of the A ring. An additional internal structure model, in which a particle has a high density core surrounded by a fluffy regolith mantle, shows that the core radius relative to the particle radius is about 0.9 for the middle A ring and is much less for the inner and outer regions of the A ring. This means that the radial variation of the internal density of ring particles exists across the A ring. Some mechanisms may be confining dense particles in the middle A ring against viscous diffusion. Alternatively, the (middle) A ring might have recently formed ($<$ 10$^{8}$ yr) by destruction of an icy satellite, so that dense particles have not yet diffused over the A ring and regolith mantles of particles have not grown thick. Our model results also indicate that the composition of the core is predominantly water ice, not rock.

preprint2014arXiv

Terrestrial Planet Formation at Home and Abroad

We review the state of the field of terrestrial planet formation with the goal of understanding the formation of the inner Solar System and low-mass exoplanets. We review the dynamics and timescales of accretion from planetesimals to planetary embryos and from embryos to terrestrial planets. We discuss radial mixing and water delivery, planetary spins and the importance of parameters regarding the disk and embryo properties. Next, we connect accretion models to exoplanets. We first explain why the observed hot Super Earths probably formed by in situ accretion or inward migration. We show how terrestrial planet formation is altered in systems with gas giants by the mechanisms of giant planet migration and dynamical instabilities. Standard models of terrestrial accretion fail to reproduce the inner Solar System. The "Grand Tack" model solves this problem using ideas first developed to explain the giant exoplanets. Finally, we discuss whether most terrestrial planet systems form in the same way as ours, and highlight the key ingredients missing in the current generation of simulations.

preprint2013arXiv

N-body simulations of oligarchic growth of Mars: Implications for Hf-W chronology

Dauphas and Pourmand (2011) [Nature 473, 489--492] estimated the accretion timescale of Mars to be 1.8 $^{+0.9}_{-1.0}$ Myr from the W isotopes of martian meteorites. This timescale was derived assuming perfect metal-silicate equilibration between the impactor and the target's mantle. However, in the case of a small impactor most likely only a fraction of the target's mantle is involved in the equilibration, while only a small part of the impactor's core equilibrates in the case of a giant impact. We examined the effects of imperfect equilibration using results of high-resolution $N$-body simulations for the oligarchic growth stage. These effects were found to be small as long as a planetary embryo has a deep liquid magma ocean during its accretion. The effect due to partial involvement of the target's mantle in equilibration is small due to the low metal-silicate partition coefficient for W suggested from the low Hf/W ratio of the martian mantle. The effect due to partial involvement of the impactor's core is also small because a large fraction of the embryo mass is delivered from small planetesimals, which are likely to fully equilibrate in the deep magma ocean on the embryo. The accretion timescale of Mars estimated by the Hf-W chronology is shorter than that expected for the minimum mass solar nebula model as long as more than 10% of each impactor's core re-equilibrates with the martian mantle and the final stages of accretion are prolonged. This probably indicates that accretion of Mars proceeded rapidly due to solid and gas surface densities significantly larger than those for the minimum mass solar nebula or due to accretion of small fragments or pebbles.

preprint2012arXiv

A multilayer model for thermal infrared emission of Saturn's rings. III: Thermal inertia inferred from Cassini CIRS

The thermal inertia values of Saturn's main rings (the A, B, and C rings and the Cassini division) are derived by applying our thermal model to azimuthally scanned spectra taken by the Cassini Composite Infrared Spectrometer (CIRS). Model fits show the thermal inertia of ring particles to be 16, 13, 20, and 11 Jm$^{-2}$K$^{-1}$s$^{-1/2}$ for the A, B, and C rings, and the Cassini division, respectively. However, there are systematic deviations between modeled and observed temperatures in Saturn's shadow depending on solar phase angle, and these deviations indicate that the apparent thermal inertia increases with solar phase angle. This dependence is likely to be explained if large slowly spinning particles have lower thermal inertia values than those for small fast spinning particles because the thermal emission of slow rotators is relatively stronger than that of fast rotators at low phase and vise versa. Additional parameter fits, which assume that slow and fast rotators have different thermal inertia values, show the derived thermal inertia values of slow (fast) rotators to be 8 (77), 8 (27), 9 (34), 5 (55) Jm$^{-2}$K$^{-1}$s$^{-1/2}$ for the A, B, and C rings, and the Cassini division, respectively. The values for fast rotators are still much smaller than those for solid ice with no porosity. Thus, fast rotators are likely to have surface regolith layers, but these may not be as fluffy as those for slow rotators, probably because the capability of holding regolith particles is limited for fast rotators due to the strong centrifugal force on surfaces of fast rotators. Other additional parameter fits, in which radii of fast rotators are varied, indicate that particles less than $\sim$ 1 cm should not occupy more than a half of the cross section for the A, B, and C rings.

preprint2012arXiv

Regolith grain sizes of Saturn's rings inferred from Cassini-CIRS far-infrared spectra

We analyze far-infrared (10-650 cm$^{-1}$) emissivity spectra of Saturn's main rings obtained by the Cassini Composite Infrared Spectrometer (CIRS). In modeling of the spectra, the single scattering albedos of regolith grains are calculated using the Mie theory, diffraction is removed with the delta-Eddington approximation, and the hemispherical emissivities of macroscopic free-floating ring particles are calculated using the Hapke's isotropic scattering model. Only pure crystalline water ice is considered and the size distribution of regolith grains is estimated. We find that good fits are obtained if the size distribution is broad ranging from 1 $μ$m to 1-10 cm with a power law index of $ \sim 3$. This means that the largest regolith grains are comparable to the smallest free-floating particles in size and that the power law indices for both free-floating particles and regolith grains are similar to each other. The apparent relative abundance of small grains increases with decreasing solar phase angle (or increasing mean temperature). This trend is particularly strong for the C ring and is probably caused by eclipse cooling in Saturn's shadow, which relatively suppresses warming up of grains larger than the thermal skin depth ($\sim$ 1 mm) under subsequent solar illumination.

preprint2011arXiv

Gap opening beyond dead zones by photoevaporation

We propose a new hypothesis for the origin of protoplanetary discs with large inner holes (or gaps), so-called transition discs. Our gas disc model takes into account layered accretion, in which poorly-ionized low-viscosity dead zones are sandwiched by high-viscosity surface layers, and photoevaporative winds induced by X-rays from the central stars. We find that a gap opens at a radius outside a dead zone, if the mass loss rate due to photoevaporative winds exceeds the mass accretion rate in the dead zone region. Since the dead zone survives even after the gap opens, mass accretion onto the central star continues for a long time. This model can reproduce large gap sizes and high mass accretion rates seen in observed transition discs.

preprint2011arXiv

How common are Earth-Moon planetary systems?

The Earth's comparatively massive moon, formed via a giant impact on the proto-Earth, has played an important role in the development of life on our planet, both in the history and strength of the ocean tides and in stabilizing the chaotic spin of our planet. Here we show that massive moons orbiting terrestrial planets are not rare. A large set of simulations by Morishima et al., 2010, where Earth-like planets in the habitable zone form, provides the raw simulation data for our study. We use limits on the collision parameters that may guarantee the formation of a circumplanetary disk after a protoplanet collision that could form a satellite and study the collision history and the long term evolution of the satellites qualitatively. In addition, we estimate and quantify the uncertainties in each step of our study. We find that giant impacts with the required energy and orbital parameters for producing a binary planetary system do occur with more than 1 in 12 terrestrial planets hosting a massive moon, with a low-end estimate of 1 in 45 and a high-end estimate of 1 in 4.

preprint2010arXiv

From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5--4.0 AU and the collisional growth is followed until the completion of planetary accretion (> 100 Myr). All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas, and orbits of Jupiter and Saturn. We end up with one to five planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today's eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts.