One of the key parameters that characterizes spiral arms in disk galaxies is a pitch angle that measures the inclination of a spiral arm to the direction of galactic rotation. The pitch angle differs from galaxy to galaxy, which suggests that the rotation law of galactic disks determines it. In order to investigate the relation between the pitch angle of spiral arms and the shear rate of galactic differential rotation, we perform local N-body simulations of pure stellar disks. We find that the pitch angle increases with the epicycle frequency and decreases with the shear rate and obtain the fitting formula. This dependence is explained by the swing amplification mechanism.

Abstract.We consider the scattering of gravitational waves by the weak gravitational fields of lens objects. We obtain the scattered gravitational waveform by treating the gravitational potential of the lens to first order, i.e. using the Born approximation. We find that the effect of scattering on the waveform is roughly given by the Schwarzschild radius of the lens divided by the wavelength of gravitational wave for a compact lens object. If the lenses are smoothly distributed, the effect of scattering is of the order of the convergence field κ along the line of sight to the source. In the short wavelength limit, the amplitude is magnified by 1 + κ, which is consistent with the result in weak gravitational lensing.

We analyze the stability of dust layer in protoplanetary disk to understand the effect of the relative motion between gas and dust. The previous analyses not including the effect of relative motion between gas and dust show that the shear-induced turbulence may prevent the dust grains from settling sufficiently to be gravitationally unstable. We determine the growth rate of Kelvin-Helmholtz instability in wide range of parameter space, and propose a possible path toward the planetesimal formation through the gravitational instability. We expect the density of dust layer becomes ρ d /ρ g ∼ 100 if the dust grains can grow up to 10m.

We perform a linear stability analysis of a dust layer in a turbulent gas disk. Youdin (2011) investigated the secular gravitational instability of a dust layer using hydrodynamic equations with a turbulent diffusion term. We obtain essentially the same result independently of Youdin (2011). In the present analysis, we restrict the area of interest to small dust particles, while investigating the secular gravitational instability in a more rigorous manner. We discuss the time evolution of the dust surface density distribution using a stochastic model and derive the advection-diffusion equation. The validity of the analysis by Youdin (2011) is confirmed in the strong drag limit. We demonstrate quantitatively that the finite thickness of a dust layer weakens the secular gravitational instability and that the density-dependent diffusion coefficient changes the growth rate. We apply the obtained results to the turbulence driven by the shear instability and find that the secular gravitational instability is faster than the radial drift when the gas density is three times as large as that in the minimum-mass disk model. If the dust particles are larger than chondrules, the secular gravitational instability grows within the lifetime of a protoplanetary disk.

We performed N-body simulations of a dust layer without a gas component and examined the formation process of planetesimals. We found that the formation process of planetesimals can be divided into three stages: the formation of non-axisymmetric wake-like structures, the creation of aggregates, and the collisional growth of the aggregates. Finally, a few large aggregates and many small aggregates are formed. The mass of the largest aggregate is larger than the mass predicted by the linear perturbation theory. We examined the dependence of system parameters on the planetesimal formation. We found that the mass of the largest aggregates increase as the size of the computational domain increases. However the ratio of the aggregate mass to the total mass M aggr /M total is almost constant 0.8 − 0.9. The mass of the largest aggregate increases with the optical depth and the Hill radius of particles.

Wave effects can be important for the gravitational lensing of gravitational waves. In such a case, wave optics must be used in stead of geometric optics. We consider a plane wave entering a lens object and solve numerically the wave equation for three lens models: the uniform density sphere, the singular isothermal sphere, and the Hernquist model. By comparing our numerical solutions with the analytical solutions under the thin lens approximation, we evaluate the error of this approximation. The results show that the relative error of the thin lens approximation is small if the geometrical thickness of the lens is much smaller than the distance between the lens and the observer.

We investigate the formation process of planetesimals from the dust layer by the gravitational instability in the gas disk using local N -body simulations. The gas is modeled as a background laminar flow. We study the formation process of planetesimals and its dependence on the strength of the gas drag. Our simulation results show that the formation process is divided into three stages qualitatively: the formation of wake-like density structures, the creation of planetesimal seeds, and their collisional growth. The linear analysis of the dissipative gravitational instability shows that the dust layer is secularly unstable although Toomre's Q value is larger than unity. However, in the initial stage, the growth time of the gravitational instability is longer than that of the dust sedimentation and the decrease in the velocity dispersion. Thus, the velocity dispersion decreases and the disk shrinks vertically. As the velocity dispersion becomes sufficiently small, the gravitational instability finally becomes dominant. Then wake-like density structures are formed by the gravitational instability. These structures fragment into planetesimal seeds. The seeds grow rapidly owing to mutual collisions.

The gravitational instability of a dust layer is one of the scenarios for planetesimal formation. If the density of a dust layer becomes sufficiently high as a result of the sedimentation of dust grains toward the midplane of a protoplanetary disk, the layer becomes gravitationally unstable and spontaneously fragments into planetesimals. Using a shearing box method, we performed local N-body simulations of gravitational instability of a dust layer and subsequent coagulation without gas and investigated the basic formation process of planetesimals. In this paper, we adopted the accretion model as a collision model. A gravitationally bound pair of particles is replaced by a single particle with the total mass of the pair. This accretion model enables us to perform long-term and large-scale calculations. We confirmed that the formation process of planetesimals is the same as that in the previous paper with the rubble pile models. The formation process is divided into three stages: the formation of nonaxisymmetric structures; the creation of planetesimal seeds; and their collisional growth. We investigated the dependence of the planetesimal mass on the simulation domain size. We found that the mean mass of planetesimals formed in simulations is proportional to L 3/2 y , where L y is the size of the computational domain in the direction of rotation. However, the mean mass of planetesimals is independent of L x , where L x is the size of the computational domain in the radial direction if L x is sufficiently large. We presented the estimation formula of the planetesimal mass taking into account the simulation domain size.