Bounds for the product of modified Bessel functions