Let $n \geq 2$ and $p$ be a prime. Let $K$ be a number field and consider two Galois representations $Ï_1, Ï_2 : \operatorname{Gal}(\overline{K} / K) \to \operatorname{GL}_n(\mathbb{Z}_p)$ having residual image a $p$-group. We explain and implement an algorithm that makes effective a result of Loïc Grenié to decide wether the semisimplifications of $Ï_1$ and $Ï_2$ are isomorphic.
As an application, we show that an irreducible representation $Ï: G_{\mathbb{Q}(\sqrt{-3})} \to \operatorname{GL}_2(\mathbb{Z}_3)$ unramified outside 3 is determined by the characteristic polynomials of Frobenius elements at five primes of small norm. As an additional check, we apply it to a 2-adic example studied by Grenié, recovering Grenié's result in a fully automated way.
Paper ID : 2510.12956Title : Comparing Galois representations in the residually reducible caseAuthors : Nuno Freitas, Ignasi Sánchez-RodríguezClassification : math.NT (Number Theory)Publication Date : October 14, 2025Paper Link : https://arxiv.org/abs/2510.12956v1 Let n ≥ 2 n \geq 2 n ≥ 2 and p p p be a prime, with K K K a number field. Consider two Galois representations ρ 1 , ρ 2 : Gal ( K ‾ / K ) → GL n ( Z p ) \rho_1, \rho_2 : \text{Gal}(\overline{K}/K) \to \text{GL}_n(\mathbb{Z}_p) ρ 1 , ρ 2 : Gal ( K / K ) → GL n ( Z p ) whose residual images are p p p -groups. This paper explains and implements an algorithm that makes Loïc Grenié's results effective for determining whether the semisimplifications of ρ 1 \rho_1 ρ 1 and ρ 2 \rho_2 ρ 2 are isomorphic. As an application, it is proved that an irreducible representation ρ : G Q ( − 3 ) → GL 2 ( Z 3 ) \rho: G_{\mathbb{Q}(\sqrt{-3})} \to \text{GL}_2(\mathbb{Z}_3) ρ : G Q ( − 3 ) → GL 2 ( Z 3 ) unramified outside 3 is completely determined by the characteristic polynomials of Frobenius elements at five primes of small norm. As additional verification, the algorithm is applied to the 2-adic example studied by Grenié, recovering Grenié's results in a completely automated manner.
This research concerns the comparison problem for Galois representations, a central problem in algebraic number theory. According to Faltings' foundational work, two continuous representations ρ 1 , ρ 2 : G K → GL n ( E ) \rho_1, \rho_2 : G_K \to \text{GL}_n(E) ρ 1 , ρ 2 : G K → GL n ( E ) unramified outside a set S S S are isomorphic if and only if they have the same traces at Frobenius elements in a finite set T T T .
Faltings-Serre method : Applies to representations with values in GL 2 ( Q 2 ) \text{GL}_2(\mathbb{Q}_2) GL 2 ( Q 2 ) and absolutely irreducible residual imageLivné's extension : Handles 2-adic 2-dimensional representations with 2-group residual imageGrenié's generalization : Extends to arbitrary p p p -adic n n n -dimensional representations, but still requires p p p -group residual imageThe main limitation of these methods is that even for small n n n , p p p , and S S S , they require constructing large number fields, making computation infeasible.
The core motivations of this paper are:
Algorithm Implementation : Transform Grenié's theoretical results into computable algorithmsGenerality : Design algorithms applicable to general cases without relying on additional information about specific representationsDatabase Potential : Run the algorithm once for each parameter set, potentially creating databasesPractical Applications : Prove modularity of abelian surfacesAlgorithm Implementation : Implements a variant of Grenié's theorem (Theorem 4.1), providing a fully automated algorithmTheoretical Improvement : Proposes a method for constructing smaller extension fields K S K_S K S , significantly reducing computational complexityConcrete Application Results :
Proves specific results in the 3-adic case (Corollary 1.1) Recovers Grenié's 2-adic results (Corollary 1.2) Modularity Proof : First successful application of 3-adic Faltings-Serre type methods to prove modularity in the residually reducible caseOpen Source Implementation : Provides complete Magma code implementationInput : Two continuous Galois representations ρ 1 , ρ 2 : G K → GL n ( E ) \rho_1, \rho_2 : G_K \to \text{GL}_n(E) ρ 1 , ρ 2 : G K → GL n ( E ) unramified outside a set S S S , with p p p -group residual images
Output : Determine whether the semisimplifications of ρ 1 \rho_1 ρ 1 and ρ 2 \rho_2 ρ 2 are isomorphic
Constraints : n ≥ 2 n \geq 2 n ≥ 2 , p p p prime, K K K a number field
The traditional K S K_S K S construction requires λ + ε + m \lambda + \varepsilon + m λ + ε + m layers, where:
λ \lambda λ satisfies 2 λ ≥ r = N 2 ( 1 + ε ) N ( N − 1 ) / 2 2^\lambda \geq r = N^2(1+\varepsilon)^{N(N-1)/2} 2 λ ≥ r = N 2 ( 1 + ε ) N ( N − 1 ) /2 ε = 0 \varepsilon = 0 ε = 0 (if p ≠ 2 p \neq 2 p = 2 ) or ε = 1 \varepsilon = 1 ε = 1 (if p = 2 p = 2 p = 2 )m m m satisfies p m ≥ n p^m \geq n p m ≥ n The improved method proposed in this paper:
K₀ = K
For i ≥ 0:
Kᵢ₊₁ = maximal p-elementary extension of Kᵢ satisfying:
1. Unramified outside S
2. Kᵢ₊₁/K is Galois
3. Exponent of Gal(Kᵢ₊₁/K) divides pᵐ
If Kᵢ₊₁ = Kᵢ for some i < λ + ε + m, stop and set K_S = Kᵢ
Let Σ ⊂ G K \Sigma \subset G_K Σ ⊂ G K be a set, and define Σ ′ = { σ k : σ ∈ Σ , k ∈ Z ≥ 0 } \Sigma' = \{\sigma^k : \sigma \in \Sigma, k \in \mathbb{Z}_{\geq 0}\} Σ ′ = { σ k : σ ∈ Σ , k ∈ Z ≥ 0 } .
If:
Σ ′ \Sigma' Σ ′ covers generators of each maximal cyclic subgroup of Gal ( K S / K ) \text{Gal}(K_S/K) Gal ( K S / K ) For all σ ∈ Σ \sigma \in \Sigma σ ∈ Σ , CharPoly ( ρ 1 ( σ ) ) = CharPoly ( ρ 2 ( σ ) ) \text{CharPoly}(\rho_1(\sigma)) = \text{CharPoly}(\rho_2(\sigma)) CharPoly ( ρ 1 ( σ )) = CharPoly ( ρ 2 ( σ )) Then the semisimplifications of ρ 1 \rho_1 ρ 1 and ρ 2 \rho_2 ρ 2 are isomorphic.
Part A: Construction of first m layers
K₀ = K
For 0 ≤ i ≤ m-1:
Kᵢ₊₁ = Kᵢ(ᵖ√α₁, ..., ᵖ√αₐ)
where αᵢ are generators of Selₚ(Kᵢ, S)
Part B: Subsequent layer construction (i ≥ m)
Compute Gal(Kᵢ/K)-submodules : Calculate all Galois submodules of Selₚ(Kᵢ, S)Eliminate 1-dimensional submodules : Check residue degree condition f ( Q ∣ Q ∩ O K ) ≤ p m f(Q|Q \cap O_K) \leq p^m f ( Q ∣ Q ∩ O K ) ≤ p m Eliminate higher-dimensional submodules : Apply Lemma 3.4 recursivelyConstruct field Kᵢ₊₁ : Select the unique submodule satisfying the conditionsMemory Optimization : Avoid storing all submodules by using inductive computation of relevant submodulesResidue Degree Bound : Limit extension field size by controlling residue degree ≤ p m \leq p^m ≤ p m Early Stopping Criterion : Theorem 8.1 provides natural stopping conditions for the algorithmGRH Avoidance : Under specific conditions, can avoid the Generalized Riemann HypothesisSoftware : Magma V2.28-20Hardware : 4-thread configurationCode : Complete implementation available on GitHubParameters : n = 3 n = 3 n = 3 , p = 2 p = 2 p = 2 , K = Q K = \mathbb{Q} K = Q , S = { 2 } S = \{2\} S = { 2 } Extension field layers : K S = K 3 K_S = K_3 K S = K 3 (degree 64 extension)Runtime : Approximately 30 hours (assuming GRH)Parameters : n = 2 n = 2 n = 2 , p = 3 p = 3 p = 3 , K = Q ( − 3 ) K = \mathbb{Q}(\sqrt{-3}) K = Q ( − 3 ) , S = { q 3 } S = \{q_3\} S = { q 3 } Extension field layers : K S = K 2 K_S = K_2 K S = K 2 (degree 54 extension)Runtime : Approximately 17 hours (assuming GRH)Correctness : Consistency with known theoretical resultsEfficiency : Computation time and memory usagePracticality : Ability to handle real modularity problemsLet K = Q ( − 3 ) K = \mathbb{Q}(\sqrt{-3}) K = Q ( − 3 ) , and ρ 1 , ρ 2 : G K → GL 2 ( Z 3 ) \rho_1, \rho_2 : G_K \to \text{GL}_2(\mathbb{Z}_3) ρ 1 , ρ 2 : G K → GL 2 ( Z 3 ) be continuous representations unramified outside q 3 q_3 q 3 with the same determinant, which is trivial modulo 3. Then ρ 1 \rho_1 ρ 1 and ρ 2 \rho_2 ρ 2 have isomorphic semisimplifications if and only if for all t ∈ T t \in T t ∈ T , where
T = { 2 O K , p 7 1 , p 7 2 , p 19 ∗ , p 73 ∗ } T = \{2O_K, p_7^1, p_7^2, p_{19}^*, p_{73}^*\} T = { 2 O K , p 7 1 , p 7 2 , p 19 ∗ , p 73 ∗ } ρ 1 ( Frob t ) \rho_1(\text{Frob}_t) ρ 1 ( Frob t ) and ρ 2 ( Frob t ) \rho_2(\text{Frob}_t) ρ 2 ( Frob t ) have the same trace.
For n = 3 n = 3 n = 3 or n = 4 n = 4 n = 4 , let ρ 1 , ρ 2 : G Q → GL n ( Z 2 ) \rho_1, \rho_2 : G_{\mathbb{Q}} \to \text{GL}_n(\mathbb{Z}_2) ρ 1 , ρ 2 : G Q → GL n ( Z 2 ) be continuous representations unramified outside 2 with 2-group residual image. The prime sets
T = { 5 , 7 , 11 , 17 , 23 , 31 } , T ′ = { 5 , 7 , 11 , 17 , 19 , 23 , 31 , 73 , 137 , 257 , 337 } T = \{5,7,11,17,23,31\}, \quad T' = \{5,7,11,17,19,23,31,73,137,257,337\} T = { 5 , 7 , 11 , 17 , 23 , 31 } , T ′ = { 5 , 7 , 11 , 17 , 19 , 23 , 31 , 73 , 137 , 257 , 337 }
can be used to determine isomorphism.
The algorithm verifies the modularity of the Jacobian of the genus 2 curve with LMFDB label 2187.a.6561.1.
Verifies the modularity of the Jacobian of the genus 2 curve with label 59049.a.177147.1.
Memory Usage : Reduced from 52GB to 600MB (through optimized submodule computation)Extension Field Size : Significantly smaller than theoretical boundsPractical Feasibility : First successful automated application of the 3-adic Faltings-Serre methodFaltings (1983) : Establishes foundational theoretical frameworkSerre : Proposes computational methods for the GL₂(ℚ₂) caseLivné (1987) : Extends to 2-adic 2-dimensional residually reducible caseGrenié (2007) : Generalizes to arbitrary p p p -adic n n n -dimensional caseDuan (2021) : Variant methods exploiting special structureThis paper is the first complete automated implementation of Grenié's general theory and successfully applies it to real modularity problems.
Theoretical Contribution : Successfully implements an automated algorithm for comparing Galois representations in the residually reducible caseComputational Breakthrough : Significantly reduces computational complexity through improved extension field constructionPractical Application : First successful application of 3-adic methods to prove modularity of abelian surfacesGRH Dependence : Most computations still require assuming the Generalized Riemann HypothesisParameter Restrictions : Only applies to cases where residual image is a p p p -groupComputational Complexity : Still computationally intensive for large parametersGRH Avoidance : Find more cases that do not depend on GRHAlgorithm Optimization : Further reduce computational complexityApplication Extension : Apply to more modularity problemsDatabase Construction : Build precomputed databases for commonly used parametersTheory-Practice Integration : Transforms deep theoretical results into practical algorithmsTechnical Innovation : Important innovations in memory optimization and algorithm designPractical Value : Solves concrete modularity problems, demonstrating method practicalityOpen Source Contribution : Provides complete code implementation, promoting research reproducibilityComputational Limitations : Still infeasible for larger parametersTheoretical Assumptions : Heavily relies on unproven assumptions like GRHLimited Scope : Restricted to specific types of Galois representationsAcademic Value : Provides new tools for computational problems in algebraic number theoryPractical Value : Can be used to verify modularity of elliptic curves and abelian varietiesMethodological Contribution : Demonstrates how to algorithmize theoretical resultsModularity Verification : Modularity proofs for elliptic curves and abelian varietiesGalois Representation Comparison : Determining representation isomorphism for specific parametersTheoretical Verification : Computational verification of known resultsThe paper cites 17 important references, including Faltings' foundational work, Grenié's theoretical contributions, and related computational method developments. Particularly important are Grenié (2007) for the original theory and van Geemen-Top (1994) for providing test examples.