Send in a pulse of light, not a pulse of voltage

[apt1002]

I confess I haven't run the code, but I did read it, and I noticed this:

def original_electric_field(t, y):
    #return plane_wave(t, y, angular_frequency)
    return gaussian(t, y)

This describes a short voltage pulse, centred on frequency zero. However, what you actually want to model is a pulse of light, centred on a non-zero frequency. I think you need to multiply the Gaussian by the plane wave.

I haven't tried it, sorry. Maybe I have misunderstood.