---
title: 'Shifted Pell Equation: Affine and ζ₃ Variants'
url: https://www.emergentmind.com/topics/shifted-pell-equation
type: topic
---

# Shifted Pell Equation: Affine and ζ₃ Variants

Searching arXiv for recent and foundational papers on shifted Pell equations and related norm-equation generalizations.
Use the arXiv search tool with queries:
1. "shifted Pell equation"
2. "\"zeta_3\" Pell equation"
3. "negative Pell Stevenhagen conjecture"

The shifted Pell equation denotes, in recent arXiv usage, two distinct generalizations of the classical Pell equation. In one sense it is the affine translate
\[
(x-a)^2-D(y-b)^2=1,
\]
with fixed integers \(a,b\) and fixed non-square \(D>1\); in another it is a twisted norm equation over \(K=\mathbb Q(\zeta_3)\),
\[
N_{K(\sqrt[3]{\alpha})/K}(x)=\zeta_3,
\]
which functions as a cubic analogue of the negative Pell equation. In both settings the central structure is unit-theoretic: integer or algebraic-integer solutions are controlled by units in a quadratic or cyclic cubic extension, and the resulting theory combines explicit parametrization with class-group, cohomological, and analytic techniques [2509.17882] [1910.14097].

## 1. Classical norm equations and the meaning of “shift”

The classical Pell equation is
\[
x^2-Dy^2=1,
\]
equivalently the norm equation \(N_{\mathbb Q(\sqrt D)/\mathbb Q}(x)=+1\), while the negative Pell equation is
\[
x^2-Dy^2=-1,
\]
equivalently \(N_{\mathbb Q(\sqrt D)/\mathbb Q}(x)=-1\). The affine shifted Pell equation replaces the quadratic form by a translate,
\[
(x-a)^2-D(y-b)^2=1,
\]
whereas the \(\zeta_3\)-Pell equation keeps the norm-form paradigm but shifts the target norm from \(1\) or \(-1\) to the nontrivial cube root of unity \(\zeta_3\) [2509.17882] [1910.14097].

These two usages encode different kinds of deformation. The affine shift moves the Pell conic inside \(\mathbb Z^2\) without changing the underlying real quadratic unit group. The \(\zeta_3\)-shift changes the codomain of the norm map and places the problem in the Kummer-theoretic setting of cyclic cubic extensions of \(K=\mathbb Q(\zeta_3)\). A common misconception is that “shifted Pell equation” has a unique canonical meaning; the cited literature shows instead that the term is currently used for both an affine translation and a target-norm twist.

## 2. Affine shifted Pell equations over \(\mathbb Z\)

For fixed integers \(a,b\) and fixed non-square \(D>1\), the solution set is
\[
F_D^{(a,b)}=\{(x,y)\in \mathbb Z^2:(x-a)^2-D(y-b)^2=1\}.
\]
Let \((\alpha,\beta)\) be the fundamental solution of
\[
u^2-Dv^2=1,\qquad \alpha+\beta\sqrt D>1,
\]
and define sequences \(\{u_n\}\), \(\{v_n\}\) by
\[
u_n+v_n\sqrt D=(\alpha+\beta\sqrt D)^n,\qquad n=0,1,2,\dots.
\]
Then all integer solutions of the shifted equation arise from the unit group of \(\mathbb Z[\sqrt D]\), concretely through
\[
x-a=\pm u_n,\qquad y-b=\pm v_n.
\]
This yields the complete parametrization of \(F_D^{(a,b)}\) by the powers of the fundamental unit [2509.17882].

The same paper records two equivalent recurrence descriptions. First,
\[
u_0=1,\;v_0=0,\qquad
u_{n+1}=\alpha u_n+D\beta v_n,\qquad
v_{n+1}=\beta u_n+\alpha v_n.
\]
Second,
\[
u_0=1,\;u_1=\alpha,\qquad v_0=0,\;v_1=\beta,
\]
and for \(n\ge 0\),
\[
u_{n+2}=2\alpha u_{n+1}-u_n,\qquad
v_{n+2}=2\alpha v_{n+1}-v_n.
\]
Thus the affine shifted Pell equation is not solved by a new Diophantine mechanism; rather, its integer points are inherited from the classical Pell recurrence and then translated by \((a,b)\).

