Con la reciente edición es posible algo como lo siguiente, pero puede que no sea tan bonito como esperabas; lo pongo por si alguien puede detectar cómo mejorarlo.
Supongo que α∈Z ya que aparece como límite superior en la suma, pero si no es así lo siguiente no se cumple. También restrinjo 0<u para poder utilizar una serie geométrica.
A simple vista esto parece el producto de dos términos que se asemejan a la función generadora, así que empezamos diciendo que su expresión es el producto de: ueu(αx+y)−1+eu=∞∑k=0Bk(αx+y)ukk! y: 1αuα−1∑k=11eu−e2πik/α=u1αα−1∑k=1∞∑n=1e−une2iπk(n−1)α=u∞∑n=1e−un1αα−1∑k=1e2iπk(n−1)α
donde en la primera línea simplemente usamos la fórmula y en la segunda usamos una serie geométrica y luego invertimos el orden de la suma ya que la suma interna es absolutamente convergente para u>0 . Entonces, ten en cuenta que: 1αα−1∑k=1e2iπk(n−1)α=δαl,n−1−1α=δαl+1,n−1α
donde he introducido l para significar que el delta de Kronecker será 1 sólo cuando n−1 es un múltiplo de α y 0 en caso contrario; nótese también que como n comienza en 1 , l comienza en 0 . Teniendo en cuenta esto tenemos: u∞∑n=1e−un1αα−1∑k=1e2iπk(n−1)α=u∞∑l=0e−u(αl+1)−u1α∞∑n=1e−un =−e−uu−1+e−uα−uα(−1+eu)=1α∞∑j=0(Bj(1α)(−α)j−Bj(0))ujj!
donde de nuevo sumamos la serie geométrica y luego utilizamos la fórmula; nota Bj(0) es un número de Bernoulli. Entonces tienes eso: 1αu2eu(αx+y)eu−1α−1∑k=11eu−e2πik/α=(∞∑k=0Bk(αx+y)ukk!)(1α∞∑j=0(Bj(1α)(−α)j−Bj(0))ujj!) y después de tomar el producto de Cauchy de las dos series: 1αu2eu(αx+y)eu−1α−1∑k=11eu−e2πik/α=∞∑n=0cnunn! donde cn es la convolución: cn=1αn∑m=0n!(Bn−m(1α)(−α)n−m−Bn−m(0))Bm(αx+y)(n−m)!m!
cn=1αn∑m=0(nm)(Bn−m(1α)(−α)n−m−Bn−m(0))Bm(αx+y)
Cuando la serie converge, esta expresión se comprueba numéricamente para los valores que he probado . Confieso que no he encontrado el límite superior de u todavía pero conjeturo que la serie se mantiene para 0<u<1 aunque posiblemente tenga una dependencia de α que aún no he encontrado.