Hi Chris,
I believe your problem is similar to what is discussed here.
In your case, a workaround is to use the Field version of propagation, more precisely with the Gradient class. You’ll be able to retrieve in your ephemerides the first order partial derivatives w.r.t. the parameters you defined.
Romain.