
We introduce a new method for performing robust Bayesian estimation of the
threedimensional spatial power spectrum at the Epoch of Reionization (EoR),
from interferometric observations. The versatility of this technique allows us
to present two approaches. First, when the observations span only a small
number of independent spatial frequencies ($k$modes) we sample directly from
the spherical power spectrum coefficients that describe the EoR signal
realisation. Second, when the number of $k$modes to be included in the model
becomes large, we sample from the joint probability density of the spherical
power spectrum and the signal coefficients, using Hamiltonian Monte Carlo
methods to explore this high dimensional ($\sim$ 20000) space efficiently. This
approach has been successfully applied to simulated observations that include
astrophysically realistic foregrounds in a companion publication (Sims et al.
2016). Here we focus on explaining the methodology in detail, and use simple
foreground models to both demonstrate its efficacy, and highlight salient
features. In particular, we show that including an arbitrary flat spectrum
continuum foreground that is $10^8$ times greater in power than the EoR signal
has no detectable impact on our parameter estimates of the EoR power spectrum
recovered from the data.

We present a comparison between three cluster models applied to data obtained
by the Arcminute Microkelvin Imager radio interferometer system. The physical
model (PM) parameterises a cluster in terms of its physical quantities to model
the dark matter and baryonic components of the cluster using NFW and GNFW
profiles respectively. The observational models (OM I and OM II) model only the
gas content of the cluster. The two OMs vary only in the priors they use in
Bayesian inference: OM I has a joint prior on angular radius $\theta$ and
integrated Comptonisation $Y$, derived from simulations, while OM II uses
separable priors on $\theta$ and $Y$ which are based on calculations of the
physical model. For the comparison we consider a sample of $54$ clusters which
are a subsample of the second Planck catalogue of SunyaevZel'dovich sources.
We first compare the $Y$ estimates of the three models, and find that the PM
generally yields lower estimates relative to the OMs. We then compute the Earth
Mover's Distance between the $\theta$  $Y$ posterior distributions obtained
from each model for each cluster, and find that the two models which are most
discrepant are PM and OM I. Finally, we compare the Bayesian evidence values
obtained from each model for each cluster. OM I generally provides the best fit
to the data but not at a statistically significant level, according to the
Jeffreys scale. The highest evidence ratio obtained is actually in favour of
the PM over OM I.

