66 votos

Proyecto abierto: Calculemos la expansión de Fourier de una forma algebraica de Maass no resoluble.

Bien, veamos si puedo usar MO para calcular explícitamente un ejemplo de algo, haciendo que otras personas se unan. Una especie de "un nivel superior" - a menudo la gente responde a las preguntas aquí, pero voy a ver si puedo hacer que la gente haga un proyecto más sustancial. Antes de empezar, tened en cuenta

(1) el cómputo podría haberse realizado ya [me encantaría conocer una referencia, si es así]

(2) el cómputo puede o no ser digno de publicarse

(3) Si vale la pena publicarlo, puede haber o no un debate sobre quiénes son los autores.

A mí, personalmente, no me importan (2) ni (3) en este momento, pero a otros sí. Permítanme pasar a las matemáticas. Oh---sólo un par de cosas más antes de empezar---este proyecto está relacionado con las matemáticas en esta pregunta pero tal vez lo lleve un poco más lejos (si conseguimos que funcione). Inicialmente había pensado en estos temas porque iba a dárselos a un estudiante de licenciatura, pero el estudiante me dice hoy que ha decidido hacer su proyecto sobre el caso holomorfo, y me pareció un poco tonto dejar que mi inversión inicial en el problema se desperdiciara, así que pensé en decírselo a quien estuviera interesado. Si nadie pica el anzuelo aquí, probablemente haré otro proyecto UG.

Bien, este es el trato. Digamos que $K/\mathbf{Q}$ es una extensión de Galois finita, y $\rho:Gal(K/\mathbf{Q})\to GL(2,\mathbf{C})$ es una representación bidimensional irreducible. Las conjeturas generales de la filosofía de Langlands predicen que $\rho$ proviene de una forma automórfica en $GL(2)$ en $\mathbf{Q}$ . La idea es que vamos a "ver" esta forma en un ejemplo explícito en el que la teoría general aún no demuestra que exista.

Ahora el determinante de $\rho$ es una representación de Galois unidimensional, y tiene sentido preguntarse si $det(\rho(c))$ es $+1$ ou $-1$ , donde $c$ es la conjugación compleja. Tiene que ser una de estas, porque $c^2=1$ . La naturaleza de la forma automórfica que se predice depende del signo.

Si el determinante es $-1$ entonces la forma debe ser holomorfa, y una forma de cúspide clásica de peso 1. En este caso se conoce la existencia de la forma, porque está implícita en la conjetura de Serre, que ahora es un teorema de Khare y Wintenberger.

Si el determinante es $+1$ entonces la forma conjetural debe ser una función real-analítica en el semiplano superior, invariante bajo un subgrupo de congruencia, y que satisface una cierta ecuación diferencial (que no son las ecuaciones de Cauchy-Riemann en este caso). Si la imagen de $\rho$ es un grupo soluble, entonces la existencia de esta forma es conocida por antiguos trabajos de Langlands y Tunnell.

Así que, en resumen, el único caso en el que no se conoce la forma es cuando el determinante de la conjugación compleja es $+1$ y la imagen de Galois no es resoluble.

He aquí un ejemplo explícito. El polinomio

g5=344 + 3106*x - 1795*x^2 - 780*x^3 - x^4 + x^5

tiene campo de división $L$ , un $A_5$ -extensión de $\mathbf{Q}$ ramificado sólo en el primer 1951. Ahora $A_5$ es isomorfo al cociente de $SL(2,\mathbf{F}_5)$ por su centro $\pm 1$ y $L$ tiene una extensión de grado dos $K_0$ también sin ramificar fuera de 1951, con $Gal(K_0/\mathbf{Q})$ ser $SL(2,\mathbf{F}_5)$ . Resulta que $K_0$ puede tomarse como el campo de división del polinomio más complicado

