## Particle and photon orbits in McVittie spacetimes

Brien C Nolan

School of Mathematical Sciences, Dublin City University, Glasnevin, Dublin 9, Ireland. (Dated: August 4, 2014)

McVittie spacetimes represent an embedding of the Schwarzschild field in isotropic cosmological backgrounds. Depending on the scale factor of the background, the resulting spacetime may contain black and white hole horizons, as well as other interesting boundary features. In order to further clarify the nature of these spacetimes, we address this question: do there exist bound particle and photon orbits in McVittie spacetimes? Considering first circular photon orbits, we obtain an explicit characterization of all McVittie spacetimes for which such orbits exist: there is a 2-parameter class of such spacetimes, and so the existence of a circular photon orbit is a highly specialised feature of a McVittie spacetime. However, we prove that in two large classes of McVittie spacetimes, there are bound particle and photon orbits: future-complete non-radial timelike and null geodesics along which the areal radius r has a finite upper bound. These geodesics are asymptotic at large times to circular orbits of a corresponding Schwarzschild or Schwarzschild-de Sitter spacetime. The existence of these geodesics lays the foundations for and shows the theoretical possibility of the formation of accretion disks in McVittie spacetimes. We also summarize and extend some previous results on the global structure of McVittie spacetimes. The results on bound orbits are established using centre manifold and other techniques from the theory of dynamical systems.

## I. INTRODUCTION AND SUMMARY.

Perhaps one of the most striking features of black holes is their ability to create circular photon orbits: by travelling to the vicinity of a Schwarzschild black hole and settling in an orbit at a radius r = 3 m , where m is the mass of the black hole, a physicist can exploit the existence of circular photon orbits at this radius to look at the back of their own head. If a Schwarzschild black hole is not available, a charged or rotating black hole can be used for the same purpose. See e.g. Sections 20, 40 and 61 of [1]. Indeed circular photon orbits exist in Kerr-Newman-(anti)de Sitter spacetime for all parameter values corresponding to black holes [2]. Thus it seems reasonable to say that the existence of circular photon orbits (CPOs) is a generic feature of black holes.

Similarly, in the case of massive particles, the existence of an innermost stable circular orbit (ISCO) is a feature of black holes that is not generally present in Newtonian gravity (but see [3]). The existence of the ISCO is central to the formation of thin accrection disks around the black hole which, in turn, encodes useful information about the black hole. See [4].

The question arises as to whether or not CPOs and ISCOs arise in black holes in more general, non-vacuum settings. A first step here is to consider how such black holes can be modelled, and this question has arisen in the broader discussion of how the background expansion of the universe affects local systems [5] and in the more wide-ranging study of inhomogeneous cosmological models (see e.g. [6]). An early contribution to this discussion was McVittie's discovery of an intriguing solution of Einstein's equations that the author himself referred to as a 'mass-particle in an expanding universe' [7]. McVittie's solution represents an embedding of the Schwarzschild field in an isotropic cosmological background. In the form presented in [7], three families of line element were given, corresponding to the curvature index k = 0 , ± 1 of the isotropic background. There appear to be clear reasons to dispute the interpretation that McVittie's solutions with k = ± 1 represent some form of isolated system in an isotropic background [8, 9], and so we will focus on the spatially flat case k = 0 [10]. In this case, the line element may be written in the form [11]

$$\begin{smallmatrix} \bullet & \text{look} \\ \text{black} \\ \kappa & \text{hole} \end{smallmatrix} \quad d s ^ { 2 } = - ( f - r ^ { 2 } H ^ { 2 } ) d t ^ { 2 } - 2 r H f ^ { - 1 / 2 } d t d r + f ^ { - 1 } d r ^ { 2 } + r ^ { 2 } d \Omega ^ { 2 }, \\ \kappa & \text{hole} \quad. \quad \kappa \quad \kappa \quad \kappa \quad \kappa \quad \kappa \,.$$

where m&gt; 0 is a constant, f ( r ) = 1 -2 m/r &gt; 0 (i.e. we restrict to r &gt; 2 m ) and H = H t ( ) is the Hubble function of the isotropic background. That is,

$$H ( t ) = S ^ { \prime } ( t ) / S ( t ), \text{ \quad \ \ } ( 2 )$$

where S t ( ) is the scale factor S t ( ) of the spacetime obtained by setting m = 0 in (1):

$$\overset { \sim } { \text{n} } \quad d s ^ { 2 } = - ( 1 - r ^ { 2 } H ^ { 2 } ) d t ^ { 2 } - 2 r H d t d r + d r ^ { 2 } + r ^ { 2 } d \Omega ^ { 2 } \ ( 3 )$$

$$= - d t ^ { 2 } + S ^ { 2 } ( d x ^ { 2 } + x ^ { 2 } d \Omega ^ { 2 } ).$$

The spacetime with this line element will be referred to as the background , and we will use terms such as background metric in the obvious way. We note that t is a global time coordinate on the spacetime with line element (1): g ab t ,a t ,b = -f -1 &lt; 0. We set the time orientation of the spacetime by taking t to increase into the future. In (4), x = rS -1 is the usual co-moving radial coordinate of the isotropic spacetime. (We note that in (2) and throughout the paper, a prime means derivative with respect to argument.)

The limit m = 0 of (1) is well-defined, as it corresponds to setting the Weyl tensor of the spacetime to zero. The invariantly defined Newman-Penrose Weyl scalar Ψ 2 of the spherically symmetric line element (1) is given by

$$\Psi _ { 2 } = - \frac { m } { r ^ { 3 } }.$$

With H = 0, (1) is the line element of the Schwarzschild exterior. With H = H 0 = constant, the line element is that of Schwarzschild-de Sitter spacetime with cosmological constant Λ = 3 H 2 0 . These limits are likewise well-defined, being respectively the vacuum and Einstein-space limits of (1). In the case H = H 0 , the coordinate transformation T = t + u r ( ) with u ′ ( r ) = H r/ 0 ( √ f ( f -r 2 H 2 0 )) yields the more familiar form

$$d s ^ { 2 } = - ( 1 - \frac { 2 m } { r } - \frac { \Lambda } { 3 } r ^ { 2 } ) d T ^ { 2 } - ( 1 - \frac { 2 m } { r } - \frac { \Lambda } { 3 } r ^ { 2 } ) ^ { - 1 } d r ^ { 2 } + r ^ { 2 } d \Omega ^ { 2 }. \quad & \text{from $t$} \\ \mathfrak { r } _ { \dots \dots \dots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots$$

Imposing the Einstein equations with a cosmological constant, and with a perfect fluid energy-momentum tensor, yields expressions for the density ρ and pressure p :

$$8 \pi \rho = 3 H ^ { 2 } - \Lambda, & & ( 7 ) & \quad \text{$\kappa$}$$

$$8 \pi p = - 2 H ^ { \prime } f ^ { - 1 / 2 } - 3 H ^ { 2 } + \Lambda. \quad \ \ ( 8 ) \quad \text{fect}$$

Setting m = 0 yields the background density and pressure:

$$8 \pi \rho _ { 0 } = 3 H ^ { 2 } - \Lambda, & & ( 9 ) & & \text{$\mathbb{Psa}$} \\ \text{clu}$$

$$8 \pi p _ { 0 } = - 2 H ^ { \prime } - 3 H ^ { 2 } + \Lambda. \quad \ \ ( 1 0 ) \quad \ \ ( \text{ra} )$$

glyph[negationslash]

Assuming that H t ′ ( ) = 0, we see that McVittie spacetimes always have a scalar curvature singularity at r = 2 m .

Comment 1 It should be noted that since ρ = ρ t ( ) and p = p t, r ( ), there is no equation of state of the form p = g ρ ( ) in the McVittie spacetime. Thus the term 'perfect fluid' is not fully appropriate: we use it in the not uncommon sense of a fluid with isotropic pressure - the radial and tangential pressures are equal. On the other hand, as the background m = 0 spacetime is homogeneous and isotropic, one can appeal to an equation of state to close the Einstein equations and so govern the evolution. This is the perspective that we take: the insertion of the mass parameter m is done post hoc , after the Hubble function has been determined. In theory, one could then consider back reaction effects, redoing CMB and other calculations to take account of the mass parameter m . As we will argue below, this mass parameter provides a model of a (highly) localized inhomogeneity in an otherwise isotropic universe (we use the qualification to reflect the fact that the energy density of the universe is unaffected by m ρ : = ρ t ( ) = ρ 0 ( ) - there is no 'addit tional' matter in the universe).

We note that the transformation ( t, r ) → ( T, x ) with T = t and

$$r ( T, x ) = u ( 1 + \frac { m } { 2 u } ) ^ { 2 }, \ \ u ( T, x ) = x a ( T ) \quad \ ( 1 1 ) \quad \ \text{Tr} \quad \text{th} \\ \text{have}$$

can be used to write the line element (1) in comoving coordinates, which is the form originally given by McVittie [7]. In these coordinates, the curvature singularity at r = 2 m arises at u = m/ 2. The coordinate transformation is a diffeomorphism only if we restrict u to either (0 , m/ 2) or ( m/ , 2 ∞ ). Either interval provides a full cover of r ∈ (2 m, ∞ ). Note also that at fixed T , lim x →∞ r T, x ( ) = ∞ .

With these properties, it is tempting to conclude that the line element (1) corresponds to a spacetime that contains a black hole embedded in an isotropic universe. However, this interpretation is too simplistic. A correct interpretation requires the thorough study of the global properties of the spacetime, based on an analysis of its geodesics. Perhaps the dominant theme that has emerged from the various studies along these lines is that different outcomes emerge depending on the background scale factor S t ( ) and corresponding Hubble function H t ( ).

As far as the author is aware, the first study of the global structure of McVittie spacetimes was undertaken by Sussman as part of a comprehensive study of the global properties of spherically symmetric, shear-free perfect fluid spacetimes [12]. In this work, a variety of possible global structures was identified. However, the interpretation of some of the results must be questioned as Sussman takes limits x → 0 and x →∞ within the same spacetimes. This does not seem possible without the inclusion of the singularity r = 2 m as part of the spacetime (rather than part of its boundary).

For expanding universe models with a big bang at a finite time in the past - that is, when the isotropic background (3) has these features - the singularity r = 2 m forms a past boundary of the spacetime. This was first pointed out in [11] for backgrounds satisfying a linear equation of state p 0 = κρ 0 and with Λ = 0, and the result was generalised in [13] and [14]. As we will see below, this feature of McVittie spacetimes is not universal, but does arise if the big bang condition lim t → 0 + S t ( ) = 0 holds. When Λ = 0, the singularity r = 2 m also arises as a future boundary of the spacetime: future-pointing ingoing radial null geodesics run into this singularity in finite affine parameter time. However, as pointed out in [13], when there is a positive cosmological constant (and when some other technical conditions on H t ( ) hold), ingoing radial null geodesics meet a horizon rather that this singularity. Subsequently, [14] showed how the spacetime can extend through the horizon to a Schwarzschild-de Sitter spacetime. Together, these points indicate that with a background Λ-CDM model with Λ &gt; 0, McVittie spacetimes can indeed model black holes in expanding universes. In fact it was shown in [14] that the boundary includes a black hole horizon and a white hole horizon. However, when Λ = 0, it is very difficult to see how this interpretation can be given: all ingoing radial null geodesics either escape to infinity, or terminate at a scalar curvature singularity rather than reaching a horizon.

Thus, as emphasized in [15], McVittie spacetimes can have a variety of global structures depending on the scale factor of the background. In some cases, including the Λ-CDM model, black and white hole horizons arise. Returning to the question of how local systems are affected by cosmological expansion, it is clear that these McVittie spacetimes provide an interesting testing ground for such questions. On the other hand, studying e.g. the

existence of CPOs and ISCOs in McVittie spacetimes adds further to our understanding and correct interpretation of these spacetimes. So in this paper, we begin the study of particle and photon orbits in McVittie spacetimes by addressing these questions: do there exist CPOs in McVittie spacetimes? Do there exist bound particle and photon orbits in McVittie spacetimes? We obtain a complete answer to the first of these questions. That is, we explicitly determine all McVittie spacetimes that admit CPOs. Neither the Λ-CDM models of [13] and [14] nor the Λ = 0 models of [11] mentioned above are among these. Such McVittie spacetimes do not possess this characteristic feature of black hole spacetimes. However, both models do have the following feature: both admit bound particle and photon orbits. That is, there are futurecomplete, non-radial null and timelike geodesics γ in these spacetimes with the property that r | γ &lt; + ∞ , and lim s →∞ r | γ = r c where s is an affine parameter (or proper time) along the geodesic. The constant r c corresponds to the radius of a circular orbit in either a Schwarzschild or Schwarzschild-de Sitter background. This result is established for two classes of McVittie spacetimes which we define below. These classes include, respectively, the spacetimes of [13] and [14] and those of [11]. Thus these two classes of McVittie spacetimes both have this characteristic black hole property: particles and photons can be confined to a spatially compact region of spacetime by means of the spacetime geometry.

Before proceeding, we give a brief but detailed summary of the results of this paper, and note their relation to previous work.

## A. Summary of Section II

In Section II, we write down the relevant geodesic equations and state relevant dynamical systems results. We point out useful bounds relating ˙ t and ˙ along individual r geodesics - see (20).

## B. Summary of Section III

In Section III, we derive and solve an equation that must hold for the Hubble function H t ( ) of a McVittie spacetime that admits a CPO. We briefly discuss global properties of the associated spacetime: a fuller discussion is given in Appendix A. As we see in (1), McVittie spacetimes are determined by a parameter m and a function H . We show in Section III that there is (only) a 2-parameter family of McVittie spacetimes that admit a CPO. Thus this is a 'no-go' result: CPOs are absent from nearly all McVittie spacetimes. We note that the same conclusion holds in relation to circular timelike orbits.

## C. Summary of Section IV

In Section IV we define a class of McVittie spacetimes that we will refer to as expanding McVittie spacetimes with a big bang background . These are the focus of the remainder of the paper: this class appears to us to be the class of McVittie spacetimes of most physical interest. See Definition 1. We prove that in this class, all causal geodesics (timelike and null, radial and non-radial) originate at the singular boundary r = 2 m at finite affine distance (proper time) in the past. This generalises previous results relating to radial null geodesics [11, 13, 14].

## D. Summary of Section V

Here, we deal with the future evolution of causal geodesics. We define two subclasses of expanding McVittie spacetimes with a big bang background that are distinguished by the (future) asymptotic value of the Hubble function H t ( ). For Class 1, H t ( ) → H &gt; 0 0 as t → + ∞ and for Class 2, H t ( ) → 0. See Definitions 2 and 3. We derive rigorously and generalise results that have been presented previously either in a heuristic manner, or for special cases. In particular, we prove that the global structure (in the sense of a Penrose-Carter conformal diagram) obtained in [14] for a special case - i.e. H t ( ) specified -holds generally for Class 1 (we note that several of the results of [14] were proven in general). Included under this heading are existence proofs relating to important radial null geodesics that delineate the global structure (see Propositions 4 - 7). The crucial role of a non-zero value of H 0 was first identified in [13]: this has far-reaching consequences in that it changes radically the interpretation of the corresponding McVittie spacetime, as discussed above. In Section V, our intention, in part, is to put the physical insights of this paper on a sounder mathematical footing. We note also a generalisation and some corrections (Proposition 7 (b) - noted in [14]) to the results of [13]. For example, the linear equation of state used in Appendix A of that paper is not required. In addition, we generalise some previous results for the case H 0 = 0 [11].

## E. Summary of Section VI

The main results of the paper are given here. Using various techniques from dynamical systems, we prove the existence of bound photon and particle orbits in both Class 1 and Class 2 spacetimes. That is, we prove the existence of future-complete non-radial timelike and null geodesics γ with the property that r | γ is finite along the whole history of the geodesic. Taking s to be the parameter (affine parameter or proper time) along the geodesic, we show that lim s → ∞ + r s ( ) = r c , where r c is the radius of a circular orbit in the corresponding Schwarzschild-de Sitter (Class 1) or Schwarzschild (Class 2) spacetime. For

photon orbits, this forces r c = 3 m and for particle orbits, r c corresponds to a stable circular orbit. See Propositions 10 and 12-14.

We use units with G = c = 1, and we use a glyph[square] to indicate the end of a proof.

## II. GENERAL PROPERTIES OF CAUSAL GEODESICS

We begin with the line element (1) and define

$$\chi ( t, r ) = f - r ^ { 2 } H ^ { 2 }, \quad \alpha ( t, r ) = r f ^ { - 1 / 2 } H, \quad \lambda ( r ) = f ^ { - 1 }. \quad \text{ of } ( 1 7 \\ \psi \quad.. \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ddots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

Recall that f ( r ) = 1 -2 m/r , and note that χλ + α 2 = 1. The geodesic equations of the spacetime may be written in the following form:

$$\ddot { r } = r f ^ { 1 / 2 } H ^ { \prime } \dot { t } ^ { 2 } + ( 1 - \frac { 3 m } { r } ) \frac { \ell ^ { 2 } } { r ^ { 3 } } + \epsilon ( \frac { m } { r ^ { 2 } } - r H ^ { 2 } ), \quad ( 1 3 ) \quad \text{so th} i$$

$$\ddot { t } = - ( 1 - \frac { 3 m } { r } ) f ^ { - 1 / 2 } H \dot { t } ^ { 2 } - \frac { 2 m } { r ^ { 2 } } f ^ { - 1 } \dot { t } \dot { r } - \epsilon f ^ { - 1 / 2 } H ( 1 4 )$$

where glyph[lscript] is the conserved angular momentum of the geodesic, glyph[epsilon1] = 0 for null geodesics and glyph[epsilon1] = -1 for timelike geodesics. The overdot represents the derivative with respect to the parameter along the geodesic: an affine parameter for null geodesics and proper time for timelike geodesics. Throughout the remainder of this paper, we use s to represent this parameter. We also have the first integral of (13)-(14):

$$- \, \chi \dot { t } ^ { 2 } - 2 \alpha \dot { t } \dot { r } + f ^ { - 1 } \dot { r } ^ { 2 } + \frac { \ell ^ { 2 } } { r ^ { 2 } } = \epsilon. \quad \ \ ( 1 5 )$$

glyph[negationslash]

It follows that ˙ t = 0 everywhere on a causal geodesic. Then by a choice of parameter orientation, we have ˙ t &gt; 0 everywhere along all causal geodesics.

We consider the geodesic equations on the region Ω = { ( t, r ) : t &gt; 0 , r &gt; 2 m } . This region has boundaries ∂ Ω 0 = { (0 , r ) : r ≥ 2 m } and ∂ Ω 2 m = { ( t, 2 m ) : t &gt; 0 } . It will be convenient to consider the geodesic equations as a first order dynamical system. To this end, we define glyph[vector] x = ( r, t, r, t ˙ ˙ ) and write the geodesic equations in the form

$$\frac { d \vec { x } } { d s } = \vec { F } ( \vec { x } ), \quad \quad ( 1 6 ) \quad w i$$

where F 1 ( glyph[vector] x ) = ˙ = r x 3 , F 2 ( glyph[vector] x ) = ˙ = t x 4 and F 3 ( glyph[vector] x , F ) 4 ( glyph[vector] x ) can be read-off the right hand sides of (13) and (14). We note that glyph[vector] F ∈ C 1 ( E ), where E = Ω × R 2 . In this context, (15) is a zero-order constraint that is propagated by the equations. That is, if (15) holds for glyph[vector] x = glyph[vector] x 0 , then this equation is valid everywhere along the solution glyph[vector] x of (16) with the initial condition glyph[vector] x (0) = glyph[vector] x 0 .

Local existence of C 1 solutions of (16) follows from standard results of dynamical systems (see e.g. [16]). For any glyph[vector] x 0 ∈ E , there is a maximal interval of existence ( s α , s ω ) = ( s α ( glyph[vector] x 0 ) , s ω ( glyph[vector] x 0 )) of the initial value problem

$$\frac { d \vec { x } } { d s } = \vec { F } ( \vec { x } ), \ \vec { x } ( 0 ) = \vec { x } _ { 0 }. \quad \quad ( 1 7 ) \quad \text{Th}$$

The intervals ( s α , 0] and [0 , s ω ) are respectively the leftmaximal and right-maximal intervals of existence. The following theorem plays an important role in the next section:

glyph[negationslash]

Theorem 1 (Perko,[16]) Let E be an open subset of R n containing glyph[vector] x 0 , let glyph[vector] F ∈ C 1 ( E ) and let ( s α , 0] be the leftmaximal interval of existence of the IVP (17). If s α &gt; -∞ , then given any compact set K ⊂ E , there exists s ∈ ( s α , 0) such that glyph[vector] x s ( ) ∈ K .

Roughly speaking, this theorem tells us that solutions of (17) continue to exist while glyph[vector] x remains bounded. We will exploit this theorem to prove extension results for causal geodesics.

By completing a square, we can write (15) in the form

$$\dot { t } ^ { 2 } - ( \alpha \dot { t } - \lambda \dot { r } ) ^ { 2 } = \lambda ( \ell ^ { 2 } r ^ { - 2 } + | \epsilon | ), \quad \ \ ( 1 8 )$$

glyph[negationslash]

so that for non-radial geodesics ( glyph[lscript] = 0),

$$( ( 1 - \alpha ) \dot { t } + \lambda \dot { r } ) ( ( 1 + \alpha ) \dot { t } - \lambda \dot { r } ) > 0. \quad ( 1 9 )$$

Noting that λ and ˙ t are both positive, and that α &gt; 0 in an expanding ( H &gt; 0) McVittie spacetime, we obtain the following:

Lemma 1 At every point on a non-radial causal geodesic in an expanding McVittie spacetime,

$$( \alpha - 1 ) \dot { t } < \lambda \dot { r } < ( \alpha + 1 ) \dot { t }. \text{ \quad \ \ } ( 2 0 )$$

glyph[square]

## III. CIRCULAR PHOTON ORBITS

glyph[negationslash]

For circular photon orbits, we have r = a = constant along the null geodesic. We take a &gt; 2 m . Imposing this condition in (13) and (15), eliminating glyph[lscript] and noting that ˙ t = 0, we obtain the following ODE:

$$H ^ { \prime } ( t ) = A ( H ^ { 2 } - H _ { 0 } ^ { 2 } ), \text{ \quad \ \ } ( 2 1 )$$

where

$$A = \frac { ( 1 - \frac { 3 m } { a } ) } { \sqrt { 1 - \frac { 2 m } { a } } }, \ \ H _ { 0 } = a ^ { - 1 } \sqrt { 1 - \frac { 2 m } { a } }. \quad ( 2 2 )$$

We note that if the circular photon orbit has radius a = 3 m , then A = 0 and we are in Schwarzschild-de Sitter spacetime. So for a non-trivial solution, we exclude a = 3 m . We note then that mH &lt; 0 1 3 √ 3 (equality holds for a = 3 m ). The only remaining condition of (13)-(15) is

$$\chi ( t, a ) \dot { t } ^ { 2 } = \frac { \ell ^ { 2 } } { a ^ { 2 } }.$$

This yields the following result.

Proposition 1 The McVittie spacetime with line element (1) admits circular photon orbits if and only if the Hubble function H t ( ) satisfies (21), with the constants A,H 0 determined by the mass parameter m and the orbital radius a via (22), and χ t, a ( ) &gt; 0 for all t along the orbit. glyph[square]

The function χ , which is invariantly defined by χ = g ab ∇ a r ∇ b r , plays an important role in spherical symmetry: this leads to an interesting corollary to Proposition 1. As is well known, χ is proportional to the product of the expansions of the ingoing and outgoing radial null geodesics of a spherically symmetric spacetime. Hence the surface (or surfaces) χ = 0 corresponds to a horizon, where (at least) one of the null expansions vanishes. The regular (or untrapped) region of the spacetime corresponds to χ &gt; 0, where one null expansion is positive and one is negative, while χ &lt; 0 corresponds to either a trapped region (two negative null expansions) or an antitrapped region (two positive null expansions). We immediately have the following, which is a particular case of a more general result:

Corollary 1 A circular photon orbit of a McVittie spacetime is confined to the regular region of the spacetime. glyph[square]

We can also see from (15) that any turning points (periapsis or apasis, whereat ˙ r = 0) of a photon orbit in a McVittie spacetime must lie in the regular region of the spacetime. The same conclusion holds for particle orbits (timelike geodesics), where (15) holds but with -1 on the right hand side. In fact it is readily seen that Corollary 1 holds much more generally in spherically symmetric spacetimes (this is the more general result referred to above). Using double null coordinates ( u, v ), both taken to increase into the future, the null expansion θ u along the future pointing null direction ∂ ∂u satisfies θ u = 1 κ 2 u ∂r ∂u for some positive function κ 2 u . A corresponding statement holds with u replaced by v . Then along any future pointing null geodesic,

$$\dot { r } = \kappa _ { u } ^ { 2 } \theta _ { u } \dot { u } + \kappa _ { v } ^ { 2 } \theta _ { v } \dot { v }.$$

With u, v increasing into the future, so that ˙ u &gt; 0 and ˙ v &gt; 0, we see that ˙ r can vanish only when either the null expansions θ u , θ v have opposite sign - that is, when the geodesic is in a regular region of the spacetime - or when both θ u and θ v vanish. This latter situation is non-generic: it occurs for example at the bifurcation 2sphere of the extended Schwarzschild (Kruskal-Szekeres) spacetime. The former conclusion should hold in general.

Corollary 1 is almost enough to rule out CPO's in any McVittie spacetime for which the isotropic background has a big bang - i.e. S t ( 0 ) = 0 at some time t 0 in the past. From (2), we see that H t ( ) will typically diverge as t → t + 0 . But then χ t, a ( ) = 1 -2 m a -a H 2 2 becomes negative at early times. It is very difficult to see how the CPOcan exit the trapped region, or terminate in the past at a finite positive value of t , without there being some

@ D

FIG. 1. The graph of (i) the Hubble function H t ( ) and (ii) the scale factor S t ( ) for the McVittie spacetime admitting a circular photon orbit at radius a &lt; 3 m . Here m = 1 and a = 2 5, giving . A = -0 447214 . &lt; 0 and H 0 = 0 178885. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0044v1-figure-aee167cb-d274-4b31-999a-9674ea581ac6.jpg)

serious pathology of the spacetime. Thus the existence of a CPO is a very strong restriction on McVittie spacetimes with big bang backgrounds.

We turn now to determining which line elements of the form (1) are admitted by Proposition 1. This amounts to solving (21), and ensuring that the resulting χ t, a ( ) is positive. There is a unique non-trivial solution (where non-trivial means that H is not constant: recall that this corresponds to Schwarzschild-de Sitter spacetime):

$$H ( t ) = - H _ { 0 } \tanh ( A H _ { 0 } t ), \ \ t \in \mathbb { R }. \quad \ ( 2 5 )$$

We have used time translation freedom to set t = 0 at the zero of H . The corresponding scale factor is

$$S ( t ) = ( \cosh ( A H _ { 0 } t ) ) ^ { - 1 / A }, \text{ \quad \ \ } ( 2 6 )$$

and we find

$$\chi ( t, a ) = a ^ { 2 } H _ { 0 } ^ { 2 } / \cosh ^ { 2 } ( A H _ { 0 } t ). \quad \quad ( 2 7 )$$

There are two inequivalent spacetimes, depending on the sign of A . We see from (22) that A is positive (respectively negative) if the CPO radius a is greater than (respectively less than) 3 m . The expansion histories of the background universes corresponding to the two choices are shown in Figures 1 and 2, and their global structure is analysed in the Appendix A.

We find it curious that (21) has a solution which does satisfy the Λ-CDM conditions of [14] (see also Definition 1 below). This solution is H t ( ) = -H 0 coth( AH t 0 ). However this yields χ t, a ( ) = -a H / 2 2 0 sinh 2 ( AH t 0 ), and so the second condition (23) for a CPO is violated.

@ D

FIG. 2. The graph of (i) the Hubble function H t ( ) and (ii) the scale factor S t ( ) for the McVittie spacetime admitting a circular photon orbit at radius a &gt; 3 m . Here m = 1 and a = 3 5,giving . A = 0 218218 . &gt; 0 and H 0 = 0 187044. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0044v1-figure-c2508ace-f418-4b65-9772-8f050911510f.jpg)