We present a comparison of mass estimates for $54$ galaxy cluster candidates
from the second Planck catalogue (PSZ2) of SunyaevZel'dovich sources. We
compare the mass values obtained with data taken from the Arcminute Microkelvin
Imager (AMI) radio interferometer system and from the Planck satellite. The
former of these uses a Bayesian analysis pipeline that parameterises a cluster
in terms of its physical quantities, and models the dark matter & baryonic
components of a cluster using NFW and GNFW profiles respectively. Our mass
estimates derived from Planck data are obtained from the results of the
Bayesian detection algorithm PowellSnakes (PwS), are based on the methodology
detailed in the PSZ2 paper, and produce two sets of mass estimates; one
estimate is calculated directly from the angular radius $\theta$  integrated
Comptonisation parameter $Y$ posterior distributions, and the other uses a
`slicing function' to provide information on $\theta$ based on Xray
measurements and previous Planck mission samples. We find that for $37$ of the
clusters, the AMI mass estimates are lower than both values obtained from
Planck data. However the AMI and slicing function estimates are within one
combined standard deviation of each other for $31$ clusters. We also generate
cluster simulations based on the slicingfunction mass estimates, and analyse
them in the same way as we did the real AMI data. We find that inclusion in the
simulations of radiosource confusion & CMB noise and measurable radiosources
causes AMI mass estimates to be systematically low.

Extrememassratioinspiral observations from future spacebased
gravitationalwave detectors such as LISA will enable strongfield tests of
general relativity with unprecedented precision, but at prohibitive
computational cost if existing statistical techniques are used. In one such
test that is currently employed for LIGO blackhole binary mergers, generic
deviations from relativity are represented by $N$ deformation parameters in a
generalised waveform model; the Bayesian evidence for each of its $2^N$
combinatorial submodels is then combined into a posterior odds ratio for
modified gravity over relativity in a nullhypothesis test. We adapt and apply
this test to a generalised model for extrememassratio inspirals constructed
on deformed blackhole spacetimes, and focus our investigation on how
computational efficiency can be increased through an evidencefree method of
model selection. This method is akin to the algorithm known as productspace
Markov chain Monte Carlo, but uses nested sampling and improved error estimates
from a rethreading technique. We perform benchmarking and robustness checks for
the method, and find orderofmagnitude computational gains over regular nested
sampling in the case of synthetic data generated from the null model.

We apply our tetradbased approach for constructing sphericallysymmetric
solutions in general relativity to modelling a void, and compare it with the
standard Lema\^itreTolmanBondi (LTB) formalism. In particular, we construct
models for the void observed in the direction of Draco in the WISE2MASS galaxy
survey, and a corresponding cosmic microwave background (CMB) temperature
decrement in the Planck data in the same direction. We find that the
presentday density and velocity profiles of the void are not well constrained
by the existing data, so that void models produced from the two approaches can
differ substantially while remaining broadly consistent with the observations.
We highlight the importance of considering the velocity as well as the density
profile in constraining voids.

The Arcminute Microkelvin Imager (AMI) carried out a blind survey for galaxy
clusters via their SunyaevZel'dovich effect decrements between 2008 and 2011.
The first detection, known as AMICL J0300+2613, has been reobserved with AMI
equipped with a new digital correlator with high dynamic range. The combination
of the new AMI data and more recent highresolution submm and infrared maps
now shows the feature in fact to be a ring of positive dustcorrelated Galactic
emission, which is likely to be anomalous microwave emission (AME). If so, this
is the first completely blind detection of AME at arcminute scales.

AMI observations towards CIZA J2242+5301, in comparison with observations of
weak gravitational lensing and Xray emission from the literature, are used to
investigate the behaviour of nonbaryonic dark matter (NBDM) and gas during the
merger. Analysis of the SunyaevZel'dovich (SZ) signal indicates the presence
of high pressure gas elongated perpendicularly to the Xray and weaklensing
morphologies which, given the mergeraxis constraints in the literature,
implies that high pressure gas is pushed out into a linear structure during
core passing. Simulations in the literature closely matching the inferred
merger scenario show the formation of gas density and temperature structures
perpendicular to the merger axis. These SZ observations are challenging for
modified gravity theories in which NBDM is not the dominant contributor to
galaxycluster gravity.

A new method is presented for modelling the physical properties of galaxy
clusters. Our technique moves away from the traditional approach of assuming
specific parameterised functional forms for the variation of physical
quantities within the cluster, and instead allows for a 'freeform'
reconstruction, but one for which the level of complexity is determined
automatically by the observational data and may depend on position within the
cluster. This is achieved by representing each independent cluster property as
some interpolating or approximating function that is specified by a set of
control points, or 'nodes', for which the number of nodes, together with their
positions and amplitudes, are allowed to vary and are inferred in a Bayesian
manner from the data. We illustrate our nodal approach in the case of a
spherical cluster by modelling the electron pressure profile Pe(r) in analyses
both of simulated SunyaevZel'dovich (SZ) data from the Arcminute MicroKelvin
Imager (AMI) and of real AMI observations of the cluster MACS J0744+3927 in the
CLASH sample. We demonstrate that one may indeed determine the complexity
supported by the data in the reconstructed Pe(r), and that one may constrain
two very important quantities in such an analysis: the cluster total volume
integrated Comptonisation parameter (Ytot) and the extent of the gas
distribution in the cluster (rmax). The approach is also wellsuited to
detecting clusters in blind SZ surveys.

A novel approach is presented for fast generation of synthetic seismograms
due to microseismic events, using heterogeneous marine velocity models. The
partial differential equations (PDEs) for the 3D elastic wave equation have
been numerically solved using the Fourier domain pseudospectral method which
is parallelizable on the graphics processing unit (GPU) cards, thus making it
faster compared to traditional CPU based computing platforms. Due to
computationally expensive forward simulation of large geological models,
several combinations of individual synthetic seismic traces are used for
specified microseismic event locations, in order to simulate the effect of
realistic microseismic activity patterns in the subsurface. We here explore the
patterns generated by few hundreds of microseismic events with different source
mechanisms using various combinations, both in event amplitudes and origin
times, using the simulated pressure and three component particle velocity
fields via 1D, 2D and 3D seismic visualizations.

The $R_h = ct$ cosmological model has received considerable attention in
recent years owing to claims that it is favoured over the standard $\Lambda$
cold dark matter ($\Lambda$CDM) model by most observational data. A key feature
of the $R_h = ct$ model is that the zero active mass condition $\rho + 3p = 0$
holds at all epochs. Most recently, Melia has claimed that this condition is a
requirement of the symmetries of the FriedmannRobertsonWalker (FRW)
spacetime. We demonstrate that this claim is false and results from a flaw in
the logic of Melia's argument.

We present a tetradbased method for solving the Einstein field equations for
sphericallysymmetric systems and compare it with the widelyused
Lema\^itreTolmanBondi (LTB) model. In particular, we focus on the issues of
gauge ambiguity and the use of comoving versus 'physical' coordinate systems.
We also clarify the correspondences between the two approaches, and illustrate
their differences by applying them to the classic examples of the Schwarzschild
and FriedmannRobertsonWalker spacetimes. We demonstrate that the tetradbased
method does not suffer from the gauge freedoms inherent to the LTB model,
naturally accommodates nonzero pressure and has a more transparent physical
interpretation. We further apply our tetradbased method to a generalised form
of 'Swiss cheese' model, which consists of an interior spherical region
surrounded by a spherical shell of vacuum that is embedded in an exterior
background universe. In general, we allow the fluid in the interior and
exterior regions to support pressure, and do not demand that the interior
region be compensated. We pay particular attention to the form of the solution
in the intervening vacuum region and verify the validity of Birkhoff's theorem
at both the metric and tetrad level. We then reconsider critically the original
theoretical arguments underlying the socalled $R_h = ct$ cosmological model,
which has recently received considerable attention. These considerations in
turn illustrate the interesting behaviour of a number of 'horizons' in general
cosmological models.

Using Arcminute Microkelvin Imager (AMI) SZ observations towards ten CLASH
clusters we investigate the influence of cluster mergers on observational
galaxy cluster studies. Although selected to be largely relaxed, there is
disagreement in the literature on the dynamical states of CLASH sample members.
We analyse our AMI data in a fully Bayesian way to produce estimated cluster
parameters and consider the intrinsic correlations in our NFW/GNFWbased model.
Varying pressure profile shape parameters, illustrating an influence of mergers
on scaling relations, induces small deviations from the canonical selfsimilar
predictions  in agreement with simulations of Poole et al. 2007 who found
that merger activity causes only small scatter perpendicular to the relations.
We demonstrate this effect observationally using the different dependencies of
SZ and Xray signals to $n_{\rm e}$ that cause different sensitivities to the
shocking and/or fractionation produced by mergers. Plotting $Y_{\rm
X}$$M_{\rm gas}$ relations (where $Y_{\rm X}=M_{\rm gas}T$) derived from AMI
SZ and from $Chandra$ Xray gives ratios of AMI and $Chandra$ $Y_{\rm X}$ and
$M_{\rm gas}$ estimates that indicate movement of clusters \textit{along} the
scaling relation, as predicted by Poole et al. 2007. Clusters that have moved
most along the relation have the most discrepant $T_{\rm SZ}$ and $T_{\rm X}$
estimates: all the other clusters (apart from one) have SZ and Xray estimates
of $M_{\rm gas}$, $T$ and $Y_{\rm X}$ that agree within $r_{500}$. We use SZ vs
Xray discrepancies in conjunction with $Chandra$ maps and $T_{\rm X}$
profiles, making comparisons with simulated cluster merger maps in Poole et al.
2006, to identify disturbed members of our sample and estimate merger stages.

We introduce a method for performing a robust Bayesian analysis of
nonGaussianity present in pulsar timing data, simultaneously with the pulsar
timing model, and additional stochastic parameters such as those describing red
spin noise and dispersion measure variations. The parameters used to define the
presence of nonGaussianity are zero for Gaussian processes, giving a simple
method of defining the strength of nonGaussian behaviour. We use simulations
to show that assuming Gaussian statistics when the noise in the data is drawn
from a nonGaussian distribution can significantly increase the uncertainties
associated with the pulsar timing model parameters. We then apply the method to
the publicly available 15 year Parkes Pulsar Timing Array data release 1
dataset for the binary pulsar J0437$$4715. In this analysis we present a
significant detection of nonGaussianity in the uncorrelated nonthermal noise,
but we find that it does not yet impact the timing model or stochastic
parameter estimates significantly compared to analysis performed assuming
Gaussian statistics. The methods presented are, however, shown to be of
immediate practical use for current European Pulsar Timing Array (EPTA) and
International Pulsar Timing Array (IPTA) datasets.

We present the first public release of our generic neural network training
algorithm, called SkyNet. This efficient and robust machine learning tool is
able to train large and deep feedforward neural networks, including
autoencoders, for use in a wide range of supervised and unsupervised learning
applications, such as regression, classification, density estimation,
clustering and dimensionality reduction. SkyNet uses a `pretraining' method to
obtain a set of network parameters that has empirically been shown to be close
to a good solution, followed by further optimisation using a regularised
variant of Newton's method, where the level of regularisation is determined and
adjusted automatically; the latter uses secondorder derivative information to
improve convergence, but without the need to evaluate or store the full Hessian
matrix, by using a fast approximate method to calculate Hessianvector
products. This combination of methods allows for the training of complicated
networks that are difficult to optimise using standard backpropagation
techniques. SkyNet employs convergence criteria that naturally prevent
overfitting, and also includes a fast algorithm for estimating the accuracy of
network outputs. The utility and flexibility of SkyNet are demonstrated by
application to a number of toy problems, and to astronomical problems focusing
on the recovery of structure from blurred and noisy images, the identification
of gammaray bursters, and the compression and denoising of galaxy images. The
SkyNet software, which is implemented in standard ANSI C and fully parallelised
using MPI, is available at http://www.mrao.cam.ac.uk/software/skynet/.

