Besides the change in quadrupole splitting for 1, the changes observed upon increasing the system size are within the trust regions. It would thus not be possible to determine which system size better represents experimental data. In other words, the Mössbauer parameters are converged at the smaller system size, and hence we will focus on the smaller models in the following.