The law of iterated logarithm asserts that if $x_1,x_2,\dots$ are i.i.d $\cal N(0,1)$ random variables and $S_n=x_1+x_2+\cdots+x_n$, then
$$\limsup_{n \to \infty} S_n/\sqrt {n \log \log n} = \sqrt 2, $$ almost surely.
1. Random matrices
Let $(x_{ij})_{1 \le i < \infty, 1 \le j \le i }$ be an array of i.i.d random $\cal N(0,1)$ variables, define $x_{ji}=x_{ij}$ and let $M_n$ be the $n$ by $n$ matrix obtained by restricting the array to $1 \le i \le n, 1 \le j \le n$.
Let $\lambda(n)$ denote the maximum eigenvalue of $M_n$. The variance and limiting distribution of the maximum eigenvalue (for these and related classes of matrices) were discovered by Tracy and Widom.
Question 1: What is the law of iterated logarithm for this scenario? Namely, can one identifies functions $F(n)$ and $G(n)$ such that
$\limsup_{n \to \infty} (\lambda (n)- {\bf E}(\lambda(n))/ F(n) = 1,$ almost surely, and
$\liminf_{n \to \infty} (\lambda (n)- {\bf E}(\lambda(n))/ G(n) = -1$ almost surely.
Update (May, 24, 2015): A full answer for the limsup (including constants) and a partial answer for the liminf was achieved by Elliot Paquette and Ofer Zeitouni: arxiv.org/abs/1505.05627 !
2. Random permutations
Choose a sequence of real numbers in the unit interval at random independently.
Use the first $n$ numbers in the sequence to describe a random permutation $\pi$ in $S_n$. Look at the length $L_n$ of the maximum increasing sub-sequence of $\pi$. The expectation, variance and limiting distributions of the maximal increasing subsequence are known by the works of Logan-Shepp, Kerov-Vershik, (expectation) Baik, Deift & Johansson (variance and limiting distribution) and many subsequent works.
Question 2: What is the law of iterated logarithm for this scenario?
3. Motivation.
My initial motivation is to try to understand the law of iterated logarithm itself, and more precisely, to understand why do we have the $\sqrt {\log \log n}$ factor. (I cannot rationally expect to understand the $\sqrt 2$.)
Best Answer
My initial answer was wrong. Instead of completely deleting it, I left it at the bottom of the post. I use notation now that slightly differ from my original post.
As before, I give only a partial answer related to the max eigenvalue question and the limsup. Recall that the tail estimates for the max eigenvalue is $$(P(\lambda_1(n) >2\sqrt{n}+tn^{-1/6})\sim e^{-Ct^{3/2}}$$ (see e.g. Ledoux's reviews or http://math.univ-lyon1.fr/~aubrun/recherche/smalldeviations/smalldeviations.pdf).
Now, the question is what is the scale in which the recentered and rescaled maximal eigenvalues start behaving like independent random variables. In my original post I essentially claimed this starts happening at scale roughly $n$. This is wrong: the correct scale is given by the so-called GUE minor process (described by Baryshnikov in 2001). Relevant to Gil's problem is the scaling limit computed by Forrester and Nagao: http://arxiv.org/pdf/0801.0100.pdf. They show that the correct scaling in $n^{2/3}$, that is that the (centered and rescaled) maximal eigenvalues of the GUE(n) and GUE(m) decorelate if $m-n>>n^{2/3}$ and are strongly correlated in $m-n<<n^{2/3}$. Using that, and Borel Cantelli, one gets the correct scale. The limsup can be taken over the sequence $a_n=n^3$ instead of $n$. If you go over that sequence, you get from Borel-Cantelli that you cannot exceed $t_{a_n}a_n^{-1/6}$ infinitely often as long as $$\sum_n e^{-C t_{a_n}^{3/2}}<\infty.$$ This will happen as soon as $t_{a_n}^{3/2}>C'\log n$, i.e. $t_{a_n}\sim C'' (\log n)^{2/3}$. Unravelling the definitions, this will give you an upper bound on the limsup that is $t_n\sim (\log n)^{2/3}$; that is, this will show that the function $F(n)$ is at most $C (\log n)^{2/3}/n^{1/6}$. The lower bound on the limsup will follows a similar pattern, using the asymptotic independence stated above, and will give the same function $F(n)$. One can also compute constants (for the a.s. convergence) in the same manner.
Filling in the details in the above is beyond a mathoverflow answer, I may return to it in detail later.
I speculate that the liminf will use the lower tail estimates, which are of the form $e^{-Ct^{3})}$. Working with the same argument would then give $(\log n)^{1/3}/n^{1/6}$ as the correct scaling.
OLD (WRONG) POST:
The following is a partial answer with some hints on how to complete it; I may revisit it later with more details. The answer below refers to the max eigenvalue question, to the limsup, and does not attempt to get the constants, only the scale.
First, the upper bound: recall that the tail estimates for the max eigenvalue is $$(P(\lambda_1(n) >2\sqrt{n}+t\sqrt{n})\sim e^{-Cnt^{3/2}}$$ (see e.g. Ledoux's reviews or http://math.univ-lyon1.fr/~aubrun/recherche/smalldeviations/smalldeviations.pdf). Also recall (a kindly neighbor next door, A. Dembo, pointed out the usefulness of this observation in the current context) that the sequence $\lambda_1(n)$ is monotone in $n$. So, the limsup can be taken over a sequence $a_n=(1+\epsilon)^n$ instead of $n$. If you go over that sequence, you get from Borel-Cantelli that you cannot exceed $t_{a_n}\sqrt{a_n}$ infinitely often as long as $$\sum_n e^{-C t_{a_n}^{3/2} a_n}<\infty.$$ This will happen as soon as $a_nt_{a_n}^{3/2}>C'\log n$, i.e. $t_{a_n}\sim C'' (\log n/a_n)^{3/2}$. Unravelling the definitions, this will give you an upper bound on the limsup that is $t_n\sim (\log \log n /n)^{2/3}$; that is, this will show that the function $F(n)$ is at most $C (\log \log n)^{2/3}/n^{1/6}$.
(Remark to @Carlo Beenakker: the difference with TSAW is that the tail estimates there are proportional to $t^3$, not $t^{3/2}$; This will be in line with the liminf computation in the RM case, in which case I expect the same discrepency in exponents of the $\log \log n$ as you noted, see below.)
For the lower bound (on the limsup), you need to take a sequence that is sparse enough so that you have approximate independence. Here is a leap of faith, that would require a separate proof: if you take an $N\times N$ GUE and resample the top $\epsilon N\times \epsilon N$ entries, the maximum eigenvalue does not change much if $\epsilon$ is small. I believe this follows from delocalization results for the eigenvectors of GUE, but have not checked details. If that is true, then for an exponentially growing sequence $a_n=C^n$, the sequence $\lambda_1(a_n)$ will behave like an i.i.d. sequence, to which you now can again apply the independent Borel-Cantelli to match the previous upper bound.
I speculate that the liminf will use the lower tail estimates, which are of the form $e^{-C(nt^{3/2})^2}$. Working with the same argument would then give $(\log \log n/n^2)^{1/3}$ as the correct scaling.