Here we present a Bayesian method of including discrete measurements of
dispersion measure due to the interstellar medium in the direction of a pulsar
as prior information in the analysis of that pulsar. We use a simple simulation
to show the efficacy of this method, where the inclusion of the additional
measurements results in both a significant increase in the precision with which
the timing model parameters can be obtained, and an improved upper limit on the
amplitude of any red noise in the dataset. We show that this method can be
applied where no multifrequency data exists across much of the dataset, and
where there is no simultaneous multifrequency data for any given observing
epoch. Including such information in the analysis of upcoming International
Pulsar Timing Array (IPTA) and European Pulsar Timing Array (EPTA) data
releases could therefore prove invaluable in obtaining the most constraining
limits on gravitational wave signals within those datasets.

A new Bayesian software package for the analysis of pulsar timing data is
presented in the form of TempoNest which allows for the robust determination of
the nonlinear pulsar timing solution simultaneously with a range of additional
stochastic parameters. This includes both red spin noise and dispersion measure
variations using either power law descriptions of the noise, or through a
modelindependent method that parameterises the power at individual frequencies
in the signal. We use TempoNest to show that at noise levels representative of
current datasets in the European Pulsar Timing Array (EPTA) and International
Pulsar Timing Array (IPTA) the linear timing model can underestimate the
uncertainties of the timing solution by up to an order of magnitude. We also
show how to perform Bayesian model selection between different sets of timing
model and stochastic parameters, for example, by demonstrating that in the
pulsar B1937+21 both the dispersion measure variations and spin noise in the
data are optimally modelled by simple power laws. Finally we show that not
including the stochastic parameters simultaneously with the timing model can
lead to unpredictable variation in the estimated uncertainties, compromising
the robustness of the scientific results extracted from such analysis.

