In models of planetary accretion, pebbles form by dust coagulation and rapidly migrate toward the central star. Planetesimals may continuously form from pebbles over the age of the protoplanetary disk by yet uncertain mechanisms. Meanwhile, large planetary embryos grow efficiently by accumulation of leftover pebbles that are not incorporated in planetesimals. Although this process, called pebble accretion, is offering a new promising pathway for formation of giant planets' cores, architectures of planetary systems formed through the process remain elusive. In the present paper, we perform simulations of formation of planetary systems using a particle-based hybrid code, to which we implement most of the key physical effects as precisely as possible. We vary the size of a protoplanetary disk, the turbulent viscosity, the pebble size, the planetesimal formation efficiency, and the initial mass distribution of planetesimals. Our simulations show that planetesimals first grow by mutual collisions if their initial size is the order of $100$ km or less. Once planetesimals reach $\sim$ 1000 km in size, they efficiently grow by pebble accretion. If pebble supply from the outer region continues for a long period of time in a large protoplanetary disk, planetary embryos become massive enough to commence runaway gas accretion, resulting in gas giant planets. Our simulations suggest that planetary systems like ours form from protoplanetary disks with moderately high turbulent viscosities. If the disk turbulent viscosity is low enough, a planet opens up a gap in the gas disk and halts accretion of pebbles even before the onset of runway gas accretion. Such a disk produces a planetary system with several Neptune-size planets.

We study the B ring's complex optical depth structure. The source of this structure may be the complex dynamics of the Keplerian shear and the self-gravity of the ring particles. The outcome of these dynamic effects depends sensitively on the collisional and physical properties of the particles. Two mechanisms can emerge that dominate the macroscopic physical structure of the ring: self-gravity wakes and viscous overstability. Here we study the interplay between these two mechanisms by using our recently developed particle collision method that allows us to better model the inter-particle contact physics. We find that for a constant ring surface density and particle internal density, particles with rough surfaces tend to produce axisymmetric ring features associated with the viscous overstability, while particles with smoother surfaces produce self-gravity wakes.

Aims. To investigate the mid-infrared (MIR) characteristics of Saturn's rings. Methods. We collected and analyzed MIR high spatial resolution images of Saturn's rings obtained in January 2008 and April 2005 with COMICS mounted on Subaru Telescope, and investigated the spatial variation in the surface brightness of the rings in multiple bands in the MIR. We also composed the spectral energy distributions (SEDs) of the C, B, and A rings and the Cassini Division, and estimated the temperatures of the rings from the SEDs assuming the optical depths. Results. We find that the C ring and the Cassini Division were warmer than the B and A rings in 2008, which could be accounted for by their lower albedos, lower optical depths, and smaller self-shadowing effect. We also find that the C ring and the Cassini Division were considerably brighter than the B and A rings in the MIR in 2008 and the radial contrast of the ring brightness is the inverse of that in 2005, which is interpreted as a result of a seasonal effect with changing elevations of the sun and observer above the ring plane.

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

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

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

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

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

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

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

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

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

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

Remnant planetesimals might have played an important role in reducing the orbital eccentricities of the terrestrial planets after their formation via giant impacts. However, the population and the size distribution of remnant planetesimals during and after the giant impact stage are unknown, because simulations of planetary accretion in the runaway growth and giant impact stages have been conducted independently. Here we report results of direct N-body simulations of the formation of terrestrial planets beginning with a compact planetesimal disk. The initial planetesimal disk has a total mass and angular momentum as observed for the terrestrial planets, and we vary the width (0.3 and 0.5AU) and the number of planetesimals (1000-5000). This initial configuration generally gives rise to three final planets of similar size, and sometimes a fourth small planet forms near the location of Mars. Since a sufficient number of planetesimals remains, even after the giant impact phase, the final orbital eccentricities are as small as those of the Earth and Venus.

We have undertaken a thorough dynamical investigation of five extrasolar planetary systems using extensive numerical experiments. The systems Gl 777 A, HD 72659, Gl 614, 47 Uma and HD 4208 were examined concerning the question of whether they could host terrestrial like planets in their habitable zones (=HZ). First we investigated the mean motion resonances between fictitious terrestrial planets and the existing gas giants in these five extrasolar systems. Then a fine grid of initial conditions for a potential terrestrial planet within the HZ was chosen for each system, from which the stability of orbits was then assessed by direct integrations over a time interval of 1 million years. The computations were carried out using a Lie-series integration method with an adaptive step size control. This integration method achieves machine precision accuracy in a highly efficient and robust way, requiring no special adjustments when the orbits have large eccentricities. The stability of orbits was examined with a determination of the Renyi entropy, estimated from recurrence plots, and with a more straight forward method based on the maximum eccentricity achieved by the planet over the 1 million year integration. Additionally, the eccentricity is an indication of the habitability of a terrestrial planet in the HZ; any value of e>0.2 produces a significant temperature difference on a planet's surface between apoapse and periapse. The results for possible stable orbits for terrestrial planets in habitable zones for the five systems are summarized as follows: for Gl 777 A nearly the entire HZ is stable, for 47 Uma, HD 72659 and HD 4208 terrestrial planets can survive for a sufficiently long time, while for Gl 614 our results exclude terrestrial planets moving in stable orbits within the HZ.