Produit harmonique, sommation de

Ramanujan et fonctions zêta

d’Arakawa-Kaneko
par

Marc-Antoine COPPO

Mémoire d’Habilitation à Diriger des Recherches

Spécialité : Mathématiques

Présenté et soutenu le 21 juin 2016

au

Laboratoire Jean Alexandre Dieudonné

UMR CNRS 7351

devant le jury

M. Eric DELABAERE, Président du jury

M. Abdelmejid BAYAD, Rapporteur du jury
M. Driss ESSOUABRI, Rapporteur du jury
M. Dominique MANCHON, Rapporteur du jury
M. Bernard CANDELPERGHER, Membre du jury
M. Frédéric PATRAS, Membre du jury
M. Paul Thomas YOUNG, Membre du jury
À la mémoire d’André Cerezo (1945-2003)
Remerciements
Je suis infiniment reconnaissant à mes collègues et amis Bernard Candelpergher
et Frédéric Patras de l’aide précieuse qu’ils m’ont apportée. Sans l’inspiration et
la technicité de Bernard, sans l’ouverture d’esprit et le savoir-faire de Frédéric, je
n’aurais pu mener à bien ce projet d’habilitation.

Je remercie chaleureusement les rapporteurs, Driss Essouabri, Abdelmejid Bayad

et Dominique Manchon, pour l’attention qu’ils ont bien voulu accorder à ce mé-
moire, les appréciations qu’ils ont portées et les suggestions qu’ils ont exprimées.

Professeur dans l’une des plus anciennes universités américaines, Paul Thomas
Young me fait l’honneur de sa présence ; je l’en remercie vivement ainsi que pour
la collaboration fructueuse que nous avons entamée.

Je remercie “mon petit camarade” de l’Université de Nice (promotion 1981) Eric

Delabaere pour son amitié fidèle, son hospitalité et son soutien sans faille.

Je remercie Yasuo Ohno, Iaroslav Blagouchine, Clemens Berger, Vladimir Kos-

tov et Didier Clamond pour l’intérêt qu’ils ont manifesté pour mes travaux lors
de discussions ou à l’occasion de mes exposés, ainsi que Michel Merle pour les
encouragements et les conseils avisés qu’il m’a toujours prodigués.

Enfin, j’adresse mes remerciements sincères à Sorin Dumitrescu et à Clotilde Fer-

manian pour la confiance qu’ils m’ont témoignée.

5
Sommaire

Introduction 9

I Le produit harmonique dans l’espace des suites 13

1 Opérateurs dans l’espace des suites 13
1.1 Les opérateurs L et R . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.2 Les opérateurs D et S . . . . . . . . . . . . . . . . . . . . . . . . . 15

2 Le produit harmonique des suites 17

2.1 L’algèbre H = (E ∗ , o
n) . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.2 Expression explicite du produit harmonique . . . . . . . . . . . . . 19
2.3 Puissances harmoniques k-ièmes . . . . . . . . . . . . . . . . . . . . 20
2.4 Propriété d’harmonicité . . . . . . . . . . . . . . . . . . . . . . . . 20
2.5 Les sommes harmoniques . . . . . . . . . . . . . . . . . . . . . . . . 21

II La sommation de Ramanujan 23
3 L’opérateur D dans l’espace des fonctions 23
3.1 Transformation de Laplace-Borel . . . . . . . . . . . . . . . . . . . 23
3.2 L’opérateur D et le difféomorphisme Λ . . . . . . . . . . . . . . . . 24
3.3 Sommation de Ramanujan . . . . . . . . . . . . . . . . . . . . . . . 25

4 Le produit harmonique des fonctions 28

4.1 Le Λ-produit de convolution . . . . . . . . . . . . . . . . . . . . . . 28
4.2 Propriété d’harmonicité . . . . . . . . . . . . . . . . . . . . . . . . 29
4.3 La fonction zêta modifiée . . . . . . . . . . . . . . . . . . . . . . . . 30

5 Somme de Ramanujan des nombres hyperharmoniques 33

5.1 Nombres hyperharmoniques . . . . . . . . . . . . . . . . . . . . . . 33
5.2 Séries de Mascheroni décalées à droite . . . . . . . . . . . . . . . . . 33

III La fonction zêta d’Arakawa-Kaneko 37

6 Valeurs spéciales de la fonction ξk 38
6.1 Valeurs sur les entier positifs . . . . . . . . . . . . . . . . . . . . . . 38
6.2 Valeurs sur les entiers négatifs . . . . . . . . . . . . . . . . . . . . . 41

7
7 La fonction zêta d’Arakawa-Kaneko alternée 42
7.1 Valeurs spéciales de la fonction βk . . . . . . . . . . . . . . . . . . . 42

Références et publications 45

8
Introduction
Ce mémoire s’articule autour des trois principaux thèmes auxquels j’ai consa-
cré mes recherches au cours de ces 15 dernières années, thèmes qui sont assez
étroitement reliés entre-eux comme je m’attacherai à le montrer. Il s’agit du pro-
duit harmonique, du procédé de sommation de Ramanujan et de la fonction zêta
d’Arakawa-Kaneko qui ont fait l’objet des articles [1] à [10] de la liste de publica-
tions 1 . Ces travaux se situent au confluent de la combinatoire, de l’algèbre et de
la théorie des nombres.

Les principales propriétés du produit harmonique dans l’espace des suites à

valeurs complexes ont été exposées dans [3] et sont rappelées dans le chapitre I.
Ce produit est défini au moyen d’un opérateur de différence finie noté D (Défi-
nitions 5 et 8) qui est une transformation binomiale involutive laissant la suite
harmonique invariante (Théorème 1). Les suites invariantes par l’opérateur D sont
joliment caractérisées à l’aide du produit harmonique, conduisant à une élégante
reformulation algébrique d’un critère d’invariance de Sun (Théorème 4). Le pro-
duit harmonique admet une expression explicite qui le fait apparaître comme une
déformation du produit de convolution discret (Théorème 5).
Le produit harmonique possède de remarquables propriétés vis-à-vis des som-
mes harmoniques ; il permet notamment de généraliser les nombres harmoniques
de Rota (Théorème 7 et Exemple 12) et de donner une extension naturelle d’une
formule de Dilcher (Théorème 8 et Exemple 13).

Le produit harmonique réapparaît dans le chapitre II en relation avec la som-

mation de Ramanujan. Les propriétés fondamentales du procédé de sommation
de Ramanujan des séries divergentes ([Be], chapitre VI) ont été étudiées en détail
dans [10] ; récemment, cette étude a été profondément renouvelée et enrichie dans
[Ca]. Dans ce chapitre, on se place dans le cadre d’un espace de fonctions analy-
tiques dans le demi-plan { <(x) > 0} qui peuvent s’écrire comme des transformées
de Laplace ; dans ce cadre, la somme de Ramanujan d’une série peut être définie
à l’aide de l’opérateur D et des nombres de Bernoulli de deuxième espèce (Défi-
nition 17), et le produit harmonique peut être construit au moyen du produit de
convolution des fonctions (Définition 19).
La combinaison du produit harmonique et de la sommation de Ramanujan
permet d’introduire d’une manière algébrique une intéressante famille de fonctions
analytiques Fk de la variable complexe s, indexée par les entiers naturels, qu’on
1. Les références numérotées [11] à [13] se rattachent à une autre branche des mathématiques
et sont sans rapport avec le sujet de ce mémoire.

9
appelle les fonctions zêta modifiées. L’étude des propriétés de ces fonctions a été
initiée dans [4], puis poursuivie par Young (cf. [Y3] et [Y4]). Dans le cas où k = 0,
F0 (s) n’est autre que ζ(s) − s−1
1
où ζ(s) est la valeur en s de la classique fonction
zêta de Riemann. Les valeurs spéciales sur les entiers positifs des fonctions zêta
modifiées Fk s’expriment à l’aide de sommes infinies qui font intervenir les nombres
de Bernoulli de deuxième espèce et les polynômes de Bell modifiés évalués sur les
nombres harmoniques (Théorème 18). Lorsque k = 0, on retrouve une formule
classique due à Hermite 2 (Corollaire 3).
Enfin, on conclut ce chapitre II en présentant dans le paragraphe 5 une élégante
expression de la somme de Ramanujan de la série des nombres hyperharmoniques
de Conway ([CG] p. 258) qui a été démontrée dans [1]. Plus précisément, on montre
que la somme de Ramanujan des nombres hyperharmoniques d’ordre q est une
combinaison Q-linéaire de {1, γ, ln(2π)} ∪ {ζ 0 (−k)}k=1
q−1
, où γ désigne la constante
d’Euler (Corollaire 5 et Exemple 22).

