The numerology of aliquot sums and perfect numbers
The numerology of the Pythagorean sages among the old yavana-s is one of the foundations of science and mathematics as we know it. One remarkable class of numbers which they discovered were the perfect numbers — teleios as they termed it. What are these numbers? Let be a number and be all its proper divisors, i.e. those divisors of , which are less than . Then we can define an arithmetic function known as the aliquot sum thus,
For example, let us consider the number 10. Its divisors are 1, 2,5,10. Its proper divisors are 1, 2, 5. Hence, . Now, if then is called a perfect number (termed pūrṇāṅka in Jagannātha’s Sanskrit edition of Euclid). From Figure 1, we can see that the numbers 6 and 28 are perfect numbers. The Pythagorean interest in them also becomes apparent from the fact that these numbers are associated with certain natural periodicities that have an old Indo-European significance. This becomes clear from their occurrence even in old Hindu tradition. The number 6 is associated with the 6 seasons in brāhmaṇa-s like the Śatapatha-brāhmaṇa and encoded into the śrauta altar. Similarly 28 is also encoded into the altar according to the Śatapatha-brāhmaṇa and corresponds to the count of the nakṣatra-s or days of the lunar month. The number 28 is also used in Vaidika tradition as one of the prescribed counts for the japa of the Savitṛ gāyatrī
Returning to the arithmetic, with our above example of , we can see that . Such numbers are called abundant numbers. Now, a subset of the proper divisors of 20 add up to 20 (Figure 1). Hence, it is also a semi-perfect or a pseudo-perfect number. There is a rare set of abundant numbers for which no subset of their proper divisors add up to them; they are known as weird numbers. For example, ; hence, it is an abundant number. However, no subset of its proper divisors, 1, 2, 5, 7, 10, 14, 35 add up to 70. Hence, 70 is a weird number. The next weird number is 836. Thus, due their rarity majority of abundant numbers are semiperfect numbers.
One of the high points of yavana brilliance was the discovery of a general formula for perfect numbers. In Book 9, proposition 36 Euclid states:
“ean apo monados hoposoioun arithmoi hexēs ektethōsin en tē diplasioni analogia, heōs hou ho sumpas suntetheis prōtos genētai, kai ho sumpas epi ton eskhaton pollaplasiastheis poiē tina, ho genomenos teleios estai.”
If as many numbers as we please beginning from an unit be set out continuously in double proportion, until the sum of all becomes prime, and if the sum multiplied into the last make some number, the product will be perfect.
In modern terms we can lay it out thus: If the sum of the series of the powers of 2, where the power are integers , is a prime number then the product of that prime with the last power of 2 in the series is a perfect number:
If is a prime then is a perfect number. Let us take the first few examples.
Not a prime
Not a prime
After this point the perfect numbers become much less frequent and large.
Another way of expressing this is thus: Given a prime number if is also a prime then is a perfect number. The corresponding prime numbers are today famous as the Mersenne primes. For all , the following yield : 2, 3, 5, 7, 13, 17, 19, 31, 61, 89
The corresponding are:
The corresponding perfect numbers are:
These primes start getting huge and rare rapidly. We computed the above in a few seconds with a modern programming language and computer. However, some of them mark historic feats of arithmetic computation by the unaided human brain in the pre-computer era. For instance, Leonhard Euler computed the 8th (10 digits) and the corresponding perfect number. The Russian village mathematician I.M. Pervushin went further by computing the 9th of these numbers ( digits). With the Great Internet Mersenne Prime Search we now have 50 of them and the corresponding perfect numbers. This also helps us to see why the always end in 1 or 7 and the corresponding perfect numbers end in 6 or 28. That the above are the only even perfect numbers was established by Euler. Unlikely as they seem, there has been no formal proof to show that no odd perfect numbers exist. The modern efforts also show that the perfect numbers are rarer than what the yavana sages thought them to be. The old Pythagorean numerologist Nicomachus states:
“It comes about that even as fair and excellent things are few and easily numerated, while ugly and vile ones are widespread, so also the abundant and deficient numbers are found in great multitude and irregularly placed – for the method of their discovery is irregular – but the perfect numbers are easily enumerated and arranged with suitable order; for only one is found among the units, 6, only one among the tens, 28, and a third in the rank of the hundreds, 496 alone, and a fourth within the limits of the thousands, that is, below ten thousand, 8128.”
The above statement suggests that Nicomachus might have thought that for each decade there is one perfect number. Alternatively, he might have simply stopped counting with the greatest perfect number he knew. That it was the former is strengthened by Platonic siddha (to use Gregory Shaw’s term), Iamblichus, even more explicitly claiming that there is one perfect number per decade. That, however, is clearly wrong as the next perfect number 33550336 is in the crores. In any case the modern confirmation of the real rarity of perfect numbers vindicates the philosophical analogy drawn from these numbers by the yavana siddha — the scarcity of truly perfect things as opposed to the profusion of the supernumerary and the deficient. With respect to the easy enumeration of the perfect numbers Nicomachus and Iamblichus likely meant Euclid’s proposition. This rarity of perfect numbers indicates that the sum of the reciprocals should converge to a constant:
It would be immensely remarkable if this constant turns up somewhere in nature.
In any case, this provides a bound for the distribution of abundant numbers and thus, brings us to the point of whether Nicomachus’ statement regarding the abundant numbers is really so? This was what we set out to investigate in our youth armed with a rather meager arithmetic knowledge and a computer. Computing the first few abundant numbers yields a sequence like below:
12, 18, 20, 24, 30, 36, 40, 42, 48, 54, 56, 60, 66, 70, 72, 78, 80, 84, 88, 90, 96, 100, 102, 104, 108, 112, 114, 120, 126, 132, 138, 140, 144, 150, 156, 160, 162, 168, 174, 176, 180, 186, 192, 196, 198, 200, 204, 208, 210, 216…
A few rules become immediately apparent:
1) if is a perfect number then is an abundant number for all . Thus, initiates a line of abundant numbers 12, 18, 24, 30…. initiates 56, 84, 112…
2) We notice that beyond these a few abundant numbers emerge in the interstices e.g. 20, 88 etc. A notable class of these emergent abundant numbers are of the form , where is the odd prime and . The first few of these are tabulated below.
These are all shown as red points in panel 1 of Figure 2. One notices that for each jump of we get a jump in these abundant numbers. Of these, some like 12, 56 and 992 are already accounted for as doubles of the perfect numbers 6, 28 and 496 but the rest are distinct. Their multiples, as with the perfect numbers, are also further abundant numbers.
3) A few further abundant numbers emerge in the interstices, which have a more complex description. They found new lineages of abundant numbers as their multiples continue to be abundant numbers. The first of these is . 70 is the product of successive odd primes 5, 7 with 2 leaving out 3, which would yield already accounted-for abundant numbers in a product with 2 via the perfect number 6. Another such is . Here 5 is left out because that will again result in already accounted-for abundants in a product with 7. Likewise, leaving out 7, we have . For further numbers in this category we need the next power of 2, e.g. .
4) Till odd abundant numbers are very rare. The is the first odd abundant number and till there are only 391 odd abundant numbers as opposed to 49090 even abundant numbers. Of these first 391 odd ones, 387 are divisible by 15. The remaining 4 which are not namely 81081, 153153, 171171, 189189 are divisible by 9009. All odd abundant numbers necessarily need to have the product of at least 3 distinct odd primes as a divisors. As with other abundant numbers, the multiples of odd abundant numbers are also abundant numbers. Based on the separation between the odd abundant numbers, we observed that the first of them 945 is the first of a lineage of odd abundant numbers, which are generated by the formula: , where . Similarly, we observed that a related formula , where produces a continuous run of 193 odd abundant numbers starting from 3465. After and respectively these formulae do not necessarily produce abundant numbers; nevertheless numbers emerging from these rules continue to remain enriched in odd abundant numbers. Of course, the other emergent odd abundant numbers might have further less-apparent rules. In any case, it appears the difference between two successive odd abundant numbers is most of the times divisible by and always divisible by the perfect number 6.
Thus, unlike what the yavana siddha-s thought, the abundant numbers are not entirely unruly. They have complex patterns, but can be described by some rules. The biggest bulk of them are children of perfect numbers — like mortal descendants of the immortal gods. This leads us to the growth of the sequence of abundant numbers, which exhibits striking regularity (Figure 2).
Given that the sum of the reciprocals of the perfect numbers converges to a constant and given what we have seen above regarding the abundant numbers, they should grow at the upper bound by the equation (the dark red line panel 1 of Figure 2). However, since the progeny of perfect numbers are not the only abundant numbers and several emerge by the other rules mentioned above should grow at a lower rate. During our initial foray into abundant numbers, by empirical examination we thought this growth rate might be exactly 4 (Figure 2, panel 1, light blue line). However, with better computers and computation of 49481 abundant numbers, we later realized that the rate was actually slightly more than 4 (Figure 2, panel 1, dark green line). Another way of visualizing the same is by defining the density of these numbers as the ratio of the number of abundant numbers to a certain number to ,
This is shown in Figure 2, panel 2, where we see converging to a value between 0.248 and 0.247 (red and blue lines). This has been a topic of deep investigation in modern mathematics, starting with the likes of Sarvadaman Chowla and Paul Erdős, and indeed the above has proven to be the approximate bound of the density of abundant numbers. We were quite pleased to have semi-empirically arrived at it for ourselves. Thus, even as the odd and even numbers are distributed in a 1:1 ratio, after initial jitters the ratio of deficient to abundant numbers converges to a constant of with new lineages of abundant numbers constantly emerging in the interstices of the established ones to maintain a linear growth at this rate.
The basic features of the aliquot sequences
The process of obtaining aliquot sums of a number can be applied iteratively (the aliquot map): . The sequence is called the aliquot sequence of . Thus, if we start with we get the sequence: 20, 22, 14, 10, 8, 7, 1, 0. More generally, by computing aliquot sequences for all numbers from 1:1000, we observe the following:
1) The number 1 converges in two steps to 0.
2) All primes converge in 3 steps, .
3) Perfect numbers converge in 1 step to themselves.
4) Certain numbers converge to a perfect number. E.g. 25, 95, 119, 143, 417, 445, 565, 675, 685, 783, 909, 913… converge to 6; 608, 650, 652, 790 converge 496. We observe that these numbers have given the appellation “aspiring numbers”.
5) Notably, 220 converges to 284 and 284 converges to 220. . Thus, if from whatever starting point if one reaches 284 or 220 one cycles between them. Such a pair is termed a pair of amicable numbers and was probably known to Platonists as indicated by the heathen Sabian from Harran, Thabit ibn Kurra’s knowledge of these numbers. Thabit discovered a rule to produce a lineage of these numbers, likely drawing on the rule to generate even perfect numbers: Let and . If are primes then we can compose the amicable pair as . For we get the pair 220, 284; produces 17296, 18416; produces the pair 9363584, 9437056. However, there are many more amicable numbers between these pairs which cannot be captured by this rule. For example, one can computationally show that 1184, 1210 are an amicable pair. Long after the days of Thabit, starting from Euler down to our times further rules have been found to capture more amicable pairs.
6) If perfect numbers are auto-cycles under the aliquot map, the amicable numbers can be considered 2-cycles. While we do not know if all higher cycles exist under the aliquot map, we can computationally find some of them. For example the pentad 12496, 14288, 15472, 14536, 14264 constitute a 5-cycle. Any number in this pentad will cycle through these values. Even more remarkable is this sequence 14316, 19116, 31704, 47616, 83328, 177792, 295488, 629072, 589786, 294896, 358336, 418904, 366556, 274924, 275444, 243760, 376736, 381028, 285778, 152990, 122410, 97946, 48976, 45946, 22976, 22744, 19916, 17716. These constitute a 28-cycle and any number in this 28-ad will cycle between these values. The lowest member of a 2- and greater cycle is always an abundant number.
7) Still other numbers will eventually converge to an odd prime and via that prime reach 0. There is a remarkable pattern in terms of the preference of the prime via which convergence occurs (Figure 3). 43 is the most preferred prime for convergence for (any links to Heegner numbers or twin primes?). We are not clear as to whether this trend survives with increasing .
8) If we leave out the starting will every other integer be reached by rest of the aliquot trajectory of some ? We can easily provide the answer to this question as no. For example, . You cannot have two 1s as proper divisors, hence 2 can never be reached from any other number under the aliquot map. Similarly, . But each of these sums leading to 5 cannot be a because we would leave out divisors 2 and 1. Thus, 5 can never be reached from any other number under the aliquot map. On the other hand, a number (where is a prime number) is never untouchable because . Similarly, a number (where is an odd prime) is never untouchable because . Nevertheless, some numbers are not reached easily until one of the above configurations is encountered for the first time. For example, contains 38=37+1 for the first time. Similarly, contains 48=47+1 for the first time. But the number 52 lies a in sweet spot as as neither 51=52-1 nor 49=52-3 are primes and is untouchable from any other number under the aliquot map. These are some obvious examples of untouchability and touchability; Erdős has shown that ultimately there are an of cases.
Patterns in the convergence length of aliquot sequences
Are these the only convergence patterns seen in the aliquot sequences? Is there any further pattern to the number of iterations needed to converge? To address these questions we can define a further sequence where each element is the number of iterations taken by the integer to converge. The plot of for the first 1000 elements is shown in Figure 4.
The picture is quite remarkable — the majority of values are rather pedestrian but from time to time there are huge eruptive values. The perfect numbers ( ), 1 and the primitive aspiring numbers, e.g. 25 ( ), and primes ( ) form the lowest values of . might also be reached by certain secondary aspiring numbers like 95 or 119. Of these the primes are the most common. It is quite obvious that even numbers have significantly longer convergence paths than odd numbers ( ; Figure 5). We also observed that on an average the abundant numbers have significantly longer convergence paths than the remaining numbers ( for ). This holds even after the primes are removed from the non-abundant numbers ( for ). Similarly, if we compare even abundant numbers and even non-abundant numbers, which occur in the roughly similar counts, the abundant even numbers still have significantly longer trajectories the non-abundant even numbers under the aliquot map ( ). Due to the relative rarity of odd abundant numbers, these effects are likely to persist beyond . The basic statistics for the convergence trajectory lengths for different types of numbers from are tabulated below and summed up in Figure 5.
During our initial foray into the first 200 terms of this sequence it appeared that certain integers, like 138, never converged under the aliquot map. Going on till more such terms appeared; however, better computation showed us that the aliquot sequences for some of these indeed converge after a large number of steps after growing to monstrous values. It was then that we learned that the mathematicians have been studying this problem quite a bit computationally starting with Derrick Lehmer-II (the son of the father-son pair of arithmeticians), one of the doyens of computational mathematics. He was the first who, back in his days after what was a tough fight, showed that indeed converges. For us, the hardest nut to crack for was . It simply would not complete in our laptop at home even after running it for over a day with an efficient divisor finding function. Nor did it complete on our primary work station used for most of routine computations. Then we had to bring out our Indrāstra, a 120 core machine that we use for big things. For that we had to first write a multi-threaded version, which was run on 100 cores of this machine and it completed overnight. We term these large as the monster numbers. For there are 18 of them, which come in the below-tabulated families, all of which are abundant numbers descending from 6. After the first few terms the members of each family follow the same trajectory to convergence. The final prime via which they converge is termed
One notices that families and have a higher-order relationship because they converge via 59. However, we keep them as different families for, barring the last 8 elements of , they followed very different trajectories (Figure 6). Their is the 6th most common for , with 45 numbers converging via it (Figure 3). However, 12 of them being monster values it is one over-represented in the aliquot sequences attaining monster values. In contrast, the most common in this range 43 (Figure 3) has only a single monster value (936) in this range. shows the most monstrous behavior in this range (Figure 6). It reaches a maximum value of the behemoth number: 3463982260143725017429794136098072146586526240388. I felt please on reaching this number and experiencing the magnificence of for myself. After that it remains in the high range for while with two prominent peaks before converging after 749 steps. Such a behavior, with a multiplicity of peaks before convergence, is also seen in the other monstrous families (Figure 6).
We wondered if there was any features of these sequences, which allowed their eruptive growth before finally converging. We observed the following:
1) They start as even abundant numbers, which means that they are statistically in the general group of numbers which are going to have longer trajectories to convergence.
2) Their trajectories tend to remain even, which again increases their probability of having a long convergence path. By specifically studying the first major eruptive case , we see that it has a convergence trajectory of length 179. It remains even from to . Thus, remaining even is a key factor for a long convergence path.
3) The next key feature observation regarding their growth came from examining the behavior of the first two cases which show above average behavior. These are and . First, both are even abundant numbers satisfying the above condition. Then we note their behavior under the aliquot map:
We see that both of them have a run of 7 continuous abundant numbers, which are multiples of the perfect number 6. This allows a monotonic increase. When they lose this divisor 6 in the 8th iteration they become odd and start falling rapidly and converge. We then looked if this pattern might play out in a truly eruptive example by taking the case of . We see that from to , then to and again to each number is a multiple of the perfect number 6. The repeated generation of a multiple of a perfect number under the aliquot map means there can be stretches of unhindered growth with one abundant number succeeding the previous one until that perfect number is lost among the divisors. It is the final loss of 6 and then 2 among the divisors that allows to converge. In making these observations we had recapitulated Lehmer’s key findings in the context of the growth of the most mysterious class of aliquot sequences we shall talk about next.
There are 12 , for , which could not be determined because the corresponding never converged in our computation (the blue dots in Figure 4). A search on the Internet reveals that deep computational efforts by various investigators have not seen an end to these sequences. They belong to 5 families which are known after Lehmer-II who first studied them. Within each family (listed below), after the first few terms the evolution converges to the same trajectory (Figure 7). These families are:
Figure 7. : dark red; : dark blue; : dark green; : hot pink; : cyan
All these numbers are abundant numbers descending from 6, which is a tendency they share with the above monster numbers that eventually converge. The frequencies of their largest prime factors are tabulated below:
Notably, in this set except , 4 of the five families have a member with 23 as the highest prime factor. It is now conjectured by mathematicians that these define the founding members of a class of aliquot sequences that would never converge but eventually shoot off to infinity. Examining the eruptive growth of these sequences, we found that the same features that cause the eruptive growth of the converging but monster hold true here too. Thus, if we consider the first case , which starts with a double of the first monster , we notice that it starts off well with the perfect number 6 as the divisor for the first 5 terms successive terms — to . But then the sequence loses 6 in the 6th term a stutters for a few terms even going down from the high point it reached. But then in it acquires a new perfect number divisor 28, which drives the abundance of successive terms until the it reaches the giant 28 digit . Lehmer’s study was the first which noted the role of such persistent divisors in driving up the aliquot sequences. Accordingly, the grand old man of arithmeticians, Richard Guy, called these ‘drivers’ and in addition to 2 and the perfect numbers defined a few more of such. What seems to happen is that in these apparently non-converging aliquot sequences the drivers drive the growth to some large number after which they are lost. But before the sequence can crash by attaining an odd number, a new driver emerges and pushes up its growth again. Thus, these sequences might not converge. On the other hand the case of the monstrous with tells us that even after many such runs there is some chance of hitting a number, which brings you crashing down though its even almost all the way through. Whatever the case, the aliquot sequence of these non-converging numbers tell us that there are several more numbers in their trajectories, which would similarly not converge. Hence, non-convergence will be encountered frequently as the starting numbers get large.
In conclusion, following along the lines of the ancients, we find that the aliquot sequences offer us some philosophical lessons and mysteries:
1) There is the predictable: The primes and perfects display very predictable and uniform behavior.
2) There is the somewhat predictable: like the amicable numbers, for some of which a rule was found as early as Thabit ibn Kurra. But not all of them can be captured by an easy rule.
3) There is the statistically tendency: We know that an even or abundant number is more likely to have a longer path to convergence though it is not clear at all if any rule can say which of them might be monstrous and which rather common place.
4) There is the unpredictable and mysterious: Why are some numbers “aspiring numbers” which converge to a perfect number? Is there any rule to describe them? Why is 43 the preferred prime for convergence at least for . Why do some numbers not converge? Can we predict which numbers will not converge? Do they really not converge at all? As far as I know these remain unanswered. Thus, a very simple arithmetic process gives rise to a whole range of unexpected behavior.