Improve fitting of broken power-law models

[astrofrog]

As discussed in #6198, the broken power-law models may behave better when fitting if parametrized with log_x_break instead of x_break. There are several options:

1. Change both models to log_x_break and deprecate x_break
2. Make it possible to describe the model with x_break and use log_x_break internally
3. Implement two versions of the model, one with x_break and one with log_x_break (this could be done with minimal code duplication).
4. Make it so that Parameters can indicate whether to be fit in log space (or according to other functions)

@gcalderone suggested that 3 might be the easiest to do. Just opening this issue so that we remember to do it for the next version.

(I'd also like to see examples where the log version would work and the non-log version fails, to make sure we really need this)