We consider the propagation of high frequency electromagnetic pulses in complex materials with nonlinear forced polarization. The physical problem is modeled by Maxwell’s equations in variational form, and well-posedness results are established with respect to probability distributions on the polarization parameters (in a Prohorov metric sense).