Let $p$ be a prime number and $N$ a positive integer not divisible by $p$. 


For some easy choices of $p$ and $N$, can anybody provide me with explicit examples of collections $$\{f\_k,\quad 2\leq k \leq k\_0 \}$$ of q-expansions such that, for each $k$, $f\_k$ is the $q$-series of a classical eigenform of weight $k$ and level $Np$, all being members of an ordinary $p$-adic Hida family of tame level $N$? Here $k\_0$ is some reasonable upper bound for the weight, and I leave to the reader the meaning of the term "reasonable". Even if $k\_0=4$ I will already be happy, and it is also fine if $k$ only runs among some of the integers between $2$ and $k\_0$. But of course, the cardinal of $\{ f\_k\}$ should be greater than $1$!


I am interested in Hida families which are neither Eisenstein not CM, as explicit examples of those are already well-known.


The simplest example seems to arise when one takes $N=1$ and $p=11$. In this case the associated Hida Hecke algebra is $\mathbf{T}=\Lambda:=\mathbb{Z}\_p[[T]]$. I guess that people like Kevin Buzzard, Robert Pollack, William Stein and many others have computed the $q$-expansions $f\_k$ in this case, up to some reasonable $k\_0$, and I would be happy to have access to these computations.


Yet another interesting example: reading the paper "How can we construct abelian Galois extensions of basic number fields?" by Barry Mazur, he reports on an example computed by Citro and Stein, where $N=1$ and $p=691$. They compute explicitly $\{f\_2,f\_{12}\}$ where $f\_{12}$ is the $p$-stabilization of the discriminant modular form $\Delta$ and $f\_2$ turns out to be the only weight 2 newform of level $p$ and character $\omega^{-10}$, where $\omega$ is the Teichmuller character. I would like to know the $q$-expansions $f\_k$ in this case for some other integers in between $2$ and $12$. (For integers $k$ larger than $12$ one can possibly compute $f\_k$ by multiplying $f\_{12}$ by a suitable Eisenstein series of weight $k-12$ and applying Hida's ordinary projector.) 


And of course, any other collection of examples is welcome!



Examples of q-expansions in a Hida family
In the case when the Hida algebra is simply $\Lambda$, one can use *families* of overconvergent modular symbols to compute the $q$-expansion of the Hida family where the coefficients are functions of the weight. 


The idea is the following: form a "random" family of overconvergent modular symbols -- that is modular symbols with values in a distribution module tensor power series in the weight. Then iterate $U\_p$ on this family until you are in the ordinary subspace (modulo the accuracy at which you are working). Since the Hida algebra is just $\Lambda$, the result will be an eigensymbol. Its Hecke-eigenvalues will then be functions of the weight and these are the coefficients of the formal $q$-expansion I mentioned above.


With this $q$-expansion in hand, just plug in your favorite weight and you'll get the form in the Hida family of that weight (of course, computed modulo your fixed accuracy). 


This approach actually came out of an Arizona Winter School 2011 student project. Currently three graduate students at University of Madison (Lalit Jain, Marton Hablicsek, and Daniel Ross), together with Rob Harron and myself, are writing up a paper on this. 


Our code is still very much in beta form, but hopefully this example works. Take $p=3$ and $N=5$. Then the first few Hecke-eigenvalues (mod $3^3$) of the associated Hida family are:


$$
a\_2(w) = 12w^2 + 7w - 1,
$$
$$
a\_3(w) = 21w^2 + 11w - 1,
$$
and
$$
a\_5(w) = 16w^2 + 17w + 1.
$$
Here $w$ is not exactly the weight variable. Instead, to specialize these eigenvalues to weight $k$, one simply sets $w = 4^{k-2} - 1$. (Here is 4 is a topological generator of $1+3{\mathbb Z}\_3$.)


Here's the example with $N=1$ and $p=11$ -- below I'm working modulo $11^8$:


$$
a\_2(w) = 37w^7 + 880w^6 + 12388w^5 + 151975w^4 + 385840w^3 + 10344442w^2 + 40094463w + 857435522,
$$
$$
a\_3(w) = 71w^7 + 681w^6 + 8325w^5 + 94314w^4 + 220797w^3 + 2758794w^2 + 16210985w + 428717761,
$$
and
$$
a\_5(w) = 12w^7 + 472w^6 + 2245w^5 + 44445w^4 + 30443w^3 + 14355835w^2 + 180294533w + 1.
$$
As before, setting $w = (1+11)^{k-2}-1$ specializes to weight $k$. For instance, taking $k=12$ one should recover the first few Fourier coefficients of $\Delta$ modulo $11^8$.


