Temporal evolution of the gas surface density for our main simulation (solid curves). The dashed curves are for our additional simulation with lower turbulent viscosity. Initially, the distribution has a power-law index of −3/2, but the index changes with time in the inner region due to the influence of disk winds. As in previous studies, the inner disk edge is placed at r = 0.015 au. The start time of N-body simulations is delayed by t₀ from t_disk.