Source author record

Masahiro Ogihara

Masahiro Ogihara 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

20works
3topics
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

20 published item(s)

preprint2022arXiv

3-D selection of 167 sub-stellar companions to nearby stars

We analyze 5108 AFGKM stars with at least five high precision radial velocity points as well as Gaia and Hipparcos astrometric data utilizing a novel pipeline developed in previous work. We find 914 radial velocity signals with periods longer than 1000\,d. Around these signals, 167 cold giants and 68 other types of companions are identified by combined analyses of radial velocity, astrometry, and imaging data. Without correcting for detection bias, we estimate the minimum occurrence rate of the wide-orbit brown dwarfs to be 1.3\%, and find a significant brown dwarf valley around 40 $M_{\rm Jup}$. We also find a power-law distribution in the host binary fraction beyond 3 au similar to that found for single stars, indicating no preference of multiplicity for brown dwarfs. Our work also reveals nine sub-stellar systems (GJ 234 B, GJ 494 B, HD 13724 b, HD 182488 b, HD 39060 b and c, HD 4113 C, HD 42581 d, HD 7449 B, and HD 984 b) that have previously been directly imaged, and many others that are observable at existing facilities. Depending on their ages we estimate that an additional 10-57 sub-stellar objects within our sample can be detected with current imaging facilities, extending the imaged cold (or old) giants by an order of magnitude.

preprint2022arXiv

A Super-Earth Orbiting Near the Inner Edge of the Habitable Zone around the M4.5-dwarf Ross 508

We report the near-infrared radial-velocity (RV) discovery of a super-Earth planet on a 10.77-day orbit around the M4.5 dwarf Ross 508 ($J_\mathrm{mag}=9.1$). Using precision RVs from the Subaru Telescope IRD (InfraRed Doppler) instrument, we derive a semi-amplitude of $3.92^{+0.60}_{-0.58}$ ${\rm m\,s}^{-1}$, corresponding to a planet with a minimum mass $m \sin i = 4.00^{+0.53}_{-0.55}\ M_{\oplus}$. We find no evidence of significant signals at the detected period in spectroscopic stellar activity indicators or MEarth photometry. The planet, Ross 508 b, has a semimajor-axis of $0.05366^{+0.00056}_{-0.00049}$ au. This gives an orbit-averaged insolation of $\approx$1.4 times the Earth's value, placing Ross 508 b near the inner edge of its star's habitable zone. We have explored the possibility that the planet has a high eccentricity and its host is accompanied by an additional unconfirmed companion on a wide orbit. Our discovery demonstrates that the near-infrared RV search can play a crucial role to find a low-mass planet around cool M dwarfs like Ross 508.

preprint2022arXiv

Rapid-then-slow migration reproduces mass distribution of TRAPPIST-1 system

The TRAPPIST-1 system is an iconic planetary system in various aspects (e.g., habitability, resonant relation, and multiplicity) and hence has attracted considerable attention. The mass distribution of the TRAPPIST-1 planets is characterized by two features: the two inner planets are large, and the masses of the four planets in the outer orbit increase with orbital distance. The origin of these features cannot be explained by previous formation models. We investigate whether the mass distribution of the TRAPPIST-1 system can be reproduced by a planet formation model using N-body simulations. We used a gas disk evolution model around a low-mass star constructed by considering disk winds and followed the growth and orbital migration from planetary embryos with the isolation mass, which increases with orbital distance. As a result, we find that from the initial phase, planets in inner orbits undergo rapid orbital migration, and the coalescence growth near the inner disk edge is enhanced. This allows the inner planets to grow larger. Meanwhile, compared with the inner planets, planets in outer orbits migrate more slowly and do not frequently collide with neighboring planets. Therefore, the trend of increasing mass toward the outer orbit, called reversed mass ranking, is maintained. The final mass distribution approximately agrees with the two features of the mass distribution in the TRAPPIST-1 system. We discover that the mass distribution in the TRAPPIST-1 system can be reproduced when embryos experience rapid migration and become trapped near the disk inner edge, and then more massive embryos undergo slower migration. This migration transition can be achieved naturally in a disk evolution model with disk winds.

preprint2021arXiv

Elemental abundances of nearby M dwarfs based on high-resolution near-infrared spectra obtained by the Subaru/IRD survey: Proof of concept

