- Abstract. We present a new formula for the Bernoulli numbers as the following integral $${B}_{2m}=\frac{(-1{)}^{m-1}}{{2}^{2m+1}}{\int}_{-\infty}^{+\infty}(\frac{{d}^{m-1}}{d{x}^{m-1}}{\text{sech}}^{2}x{)}^{2}dx.$$ This formula is motivated by the results of Fairlie and Veselov, who discovered the relation of Bernoulli polynomials with soliton theory.

1 Introduction

In the paper [1] D. Fairlie and one of the authors discovered an interesting relation of the Bernoulli polynomials with the theory of the Korteweg-de Vries (KdV) equation
$${u}_{t}-6u{u}_{x}+{u}_{xxx}=0.$$
It is known since 1967 due to Gardner, Green, Kruskal and Miura [2] that this equation has infinitely many conservation laws of the form
$${I}_{m}\left[u\right]=\int {P}_{m}(u,{u}_{x},{u}_{xx},...,{u}_{m})dx,$$
where
${P}_{m}$
are some polynomials of the function
$u$
and its
$x$
-derivatives up to order
$m$
. They are uniquely defined by some homogeneity property modulo adding a total derivative and multiplication by a constant. This constant can be fixed by demanding that
${P}_{m}(u,{u}_{x},{u}_{xx},...,{u}_{m})={u}_{m}^{2}$
plus a function of derivatives of order less than
$m.$
The KdV equation is famous for its remarkable family of solutions, known as solitons, the simplest of which is a one-soliton solution
$$u=-2{\text{sech}}^{2}(x-4t),$$
corresponding to the initial profile
$u(x,0)=-2{\text{sech}}^{2}x.$
The main result of [1] is the following formula relating the Faulhaber polynomials
${F}_{m}$
with the integrals of the KdV equation:

Recall that the Faulhaber polynomials are directly related to the Bernoulli polynomials through the formula

where
${B}_{k}\left(x\right)$
and
${B}_{k}$
are Bernoulli polynomials and Bernoulli numbers respectively (see [3, 4] ). The Bernoulli numbers have the following generating function:

$$\begin{array}{c}{I}_{m-1}[-2\lambda {\text{sech}}^{2}x]=(-1{)}^{m-1}\frac{{2}^{2m+2}}{2m+1}{F}_{m}(\lambda ).\end{array}$$ | (1) |

$$\begin{array}{c}{B}_{2m+2}(x+1)=(2m+2){F}_{m}\left(\frac{{x}^{2}+x}{2}\right)+{B}_{2m+2},\end{array}$$ | (2) |

$$\frac{z}{{e}^{z}-1}={\sum}_{k=0}^{\infty}\frac{{B}_{k}}{k!}{z}^{k}.$$
All odd Bernoulli numbers except
${B}_{1}=-\frac{1}{2}$
are zero and the first even Bernoulli numbers are
$${B}_{0}=1,{B}_{2}=\frac{1}{6},{B}_{4}=-\frac{1}{30},{B}_{6}=\frac{1}{42},{B}_{8}=-\frac{1}{30},{B}_{10}=\frac{5}{66},{B}_{12}=-\frac{691}{2730},...$$
They play an important role in analysis, number theory, algebraic topology and many other areas of mathematics.

In this note we show that the relationship with soliton theory brings the following formula for the Bernoulli numbers:

where
$({\text{sech}}^{2}x{)}^{(m-1)}$
denotes the
$(m-1)$
th derivative of
${\text{sech}}^{2}x.$
We present also a direct proof of this formula based on the elegant arguments due to Logan, which we found in the book [5] by Graham, Knuth and Patashnik.

$$\begin{array}{c}{B}_{2m}=\frac{(-1{)}^{m-1}}{{2}^{2m+1}}{\int}_{-\infty}^{+\infty}\left(\right({\text{sech}}^{2}x{)}^{(m-1)}{)}^{2}dx,m\ge 1,\end{array}$$ | (3) |

2 Proofs of the main formula.

