Calculating the Riemann tensor for a 3-Sphere I have worked out all the connection symbols for the 3-sphere using calculus of variations, cf. this Phys.SE post. So to find the Riemann tensor I am trying to find all the nonzero components of:
\begin{equation}
R^\rho_{\sigma \mu \nu} = \partial_\mu\Gamma^\rho_{\nu\sigma}-\partial_\nu\Gamma^\rho_{\mu\sigma}+\Gamma^e_{\mu\lambda}\Gamma^\lambda_{\nu\sigma}-\Gamma^\rho_{\nu\lambda}\Gamma^\lambda_{\mu\sigma}
\end{equation}
The only strategy I can see is trying every possible combination of Christoffel symbols that  produces a nonzero component (Which is of course possible). Surely there must be a better approach than this? What is the general strategy for finding the Riemann tensor?
 A: Here's a way to find the Riemann tensor of the 3-sphere with a lot of intelligence but no calculations.
At any point $p$ on a sphere, all directions look the same. Therefore there can be no privileged vector at a point $p$. Now consider the eigenvalue problem for the Ricci tensor, $$R^\alpha{}_\beta x^\beta = \lambda x^\alpha.$$
Since no vector is better than any other vector, either no vector is an eigenvector, or every vector is. Since the Ricci tensor is symmetric, it must be the latter. But the only way for every vector to be an eigenvector is if the tensor is proportional to the identity tensor. Hence, $$R^\alpha{}_\beta = \lambda \delta^\alpha{}_\beta \Rightarrow R_{\alpha\beta} = \lambda g_{\alpha\beta}.$$
Now take the trace of both sides to conclude that $\lambda = R/3$ where $R$ is the curvature scalar.
Now we need to be even more clever. A diagonal element $R_{\alpha\beta}x^\alpha x^\beta$ of the Ricci tensor in $n$ dimensions is $(n-1)$ times the average of the sectional curvature over planes containing $x^\alpha$, for $x_\alpha x^\alpha = 1$. The important thing is that the sectional curvature is the Gaussian curvature of the 2-surface generated by geodesics starting at $p$. What can this 2-surface be? It's a 2-surface characterized by constant curvature (because of all the symmetry we have) -- that's (some part of) a 2-sphere! The Gaussian curvature of a 2-sphere is $1/r^2$. Thus $2/r^2 = \lambda g_{\alpha\beta}x^\alpha x^\beta = \lambda$ and we can conclude $$R_{\alpha\beta} = \frac{2}{r^2} g_{\alpha\beta}.$$
A: 
Surely there must be a better approach than this? 

Certainly; it's the Cartan formalism which employs differential forms. Consider your case of a sphere, with a metric tensor,
$$ds^2 = r^2dr^2+r^2\sin^2 \theta \, d\theta^2$$
We can choose an orthonormal basis $e^a$ such that $ds^2 = \eta_{ab}e^{a}e^{b}$, and $e^a = e^a_\mu dx^\mu$. For our case:
$$e^r = r \, dr, \, \, e^\theta = r\sin \theta \, d\theta$$
We now compute the exterior derivative of each, and re-express them in terms of the basis:
$$de^r = 0, \quad de^\theta = \sin \theta \, dr \wedge d\theta = \frac{1}{r^2} \, e^r \wedge e^\theta$$
We can now use the first structure equation,$^{1}$ namely $de^a + \omega^a_{b} \wedge e^b = 0$, to deduce the connection components $\omega^a_b$. With some experience, you'll be able to read them off:
$$\omega^r_r = \omega^\theta_\theta = 0, \quad \omega^\theta_r = \frac{1}{r^2}e^\theta = \frac{1}{r}\sin \theta \, d\theta$$
The second structure equation allows us to directly compute components of the Ricci tensor $R^{a}_b$ in the orthonormal basis, namely, $R^a_b = d\omega^a_b + \omega^a_c \wedge \omega^c_b$. We find that,
$$R^\theta_\theta = R^r_r = 0, \quad R^\theta_r = \frac{1}{r^4} e^\theta \wedge e^r$$
Given that $R^a_{b} = R^a_{bcd} e^c \wedge e^d$, we can compute the terms of the full Riemann tensor. To convert back from our orthonormal basis, we can use the relation
$$R^\lambda_{\mu \nu \sigma} = ( e^{-1})^{\lambda}_a \, R^a_{bcd} \, e^b_\mu \, e^c_\nu \, e^d_\sigma$$
By this method, we find for example, $R^\theta_{r\theta r} = \sin^2 \theta$ as expected for the sphere. Generally, the Cartan formalism is much easier than direct computation, but be careful in particular when,


*

*Deducing the connections from the first structure equation; it's not always straightforward;

*Converting back to the orthonormal basis.


For an excellent exposition of the method, see the gravitational physics lectures by Professor Ruth Gregory available from the Perimeter Scholars website.

$1$ In fact, the structure equation states $de^a + \omega^a_b + e^b = T^a$, where $T^a$ is the torsion, but in general relativity, we may assume $T^a = 0$. See my own question: Why can we assume torsion is zero in GR?
A: I have just now finished an article, "Geometry of the 3-sphere", in which at the end of the paper I give a simple derivation of the Riemann curvature bivector for the unit 3-sphere, using (Clifford) geometric algebra. I also discuss the Lie group $SU(2)$ and Lie algebra $SU(2)$ on the unit 3-sphere, using the powerful, but still rather unknown geometric algebra $G(4)$, the Pauli geometric algebra $G(3)$, and the quaternions $H$, isomorphic to the even subalgebra of $G(3)$.
A: In 3 dimensions, all the data in the Riemman tensor is contained in the Ricci tensor as the Weyl tensor vanishes. So it suffices to compute the Ricci tensor and then using the decomposition given here.
