Lambert series#Factorization theorems

{{short description|Mathematical term}}

{{For|generalized Lambert series|Appell–Lerch sum}}

File:Cplot Lambert series.png plot, using a version of the domain coloring method{{Cite web | url=http://nbviewer.ipython.org/github/empet/Math/blob/master/DomainColoring.ipynb | title=Jupyter Notebook Viewer}}]]

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

:S(q)=\sum_{n=1}^\infty a_n \frac {q^n}{1-q^n}.

It can be resummed formally by expanding the denominator:

:S(q)=\sum_{n=1}^\infty a_n \sum_{k=1}^\infty q^{nk} = \sum_{m=1}^\infty b_m q^m

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

:b_m = (a*1)(m) = \sum_{n\mid m} 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_{n=1}^\infty q^n \sigma_0(n) = \sum_{n=1}^\infty \frac{q^n}{1-q^n}

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_{n=1}^\infty q^n \sigma_\alpha(n) = \sum_{n=1}^\infty \frac{n^\alpha q^n}{1-q^n} = \sum_{n=1}^\infty \operatorname{Li}_{-\alpha}(q^n)

where \alpha is any complex number, \operatorname{Li} is the polylogarithm, and

:\sigma_\alpha(n) = (\textrm{Id}_\alpha*1)(n) = \sum_{d\mid n} d^\alpha \,

is the divisor function. In particular, for \alpha = 1, the Lambert series one gets is

:q \frac{F'(q)}{F(q)}

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

:F(q) := \frac{1}{\phi(q)} = \sum_{k=0}^\infty p(k) q^k = \prod_{n=1}^\infty \frac{1}{1-q^n}.

Additional Lambert series related to the previous identity include those for the variants of the

Möbius function given below \mu(n)

:See the forum post [https://mathoverflow.net/q/98174 here] (or the article {{arXiv|1112.4911}}) and the conclusions section of {{ArXiv|1712.00611}} by Merca and Schmidt (2018) for usage of these two less standard Lambert series for the Moebius function in practical applications.

:\sum_{n=1}^\infty \mu(n)\,\frac{q^n}{1-q^n} = q.

Related Lambert series over the Möbius function include the following identities for any

prime \alpha \in \mathbb{Z}^{+}:

:

\sum_{n \geq 1} \frac{\mu(n) q^n}{1+q^n} = q-2q^2

:

\sum_{n \geq 1} \frac{\mu(\alpha n) q^n}{1-q^n} = -\sum_{n \geq 0} q^{\alpha^n}.

{{Citation needed|date=May 2024}}

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_{f}(q) := q to be the Lambert series generating function of the arithmetic function f:

:

\begin{align}

\sum_{n \geq 1} \frac{f(n) q^n}{1+q^n} & = \sum_{n \geq 1} \frac{f(n)q^{n}}{1-q^{n}} -

\sum_{n \geq 1} \frac{2 f(n) q^{2n}}{1-q^{2n}} \\

& =

L_f(q) - 2 \cdot L_f(q^2).

\end{align}

For Euler's totient function \varphi(n):

:\sum_{n=1}^\infty \varphi(n)\,\frac{q^n}{1-q^n} = \frac{q}{(1-q)^2}.

For Von Mangoldt function \Lambda(n):

:\sum_{n=1}^\infty \Lambda(n)\,\frac{q^n}{1-q^n} = \sum_{n=1}^{\infty} \log(n)q^n

For Liouville's function \lambda(n):

:\sum_{n=1}^\infty \lambda(n)\,\frac{q^n}{1-q^n} =

\sum_{n=1}^\infty q^{n^2}

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 an are trigonometric functions, for example, an = sin(2n 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^{th} powers, n = k^m \in \mathbb{Z}^{+}, 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_{d^m|n} \mu\left(\frac{n}{d^m}\right), which in turn shows that

:\sum_{n \geq 1} \frac{\lambda_m(n) q^n}{1-q^n} = \sum_{n \geq 1} q^{n^m},\ \text{ for } 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

{{cite web|last1=Weisstein|first1=Eric W.|title=Lambert Series|url=http://mathworld.wolfram.com/LambertSeries.html|website=MathWorld|access-date=22 April 2018}}

:\sum_{n=1}^{\infty} \frac{4 \cdot (-1)^{n+1} q^{2n+1}}{1-q^{2n+1}} = \sum_{m=1}^{\infty} 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{Id}_1), (\lambda, \chi_{\operatorname{sq}}), (\Lambda, \log),

(|\mu|, 2^{\omega}), (J_t, \operatorname{Id}_t), (d^3, (d \ast 1)^2),

where \varepsilon(n) = \delta_{n,1} is the multiplicative identity for Dirichlet convolutions, \operatorname{Id}_k(n) = n^k is the identity function for k^{th} powers, \chi_{\operatorname{sq}} 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^{-z} one obtains another common form for the series, as

:\sum_{n=1}^\infty \frac {a_n}{e^{zn}-1}= \sum_{m=1}^\infty b_m e^{-mz}

where

:b_m = (a*1)(m) = \sum_{d\mid m} 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{Li}_0(q^{n}) is a polylogarithm function, we may refer to any sum of the form

:\sum_{n=1}^{\infty} \frac{\xi^n \,\mathrm{Li}_u (\alpha q^n)}{n^s} = \sum_{n=1}^{\infty} \frac{\alpha^n \,\mathrm{Li}_s(\xi q^n)}{n^u}

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

:12\left(\sum_{n=1}^{\infty} n^2 \, \mathrm{Li}_{-1}(q^n)\right)^{\!2} = \sum_{n=1}^{\infty}

n^2 \,\mathrm{Li}_{-5}(q^n) -

\sum_{n=1}^{\infty} n^4 \, \mathrm{Li}_{-3}(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{{cite journal|last1=Merca|first1=Mircea|title=The Lambert series factorization theorem|journal=The Ramanujan Journal|date=13 January 2017|volume=44|issue=2|pages=417–435|doi=10.1007/s11139-016-9856-3|s2cid=125286799}}

:\sum_{n \geq 1} \frac{a_n q^n}{1\pm q^n} = \frac{1}{(\mp q; q)_{\infty}} \sum_{n \geq 1} \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) = [q^n] (\mp q; q)_{\infty} \frac{q^k}{1 \pm q^k} is the respective sum or difference of the

restricted partition functions s_{e/o}(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_{n,k} := s_e(n, k) - s_o(n, k) = [q^n] (q; q)_{\infty} \frac{q^k}{1-q^k} denote the invertible lower triangular sequence whose first few values are shown in the table below.

class="wikitable"
n \ k12345678
110000000
201000000
3-1-1100000
4-10-110000
5-1-1-1-11000
6001-1-1100
700-10-1-110
810010-1-11

Another characteristic form of the Lambert series factorization theorem expansions is given by{{cite journal|author=Merca, M.|author2=Schmidt, M. D.|name-list-style=amp|title=Generating Special Arithmetic Functions by Lambert Series Factorizations|journal=Contributions to Discrete Mathematics|date=2019|volume=14|issue=1|pages=31–45|doi=10.11575/cdm.v14i1.62425|doi-access=free|arxiv=1706.00393|bibcode=2017arXiv170600393M}}

:L_f(q) := \sum_{n \geq 1} \frac{f(n) q^n}{1-q^n} = \frac{1}{(q; q)_{\infty}} \sum_{n \geq 1} \left(s_{n,k} f(k)\right) q^n,

where (q; q)_{\infty} 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_{n,k}^{(-1)} = \sum_{d|n} p(d-k) \mu\left(\frac{n}{d}\right)

The next table lists the first several rows of these corresponding inverse matrices.{{cite web|title=A133732|url=https://oeis.org/A133732|website=Online Encyclopedia of Integer Sequences|access-date=22 April 2018}}

class="wikitable"
n \ k12345678
110000000
201000000
311100000
421110000
543211000
653221100
7107532110
8129643211

We let G_j := \frac{1}{2} \left\lceil \frac{j}{2} \right\rceil \left\lceil \frac{3j+1}{2} \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)_{\infty} = \sum_{n \geq 0} (-1)^{\left\lceil \frac{n}{2} \right\rceil} q^{G_n}.

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{{cite journal|last1=Schmidt|first1=Maxie D.|title=New Recurrence Relations and Matrix Equations for Arithmetic Functions Generated by Lambert Series|journal=Acta Arithmetica|date=8 December 2017|volume=181|issue=4|pages=355–367|doi=10.4064/aa170217-4-8|arxiv=1701.06257|bibcode=2017arXiv170106257S|s2cid=119130467}}

:f(n) = \sum_{k=1}^n \sum_{d|n} p(d-k) \mu(n/d) \times \sum_{j: k-G_j > 0} (-1)^{\left\lceil \frac{j}{2} \right\rceil} b(k-G_j).

This work on Lambert series factorization theorems is extended in{{cite arXiv|author=M. Merca|author2=Schmidt, M. D.|name-list-style=amp|title=New Factor Pairs for Factorizations of Lambert Series Generating Functions|eprint=1706.02359|class=math.CO|year=2017}} to more general expansions of the form

:\sum_{n \geq 1} \frac{a_n q^n}{1-q^n} = \frac{1}{C(q)} \sum_{n \geq 1} \left(\sum_{k=1}^n s_{n,k}(\gamma) \widetilde{a}_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{a}_k(\gamma) = \sum_{d|k} \sum_{r| \frac{k}{d}} a_d \gamma(r).

The corresponding inverse matrices in the above expansion satisfy

:s_{n,k}^{(-1)}(\gamma) = \sum_{d|n} [q^{d-k}] \frac{1}{C(q)} \gamma\left(\frac{n}{d}\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{a}_k(\gamma) = \sum_{k=1}^{n} s_{n,k}^{(-1)}(\gamma) \times [q^k]\left(\sum_{d=1}^k \frac{a_d q^d}{1-q^d} 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_{1 \leq n \leq x} g_f(n).

We also adopt the notation from the previous section that

:s_{n,k} = [q^n] (q; q)_{\infty} \frac{q^k}{1-q^k},

where (q; q)_{\infty} 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_{b = \pm 1} \sum_{k=1}^{\left\lfloor \frac{\sqrt{24n+1}-b}{6}\right\rfloor}

(-1)^{k+1} g_f\left(n+1-\frac{k(3k+b)}{2}\right) +

\sum_{k=1}^{n+1} s_{n+1,k} f(k),

:\Sigma_f(x+1) = \sum_{b = \pm 1} \sum_{k=1}^{\left\lfloor \frac{\sqrt{24x+1}-b}{6}\right\rfloor}

(-1)^{k+1} \Sigma_f\left(n+1-\frac{k(3k+b)}{2}\right) +

\sum_{n=0}^x \sum_{k=1}^{n+1} s_{n+1,k} 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^{th} derivatives of a Lambert series for any s \geq 1{{cite arXiv|last1=Schmidt|first1=Maxie D.|title=Combinatorial Sums and Identities Involving Generalized Divisor Functions with Bounded Divisors|date=2017|eprint=1704.05595|class=math.NT}}{{cite arXiv|last1=Schmidt|first1=Maxie D.|title=Factorization Theorems for Hadamard Products and Higher-Order Derivatives of Lambert Series Generating Functions|date=2017|eprint=1712.00608|class=math.NT}}

:q^s \cdot D^{(s)}\left[\frac{q^i}{1-q^i}\right] = \sum_{m=0}^s \sum_{k=0}^m \left[\begin{matrix} s \\ m\end{matrix}\right]

\left\{\begin{matrix} m \\ k\end{matrix}\right\} \frac{(-1)^{s-k} k! i^m}{(1-q^i)^{k+1}}

:q^s \cdot D^{(s)}\left[\frac{q^i}{1-q^i}\right] = \sum_{r=0}^s\left[\sum_{m=0}^s \sum_{k=0}^m \left[\begin{matrix} s \\ m\end{matrix}\right]

\left\{\begin{matrix} m \\ k\end{matrix}\right\} \binom{s-k}{r} \frac{(-1)^{s-k-r} k! i^m}{(1-q^i)^{k+1}}\right] q^{(r+1)i},

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

:[q^n]\left(\sum_{i \geq t} \frac{a_i q^{mi}}{(1-q^i)^{k+1}}\right) = \sum_{\begin{matrix} d|n \\ t \leq d \leq \left\lfloor \frac{n}{m} \right\rfloor\end{matrix}}

\binom{\frac{n}{d}-m+k}{k} a_d.

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

:A_t(n) := \sum_{\begin{matrix} 0 \leq k \leq m \leq t \\ 0 \leq r \leq t\end{matrix}} \sum_{d|n} \left[\begin{matrix} t \\ m\end{matrix}\right]

\left\{\begin{matrix} m \\ k\end{matrix}\right\} \binom{t-k}{r} \binom{\frac{n}{d}-1-r+k}{k} (-1)^{t-k-r} k! d^m \cdot a_d \cdot

\left[t \leq d \leq \left\lfloor \frac{n}{r+1} \right\rfloor\right]_{\delta},

where [\cdot]_{\delta} denotes Iverson's convention, then we have the coefficients for the t^{th} derivatives of a Lambert series

given by

:\begin{align}

A_t(n) & = [q^n]\left(q^t \cdot D^{(t)}\left[\sum_{i \geq t} \frac{a_i q^i}{1-q^i}\right]\right) \\

& = [q^n]\left(\sum_{n \geq 1} \frac{(A_t \ast \mu)(n) q^n}{1-q^n}\right).

\end{align}

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

:[q^n]\left(q^t \cdot D^{(t)}\left[\sum_{i \geq 1} \frac{f(i) q^i}{1-q^i}\right]\right) = \frac{n!}{(n-t)!} \cdot (f \ast 1)(n).

See also

References

{{Reflist}}

  • {{cite book|last=Berry|first=Michael V. |title=Functions of Number Theory|year=2010|publisher=CAMBRIDGE UNIVERSITY PRESS|isbn=978-0-521-19225-5|pages=637–641|url=http://dlmf.nist.gov}}
  • {{cite journal|first1=Preston A. | last1=Lambert

|title=Expansions of algebraic functions at singular points

|journal=Proc. Am. Philos. Soc.

|year=1904

|volume=43 | number=176

|jstor=983503

|pages=164–172

}}

  • {{Apostol IANT}}
  • {{springer|title=Lambert series|id=p/l057340}}
  • {{MathWorld|urlname=LambertSeries|title=Lambert Series}}
  • {{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}}

Category:Analytic number theory

Category:Q-analogs

Category:Series (mathematics)