Figure shows the correlation between the polarizable embedding schemes compared to the reference electrostatic embedding. The analysis is performed only for the MBIS charges, since this model exhibited the lower MRAD in the previous analysis. These results are obtained evaluating equation ( ), using MBIS charges and two types of polarizabilities. A substantial improvement in the estimation of the coupling energy is observed when atomic polarizabilities are incorporated for each element. The relative difference between each approximation and the reference energy obtained from DFT calculations is showcased by the histograms in the bottom panels from Figure ). Data points corresponding to stronger ML/MM interactions (i.e. more negative coupling energies energies) exhibit larger deviations in the correlation plots. The polarization correction is able to correct the coupling energy even for data points across the whole sampled range of E EE QM/MM values, reflecting the polarization energy becomes larger in module as the coupling energy becomes larger. This feature outlines that structures in which the ML and the MM subsystem interact strongly, are expected to exhibit stronger polarization effects.