A new model independent method is presented for the analysis of pulsar timing
data and the estimation of the spectral properties of an isotropic
gravitational wave background (GWB). We show that by rephrasing the likelihood
we are able to eliminate the most costly aspects of computation normally
associated with this type of data analysis. When applied to the International
Pulsar Timing Array Mock Data Challenge data sets this results in speedups of
approximately 2 to 3 orders of magnitude compared to established methods. We
present three applications of the new likelihood. In the low signal to noise
regime we sample directly from the power spectrum coefficients of the GWB
signal realization. In the high signal to noise regime, where the data can
support a large number of coefficients, we sample from the joint probability
density of the power spectrum coefficients for the individual pulsars and the
GWB signal realization. Critically in both these cases we need make no
assumptions about the form of the power spectrum of the GWB, or the individual
pulsars. Finally we present a method for characterizing the spatial correlation
between pulsars on the sky, making no assumptions about the form of that
correlation, and therefore providing the only truly general Bayesian method of
confirming a GWB detection from pulsar timing data.

We have obtained deep SZ observations towards 15 of the apparently hottest
XMM Cluster Survey (XCS) clusters that can be observed with the Arcminute
Microkelvin Imager (AMI). We use a Bayesian analysis to quantify the
significance of our SZ detections. We detect the SZ effect at high significance
towards three of the clusters and at lower significance for a further two
clusters. Towards the remaining ten clusters, no clear SZ signal was measured.
We derive cluster parameters using the XCS mass estimates as a prior in our
Bayesian analysis. For all AMIdetected clusters, we calculate largescale mass
and temperature estimates while for all undetected clusters we determine upper
limits on these parameters. We find that the large scale mean temperatures
derived from our AMI SZ measurements (and the upper limits from null
detections) are substantially lower than the XCSbased coretemperature
estimates. For clusters detected in the SZ, the mean temperature is, on
average, a factor of 1.4 lower than temperatures from the XCS. For clusters
undetected in SZ, the average 68% upper limit on the mean temperature is a
factor of 1.9 below the XCS temperature.

