Circumplanetary particle disks would be created in the late stage of planetary formation either by impacts of planetary bodies or disruption of satellites or passing bodies, and satellites can be formed by accretion of disk particles spreading across the Roche limit. Previous N-body simulation of lunar accretion focused on the formation of single-satellite systems from disks with large disk-to-planet mass ratios, while recent models of the formation of multiple-satellite systems from disks with smaller mass ratios do not take account of gravitational interaction between formed satellites. In the present work, we investigate satellite accretion from particle disks with various masses, using N-body simulation. In the case of accretion from somewhat less massive disks than the case of lunar accretion, formed satellites are not massive enough to clear out the disk, but can become massive enough to gravitationally shepherd the disk outer edge and start outward migration due to gravitational interaction with the disk. When the radial location of the 2:1 mean motion resonance of the satellite reaches outside the Roche limit, the second satellite can be formed near the disk outer edge, and then the two satellites continue outward migration while being locked in the resonance. Co-orbital satellites are found to be occasionally formed on the orbit of the first satellite. Our simulations also show that stochastic nature involved in gravitational interaction and collision between aggregates in the tidal environment can lead to diversity in the final mass and orbital architecture, which would be expected in satellite systems of exoplanets.

A compact gas cloud G2 is predicted to reach the pericenter of its orbit around the super massive black hole (SMBH) of our galaxy, Sagittarius A* (Sgr A*). This event will give us a rare opportunity to observe the interaction between SMBH and gas around it. We report the result of the fully three-dimensional simulation of the evolution of G2 during the first pericenter passage. The strong tidal force by the SMBH stretches the cloud along its orbit, and compresses it strongly in the vertical direction, resulting in the heating up and flaring up of the cloud. The bolometric luminosity will reach the maximum of $\sim100 L_{\odot}$. This flare should be easily observed in the near infrared.

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