Le chapitre III est consacré aux valeurs spéciales de la fonction zêta d’Arakawa-
Kaneko. Ces valeurs ont d’abord été étudiées par Arakawa et Kaneko ([AK]) puis
par Ohno qui a donné une remarquable expression des valeurs spéciales de la
fonction zêta d’Arakawa-Kaneko sur les entiers positifs en terme de “multiple zeta-
star values” ([O1], Théorème 2). Cette famille de fonctions analytiques ξk de la
variable complexe s, indexée par les entiers positifs, généralise la classique fonction
zêta de Riemann au travers de la relation ξ1 (s) = sζ(s + 1).
Dans [5], la fonction zêta d’Arakawa-Kaneko a été reconsidérée d’un point de
vue plus général en s’inspirant du modèle suivant lequel la fonction zêta d’Hurwitz
englobe la fonction zêta de Riemann comme un cas particulier (Définition 24). Pour
ce faire, on introduit une seconde variable x de partie réelle positive de telle sorte
que ξk (s, x) n’est autre que ξk (s) lorsque x = 1. Sur les valeurs entières et négatives
de s, la fonction ξk (s, x) interpole les polynômes de poly-Bernoulli (Théorème 23),
tandis que ses valeurs sur les entiers positifs peuvent s’exprimer au moyen de
sommes infinies qui font intervenir les polynômes de Bell modifiés évalués sur les
nombres harmoniques généralisés (Théorème 21) ; cette expression est similaire à
celle donnée pour la fonction zêta modifiée Fk au chapitre II précédent, l’explication
de cette analogie étant contenue dans la Remarque 12.
Ces remarquables propriétés font de la fonction zêta d’Arakawa-Kaneko géné-
ralisée ξk (s, x) (ainsi que de sa variante ”alternée” ξk∗ (s, x)) un puissant outil pour
l’étude des sommes d’Euler 3 (cas où x prend la valeur 1) et des sommes binomiales
2. Cette formule apparaît pour la première fois dans une lettre de Hermite à Pincherle datée
du 10 août 1900.
3. Introduites par Euler et Goldbach au milieu du 18ème siècle, les sommes d’Euler sont aussi
appelées valeurs zêta multiples ou nombres polyzêtas.

10
inverses 4 (cas où x prend la valeur 1/2), comme l’ont montré les belles identités
obtenues dans [2] (Exemples 23, 26, et 27). Les valeurs spéciales de la fonction
βk (s) = 2−s ξk∗ (s, 12 ) qui généralise la classique fonction bêta de Dirichlet (β0 n’est
autre que la fonction β) sont des périodes au sens de Kontsevich et Zagier ([KZ]).
Les sommes harmoniques et binomiales inverses qui interviennent dans le cadre
de cette étude sont apparues pour la première fois dans la littérature scientifique en
relation avec les diagrammes de Feynman ([DK]), mettant en évidence l’existence
de profondes connexions entre certaines branches de la théorie des nombres et de
la physique quantique ([Br], [Zh]).

4. Ces sommes binomiales inverses sont des cas particuliers de séries factorielles inverses consi-
dérées par Stirling dans son Methodus Differentialis (1730).

11
Chapitre I
Le produit harmonique dans
l’espace des suites
Tous les résultats énoncés dans ce chapitre ont été démontrés dans [3]. On
renvoie à cet article pour le détail des preuves.

1 Opérateurs dans l’espace des suites

On commence par introduire quelques opérateurs dans l’espace des suites à
valeurs complexes qui jouent un rôle crucial dans la construction du produit har-
monique, ainsi que leurs images respectives dans l’espace isomorphe des séries
formelles.
∗
Notation. Le C-espace vectoriel CN des suites a = (a(1), a(2), a(3), . . . , a(n), · · · )
à valeurs dans C est noté E ∗ .

Définition 1. Si C[[z]] désigne l’espace des séries formelles, on a un isomorphisme

naturel :
Φ : E ∗ −→ C[[z]]
défini par

zn z2 z3
Φ(a)(z) = a(n + 1) = a(1) + a(2)z + a(3) + a(4) + · · ·
X

n≥0 n! 2 6

Exemple 1. a) La suite δm définie pour tout m ≥ 0 et n ≥ 1 par


1 si n = m + 1
δm (n) =
0 sinon

vérifie la relation
zm
Φ(δm )(z) = .
m!
b) La suite 1 := (1, 1, 1, . . . ) vérifie Φ(1)(z) = ez .
c) La suite N := (1, 2, 3, . . . ) vérifie Φ(N )(z) = (1 + z)ez .
d) La suite géométrique de raison α ∈ C : αN −1 := (1, α, α2 , α3 , . . . ) vérifie la
relation
Φ(αN −1 )(z) = eαz .

13
 1 1
e) On appelle suite harmonique et on note H0 la suite (1, , , . . . ) des inverses
2 3
des entiers naturels. Elle vérifie la relation
1
Φ(H0 )(z) = (ez − 1) .
z
Définition 2. Les opérateurs sur E ∗ se transforment en opérateurs sur C[[z]]
via l’isomorphisme Φ. Plus précisément, si U désigne un opérateur sur E ∗ , il lui
correspond l’opérateur u sur C[[z]] défini par la relation
ΦU = uΦ ⇔ u = ΦU Φ−1
que l’on appellera l’image de U . On a donc le diagramme :
U
E∗ −−−→ E∗
x 
 
Φ−1 Φ
 y
u
C[[z]] −−−→ C[[z]]
L’image de l’opérateur I d’identité sur E ∗ sera notée Id.

1.1 Les opérateurs L et R

Définition 3. L’opérateur L de décalage à gauche (left) sur E ∗ est défini par
L(a)(n) = a(n + 1),
c’est à dire :
L
(a(1), a(2), a(3), . . . ) 7−→ (a(2), a(3), a(4), . . . ) .
L’image de L est l’opérateur de dérivation formelle ∂ :
zn z2
Φ(L(a))(z) = a(n + 2) = a(2) + a(3)z + a(4) + · · · = ∂Φ(a)(z).
X

n≥0 n! 2!

Définition 4. L’opérateur R de décalage à droite (right) sur E ∗ est défini par


a(n − 1) si n > 1
R(a)(n) = 
0 si n = 1,
c’est à dire :
R
(a(1), a(2), a(3), . . . ) 7−→ (0, a(1), a(2), a(3), . . . ) .
L’image de R est l’opérateur d’intégration formelle :
R

z n+1 z2 Z z
Φ(R(a))(z) = a(n + 1) = a(1)z + a(2) + · · · = Φ(a)(t) dt.
X

n≥0 (n + 1)! 2! 0

14
Remarque 1. Les opérateurs L et R ne sont pas inverses l’un de l’autre. On a la
relation LR = I, mais on notera que l’opérateur RL n’est pas l’identité :
RL
(a(1), a(2), a(3), . . . ) 7−→ (0, a(2), a(3), a(4), . . . ) .

1.2 Les opérateurs D et S

Définition 5. Soit V : E ∗ −→ C le morphisme d’évaluation défini par

V (a) = a(1).

L’image de V est l’application v : C[[z]] −→ C telle que v(Φ(a)) = Φ(a)(0).

L’opérateur de différence finie D : E ∗ −→ E ∗ est défini par
   
D(a)(n) = V (I − L)n−1 a = v (Id − ∂)n−1 Φ(a) ,

ce qui se traduit par

n
!
n
D(a)(n + 1) = (−1) j
a(j + 1) pour tout n ≥ 0.
X

j=0 j

On obtient ainsi pour de petites valeurs de n,

D(a)(1) = a(1),
D(a)(2) = a(1) − a(2),
D(a)(3) = a(1) − 2a(2) + a(3), etc.

On a la relation :
Φ(D(a))(z) = ez Φ(a)(−z) ,
ce qui signifie que l’image d de l’opérateur D vérifie pour tout f ∈ C[[z]],

d(f )(z) = ez f (−z) .

Théorème 1. L’opérateur D est un automorphisme auto-inverse qui laisse la suite

harmonique H0 invariante ; autrement dit, pour toute suite a ∈ E ∗ , on a

D(D(a)) = a

et

D(H0 ) = H0 .

Exemple 2. a) On a D(1) = δ0 et D(N ) = δ0 − δ1 .

15
b) Si αN −1 est la suite géométrique de raison α ∈ C, alors D(αN −1 ) = (1 − α)N −1 .
1
En particulier la suite ( )N −1 est invariante par D.
2
Définition 6. L’opérateur de sommation S : E ∗ −→ E ∗ est défini par
n
S(a) (n) =
X
a(j) .
j=1

L’opérateurR S est un automorphisme d’inverse I − R. L’image de S est l’opérateur

s = Id − d d.

Exemple 3. 1) S(δ0 ) = 1, S(1) = N .

1
2) S(αN −1 ) = (1 − αN ) pour α 6= 1. En particulier,
1−α
1
S((−1)N −1 ) = (1 + (−1)N −1 ) = (1, 0, 1, 0, . . . ) .
2
Théorème 2 (Relation entre les quatre opérateurs précédents). Pour toute suite
a ∈ E ∗ , on a les relations suivantes :

DL(a) = (I − L)D(a),
DS(a) = (I − R)D(a),
DR(a) = (I − S)D(a) .

1
Exemple 4. On rappelle que la suite harmonique H0 est définie par H0 (n) = .
n
Pour tout entier q ≥ 0, on considère la suite Hq des nombres hyperharmoniques
d’ordre q définie par
Hq = S q (H0 ) ,
ce qui se traduit récursivement par la relation
n
Hq (n) = Hq−1 (j) pour q ≥ 1.
X

j=1

En particulier, la suite H1 des nombres hyperharmoniques d’ordre 1 est la suite

des nombres harmoniques :
n
1
H1 (n) =
X
.
j=1 j

D’après la relation DR = (I − S)D et l’invariance de la suite H0 par D, on a

DR(H0 ) = (I − S)(H0 ) = H0 − H1 ,

16
ou, de manière équivalente,
D(H1 ) = H0 − R(H0 ) .
D’une manière plus explicite :
1

−

si n > 1
D(H1 )(n) = n(n − 1)
1

si n = 1.

2 Le produit harmonique des suites

On construit à présent le produit harmonique dans l’espace des suites au moyen
de l’opérateur D introduit précédemment.

