Consider a Lorentzian manifold completely covered by the following coordinate system:
$r\in[-\frac{L}{2},\infty) \\ \theta\in[0,\pi]\text{ when }r\ge\frac{L}{2} \\ \theta\in[0,\frac{\pi}{2}]\text{ when }r\lt\frac{L}{2} \\ \phi\in[0,2\pi) \\ t\in(-\infty,\infty)$
Now, we identify each point $(-\frac{L}{2},\theta,\phi,t)$ with the point $(\frac{L}{2},\pi-\theta,\phi,t)$. (Remember, this only includes $\theta\in[0,\frac{\pi}{2}]$.)
Then, we impose the following metric. This metric is just an ansatz to come up with the basic shape, so I'm not too worried about the particular formulas -- in a second I'll explain what's really going on.
$g_{tt} = -(1-\frac{r^2}{\alpha^2}) \\ g_{rr} = (1-\frac{r^2}{\alpha^2})^{-1} \\ g_{\theta\theta} = e^{-kr^2}\sin(\frac{5\pi}{2}e^{-r^2}) + r^2 - C \\ g_{\phi\phi} = [a(r)\sin^2{2\theta} + b(r)\sin^2{\theta}]\ g_{\theta\theta}$
where $k$ and $C$ are chosen to give $g_{\theta\theta}$ two zeros at $r=\pm\frac{L}{2}$, and
$a(r) = \begin{cases} 1 & r\le \frac{L}{2} \\ 1 - e^{-\frac{1}{(r-\frac{L}{2})^2}} & r\gt \frac{L}{2} \end{cases} \\ b(r) = \begin{cases} 0 & r\le \frac{L}{2} \\ e^{-\frac{1}{(r-\frac{L}{2})^2}} & r\gt \frac{L}{2} \end{cases} = 1 - a(r)$
You can see the shape of these functions here, where I have renamed $g_{\theta\theta}$ to $g_{hh}$.
To understand what this manifold is, first note that, because $g_{rr}$ is only a function of $r$, the $r$ coordinate lines are also geodesics, at least with respect to a spatial slice (hypersurface of constant $t$). Now consider the cross-sections of these radial geodesics (surfaces of constant $r$ and $t$). As $r \to \infty$, such a surface becomes a sphere. But as $r \to \frac{L}{2}^+$, it becomes a pair of spheres, because of the $\sin^2{2\theta}$ in $g_{\phi\phi}$ -- one for $\theta\in [0,\frac{\pi}{2}]$ (call this $S_1$) and another for $\theta\in [\frac{\pi}{2},\pi]$ (call it $S_2$). And when $r\lt\frac{L}{2}$, it is again a single sphere, because we lose half the range of $\theta$. These single spheres are the natural continuation of the $S_1$ across $r=\frac{L}{2}$, but they are also the continuation of the $S_2$ across the identification mentioned earlier!
In other words, in spite of its strange coordinate system, this manifold is both smooth and closed. What's really going on is that the region $r\in[-\frac{L}{2}, \frac{L}{2}]$ is a "bridge" that connects $\theta\in[0,\frac{\pi}{2}]$ to $\theta\in[\frac{\pi}{2},\pi]$ via the point at $r=\frac{L}{2}$ -- which is a single point because $g_{\theta\theta}=0$ here -- this point is the origin of a spherical coordinate system where $r\in[\frac{L}{2},\infty)$. Since both ends of the bridge are at the same point, it can also be thought of as a loop.
Note that this manifold is asymptotically a de Sitter space, according to this formula. (Except for the offset of $C$ in $g_{\theta\theta}$, but I don't think that's really an issue as long as it varies like $r^2$.)
Now here's what I want to do to this manifold. In spite of the loop that it contains, I want to distort the metric in such a way as to make it "topologically neutral", in the sense that, when you take its Chern-Gauss-Bonnet integral, you get the same value as de Sitter space alone, if it didn't have this loop structure in it. And I want this spacetime to remain both smooth and stationary, ie. $g_{\mu\nu}$ should be independent of $t$, though I presume it will now contain off-diagonal terms. While doing this, we can allow $L$ to vary, but $\alpha$ should be taken as given.
This is inspired by the idea that you can induce a loop in a flat belt without lifting its ends, as long as there is a compensating $360^{\circ}$ twist along the belt.
So how might the "twist" manifest, if the analogy holds? I was imagining that the eigenvectors of the Ricci tensor could twist about those radial geodesics, relative to the Levi-Civita connection (eg. if $\vec v$ is a Ricci eigenvector field, then $\nabla_r{\vec v}$ is roughly perpendicular to both $\vec v$ and $\hat r$).
But aside from such vague notions, I haven't the slightest clue how to try modifying this metric to achieve the desired "topological neutrality". As far as I'm concerned, any modification is as good a guess as any other. Yet I cannot possibly get anywhere without narrowing down the field of candidates.
I don't need a closed form solution, either; I assume this is more likely to be resolved by a numerical simulation. Also, I'm very interested in whether the condition of "topological neutrality" limits the possible solutions to a discrete set. But I also want to find the solution(s) themselves, approximate or not.
So any insight or general thoughts are welcome and greatly appreciated! And thank you for taking the time to read this!