Gas flow structure averaged in the y direction within the calculation domain of hydrodynamical simulation around a planet (panel a), the difference of the azimuthally and vertically averaged velocity of dust particles from the unperturbed steady-state drift velocity, (panel b), advective (solid lines) and diffusive (dashed lines) flux of dust with St = 10^‒3 in the steady state (panel c), and the dust surface density (panel d). Panels a, b, and d: different colors correspond to different Stokes numbers. Panel a: result obtained from m81-hw881. The solid gray lines correspond to the streamlines. The dashed circle is the Hill radius. The regions enclosed by blue curves correspond to the forbidden regions. The horizontal dashed lines in panels b and d correspond to and Σ_d/Σ_d,₀ = 1. The figure displayed at the lower right of panel b shows . Panels b-d: we adopt α_diff = 10^‒5. To improve legibility, we only plot the region where |x| ≤ 2.