We perform local $N$-body simulations of disk galaxies and investigate the evolution of spiral arms. We calculate the time autocorrelation of the surface density of spiral arms and find that the typical evolution timescale is described by the epicycle period. We investigate the distribution of the orbital elements of stars and find that in spiral arms the epicycle motions of stars are in phase while the spatial distribution of the guiding center is nearly uniform. These facts clearly show that the phase synchronization of the epicycle motion takes place, which is theoretically predicted by the swing amplification.

The origin of galactic spiral arms is one of fundamental problems in astrophysics. Based on the local analysis Toomre (1981) proposed the swing amplification mechanism in which the self-gravity forms spiral arms as leading waves of stars rotate to trailing ones due to galactic shear. The structure of spiral arms is characterized by their number and pitch angle. We perform global $N$-body simulations of spiral galaxies to investigate the dependence of the spiral structure on disk parameters and compare the simulation results with the swing amplification model. We find that the spiral structure in the $N$-body simulations agrees well with that predicted by the swing amplification for the wide range of parameters. The pitch angle decreases with increasing the shear rate and is independent of the disk mass fraction. The number of spiral arms decreases with both increasing the shear rate and the disk mass fraction. If the disk mass fraction is fixed, the pitch angle increases with the number of spiral arms.

Planetesimal formation is one of the most important unsolved problems in planet formation theory. In particular, rocky planetesimal formation is difficult because silicate dust grains are easily broken when they collide. Recently, it has been proposed that they can grow as porous aggregates when their monomer radius is smaller than $\sim$ 10 nm, which can also avoid the radial drift toward the central star. However, the stability of a layer composed of such porous silicate dust aggregates has not been investigated. Therefore, we investigate the gravitational instability of this dust layer. To evaluate the disk stability, we calculate Toomre's stability parameter $Q$, for which we need to evaluate the equilibrium random velocity of dust aggregates. We calculate the equilibrium random velocity considering gravitational scattering and collisions between dust aggregates, drag by mean flow of gas, stirring by gas turbulence, and gravitational scattering by gas density fluctuation due to turbulence. We derive the condition of the gravitational instability using the disk mass, dust-to-gas ratio, turbulent strength, orbital radius, and dust monomer radius. We find that, for the minimum mass solar nebula model at 1 au, the dust layer becomes gravitationally unstable when the turbulent strength $\alpha\lesssim10^{-5}$. If the dust-to-gas ratio is increased twice, the gravitational instability occurs for $\alpha\lesssim10^{-4}$. We also find that the dust layer is more unstable in disks with larger mass, higher dust-to-gas ratio, and weaker turbulent strength, at larger orbital radius, and with a larger monomer radius.

We consider the dynamics of porous icy dust aggregates in a turbulent gas disk and investigate the stability of the disk. We evaluate the random velocity of porous dust aggregates by considering their self-gravity, collisions, aerodynamic drag, turbulent stirring and scattering due to gas. We extend our previous work by introducing the anisotropic velocity dispersion and the relaxation time of the random velocity. We find the minimum mass solar nebular model to be gravitationally unstable if the turbulent viscosity parameter $\alpha$ is less than about $4 \times 10^{-3}$. The upper limit of $\alpha$ for the onset of gravitational instability is derived as a function of the disk parameters. We discuss the implications of the gravitational instability for planetesimal formation.

A ring system consisting of two dense narrow rings has been discovered around Centaur Chariklo. The existence of these rings around a small object poses various questions, such as their origin, stability, and lifetime. In order to understand the nature of Chariklo's rings, we perform global $N$-body simulations of the self-gravitating collisional particle rings for the first time. We find that Chariklo should be denser than the ring material to avoid the rapid diffusion of the rings. If Chariklo is denser than the ring material, fine spiral structures called self-gravity wakes occur in the inner ring. These wakes accelerate the viscous spreading of the ring significantly and they typically occur on timescales of about $100\,\mathrm{years}$ for m-sized ring particles, which is considerably shorter than the timescales suggested in previous studies. The existence of these narrow rings implies smaller ring particles or the existence of shepherding satellites.