Let us first present the proof based on the results of [1] . Recall that the Faulhaber polynomial
${F}_{m}\left(\lambda \right)$
with
$m\ge 1$
has the form
$${F}_{m}\left(\lambda \right)={\alpha}_{2}^{m}{\lambda}^{2}+{\alpha}_{3}^{m}{\lambda}^{3}+...+{\alpha}_{m+1}^{m}{\lambda}^{m+1}$$
with some rational coefficients
${\alpha}_{2}^{m},{\alpha}_{3}^{m},...,{\alpha}_{m+1}^{m}$
(see [3] ). Since the only quadratic term in the density of the KdV integral
${I}_{m-1}$
is
$\int {u}_{m-1}^{2}dx,$
we have from the relation ( 1 ) that

Differentiating the formula ( 2 ) twice with respect to
$x$
gives
$$\left(\frac{2x+1}{2}{)}^{2}{F}_{m}^{\prime \prime}\right(\frac{x(x+1)}{2})+{F}_{m}^{\prime}(\frac{x(x+1)}{2})=\frac{1}{2m+2}{B}_{2m+2}^{\prime \prime}(1+x)=(2m+1\left){B}_{2m}\right(1+x)$$
because
${B}_{2m+2}(x{)}^{\prime \prime}=(2m+2\left)\right(2m+1\left){B}_{2m}\right(x).$
Since
${F}_{m}^{\prime}\left(0\right)=0$
this reduces to
$$\frac{1}{4}{F}_{m}^{\prime \prime}\left(0\right)=(2m+1){B}_{2m}\left(1\right).$$
Now using the well-known symmetry
${B}_{k}(1-x)=(-1{)}^{k}{B}_{k}(x),$
we have
${B}_{2m}\left(1\right)={B}_{2m}\left(0\right)={B}_{2m}$
and thus
$${\alpha}_{2}^{m}=\frac{1}{2}{F}_{m}^{\prime \prime}\left(0\right)=2(2m+1){B}_{2m}.$$
Substituting this into ( 4 ) we come to the formula ( 3 ).

$$\begin{array}{c}{\int}_{-\infty}^{+\infty}\left(\right({\text{sech}}^{2}x{)}^{(m-1)}{)}^{2}dx=(-1{)}^{m-1}\frac{{2}^{2m}}{2m+1}{\alpha}_{2}^{m}.\end{array}$$ | (4) |

We are now going to prove the formula ( 3 ) directly without reference to soliton theory. We borrow the main idea from the book [5] , where it is attributed to Logan. Consider the integral
${J}_{m}=\frac{(-1{)}^{m-1}}{{2}^{2m+1}}{\int}_{-\infty}^{+\infty}\left(\right({\text{sech}}^{2}x{)}^{(m-1)}{)}^{2}dx.$
Integrating
${J}_{m}$
by parts
$m-1$
times gives
$${J}_{m}=\frac{1}{{2}^{2m+1}}{\int}_{-\infty}^{+\infty}({\text{sech}}^{2}x{)}^{(2m-2)}{\text{sech}}^{2}xdx=\frac{1}{{2}^{2m+1}}{\int}_{-\infty}^{+\infty}tanh{x}^{(2m-1)}tanh{x}^{\left(1\right)}dx.$$
Let
$y=tanhx$
then
$${J}_{m}=\frac{1}{{2}^{2m+1}}{\int}_{-1}^{+1}{T}_{2m-1}\left(y\right)dy,$$
where the polynomial
${T}_{k}\left(y\right)$
is the
$k$
-th derivative of
$y=tanhx$
rewritten in terms of
$y$
:
$${T}_{1}={y}^{\left(1\right)}=1-{y}^{2}$$
$${T}_{2}={y}^{\left(2\right)}=-2y{y}^{\left(1\right)}=-2y(1-{y}^{2})=-2y+2{y}^{3}$$
$${T}_{3}={y}^{\left(3\right)}=-2{y}^{\left(1\right)}+6{y}^{2}{y}^{\left(1\right)}=(1-{y}^{2})(-2+6{y}^{2})=-2+8{y}^{2}-6{y}^{4}$$
$${T}_{4}={y}^{\left(4\right)}=16y{y}^{\left(1\right)}-24{y}^{3}{y}^{\left(1\right)}=(1-{y}^{2})(16y-24{y}^{3})=16y-40{y}^{3}+24{y}^{5}$$
and so on. These polynomials can be determined by the recurrence formula