Detailed chemical analyses of M dwarfs are scarce but necessary to constrain the formation environment and internal structure of planets being found around them. We present elemental abundances of 13 M dwarfs (2900 < Teff < 3500 K) observed in the Subaru/IRD planet search project. They are mid-to-late M dwarfs whose abundance of individual elements has not been well studied. We use the high-resolution (~70,000) near-infrared (970-1750 nm) spectra to measure the abundances of Na, Mg, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Sr by the line-by-line analysis based on model atmospheres, with typical errors ranging from 0.2 dex for [Fe/H] to 0.3-0.4 dex for other [X/H]. We measure radial velocities from the spectra and combine them with Gaia astrometry to calculate the Galactocentric space velocities UVW. The resulting [Fe/H] values agree with previous estimates based on medium-resolution K-band spectroscopy, showing a wide distribution of metallicity (-0.6 < [Fe/H] < +0.4). The abundance ratios of individual elements [X/Fe] are generally aligned with the solar values in all targets. While the [X/Fe] distributions are comparable to those of nearby FGK stars, most of which belong to the thin disk population, the most metal-poor object, GJ 699, could be a thick disk star. The UVW velocities also support this. The results raise the prospect that near-infrared spectra of M dwarfs obtained in the planet search projects can be used to grasp the trend of elemental abundances and Galactic stellar population of nearby M dwarfs.

preprint2020arXiv

Breaking Resonant Chains: Destabilization of Resonant Planets due to Long-term Mass Evolution

Recent exoplanet observations reported a large number of multiple-planet systems, in which some of the planets are in a chain of resonances. The fraction of resonant systems to non-resonant systems provides clues about their formation history. We investigated the orbital stability of planets in resonant chains by considering the long-term evolution of planetary mass and stellar mass and using orbital calculations. We found that while resonant chains were stable, they can be destabilized by a change of $\sim$10% in planetary mass. Such a mass evolution can occur by atmospheric escape due to photoevaporation. We also found that resonant chains can be broken by a stellar mass loss of $\lesssim1$%, which would be explained by stellar winds or coronal mass ejections. The long-term mass change of planets and stars plays an important role in the orbital evolutions of planetary systems including super-Earths.

preprint2020arXiv

Do the TRAPPIST-1 Planets Have Hydrogen-rich Atmospheres?

Recently, transmission spectroscopy in the atmospheres of the TRAPPIST-1 planets revealed flat and featureless absorption spectra, which rule out cloud-free hydrogen-dominated atmospheres. Earth-sized planets orbiting TRAPPIST-1 likely have either a clear or a cloudy/hazy hydrogen-poor atmosphere. In this paper, we investigate whether a proposed formation scenario is consistent with expected atmospheric compositions of the TRAPPIST-1 planets. We examine the amount of a hydrogen-rich gas that TRAPPIST-1-like planets accreted from the ambient disk until disk dispersal. Since TRAPPIST-1 planets are trapped into a resonant chain, we simulate disk gas accretion onto a migrating TRAPPIST-1-like planet. We find that the amount of an accreted hydrogen-rich gas is as small as 10$^{-2}$ wt% and 0.1 wt% for TRAPPIST-1b and 1c, 10$^{-2}$ wt% for 1d, 1 wt% for 1e, a few wt% for 1f and 1g and 1 wt% for 1h, respectively. We also calculate a long-term thermal evolution of TRAPPIST-1-like planets after disk dissipation and estimate the mass loss of their hydrogen-rich atmospheres driven by a stellar X-ray and UV irradiation. We find that all the accreted hydrogen-rich atmospheres can be lost via hydrodynamic escape. Therefore, we conclude that TRAPPIST-1 planets should have no primordial hydrogen-rich gases but secondary atmospheres such as a Venus-like one and water vapor, if they currently retain atmospheres.

preprint2020arXiv

Formation of single-moon systems around gas giants