2.1 L’algèbre H = (E ∗ , o
n)
Définition 7. Si a et b sont deux suites dans E ∗ , on note ab la suite définie par
(ab)(n) = a(n)b(n) .
On a en particulier : 1 a = a et δm a = a(m + 1)δm pour tout m ≥ 0. Muni de
ce produit, appelé produit de Hadamard des suites, l’espace E ∗ est une algèbre
commutative, associative et unitaire notée A. L’élément unité de A est la suite 1.
Définition 8. On définit le produit harmonique a o
n b de deux suites a et b dans
E par
∗

aon b := D(D(a)D(b)) .
Comme D = D−1 , on déduit immédiatement de la définition précédente les deux
relations fondamentales suivantes :

n b) = D(a)D(b) ,
D(a o
et
D(ab) = D(a) o
n D(b) .
Exemple 5. 1) 1 o
n a = a(1)1, car
n a) = D(1)D(a) = δ0 D(a) = D(a)(1)δ0 = a(1)δ0 = a(1)D(1).
D(1 o

2) N o
n a = a(2)1 + (a(1) − a(2))N , car
D(N )D(a) = (δ0 − δ1 )D(a)
= D(a)(1)δ0 − D(a)(2)δ1
= a(1)D(N ) + a(2)D(1 − N ).

17
3) αN −1 o
n β N −1 = (α + β − αβ)N −1 , car

n β N −1 ) = (1 − α)N −1 (1 − β)N −1
D(αN −1 o
= (1 − (α + β − αβ))N −1
= D((α + β − αβ)N −1 ) .

Théorème 3. L’espace (E ∗ , o
n) est une C-algèbre commutative, associative et uni-
taire notée H, isomorphe à l’algèbre A. L’élément unité dans H est la suite δ0 .

Corollaire 1. Une suite a est inversible dans H si et seulement si la suite D(a) est
inversible dans A (i.e. D(a)(n) 6= 0 pour tout n). Dans ce cas, l’inverse harmonique
de a est donné par la formule

1
!
a o(−1)
n
=D .
D(a)

Exemple 6. a)
(H0 )no(−1) = D(N ) = δ0 − δ1 ,
b)
N −1
α
o(−1) 
N −1 n

α = .
α−1
Remarque 2. On notera que l’algèbre H contient des diviseurs de zéro. On a par
exemple
1on δ1 = 0 .

On rappelle que la suite géométrique de raison 1/2 est invariante par D. Elle
peut s’écrire
1
 N −1
n δ0 .
o
2
Plus généralement, on donne la caractérisation suivante qui est une reformulation
algébrique du critère de Sun ([Su]).

Théorème 4 (Caractérisation des suites invariantes par D). Une suite a ∈ E ∗ est
invariante par D si et seulement si elle peut s’écrire sous la forme

1
 N −1
a= nb
o
2
où la suite b ∈ E ∗ est telle que b(2k) = 0 pour tout k ≥ 1.

18
Exemple 7. a) La suite harmonique H0 est invariante par D. Elle peut s’écrire
1
H0 = ( )N −1 o
nb
2
1 1
avec b = H0 on (−1)N −1 = (1, 0, , 0, , . . . ).
3 5
b) La suite
1 1 1 1
a = (δ0 + 1) = (1, , , , · · · )
2 2 2 2
est invariante par D. Elle peut s’écrire
1
a = ( )N −1 o
n (1, 0, 1, 0, · · · ) .
2

2.2 Expression explicite du produit harmonique

Le produit harmonique admet l’expression explicite suivante :
Théorème 5. Pour toutes suites a et b ∈ E ∗ et tout entier n ≥ 0,

(a o
n b)(n + 1) = Cni,j a(i + 1)b(j + 1)
X

0≤i≤n
0≤j≤n

où les nombres Cni,j sont définis par l’identité

(X + Y − XY )n = Cni,j X i Y j .
X

0≤i≤n
0≤j≤n

Remarque 3. En développant (X + Y − XY )n par la formule du binôme et en

identifiant le coefficient de X i Y j , ceci peut encore s’écrire
! !
n k
(a o
n b)(n + 1) = (−1) k−j
a(k + 1)b(n + 1 − j) (n ≥ 0) .
X

0≤j≤k≤n k j

Exemple 8. Pour de petites valeurs de n, on obtient ainsi

(a o
n b)(1) = a(1)b(1),
(a o
n b)(2) = a(2)b(1) + a(1)b(2) − a(2)b(2),
(a o
n b)(3) = a(3)b(1) + a(1)b(3) + 2a(2)b(2) − 2a(3)b(2) − 2a(2)b(3) + a(3)b(3),

etc.

19
2.3 Puissances harmoniques k-ièmes
Définition 9. Pour toute suite a ∈ E ∗ , on définit pour tout entier k ≥ 0, la
puissance harmonique k-ième de a notée anok par
ano0 = δ0 et ano(k+1) = anok o
n a.
Par récurrence sur k, on en déduit immédiatement la formule suivante :
 
anok = D(D(a) . . . D(a)) = D (D(a))k .
| {z }
k

En particulier, si a est une suite invariante par D, alors on a

anok = D(ak ) .
Exemple 9. a)

N si k est impair
N nok = D((δ0 − δ1 )k ) = 1 + (−1)k (1 − N ) = 
2−N si k est pair.

b)
k
(δ1 )nok = m!S2 (k, m)δm ,
X

m=0
où les S2 (k, m) sont les nombres de Stirling de deuxième espèce définis par
1 X m
!
m k
S2 (k, m) = (−1)m−j j .
m! j=0 j

c) !
k
D(N ) = (D(N ))
k ok
= (δ0 − δ1 ) ok
= (−1) j!S2 (i, j)δj .
i
X
n n

0≤j≤i≤k i

2.4 Propriété d’harmonicité

Notation. Pour p entier naturel, on pose (N )0 = 1, et pour p ≥ 1
(N )p = N (N + 1) · · · (N + p − 1) .
Théorème 6 (Propriété d’harmonicité). Pour toute suite a ∈ E ∗ et tout entier
p ≥ 0, on a la relation
(N )p
!
p! p!
na=
o S a ,
(N )p+1 (N )p+1 p!
En particulier, pour p = 0, on en déduit l’important corollaire :
n a = H0 S(a) .
H0 o

20
Exemple 10.
 
n H0 = H0 S(H0 ) = H0 H1 = D (H0 )2 ,
H0 o
 
(H0 )no3 = H0 o
n (H0 o
n H0 ) = H0 S(H0 H1 ) = D (H0 )3 .

Exemple 11. On a
1 1
na=
o S(N a)
N (N + 1) N (N + 1)
c’est à dire
1 1 n
( n a)(n) = ka(k) (n ≥ 1).
X
o
N (N + 1) n(n + 1) k=1

2.5 Les sommes harmoniques

On rappelle la remarquable propriété de la suite harmonique H0 vis-à-vis du
produit harmonique énoncée au paragraphe précédent :

n a = H0 S(a) ,
H0 o

ce qui peut se traduire par D (H0 D(a)) = H0 S(a) pour toute suite a.
Plus généralement, on introduit à présent une notion de somme harmonique
de la manière suivante :
Définition 10. Soit une suite a ∈ E ∗ , on définit pour tout entier naturel k, la
somme harmonique k-ième de a notée S (k) (a) par la formule

(H0 )nok o
n a = H0 S (k) (a) .

Théorème 7. Pour toute suite a ∈ E ∗ , on a S (1) (a) = S(a) et la relation de

récurrence : n
1 (k)
S (k+1) (a)(n) = S (a)(m) pour k ≥ 1.
X

m=1 m

Il en résulte que pour k ≥ 1,

1
S (k) (a)(n) = a(nk ) .
X

n≥n1 ≥···≥nk ≥1 n1 . . . nk−1

Exemple 12. Dans le cas où a est la suite harmonique H0 , les nombres S (k) (a)(n)
ne sont autres que les nombres harmoniques de Rota ([Ro]).
1
S (k) (H0 )(n) =
X
.
n≥n1 ≥···≥nk ≥1 n1 . . . nk

21
Pour les petites valeurs de k, on a ainsi :

S (0) (H0 )(n) = 1 ,

1 1
S (1) (H0 )(n) = 1 +
+ · · · + = H1 (n) ,
2 n
1 1 1 1 1 n
H1 (m)
S (H0 )(n) = 1 + (1 + ) + · · · + (1 + + · · · + ) =
(2)
X
.
2 2 n 2 n m=1 m

Théorème 8 (Formule de Dilcher étendue). Pour toute suite a ∈ E ∗ , et pour

k ≥ 1, on a l’identité
n
1
!
n
(k)
(a)(n) = (−1) m−1
X
S D(a)(m) .
m=1 m mk−1

Il en résulte que, pour tout entier naturel r,

1 n
n 1 (r)
!
= (−1)m−1 S (H0 )(m).
X X
r+1
n≥n1 ≥···≥nk ≥1 n1 . . . nk−1 nk m=1 m mk

Exemple 13. a) Pour r = 0, on retrouve la forme classique de la formule de

Dilcher ([D]) :

1 n
n 1
!
= (−1)m−1
X X
.
n≥n1 ≥···≥nk ≥1 n1 . . . nk m=1 m mk

b) Pour r = 1,

1 n
n 1
!
= (−1)m−1 H1 (m) .
X X
2
n≥n1 ≥···≥nk ≥1 n1 . . . nk−1 nk m=1 m mk