Recently it is proposed that porous icy dust aggregates are formed by pairwise accretion of dust aggregates beyond the snowline. We calculate the equilibrium random velocity of porous dust aggregates taking into account mutual gravitational scattering, collisions, gas drag, and turbulent stirring and scattering. We find that the disk of porous dust aggregates becomes gravitationally unstable as they evolve through gravitational compression in the minimum-mass solar nebula model for a reasonable range of turbulence strength, which leads to rapid formation of planetesimals.

We revisit the swing amplification model of galactic spiral arms proposed by Toomre (1981). We describe the derivation of the perturbation equation in detail and investigate the amplification process of stellar spirals. We find that the elementary process of the swing amplification is the phase synchronization of the stellar epicycle motion. Regardless of the initial epicycle phase, the epicycle phases of stars in a spiral are synchronized during the amplification. Based on the phase synchronization, we explain the dependence of the pitch angle of spirals on the epicycle frequency. We find the most amplified spiral mode and calculate its pitch angle, wavelengths, and amplification factor, which are consistent with those obtained by the more rigorous model based on the Boltzmann equation by Julian and Toomre (1966).

Based on the swing amplification model of Julian and Toomre (1966), we investigate the formation and structure of stellar spirals in disk galaxies. We calculate the pitch angle, wavelengths, and amplification factor of the most amplified mode. We also obtain the fitting formulae of these quantities as a function of the epicycle frequency and Toomre's $Q$. As the epicycle frequency increases, the pitch angle and radial wavelength increases, while the azimuthal wavelength decreases. The pitch angle and radial wavelength increases with $Q$, while the azimuthal wavelength weakly depends on $Q$. The amplification factor decreases with $Q$ rapidly. In order to confirm the swing amplification model, we perform local $N$-body simulations. The wavelengths and pitch angle by the swing amplification model are in good agreement with those by $N$-body simulations. The dependence of the amplification factor on the epicycle frequency in $N$-body simulations is generally consistent with that in the swing amplification model. Using these results, we estimate the number of spiral arms as a function of the shear rate. The number of spiral arms increases with the shear rate if the disk to halo mass ratio is fixed.

We investigate the dynamics of self-gravity wakes in dense planetary rings. In particular, we examine how the pitch angle of self-gravity wakes depend on ring parameters using N-body simulations. We calculate the pitch angles using the two-dimensional autocorrelation function of the ring surface density. We obtain the pitch angles for the inner and outer parts of the autocorrelation function separately. We confirm that the pitch angles are 15 to 30 degrees for reasonable ring parameters, which are consistent with previous studies. We find that the inner pitch angle increases with the Saturnicentric distance, while it barely depends on the optical depth and the restitution coefficient of ring particles. The increase of the inner pitch angle with the Saturnicentric distance is consistent with the observations of the A ring. The outer pitch angle does not have the clear dependence on any ring parameters and is about 10 - 15 degrees. This value is consistent with the pitch angle of spiral arms in collisionless systems.

One of the key parameters that characterize 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.

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.

The Cassini spacecraft discovered a propeller-shaped structure in Saturn's A ring. This propeller structure is thought to be formed by gravitational scattering of ring particles by an unseen embedded moonlet. Self-gravity wakes are prevalent in dense rings due to gravitational instability. Strong gravitational wakes affect the propeller structure. Here, we derive the condition for formation of a propeller structure by a moonlet embedded in a dense ring with gravitational wakes. We find that a propeller structure is formed when the wavelength of the gravitational wakes is smaller than the Hill radius of the moonlet. We confirm this formation condition by performing numerical simulations. This condition is consistent with observations of propeller structures in Saturn's A ring.

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 non-axisymmetric 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_y^{3/2}$, 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.

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_\mathrm{aggr}/M_\mathrm{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.

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 $\rho_\mathrm{d}/\rho_\mathrm{g} \sim 200$ if the dust grains can grow up to 5m.

We consider the scattering of the 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 $\kappa$ along the line of sight to the source. In the short wavelength limit, the amplitude is magnified by $1+\kappa$, which is consistent with the result in weak gravitational lensing.

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.