def clumping_factor_HB(z, beta=2):
    """Clumping factor as a function of redshift used by Haiman & Bryan (2006).

    See Haiman & Bryan (2006ApJ...650....7H).
    """
    return 1 + 9 * ((1 + z) / 7) ** -beta
