Use BoundedSecant solver in density_Ts

- constrain density to positive values