We calculate the skewness (the third moment $S_3$) of matter distribution in dynamical dark energy cosmologies. We pay particular attention to the impact of dark energy perturbations on this quantity. There is indeed a clear signature of dark energy perturbations on this quantity. By properly allowing dark energy perturbations we show that their impact on $S_3$ is strong enough (a factor $\sim 3$ greater) to easily discriminate between clustering and non-clustering dark energy cosmologies. This indicates that high order statistics of the cosmic density field are useful to the study of dark energy models and are potentially able to rule out clustering dark energy cosmologies.

As recently pointed out in Ref. [Phys.\ Rev.\ D {\bf 96}, 8, 083502 (2017)] the evolution of the linear matter perturbations in nonadiabatic dynamical dark energy models is almost indistinguishable (quasi-degenerated) to the standard $\Lambda$CDM scenario. In this work we extend this analysis to CMB observables in particular the integrated Sachs-Wolfe effect and its cross-correlation with large scale structure. We find that this feature persists for such CMB related observable reinforcing that new probes and analysis are necessary to reveal the nonadiabatic features in the dark energy sector.

We study structure formation using relativistic cosmological linear perturbation theory in the presence of intrinsic and relative (with respect to matter) non-adiabatic dark energy perturbations. For different dark energy models we assess the impact of non-adiabaticity on the matter growth promoting a comparison with growth rate data. The dark energy models studied lead to peculiar signatures of the (non)adiabatic nature of dark energy perturbations in the evolution of the $f \sigma_{8}(z)$ observable. We show that non-adiabatic DE models become close to be degenerated with respect to the $\Lambda$CDM model at first order in linear perturbations. This would avoid the identification of the non-adiabatic nature of dark energy using current available data. Therefore, such evidence indicates that new probes are necessary to reveal the non-adiabatic features in the dark energy sector.

The concordance particle creation model - a class of $\Lambda(t)$CDM cosmologies - is studied using large scale structure (LSS) formation, with particular attention to the integrated Sachs-Wolfe (ISW) effect. The evolution of the gravitational potential and the amplitude of the cross-correlation of the cosmic microwave background (CMB) signal with LSS surveys are calculated in detail. We properly include in our analysis the peculiarities involving the baryonic dynamics of the $\Lambda(t)$CDM model which were not included in previous works. Although both the $\Lambda(t)$CDM and the standard cosmology are in agreement with available data for the CMB-LSS correlation, the former presents a slightly higher signal which can be identified with future data.