General formulas for a class of Euler sums

Titel: Algemene formules voor een klasse van Euler-sommen

Auteurs: David H. Bailey, Ross McPhedran en Bruno Salvy
Datum: 6 april 2026

1. Het Probleem

Euler-sommen (ook wel Euler-Zagier-sommen of harmonische sommen genoemd) zijn oneindige reeksen die harmonische getallen $H_k = \sum_{j=1}^k 1/j$ bevatten. Het paper richt zich op de analytische evaluatie van sommen van de vorm:
$$\sum_{k=1}^{\infty} R(k)H_k$$
waarbij $R(k)$ een rationale functie is (een quotiënt van twee polynomen) met de voorwaarde dat de graad van de noemer ten minste twee hoger is dan de graad van de teller (voor convergentie).

Specifiek behandelen de auteurs gevallen met gemengde noemers, zoals:
$$\sum_{k=1}^{\infty} \frac{H_k}{(m_1k + n_1)^{p_1}(m_2k + n_2)^{p_2} \cdots (m_rk + n_r)^{p_r}}$$
en varianten met een teller $k^q$. Hoewel veel van deze sommen elegante analytische oplossingen hebben (vaak uitgedrukt in termen van de Riemann-zètafunctie), ontbrak er tot nu toe een systematische, algemene methode om gesloten vormen af te leiden voor willekeurige rationale functies $R(k)$, vooral bij complexe noemers.

2. Methodologie

De kern van de oplossing ligt in de residumethode uit de complexe analyse, gekoppeld aan eigenschappen van de digamma- ($\psi$) en polygammafuncties.

• Residumethode: De auteurs gebruiken een stelling (Propositie 1) die de som relateert aan de som van de residuen van een specifieke complexe functie. Voor een rationale functie $R(s)$ geldt:
  $$\sum_{k=1}^{\infty} R(k)H_k = \frac{1}{2} \sum_{\alpha \in \text{polen van } R} \text{Res}_{s=\alpha} \left[ R'(s)(\psi(-s) + \gamma) - R(s)(\psi(-s) + \gamma)^2 \right]$$
  Hierbij is $\gamma$ de constante van Euler-Mascheroni en $\psi(z) = \Gamma'(z)/\Gamma(z)$.
• Symbolische Berekening: In plaats van de polen numeriek te benaderen, gebruiken de auteurs een symbolische aanpak. Ze ontleden de noemer in irreducibele factoren en berekenen de Taylor- en Laurent-reeksen rond de polen (die wortels van polynomen kunnen zijn). Dit maakt het mogelijk om exacte formules af te leiden zonder expliciete numerieke waarden van de polen te hoeven kennen.
• Numerieke Verificatie: Om de afgeleide formules te verifiëren, hebben de auteurs een hoogprecisie numeriek algoritme ontwikkeld. Dit omvat het berekenen van de eerste $10^8$ termen van de reeks en het gebruik van de Euler-Maclaurin-somformule voor de staart. Vervolgens wordt het PSLQ-algoritme (Integer Relation Detection) gebruikt om de rationele coëfficiënten van de gesloten vorm te identificeren op basis van de numerieke waarde.

3. Belangrijkste Bijdragen

1. Algoritme HarmonicSum:
   De auteurs presenteren een implementeerbaar algoritme (in Maple en Mathematica) dat automatisch de gesloten vorm berekent voor elke convergente som $\sum R(k)H_k$. Het algoritme verwerkt de partial fraction decomposition en berekent de residuen symbolisch.

2. Algemene Formules:
   Er worden expliciete formules afgeleid voor:

   • Enkele polen: Formules voor $\sum \frac{H_k}{(k+t)^p}$ uitgedrukt in termen van $\psi$ en $\gamma$.
   • Gemengde noemers: Formules voor producten van lineaire factoren in de noemer (bijv. $(2k+1)^p(3k+1)^q$).
   • Teller met macht: Uitbreidingen voor sommen met een teller $k^q$.

3. Uitgebreide Catalogus:
   Het paper bevat een uitgebreide catalogus (Appendix B) met honderden specifieke formules voor gevallen waar de machten $p_i$ tot 3 lopen en de multipliers $m_i$ tot 4. Deze formules zijn uitgedrukt in termen van $\psi^{(n)}(z)$ (polygammafuncties) en $\gamma$.

4. Efficiëntie-analyse:
   De auteurs tonen aan dat het evalueren van deze gesloten formules (via snelle algoritmen voor de digamma/polygammafuncties) vele malen sneller is dan het direct numeriek sommeren van de reeks. In een voorbeeld bleek de formule-gebaseerde methode ongeveer 30.000 keer sneller te zijn dan directe sommatie voor 280-cijferige precisie.

4. Resultaten

• Analytische Oplossingen: De paper levert bewijzen en formules voor een breed scala aan Euler-sommen die eerder alleen numeriek of in specifieke gevallen bekend waren.
• Voorbeelden:
  • Voor $\sum \frac{H_k}{(2k-1)^2}$ wordt een exacte uitdrukking gevonden die $\pi^2$, $\ln(2)$ en $\zeta(3)$ bevat.
  • Voor complexere gevallen zoals $\sum \frac{H_k}{(2k+1)^3(3k+1)^2}$ worden formules afgeleid die $\psi$-functies met argumenten als $1/3$ en $1/2$ bevatten.
• Validatie: Alle afgeleide formules zijn geverifieerd door de numerieke waarde van de som (via het directe algoritme) te vergelijken met de numerieke waarde van de gesloten formule (via de polygamma-berekeningen) tot op honderden decimalen.

5. Betekenis en Toekomstperspectief

• Wiskundige Elegantie en Efficiëntie: De formules vervangen oneindige sommaties door een eindig aantal termen met bekende speciale functies. Dit is niet alleen wiskundig elegant, maar ook cruciaal voor efficiënte numerieke berekeningen in toepassingen zoals kwantumveldtheorie en getaltheorie.
• Generalisatie: De methode is flexibel en kan worden uitgebreid naar andere soorten sommen, zoals die met $H_k^2$, alternerende tekens $(-1)^k$, of even harmonische getallen $H_{2k}$.
• Toepasbaarheid: De beschikbaarheid van deze algemene formules en het bijbehorende algoritme stelt onderzoekers in staat om complexe sommen in de natuurkunde en wiskunde sneller en nauwkeuriger te evalueren zonder afhankelijk te zijn van brute-force numerieke sommatie.

Samenvattend biedt dit paper een robuust theoretisch raamwerk en praktische tools om een grote klasse van Euler-sommen volledig analytisch op te lossen, waarbij de brug wordt geslagen tussen complexe analyse, symbolische berekening en hoogprecisie numerieke methoden.