During planet formation gravitational interaction between a planetary embryo and the protoplanetary gas disc causes orbital migration of the planetary embryo, which plays an important role in shaping the final planetary system. While migration sometimes occurs in the supersonic regime, wherein the relative velocity between the planetary embryo and the gas is higher than the sound speed, migration prescriptions proposed thus far describing the planet-disc interaction force and the timescales of orbital change in the supersonic regime are inconsistent with one another. Here we discuss the details of existing prescriptions in the literature and derive a new simple and intuitive formulation for planet-disc interactions based on dynamical friction that can be applied in both supersonic and subsonic cases. While the existing prescriptions assume particular disc models, ours include the explicit dependence on the disc parameters; hence it can be applied to discs with any radial surface density and temperature dependence (except for the local variations with radial scales less than the disc scale height). Our prescription will reduce the uncertainty originating from different literature formulations of planet migration and will be an important tool to study planet accretion processes, especially when studying the formation of close-in low-mass planets that are commonly found in exoplanetary systems.

The ice-giant planet Uranus likely underwent a giant impact, given that its spin axis is tilted by 98 degrees. That its satellite system is equally inclined and prograde suggests that it was formed as a consequence of the impact. However, the disks predicted by the impact simulations generally have sizes one order smaller and masses two orders larger than those of the observed system at present. Here we show, by means of a theoretical model, that the Uranian satellite formation is regulated by the evolution of the impact-generated disk. Because the vaporization temperature of water ice is low and both Uranus and the impactor are assumed to be ice-dominated, we can conclude that the impact-generated disk has mostly vaporized. We predict that the disk lost a significant amount of water vapour mass and spread to the levels of the current system until the disk cooled down enough for ice condensation and accretion of icy particles to begin. From the predicted distribution of condensed ices, our N-body simulation is able to reproduce the observed mass-orbit configuration of Uranian satellites. This scenario contrasts with the giant-impact model for the Earth's Moon, in which about half of the compact, impact-generated, solid or liquid disk is immediately incorporated into the Moon on impact.

The detection of Earth-size exoplanets around low-mass stars -- in stars such as Proxima Centauri and TRAPPIST-1 -- provide an exceptional chance to improve our understanding of the formation of planets around M stars and brown dwarfs. We explore the formation of such planets with a population synthesis code based on a planetesimal-driven model previously used to study the formation of the Jovian satellites. Because the discs have low mass and the stars are cool, the formation is an inefficient process that happens at short periods, generating compact planetary systems. Planets can be trapped in resonances and we follow the evolution of the planets after the gas has dissipated and they undergo orbit crossings and possible mergers. We find that formation of planets above Mars mass and in the planetesimal accretion scenario, is only possible around stars with masses $M_{\star} \ge 0.07 M_{sun}$ and discs of $M_{disc} \ge 10^{-2}~M_{sun} $. We find that planets above Earth-mass form around stars with masses larger than $0.15 M_{\oplus}$ while planets larger than $5 M_{\oplus}$do not form in our model, even not under the most optimal conditions (massive disc), showing that planets such as GJ 3512b form with another, more efficient mechanism. Our results show that the majority of planets form with a significant water fraction; that most of our synthetic planetary systems have 1, 2, or 3 planets, but those with 4, 5, 6, and 7 planets are also common, confirming that compact planetary systems with many planets should be a relatively common outcome of planet formation around small stars.

It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter's formation. However, understanding the specifics of the formation process is challenging as both small particles (pebbles) as well as the satellites are subject to fast migration processes. Here, we hypothesize a new scenario for the origin of the Galilean system, based on the capture of several planetesimal seeds and subsequent slow accretion of pebbles. To halt migration, we invoke an inner disk truncation radius, and other parameters are tuned for the model to match physical, dynamical, compositional, and structural constraints. In our scenario it is natural that Ganymede's mass is determined by pebble isolation. Our slow-pebble-accretion scenario then reproduces the following characteristics: (1) the mass of all the Galilean satellites; (2) the orbits of Io, Europa, and Ganymede captured in mutual 2:1 mean motion resonances; (3) the ice mass fractions of all the Galilean satellites; (4) the unique ice-rock partially differentiated Callisto and the complete differentiation of the other satellites. Our scenario is unique to simultaneously reproduce these disparate properties.

The Earth's ocean mass is only 2.3 x 10^{-4} of the whole planet mass. Even including water in the interior, it would be at most 10^{-3}-10^{-2}. Ancient Mars may have had a similar or slightly smaller water fraction. It is important to clarify the water delivery mechanism to rocky planets in habitable zones in exoplanetary systems, as well as that to the Earth and Mars. Here, we consider water delivery to planets by icy pebbles after the snowline inwardly passes the planetary orbits and derive the water mass fraction (f_{water}) of the final planet as a function of disk parameters and discuss the parameters that reproduce f_{water} comparable to that inferred for the Earth and ancient Mars. We calculate the growth of icy pebbles and their radial drift with a 1D model, and accretion of icy pebbles onto planets, by simultaneously solving the snowline migration and the disk dissipation, to evaluate f_{water} of the planets. We find that f_{water} is regulated by the total mass (M_{res}) of icy dust materials preserved in the outer disk regions at the timing (t = t_{snow}) of the snowline passage of the planetary orbit. Because M_{res} decays rapidly after the pebble formation front reaches the disk outer edge (at t = t_{pff}), f_{water} is sensitive to the ratio t_{snow}/t_{pff}, which is determined by the disk parameters. We find t_{snow}/t_{pff} < 10 or > 10 is important. Deriving an analytical formula for f_{water} that reproduces the numerical results, we find that f_{water} of a rocky planet near 1 au is ~ 10^{-4}-10^{-2}, in the disks with initial disk size ~ 30-50 au and the initial disk mass accretion rate ~ (10^{-8}-10^{-7}) M_sun/y. Because these disks may be median or slightly compact/massive disks, the water fraction of rocky planets in habitable zones may be often similar to that of the Earth, if the icy pebble accretion is responsible for the water delivery.

We report the detection of an excess in dust continuum emission at 233~GHz (1.3~mm in wavelength) in the protoplanetary disk around TW~Hya revealed through high-sensitivity observations at $\sim$3~au resolution with the Atacama Large Millimeter/submillimeter Array (ALMA). The sensitivity of the 233~GHz image has been improved by a factor of 3 with regard to that of our previous cycle 3 observations. The overall structure is mostly axisymmetric, and there are apparent gaps at 25 and 41 au as previously reported. The most remarkable new finding is a few au-scale excess emission in the south-west part of the protoplanetary disk. The excess emission is located at 52 au from the disk center and is 1.5 times brighter than the surrounding protoplanetary disk at a significance of 12$\sigma$. We performed a visibility fitting to the extracted emission after subtracting the axisymmetric protoplanetary disk emission and found that the inferred size and the total flux density of the excess emission are 4.4$\times$1.0~au and 250~$\mu$Jy, respectively. The dust mass of the excess emission corresponds to 0.03~$M_\oplus$ if a dust temperature of 18~K is assumed. Since the excess emission can also be marginally identified in the Band 7 image at almost the same position, the feature is unlikely to be a background source. The excess emission can be explained by a dust clump accumulated in a small elongated vortex or a massive circumplanetary disk around a Neptune mass forming-planet.

Most major planetary bodies in the solar system rotate in the same direction as their orbital motion: their spin is prograde. Theoretical studies to explain the direction as well as the magnitude of the spin vector have had mixed success. When the accreting building blocks are $\sim$ km-size planetesimals -- as predicted by the classical model -- the accretion process is so symmetric that it cancels out prograde with retrograde spin contributions, rendering the net spin minute. For this reason, the currently-favored model for the origin of planetary rotation is the giant impact model, in which a single collision suffices to deliver a spin, which magnitude is close to the breakup rotation rate. However, the giant impact model does not naturally explain the preference for prograde spin. Similarly, an increasing number of spin-vector measurement of asteroids also shows that the spin vector of large (primordial) asteroids is not isotropic. Here, we re-assess the viability of smaller particles to bestow planetary bodies with a net spin, focusing on the pebble accretion model in which gas drag and gravity join forces to accrete small particles at a large cross section. Similar to the classical calculation for planetesimals, we integrate the pebble equation of motion and measure the angular momentum transfer at impact. We consider a variety of disk conditions and pebble properties and conduct our calculations in the limits of 2D (planar) and 3D (homogeneous) pebble distributions. We find that in certain regions of the parameter space the angular momentum transfer is significant, much larger than with planetesimals and on par with or exceeding the current spin of planetary bodies.

It is important to clarify where and when rocky and icy planetesimals are formed in a viscously evolving disk. We wish to understand how local runaway pile-up of solids occurs inside or outside the snow line. We assume an icy pebble contains micron-sized silicate grains that are uniformly mixed with ice and are released during the ice sublimation. Using a local one-dimensional code, we solve the radial drift and the turbulent diffusion of solids and the water vapor, taking account of their sublimation/condensation around the snow line. We systematically investigate effects of back-reactions of the solids to gas on the radial drift and diffusion of solids, scale height evolution of the released silicate particles, and possible difference in effective viscous parameters between that for turbulent diffusion ($\alpha_{\rm tur}$) and that for the gas accretion rate onto the central star ($\alpha_{\rm acc}$). We study the dependence on the ratio of the solid mass flux to the gas ($F_{\rm p/g}$). We show that the favorable locations for the pile-up of silicate grains and icy pebbles are the regions in the proximity of the water snow line inside and outside it, respectively. We found that runaway pile-ups occur when both the back-reactions for radial drift and diffusion are included. In the case with only the back-reaction for the radial drift, no runaway pile-up is found except for extremely high pebble flux, while the condition of streaming instability can be satisfied for relatively large $F_{\rm p/g}$ as found in the past literatures. If the back-reactions for radial diffusion is considered, the runaway pile-up occurs for reasonable value of pebble flux. The runaway pile-up of silicate grains that would lead to formation of rocky planetesimals occurs for $\alpha_{\rm tur} \ll \alpha_{\rm acc}$, while the runaway pile-up of icy pebbles is favored for $\alpha_{\rm tur} \sim \alpha_{\rm acc}$.

Planets in the "Habitable Zones" around M-type stars are important targets for characterization in future observations. Due to tidal-locking in synchronous spin-orbit rotations, the planets tend to have a hot dayside and a cold nightside. On the cold nightside, water vapor transferred from the dayside can be frozen in ("cold trap") or the major atmospheric constituent could also condense ("atmospheric collapse") if the atmosphere is so thin that the heat re-distribution is not efficient, in the case of a single M-type star. Motivated by the abundance of binary star systems, we investigate the effects of irradiation from a G-type companion star on the climate of a tidally locked planet around an M-type star using the 2D energy balance model. We find that the irradiation from the G-type star is more effective at warming up the nightside of the planet than the dayside. This contributes to the prevention of the irreversible trapping of water and atmosphere on the cold nightside, broadening the parameter space where tidally locked planets can maintain surface liquid water. Tidally locked ocean planets with < ~ 0.3 bar atmospheres or land planets with < ~ 3 bar atmospheres can realize temperate climate with surface liquid water when they are also irradiated by a companion star with a separation of 1 - 4 au. We also demonstrate that planets with given properties can be in the Earth-like temperate climate regime or in a completely frozen state under the same total irradiation.

Understanding the fate of planetary systems through white dwarfs which accrete debris crucially relies on tracing the orbital and physical properties of exo-asteroids during the giant branch phase of stellar evolution. Giant branch luminosities exceed the Sun's by over three orders of magnitude, leading to significantly enhanced Yarkovsky and YORP effects on minor planets. Here, we place bounds on Yarkovsky-induced differential migration between asteroids and planets during giant branch mass loss by modelling one exo-Neptune with inner and outer exo-Kuiper belts. In our bounding models, the asteroids move too quickly past the planet to be diverted from their eventual fate, which can range from: (i) populating the outer regions of systems out to 10^4-10^5 au, (ii) being engulfed within the host star, or (iii) experiencing Yarkovsky-induced orbital inclination flipping without any Yarkovsky-induced semimajor axis drift. In these violent limiting cases, temporary resonant trapping of asteroids with radii of under about 10 km by the planet is insignificant, and capture within the planet's Hill sphere requires fine-tuned dissipation. The wide variety of outcomes presented here demonstrates the need to employ sophisticated structure and radiative exo-asteroid models in future studies. Determining where metal-polluting asteroids reside around a white dwarf depends on understanding extreme Yarkovsky physics.

The three-dimensional structure of the gas flow around a planet is thought to influence the accretion of both gas and solid materials. In particular, the outflow in the mid-plane region may prevent the accretion of the solid materials and delay the formation of super-Earths' cores. However, it is not yet understood how the nature of the flow field and outflow speed change as a function of the planetary mass. In this study, we investigate the dependence of gas flow around a planet embedded in a protoplanetary disc on the planetary mass. Assuming an isothermal, inviscid gas disc, we perform three-dimensional hydrodynamical simulations on the spherical polar grid, which has a planet located at its centre. We find that gas enters the Bondi or Hill sphere at high latitudes and exits through the mid-plane region of the disc regardless of the assumed dimensionless planetary mass $m=R_{\rm Bondi}/H$, where $R_{\rm Bondi}$ and $H$ are the Bondi radius of the planet and disc scale height, respectively. The altitude from where gas predominantly enters the envelope varies with the planetary mass. The outflow speed can be expressed as $|u_{\rm out}|=\sqrt{3/2}mc_{\rm s}$ $(R_{\rm Bondi}\leq R_{\rm Hill})$ or $|u_{\rm out}|=\sqrt{3/2}(m/3)^{1/3} c_{\rm s}$ ($R_{\rm Bondi}\geq R_{\rm Hill}$), where $c_{\rm s}$ is the isothermal sound speed and $R_{\rm Hill}$ is the Hill radius. The outflow around a planet may reduce the accretion of dust and pebbles onto the planet when $m\gtrsim\sqrt{\rm St}$, where St is the Stokes number. Our results suggest that the flow around proto-cores of super-Earths may delay their growth and, consequently, help them to avoid runaway gas accretion within the lifetime of the gas disc.

