We develop estimators of agreement and disagreement between correlated cosmological data sets. These account for data correlations when computing the significance of both tensions and excess confirmation while remaining statistically optimal. We discuss and thoroughly characterize different approaches commenting on the ones that have the best behavior in practical applications. We complement the calculation of their statistical distribution within the Gaussian model with one estimator that takes non-Gaussianities fully into account. To illustrate the use of our techniques, we apply these estimators to supernovae measurements of the distance-redshift relation, absolutely calibrated by the local distance ladder. The suite of best estimators that we discuss finds results that are in excellent agreement between estimators and find no indications of significant internal inconsistencies in this data set above the $1\%$ probability threshold. This shows the robustness of local determinations of the Hubble constant to features in the distance-redshift relation.

The nonresonant cosmic ray instability, predicted by Bell (2004), is thought to play an important role in the acceleration and confinement of cosmic rays (CRs) close to supernova remnants. Despite its importance, the exact mechanism responsible for the saturation of the instability has not been determined, and there is no first-principle prediction for the amplitude of the saturated magnetic field. Using a survey of self-consistent kinetic hybrid simulations (with kinetic ions and fluid electrons), we study the saturation of the non-resonant streaming instability as a function of the parameters of both the thermal background plasma and the CR population. The strength of the saturated magnetic field has important implications for both CR acceleration in supernova remnants and CR diffusion in the Galaxy.

The limitation of the Quasilinear Theory (QLT) to describe the diffusion of electrons and ions in velocity space when interacting with a spectrum of large amplitude electrostatic Langmuir, Upper and Lower hybrid waves, is analyzed. We analytically and numerically estimate the threshold for the amplitude of the waves above which the QLT breaks down, using a test particle code. The evolution of the velocity distribution, the velocity-space diffusion coefficients, the driven current, and the heating of the particles are investigated, for the interaction with small and large amplitude electrostatic waves, i.e. in both regimes, there where QLT is valid and there where it clearly breaks down.