Posts: 1,631
Threads: 107
Joined: Aug 2007
11/19/2008, 02:20 PM
(This post was last modified: 11/19/2008, 02:22 PM by bo198214.)
In order to understand it myself and as a service for the non-German community, in the following posts I will provide the construction of Kneser (for an analytic super logarithm) which he gives in his paper:
Hellmuth Kneser, Reelle analytische Lösungen der Gleichung \( \phi(\phi(x))=e^x \) und verwandter Funktionalgleichungen. J.Reine Angew. Math. 187, (1949), 56-67
I will proceed in several posts because my time is currently quite limited, and so I will give account of succeeding portions of his construction.
Posts: 1,631
Threads: 107
Joined: Aug 2007
11/19/2008, 03:25 PM
(This post was last modified: 11/19/2008, 03:30 PM by bo198214.)
As mentioned in the title of his paper Kneser originally seeks for an analytic solution \( \phi \) of the functional equation \( \phi(\phi(x))=e^x \) (This functional equation was considerably discussed at the conference of the German mathematician's association, October 1941 in Jena. So there the idea was born to construct an analytic solution.)
As is very well-known the construction of a half-iterate can be reduced to the solution of the Abel equation, which Kneser gives here in a bit more generality:
\( \psi(f(x))=\psi(x)+\beta \).
If one has such a solution \( \psi \) then it is easy to construct the half iterate by:
\( \phi(x)=\psi^{-1}(\frac{\beta}{2}+\psi(x)) \)
because then
\( \begin{align*}
\phi(\phi(x))
&=\psi^{-1}(\frac{\beta}{2}+\psi(\psi^{-1}(\frac{\beta}{2}+\psi(x)))\\
&=\psi^{-1}(\frac{\beta}{2}+\frac{\beta}{2} +\psi(x))\\
&=\psi^{-1}(\beta+\psi(x))\\
&= \psi^{-1}(\psi(f(x))\\
&=f(x)
\end{align*} \)
Closely related to the Abel equation is the Schroeder equation:
\( \chi(f(x))=\gamma \chi(x) \)
If one has a solution of the Schroeder equation, then \( \psi(x)=\log(\chi(x)) \) is a solution of the Abel equation with \( \beta=\log(\gamma) \), because:
\( \psi(f(x))=\log(\chi(f(x))=\log(\gamma)+\log(\chi(x))=\log(\gamma)+\psi(x) \)
And one can also derive the fractional iterates directly from the Schroeder function by:
\( f^{t}(x)=\chi^{-1}(\gamma^t\chi(x)) \).
The good thing is now that for holomorphic functions with fixed point \( c \) and \( a:=f'( c)\neq 0,1 \) there is always the so called principal Schroeder function (holomorphic in a vicinity of \( c \)) which satisfies the Schroeder equation for \( \gamma=a \) and which satisfies \( \chi'( c)=1 \).
There are two possibilities to construct this principal Schroeder function. First by development of the powerseries in \( c \) and chosing the coefficients of \( \chi \) such that they satisfy the Schroeder equation.
And second by the limit
\( \chi(x)=\lim_{n\to\infty} \frac{f^n(x)-c}{a^n} \)
As you can easily verify if \( \chi \) is a solution of the Schroeder equation then \( b\chi(z) \) is also a solution for any constant \( b \). They are called "regular" (by Szekeres) if \( \chi \) is the principal Schroeder function.
Those regular solutions are characterized by that \( \chi^{-1}(\gamma^t \chi(x)) \) is holomorphic in \( c \).
As a first step Kneser computes this principal Schroeder function \( \chi \) of \( \exp(z) \) at \( \exp \)'s first fixed point \( c\approx0.318+1.337i \) in the upper half plane, the fixed point nearest to the real axis. \( \chi(z) \) is however not real on the real axis, so he determines some mapping properties of \( \chi \) and later manipulates \( \chi \) to become real and analytic on the real axis.
Posts: 1,631
Threads: 107
Joined: Aug 2007
11/22/2008, 06:11 PM
(This post was last modified: 11/23/2008, 09:55 AM by bo198214.)
So up to now we have a Schroeder function \( \chi \) in a vicinity \( H_0 \) of the fixed point \( c \). This function satisfies:
\( \chi(\log(z))=\chi(z)/c \)
\( \chi'( c)=1 \)
The next thing Kneser does is to analytically continue \( \chi \) from this small vicinity \( H_0 \) to the whole upper halfplane without the points \( \exp^n(0) \), that is:
\( H=\{z\in\mathbb{C}: \Im(z)\ge 0, z\neq \exp^n(0), n= 0,1,\dots\} \).
He first verifies that \( \lim_{n\to\infty} \log^n(z) = c \) for each \( z\in H \) (where the cut of the logarithm is the usual one, i.e. \( -\pi<\Im(\log(z))\le \pi \).)
And then he continues the function along the increasing sets \( H_n \), which contain all the points \( z \) such that \( \log^n(z) \) is contained in that initial vicinity \( H_0 \) of \( c \).
The inverse function (in a vicinity of \( c \)) \( \chi^{-1} \) satisifies:
\( \chi^{-1}(cz)=\exp(\chi^{-1}(z)) \) and hence can be continued to the whole complex plane, i.e. is an entire function.
To see the properties of \( \chi \) on \( H \) he considers the following lines and areas. Each increasing index number indicates the application of exponentiation, for example \( A_1=\exp(A_0) \), \( H_{-1}=\log(H_0) \). The lines are without end points, the areas are without boundaries.
These areas are mapped by \( \chi \) to (the letters indicate the source area):
Now \( \chi(H_{-1})\cup \chi(H_0)\cup \chi(H_1) \) is simply connected and does not contain 0 (only in the boundary). Hence it is possible to define a holomorphic logarithm \( \underline{\log} \) on that domain, he defines
\( \psi(z)=\underline{\log}(\chi(z)) \) on \( H_{-1}\cup H_0\cup H_1 \). Which then satisfies
\( \psi(e^z)=\psi(z)+\underline{\log}( c)=\psi(z)+c \).
And this is the image under \( \psi \), considering that \( \underline{\log} \) is biholomorphic:
Define \( L_0=\psi(H_0) \) and define \( L_n=L_0+nc \), then we have \( L_{-1}=\psi(H_{-1}) \) and \( L_1=\psi(H_1) \). Define \( L=\bigcup_{n=-\infty}^\infty L_n \). So one can see that the boundary of \( L \) consists only of the cyan and violet arcs, hence the points of the real axis in \( H \) are mapped to the boundary of \( L \). This property is used in the final step of Kneser's construction, which follows in my next post.
Posts: 1,631
Threads: 107
Joined: Aug 2007
11/23/2008, 01:00 PM
(This post was last modified: 11/23/2008, 01:03 PM by bo198214.)
As the last step we use the Riemann mapping theorem, to map the area \( L \) biholomorphically to the upper halfplane, which maps the boundary to the real axis. By some reason (to be explained) the corresponding function \( \Psi \) indeed satisfies the Abel equation.
But we do that in two steps first we map \( L \) to the unit disk \( \mathbb{E} \) via \( \alpha: L\leftrightarrow\mathbb{E} \) and then we map \( \mathbb{E} \) to the upper (open) halfplane \( \mathbb{H} \) via \( \beta:\mathbb{E}\leftrightarrow\mathbb{H} \).
The existence of \( \alpha \) is guarantied by the Riemann mapping theorem. The translation \( z\mapsto z+c \) maps \( L \) to \( L \) without having a fixed point, hence the corresponding mapping in \( \mathbb{E} \), \( \tau:\alpha(z)\mapsto \alpha(z+c) \) maps the unit disk into the unit disk without a fixed point. It is well known that such a function needs to be linear fractional, i.e. of the form \( \tau(z)=\frac{d_1 z+d_2}{d_3 z +d_4} \). Kneser shows that must be parabolic, that means it has only one fixed point \( p \) at the boundary of \( \mathbb{E} \).
To map the unit disk to the upper half plane we can again use a linear fractional transformation \( \beta(z) \). There are enough parameters to chose it such that \( p \) is mapped to infinity and such that \( \beta(\tau(z))=z+1 \).
Now we can define
\( \Psi(z)=\beta(\alpha(\psi(z)) \) which has the property
\( \Psi(e^z)=\beta(\alpha(\psi(z)+c))=\beta(\tau(\alpha(\psi(z)))=1+\Psi(z) \).
and is biholomorphic on the interior of \( H_{-1}\cup H_0\cup H_1 \).
The interior of \( H_{-1}\cup H_0\cup H_1 \) is mapped to some area bordering on the real line. The boundary on the real axis is hence mapped to the real axis. More precisely the interval \( (\log(\Re( c)),e^{|c|}) \) is mapped to the real axis. By the Schwarz reflection principle it can be continued to the complex conjugate of \( H_{-1}\cup H_0 \cup H_1 \) especially it is analytic on \( (\log(\Re( c)),e^{|c|}) \) and can from there be continued to a vicinity of the whole real axis by \( \Psi(e^x)=\Psi(x)+1 \) and \( \Psi(\ln(x))=\Psi(x)-1 \).
Posts: 1,631
Threads: 107
Joined: Aug 2007
(06/02/2009, 10:04 AM)Kouznetsov Wrote: (06/02/2009, 09:24 AM)bo198214 Wrote: (06/02/2009, 03:10 AM)Kouznetsov Wrote: About Kneser's expansion: it would be good to check, that we evaluate the same function expanding it at fixed point L and doing that at L^*.... Yes, there are singularities at the real axis at \( \exp^{[n]}(0) \). But thats only an intermediate step of Kneser's construction. It is sufficient to show that the two expansions coincide along any segment of positive length. Then the two functions coincide in the whole range of the holomorphizm.
They dont coincide, because they are not real at the real axis.
Kneser's trick is to prepend a conformal mapping, that makes it real at the real axis.
Please have a look into my previous posts in this thread.
Moderator's note: I moved the two previous posts from thread overview paper co-author inviation to here.
Posts: 684
Threads: 24
Joined: Oct 2008
01/23/2010, 01:01 PM
(This post was last modified: 01/27/2010, 08:36 PM by sheldonison.)
(corrections and additions, corrected notation)
It occured to me that complex fourier analysis is at the heart of a lot of key topics for tetration. I finally have gotten to the point where I understand Kouznetsov's Cauchy integral, and implemented parts of it. I also think I might understand Kneser's approach, based on the Riemann mapping, starting with the entire analytic super-exponential, developed from the fixed point, c~=0.318+i*1.337. This is the inverse Abel function, where the Abel function is the slog developed from the fixed point c. \( \psi^{-1}(z)=\lim_{n\to\infty}( \exp^{\circ n}(c+c^{(z-n)}) \), at least I think that is what is meant by \( \psi^{-1}(z) \). I tried really hard to understand Jay's post, and this post, but there are gaps.... I tried to fill in the gaps, but I have a feeling I'm going down a slightly different approach here.
Is the real valued sexp at the real axis connected to the Schroeder equation, chi, by a complex 1-cyclic function, \( \theta(z)-z \)
\( \text{sexp}_e(z)=\psi^{-1}(\theta(z)) \)
Where \( \theta(z) \) is the Riemann mapping of contour of the limit of iterating the natural logarithm, starting with the real interval, 0..1 (which corresponds to sexp(z) from z=0 to z=1). Then
subtracting c, which is the fixed point, ~= 0.318+i*1.337? Then the Riemann mapping unit circle is scaled by 2*pi and converted to a repeating unit length at the real axis. Is this equation correct, and is this Kneser's approach?
\( \psi(z)=\lim_{n\to\infty}(
\log_c(\ln^{\circ n}(z)-c)+{n}) \)
\(
\theta(z)=\operatorname{RiemannMapping} \lim_{n\to\infty}(
\log_c(\ln^{\circ n}[0..1]-c)+{n}) \)
This equation is developed from the inverse sexp function, or the Abel function, \( \psi(z) \). Assuming the equation is correct, here is my rough calculation of the resulting limit, whose contour is what I think is to be Riemann mapped, and then used to generate \( \theta(z) \). For double precision numbers >0 and <1, the real part of the contour never exceeds 1.78, so most of the contour graph is only for values of z super-exponentially close to zero or one, but I thought I would include the entire contour for completeness. Also, I removed the scaling, and added the next iteration, which is developed from sexp(z), from z=1 to z=2, which shows that \( \theta(z+1)=\theta(z)+1 \), and that \( \theta(z)-z \) is 1-cyclic. Contour values from this graph should have imaginary values of zero for \( \psi^{-1}(z) \). Also, the outermost \( \log_c \) for the contour is multi-valued, and as n increases, sometimes multiples of -2pi*i have to be added so that the contour equation converges. The other values correspond to other identical parts of the \( \psi{-1}(z) \) function, which repeats infinitely many times in the complex plane.
Once the Riemann mapping is generated, and converted to a unit length at the real axis, to generate \( \theta(z) \), I assume that the \( \theta(z) \) goes to a small constant (the center of the Riemann mapping unit circle), as the imaginary axis grows, and has a singularities at unit intervals at the real axis. Presumably, the singularity in the Riemann mapping is at one point on the edge of the unit circle, but not inside the unit circle.
The contour is valid before and after the Riemann mapping. But after the Riemann mapping, it would seem that \( \theta(z)-z \) is a 1-cyclic function, and it is the only analytic 1-cyclic function that converts from \( \psi^{-1}(z) \) to the sexp_e such that \( \theta(z)-z \) will be a constant plus the sum of a complex fourier series, with complex coefficients, of \( A_n*e^{(+2z\pi n {i})} \). Because of the unique Riemann mapping, all Fourier terms of the form \( B_n*e^{(-2z\pi n {i})} \) are zero. Because the Riemann mapping to drive these terms to zero is unique, there is only one sexp solution approaching the fixed point as i increases, with the complex fourier series decaying to zero as i increases.
Given that Kouznetsov's solution also assumes the limiting value of f as imaginary increases, it would seem that complex Fourier analysis would allow one to prove why the solution is unique, why Kouznetsov's solution is equal to Kneser's solution, and perhaps even why Kouznetsov's iterated Cauchy integral converges.
- Sheldon
Posts: 368
Threads: 44
Joined: Sep 2009
What do you mean by \( \ln_f \)?
Posts: 684
Threads: 24
Joined: Oct 2008
01/25/2010, 07:42 AM
(This post was last modified: 01/25/2010, 07:50 AM by sheldonison.)
(01/25/2010, 06:35 AM)mike3 Wrote: What do you mean by \( \ln_f \)? Typo, that should've been \( \log_f \)
f=fixed point of "e", 0.318131505+1.337235701i, \( \ln(f)=f \)
\( \log_f(z)=\ln(z)/\ln(f))=(1/f)*\ln(z) \)
For clairty, I'm also changing this in the original reply I made.
- Sheldon
Posts: 368
Threads: 44
Joined: Sep 2009
(01/25/2010, 07:42 AM)sheldonison Wrote: (01/25/2010, 06:35 AM)mike3 Wrote: What do you mean by \( \ln_f \)? Typo, that should've been \( \log_f \)
f=fixed point of "e", 0.318131505+1.337235701i, \( \ln(f)=f \)
\( \log_f(z)=\ln(z)/\ln(f))=(1/f)*\ln(z) \)
For clairty, I'm also changing this in the original reply I made.
- Sheldon
Ah. So how would you use the given equation to calculate \( \theta(z) \)?
Posts: 684
Threads: 24
Joined: Oct 2008
01/26/2010, 01:22 PM
(This post was last modified: 01/28/2010, 11:25 AM by sheldonison.)
(01/26/2010, 06:24 AM)mike3 Wrote: Ah. So how would you use the given equation to calculate \( \theta(z) \)? Update: fixed notation to use the Abel function, \( \psi \) which is the inverse of the super-exponential developed from the fixed point "c".
Well, again, I was sort of trying to follow the graphs that Jay made, and Henryk's description of Kneser's approach, and it was not quite accessible to me. It appears they were iterating the logarithm on a unit length, [0..1], which is the first step in the limit towards the inverse super-exponential, or the Abel function, \( \psi \). The sequence of logarithms converges towards the fixed point, so the logical next step is the equation for \( \psi \), or the inverse super-exponential developed from the complex fixed point c. And the equation I got for \( \psi \) is
\( \psi(z)=\lim_{n\to\infty}(\log_c(\ln^{\circ n}(z)-c)+n) \)
Of course, I'm assuming that \( \psi^{-1} \) is the entire sexp_e developed from the fixed point, developed from the approximation for small (large negative) values of z, \( \psi{-1}(z)=c+c^{z} \), which has complex values at the real axis. But I don't think that matches what Kneser is doing, but I think it must be closely related. Then I graphed the contour of \( \psi([0..1]) \), and I was wondering if that was indeed the contour that Kneser Reimann mapped, \( \psi[0..1] \). But I couldn't tell from Jay's posts, or Henryk's post. It seemed more complicated than that. Anyway, If you already had the sexp_e to start with than
\( \theta(z)=\psi(\text{sexp}_e(z)) \)
The 1-cyclic function is \( \theta(z)-z \). Plugging into the equation for \( \psi(z) \).
\( \theta(z)=\lim_{n\to\infty}(
\log_c(\ln^{\circ n}(\text{sexp}_e(z))-c)+{n}) \)
Of course, we don't know what sexp_e is, so we start with the contour generated over a unit length from [0..1]. A linear approximation for sexp_e is probably a reasonbly starting point, since it has a continuous first derivative. Then do a Riemann mapping of the contour of \( \psi^{-1} \) of the unit length, to a unit circle. After the Riemann mapping, unroll the the unit circle, divide by 2*pi, and we're back to a unit length. Only now \( \theta(z)-z \) is a 1-cyclic function that can be expressed as the unique sum of a complex fourier series exponentially decaying to a constant as i increases. I gather the computational difficulty comes in actually doing the Riemann mapping...
\( \theta(z)=\text{Riemann Mapping} \lim_{n\to\infty}(
\log_c(\ln^{\circ n}([0..1])-c)+{n}) \)
- Shel
|