76

Earlier today, I was talking with my friend about some "cool" infinite series and the value they converge to like the Basel problem, Madhava-Leibniz formula for $\pi/4, \log 2$ and similar alternating series etc.

One series that popped into our discussion was $\sum\limits_{n=1}^{\infty} \frac{1}{n^n}$.

Proving the convergence of this series is trivial but finding the value to which converges has defied me so far. Mathematica says this series converges to $\approx 1.29129$.

I tried Googling about this series and found very little information about this series (which is actually surprising since the series looks cool enough to arise in some context).

We were joking that it should have something to do with $\pi,e,\phi,\gamma$ or at the least it must be a transcendental number :-).

My questions are:

  1. What does this series converge to?
  2. Does this series arise in any context and are there interesting trivia to be known about this series?

I am actually slightly puzzled that I have not been able to find much about this series on the Internet. (At least my Google search did not yield any interesting results).

  • 5
  • 2
    I asked a related question on MO: http://mathoverflow.net/questions/33397/is-there-a-standard-name-for-mz1-sum-n1-infty-fracznnn . It is possible to make "deformed" trigonometric functions by replacing $e^{z}$ with $M(z)$, but so far I haven't found references to them in lists of special functions. – graveolensa Feb 10 '11 at 09:32
  • @deoxygerbe: Very Interesting! I was thinking of a function on similar lines and study its behavior. The function actually looks a cool function. It converges for all $z$ and is holomorphic. An interesting function to play around with! –  Feb 10 '11 at 09:40

4 Answers4

57

Certainly you need to check out Sophomore's Dream.

27

This series turns out to be equal to the integral of $1/x^x$ from $0$ to $1$. Sorry I am bad with the notation needed for math symbols here; but if you rewrite $x^{-x}$ as the exponential of $x\log x$, you can expand it as a taylor series and integrate term by term; to find it equals your sum. This doesn't give an explicit value; but I think it's a pretty cool identity.

EDIT: Just putting what you said into symbols, probably for the fun of it.

$$\int_0^1 x^{-x} dx=\int_0^1 e^{-x \log x} dx$$

$$ e^{-x \log x} =\sum_{k=0}^\infty (-1)^k \frac{x^k\log^kx}{k!}$$

But since

$$\int_0^1 x^k \log^k x dx =(-1)^k \frac{k!}{(k+1)^{k+1}}$$

we get

$$\int_0^1 x^{-x} dx= \sum_{k=0}^\infty \frac{1}{(k+1)^{k+1}}=\sum_{k=1}^\infty \frac{1}{k^k}$$

Pedro
  • 122,002
Steve
  • 919
1

Liouville's theorem implies that if your constant is irrational then it is in fact transcendental. Unfortunately, the trivial way of proving irrationality doesn't work. On second thought, the same problem also prevents us from applying Liouville's theorem!

Yuval Filmus
  • 57,157
  • I am unable to immediately see that if the constant is irrational then it has to be a Liouville number. Could you throw light on the proof? –  Feb 11 '11 at 02:57
  • @Sivaram: You have a point... – Yuval Filmus Feb 11 '11 at 04:45
  • 4
    To apply Liouville's theorem (or the much stronger Roth version) you need an infinite sequence of rational approximations. The smallest upper bound I get on the denominators of the obvious sequence of rational approximations ($\sum_{n=1}^l \frac{1}{n^n}$) is $q = \prod_{i=1}^n i^i$, whereas the best bound I get for $\left| \frac{p}{q} - \alpha \right|$ is $\frac{1}{n^n} \sum_{k=1}^\infty \frac{1}{k^k}$. That's huge compared to $1/q^r$ so the approximation theory doesn't apply at all - unless you produce a much better sequence of rationals and bound them very tightly. – quanta May 07 '11 at 21:38
1

Since the question asks for interesting trivia, I believe the following is appropriate as an answer.

Digits

The digit sequence is A073009 on OEIS.

This simple Python Script may be used to compute the digits of the sequence.

from decimal import Decimal, getcontext

getcontext().prec = 1003

total = Decimal(0) for n in range(1, 10000): term = Decimal(1) / Decimal(n) ** n total += term

print(total)

With this I've obtained the following:

1.291285997062663540407282590595600541498619368274522317310002445136944538765234455558817041129429708984995070924815430548410487419284864197579163555947913696496974156878020799729177948273009025649230550720966638128467012053685745978703001277894129288253551770222383375319345749259967779648300849549111066964975501051975742911621097021561669532897689242789005809390814788094036799305589535200633716110465094638606808864998606531021853412479159737305271068682465224677033686046987023420196583143133968738817295689355368517985214206662641654380612245699409663560438852399693813044840101532338556989547899226146597068180753342912289091004995136410358472374167966099403742887228090823947240301242337506966587431476835029834700965969301980712205941547423918884954889204314784037389693592832744937301860181757952468190913559650620576842700890732654713723383484718562324804417342338565270511374482208606983811697064478963155480311086868468078070105703423000095477662829927022264266182213029160934485049255679992

Continued Fraction

Using the first few of the digits found, Wolfram Alpha computes the first terms of its continued fraction's sequence.

[1; 3, 2, 3, 4, 3, 1, 2, 1, 1, 6, 7, 2, 5, 3, 1, 2, 1, 8, 1, 2, 4, 1, 9, 3, 1, 1, 18, 1, 1, 29, 4, 1, 5, 2, 167, 1, 62, 4, 2, 1, 3, 3, 27, 1, 9, 1, 46, 1, 3, 2, 2, 1, 1, 3, 2, 10, 73, 1, 11, 1, 2, 1, 1, 18, 1, 4, 1, 4, 6, 1, 4, 4, 1, 6, 1, 1, 1, 2, 1, 7, 8, 4, 1, 3, 1, 4, 28, 2, 1, 6, 2, 10, 3, 1, 2, 2, 1, 4, 1, 1, 10, 1, 2, 1, 26, 5, 2, 3, 1, 4, 2, 1, 1, 1, 3, 1, 4, 7, 45, 2, 2, 6, 7, 1, 1, 5, 1]

Approximations

$$\begin{aligned} \ln(3)^e=1.2912\color{red}{98757}\dots~&(E\sim 10^{-5})\\ \sqrt[\pi+\frac14]{\frac{\phi+\pi}{2}}=1.29128599\color{red}{20365}\dots~&(E\sim10^{-9})\\ \frac38\pi\cot^2\left(\frac{1830385}{2400594}\right)=1.29128599706266354\color{red}{88534}\dots~&(E\sim10^{-18}) \end{aligned}$$

Alma Arjuna
  • 3,759