A useful auxiliary quantity is the “norm-sum” function
\[
f(n)=u_n+v_n,\qquad n\ge 0,
\]
with strictly increasing inverse \(g\) on \([1,\infty)\). The asymptotic relation
\[
f(n)\sim \frac{(\alpha+\beta\sqrt D)^n}{2\sqrt D}(1+\sqrt D)
\]
implies
\[
g(x)=\frac{\log x + C}{\log(\alpha+\beta\sqrt D)}+\mu(x),\qquad
C=\log\!\Bigl(\frac{2\sqrt D}{1+\sqrt D}\Bigr),
\]
where \(\mu(x)\) is a bounded “saw-tooth” function vanishing exactly at the integer points \(x=f(n)\). This gives a quantitative inversion of the exponential growth of Pell solutions.

## 3. Exact enumeration in bounded diamonds

The bounded-region problem in [2509.17882] asks for the points of \(F_D^{(a,b)}\) inside the diamond
\[
|x|+|y|\le \lambda.
\]
Writing
\[
Q_{F_D^{(a,b)}}(\lambda)=\{(x,y)\in F_D^{(a,b)}:|x|+|y|\le \lambda\},
\]
and setting
\[
K=|a|+\max\{|b|+1,\sqrt{1+Db^2}\},
\]
the main enumeration theorem states that if \(\lambda\ge K\), then
\[
Q_{F_D^{(a,b)}}(\lambda)
=
\bigcup_{i,j\in\{\pm1\}}
\{(i\,u_n+a,\;j\,v_n+b):0\le n\le g(\lambda-i\,a-j\,b)\}.
\]
In particular,
\[
\bigl|Q_{F_D^{(a,b)}}(\lambda)\bigr|
=
2+\sum_{i,j=\pm1}\Bigl\lfloor g(\lambda-i\,a-j\,b)\Bigr\rfloor.
\]

The proof strategy is explicit. One first uses the unit-group parametrization to show that every solution has the form \((x-a,y-b)=(\pm u_n,\pm v_n)\). One then rewrites the diamond condition as
\[
|\,i\,u_n+a\,|+|\,j\,v_n+b\,|\le \lambda,
\]
which, once the signs stabilize by quadrant for large \(n\), becomes the inequality
\[
u_n+v_n\le \lambda-i\,a-j\,b.
\]
Since \(f(n)=u_n+v_n\) is strictly increasing, the admissible indices are exactly those with \(n\le g(\lambda-i\,a-j\,b)\). The four sign choices therefore produce four disjoint monotone families of solutions, and for sufficiently large \(\lambda\) no extra points occur.

The resulting algorithm is equally explicit: compute the fundamental solution \((\alpha,\beta)\), for example via the continued fraction of \(\sqrt D\); generate \((u_n,v_n)\) recursively; set \(\Lambda_{i,j}=\lambda-i\,a-j\,b\); compute \(N_{i,j}=\lfloor g(\Lambda_{i,j})\rfloor\); and list all points
\[
(x_n^{\,i,j},y_n^{\,i,j})=(i\,u_n+a,\;j\,v_n+b),\qquad 0\le n\le N_{i,j}.
\]
For the classical case \(D=2\), \((a,b)=(0,0)\), and \(\lambda=20\), the paper gives \(g(20)=2\) and the ten solutions
\[
(\pm1,0),\;(\pm3,\pm2),\;(\pm17,\pm12).
\]
For the shifted example \((x-1)^2-2(y-0)^2=1\) with \(\lambda=10\), the same construction generates candidate families and the diamond constraint filters them to the predicted subset.

The significance of this result lies in its exactness for fixed \(D\). The abstract emphasizes that much of the Pell literature varies \(D\) and studies asymptotics of the fundamental unit, whereas the bounded-region distribution for a fixed \(D\) had received comparatively little attention.

