Somos sequence

In mathematics, a Somos sequence is a sequence of numbers defined by a certain recurrence relation, described below. They were discovered by mathematician Michael Somos. From the form of their defining recurrence (which involves division), one would expect the terms of the sequence to be fractions, but surprisingly, a few Somos sequences have the property that all of their members are integers.

Recurrence equations

For an integer number ${\displaystyle k}$ larger than ${\displaystyle 1}$, a Somos-${\displaystyle k}$ sequence ${\displaystyle \{a_{n}\}_{n\in \mathbb {Z} }}$ is a solution of the equation $${\displaystyle a_{n}a_{n-k}=\sum _{i=1}^{\lfloor k/2\rfloor }\alpha _{i}a_{n-i}a_{n-k+i}}$$ where ${\displaystyle \alpha _{1}}$, ${\displaystyle \alpha _{2}}$, ..., ${\displaystyle \alpha _{\lfloor k/2\rfloor }}$ are fixed parameters. It can be rearranged into the form of an order ${\displaystyle k}$ recurrence relation $${\displaystyle a_{n}={\frac {\sum _{i=1}^{\lfloor k/2\rfloor }\alpha _{i}a_{n-i}a_{n-k+i}}{a_{n-k}}}.}$$ Hence a non-degenerate solution is determined by a choice of ${\displaystyle k}$ initial values ${\displaystyle a_{0}}$, ${\displaystyle a_{1}}$, ..., ${\displaystyle a_{k-1}}$. The sequence yielded by setting ${\displaystyle \alpha _{1}=\alpha _{2}=\dots =\alpha _{\lfloor k/2\rfloor }=1}$ and ${\displaystyle a_{0}=a_{1}=\dots =a_{k-1}=1}$ is referred to as the Somos-${\displaystyle k}$ sequence. The Somos-${\displaystyle k}$ sequence is symmetric, i.e., ${\displaystyle s_{-n}=s_{n+k-1}}$.

For ${\displaystyle k=2}$ or ${\displaystyle 3}$, the defining relations are very simple (there is no addition on the right-hand side).

In the first nontrivial case, ${\displaystyle k=4}$, the relation is $${\displaystyle a_{n}={\frac {\alpha a_{n-1}a_{n-3}+\beta a_{n-2}^{2}}{a_{n-4}}}}$$ where ${\displaystyle \alpha }$ and ${\displaystyle \beta }$ are constant parameters. In the case ${\displaystyle k=5}$ the relation is $${\displaystyle a_{n}={\frac {\alpha a_{n-1}a_{n-4}+\beta a_{n-2}a_{n-3}}{a_{n-5}}}.}$$

Sequence values

The Somos-${\displaystyle 2}$ and Somos-${\displaystyle 3}$ sequences are the all-ones sequences (..., ${\displaystyle 1}$, ${\displaystyle 1}$, ${\displaystyle 1}$, ${\displaystyle 1}$, ${\displaystyle 1}$, ...).

The values in the Somos-${\displaystyle 4}$ sequence are

1, 1, 1, 1, 2, 3, 7, 23, 59, 314, 1529, 8209, 83313, 620297, 7869898, ... (sequence A006720 in the OEIS).

The values in the Somos-${\displaystyle 5}$ sequence are

1, 1, 1, 1, 1, 2, 3, 5, 11, 37, 83, 274, 1217, 6161, 22833, 165713, ... (sequence A006721 in the OEIS).

The values in the Somos-${\displaystyle 6}$ sequence are

1, 1, 1, 1, 1, 1, 3, 5, 9, 23, 75, 421, 1103, 5047, 41783, 281527, ... (sequence A006722 in the OEIS).

The values in the Somos-${\displaystyle 7}$ sequence are

1, 1, 1, 1, 1, 1, 1, 3, 5, 9, 17, 41, 137, 769, 1925, 7203, 34081, ... (sequence A006723 in the OEIS).

The first 17 values in the Somos-${\displaystyle 8}$ sequence are

1, 1, 1, 1, 1, 1, 1, 1, 4, 7, 13, 25, 61, 187, 775, 5827, 14815 [the next value is fractional].

Integrality and the Laurent phenomenon

The form of the recurrences describing the Somos sequences involves divisions, making it appear likely that the sequences defined by these recurrence will contain fractional values. Nevertheless, for ${\displaystyle k\leq 7}$ the Somos sequences contain only integer values. Generally, for ${\displaystyle k\leq 7}$, a Somos-${\displaystyle k}$ sequence satisfies the so-called Laurent property. That is, as a function of the initial terms ${\displaystyle a_{0}}$, ..., ${\displaystyle a_{k-1}}$, every term ${\displaystyle a_{n}}$ is a multivariate Laurent polynomial with coefficients in ${\displaystyle \mathbb {Z} [\alpha _{1},\dots ,\alpha _{\lfloor k/2\rfloor }]}$. Several mathematicians have studied the problem of proving and explaining this property of the Somos sequences; it is closely related to the combinatorics of cluster algebras.

For ${\displaystyle k\geq 8}$ the Somos-${\displaystyle k}$ sequences eventually contain fractional values. For Somos-${\displaystyle 8}$ the first fractional value is the 18th term with value ${\displaystyle 420514/7}$.

Closed formula

Consider a Somos-${\displaystyle 4}$ relation $${\displaystyle a_{n}a_{n-4}=\alpha a_{n-1}a_{n-3}+\beta a_{n-2}^{2}}$$ with ${\displaystyle \alpha \neq 0}$. Then a complex-valued Somos-${\displaystyle 4}$ sequence corresponds to an arithmetic sequence of points ${\displaystyle Q+nP}$ on an elliptic curve ${\displaystyle E}$ (see Elliptic curve#Group law). The general term is given by the formula $${\displaystyle a_{n}=AB^{n}{\frac {\sigma (z_{0}+n\kappa )}{\sigma (\kappa )^{n^{2}}}}}$$ where ${\displaystyle \sigma (z)=\sigma (z;g_{2},g_{3})}$ denotes the Weierstrass sigma function associated with the curve ${\displaystyle E}$ written in the canonical form $${\displaystyle E:\quad y^{2}=4x^{3}-g_{2}x-g_{3}.}$$ The six parameters ${\displaystyle g_{2},g_{3},A,B,z_{0},\kappa \in \mathbb {C} }$ are determined uniquely (up to a choice of sign) by the coefficients ${\displaystyle \alpha ,\beta }$ and the initial terms ${\displaystyle a_{0},a_{1},a_{2},a_{3}}$. The ${\displaystyle (x,y)}$ coordinates of the sequence ${\displaystyle Q+nP}$ are given by ${\displaystyle (\wp (z_{0}+n\kappa ),\wp '(z_{0}+n\kappa ))}$, where ${\displaystyle \wp (z)=\wp (z;g_{2},g_{3})}$ denotes the Weierstrass elliptic function. The coefficients ${\displaystyle \alpha }$ and ${\displaystyle \beta }$ are given as elliptic functions of ${\displaystyle \kappa }$ by $${\displaystyle \alpha =\wp (\kappa )^{2},\quad \beta =\wp '(\kappa )^{2}(\wp (2\kappa )-\wp (\kappa )).}$$

See also

• Göbel's sequence

References

External links

• Jim Propp's Somos Sequence Site
• Weisstein, Eric W., "Somos Sequence", MathWorld
• The Troublemaker Number. Numberphile video on the Somos sequences