Coagulation in the Vicinity of a Gap-Opening Jupiter-Mass Planet

We analyze the coagulation of dust in and around a gap opened by a Jupiter-mass planet. To this end, we carry out a high-resolution magnetohydrodynamic (MHD) simulation of the gap environment, which is turbulent due to the magnetorotational instability. From the MHD simulation, we obtain values of the gas velocities, densities, and turbulent stresses (a) close to the gap edge, (b) in one of the two gas streams that accrete onto the planet, (c) inside the low-density gap, and (d) outside the gap. The MHD values are then input into a Monte Carlo dust-coagulation algorithm which models grain sticking and compaction. We also introduce a simple implementation for bouncing, for comparison purposes. We consider two dust populations for each region: one whose initial size distribution is monodisperse, with monomer radius equal to 1 μm, and another one whose initial size distribution follows the Mathis–Rumpl–Nordsieck distribution for interstellar dust grains, with an initial range of monomer radii between 0.5 and 10 μm. Without bouncing, our Monte Carlo calculations show steady growth of dust aggregates in all regions, and the mass-weighted (m-w) average porosity of the initially monodisperse population reaches xtremely high final values of 98%. The final m-w porosities in all other cases without bouncing range between 30% and 82%. The efficiency of compaction is due to high turbulent relative speeds between dust particles. When bouncing is introduced, growth is slowed down in the planetary wake and inside the gap. Future studies will need to explore the effect of different planet masses and electric charge on grains.