def METHOD_NAME(xi=None, d=None, bjerrum_length=None, sigma=None, valency=None):
    el_char = 1.0
    return xi * np.tan(xi * d / 2.0) + 2.0 * np.pi * \
        bjerrum_length * sigma / (valency * el_char)