Bernoulli umbra

From HandWiki


In Umbral calculus, the Bernoulli umbra B is an umbra, a formal symbol, defined by the relation evalBn=Bn, where eval is the index-lowering operator,[1] also known as evaluation operator [2] and Bn are Bernoulli numbers, called moments of the umbra.[3] A similar umbra, defined as evalB+n=Bn+, where B1+=1/2 is also often used and sometimes called Bernoulli umbra as well. They are related by equality B+=B+1. Along with the Euler umbra, Bernoulli umbra is one of the most important umbras.

In Levi-Civita field, Bernoulli umbras can be represented by elements with power series B=ε112ε24+3ε36401525ε5580608+ and B+=ε1+12ε24+3ε36401525ε5580608+, with lowering index operator corresponding to taking the coefficient of 1=ε0 of the power series. The numerators of the terms are given in OEIS A118050[4] and the denominators are in OEIS A118051.[5] Since the coefficients of ε1 are non-zero, the both are infinitely large numbers, B being infinitely close (but not equal, a bit smaller) to ε11/2 and B+ being infinitely close (a bit smaller) to ε1+1/2.

In Hardy fields (which are generalizations of Levi-Civita field) umbra B+ corresponds to the germ at infinity of the function ψ1(lnx) while B corresponds to the germ at infinity of ψ1(lnx)1, where ψ1(x) is inverse digamma function.

Plot of the function ψ1(ln(x)), whose germ at positive infinity corresponds to B+.

Exponentiation

Since Bernoulli polynomials is a generalization of Bernoulli numbers, exponentiation of Bernoulli umbra can be expressed via Bernoulli polynomials:

eval(B+a)n=Bn(a),

where a is a real or complex number. This can be further generalized using Hurwitz Zeta function:

eval(B+a)p=pζ(1p,a).

From the Riemann functional equation for Zeta function it follows that

evalB+p=evalB+p+12pπp+1sin(πp/2)Γ(p)(p+1)

Derivative rule

Since B1+=1/2 and B1=1/2 are the only two members of the sequences Bn+ and Bn that differ, the following rule follows for any analytic function f(x):

f(x)=eval(f(B++x)f(B+x))=evalΔf(B+x)

Elementary functions of Bernoulli umbra

As a general rule, the following formula holds for any analytic function f(x):

evalf(B+x)=DeD1f(x).

This allows to derive expressions for elementary functions of Bernoulli umbra.

evalcos(zB)=evalcos(zB+)=z2cot(z2)
evalcosh(zB)=evalcosh(zB+)=z2coth(z2)
evalezB=zez1
evalln(B+z)=ψ(z)

Particularly,

evallnB+=γ [6]
eval1πln(B+zπB+zπ)=cotz
eval1πln(B+1/2+zπB+1/2zπ)=tanz
evalcos(aB+x)=a2csc(a2)cos(a2x)
evalsin(aB+x)=a2cot(a2)sinxa2cosx

Particularly,

evalsinB=1/2,
evalsinB+=1/2,

Relations between exponential and logarithmic functions

Bernoulli umbra allows to establish relations between exponential, trigonometric and hyperbolic functions on one side and logarithms, inverse trigonometric and inverse hyperbolic functions on the other side in closed form:

eval(cosh(2xB±)1)=evalxπartanh(xπB±)=evalxπarcoth(πB±x)=xcoth(x)1
evalz2πln(B+z2πB+z2π)=evalcos(zB)=evalcos(zB+)=z2cot(z2)

References

  1. Taylor, Brian D. (1998). "Difference Equations via the Classical Umbral Calculus". Mathematical Essays in honor of Gian-Carlo Rota. pp. 397–411. doi:10.1007/978-1-4612-4108-9_21. ISBN 978-1-4612-8656-1. 
  2. Di Nardo, E. (February 14, 2022). "A new approach to Sheppard's corrections". arXiv:1004.4989 [math.ST].
  3. "The classical umbral calculus: Sheffer sequences". Lecture Notes of Seminario Interdisciplinare di Matematica 8: 101–130. 2009. http://math.fau.edu/Niederhausen/HTML/Papers/The%20classical%20umbral%20calculus%20dinardo_09.pdf. 
  4. Sloane, N. J. A., ed. "Sequence A118050". OEIS Foundation. https://oeis.org/A118050. 
  5. Sloane, N. J. A., ed. "Sequence A118051". OEIS Foundation. https://oeis.org/A118051. 
  6. Yu, Yiping (2010). "Bernoulli Operator and Riemann's Zeta Function". arXiv:1011.3352 [math.NT].