Several mechanisms have been proposed to explain the formation process of satellite systems, and relatively large moons are thought to be born in circumplanetary disks. Making a single-moon system is known to be more difficult than multiple-moon or moonless systems. We aim to find a way to form a system with a single large moon, such as Titan around Saturn. We examine the orbital migration of moons, which change their direction and speed depending on the properties of circumplanetary disks. We modeled dissipating circumplanetary disks with taking the effect of temperature structures into account and calculated the orbital evolution of Titan-mass satellites in the final evolution stage of various circumplanetary disks. We also performed N-body simulations of systems that initially had multiple satellites to see whether single-moon systems remained at the end. The radial slope of the disk-temperature structure characterized by the dust opacity produces a patch of orbits in which the Titan-mass moons cease inward migration and even migrate outward in a certain range of the disk viscosity. The patch assists moons initially located in the outer orbits to remain in the disk, while those in the inner orbits fall onto the planet. We demonstrate for the first time that systems can form that have only one large moon around giant planet. Our N-body simulations suggest satellite formation was not efficient in the outer radii of circumplanetary disks.

preprint2020arXiv

Promoted Mass Growth of Multiple, Distant Giant Planets through Pebble Accretion and Planet-Planet Collision

We propose a pebble-driven planet formation scenario to form giant planets with high multiplicity and large orbital distances in the early gas disk phase. We perform N-body simulations to investigate the growth and migration of low-mass protoplanets in the disk with inner viscously heated and outer stellar irradiated regions. The key feature of this model is that the giant planet cores grow rapidly by a combination of pebble accretion and planet-planet collisions. This consequently speeds up their gas accretion. Because of efficient growth, the planet transitions from rapid type I migration to slow type II migration early, reducing the inward migration substantially. Multiple giant planets can sequentially form in this way with increasing semimajor axes. Both mass growth and orbital retention are more pronounced when a large number of protoplanets are taken into account compared to the case of single planet growth. Eventually, a few numbers of giant planets form with orbital distances of a few to a few tens of AUs within $1.5{-}3$ Myr after the birth of the protoplanets. The resulting simulated planet populations could be linked to the substructures exhibited in disk observations as well as large orbital distance exoplanets observed in radial velocity and microlensing surveys.

preprint2020arXiv

Unified simulations of planetary formation and atmospheric evolution II: Rapid disk clearing by photoevaporation yields low-mass super-Earth atmospheres

Super-Earths possess low-mass H$_2$/He atmospheres (typically less than 10% by mass). However, the origins of super-Earth atmospheres have not yet been ascertained. We investigate the role of rapid disk clearing by photoevaporation during the formation of super-Earths and their atmospheres. We perform unified simulations of super-Earth formation and atmospheric evolution in evolving disks that consider both photoevaporative winds and magnetically driven disk winds. For the growth mode of planetary cores, we consider two cases in which planetary embryos grow with and without pebble accretion. Our main findings are summarized as follows. (i) The time span of atmospheric accretion is shortened by rapid disk dissipation due to photoevaporation, which prevents super-Earth cores from accreting massive atmospheres. (ii) Even if planetary cores grow rapidly by embryo accretion in the case without pebble accretion, the onset of runaway gas accretion is delayed because the isolation mass for embryo accretion is small. Together with rapid disk clearing, the accretion of massive atmospheres can be avoided. (iii) After rapid disk clearing, a number of high-eccentricity embryos can remain in outer orbits. Thereafter, such embryos may collide with the super-Earths, leading to efficient impact erosion of accreted atmospheres. We, therefore, find that super-Earths with low-mass H$_2$/He atmospheres are naturally produced by N-body simulations that consider realistic disk evolution.

preprint2020arXiv

Unified simulations of planetary formation and atmospheric evolution: Effects of pebble accretion, giant impacts, and stellar irradiations on super-Earth formation

A substantial number of super-Earths have been discovered, and atmospheres of transiting super-Earths have also been observed by transmission spectroscopy. Several lines of observational evidence indicate that most super-Earths do not possess massive H$_2$/He atmospheres. However, accretion and retention of less massive atmospheres on super-Earths challenge planet formation theory. We consider the following three mechanisms: (i) envelope heating by pebble accretion, (ii) mass loss during giant impacts, and (iii) atmospheric loss by stellar X-ray and EUV photoevaporation. We investigate whether these mechanisms influence the amount of the atmospheres that form around super-Earths. We develop a code combining an N-body simulation of pebble-driven planetary formation and an atmospheric evolution simulation. We demonstrate that the observed orbital properties of super-Earths are well reproduced by the results of our simulations. However, (i) heating by pebble accretion ceases prior to disk dispersal, (ii) the frequency of giant impact events is too low to sculpt massive atmospheres, and (iii) many super-Earths having H$_2$/He atmospheres of $\gtrsim 10$ wt% survive against stellar irradiations for 1 Gyr. Therefore, it is likely that other mechanisms such as suppression of gas accretion are required to explain less massive atmospheres ($\lesssim 10$ wt%) of super-Earths.

