Fix bug on bracketing in calculate_ref_properties

Signed-off-by: Howard Cheung <howard.at@gmail.com>
Add a function to estimate coefficients of model by Limited BFGS-B method.