The formation of planets remains one of the most challenging problems of contemporary astrophysics. Starting with micron-sized dust grains, coagulation models predict growth up to centimeter (pebbles), but growth beyond this size is difficult because of fragmentation and drift. Ways to bypass this problem have focused on inhomogeneities in the flow, be that zonal flows, streaming instability, or vortices. Because vortices are in equilibrium between the Coriolis and the pressure force, the pressureless grains will orbit along a vortex streamline experiencing a drag force. This is a very effective mechanism to concentrate pebbles as also seen in numerical simulations and possibly in ALMA observations. Yet, a high pebble load is dangerous for the vortex, and we showed that in two-dimensional simulations the backreaction eventually leads to vortex disruption. We investigate whether the same happens in three dimensions. We perform 3D simulations with pebbles in a local box finding that, although the pebbles disturb the vortex around the midplane, the column does not get destroyed. This result is important because, based on the previous 2D result suggesting complete disruption, the vortex interpretation of ALMA observations has been called into question. We show instead that the vortex behaves like a Taylor column, and the pebbles as obstacles to the flow. Pebble accumulation in the center of the vortices proceeds to roughly the same concentration as in the control run without backreaction.

We analyse the concentration of solid particles in vortices created and sustained by radial buoyancy in protoplanetary disks, i.e. baroclinic vortex growth. Besides the gas drag acting on particles we also allow for back-reaction from dust onto the gas. This becomes important when the local dust-to-gas ratio approaches unity. In our 2D, local, shearing sheet simulations we see high concentrations of grains inside the vortices for a broad range of Stokes numbers, ${\rm St}$. An initial dust-to-gas ratio of 1:100 can easily be reversed to 100:1 for ${\rm St}=1$. The increased dust-to-gas ratio triggers the streaming instability, thus counter-intuitively limiting the maximal achievable overdensities. We find that particle trapping inside vortices opens the possibility for gravity-assisted planetesimal formation even for small particles ($\rm{St}=0.01$) and low initial dust-to-gas ratios (1:$10^4$)

Recent studies have shown that baroclinic vortex amplification is strongly dependent on certain factors, namely, the global entropy gradient, the efficiency of thermal diffusion and/or relaxation as well as numerical resolution. We conduct a comprehensive study of a broad range and combination of various entropy gradients, thermal diffusion and thermal relaxation time-scales via local shearing sheet simulations covering the parameter space relevant for protoplanetary disks. We measure the Reynolds stresses as a function of our control parameters and see that there is angular momentum transport even for entropy gradients as low as $\beta=-{d\ln s}/{d\ln r}={1}/{2}$, which corresponds to values observed in protoplanetary accretion disks. The amplification-rate of the perturbations, $\Gamma$, appears to be proportional to $\beta^2$ and thus proportional to the square of the \BV ($\Gamma \propto \beta^2 \propto N^2$). The saturation level of Reynolds stresses on the other hand seems to be proportional to $\beta^{1/2}$. This highlights the importance of baroclinic effects even for the low entropy gradients expected in protoplanetary disks.