The pebble accretion model has the potential to explain the formation of various types of planets. A growing planet embedded in a disk induces three-dimensional (3D) gas flow, which may influence pebble accretion. In this study, we investigate the influence of the 3D planet-induced gas flow on pebble accretion. Assuming a non-isothermal, inviscid gas disk, we perform 3D hydrodynamical simulations on the spherical polar grid. Then we numerically integrate the equation of motion of pebbles in 3D using hydrodynamical simulations data. We find that the trajectories of pebbles in the planet-induced gas flow differ significantly from those in the unperturbed shear flow for a wide range of pebble sizes investigated (${\rm St}=10^{-3}$--$10^{0}$, where ${\rm St}$ is the Stokes number). The horseshoe flow and outflow of the gas alter the motion of the pebbles, which leads to the reduction of the width of the accretion window, $w_{\rm acc}$, and the accretion cross section, $A_{\rm acc}$. On the other hand, the changes in trajectories also cause an increase in relative velocity of pebbles to the planet, which offsets the reduction of $w_{\rm acc}$ and $A_{\rm acc}$. As a consequence, in the Stokes regime, the accretion probability of pebbles, $P_{\rm acc}$, in the planet-induced gas flow is comparable to that in the unperturbed shear flow except when the Stokes number is small, ${\rm St}\sim10^{-3}$, in 2D accretion, or when the thermal mass of the planet is small, $m=0.03$ in 3D accretion. In contrast, in the Epstein regime, $P_{\rm acc}$ in the planet-induced gas flow becomes smaller than that in the shear flow in the Stokes regime in both 2D and 3D accretion, regardless of assumed ${\rm St}$ and $m$. Our results suggest that the 3D planet-induced gas flow may be helpful to explain the distribution of exoplanets as well as the architecture of the solar system.

The three-dimensional structure of the gas flow around a planet is thought to influence the accretion of both gas and solid materials. In particular, the outflow in the mid-plane region may prevent the accretion of the solid materials and delay the formation of super-Earths' cores. However, it is not yet understood how the nature of the flow field and outflow speed change as a function of the planetary mass. In this study, we investigate the dependence of gas flow around a planet embedded in a protoplanetary disc on the planetary mass. Assuming an isothermal, inviscid gas disc, we perform three-dimensional hydrodynamical simulations on the spherical polar grid, which has a planet located at its centre. We find that gas enters the Bondi or Hill sphere at high latitudes and exits through the mid-plane region of the disc regardless of the assumed dimensionless planetary mass $m=R_{\rm Bondi}/H$, where $R_{\rm Bondi}$ and $H$ are the Bondi radius of the planet and disc scale height, respectively. The altitude from where gas predominantly enters the envelope varies with the planetary mass. The outflow speed can be expressed as $|u_{\rm out}|=\sqrt{3/2}mc_{\rm s}$ $(R_{\rm Bondi}\leq R_{\rm Hill})$ or $|u_{\rm out}|=\sqrt{3/2}(m/3)^{1/3} c_{\rm s}$ ($R_{\rm Bondi}\geq R_{\rm Hill}$), where $c_{\rm s}$ is the isothermal sound speed and $R_{\rm Hill}$ is the Hill radius. The outflow around a planet may reduce the accretion of dust and pebbles onto the planet when $m\gtrsim\sqrt{\rm St}$, where St is the Stokes number. Our results suggest that the flow around proto-cores of super-Earths may delay their growth and, consequently, help them to avoid runaway gas accretion within the lifetime of the gas disc.