/*hmirat compute radiative association rate for H- */
double hmirat(double te);