## 4. The \(\zeta_3\)-Pell equation as a shifted norm problem

Let \(\zeta_3=e^{2\pi i/3}\) and \(K=\mathbb Q(\zeta_3)\), with ring of integers \(\mathcal O_K=\mathbb Z[\zeta_3]\). By Kummer theory, every cyclic cubic extension of \(K\) has the form
\[
L=K(\sqrt[3]{\alpha}),\qquad \alpha\in K^\times/(K^\times)^3.
\]
The \(\zeta_3\)-Pell equation is the twisted norm equation
\[
N_{L/K}(x)=\zeta_3,\qquad x\in \mathcal O_L,
\]
equivalently the question whether there exists a unit \(u\in \mathcal O_L^\times\) with relative norm \(N_{L/K}(u)=\zeta_3\) [1910.14097].

This is the direct cubic analogue of the negative Pell equation. In the quadratic case, the local obstruction to
\[
x^2-Dy^2=-1
\]
is that every odd prime \(p\mid D\) must satisfy \(p\equiv 1\pmod 4\). In the \(\zeta_3\)-setting, the local condition is formulated prime by prime for \(\alpha\). For every finite prime \(\mathfrak p\mid \alpha\), one requires that the local extension \(K_{\mathfrak p}(\sqrt[3]{\alpha})/K_{\mathfrak p}\) admit a local norm surjection onto all of \(\mathcal O_{K_{\mathfrak p}}^\times\). Concretely, if \(\mathfrak p\nmid (1-\zeta_3)\), one needs
\[
N_{K/\mathbb Q}(\mathfrak p)\equiv 1\pmod 9,
\]
whereas if \(\mathfrak p\mid (1-\zeta_3)\), the prime above \(3\), there is no local obstruction.

Hence a necessary condition is that in the factorization
\[
\alpha=\prod_i \mathfrak p_i^{e_i},
\]
every prime \(\mathfrak p_i\nmid (1-\zeta_3)\) satisfy
\[
N_{K/\mathbb Q}(\mathfrak p_i)\equiv 1\pmod 9.
\]
The role of this condition parallels the congruence restrictions \(p\equiv 1\pmod 4\) in the negative Pell problem, but the cubic situation is more intricate because local compatibility is not the whole story.

## 5. Cohomology, genus theory, and the Rédei matrix criterion

Knight–Xiao formulate a precise algebraic solvability criterion. Suppose
\[
\alpha=\prod_{i=1}^n T_i,
\]
with each \(T_i\) a prime of \(K\) of norm congruent to \(1\pmod 9\), unramified at \((1-\zeta_3)\). Then there exists \(u\in \mathcal O_L^\times\) with \(N_{L/K}(u)=\zeta_3\) if and only if two conditions hold: first, every \(T_i\equiv 1\pmod{1-\zeta_3}\), equivalently \(N_{K/\mathbb Q}(T_i)\equiv 1\pmod 9\); second, the ramified primes \(\mathfrak P_i\mid T_i\) in \(L/K\) generate the \((1-\zeta_3)\)-torsion subgroup
\[
\mathrm{Cl}(L)[\,1-\zeta_3\,]
\]
of the class group. By Galois cohomology and Hilbert Theorem 90, this is equivalent to the vanishing of
\[
H^2\bigl(\mathrm{Gal}(L/K),\mathcal O_L^\times\bigr)
\]
[1910.14097].

The criterion shows that local congruence conditions are necessary but not sufficient. The global obstruction is encoded in class-group torsion, exactly as in genus-theoretic refinements of the negative Pell problem. This is the main structural difference between merely satisfying local cubic reciprocity conditions and actually solving the twisted norm equation.