c) Pour r = 2,

1 n
1 X m
H1 (j)
!
m−1 n
= (−1)
X X
3
.
n≥n1 ≥···≥nk ≥1 n1 . . . nk−1 nk m=1 m mk j=1 j

d) Enfin, pour r = 3,

1 n
1 Xm
1X j
H1 (i)
!
m−1 n
= (−1)
X X
4 k
.
n≥n1 ≥···≥nk ≥1 n1 . . . nk−1 nk m=1 m m j=1 j i=1 i

22
Chapitre II
La sommation de Ramanujan
Tous les résultats énoncés dans ce chapitre ont été démontrés dans [1], [4], et
[9]. On renvoie à ces articles pour le détail des preuves.

3 L’opérateur D dans l’espace des fonctions

On commence par introduire un cadre analytique qui permet d’étendre les
opérateurs de différence finie D et de sommation S définis au chapitre I.

3.1 Transformation de Laplace-Borel

Définition 11. On considère l’espace vectoriel E des fonctions f ∈ C 1 (]0, +∞[) à
valeurs complexes vérifiant la propriété suivante :

Pour tout ε > 0, il existe Cε > 0 tel que |f (t)| ≤ Cε eεt pour tout t ∈ ]0, +∞[ .

Définition 12. Soit f une fonction dans l’espace E. La transformée de Laplace

L(f ) de f est définie par
Z +∞
L(f )(x) = e−xt f (t) dt pour <(x) > 0 .
0

Notation. Dans la suite, on note E = L(E) l’image de E par L.

Théorème 9. Si a est une fonction dans E, alors elle vérifie les propriétés sui-
vantes :
a) a est une fonction analytique dans le demi-plan { <(x) > 0 } ,
b) a(x) → 0 quand <(x) → +∞ ,
c) L : E → E est un isomorphisme.

Définition 13. Soit a ∈ E. La transformée de Borel de a est l’unique fonction

ab ∈ E telle que a = L(ab). On a les formules réciproques :

1 Z c+i∞ zt
ab(t) = e a(z) dz pour c > 0 et t > 0 ,
2iπ c−i∞
et Z +∞
a(x) = e−xt ab(t) dt pour <(x) > 0 .
0

23
Définition 14. Soient f et g deux fonctions dans E. Le produit de convolution
f ∗ g de f et g est la fonction définie pour tout t > 0 par
Z t
(f ∗ g)(t) = f (u)g(t − u) du .
0

Théorème 10. Pour tout f ∈ E et g ∈ E, alors f ∗ g ∈ E et

L(f ∗ g) = L(f ) L(g) .

Il en résulte que si a ∈ E et b ∈ E alors le produit ab ∈ E car ab = L(ab ∗ bb).

3.2 L’opérateur D et le difféomorphisme Λ

Définition 15. Soit a une fonction de E alors l’intégrale
Z +∞
e−t (1 − e−t )x−1 ab(t)dt
0

converge pour tout x vérifiant <(x) > 0. On appelle D(a) la fonction définie pour
tout x tel que <(x) > 0 par
Z +∞
D(a)(x) = e−t (1 − e−t )x−1 ab(t)dt .
0

Remarque 4.
Les valeurs de D(a) sur les entiers positifs peuvent être calculées directement sans
recourir à ab. Le développement de (1 − e−t )n par la formule du binôme conduit à
l’expression :
n
!
n
D(a)(n + 1) = (−1)j a(j + 1) pour tout entier n ≥ 0.
X

j=0 j

Autrement dit, l’opérateur D dans l’espace E des fonctions étend l’opérateur D

défini au Chapitre I dans l’espace E ∗ des suites.

Notation. On appelle Λ le C 1 -diffeomorphisme de R+ défini par

Λ(u) = − ln(1 − e−u ).

En particulier, il est important de noter que Λ est involutif :

Λ−1 = Λ .

24
Théorème 11. Soit a une fonction dans E. Alors la fonction D(a) ∈ E et, de
plus, verifie la relation
[ = ab(Λ) ,
D(a)
où ab(Λ) désigne ab ◦ Λ.
1 ts−1
Exemple 14. Soit a(x) = avec <(s) ≥ 1. Alors b (t) =
a . Il en résulte, par
xs Γ(s)
le changement de variable t = Λ(u), que

1 Λs−1
!
D( s ) = L ,
x Γ(s)

Remarque 5. Le Théorème 11 peut être visualisé par le diagramme suivant

D
E −−−→ E
 x
 
L−1 L
y 
Λ?
E −−−→ E
où Λ? (ab) = ab(Λ). Les propriétés algébriques de D sont résumées dans le théorème
suivant qui est l’analogue continu du Théorème 1.

Théorème 12. L’opérateur D est un automorphisme de E qui verifie D = D−1

1
et laisse invariante la fonction x 7→ .
x

3.3 Sommation de Ramanujan

Définition 16. La suite des nombres de Bernoulli de seconde espèce (bn ) ([J],
[Y1]), encore appelés coefficients de Gregory 5 ([Bl]), est définie par la fonction
génératrice
z
= bn z n .
X
ln(1 + z) n≥0

Exemple 15. Les bn sont des nombres rationnels qui peuvent se calculer au moyen
de la relation de récurrence
n−1
(−1)k bk
b0 = 1 et =0 pour n ≥ 2.
X

k=0 n − k

5. Cette famille de nombres apparaît en effet pour la première fois dans une lettre de James
Gregory datant de 1670.

25
Une représentation intégrale des nombres bn est donnée par
Z 1
n!bn = x(x − 1)(x − 2) . . . (x − n + 1) dx pour n ≥ 1.
0

On a ainsi pour les petites valeurs de n,

1 1 1 19 3
b1 = , b 2 = − , b 3 = , b 4 = − , b5 = , etc.
2 12 24 720 160
On notera que |bn | = (−1)n−1 bn pour tout n ≥ 1.
Théorème 13. Soit a une fonction dans E. La serie
Z +∞
e−t (1 − e−t )n−1 ab(t) dt
X
|bn |
n≥1 0

converge et
∞ Z +∞ Z +∞
1 1
e−t (1 − e−t )n−1 ab(t)dt = e−t ( − ) ab(t) dt .
X
|bn |
n=1 0 0 1−e −t t
Remarque 6. Si a ∈ E, alors la série a(n) peut s’écrire
P
n≥1

X Z +∞
a(n) = e−nt ab(t) dt ,
X

n≥1 n≥1 0

et une permutation formelle de et conduirait à écrire

P R +∞
n≥1 0
Z +∞
e−t
a(n) = ab(t) dt.
X

n≥1 0 1 − e−t
Cependant, cette dernière intégrale peut diverger en 0. On peut la renormaliser
en retirant la singularité en 0 et ceci peut être fait simplement en soustrayant la
partie polaire 1t de 1−e1 −t . Ceci conduit à la définition suivante :
Définition 17. Soit a une fonction dans E = L(E). La somme de Ramanujan de
la série n≥1 a(n) est définie par
P

R Z +∞
1 1 ∞
a(n) = e−t ( ) b (t)dt = |bn |D(a) (n) .
X X
− a
n≥1 0 1 − e−t t n=1

1 ts−1
Exemple 16. Soit a(x) = avec <(s) > 1. Alors, a ∈ E et b (t) =
a . Par
xs Γ(s)
conséquent
R
1 1
= ζ(s) −
X

n≥1 n
s s−1
où ζ(s) est la valeur en s de la fonction zêta de Riemann.

26
Remarque 7. a) La définition de la somme de Ramanujan donnée dans ce cha-
pitre coïncide avec celle donnée dans un cadre plus général dans [10] et [Ca] ;
elle possède donc les mêmes propriétés fondamentales. En particulier, il résulte
de la définition que l’application a 7→ Rn≥1 a(n) est linéaire.
P

b) Soit m un entier positif et a ∈ E. La somme de Ramanujan de la série translatée

n≥1 a(n + m) s’exprime par la formule :
P

R Z +∞
1 1
a(n + m) = e−t e−mt ( − ) ab(t)dt .
X

n≥1 0 1−e −t t
Cette somme ne vérifie cependant pas la propriété de décalage usuelle :
R R m
a(n + m) =
X X X
a(n) − a(j)
n≥1 n≥1 j=1

mais seulement la relation inhabituelle :

R R m Z m+1
a(n + m) = a(j) +
X X X
a(n) − a(x) dx
n≥1 n≥1 j=1 1

Dans le Chapitre V de [Ca] est développé un formalisme algébrique adéquat

qui explique pourquoi cette relation est “naturelle”.
c) La formule
R ∞ ∞ n
!
n
a(n) = |bn+1 |D(a) (n + 1) = (−1) j
a(j + 1)
X X X X
|bn+1 |
n≥1 n=0 n=0 j=0 j
est vraie dans un cadre plus général (cf. [Ca], Théorème 15).
Exemple 17 (Constantes de Stieltjes). Les constantes de Stieltjes γk sont définies
par le développement en série de Laurent de ζ :
1 ∞
(−1)k
ζ(s) = +γ+ γk (s − 1)k (s 6= 1) .
X
s−1 k=1 k!
Ces constantes peuvent se représenter comme des sommes de Ramanujan (cf. [9]),
plus précisément, on a l’expression suivante :
R
lnk (n) X ∞ n
n lnk (j + 1)
!
γk = = |bn+1 | (−1)j
X X

n≥1 n n=0 j=0 j j+1

∞ n X m
!
|bn+1 | X m
= (−1)j lnk (j + 1) (par la propriété d’harmonicité) .
X

