Galvanized by comments in the recent post Bernoulli Numbers and the Harmonic Oscillator by John Baez on the blog site The n-Category Café, I’ve decided to post a fraction of some sets of notes I’ve developed on the appearance of the Bernoulli numbers, polynomials, and function in diverse areas of mathematics. I find the umbral Sheffer calculus to often be the most comprehensive and elegant way to derive and view the various relationships among the formulas involved in the related concepts, so I pay homage to Blissard in the title who was one of the first to use umbral characters in deriving Bernoulli number identities.
First, I’ll reiterate the relationship between a formal group law and the Bernoulli numbers:
The refs attached to my formulas dated Sep 18 2014 in OEIS A008292 are relevant to the connections between formal group laws and the Bernoulli numbers. As I commented in the Café, the shifted reciprocal of the bivariate e.g.f. for the Eulerian polynomials, the
for the hyperbolic formal group law, generates the Bernoulli numbers. In more detail, given the compositional-inverse pair of functions


and


the Bernoulli numbers appear in the shifted reciprocal


Relevant refs are “Elliptic formal group laws, integral Hirzebruch genera and Krichever genera” by Victor Buchstaber and Elena Bunkova (the Bernoulli numbers occur on pgs. 36 & 37), “Towards generalized cohomology Schubert calculus via formal root polynomials” by Cristian Lenart and Kirill Zainoulline, and the more recent papers “Theta divisors and permutohedra” and “Chern-Dold character in complex cobordisms and theta divisors” by Buchstaber and A.P. Veselov (the Eulerian polynomials are the h-polynomials of the permutohedra).
satisfies

a Ricatti equation, and can be written in terms of a Weierstrass elliptic function (see Buchstaber & Bunkova), so

and the formal group law is

called the hyperbolic formal group law and related to a generalized cohomology theory by Lenart and Zainoulline.
Note
is a linear combination of the infinitesimal generators for
.
In my pdf in my post “The Elliptic Lie Triad: Riccati and KdV Equations, Infinigens, and Elliptic Genera” at my Web blog Shadows of Simplicity, I note in detail the relations of the compositional inverse pair
and
to the KdV equation and, therefore, to the Burgers’ equation, the heat equation, the Schwarzian derivative, and more.
The coefficients of
containing the Bernoulli numbers can be used to generate
via the partition polynomials of OEIS A134264, the refined Narayana partition polynomials, enumerating / labeling Dyck paths, polygon dissections, forests, and noncrossing partitions; therefore, the Bernoulli numbers can be generated from
by the inverse set of partition polynomials of A350499. This pair of sets of partition polynomials allows conversions between the free cumulants and the free moments of free probability theory and between the coefficients of a compositional-inverse pair of Laurent series. (In fact, at the bottom of p. 11 of “Todd polynomials and Hirzebruch numbers” by Buchstaber and Veselov is an example of the use of A350499 to calculate the Bernoulli numbers, of which they were unaware until I notified them in June.)
For details of the relationship between the Hirzebruch-Todd criterion and the K-P equation mentioned in comments in the Café, see my pdf
The Khovanskii-Pukhliko formula, the Hirzebruch-Todd criterion, the Bernoulli op, and umbral witchcraft.
(See also the answers to the MO-Q “Hirzebruch’s motivation of the Todd class” and “Formal Groups, Witt vectors and Free Probability” by R. Friedrich and J. McKay.)
For the links between the Bernoulli polynomials, Newton series, and the derivative / tangent space of functions, see my pdf
Conjuring the continuous from the discrete with umbral witchcraft: the derivative and the Bernoulli function.
For a review of a natural extension of the Bernoulli polynomials to the Bernoulli function of Milnor, see my pdf
Résumé of the Bernoulli function.
(See also my answer to the MO-Q “Transformation converting power series to Bernoulli polynomial series” and my posts “The Kervaire-Milnor Formula and Bernoulli Numbers“, “Differintegral Ops and the Bernoulli and Reciprocal Polynomials“, and “Appell polynomials, cumulants, noncrossing partitions, Dyck lattice paths, and inversion“.
For some relationships among the quantum zeta function, the quantum partition function, heat kernels, statistical mechanics, the harmonic oscillator, Hermite and Bernoulli polynomials, the Bernoulli function and Hurwitz zeta function, and Bose-Einstein condensates, see my pdf
(forthcoming)