
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.

A method is presented for Bayesian model selection without explicitly
computing evidences, by using a combined likelihood and introducing an integer
model selection parameter $n$ so that Bayes factors, or more generally
posterior odds ratios, may be read off directly from the posterior of $n$. If
the total number of models under consideration is specified a priori, the full
joint parameter space $(\theta, n)$ of the models is of fixed dimensionality
and can be explored using standard Markov chain Monte Carlo (MCMC) or nested
sampling methods, without the need for reversible jump MCMC techniques. The
posterior on $n$ is then obtained by straightforward marginalisation. We
demonstrate the efficacy of our approach by application to several toy models.
We then apply it to constraining the dark energy equationofstate using a
freeform reconstruction technique. We show that $\Lambda$CDM is significantly
favoured over all extensions, including the simple $w(z){=}{\rm constant}$
model.

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.

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/.

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.

We present 16 GHz (1.9 cm) deep radio continuum observations made with the
Arcminute Microkelvin Imager (AMI) of a sample of lowmass young stars driving
jets. We combine these new data with archival information from an extensive
literature search to examine spectral energy distributions (SEDs) for each
source and calculate both the radio and submm spectral indices in two
different scenarios: (1) fixing the dust temperature (Td) according to
evolutionary class; (2) allowing Td to vary. We use the results of this
analysis to place constraints on the physical mechanisms responsible for the
radio emission. From AMI data alone, as well as from model fitting to the full
SED in both scenarios, we find that 80 per cent of the objects in this sample
have spectral indices consistent with freefree emission. We find an average
spectral index in both Td scenarios consistent with freefree emission. We
examine correlations of the radio luminosity with bolometric luminosity,
envelope mass, and outflow force and find that these data are consistent with
the strong correlation with envelope mass seen in lower luminosity samples. We
examine the errors associated with determining the radio luminosity and find
that the dominant source of error is the uncertainty on the opacity index,
beta. We examine the SEDs for variability in these young objects, and find
evidence for possible radio flare events in the histories of L1551 IRS 5 and
Serpens SMM 1.

Most SunyaevZel'dovich (SZ) and Xray analyses of galaxy clusters try to
constrain the cluster total mass and/or gas mass using parameterised models and
assumptions of spherical symmetry and hydrostatic equilibrium. By numerically
exploring the probability distributions of the cluster parameters given the
simulated interferometric SZ data in the context of Bayesian methods, and
assuming a betamodel for the electron number density we investigate the
capability of this model and analysis to return the simulated cluster input
quantities via three rameterisations. In parameterisation I we assume that the
T is an input parameter. We find that parameterisation I can hardly constrain
the cluster parameters. We then investigate parameterisations II and III in
which fg(r200) replaces temperature as a main variable. In parameterisation II
we relate M_T(r200) and T assuming hydrostatic equilibrium. We find that
parameterisation II can constrain the cluster physical parameters but the
temperature estimate is biased low. In parameterisation III, the virial theorem
replaces the hydrostatic equilibrium assumption. We find that parameterisation
III results in unbiased estimates of the cluster properties. We generate a
second simulated cluster using a generalised NFW (GNFW) pressure profile and
analyse it with an entropy based model to take into account the temperature
gradient in our analysis and improve the cluster gas density distribution. This
model also constrains the cluster physical parameters and the results show a
radial decline in the gas temperature as expected. The mean cluster total mass
estimates are also within 1 sigma from the simulated cluster true values.
However, we find that for at least interferometric SZ analysis in practice at
the present time, there is no differences in the AMI visibilities between the
two models. This may of course change as the instruments improve.

We present deep 1.8 cm (16 GHz) radio continuum imaging of seven young
stellar objects in the Taurus molecular cloud. These objects have previously
been extensively studied in the submm to NIR range and their SEDs modelled to
provide reliable physical and geometrical parametres.We use this new data to
constrain the properties of the longwavelength tail of the greybody spectrum,
which is expected to be dominated by emission from large dust grains in the
protostellar disk. We find spectra consistent with the opacity indices expected
for such a population, with an average opacity index of beta = 0.26+/0.22
indicating grain growth within the disks. We use spectra fitted jointly to
radio and submm data to separate the contributions from thermal dust and radio
emission at 1.8 cm and derive disk masses directly from the cmwave dust
contribution. We find that disk masses derived from these flux densities under
assumptions consistent with the literature are systematically higher than those
calculated from submm data, and meet the criteria for giant planet formation
in a number of cases.

