scipy.special.lambertw
scipy.special.zetac
Riemann zeta function.
The two-argument version is the Hurwitz zeta function:
Riemann zeta function corresponds to q = 1.
See also
zetac