Planetary migration is a major challenge for planet formation theories. The speed of Type I migration is proportional to the mass of a protoplanet, while the final decade of growth of a pebble-accreting planetary core takes place at a rate that scales with the mass to the two-thirds power. This results in planetary growth tracks (i.e., the evolution of a protoplanet's mass versus its distance from the star) that become increasingly horizontal (migration-dominated) with rising mass of the protoplanet. It has been shown recently that the migration torque on a protoplanet is reduced proportional to the relative height of the gas gap carved by the growing planet. Here we show from 1-D simulations of planet-disc interaction that the mass at which a planet carves a 50% gap is approximately 2.3 times the pebble isolation mass. Our measurements of the pebble isolation mass from 1-D simulations match published 3-D results relatively well, except at very low viscosities where the 3-D pebble isolation mass is significantly higher, possibly due to gap edge instabilities not captured in 1-D. The pebble isolation mass demarks the transition from pebble accretion to gas accretion. Gas accretion to form gas-giant planets therefore takes place over a few astronomical units of migration after reaching first the pebble isolation mass and, shortly after, the 50% gap mass. Our results demonstrate how planetary growth can outperform migration, both during core accretion and during gas accretion, even when the Stokes number of the pebbles is small, St~0.01, and the pebble-to-gas flux ratio in the protoplanetary disc is in the nominal range of 0.01-0.02. We find that planetary growth is very rapid in the first million years of the protoplanetary disc and that the probability for forming gas-giant planets increases with the initial size of the protoplanetary disc and with decreasing turbulent diffusion.

We compare the planet-to-star mass-ratio distribution measured by gravitational microlensing to core accretion theory predictions from population synthesis models. The core accretion theory's runaway gas accretion process predicts a dearth of intermediate-mass giant planets that is not seen in the microlensing results. In particular, the models predict $\sim10\,\times$ fewer planets at mass ratios of $10^{-4} \leq q \leq 4 \times 10^{-4}$ than inferred from microlensing observations. This tension implies that gas giant formation may involve processes that have hitherto been overlooked by existing core accretion models or that the planet-forming environment varies considerably as a function of host-star mass. Variation from the usual assumptions for the protoplanetary disk viscosity and thickness could reduce this discrepancy, but such changes might conflict with microlensing results at larger or smaller mass ratios, or with other observations. The resolution of this discrepancy may have important implications for planetary habitability because it has been suggested that the runaway gas accretion process may have triggered the delivery of water to our inner solar system. So, an understanding of giant planet formation may help us to determine the occurrence rate of habitable planets.

The mass and semimajor axis distribution of gas giants in exoplanetary systems obtained by radial velocity surveys shows that super-jupiter-mass planets are piled up at > 1 au, while jupiter/sub-jupiter-mass planets are broadly distributed from ~0.03 au to beyond 1 au. This feature has not been explained by theoretical predictions. In order to reconcile this inconsistency, we investigate evolution of gas giants with a new type II migration formula by Kanagawa et al. (2018), by comparing the migration, growth timescales of gas giants, and disk lifetime and by population synthesis simulation. While the classical migration model assumes that a gas giant opens up a clear gap in the protoplanetary disk and the planet migration tied to the disk gas accretion, recent high-resolution simulations show that the migration of gap-opening planets is decoupled from the disk gas accretion and Kanagawa et al. (2018) proposed that type II migration speed is no other than type I migration speed with the reduced disk gas surface density in the gap. We show that with this new formula, type II migration is significantly reduced for super-jupiter-mass planets, if the disk accretion is driven by the disk wind as suggested by recent MHD simulations. Population synthesis simulations show that super-jupiter-mass planets remain at > 1 au without any additional ingredient such as disk photoevaporation. Therefore, the mystery of the pile-up of gas giants at > 1 au would be theoretically solved, if the new formula is confirmed and wind-driven disk accretion dominates.

The formation and orbital evolution of Saturn's inner mid-sized moons are still debated. The most puzzling aspects are 1) how the Tethys-Dione pair and the Mimas-Enceladus pair passed through their strong 3:2 mean-motion resonances during the tidal orbital evolution, and 2) the current strong heat flow from Enceladus, which is a few orders of magnitude higher than the tidal energy dissipation caused by the present orbital eccentricity of Enceladus. Here we perform N-body simulations of the moons' orbital evolution from various initial conditions -- assuming that the moons were formed from Saturn's hypothetical massive ring -- and investigate possible paths to solve the above difficulties. If the influence of the rings is neglected, we find that the Tethys-Dione pair cannot avoid becoming trapped in the mean-motion resonances as they recede from Saturn, and that the Tethys-Enceladus pair cannot avoid collisions after the resonance trapping, in case Saturn's quality factor is smaller than 15,000. These findings are inconsistent with the current orbital configuration. However, taking into account both the eccentricity excitation and the orbital expansion caused by the ring torque, we find that these resonance captures are avoided. With the high eccentricity pumped up by the torque, Enceladus passes through the resonances with Tethys, and the Dione-Tethys pair passes through their 2:1 and possibly the 3:2 resonance as well. After Enceladus resides beyond the 2:1 resonance with the outer ring edge, the eccentricity can be tidally damped. While this is a promising path of evolution, in most runs, Enceladus collides with Tethys by the excited eccentricity. There is a hint that a ring mass decrease could avoid the collision. The tidal heat due to the eccentricity excitation by the ring torque is stored in the moons and slowly radiated away. It may account for the current high heat flow.

We want to investigate how planet formation is imprinted on stellar surface composition using up-to-date stellar evolution models. We simulate the evolution of pre-main-sequence stars as a function of the efficiency of heat injection during accretion, the deuterium mass fraction, and the stellar mass. For simplicity, we assume that planet formation leads to the late accretion of zero-metallicity gas, diluting the surface stellar composition as a function of the mass of the stellar outer convective zone. We adopt $150\,{\mathrm{M}_\oplus}(M_\star/\mathrm{M}_\odot)(Z/\mathrm{Z}_\odot)$ as an uncertain but plausible estimate of the mass of heavy elements that is not accreted by stars with giant planets, including our Sun. By combining our stellar evolution models to these estimates, we evaluate the consequences of planet formation on stellar surface composition. We show that after the first $\sim0.1$ Myr, the evolution of the convective zone follows classical evolutionary tracks within a factor of two in age. We find that planet formation should lead to a scatter in stellar surface composition that is larger for high-mass stars than for low-mass stars. We predict a spread in [Fe/H] of approximately $0.02$ dex for stars with $T_\mathrm{eff}\sim 5500\,$K, marginally compatible with differences in metallicities observed in some binary stars with planets. Stars with $T_\mathrm{eff}\geq 7000\,$K may show much larger [Fe/H] deficits, by 0.6 dex or more, compatible with the existence of refractory-poor $\lambda$ Boo stars. We also find that planet formation may explain the lack of refractory elements seen in the Sun as compared to solar twins, but only if the ice-to-rock ratio in the solar-system planets is less than $\approx0.4$ and planet formation began less than $\approx1.3$ Myr after the beginning of the formation of the Sun. (abbreviated)

The formation of multiple close-in low-mass exoplanets is still a mystery. The challenge is to build a system wherein the outermost planet is beyond 0.2 AU from the star. Here we investigate how the prescription for type I planet migration affects the ability to trap multiple planets in a resonant chain near the inner edge of the protostellar disc. A sharp edge modelled as a hyperbolic tangent function coupled with supersonic corrections to the classical type I migration torques results in the innermost planets being pushed inside the cavity through resonant interaction with farther planets because migration is starward at slightly supersonic eccentricities. Planets below a few Earth masses are generally trapped in a resonant chain with the outermost planet near the disc edge, but long-term stability is not guaranteed. For more massive planets the migration is so fast that the eccentricity of the innermost resonant pair is excited to highly supersonic levels due to decreased damping on the innermost planet as it is pushed inside the cavity; collisions frequently occur and the system consists one or two intermediate-mass planets residing closer to the star than the disc's inner edge. We found a neat pileup of resonant planets outside the disc edge only if the corotation torque does not rapidly diminish at high eccentricity. We call for detailed studies on planet migration near the disc's inner edge, which is still uncertain, and for an improved understanding of eccentricity damping and disc torques in the supersonic regime.

Dynamical models of planet formation coupled with cosmochemical data from martian meteorites show that Mars' isotopic composition is distinct from that of Earth. Reconciliation of formation models with meteorite data require that Mars grew further from the Sun than its present position. Here, we evaluate this compositional difference in more detail by comparing output from two $N$-body planet formation models. The first of these planet formation models simulates what is termed the "Classical" case wherein Jupiter and Saturn are kept in their current orbits. We compare these results with another model based on the "Grand Tack", in which Jupiter and Saturn migrate through the primordial asteroid belt. Our estimate of the average fraction of chondrite assembled into Earth and Mars assumes that the initial solid disk consists of only sources of enstatite chondrite composition in the inner region, and ordinary chondrite in the outer region. Results of these analyses show that both models tend to yield Earth and Mars analogues whose accretion zones overlap. The Classical case fares better in forming Mars with its documented composition (29% to 68% enstatite chondrite plus 32% to 67% ordinary chondrite) though the Mars analogues are generally too massive. However, if we include the restriction of mass on the Mars analogues, the Classical model does not work better. We also further calculate the isotopic composition of $^{17} \rm O$, $^{50} \rm Ti$, $^{54} \rm Cr$, $^{142} \rm Nd$, $^{64} \rm Ni$, and $^{92} \rm Mo$ in the martian mantle from the Grand Tack simulations. We find that it is possible to match the calculated isotopic composition of all the above elements in Mars' mantle with their measured values, but the resulting uncertainties are too large to place good restriction on the early dynamical evolution and birth place of Mars.

We report the detection of planets around two evolved giant stars from radial velocity measurements at Okayama Astrophysical observatory. 24 Boo (G3IV) has a mass of $0.99\,M_{\odot}$, a radius of $10.64\,R_{\odot}$, and a metallicity of ${\rm [Fe/H]}=-0.77$. The star hosts one planet with a minimum mass of $0.91\,M_{\rm Jup}$ and an orbital period of $30.35{\rm d}$. The planet has one of the shortest orbital periods among those ever found around evolved stars by radial-veloocity methods. The stellar radial velocities show additional periodicity with $150{\rm d}$, which are probably attributed to stellar activity. The star is one of the lowest-metallicity stars orbited by planets currently known. $\gamma$ Lib (K0III) is also a metal-poor giant with a mass of $1.47\,M_{\odot}$, a radius of $11.1\,R_{\odot}$, and ${\rm [Fe/H]}=-0.30$. The star hosts two planets with minimum masses of $1.02M_{\rm Jup}$ and $4.58\,M_{\rm Jup}$, and periods of $415{\rm d}$ and $964{\rm d}$, respectively. The star has the second lowest metallicity among the giant stars hosting more than two planets. Dynamical stability analysis for the $\gamma$ Lib system sets a minimum orbital inclination angle to be about $70^{\circ}$ and suggests that the planets are in 7:3 mean-motion resonance, though the current best-fitted orbits to the radial-velocity data are not totally regular.

We have investigated the formation of a circumstellar wide-orbit gas giant planet in a multiple stellar system. We consider a model of orbital circularization for the core of a giant planet after it is scattered from an inner disk region by a more massive planet, which was proposed by Kikuchi et al (2014). We extend their model for single star systems to binary (multiple) star systems, by taking into account tidal truncation of the protoplanetary gas disk by a binary companion. As an example, we consider a wide-orbit gas giant in a hierarchical triple system, HD131399Ab. The best-fit orbit of the planet is that with semimajor axis $\sim 80$ au and eccentricity $\sim 0.35$. As the binary separation is $\sim 350$ au, it is very close to the stability limit, which is puzzling. With the original core location $\sim 20$-30 au, the core (planet) mass $\sim 50 M_{\rm E}$ and the disk truncation radius $\sim 150$ au, our model reproduces the best-fit orbit of HD131399Ab. We find that the orbit after the circularization is usually close to the stability limit against the perturbations from the binary companion, because the scattered core accretes gas from the truncated disk. Our conclusion can also be applied to wider or more compact binary systems if the separation is not too large and another planet with $> \sim$ 20-30 Earth masses that scattered the core existed in inner region of the system.

The icy satellites around Jupiter are considered to have formed in a circumplanetary disk. While previous models focused on the formation of satellites starting from satellitesimals, the question of how satellitesimals form from smaller dust particles has not been addressed so far. In this work, we study the possibility that satellitesimals form in situ in a circumplanetary disk. We calculate the radial distribution of the surface density and representative size of icy dust particles that grow by colliding with each other and drift toward the central planet in a steady circumplanetary disk with a continuous supply of gas and dust from the parent protoplanetary disk. The radial drift barrier is overcome if the ratio of the dust to gas accretion rates onto the circumplanetary disk, $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$, is high and the strength of turbulence, $\alpha$, is not too low. The collision velocity is lower than the critical velocity of fragmentation when $\alpha$ is low. Taken together, we find that the conditions for satellitesimal formation via dust coagulation are given by $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}\ge1$ and $10^{-4}\le\alpha<10^{-2}$. The former condition is generally difficult to achieve, suggesting that the in-situ satellitesimal formation via particle sticking is viable only under an extreme condition. We also show that neither satellitesimal formation via the collisional growth of porous aggregates nor via streaming instability is viable as long as $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$ is low.

Context. Planet formation with pebbles has been proposed to solve a couple of long-standing issues in the classical formation model. Some sophisticated simulations have been done to confirm the efficiency of pebble accretion. However, there has not been any global N-body simulations that compare the outcomes of planet formation via pebble accretion with observed extrasolar planetary systems. Aims. In this paper, we study the effects of a range of initial parameters of planet formation via pebble accretion, and present the first results of our simulations. Methods. We incorporate the pebble accretion model by Ida et al. (2016) in the N-body code SyMBA (Duncan et al. 1998), along with the effects of gas accretion, eccentricity and inclination damping and planet migration in the disc. Results. We confirm that pebble accretion leads to a variety of planetary systems, but have difficulty in reproducing observed properties of exoplanetary systems, such as planetary mass, semimajor axis, and eccentricity distributions. The main reason behind this is a too-efficient type I migration, which sensitively depends on the disc model. However, our simulations also lead to a few interesting predictions. First, we find that formation efficiencies of planets depend on the stellar metallicities, not only for giant planets, but also for Earths (Es) and Super-Earths (SEs). The dependency for Es/SEs is subtle. Although higher metallicity environments lead to faster formation of a larger number of Es/SEs, they also tend to be lost later via dynamical instability. Second, our results indicate that a wide range of bulk densities observed for Es and SEs is a natural consequence of dynamical evolution of planetary systems. Third, the ejection trend of our simulations suggest that one free-floating E/SE may be expected for two smaller-mass planets.

