PeterDonis said:
For the case under discussion, yes, I believe Raychaudhuri's equation gives ##\dot{\Theta} \neq 0##.
Just to confirm this, here's a quick computation of the expansion scalar ##\theta## (I'll stick to the lower-case Greek letter here since that's the standard symbol) and the Raychaudhuri equation for the Bell congruence.
We'll work in the inertial frame in which the ships are initially at rest. In that frame we have the following for the coordinates ##(t, x)## of a given ship [Edit: I've normalized so that the proper acceleration of each ship is ##1##; that means ##t## and ##x## are essentially dimensionless coordinates, in ordinary units what I'm writing as ##t## and ##x## would be ##at## and ##ax##, where ##a## is the proper acceleration]:
$$
(t, x) = (\gamma v, x_0 - 1 + \gamma)
$$
where ##\gamma = \cosh \tau## and ##\gamma v = \sinh \tau## (##\tau## is the proper time along any ship's worldline [Edit: again, this is normalized, in ordinary units what I'm writing as ##\tau## would be ##a \tau##]), and ##x_0## is the ##x## coordinate at which the ship starts (i.e., its ##x## coordinate before it turns on its engines and begins accelerating). Taking derivatives with respect to ##\tau## gives us the 4-velocity ##(u^t, u^x## and the 4-acceleration ##(a^t, a^x)##:
$$
(u^t, u^x) = (\gamma, \gamma v)
$$
$$
(a^t, a^x) = (\gamma v, \gamma)
$$
It's straightforward to eliminate ##\tau##, by observing that ##\tau = \sinh^{-1} t## and working through the math accordingly; we find that ##\gamma## and ##\gamma v## are functions of ##t## only, such that ##\partial_t \gamma = v## and ##\partial_t \left( \gamma v \right) = 1##. (It's actually evident that these relations must hold from looking at the ##(t, x)## vector above and realizing that ##v = dx / dt##.)
Armed with all this, computing the expansion scalar is simple, because we're in a global inertial frame so all the connection coefficients are zero; we get
$$
\theta = \partial_a u^a = \partial_t u^t = v
$$
So ##\theta > 0## as soon as the ships start moving. The Raychaudhuri equation (with terms that are identically zero not shown) is:
$$
\dot{\theta} = - \frac{\theta^2}{3} + \partial_a a^a = - \frac{v^2}{3} + \partial_t a^t = 1 - \frac{v^2}{3} = \frac{1}{\gamma^2} + \frac{2}{3} v^2
$$
If this is correct, then ##\dot{\theta}## starts at ##1## but approaches a limit of ##2/3## as ##t \rightarrow \infty##.
The only thing I'm not sure about is the factor of ##1/3## in the first term; it's possible that this is due to a different definition of ##\theta^2## than just the square of the scalar ##\theta## that I calculated above. It would be nice and neat if that factor of ##1/3## were not there, since that would make ##\dot{\theta} = 1 / \gamma^2##, which looks nicer (and would make it approach ##0## as ##t \rightarrow \infty##, which is also neater). But as far as I can tell the factor of ##1/3## should be there given the definitions I'm using.