The Planck Early Release Compact Source Catalogue includes 26 sources with no
obvious matches in other radio catalogues (of primarily extragalactic sources).
Here we present observations made with the Arcminute Microkelvin Imager Small
Array (AMI SA) at 15.75 GHz of the eight of the unmatched sources at
declination > +10 degrees. Of the eight, four are detected and are associated
with known objects. The other four are not detected with the AMI SA, and are
thought to be spurious.

We present deep radio continuum observations of the cores identified as
deeply embedded young stellar objects in the Serpens molecular cloud by the
Spitzer c2d programme at a wavelength of 1.8cm with the Arcminute Microkelvin
Imager Large Array (AMILA). These observations have a resolution of ~30arcsec
and an average sensitivity of 19microJy/beam. The targets are predominantly
Class I sources, and we find the detection rate for Class I objects in this
sample to be low (18%) compared to that of Class 0 objects (67%), consistent
with previous works. For detected objects we examine correlations of radio
luminosity with bolometric luminosity and envelope mass and find that these
data support correlations found by previous samples, but do not show any
indiction of the evolutionary divide hinted at by similar data from the Perseus
molecular cloud when comparing radio luminosity with envelope mass. We conclude
that envelope mass provides a better indicator for radio luminosity than
bolometric luminosity, based on the distribution of deviations from the two
correlations. Combining these new data with archival 3.6cm flux densities we
also examine the spectral indices of these objects and find an average spectral
index of 0.53+/1.14, consistent with the canonical value for a partially
optically thick spherical or collimated stellar wind. However, we caution that
possible interepoch variability limits the usefulness of this value, and such
variability is supported by our identification of a possible flare in the radio
history of Serpens SMM 1.

We present deep radio continuum observations of the cores identified as
deeply embedded young stellar objects in the Perseus molecular cloud by the
Spitzer c2d programme at a wavelength of 1.8 cm with the Arcminute Microkelvin
Imager Large Array (AMILA). We detect 72% of Class 0 objects from this sample
and 31% of Class I objects. No starless cores are detected. We use the flux
densities measured from these data to improve constraints on the correlations
between radio luminosity and bolometric luminosity, infrared luminosity and,
where measured, outflow force. We discuss the differing behaviour of these
objects as a function of protostellar class and investigate the differences in
radio emission as a function of core mass. Two of four possible very low
luminosity objects (VeLLOs) are detected at 1.8 cm.

We present 13.918.2 GHz observations of the SunyaevZel'dovich (SZ) effect
towards Abell 2146 using the Arcminute Microkelvin Imager (AMI). The cluster is
detected with a peak SNR ratio of 13 sigma in the radio source subtracted map.
Comparison of the SZ and Xray images suggests that they both have extended
regions which lie approximately perpendicular to one another, with their
emission peaks significantly displaced. These features indicate
nonuniformities in the distributions of the gas temperature and pressure,
indicative of a cluster merger. We use a Bayesian cluster analysis to explore
the highdimensional parameter space of the clusterplussources model to
obtain cluster parameter estimates in the presence of radio point sources,
receiver noise and primordial CMB anisotropy; the probability of SZ + CMB
primordial structure + radio sources + receiver noise to CMB + radio sources +
receiver noise is 3 x 10^{6}:1. We compare the results from three different
cluster models. Our preferred model exploits the observation that the gas
fractions do not appear to vary greatly between clusters. Given the relative
masses of the two merging systems in Abell 2146, the mean gas temperature can
be deduced from the virial theorem (assuming all of the kinetic energy is in
the form of internal gas energy) without being affected significantly by the
merger event, provided the primary cluster was virialized before the merger. In
this model we fit a simple spherical isothermal betamodel, despite the
inadequacy of this model for a merging system like Abell 2146, and assume the
cluster follows the masstemperature relation of a virialized, singular,
isothermal sphere. We note that this model avoids inferring largescale cluster
parameters internal to r_200 under the widely used assumption of hydrostatic
equilibrium. We find that at r_200 M_T= 4.1 \pm 0.5 x 10^{14} h^{1}M_sun and
T=4.5 \pm 0.5 keV.