Wesummarise the main result of this section as follows.

Proposition 2 The McVittie spacetime with line element (1) admits a circular photon orbit if and only if the Hubble function is given by (25), with A,H 0 determined by the mass parameter m and photon orbital radius a as in (22). glyph[square]

Comment 2 It is straightforward to derive and solve the equation corresponding to (21) in the case of circular particle orbits. This yields a different McVittie spacetime to that of Proposition 2 - but again, only a 2-parameter family of spacetimes arises. We will not pursue this further: the main point is that circular orbits almost never exist in McVittie spacetimes.

## IV. PAST EVOLUTION

We begin the dicsussion by considering the past evolution of causal geodesics in the class of McVittie spacetimes of interest. These correspond to McVittie spacetimes for which the isotropic background is an expanding cosmological model - an expanding isotropic spacetime with a big bang at a finite time in the past.

Definition 1 The spacetime with line element (1) is said to be an expanding McVittie spacetime with a big bang background if the following conditions on the Hubble function hold:

$$( i ) \ H \in C ^ { 2 } ( ( 0, + \infty ), \mathbb { R } _ { + } ) ;$$

$$( i i ) \, \lim _ { t \rightarrow 0 ^ { + } } \int _ { t } ^ { t _ { 0 } } H ( u ) d u = + \infty \text{ for all } t _ { 0 } > 0 ;$$

$$( i i i ) \ H ^ { \prime } ( t ) < 0 \ f o r \ a l l \ t > 0.$$

The first condition here includes a technical differentiability condition and the condition for expansion: H t ( ) &gt; 0 for all t &gt; 0. The second incorporates the big bang condition. This is equivalent to the existence of some t 0 ∈ R -which we set to zero by a translation -for which the scale factor S satisfies S t ( 0 ) = 0. The third condition is equivalent to the weak energy condition ρ + p &gt; 0 in both the background and the McVittie spacetime - see (7) and (8).

In this section, we establish the fact that ∂ Ω 2 m = { ( t, r ) : t &gt; 0 , r = 2 m } -which, recall, is a curvature singularity - forms the past boundary of the spacetime. That is, all causal geodesics originate in the past at ∂ Ω 2 m at finite affine distance (null geodesics) or finite proper time (timelike geodesics). All such geodesics extend back to ∂ Ω 2 m at a positive value of t : the background big bang surface ∂ Ω 0 = { ( t, r ) : t = 0 , r ≥ 2 m } is cut off by ∂ Ω 2 m . These statements are non-trivial and so proofs are given below. We note that it is straightforward to see that r = 2 m forms a spacelike portion of the past boundary of the spacetime [13]: what is not obvious is that all causal geodesics originate here. This is stated formally as Proposition 3 below.

The radial null geodesics of (1) are of particular importance in determining the global structure of the spacetime. These are the geodesics satisfying (13)-(15) with glyph[lscript] = glyph[epsilon1] = 0. Noting that t is a global time coordinate of (1) that increases into the future by a choice made in Section I above, we have ˙ t &gt; 0 everywhere along a causal geodesic, and so from (15) we can write down

$$\frac { d r } { d t } = H r ( 1 - \frac { 2 m } { r } ) ^ { 1 / 2 } \pm ( 1 - \frac { 2 m } { r } ) \quad \ ( 2 8 )$$

for radial null geodesics (RNG). There are two families: outgoing, corresponding to the upper sign, and ingoing, corresponding to the lower. We introduce the radial coordinate z as defined in [11]:

$$z = \sqrt { 1 - \frac { 2 m } { r } }. \quad \quad \quad ( 2 9 )$$

Notice then that z ∈ [0 , 1) with z = 0 ⇔ r = 2 m and z → 1 as r →∞ . In terms of z , Ω = { ( t, z ) : t ≥ 0 , z ∈ (0 1) , } and the RNGs satisfy

$$\frac { d z } { d t } = \frac { H } { 2 } ( 1 - z ^ { 2 } ) \pm \frac { z } { 4 m } ( 1 - z ^ { 2 } ) ^ { 2 }. \quad \ ( 3 0 )$$

Our first step is to establish the following: given any point P : ( t, z ) = ( t 1 , z 1 ) ∈ Ω (so that t 1 &gt; 0 and z 1 ∈ (0 1)), , the ingoing and outgoing radial null geodesics (IRNGs and ORNGs) through P both originate on the surface ∂ Ω 2 m = { r = 2 m, 0 &lt; t &lt; ∞} =

{ z = 0 0 , &lt; t &lt; ∞} . We then show that the same holds for all causal geodesics passing through P , and thereby show that the set ∂ Ω 0 does not form part of the past boundary of the spacetime.

Lemma 2 In an expanding McVittie spacetime with a big bang background, all outgoing radial null geodesics of the spacetime originate at ∂ Ω 2 m at finite affine distance in the past. That is, for each ORNG γ , there exists s 0 &gt; -∞ such that lim s → s + 0 ( ( t s , r ) ( s )) | γ = ( t 0 , 2 m ) for some t 0 &gt; 0 .

Proof: Consider the ORNG through ( t, z ) = ( t 1 , z 1 ) with t 1 &gt; 0 and z 1 ∈ (0 , 1). This satisfies

$$\frac { d z } { d t } = \frac { H } { 2 } ( 1 - z ^ { 2 } ) + \frac { z } { 4 m } ( 1 - z ^ { 2 } ) ^ { 2 } > 0 \quad ( 3 1 ) \quad \gamma ( 0 ) \quad \ s o i n$$

for t &gt; 0 0 , ≤ z &lt; 1. So for t &lt; t 1 , we have z ≤ z 1 , and hence

$$\frac { d z } { d t } \geq \frac { H } { 2 } ( 1 - z _ { 1 } ^ { 2 } ) + \lambda _ { 1 } z, \quad \quad \ ( 3 2 ) \quad i s \ \tau \\ \ t h e$$

where λ 1 = (1 -z 2 1 ) 2 / 4 m &gt; 0, and we recall from part (i) of Definition 1 that H &gt; 0. Integrating this linear differential inequality from t to t 1 yields

$$z ( t ) & \leq z _ { 1 } e ^ { \lambda _ { 1 } ( t - t _ { 1 } ) } - \frac { ( 1 - z _ { 1 } ) ^ { 2 } } { 2 } \int _ { t } ^ { t _ { 1 } } H ( u ) e ^ { \lambda _ { 1 } ( t - u ) } d u \\ & \leq z _ { 1 } e ^ { \lambda _ { 1 } ( t - t _ { 1 } ) } - \frac { ( 1 - z _ { 1 } ) ^ { 2 } } { 2 } e ^ { - \lambda _ { 1 } t _ { 1 } } \int _ { t } ^ { t _ { 1 } } H ( u ) d u. \quad \ \text{T} \quad \mathfrak { d }$$

These inequalities hold for t &lt; t 1 . From part (ii) of Definition 1, we see that z t ( ) must reach 0 (and so r reaches 2 m ) at some positive value t 0 of t . From (14) with glyph[lscript] = 0 and part (iii) of Definition 1, we see that ¨ r &lt; 0 along the geodesic. This guarantees that r reaches 2 m at a finite value s 0 of the affine parameter s . glyph[square]

Lemma 3 In an expanding McVittie spacetime with a big bang background, all ingoing radial null geodesics of the spacetime originate at ∂ Ω 2 m at finite affine distance in the past. That is, for each IRNG γ , there exists s 0 &gt; -∞ such that lim s → s + 0 ( ( t s , r ) ( s )) | γ = ( t 0 , 2 m ) for some t 0 &gt; 0 .

Proof: For IRNGs, we have

$$\frac { d z } { d t } = \frac { H } { 2 } ( 1 - z ^ { 2 } ) - \frac { z } { 4 m } ( 1 - z ^ { 2 } ) ^ { 2 }. \quad \ \ ( 3 3 ) \quad \ \ g r o \iota ^ { \nu \ u } _ { \ } t h a t$$

We note that 0 ≤ z (1 -z 2 ) ≤ 2 / (3 √ 3) for all z ∈ [0 , 1). Then defining H c = 1 3 √ 3 m , we find

$$\frac { d z } { d t } \geq \frac { 1 } { 2 } ( H - H _ { c } ) ( 1 - z ^ { 2 } ), \quad t > 0, \ \ 0 \leq z < 1. \ ( 3 4 ) \quad \text{$groun$}$$

Integrating from ( t, z ) to ( t 1 , z 1 ) where t 1 &gt; t eventually yields

$$z ( t ) \leq \frac { ( 1 + z _ { 1 } ) - ( 1 - z _ { 1 } ) e ^ { I ( t ) } } { ( 1 + z _ { 1 } ) + ( 1 - z _ { 1 } ) e ^ { I ( t ) } }, \quad t < t _ { 1 } \quad ( 3 5 )$$

where

$$I ( t ) = \int _ { t } ^ { t _ { 1 } } ( H ( u ) - H _ { c } ) d u.$$

Divergence of I in the limit t → 0 + (which comes from part (ii) of Definition 1) shows that z reaches 0 at a positive value of t . The proof that this occurs at a finite value of the affine parameter is identical to the corresponding proof in Lemma 2. glyph[square]

The following results establish that the results of these two lemmas apply to all causal geodesics. We first establish a somewhat obvious 'radial confinement' result.

Lemma 4 Let p ∈ Ω and let γ be a causal geodesic with γ (0) = p . Then the past branch of γ -that is, the set of points

$$\gamma ^ { - } = \{ \gamma ( s ) \, \colon s _ { \alpha } < s < 0 \} \, \text{ \quad \ \ } ( 3 7 )$$

is contained in the interior of the region of Ω bounded by the point p , by the paths of the past-directed radial null geodesics through p and by the boundary ∂ Ω 2 m .

Proof: Noting that ˙ t &gt; 0 along γ , we find from (20) that

$$\lambda ^ { - 1 } ( \alpha - 1 ) < \frac { d r } { d t } < \lambda ^ { - 1 } ( \alpha + 1 ). \quad \ \ ( 3 8 )$$

