The growth process of proto-planets can be sped-up by accreting a large number of solid, pebble-sized objects that are still present in the protoplanetary disc. It is still an open question on how efficient this process works in realistic turbulent discs. Here, we investigate the accretion of pebbles in turbulent discs that are driven by the purely hydrodynamical vertical shear instability (VSI). For this purpose, we perform global three-dimensional simulations of locally isothermal, VSI turbulent discs with embedded protoplanetary cores from 5 to 100 $M_\oplus$ that are placed at 5.2 au distance from the star. In addition, we follow the evolution of a swarm of embedded pebbles of different size under the action of drag forces between gas and particles in this turbulent flow. Simultaneously, we perform a set of comparison simulations for laminar viscous discs where the particles experience stochastic kicks. For both cases, we measure the accretion rate onto the cores as a function of core mass and Stokes number ($\tau_s$) of the particles and compare it to recent MRI turbulence simulations. Overall the dynamic is very similar for the particles in the VSI turbulent disc and the laminar case with stochastic kicks. For the small mass planets (i.e. 5 and 10 $M_\oplus$), well-coupled particles with $\tau_s = 1$, which have a size of about one meter at this location, we find an accretion efficiency (rate of particles accreted over drifting inward) of about 1.6-3%. For smaller and larger particles this efficiency is higher. However, the fast inward drift for $\tau_s = 1$ particles makes them the most effective for rapid growth, leading to mass doubling times of about 20,000 yr. For masses between 10 and 30 $M_\oplus$ the core reaches the pebble isolation mass and the particles are trapped at the pressure maximum just outside of the planet, shutting off further particle accretion.

In weakly ionized discs turbulence can be generated through the vertical shear instability (VSI). Embedded planets feel a stochastic component in the torques acting on them which can impact their migration. In this work we study the interplay between a growing planet embedded in a protoplanetary disc and the VSI-turbulence. We performed a series of three-dimensional hydrodynamical simulations for locally isothermal discs with embedded planets in the mass range from 5 to 100 Earth masses. We study planets embedded in an inviscid disc that is VSI unstable, becomes turbulent and generates angular momentum transport with an effective $\alpha = 5 \cdot 10^{-4}$. This is compared to the corresponding viscous disc using exactly this $\alpha$-value. In general we find that the planets have only a weak impact on the disc turbulence. Only for the largest planet ($100 M_\oplus$) the turbulent activity becomes enhanced inside of the planet. The depth and width of a gap created by the more massive planets ($30, 100 M_\oplus$) in the turbulent disc equal exactly that of the corresponding viscous case, leading to very similar torque strengths acting on the planet, with small stochastic fluctuations for the VSI disc. At the gap edges vortices are generated that are stronger and longer lived in the VSI disc. Low mass planets (with $M_p \leq 10 M_\oplus$) do not open gaps in the disc in both cases but generate for the turbulent disc an over-density behind the planet that exerts a significant negative torque. This can boost the inward migration in VSI turbulent discs well above the Type I rate. Due to the finite turbulence level in realistic three-dimensional discs the gap depth will always be limited and migration will not stall in inviscid discs.

Recently, the vertical shear instability (VSI) has become an attractive purely hydrodynamic candidate for the anomalous angular momentum transport required for weakly ionized accretion disks. In direct three-dimensional numerical simulations of VSI turbulence in disks, a meridional circulation pattern was observed that is opposite to the usual viscous flow behavior. Here, we investigate whether this feature can possibly be explained by an anisotropy of the VSI turbulence. Using three-dimensional hydrodynamical simulations, we calculate the turbulent Reynolds stresses relevant for angular momentum transport for a representative section of a disk. We find that the vertical stress is significantly stronger than the radial stress. Using our results in viscous disk simulations with different viscosity coefficients for the radial and vertical direction, we find good agreement with the VSI turbulence for the stresses and meridional flow; this provides additional evidence for the anisotropy. The results are important with respect to the transport of small embedded particles in disks.

Among the candidates for generating turbulence in accretion discs in situations with low intrinsic ionization the vertical shear instability (VSI) has become an interesting candidate, as it relies purely on a vertical gradient in the angular velocity. Existing simulations have shown that $\alpha$-values a few times $10^{-4}$ can be generated. The particle growth in the early planet formation phase is determined by the dynamics of dust particles. Here, we address in particular the efficiency of VSI-turbulence in concentrating particles in order to generate overdensities and low collision velocities. We perform 3D numerical hydrodynamical simulations of accretion discs around young stars that include radiative transport and irradiation from the central star. The motion of particles within a size range of a fraction of mm up to several m is followed using standard drag formula. We confirm that under realistic conditions the VSI is able to generate turbulence in full 3D protoplanetary discs. The irradiated disc shows turbulence within 10 to 60au. The mean radial motion of the gas is such that it is directed inward near the midplane and outward in the surface layers. We find that large particles drift inward with the expected speed, while small particles can experience phases of outward drift. Additionally, the particles show bunching behaviour with overdensities reaching 5 times the average value, which is strongest for dimensionless stopping times around unity. Particles in a VSI-turbulent discs are concentrated in large scale turbulent eddies and show low relative speeds that allow for growing collisions. The reached overdensities will also allow for the onset streaming instabilities further enhancing particle growth. The outward drift for small particles at higher disk elevations allows for the transport of processed high temperature material in the Solar System to larger distances.

The origin of turbulence in accretion discs is still not fully understood. While the magneto-rotational instability is considered to operate in sufficiently ionized discs, its role in the poorly ionized protoplanetary disc is questionable. Recently, the vertical shear instability (VSI) has been suggested as a possible alternative. Our goal is to study the characteristics of this instability and the efficiency of angular momentum transport, in extended discs, under the influence of radiative transport and irradiation from the central star. We use multi-dimensional hydrodynamic simulations to model a larger section of an accretion disc. First we study inviscid and weakly viscous discs using a fixed radial temperature profile in two and three spatial dimensions. The simulations are then extended to include radiative transport and irradiation from the central star. In agreement with previous studies we find for the isothermal disc a sustained unstable state with a weak positive angular momentum transport of the order of $\alpha \approx 10^{-4}$. Under the inclusion of radiative transport the disc cools off and the turbulence terminates. For discs irradiated from the central star we find again a persistent instability with a similar $\alpha$ value as for the isothermal case. We find that the VSI can indeed generate sustained turbulence in discs albeit at a relatively low level with $\alpha$ about few times $10^{-4}$