To illustrate this difference in predicted behaviour between a simple charge-compensation model and our full self-consistent model, we consider the doping response of Li 3 OCl as a function of effective dopant concentration, r[M r ], for both a simple "lithium-only" charge compensation model, where only charge-carrying lithium defects change concentration in response to doping, and for our full model, where all defects re-equilibrate. Figure shows the predicted doping response for these two models, for supervalent doping and subvalent doping at a range of dopant concentrations, under Li-rich and Lipoor synthesis conditions. Here, we have considered a maximum dopant concentration equivalent to ∼ 4 % of lithium sites, which is a similar dopant/defect concentration to that used in previous ab initio molecular dynamics studies of doping in Li 3 OCl . Because our model assumes that defects behave ideally (via our use of eq. ( )) our results are expected to be quantitatively accurate only at low dopant and defect concentrations. At high dopant or high defect concentrations, non-ideal defect-defect interactions may make a significant contribution to defect chemical potentials, causing defect concentrations to deviate from the dilute-limit behaviour modelled here . We expect, however, that general trends regarding the relative response of different defects to doping will be qualitatively correct even at these high dopant concentrations.