To test our hypothesis and analyse the correlations of the node disorder experimentally, we turn to diffuse scattering (DS). A full description of the data collection, treatment, and 3D-∆PDF extraction can be found in the S.I. We observe significant DS in the 0kl and h0l layers, structured into sharp planes perpendicular to l [Fig. ]. The linewidth of these planes is resolution limited (FWHM < 0.03 r.l.u.), implying a correlation length along the c direction of ζx > 14 nm. In addition, there are maxima around [h, k, l + ! " ] : h, k, l ∈ Z within these planes of DS, corresponding to a shorter-range correlation of ζx = 3.9 nm. Accordingly, we observe strong features in the 0yz and x0z layers of the 3D-∆PDF, with alternating signs along the z direction and the same sign along the x-and y directions [Fig. ]. We first note that the features of alternating signs along z correspond to the observed planes of DS, using the condition ℓ = n + ! " : n ∈ Z. Therefore, antiferrohexapolar correlations apply along the z direction. Second, we derive the shorter-range correlations in the xy plane to be ferrohexapolar, based on the same sign of the features within these planes [Fig. )]. Thus, while our electrostatic arguments correctly predicted antiferrohexapolar correlations along z, they cannot explain the ferrohexapolar interactions correlations present in the xy plane. To better understand the correlated disorder governing the system, we turn to a microscopic model of the SBU orientations using a Monte Carlo (MC) approach. This method involves generating a representative sample of the statistical distribution of disordered super-cells. It is evident from the 3D-∆PDF that the node orientations are not random, hence a method of including these correlations in our simulation is necessary. Given that the correlation length along c is long, we make the approximation that the disorder is present only in the ab plane, with longrange antiferrohexapolar correlations persisting in the c direction. Within the plane, the tendency for ferrohexapolar interactions can be encoded in the Hamiltonian,