with
${T}_{0}\left(x\right)=x.$
^{1 }
They have integer coefficients with the highest one being equal to
$(-1{)}^{n}n!$
and have the symmetry
${T}_{n}(-x)=(-1{)}^{n-1}{T}_{n}(x).$
Note that the non-zero coefficients have the alternating signs and their total sum is zero. The integers
$\left|{T}_{2m-1}\right(0\left)\right|$
are called tangent numbers (see e.g. [5] ).

$$\begin{array}{c}{T}_{n}\left(x\right)=(1-{x}^{2}){T}_{n-1}(x{)}^{\prime}\end{array}$$ | (5) |

Let us consider now the generating function
$$T(x,z)={\sum}_{n\ge 0}{T}_{n}\left(x\right)\frac{{z}^{n}}{n!}.$$
Lemma 1. The generating function for the polynomials
${T}_{n}\left(x\right)$
is

Indeed it is easy to see that when
$x=tanhw$
the function
$T(x,z)$
becomes
$tanh(z+w)$
. Now the claim follows from the Taylor formula and the definition of the polynomials
${T}_{n}\left(x\right).$
From this one can derive an interesting relation between Bernoulli and tangent numbers, see [5] , formula (6.93) and discussion after that. We will need however a slightly different result.

$$\begin{array}{c}T(x,z)=\frac{sinhz+xcoshz}{coshz+xsinhz}.\end{array}$$ | (6) |

Lemma 2. Bernoulli numbers can be written as

for all
$m>1.$
Indeed, this formula is obvious for odd
$m$
bigger than 1 since
${T}_{m-1}\left(x\right)$
is an odd function and both sides of ( 7 ) are then equal to
$0.$
In order to prove the formula for even
$m$
, let us consider first the left hand side of ( 6 ). The function
$T(x,z)$
can be rewritten as
$$T(x,z)=\frac{coshz(x+\frac{sinhz}{coshz})}{sinhz(x+\frac{coshz}{sinhz})}=cothz-\frac{1}{{sinh}^{2}z(x+cothz)}.$$
Integrating with respect to
$x$
gives
$$\int T(x,z)dx=xcothz-\frac{1}{{sinh}^{2}z}ln|x+cothz|,$$
which leads to

The expansion of the function
$cothz$
can be written in terms of the Bernoulli numbers:

This easily follows from the identity
$$\frac{z}{2}coth\frac{z}{2}=\frac{z}{2}+\frac{z}{{e}^{z}-1}.$$
From ( 8 ) and ( 9 ) it follows that

On the other hand from ( 6 ) we have

since
${T}_{2n}\left(x\right)$
are odd polynomials. Comparing ( 10 ) and ( 11 ) and replacing
${a}_{2n-1}$
by
$\frac{{2}^{2n}{B}_{2n}}{\left(2n\right)!}$
we obtain
$${B}_{2n}=\frac{1}{{2}^{2n+1}}{\int}_{-1}^{1}{T}_{2n-1}\left(x\right)dx.$$
This proves Lemma 2.

$$\begin{array}{c}{B}_{m}=\frac{1}{{2}^{m+1}}{\int}_{-1}^{1}{T}_{m-1}\left(x\right)dx\end{array}$$ | (7) |

$$\begin{array}{c}{\int}_{-1}^{1}T(x,z)dx=2cothz-\frac{2z}{{sinh}^{2}z}=2cothz+2z{coth}^{\prime}z.\end{array}$$ | (8) |

$$\begin{array}{c}cothz=\frac{1}{z}+{\sum}_{n=1}{a}_{2n-1}{z}^{2n-1},{a}_{2n-1}=\frac{{2}^{2n}{B}_{2n}}{\left(2n\right)!}.\end{array}$$ | (9) |

$$\begin{array}{c}{\int}_{-1}^{1}T(x,z)dx=4{\sum}_{n=1}n{a}_{2n-1}{z}^{2n-1}.\end{array}$$ | (10) |

$$\begin{array}{c}{\int}_{-1}^{1}T(x,z)dx={\sum}_{n\ge 1}\left({\int}_{-1}^{1}{T}_{2n-1}\right(x\left)dx\right)\frac{{z}^{2n-1}}{(2n-1)!},\end{array}$$ | (11) |