With approximately one ninth of Earth's mass, Mars is widely considered to be a stranded planetary embryo that never became a fully-grown planet. A currently popular planet formation theory predicts that Mars formed near Earth and Venus and was subsequently scattered outwards to its present location. In such a scenario, the compositions of the three planets are expected to be similar to each other. However, bulk elemental and isotopic data for martian meteorites demonstrate that key aspects of Mars' composition are markedly different from that of Earth. This suggests that Mars formed outside of the terrestrial feeding zone during primary accretion. It is therefore probable that Mars always remained significantly farther from the Sun than Earth; its growth was stunted early and its mass remained relatively low. Here we identify a potential dynamical pathway that forms Mars in the asteroid belt and keeps it outside of Earth's accretion zone while at the same time accounting for strict age and compositional constraints, as well as mass differences. Our uncommon pathway (approximately 2% probability) is based on the Grand Tack scenario of terrestrial planet formation, in which the radial migration by Jupiter gravitationally sculpts the planetesimal disc at Mars' current location. We conclude that Mars' formation requires a specific dynamical pathway, while this is less valid for Earth and Venus. We further predict that} Mars' volatile budget is most likely different from Earth's and that Venus formed close enough to our planet that it is expected to have a nearly identical composition from common building blocks.

Recent theoretical work has shown that the pre-main-sequence (PMS) evolution of stars is much more complex than previously envisioned. Instead of the traditional steady, one-dimensional solution, accretion may be episodic and not necessarily symmetrical, thereby affecting the energy deposited inside the star and its interior structure. Given this new framework, we want to understand what controls the evolution of accreting stars. We use the MESA stellar evolution code with various sets of conditions. In particular, we account for the (unknown) efficiency of accretion in burying gravitational energy into the protostar through a parameter, $\xi$, and we vary the amount of deuterium present. We confirm the findings of previous works that the evolution changes significantly with the amount of energy that is lost during accretion. We find that deuterium burning also regulates the PMS evolution. In the low-entropy accretion scenario, the evolutionary tracks in the H-R diagram are significantly different from the classical tracks and are sensitive to the deuterium content. A comparison of theoretical evolutionary tracks and observations allows us to exclude some cold accretion models ($\xi\sim 0$) with low deuterium abundances. We confirm that the luminosity spread seen in clusters can be explained by models with a somewhat inefficient injection of accretion heat. The resulting evolutionary tracks then become sensitive to the accretion heat efficiency, initial core entropy, and deuterium content. In this context, we predict that clusters with a higher D/H ratio should have less scatter in luminosity than clusters with a smaller D/H. Future work on this issue should include radiation-hydrodynamic simulations to determine the efficiency of accretion heating and further observations to investigate the deuterium content in star-forming regions. (abbrev.)

We have investigated the probability of temporary capture of asteroids in eccentric orbits by a planet in a circular or an eccentric orbit through analytical and numerical calculations. We found that in the limit of the circular orbit, the capture probability is $\sim 0.1\%$ of encounters to the planet's Hill sphere, independent of planetary mass and semimajor axis. In general, the temporary capture becomes more difficult as the planet's eccentricity ($e_{\rm p}$) increases. We found that the capture probability is almost independent of $e_{\rm p}$ until a critical value ($e_{\rm p}^{\rm c}$) that is given by $\simeq$ 5 times the Hill radius scaled by the planet's semimajor axis. For $e_{\rm p} > e_{\rm p}^{\rm c}$, the probability decreases approximately in proportion to $e_{\rm p}^{-1}$. The current orbital eccentricity of Mars is several times larger than $e_{\rm p}^{\rm c}$. However, since the range of secular change in Martian eccentricity overlaps $e_{\rm p}^{\rm c}$, the capture of minor bodies by the past Mars is not ruled out.

Content: For up to a few millions of years, pebbles must provide a quasi-steady inflow of solids from the outer parts of protoplanetary disks to their inner regions. Aims: We wish to understand how a significant fraction of the pebbles grows into planetesimals instead of being lost to the host star. Methods:We examined analytically how the inward flow of pebbles is affected by the snow line and under which conditions dust-rich (rocky) planetesimals form. When calculating the inward drift of solids that is due to gas drag, we included the back-reaction of the gas to the motion of the solids. Results: We show that in low-viscosity protoplanetary disks (with a monotonous surface density similar to that of the minimum-mass solar nebula), the flow of pebbles does not usually reach the required surface density to form planetesimals by streaming instability. We show, however, that if the pebble-to-gas-mass flux exceeds a critical value, no steady solution can be found for the solid-to-gas ratio. This is particularly important for low-viscosity disks (alpha < 10^(-3)) where we show that inside of the snow line, silicate-dust grains ejected from sublimating pebbles can accumulate, eventually leading to the formation of dust-rich planetesimals directly by gravitational instability. Conclusions: This formation of dust-rich planetesimals may occur for extended periods of time, while the snow line sweeps from several au to inside of 1 au. The rock-to-ice ratio may thus be globally significantly higher in planetesimals and planets than in the central star.

