calculate sum_a I_a sqrt(tau_a)