preprint2016arXiv

Evolution of Protoplanetary Discs with Magnetically Driven Disc Winds

Aims: We investigate the evolution of protoplanetary discs (PPDs hereafter) with magnetically driven disc winds and viscous heating. Methods: We consider an initially massive disc with ~0.1 Msun to track the evolution from the early stage of PPDs. We solve the time evolution of surface density and temperature by taking into account viscous heating and the loss of the mass and the angular momentum by the disc winds within the framework of a standard alpha model for accretion discs. Our model parameters, turbulent viscosity, disc wind mass loss, and disc wind torque, which are adopted from local magnetohydrodynamical simulations and constrained by the global energetics of the gravitational accretion, largely depends on the physical condition of PPDs, particularly on the evolution of the vertical magnetic flux in weakly ionized PPDs. Results: Although there are still uncertainties concerning the evolution of the vertical magnetic flux remaining, surface densities show a large variety, depending on the combination of these three parameters, some of which are very different from the surface density expected from the standard accretion. When a PPD is in a "wind-driven accretion" state with the preserved vertical magnetic field, the radial dependence of the surface density can be positive in the inner region <1-10 au. The mass accretion rates are consistent with observations, even in the very low level of magnetohydrodynamical turbulence. Such a positive radial slope of the surface density gives a great impact on planet formation because (i)it inhibits the inward drift or even results in the outward drift of pebble/boulder-sized solid bodies, and (ii) it also makes the inward type-I migration of proto-planets slower or even reversed. Conclusions: The variety of our calculated PPDs should yield a wide variety of exoplanet systems.

preprint2015arXiv

A reassessment of the in situ formation of close-in super-Earths

A large fraction of stars host one or multiple close-in super-Earth planets. There is an active debate about whether these planets formed in situ or at greater distances from the central star and migrated to their current position. It has been shown that part of their observed properties (e.g., eccentricity distribution) can be reproduced by N-body simulations of in situ formation starting with a population of protoplanets of high masses and neglecting the effects of the disk gas. We plan to reassess the in situ formation of close-in super-Earths through more complete simulations. We performed N-body simulations of a population of small planetary embryos and planetesimals that include the effects of disk-planet interactions (e.g., eccentricity damping, type I migration). In addition, we also consider the accretion of a primitive atmosphere from a protoplanetary disk. We find that planetary embryos grow very quickly well before the gas dispersal, and thus undergo rapid inward migration, which means that one cannot neglect the effects of a gas disk when considering the in-situ formation of close-in super-Earths. Owing to their rapid inward migration, super-Earths reach a compact configuration near the disk's inner edge whose distribution of orbital parameters matches the observed close-in super-Earths population poorly. On the other hand, simulations including eccentricity damping, but no type I migration, reproduce the observed distributions better. Including the accretion of an atmosphere does not help reproduce the bulk architecture of observations. Interestingly, we find that the massive embryos can migrate inside the disk edge while capturing only a moderately massive hydrogen/helium atmosphere. By this process they avoid becoming giant planets. The bulk of close-in super-Earths cannot form in situ, unless type I migration is suppressed in the entire disk inside 1 AU.

preprint2015arXiv

Formation of terrestrial planets in disks evolving via disk winds and implications for the origin of the solar system's terrestrial planets

Recent three-dimensional magnetohydrodynamical simulations have identified a disk wind by which gas materials are lost from the surface of a protoplanetary disk, which can significantly alter the evolution of the inner disk and the formation of terrestrial planets. A simultaneous description of the realistic evolution of the gaseous and solid components in a disk may provide a clue for solving the problem of the mass concentration of the terrestrial planets in the solar system. We simulate the formation of terrestrial planets from planetary embryos in a disk that evolves via magnetorotational instability and a disk wind. The aim is to examine the effects of a disk wind on the orbital evolution and final configuration of planetary systems. We perform N-body simulations of sixty 0.1 Earth-mass embryos in an evolving disk. The evolution of the gas surface density of the disk is tracked by solving a one-dimensional diffusion equation with a sink term that accounts for the disk wind. We find that even in the case of a weak disk wind, the radial slope of the gas surface density of the inner disk becomes shallower, which slows or halts the type I migration of embryos. If the effect of the disk wind is strong, the disk profile is significantly altered (e.g., positive surface density gradient, inside-out evacuation), leading to outward migration of embryos inside ~ 1 AU. Disk winds play an essential role in terrestrial planet formation inside a few AU by changing the disk profile. In addition, embryos can undergo convergent migration to ~ 1 AU in certainly probable conditions. In such a case, the characteristic features of the solar system's terrestrial planets (e.g., mass concentration around 1 AU, late giant impact) may be reproduced.

