conditioning of the monomial basis It is well-known that the Vandermonde matrix is ill-conditioned when all nodes, i.e., the points that generate it, are real. I have read the following on a paper:

The ill-
conditioning of the Vandermonde basis has an elementary explanation
in cases where the points xj are unequal in size. If $|x_j |$ varies with j, then the powers
$|(x_j)^k|$ vary exponentially. This means that function information associated with smaller
values of $x_j$ will only be resolvable through exponentially large expansion coefficients,
and accuracy will quickly be lost in floating point arithmetic

I have not understood the explanation because he terms "function information" and "expansion coefficients" are not defined in this context on the paper. Could somebody provide an explanation on why the Vandermonde matrix is ill-conditioned?.

The Vandermonde matrices are not necessarily ill-conditioned when the nodes are real.
Example: In the case of $n=2$ nodes, the Vandermonde matrix is
$$A = \begin{bmatrix} 1 & x_0 \\ 1 & x_1 \end{bmatrix}.$$
If $x_0 = -1$ and $x_1 = 1$, then columns are orthogonal and the 2-norm condition number of $A$ is given by
$$\kappa_2(A) = \|A\|_2 \|A^{-1}\|_2 = 1.$$
This is the smallest possible value.
In order to understand exactly why the Vandermonde matrices can be arbitrarily ill-conditioned we investigate the singular values. Let $$\sigma_1 \ge \sigma_2 \dots \ge \sigma_n \ge 0$$
denote the singular values of $A$. Then $\{\sigma_i^2\}_{i=1}^n$ are the eigenvalues of $A^TA$ and $$\kappa_2(A) = \frac{\sigma_1}{\sigma_n}.$$ We can establish simple bounds on the singular values as follows. For a general symmetric matrix $B = [b_{ij}]$, the smallest eigenvalue satisfies
$$ \lambda_{\min}(B) = \min \{ x^T B x \: : \: \|x\|_2 = 1 \} \leq b_{ii} $$
for all $i$ and the largest eigenvalue satisfies
$$ \lambda_{\max}(B) = \max \{ x^T B x \: : \: \|x\|_2 = 1 \} \ge b_{ii} $$
for all $i$.
Example: In the case of $n=2$ nodes, we have
$$ A^T A = \begin{bmatrix} 2 & x_0+x_1 \\ x_0 + x_1 & x_0^2 + x_1^2 \end{bmatrix}.$$
It follows that $$ \sigma_2^2 \leq 2, \quad x_0^2 + x_1^2 \leq \sigma_1^2$$
We conclude that $$\kappa_2(A) \ge \sqrt{ \frac{x_0^2 + x_1^2}{2} }.$$
The case of $x_1 = -x_0 = -x$ is particular clear. In this case, the relevant linear system is
$$ \begin{bmatrix} 1 & x \\ 1 & -x \end{bmatrix} \begin{bmatrix} c_0 \\ c_1 \end{bmatrix} = \begin{bmatrix} y_0 \\ y_1 \end{bmatrix}. $$
The solution is
$$ \begin{bmatrix} c_0 \\ c_1 \end{bmatrix} = \frac{1}{2x} \begin{bmatrix} x & x  \\ 1 & -1\end{bmatrix} \begin{bmatrix} y_0 \\ y_1 \end{bmatrix} = \begin{bmatrix} \frac{y_0+y_1}{2} \\ \frac{y_1-y_0}{2x} \end{bmatrix}.$$
Here we note that while $c_0 = \frac{y_0+y_1}{2}$ is insensitive to small relative changes in the right-hand side when $y_1$ and $y_2$ are nearly equal, the divided difference
$$c_1 = \frac{y_1 - y_0}{2x}$$ is extremely sensitive to small relative changes in the right-hand side when $y_1$ and $y_2$ are nearly equal. Specifically, a tiny relative change can flip the sign of $c_1$. The situation is reversed when $y_1 + y_2 \approx 0$. In both cases, we see that a small componentwise relative change in the right hand-side $y$, can cause a large componentwise relative change in the solution $c$.

The analysis presented here extends to the general case of $n>2$. Here we have $$\sigma_n^2 \leq n, \quad \sum_{j=0}^{n-1} (x_j^{n-1})^2 \leq \sigma_1^2$$
and 
$$ \kappa_2(A) \ge \sqrt{\frac{\sum_{j=0}^{n-1} (x_j^{n-1})^2}{n }}.$$
