Lambert Series

In mathematics, a Lambert series, named for Johann Heinrich Lambert, is a series taking the form

:$S(q)=\sum_^\infty a_n \frac .$

It can be resumed formally by expanding the denominator:

:$S(q)=\sum_^\infty a_n \sum_^\infty q^ = \sum_^\infty b_m q^m$

where the coefficients of the new series are given by the Dirichlet convolution of ''a''_''n'' with the constant function 1(''n'') = 1:

:$b_m = (a*1)(m) = \sum_ a_n. \,$

This series may be inverted by means of the Möbius inversion formula, and is an example of a Möbius transform.

Examples

Since this last sum is a typical number-theoretic sum, almost any natural multiplicative function will be exactly summable when used in a Lambert series. Thus, for example, one has

:$\sum_^\infty q^n \sigma_0(n) = \sum_^\infty \frac$

where $\sigma_0(n)=d(n)$ is the number of positive divisors of the number ''n''.

For the higher order sum-of-divisor functions, one has

:$\sum_^\infty q^n \sigma_\alpha(n) = \sum_^\infty \frac$

where $\alpha$ is any complex number and

:$\sigma_\alpha(n) = (\textrm_\alpha*1)(n) = \sum_ d^\alpha \,$
is the divisor function. In particular, for $\alpha = 1$, the Lambert series one gets is

:$q \frac$

which is (up to the factor of $q$) the logarithmic derivative of the usual generating function for partition numbers

:$F(q) := \frac = \sum_^\infty p(k) q^k = \prod_^\infty \frac.$

Additional Lambert series related to the previous identity include those for the variants of the
Möbius function given below $\mu(n)$
:

:$\sum_^\infty \mu(n)\,\frac = q.$
Related Lambert series over the Moebius function include the following identities for any
prime $\alpha \in \mathbb^$:
:$\begin \sum_ \frac & = q-2q^2 \\ \sum_ \frac & = -\sum_ q^. \end$
The proof of the first identity above follows from a multi-section (or bisection) identity of these
Lambert series generating functions in the following form where we denote
$L_(q) := q$ to be the Lambert series generating function of the arithmetic function ''f'':
:$\begin \sum_ \frac & = \sum_ \frac - \sum_ \frac \\ & = L_f(q) - 2 \cdot L_f(q^2). \end$
The second identity in the previous equations follows from the fact that the coefficients of the left-hand-side sum are given by
:$\begin \sum_ \mu(\alpha d) = \sum_ \mu(d) = \varepsilon\left(\frac\right) = 1 \iff n = (n, \alpha) \iff n = \alpha^k,\ \text k \geq 1, \end$
where the function $\varepsilon(n) = \delta_$ is the multiplicative identity with respect to the operation of Dirichlet convolution of arithmetic functions.

For Euler's totient function $\varphi(n)$:
:$\sum_^\infty \varphi(n)\,\frac = \frac.$

For Von Mangoldt function $\Lambda(n)$:

:$\sum_^\infty \Lambda(n)\,\frac = \sum_^ \log(n)q^n$

For Liouville's function $\lambda(n)$:

:$\sum_^\infty \lambda(n)\,\frac = \sum_^\infty q^$

with the sum on the right similar to the Ramanujan theta function, or Jacobi theta function $\vartheta_3(q)$. Note that Lambert series in which the ''a''_''n'' are trigonometric functions, for example, ''a''_''n'' = sin(2''n'' ''x''), can be evaluated by various combinations of the logarithmic derivatives of Jacobi theta functions.

Generally speaking, we can extend the previous generating function expansion by letting $\chi_m(n)$ denote the characteristic function of the $m^$ powers, $n = k^m \in \mathbb^$, for positive natural numbers $m > 2$ and defining the generalized ''m''-Liouville lambda function to be the arithmetic function satisfying $\chi_m(n) := (1 \ast \lambda_m)(n)$. This definition of $\lambda_m(n)$ clearly implies that $\lambda_m(n) = \sum_ \mu\left(\frac\right)$, which in turn shows that

:$\sum_ \frac = \sum_ q^,\ \text m \geq 2.$

