Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^\circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$\mu$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $\alpha \simeq 10^{-4}$. The value of $\alpha$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $\alpha$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $\sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation.

Several protoplanetary disks observed by ALMA show dust concentrations consistent with particle trapping in giant vortices. The formation and survival of vortices is of major importance for planet formation, because vortices act as particle traps and are therefore preferred locations of planetesimal formation. Recent studies showed that the vertical shear instability (VSI) is capable of generating turbulence and small vortices in protoplanetary disks that have the proper radial and vertical stratification and thermally relax on sufficiently short time scales. But the effect of the azimuthal extend of the disk is often neglected as the disks azimuth is limited to $\Delta \phi \leq \pi/2$. We aim to investigate the influence of the azimuthal extent of the disk on the long-term evolution of a protoplanetary disk and the possibility of large vortices forming. To this end, we perform 3-dimensional simulations for up to 1000 local orbits using different values of $\Delta \phi = \pi/2 $ to $2\pi$ for VSI in disks with a prescribed radial density and temperature gradient cooling on short timescales. We find the VSI capable of forming large vortices which can exist at least several hundred orbits in simulations covering a disk with $\Delta \phi \geq \pi$. This suggests the VSI to be capable to form vortices or at least to trigger vortex formation via a secondary instability, e.g. Rossby Wave Instability or Kelvin Helmholtz Instability.