We present CARMA observations of a massive galaxy cluster discovered in the
AMI blind SZ survey. Without knowledge of the cluster redshift a Bayesian
analysis of the AMI, CARMA and joint AMI & CARMA uvdata is used to quantify
the detection significance and parameterise both the physical and observational
properties of the cluster whilst accounting for the statistics of primary CMB
anisotropies, receiver noise and radio sources. The joint analysis of the AMI &
CARMA uvdata was performed with two parametric physical cluster models: the
{\beta}model; and the model described in Olamaie et al. 2012 with the pressure
profile fixed according to Arnaud et al. 2010. The cluster mass derived from
these different models is comparable but our Bayesian evidences indicate a
preference for the {\beta}profile which we, therefore, use throughout our
analysis. From the CARMA data alone we obtain a Bayesian probability of
detection ratio of 12.8:1 when assuming that a cluster exists within our search
area; alternatively assuming that Jenkins et al. 2001 accurately predicts the
number of clusters as a function of mass and redshift, the Bayesian probability
of detection is 0.29:1. From the analysis of the AMI or AMI & CARMA data the
probability of detection ratio exceeds 4.5x10^3:1. Performing a joint analysis
of the AMI & CARMA data with a physical cluster model we derive the total mass
internal to r200 as MT,200 = 4.1x10^14Msun. Using a phenomenological
{\beta}model to quantify the temperature decrement as a function of angular
distance we find a central SZ temperature decrement of 170{\mu}K in the AMI &
CARMA data. The SZ decrement in the CARMA data is weaker than expected and we
speculate that this is a consequence of the cluster morphology. In a
forthcoming study we will assess the impact of cluster morphology on the SZ
decrements that are observed with interferometers such as AMI and CARMA.