We also have a slightly more generalized Lambert series expansion generating the sum of squares function $r_2(n)$ in the form of

:$\sum_^ \frac = \sum_^ r_2(m) q^m.$

In general, if we write the Lambert series over $f(n)$ which generates the arithmetic functions $g(m) = (f \ast 1)(m)$, the next pairs of functions correspond to other well-known convolutions expressed by their Lambert series generating functions in the forms of

:$(f, g) = (\mu, \varepsilon), (\varphi, \operatorname_1), (\lambda, \chi_), (\Lambda, \log), (, \mu, , 2^), (J_t, \operatorname_t), (d^3, (d \ast 1)^2),$

where $\varepsilon(n) = \delta_$ is the multiplicative identity for Dirichlet convolutions, $\operatorname_k(n) = n^k$ is the identity function for $k^$ powers, $\chi_$ denotes the characteristic function for the squares, $\omega(n)$ which counts the number of distinct prime factors of $n$ (see prime omega function), $J_t$ is Jordan's totient function, and $d(n) = \sigma_0(n)$ is the divisor function (see Dirichlet convolutions).

The conventional use of the letter ''q'' in the summations is a historical usage, referring to its origins in the theory of elliptic curves and theta functions, as the nome.

Alternate form

Substituting $q=e^$ one obtains another common form for the series, as

:$\sum_^\infty \frac = \sum_^\infty b_m e^$

where
:$b_m = (a*1)(m) = \sum_ a_d\,$

as before. Examples of Lambert series in this form, with $z=2\pi$, occur in expressions for the Riemann zeta function for odd integer values; see Zeta constants for details.

Current usage

In the literature we find ''Lambert series'' applied to a wide variety of sums. For example, since $q^n/(1 - q^n ) = \mathrm_0(q^)$ is a polylogarithm function, we may refer to any sum of the form

:$\sum_^ \frac = \sum_^ \frac$

as a Lambert series, assuming that the parameters are suitably restricted. Thus

:$12\left(\sum_^ n^2 \, \mathrm_(q^n)\right)^ = \sum_^ n^2 \,\mathrm_(q^n) - \sum_^ n^4 \, \mathrm_(q^n),$

which holds for all complex ''q'' not on the unit circle, would be considered a Lambert series identity. This identity follows in a straightforward fashion from some identities published by the Indian mathematician S. Ramanujan. A very thorough exploration of Ramanujan's works can be found in the works by Bruce Berndt.

Factorization theorems

A somewhat newer construction recently published over 2017–2018 relates to so-termed ''Lambert series factorization theorems'' of the form

:$\sum_ \frac = \frac \sum_ \left((s_o(n, k) \pm s_e(n, k)) a_k\right) q^n,$

where s_o(n, k) \pm s_e(n, k) = ^n(\mp q; q)_ \frac is the respective sum or difference of the
restricted partition functions $s_(n, k)$ which denote the number of $k$'s in all partitions of $n$ into an ''even'' (respectively, ''odd'') number of distinct parts. Let s_ := s_e(n, k) - s_o(n, k) = ^n(q; q)_ \frac denote the invertible lower triangular sequence whose first few values are shown in the table below.

Another characteristic form of the Lambert series factorization theorem expansions is given by

:$L_f(q) := \sum_ \frac = \frac \sum_ \left(s_ f(k)\right) q^n,$

where $(q; q)_$ is the (infinite) q-Pochhammer symbol. The invertible matrix products on the right-hand-side of the previous equation correspond to inverse matrix products whose lower triangular entries are given in terms of the partition function and the Möbius function by the divisor sums

:$s_^ = \sum_ p(d-k) \mu\left(\frac\right)$

The next table lists the first several rows of these corresponding inverse matrices.

We let $G_j := \frac \left\lceil \frac \right\rceil \left\lceil \frac \right\rceil$ denote the sequence of interleaved pentagonal numbers, i.e., so that the pentagonal number theorem is expanded in the form of

:$(q; q)_ = \sum_ (-1)^ q^.$

