Perhaps just `1.0 - Gaussian1D.eval(x, amplitude, mean, stddev)`?  Or would that be too cheeky?  It should also have a `fit_deriv` method which should just be -1 times the result of `Gaussian1D.fit_deriv`.