Functional Analysis behind a Family of Multidimensional Continued Fractions: Part I
Ilya Amburg, Thomas Garrity

TL;DR
This paper explores the functional analysis of transfer operators linked to triangle partition maps, a family of multidimensional continued fraction algorithms, establishing their classes and properties.
Contribution
It introduces the classification of transfer operators for triangle partition maps and analyzes their origins and compositions, laying groundwork for further spectral analysis.
Findings
Two classes of transfer operators identified
Theorems on the origin of these classes established
Foundations for spectral analysis in Part II laid out
Abstract
Triangle partition maps form a family that includes many, if not most, well-known multidimensional continued fraction algorithms. This paper begins the exploration of the functional analysis behind the transfer operator of each of these maps. We show that triangle partition maps give rise to two classes of transfer operators and present theorems regarding the origin of these classes; we also present related theorems on the form of transfer operators arising from compositions of triangle partition maps. In the next paper, Part II, we will find eigenfunctions of eigenvalue 1 for transfer operators associated with select triangle partition maps on specified Banach spaces and then proceed to prove that the transfer operators, viewed as acting on one-dimensional families of Hilbert spaces, associated with select triangle partition maps are nuclear of trace class zero. We will finish in part…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsMathematical Dynamics and Fractals · Fractional Differential Equations Solutions · Functional Equations Stability Results
Functional Analysis behind a Family of Multidimensional Continued Fractions:
Part I
Ilya Amburg111Center for Applied Mathematics, Cornell University, Ithaca, NY 14853; [email protected] and Thomas Garrity222Williams College Department of Mathematics, Williamstown, MA 01267; [email protected]
Abstract
Triangle partition maps form a family that includes many, if not most, well-known multidimensional continued fraction algorithms. This paper begins the exploration of the functional analysis behind the transfer operator of each of these maps. We show that triangle partition maps give rise to two classes of transfer operators and present theorems regarding the origin of these classes; we also present related theorems on the form of transfer operators arising from compositions of triangle partition maps. In the next paper, Part II, we will find eigenfunctions of eigenvalue 1 for transfer operators associated with select triangle partition maps on specified Banach spaces and then proceed to prove that the transfer operators, viewed as acting on one-dimensional families of Hilbert spaces, associated with select triangle partition maps are nuclear of trace class zero. We will finish in part II by deriving Gauss-Kuzmin distributions associated with select triangle partition maps.
1 Introduction
Using continued fractions to understand the structure of real numbers has been going on for many centuries. One of the most interesting properties of continued fractions goes back to work of Lagrange in 1789, who showed that a real number has an eventually periodic continued fraction expansion if and only if the number is a quadratic irrational. This is in analog to the fact that a real number has an eventually periodic decimal expansion if and only if the number is rational. As most numbers are neither rational nor quadratics, this leads to the following first of a whole slew of open questions:
The Hermite Problem: Find methods for writing numbers that reflect special algebraic properties.
Even for cubic irrationals, this problem is still quite open. Attempts to solve it are usually called multidimensional continued fractions. There are many, seemingly too many. (Multidimensional continued fractions have many other uses.)
Almost all multidimensional continued fraction algorithms associate to a point in a triangle a sequence of integers. For most, eventual periodicity of this sequence means that the and are in a number field of degree at most three. What is difficult is the converse.
Most multidimensional continued fraction algorithms can be interpreted as iterative systems on a triangle or, in higher dimensions, on an appropriate simplex. For a review of many types of multidimensional continued fraction algorithms, see Schweiger’s Multidimensional Continued Fractions [38] and Karpenkov’s Geometry of Continued Fractions [25]. Thinking of these as iterative systems, it is natural to investigate them as dynamical systems, which in turn will lead to the study of transfer operators (for background on transfer operators, see Baladi [7]). This has long been done for traditional continued fractions, using the Gauss map acting on the unit interval. For background on this, see Hensley[18], Iosifescu and Kraaikamp [21], Kesseböhmer, Munday and Stratmann [26], Khinchin [27], Rockett and Szusz [36] and Schweiger [37]. In [17], this was done for the triangle map, where it was shown that much of the work of Mayer and Roepstorff [32, 33] has nontrivial analogs for the triangle map. Building on an earlier version of [17] , there is also the interesting work of Bonanno, Del Vigna and Munday [9]. Further, there is the work of Fourgeron and Skripchenko [15] on the Lyapunov exponents of the triangle map. The goal of this paper is to see which of these analogs hold for various members of the family of triangle partition maps.
Triangle partition maps [12, 13, 25] are a family of multidimensional continued fraction algorithms that include (when combinations are allowed) many, if not most, well-known multidimensional continued fraction algorithms, which is why they are a natural class of algorithms to study. These maps are reviewed in section 2.
In section 3, we will start to look at transfer operators associated with triangle partition maps. For the Gauss map, we know that the transfer operator is
[TABLE]
and for the triangle map, the transfer operator is the somewhat similar looking,
[TABLE]
In section 4, we will show for precisely half of the triangle partition maps that the transfer operator will have the form
[TABLE]
where will be a polynomial linear in (with the possibility of a term) and linear in and . We will say that such transfer operators have polynomial behavior. We will then show for the remaining half of triangle partition maps that their transfer operators have drastically different behavior. We will see that their transfer operators are of the general form
[TABLE]
where is exponential in . Such triangle partition algorithms are naturally enough said to have non-polynomial behavior. It is here where this paper (Part I) ends.
But this paper sets up, more importantly, the needed machinery for the next paper “Functional Analysis behind a Family of Multidimensional Continued Fractions: Part II” [4] In this next paper we turn to the spectrum of the transfer operators. Here we need to be concerned with which vector spaces of functions we are considering. We will see in section two of that paper that for some triangle partition maps there are natural Banach spaces of functions, that for of these the largest eigenvalue of the associated transfer operator is one, and that for some of these the dimension of the corresponding eigenspace is one. Then we will see in section three that a number of the triangle partition maps’ transfer operators are nuclear operator of trace class zero when viewed as acting on naturally occurring one-dimensional families of Hilbert spaces. In section four we present the Gauss-Kuzmin statistics for select triangle partition maps. In the conclusion of this next paper we discuss future directions for work.
We would like to thank the referees for their suggestions.
2 Background on Triangle Partition Maps
As this is background, much of this section is similar to certain sections in [3, 6, 12, 13, 16, 25, 39].
2.1 The Triangle Map
We start with the triangle map, from which, in the next section, the 216 triangle partition maps are constructed. For a further explanation of the triangle map, see [16, 6].
Partition the triangle
[TABLE]
into subtriangles
[TABLE]
for each non-negative integer . The partitioning is represented in the following diagram:
(0,0)(1,0)(1,1)\triangle_{0}$$\triangle_{1}$$\triangle_{2}$$\triangle_{3}$$\triangle_{4}
Each subtriangle has vertices The triangle map is defined by setting, for any ,
[TABLE]
This is an analog of the traditional Gauss map for continued fractions.
And as with the Gauss map and continued fractions, it is useful to translate this into matrices. We have the projection map
[TABLE]
mapping rays in to points in . This allows us to associate to the triangle in the cone in (which we will also denote by )
[TABLE]
The cone is the positive span of the column vectors of the matrix
[TABLE]
Define matrices
[TABLE]
These matrices occur naturally since the cone corresponding to the subtriangles is the positive span of the columns of , which by an abuse of notation, we will often write as
[TABLE]
Then the triangle map can be shown to be equal to:
[TABLE]
if
This enables us to associate the triangle sequence to a point by letting equal if . If for any we have , the sequence terminates.
2.2 Triangle Partition Maps: TRIP maps
The above triangle map depends on an initial choice of three vertices for . If we permute these vertices both before and after applying the matrices and , we create different multidimensional continued fraction algorithms. This is at the heart of the definition for the 216 triangle partition maps.
More precisely, we will allow permutations of the initial vertices by some by some after applying and by some after applying
This leads us to define the matrices
[TABLE]
for every
In particular, we denote the permutation matrices as follows:
e=\left(\begin{array}[]{ccc}1&0&0\\ 0&1&0\\ 0&0&1\\ \end{array}\right), (12)=\left(\begin{array}[]{ccc}0&1&0\\ 1&0&0\\ 0&0&1\\ \end{array}\right), (13)=\left(\begin{array}[]{ccc}0&0&1\\ 0&1&0\\ 1&0&0\\ \end{array}\right),
(23)=\left(\begin{array}[]{ccc}1&0&0\\ 0&0&1\\ 0&1&0\\ \end{array}\right), (123)=\left(\begin{array}[]{ccc}0&1&0\\ 0&0&1\\ 1&0&0\\ \end{array}\right), and (132)=\left(\begin{array}[]{ccc}0&0&1\\ 1&0&0\\ 0&1&0\\ \end{array}\right).
Let be the image of under the action of Thus, with an abuse of notation, we have
[TABLE]
Then define a triangle partition map as where
[TABLE]
if
As has six elements, we have different triangle partition maps. Of course, the original triangle map corresponds to Appendix A in Part II gives the explicit forms of for of the maps. One of the goals of this paper is to explain why we do not list in an easy-to-read table the forms of for the remaining maps.
We obtain an even larger family of maps by allowing compositions of triangle partition maps [12]. As an example, we might perform the first subdivision of using the second subdivision using and so on. We can represent such compositions as where, of course, each subscript is short for a permutation triplet.
Again, many, if not most, multidimensional continued fraction algorithms can be put into this language [12], which is why these are natural maps to consider. To get a feel of these maps, we consider two examples in a bit more detail.
2.3 The Case
This is simply an example of what we will later call a polynomial-behavior triangle partition map. It is in and of itself not that significant. We have
[TABLE]
We want to see what the subtriangles look like. We have
[TABLE]
Pictorially, we have
(0,0)(1,0)(1,1)\triangle_{0}((23),(23),(23))$$\triangle_{1}((23),(23),(23))$$\triangle_{2}((23),(23),(23))$$\triangle_{3}((23),(23),(23))$$\triangle_{4}((23),(23),(23))
In the plane , the vertices of each are
[TABLE]
Now, for any , we have
[TABLE]
in which case,
[TABLE]
a formula that we will need later.
2.4 The Case
Here we present another example of a TRIP map. We will see that its behavior is quite different than that of TRIP maps and a difference that will be explained in section 4.
We first need to calculate and . is easy, as
[TABLE]
For the other matrix, we have
[TABLE]
We need to find a clean formula for We will need the Fibonacci numbers, which are given by with initial terms Note that this means that
Proposition 1
[TABLE]
Proof By induction. For , we have that is the identity. But for , we have
[TABLE]
Then in general we have
[TABLE]
as desired.
Now to find the subtriangles .
Proposition 2
[TABLE]
Proof By calculation this is
[TABLE]
Thus the subtriangle in the plane has vertices
[TABLE]
Pictorially, we have
(0,0)(1,0)(1,1)\triangle_{0}$$\triangle_{1}$$\triangle_{2}$$\triangle_{3}
We can now calculate the TRIP map, whose proof is just a somewhat painful calculation.
Proposition 3
For , the TRIP map
[TABLE]
has coordinate
[TABLE]
and -coordinate
[TABLE]
The fact that we had to use two lines to write out the TRIP map is not just chance, as we will see in section 4.
3 Transfer Operators
For general background on transfer operators and their importance, see Baladi [7].
In general, for a dynamical system corresponding to the map a transfer operator linearly maps functions from a vector space of functions on to another (not necessarily different) vector space of functions on For a more concrete definition, define a function Then the transfer operator acting on is defined by
[TABLE]
If is piece-wise differentiable, as for our maps, we choose , as is traditional and natural.
The spectrum for each transfer operator yields much information about the algorithm. One of the main motivations is that if the transfer operator has a leading eigenfunction with eigenvalue one with multiplicity one, we expect that this eigenfunction will give rise to the natural invariant measure for the map . This is exactly what Gauss conjectured for the transfer operator for the traditional continued fraction, conjectures only proven in the 1920s by Kuzmin. (Again, the goal of this paper is to start such investigations for TRIP maps, and hence for almost all multi-dimensional continued fraction algorithms.)
We define our transfer operators to be
[TABLE]
While a long calculation is required to arrive at the transfer operator corresponding to each the calculations are not difficult theoretically and have been carried out for each of the 216 triangle partition maps. As an example, as shown in [17],
[TABLE]
We have calculated the explicit form of
[TABLE]
for all As mentioned in the introduction, one of the main goals for this paper is to show that the family of 216 triangle partition maps naturally splits into two classes: 108 maps exhibiting what we term polynomial behavior, and 108 maps exhibiting non-polynomial behavior. Polynomial-behavior maps’ transfer operators have compact explicit forms while non-polynomial-behavior maps’ transfer operators have very lengthy and complex forms. The explicit forms for all of these maps are available online (see [2]).
As another example, let us return to the TRIP map of subsection 2.3. This has to have polynomial behavior, as otherwise it would be too cumbersome to express. We need to find the Jacobian of with any replaced by the first component of and any replaced by the second component of The Jacobian is:
[TABLE]
Substituting the first component of for the Jacobian becomes so that
[TABLE]
Hence
[TABLE]
Of course, we would now like to calculate the transfer operator for the TRIP map of subsection 2.4, but the calculations are a bit too cumbersome to actually write out. However, note that the transfer operator will be in terms of Fibonacci numbers, and hence exponential in terms of the parameter .
4 Polynomial and Non-Polynomial Behavior in Transfer Operators
4.1 A Clean Criterion
The transfer operator of the Gauss map is
[TABLE]
Note that the denominator is a polynomial that is the square of a polynomial linear in and linear in . The transfer operator for the original triangle map is
[TABLE]
Now the denominator is a polynomial that is the cube of a polynomial linear in and linear in and . In the previous section, we have shown that the transfer operator for the TRIP map also has this behavior. The goal of this section is to show that similar analogs hold for precisely 106 of the other triangle partition maps (and thus giving us polynomials behavior for 108 fo the maps), while the behavior of the transfer operator for the remaining 108 triangle maps is drastically different.
Definition 4
The transfer operator for a triangle partition map is said to have polynomial behavior if its Jacobian is a polynomial in the variables , and , with possibly a term appearing. Otherwise, we say that the transfer operator has non-polynomial behavior.
We will show that those transfer operators that have polynomial behavior actually have their Jacobians equal to the cube of a polynomial that is linear in and linear in and , with possibly a term appearing. In Appendix B of Part II, we explicitly list the 108 transfer operators that have polynomial behavior. Those without this behavior are far too messy to list, a fact that we find interesting.
We will need the following well-known lemma, whose proof is a calculation:
Lemma 5
Consider a vector , a matrix and the map
[TABLE]
Then
[TABLE]
Theorem 6
A triangle partition map has polynomial behavior if and only if the eigenvalues of the associated all have magnitude 1; it has non-polynomial-behavior otherwise.
Proof Our maps are all infinite-to-one, but, for each nonegative integer , the map is one-to-one and onto when restricted to the subtriangle :
[TABLE]
For the rest of this proof, we will fix , allowing us to consider the inverse
By the previous lemma and the definition of it follows that
[TABLE]
where
[TABLE]
Since is real and has magnitude (since all matrices involved in its calculation also have real determinants of magnitude ), this implies that a given is polynomial-behavior if and only if the first column of contains terms polynomial in
For simplicity, write where are the columns of Note that the inverse of a triangle partition map, can be written as
[TABLE]
Note that for any must be bounded as it must land back inside Since the inverse of each triangle partition map is bijective and since the choice of depends on which the original lies in, the first column of must depend on Thus, to show that the first column of depends on polynomially, it is sufficient to show that exhibits only polynomial dependence on – for then, by the above argument, the first row of must necessarily exhibit only polynomial dependence on
We have (suppressing the in some of the matrices)
[TABLE]
Let We can find the Jordan decomposition of ; i.e. we can write it as
[TABLE]
where is some invertible matrix of dimensions identical to those of and is the Jordan normal form of . Performing the Jordan decomposition for all 36 unique s we see that there exist only six unique s, namely
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
and
[TABLE]
where corresponds to a three-by-three square matrix with diagonal entries defined by the roots of ; similarly for
It can be shown that
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
and
[TABLE]
It is well-known that the diagonal elements of each are precisely the eigenvalues of the matrix from which it originated; further, through are diagonal and each contain at least one entry of magnitude greater than 1 on their respective diagonals. From this, it is clear that if (and hence ) has eigenvalues that all have magnitude 1, then the first column of depends polynomially on and hence the associated has polynomial behavior. Otherwise, the first column of depends exponentially on and hence the associated has non-polynomial behavior.
Assume has polynomial behavior in Then the first column of must depend strictly polynomially on From the above decomposition of we see that the only place where -dependance may enter the first column of is through ; hence, by the above argument the eigenvalues of all have magnitude 1. Now assume the eigenvalues of all have magnitude 1. Running above argument in reverse, we see that must necessarily have polynomial behavior.
Corollary 7
A triangle partition map is non-polynomial behavior if and only if the associated depends exponentially on
Corollary 8
A triangle partition map has polynomial behavior if and only if the associated is not diagonalizable.
Proof This follows immediately from the form of through in Theorem 6.
Corollary 9
Let and be three-by-three permutation matrices. If a triangle partition map has polynomial behavior, then so does the triangle partition map ; similarly, if have non-polynomial behavior, then so does .
Proof We need to compare and . We have
[TABLE]
Since conjugation will not change the Jordan canonical form of the relevant matrices, we are done.
Note that transfer operators with polynomial behavior are basically traditional zeta functions (though not dynamical zeta functions). In fact, we were initially tempted to call such transfer operators “zeta”-like.
Finally, what determines if a TRIP map has polynomial behavior depends only on . By calculation, the 108 polynomial behavior TRIP maps are
[TABLE]
5 On Polynomial Behavior in Combination Triangle Partition Maps
We are concerned with the form of the transfer operator where
[TABLE]
and is a finite composition of triangle partition maps defined by
[TABLE]
where each is a TRIP map for a permutation triplet .
If the denominator of is (non-trivially) polynomial in (also allowing for factors of ) then gives rise to a transfer operator that has polynomial behavior in and is itself polynomial behavior in otherwise, gives rise to a transfer operator that has non-polynomial behavior in and has itself non-polynomial behavior. Here we state and prove a theorem regarding the polynomial behavior of combo triangle partition maps, and present some corollaries. To avoid redundancy, we will refer to polynomial dependence on aside from factors of as polynomial dependence on or strictly polynomial dependence on
Theorem 10
A combination triangle partition map has polynomial behavior in if and only if the eigenvalues of the associated (the matrix corresponding to ) all have magnitude 1; it has non-polynomial behavior in otherwise.
Proof We have already shown explicitly that where
[TABLE]
and
[TABLE]
Since is real and has magnitude , this implies that has polynomial behavior in if and only if the first column of contains terms polynomial in . But now the argument will mirror the proof of theorem 6, since what matters is the nature of each , right down to there only being six unique Jordan canonical forms.
We have the following corollaries, whose proofs mimic the corollaries of theorem 6.
Corollary 11
A combination triangle partition map as defined in the above theorem is non-polynomial-behavior in if and only if the associated depends exponentially on
Corollary 12
A combination triangle partition map as defined in the above theorem is polynomial-behavior in if and only if the associated is not diagonalizable; it is non-polynomial-behavior otherwise.
6 Conclusion
While we find the distinction between polynomial versus non-polynomial behavior interesting, the real goal of this paper is to set-up the machinery for our paper “Functional Analysis behind a Family of Multidimensional Continued Fractions: Part II.” In this next paper we will find eigenfunctions of eigenvalue 1 for transfer operators associated with 18 polynomial-behavior triangle partition maps The formidably complex form of the non-polynomial-behavior transfer operators and lack of appropriate techniques makes finding leading eigenfunctions prohibitively difficult; however, finding leading eigenfunctions for the remaining 90 polynomial-behavior transfer operators appears a doable, though computationally-intensive, problem. We will then show the nuclearity of transfer operators, thought of as acting on one-dimensional families of Hilbert spaces, associated with 36 polynomial-behavior maps.
In the conclusion of the next paper, we further discuss many of natural and important questions that remain open.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] G. Alkauskas, Transfer operator for the Gauss’ continued fraction map. I. Structure of the eigenvalues and trace formulas, http://arxiv.org/pdf/1210.4083 v 6.pdf
- 2[2] I. Amburg, Explicit Forms for and Some Functional Analysis behind a Family of Multidimensional Continued Fractions – Triangle Partition Maps – and Their Associated Transfer Operators , Senior Thesis, Williams, 2014, available at https://drive.google.com/open?id=0B 7d Xp Xfy 8ngh V Hp WTGV 0bj M 3c Tg
- 3[3] I. Amburg, K. Dasaratha, L. Flapan, T. Garrity, C. Lee, C. Mihaila, N. Neumann-Chun, S. Peluse, M. Stoffregen, Stern Sequences for a Family of Multidimensional Continued Fractions: TRIP-Stern Sequences, Journal of Integer Sequences , vol. 20 (2017), Article 17.1.7.
- 4[4] I. Amburg and T. Garrity, Functional Analysis behind a Family of Multidimensional Continued Fractions: Part II, Publicationes Mathematicae Debrecen , to appear.
- 5[5] I. Antoniou and S. A. Shkarin, Analyticity of smooth eigenfunctions and spectral analysis of the Gauss map, Journal of Statistical Physics , 111, no. 1-2, (2003), pp. 355-369.
- 6[6] S. Assaf, L. Chen, T. Cheslack-Postava, B. Cooper, A. Diesl, T. Garrity, M. Lepinski and A. Schuyler, Dual Approach to Triangle Sequences: A Multidimensional Continued Fraction Algorithm, Integers , Vol. 5, (2005), A 08.
- 7[7] V. Baladi, Positive Transfer Operators and Decay of Correlation (Advanced Series in Nonlinear Dynamics), World Scientific, 2000.
- 8[8] S. Ben Ammou, C. Bonanno, I. Chouari, S. Isola, On the spectrum of the transfer operators of a one-parameter family with intermittency transition, http://lanl.arxiv.org/pdf/1506.02573 v 1.pdf