Remark. As a corollary we have the following representation of the Bernoulli number
${B}_{2n}$
as a sum of the fractions of the type
$${B}_{2n}=\frac{1}{{2}^{2n}}\left(\frac{{k}_{0}}{1}+\frac{{k}_{1}}{3}+\frac{{k}_{2}}{5}+\cdot \cdot \cdot +\frac{{k}_{n}}{2n+1}\right),$$
where the denominators are consecutive odd numbers and the numerators are the coefficients of the polynomials
${T}_{2n-1}\left(x\right).$
Comparison with the formula (6.93) from the book [5] leads to the following relation with the tangent numbers:

Now combining Lemma 2 with the definition of the polynomials
${T}_{k}$
we have our main Theorem. Bernoulli numbers
${B}_{2m}$
with
$m\ge 1$
have the following integral representation
$${B}_{2m}=\frac{(-1{)}^{m-1}}{{2}^{2m+1}}{\int}_{-\infty}^{+\infty}(\frac{{d}^{m-1}}{d{x}^{m-1}}{\text{sech}}^{2}x{)}^{2}dx.$$

$$\begin{array}{c}{\int}_{-1}^{1}{T}_{n}\left(x\right)dx=\frac{2n+2}{{2}^{n+1}-1}{T}_{n}\left(0\right).\end{array}$$ | (12) |

^{1 }
These polynomials should not be confused with the classical Chebyshev (Tchebycheff ) polynomials, which are also denoted as
${T}_{n}\left(x\right)$
. We follow the notations of the book [5] , where trigonometric version of these polynomials was considered.

3 Concluding remark.

There exists a generalisation of the Bernoulli numbers, known as Bernoulli-Hurwitz numbers
$B{H}_{k}$
and related to the coefficients of the Laurent series of the Weierstrass elliptic function
$\mathcal{\wp}\left(z\right)$
at zero:

(see [6, 7] ).

$$\begin{array}{c}\mathcal{\wp}\left(z\right)={z}^{-2}+{\sum}_{k=1}^{\infty}\frac{B{H}_{2k+2}}{(2k+2)\left(2k\right)!}{z}^{2k}\end{array}$$ | (13) |

Our formula suggests a different generalisation of Bernoulli numbers related to an elliptic curve, namely

where the integral is taking over one of the cycles on the corresponding elliptic curve. Obviously these numbers depend also on the choice of such a cycle, which is an element of the first homology group of the curve. They naturally appear in relation with the spectral density for the classical Lamé operator [8] .

$$\begin{array}{c}{B}_{2m}^{ell}=\frac{(-1{)}^{m-1}}{{2}^{2m+1}}\int (\mathcal{\wp}(z{)}^{(m-1)}{)}^{2}dz,m>1,\end{array}$$ | (14) |

4 Acknowledgements.

We are grateful to Victor Enolski and Nick Trefethen for very useful and stimulating discussions.

References
Department of Mathematical Sciences, Loughborough University, Loughborough, Leicestershire, LE11 3TU, UK E-mail address : M.Grosset@lboro.ac.uk Department of Mathematical Sciences, Loughborough University, Loughborough, Leicestershire, LE11 3TU, UK and Landau Institute for Theoretical Physics, Moscow, Russia E-mail address : A.P.Veselov@lboro.ac.uk

- D. B. Fairlie and A. P. Veselov, Faulhaber and Bernoulli polynomials and solitons, Physica D 152-153 (2001) 47-50.
- R.M.Miura, C.S.Gardner and M.D.Kruskal, Korteweg-de Vries equation and generalisations.II. J. Math.Phys., v.9, n.8 (1968),1204-1209.
- D.E.Knuth, Johann Faulhaber and the sums of powers. Math. of Computation, v.61, n.203 (1993), 277-294.
- M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables, 1968.
- R. L. Graham, D. E. Knuth, O. Patashnik, Concrete mathematics: a foundation for computer science, Addison-Wesley Publishing Company Inc, 1994, p. 283-290.
- N.M. Katz, The congruences of Clausen-von Staudt and Kummer for Bernoulli-Hurwitz numbers. Math. Ann. 216 (1975), 1-4.
- Y. Onishi, Theory of the generalized Bernoulli-Hurwitz numbers for the algebraic functions of cyclotomic type and the universal Bernoulli numbers. math.NT/0406096.
- M-P. Grosset and A.P. Veselov, Lamé equation, quantum top and elliptic Bernoulli polynomials. In preparation.