g24 = 14488688572801 - 2922378139308818*x^2 + 134981448876235615*x^4 -   1381768039105642956*x^6 + 4291028045077743465*x^8 -   2050038614413776542*x^10 + 287094814384960835*x^12 - 9040633522810414*x^14 + 63787035668165*x^16 - 158664037068*x^18 +   152929135*x^20 - 50726*x^22 + x^24

Tal vez debería decir que David Roberts me contó estos polinomios en enero de 2008; están en un artículo suyo y de John Jones, pero los aprendieron de un artículo de Doud y Moore.

Ahora $SL(2,\mathbf{F}_5)$ tiene dos representaciones complejas bidimensionales fieles; las trazas de cada representación toman valores en $\mathbf{Q}(\sqrt{5})$ y uno es, por supuesto, el conjugado del otro. El determinante de ambas representaciones es trivial así que de hecho son $SL(2,\mathbf{C})$ -valorado. Oh también, todas las raíces de g24 son reales y por lo tanto $K_0$ es totalmente real. Así que lo que tenemos aquí es una representación $$\rho_0:Gal(K_0/\mathbf{Q})\to GL(2,\mathbf{C})$$ que se conjetura que proviene de las formas automórficas, pero, hasta donde yo sé, la conjetura no se conoce en este caso.

Por desgracia, el conductor de $\rho_0$ es $1951^2$ que es un poco grande. De hecho, permítanme decir algo más sobre lo que ocurre en 1951. En el $A_5$ los grupos de descomposición e inercia en 1951 son ambos cíclicos de orden 5. Si mi comprensión de lo que me dijo David Roberts es correcta, en el $SL(2,\mathbf{F}_5)$ extensión los grupos de descomposición e inercia son ambos cíclicos de orden 10 (de hecho acabo de conseguir que magma lo compruebe). Pero el resultado es que $\rho_0$ restringido a un grupo de descomposición en 1951 es de la forma $\psi+\psi^{-1}$ avec $\psi$ de orden 10. ¿Qué carácter de orden 10? Pues hay cuatro caracteres $(\mathbf{Z}/1951\mathbf{Z})^\times\to\mathbf{C}^\times$ de orden 10, y dos de ellos servirán, y dos no, y cuáles servirán depende de qué representación bidimensional de $SL(2,\mathbf{F}_5)$ que usted eligió.

Sin embargo, el punto clave es que si consigues $\psi$ derecho, entonces el giro $\rho:=\rho_0\otimes\psi$ tendrá el conductor 1951, que es pequeño para estos fines.

Ahora, como Marty hizo en un $A_4$ ejemplo y como hice y Junkie hizo en un ejemplo de diedro en la pregunta sobre la forma de Maass citada anteriormente, es posible averiguar explícitamente los números $b_1$ , $b_2$ , $b_3$ ..., con la propiedad de que

$$L(\rho,s)=\sum_{n\geq1}b_n/n^s.$$

Si uno tuviera un programa de ordenador que pudiera calcular $b_n$ para $n\geq1$ Entonces no hay una, sino dos formas en las que se podría intentar dar evidencia computacional a las predicciones dadas por la filosofía de Langlands:

(A) se podrían utilizar las técnicas que Fernando Rodríguez-Villegas me explicó hace unos meses para intentar conseguir pruebas computacionales que $L(\rho,s)$ tenía continuación analítica en el plano complejo y satisfacía la ecuación funcional correcta, y

(B) se podría calcular la correspondiente función analítica real en el semiplano superior, evaluarla en varios lugares con 30 decimales, y ver si la función es invariante bajo el grupo $\Gamma_1(1951)$ .

No sé mucho sobre (A) pero una vez intenté, y fracasé, hacer (B), y mi intuición es que mi error está en el programa de ordenador que escribí para calcular el $b_n$ . Pero como indica la respuesta de junkie en la pregunta anterior, parece que ahora hay varias formas de calcular el $b_n$ y una cosa que me pregunto es si podemos utilizar los métodos que él/ella indicó en esta pregunta.

