We study the torque acting on a planet embedded in an optically thick accretion disc, using global two-dimensional hydrodynamic simulations. The temperature of an optically thick accretion disc is determined by the energy balance between the viscous heating and the radiative cooling. The radiative cooling rate depends on the opacity of the disc. The opacity is expressed as a function of the temperature. We find the disc is divided into three regions that have different temperature distributions. The slope of the entropy distribution becomes steep in the inner region of the disc with the high temperature and the outer region of the disc with the low temperature, while it becomes shallow in the middle region with the intermediate temperature. Planets in the inner and outer regions move outward owing to the large positive corotation torque exerted on the planet by an adiabatic disc, on the other hand, a planet in the middle region moves inward toward the central star. Planets are expected to accumulate at the boundary between the inner and middle regions of the adiabatic disc. The positive corotation torque decreases with an increase in the viscosity of the disc. We find that the positive corotation torque acting on the planet in the inner region becomes too small to cancel the negative Lindblad torque when we include the large viscosity, which destroys the enhancement of the density in the horseshoe orbit of the planet. This leads to the inward migration of the planet in the inner region of the disc. A planet with 5 Earth masses in the inner region can move outward in a disc with the surface density of 100 g/cm^2 at 1 AU when the accretion rate of a disc is smaller than 2x10^{-8} solar mass/yr.

We study Type I migration of a planet in a radiatively efficient disk using global two dimensional hydrodynamic simulations. The large positive corotation torque is exerted on a planet by an adiabatic disk at early times when the disk has the steep negative entropy gradient. The gas on the horseshoe orbit of the planet is compressed adiabatically during the change of the orbit from the slow orbit to the fast orbit, increasing its density and exerting the positive torque on the planet. The planet would migrate outward in the adiabatic disk before saturation sets in. We further study the effect of energy dissipation by radiation on Type I migration of the planet. The corotation torque decreases when the energy dissipates effectively because the density of the gas on the horseshoe orbit does not increase by the compression compared with the gas of the adiabatic disk. The total torque is mainly determined by the negative Lindblad torque and becomes negative. The planet migrates inward toward the central star in the radiatively efficient disk. The migration velocity is dependent on the radiative efficiency and greatly reduced if the radiative cooling works inefficiently.

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct $N$-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about $\sim 0.1$ Earth mass within $10^7$ years. This brings another difficulty in formation of gas giant planets, which requires cores with $\sim 10$ Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos' collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the solar system.