A second criterion packages the obstruction in a Rédei matrix. If \(\mathfrak P_i\) denotes the prime above \(T_i\) in \(L\), one forms an \(n\times n\) symmetric matrix
\[
\mathsf R_L=(r_{ij})
\]
over \(\mathbb F_3\). For \(i\ne j\), the off-diagonal entries record the cubic residue symbol,
\[
r_{ij}=\log_{\zeta_3}\!\bigl((T_j)^{\tfrac{N(T_i)-1}{3}}\bigr),
\]
and the diagonal entries are chosen so that each row sums to zero in \(\mathbb F_3\). The rank of \(\mathsf R_L\) detects solvability: if \(\mathsf R_L\) has full rank \(n-1\), then the ramified primes generate \(\mathrm{Cl}(L)[1-\zeta_3]\), hence there is a solution \(u\) with \(N(u)=\zeta_3\); if \(\mathrm{rank}(\mathsf R_L)<n-1\), then no such \(u\) exists.

The Rédei matrix criterion is the computationally effective form of the cohomological theorem. It turns the existence of norm-\(\zeta_3\) units into a finite-field linear-algebra problem whose entries are cubic residue symbols.

## 6. Density theorems, heuristics, and analytic methods

To quantify solvability, the paper defines \(S(X)\) as the set of \(\alpha\in \mathcal O_K\) with \(N_{K/\mathbb Q}(\alpha)\le X\) and no prime divisors except those of norm \(1\pmod 9\), and \(S_3(X)\subset S(X)\) as the subset for which the \(\zeta_3\)-Pell equation is solvable. The key constant is
\[
B=\prod_{i=0}^\infty (1-3^{-2i-1})\approx 0.7112\ldots,
\]
which is exactly the probability that a large random symmetric matrix over \(\mathbb F_3\) is nonsingular. Knight–Xiao prove
\[
B\le
\liminf_{X\to\infty}\frac{|S_3(X)|}{|S(X)|}
\le
\limsup_{X\to\infty}\frac{|S_3(X)|}{|S(X)|}
\le 4B
\]
[1910.14097].

The same work discusses a Stevenhagen-type heuristic. By analogy with the negative Pell equation, one might expect the exact density to equal the naive prediction
\[
2^{-2/3}\approx 0.6300\ldots,
\]
but the paper shows, both theoretically and numerically, that this precise value fails in the \(\zeta_3\)-case. Numerical data up to \(N(\alpha)\approx 10^8\) suggest instead that
\[
\frac{|S_3(X)|}{|S(X)|}\approx 0.75-0.80,
\]
which lies well above \(2^{-2/3}\) and within the interval \([B,4B]\).

The proof architecture has three components. The algebraic component uses classical genus theory together with long exact sequences in Galois cohomology for the unit–ideal exact sequences to derive the local-global criterion and to relate solvability to the \((1-\zeta_3)\)-torsion of the class group. The governing-field component shows that when \(\mathsf R_L\) has \(\mathbb F_3\)-kernel of dimension at least \(d\), one can build a governing field \(F\supset K(\sqrt[3]{\alpha})\) of controlled degree whose Frobenius conditions force failure of full rank; a Chebotarev argument then bounds the frequency of such fields. The analytic component proves that the Rédei matrices behave like random symmetric matrices over \(\mathbb F_3\). This reduces to estimating sums of cubic characters
\[
\sum_{N(\mathfrak p_i)\le X}\Bigl(\frac{\mathfrak p_j}{\mathfrak p_i}\Bigr)_3
\]
and establishing sufficient cancellation via a large-sieve or Heath-Brown style mean-value bound for cubic residue symbols. Combined with a refined count of integers in \(\mathcal O_K\) having exactly \(r\) prime factors, this yields limiting proportions for fixed kernel dimension and the constants
\[
B_r=\prod_{i=0}^{r-1}(1-3^{-2i-1}).
\]

Taken together, the two branches of the subject show that Pell-type shifting can act either on coordinates or on the target of the norm map. In the affine case, the principal achievement is an exact finite-region enumeration in terms of powers of a fundamental unit. In the \(\zeta_3\)-case, the emphasis is instead on solvability criteria, random-matrix statistics, and the failure of the naive Stevenhagen analogue.

Source: https://www.emergentmind.com/topics/shifted-pell-equation