n=0 n + 1 m=0 j=0 j

En particulier, la première constante de Stieltjes γ1 peut s’écrire
R
ln (n) 3 π2 X ∞
|bn | n−1
" #
X |bn−k+1 bk |
γ1 = = − + + (H1 (n) − H1 (k)) .
X

n≥1 n 2 6 n=2 n
2
k=1 n − k + 1

27
4 Le produit harmonique des fonctions
On donne à présent une construction du produit harmonique dans ce cadre
analytique au moyen du produit de convolution.

4.1 Le Λ-produit de convolution

Définition 18. Si a et b sont deux fonctions dans E, alors le Λ-produit de convo-
lution ab ~ bb de ab et bb est défini par

ab ~ bb = Λ? (Λ? (ab) ∗ Λ? (bb)) ,

ou de manière équivalente (puisque Λ? = (Λ? )−1 ) par

(ab ~ bb)(Λ) = ab(Λ) ∗ bb(Λ) .

Le Λ-produit de convolution hérite des propriétés algébriques du produit de convo-

lution ordinaire : il est bilinéaire, commutatif et associatif.

Définition 19. Soient a et b deux fonctions dans E. Le produit harmonique a o

nb
de a et b est défini par
aon b = L(ab ~ bb) ∈ E .
Cette construction peut être synthétisée dans le diagramme suivant :

(a, b) −−−→ (ab, bb) −−−→ (ab(Λ), bb(Λ))

  
  
  
y y y

n b ←−−− ab ~ bb ←−−− ab(Λ) ∗ bb(Λ)

ao

Le produit harmonique dans E hérite des propriétés de bilinéarité, de commutati-

vité et d’associativité du Λ-produit de convolution.

Théorème 14. Si a et b sont deux fonctions dans E alors on a les relations

fondamentales
D(a on b) = D(a) D(b) ,
et
D(ab) = D(a) o
n D(b) .

Du point de vue de la sommation de Ramanujan, les deux relations précédentes

s’interprètent de la façon suivante :

28
Corollaire 2. Soient a et b dans E, on a les identités :
R ∞
(a o
n b)(n) =
X X
|bn |D(a)(n) D(b)(n) .
n≥1 n=1

et
R ∞
(ab)(n) =
X X
|bn |(D(a) o
n D(b))(n) .
n≥1 n=1

Remarque 8. Les valeurs du produit a o n b sur les entiers positifs peuvent être
évaluées sans recourir à ab ni à bb. Par des transformations élémentaires, on peut en
effet montrer que
Z +∞ Z +∞
(a o
n b)(n + 1) = (e−t−s )(e−t + e−s − e−t e−s )n ab(t)bb(s)dtds .
0 0

Par conséquent, si les nombres Cni,j sont définis par

(X + Y − XY )n = Cni,j X i Y j ,
X

0≤i≤n
0≤j≤n

on a l’expression explicite

(a o
n b)(n + 1) = Cni,j a(i + 1)b(j + 1) ,
X

0≤i≤n
0≤j≤n

Autrement dit, le produit harmonique dans l’espace E des fonctions étend natu-
rellement le produit harmonique défini au Chapitre I dans l’espace E ∗ des suites.

4.2 Propriété d’harmonicité

On énonce à présent une propriété d’harmonicité dans l’espace E des fonctions
analogue à celle déjà vue dans l’espace des suites (cf. Théorème 6).

Théorème 15. Soit a ∈ E. Alors

1 S(a)(x)
na=
o ,
x x
où S(a) désigne la fonction définie pour <(x) > 0 par
Z +∞
e−t
S(a)(x) = (1 − e−xt )ab(t)dt .
0 1 − e−t

29
Remarque 9. L’opérateur S étend l’opérateur de sommation introduit au Cha-
pitre I. Pour tout entier n ≥ 1, on a en effet
n
S(a)(n) =
X
a(j) .
j=1

Exemple 18.
1 1 ψ(x + 1) + γ
n = L(Λ) =
o ,
x x x
où ψ désigne la fonction digamma (dérivée logarithmique de la fonction Γ).

4.3 La fonction zêta modifiée

Définition 20. Pour tout entier k ≥ 1, on considère la puissance harmonique
1 1
k-ième de la fonction x 7→ notée (abusivement) .
x x
1 ok 1 1 1
 n
= o
n o n ... o
n (k = 1, 2, 3, · · · ) ,
x |x x {z x}
k

Pour tout entier k ≥ 0 et <(s) ≥ 1, on définit la fonction zêta modifiée d’ordre k

par la formule suivante :
P
 R
 1
n≥1 
ns
si k = 0,
Fk (s) :=  nok

(n) si k ≥ 1
PR 1 1

 n≥1 x
o
n xs

Théorème 16. Pour tout entier k ≥ 0,

∞
|bn | 1
 
Fk (s) = (n)
X
k
D s
n=1 n x
1
où désigne la fonction x 7→ x−s . En particulier,
xs
∞
1 1
 
F0 (s) = (n) = ζ(s) −
X
|bn |D s .
n=1 x s−1

Théorème 17. La fonction Fk se prolonge analytiquement dans C en une fonction

entière. Les valeurs aux entiers négatifs de la fonction zêta modifiée sont données
par :
n m
!
|bm+1 | X m
Fk (−n) = (−1)j (j + 1)n . (n = 0, 1, 2, . . . )
X

m=0 (m + 1) k
j=0 j

30
 En particulier, pour k = 0, on déduit la relation
1 − Bn+1 n m
!
m
= |bm+1 | (−1)j (j + 1)n
X X
n+1 m=0 j=0 j
où les Bn sont les nombres de Bernoulli ([AIK]).
Définition 21 (Polynômes de Bell modifiés). La suite des polynômes de Bell
modifiés (Pn ) est définie par la fonction génératrice
∞ ∞
zk
!
exp = Pn (x1 , · · · , xn ) z n
X X
xk
k=1 k n=0

ou par la relation de récurrence équivalente :

n
P0 = 1, nPn (x1 , . . . , xn ) = xj Pn−j (x1 , . . . , xn−j ) (n ≥ 1) ,
X

j=1

ou bien encore par la représentation explicite :

1
k 1  k2 kn
x1 x2 xn
 
Pn (x1 , ..., xn ) =
X
...
k1 +2k2 +···+nkn =n k1 !k2 ! . . . kn ! 1 2 n
Pour de petites valeurs de n, on obtient ainsi les expressions suivantes :
P0 = 1 ,
P1 (x1 ) = x1 ,
1 1
P2 (x1 , x2 ) = x21 + x2 ,
2 2
1 3 1 1
P3 (x1 , x2 , x3 ) = x1 + x1 x2 + x3 ,
6 2 3
1 4 1 2 1 1 1
P4 (x1 , x2 , x3 , x4 ) = x1 + x1 x2 + x22 + x1 x3 + x4 .
24 4 8 3 4

Théorème 18 (valeurs spéciales de Fk sur les entiers positifs). Pour tout entier
m ≥ 0, on a
1 Pm (Hn(1) , Hn(2) , . . . , Hn(m) )
 
D m+1 (n) =
x n
d’où il résulte que
∞
|bn |
Fk (m + 1) = Pm (Hn(1) , Hn(2) , . . . , Hn(m) ) ,
X
k+1
n=1 n

où bn désigne le n-ième nombre de Bernoulli de seconde espèce, Pm le m-ième

n
1
polynôme de Bell modifié et Hn(m) =
X
m
.
j=1 j

31
Remarque 10. Avec les notations du Chapitre I, on a l’identité Hn(1) = H1 (n).
Dans le cas où k = 0, on obtient une reformulation de la formule d’Hermite :
Corollaire 3 (formule d’Hermite revisitée). Pour m ≥ 1,
1 ∞
|bn |
F0 (m + 1) = ζ(m + 1) − = Pm (Hn(1) , Hn(2) , . . . , Hn(m) ) ,
X
m n=1 n
Corollaire 4 (formule de dualité). Pour tout entier k ≥ 0,
∞
|bn | R
Pk (Hn(1) , Hn(2) , . . . , Hn(k) ) X R
1 1
Fk (1) = = =
X X X
k+1
.
n=1 n n≥1 n n≥1 n n≥n1 ≥···≥nk ≥1 n1 . . . nk

Exemple 19.
∞
|bn |
F0 (1) = =γ,
X

n=1 n
∞
|bn |Hn(1)
F0 (2) = = ζ(2) − 1 ,
X

n=1 n
1X ∞
|bn |(Hn(1) )2 1 X ∞
|bn |Hn(2) 1
F0 (3) = + = ζ(3) − ,
2 n=1 n 2 n=1 n 2
∞ R
|bn | X Hn(1)
F1 (1) = =
X
2
,
n=1 n n≥1 n
∞
|bn | X R
1 Xn (1)
Hm 1X R
(Hn(1) )2 1 X R
Hn(2)
F2 (1) = = = +
X
.
n=1 n
3
n≥1 n m=1 m 2 n≥1 n 2 n≥1 n
Les valeurs spéciales de la fonction F1 sur les entiers positifs admettent une
expression particulière :
Théorème 19. Pour tout entier q ≥ 2,
∞ X 1 X
Hn(1) q−1 ∞
1 ∞
ln(n + 1)
F1 (q) = γζ(q) + ζ(q + 1) − +
X X
− .
k=1 k n=1 (n + 1) n
q k q−k nq
n=1 n n=1

