Allow callback for l_bfgs_b. Allow kwargs to optimization functions.

If you decide to bumpy minimum scipy version, might even be worth using the new minimize API. We haven't done it in statsmodels yet though FWIW.