Is the ΛCDM model in crisis?

We present strong evidences for dynamical dark energy challenging the ΛCDM model. Several dark energy models are explored, including ω₀ω_aCDM, logarithmic, Exponential, JBP, and BA, along with non-flat cosmological models accounting also for potential spatial curvature different from zero. Through our analysis, we find evidences supporting a flat Universe (Ω_k≈0). Using the Metropolis-Hastings Markov Chain Monte Carlo algorithm, we analyze observational data from Baryon Acoustic Oscillations of DESI DR2, Type Ia Supernovae, and Compressed CMB likelihood to constrain the parameters of these models. Our findings provide strong evidences that ω ≠ −1, with deviations from the ΛCDM model favoring dynamical dark energy models characterized by the Quintom-B scenario (ω₀ > −1, ω_a <0, and ω₀ + ω_a <−1). We also derive the upper bounds on ∑m_ν using the combination of CMB and DESI DR2 data. For the ΛCDM model, we find ∑m_ν < 0.066eV, while for ωCDM, it is ∑m_ν<0.075eV. In the oΛCDM and o ω CDM models, the limits are ∑m_ν < 0.263eV and ∑m_ν < 0.520 eV, respectively. For other models, including ω₀ω_aCDM, Logarithmic, Exponential, JBP, BA, and GEDE, the upper limits range from <0.043 eV to < 0.127 eV, depending on the model. Constraints on the effective number of relativistic species, N_eff, show that our results remain consistent with the standard value of N_eff = 3.044 for each dark energy model. Bayesian evidence shows that combining the DES-SN5Y and Union3 SNe Ia samples with CMB + DESI DR2 reveals a deviation from the ΛCDM model. Finally, we found that none of the models reach the 5 σ threshold of deviation from ΛCDM, but some models show tensions exceeding 3 σ with DES-SN5YR or Union3, indicating we are starting to see the cracks in the cosmological constant Λ.