We present observations from the Small Array of the Arcminute Microkelvin
Imager (AMI) of eight high Xray luminosity galaxy cluster systems selected
from the Local Cluster Substructure Survey (LoCuSS) sample.We detect the
SunyaevZel'dovich (SZ) effect in seven of these clusters. With the assumptions
that galaxy clusters are isothermal, have a density profile described by a
spherical b model and obey the theoretical MT relation, we are able to derive
cluster parameters at r200 from our SZ data. With the additional assumption of
hydrostatic equilibrium we are able to derive parameters at r500. We present
posterior probability distributions for cluster parameters such as mass, radius
and temperature (TSZ, MT). Combining our sample with that of AMI Consortium:
Rodr'\iguezGonz\'alvez et al. (2011) and using largeradius Xray temperature
estimates (TX) from Chandra and Suzaku observations, we find that there is
reasonable correspondence between TX and TSZ,MT values at low TX, but that for
clusters with TX above around 6keV the correspondence breaks down with TX
exceeding TSZ, MT; we stress that this finding is based on just ten clusters.

We perform deep 1.8 cm radio continuum imaging towards thirteen protostellar
regions selected from the Spitzer c2d small clouds and cores programme at high
resolution (25") in order to detect and quantify the cmwave emission from
deeply embedded young protostars. Within these regions we detect fifteen
compact radio sources which we identify as radio protostars including two
probable new detections. The sample is in general of low bolometric luminosity
and contains several of the newly detected VeLLO sources. We determine the 1.8
cm radio luminosity to bolometric luminosity correlation, L_rad L_bol, for the
sample and discuss the nature of the radio emission in terms of the available
sources of ionized gas. We also investigate the L_radL_IR correlation and
suggest that radio flux density may be used as a proxy for the internal
luminosity of low luminosity protostars.

We present observations using the Small Array of the Arcminute Microkelvin
Imager (AMI; 1418 GHz) of four Abell and three MACS clusters spanning
0.1710.686 in redshift. We detect SunyaevZel'dovich (SZ) signals in five of
these without any attempt at source subtraction, although strong source
contamination is present. With radiosource measurements from highresolution
observations, and under the assumptions of spherical $\beta$model,
isothermality and hydrostatic equilibrium, a Bayesian analysis of the data in
the visibility plane detects extended SZ decrements in all seven clusters over
and above receiver noise, radio sources and primary CMB imprints. Bayesian
evidence ratios range from 10^{11}:1 to 10^{43}:1 for six of the clusters and
3000:1 for one with substantially less data than the others. We present
posterior probability distributions for, e.g., total mass and gas fraction
averaged over radii internal to which the mean overdensity is 1000, 500 and
200, r_200 being the virial radius. Reaching r_200 involves some extrapolation
for the nearer clusters but not for the moredistant ones. We find that our
estimates of gas fraction are low (compared with most in the literature) and
decrease with increasing radius. These results appear to be consistent with the
notion that gas temperature in fact falls with distance (away from near the
cluster centre) out to the virial radius.

A large submm survey with Herschel will enable many exciting science
opportunities, especially in an era of widefield optical and radio surveys and
high resolution cosmic microwave background experiments. The HerschelSPIRE
Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250,
350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3
million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion
noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere
(3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which
have extensive multiwavelength coverage and are easily accessible from ALMA.
Two thirds of the of the sources are expected to be at z > 1, one third at z >
2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary
anisotropy studies with Planck and groundbased CMB data. (c) Find at least
1200 strongly lensed bright submm sources leading to a 2% test of general
relativity. (d) Identify 200 protocluster regions at z of 2 and perform an
unbiased study of the environmental dependence of star formation. (e) Perform
an unbiased survey for star formation and dust at high Galactic latitude and
make a census of debris disks and dust around AGB stars and white dwarfs.