Then for any Lambert series $L_f(q)$ generating the sequence of $g(n) = (f \ast 1)(n)$, we have the corresponding inversion relation of the factorization theorem expanded above given by

:$f(n) = \sum_^n \sum_ p(d-k) \mu(n/d) \times \sum_ (-1)^ b(k-G_j).$

This work on Lambert series factorization theorems is extended in to more general expansions of the form

:$\sum_ \frac = \frac \sum_ \left(\sum_^n s_(\gamma) \widetilde_k(\gamma)\right) q^n,$

where $C(q)$ is any (partition-related) reciprocal generating function, $\gamma(n)$ is any arithmetic function, and where the
modified coefficients are expanded by

:$\widetilde_k(\gamma) = \sum_ \sum_ a_d \gamma(r).$

The corresponding inverse matrices in the above expansion satisfy

:s_^(\gamma) = \sum_ ^\frac \gamma\left(\frac\right),

so that as in the first variant of the Lambert factorization theorem above we obtain an inversion relation for the right-hand-side coefficients of the form

:\widetilde_k(\gamma) = \sum_^ s_^(\gamma) \times ^kleft(\sum_^k \frac C(q)\right).

Recurrence relations

Within this section we define the following functions for natural numbers $n,x \geq 1$:
:$g_f(n) := (f \ast 1)(n),$
:$\Sigma_f(x) := \sum_ g_f(n).$

We also adopt the notation from the previous section that

:s_ = ^n(q; q)_ \frac,

where $(q; q)_$ is the infinite q-Pochhammer symbol. Then we have the following recurrence relations for involving these functions and the pentagonal numbers proved in:

:$g_f(n+1) = \sum_ \sum_^ (-1)^ g_f\left(n+1-\frac\right) + \sum_^ s_ f(k),$
:$\Sigma_f(x+1) = \sum_ \sum_^ (-1)^ \Sigma_f\left(n+1-\frac\right) + \sum_^x \sum_^ s_ f(k).$

Derivatives

Derivatives of a Lambert series can be obtained by differentiation of the series termwise with respect to $q$. We have the following identities for the termwise $s^$ derivatives of a Lambert series for any $s \geq 1$

:q^s \cdot D^\left frac\right= \sum_^s \sum_^m \left begin s \\ m\end\right
\left\ \frac
:q^s \cdot D^\left frac\right= \sum_^s\left begin_s_\\_m\end\right
_______\left\_\binom_\frac\right.html" ;"title="sum_^s \sum_^m \left begin s \\ m\end\right
\left\ \binom \frac\right">sum_^s \sum_^m \left begin s \\ m\end\right
\left\ \binom \frac\rightq^,

where the bracketed triangular coefficients in the previous equations denote the Stirling numbers of the first and second kinds.
We also have the next identity for extracting the individual coefficients of the terms implicit to the previous expansions given in the form of

: ^nleft(\sum_ \frac\right) = \sum_
\binom a_d.

Now if we define the functions $A_t(n)$ for any $n,t \geq 1$ by

:A_t(n) := \sum_ \sum_ \left begin t \\ m\end\right
\left\ \binom \binom (-1)^ k! d^m \cdot a_d \cdot
\left \leq d \leq \left\lfloor \frac \right\rfloor\right,

where $$cdot denotes Iverson's convention, then we have the coefficients for the $t^$ derivatives of a Lambert series
given by

:\begin
A_t(n) & = ^nleft(q^t \cdot D^\left sum_ \frac\rightright) \\
& = ^nleft(\sum_ \frac\right).
\end

Of course, by a typical argument purely by operations on formal power series we also have that

: ^nleft(q^t \cdot D^\left sum_ \frac\rightright) = \frac \cdot (f \ast 1)(n).

See also

* Erdős–Borwein constant
* Arithmetic function
* Dirichlet convolution

References

*
*
*
*
*
* {{cite arXiv, last=Schmidt, first=Maxie Dion, date=2020-04-06, title=A catalog of interesting and useful Lambert series identities, class=math.NT, eprint=2004.02976

Analytic number theory
Q-analogs
Mathematical series