Exemple 20.
∞
|bn |Hn(1) π2 ∞
ln (n + 1)
F1 (2) = = 1 +
X X
γ − ζ(3) − ,
n=1 n 2 6 n=1 n2
π4 π2 X∞
ln (n + 1)
F1 (3) = γζ(3) − − + ,
360 12 n=1 n3
π4 π2 2 π2 1 X∞
ln (n + 1)
F1 (4) = γ − 2ζ(5) + ζ(3) − ζ(3) + − + .
90 6 3 18 2 n=1 n4

32
5 Somme de Ramanujan des nombres hyperhar-
moniques
5.1 Nombres hyperharmoniques
Définition 22. On rappelle que H0 désigne la suite définie par H0 (n) = n1 . Pour
tout entier q ≥ 0, on considère la suite Hq des nombres hyperharmoniques d’ordre
q ([CG], [DB]) définie par
Hq = S q (H0 )
où S est l’opérateur de sommation défini au chapitre I. Ceci se traduit par la
relation de récurrence : n
Hq (n) = Hq−1 (j) .
X

j=1

Les nombres hyperharmoniques d’ordre 1 sont simplement les nombres harmo-

niques
n
1
H1 (n) =
X
,
j=1 j

et pour q ≥ 2,

n+q−1
!
Hq (n) = (H1 (n + q − 1) − H1 (q − 1)) .
q−1

En particulier,
n
H2 (n) = H1 (j) = (n + 1)H1 (n) − n .
X

j=1

5.2 Séries de Mascheroni décalées à droite

|bn |
Définition 23. Si r désigne entier naturel, on considère la série décalée
P
n≥1 n
r fois à droite, c’est à dire la série
∞
|bn |
σr :=
X
.
n=r+1 n − r

Théorème 20. Pour tout entier naturel r, les sommes σr admettent l’expression
suivante :
1 1 1
σ0 = γ , σ1 = −ζ 0 (0) − γ − = (ln(2π) − γ − 1) ,
2 2 2

33
et pour r ≥ 2,

1 r−1
σr = − S1 (r − 1, k)ζ 0 (−k) + (−1)r br γ + tr
X
(r − 1)! k=1
avec
1 r
Bk
tr = − S1 (r − 1, k − 1) 2 ,
X
(r − 1)! k=1 k

où les S1 (n, k) sont les nombres de Stirling de première espèce (sans signe) définis
comme les coefficients de la factorielle ascendante
n
(x)n := x(x + 1) . . . (x + n − 1) = S1 (n, k)xk et S1 (n, k) = 0 pour n < k.
X

k=0

Inversement, la formule pour σr précédente admet la formule duale :

k+1
Bk+1 Bk+1
ζ 0 (−k) = (−1)k−r (r − 1)!S2 (k, r − 1)σr − pour k ≥ 1 ,
X
γ−
r=2 k+1 (k + 1)2
où les S2 (k, m) sont les nombres de Stirling de deuxième espèce.
Exemple 21.
1 1
σ2 = −ζ 0 (−1) − γ− ,
12 24
1 1 1 1
σ3 = − ζ 0 (−1) − ζ 0 (−2) − γ − ,
2 2 24 48
1 1 1 19 13
σ4 = − ζ 0 (−1) − ζ 0 (−2) − ζ 0 (−3) − γ− ,
3 2 6 720 960
1 11 1 1 3 19
σ5 = − ζ 0 (−1) − ζ 0 (−2) − ζ 0 (−3) − ζ 0 (−4) − γ− .
4 24 4 24 160 1920
Et inversement,
1 1
ζ 0 (−1) = −σ2 − γ− ,
12 24
ζ 0 (−2) = σ2 − 2σ3 ,
1 1
ζ 0 (−3) = −σ2 + 6σ3 − 6σ4 + γ+ ,
120 480
ζ 0 (−4) = σ2 − 14σ3 + 36σ4 − 24σ5 .
L’évaluation des sommes σr permet d’exprimer la somme de Ramanujan de la
série des nombres hyperharmoniques d’ordre q comme une combinaison Q-linéaire
de {1, γ, ln(2π)} ∪ {ζ 0 (−k)}q−1
k=1 . Plus précisément, on a la formule suivante.

34
Corollaire 5.
R q !
q
Hq (n) = (−1)r
X X
σr
n≥1 r=0 r
q " ! #
q q q
= − ln(2π) + γ br +
X
2 r=0 r 2
1
!
q
(−1) r
S1 (r − 1, k)ζ 0 (−k)
X
−
1≤k<r≤q r (r − 1)!
1
!
q Bk
(−1)r S1 (r − 1, k − 1) 2 .
X
−
2≤k≤r≤q r (r − 1)! k

Exemple 22.
R
1 3 1
H1 (n) = σ0 − σ1 = − ln(2π) + γ + ,
X

n≥1 2 2 2
R
23 23
H2 (n) = − ln(2π) + γ − ζ 0 (−1) +
X
,
n≥1 12 24
R
3 55 5 1 67
H3 (n) = − ln(2π) + γ − ζ 0 (−1) + ζ 0 (−2) +
X
,
n≥1 2 24 2 2 48
R
1901 13 3 1 1747
H4 (n) = −2 ln(2π) + γ − ζ 0 (−1) + ζ 0 (−2) − ζ 0 (−3) +
X
.
n≥1 720 3 2 6 960

Remarque 11. L’introduction des polynômes de Bernoulli de deuxième espèce

bn (x) qui sont définis par la fonction génératrice
∞
z
(1 + z)x = bn (x)z n
X
ln(1 + z) n=0

permet d’exprimer plus simplement le coefficient de γ dans l’expression précédente

de la somme de Ramanujan de la série des nombres hyperharmoniques d’ordre q.
On a q !
q
coefficient de γ = br = bq (q) .
X

r=0 r

35
Chapitre III
La fonction zêta
d’Arakawa-Kaneko
Tous les résultats énoncés dans ce chapitre ont été démontrés dans [2] et [5]. On
renvoie à ces deux articles pour le détail des preuves. On commence par introduire
la fonction zêta d’Arakawa-Kaneko généralisée ξk (s, x) puis sa variante alternée
ξk∗ (s, x). On s’intéresse tout particulièrement à leurs valeurs spéciales sur les entiers
qui sont des périodes au sens de Kontsevich et Zagier ([KZ]).
Définition 24. Pour tout entier k ≥ 1, la fonction zêta d’Arakawa-Kaneko ξk (s, x)
est définie pour <(s) > 0 et <(x) > 0 par la transformée de Mellin normalisée
1 Z +∞ −xt Lik (1 − e−t ) s−1
ξk (s, x) = e t dt
Γ(s) 0 1 − e−t
avec ∞
zn
Lik (z) =
X
k
.
n=1 n
En particulier, comme Li1 (1 − e−t ) = t, on notera que ξ1 (s, x) n’est autre que
sζ(s + 1, x) où ζ(s, x) est la fonction zêta d’Hurwitz. Dans la suite, on s’intéresse
tout particulièrement aux valeurs x = 1 et x = 1/2. Pour cela, on pose
1
ξk (s) := ξk (s, 1) et αk (s) := 2−s ξk (s, ) .
2
Pour k = 1, il existe une relation simple entre les valeurs des fonctions ζ, ξk et αk :
1
ξ1 (s) = s ζ(s + 1) et α1 (s) = 2−s sζ(s + 1, ) = (2 − 2−s )sζ(s + 1) .
2
Remarque 12 (Lien avec la fonction zêta modifiée). Pour <(s) ≥ 1, on peut
donner au moyen de l’opérateur D et du produit harmonique une définition plus
algébrique de la fonction ξk d’Arakawa-Kaneko comme suit :
∞
1 ok 1
 n !
ξk (s) = (n)
X
D o
n s
n=1 x x
∞
1 1
 
= (n)
X
k
D
n=1 n xs
1
où désigne la fonction x 7→ x−s . En particulier, on remarque l’analogie formelle
xs
avec l’expression de la fonction Fk vue au chapitre II (cf. Théorème 16) :
∞
|bn | 1
 
Fk (s) = (n) .
X
k
D
n=1 n xs

37
6 Valeurs spéciales de la fonction ξk
6.1 Valeurs sur les entier positifs
Théorème 21 (valeurs aux entiers positifs). Pour tout entier m ≥ 0 et <(x) > 0,
∞
n!
ξk (m + 1, x) = Pm (h(1)
n (x), . . . , hn (x))
(m)
X

n=0 (n + 1)k x(x + 1) . . . (x + n)

où Pm (x1 , . . . , xm ) désigne le m-ième polynôme de Bell modifié évalué sur les

nombres harmoniques généralisés
n
1
n (x) =
h(m)
X
.
j=0 (j + x)
m

En spécialisant cette relation en x = 1 et x = 1/2, on déduit le corollaire

suivant :

Corollaire 6. Pour tout entier m ≥ 0,

∞
1
ξk (m + 1) = Pm (Hn(1) , . . . , Hn(m) ) ,
X

n=1 nk+1

et
∞
22n−1
αk (m + 1) = Pm (On(1) , . . . , On(m) ) ,
X
 
2n
n=1 n
nk+1
avec
(m)
n
1
On(m) = 2−m hn−1 (1/2) =
X
,
j=1 (2j − 1)
m

et

(m)
n
1
Hn(m) = hn−1 (1) =
X
m
.
j=1 j

En particulier,