preprint2015arXiv

Suppression of type I migration by disk winds

Planets less massive than Saturn tend to rapidly migrate inward in protoplanetary disks. This is the so-called type I migration. Simulations attempting to reproduce the observed properties of exoplanets show that type I migration needs to be significantly reduced over a wide region of the disk for a long time. However, the mechanism capable of suppressing type I migration over a wide region has remained elusive. The recently found turbulence-driven disk winds offer new possibilities. We investigate the effects of disk winds on the disk profile and type I migration for a range of parameters that describe the strength of disk winds. We also examine the in situ formation of close-in super-Earths in disks that evolve through disk winds. The disk profile, which is regulated by viscous diffusion and disk winds, was derived by solving the diffusion equation. We carried out a number of simulations and plot here migration maps that indicate the type I migration rate. We also performed N-body simulations of the formation of close-in super-Earths from a population of planetesimals and planetary embryos. We define a key parameter, Kw, which determines the ratio of strengths between the viscous diffusion and disk winds. For a wide range of Kw, the type I migration rate is presented in migration maps. These maps show that type I migration is suppressed over the whole close-in region when the effects of disk winds are relatively strong (Kw < 100). From the results of N-body simulations, we see that type I migration is significantly slowed down assuming Kw = 40. We also show that the results of N-body simulations match statistical orbital distributions of close-in super-Earths.

preprint2014arXiv

N-body Simulations of Terrestrial Planet Formation under the Influence of a Hot Jupiter

We investigate the formation of multiple-planet systems in the presence of a hot Jupiter using extended N-body simulations that are performed simultaneously with semi-analytic calculations. Our primary aims are to describe the planet formation process starting from planetesimals using high-resolution simulations, and to examine the dependences of the architecture of planetary systems on input parameters (e.g., disk mass, disk viscosity). We observe that protoplanets that arise from oligarchic growth and undergo type I migration stop migrating when they join a chain of resonant planets outside the orbit of a hot Jupiter. The formation of a resonant chain is almost independent of our model parameters, and is thus a robust process. At the end of our simulations, several terrestrial planets remain at around 0.1 AU. The formed planets are not equal-mass; the largest planet constitutes more than 50 percent of the total mass in the close-in region, which is also less dependent on parameters. In the previous work of this paper (Ogihara et al. 2013), we have found a new physical mechanism of induced migration of the hot Jupiter, which is called a crowding-out. If the hot Jupiter opens up a wide gap in the disk (e.g., owing to low disk viscosity), crowding-out becomes less efficient and the hot Jupiter remains. We also discuss angular momentum transfer between the planets and disk.

preprint2013arXiv

Condition for Capture into First-order Mean Motion Resonances and Application to Constraints on Origin of Resonant Systems

We investigate the condition for capture into first-order mean motion resonances using numerical simulations with a wide range of various parameters. In particular, we focus on deriving the critical migration timescale for capture into the 2:1 resonance; additional numerical experiments for closely spaced resonances (e.g., 3:2) are also performed. We find that the critical migration timescale is determined by the planet-to-stellar mass ratio, and its dependence exhibits power-law behavior with index -4/3. This dependence is also supported by simple analytic arguments. We also find that the critical migration timescale for systems with equal-mass bodies is shorter than that in the restricted problem; for instance, for the 2:1 resonance between two equal-mass bodies, the critical timescale decreases by a factor of 10. In addition, using the obtained formula, the origin of observed systems that include first-order commensurabilities is constrained. Assuming that pairs of planets originally form well separated from each other and then undergo convergent migration and are captured in resonances, it is possible that a number of exoplanets experienced rapid orbital migration. For systems in closely spaced resonances, the differential migration timescale between the resonant pair can be constrained well; it is further suggested that several exoplanets underwent migration that can equal or even exceed the type I migration rate predicted by the linear theory. This implies that some of them may have formed in situ. Future observations and the use of our model will allow us to statistically determine the typical migration speed in a protoplanetary disk.

