Thanks, Kaveh, for wishing to look at chapters on proof complexity!
Concerning Robin's question, first that note $AC^0$ contains functions requiring formulas (and even circuits) of size $n^k$ for any constant $k$. This follows, say, from a simple fact that $AC^0$ contains all DNFs with constantly long monomials. Thus, $AC^0$ contains at least $\exp(n^k)$ distinct functions, for any $k$. On the other hand, we have at most about $\exp(t\log n)$ functions computable by formulas of size $t$.
I shortly discussed the issue of getting explicit lower bounds of $n^2$ or larger with Igor Sergeev (from Moscow university). One possibility could be to use Andreev's method, but applied to some another, easier computable function instead of Parity. That is, consider a function of $n$ variables of the form $F(X)=f(g(X_1),\ldots,g(X_b))$ where $b=\log n$ and $g$ is a function in $AC^0$ of $n/b$ variables; $f$ is some most complex function of $b$ variables (mere existence of $f$ is enough). We only need that the function $g$ cannot be "killed" in the following sense: if we fix all but $k$ variables in $X$, then it must be possible to fix all but one of the remaining variables of $g$ so that the obtained subfunction of $g$ is a single variable. Then applying Andreev's argument and using Hastad's result that the shrinking constant is at least $2$ (not just $3/2$ as earlier shown by Sybbotovskaya), the resulting lower bound for $F(X)$ will be about $n^3/k^2$. Of course, we know that every function in $AC^0$ can be killed by fixing all but $n^{1/d}$ variables, for some constant $d\geq 2$. But to get a $n^2$ lower bound it would be enough to find an explicit function in $AC^0$ which cannot be killed by fixing all but, say, $n^{1/2}$ variables. One should search for such a function in depth larger than two.
Actually, for the function $F(X)$ as above, one can obtain lower bounds about $n^2/\log n$ via simple greedy argument, no Nechiporuk, no Subbotovskaya and no random restrictions! For this, it is merely enough that the "inner function" g(Y) is non-trivial (depends on all its $n/b$ variables). Moreover, the bound holds for any basis of constant fanin-gates, not just for De Morgan formulas.
Proof: Given a formula for $F(X)$ with $s$ leaves, select in each block $X_i$ a variable which appears the smallest number of times as a leaf. Then set all remaining variables to the corresponding constants so that each $g(X_i)$ turns to a variable or its negation. The obtained formula will then be at least $n/b$ times smaller than the original formula. Thus, $s$ is at least $n/b=n/\log n$ times the formula size $2^b/\log b=n/\log\log n$ of $f$, that is, $s\geq n^{2-o(1)}$. Q.E.D.
To get $n^2$ or more, one has to incorporate Subbotovskaya-Hastad shrinking effect under random restrictions. A possible candidate could be some version of Sipser's function used by Hastad to show that depth-$(d+1)$ circuits are more powerful than those of depth $d$.