ξk (1) = ζ(k + 1) ,
∞
Hn(1) 1 1X k
ξk (2) = = (k + 3)ζ(k + 2) ζ(j)ζ(k + 2 − j) ,
X
−
n=1 n
k+1 2 2 j=2

38
et pour tout entier m ≥ 0,
∞
1
ξ1 (m + 1) = (m + 1)ζ(m + 2) = P (Hn(1) , . . . , Hn(m) ) ,
X
n 2 m
n=1
∞
22n−1
α1 (m + 1) = (2 − 2−m−1 )(m + 1)ζ(m + 2) = Pm (On(1) , . . . , On(m) ) .
X
 
2n
n=1 n
n2

Remarque 13. Pour k = 0, la fonction ξk n’est pas définie car l’intégrale

1 Z +∞ −t Li0 (1 − e−t ) s−1 1 Z +∞ s−1

e t dt = t dt
Γ(s) 0 1 − e−t Γ(s) 0

diverge pour tout s. Cependant, au sens de la sommation de Ramanujan, d’après

la première formule du Corollaire 6, on déduit du Corollaire 4 l’identité
R
Pm (Hn(1) , Hn(2) , . . . , Hn(m) ) X ∞
|bn |
ξ0 (m + 1) = = pour tout entier naturel m.
X
n m+1
n≥1 n=1 n

On déduit également du Corollaire 6 une intéressante réécriture de la formule

sommatoire d’Ohno ([O2], Théorème 8).

Corollaire 7. Pour tout entier k ≥ 2,

k−2 k−2 ∞
1
ξk−m−1 (m + 1) = Pm (Hn(1) , . . . , Hn(m) ) = (2 − 22−k )(k − 1)ζ(k) .
X X X

m=0 m=0 n=1 nk−m

En particulier, il en résulte que

∞
(Hn(1) )2 X ∞
Hn(2)
+ = 2ζ(2)ζ(3) − ζ(5) .
X

n=1 n3 n=1 n
3

Théorème 22. Pour tout entier k ≥ 1,

∞
22n 1 k
(ln 2)k−j
2αk (1) = = 2k−1 (−1)j−1
X X
Lj ,
(j − 1)!(k − j)!
 
n=1
2n nk+1 j=1
n

avec Z π
u
 
Lk = u lnk−1 2 sin du .
0 2
Remarque 14. Avec les notations introduites par Lewin (cf. [L]), l’intégrale Lk
(1)
n’est autre que −Lsk+1 (π).

39
Exemple 23.
∞
Hn(1)
ξ1 (2) = = 2ζ(3) ,
X
2
n=1 n
∞
22n−1 On(1) 7
α1 (2) = = ζ(3) ,
X
2
 
n=1
2n n2
n
∞
Hn(1) 5 π4
ξ2 (2) = = ζ(4) =
X
,
n=1 n
3 4 72
∞
2
2n−1
12
π 7
α2 (1) = = ln 2 − ζ(3)
X
2 4
 
2n n3
n=1 n
∞
2
2n−1
On(1) π4
α2 (2) = = 7ζ(3) ln 2 − 8G(1) ,
X
−
32
 
n=1
2n n3
n
∞
2
2n−1
1 π2 7 π4
α3 (1) = = (ln 2) 2
ln 2 + + 4G(1)
X
− ζ(3)
2 2 96
 
n=1
2n n4
n

∞
On(1)
où G(1) = désigne la constante de Ramanujan ([Be] p. 257, [Si]).
X

n=1 (2n)
3

Remarque 15. Comme

7 1 π3 X (−1)n
ζ(3) = et =
X
,
8 n≥0 (2n + 1)
3 32 n≥0 (2n + 1)3

la formule
7 π4 1X∞
22n On
G(1) = ζ(3) ln 2 − −
8 256 2 n=1 2n (2n)3
 
n

déduite de l’évaluation de α2 (2), peut encore se réécrire

∞
1 ∞
(−1)n 1X ∞
22n On
G(1) = 2 ln(2) − (π/8 + ln(2))
X X
− ,
(4n + 1) (2n + 1) 2 2n (2n)3
 
3 3
n=0 n=0 n=1 n

qui est l’expression exacte la plus proche de la formule

πX ∞
(−1)n π X ∞
1
G(1) = − √
4 n=0 (4n + 1) 3
3 3 n=0 (2n + 1)3

donnée par Ramanujan dans le chapitre IX de ses Notebooks et qui s’est avérée
erronée (cf. [Be] p. 257, [Si]).

40
6.2 Valeurs sur les entiers négatifs
Théorème 23 (valeurs aux entiers négatifs). Pour tout entier k ≥ 1 et pour
<(x) > 0, la fonction s 7→ ξk (s, x) se prolonge analytiquement dans C en une
fonction entière. Les valeurs aux entiers négatifs de la fonction d’Arakawa-Kaneko
sont données par :

ξk (−n, x) = (−1)n Bn(k) (x) (n = 0, 1, 2, . . . )

où les Bn(k) (x) sont les polynômes de poly-Bernoulli définis par la fonction généra-
trice :
−xt Lik (1 − e )
−t ∞
tn
= Bn (x) .
(k)
X
e
1 − e−t n=0 n!

Remarque 16. Les polynômes Bn(k) (x) sont des polynômes de degré n en x. Pour
k = 1, on retrouve (au signe près) les polynômes de Bernoulli classiques ([J]). Les
nombres Bn(k) := Bn(k) (0) sont les nombres de poly-Bernoulli introduits par Kaneko
([AIK]). On a l’expression suivante :
n
!
n (k) n−j
Bn(k) (x) = (−1) n−j
X
Bj x .
j=0 j

Les polynômes de poly-Bernoulli admettent l’expression explicite

n
1 m
!
m
(−1) n
Bn(k) (x) = (−1)j (x + j)n .
X X

m=0 (m + 1) j=0
k j

Exemple 24. Pour de petites valeurs de k et n, on obtient ainsi

(2) 1
ξ2 (−1, x) = −B1 (x) = x − ,
4
(2) 1 1
ξ2 (−2, x) = B2 (x) = x2 − x − ,
2 36
(2) 3 1 1
ξ2 (−3, x) = −B3 (x) = x3 − x2 − x + ,
4 12 24
(3) 1
ξ3 (−1, x) = −B1 (x) = x − ,
8
(3) 1 11
ξ3 (−2, x) = B2 (x) = x2 − x − ,
4 216
(3) 3 11 1
ξ3 (−3, x) = −B3 (x) = x3 − x2 − x + .
8 72 288

41
7 La fonction zêta d’Arakawa-Kaneko alternée
Définition 25. On considère à présent la fonction zêta d’Arakawa-Kaneko alter-
née ξk∗ (s, x) définie pour <(s) > 0, <(x) > 0 et k ≥ 0 par

1 Z +∞ e−xt 1 − e−t s−1

ξk∗ (s, x) = Lik ( ) t dt
Γ(s) 0 1 − e−t 2
On s’intéresse tout particulièrement au cas où x = 1/2. Pour cela, on définit
1
βk (s) := 2−s ξk∗ (s, ) .
2
En particulier, pour k = 0,
∞
(−1)n−1
β0 (s) = = β(s)
X

n=1 (2n − 1)
s

où β désigne la fonction bêta de Dirichlet ([RZ]). Autrement dit, la fonction βk est

une généralisation assez naturelle de la fonction bêta de Dirichlet .

7.1 Valeurs spéciales de la fonction βk

Théorème 24. Pour tout entier m ≥ 0,
∞
2n−1
βk (m + 1) = Pm (On(1) , . . . , On(m) ) (k ≥ 0) .
X
 
2n
n=1 n
nk+1
Exemple 25. Comme β0 (s) = β(s), on en déduit que pour tout entier m ≥ 1,
∞
2n
2β(2m) = P2m−1 (On(1) , . . . , On(2m−1) ) ,
X
 
2n
n=1 n
n
et
∞
2n
 2m+1
E2m π
2β(2m + 1) = P2m (On(1) , On(2) , . . . , On(2m) ) = (−1) m
X
(2m)! 2
 
2n
n=1 n
n
où les E2m sont les nombres d’Euler. Si G = β(2) désigne la constante de Catalan,
il en résulte en particulier que
∞
2n On(1)
2G =
X
 
2n
,
n=1 n
n

π3 ∞
2 (On(1) )2 X
n ∞
2n O(2)
= +   n ,
X
8
 
2n n 2n n
n=1 n n=1 n
∞
2 n
(On(1) )3 ∞
2n On(1) On(2) ∞
2n O(3)
12β(4) = +3 + 2   n .
X X X
   
2n n 2n n 2n n
n=1 n n=1 n n=1 n

42
 Les valeurs de la fonction β1 admettent une expression particulière.
Théorème 25 (Transformation d’Euler). Pour tout entier m ≥ 1,
∞ m
On(m)
β1 (m) = (−1)n−1 = (1 − 2−m−1 )mζ(m + 1) − β(j)β(m + 1 − j) .
X X

n=1 n j=1

Exemple 26.
∞
On(1) ∞
2n−1 1 π2
(−1)n−1 = =
X X
,
16
 
n 2n n2
n=1 n=1 n
∞ (2)
n−1 On
∞
2 n−1
On(1) 7 π
(−1) = = ζ(3) − G ,
X X
4 2
 
n 2n n 2
n=1 n=1 n
∞ (3)
n−1 On
∞
2n
(On(1) )2 X ∞
2n O(2) π4
(−1) = +   n = − G2 ,
X X
(2n) (2n) 64
 
n 2n 2 2n 2
n=1 n=1 n n=1 n