The lower and upper bounds correspond to the uniquely defined value at each point of Ω of the slopes of the ingoing and outgoing radial null geodesics. Thus at any given point of Ω, the geodesic γ crosses the ingoing (respectively outgoing) radial null geodesic from below (respectively above). It follows that to the past of p (i.e. for s &lt; 0), γ remains below (respectively above) the ingoing (respectively outgoing) radial null geodesic through p . The conclusion follows. glyph[square]

Comment 3 It is tempting to conclude on the basis of this lemma that any causal γ must extend back to r = 2 m . However it remains to prove that γ extends sufficiently far into the past in order that this happens. We now prove this extension result.

Proposition 3 Let γ be a future-pointing causal geodesic of an expanding McVittie spacetime with a big bang background. Then there exists t α &gt; 0 and s α &gt; -∞ such that lim s → s + α ( ( t s , r ) ( s )) | γ = ( t α , 2 m ) where s is an affine parameter (respectively proper time) along a null (respectively timelike) geodesic. That is, all causal geodesics of an expanding McVittie spacetime with a big bang background originate at ∂ Ω 2 m .

Proof: Let γ (0) = p ∈ Ω, and as before, let ( s α , 0] be the left maximal interval of existence for γ . Applying the radial confinement result of the previous lemma and recalling that ˙ t &gt; 0 along γ , we have

$$t _ { i } ( p ) < t ( s ) < t ( 0 ), \ \ s \in ( s _ { \alpha }, 0 ), \quad \ \ ( 3 9 )$$

where t i ( p ) is the value of t at which the ingoing radial null geodesic through p meets ∂ Ω 2 m . Likewise,

$$r ( s ) < r _ { + } ( p ), \ \ s \in ( s _ { \alpha }, 0 ), \quad \ \ ( 4 0 ) \quad \ \stackrel { \ll } { F }$$

where r + ( p ) is the maximum of the value of r along the past branches of the radial null geodesics through p .

From Lemma 1, and recalling that λ = f -1 and that ˙ t &gt; 0, we have

$$- \, \frac { 2 m } { r } f ^ { - 1 } \dot { t } \dot { r } > - \frac { 2 m } { r } ( \alpha + 1 ) \dot { t } ^ { 2 }. \quad \ \ ( 4 1 ) \quad I t \ f i$$

Then from (14), we can write

$$\ddot { t } > - \left ( ( 1 - \frac { m } { r } ) f ^ { - 1 / 2 } H + \frac { 2 m } { r ^ { 2 } } \right ) \dot { t } ^ { 2 }, \quad ( 4 2 ) \quad \text{Usir}$$

where we note that the coefficient of ˙ t 2 on the right hand side is strictly negative on Ω.

For p ∈ Ω and glyph[epsilon1] &gt; 0, define the compact sets

$$K _ { p, \epsilon } & = \{ ( t, r ) \, \colon t _ { i } ( p ) - \epsilon \leq t \leq t ( 0 ), 2 m + \epsilon \leq r \leq r _ { + } ( p ) \}. \quad \text{and } E \\ \underset { \longrightarrow } \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,.$$

Then there exists a positive constant C p,glyph[epsilon1] such that

$$( 1 - \frac { m } { r } ) f ^ { - 1 / 2 } H + \frac { 2 m } { r ^ { 2 } } \leq C _ { p, \epsilon }, \quad \text{for all $(t,r)\in K_{p,\epsilon}$}. \quad \text{may $\mathfrak{1}$} \\ \mathfrak { n } _ { \dots } \dots \dots \dots \lambda \dots \lambda \dots \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \lambda \dots \varphi$}$$

Now suppose that the past branch of γ is confined to K p,glyph[epsilon1] - that is, that ( ( t s , r ) ( s )) ∈ K p,glyph[epsilon1] for all s ∈ ( s α , 0). Then along γ ,

$$\ddot { t } \geq - C _ { p, \epsilon } \dot { t } ^ { 2 }, \ \ s \in ( s _ { \alpha }, 0 ). \quad \quad ( 4 5 )$$

Integrating yields

$$\dot { t } ( s ) < \dot { t } ( 0 ) \exp [ ( t ( 0 ) - t ( s ) ) C _ { p, \epsilon } ], \ \ s \in ( s _ { \alpha }, 0 ), \ \ ( 4 6 ) \quad \Box$$

and so

$$0 \leq \dot { t } ( s ) \leq \dot { t } ( 0 ) \exp [ t ( 0 ) C _ { p, \epsilon } ], \ \ s \in ( s _ { \alpha }, 0 ). \quad ( 4 7 )$$

Using (20), this yields finite upper and lower bounds for ˙: r

$$C _ { p, \epsilon } ^ { - } \leq \dot { r } ( s ) \leq C _ { p, \epsilon } ^ { + }, \ \ s \in ( s _ { \alpha }, 0 ), \quad \ ( 4 8 ) \quad \ r a d i$$

where the bounding constants C ± p,glyph[epsilon1] depend on the parameter glyph[epsilon1] , the initial point p and the initial value and derivative of t .

From (47) and (48), it follows that the corresponding solution glyph[vector] x s ( ) of the dynamical system (16) is confined to the compact subset K of E = Ω × R 2 , where

$$K = K _ { p, \epsilon } \times [ C _ { p, \epsilon } ^ { - }, C _ { p, \epsilon } ^ { + } ] \times [ 0, \dot { t } ( 0 ) \exp [ t ( 0 ) C _ { p, \epsilon } ]. \ \ ( 4 9 ) \quad \text{result}$$

This contradicts Theorem 1, and so γ must exit K p,glyph[epsilon1] . The radial confinement result indicates that γ must exit the set via the lower boundary r = 2 m + . That is, for glyph[epsilon1] every glyph[epsilon1] &gt; 0, there exists s 0 ∈ ( s α , 0) so that r s ( 0 ) &lt; 2 m + . It follows that there exists glyph[epsilon1] s 2 m ∈ -∞ [ , 0) such that lim s → s + 2 m r s ( ) = 2 m . We now show that s 2 m &gt; -∞ , completing the proof. (Note that s 2 m = -∞ corresponds to s α = -∞ .)

Recall that f = 1 -2 m/r and that H t ( ) &gt; 0 with H t ′ ( ) &lt; 0. It follows that H t s ( ( )) ≤ H t ( i ( p )) for all s ∈ ( s α , 0). Recall that we may write (15) in the form

$$\dot { t } ^ { 2 } - ( \alpha \dot { t } - \lambda \dot { r } ) ^ { 2 } = f ^ { - 1 } ( | \epsilon | + \frac { \ell ^ { 2 } } { r ^ { 2 } } ). \quad \ \ ( 5 0 )$$

It follows that

$$\lim _ { s \to s _ { 2 m } ^ { + } } \dot { t } ( s ) = + \infty. \quad \quad \quad ( 5 1 )$$

Using the left-hand inequality of (20) in (14), we can write

$$\ddot { t } < - ( 1 - \frac { m } { r } ) f ^ { - 1 / 2 } H \dot { t } ^ { 2 } + \frac { 2 m } { r ^ { 2 } } \dot { t } ^ { 2 } + f ^ { - 1 / 2 } H. \quad ( 5 2 )$$

On the right hand side, the terms (1 -m r ) f -1 / 2 H , 2 m r 2 and H all have finite positive limits as the geodesic approaches r = 2 m . As ˙ and t f -1 / 2 both diverge to + ∞ in the limit, it follows that the first term on the right hand side (rhs) of (52) dominates the other two, so that we may write

$$\stackrel { \mu } { ( 4 4 ) } _ { K _ { p, \epsilon } } \quad \stackrel { \dots } { t } < [ \text{rhs of } ( 5 2 ) ] \sim - ( 1 - \frac { m } { r } ) f ^ { - 1 / 2 } \stackrel { \dots } { t } ^ { 2 }, \ \ s \to s _ { 2 m } ^ { + }. \ \ ( 5 3 )$$

Thus there exists s ∗ ∈ ( s 2 m , 0) such that

$$\ddot { t } ( s ) < 0, \ \ s < s _ { * }. \quad \quad \quad ( 5 4 )$$

It follows by integrating twice that s 2 m &gt; -∞ : the geodesic reaches r = 2 m in finite affine (or proper) time. glyph[square]

## V. FUTURE EVOLUTION

The main aim of this section is to build a comprehensive picture of the future evolution of radial null geodesics in expanding McVittie spacetimes with a big bang. We will also, where possible, draw conclusions about nonradial geodesics. The results we obtain both generalise and (attempt to) clarify previous results. The intention is to clarify the connection between the different global structures derived in previous work and the features of the corresponding background spacetimes (as encoded in the Hubble function H ). The main results of the paper, which relate to bound photon and particle orbits, are presented in the following section and many of the results of this section are not required for those results. However this section is not wholly an aside: the reader interested in the results of Section VI should review Definitions 2 and 3, Lemmas 5 and 6 and the paragraphs between Propositions 4 and 5.

As is evident from the contrasting results of [11] and [13], the asymptotic value as t → + ∞ of the Hubble

function H t ( ) has a significant influence on the global structure of the spacetime. This is reflected in the results below, and to allow us to present those results in a clear manner, we define two subclasses of expanding McVittie spacetimes with a big bang.

The first class comprises McVittie spacetimes where the scale factor of the isotropic background is that of an expanding, k = 0, Λ-CDM isotropic universe. For the purposes of this paper, the defining properties are these:

Definition 2 The spacetime with line element (1) is a Class 1 McVittie spacetime if the Hubble function H satisfies (i)-(iii) of Definition 1, plus the following conditions:

$$( i v ) \lim _ { \substack { t \to \infty \\ H _ { 0 } > 0. } } ( H ( t ), H ^ { \prime } ( t ), H ^ { \prime \prime } ( t ) ) \ = \ ( H _ { 0 }, 0, 0 ) \ \text{where}$$

- (v) The background density and pressure satisfy (9) and (10) with Λ = 3 H 2 0 , and there is an equation of state p 0 = g ρ ( 0 ) satisfying g ∈ C 1 [0 , ∞ ) with the condition

$$\kappa \coloneqq \frac { 3 } { 2 } ( 1 + g ^ { \prime } ( 0 ) ) > 0 \quad \quad ( 5 5 ) \quad \text{ with } \quad \text{$w$} _ { \mathcal { L } \,. }$$

on the sound speed at zero density.

Comment 4 We note that the technical condition on the sound speed corresponds to the existence and positivity of the limit lim ρ 0 → 0 1 + p 0 ρ 0 and so expresses a physically motivated energy condition, as well as a differentiability condition on g .

The second class we consider corresponds to isotropic backgrounds that share features with spatially flat Robertson-Walker universes with equation of state p = constant × ρ and zero cosmological constant. The defining properties are these:

Definition 3 The spacetime with line element (1) is a Class 2 McVittie spacetime if the Hubble function H satisfies (i)-(iii) of Definition 1, plus the following conditions:

$$( i v ) \lim _ { t \rightarrow \infty } ( H ( t ), H ^ { \prime } ( t ), H ^ { \prime \prime } ( t ) ) = ( 0, 0, 0 ).$$

- (v) The background density and pressure satisfy (9) and (10) with Λ = 0 , and there is an equation of state p 0 = g ρ ( 0 ) satisfying g ∈ C 2 [0 , ∞ ) with the condition

$$\kappa \coloneqq \frac { 3 } { 2 } ( 1 + g ^ { \prime } ( 0 ) ) > 0 \quad \quad ( 5 6 ) \quad \text{$\text{Cl}$} \\ \Omega _ { R }$$

on the sound speed at zero density.

↦

Comment 5 The additional differentiability requirement on the equation of state function ρ 0 → g ρ ( 0 ) is a technical condition required for some of the proofs of Section VI below.

Our first result proves the existence of outgoing photon orbits that extend to infinity. This demonstrates the existence of an asymptotic region of every expanding McVittie spacetime with a big bang background where light rays extend to arbitrarily large radii.

Proposition 4 Let γ be an outgoing radial null geodesic of an expanding McVittie spacetime with a big bang background. Then γ is future complete, and lim s →∞ ( t, r ) | γ = (+ ∞ , + ∞ ) .

Proof: Along an outgoing radial null geodesic, we have

$$\frac { d z } { d t } = \frac { H } { 2 } ( 1 - z ^ { 2 } ) + \frac { z } { 4 m } ( 1 - z ^ { 2 } ) ^ { 2 }, \quad \ \ ( 5 7 )$$

which is positive so that z increases with t . Since H &gt; 0 (part (i) of Definition 1), we see that

$$\Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \PhI _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ { 0 } \quad \Phi _ }$$

where z 0 = z t ( 0 ) and t 0 &gt; 0 is arbitrary. Integrating shows that z → 1 (i.e. r → + ∞ ) as t → + ∞ . From (13) and part (iii) of Definition 1, we see that ¨ r &lt; 0 along the geodesic, and so the affine parameter s must extend to + ∞ as t, r do. glyph[square]

Next, we show that in Class 2 McVittie spacetimes, all photon and particle orbits that originate outside a certain radius extend to infinity and are future complete. This and some subsequent proofs require some details of the horizon structure of McVittie spacetimes [11, 13, 14].

Recall that the horizon of the McVittie spacetime with line element (1) is the set of points with χ = g ab ∇ a r ∇ b r = 0, the regular region is the set of points with χ &gt; 0 and the anti-trapped region is the set of points with χ &lt; 0. It is straightforward to show that the expanding condition H &gt; 0 leads to the region χ &lt; 0 being anti-trapped rather than trapped. We will denote the regular region by Ω R , the anti-trapped region by Ω A and the horizon by H . We note that Ω is the disjoint union of Ω R , Ω A and H and that

$$\Omega _ { R } = \{ ( t, r ) \in \Omega \, \colon H ( t ) < \frac { 1 } { r } ( 1 - \frac { 2 m } { r } ) ^ { 1 / 2 } \}, \quad ( 5 9 )$$

$$\Omega _ { A } = \{ ( t, r ) \in \Omega \, \colon H ( t ) > \frac { \cdot } { r } ( 1 - \frac { \cdot } { r } ) ^ { 1 / 2 } \}, \quad ( 6 0 )$$

$$\mathcal { H } = \{ ( t, r ) \in \Omega \colon H ( t ) = \frac { \cdot } { r } ( 1 - \frac { \cdot } { r } ) ^ { 1 / 2 } \}. \quad ( 6 1 )$$

There are crucial differences between the horizons of Class 1 and Class 2 McVittie spacetimes. First, we note that in a Class 2 McVittie spacetime, all three sets Ω A , Ω R and H are non-empty. However, in Class 1, this is not necessarily the case. Noting that

$$H ( t ) > H _ { 0 } \text{ \ for all } \ t > 0 \text{ \quad \ \ } ( 6 2 )$$

and that

$$\frac { 1 } { r } ( 1 - \frac { 2 m } { r } ) ^ { 1 / 2 } \leq \frac { 1 } { 3 \sqrt { 3 } m },$$

$$( 6 3 )$$

we see that a necessary and sufficient condition for the existence of a horizon and a regular region in a Class 1 McVittie spacetime is that ([13, 14])

$$\ m H _ { 0 } < \frac { 1 } { 3 \sqrt { 3 } }. \quad \quad \quad ( 6 4 )$$

We note that this is identically the necessary and sufficient condition for the existence of horizons in the Schwarzschild-de Sitter spacetime with mass parameter m and cosmological constant Λ = 3 H 2 0 . For Class 1 McVittie spacetimes, we restrict our attention to those cases where a horizon exists:

## Definition 2 [continued]

(vi) The horizon existence condition (64) is satisfied.

↦

From (61) and monotonicity of t → H t ( ), we can describe the horizon by

$$t = t _ { h } ( r ) = H _ { \text{inv} } ( \frac { 1 } { r } \sqrt { 1 - \frac { 2 m } { r } } ), \quad \quad ( 6 5 )$$

↦

↦

where H inv is the inverse of the function t → H t ( ): H inv ( H t ( )) = t for all t &gt; 0. The function r → t h ( r ) has a positive global minimum at r = 3 m , is decreasing on (2 m, 3 m ) and is increasing on (3 m, + ∞ ). We denote the global minimum by t h, 3 m . Then we may also describe the two branches of the horizon by functions

$$r _ { h } ^ { + } \colon [ t _ { h, 3 m }, + \infty ) \to [ 3 m, + \infty ) \quad \ \ ( 6 6 ) \quad \text{fig} \quad \text{ and}$$

and