We present 16GHz SunyaevZel'dovich observations using the Arcminute
Microkelvin Imager (AMI) and subsequent Bayesian analysis of six galaxy
clusters at redshift ($z \approx 1$) chosen from an Xray and Infrared selected
sample from Culverhouse et al. (2010). In the subsequent analysis we use two
cluster models, an isothermal \betamodel and a Dark Matter GNFW (DMGNFW)
model in order to derive a formal detection probability and the cluster
parameters. We detect two clusters (CLJ1415+3612 & XMJ0830+5241) and measure
their total masses out to a radius of 200 $\times$ the critical density at the
respective cluster's redshift. For CLJ1415+3612 and XMJ0830+5241, we find
M_{\mathrm{T},200} for each model, which agree with each other for each
cluster. We also present maps before and after source subtraction of the entire
sample and provide 1D and 2D posterior marginalised probability distributions
for each fitted cluster profile parameter of the detected clusters. Using
simulations which take into account the measured source environment from the
AMI Large Array (LA), source confusion noise, CMB primordials, instrument
noise, we estimate from lowradius Xray data from Culverhouse et al. (2010),
the detectability of each cluster in the sample and compare it with the result
from the Small Array (SA) data. Furthermore, we discuss the validity of the
assumptions of isothermality and constant gas mass fraction. We comment on the
bias that these smallradius estimates introduce to largeradius SZ
predictions. In addition, we followup the two detections with deep,
singlepointed LA observations. We find a 3 sigma tentative decrement toward
CLJ1415+3612 at highresolution and a 5 sigma highresolution decrement towards
XMJ0830+5241.

The AMI Galactic Plane Survey (AMIGPS) is a large area survey of the outer
Galactic plane to provide arcminute resolution images at milliJansky
sensitivity in the centimetrewave band. Here we present the first data release
of the survey, consisting of 868 deg^2 of the Galactic plane, covering the area
76 deg \lessapprox l \lessapprox 170 deg between latitudes of b \lessapprox 5
deg, at a central frequency of 15.75 GHz (1.9 cm). We describe in detail the
drift scan observations which have been used to construct the maps, including
the techniques used for observing, mapping and source extraction, and summarise
the properties of the finalized datasets. These observations constitute the
most sensitive Galactic plane survey of large extent at centimetrewave
frequencies greater than 1.4 GHz.

We present 16GHz AMI SZ observations of 19 clusters with L_X >7x10^37 W
(h50=1) selected from the LoCuS survey (0.142<z<0.295) and of A1758b, in the
FoV of A1758a. We detect 17 clusters with 523sigma peak surface brightnesses.
Cluster parameters are obtained using a Bayesian cluster analysis. We fit
isothermal betamodels to our data and assume the clusters are virialized (with
all the kinetic energy in gas internal energy). Our gas temperature, T_AMI, is
derived from AMI SZ data, not from Xray spectroscopy. Cluster parameters
internal to r500 are derived assuming HSE. We find: (i) Different gNFW
parameterizations yield significantly different parameter degeneracies. (ii)
For h70 = 1, we find the virial radius r200 to be typically 1.6+/0.1 Mpc and
the total mass M_T(r200) typically to be 2.02.5xM_T(r500).(iii) Where we have
found M_T Xray (X) and weaklensing (WL) values in the literature, there is
good agreement between WL and AMI estimates (with M_{T,AMI}/M_{T,WL}
=1.2^{+0.2}_{0.3} and =1.0+/0.1 for r500 and r200, respectively). In
comparison, most Suzaku/Chandra estimates are higher than for AMI (with
M_{T,X}/M_{T,AMI}=1.7+/0.2 within r500), particularly for the stronger
mergers.(iv) Comparison of T_AMI to T_X sheds light on high Xray masses: even
at large r, T_X can substantially exceed T_AMI in mergers. The use of these
higher T_X values will give higher Xray masses. We stress that larger T_SZ
and T_X data are scarce and must be increased. (v) Despite the paucity of data,
there is an indication of a relation between merger activity and SZ
ellipticity. (vi) At small radius (but away from any cooling flow) the SZ
signal (and T_AMI) is less sensitive to ICM disturbance than the Xray signal
(and T_X) and, even at high r, mergers affect n^2weighted Xray data more than
nweighted SZ, implying significant shocking or clumping or both occur even in
the outer parts of mergers.

