Suppose a matrix $A$ has linearly independent columns and QR factorization $A = QR$.
I was able to prove that $A^\top A = R^\top R$.
$$G = A^\top A = (QR)^\top(QR) = R^\top Q^\top QR = R^\top (Q^\top Q)R = R^\top R.$$
Apparently it is true that the angles between columns $i,j$ of $A$ are equivalent to the angles between columns $i,j$ of $R$, but I've double checked my reading and can't seem to figure out why.
I'm guessing it may have something to do with norms. I recall that angle $\theta$ between two vectors $x$ and $y$ is $$\theta = \arccos\frac{x^\top y}{\|x\|\|y\|}$$ and know how the Gram matrix is formed simply using inner-products, but haven't been able to find a connection yet.
You are very close. Note that, for any matrix $B$, the Gram matrix for the columns of $B$ is just $B^\top B$; you can verify this using the definition of matrix multiplication.
If $v_1, \ldots v_n$ are the column vectors of $A$, and $\theta_{ij}$ is the angle between $v_i$ and $v_j$, then according to your formula, $$\theta_{ij} = \arccos\frac{v_i^\top v_j}{\sqrt{v_i^\top v_i} \sqrt{v_j^\top v_j}} = \arccos\frac{(A^\top A)_{ij}}{\sqrt{(A^\top A)_{ii}} \sqrt{(A^\top A)_{jj}}}.$$ This is true, since $A^\top A$ is the gram matrix for $v_1, \ldots, v_n$.
Correspondingly, if $w_1, \ldots, w_n$ are the columns of $R$, with $\psi_{ij}$ being the angle between $w_i$ and $w_j$, we similarly see
$$\psi_{ij} = \arccos\frac{(R^\top R)_{ij}}{\sqrt{(R^\top R)_{ii}} \sqrt{(R^\top R)_{jj}}} = \arccos\frac{(A^\top A)_{ij}}{\sqrt{(A^\top A)_{ii}} \sqrt{(A^\top A)_{jj}}} = \theta_{ij}.$$