It is generally accepted that silicate-metal (`rocky') planet formation relies on coagulation from a mixture of sub-Mars sized planetary embryos and (smaller) planetesimals that dynamically emerge from the evolving circum-solar disc in the first few million years of our Solar System. Once the planets have, for the most part, assembled after a giant impact phase, they continue to be bombarded by a multitude of planetesimals left over from accretion. Here we place limits on the mass and evolution of these planetesimals based on constraints from the highly siderophile element (HSE) budget of the Moon. Outcomes from a combination of N-body and Monte Carlo simulations of planet formation lead us to four key conclusions about the nature of this early epoch. First, matching the terrestrial to lunar HSE ratio requires either that the late veneer on Earth consisted of a single lunar-size impactor striking the Earth before 4.45 Ga, or that it originated from the impact that created the Moon. An added complication is that analysis of lunar samples indicates the Moon does not preserve convincing evidence for a late veneer like Earth. Second, the expected chondritic veneer component on Mars is 0.06 weight percent. Third, the flux of terrestrial impactors must have been low ( <=10^(-6) M_earth/Myr) to avoid wholesale melting of Earth's crust after 4.4~Ga, and to simultaneously match the number of observed lunar basins. This conclusion leads to an Hadean eon which is more clement than assumed previously. Last, after the terrestrial planets had fully formed, the mass in remnant planetesimals was ~10^(-3) M_earth, lower by at least an order of magnitude than most previous models suggest. Our dynamically and geochemically self-consistent scenario requires that future N-body simulations of rocky planet formation either directly incorporate collisional grinding or rely on pebble accretion.

We report $\sim$3 au resolution imaging observations of the protoplanetary disk around TW Hya at 145 and 233 GHz with the Atacama Large Millimeter/Submillimeter Array. Our observations revealed two deep gaps ($\sim$25--50 \%) at 22 and 37~au and shallower gaps (a few \%) at 6, 28, and 44~au, as recently reported by Andrews et al. (2016). The central hole with a radius of $\sim3$~au was also marginally resolved. The most remarkable finding is that the spectral index $\alpha (R)$ between bands 4 and 6 peaks at the 22~au gap. The derived power-law index of the dust opacity $\beta (R)$ is $\sim1.7$ at the 22~au gap and decreases toward the disk center to $\sim0$. The most prominent gap at 22~au could be caused by the gravitational interaction between the disk and an unseen planet with a mass of $\lesssim$1.5 $M_\mathrm{Neptune}$ although other origins may be possible. The planet-induced gap is supported by the fact that $\beta (R)$ is enhanced at the 22~au gap, indicating a deficit of $\sim$mm-sized grains within the gap due to dust filtration by a planet.

We calculate the microlensing event rate and typical time-scales for the free-floating planet (FFP) population that is predicted by the core accretion theory of planet formation. The event rate is found to be ~$1.8\times 10^{-3}$ of that for the stellar population. While the stellar microlensing event time-scale peaks at around 20 days, the median time-scale for FFP events (~0.1 day) is much shorter. Our values for the event rate and the median time-scale are significantly smaller than those required to explain the \cite{Sum+11} result, by factors of ~13 and ~16, respectively. The inclusion of planets at wide separations does not change the results significantly. This discrepancy may be too significant for standard versions of both the core accretion theory and the gravitational instability model to explain satisfactorily. Therefore, either a modification to the planet formation theory is required, or other explanations to the excess of short-time-scale microlensing events are needed. Our predictions can be tested by ongoing microlensing experiment such as KMTNet, and by future satellite missions such as WFIRST and Euclid.

We investigate the simultaneous evolution of dust and gas density profiles at a radial pressure bump located in a protoplanetary disk. If dust particles are treated as test particles, a radial pressure bump traps dust particles that drift radially inward. As the dust particles become more concentrated at the gas pressure bump, however, the drag force from dust to gas (back-reaction), which is ignored in a test-particle approach, deforms the pressure bump. We find that the pressure bump is completely deformed by the back-reaction when the dust-to-gas mass ratio reaches $\sim 1$ for a slower bump restoration. The direct gravitational instability of dust particles is inhibited by the bump destruction. In the dust-enriched region, the radial pressure support becomes $\sim 10-100$ times lower than the global value set initially. Although the pressure bump is a favorable place for streaming instability (SI), the flattened pressure gradient inhibits SI from forming large particle clumps corresponding to $100-1000$ km sized bodies, which has been previously proposed. If SI occurs there, the dust clumps formed would be $10-100$ times smaller, that is, of about $1 - 100$ km.

Context. The classical "planetesimal" accretion scenario for the formation of planets has recently evolved with the idea that "pebbles", centimeter- to meter-sized icy grains migrating in protoplanetary disks, can control planetesimal and/or planetary growth. Aims. We investigate how pebble accretion depends on disk properties and affects the formation of planetary systems Methods. We construct analytical models of pebble accretion onto planetary embryos that consistently account for the mass and orbital evolution of the pebble flow and reflect disk structure. Results. We derive simple formulas for pebble accretion rates in the so-called "settling" regime for planetary embryos with more than 100 km in size. For relatively smaller embryos or in outer disk regions, the accretion mode is 3D, meaning that the thickness of the pebble flow must be taken into account, and resulting in an accretion rate that is independent of the embryo mass. For larger embryos or in inner regions, the accretion is in a 2D mode, i.e., the pebble disk may be considered to be infinitely thin. We show that the radial dependence of the pebble accretion rata is different (even the sign of the power-law exponent changes) for different disk conditions such as the disk heating source (viscous heating or stellar irradiation), the drag law (Stokes or Epstein, and weak or strong coupling), and in the 2D or 3D accretion modes. We also discuss the effect of the sublimation/destruction of icy pebbles inside the snow line. Conclusions. Pebble accretion easily produces a large diversity of planetary systems. In other words, to correctly infer the results of planet formation through pebble accretion, detailed prescriptions of disk evolution and pebble growth, sublimation/destruction and migration are required.

A radial velocity (RV) survey for intermediate-mass giants has been operated for over a decade at Okayama Astrophysical Observatory (OAO). The OAO survey has revealed that some giants show long-term linear RV accelerations (RV trends), indicating the presence of outer companions. Direct imaging observations can help clarify what objects generate these RV trends. We present the results of high-contrast imaging observations or six intermediate-mass giants with long-term RV trends using the Subaru Telescope and HiCIAO camera. We detected co-moving companions to $\gamma$ Hya B ($0.61^{+0.12}_{-0.14} M_\odot$), HD 5608 B ($0.10 \pm 0.01 M_\odot$), and HD 109272 B ($0.28 \pm 0.06 M_\odot$). For the remaining targets($\iota$ Dra, 18 Del, and HD 14067) we exclude companions more massive than 30-60 $M_\mathrm{Jup}$ at projected separations of 1arcsec-7arcsec. We examine whether these directly imaged companions or unidentified long-period companions can account for the RV trends observed around the six giants. We find that the Kozai mechanism can explain the high eccentricity of the inner planets $\iota$ Dra b, HD 5608 b, and HD 14067 b.

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.

We report the first detection of a gap and a ring in 336 GHz dust continuum emission from the protoplanetary disk around TW Hya, using the Atacama Large Millimeter/Submillimeter Array (ALMA). The gap and ring are located at around 25 and 41 au from the central star, respectively, and are associated with the CO snow line at ~30 au. The gap has a radial width of less than 15 au and a mass deficit of more than 23%, taking into account that the observations are limited to an angular resolution of ~15 au. In addition, the 13CO and C18O J = 3 - 2 lines show a decrement in CO line emission throughout the disk, down to ~10 au, indicating a freeze-out of gas-phase CO onto grain surfaces and possible subsequent surface reactions to form larger molecules. The observed gap could be caused by gravitational interaction between the disk gas and a planet with a mass less than super-Neptune (2M_{Neptune}), or could be the result of the destruction of large dust aggregates due to the sintering of CO ice.

Standard accretion disk models suggest that the snow line in the solar nebula migrated interior to the Earth's orbit in a late stage of nebula evolution. In this late stage, a significant amount of ice could have been delivered to 1 AU from outer regions in the form of mm to dm-sized pebbles. This raises the question why the present Earth is so depleted of water (with the ocean mass being as small as 0.023% of the Earth mass). Here we quantify the amount of icy pebbles accreted by terrestrial embryos after the migration of the snow line assuming that no mechanism halts the pebble flow in outer disk regions. We use a simplified version of the coagulation equation to calculate the formation and radial inward drift of icy pebbles in a protoplanetary disk. The pebble accretion cross section of an embryo is calculated using analytic expressions presented by recent studies. We find that the final mass and water content of terrestrial embryos strongly depends on the radial extent of the gas disk, the strength of disk turbulence, and the time at which the snow lines arrives at 1 AU. The disk's radial extent sets the lifetime of the pebble flow, while turbulence determines the density of pebbles at the midplane where the embryos reside. We find that the final water content of the embryos falls below 0.023 wt% only if the disk is compact (< 100 AU), turbulence is strong at 1 AU, and the snow line arrives at 1 AU later than 2-4 Myr after disk formation. If the solar nebula extended to 300 AU, initially rocky embryos would have evolved into icy planets of 1-10 Earth masses unless the snow-line migration was slow. If the proto-Earth contained water of ~ 1 wt% as might be suggested by the density deficit of the Earth's outer core, the formation of the proto-Earth was possible with weaker turbulence and with earlier (> 0.5-2 Myr) snow-line migration.

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.

At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free surface such as a planetary surface. This difficulty comes from the assumption of differentiability of density in SSPH. We have developed an alternative formulation of SPH, density independent SPH (DISPH), which is based on differentiability of pressure instead of density to solve the problem of a contact discontinuity. In this paper, we report the results of the GI simulations with DISPH and compare them with those obtained with SSPH. We found that the disk properties, such as mass and angular momentum produced by DISPH is different from that of SSPH. In general, the disks formed by DISPH are more compact: while formation of a smaller mass moon for low-oblique impacts is expected with DISPH, inhibition of ejection would promote formation of a larger mass moon for high-oblique impacts. Since only the improvement of core-mantle boundary significantly affects the properties of circumplanetary disks generated by GI and DISPH has not been significantly improved from SSPH for a free surface, we should be very careful when some conclusions are drawn from the numerical simulations for GI. And it is necessary to develop the numerical hydrodynamical scheme for GI that can properly treat the free surface as well as the contact discontinuity.

We report the detection of a double planetary system around the evolved intermediate-mass star HD 47366 from precise radial-velocity measurements at Okayama Astrophysical Observatory, Xinglong Station, and Australian Astronomical Observatory. The star is a K1 giant with a mass of 1.81+-0.13M_sun, a radius of 7.30+-0.33R_sun, and solar metallicity. The planetary system is composed of two giant planets with minimum mass of 1.75^{+0.20}_{-0.17}Mjup and 1.86^{+0.16}_{-0.15}Mjup, orbital period of 363.3^{+2.5}_{-2.4} d and 684.7^{+5.0}_{-4.9} d, and eccentricity of 0.089^{+0.079}_{-0.060} and 0.278^{+0.067}_{-0.094}, respectively, which are derived by a double Keplerian orbital fit to the radial-velocity data. The system adds to the population of multi-giant-planet systems with relatively small orbital separations, which are preferentially found around evolved intermediate-mass stars. Dynamical stability analysis for the system revealed, however, that the best-fit orbits are unstable in the case of a prograde configuration. The system could be stable if the planets were in 2:1 mean-motion resonance, but this is less likely considering the observed period ratio and eccentricity. A present possible scenario for the system is that both of the planets have nearly circular orbits, namely the eccentricity of the outer planet is less than ~0.15, which is just within 1.4sigma of the best-fit value, or the planets are in a mutually retrograde configuration with a mutual orbital inclination larger than 160 degree.

Recent observations started revealing the compositions of protostellar discs and planets beyond the Solar System. In this paper, we explore how the compositions of terrestrial planets are affected by dynamical evolution of giant planets. We estimate the initial compositions of building blocks of these rocky planets by using a simple condensation model, and numerically study the compositions of planets formed in a few different formation models of the Solar System. We find that the abundances of refractory and moderately volatile elements are nearly independent of formation models, and that all the models could reproduce the abundances of these elements of the Earth. The abundances of atmophile elements, on the other hand, depend on the scattering rate of icy planetesimals into the inner disc as well as the mixing rate of the inner planetesimal disc. For the classical formation model, neither of these mechanisms are efficient and the accretion of atmophile elements during the final assembly of terrestrial planets appears to be difficult. For the Grand Tack model, both of these mechanisms are efficient, which leads to a relatively uniform accretion of atmophile elements in the inner disc. It is also possible to have a "hybrid" scenario where the mixing is not very efficient but the scattering is efficient. The abundances of atmophile elements in this case increases with orbital radii. Such a scenario may occur in some of the extrasolar planetary systems which are not accompanied by giant planets or those without strong perturbations from giants. We also confirm that the Grand Tack scenario leads to the distribution of asteroid analogues where rocky planetesimals tend to exist interior to icy ones, and show that their overall compositions are consistent with S-type and C-type chondrites, respectively.

We have investigated the dependence of the prograde/retrograde temporary capture of asteroids by a planet on their original heliocentric semimajor axes through analytical arguments and numerical orbital integrations in order to discuss the origins of irregular satellites of giant planets. We found that capture is mostly retrograde for the asteroids near the planetary orbit and is prograde for those from further orbits. An analytical investigation reveals the intrinsic dynamics of these dependences and gives boundary semimajor axes for the change in prograde/retrograde capture. The numerical calculations support the idea of deriving the analytical formulae and confirm their dependence. Our numerical results show that the capture probability is much higher for bodies from the inner region than for outer ones. These results imply that retrograde irregular satellites of Jupiter are most likely to be captured bodies from the nearby orbits of Jupiter that may have the same origin as Trojan asteroids, while prograde irregular satellites originate from far inner regions such as the main-belt asteroid region.

A better knowledge of Jovian satellites' origins will bring light on the environment that surrounded Jupiter during its formation and can help us to understand the characteristics of this unique satellite system. We developed a semi-analytical model to investigate Jupiter's regular satellite formation and present the results of our population synthesis calculations. We performed simulations adopting a massive, static, low-viscosity circumplanetary disk model, in agreement with a current study of magnetorotational instability in a circum-planetary disk. We find that the high gas density leads to very rapid migration of satellitesimals due to gas drag and type II migration of satellites in a faster disk-dominated mode. A large concentration of solids, large building blocks and longer type II migration time-scales favor formation and survival of large satellites. However, bodies as massive as Ganymede and those located far away from Jupiter, such as Callisto, are difficult to form with this scenario.

No Hadean rocks have ever been found on Earth's surface except for zircons---evidence of continental crust, suggesting that Hadean continental crust existed but later disappeared. One hypothesis for the disappearance of the continental crust is excavation/melting by the Late Heavy Bombardment (LHB), a concentration of impacts in the last phase of the Hadean eon. In this paper, we calculate the effects of LHB on Hadean continental crust in order to investigate this hypothesis. Approximating the size-frequency distribution of the impacts by a power-law scaling with an exponent {\alpha} as a parameter, we have derived semi-analytical expressions for the effects of LHB impacts. We calculated the total excavation/melting volume and area affected by the LHB from two constraints of LHB on the moon, the size of the largest basin during LHB, and the density of craters larger than 20 km. We also investigated the effects of the value of {\alpha}. Our results show that LHB does not excavate/melt all of Hadean continental crust directly, but over 70% of the Earth's surface area can be covered by subsequent melts in a broad range of {\alpha}. If there have been no overturns of the continental crust until today, LHB could be responsible for the absence of Hadean rocks because most of Hadean continental crust is not be exposed on the Earth's surface in this case.

Most of super-Earths detected by the radial velocity (RV) method have significantly smaller eccentricities than the eccentricities corresponding to velocity dispersion equal to their surface escape velocity ("escape eccentricities"). If orbital instability followed by giant impacts among protoplanets that have migrated from outer region is considered, it is usually considered that eccentricities of the merged bodies become comparable to those of orbital crossing bodies, which are excited up to their escape eccentricities by close scattering. However, the eccentricity evolution in the {\it in situ} accretion model has not been studied in detail. Here, we investigate the eccentricity evolution through {\it N}-body simulations. We have found that the merged planets tend to have much smaller eccentricities than the escape eccentricities due to very efficient collision damping. If the protoplanet orbits are initially well separated and their eccentricities are securely increased, an inner protoplanet collides at its apocenter with an outer protoplanet at its pericenter. The eccentricity of the merged body is the smallest for such configuration. Orbital inclinations are also damped by this mechanism and planets tend to share a same orbital plane, which is consistent with {\it Kepler} data. Such efficient collision damping is not found when we start calculations from densely packed orbits of the protoplanets. If the protoplanets are initially in the mean-motion resonances, which corresponds to well separated orbits, the {\it in situ} accretion model well reproduces the features of eccentricities and inclinations of multiple super-Earths/Earth systems discovered by RV and {\it Kepler} surveys.

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.

We report detections of new exoplanets from a radial velocity (RV) survey of metal-rich FGK stars by using three telescopes. By optimizing our RV analysis method to long time-baseline observations, we have succeeded in detecting five new Jovian-planets around three metal-rich stars HD 1605, HD 1666, and HD 67087 with the masses of $1.3 M_{\odot}$, $1.5 M_{\odot}$, and $1.4 M_{\odot}$, respectively. A K1 subgiant star HD 1605 hosts two planetary companions with the minimum masses of $ M_p \sin i = 0.96 M_{\mathrm{JUP}}$ and $3.5 M_{\mathrm{JUP}}$ in circular orbits with the planets' periods $P = 577.9$ days and $2111$ days, respectively. HD 1605 shows a significant linear trend in RVs. Such a system consisting of Jovian planets in circular orbits has rarely been found and thus HD 1605 should be an important example of a multi-planetary system that is likely unperturbed by planet-planet interactions. HD 1666 is a F7 main sequence star which hosts an eccentric and massive planet of $ M_p \sin i = 6.4 M_{\mathrm{JUP}}$ in the orbit with $a_{\rm p} = 0.94$ AU and an eccentricity $e=0.63$. Such an eccentric and massive planet can be explained as a result of planet-planet interactions among Jovian planets. While we have found the large residuals of $\mathrm{rms} = 35.6\ \mathrm{m\ s^{-1}}$, the periodogram analysis does not support any additional periodicities. Finally, HD 67087 hosts two planets of $ M_p \sin i = 3.1 M_{\mathrm{JUP}}$ and $4.9 M_{\mathrm{JUP}}$ in orbits with $P=352.2$ days and $2374$ days, and $e=0.17$ and $0.76$, respectively. Although the current RVs do not lead to accurate determinations of its orbit and mass, HD 67087 c can be one of the most eccentric planets ever discovered in multiple systems.

Although the two moons of Mars, Phobos and Deimos, have long been thought to be captured asteroids, recent observations of their compositions and orbits suggest that they may have formed from debris generated by one or more giant impacts of bodies with ~ 0.01 x target mass. Recent studies have both analytically estimated debris produced by giant impacts on Mars and numerically examined the evolution of circum-Mars debris disks. We perform a numerical study (Smoothed Particle Hydrodynamics simulation) of debris retention from giant impacts onto Mars, particularly in relation to a Borealis-scale giant impact (E ~ 3 x 10^29 J) capable of producing the Borealis basin. We find that a Borealis-scale impact is capable of producing a disk of mass ~ 5 x 10^20 kg (~ 1 - 4 % of the impactor mass), sufficient debris to form at least one of the martian moons according to recent numerical studies of martian debris disk evolution. While a Borealis-scale impact may generate sufficient debris to form both Phobos and Deimos, further studies of the debris disk evolution are necessary. Our results can serve as inputs for future studies of martian debris disk evolution.

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.

The Transiting Exoplanet Survey Satellite (TESS) will search for planets transiting bright and nearby stars. TESS has been selected by NASA for launch in 2017 as an Astrophysics Explorer mission. The spacecraft will be placed into a highly elliptical 13.7-day orbit around the Earth. During its two-year mission, TESS will employ four wide-field optical CCD cameras to monitor at least 200,000 main-sequence dwarf stars with I = 4-13 for temporary drops in brightness caused by planetary transits. Each star will be observed for an interval ranging from one month to one year, depending mainly on the star's ecliptic latitude. The longest observing intervals will be for stars near the ecliptic poles, which are the optimal locations for follow-up observations with the James Webb Space Telescope. Brightness measurements of preselected target stars will be recorded every 2 min, and full frame images will be recorded every 30 min. TESS stars will be 10-100 times brighter than those surveyed by the pioneering Kepler mission. This will make TESS planets easier to characterize with follow-up observations. TESS is expected to find more than a thousand planets smaller than Neptune, including dozens that are comparable in size to the Earth. Public data releases will occur every four months, inviting immediate community-wide efforts to study the new planets. The TESS legacy will be a catalog of the nearest and brightest stars hosting transiting planets, which will endure as highly favorable targets for detailed investigations.

Context. Circumstellar disks are known to contain a significant mass in dust ranging from micron to centimeter size. Meteorites are evidence that individual grains of those sizes were collected and assembled into planetesimals in the young solar system. Aims. We assess the efficiency of dust collection of a swarm of non-drifting planetesimals {\rev with radii ranging from 1 to $10^3$\,km and beyond. Methods. We calculate the collision probability of dust drifting in the disk due to gas drag by planetesimal accounting for several regimes depending on the size of the planetesimal, dust, and orbital distance: the geometric, Safronov, settling, and three-body regimes. We also include a hydrodynamical regime to account for the fact that small grains tend to be carried by the gas flow around planetesimals. Results. We provide expressions for the collision probability of dust by planetesimals and for the filtering efficiency by a swarm of planetesimals. For standard turbulence conditions (i.e., a turbulence parameter $\alpha=10^{-2}$), filtering is found to be inefficient, meaning that when crossing a minimum-mass solar nebula (MMSN) belt of planetesimals extending between 0.1 AU and 35 AU most dust particles are eventually accreted by the central star rather than colliding with planetesimals. However, if the disk is weakly turbulent ($\alpha=10^{-4}$) filtering becomes efficient in two regimes: (i) when planetesimals are all smaller than about 10 km in size, in which case collisions mostly take place in the geometric regime; and (ii) when planetary embryos larger than about 1000 km in size dominate the distribution, have a scale height smaller than one tenth of the gas scale height, and dust is of millimeter size or larger in which case most collisions take place in the settling regime. These two regimes have very different properties: we find that the local filtering efficiency $x_{filter,MMSN}$ scales with $r^{-7/4}$ (where $r$ is the orbital distance) in the geometric regime, but with $r^{-1/4}$ to $r^{1/4}$ in the settling regime. This implies that the filtering of dust by small planetesimals should occur close to the central star and with a short spread in orbital distances. On the other hand, the filtering by embryos in the settling regime is expected to be more gradual and determined by the extent of the disk of embryos. Dust particles much smaller than millimeter size tend only to be captured by the smallest planetesimals because they otherwise move on gas streamlines and their collisions take place in the hydrodynamical regime. Conclusions. Our results hint at an inside-out formation of planetesimals in the infant solar system because small planetesimals in the geometrical limit can filter dust much more efficiently close to the central star. However, even a fully-formed belt of planetesimals such as the MMSN only marginally captures inward-drifting dust and this seems to imply that dust in the protosolar disk has been filtered by planetesimals even smaller than 1 km (not included in this study) or that it has been assembled into planetesimals by other mechanisms (e.g., orderly growth, capture into vortexes). Further refinement of our work concerns, among other things: a quantitative description of the transition region between the hydro and settling regimes; an assessment of the role of disk turbulence for collisions, in particular in the hydro regime; and the coupling of our model to a planetesimal formation model.

We report the detection of a substellar companion orbiting an evolved intermediate-mass ($M_\star=2.4\,M_\odot$) star HD 14067 (G9 III) using precise Doppler technique. Radial velocities of this star can be well fitted either by a periodic Keplerian variation with a decreasing linear velocity trend (P=1455 days, $K_1=92.2$ m s$^{-1}$, $e=0.533$, and $\dot{\gamma}=-22.4$ m s$^{-1}$ yr$^{-1}$) or a single Keplerian orbit without linear trend (P=2850 days, $K_1=100.1$ m s$^{-1}$, and $e=0.697$). The minimum mass ($m_2\sin{i}=7.8\,M_{\rm J}$ for the model with a linear trend, or $m_2\sin{i}=9.0\,M_{\rm J}$ for the model without a linear trend) suggests a long-period giant planet around an evolved intermediate-mass star. The eccentricity of the orbit is among the highest known for planets ever detected around evolved stars.

We have investigated i) the formation of gravitationally bounded pairs of gas-giant planets (which we call "binary planets") 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.

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'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.

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.

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.

Planetary migration is one of the most serious problems to systematically understand the observations of exoplanets. We clarify that the theoretically predicted type II migration is too fast, as well as type I migration, by developing detailed analytical arguments in which the timescale of type II migration is compared with the disk lifetime. In the disk-dominated regime, the type II migration timescale is characterized by a {\it local} viscous diffusion timescale, while the disk lifetime characterized by a {\it global} diffusion timescale that is much longer than the local one. Even in the planet-dominated regime where the inertia of the planet mass reduces the migration speed, the timescale is still shorter than the disk lifetime except in the final disk evolution stage where the total disk mass decays below the planet mass. This suggests that most giant planets plunge into the central stars within the disk lifetime, and it contradicts the exoplanet observations that gas giants are piled up at $r \ga 1$AU. We examine additional processes that may arise in protoplanetary disks: dead zones, photoevaporation of gas, and gas flow across a gap formed by a type II migrator. Although they make the type II migration timescale closer to the disk lifetime, we show that none of them can act as an effective barrier for rapid type II migration with the current knowledge of these processes. We point out that gas flow across a gap and the fraction of the flow accreted onto the planets are uncertain and they may have a potential to solve the problem. Much more detailed investigation for each process may be needed to explain the observed distribution of gas giants in extrasolar planetary systems.

We have carried out high contrast imaging of 70 young, nearby B and A stars to search for brown dwarf and planetary companions as part of the Gemini NICI Planet-Finding Campaign. Our survey represents the largest, deepest survey for planets around high-mass stars (~1.5-2.5 M_sun) conducted to date and includes the planet hosts beta Pic and Fomalhaut. We obtained follow-up astrometry of all candidate companions within 400 AU projected separation for stars in uncrowded fields and identified new low-mass companions to HD 1160 and HIP 79797. We have found that the previously known young brown dwarf companion to HIP 79797 is itself a tight (3 AU) binary, composed of brown dwarfs with masses 58 (+21, -20) M_Jup and 55 (+20, -19) M_Jup, making this system one of the rare substellar binaries in orbit around a star. Considering the contrast limits of our NICI data and the fact that we did not detect any planets, we use high-fidelity Monte Carlo simulations to show that fewer than 20% of 2 M_sun stars can have giant planets greater than 4 M_Jup between 59 and 460 AU at 95% confidence, and fewer than 10% of these stars can have a planet more massive than 10 M_Jup between 38 and 650 AU. Overall, we find that large-separation giant planets are not common around B and A stars: fewer than 10% of B and A stars can have an analog to the HR 8799 b (7 M_Jup, 68 AU) planet at 95% confidence. We also describe a new Bayesian technique for determining the ages of field B and A stars from photometry and theoretical isochrones. Our method produces more plausible ages for high-mass stars than previous age-dating techniques, which tend to underestimate stellar ages and their uncertainties.

The orbital distributions of currently observed extrasolar giant planets allow marginally stable orbits for hypothetical, terrestrial planets. In this paper, we propose that many of these systems may not have additional planets on these "stable" orbits, since past dynamical instability among giant planets could have removed them. We numerically investigate the effects of early evolution of multiple giant planets on the orbital stability of the inner, sub-Neptune-like planets which are modeled as test particles, and determine their dynamically unstable region. Previous studies have shown that the majority of such test particles are ejected out of the system as a result of close encounters with giant planets. Here, we show that secular perturbations from giant planets can remove test particles at least down to 10 times smaller than their minimum pericenter distance. Our results indicate that, unless the dynamical instability among giant planets is either absent or quiet like planet-planet collisions, most test particles down to ~0.1 AU within the orbits of giant planets at a few AU may be gone. In fact, out of ~30 % of survived test particles, about three quarters belong to the planet-planet collision cases. We find a good agreement between our numerical results and the secular theory, and present a semi-analytical formula which estimates the dynamically unstable region of the test particles just from the evolution of giant planets. Finally, our numerical results agree well with the observations, and also predict the existence of hot rocky planets in eccentric giant planet systems.

We report the detections of planetary companions orbiting around three evolved intermediate-mass stars from precise radial velocity measurements at Okayama Astrophysical Observatory. HD 2952 (K0III, 2.5 M_sun) and omega Ser (G8III, 2.2 M_sun) host a relatively low mass planet with minimum mass of m_2sin i=1.6 M_J and 1.7 M_J in nearly circular orbits with period of P=312 and 277 d, respectively. HD 120084 (G7 III, 2.4 M_sun) hosts an eccentric planet with m_2sin i=4.5 M_J in an orbit with P=2082 d and eccentricity of e=0.66. The planet has one of the largest eccentricities among those ever discovered around evolved intermediate-mass stars, almost all of which have eccentricity smaller than 0.4. We also show that radial velocity variations of stellar oscillations for G giants can be averaged out below a level of a few m/s at least in timescale of a week by high cadence observations, which enables us to detect a super-Earth and a Neptune-mass planet in short-period orbits even around such giant stars.

We review the basic dynamics and accretion of planetesimals by showing N-body simulations. The orbits of planetesimals evolve through two-body gravitational relaxation: viscous stirring increases the random velocity and dynamical friction realizes the equiparation of the random energy. In the early stage of planetesimal accretion the growth mode of planetesimals is runaway growth where larger planetesimals grow faster than smaller ones. When a protoplanet (runaway-growing planetesimal) exceeds a critical mass the growth mode shifts to oligarchic growth where similar-sized protoplanets grow keeping a certain orbital separation. The final stage of terrestrial planet formation is collision among protoplanets known as giant impacts. We also summarize the dynamical effects of disk gas on planets and the core accretion model for formation of gas giants and discuss the diversity of planetary systems.

We report the detection of a double planetary system orbiting around the evolved intermediate-mass star HD 4732 from precise Doppler measurements at Okayama Astrophysical Observatory (OAO) and Anglo-Australian Observatory (AAO). The star is a K0 subgiant with a mass of 1.7 M_sun and solar metallicity. The planetary system is composed of two giant planets with minimum mass of msini=2.4 M_J, orbital period of 360.2 d and 2732 d, and eccentricity of 0.13 and 0.23, respectively. Based on dynamical stability analysis for the system, we set the upper limit on the mass of the planets to be about 28 M_J (i>5 deg) in the case of coplanar prograde configuration.

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%.

We report that HAT-P-7 has a common proper motion stellar companion. The companion is located at $\sim3.9$ arcsec to the east and estimated as an M5.5V dwarf based on its colors. We also confirm the presence of the third companion, which was first reported by Winn et al. (2009), based on long-term radial velocity measurements. We revisit the migration mechanism of HAT-P-7b given the presence of those companions, and propose sequential Kozai migration as a likely scenario in this system. This scenario may explain the reason for an outlier in the discussion of the spin-orbit alignment timescale for HAT-P-7b by Albrecht et al. (2012).

Many extrasolar planetary systems containing multiple super-Earths have been discovered. N-body simulations taking into account standard type-I planetary migration suggest that protoplanets are captured into mean-motion resonant orbits near the inner disk edge at which the migration is halted. Previous N-body simulations suggested that orbital stability of the resonant systems depends on number of the captured planets. In the unstable case, through close scattering and merging between planets, non-resonant multiple systems are finally formed. In this paper, we investigate the critical number of the resonantly trapped planets beyond which orbital instability occurs after disk gas depletion. We find that when the total number of planets ($N$) is larger than the critical number ($N_{\rm crit}$), crossing time that is a timescale of initiation of the orbital instability is similar to non-resonant cases, while the orbital instability never occurs within the orbital calculation time ($10^8$ Kepler time) for $N\leq N_{\rm crit}$. Thus, the transition of crossing time across the critical number is drastic. When all the planets are trapped in 7:6 resonance of adjacent pairs, $N_{\rm crit} = 4$. We examine the dependence of the critical number of 4:3, 6:5 and 8:7 resonance by changing the orbital separation in mutual Hill radii and planetary mass. The critical number increases with increasing the orbital separation in mutual Hill radii with fixed planetary mass and increases with increasing planetary mass with fixed the orbital separation in mutual Hill radii. We also calculate the case of a system which is not composed of the same resonance. The sharp transition of the stability can be responsible for the diversity of multiple super-Earths (non-resonant or resonant), that is being revealed by $Kepler$ mission.

Planets migrate due to the recoil they experience from scattering solid (planetesimal) bodies. To first order, the torques exerted by the interior and exterior disks cancel, analogous to the cancellation of the torques from the gravitational interaction with the gas (type I migration). Assuming the dispersion-dominated regime and power-laws characterized by indices {\alpha} and {\beta} for the surface density and eccentricity profiles, we calculate the net torque on the planet. We consider both distant encounters and close (orbit-crossing) encounters. We find that the close and distant encounter torques have opposite signs with respect to their {\alpha} and {\beta} dependences; and that the torque is especially sensitive to the eccentricity gradient ({\beta}). Compared to type-I migration due to excitation of density waves, the planetesimal-driven migration rate is generally lower due to the lower surface density of solids in gas-rich disk, although this may be partially or fully offset when their eccentricity and inclination are small. Allowing for the feedback of the planet on the planetesimal disk through viscous stirring, we find that under certain conditions a self-regulated migration scenario emerges, in which the planet migrates at a steady pace that approaches the rate corresponding to the one-sided torque. If the local planetesimal disk mass to planet mass ratio is low, however, migration stalls. We quantify the boundaries separating the three migration regimes.

We report the detections of substellar companions orbiting around seven evolved intermediate-mass stars from precise Doppler measurements at Okayama Astrophysical Observatory. o UMa (G4 II-III) is a giant with a mass of 3.1 M_sun and hosts a planet with minimum mass of m_2sini=4.1 M_J in an orbit with a period P=1630 d and an eccentricity e=0.13. This is the first planet candidate (< 13 M_J) ever discovered around stars more massive than 3 M_sun. o CrB (K0 III) is a 2.1 M_sun giant and has a planet of m_2sini=1.5 M_J in a 187.8 d orbit with e=0.19. This is one of the least massive planets ever discovered around ~2 M_sun stars. HD 5608 (K0 IV) is an 1.6 M_sun subgiant hosting a planet of m_2sini=1.4 M_J in a 793 d orbit with e=0.19. The star also exhibits a linear velocity trend suggesting the existence of an outer, more massive companion. 75 Cet (G3 III:) is a 2.5 M_sun giant hosting a planet of m_2sini=3.0 M_J in a 692 d orbit with e=0.12. The star also shows possible additional periodicity of about 200 d and 1880 d with velocity amplitude of ~7--10 m/s, although these are not significant at this stage. nu Oph (K0 III) is a 3.0 M_sun giant and has two brown-dwarf companions of m_2sini= 24 M_J and 27 M_J, in orbits with P=530.3 d and 3190 d, and e=0.126 and 0.17, respectively, which were independently announced by Quirrenbach et al. (2011). The ratio of the periods is close to 1:6, suggesting that the companions are in mean motion resonance. We also independently confirmed planets around k CrB (K0 III-IV) and HD 210702 (K1 IV), which had been announced by Johnson et al. (2008) and Johnson et al. (2007a), respectively. All of the orbital parameters we obtained are consistent with the previous results.

We have investigated Saturn's core formation at a radial pressure maximum in a protoplanetary disk, which is created by gap opening by Jupiter. A core formed via planetesimal accretion induces the fragmentation of surrounding planetesimals, which generally inhibits further growth of the core by removal of the resulting fragments due to radial drift caused by gas drag. However, the emergence of the pressure maximum halts the drift of the fragments, while their orbital eccentricities and inclinations are efficiently damped by gas drag. As a result, the core of Saturn rapidly grows via accretion of the fragments near the pressure maximum. We have found that in the minimum-mass solar nebula, kilometer sized planetesimals can produce a core exceeding 10 Earth masses within two million years. Since Jupiter may not have undergone significant type II inward migration, it is likely that Jupiter's formation was completed when the local disk mass has already decayed to a value comparable to or less than Jovian mass. The expected rapid growth of Saturn's core on a timescale comparable to or shorter than observationally inferred disk lifetime enables Saturn to acquire the current amount of envelope gas before the disk gas is completely depleted. The high heat energy release rate onto the core surface due to the rapid accretion of the fragments delays onset of runaway gas accretion until the core mass becomes somewhat larger than that of Jupiter, which is consistent with the estimate based on interior modeling. Therefore, the rapid formation of Saturn induced by gap opening of Jupiter can account for the formation of multiple gas giants (Jupiter and Saturn) without significant inward migration and larger core mass of Saturn than that of Jupiter.

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.

We estimate minimum dust abundances required for secular gravitational instability (SGI) to operate at the midplane dust layer of protoplanetary disks. For SGI to be a viable process, the growth time of the instability T_grow must be shorter than the radial drift time of the dust T_drift. The growth time depends on the turbulent diffusion parameter alpha, because the modes with short wavelengths are stabilized by turbulent diffusion. Assuming that turbulence is excited via the Kelvin-Helmholtz or streaming instabilities in the dust layer, and that its strength is controlled by the energy supply rate from dust accretion, we estimate the diffusion parameter and the growth time of the instability. The condition T_grow < T_drift requires that the dust abundance must be greater than a critical abundance Z_min, which is a function of the Toomre parameter Q_g and aspect ratio h_g / r of the gas disk. For a wide range of parameter space, the required dust abundance is less than 0.1. A slight increase in dust abundance opens a possible route for the dust to directly collapse to planetesimals.

We have studied formation of planetesimals at a radial pressure bump in a protoplanetary disk created by radially inhomogeneous magnetorotational instability (MRI), through three-dimensional resistive MHD simulations including dust particles. In our previous papers, we showed that the inhomogeneous MRI developing in non-uniform structure of magnetic field or magnetic resistivity can transform the local gas flow in the disk to a quasi-steady state with local rigid rotation that is no more unstable against the MRI. Since the outer part of the rigid rotation is super-Keplerian flow, a quasi-static pressure bump is created and dust concentration is expected there. In this paper, we perform simulations of the same systems, adding dust particles that suffer gas drag and modulate gas flow via the back-reaction of the gas drag (dust drag). We use O(10^7) super-particles, each of which represents many dust particles with sizes of centimeter to meter. We have found that the dust drag suppresses turbulent motion to decrease the velocity dispersion of the dust particles while it broadens the dust concentrated regions to limit peaky dust concentration, compared with the simulation without the dust drag. We found that reduction in the velocity dispersion) is dominated over the suppression in particle concentration. For meter-size particles with the friction time ~1/Omega, where Omega is Keplerian frequency, the gravitational instability of the dust particles that may lead to planetesimal formation is expected. For such a situation, we further introduced the self-gravity of dust particles to the simulation to demonstrate that several gravitationally bound clumps are actually formed. Through analytical arguments, we found that the planetesimal formation from meter-sized dust particles can be possible at ~5AU, if dust spatial density is a few times larger than that in the minimum mass solar nebula.

We study the turbulence induced in the dust layer of a protoplanetary disk based on the energetics of dust accretion due to gas drag. We estimate turbulence strength from the energy supplied by dust accretion, using the radial drift velocity of the dust particles in a laminar disk. Our estimate of the turbulence strength agrees with previous analytical and numerical research on the turbulence induced by Kelvin-Helmholtz and/or streaming instabilities for particles whose stopping time is less than the Keplerian time. For such small particles, the strongest turbulence is expected to occur when the dust-to-gas ratio of the disk is ~C_eff^(1/2) (h_g / r) ~ 10^(-2), where C_eff ~ 0.2 represents the energy supply efficiency to turbulence and h_g / r ~ 5 x 10^(-2) is the aspect ratio of the gas disk. The maximum viscosity parameter is alpha_max ~ C_eff T_s (h_g / r)^2 ~ 10^(-4) T_s, where T_s (<1) is the non-dimensional stopping time of the dust particles. Modification in the dust-to-gas ratio from the standard value, 10^(-2), by any process, results in weaker turbulence and a thinner dust layer, and consequently may accelerate the growth process of the dust particles.

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, $\gamma$, the total mass of two protoplanets, $M_{\rm T}$, the impact angle, $\theta$, 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 $\gamma$ and $\theta$, 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 $\gamma$ and $\theta$, 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.

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$\mu$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.

Evolution of a snow line in an optically-thick protoplanetary disk is investigated with numerical simulations. The ice-condensing region in the disk is obtained by calculating the temperature and the density with the 1+1D approach. The snow line migrates as the mass accretion rate (\dot{M}) in the disk decreases with time. Calculations are carried out from an early phase with high disk accretion rates (\dot{M} \sim 10^{-7} M_sun/yr) to a later phase with low disk accretion rates (\dot{M} \sim 10^{-12} M_sun/yr) using the same numerical method. It is found that the snow line moves inward for \dot{M} > 10^{-10} M_sun/yr, while it gradually moves outward in the later evolution phase with \dot{M} < 10^{-10} M_sun/yr. In addition to the silicate opacity, the ice opacity is taken into consideration. In the inward migration phase, the additional ice opacity increases the distance of the snow line from the central star by a factor of 1.3 for dust grains < 10 micro-meter in size and 1.6 for > 100 micro-meter. It is inevitable that the snow line comes inside the Earth's orbit in the course of the disk evolution, if the alpha viscosity parameter is in a range 0.001-0.1, the dust-to-gas mass ratio is higher than a tenth of the solar abundance value, and the dust grains are smaller than 1 mm. The formation of water-devoid planetesimals in the terrestrial planet region seems to be difficult throughout the disk evolution, which imposes a new challenge to planet formation theory.

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.

We present a new analytic approach to the disk-planet interaction that is especially useful for planets with eccentricity larger than the disk aspect ratio. We make use of the dynamical friction formula to calculate the force exerted on the planet by the disk, and the force is averaged over the period of the planet. The resulting migration and eccentricity damping timescale agrees very well with the previous works in which the planet eccentricity is moderately larger than the disk aspect ratio. The advantage of this approach is that it is possible to apply this formulation to arbitrary large eccentricity. We have found that the timescale of the orbital evolution depends largely on the adopted disk model in the case of highly eccentric planets. We discuss the possible implication of our results to the theory of planet formation.

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.

We present the discovery of a wide (67 AU) substellar companion to the nearby (21 pc) young solar-metallicity M1 dwarf CD-35 2722, a member of the ~100 Myr AB Doradus association. Two epochs of astrometry from the NICI Planet-Finding Campaign confirm that CD-35 2722 B is physically associated with the primary star. Near-IR spectra indicate a spectral type of L4\pm1 with a moderately low surface gravity, making it one of the coolest young companions found to date. The absorption lines and near-IR continuum shape of CD-35 2722 B agree especially well the dusty field L4.5 dwarf 2MASS J22244381-0158521, while the near-IR colors and absolute magnitudes match those of the 5 Myr old L4 planetary-mass companion, 1RXS J160929.1-210524 b. Overall, CD-35 2722 B appears to be an intermediate-age benchmark for L-dwarfs, with a less peaked H-band continuum than the youngest objects and near-IR absorption lines comparable to field objects. We fit Ames-Dusty model atmospheres to the near-IR spectra and find T=1700-1900 K and log(g) =4.5\pm0.5. The spectra also show that the radial velocities of components A and B agree to within \pm10 km/s, further confirming their physical association. Using the age and bolometric luminosity of CD-35 2722 B, we derive a mass of 31\pm8 Mjup from the Lyon/Dusty evolutionary models. Altogether, young late-M to mid-L type companions appear to be over-luminous for their near-IR spectral type compared to field objects, in contrast to the under-luminosity of young late-L and early-T dwarfs.

The questions of how planets form and how common Earth-like planets are can be addressed by measuring the distribution of exoplanet masses and orbital periods. We report the occurrence rate of close-in planets (with orbital periods less than 50 days) based on precise Doppler measurements of 166 Sun-like stars. We measured increasing planet occurrence with decreasing planet mass (M). Extrapolation of a power law mass distribution fitted to our measurements, df/dlogM = 0.39M^-0.48, predicts that 23% of stars harbor a close-in Earth-mass planet (ranging from 0.5 to 2.0 Earth masses). Theoretical models of planet formation predict a deficit of planets in the domain from 5 to 30 Earth masses and with orbital periods less than 50 days. This region of parameter space is in fact well populated, implying that such models need substantial revision.

We have numerically studied the thermal evolution of various-mass terrestrial planets in habitable zones, focusing on duration of dynamo activity to generate their intrinsic magnetic fields, which may be one of key factors in habitability on the planets. In particular, we are concerned with super-Earths, observations of which are rapidly developing. We calculated evolution of temperature distributions in planetary interior, using Vinet equations of state, Arrhenius-type formula for mantle viscosity, and the astrophysical mixing length theory for convective heat transfer modified for mantle convection. After calibrating the model with terrestrial planets in the Solar system, we apply it for 0.1--$10M_{\oplus}$ rocky planets with surface temperature of $300~\mbox{K}$ (in habitable zones) and the Earth-like compositions. With the criterion for heat flux at the CMB (core-mantle boundary), the lifetime of the magnetic fields is evaluated from the calculated thermal evolution. We found that the lifetime slowly increases with the planetary mass ($M_p$) independent of initial temperature gap at the core-mantle boundary ($\Delta T_{\rm CMB}$) but beyond a critical value $M_{c,p}$ ($\sim O(1)M_{\oplus}$) it abruptly declines by the mantle viscosity enhancement due to the pressure effect. We derived $M_{c,p}$ as a function of $\Delta T_{\rm CMB}$ and a rheological parameter (activation volume, $V^*$). Thus, the magnetic field lifetime of super-Earths with $M_p > M_{p,c}$ sensitively depends on $\Delta T_{\rm CMB}$, which reflects planetary accretion, and $V^*$, which has uncertainty at very high pressure. More advanced high-pressure experiments and first-principle simulation as well as planetary accretion simulation are needed to discuss habitability of super-Earths.

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.

Our team is carrying out a multi-year observing program to directly image and characterize young extrasolar planets using the Near-Infrared Coronagraphic Imager (NICI) on the Gemini-South 8.1-meter telescope. NICI is the first instrument on a large telescope designed from the outset for high-contrast imaging, comprising a high-performance curvature adaptive optics system with a simultaneous dual-channel coronagraphic imager. Combined with state-of-the-art observing methods and data processing, NICI typically achieves ~2 magnitudes better contrast compared to previous ground-based or space-based programs, at separations inside of ~2 arcsec. In preparation for the Campaign, we carried out efforts to identify previously unrecognized young stars, to rigorously construct our observing strategy, and to optimize the combination of angular and spectral differential imaging. The Planet-Finding Campaign is in its second year, with first-epoch imaging of 174 stars already obtained out of a total sample of 300 stars. We describe the Campaign's goals, design, implementation, performance, and preliminary results. The NICI Campaign represents the largest and most sensitive imaging survey to date for massive (~1 Mjup) planets around other stars. Upon completion, the Campaign will establish the best measurements to date on the properties of young gas-giant planets at ~5-10 AU separations. Finally, Campaign discoveries will be well-suited to long-term orbital monitoring and detailed spectrophotometric followup with next-generation planet-finding instruments.

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 report the discovery of a tight substellar companion to the young solar analog PZ Tel, a member of the Beta Pictoris moving group observed with high contrast adaptive optics imaging as part of the Gemini NICI Planet-Finding Campaign. The companion was detected at a projected separation of 16.4 +/- 1.0 AU (0.33 +/- 0.01") in April 2009. Second-epoch observations in May 2010 demonstrate that the companion is physically associated and shows significant orbital motion. Monte Carlo modeling constrains the orbit of PZ Tel B to eccentricities > 0.6. The near-IR colors of PZ Tel B indicate a spectral type of M7+/-2 and thus this object will be a new benchmark companion for studies of ultracool, low-gravity photospheres. Adopting an age of 12 +8 -4 Myr for the system, we estimate a mass of 36 +/- 6 Mjup based on the Lyon/DUSTY evolutionary models. PZ Tel B is one of few young substellar companions directly imaged at orbital separations similar to those of giant planets in our own solar system. Additionally, the primary star PZ Tel A shows a 70 um emission excess, evidence for a significant quantity of circumstellar dust that has not been disrupted by the orbital motion of the companion.

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' 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.

We report the detections of two substellar companions orbiting around evolved intermediate-mass stars from precise Doppler measurements at Subaru Telescope and Okayama Astrophysical Observatory. HD 145457 is a K0 giant with a mass of 1.9 M_sun and has a planet of minimum mass m_2sini=2.9 M_J orbiting with period of P=176 d and eccentricity of e=0.11. HD 180314 is also a K0 giant with 2.6 M_sun and hosts a substellar companion of m_2sin i=22 M_J, which falls in brown-dwarf mass regime, in an orbit with P=396 d and e=0.26. HD 145457 b is one of the innermost planets and HD 180314 b is the seventh candidate of brown-dwarf-mass companion found around intermediate-mass evolved stars.

Planet migration is the process by which a planet's orbital radius changes in time. The main agent for causing gas giant planet migration is the gravitational interaction of the young planet with the gaseous disk from which it forms. We describe the migration rates resulting from these interactions based on a simple model for disk properties. These migration rates are higher than is reasonable for planet survival. We discuss some proposed models for which the migration rates are lower. There are major uncertainties in migration rates due to a lack of knowledge about the detailed physical properties of disks. We also describe some additional forms of migration.

We present results of direct imaging observations for HAT-P-7 taken with the Subaru HiCIAO and the Calar Alto AstraLux. Since the close-in transiting planet HAT-P-7b was reported to have a highly tilted orbit, massive bodies such as giant planets, brown dwarfs, or a binary star are expected to exist in the outer region of this system. We show that there are indeed two candidates for distant faint stellar companions around HAT-P-7. We discuss possible roles played by such companions on the orbital evolution of HAT-P-7b. We conclude that as there is a third body in the system as reported by Winn et al. (2009, ApJL, 763, L99), the Kozai migration is less likely while planet-planet scattering is possible.

We report the detection of a large mass planet orbiting around the K0 metal-rich subgiant HD38801 ($V=8.26$) by precise radial velocity (RV) measurements from the Subaru Telescope and the Keck Telescope. The star has a mass of $1.36M_{\odot}$ and metallicity of [Fe/H]= +0.26. The RV variations are consistent with a circular orbit with a period of 696.0 days and a velocity semiamplitude of 200.0\mps, which yield a minimum-mass for the companion of $10.7\mjup$ and semimajor axis of 1.71 AU. Such super-massive objects with very low-eccentricities and hundreds of days period are uncommon among the ensemble of known exoplanets.

The Jovian regular satellite system mainly consists of four Galilean satellites that have similar masses and are trapped in mutual mean motion resonances except for the outer satellite, Callisto. On the other hand, the Saturnian regular satellite system has only one big icy body, Titan, and a population of much smaller icy moons. We have investigated the origin of these major differences between the Jovian and Saturnian satellite systems by semi-analytically simulating the growth and orbital migration of proto-satellites in an accreting proto-satellite disk. We set up two different disk evolution/structure models that correspond to Jovian and Saturnian systems, by building upon previously developed models of an actively-supplied proto-satellite disk, the formation of gas giants, and observations of young stars. Our simulations extend previous models by including the (1) different termination timescales of gas infall onto the proto-satellite disk and (2) different evolution of a cavity in the disk, between the Jovian and Saturnian systems. We have performed Monte Carlo simulations and show that in the case of the Jovian systems, four to five similar-mass satellites are likely to remain trapped in mean motion resonances. This orbital configuration is formed by type I migration, temporal stopping of the migration near the disk inner edge, and quick truncation of gas infall caused by Jupiter opening a gap in the Solar nebula. The Saturnian systems tend to end up with one dominant body in the outer regions caused by the slower decay of gas infall associated with global depletion of the Solar nebula. The total mass and compositional zoning of the predicted Jovian and Saturnian satellite systems are consistent with the observed satellite systems.

How to create planetesimals from tiny dust particles in a proto-planetary disk before the dust particles spiral to the central star is one of the most challenging problems in the theory of planetary system formation. In our previous paper Kato et al. (2009), we have shown that a steady angular velocity profile that consists of both super and sub-Keplerian regions is created in the disk through non-uniform excitation of Magneto-Rotational Instability (MRI). Such non-uniform MRI excitation is reasonably expected in a part of disks with relatively low ionization degree. In this paper, we show through three-dimensional resistive MHD simulations with test particles that this radial structure of the angular velocity indeed leads to prevention of spiral-in of dust particles and furthermore to their accumulation at the boundary of super-Keplerian and sub-Keplerian regions. Treating dust particles as test particles, their motions under the influence of the non-uniform MRI through gas drag are simulated. In the most favorable cases (meter-size dust particles in the disk region with a relatively large fraction of MRI-stable region), we found that the dust concentration is peaked around the super/sub-Keplerian flow boundary and the peak dust density is 10,000 times as high as the initial value. The peak density is high enough for the subsequent gravitational instability to set in, suggesting a possible route to planetesimal formation via non-uniformly excited MRI in weakly ionized regions of a disk.

We report the detection of a substellar companion orbiting the G5 dwarf HD 16760 from the N2K sample. Precise Doppler measurements of the star from Subaru and Keck revealed a Keplerian velocity variation with a period of 466.47+-0.35 d, a semiamplitude of 407.71+-0.84 m/s, and an eccentricity of 0.084+-0.003. Adopting a stellar mass of 0.78+-0.05 M_Sun, we obtain a minimum mass for the companion of 13.13+-0.56 M_JUP, which is close to the planet/brown-dwarf transition, and the semimajor axis of 1.084+-0.023 AU. The nearly circular orbit despite the large mass and intermediate orbital period makes this companion unique among known substellar companions.

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 ask if Earth-like planets (terrestrial mass and habitable-zone orbit) can be detected in multi-planet systems, using astrometric and radial velocity observations. We report here the preliminary results of double-blind calculations designed to answer this question.

The extrasolar planets (EPs) so far detected are very different to the planets in our own Solar System. Many of them have Jupiter-like masses and close-in orbits (the so-called hot planets, HPs), with orbital periods of only a few days. In this paper, we present a new statistical analysis of the observed EPs, focusing on the origin of the HPs. Among the several HP formation mechanisms proposed so far, the two main formation mechanisms are type II migration and scattering. In both cases, planets form beyond the so-called snow-line of the protoplanetary disk and then migrate inward due to angular momentum and energy exchange with either the protoplanetary disk or with companion planets. Although theoretical studies produce a range of observed features, no firm correspondence between the observed EPs and models has yet been established. In our analysis, by means of principal component analysis and hierarchical cluster analysis, we find convincing indications for the existence of two types of HPs, whose parameters reflect physical mechanisms of type II migration and scattering.

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.

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 "desert" in the Mp-a distribution centered around Mp ~ 30-50 M_Earth and a < 1 AU. We show that the expected "desert" 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.

We study the conditions for collisions between planetesimals to be accretional or disruptive in turbulent disks, through analytical arguments based on fluid dynamical simulations and orbital integrations. In turbulent disks, the velocity dispersion of planetesimals is pumped up by random gravitational perturbations from density fluctuations of the disk gas. When the velocity dispersion is larger than the planetesimals' surface escape velocity, collisions between planetesimals do not result in accretion, and may even lead to their destruction. In disks with a surface density equal to that of the ``minimum mass solar nebula'' and with nominal MRI turbulence, we find that accretion proceeds only for planetesimals with sizes above $\sim 300$ km at 1AU and $\sim 1000$ km at 5AU. We find that accretion is facilitated in disks with smaller masses. However, at 5AU and for nominal turbulence strength, km-sized planetesimals are in a highly erosive regime even for a disk mass as small as a fraction of the mass of Jupiter. The existence of giant planets implies that either turbulence was weaker than calculated by standard MRI models or some mechanism was capable of producing Ceres-mass planetesimals in very short timescales. In any case, our results show that in the presence of turbulence planetesimal accretion is most difficult in massive disks and at large orbital distances.

We report on the detection of four extrasolar planets orbiting evolved intermediate-mass stars from a precise Doppler survey of G and K giants at Okayama Astrophysical Observatory. All of the host stars are considered to be formerly early F-type or A-type dwarfs when they were on the main sequence. 14 And (K0 III) is a clump giant with a mass of 2.2 M_solar and has a planet of minimum mass m_2sin i=4.8 M_Jup in a nearly circular orbit with a 186 day period. This is one of the innermost planets around evolved intermediate-mass stars and such planets have only been discovered in clump giants. 81 Cet (G5 III) is a clump giant with 2.4 M_solar hosting a planet of m_2sin i=5.3 M_Jup in a 953 day orbit with an eccentricity of e=0.21. 6 Lyn (K0 IV) is a less evolved subgiant with 1.7 M_solar and has a planet of m_2sin i=2.4 M_Jup in a 899 day orbit with e=0.13. HD 167042 (K1 IV) is also a less evolved star with 1.5 M_solar hosting a planet of m_2sin i=1.6 M_Jup in a 418 day orbit with e=0.10. This planet was independently announced by Johnson et al. (2008, ApJ, 675, 784). All of the host stars have solar or sub-solar metallicity, which supports the lack of metal-rich tendency in planet-harboring giants in contrast to the case of dwarfs.

We address two outstanding issues in the sequential accretion scenario for gas giant planet formation, the retention of dust grains in the presence of gas drag and that of cores despite type I migration. The efficiency of these processes is determined by the disk structure. Theoretical models suggest that planets form in protostellar disk regions with an inactive neutral ``dead zone'' near the mid plane, sandwiched together by partially ionized surface layers where magnetorotational instability is active. Due to a transition in the abundance of dust grains, the active layer's thickness decreases abruptly near the ice line. Over a range of modest accretion rates ($\sim 10^{-9}-10^{-8} M_\odot$ yr$^{-1}$), the change in the angular momentum transfer rate leads to local surface density and pressure distribution maxima near the ice line. The azimuthal velocity becomes super-Keplerian and the grains accumulate in this transition zone. This barrier locally retains protoplanetary cores and enhances the heavy element surface density to the critical value needed to initiate efficient gas accretion. It leads to a preferred location and epoch of gas giant formation. We simulate and reproduce the observed frequency and mass-period distribution of gas giants around solar type stars without having to greatly reduce the type I migration strength. The mass function of the short-period planets can be utilized to calibrate the efficiency of type I migration and to extrapolate the fraction of stars with habitable terrestrial planets.

We have investigated the planetesimal accretion rate onto giant planets that are growing through gas accretion, using numerical simulations and analytical arguments. We derived the condition for gap opening in the planetesimal disk, which is determined by a competition between the expansion of the planet's Hill radius due to the planet growth and the damping of planetesimal eccentricity due to gas drag. We also derived the semi-analytical formula for the planetesimal accretion rate as a function of ratios of the rates of the Hill radius expansion, the damping, and planetesimal scattering by the planet. The predicted low planetesimal accretion rate due to gap opening in early gas accretion stages quantitatively shows that "phase 2," which is a long slow gas accretion phase before onset of runaway gas accretion, is not likely to occur. In late stages, rapid Hill radius expansion fills the gap, resulting in significant planetesimal accretion, which is as large as several $M_{\oplus}$ for Jupiter and Saturn. The efficient onset of runaway gas accretion and the late pollution may reconcile the ubiquity of extrasolar giant planets with metal-rich envelopes of Jupiter and Saturn inferred from interior structure models. These formulae will give deep insights into formation of extrasolar gas giants and the diversity in metallicity of transiting gas giants.

We report the detection of 3 new extrasolar planets from the precise Doppler survey of G and K giants at Okayama Astrophysical Observatory. The host stars, namely, 18 Del (G6 III), xi Aql (K0 III) and HD 81688 (K0 III-IV), are located at the clump region on the HR diagram with estimated masses of 2.1-2.3 M_solar. 18 Del b has a minimum mass of 10.3 M_Jup and resides in a nearly circular orbit with period of 993 days, which is the longest one ever discovered around evolved stars. xi Aql b and HD 81688 b have minimum masses of 2.8 and 2.7 M_Jup, and reside in nearly circular orbits with periods of 137 and 184 days, respectively, which are the shortest ones among planets around evolved stars. All of the substellar companions ever discovered around possible intermediate-mass (1.7-3.9 M_solar) clump giants have semimajor axes larger than 0.68 AU, suggesting the lack of short-period planets. Our numerical calculations suggest that Jupiter-mass planets within about 0.5 AU (even up to 1 AU depending on the metallicity and adopted models) around 2-3 M_solar stars could be engulfed by the central stars at the tip of RGB due to tidal torque from the central stars. Assuming that most of the clump giants are post-RGB stars, we can not distinguish whether the lack of short-period planets is primordial or due to engulfment by central stars. Deriving reliable mass and evolutionary status for evolved stars is highly required for further investigation of formation and evolution of planetary systems around intermediate-mass stars.

We have investigated the formation of close-in extrasolar giant planets through a coupling effect of mutual scattering, Kozai mechanism, and tidal circularization, by orbital integrations. We have carried out orbital integrations of three planets with Jupiter-mass, directly including the effect of tidal circularization. We have found that in about 30% runs close-in planets are formed, which is much higher than suggested by previous studies. We have found that Kozai mechanism by outer planets is responsible for the formation of close-in planets. During the three-planet orbital crossing, the Kozai excitation is repeated and the eccentricity is often increased secularly to values close enough to unity for tidal circularization to transform the inner planet to a close-in planet. Since a moderate eccentricity can remain for the close-in planet, this mechanism may account for the observed close-in planets with moderate eccentricities and without nearby secondary planets. Since these planets also remain a broad range of orbital inclinations (even retrograde ones), the contribution of this process would be clarified by more observations of Rossiter-McLaughlin effects for transiting planets.

In a further development of a deterministic planet-formation model (Ida & Lin 2004), we consider the effect of type-I migration of protoplanetary embryos due to their tidal interaction with their nascent disks. During the early embedded phase of protostellar disks, although embryos rapidly emerge in regions interior to the ice line, uninhibited type-I migration leads to their efficient self-clearing. But, embryos continue to form from residual planetesimals at increasingly large radii, repeatedly migrate inward, and provide a main channel of heavy element accretion onto their host stars. During the advanced stages of disk evolution (a few Myr), the gas surface density declines to values comparable to or smaller than that of the minimum mass nebula model and type-I migration is no longer an effective disruption mechanism for mars-mass embryos. Over wide ranges of initial disk surface densities and type-I migration efficiency, the surviving population of embryos interior to the ice line has a total mass several times that of the Earth. With this reservoir, there is an adequate inventory of residual embryos to subsequently assemble into rocky planets similar to those around the Sun. But, the onset of efficient gas accretion requires the emergence and retention of cores, more massive than a few M_earth, prior to the severe depletion of the disk gas. The formation probability of gas giant planets and hence the predicted mass and semimajor axis distributions of extrasolar gas giants are sensitively determined by the strength of type-I migration. We suggest that the observed fraction of solar-type stars with gas giant planets can be reproduced only if the actual type-I migration time scale is an order of magnitude longer than that deduced from linear theories.

We report the detection of a substellar companion orbiting the intermediate-mass giant star 11 Com (G8 III). Precise Doppler measurements of the star from Xinglong station and Okayama Astrophysical Observatory (OAO) revealed Keplerian velocity variations with an orbital period of 326.03 +/- 0.32 days, a semiamplitude of 302.8 +/- 2.6 m/s, and an eccentricity of 0.231 +/- 0.005. Adopting a stellar mass of 2.7 +/- 0.3 M_solar, the minimum mass of the companion is 19.4 +/- 1.5 M_Jup, well above the deuterium burning limit, and the semimajor axis is 1.29 +/- 0.05 AU. This is the first result from the joint planet search program between China and Japan aiming at revealing statistics of substellar companions around intermediate-mass giants. 11 Com b emerged from 300 targets of the planet search program at OAO. The current detection rate of a brown dwarf candidate seems to be comparable to that around solar-type stars within orbital separations of $\sim$3 AU.

We report the detection of five Jovian mass planets orbiting high metallicity stars. Four of these stars were first observed as part of the N2K program and exhibited low RMS velocity scatter after three consecutive observations. However, follow-up observations over the last three years now reveal the presence of longer period planets with orbital periods ranging from 21 days to a few years. HD 11506 is a G0V star with a planet of \msini = 4.74 \mjup in a 3.85 year orbit. HD 17156 is a G0V star with a 3.12 \mjup planet in a 21.2 day orbit. The eccentricity of this orbit is 0.67, one of the highest known for a planet with a relatively short period. The orbital period for this planet places it in a region of parameter space where relatively few planets have been detected. HD 125612 is a G3V star with a planet of \msini = 3.5 \mjup in a 1.4 year orbit. HD 170469 is a G5IV star with a planet of \msini = 0.67 \mjup in a 3.13 year orbit. HD 231701 is an F8V star with planet of 1.08 \mjup in a 142 day orbit. All of these stars have supersolar metallicity. Three of the five stars were observed photometrically but showed no evidence of brightness variability. A transit search conducted for HD 17156 was negative but covered only 25% of the search space and so is not conclusive.

We investigate both analytically and numerically the motion of massless particles orbiting primary star in a close circular binary system with particular focus on the gas drag effects. These are the first calculations with particles ranging in size from 1m to 10km, which account for the presence of a tidally perturbed gaseous disk. We have found numerically that the radial mass transport by the tidal waves plays a crucial role in the orbital evolution of particles. Numerical results are confirmed analytical calculations that do not assume anything about origin of the radial flow in the disk. We demonstrate that the migration rate of a particle in a disk out of radial equilibrium is enhanced due to the enhanced mass flux of gas colliding with the particle and the migration is always directed inward regardless of the sign of the radial gas velocity. Within the framework of the perturbation theory we derive general, formulae for short-term variations of the particle semi-major axis, eccentricity and inclination in such disk. The formulae account for departures from axial symmetry by introducing effective components of the gas velocity. They agree with numerical results within several percent. We have also found in numerical simulations that the tidal waves introduce coherence in periastron longitude and eccentricity for particles on neighbouring orbits. The degree of the coherence depends on the particle size and on the distance from the primary star, being most prominent for particles with 10m radius. The results are important mainly in the context of planetesimal formation and, to a lesser degree, during the early planetesimal accretion stage.

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).

We have investigated the obliquity evolution of terrestrial planets in habitable zones (at ~ 1AU) in extrasolar planetary systems, due to tidal interactions with their satellite and host star with wide varieties of satellite-to-planet mass ratio and initial obliquity, through numerical calculations and analytical arguments. The obliquity, the angle between planetary spin axis and its orbit normal, of a terrestrial planet is one of the key factors in determining the planetary surface environments. A recent scenario of terrestrial planet accretion implies that giant impacts of Mars-sized bodies determine the planetary spin and form satellites. With isotropic giant impacts, tilted spins are more likely to be produced than straight ones and satellites with various mass are formed. However, most of previous studies have focused on a particular case of the Earth-Moon systems or the two-body planar problem. We numerically integrated the evolution of planetary spin and a satellite orbit with various satellite mass and initial obliquity. We found that in the case of initially tiled spins, the satellite's orbit migrates outward until the orbit reaches the critical radius ~ 10-20 planetary radii, but then the migration is reversed to inward one with large-amplitude oscillation. The satellite eventually falls onto the planetary surface or it is captured at the synchronous state at several planetary radii. With the results of numerical integration and analytical arguments, we identified the parameter regions of qualitatively different evolution.

The evidence for a shortage of exosolar planets with semimajor axes -1.1 <= log (a/AU) <= -0.2 is investigated. It is shown that this valley results from a gap in the radial distribution of planets, orbiting stars with masses M >= 1.2 Msolar (the high-mass sample, HMS). No underabundance is found for planets orbiting stars with smaller masses. The observational data also indicate that within the HMS population it is preferentially the more massive planets with M sin(i) >= 0.8 M_J that are missing. Monte-Carlo simulations of planet formation and migration are presented that reproduce the observed shortage of planets in the observed radius regime. A dependence on the disk depletion timescale tau_dep is found. The gap is more pronounced for tau_dep = 10^6 - 10^7 yrs than for tau_dep = 3*10^6 - 3*10^7 yrs. This might explain the observed trend with stellar mass if disks around stars with masses M* >= 1.2 Msolar have shorter depletion timescales than those around less massive stars. Possible reasons for such a dependence are a decrease of disk size and an increase of stellar EUV flux with stellar mass.

The high density of the close-in extrasolar planet HD149026b suggests the presence of a huge core in the planet, which challenges planet formation theory. We first derive constraints on the amount of heavy elements and hydrogen/helium present in the planet: We find that preferred values of the core mass are between 50 and 80 M_E. We then investigate the possibility of subcritical core accretion as envisioned for Uranus and Neptune and find that the subcritical accretion scenario is unlikely in the case of HD149026b for at least two reasons: (i) Subcritical planets are such that the ratio of their core mass to their total mass is above ~0.7, in contradiction with constraints for all but the most extreme interior models of HD149026b; (ii) High accretion rates and large isolation mass required for the formation of a subcritical core of 30 M_E are possible only at specific orbital distances in a disk with a surface density of dust equal to at least 10 times that of the minimum mass solar nebula. This value climbs to 30 when considering a 50 M_E core. These facts point toward two main routes for the formation of this planet: (i) Gas accretion that is limited by a slow viscous inflow of gas in an evaporating disk; (ii) A significant modification of the composition of the planet after as accretion has stopped. These two routes are not mutually exclusive. Illustrating the second route, we show that for a wide range of impact parameters, giant impacts lead to a loss of the gas component of the planet and thus may lead to planets that are highly enriched in heavy elements. In the giant impact scenario, we expect an outer giant planet to be present. Observational studies by imaging, astrometry and long term interferometry of this system are needed to better narrow down the ensemble of possibilities.

Doppler measurements from Subaru and Keck have revealed radial velocity variations in the V=8.15, G0IV star HD 149026 consistent with a Saturn-Mass planet in a 2.8766 day orbit. Photometric observations at Fairborn Observatory have detected three complete transit events with depths of 0.003 mag at the predicted times of conjunction. HD 149026 is now the second brightest star with a transiting extrasolar planet. The mass of the star, based on interpolation of stellar evolutionary models, is 1.3 +/- 0.1 solar masses; together with the Doppler amplitude, K=43.3 m s^-1, we derive a planet mass Msin(i)=0.36 Mjup, and orbital radius of 0.042 AU. HD 149026 is chromospherically inactive and metal-rich with spectroscopically derived [Fe/H]=+0.36, Teff=6147 K, log g=4.26 and vsin(i)=6.0 km s^-1. Based on Teff and the stellar luminosity of 2.72 Lsun, we derive a stellar radius of 1.45 Rsun. Modeling of the three photometric transits provides an orbital inclination of 85.3 +/- 1.0 degrees and (including the uncertainty in the stellar radius) a planet radius of 0.725 +/- 0.05 Rjup. Models for this planet mass and radius suggest the presence of a ~67 Mearth core composed of elements heavier than hydrogen and helium. This substantial planet core would be difficult to construct by gravitational instability.

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' 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.

Dust growth and settling considerably affect the spectral energy distributions (SEDs) of protoplanetary disks. We investigated dust growth and settling in protoplanetary disks through numerical simulations to examine time-evolution of the disk optical thickness and SEDs. In the present paper, we considered laminar disks as a first step of a series of papers. As a result of dust growth and settling, a dust layer forms around the mid-plane of a gaseous disk. After the formation of the dust layer, small dust grains remain floating above the layer. Although the surface density of the floating small grains is much less than that of the dust layer, they govern the disk optical thickness and the emission. The floating small grains settle onto the dust layer in a long time scale compared with the formation of the dust layer. Rapid grain growth in the inner part of disks makes the radial distribution of the disk optical thickness less steep than that of the disk surface density. At t > 10^6yr, the optical thickness of the inner disk (> a few AU) almost vanishes, which may correspond to disk inner holes observed by Spitzer Space Telescope. Furthermore, we examined time-evolution of disk SEDs, using our numerical results and the two-layer model. The grain growth and settling decrease the magnitude of the SEDs especially. Our results indicate that grain growth and settling can explain the decrease in observed energy fluxes at millimeter/sub-millimeter wavelengths with time scales of 10^{6-7}yr without depletion of the disks.

The N2K consortium is carrying out a distributed observing campaign with the Keck, Magellan and Subaru telescopes, as well as the automatic photometric telescopes of Fairborn Observatory, in order to search for short-period gas giant planets around metal-rich stars. We have established a reservoir of more than 14,000 main sequence and subgiant stars, closer than 110 pc, brighter than V=10.5 and with $0.4 < B-V < 1.2$. Because the fraction of stars with planets is a sensitive function of stellar metallicity, a broadband photometric calibration has been developed to identify a subset of 2000 stars with [Fe/H] $> 0.1$ dex for this survey. We outline the strategy and report the detection of a planet orbiting the metal-rich G5IV star HD 88133 with a period of 3.41 days, semi-velocity amplitude, K=35.7 \mse and \msini = 0.29 \mjupe. Photometric observations reveal that HD 88133 is constant on the 3.415-day radial velocity period to a limit of 0.0005 mag. Despite a transit probability of 15.6%, our photometry rules out the shallow transits predicted by the large stellar radius.

The apparent dependence of detection frequency of extrasolar planets on the metallicity of their host stars is investigated with Monte Carlo simulations using a deterministic core-accretion planet formation model. According to this model, gas giants formed and acquired their mass $M_{\rm p}$ through planetesimal coagulation followed by the emergence of cores onto which gas is accreted. These protoplanets migrate and attain their asymptotic semi-major axis $a$ through their tidal interaction with their nascent disk. Based on the observed properties of protostellar disks, we generate $M_{\rm p}$-$a$ distribution. Our results reproduce the observed lack of planets with intermediate mass $M_{\rm p} = 10$--100$M_{\oplus}$ and $a \la 3$AU and with large mass $M_{\rm p} \ga 10^3 M_{\oplus}$ and $a \la 0.2$AU. Based on the simulated $M_{\rm p}$-$a$ distributions, we also evaluate the metallicity dependence of fraction of stars harboring planets that are detectable with current radial velocity survey. If protostellar disks attain the same fraction of heavy elements which are contained in their host stars, the detection probability around metal-rich stars would be greatly enhanced because protoplanetary cores formed in them can grow to several Earth masses prior to their depletion. These large masses are required for the cores to initiate rapid gas accretion and to transform into giant planets. The theoretically extrapolated metallicity dependence is consistent with the observation. This correlation does not arise naturally in the gravitational-instability scenario. We also suggest other metallicity dependence of the planet distributions that can be tested by on-going observations.

We examine the accretion of cores of giant planets from planetesimals, gas accretion onto the cores, and their orbital migration. We adopt a working model for nascent protostellar disks with a wide variety of surface density distributions in order to explore the range of diversity among extra solar planetary systems. If some cores can acquire more mass than a critical value of several Earth masses during the persistence of the disk gas, they would be able to rapidly accrete gas and evolve into gas giant planets. The gas accretion process is initially regulated by the Kelvin-Helmholtz contraction of the planets' gas envelope. Based on the assumption that the exponential decay of the disk-gas mass occurs on the time scales $\sim 10^{6}-10^{7}$ years and that the disk mass distribution is comparable to those inferred from the observations of circumstellar disks of T Tauri stars, we carry out simulations to predict the distributions of masses and semi major axes of extra solar planets. Since planets' masses grow rapidly from $10 M_{\oplus}$ to $100 M_{\oplus}$, the gas giant planets rarely form with asymptotic masses in this intermediate range. Our model predicts a paucity of extra solar planets with mass in the range 10-$100 M_{\oplus}$ and semi major axis less than 3AU. We refer to this deficit as a ``planet desert''. The effect of migration is to sharpen the boundaries and to enhance the contrast of the planet desert. The mass and semi major axis distributions generated in our simulations for the gas giants are consistent with those of the known extra solar planets.

Most extrasolar planets are observed to have eccentricities much larger than those in the solar system. Some of these planets have sibling planets, with comparable masses, orbiting around the same host stars. In these multiple planetary systems, eccentricity is modulated by the planets' mutual secular interaction as a consequence of angular momentum exchange between them. For mature planets, the eigenfrequencies of this modulation are determined by their mass and semi-major axis ratios. But, prior to the disk depletion, self gravity of the planets' nascent disks dominates the precession eigenfrequencies. We examine here the initial evolution of young planets' eccentricity due to the apsidal libration or circulation induced by both the secular interaction between them and the self gravity of their nascent disks. We show that as the latter effect declines adiabatically with disk depletion, the modulation amplitude of the planets' relative phase of periapse is approximately invariant despite the time-asymmetrical exchange of angular momentum between planets. However, as the young planets' orbits pass through a state of secular resonance, their mean eccentricities undergo systematic quantitative changes. For applications, we analyze the eccentricity evolution of planets around Upsilon Andromedae and HD168443 during the epoch of protostellar disk depletion. We find that the disk depletion can change the planets' eccentricity ratio. However, the relatively large amplitude of the planets' eccentricity cannot be excited if all the planets had small initial eccentricities.

We numerically calculated angular momentum transfer processes in a dense particulate disk within Roche limit by global $N$-body simulations, up to $N=10^5$, for parameters corresponding to a protolunar disk generated by a giant impact on a proto-Earth. In the simulations, both self-gravity and inelastic physical collisions are included. We first formalized expressions for angular momentum transfer rate including self-gravity and calculated the transfer rate with the results of our $N$-body simulations. Spiral structure is formed within Roche limit by self-gravity and energy dissipation of inelastic collisions, and angular momentum is effectively transfered outward. Angular momentum transfer is dominated by both gravitational torque due to the spiral structure and particles' collective motion associated with the structure. Since formation and evolution of the spiral structure is regulated by the disk surface density, the angular momentum transfer rate depends on surface density, but not on particle size or number, so that the time scale of evolution of a particulate disk is independent on the number of particles ($N$) that is used to represent the disk, if $N$ is large enough to represent the spiral structure. With $N=10^5$, the detailed spiral structure is resolved while it is only poorly resolved with $N=10^3$, however, we found that calculated angular momentum transfer does not change as long as $N \gtrsim 10^3$.

We investigate the effects of a passing stellar encounter on a planetesimal disk through analytical calculations and numerical simulations, and derive the boundary radius ($a_{\rm planet}$) outside which planet formation is inhibited by disruptive collisions with high relative velocities.

We investigate the formation of spatial structure in dense, self-gravitating particle systems such as Saturn'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{\lambda}{cr} = 4\pi^2 G\Sigma/\kappa^2$ in the linear theory of axisymmetric gravitational instability in a thin disk, where $G$, $\Sigma$, and $\kappa$ are the gravitational constant, surface density, and a epicyclic frequency.

We have investigated effects of early stellar encounters on a protoplanetary disk (planetesimal disk) and found that they can explain the high eccentricities and inclinations observed in the outer part ($>42$AU) of the Edgeworth-Kuiper Belt (EKB). The proto-sun is considered as a member of a stellar aggregation that undergoes dissolution on a timescale $\sim 10^8$yrs, such that a planetesimal disk experiences a flyby encounter at pericenter distance ($q$) on the order of 100AU. We simulated dynamical evolution of a planetesimal (test particle) disk perturbed by a passing star. We show that the stellar encounter pumps the velocity dispersion in the disk in the outer parts ($> 0.25q$). Planet formation is forestalled in that region. We also find that a stellar encounter with pericenter distance $q \sim 100-200$AU could have pumped up the velocity dispersion of EKB objects outside 42AU to the observed magnitude while preserving that inside Neptune's 3:2 mean-motion resonance (located at 39.5AU), which allows for the efficient capture of objects by the sweeping of the 3:2 resonance during orbital migration by proto-Neptune.

The ratio of the vertical velocity dispersion to radial one (sigma_z / sigma_R) of self-gravitating bodies in various disc potentials is investigated through two different numerical methods (statistical compilation of two-body encounters and N-body simulations). The velocity dispersion generated by two-body relaxation is considered. The ratio is given as a function of a disc potential parameter, kappa/Omega, where kappa and Omega are the epicycle and circular frequencies (the parameters kappa/Omega=1 and 2 correspond to Kepler rotation and solid-body rotation). For 1 <= kappa/Omega <~ 1.5, the velocity dispersion increases keeping some anisotropy (sigma_z / sigma_R ~ 0.5-0.7) if the amplitude of radial excursion is larger than tidal radius, while sigma_z / sigma_R << 1 for smaller amplitude. On the other hand, for 1.5 <~ kappa/Omega <= 2.0, we found isotropic state (sigma_z / sigma_R ~ 1) in the intermediate velocity regime, while anisotropic state (sigma_z / sigma_R < 1) still exists for higher and lower velocity regimes. The range of the intermediate velocity regime expands with kappa/Omega. In the limit of solid-body rotation, the regime covers all over the velocity space. Thus, the velocity dispersion generally has two different anisotropic states for each disc potential (1 <= kappa/Omega < 2) and one isotropic state for 1.5 <~ kappa/Omega < 2 where the individual states correspond to different amplitude of velocity dispersion, while in the limit of solid-body rotation (kappa/Omega = 2.0), entire velocity space is covered by the isotropic state.

We consider several processes operating during the late stages of planet formation that can affect observed orbital elements. Disk-planet interactions, tidal interactions with the central star, long term orbital instability and the Kozai mechanism are discussed.

Non-Gaussian velocity distribution in star forming region is reproduced by inelastic clump collision model. We numerically calculated the evolution of inelastic hard spheres in sheared flow, which corresponds to cloud clumps in differential galactic rotation. This system fluctuates largely around equilibrium state, creating clusters with inelastic collisions and destroying them with shear motion. The fluctuation makes spheres have non-Gaussian velocity distribution with nearly exponential tail. How far from Gaussian distribution depends upon coefficient of restitution, which can produce the variety of degree of deviation from Gaussian among regions.