As we know, to derive the Schwarzschild solution in GR we want to determine the unknown functions $a,b$ (sometimes $c$) of $$(ds)^2 = a(r,t)\;dt^2 - b(r,t)\; dr^2 - r^2d\Omega^2 $$ or $$ (ds)^2 = a(r,t)\;dt^2 + 2b(r,t)\;drdt - c(r,t)\; dr^2 - r^2d\Omega^2 $$ where $d\Omega^2 = d\theta^2 + \sin^2 \theta \; d\phi^2 \; $ is the usual line element of the 2-sphere. When applied to the empty space field equations where $G_{ab} \equiv 0$, Birkhoff's theorem says any static, spherically symmetric solution to the vacuum equations of general relativity is necessarily the Schwarzschild solution up to a coordinate transformation. That basically wraps up looking for exterior spherically symmetric solutions in GR. The next logical step is then to look for interior solutions which is a whole different ball game.
I'm interested in considering axially symmetric spacetimes in GR. The usual approach is to look at Weyl metrics (as mentioned in the comments below) and go for broke.
What I am confused about is trying to consider one of the following line element $$ (ds)^2 = a(r,\theta)\;dt^2 + 2b(r,\theta)\;drdt - c(r,\theta)\; dr^2 - r^2d\Omega^2. $$
The above line element is axially symmetric also correct? In terms of polar spherical coordinates it is symmetric about the polar angle coordinate.
In terms of a spacetime that models a non-spherically symmetric body for e.g. a planet, I would have thought that this line element is more suitable than the general Weyl metric. However, I have never seen a method beginning with the above line element.
Question Is the Weyl solution to axially symmetric gravitational fields as Schwarzschild solutions are to spherically symmetric gravitational fields? In other words, is there a Birkhoff-type theorem for axially symmetric solutions to the field equations of GR?