Interaction of fast and slow varying electromagnetic waves propagating in paraelectric or ferroelectric material

Once, a referee asked how one can write an oscillator model for a ferroelectric and the total Maxwell equation for a rapidly varying electric field. From what is known about polarization, it is a slow function of time and coordinates, but the optical wave is a fast function. However, there are examples for the interaction of high frequency and low frequency waves in nonlinear wave theory. This means that similar equations can be written for ferroelectric polarization and electromagnetic waves.