We present 16GHz observations using the Arcminute Microkelvin Imager (AMI)
of 11 clusters with 7 x 10^{37}W < L_X < 11 x 10^{37}W (h_{50}=1.0) selected
from the Local Cluster Substructure Survey (LoCuSS) and compare them to Xray
data. We use a fast, Bayesian cluster analysis to explore the highdimensional
parameter space of the clusterplussources model and obtain robust cluster
parameter estimates in the presence of radio point sources, receiver noise and
primordial CMB anisotropy. Our analysis fits a spherical, isothermal betamodel
to our data and assumes the cluster follows the theoretical masstemperature
relation. Largescale cluster parameters internal to r_{500} are derived under
the assumption of hydrostatic equilibrium. Posterior distributions for the
largescale parameters of 8 of our clusters are given; SZ effects towards Abell
1704 and Zw0857.9+2107 were not detected and our spherical betaprofile was
found to be an inadequate fit to the decrement on our map for Abell 2409.

A tetradbased procedure is presented for solving Einstein's field equations
for sphericallysymmetric systems; this approach was first discussed by Lasenby
et al. in the language of geometric algebra. The method is used to derive
metrics describing a point mass in a spatiallyflat, open and closed expanding
universe respectively. In the spatiallyflat case, a simple coordinate
transformation relates the metric to the corresponding one derived by McVittie.
Nonetheless, our use of noncomoving (`physical') coordinates greatly
facilitates physical interpretation. For the open and closed universes, our
metrics describe different spacetimes to the corresponding McVittie metrics and
we believe the latter to be incorrect. In the closed case, our metric possesses
an image mass at the antipodal point of the universe. We calculate the geodesic
equations for the spatiallyflat metric and interpret them. For radial motion
in the Newtonian limit, the force acting on a test particle consists of the
usual $1/r^2$ inwards component due to the central mass and a cosmological
component proportional to $r$ that is directed outwards (inwards) when the
expansion of the universe is accelerating (decelerating). For the standard
$\Lambda$CDM concordance cosmology, the cosmological force reverses direction
at about $z\approx 0.67$. We also derive an invariant fully
generalrelativistic expression, valid for arbitrary sphericallysymmetric
systems, for the force required to hold a test particle at rest relative to the
central point mass.

We present some astrophysical consequences of the metric for a point mass in
an expanding universe derived in Nandra, Lasenby & Hobson, and of the
associated invariant expression for the force required to keep a test particle
at rest relative to the central mass. We focus on the effect of an expanding
universe on massive objects on the scale of galaxies and clusters. Using
Newtonian and generalrelativistic approaches, we identify two important
timedependent physical radii for such objects when the cosmological expansion
is accelerating. The first radius, $r_F$, is that at which the total radial
force on a test particle is zero, which is also the radius of the largest
possible circular orbit about the central mass $m$ and where the gas pressure
and its gradient vanish. The second radius, $r_S$, which is \approx r_F/1.6$,
is that of the largest possible stable circular orbit, which we interpret as
the theoretical maximum size for an object of mass $m$. In contrast, for a
decelerating cosmological expansion, no such finite radii exist. Assuming a
cosmological expansion consistent with a $\Lambda$CDM concordance model, at the
present epoch we find that these radii put a sensible constraint on the typical
sizes of both galaxies and clusters at low redshift. For galaxies, we also find
that these radii agree closely with zeroes in the radial velocity field in the
neighbourhood of nearby galaxies, as inferred by Peirani & Pacheco from recent
observations of stellar velocities. We then consider the future effect on
massive objects of an accelerating cosmological expansion driven by phantom
energy, for which the universe is predicted to end in a `Big Rip' at a finite
time in the future at which the scale factor becomes singular. In particular,
we present a novel calculation of the time prior to the Big Rip that an object
of a given mass and size will become gravitationally unbound.