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.

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.

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.

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.

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.

The core-accretion model predicts that planetary cores as massive as super-Earths undergo runaway gas accretion to become gas giants. However, the exoplanet census revealed the prevalence of super-Earths close to their host stars, which should have avoided runaway gas accretion. In fact, mass-radius relationships of transiting planets suggest that some close-in super-Earths possess H$_2$/He atmospheres of ~ 0.1-10% by mass. Previous studies indicated that properties of a disk gas such as metallicity and the inflow/outflow cycle of a disk gas around a super-Earth can regulate accumulation of a H$_2$/He atmosphere onto itself. In this paper, we propose a new mechanism that radial mass accretion in a disk can limit the gas accretion onto super-Earth cores. Recent magneto-hydrodynamic simulations found that magnetically driven disk winds can drive a rapid gas flow near the disk surface. Such a rapid gas flow may slip out of a planetary core and regulate gas supply to an accreting gas onto the core. We performed N-body simulations for formation of super-Earths with accretion of atmospheres in a viscous accretion disk including effects of wind-driven accretion. We found that even super-Earth cores can avoid triggering runaway gas accretion if the inflow of a disk gas toward the cores is limited by viscous accretion. Our model predicts that super-Earths having H$_2$/He atmosphere of ~ 0.1-10 wt % form within $\lesssim$ 1 au of the central star, whereas gas giants are born in the outer region. This mechanism can explain the radial dependence of observed giant planets beyond the solar system.

No planets exist inside the orbit of Mercury and the terrestrial planets of the solar system exhibit a localized configuration. According to thermal structure calculation of protoplanetary disks, a silicate condensation line (~ 1300 K) is located around 0.1 au from the Sun except for the early phase of disk evolution, and planetesimals could have formed inside the orbit of Mercury. A recent study of disk evolution that includes magnetically driven disk winds showed that the gas disk obtains a positive surface density slope inside ~ 1 au from the central star. In a region with positive midplane pressure gradient, planetesimals undergo outward radial drift. We investigate the radial drift of planetesimals and type I migration of planetary embryos in a disk that viscously evolves with magnetically driven disk winds. We show a case in which no planets remain in the close-in region. Radial drifts of planetesimals are simulated using a recent disk evolution model that includes effects of disk winds. The late stage of planet formation is also examined by performing N-body simulations of planetary embryos. We demonstrate that in the middle stage of disk evolution, planetesimals can undergo convergent radial drift in a magnetorotational instability (MRI)-inactive disk, in which the pressure maximum is created, and accumulate in a narrow ring-like region with an inner edge at ~ 0.7 au from the Sun. We also show that planetary embryos that may grow from the narrow planetesimal ring do not exhibit significant type I migration in the late stage of disk evolution. The origin of the localized configuration of the terrestrial planets of the solar system, in particular the deficit of close-in planets, can be explained by the convergent radial drift of planetesimals in disks with a positive pressure gradient in the close-in region.

Planets with masses larger than about 0.1 Earth-masses undergo rapid inward migration (type I migration) in a standard protoplanetary disk. Recent magnetohydrodynamical simulations revealed the presence of magnetically driven disk winds, which would alter the disk profile and the type I migration in the close-in region. We investigate orbital evolution of planetary embryos in disks that viscously evolve under the effects of disk winds. The aim is to discuss effects of altered disk profiles on type I migration. In addition, we aim to examine whether observed distributions of close-in super-Earths can be reproduced by simulations that include effects of disk winds. We perform N-body simulations of super-Earth formation from planetary embryos, in which a recent model for disk evolution is used. We explore a wide range of parameters and draw general trends. We also carry out N-body simulations of close-in super-Earth formation from embryos in such disks under various conditions. We find that the type I migration is significantly suppressed in many cases. Even in cases in which inward migration occurs, the migration timescale is lengthened to 1 Myr, which mitigates the type I migration problem. This is because the gas surface density is decreased and has a flatter profile in the close-in region due to disk winds. We find that when the type I migration is significantly suppressed, planets undergo late orbital instability during the gas depletion, leading to a non-resonant configuration. We also find that observed distributions of close-in super-Earths (e.g., period ratio, mass ratio) can be reproduced. In addition, we show that in some results of simulations, systems with a chain of resonant planets, like the TRAPPIST-1 system, form.

