  We study the distribution of multivalent counterions next to a dielectric
slab, bearing a quenched, random distribution of charges on one of its solution
interfaces, with a given mean and variance, both in the absence and in the
presence of a bathing monovalent salt solution. We use the previously derived
approach based on the dressed multivalent-ion theory that combines aspects of
the strong and weak coupling of multivalent and monovalent ions in a single
framework. The presence of quenched charge disorder on the charged surface of
the dielectric slab is shown to substantially increase the density of
multivalent counterions in its vicinity. In the counterion-only model (with no
monovalent salt ions), the surface disorder generates an additional logarithmic
attraction potential and thus an algebraically singular counterion density
profile at the surface. This behavior persists also in the presence of a
monovalent salt bath and results in significant violation of the contact-value
theorem, reflecting the anti-fragility effects of the disorder that drive the
system towards a more ordered state. In the presence of an interfacial
dielectric discontinuity, depleting the counterion layer at the surface, the
charge disorder still generates a much enhanced counterion density further away
from the surface. Likewise, the charge inversion and/or overcharging of the
surface occur more strongly and at smaller bulk concentrations of multivalent
counterions when the surface carries quenched charge disorder. Overall, the
presence of quenched surface charge disorder leads to sizable effects in the
distribution of multivalent counterions in a wide range of realistic parameters
and typically within a distance of a few nanometers from the charged surface.
