C_n_numu in thermo has a stray T_nue where it should have a T_numt:
C_n_numu = 8 * f_n * const.GF**2 / jnp.pi**5 * (
4 * (gmuL**2 + gmuR**2)
* (T_g**8 - T_nue**8 * jnp.exp(2 * mu_numt / T_numt))
- (
T_numt**8 * jnp.exp(2 * mu_numt / T_numt)
- T_nue**8 * jnp.exp(2 * mu_nue / T_nue)
)
) / const.hbar
should instead be
C_n_numu = 8 * f_n * const.GF**2 / jnp.pi**5 * (
4 * (gmuL**2 + gmuR**2)
* (T_g**8 - T_numt**8 * jnp.exp(2 * mu_numt / T_numt))
- (
T_numt**8 * jnp.exp(2 * mu_numt / T_numt)
- T_nue**8 * jnp.exp(2 * mu_nue / T_nue)
)
) / const.hbar
C_n_numuinthermohas a strayT_nuewhere it should have aT_numt:should instead be