Permítanme hablar más sobre cómo traté de calcular el $b_n$ . El carácter de la forma Maass es $\psi^2$ el determinante de $\rho$ . La teoría general nos dice que $b_n$ es una función multiplicativa de $n$ por lo que sólo necesitamos calcular $b_n$ para $n$ una potencia principal. De nuevo, la teoría general (considere la $L$ -) dice que para $p\not=1951$ se puede calcular $b_{p^n}$ de $b_p$ . Si $p=1951$ entonces $b_{p^n}=1$ para todos $n$ porque los grupos de descomposición e inercia coinciden para 1951 en $K_0$ [EDIT: Esta parte del argumento es errónea, y explica por qué mis programas no funcionaban. Finalmente descubrí mi error tras comparar la salida de mis programas y los de Junkie y ver en qué se diferenciaban. Es cierto que la descomposición y la inercia coinciden en $K_0$ pero cuando se tuerce por el carácter de orden 10 esto deja de ser cierto. De hecho $b_{1951}$ es una raíz 5ª primitiva de la unidad que no sé cómo calcular con mi método más que por ensayo y error]. Por último, si el $L$ -función de $\rho_0$ es $\sum_n a_n/n^s$ entonces $b_n=\psi(n)a_n$ para todos $n$ primo a 1951, por lo que basta con calcular $a_p$ para $p\not=1951$ .

Para calcular $a_p$ Voy a calcular la traza de $\rho_0(Frob_p)$ . Primero calculo el GCD de los grados de los factores irreducibles de g24 mod $p$ . Si $p$ no divide el discriminante de g24 entonces este GCD es el orden de $Frob_p$ en $SL(2,\mathbf{F}_5)$ . Si $p$ divide el discriminante de g24 entonces mala suerte, necesito factorizar $p$ en el anillo de enteros del campo numérico generado por una raíz de g24. He hecho esto usando magma. Aquí están los resultados:

                                                   prime   order
                                                     2       6
                                                     3      10
                                                     5       4
                                                   163       5
                                                 16061       1
                                                889289      10
                                          451400586583       2
                             1188493301983785760551727       2
120450513180827412314298160097013390669723824832697847       1

Ahora bien, desgraciadamente el simple cálculo del orden de la clase conj de Frobenius no es suficiente para determinar la traza de la representación de Galois, porque $SL(2,\mathbf{F}_5)$ contiene dos clases conjuntas de elementos de orden 5, y dos de orden 10. Sin embargo, en ambos casos, las clases conj se mantienen distintas en $A_5$ por lo que basta con tener un algoritmo que pueda distinguir entre las dos clases de conjugación en $A_5$ . Más concretamente, tenemos que resolver el siguiente problema: etiquetamos las clases conjuntas de elementos de orden 5 en $A_5$ como C1 y C2, y queremos un algoritmo que, dado un primo $p$ para el que g5 es irreducible mod $p$ queremos que el algoritmo devuelva "C1" o "C2" dependiendo de la clase $Frob_p$ está en. He aquí una bonita manera de hacerlo, que me explicó Bjorn Poonen: si g5 es irreducible mod $p$ entonces sus raíces en un cierre alg de $\mathbf{F}_p$ son $x,x^p,x^{p^2},...$ . Establecer $x_i=x^{p^i}$ y calcular $\prod_{i<j}(x_i-x_j)$ . Este producto es una raíz cuadrada del discriminante de g5. Elegir de una vez por todas una raíz cuadrada del discriminante de g5 en los enteros; si el producto es congruente con este mod $p$ devuelve "C1", si no devuelve "C2".

¡Eso es! He implementado esto. Tenía un programa que devolvía un montón de $a_n$ y, por lo tanto, un montón de $b_n$ 's. Construí la función en el semiplano superior como la suma habitual que implica funciones de Bessel y demás, pero computacionalmente no resultó ser invariante bajo $\Gamma_1(1951)$ .

