where symmetric breakup n x = n y = n/2 corresponds to ξ = 1/2, and maximally asymmetric breakup n x = n -1, n y = 1 is denoted by ξ = 1/n. The results of Equation 32 with parameters D, ∆χ, and ∆η computed from parametrized Equations 26 and 28 are also included in Figure . If we compare our results with the experimental data, we see very similar trends of increasing stability of the poly-charged helium clusters with the increase of cluster size. It is worth noting that, unlike the polyenes and annulenes, helium clusters are more stable for the symmetric breakup (ξ = 1/2) than for the small cluster removal (ξ → 0). This is in agreement with the experimental observations, where the strong evaporation of charged helium clusters was observed via the removal of small charged fragments and not through symmetric Coulomb explosions. Despite the qualitative agreement of our theoretically predicted stability trends with the experimental observations, the quantitative predictions differ from the actual observed trends. The main reason for this is that our theoretical methodology has a limited quality of the quantum-chemical calculations, in particular for the sampling of the conformational space, poor harmonic approximation of the ZPVE for the superfluid vibrational state, limited size of the basis set for calculation of the electronic energies, and significant failures of the conventional basis sets to describe the EA of helium (see SI for details). These defects of the theoretical model, combined with extrapolation of the parameters from cluster size ranges of 1 ≤ n < 1000 up to n ∼ 10 8 , can lead to significant deviations in the observed trends.