BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

S. Ida

S. Ida contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.EPastro-phastro-ph.SRastro-ph.IMphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
5topics
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

19 published item(s)

preprint2016arXiv

Analysis of terrestrial planet formation by the Grand Tack model: System architecture and tack location

The Grand Tack model of terrestrial planet formation has emerged in recent years as the premier scenario used to account for several observed features of the inner solar system. It relies on early migration of the giant planets to gravitationally sculpt and mix the planetesimal disc down to ~1 AU, after which the terrestrial planets accrete from material left in a narrow circum-solar annulus. Here we have investigated how the model fares under a range of initial conditions and migration course-change (`tack') locations. We have run a large number of N-body simulations with a tack location of 1.5 AU and 2 AU and tested initial conditions using equal mass planetary embryos and a semi-analytical approach to oligarchic growth. We make use of a recent model of the protosolar disc that takes account of viscous heating, include the full effect of type 1 migration, and employ a realistic mass-radius relation for the growing terrestrial planets. Results show that the canonical tack location of Jupiter at 1.5 AU is inconsistent with the most massive planet residing at 1 AU at greater than 95% confidence. This favours a tack farther out at 2 AU for the disc model and parameters employed. Of the different initial conditions, we find that the oligarchic case is capable of statistically reproducing the orbital architecture and mass distribution of the terrestrial planets, while the equal mass embryo case is not.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effect of turbulence on collisions of dust particles with planetesimals in protoplanetary disks

Planetesimals in gaseous protoplanetary disks may grow by collecting dust particles. Hydrodynamical studies show that small particles generally avoid collisions with the planetesimals because they are entrained by the flow around them. This occurs when $St$, the Stokes number, defined as the ratio of the dust stopping time to the planetesimal crossing time, becomes much smaller than unity. However, these studies have been limited to the laminar case, whereas these disks are believed to be turbulent. We want to estimate the influence of gas turbulence on the dust-planetesimal collision rate and on the impact speeds. We used three-dimensional direct numerical simulations of a fixed sphere (planetesimal) facing a laminar and turbulent flow seeded with small inertial particles (dust) subject to a Stokes drag. A no-slip boundary condition on the planetesimal surface is modeled via a penalty method. We find that turbulence can significantly increase the collision rate of dust particles with planetesimals. For a high turbulence case (when the amplitude of turbulent fluctuations is similar to the headwind velocity), we find that the collision probability remains equal to the geometrical rate or even higher for $St\geq 0.1$, i.e., for dust sizes an order of magnitude smaller than in the laminar case. We derive expressions to calculate impact probabilities as a function of dust and planetesimal size and turbulent intensity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Extrasolar Binary Planets II: Detectability by Transit Observations

We discuss the detectability of gravitationally bounded pairs of gas-giant planets (which we call "binary planets") in extrasolar planetary systems that are formed through orbital instability followed by planet-planet dynamical tides during their close encounters, based on the results of N-body simulations by Ochiai, Nagasawa and Ida (Paper I). Paper I showed that the formation probability of a binary is as much as $\sim 10\%$ for three giant planet systems that undergo orbital instability, and after post-capture long-term tidal evolution, the typical binary separation is 3--5 times the sum of physical radii of the planets. The binary planets are stable during main sequence lifetime of solar-type stars, if the stellarcentric semimajor axis of the binary is larger than 0.3 AU. We show that detecting modulations of transit light curves is the most promising observational method to detect binary planets. Since the likely binary separations are comparable to the stellar diameter, the shape of the transit light curve is different from transit to transit, depending on the phase of the binary's orbit. The transit durations and depth for binary planet transits are generally longer and deeper than those for the single planet case. We point out that binary planets could exist among the known inflated gas giant planets or objects classified as false positive detections at orbital radii > 0.3 AU, propose a binary planet explanation for the CoRoT candidate SRc01 E2 1066, and show that binary planets are likely to be present in, and could be detected using Kepler-quality data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A dynamical study on the habitability of terrestrial exoplanets II: The super Earth HD 40307 g

HARPS and it Kepler results indicate that half of solar-type stars host planets with periods P<100 d and masses M < 30 M_E. These super Earth systems are compact and dynamically cold. Here we investigate the stability of the super Earth system around the K-dwarf HD40307. It could host up to six planets, with one in the habitable zone. We analyse the system&#39;s stability using numerical simulations from initial conditions within the observational uncertainties. The most stable solution deviates 3.1 sigma from the published value, with planets e and f not in resonance and planets b and c apsidally aligned. We study the habitability of the outer planet through the yearly-averaged insolation and black-body temperature at the pole. Both undergo large variations because of its high eccentricity and are much more intense than on Earth. The insolation variations are precession dominated with periods of 40 kyr and 102 kyr for precession and obliquity if the rotation period is 3 d. A rotation period of about 1.5 d could cause extreme obliquity variations because of capture in a Cassini state. For faster rotation rates the periods converge to 10 kyr and 20 kyr. The large uncertainty in the precession period does not change the overall outcome.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Extrasolar Binary Planets I: Formation by tidal capture during planet-planet scattering

We have investigated i) the formation of gravitationally bounded pairs of gas-giant planets (which we call &#34;binary planets&#34;) from capturing each other through planet-planet dynamical tide during their close encounters and ii) the following long-term orbital evolution due to planet-planet and planet-star {\it quasi-static} tides. For the initial evolution in phase i), we carried out N-body simulations of the systems consisting of three jupiter-mass planets taking into account the dynamical tide. The formation rate of the binary planets is as much as 10% of the systems that undergo orbital crossing and this fraction is almost independent of the initial stellarcentric semi-major axes of the planets, while ejection and merging rates sensitively depend on the semi-major axes. As a result of circularization by the planet-planet dynamical tide, typical binary separations are a few times the sum of the physical radii of the planets. After the orbital circularization, the evolution of the binary system is governed by long-term quasi-static tide. We analytically calculated the quasi-static tidal evolution in later phase ii). The binary planets first enter the spin-orbit synchronous state by the planet-planet tide. The planet-star tide removes angular momentum of the binary motion, eventually resulting in a collision between the planets. However, we found that the binary planets survive the tidal decay for main-sequence life time of solar-type stars (~10Gyrs), if the binary planets are beyond ~0.3 AU from the central stars. These results suggest that the binary planets can be detected by transit observations at >0.3AU.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Orbital Circularization of a Planet Accreting Disk Gas: Formation of Distant Jupiters in Circular Orbits based on Core Accretion Model

Recently, gas giant planets in nearly circular orbits with large semimajor axes ($a \sim$ 30--1000AU) have been detected by direct imaging. We have investigated orbital evolution in a formation scenario for such planets, based on core accretion model: i) Icy cores accrete from planetesimals at $\lesssim$ 30AU, ii) they are scattered outward by an emerging nearby gas giant to acquire highly eccentric orbits, and iii) their orbits are circularized through accretion of disk gas in outer regions, where they spend most of time. We analytically derived equations to describe the orbital circularization through the gas accretion. Numerical integrations of these equations show that the eccentricity decreases by a factor of more than 5 during the planetary mass increases by a factor of 10. Because runaway gas accretion increases planetary mass by $\sim$ 10--300, the orbits are sufficiently circularized. On the other hand, $a$ is reduced at most only by a factor of 2, leaving the planets in outer regions. If the relative velocity damping by shock is considered, the circularization is slowed down, but still efficient enough. Therefore, this scenario potentially accounts for the formation of observed distant jupiters in nearly circular orbits. If the apocenter distances of the scattered cores are larger than the disk sizes, their $a$ shrink to a quarter of the disk sizes; the $a$-distribution of distant giants could reflect outer edges of the disks in a similar way that those of hot jupiters may reflect inner edges.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Planet Population Synthesis

With the increasing number of exoplanets discovered, statistical properties of the population as a whole become unique constraints on planet formation models provided a link between the description of the detailed processes playing a role in this formation and the observed population can be established. Planet population synthesis provides such a link. The approach allows to study how different physical models of individual processes (e.g., proto-planetary disc structure and evolution, planetesimal formation, gas accretion, migration, etc.) affect the overall properties of the population of emerging planets. By necessity, planet population synthesis relies on simplified descriptions of complex processes. These descriptions can be obtained from more detailed specialised simulations of these processes. The objective of this chapter is twofold: 1) provide an overview of the physics entering in the two main approaches to planet population synthesis and 2) present some of the results achieved as well as illustrate how it can be used to extract constraints on the models and to help interpret observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Toward a Deterministic Model of Planetary Formation VII: Eccentricity Distribution of Gas Giants

The ubiquity of planets and diversity of planetary systems reveal planet formation encompass many complex and competing processes. In this series of papers, we develop and upgrade a population synthesis model as a tool to identify the dominant physical effects and to calibrate the range of physical conditions. Recent planet searches leads to the discovery of many multiple-planet systems. Any theoretical models of their origins must take into account dynamical interaction between emerging protoplanets. Here, we introduce a prescription to approximate the close encounters between multiple planets. We apply this method to simulate the growth, migration, and dynamical interaction of planetary systems. Our models show that in relatively massive disks, several gas giants and rocky/icy planets emerge, migrate, and undergo dynamical instability. Secular perturbation between planets leads to orbital crossings, eccentricity excitation, and planetary ejection. In disks with modest masses, two or less gas giants form with multiple super-Earths. Orbital stability in these systems is generally maintained and they retain the kinematic structure after gas in their natal disks is depleted. These results reproduce the observed planetary mass-eccentricity and semimajor axis-eccentricity correlations. They also suggest that emerging gas giants can scatter residual cores to the outer disk regions. Subsequent in situ gas accretion onto these cores can lead to the formation of distant (> 30AU) gas giants with nearly circular orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A dynamical study on the habitability of terrestrial exoplanets I: Tidally evolved planet-satellite pairs

We investigate the obliquity and spin period of Earth-Moon like systems after 4.5 Gyr of tidal evolution with various satellite masses and initial planetary obliquity and discuss their relations to the habitability of the planet. We find three possible outcomes: either i) the system is still evolving, ii) the system is double synchronous or iii) the satellite has collided with the planet. The transition between case i) and ii) is abrupt and occurs at slightly larger satellite mass ($m_s \sim 0.02m_p$) than the lunar mass. We suggest that cases ii) and iii) are less habitable than case i). Using results from models of giant impacts and satellite accretion, we found that the systems that mimic our own with rotation period $12 < P_p < 48$ h and current planetary obliquity $\varepsilon_p < 40^\circ$ or $\varepsilon_p > 140^\circ$ only represent 14% of the possible outcomes. Elser et al. (2011) conclude that the probability of a terrestrial planet having a heavy satellite is 13%. Combining these results suggests that the probability of ending up with a system such as our own is of the order of 2%.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Merging Criteria for Giant Impacts of Protoplanets

At the final stage of terrestrial planet formation, known as the giant impact stage, a few tens of Mars-sized protoplanets collide with one another to form terrestrial planets. Almost all previous studies on the orbital and accretional evolution of protoplanets in this stage have been based on the assumption of perfect accretion, where two colliding protoplanets always merge. However, recent impact simulations have shown that collisions among protoplanets are not always merging events, that is, two colliding protoplanets sometimes move apart after the collision (hit-and-run collision). As a first step towards studying the effects of such imperfect accretion of protoplanets on terrestrial planet formation, we investigated the merging criteria for collisions of rocky protoplanets. Using the smoothed particle hydrodynamic (SPH) method, we performed more than 1000 simulations of giant impacts with various parameter sets, such as the mass ratio of protoplanets, $γ$, the total mass of two protoplanets, $M_{\rm T}$, the impact angle, $θ$, and the impact velocity, $v_{\rm imp}$. We investigated the critical impact velocity, $v_{\rm cr}$, at the transition between merging and hit-and-run collisions. We found that the normalized critical impact velocity, $v_{\rm cr}/v_{\rm esc}$, depends on $γ$ and $θ$, but does not depend on $M_{\rm T}$, where $v_{\rm esc}$ is the two-body escape velocity. We derived a simple formula for $v_{\rm cr}/v_{\rm esc}$ as a function of $γ$ and $θ$, and applied it to the giant impact events obtained by \textit{N}-body calculations in the previous studies. We found that 40% of these events should not be merging events.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Orbital Evolution of Planets around Intermediate-Mass Giants

Around low- and intermediate-mass (1.5-3 M_sun) red giants, no planets have been found inside 0.6 AU. Such a paucity is not seen in the case of 1 M_sun main sequence stars. In this study, we examine the possibility that short-period planets were engulfed by their host star evolving off the main sequence. To do so, we have simulated the orbital evolution of planets, including the effects of stellar tide and mass loss, to determine the critical semimajor axis, a_crit, beyond which planets survive the RGB expansion of their host star. We have found that a_crit changes drastically around 2 M_sun: In the lower-mass range, a_crit is more than 1 AU, while a_crit is as small as about 0.2 AU in the higher-mass range. Comparison with measured semimajor axes of known planets suggests that there is a lack of planets that only planet engulfment never accounts for in the higher-mass range. Whether the lack is real affects our understanding of planet formation. Therefore, increasing the number of planet samples around evolved intermediate-mass stars is quite meaningful to confirm robustness of the lack of planets.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Planet Engulfment by ~1.5-3 Solar-Mass Red Giants

Recent radial-velocity surveys for GK clump giants have revealed that planets also exist around ~1.5-3 Msun stars. However, no planets have been found inside 0.6 AU around clump giants, in contrast to solar-type main-sequence stars, many of which harbor short-period planets such as hot Jupiters. In this study we examine the possibility that planets were engulfed by host stars evolving on the red-giant branch (RGB). We integrate the orbital evolution of planets in the RGB and helium burning (HeB) phases of host stars, including the effects of stellar tide and stellar mass loss. Then we derive the critical semimajor axis (or the survival limit) inside which planets are eventually engulfed by their host stars after tidal decay of their orbits. Especially, we investigate the impact of stellar mass and other stellar parameters on the survival limit in more detail than previous studies. In addition, we make detailed comparison with measured semimajor axes of planets detected so far, which no previous study did. We find that the critical semimajor axis is quite sensitive to stellar mass in the range between 1.7 and 2.1 Msun, which suggests a need for careful comparison between theoretical and observational limits of existence of planets. Our comparison demonstrates that all those planets are beyond the survival limit, which is consistent with the planet-engulfment hypothesis. However, on the high-mass side (> 2.1 Msun), the detected planets are orbiting significantly far from the survival limit, which suggests that engulfment by host stars may not be the main reason for the observed lack of short-period giant planets. To confirm our conclusion, the detection of more planets around clump giants, especially with masses > 2.5 Msun, is required.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The SPICA coronagraphic instrument (SCI) for the study of exoplanets

We present the SPICA Coronagraphic Instrument (SCI), which has been designed for a concentrated study of extra-solar planets (exoplanets). SPICA mission provides us with a unique opportunity to make high contrast observations because of its large telescope aperture, the simple pupil shape, and the capability for making infrared observations from space. The primary objectives for the SCI are the direct coronagraphic detection and spectroscopy of Jovian exoplanets in infrared, while the monitoring of transiting planets is another important target. The specification and an overview of the design of the instrument are shown. In the SCI, coronagraphic and non-coronagraphic modes are applicable for both an imaging and a spectroscopy. The core wavelength range and the goal contrast of the coronagraphic mode are 3.5--27$μ$m, and 10$^{-6}$, respectively. Two complemental designs of binary shaped pupil mask coronagraph are presented. The SCI has capability of simultaneous observations of one target using two channels, a short channel with an InSb detector and a long wavelength channel with a Si:As detector. We also give a report on the current progress in the development of key technologies for the SCI.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Population of Very-Hot Super-Earths in Multiple-Planet Systems Should be Uncovered by Kepler

We simulate a Kepler-like observation of a theoretical exoplanet population and we show that the observed orbital period distribution of the Kepler giant planet candidates is best matched by an average stellar specific dissipation function Q_* in the interval 10^6 ~< Q_* ~< 10^7. In that situation, the few super-Earths that are driven to orbital periods P < 1 day by dynamical interactions in multiple-planet systems will survive tidal disruption for a significant fraction of the main-sequence lifetimes of their stellar hosts. Consequently, though these very-hot super-Earths are not characteristic of the overall super-Earth population, their substantial transit probability implies that they should be significant contributors to the full super-Earth population uncovered by Kepler. As a result, the CoRoT-7 system may be the first representative of a population of very-hot super-Earths that we suggest should be found in multiple-planet systems preferentially orbiting the least-dissipative stellar hosts in the Kepler sample.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Toward a Deterministic Model of Planetary Formation VI: Dynamical Interaction and Coagulation of Multiple Rocky Embryos and Super-Earth Systems around Solar Type Stars

Radial velocity and transit surveys indicate that solar-type stars bear super-Earths, with mass and period up to ~ 20 M_E and a few months, are more common than those with Jupiter-mass gas giants. In many cases, these super-Earths are members of multiple-planet systems in which their mutual dynamical interaction has influenced their formation and evolution. In this paper, we modify an existing numerical population synthesis scheme to take into account protoplanetary embryos&#39; interaction with their evolving natal gaseous disk, as well as their close scatterings and resonant interaction with each other. We show that it is possible for a group of compact embryos to emerge interior to the ice line, grow, migrate, and congregate into closely-packed convoys which stall in the proximity of their host stars. After the disk-gas depletion, they undergo orbit crossing, close scattering, and giant impacts to form multiple rocky Earths or super-Earths in non-resonant orbits around ~ 0.1AU with moderate eccentricities of ~0.01-0.1. We suggest that most refractory super-Earths with period in the range of a few days to weeks may have formed through this process. These super-Earths differ from Neptune-like ice giants by their compact sizes and lack of a substantial gaseous envelope.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Modification of Angular Velocity by Inhomogeneous MRI Growth in Protoplanetary Disks

We have investigated evolution of magneto-rotational instability (MRI) in protoplanetary disks that have radially non-uniform magnetic field such that stable and unstable regions coexist initially, and found that a zone in which the disk gas rotates with a super-Keplerian velocity emerges as a result of the non-uniformly growing MRI turbulence. We have carried out two-dimensional resistive MHD simulations with a shearing box model. We found that if the spatially averaged magnetic Reynolds number, which is determined by widths of the stable and unstable regions in the initial conditions and values of the resistivity, is smaller than unity, the original Keplerian shear flow is transformed to the quasi-steady flow such that more flattened (rigid-rotation in extreme cases) velocity profile emerges locally and the outer part of the profile tends to be super-Keplerian. Angular momentum and mass transfer due to temporally generated MRI turbulence in the initially unstable region is responsible for the transformation. In the local super-Keplerian region, migrations due to aerodynamic gas drag and tidal interaction with disk gas are reversed. The simulation setting corresponds to the regions near the outer and inner edges of a global MRI dead zone in a disk. Therefore, the outer edge of dead zone, as well as the inner edge, would be a favorable site to accumulate dust particles to form planetesimals and retain planetary embryos against type I migration.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Signature of the Ice Line and Modest Type I Migration in the Observed Exoplanet Mass-Semimajor Axis Distribution

Existing exoplanet radial velocity surveys are complete in the planetary mass-semimajor axis (Mp-a) plane over the range 0.1 AU < a < 2.0 AU where Mp >~ 100 M_Earth. We marginalize over mass in this complete domain of parameter space and demonstrate that the observed semimajor axis distribution is inconsistent with models of planet formation that use the full Type I migration rate derived from a linear theory and that do not include the effect of the ice line on the disk surface density profile. However, the efficiency of Type I migration can be suppressed by both nonlinear feedback and the barriers introduced by local maxima in the disk pressure distribution, and we confirm that the synthesized Mp-a distribution is compatible with the observed data if we account for both retention of protoplanetary embryos near the ice line and an order-of-magnitude reduction in the efficiency of Type I migration. The validity of these assumption can be checked because they also predict a population of short-period rocky planets with a range of masses comparable to that of the Earth as well as a &#34;desert&#34; in the Mp-a distribution centered around Mp ~ 30-50 M_Earth and a < 1 AU. We show that the expected &#34;desert&#34; in the Mp-a plane will be discernible by a radial velocity survey with 1 m/s precision and n ~ 700 radial velocity observations of program stars.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Towards a deterministic model of planetary formation. III. Mass distribution of short-period planets around stars of various masses

The origin of a recently discovered close-in Neptune-mass planet around GJ436 poses a challenge to the current theories of planet formation. Based on the sequential accretion hypothesis and the standard theory of gap formation and orbital migration, we show that around M dwarf stars, close-in Neptune-mass ice-giant planets may be relatively common, while close-in Jupiter-mass gas-giant planets are relatively rare. The mass distribution of close-in planets generally has two peaks at about Neptune mass and Jupiter mass. The lower-mass peak takes the maximum frequency for M dwarfs. Around more massive solar-type stars (G dwarfs), the higher-mass peak is much more pronounced. These are because planets tend to undergo type II migration after fully accreting gas around G dwarfs while they tend to migrate faster than gas accretion around M stars. Close-in Neptune-mass planets may also exist around G dwarfs, though they tend to be mostly composed of silicates and iron cores and their frequency is expected to be much smaller than that of Neptune-mass planets around M dwarfs and that of gas giants around G dwarfs. We also show that the conditions for planets&#39; migration due to their tidal interaction with the disk and the stellar-mass dependence in the disk-mass distribution can be calibrated by the mass distribution of short-period planets around host stars with various masses.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Spatial Structure and Coherent Motion in Dense Planetary Rings Induced by Self-Gravitational Instability

We investigate the formation of spatial structure in dense, self-gravitating particle systems such as Saturn&#39;s B-ring through local $N$-body simulations to clarify the intrinsic physics based on individual particle motion. In such a system, Salo (1995) showed that the formation of spatial structure such as wake-like structure and particle grouping (clump) arises spontaneously due to gravitational instability and the radial velocity dispersion increases as the formation of the wake structure. However, intrinsic physics of the phenomena has not been clarified. We performed local $N$-body simulations including mutual gravitational forces between ring particles as well as direct (inelastic) collisions with identical (up to $N\sim40000$) particles. In the wake structure particles no longer move randomly but coherently. We found that particle motion was similar to Keplerian motion even in the wake structure and that the coherent motion was produced since the particles in a clump had similar eccentricity and longitude of perihelion. This coherent motion causes the increase and oscillation in the radial velocity dispersion. The mean velocity dispersion is rather larger in a more dissipative case with a smaller restitution coefficient and/or a larger surface density since the coherence is stronger in the more dissipative case. Our simulations showed that the wavelength of the wake structure was approximately given by the longest wavelength $\hsλ{cr} = 4π^2 GΣ/κ^2$ in the linear theory of axisymmetric gravitational instability in a thin disk, where $G$, $Σ$, and $κ$ are the gravitational constant, surface density, and a epicyclic frequency.

0 citations0 reviewsNo code linkedOpen access