preprint2013arXiv

Crowding Out of Giants by Dwarfs: An Origin for the Lack of Companion Planets in Hot Jupiter Systems

We investigate formation of close-in terrestrial planets from planetary embryos under the influence of a hot Jupiter (HJ) using gravitational N-body simulations that include gravitational interactions between the gas disk and the terrestrial planet (e.g., type I migration). Our simulations show that several terrestrial planets efficiently form outside the orbit of the HJ, making a chain of planets, and all of them gravitationally interact directly or indirectly with the HJ through resonance, which leads to inward migration of the HJ. We call this mechanism of induced migration of the HJ as "crowding out." The HJ is eventually lost by collision with the central star, and only several terrestrial planets remain. We also find that the efficiency of the crowding-out effect depends on model parameters; for example, the heavier the disk is, the more efficient the crowding out is. When planet formation occurs in a massive disk, the HJ can be lost to the central star and is never observed. On the other hand, for a less massive disk, the HJ and terrestrial planets can coexist; however, the companion planets can be below the detection limit of current observations. In both cases, systems with the HJ and terrestrial planets have little chance for detection. Therefore, our model naturally explains the lack of companion planets in HJ systems regardless of the disk mass. In effect, our model provide a theoretical prediction for future observations; additional planets can be discovered just outside the HJ, and their masses should generally be small.

preprint2012arXiv

N-body Simulations of Satellite Formation around Giant Planets: Origin of Orbital Configuration of the Galilean Moons

As the number of discovered extrasolar planets has been increasing, diversity of planetary systems requires studies of new formation scenarios. It is important to study satellite formation in circumplanetary disks, which is often viewed as analogous to formation of rocky planets in protoplanetary disks. We investigated satellite formation from satellitesimals around giant planets through N-body simulations that include gravitational interactions with a circumplanetary gas disk. Our main aim is to reproduce the observable properties of the Galilean satellites around Jupiter through numerical simulations, as previous N-body simulations have not explained the origin of the resonant configuration. We performed accretion simulations based on the work of Sasaki et al. (2010), in which an inner cavity is added to the model of Canup & Ward (2002, 2006). We found that several satellites are formed and captured in mutual mean motion resonances outside the disk inner edge and are stable after rapid disk gas dissipation, which explains the characteristics of the Galilean satellites. In addition, owing to the existence of the disk edge, a radial compositional gradient of the Galilean satellites can also be reproduced. An additional objective of this study is to discuss orbital properties of formed satellites for a wide range of conditions by considering large uncertainties in model parameters. Through numerical experiments and semianalytical arguments, we determined that if the inner edge of a disk is introduced, a Galilean-like configuration in which several satellites are captured into a 2:1 resonance outside the disk inner cavity is almost universal. In fact, such a configuration is produced even for a massive disk and rapid type I migration. This result implies the inevitability of a Galilean satellite formation in addition to providing theoretical predictions for extrasolar satellites.

preprint2010arXiv

Eccentricity Trap: Trapping of Resonantly Interacting Planets near the Disk Inner Edge

Using orbital integration and analytical arguments, we have found a new mechanism (an "eccentricity trap") to halt type I migration of planets near the inner edge of a protoplanetary disk. Because asymmetric eccentricity damping due to disk-planet interaction on the innermost planet at the disk edge plays a crucial role in the trap, this mechanism requires continuous eccentricity excitation and hence works for a resonantly interacting convoy of planets. This trap is so strong that the edge torque exerted on the innermost planet can completely halt type I migrations of many outer planets through mutual resonant perturbations. Consequently, the convoy stays outside the disk edge, as a whole. We have derived semi-analytical formula for the condition for the eccentricity trap and predict how many planets are likely to be trapped. We found that several planets or more should be trapped by this mechanism in protoplanetary disks that have cavities. It can be responsible for the formation of non-resonant, multiple, close-in super-Earth systems extending beyond 0.1AU. Such systems are being revealed by radial velocity observations to be quite common around solar-type stars.

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.