Magnetically-driven disk winds would alter the surface density slope of gas in the inner region of a protoplanetary disk $(r \lesssim 1 {\rm au})$. This in turn affects planet formation. Recently, the effect of disk wind torque has been considered with the suggestion that it would carve out the surface density of the disk from inside and would induce global gas flows (wind-driven accretion). We aim to investigate effects of global gas flows on type I migration and also examine planet formation. A simplified approach was taken to address this issue, and N-body simulations with isolation-mass planets were also performed. In previous studies, the effect of gas flow induced by turbulence-driven accretion has been taken into account for its desaturation effect of the corotation torque. If more rapid gas flows (e.g., wind-driven accretion) are considered, the desaturation effect can be modified. In MRI-inactive disks, in which the wind-driven accretion dominates the disk evolution, the gas flow at the midplane plays an important role. If this flow is fast, the corotation torque is efficiently desaturated. Then, the fact that the surface density slope can be positive in the inner region due to the wind torque can generate an outward migration region extended to super-Earth mass planets. In this case, we observe that no planets fall onto the central star in N-body simulations with migration forces imposed to reproduce such migration pattern. We also see that super-Earth mass planets can undergo outward migration. Relatively rapid gas flows affects type I migration and thus the formation of close-in planets.

"Hot super-Earths" (or "Mini-Neptunes") between 1 and 4 times Earth's size with period shorter than 100 days orbit 30-50\% of Sun-like type stars. Their orbital configuration -- measured as the period ratio distribution of adjacent planets in multi-planet systems -- is a strong constraint for formation models. Here we use N-body simulations with synthetic forces from an underlying evolving gaseous disk to model the formation and long-term dynamical evolution of super-Earth systems. While the gas disk is present, planetary embryos grow and migrate inward to form a resonant chain anchored at the inner edge of the disk. These resonant chains are far more compact than the observed super-Earth systems. Once the gas dissipates resonant chains may become dynamically unstable. They undergo a phase of giant impacts that spreads the systems out. Disk turbulence has no measurable effect on the outcome. Our simulations match observations if a small fraction of resonant chains remain stable, while most super-Earths undergo a late dynamical instability. Our statistical analysis restricts the contribution of stable systems to less than $25\%$. Our results also suggest that the large fraction of observed single planet systems does not necessarily imply any dichotomy in the architecture of planetary systems. Finally, we use the low abundance of resonances in Kepler data to argue that, in reality, the survival of resonant chains happens likely only in $\sim 5\%$ of the cases. This leads to a mystery: in our simulations only 50-60\% of resonant chains became unstable whereas at least 75\% (and probably 90-95\%) must be unstable to match observations.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We have investigated the final accretion stage of terrestrial planets from Mars-mass protoplanets that formed through oligarchic growth in a disk comparable to the minimum mass solar nebula (MMSN), through N-body simulation including random torques exerted by disk turbulence due to Magneto-Rotational-Instability. For the torques, we used the semi-analytical formula developed by Laughlin et al.(2004). The damping of orbital eccentricities (in all runs) and type-I migration (in some runs) due to the tidal interactions with disk gas are also included. We found that the orbital eccentricities pumped up by the turbulent torques and associated random walks in semimajor axes tend to delay isolation of planets, resulting in more coagulation of planets than in the case without turbulence. The eccentricities are still damped after planets become isolated. As a result, the number of final planets decreases with increase in strength of the turbulence, while Earth-mass planets with small eccentricities are still formed. In the case of relatively strong turbulence, the number of final planets are 4-5 at 0.5-2AU, which is consistent with Solar system, for relatively wide range of disk surface density (~10^{-4}-10^{-2} times MMSN).