$$r _ { h } ^ { - } \colon [ t _ { h, 3 m }, + \infty ) \to ( 2 m, 3 m ] \quad \ \ ( 6 7 ) \quad ^ { ( t, r } _ { \text{am} }$$

where by implicit differentiation we find

$$\frac { d r _ { h } ^ { \pm } } { d t } = - r ^ { 2 } ( 1 - \frac { 3 m } { r } ) ^ { - 1 } ( 1 - \frac { 2 m } { r } ) ^ { 1 / 2 } H ^ { \prime } ( t ) \Big | _ { r = r _ { h } ^ { \pm } }. \ \ ( 6 8 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

In a Class 1 McVittie spacetime, H t ( ) → H 0 &gt; 0 as t →∞ . Then on the horizon, r is restricted to the interval ( r -, r + ) where r ± are the larger ( r + ) and the smaller ( r -) of the two positive roots of 1 -2 m r -r 2 H 2 0 = 0, and we have t → + ∞ as r → r ± along the horizon. It is straightforward to prove that

$$2 m < r _ { - } < 3 m < r _ { + }. \quad \quad ( 6 9 ) \quad \text{col} _ { h \infty }$$

From (68), we see that on the inner branch , r -h decreases with t from r = 3 m to r = r -, with lim t → ∞ + r -h = r -. On the outer branch , r + h increases with t from r = 3 m to r = r + , and lim t → ∞ + r + h = r + . Note in particular that χ t, r ( ) &lt; 0 for all r &gt; r + . See Figure 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0044v1-figure-d594b486-5964-4065-8e6a-79d26270dd4e.jpg)













FIG. 3. The horizon of (i) a Class 1 McVittie spacetime (top) and (ii) a Class 2 McVittie spacetime (bottom). In both cases, t and r are shown on the horizontal and vertical axes respectively, and m = 1. The point at which the horizon first forms, ( t, r ) = ( t h, 3 m , 3 m ), is also shown in both. For the Class 1 example, we have taken H 0 = 0 15 and . H t ( ) = H 0 coth(3 H t 0 ). This puts the asymptotes of the horizon at r = r -= 2 26 (in-. ner branch) and r = r + = 5 24 (outer branch). For the Class . 2 example, we have taken H t ( ) = 2 5 . t -1 . The asymptote of the inner branch of the horizon is r = 2 m , and r →∞ on the outer branch. In both cases, the regular region of spacetime ( χ t, r ( ) &gt; 0) is the convex region bounded above and below by the horizon. The exterior of this region is the anti-trapped region ( χ t, r ( ) &lt; 0).

In a Class 2 McVittie spacetime, we have H 0 = lim t →∞ H t ( ) = 0. There is no restriction of r - except of course that r &gt; 2 m . From (68), we see that on the inner branch , r -h decreases with t from r = 3 m to r = 2 m , and lim t → ∞ + r -h = 2 m . On the outer branch , r -h increases with t from r = 3 m to r = + ∞ , and lim t → ∞ + r + h = + ∞ . See Figure 3.

The behaviour of ingoing radial null geodesics that are initially in the regular region of the spacetime is essen-

tially the same in Class 1 and Class 2 McVittie spacetimes, as the next result shows.

Proposition 5 Let γ be an IRNG of a Class 1 or a Class 2 McVittie spacetime with initial point γ (0) . If γ (0) ∈ Ω R ∪ H = { ( t, r ) ∈ Ω : χ t, r ( ) ≥ 0 } , then there exists s 0 ∈ (0 , ∞ ) such that lim s → s -0 ( ( t s , r ) ( s )) = (+ ∞ , r -) in Class 1, and lim s → s -0 ( ( t s , r ) ( s )) = (+ ∞ , 2 m ) in Class 2.

Proof: We use r in ( ) to denote the solution of the IRNG t equations considered as a curve in the t -r plane and we have

$$r ^ { \prime } _ { i n } ( t ) = r f ^ { 1 / 2 } H ( t ) - f. \text{ \quad \ \ } ( 7 0 ) \text{ \ \ } \text{ has}$$

Essentially by definition,

$$r _ { i n } ^ { \prime } ( t ) | _ { \mathcal { H } } = 0. \quad \quad \ \ ( 7 1 ) \quad \mathbf P _ { 1 }$$

It follows from the description given above that the horizon acts as a one-way membrane with respect to IRNGs, injecting the trajectories into the regular region. Thus an IRNG with initial point in Ω R remains in Ω R on the maximal interval of existence. Furthermore, it follows that the geodesic cannot meet H again at some finite value of t (this would violate uniqueness, as there is a unique IRNG at each point of the horizon which must enter rather than exit the regular region). It follows that r ′ in ( ) t &lt; 0 on the maximal interval of existence and that r &gt; r -along the geodesic. (We consider Class 1 explicitly. The proof is identical for Class 2, with r -replaced by 2 m .) Furthermore, equation (70) and the conditions on H t ( ) show that this slope is bounded below, and so global existence in t follows. By monotonicity, lim t →∞ r in ( ) t ≥ r -exists, and it is clear that the limit must in fact equal r -(any other limit would yield a contradiction when t → ∞ in (70)).

We prove that these geodesics have finite affine length to the future as follows. Using (14) with the null condition glyph[epsilon1] = 0 and the ingoing condition fr ˙ = ( rf -1 / 2 H -1) ˙ t we obtain

$$\ddot { t } = \left [ - ( 1 - \frac { m } { r } ) f ^ { - 1 / 2 } H + \frac { 2 m } { r ^ { 2 } } \right ] \dot { t } ^ { 2 } \quad \ \ ( 7 2 ) \quad \ \ ( \text{th} { 8 } \\ \lim.$$

along IRNGs. Given that the geodesic is (eventually) confined to Ω R , we have χ &gt; 0 or equivalently f -1 / 2 H &gt; r -1 (see (59)). This yields

$$\ddot { t } > u ( r ) \dot { t } ^ { 2 }, \quad u ( r ) = \frac { 3 m } { r ^ { 2 } } - \frac { 1 } { r }. \quad \ \ ( 7 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

↦

It is straightforward to show that r → u r ( ) is a strictly decreasing function for 0 &lt; r &lt; 6 m : this condition holds along the geodesics in question for sufficiently late parameter times as r → r -&lt; 3 m . Thus there exists s ∗ &lt; s ω such that (recalling ˙ r &lt; 0 along these geodesics)

$$u ( r ( s ) ) > u ( r ( s _ { * } ) ) > 0, \ \ s \in ( s _ { * }, s _ { \omega } ). \quad ( 7 4 ) \quad \text{RN} \quad.$$

Then

$$\ddot { t } > u ( r ( s _ { * } ) ) t ^ { 2 }, \ \ s \in ( s _ { * }, s _ { \omega } ). \quad \quad ( 7 5 )$$

Integrating twice shows that t → + ∞ in finite affine time - that is, s ω &lt; + ∞ . glyph[square]

Comment 6 In [11], which dealt with a subset of Class 2 McVittie spacetimes, it was erroneously claimed that the IRNGs of the previous result have infinite affine length to the future. A heuristic argument was put forth in [13] correcting this claim, and generalising to Class 1 McVittie spacetimes. The result finite affine length to the future - was proven rigorously in Appendix A of [14], using the argument repeated above. We have included the proof of finite affine length for purposes of continuity and clarity.

In Class 1, we can identify the following universal behaviour for all causal geodesics that are initially at sufficiently large radius.

Proposition 6 Let γ be a causal geodesic of a Class 1 McVittie spacetime with initial point γ (0) . If r (0) = r | γ (0) &gt; r + , then γ is future complete and r → ∞ as s →∞ .

Proof: We note first that χ t ( (0) , r (0)) &lt; 0 where t (0) = t | γ (0) . From (15), we see that ˙(0) r &gt; 0, and so the conditions r | γ s ( ) &gt; r + and ˙( r s ) &gt; 0 hold for s ∈ [0 , s ω ), the right-maximal interval of existence for the geodesic. Hence r s ( ) &gt; r (0) &gt; r + &gt; 3 m for all s ∈ [0 , s ω ). It follows from (14) that

$$\ddot { t } < \left ( 1 - \frac { 2 m } { r ( 0 ) } \right ) ^ { - 1 / 2 } H ( t ( 0 ) ), \text{ \quad \ \ } ( 7 6 )$$

where we have used H t ′ ( ) &lt; 0. Integrating shows that ˙ ( t s ) and t ( s ) are bounded on finite intervals. From (20) it follows that ˙( r s ) is bounded on finite intervals. Then a straightforward application of Theorem 1 (applied to the right-maximal rather than the left-maximal interval of existence) shows that s ω = + ∞ , i.e. the geodesic γ is future-complete. Since ˙( r s ) &gt; 0 for all s ∈ [0 , ∞ ), the limit r ∞ = lim s →∞ r s ( ) ≤ + ∞ exists. If r ∞ &lt; + ∞ (that is, if r has a finite limit along the geodesic), then lim s →∞ ˙ r = 0. From (20), we obtain

$$- \chi \dot { t } < ( \alpha + 1 ) \dot { r }.$$

Since r s ( ) &gt; r (0) &gt; r + , we know that χ is negative and bounded away from zero on the geodesic, and so ˙ r → 0 implies ˙ t → 0 as s → ∞ . But then taking the limit s →∞ in (15) yields glyph[lscript] 2 r 2 ∞ = . This yields a contradiction glyph[epsilon1] -and hence the conclusion that lim s →∞ r s ( ) = + ∞ -unless glyph[lscript] = glyph[epsilon1] = 0. So the theorem is proven except for the case of radial null geodesics. The case of outgoing RNGs is dealt with in Proposition 4. For ingoing RNGs satisfying the hypothesis of the theorem, we have the inequality

$$z ( s ) > z _ { + } = \sqrt { 1 - \frac { 2 m } { r _ { + } } }. \quad \quad \quad ( 7 8 )$$

Then

$$\frac { d z } { d t } & > \frac { ( 1 - z ^ { 2 } ) } { 2 } ( H - \frac { z } { 2 m } ( 1 - z _ { + } ^ { 2 } ) ) \\ & = \frac { ( 1 - z ^ { 2 } ) } { 2 } ( H - \frac { z } { z _ { + } } H _ { 0 } ) \\ & > \frac { H _ { 0 } } { 2 } ( 1 - z ^ { 2 } ) ( 1 - \frac { z } { z _ { + } } ), \\ \mathcal { U } ^ { \prime } / \omega \succ \partial _ { t \infty } + \mathfrak { h } \mathcal { H } / \omega \succ \mathcal { U } \mathcal { U } \mathcal { H } \mathcal { H } \mathcal { H } ^ { 2 }.$$

using H t ′ ( ) &lt; 0 (so that H t ( ) &gt; H 0 ) and H 0 = z + 2 m (1 -z 2 + ) which follows from the definition of r + . Integrating (and using the initial condition z t ( 0 ) &gt; z + ) yields z → 1 (i.e. r → + ∞ ) as t →∞ . glyph[square]

We can identify two further special cases of IRNGs in Class 1 McVittie spacetimes.

## Proposition 7 In a Class 1 McVittie spacetime,

- (a) there exists a 1-parameter family of IRNGs with ( t, r ) → (+ ∞ , r -) and r t ( ) &lt; r -for all t &lt; + ∞ . These geodesics satisfy γ s ( ) ∈ Ω A throughout their interval of existence and have finite affine length.
- (b) There exists a unique IRNG with ( t, r ) → (+ ∞ , r + ) ; this geodesic is future-complete.

To prove these results, we use a dynamical systems approach. It is useful towards this end to consider as independent variables z t ( ) (i.e. the coordinate z considered as a function of t along the geodesic) and ξ = H t ( ). The evolution of z in terms of the state variables ( z, ξ ) is given by (30) (lower sign for IRNGs). The corresponding equation for the evolution of H is furnished by the following lemma. This lemma is central to the analysis of Section VI, and so results for both Classes 1 and 2 are given here.

Lemma 5 In any expanding McVittie spacetime with a big bang background, there exists a function w ∈ C 1 ([ H , 0 + ∞ ) , R -) such that

$$H ^ { \prime } ( t ) = w ( H ( t ) ). \text{ \quad \ \ } ( 7 9 )$$

In a Class 1 McVittie spacetime, this function satisfies

$$w ( H _ { 0 } ) = 0, \ \ w ^ { \prime } ( H _ { 0 } ) = - 2 H _ { 0 } \kappa, \quad \ ( 8 0 ) \quad \text{that} _ { - }$$

and in a Class 2 McVittie spacetime,

$$w ( 0 ) = w ^ { \prime } ( 0 ) = 0, \ \ w ^ { \prime \prime } ( 0 ) = - 2 \kappa, \quad \ ( 8 1 ) \quad \stackrel { \Theta \, \Theta _ { \Gamma } } { \text{rem} }$$

where κ is the constant defined in part (v) of Definition 2 and Definition 3. Furthermore, this function extends to a function w ∈ C 1 ( R R , ) .

↦

Proof: Since t → H t ( ) is monotone decreasing and C 1 , this function has a monotone C 1 inverse H inv . Thus we can write

$$\xi = H ( t ) \Leftrightarrow t = H _ { \text{inv} } ( \xi ), \text{ \quad \ \ } ( 8 2 ) \quad \text{def}$$

where H inv ∈ C 1 (( H , 0 + ∞ ) , R ) and satisfies lim ξ → H + 0 H inv ( ξ ) = + ∞ and lim ξ → ∞ + H inv ( ξ ) = 0. Then we can write

$$H ^ { \prime } ( t ) = w ( H ( t ) )$$

for the C 1 function w = H ′ ◦ H inv . This proves existence of the required function w and shows that w ∈ C 1 (( H , 0 + ∞ ) , R -), since we know by hypothesis that H t ′ ( ) &lt; 0 for all t &gt; 0. To prove the second part, we establish the equalities in limiting form (e.g. lim ξ → H + 0 w ξ ( ) = 0): this shows that w extends to a function with the required properties. For convenience, we use the same name w for this extension.

To derive the required limits, we take the derivative of the background field equations (9) and (10) with respect to t (with p 0 = g ρ ( 0 )) and eliminate

$$\rho _ { 0 } ^ { \prime } ( t ) = \frac { 3 } { 4 \pi } H ( t ) H ^ { \prime } ( t )$$

to obtain

$$\frac { H ^ { \prime \prime } ( t ) } { H ^ { \prime } ( t ) } = - 3 H ( t ) ( 1 + g ^ { \prime } ( \rho _ { 0 } ) ). \quad \quad ( 8 5 )$$

Taking the limit t →∞ and recalling that ρ 0 → 0 in this limit yields

$$\lim _ { t \to \infty } \frac { H ^ { \prime \prime } ( t ) } { H ^ { \prime } ( t ) } = - 2 H _ { 0 } \kappa,$$

using Definition 2(v). On the other hand, we note that

$$w ^ { \prime } ( \xi ) & = \frac { d } { d t } \left \{ H ^ { \prime } ( t ) \right \} \frac { d t } { d \xi } \\ & = H ^ { \prime \prime } ( t ) \frac { d t } { d H } = \frac { H ^ { \prime \prime } ( t ) } { H ^ { \prime } ( t ) }.$$

The result regarding the limit follows by noting that ξ → H + 0 ⇔ → t + ∞ . Using Definition 2(iv), we have

$$\lim _ { \xi \to H _ { 0 } ^ { + } } w ( \xi ) = \lim _ { t \to \infty } H ^ { \prime } ( t ) = 0. \quad \quad ( 8 8 )$$

Thus the domain of w may be extended to include H 0 , and the required one-sided derivative exists to guarantee that w ∈ C 1 ([ H , 0 ∞ )).

This completes the proof for Class 1 McVittie spacetimes. For Class 2, the results w (0) = w ′ (0) = 0 arise as a special case ( H 0 = 0) of the Class 1 results, so it only remains to show that w ′′ (0) = -2 κ . Combining the ξ derivative of (87) with the t derivative of (85) and using (84) yields

$$w ^ { \prime \prime } ( \xi ) = - 3 ( 1 + g ^ { \prime } ( \rho _ { 0 } ) ) - \frac { 9 H ^ { 2 } } { 4 \pi } g ^ { \prime \prime } ( \rho _ { 0 } ). \quad ( 8 9 )$$

Taking the limit ξ → 0 + yields the required result, recalling the definition of κ and that g ∈ C 2 [0 , + ∞ ).

To extend the domain of w to the whole real line, we define w H ( 0 -ξ ) = -w H ( 0 + ) for ξ ξ &gt; 0. glyph[square]

Proof of Proposition 7 The equation governing IRNGs is (30) with the lower sign. Using Lemma 5 and taking ξ = H t ( ), we can write this as a 2-dimensional dynamical system:

$$\frac { d z } { \delta t } = \frac { \xi } { 2 } ( 1 - z ^ { 2 } ) - \frac { z } { 4 m } ( 1 - z ^ { 2 } ) ^ { 2 }, \quad \ \ ( 9 0 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\frac { d \xi } { d t } = w ( \xi ). \quad \begin{matrix} & \text{sum} \\ ( 9 1 ) & \text{can} \\ & \text{is tl} \end{matrix}$$

Then (90) and (91) define a C 1 dynamical system on R 2 . The points P 1 : ( z, ξ ) = ( z -, H 0 ) and P 2 : ( z, ξ ) = ( z + , H 0 ) are equilibrium points of this system, where z ± := z r ( ± ). Both are hyperbolic equilibrium points [16], with eigenvalues

$$\lambda _ { 1 } = - \frac { 1 } { 4 m } ( 1 - z _ { \pm } ^ { 2 } ) ( 1 - 3 z _ { \pm } ^ { 2 } ), \quad \lambda _ { 2 } = - 2 H _ { 0 } \kappa. \ ( 9 2 ) \quad \text{Class} \quad \text{correc} \quad \text{locio} \quad.$$

Thus λ 2 &lt; 0 in both cases, but λ 1 &lt; 0 for P 1 and λ 1 &gt; 0 for P 2 . This follows from (69) and the definition (29) of z .

It follows that P 1 is an asymptotically stable equilibrium point: there exists a neighbourhood of P 1 in R 2 such that every solution of (90) and (91) with initial point in this neighbourhood satisfies lim t → ∞ + ( z t , ξ ( ) ( )) t = ( z -, H 0 ). Since the stable manifold is 2-dimensional, there are trajectories that approach the equilibrium point along every tangent at P 1 . The 1-parameter family of the statement of the theorem correspond to trajectories of the dynamical system that approach the equilibrium point along directions with z &lt; z -. Thus these geodesics satisfy z &lt; z -at late times, and so are confined to the anti-trapped region Ω A at late times ( z -&lt; z for points in Ω R ∪ H ). But z ′ ( ) t &gt; 0 for IRNGs in Ω A , so the geodesics are confined to Ω A for all times.

To complete the proof of part (a) of the proposition, it remains to prove that these geodesics have finite affine length. Proposition 3 establishes that the geodesics have finite affine length to the past. To show finite affine length to the future, we note that the coefficient of ˙ t 2 in the right hand side of (72) satsfies

$$\lim _ { s \to s _ { \omega } } - ( 1 - \frac { m } { r } ) f ^ { - 1 / 2 } H + \frac { 2 m } { r ^ { 2 } } = u ( r _ { - } ), \quad ( 9 3 )$$

(see (73)) where ( s α , s ω ) is the maximal interval of existence of the geodesic. Since r -&lt; 3 m , this coefficient is positive, and so there exists a positive constant k 2 and s ∗ &lt; s ω such that

$$\ddot { t } > k ^ { 2 } \dot { t } ^ { 2 } \quad \text{for all } s \in ( s _ { * }, s _ { \omega } ). \quad \quad ( 9 4 ) \quad \Pr _ { \cdot \cdot \cdot }$$

Integrating twice shows that t diverges to + ∞ in finite s completing the proof of part (a).

For part (b), we note that P 2 is a saddle point of the dynamical system as λ 1 &gt; 0 , λ 2 &lt; 0. The stable manifold theorem tells us that the nonlinear system has a 1-dimensional stable manifold W s , tangent to the 1-dimensional stable space E s of the linearized system at P 2 . This stable manifold corresponds to a unique

IRNG satisfying lim t → ∞ + ( z, ξ ) = ( z + , H 0 ). Along this geodesic, we have

$$\lim _ { s \to s _ { \omega } } - ( 1 - \frac { m } { r } ) f ^ { - 1 / 2 } H + \frac { 2 m } { r ^ { 2 } } = u ( r _ { + } ) < 0, \quad ( 9 5 )$$

since r + &gt; 3 m . Thus ¨ is eventually negative, t and so cannot diverge to + ∞ in finite affine time. This geodesic is therefore future-complete. glyph[square]

Comment 7 In [14], the Penrose-Carter diagram for the McVittie spacetime with H t ( ) = H 0 coth ( 3 H t 0 2 ) was constructed using a combination of analytical and numerical approaches. The results of Propositions 3-7 show that the resulting structure applies for all McVittie spacetimes of Class 1. This reinforces the results of [13], modulo the corrections of [14], e.g. the geodesic nature of the 'cosmological horizon' r = r + ( ) t in the language and notation of [13] (this is the geodesic η 1 of [14] and is the unique IRNG of Proposition 7(b) above) and the existence of the geodesics identified in Proposition 7(a) above (also identified numerically in [14]). We note also that the results establishing the affine lengths of the geodesics in Proposition 7 we first given in [14].

For Class 2 McVittie spacetimes, the asymptotic behaviour of ingoing radial null geodesics (IRNGs) is sensitive to the value of the constant κ introduced in Definition 3. Proposition 8 generalises results of [11]: the following lemma is required for its proof.

Lemma 6 In a Class 2 McVittie spacetime,

$$\lim _ { t \to \infty } - \frac { H ^ { \prime } ( t ) } { H ( t ) ^ { 2 } } = \kappa, \text{ \quad \quad } ( 9 6 )$$

where κ is the constant defined in part (v) of Definition 3.

Proof: From the Einstein equations (9), (10) and Definitions 1 and 3, we see that lim t →∞ ρ 0 = lim t →∞ p 0 = 0, and that

$$\frac { p _ { 0 } } { \rho _ { 0 } } = - \frac { 2 } { 3 } \frac { H ^ { \prime } } { H ^ { 2 } } - 1.$$

Appealing to part (v) of Definition 3 and taking the limit t →∞ yields the result. glyph[square]

Proposition 8 Let γ be an ingoing radial null geodesic of a Class 2 McVittie spacetime, and let κ be as defined in (56).

(a) If κ &lt; 1 and γ (0) ∈ Ω A with r (0) = r | γ (0) sufficiently large, then γ remains in Ω A and is future complete with lim s →∞ ( ( t s , r ) ( s )) = (+ ∞ , + ∞ ) ;

(b) If κ &gt; 1 and γ (0) ∈ Ω A , then γ enters Ω R and terminates at ∂ Ω 2 m in finite parameter time.

Proof of part (a): We prove this result by identifying a family of curves that extends to infinity, and then show that each IRNG satisfying the hypothesis of the proposition must remain above a member of this family.

We introduce a 1-parameter family { Σ ν : ν &gt; 1 } of curves in Ω A = { ( t, r ) ∈ Ω : χ t, r ( ) &lt; 0 } . These are defined to be the curves at each point of which the unique IRNG through that point has slope r ′ in ( ) = t ν -1. Then

$$\Sigma _ { \nu } = \{ ( t, r ) \in \Omega \colon H ( t ) = \frac { 1 } { r } ( \nu - \frac { 2 m } { r } ) ( 1 - \frac { 2 m } { r } ) ^ { - 1 / 2 } \}. \quad \begin{array} { c } \Box \\ \text{Proo} \\ \text{in the} \end{array}$$

↦

Notice that the horizon corresponds to Σ 1 . It is straightforward to show that { Σ ν : ν &gt; 1 } provides a foliation of Ω , the anti-trapped region of the spacetime. A Since t → H t ( ) is monotone, we can invert (98) to obtain a 1-parameter family of functions t ν ∈ C 1 ((2 m, + ∞ ) , R + ) such that

$$( t, r ) \in \Sigma _ { \nu } \Leftrightarrow t = t _ { \nu } ( r ). \quad \quad ( 9 9 ) \quad \underset { \ m a } { \omega }$$

A straightforward calculation shows that t ′ ν ( r ) &gt; 0 for all r &gt; 3 m and all ν ≥ 1, and so there exists a 1-parameter family of functions r ν ∈ C 1 ( R + , (3 m, ∞ )) such that

$$( t, r ) \in \Sigma _ { \nu } \cap \{ ( t, r ) \in \Omega _ { A } \, \colon r > 3 m \} \Leftrightarrow r = r _ { \nu } ( t ). \ ( 1 0 0 ) \quad \text{Since}$$

This amounts to saying that the set { ( t, r ) ∈ Ω A : r &gt; 3 m } may be foliation by C 1 curves in the t -r plane, indexed by ν &gt; 1, each one of which has the property that the (different) IRNGs intersecting that curve all have slope r ′ in ( ) t = ν -1, and we note that r ν ( ) t → + ∞ as t → + ∞ (on account of H t ( ) → 0 in this limit).

We now show that r in ( ) t &gt; r ν ( ) for all t t &gt; t 0 and for an appropriate value of ν .

Implicit differentiation of the defining relation of Σ ν in (98) yields

$$\frac { d r _ { \nu } } { d t } = - \, \frac { H ^ { \prime } } { H ^ { 2 } } \left ( \nu - \frac { 2 m } { r } \right ) ^ { 2 } \left ( 1 - \frac { 2 m } { r } \right ) ^ { 1 / 2 } \left ( u _ { \nu } ( \frac { m } { r } ) \right ) ^ { - 1 } \Big | _ { r = r _ { \nu } ( t ) \ \kappa > 1, \, \text{where we} }, \, \\ \text{where } w h o r a \ w \ w \ - \, \nu \ - \, \nu \ \mp \ \Lambda \wr r \ \pm \, \kappa \wr ^ { 2 } \text{ and where two have} } \quad \ w \ h o r a \ h o r a$$

where u ν ( x ) = ν -( ν + 4) x + 6 x 2 and where we have repeatedly used (98) to make convenient substitutions. It follows from Lemma 6 that

$$\frac { d r _ { \nu } } { d t } = \kappa \left ( \nu + O ( \frac { m } { r _ { \nu } ( t ) } ) ^ { 2 } \right ), \ \ t \to \infty. \quad \ ( 1 0 2 ) \quad \text{This} _ { \iota \, - \iota \, - \iota \, - \iota }$$

Then we can calculate

$$& \frac { d } { d t } \left ( \frac { r _ { i n } } { r _ { \nu } } \right ) \Big | _ { r _ { i n } = r _ { \nu } } = \frac { 1 } { r _ { \nu } ( t ) } \left ( \nu - 1 - \kappa \nu + O ( \frac { m } { r _ { \nu } ( t ) } ) ^ { 2 } \right ), \quad t \to \frac { \nu _ { * } < \nu _ { \nu _ { \nu } } ^ { 2 } } { r _ { i n } ^ { \prime } = } \\ & \Lambda \sim \times 1 \quad \text{and} \, \text{and} \, \text{do} \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \incolon$$

As κ &lt; 1, we can choose ν &gt; 1 such that ν -1 -κν &gt; 0 (we require ν &gt; 1 / (1 -κ )). From (103), we see that there exists r ∗ &gt; m 3 such that

$$\frac { d } { d t } \left ( \frac { r _ { i n } } { r _ { \nu } } \right ) \Big | _ { r _ { i n } = r _ { \nu } } > 0 \text{ \ for all } r _ { \nu } ( t ) > r _ { * }. \text{ \ \ } ( 1 0 4 ) \text{ \ \ } \text{Since} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text
}$$

Now consider the IRNG that passes through the point ( t 0 , r 0 ) ∈ Σ ν with r 0 &gt; r ∗ . According to (104), the geodesic crosses r = r ν ( ) from below, and cannot cross t this curve again for t &gt; t 0 ( t → r ν ( ) t is increasing). Thus r in ( ) t &gt; r ν ( ) t for all t &gt; t 0 , and so r in ( ) t → + ∞ as t → ∞ along the geodesic. By uniqueness, the same conclusion holds for all IRNGs satisfying r | t 0 &gt; r 0 . The usual convexity argument based on (13) indicates that these geodesics have infinite affine length to the future. glyph[square]

Proof of part (b): Again, let r in ( ) denote the IRNG t in the t -r plane, but with the additional condition that r in ( t 0 ) = r 0 where ( t 0 , r 0 ) = ( t, r ) | γ (0) ∈ Ω . A Then we initially have r ′ in ( ) t &gt; 0. The IRNG thus either meets the horizon whereat r ′ in = 0 and then enters Ω R , or remains outside the horizon for all t &gt; t 0 . OUr aim now is to rule out the latter case. In this latter case, there exists t 1 &gt; t 0 such that lim t → t 1 r in ( ) = t ∞ (this follows from the fact that r → + ∞ on the relevant branch of the horizon). We may rule this out by the following convexity argument. First, we calculate that

$$\text{amer} \quad & r ^ { \prime \prime } _ { i n } ( t ) = r f ^ { 1 / 2 } H ^ { \prime } ( t ) + ( ( 1 - \frac { m } { r } ) H ( t ) - \frac { 2 m } { r ^ { 2 } } f ^ { 1 / 2 } ) ( r H ( t ) - f ^ { 1 / 2 } ). \\ ( 1 0 0 ) \quad & \circ \dots \dots \omega \dots \dots \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \dots \dots \omega \end{cases}$$

Since the geodesic remains in the anti-trapped region Ω A , we have H &gt; r -1 f 1 / 2 . This allows us to identify the dominant terms in (105). If t 1 &lt; + ∞ , we have

$$\frac { r _ { i n } ^ { \prime \prime } } { r _ { i n } } \sim H ^ { \prime } ( t _ { 1 } ) + H ( t _ { 1 } ) ^ { 2 }, \quad t \to t _ { 1 }. \quad ( 1 0 6 )$$

But this rules out r in →∞ in finite t . So we must have t 1 = + ∞ . In this case, we find

$$\frac { r _ { i n } ^ { \prime \prime } } { r _ { i n } } & \sim H ^ { \prime } ( t ) + H ( t ) ^ { 2 }, \quad t \to \infty \\ & \sim ( 1 - \kappa ) H ( t ) ^ { 2 }, \quad t \to \infty$$

, where we have used Lemma 6. So by the hypothesis that κ &gt; 1, we see that r ′′ in ( ) t &lt; 0 for sufficiently large t . Now choose t 2 &gt; t 0 such that r ′′ in ( ) t &lt; 0 for all t ≥ t 2 , and let ν 2 be the value of ν &gt; 1 such that ( t 2 , r in ( t 2 )) ∈ Σ ν 2 . We consider the set

$$S = \{ \nu \geq 1 \, \colon r _ { i n } ( t ) > r _ { \nu } ( t ) \text{ for all } t \geq t _ { 2 } \}. \quad ( 1 0 7 )$$

→∞ . This set is non-empty (it contains ν = ν 2 ) and is bounded below by ν = 1, since Σ 1 is the horizon which r in does not cross. Hence S has a greatest lower bound ν ∗ with 1 ≤ ν ∗ &lt; ν 2 . Then for any glyph[epsilon1] ∈ (0 , ν 2 -ν ∗ ), r in crosses Σ ν ∗ + glyph[epsilon1] from below at some time t glyph[epsilon1] . At the point of crossing, r ′ in = ν ∗ + , and by convexity, glyph[epsilon1] r ′ in ( ) t &lt; ν ∗ + glyph[epsilon1] for all t &gt; t glyph[epsilon1] . Now consider the surface Σ ν ∗ -δ (or r = r ν ∗ -δ ). From (102), r = r ν ∗ -δ ( ) has slope obeying t

$$r ^ { \prime } _ { \nu _ { * } - \delta } ( t ) \sim \kappa ( \nu _ { * } - \delta ), \ \ t \to + \infty. \quad ( 1 0 8 )$$

Since κ &gt; 1, we can choose glyph[epsilon1] and δ so that κ ν ( ∗ -δ ) &gt; ν ∗ + glyph[epsilon1] . This condition on the asymptotic slopes, along with

the convexity property r ′ in ( ) t &lt; ν ∗ + glyph[epsilon1] for all t &gt; t glyph[epsilon1] , tells us that r in crosses r ν ∗ -δ . This contradicts the infimum property of ν ∗ , and so we rule out the possibility that r in does not cross the horizon. The evolution of the geodesic is thus subject to Proposition 5. glyph[square]

Comment 8 It is evident from the previous proposition that there is no general result governing the future evolution of causal geodesics in a Class 2 McVittie spacetime as we have in Proposition 6 for Class 1. At the very least, this behaviour depends on the value of κ . However, we can conclude that if κ &lt; 1, then a result corresponding to Proposition 6 applies. For κ ≥ 1, we have not been able to draw any general conclusions, and it may be the case that no causal geodesics escape to infinity, other than the outgoing radial null geodesics covered by Proposition 4.

Proposition 9 Let γ be a causal geodesic of a Class 2 McVittie spacetime with κ &lt; 1 . If γ (0) ∈ Ω A and r (0) = r | γ (0) sufficiently large, then γ remains in Ω A and is future complete with lim s →∞ ( ( t s , r ) ( s )) = (+ ∞ , + ∞ ) .

Proof: This is a straightforward application of the 'radial confinement' result Lemma 4. As per Propositions 4 and 7, both the ingoing and outgoing radial null geodesics through γ (0) remain in Ω A , are future complete and extend to infinite values of both t and r . Thus χ is negative, and so by (15) ˙ r is positive, along the maximally extended geodesic. From (14), we then obtain

$$\ddot { t } \leq ( 1 - \frac { 2 m } { r ( 0 ) } ) ^ { - 1 / 2 } H ( t ( 0 ) ), \quad \ \ ( 1 0 9 )$$

where we have assumed (without loss of generality) that r (0) &gt; 3 m , and we have used - ≤ glyph[epsilon1] 1. Integrating shows that ˙ t remains finite for finite parameter times s , and then Lemma 1 yields the same for ˙. Applying Theorem 1 r yields global existence for the geodesic, which, returning to the radial confinement result, completes the proof. glyph[square]

## VI. BOUND PARTICLE AND PHOTON ORBITS

We now come to our main results: we prove that in any Class 1 or Class 2 McVittie spacetime, there are bound photon and particle orbits. That is, there are futurecomplete causal geodesics satisfying r s ( ) &lt; + ∞ for all s &gt; 0. For photon orbits, we prove that these orbits approach r = 3 m as s →∞ . For particle orbits, we prove that the orbits are asymptotic to a radius corresponding to a stable circular orbit of the associated Scwarzschild(de Sitter) spacetime (Schwarzschild-de Sitter for Class 1, Schwarzschild for Class 2). To prove the relevant results, we use a dynamical systems approach. The proofs run along somewhat different lines for the two classes, and so we consider them separately. However, the initial steps are the same for both classes.

To begin, we write the geodesic equations as a dynamical system - that is, as a first order autonomous system of ODEs. We introduce the variables

$$x = ( x _ { 1 }, x _ { 2 }, x _ { 3 } ) = ( H, r, \dot { r } ), \ \ y = \dot { t }. \quad ( 1 1 0 )$$

Then using (13) and (14) we can write the geodesic equations as a 4-dimensional system

$$( \dot { x }, \dot { y } ) = G ( x, y ), \text{ \quad \ \ } ( 1 1 1 )$$

for some function G : R 4 → R 4 .

The condition (15) corresponds to a zero-order constraint:

$$( 1 5 ) \Leftrightarrow g ( x, y ) = 0, \text{ \quad \ \ } ( 1 1 2 )$$

where

$$g ( x, y ) \coloneqq - \chi y ^ { 2 } - 2 \alpha x _ { 3 } y + \lambda x _ { 3 } ^ { 2 } + \frac { \ell ^ { 2 } } { x _ { 2 } ^ { 2 } } + | \epsilon | \quad ( 1 1 3 )$$

and we note that making the appropriate substitutions we have

$$\chi = f - x _ { 1 } ^ { 2 } x _ { 2 } ^ { 2 }, \text{ \quad \quad } ( 1 1 4 )$$

$$\alpha = x _ { 1 } x _ { 2 } f ^ { - 1 / 2 }, \text{ \quad \ \ } ( 1 1 5 )$$

$$\lambda = f ^ { - 1 }, \quad \quad \quad ( 1 1 6 )$$

$$f = 1 - \frac { 2 m } { x _ { 2 } }. \quad \quad \quad ( 1 1 7 )$$

glyph[negationslash]

The implicit function theorem allows us to solve g x, y ( ) = 0 locally. Provided χy + αx 3 = 0 at some point P , there is a neighbourhood of P such that

$$g ( x, y ) = 0 \Leftrightarrow y = u ( x ) \text{ \quad \quad } ( 1 1 8 )$$

for some C 1 function u . We use this to reduce the dimension of the problem. That is, we eliminate y from G x, y ( ) and write the geodesic equations as

$$\dot { x } = F ( x ), \text{ \quad \ \ } ( 1 1 9 )$$

where F : R 3 → R 3 with components read off from (13) and (14), with the substitution ˙ t → y = u x ( ).

One other technical manoeuvre is required in this process. We note that Lemma 5 allows us to write H t ′ ( ) = w H t ( ( )) for the function w introduced in that lemma, and so

$$x _ { 1 } = \frac { d H } { d s } = w ( x _ { 1 } ) \dot { t }. \text{ \quad \ \ } ( 1 2 0 )$$

It is at this point that the methods of analysis diverge for the two classes. The difference arises due to the different structures of the function w in the two cases.

## A. Bound orbits in Class 1 McVittie spacetimes

Throughout this subsection, we assume that our metric is that of a Class 1 McVittie spacetime as defined by Definitions 1 and 2. To begin, we recap some relevant results (including Lemma 5) that allow us to write the geodesic equations in dynamical systems form.

Lemma 7 The geodesic equation for x 1 can be written as

$$\dot { x _ { 1 } } = w ( x _ { 1 } ) u ( x ), \text{ \quad \ \ } ( 1 2 1 ) \text{ \ \ } a l _ { \ell }$$

where w and u are as defined in (79), (113) and (118). The function w extends to a function w ∈ C 1 ( R R , ) satisfying

$$w ( H _ { 0 } ) = 0, \ \ w ^ { \prime } ( H _ { 0 } ) = - 2 H _ { 0 } \kappa, \quad \ ( 1 2 2 )$$

where κ is defined in part (v) of Definition 2, and we have

$$w ( x _ { 1 } ) < 0, \ \ \text{ for all } x _ { 1 } > H _ { 0 }, \quad \ \ ( 1 2 3 ) \quad \text{and}$$

$$w ( H _ { 0 } - \xi ) = - w ( H _ { 0 } + \xi ), \ \ \text{ for all } \xi > 0. \quad ( 1 2 4 ) \quad \ \ F$$

glyph[square]

By making the appropriate substitutions in (14), we can write down the other equations of the dynamical system:

$$\dot { x _ { 2 } } & = x _ { 3 }, & ( 1 2 5 ) & \overset { \ a o e s \, n } { \mathbb { R } } \, \text{ wit} \\ \dot { x _ { 3 } } & = x _ { 2 } f ^ { 1 / 2 } w ( x _ { 1 } ) u ( x ) ^ { 2 } + ( 1 - \frac { 3 m } { x _ { 2 } } ) \frac { \ell ^ { 2 } } { x _ { 2 } ^ { 3 } } + \epsilon ( \frac { m } { x _ { 2 } ^ { 2 } } - x _ { 1 } ^ { 2 } x _ { 2 } ). & \overset { \ a o e s \, n } { \mathbb { R } } \, \text{ with } \\ & \text{ then }$$

## 1. Photon orbits in Class 1

We focus on the point P : x = ( H , 0 3 m, 0). We find that

$$\chi ( P ) = \chi _ { | _ { r = 3 m, H = H _ { 0 } } } = \frac { 1 } { 3 } - 9 m ^ { 2 } H _ { 0 } ^ { 2 } > 0, \quad ( 1 2 7 )$$

and so χ is positive in a neighbourhood of P (note the application here of the horizon existence condition (64)). Recalling that ˙ t &gt; 0 along all causal geodesics, we see that taking the appropriate solution of g x, y ( ) = 0 in (112) yields

$$u ( x ) = u _ { 1 } ( x ) \coloneqq - \frac { \alpha } { \chi } x _ { 3 } + \frac { 1 } { \chi } ( x _ { 3 } ^ { 2 } + \chi \frac { \ell ^ { 2 } } { x _ { 2 } ^ { 2 } } ) ^ { 1 / 2 }. \quad ( 1 2 8 ) \quad \text{by th} \\ \text{such}$$

This is a smooth function of x is a neighbourhood of P . For clarity, we collect some relevant results as follows.

Lemma 8 In a neighbourhood of P : x = ( H , 0 3 m, 0) , the equations governing null geodesics of a Class 1 McVittie spacetime may be written as (125), (126) with glyph[epsilon1] = 0 and (121) with u = u 1 and u 1 defined in (128). Furthermore, P is an equilibrium point of this system. | square

Our main result is obtained by a relatively straightforward application of standard results on dynamical systems.

Proposition 10 In a Class 1 McVittie spacetime, there exists a 3-parameter family of null geodesics that are future complete, with r → 3 m and t → + ∞ as s → ∞ along each member of the family. These geodesics correspond to bound photon orbits of the spacetime.

Proof: Our starting point is the dynamical system

$$\dot { x } = F ( x )$$

with

$$F _ { 1 } ( x ) = w ( x _ { 1 } ) u _ { 1 } ( x ), \quad F _ { 2 } ( x ) = x _ { 3 } \quad \ \ ( 1 3 0 )$$

and

$$& \downarrow ) \quad F _ { 3 } ( x ) = x _ { 2 } f ^ { 1 / 2 } w ( x _ { 1 } ) u _ { 1 } ( x ) ^ { 2 } + ( 1 - \frac { 3 m } { x _ { 2 } } ) \frac { \ell ^ { 2 } } { x _ { 2 } ^ { 3 } }. \quad ( 1 3 1 ) \\ & -$$

This is a C 1 dynamical system on a neighbourhood O P ⊆ R 3 of P . We note that the projection

$$\pi \colon \mathcal { O } _ { P } \to \mathbb { R } ^ { 2 } \colon x = ( x _ { 1 }, x _ { 2 }, x _ { 3 } ) \to ( x _ { 1 }, x _ { 2 } ) = ( H, r ) \\. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

does not project into Ω, as we must include values of H ∈ R with H &lt; H 0 . Such points are not in the spacetime: the projection is into Ω if and only if x 1 = H &gt; H 0 .

The Jacobian matrix at P is given by

$$J ( P ) = \left ( \begin{array} { c c c } a & 0 & 0 \\ 0 & 0 & 1 \\ b & c & 0 \end{array} \right ), \quad \quad \quad ( 1 3 3 )$$

where

$$a = - \kappa ^ { 2 } H _ { 0 } \frac { \ell } { m } ( \frac { 1 } { 3 } - 9 m ^ { 2 } H _ { 0 } ^ { 2 } ) ^ { - 1 / 2 }, \quad \ ( 1 3 4 )$$

$$b = - \underset { \lambda } { \sqrt { 3 } } \overset { \cdots } { H _ { 0 } } \frac { \check { \ell } ^ { 2 } } { m } ( \frac { 1 } { 3 } - 9 m ^ { 2 } H _ { 0 } ^ { 2 } ) ^ { - 1 }, \quad ( 1 3 5 )$$

$$c = \frac { \ell ^ { 2 } } { ( 3 m ) ^ { 4 } }. \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
c = \frac { \ell ^ { 2 } } { ( 3 m ) ^ { 4 } }.$$

This has eigenvalues σ = a, ± √ c , all of which are real, and so P is a hyperbolic equilibrium point of the system. Two of the eigenvalues ( a, - √ c ) are negative. Thus, by the Stable Manifold Theorem (see e.g.[16]), there exists a 2-dimensional differentiable manifold W s ⊆ O P such that any solution of (129) with x (0) ∈ W s satisfies lim s →∞ x s ( ) = P . This manifold is tangent to E s at P , where E s is the 2-dimensional stable space of the linearized system

$$\dot { \delta x } = J ( P ) \delta x. \text{ \quad \quad \ \ } ( 1 3 7 )$$

By calculating the relevant eigenvectors, we find that

$$E ^ { s } = \{ x ( P ) + \xi \vec { v } _ { 1 } + \eta \vec { v } _ { 2 } \colon \xi, \eta \in \mathbb { R } \}, \quad \ \ ( 1 3 8 )$$

where

$$\vec { v } _ { 1 } = ( ( c + a ^ { 2 } ) / b, 1, a ), \quad \vec { v } _ { 2 } = ( 0, 1, - \sqrt { c } ). \quad ( 1 3 9 )$$

So far, these results apply only to the 'unphysical' O P . To see the relevance to Ω, and hence complete the proof, we note that those trajectories of W s which have tangent at P corresponding to ξ &lt; 0 in (138) satisfy x 1 ( s ) &gt; H 0 for all sufficiently large s (note that b &lt; 0). These trajectories, of which there is a 3-parameter family (two parameters coming from the dimension of W s and glyph[lscript] providing the third) correspond to null geodesics of the spacetime which are future complete and have the asymptotic behaviour given in the statement of the proposition. glyph[square]

Comment 9 We note that the orbits described in the previous result are asymptotic to the unique circular photon orbit r = 3 m of the Schwarzschild-de Sitter spacetime with positive cosmological constant Λ = 3 H 2 0 . They are essentially as unstable as this unstable orbit at r = 3 m . The photon orbits of Schwarzschild-de Sitter spacetime are governed by the equation

$$\ddot { r } = \left ( 1 - \frac { 3 m } { r } \right ) \frac { \ell ^ { 2 } } { r ^ { 3 } }, \text{ \quad \ \ } ( 1 4 0 )$$

which is (13) with the null condition glyph[epsilon1] = 0 and the Schwarzschild-de Sitter condition H t ′ ( ) ≡ 0. In terms of the dynamical system of Proposition 10, this equation corresponds to the invariant manifold { x : x 1 = H 0 } of the flow. Then the circular photon orbit corresponds to the stationary point ( x , x 2 3 ) = (3 m, 0). This is a saddle of the 2-dimensional flow: thus, as is well know, the circular photon orbit of Schwarzschild-de Sitter spacetime is unstable. In terms of initial data, for a given r (0) sufficiently close to 3 m , there are precisely (only) two values of ˙(0) for which the corresponding solution r of (140) lies on the stable manifold of the flow, giving lim s →∞ r s ( ) = 3 m . All other values of ˙(0) yield trajecr tories which do not approach r = 3 m . In the non-static setting - that is, in the Class 1 McVittie spacetime - the stable manifold has an additional dimension. This corresponds to the flow of H t ( ) to the asymptotic value H 0 . However infinite fine tuning of the initial data is still required to ensure that the trajectory approaches r = 3 m : the stable manifold has co-dimension 1.

## 2. Particle orbits in Class 1

For timelike geodesics, we find that there are bound orbits that are asymptotic to the stable circular particle orbits of the corresponding Schwarzschild-de Sitter spacetime. By 'corresponding', we mean that the mass parameter m and limiting Hubble value H 0 of the McVittie spacetime provide the mass parameter and positive cosmological constant Λ = 3 H 2 0 of the Schwarzschild-de Sitter spacetime. We recall that r s ( ) = r c is a circular orbit of Schwarschild-de Sitter spacetime if G r ( c ) = 0, where (see (13) with glyph[epsilon1] = -1 and H t ( ) ≡ H 0 )

$$G ( r ) = H _ { 0 } ^ { 2 } r - \frac { m } { r ^ { 2 } } + ( 1 - \frac { 3 m } { r } ) \frac { \ell ^ { 2 } } { r ^ { 3 } }, \quad ( 1 4 1 ) \quad \text{$C$}$$

and that this orbit is stable if in addition

$$G ^ { \prime } ( r _ { c } ) < 0. \text{ \quad \ \ } ( 1 4 2 )$$

Stable circular orbits (or indeed circular orbits of any stability type) do not exist in Schwarzschild-de Sitter spacetime for all ranges of the parameters m and H 0 . We have the following result, which is proven in Appendix B.

Proposition 11 The Schwarzschild-de Sitter spacetime with mass parameter m and positive cosmological constant Λ = 3 H 2 0 admits stable circular particle orbits if and only if mH 0 &lt; 2 75 √ 3 . The radius of stable orbits satisfies r c ∈ (6 m,r + ) , where r + is the larger of the two positive roots of 1 -2 m/r -r 2 H 2 0 = 0 .

Comment 10 It is worth recalling what is meant by stability of these orbits. Let m H , 0 and glyph[lscript] be fixed. Write the timelike geodesic equations for Schwarzschild-de Sitter spacetime as the 2-dimensional dynamical system

$$\dot { r } = p, \quad \dot { p } = ( 1 - \frac { 3 m } { r } ) \frac { \ell ^ { 2 } } { r ^ { 3 } } - \frac { m } { r ^ { 2 } } + H _ { 0 } ^ { 2 } r. \quad ( 1 4 3 )$$

Let r c satisfy G r ( c ) = 0, so that ( r s , p ( ) ( s )) = ( r , c 0) is a circular orbit of the system. This orbit is said to be stable if for any glyph[epsilon1] &gt; 0 there exists δ &gt; 0 such that for all ( r , p 0 0 ) with | ( r , p 0 0 ) -( r , c 0) | &lt; δ , the solution of (143) with ( r (0) , p (0)) = ( r , p 0 0 ) - that is, the unique geodesic through ( r , p 0 0 ) -satisfies | ( r s , p ( ) ( s )) -( r , c 0) | &lt; glyph[epsilon1] for all s ≥ 0. By using a Lyapunov function argument, it is straightforward to see that the orbit is stable if G r ′ ( c ) &lt; 0.

Corresponding to Lemma 8, we have this result:

Lemma 9 Let r c &gt; 3 m satisfy G r ( c ) = 0 . In a neighbourhood of Q : x = ( H ,r , 0 c 0) , the equations governing timelike geodesics of a Class 1 McVittie spacetime may be written as (125), (126) with glyph[epsilon1] = -1 and (121) with

$$\cdot \quad u = u _ { 2 } \coloneqq - \frac { \alpha } { \chi } x _ { 3 } + \frac { 1 } { \chi } ( x _ { 3 } ^ { 2 } + \chi ( 1 + \frac { \ell ^ { 2 } } { x _ { 2 } ^ { 2 } } ) ) ^ { 1 / 2 }. \quad ( 1 4 4 )$$

Furthermore, Q is an equilibrium point and the hypersurface x 1 = H 0 is an invariant manifold of this system. The restriction of the dynamical system to x 1 = H 0 corresponds to the timelike geodesic equations of Schwarzschild-de Sitter spacetime.

Proof: Collecting relevant definitions establishes the first part of the proof. For the second, we use Lemma 5, and we note that a straightforward calculation yields

$$\dot { x _ { 3 } } | _ { Q } = G ( r _ { c } ). \text{ \quad \ \ } ( 1 4 5 )$$

We note also that

$$\stackrel { \text{w} } { \text{ar} } \\ 0, \quad \chi _ { Q } \coloneqq \chi | _ { r = r _ { c }, H = H _ { 0 } } = ( 1 - \frac { 3 m } { r _ { c } } ) ( 1 + \frac { \ell ^ { 2 } } { r _ { c } ^ { 2 } } ) > 0. \quad ( 1 4 6 )$$

glyph[square]

Corresponding to Proposition 10, we have this result which is established using centre manifold theory [17].

Proposition 12 In a Class 1 McVittie spacetime with line element (1), let r c satisfy G r ( c ) = 0 and G r ′ ( c ) &lt; 0 where G is defined in (141). Then Q : x = ( H ,r , 0 c 0) is a stable equilibrium point of the timelike geodesic equations. Furthermore, there exists a neighbourhood N Q ⊆ R 3 of Q : x = ( H ,r , 0 c 0) such that for any x (0) ∈ N Q , there exists a positive constant γ and a stable solution ¯( x s ) = ( r s , p ( ) ( s )) of (143) such that where κ is defined in part (v) of Definition 3, and we have

$$w ( x _ { 1 } ) < 0, \ \ \ f o r \ a l l \ x _ { 1 } > 0, \ \ \ \ ( 1 5 0 )$$

$$w ( - \xi ) = - w ( \xi ), \quad \text{for all } \xi > 0. \quad \ ( 1 5 1 )$$

glyph[square]

$$x _ { 1 } ( s ) & = O ( e ^ { - \gamma s } ), \quad ( x _ { 2 } ( s ), x _ { 3 } ( s ) ) = \bar { x } ( s ) + O ( e ^ { - \gamma s } ), \quad s \to \infty \\. \quad \, \tau \, \quad \tau \, \quad \tau \, \quad \tau \, \cdots \quad \tau \, \cdots \quad \tau \, \cdots \quad \tau \, \cdots \quad \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text}$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text`$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text#$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text'$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text.$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text<$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\textf$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text51$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{#\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$$

where ( x 1 ( s , x ) 2 ( s , x ) 3 ( s )) = x s ( ) is the unique timelike geodesic with initial data x 0 .

Proof: A straightforward calculation of the Jacobian matrix corresponding to the equilibrium point Q shows that it has one negative and two purely imagininary eigenvalues. Thus the equilibrium point has a onedimensional stable manifold and a two-dimensional centre manifold. In fact the invariant manifold x 1 = H 0 is a centre manifold of the point [17]. As noted in Lemma 9, the flow on the centre manifold is equivalent to the geodesic flow on Schwarzschild-de Sitter spacetime. Thus the equilibrium point is a stable point of the centre manifold. It follows by Theorem 2 of [17] and by the comments above that the equilibrium point is a stable point of the full flow (i.e. not restricted to the centre manifold), and that the asymptotic behaviour of the flow is as described in the statement of the proposition.

Comment 11 The interpretation of Proposition 12 is as follows. The result tells us that for sufficiently late times ( H t ( ) close to H 0 ), timelike geodesics which are sufficiently close to r = r c with sufficiently small radial velocity ( ˙ close to 0), the orbit is equivalent to the correspondr ing orbit in Schwarzschild-de Sitter spacetime, modulo an exponentially decaying correction. These orbits are future complete, and correspond to bound particle orbits of the spacetime. Their past evolution is governed by Proposition 3. Note that to obtain geodesics of the spacetime, we must choose x (0) ∈ N Q with x 1 (0) &gt; H 0 (i.e. H t ( (0)) &gt; H 0 ).

## B. Bound orbits in Class 2 McVittie spacetimes

We now consider Class 2 McVittie spacetimes as defined by Definitions 1 and 3. Again, we begin by writing down the result that gives a key geodesic equation in dynamical systems form.

Lemma 10 The geodesic equation for x 1 can be written as

$$\dot { x _ { 1 } } = w ( x _ { 1 } ) u ( x ), \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where w and u are as defined in (79), (113) and (118). The function w extends to a function w ∈ C 2 ( R R , ) satisfying

$$( w ( 0 ), w ^ { \prime } ( 0 ), w ^ { \prime \prime } ( 0 ) ) = ( 0, 0, - 2 \kappa ), \quad ( 1 4 9 ) \quad \text{the} ^ { ^ { * } }$$

As in the Class 1 case, the other geodesic equations are given by (125) and (126).

## 1. Photon orbits in Class 2

The analysis above for photon orbits in Class 1 McVittie spacetimes follows through to the present case except in one important regard. The equilibrium point P is now given by x P ( ) = (0 3 , m, 0), and the Jacobian matrix at P is given by

$$J ( P ) = \left ( \begin{array} { c c c } 0 & 0 & 0 \\ 0 & 0 & 1 \\ 0 & c & 0 \end{array} \right ), \quad c = \frac { \ell ^ { 2 } } { ( 3 m ) ^ { 4 } }. \quad \ ( 1 5 2 )$$

Thus the eigenvalues are 0 , ± √ c , and the equilibrium point in non-hyperbolic [16]. (This follows from the fact that w ′ (0) = 0, which is one of the features that distinguishes between Class 1 and Class 2.) Thus we must appeal to centre manifold theory [17]. For the linearized system at P , the centre space is given by E c = { x : x 2 = x 3 = 0 . } By Theorem 1 of [17], it follows that there is a local centre manifold W c of the nonlinear system which takes the form

$$\colon & \colon \text{$\alpha$} \\ \text{le or} \quad x _ { 2 } = h _ { 2 } ( x _ { 1 } ), \quad x _ { 3 } = h _ { 3 } ( x _ { 1 } ), \quad h _ { i } ( x _ { 1 } ) = O ( x _ { 1 } ^ { 2 } ), \quad x _ { 1 } \to 0, \\ \text{ermed} \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdots$$

where h i is C 2 in a neighbourhood of x 1 = 0. The dynamics on the local centre manifold are governed by

$$\dot { x _ { 1 } } = w ( x _ { 1 } ) u _ { 1 } ( x _ { 1 }, h _ { 2 } ( x _ { 1 } ), h _ { 3 } ( x _ { 1 } ) ). \quad ( 1 5 4 )$$

From Lemma 10 and (153), we see that this reads

$$\dot { x _ { 1 } } = - \kappa x _ { 1 } ^ { 2 } u _ { 1 } ( P ) + O ( x _ { 1 } ^ { 3 } ), \ \ x _ { 1 } \to 0. \quad \ ( 1 5 5 )$$

We note that (128) gives u 1 ( P ) = glyph[lscript] √ 3 m &gt; 0. It follows that if x 1 (0) is sufficently small and positive, and if x (0) = ( x 1 (0) , h 2 ( x 1 (0)) , h 3 ( x 1 (0))) (that is, if the initial data for the geodesic correspond to a point of the centre manifold W c ), then x 1 ( s ) → 0 on the centre manifold as s →∞ .

In addition, there are local stable and unstable manifolds of the flow, W s and W u respectively. These are tangent at the equilibrium point to the stable and unstable spaces of the linearised system. The overall stability properties of the equilibrium point may be deduced from the topological equivalence of the nonlinear flow to the

flow described by the partitioned system obtained by restricting to the local centre, stable and unstable manifolds respectively (see e.g. Theorem 2.7.2 of [18]). Each of these is one dimensional, and given the description above of the stability properties of the flow on the local centre manifold, we conclude that we have the same result as applies in Class 1. Comment 8 following Proposition 10 also holds in the present case. We note that in the language of [19], there is a two-dimensional local centre-stable manifold of the flow. This plays essentially the same role as the 2-dimensional stable manifold of Proposition 10. We summarise as follows.

Proposition 13 In a Class 2 McVittie spacetime, there exists a 3-parameter family of null geodesics that are future complete, with r → 3 m and t → + ∞ as s → ∞ along each member of the family. These geodesics correspond to bound photon orbits of the spacetime. glyph[square]

## 2. Particle orbits in Class 2

Turning to the case of timelike geodesics in Class 2 McVittie spacetimes, we encounter a technical difficulty which is not resolved by general dynamical systems theory: all of the eigenvalues of the Jacobian matrix at the equilibrium point have zero real parts. However, a useful feature of the dynamical system allows us to circumvent this problem and show that the stability properties of timelike geodesics in Class 1, as described in Proposition 12, carry over to Class 2. This feature is that we can think of the dynamical system as being a perturbation of the system describing the timelike geodesics of Schwarzschild spacetime. This allows us to exploit the existence of a conserved energy in this latter system. This energy (or more accurately, the effective potential, which plays the role of a Lyapunov function) is of course no longer conserved in the non-static McVittie spacetime, but (i) it remains positive definite and (ii) its evolution can be controlled by a Gronwall type argument. These properties allow us to conclude that the equilibrium point x = (0 , r c , 0) is stable.

We begin by recalling the form of the relevant geodesic equations in a neighbourhood of x = (0 , r c , 0), where r c corresponds to a circular particle orbit of the Schwarzschild spacetime with mass parameter m .

## Lemma 11 Define

$$V ( r ) = - \frac { m } { r } + \frac { \ell ^ { 2 } } { 2 r ^ { 2 } } - \frac { m \ell ^ { 2 } } { r ^ { 3 } }, \quad \ \ ( 1 5 6 ) \quad \text{Ler} \\ \text{sol}$$

and let r c &gt; 3 m satisfy V ′ ( r c ) = 0 . Then in a neighbourhood of Q : x = (0 , r c , 0) , the equations governing timelike geodesics of a Class 2 McVittie spacetime may be written as

$$x _ { 1 } & = w ( x _ { 1 } ) u _ { 2 } ( x ), & ( 1 5 7 )$$

$$\dot { x _ { 2 } } & = x _ { 3 }, & ( 1 5 8 )$$

$$\dot { x _ { 3 } } = - V ^ { \prime } ( x _ { 2 } ) + x _ { 2 } x _ { 1 } ^ { 2 } + x _ { 2 } f ^ { 1 / 2 } ( x _ { 2 } ) w ( x _ { 1 } ) u _ { 2 } ( x ) ^ { 2 } ( 1 5 9 )$$

where u 2 is defined in (144) and w satisfies the conditions of Lemma 10. Furthermore, Q is an equilibrium point and the hypersurface x 1 = 0 is an invariant manifold of this system. The restriction of the dynamical system to x 1 = 0 corresponds to the timelike geodesic equations of Schwarzschild spacetime. glyph[square]

Comment 12 With V ′ ( r c ) = 0, we see that r = r c corresponds to a circular orbit of the Schwarzschild spacetime obtained by setting H t ( ) ≡ 0 in (1). It is worth considering briefly how stability of such orbits may be deduced. We define the Lyapunov function

$$W ( x _ { 2 }, x _ { 3 } ) = \frac { 1 } { 2 } x _ { 3 } ^ { 2 } + V ( x _ { 2 } ) - V ( r _ { c } ). \quad ( 1 6 0 )$$

If x 2 = r c is a local minimum of V (which holds if V ′ ( r c ) = 0 , V ′′ ( r c ) &gt; 0), then there exists glyph[epsilon1] &gt; 0 such that

$$| ( x _ { 2 } - r _ { c }, x _ { 3 } ) | < \epsilon \Rightarrow \begin{cases} \begin{array} { c } W ( x _ { 2 }, x _ { 3 } ) \geq 0, \\ W ( x _ { 2 }, x _ { 3 } ) = 0 \Leftrightarrow ( x _ { 2 }, x _ { 3 } ) = ( r _ { c }, 0 ). \end{array} \\ \text{class 2} \quad \text{$\mathfrak{D.\sim\lim_{1}$ and\lim_{1}$} \end{cases}$$

By a simple calculation,

$$\dot { W } = 0, \quad \quad \quad ( 1 6 2 )$$

and so by Lyapunov's stability theorem (see e.g. [16], p. 130), ( x , x 2 3 ) = ( r , c 0) is a stable equilibrium point: given glyph[epsilon1] &gt; 0 and sufficiently small, there exists δ &gt; 0 such that initial data with

$$| ( x _ { 2 } ( 0 ), x _ { 3 } ( 0 ) ) - ( r _ { c }, 0 ) | < \delta \text{ \quad \ \ } ( 1 6 3 )$$

yield solutions with

$$\S _ { L } ^ { S } \quad | ( x _ { 2 } ( s ), x _ { 3 } ( s ) ) - ( r _ { c }, 0 ) | < \epsilon \quad \text{for all } s \geq 0. \quad ( 1 6 4 )$$

By (162), the flow is confined to the level set

$$W ( x _ { 2 }, x _ { 3 } ) = W ( x _ { 2 } ( 0 ), x _ { 3 } ( 0 ) ). \quad \ \ ( 1 6 5 )$$

We note that W is essentially (the square of) the conserved energy associated with the timelike Killing vector of the exterior Schwarzschild spacetime.

Moving away from the invariant manifold x 1 = 0, the flow no longer has a conserved energy. However, we can use the same function W to obtain useful control over the flow. The overall strategy is to introduce a mollified dynamical system that is equivalent to the system of Lemma 11 in a neighbourhood of Q . We use the properties of a mollified form of W to prove global existence of solutions and stability of Q in this system. Stability in the full system then follows as the systems coincide in a neighbourhood of Q .

We begin by noting that the geodesic equations of Lemma 11 may be written as

$$\dot { x _ { 1 } } = - k _ { 1 } x _ { 1 } ^ { 2 } A ( x ),$$

$$\dot { x _ { 2 } } & = x _ { 3 }, & ( 1 6 7 )$$

$$\dot { x _ { 3 } } = - V ^ { \prime } ( x _ { 2 } ) + k _ { 2 } x _ { 1 } ^ { 2 } B ( x ), \text{ \quad \ \ } ( 1 6 8 )$$

where k , k 1 2 are positive constants, A,B are C 2 in a neighbourhood of Q , and A Q ( ) = B Q ( ) = 1.

We introduce a mollifier ψ ∈ C ∞ 0 ( R 3 , R ) with the property that

$$\psi ( x ) = \begin{cases} \ 1 & \text{if} \quad | x - x ( Q ) | < 1, \\ 0 & \text{if} \quad | x - x ( Q ) | \geq 2, \end{cases} \quad ( 1 6 9 ) \quad \text{while}$$

and we introduce mollified versions of the coefficients of the dynamical system:

$$A _ { \epsilon } ( x ) = A ( x \psi ( \frac { x } { \epsilon } ) ), \ \ B _ { \epsilon } ( x ) = B ( x \psi ( \frac { x } { \epsilon } ) ). \quad ( 1 7 0 ) \quad \stackrel { W _ { \epsilon } ( s ) } { \text{for $A$} }$$

Then A glyph[epsilon1] and B glyph[epsilon1] coincide with A and B respectively for | x -x Q ( ) | &lt; glyph[epsilon1] , and we have A glyph[epsilon1] ( x ) = B glyph[epsilon1] ( x ) = 1 for | x -x Q ( ) | &gt; 2 . glyph[epsilon1] We also introduce a mollified form of the potential V . We introduce a smooth, bounded function V glyph[epsilon1] ∈ C ∞ 0 ( R R , ) with

$$V _ { \epsilon } ( x _ { 2 } ) = \begin{cases} \ V ( x _ { 2 } ), & | x _ { 2 } - r _ { c } | < \epsilon, \\ c, & | x _ { 2 } - r _ { c } | > 2 \epsilon, \end{cases} \quad ( 1 7 1 ) \quad.$$

with c &gt; V ( r c ). It follows that V glyph[epsilon1] -V ( r c ) is smooth, bounded and nonnegative for all x 2 , coincides with V ( x 2 ) -V ( r c ) for | x 2 -r c | &lt; glyph[epsilon1] and vanishes if and only if x 2 = r c . It follows that solutions of the system

$$\dot { x _ { 1 } } = - k _ { 1 } x _ { 1 } ^ { 2 } A _ { \epsilon } ( x ), \quad \quad ( 1 7 2 ) \quad \text{Int} \,,$$

$$\dot { x _ { 2 } } & = x _ { 3 }, & ( 1 7 3 )$$

$$\dot { x _ { 3 } } = - V _ { \epsilon } ^ { \prime } ( x _ { 2 } ) + k _ { 2 } x _ { 1 } ^ { 2 } B _ { \epsilon } ( x ) \text{ \quad \ \ } ( 1 7 4 )$$

coincide with solutions of the original system (166)-(168) when those solutions satisfy | x -x Q ( ) | &lt; glyph[epsilon1] .

Lemma 12 Let x s ( ) be a solution of (172)-(174) with | x (0) -x Q ( ) | &lt; glyph[epsilon1] and x 1 (0) &gt; 0 . Then there exists δ ∈ (0 1) , such that

$$0 < x _ { 1 } ( s ) < \frac { x _ { 1 } ( 0 ) } { 1 + k _ { 1 } ( 1 - \delta ) x _ { 1 } ( 0 ) s } \ \text{ for all } s > 0. \ ( 1 7 5 ) \ \underset { \text{that } 1 } { \text{ since } 1 }$$

Proof: From the properties of A glyph[epsilon1] , we can choose glyph[epsilon1] &gt; 0 so that there exists δ ∈ (0 1) such that ,

$$1 - \delta < A _ { \epsilon } ( x ) < 1 + \delta \quad \text{for all } x \in \mathbb { R } ^ { 3 }. \quad ( 1 7 6 ) \quad ^ { \iota \mu \nu \nu } _ { c o n s t }$$

Then (172) yields

$$\frac { d } { d s } \left ( \frac { 1 } { x _ { 1 } ( s ) } \right ) = k _ { 1 } A _ { \epsilon } ( x ( s ) ) > k _ { 1 } ( 1 - \delta ), \quad ( 1 7 7 )$$

and integrating yields the result.

glyph[square]

Corollary 2 Let x s ( ) be a solution of (172)-(174) with | x (0) -x Q ( ) | &lt; glyph[epsilon1] and x 1 (0) &gt; 0 . Then x 1 ( s ) &lt; x 1 (0) and | sx 1 ( s ) | &lt; k 1 (1 -δ ) for all s &gt; 0 . | square

Now we define

$$W _ { \epsilon } ( x _ { 2 }, x _ { 3 } ) = \frac { 1 } { 2 } x _ { 3 } ^ { 2 } + V _ { \epsilon } ( x _ { 2 } ) - V ( r _ { c } ). \quad \ ( 1 7 8 ) \quad \ \ \ s m a l { \prod \nolimits }$$

Then using (173) and (174), we have

$$\frac { d W _ { \epsilon } } { d s } = k _ { 2 } B _ { \epsilon } ( x ) x _ { 1 } ^ { 2 } x _ { 3 }, \text{ \quad \ \ } ( 1 7 9 )$$

which has the formal solution

$$\underset { 3 \text{ of} } { W _ { \epsilon } ( s ) } = W _ { \epsilon } ( 0 ) + \int _ { 0 } ^ { s } k _ { 2 } B _ { \epsilon } ( x ( \sigma ) ) x _ { 1 } ^ { 2 } ( \sigma ) x _ { 3 } ( \sigma ) d \sigma, \ \ ( 1 8 0 )$$

where for convenience we write W x glyph[epsilon1] ( 2 ( s , x ) 3 ( s )) = W s glyph[epsilon1] ( ). Introducing the same bound for B glyph[epsilon1] as we did for A glyph[epsilon1] we can write

$$\inf _ { \substack { \epsilon \\ \epsilon \colon \_ \_ } } \left | W _ { \epsilon } ( s ) \right | \leq | W _ { \epsilon } ( 0 ) | + k _ { 2 } ( 1 + \delta ) \int _ { 0 } ^ { s } x _ { 1 } ^ { 2 } ( \sigma ) | x _ { 3 } ( \sigma ) | d \sigma. \ ( 1 8 1 )$$

Since V glyph[epsilon1] ( x 2 ) -V ( r c ) ≥ 0 for all x 2 ∈ R (see (171)), we have

$$\frac { 1 } { 2 } x _ { 3 } ^ { 2 } ( s ) \leq W _ { \epsilon } ( s ) \leq Z ( s ), \text{ \quad \ \ } ( 1 8 2 )$$

where Z s ( ) is the non-negative function on the right hand side of (181). Then

$$Z ^ { \prime } ( s ) & = k _ { 2 } ( 1 + \delta ) x _ { 1 } ^ { 2 } ( s ) | x _ { 3 } ( s ) | \\ & \leq \sqrt { 2 } k _ { 1 } ( 1 + \delta ) x _ { 1 } ^ { 2 } ( s ) ( Z ( s ) ) ^ { 1 / 2 }. \quad ( 1 8 3 )$$

Integrating and using Lemma 12 yields

$$( Z ( s ) ) ^ { 1 / 2 } & \leq | W _ { \epsilon } ( 0 ) | ^ { 1 / 2 } + \frac { 1 } { \sqrt { 2 } } \frac { k _ { 2 } ( 1 + \delta ) x _ { 1 } ^ { 2 } ( 0 ) s } { 1 + k _ { 1 } ( 1 - \delta ) x _ { 1 } ( 0 ) s } \\ & \leq | W _ { \epsilon } ( 0 ) | ^ { 1 / 2 } + c _ { 1 } x _ { 1 } ( 0 ).$$

where

$$c _ { 1 } = \frac { k _ { 1 } ( 1 + \delta ) } { \sqrt { 2 } k _ { 2 } ( 1 - \delta ) } > 0. \text{ \quad \ \ } ( 1 8 5 )$$

Since W s glyph[epsilon1] ( ) is positive semi-definite, and using the fact that r = r c is a local minimum of the potential, this yields the following a priori bounds:

Lemma 13 Let x s ( ) be a solution of (172)-(174) with | x (0) -x Q ( ) | &lt; glyph[epsilon1] and x 1 (0) &gt; 0 . Then there exist positive constants c 1 , c 2 and c 3 such that

$$| x _ { 2 } ( s ) | \leq c _ { 2 } ( | W _ { \epsilon } ( 0 ) | ^ { 1 / 2 } + c _ { 1 } x _ { 1 } ( 0 ) ), \quad s \geq 0, ( 1 8 6 )$$

$$| x _ { 3 } ( s ) | \leq c _ { 3 } ( | W _ { \epsilon } ( 0 ) | ^ { 1 / 2 } + c _ { 1 } x _ { 1 } ( 0 ) ), \quad s \geq 0. \ ( 1 8 7 )$$

glyph[square]

Thus we have shown that if we choose initial data with | x (0) -x Q ( ) | &lt; glyph[epsilon1] , then Lemmas 12 and 13 show that solutions of the mollified system satisfy | x s ( ) -x Q ( ) | &lt; O glyph[epsilon1] ( ) for all s &gt; 0. In conjunction with Theorem 1 above, this is sufficient to prove global existence for solutions of the mollified system. Since solutions of the original system coincide with solutions of the mollified system for small | x -x Q ( ) , | this proves stability of the equilibrium point of the full system. We summarise as follows.

Proposition 14 In a Class 2 McVittie spacetime with line element (1), let r c satisfy V ′ ( r c ) = 0 and V ′′ ( r c ) &gt; 0 where V is defined in (156). Let x Q ( ) = (0 , r c , 0) . There exists glyph[epsilon1] &gt; 0 such that, if | x (0) -x Q ( ) | &lt; glyph[epsilon1] , with x 1 (0) &gt; 0 , then there exists a global, unique solution of the timelike geodesic equations (157)-(159) with initial data x (0) . | x s ( ) -x Q ( ) | &lt; O glyph[epsilon1] ( ) for all s &gt; 0 , and lim s → ∞ + x 1 ( s ) = 0 . glyph[square]

Comment 13 The geodesics described by this proposition are bound particle orbits of the Class 2 McVittie spacetime in question. The result shows that if a particle orbit, at sufficiently late times, is, at any instant, sufficiently close in ( r -˙) r phase space to a stable circular orbit of the corresponding Schwarzschild spacetime, then the orbit remains close for all future times.

## VII. CONCLUSIONS

McVittie spacetimes provide us with a conceptually very simple model of the gravitational field of a point mass embedded in an otherwise isotropic, spatially flat universe. By conceptual simplicity, we mean that we can write down simple conditions on a spacetime that lead uniquely to the McVittie metric (1). See for example [8] and [9]. The interpretation of the resulting spacetime is somewhat more complex. However, a clear picture has been established by this stage. In the case of an expanding spacetime, the hypersurface ∂ Ω 2 m = { ( t, r ) : t &gt; 0 , r = 2 m } forms the past boundary of the spacetime: all causal geodesics originate at this singular surface at a finite time in the past. Outgoing rays (outgoing radial null geodesics) escape to infinity, and there are families of ingoing rays (ingoing radial null geodesics) that meet either a horizon (Class 1) or a singularity (Class 2) at a finite time in their future. In addition, we have shown that large, physically relevant classes of McVittie spacetimes can 'capture' both photons and massive particles, and keep them in orbit at finite radius for all time. This reinforces the interpretation of McVittie spacetimes as representing a central massive object in an otherwise isotropic universe. One can interpret all the results on test-particle (massive and massless) motion in McVittie spacetimes in terms of a balance between the attraction of the central mass and the expansion of the cosmological background. The results presented here show that there exists a region of the spacetime where the central mass always wins out, subject to initial conditions on the motion of the test particle - as in the vacuum black hole case (with or without a cosmological constant).

As noted, these results also lay the foundations for the study of accretion disks in McVittie spacetimes. We have essentially established the existence of ISCOs in these spacetimes, and our results also make it clear that significant gravitational lensing occurs. We speculate that this may provide a useful setting for the calculation of lensing effects by galactic clusters or other over-densities in the universe. Likewise, it is clear that the central mass of a McVittie metric significantly disturbs and traps a proportion of a cloud of photons that propagates through a region of spacetime sufficiently close to that central mass. It would be of interest to determine the effect of this central mass on CMB calculations, e.g. through the recalculation of distance relations [6]. We note in this context that while the McVittie spacetime in spatially homogenous in the sense that ρ = ρ t ( ), effects of inhomogeneity would be manifest in the angular diameter distance D A , which satisfies

$$\frac { d ^ { 2 } D _ { A } } { d s ^ { 2 } } = - ( | \sigma | ^ { 2 } + 4 \pi ( \rho + p ) ( u _ { a } k ^ { a } ) ^ { 2 } ) D _ { A }, \quad ( 1 8 8 )$$

where u a is the fluid 4-velocity, k a is tangent to a bundle of null geodesics (with shear σ ) emanating from the source and the deriative is along k a . In this and other ways, we believe that McVittie spacetimes could play an interesting role in inhomogeneous cosmological modelling, by providing a conceptually simple model of an overdensity in an otherwise spatially homogeneous universe.

## ACKNOWLEDGMENTS

I thank Kayll Lake for interesting discussions of McVittie spacetimes.

## Appendix A: Global structure of McVittie spacetimes admitting a circular photon orbit.

For convenience, we recall the following facts. A McVittie spacetime has line element

$$\begin{smallmatrix} ( \text{Class} & d s ^ { 2 } = - ( f - r ^ { 2 } H ^ { 2 } ) d t ^ { 2 } - 2 r H f ^ { - 1 / 2 } d t d r + f ^ { - 1 } d r ^ { 2 } + r ^ { 2 } d \Omega ^ { 2 }, \\ \text{have} & &. &. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &... &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. &.. \end{smallmatrix}$$

where H = H t ( ) is the Hubble function of the background and f = 1 -2 m/r . As proven in Section 3 above, the spacetime with this line element admits a circular photon orbit if and only if

$$H ( t ) = - H _ { 0 } \tanh ( A H _ { 0 } t ), \ \ t \in \mathbb { R }. \quad ( A 2 )$$

where

$$A = \frac { ( 1 - \frac { 3 m } { a } ) } { \sqrt { 1 - \frac { 2 m } { a } } }, \ \ H _ { 0 } = a ^ { - 1 } \sqrt { 1 - \frac { 2 m } { a } } \quad ( A 3 )$$

and the orbital radius is r = a &gt; 2 m . The scale factor is given by

$$S ( t ) = ( \cosh ( A H _ { 0 } t ) ) ^ { - 1 / A }, \text{ \quad \ } ( A 4 )$$

which immediately shows that an important role is played by the sign of A . See Figures 1 and 2 in Section 3 above. We exclude the case a = 3 m (which corresponds to Schwarzschild-de Sitter spacetime) and note that A &lt; 0 ⇔ a &lt; 3 m .

Defining z = √ 1 -2 m/r (cf. (29)), the equations governing radial null geodesics are

$$\frac { d z } { d t } = \frac { H } { 2 } ( 1 - z ^ { 2 } ) \pm \frac { z } { 4 m } ( 1 - z ^ { 2 } ) ^ { 2 }. \quad \ ( A 5 )$$

Compactified representations of the outgoing (upper sign) and ingoing (lower sign) RNGs are shown in Figures 4 and 5. These diagrams also include the horizons of the spacetimes, given by

$$H ^ { 2 } ( t ) = \frac { z ^ { 2 } } { 4 m ^ { 2 } } ( 1 - z ^ { 2 } ) ^ { 2 }. \quad \quad ( \text{A6} ) \quad \text{Fig}$$

## 1. CPO at radius a &lt; 3 m .

In Figure 4, we see that outgoing geodesics originate in the past at either (i) ( z, t ) = (0 , t 0 ) with t 0 &gt; 0; (ii) ( z, t ) = ( z -, -∞ ); (iii) ( z, t ) = (+ ∞ -∞ , ) or (iv) ( z, t ) = ( z + , -∞ ). Note that z ± are the larger and smaller of the two roots in (0 , 1) of the asymptotic horizon equation

$$H _ { 0 } ^ { 2 } = \frac { z ^ { 2 } } { 4 m ^ { 2 } } ( 1 - z ^ { 2 } ) ^ { 2 }. \quad \quad ( \text{A7} ) \quad \begin{array} { c } \frac { \epsilon _ { 1 } } { ( \text{A} } \\ + \infty \end{array} }$$

There is a unique outgoing radial null geodesic (ORNG) in case (iv), which is a separatrix for cases (ii) and (iii). The existence of this geodesic may be established using dynamical systems, where ( z, t ) = ( z + , -∞ ) corresponds to a saddle. The ORNGs are either (v) future-complete, terminating at ( z, t ) = (+ ∞ , + ∞ ) or (vi) of finite affine length, terminate at a finite time ( t &lt; 0) in the future at z = 0 ( r = 2 m ). The symmetry H ( -t ) = -H t ( ) and the form of the RNG equations (A5) shows that the ingoing RNGs are simply the time-reversed outoing RNGs (as is manifest in Figure 4). By appealing to the geodesic equation (see (13) with glyph[lscript] = glyph[epsilon1] = 0)

$$\ddot { r } = r f ^ { 1 / 2 } H ^ { \prime } ( t ) \dot { t } ^ { 2 }, \text{ \quad \ \ } ( A 8 ) \text{ \ \ IR }$$

it is readily established that the boundary r = 2 m is at finite affine distance in cases (i) and (v). Likewise, z = z -is at a finite affine distance to the past in case (ii), and geodesics extending to r = + ∞ (cases (iii) and (iv)) have infinite affine length in the relevant temporal direction.

It is straightforward to establish the following facts in relation to the corresponding background spacetime that is, the isotropic spacetime with line element (3) with H t ( ) given by (A2).

- · The comoving world-lines x = x 0 are complete timelike geodesics, and t is proper time along these geodesics. (Recall that x = rS -1 ( ) where t S is the scale factor.)
- · Outgoing RNGs are future complete and ingoing RNGs are past complete; r and | t | become infinite in the limit of infinite affine length.
- · Outgoing RNGs either meet the regular origin at finite affine time in the past, or are past complete.
- · Ingoing RNGs either meet the regulat origin at finite affine time in the future, or are future complete.
- · The comoving coordinate x has a finite limit as the affine parameter | s | →∞ along null geodesics.

These points give rise to the conformal diagrams of Figure 5 for the McVittie spacetime and its background.

## 2. CPO at radius a &gt; 3 m .

As we see in Figure 6, all ORNGs originate at z = 0 ( r = 2 m ) at some time t &lt; 0. There is a unique ORNG which terminates at ( z, t ) = ( z -, + ∞ ), and all others terminate at z = 0 , t &gt; 0 or at ( z, t ) = ( z + , + ∞ ). The IRNGs are the time reversal of the ORNGs. Noting the exponential decay of S t ( ) as | t | → ∞ in this case (cf. (A4)), we find that the boundaries at r = 2 m and t = + ∞ are at finite affine distance.

The following properties hold in the background:

- · The comoving worldlines are complete timelike geodesics carrying proper time t , and r → 0 as t →±∞ along these geodesics.
- · Along RNGs, t becomes infinite in finite affine time, and r remains finite in this limit.
- · All ORNGs meet r = 0 at a finite t and at finite affine distance in the future; all IRNGs meet r = 0 in finite t and finite affine distance in the past.

The behaviour of the ORNGs as t → + ∞ and the IRNGs as t → -∞ suggests that the spacetime is extendible across these coordinate boundaries. We note that the scale factor is asymptotic to the scale factor of a de Sitter spacetime:

$$\text{and} \quad S ( t ) = ( \cosh ( A H _ { 0 } t ) ) ^ { - 1 / A } \sim 2 ^ { 1 / A } e ^ { - H _ { 0 } t }, \quad t \to + \infty, \\ \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

and so extending to de Sitter spacetime seems to be a natural possibility. We note that the time coordinate t here shares features of the time coordinate ˆ in t the representation

$$\L e ^ { \quad d s ^ { 2 } = - d \hat { t } ^ { 2 } + \exp ( 2 \alpha ^ { - 1 } \hat { t } ) ( d \hat { x } ^ { 2 } + d \hat { y } ^ { 2 } + d \hat { z } ^ { 2 } ) } \quad ( \mathbf A 1 0 )$$

of the de Sitter metric. These coordinates cover only half of the full de Sitter hyperboloid. See Section 5.2 of [20]. We conjecture that a similar extension to Schwarzschildde Sitter spacetime is possible in the McVittie case.

The conformal diagrams that ensue from these features are shown in Figure 7.

-

-

-

0.15

0.10

0.05

0.00

0.05

0.10

0.15

0.0

0.2

0.4

0.6

0.8

1.0

FIG. 4. Radial null geodesics and the horizon (red) of a McVittie spacetime admitting a circular photon orbit at radius a &lt; 3 m . Here m = 1 and a = 2 5, . giving A = -0 447214 . &lt; 0 and H 0 = 0 178885. . Outgoing geodesics are shown in the upper diagram, IRNGs in the lower. Due to time reversal properties of the spacetime, IRNGs are simply time reversed ORNGs: the corresponding diagram is obtained by inversion. In the diagram, the radial coordinate is z = √ 1 -2 m/r ∈ (0 1) and runs along the horizontal axis. , Thus r = 2 m at z = 0 and r →∞ at z = 1. The vertical axis corresponds to the value of H t ( ) along the geodesics.Referring to Figure 1, we note that t increases from bottom to top, giving the direction of flow along the geodesics. The critical geodesics ¯ γ O and ¯ γ I are shown in black: these are the unique ORNG/IRNG that originate/terminate at r = r + . The trapping properties of the different regions of the spacetime may be read from this diagram: for t &lt; 0, there are trapped regions (both ORNGs and IRNGs move to the left - decreasing r - with increasing t ) and regular regions (ORNGs move right, IRNGs move left). For t &gt; 0, there are anti-trapped regions (both families move right).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0044v1-figure-65469f7f-baba-442a-a403-7ae6e79b1f08.jpg)

## Appendix B: Stable circular orbits in Schwarzschild-de Sitter spacetime.

The conditions for a stable circular (particle) orbit in the Schwarzschild-de Sitter spacetime with mass parameter m and cosmological constant Λ = 3 H 2 0 are G r ( c ) = 0 , G ′ ( r c ) &lt; 0 where r c is the orbital radius and G is given in (141). We wish to answer two related questions: for what range of the parameters m,H 0

FIG. 5. Conformal diagrams for the McVittie spacetime (bottom) that admits a circular photon orbit at r = a &lt; 3 m , and the corresponding isotropic background (top). In both cases, t increases in the vertical direction and J ± corresponds to r = + ∞ , t = ±∞ . The spacetime is collapsing ( H &lt; 0) for t &lt; 0 and expanding ( H &gt; 0) for t &gt; 0. In the McVittie case, we note the existence of ORNGs that (a) originate at r = r -and terminate at r = 2 m ; (b) originate at r = r -and extend to r = + ∞ ; (c) originate at r = + ∞ , t = -∞ and extend to r = + ∞ , t = + ∞ and (d) originate at r = 2 m and extend to r = + ∞ , t = + ∞ . Similar results hold for IRNGs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0044v1-figure-ce975be6-c95b-48fc-979c-37169f355a01.jpg)

do stable circular orbits exist? What is the radius of the innermost stable circular orbit (ISCO)? The algebra is simplified by introducing variables u &gt; 2 , v &gt; 0 and a parameter h 0 &gt; 0 defined by r = mu,glyph[lscript] 2 = m v 2 and mH 0 = h 0 . Note that G r ( c ) = 0 , G ′ ( r c ) &lt; 0 if and only if ¯ ( G r c ) = 0 , G ¯ ( ′ r c ) &lt; 0 where ¯ ( G r ) = r 4 G r ( ) (recall that r &gt; 2 m ). This allows us to deal with polynomial functions only. Then the condition for a circular orbit is

$$\alpha ( u, v ) \coloneqq h _ { 0 } ^ { 2 } u ^ { 5 } - u ^ { 2 } + ( u - 3 ) v = 0, \quad ( \mathbf B 1 )$$

and the condition for stability of this orbit is

$$\beta ( u, v ) \coloneqq 3 u ^ { 2 } - 4 u v + 1 5 v < 0. \quad \text{(B2)}$$

Note that we have used (B1) to eliminate the term in h 2 0 that naturally arises in β u, v ( ).

Considering β as a parametrised quadratic in u shows that we must have v &gt; 45 4 and that (B2) holds if and only if

$$u _ { - } ( v ) < u < u _ { + } ( v ) \text{ \quad \quad } \text{(B3)}$$

FIG. 6. Radial null geodesics of a McVittie spacetime admitting a circular photon orbit at radius a &gt; 3 m . Here m = 1 and a = 3 5,giving . A = 0 218218 . &gt; 0 and H 0 = 0 187044. . Outgoing geodesics are shown in the top panel, and ingoing in the bottom. In the diagram, the radial coordinate is z = √ 1 -2 m/r ∈ (0 , 1) (see (29) below) and runs along the horizontal axis. Thus r = 2 m at z = 0 and r →∞ at z = 1. The vertical axis corresponds to the value of -H t ( ) along the geodesics. Referring to Figure 2, we note that t increases from bottom to top, giving the direction of flow along the geodesics. The expanding region of the spacetime ( t &lt; 0) contains regular and anti-trapped regions, and the collapsing region ( t &gt; 0) contains regular and trapped regions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0044v1-figure-7f2303d4-a65f-4a90-83a2-1e5c7a520741.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0044v1-figure-6ca011b6-7a95-4dd6-9b52-cfadd994c67a.jpg)

where

$$u _ { \pm } ( v ) = \frac { 2 v } { 3 } \left ( 1 \pm \sqrt { 1 - \frac { 4 5 } { 4 v } } \right ). \quad \quad ( B 4 )$$

We note that u ′ -( v ) &lt; 0 for v &gt; 45 4 and that lim v → ∞ + u -( v ) = 15 4 . This yields the (non-sharp) lower bound u &gt; 15 4 for the radius of a stable circular orbit. We also have the upper bound u -( v ) ≤ u -( 45 4 ) = 15 2 .

Suppose that for a given value of h 0 , there exists a point ( u , v 1 1 ) such that α u , v ( 1 1 ) = 0 and β u , v ( 1 1 ) &lt; 0 -that is, there is a stable circular orbit at u = u 1 with angular momentum parameter v 1 . Define

$$g ( u, v ) = \frac { u ^ { 2 } - ( u - 3 ) v } { u ^ { 5 } } \, \text{ \quad \ \ } \text{(B5)}$$

so that g u, v ( ) = h 2 0 if and only if α u, v ( ) = 0. Along the curve g u, v ( ) = h 2 0 we have have dv du &gt; 0 for β u, v ( ) &lt; 0. It follows that the minimum of u on g u, v ( ) = h 2 0 for which β &lt; 0 occurs along u = u -( v ). Furthermore any value of h 0 for which there exists points ( u, v ) with g u, v ( ) = h 2 0 gives rise to a curve with this property.

Substituting then yields

$$h _ { 0 } ^ { 2 } = \frac { u _ { i s c o } - 6 } { u _ { i s c o } ^ { 3 } ( 4 u _ { i s c o } - 1 5 ) },$$

where u isco represents the ISCO radius for a given (scaled) Hubble parameter h 0 . We now have the improved lower bound u &gt; 6. A straightforward calculation shows that h 0 and u isco increase together. If there exists stable circular orbits for a given value of h 0 , then there exists an ISCO for this value of h 0 . This yields the maximum value of h 2 0 for which a stable circular orbit exists: this corresponds to the (global) maximum ISCO value, which is the maximum of u on u -( v ), i.e. u = 15 2. / Hence

$$h _ { 0 } | _ { \max } = h _ { 0 } | _ { u _ { i s c o } = \frac { 1 5 } { 2 } } = \frac { 2 } { 7 5 \sqrt { 3 } }. \quad \quad ( B 7 )$$

- [1] S. Chandrasekhar. The mathematical theory of black holes . Oxford University Press, 1998.
- [2] Z. Stuchl´ ık and S. Hled´ ık. Equatorial photon motion in the Kerr-Newman spacetimes with a non-zero cosmological constant. Classical and Quantum Gravity , 17:4541, 2000.
- [3] P. Amsterdamski, T. Bulik, D. Gondek-Rosinska, and W. Kluzniak. Marginally stable orbits around Maclaurin spheroids and low-mass quark stars. Astronomy and Astrophysics , 381:L21-L24, 2002.
- [4] M.A. Abramowicz and C.P. Fragile. Foundations of Black Hole Accretion Disk Theory. Living Reviews in Relativity , 16(1), 2013.
- [5] V. Faraoni and A. Jacques. Cosmological expansion and local physics. Physical Review D , 76:063510, 2007.
- [6] K. Bolejko, M.-N. C´ el´rier, e and A. Krasi´ski. n Inhomogeneous cosmological models: exact solutions and their applications. Classical and Quantum Gravity , 28:164002, 2011.
- [7] G.C. McVittie. The mass-particle in an expanding universe. Monthly Notices of the Royal Astronomical Society , 93:325, 1933.
- [8] B.C. Nolan. A point mass in an isotropic universe: Existence, uniqueness, and basic properties. Physical Review D , 58:064006, 1998.
- [9] R. Nandra, A.N. Lasenby, and M.P. Hobson. The effect of a massive object on an expanding universe. Monthly Notices of the Royal Astronomical Society , 422:2931, 2012.
- [10] We note that both [8] and [9] identify solutions of the Einstein equations that do for the k = -1 isotropic spacetimes what the McVittie spacetime does for the k = 0 class. In comoving coordinates, these solutions are given in terms of elliptic integrals, and so are not readily amenable to the analytic studies that have been done in
- the k = 0 case. As far as the author is aware, the possibility that these solutions have an elementary form in area radial coordinates has not been investigated: this may make the solutions more tractable.
- [11] B.C. Nolan. A point mass in an isotropic universe: II. Global properties. Classical and Quantum Gravity , 16:1227, 1999.
- [12] R.A. Sussman. On spherically symmetric shear-free perfect fluid configurations (neutral and charged). III. Global view. Journal of Mathematical Physics , 29:1177, 1988.
- [13] N. Kaloper, M. Kleban, and D. Martin. McVittie's legacy: black holes in an expanding universe. Physical Review D , 81:104044, 2010.
- [14] K. Lake and M. Abdelqader. More on McVittie's legacy: A Schwarzschild-de Sitter black and white hole embedded in an asymptotically Λ CDM cosmology. Physical Review D , 84:044045, 2011.
- [15] A.M. da Silva, M. Fontanini, and D.C. Guariento. How the expansion of the universe determines the causal structure of McVittie spacetimes. Physical Review D , 87:064030, 2013.
- [16] L. Perko. Differential equations and dynamical systems . Springer-Verlag, New York, 1991.
- [17] J. Carr. Applications of centre manifold theory . SpringerVerlag New York, 1981.
- [18] D.K. Arrowsmith and C.M. Place. An introduction to dynamical systems . Cambridge University Press, 1990.
- [19] U. Kirchgraber and K.J. Palmer. Geometry in the neighborhood of invariant manifolds of maps and flows and linearization . Longman Scientific &amp; Technical Harlow, 1990.
- [20] S.W. Hawking and G.F.R. Ellis. The large scale structure of space-time . Cambridge University Press, 1973.

-∞

FIG. 7. Conformal diagrams for the McVittie spacetime (bottom) that admits a circular photon orbit at r = a &gt; 3 m , and the corresponding isotropic background (top). In both cases, t increases in the vertical direction and J ± corresponds to r = + ∞ , t = ±∞ . The spacetime is expanding ( H &lt; 0) for t &lt; 0 and collapsing ( H &gt; 0) for t &gt; 0. In the McVittie case, the past singular boundaries (serrated) are at r = 2 m,t &lt; 0 and the future singular boundaries are at r = 2 m,t &gt; 0. At A , ( t, r ) = (+ ∞ , r -and at B , ( t, r ) = ( -∞ , r -). We note the existence of (a) ORNGs that originate at r = 2 m and terminate at r = 2 m ; (b) a unique ORNG ¯ γ O that originates at r = 2 m and extends to r = r + with t → + ∞ and (c) ORNGs that originate at r = 2 m and extend to r = r + with t → + ∞ . Similar results hold for IRNGs: ¯ γ I is the unique IRNG that originates at r = r -. We speculate that the spacetime extends to a Schwarzschild-de Sitter spacetime across r = r + , t = ±∞ , but we have not addressed this in the paper.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0044v1-figure-abdf1bf3-35c6-49e2-8369-695c7e9ea0ea.jpg)
Brien C. Nolan
Particle and photon orbits in McVittie spacetimes
McVittie spacetimes represent an embedding of the Schwarzschild field in
isotropic cosmological backgrounds. Depending on the scale factor of the
background, the resulting spacetime may contain black and white hole horizons,
as well as other interesting boundary features. In order to further clarify the
nature of these spacetimes, we address this question: do there exist bound
particle and photon orbits in McVittie spacetimes? Considering first circular
photon orbits, we obtain an explicit characterization of all McVittie
spacetimes for which such orbits exist: there is a 2-parameter class of such
spacetimes, and so the existence of a circular photon orbit is a highly
specialised feature of a McVittie spacetime. However, we prove that in two
large classes of McVittie spacetimes, there are bound particle and photon
orbits: future-complete non-radial timelike and null geodesics along which the
areal radius $r$ has a finite upper bound. These geodesics are asymptotic at
large times to circular orbits of a corresponding Schwarzschild or
Schwarzschild-de Sitter spacetime. The existence of these geodesics lays the
foundations for and shows the theoretical possibility of the formation of
accretion disks in McVittie spacetimes. We also summarize and extend some
previous results on the global structure of McVittie spacetimes. The results on
bound orbits are established using centre manifold and other techniques from
the theory of dynamical systems.