Si alguien quiere aceptar el reto de calcular el $b_n$ eso sería genial. He explicado una forma de hacerlo más arriba, pero soy muy consciente de que puede haber otras formas de calcular el $b_n$ análogo al planteamiento del drogadicto de la pregunta anterior.

27voto

varunsrin Puntos 113

Aquí está el código Magma que te da la respuesta en pocos segundos. Hice un caso especial para los primos malos, y los hice a mano.

_<x> := PolynomialRing(Rationals());
f5 := 344 + 3106*x - 1795*x^2 - 780*x^3 - x^4 + x^5;
g24 := 14488688572801 - 2922378139308818*x^2 + 134981448876235615*x^4 -
       1381768039105642956*x^6 + 4291028045077743465*x^8 -
       2050038614413776542*x^10 + 287094814384960835*x^12 -
       9040633522810414*x^14 + 63787035668165*x^16 - 158664037068*x^18 +
       152929135*x^20 - 50726*x^22 + x^24;
K := NumberField(f5);
_,D := IsSquare(Integers()!Discriminant(f5));
prec := 30;
CHAR_TABLE := CharacterTable(GaloisGroup(g24));
chi := CHAR_TABLE[2];

BAD_FACTORS :=
 [ <2,Polynomial([1,-1,1])>,
   <3,Polynomial([1,-ComplexField(prec)!chi[9],1])>,
   <5,Polynomial([1,0,1])>,
   <7,Polynomial([1,0,1])>,
   <71,Polynomial([1,0,1])>,
   <137,Polynomial([1,1,1])>,
   <163,Polynomial([1,-ComplexField(prec)!chi[5],1])>,
   <1951,Polynomial([1])>,
   <16061,Polynomial([1,-2,1])>,
   <889289,Polynomial([1,-ComplexField(prec)!chi[8],1])> ];
BAD := [bf[1] : bf in BAD_FACTORS];
FACTORS := [bf[2] : bf in BAD_FACTORS];

function LOCAL(p,d : Precision:=prec)
  if p in BAD then return FACTORS[Position(BAD,p)]; end if;
  R := Roots(ChangeRing(f5,GF(p)));
  if #R eq 1 then return Polynomial([1,0,1]); end if;
  if #R eq 2 then
    ord := Lcm([Degree(f[1]) : f in Factorization(Polynomial(GF(p),g24))]);
    return Polynomial([1,ord eq 3 select 1 else -1,1]); end if;
  if #R eq 5 then
    ord := Lcm([Degree(f[1]) : f in Factorization(Polynomial(GF(p),g24))]);
    return Polynomial([1,ord eq 1 select -2 else 2,1]); end if;
  r := Roots(ChangeRing(f5,GF(p^5)));
  x := r[1][1];
  prod := GF(p)!&*[x^(p^i)-x^(p^j) : j in [(i+1)..4], i in [0..4]];
  wh := prod eq GF(p)!D;
  ord := Lcm([Degree(f[1]) : f in Factorization(Polynomial(GF(p),g24))]);
  if ord eq 10 then class := wh select 8 else 9; // compatible with FACTORS
  else class := wh select 6 else 5; end if;
  return Polynomial([1,-ComplexField(prec)!chi[class],1]);
 end function;

L := LSeries(1, [0,0], 1951^2, LOCAL : Precision:=prec);
// s->1-s, Gamma(s/2)^2
psi := DirichletGroup(1951, CyclotomicField(10)).1;
p1951 := Polynomial([1,-ComplexField(prec)!CyclotomicField(5).1]);
TP := TensorProduct(L, LSeries(psi : Precision:=prec), [<1951, 1, p1951>]);
CheckFunctionalEquation(TP);

Aquí están los valores especiales:

ev := Evaluate(TP,0); // 2-1.453085056...
rel := PowerRelation(ev,4 : Al:="LLL");
NF := NumberField(rel);
Q5<zeta5> := CyclotomicField(5);
assert IsIsomorphic(NF,Q5);
Q5!NF.1;