où G désigne la constante de Catalan.

Théorème 26. Pour tout entier k ≥ 1,
∞
2n 1 k
2j−1 (ln 2)k−j
2βk (1) = = (−1)j−1
X X
L̃j
(j − 1)!(k − j)!
 
n=1
2n nk+1 j=1
n

avec π
u
Z  
2
L̃k = u ln k−1
2 sin du .
0 2
Remarque 17. Avec les notations introduites par Lewin (cf. [L]), l’intégrale L̃k
(1) π
n’est autre que −Lsk+1 ( ).
2
Exemple 27. a) pour k = 2,
∞
2n 1 π2 35
2β2 (1) = = ln 2 + πG − ζ(3) ,
X
8 16
 
2n n3
n=1 n

b) pour k = 3,
∞
2n 1 π2 35
2β3 (1) = = (ln 2)2 + πG ln 2 − ζ(3) ln 2
X
16 16
 
2n n4
n=1 n
(1)
7π 5 4 ∞
H2n
+ + G(1) + 2π (−1)n
X
,
768 2 n=1 (2n + 1)2

où G et G(1) désignent respectivement les constantes de Catalan et de Ramanujan.

43
Remarque 18. Par les mêmes méthodes que précédemment, on peut montrer que
la fonction βk se prolonge analytiquement dans tout C et que ses valeurs sur les
entiers négatifs sont les nombres rationnels donnés par la formule
n
1 m
!
m
βk (−n) = (−1)j (2j + 1)n (n = 0, 1, 2, . . . )
X X

m=0 2m+1 (m + 1) k
j=0 j

En particulier, pour k = 0,
n
1 X
m
!
j m
2β0 (−n) = En = (−1) (2j + 1)n ,
X

m=0 2 j=0
m j

où les En sont les nombres d’Euler :

E0 = 1, E1 = 0, E2 = −1, E3 = 0, E4 = 5, E5 = 0, E6 = −61, etc.

On a aussi, pour k = 1,
1 1 1 1 7
β1 (0) = , β1 (−1) = , β1 (−2) = − , β1 (−3) = − , β1 (−4) = , etc.
2 4 6 2 5

44
Références

[AIK] T. Arakawa, T. Ibukiyama, et M. Kaneko, Bernoulli Numbers and Zeta

Functions, Springer Verlag, Tokyo, 2014.
[AK] T. Arakawa et M. Kaneko, Multiple zeta values, Poly-Bernoulli numbers and
related zeta functions, Nagoya Math. J., 153 (1999), 189-209.
[Be] B. C. Berndt, Ramanujan’s Notebooks Part I, Springer Verlag, New York,
1985.
[Bl] I. V. Blagouchine, Expansions of generalized Euler’s constants, J. Number
Theory 158 (2016), 365-396.
[Br] D. Broadhurst, Multiple zeta values and modular forms in quantum field
theory, in Schneider et Blümlein (eds.), Computer algebra in quantum field
theory, Springer Verlag, 33-73, 2013.
[BH] A. Bayad et Y. Hamahata, Arakawa-Kaneko L-functions and generalized
poly-Bernoulli polynomials, J. Number Theory, 131 (2011), 1020-1036.
[Ca] B. Candelpergher, Ramanujan summation of divergent series, Lecture Notes
in Math. (2016), à paraître.
[Car] P. Cartier, Fonctions polylogarithmes, nombres polyzêtas et groupes pro-
unipotents, Séminaire Bourbaki, 2000-2001, exp. n° 885, 137-173.
[Ch] K.-W. Chen, Generalized harmonic numbers and Euler sums, International
J. Number Theory (2016), à paraître.
[Cho] J. Choi, Summation formulas involving binomial coefficients, harmonic num-
bers, and generalized harmonic numbers, Abstract and Applied Analysis,
2014 (2014), Article ID 501906, 10 pages.
[Co] L. Comtet, Advanced Combinatorics, D. Reidel, Dordrecht, Holland, 1974.
[CG] J. Conway et R. Guy, The Book of Numbers, Springer Verlag, New York,
1996.
[D] K. Dilcher, Some q-series identities related to divisors functions, Discrete
Math. 145 (1995), 83-93.
[DB] A. Dil et K. Boyadzhiev, Euler sums of hyperharmonic numbers, J. Number
Theory, 147 (2015), 490-498.
[DK] A. Davydychev et M. Kalmykov, Massive Feynman diagrams and inverse
binomial sums, Nuclear Physics, B 699 (2004), 3-64.
[H] M. Hoffman, Harmonic-number summation identities, symetric functions,
and multiple zeta functions, Ramanujan J. (2016), à paraître.

45
[J] C. Jordan, Calculus of Finite Differences, Chelsea, New York, 1965.
[KZ] M. Kontsevich et D. Zagier, Periods, in Engquist (ed.) et al., Mathematics
unlimited - 2001 and beyond, Springer Verlag, 771-808, 2001.
[L] L. Lewin, Polylogarithms and Associated Functions, North-Holland, New
York, 1981.
[O1] Y. Ohno, A generalization of the duality and sum formulas on the multiple
zeta values, J. Number Theory, 74 (1999), 39-43.
[O2] Y. Ohno, Sum relations for multiple zeta values, in Aoki (ed.) et al., Zeta
functions, topology, and quantum physics, Dev. Math. 14, 131-144, Springer
Verlag, New York, 2005.
[Ra] S. Ramanujan, Notebooks, Vol. 2, Tata Institute of Fundamental Research,
Bombay, 1957.
[Ro] G.-C. Rota, Formal power series of logarithmic type, Advances in Math. 75
(1989), 1-118.
[RZ] T. Rivoal et W. Zudilin, Diophantine properties of numbers related to Ca-
talan’s constant, Math. Ann. 326 (2003), 705-721.
[Sa] Y. Sasaki, On generalized poly-Bernoulli numbers and related L-functions,
J. Number Theory, 132 (2012), 156-170.
[Si] R. Sitaramachandrarao, A formula of S. Ramanujan, J. Number Theory, 25
(1987), 1-19.
[Su] Z.-H. Sun, Invariant sequences under binomial transformation, Fibonacci
Quart. 39 (2001), 324-333.
[Y1] P. T. Young, A 2-adic formula for Bernoulli numbers of the second kind and
for the Nörlund numbers, J. Number Theory, 128 (2008), 2951-2962.
[Y2] P. T. Young, Rational series for multiple zeta and log gamma functions, J.
Number Theory, 133 (2013), 3995-4009.
[Y3] P. T. Young, Symmetries of Bernoulli polynomial series and Arakawa- Ka-
neko zeta functions, J. Number Theory, 143 (2014), 142-161.
[Y4] P. T. Young, The p-adic Arakawa-Kaneko zeta functions, J. Number Theory,
155 (2015), 13-35.
[Za] D. Zagier, The Mellin transform and other useful analytic techniques, Ap-
pendix to E. Zeidler, Quantum Field Theory I : Basics in Mathematics and
Physics, Springer Verlag, Berlin-Heidelberg-New York (2006), 305-323.
[Zh] J. Zhao, Multiple Zeta Functions, Multiple Polylogarithms and Their Special
Values, World Scientific, New Jersey, 2016.
P 2k−1 −n
[Zu] I. J. Zucker, On the series k
k and related sums, J. Number Theory,
20 (1985), 92-102.

46
Publications

[1] M-A. Coppo et P. T. Young, On shifted Mascheroni series and hyperharmonic

numbers, Journal of Number Theory 169 (2016), 1-20.
[2] M-A. Coppo et B. Candelpergher, Inverse binomial series and values of
Arakawa-Kaneko zeta functions, Journal of Number Theory 150 (2015), 98-
119.
[3] B. Candelpergher et M-A. Coppo, Le produit harmonique des suites, L’En-
seignement mathématique, 59 (2013), 39-72.
[4] B. Candelpergher et M-A. Coppo, A new class of identities involving Cauchy
numbers, harmonic numbers and zeta values, The Ramanujan Journal 27
(2012), 305-328.
[5] M-A. Coppo et B. Candelpergher, The Arakawa-Kaneko zeta function, The
Ramanujan Journal, 22 (2010), 153-162.
[6] M-A. Coppo, Nouvelles expressions des formules de Hasse et de Hermite pour
la fonction zêta d’Hurwitz, Expositiones Mathematicae 27 (2009), 79-86.
[7] M-A. Coppo, Une histoire des séries infinies, La Gazette des Mathématiciens,
120 (2009), 39-52.
[8] M-A. Coppo, Sur les sommes d’Euler divergentes, Expositiones Mathematicae,
18 (2000), 297-308.
[9] M-A. Coppo, Nouvelles expressions des constantes de Stieltjes, Expositiones
Mathematicae, 17 (1999), 349-358.
[10] B. Candelpergher, M-A. Coppo, et E. Delabaere, La sommation de Ramanu-
jan, L’Enseignement Mathématique, 43 (1997), 93-132.
[11] M-A. Coppo et C. Walter, Composante centrale du lieu de Brill-Noether du
schéma de Hilbert, Lecture Notes in Pure and Applied Mathematics, 200,
341-349, Dekker, New York, 1996
[12] M-A. Coppo, Une généralisation d’un théorème de Cayley, Comptes Rendus
de l’Académie des Sciences, Série I : Mathématiques, 314 (1992), 613-616.
[13] M-A. Coppo, Familles maximales de systèmes de points surabondants, Ma-
thematische Annalen, 291 (1991), 725-735.

47