The ground state structural and energetic properties for rocksalt and cesium chloride phases of the cesium halides were explored using the Random Phase Approximation (RPA) and beyond-RPA methods to benchmark the non-empirical SCAN meta-GGA and its empirical dispersion corrections. The importance of non-additivity and higher-order multipole moments of dispersion in these systems is discussed. RPA generally predicts the equilibrium volume for these halides within 2.4% of the experimental value, while beyond-RPA methods utilizing the renormalized adiabatic LDA (rALDA) exchange-correlation kernel are typically within 1.8%. The zero-point vibrational energy is small and shows that the stability of these halides is purely due to electronic correlation effects. The rAPBE kernel as a correction to RPA overestimates the equilibrium volume and could not predict the correct phase ordering in the case of cesium chloride, while the rALDA kernel consistently predicted results in agreement with the experiment for all of the halides. However, due to its reasonable accuracy with lower computational cost, SCAN+rVV10 proved to be a good alternative to the RPA-like methods for describing the properties of these ionic solids.Alkali halides provide a useful benchmark for new theoretical methods to test their performance in predicting equilibrium and non-equilibrium properties of ionic solids 1-7 . Among the alkali halides, cesium halides are of particular interest in terms of their phase stability and have been studied both experimentally as well as theoretically 3,8-12 . CsF is experimentally stable in the B1 structure, while the Cl, Br, and I materials exist experimentally in the B2 structure. In Strukturbericht notation, B1 corresponds to the rocksalt (NaCl) phase, whereas B2 refers to the CsCl phase 13 . This difference in phase preference for the cesium halides can only be understood through the inclusion of dispersion interactions 3,8-12 .The unexpected stability of an ionic B2 phase was explained by London 8 through the presence of relatively large van der Waals interactions between the heavy Cs + cation and the heavier halide anions (Cl − , Br − , and I − ). Since dispersion effects are proportional to the polarizability and number of electrons in the anion, they are expected to become more important as one moves down the halide column. Furthermore, the coordination number of Cs in the B2 phase is higher than that in the B1 phase so there are locally more halide anions with which to interact. Dispersion is a pure quantum mechanical effect due to instantaneous or induced electronic multipole moments and is therefore difficult to capture with classical models 14 . The simplest treatment of the dispersion interaction is modeled by simple pairwiseadditive interactions between atoms, but this type of approximation completely ignores any nonadditive, nonlocal, and collective many-body effects 14,15 , which can be important in cases where screening effects modify electron-electron interactions.Rather than rely on classical models, ab...