Así que $L(\rho,0)=-4\zeta_5(1+\zeta_5)$ para Marty. Me dan $L(\rho_0,-1)=32(48723\sqrt{5} - 778741)$ como un algebraico. Obtengo $L(\rho,-2)=8800\zeta_5^3 - 14444\zeta_5^2 + 35604\zeta_5 + 17412$ con más precisión. Determiné el factor TensorProduct en 1951 mediante prueba y error, haciendo las conjeturas obvias hasta que una funcionó (el fallo está en 100-110 dígitos). Con esto, lo llevo a 240 dígitos e incluso puedo obtener $$L(\rho,-4)=-18475535360\zeta_5^3 - 11142861380\zeta_5^2 - 12091894020\zeta_5 - 7107607296$$ y $$L(\rho,-6)=25255057273186244\zeta_5^3 - 1015274469604000\zeta_5^2 - 15695788409197884\zeta_5 + 9459547822189412$$ La precisión puede ser mayor si quieres más.

Por último, el formulario Maass:

function MaassEval(L,z)
  x:=Real(z); y:=Imaginary(z);                                                   
  printf "Using %o coefficients\n", Ceiling(11/y);                               
  C := LGetCoefficients(L,Ceiling(11/y));                                        
  pi := Pi(RealField());                                                         
  a := Sqrt(y)*&+[C[n]*KBessel(0,2*pi*n*y)*Sin(2*pi*n*x) : n in [1..#C]];        
  return a;                                                                      
end function;                                                                   
zz:=0.0001+0.0001*ComplexField().1;
MaassEval(TP,zz);   
// Using 110000 coefficients
// -1.71477211817772949974178783985E-8 + 9.01673609747756708674470686948E-9*i
MaassEval(TP,zz/(1951*zz+1));
// Using 161297 coefficients
// -1.71477211817772949974179078240E-8 + 9.01673609747756708674496293450E-9*i

8voto

Bob Somers Puntos 4186

He aceptado la respuesta de Junkie a esta pregunta. Pero al analizar su programa que funcionaba y compararlo con mi programa que no funcionaba pude encontrar fácilmente el error en mi programa: Pensé erróneamente que el coeficiente de $q^{1951}$ era 1 y no lo es, es una raíz 5ª primitiva de la unidad. Mi código funciona ahora, y está escrito en PARI-GP, así que pensé en publicarlo para que la gente que no tiene acceso a magma pueda jugar con el milagro que es una eigenforma de Maass cuspidal que (conjeturalmente) corresponde a una representación de Galois par no resoluble. Hago mi post comunitario wiki porque mi post no dice nada que el de Junkie no diga, así que no quiero ningún crédito por ello.

\\ nonsolvmaass.g: works out a maass form associated to an explicit
\\ nonsolvable finite image even irred 2-d Galois rep.
\\ Note: rho(c)=-1 so we use sines.

\\ N is number of coefficients of the L-function we will compute.

N=2000000;

\\ g5 has splitting field an A_5 field.

g5=344 + 3106*x - 1795*x^2 - 780*x^3 - x^4 + x^5

\\ g24 has splitting field a field with Galois group SL(2,5).

g24 = 14488688572801 - 2922378139308818*x^2 + 134981448876235615*x^4 -   1381768039105642956*x^6 + 4291028045077743465*x^8 -   2050038614413776542*x^10 + 287094814384960835*x^12 -   9040633522810414*x^14 + 63787035668165*x^16 - 158664037068*x^18 +   152929135*x^20 - 50726*x^22 + x^24

\\ The character table of SL(2,5) starts like this:
\\
\\
\\ -------------------------------------------
\\ Class |   1  2  3  4    5    6  7    8    9
\\ Size  |   1  1 20 30   12   12 20   12   12
\\ Order |   1  2  3  4    5    5  6   10   10
\\ -------------------------------------------
\\ X.1   +   1  1  1  1    1    1  1    1    1
\\ X.2   -   2 -2 -1  0   Z1 Z1#2  1  -Z1-Z1#2
\\ X.3   -   2 -2 -1  0 Z1#2   Z1  1-Z1#2  -Z1
\\
\\...
\\
\\ # denotes algebraic conjugation, that is,
\\ #k indicates replacing the root of unity w by w^k
\\
\\ Z1     = (CyclotomicField(5)) ! [ RationalField() | 0, 0, 1, 1 ]
\\
\\ If we let K be the splitting field of g24 then K contains a splitting
\\ field for g5. If we fix an isomorphism of Gal(K/Q) with SL(2,5) and
\\ use X.2 to define a 2-dimensional representation of Gal(K/Q), then
\\ the resulting representation has conductor 1951^2. 
\\
\\ Here is a pari script to compute the order of Frob_p for the resulting
\\ representation (for p not 1951). First a function that only works for the primes
\\ not dividing disc(g24).

ordfrob0(p)=local(m);m=factormod(g24,p);lcm(vector(#m~,i,poldegree(m[i,1])))

\\ For the primes dividing disc(g24) other than 1951 I used magma to work
\\ out the answer by hand [computing a maximal order etc].

ordfrob(p)=if(p==1951,0,if(p==2,6,if(p==3,10,if(p==5,4,if(p==163,5,if(p==16061,1,if(p==889289,10,if(p==451400586583,2,if(p==1188493301983785760551727,2,if(p==120450513180827412314298160097013390669723824832697847,1,ordfrob0(p)))))))))))

\\ If the order of Frob_p in SL(2,5) is 5 or 10 then we have to decide which
\\ of the two conjugacy classes p is in! Classes 5 and 8 reduce to the same
\\ class in A_5, and classes 6 and 9 reduce to the other. So we have to
\\ be able to distinguish between the classes in A_5. Bjorn Poonen told
\\ me how to do this. The trick is to compute the product whose square
\\ is the discriminant mod p and to see which of the square roots of the
\\ discriminant that the product itself is isomorphic to.
\\ This will not work for any p dividing the discriminant of g5.
\\ Fortunately I checked explicitly that no prime p dividing disc(g5)
\\ has Frob_p of order 5 or 10. Hooray for good fortune!
\\ The function returns 0 or 1.

whichclass(p)=local(t,conjs,pro);t=Mod(x,Pol(Mod([1, -1, -780, -1795, 3106, 344],p)));conjs=[t];for(i=1,4,conjs=concat(conjs,conjs[i]^p));pro=1;for(i=1,5,for(j=i+1,5,pro=pro*(conjs[i]-conjs[j])));pro==518348075378;

\\ Now finally an algorithm to compute a_p.

zet10=exp(2*Pi*I/10);Z1=zet10^4+zet10^6;Z1h2=zet10^2+zet10^8;

ap(p)=if(p==1951,0,orp=ordfrob(p);if(orp==1,2,if(orp==2,-2,if(orp==3,-1,if(orp==4,0,if(orp==6,1,whp=whichclass(p);if(orp==5,if(whp==0,Z1,Z1h2),if(whp==0,-Z1,-Z1h2))))))))

\\ The problem with the Maass form associated to these eigenvalues is that, although it
\\ has trivial character, it has conductor 1951^2. At a prime above 1951, the
\\ inertia subgroup has order 10, as does the decomposition group. The representation
\\ sends a generator of this inertia subgroup to something with eigenvalues zet10
\\ and zet10^9. So we need to twist by a Dirichlet character of conductor 1951
\\ and order 10 and we will get a representation of conductor 1951.
\\ But which one?! Well, there are only four, and two will work, so we have
\\ a sporting chance! If we get it wrong we can just try another one.
\\ Note that 3 is a primitive root mod 1951.

twist=vector(1950,i,0);z0=zet10^3;z=z0;n=Mod(3,1951);for(i=1,1950,twist[lift(n)]=z;n=n*3;z=z*z0);

\\ The character of the form we are after is the square of "twist".

\\ Here's the trace of Frobenius on the conductor 1951 twist.
\\ Note that the value at p=1951 was found by trial and error.
\\ The error I made, that Junkie fixed, was that I thought that "pure thought" gave
\\ that bp(1951)=1.

bp(p)=if(p==1951,zet10^6,ap(p)*twist[p%1951])

\\ Here's the character of the form.

chi(p)=if(p==1951,0,twist[p%1951]^2);

\\ Now v will be the vector of Hecke eigenvalues.

v=vector(N,i,0);
v[1]=1;
for(i=2,N,if(i%1000==0,print(N-i));fac=factor(i);k=matsize(fac)[1];\
   if(k>1,v[i]=prod(j=1,k,v[fac[j,1]^fac[j,2]]),\
   if(fac[1,2]==1,v[i]=bp(i),\
   p=fac[1,1];e=fac[1,2];v[i]=v[p]*v[p^(e-1)]-chi(p)*v[p^(e-2)]))\
);

\\ Now glue them together with the Maass form analogue of a q-expansion.
\\ The function q^n is replaced by sqrt(y)*K(2.Pi.n.y)*sin(2.Pi.n.x)
\\ where K=K_0 is a certain Bessel function.

F(z)=local(x,y,M);x=real(z);y=imag(z);M=ceil(11/y);if(M>N,error("y too small."));sqrt(y)*sum(n=1,M,if(v[n]==0,0,v[n]*besselk(1e-30*I,2*Pi*n*y)*sin(2*Pi*n*x)))

Si cortas y pegas eso en una sesión de pari, entonces se quedará sin memoria, o pasará un largo tiempo contando desde 2.000.000 hasta 0 y luego se detendrá. Si se queda sin memoria, escriba algo como "allocatemem(2000000000)" o algo parecido. Si no puede molestarse en esperar a que la cuenta atrás termine, entonces cambie la opción " $N=$ " a algo más pequeño al principio. Está computando la primera $N$ coeficientes de la forma Maass.

La forma en sí misma es la función $F$ definido en la última línea. Cuanto más grande sea la $N$ ser, la mejor precisión $F$ se definirá a, y por lo tanto, cuanto más cerca de la línea real será capaz de evaluar $F$ . La cuestión es que si $z$ está en el plano medio superior pero muy cerca de la línea real, entonces se necesitan muchos términos de $F$ para evaluar $F(z)$ con exactitud.

Como en el post de Junkie, podemos intentar una verificación aleatoria, para ver si $F$ parece que es invariable bajo $\Gamma_1(1951)$ :

(22:47) gp > zz4=0.0001+0.0001*I
%7 = 0.0001000000000000000000000000000 + 0.0001000000000000000000000000000*I
(22:47) gp > F(zz4)
%8 = 0.000000008720947535106531893030584189 - 0.00000001730002410267070596763764527*I
(22:48) gp > F(zz4/(1951*zz4+1))
%9 = 0.000000008720947535106531893415409754 - 0.00000001730002410267070596828184932*I
(22:51) gp > %8-%9
%10 = -3.848255644 E-28 + 6.442040498 E-28*I

La cuestión es que la diferencia entre los dos valores de $F$ es un número minúsculo, por lo que $F$ es, dentro de los errores de cálculo, aparentemente invariable bajo $\Gamma_1(1951)$ . Se puede jugar con otros ejemplos, pero se advierte: cuanto más pequeño $N$ es (es decir, cuanto menos tiempo se moleste en esperar al principio) más cuidado hay que tener (es decir, más probable es que el programa se niegue a computar $F$ donde quieres que lo calcule).

5voto

varunsrin Puntos 113

Un problema. En 4 horas Magma puede calcular esto:

> f:=x^24 - 50726*x^22 + 152929135*x^20 - 158664037068*x^18 + 63787035668165*x^16 - 
    9040633522810414*x^14 + 287094814384960835*x^12 - 2050038614413776542*x^10 +
    4291028045077743465*x^8 - 1381768039105642956*x^6 + 134981448876235615*x^4 -
    2922378139308818*x^2 + 14488688572801;
> L:=LSeries(NumberField(f) : Method:="Artin");

Pero el problema es que la descomposición de L`prod no contiene factores bidimensionales...

> [<Character(x[1]`parent),x[2]> : x in L`prod | #x[1]`lpoles eq 0];
[
    <( 3, 3, 0, -1, -zeta(5)_5^3 - zeta(5)_5^2, zeta(5)_5^3 + zeta(5)_5^2 + 1, 
       0, zeta(5)_5^3 + zeta(5)_5^2 + 1, -zeta(5)_5^3 - zeta(5)_5^2 ), 1>,
    <( 3, 3, 0, -1, zeta(5)_5^3 + zeta(5)_5^2 + 1, -zeta(5)_5^3 - zeta(5)_5^2, 
       0, -zeta(5)_5^3 - zeta(5)_5^2, zeta(5)_5^3 + zeta(5)_5^2 + 1 ), 1>,
    <( 5, 5, -1, 1, 0, 0, -1, 0, 0 ), 1>,
    <( 6, -6, 0, 0, 1, 1, 0, -1, -1 ), 2>
]

Así que logramos una descomposición $1+\rho_3+\bar\rho_3+\sigma_5+2\eta_6$ con esto, ¿pero no hay factor de grado 2? ¿Es esto correcto? Tal vez el campo de división es necesario aquí (cuando todos los representantes aparecerá). La suposición es que el cálculo anterior ya hizo esto en esencia, así que puedo volver a ejecutarlo después de añadir más inteligencia.

4voto

varunsrin Puntos 113

Voy a copiar un comentario mío sobre Jehanne aquí.

Su documento es: http://dx.doi.org/10.1006/jnth.2001.2656

Jehanne tiene 4 ejemplos 2-5 totalmente reales (página 353-6), incluido el de Booker. También obtiene el campo del grado 24. El método para calcular el $L$ -se enumeran en la sección 6.

Jehanne da el truco de Poonen como debido a Serre (del libro de Buhler). Booker dice lo mismo (página 332).

Como señaló Kevin Buzzard, la aplicación a la conjetura computacional de Stark era sólo para los complejos $A_5$ sin embargo, por lo que el cómputo explícito totalmente real podría no haberse realizado fuera de Booker.

No creo que Dokchitser haga nada en su código Magma que no sea en esencia lo mismo que tú haces aquí, ya que los tipos de división simplemente se realizan de forma alternativa, normalmente más complicada debido a la generalidad. Creo que su código probablemente podría hacer $A_4$ como una caja negra sin muchas modificaciones, pero $A_5$ requiere un campo numérico de grado 240.

Había olvidado dónde se había publicado la tesis de Booker. Gracias por el enlace.

3voto

varunsrin Puntos 113

En realidad, debería leer la ayuda de Magma, pues la dirección de ArtinRepresentations está ahí.

  f:=14488688572801 - 2922378139308818*x^2 + 134981448876235615*x^4 -
     1381768039105642956*x^6 + 4291028045077743465*x^8 - 
    2050038614413776542*x^10 + 287094814384960835*x^12 -
   9040633522810414*x^14 + 63787035668165*x^16 - 158664037068*x^18 + 
    152929135*x^20 - 50726*x^22 + x^24;
  K:=NumberField(f);
  A:=ArtinRepresentations(K); // Four hours for computing
  K`artinrepdata`Frob(7); // fails

Entonces A[2] tiene grado 2, aunque el conductor es $1951^2$ como tú dices. Magma hace el producto tensorial L-funciones, pero tal vez no automáticamente en los malos primos. Pero el problema es que no puede calcular el Frobenius para muchos primos, probablemente debido a la falta de distinción del tipo de ciclo. No creo que Magma sea la herramienta adecuada para un problema tan específico. Tal vez la función interna de Frob pueda ser anulada y funcione a través del truco de Serre-Poonen.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X