The excess of electron recoil events seen by the XENON1T experiment has been interpreted as a potential signal of axion-like particles (ALPs), either produced in the Sun, or constituting part of the dark matter halo of the Milky Way. It has also been explained as a consequence of trace amounts of tritium in the experiment. We consider the evidence for the solar and dark-matter ALP hypotheses from the combination of XENON1T data and multiple astrophysical probes, including horizontal branch stars, red giants, and white dwarfs. We briefly address the influence of ALP decays and supernova cooling. While the different datasets are in clear tension for the case of solar ALPs, all measurements can be simultaneously accommodated for the case of a sub-dominant fraction of dark-matter ALPs. Nevertheless, this solution requires the tuning of several a priori unknown parameters, such that for our choices of priors a Bayesian analysis shows no strong preference for the ALP interpretation of the XENON1T excess over the background hypothesis.

We present global fits of cosmologically stable axion-like particle and QCD axion models in the mass range 0.1 neV to 10 eV. We focus on the case where the Peccei-Quinn symmetry is broken before the end of inflation, such that the initial value of the axion field can be considered to be homogeneous throughout the visible Universe. We include detailed likelihood functions from light-shining-through-wall experiments, haloscopes, helioscopes, the axion relic density, horizontal branch stars, supernova 1987A, white dwarf cooling, and gamma-ray observations. We carry out both frequentist and Bayesian analyses, with and without the inclusion of white dwarf cooling. We explore the degree of fine-tuning present in different models and identify parameter regions where it is possible for QCD axion models to account for both the dark matter in the Universe and the cooling hints, comparing them to specific DFSZ- and KSVZ-type models. We find the most credible parameter regions, allowing us to set (prior-dependent) upper and lower bounds on the axion mass. Our analysis also suggests that QCD axions in this scenario most probably make up a non-negligible but sub-dominant component of the dark matter in the Universe.

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of earth-massed habitable worlds around the nearest starts, and also into distant Milky way objects up to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry ESA missions: NEAT for M3, micro-NEAT for S1 mission, and Theia for M4 and M5. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this white paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review quickly instrumentation and mission profiles.

We present a two-field model that realises inflation and the observed density of dark energy today, whilst solving the fine-tuning problems inherent in quintessence models. One field acts as the inflaton, generically driving the other to a saddle-point of the potential, from which it acts as a quintessence field following electroweak symmetry breaking. The model exhibits essentially no sensitivity to the initial value of the quintessence field, naturally suppresses its interactions with other fields, and automatically endows it with a small effective mass in the late Universe. The magnitude of dark energy today is fixed by the height of the saddle point in the potential, which is dictated entirely by the scale of electroweak symmetry breaking.

Searches for supersymmetric electroweakinos have entered a crucial phase, as the integrated luminosity of the Large Hadron Collider is now high enough to compensate for their weak production cross-sections. Working in a framework where the neutralinos and charginos are the only light sparticles in the Minimal Supersymmetric Standard Model, we use gambit to perform a detailed likelihood analysis of the electroweakino sector. We focus on the impacts of recent ATLAS and CMS searches with 36 fb$^{-1}$ of 13 TeV proton-proton collision data. We also include constraints from LEP and invisible decays of the $Z$ and Higgs bosons. Under the background-only hypothesis, we show that current LHC searches do not robustly exclude any range of neutralino or chargino masses. However, a pattern of excesses in several LHC analyses points towards a possible signal, with neutralino masses of $(m_{\tilde{\chi}_1^0}, m_{\tilde{\chi}_2^0}, m_{\tilde{\chi}_3^0}, m_{\tilde{\chi}_4^0})$ = (8-155, 103-260, 130-473, 219-502) GeV and chargino masses of $(m_{\tilde{\chi}_1^{\pm}}, m_{\tilde{\chi}_2^{\pm}})$ = (104-259, 224-507) GeV at the 95% confidence level. The lightest neutralino is mostly bino, with a possible modest Higgsino or wino component. We find that this excess has a combined local significance of $3.3\sigma$, subject to a number of cautions. If one includes LHC searches for charginos and neutralinos conducted with 8 TeV proton-proton collision data, the local significance is lowered to 2.9$\sigma$. We briefly consider the implications for dark matter, finding that the correct relic density can be obtained through the Higgs-funnel and $Z$-funnel mechanisms, even assuming that all other sparticles are decoupled. All samples, gambit input files and best-fit models from this study are available on Zenodo.

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.

We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.

We study the seven-dimensional Minimal Supersymmetric Standard Model (MSSM7) with the new GAMBIT software framework, with all parameters defined at the weak scale. Our analysis significantly extends previous weak-scale, phenomenological MSSM fits, by adding more and newer experimental analyses, improving the accuracy and detail of theoretical predictions, including dominant uncertainties from the Standard Model, the Galactic dark matter halo and the quark content of the nucleon, and employing novel and highly-efficient statistical sampling methods to scan the parameter space. We find regions of the MSSM7 that exhibit co-annihilation of neutralinos with charginos, stops and sbottoms, as well as models that undergo resonant annihilation via both light and heavy Higgs funnels. We find high-likelihood models with light charginos, stops and sbottoms that have the potential to be within the future reach of the LHC. Large parts of our preferred parameter regions will also be accessible to the next generation of direct and indirect dark matter searches, making prospects for discovery in the near future rather good.

One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.

Scalar singlet dark matter is one of the simplest and most predictive realisations of the WIMP (weakly-interacting massive particle) idea. Although the model is constrained from all directions by the latest experimental data, it still has viable regions of parameter space. Another compelling aspect of scalar singlets is their ability to stabilise the electroweak vacuum. Indeed, models of scalar dark matter are not low-energy effective theories, but can be valid all the way to the Planck scale. Using the GAMBIT framework, we present the first global fit to include both the low-energy experimental constraints and the theoretical constraints from UV physics, considering models with a scalar singlet charged under either a $\mathbb{Z}_2$ or a $\mathbb{Z}_3$ symmetry. We show that if the model is to satisfy all experimental constraints, completely stabilise the electroweak vacuum up to high scales, and also remain perturbative to those scales, one is driven to a relatively small region of parameter space. This region has a Higgs-portal coupling slightly less than 1, a dark matter mass of 1 to 2 TeV and a spin-independent nuclear scattering cross-section around 10$^{-45}$ cm$^2$.

An excess of gamma rays has been identified at the centre of the Milky Way, and annihilation of dark matter has been posited as a potential source. This hypothesis faces significant challenges: difficulty characterizing astrophysical backgrounds, the need for a non-trivial adiabatic contraction of the inner part of the Milky Way's dark matter halo, and recent observations of photon clustering, which suggest that the majority of the excess is due to unresolved point sources. Here we point out that the apparent point-like nature of the emission rules out the dark matter interpretation of the excess entirely. Attempting to model the emission with dark matter point sources either worsens the problem with the inner slope, requires an unrealistically large minihalo fraction toward the Galactic Centre, or overproduces the observed emission at higher latitudes.

The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Alberta.

Low-frequency radio continuum observations of edge-on galaxies are ideal to study cosmic-ray electrons (CREs) in halos via radio synchrotron emission and to measure magnetic field strengths. We obtained new observations of the edge-on spiral galaxy NGC 891 at 129-163 MHz with the LOw Frequency ARray (LOFAR) and at 13-18 GHz with the Arcminute Microkelvin Imager (AMI) and combine them with recent high-resolution Very Large Array (VLA) observations at 1-2 GHz, enabling us to study the radio continuum emission over two orders of magnitude in frequency. The spectrum of the integrated nonthermal flux density can be fitted by a power law with a spectral steepening towards higher frequencies or by a curved polynomial. Spectral flattening at low frequencies due to free-free absorption is detected in star-forming regions of the disk. The mean magnetic field strength in the halo is 7 +- 2 $\mu$G. The scale heights of the nonthermal halo emission at 146 MHz are larger than those at 1.5 GHz everywhere, with a mean ratio of 1.7 +- 0.3, indicating that spectral ageing of CREs is important and that diffusive propagation dominates. The halo scale heights at 146 MHz decrease with increasing magnetic field strengths which is a signature of dominating synchrotron losses of CREs. On the other hand, the spectral index between 146 MHz and 1.5 GHz linearly steepens from the disk to the halo, indicating that advection rather than diffusion is the dominating CRE transport process. This issue calls for refined modelling of CRE propagation.

The Arcminute Microkelvin Imager (AMI) telescopes located at the Mullard Radio Astronomy Observatory near Cambridge have been significantly enhanced by the implementation of a new digital correlator with 1.2 MHz spectral resolution. This system has replaced a 750-MHz resolution analogue lag-based correlator, and was designed to mitigate the effects of radio frequency interference, particularly from geostationary satellites that contaminate observations at low declinations. The upgraded instrument consists of 18 ROACH2 Field Programmable Gate Array platforms used to implement a pair of real-time FX correlators -- one for each of AMI's two arrays. The new system separates the down-converted RF baseband signal from each AMI receiver into two 2.3 GHz-wide sub-bands which are each digitized at 5-Gsps with 8 bits of precision. These digital data streams are filtered into 2048 frequency channels and cross-correlated using FPGA hardware, with a commercial 10 Gb Ethernet switch providing high-speed data interconnect. Images formed using data from the new digital correlator show over an order of magnitude improvement in dynamic range over the previous system. The ability to observe at low declinations has also been significantly improved.

QUIJOTE (Q-U-I JOint TEnerife) is an experiment designed to achieve CMB B-mode polarization detection and sensitive enough to detect a primordial gravitational-wave component if the B-mode amplitude is larger than r = 0.05. It consists in two telescopes and three instruments observing in the frequency range 10-42 GHz installed at the Teide Observatory in the Canary Islands, Spain. The observing strategy includes three raster scan deep integration fields for cosmology, a nominal wide survey covering the Northen Sky and specific raster scan deep integration observations in regions of specific interest. The main goals of the project are presented and the first scientific results obtained with the first instrument are reviewed.

We report on late time radio and X-ray observations of the tidal disruption event candidate ASASSN-14li, covering the first 1000 days of the decay phase. For the first $\sim200$ days the radio and X-ray emission fade in concert. This phase is better fit by an exponential decay at X-ray wavelengths, while the radio emission is well described by either an exponential or the canonical $t^{-5/3}$ decay assumed for tidal disruption events. The correlation between radio and X-ray emission during this period can be fit as $L_{R}\propto L_{X}^{1.9\pm0.2}$. After 400 days the radio emission at $15.5\,\textrm{GHz}$ has reached a plateau level of $244\pm8\,\mu\textrm{Jy}$ which it maintains for at least the next 600 days, while the X-ray emission continues to fade exponentially. This steady level of radio emission is likely due to relic radio lobes from the weak AGN-like activity implied by historical radio observations. We note that while most existing models are based upon the evolution of ejecta which are decoupled from the central black hole, the radio : X-ray correlation during the declining phase is also consistent with core jet emission coupled to a radiatively efficient accretion flow.

We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamlike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDcalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs, and the neutrino telescope likelihood package nuLike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for results obtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application of DarkBit's advanced direct and indirect detection routines to a simple effective dark matter model.

We present a deep survey of the SuperCLASS super-cluster - a region of sky known to contain five Abell clusters at redshift $z\sim0.2$ - performed using the Arcminute Microkelvin Imager (AMI) Large Array (LA) at 15.5$~$GHz. Our survey covers an area of approximately 0.9 square degrees. We achieve a nominal sensitivity of $32.0~\mu$Jy beam$^{-1}$ toward the field centre, finding 80 sources above a $5\sigma$ threshold. We derive the radio colour-colour distribution for sources common to three surveys that cover the field and identify three sources with strongly curved spectra - a high-frequency-peaked source and two GHz-peaked-spectrum sources. The differential source count (i) agrees well with previous deep radio source count, (ii) exhibits no evidence of an emerging population of star-forming galaxies, down to a limit of 0.24$~$mJy, and (iii) disagrees with some models of the 15$~$GHz source population. However, our source count is in agreement with recent work that provides an analytical correction to the source count from the SKADS Simulated Sky, supporting the suggestion that this discrepancy is caused by an abundance of flat-spectrum galaxy cores as-yet not included in source population models.

Ultracompact Minihalos (UCMHs) have been proposed as a type of dark matter sub-structure seeded by large-amplitude primordial perturbations and topological defects. UCMHs are expected to survive to the present era, allowing constraints to be placed on their cosmic abundance using observations within our own Galaxy. Constraints on their number density can be linked to conditions in the early universe that impact structure formation, such as increased primordial power on small scales, generic weak non-Gaussianity, and the presence of cosmic strings. We use new constraints on the abundance of UCMHs from pulsar timing to place generalised limits on the parameters of each of these cosmological scenarios. At some scales, the limits are the strongest to date, exceeding those from dark matter annihilation. Our new limits have the added advantage of being independent of the particle nature of dark matter, as they are based only on gravitational effects.

Small-scale dark matter structure within the Milky Way is expected to affect pulsar timing. The change in gravitational potential induced by a dark matter halo passing near the line of sight to a pulsar would produce a varying delay in the light travel time of photons from the pulsar. Individual transits produce an effect that would either be too rare or too weak to be detected in 30-year pulsar observations. However, a population of dark matter subhalos would be expected to produce a detectable effect on the measured properties of pulsars if the subhalos constitute a significant fraction of the total halo mass. The effect is to increase the dispersion of measured period derivatives across the pulsar population. By statistical analysis of the ATNF pulsar catalogue, we place an upper limit on this dispersion of $\log \sigma_{\dot{P}} \leq -17.05$. We use this to place strong upper limits on the number density of ultracompact minihalos within the Milky Way. These limits are completely independent of the particle nature of dark matter.

I give a brief review of some of the implications of Fermi data for theories of the identity of dark matter, and their combination with data from other complementary probes. I also preview some of the prospects for probing such models with future data.

We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to external codes, a suite of different statistical methods and parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easily-extendible than in the past. Here we give a detailed description of the framework, its design and motivation, and the current models and other specific components presently implemented in GAMBIT. Accompanying papers deal with individual modules and present first GAMBIT results. GAMBIT can be downloaded from gambit.hepforge.org.

Dark matter with momentum- or velocity-dependent interactions with nuclei has shown significant promise for explaining the so-called Solar Abundance Problem, a longstanding discrepancy between solar spectroscopy and helioseismology. The best-fit models are all rather light, typically with masses in the range of 3-5 GeV. This is exactly the mass range where dark matter evaporation from the Sun can be important, but to date no detailed calculation of the evaporation of such models has been performed. Here we carry out this calculation, for the first time including arbitrary velocity- and momentum-dependent interactions, thermal effects, and a completely general treatment valid from the optically thin limit all the way through to the optically thick regime. We find that depending on the dark matter mass, interaction strength and type, the mass below which evaporation is relevant can vary from 1 to 4 GeV. This has the effect of weakening some of the better-fitting solutions to the Solar Abundance Problem, but also improving a number of others. As a by-product, we also provide an improved derivation of the capture rate that takes into account thermal and optical depth effects, allowing the standard result to be smoothly matched to the well-known saturation limit.

We introduce ScannerBit, the statistics and sampling module of the public, open-source global fitting framework GAMBIT. ScannerBit provides a standardised interface to different sampling algorithms, enabling the use and comparison of multiple computational methods for inferring profile likelihoods, Bayesian posteriors, and other statistical quantities. The current version offers random, grid, raster, nested sampling, differential evolution, Markov Chain Monte Carlo (MCMC) and ensemble Monte Carlo samplers. We also announce the release of a new standalone differential evolution sampler, Diver, and describe its design, usage and interface to ScannerBit. We subject Diver and three other samplers (the nested sampler MultiNest, the MCMC GreAT, and the native ScannerBit implementation of the ensemble Monte Carlo algorithm T-Walk) to a battery of statistical tests. For this we use a realistic physical likelihood function, based on the scalar singlet model of dark matter. We examine the performance of each sampler as a function of its adjustable settings, and the dimensionality of the sampling problem. We evaluate performance on four metrics: optimality of the best fit found, completeness in exploring the best-fit region, number of likelihood evaluations, and total runtime. For Bayesian posterior estimation at high resolution, T-Walk provides the most accurate and timely mapping of the full parameter space. For profile likelihood analysis in less than about ten dimensions, we find that Diver and MultiNest score similarly in terms of best fit and speed, outperforming GreAT and T-Walk; in ten or more dimensions, Diver substantially outperforms the other three samplers on all metrics.

I give a brief review of a few recent developments and future directions in the search for dark matter using high-energy neutrinos from the Sun. This includes the ability to recast neutrino telescope limits on nuclear scattering of dark matter to arbitrary new theories, and new calculations of the solar atmospheric background relevant to such searches. I also touch on applications to global searches for new physics, and prospects for improving searches for asymmetric dark matter in the Sun.

The Arcminute Microkelvin Imager (AMI) carried out a blind survey for galaxy clusters via their Sunyaev-Zel'dovich effect decrements between 2008 and 2011. The first detection, known as AMI-CL 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 high-resolution sub-mm and infra-red maps now shows the feature in fact to be a ring of positive dust-correlated 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.

The candidate black hole X-ray binary Swift J1753.5-0127 faded to quiescence in 2016 November, after a prolonged outburst that was discovered in 2005. Nearly three months later the system displayed renewed activity that lasted through 2017 July. Here, we present radio and X-ray monitoring over ~3 months of the renewed activity to study the coupling between the jet and the inner regions of the disk/jet system. Our observations cover low X-ray luminosities that have not historically been well-sampled (Lx~2e33 - 1e36 erg/s; 1-10 keV), including time periods when the system was both brightening and fading. At these low luminosities Swift J1753.5-0127 occupies a parameter space in the radio/X-ray luminosity plane that is comparable to "canonical" systems (e.g., GX 339-4), regardless of whether the system was brightening or fading, even though during its >11-year outburst Swift J1753.5-0127 emitted less radio emission from its jet than expected. We discuss implications for the existence of a single radio/X-ray luminosity correlation for black hole X-ray binaries at the lowest luminosities (Lx < 1e35 erg/s), and we compare to supermassive black holes. Our campaign includes the lowest luminosity quasi-simultaneous radio/X-ray detection to date for a black hole X-ray binary during its rise out of quiescence, thanks to early notification from optical monitoring combined with fast responses from sensitive multiwavelength facilities.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

The existence of substructure in halos of annihilating dark matter would be expected to substantially boost the rate at which annihilation occurs. Ultracompact minihalos of dark matter (UCMHs) are one of the more extreme examples of this. The boosted annihilation can inject significant amounts of energy into the gas of a galaxy over its lifetime. Here we determine the impact of the boost factor from UCMH substructure on the heating of galactic gas in a Milky Way-type galaxy, by means of N-body simulation. If $1\%$ of the dark matter exists as UCMHs, the corresponding boost factor can be of order $10^5$. For reasonable values of the relevant parameters (annihilation cross section $3\times10^{-26} ~\textrm{cm}^3~ \textrm{s}^{-1}$, dark matter mass 100 GeV, 10% heating efficiency), we show that the presence of UCMHs at the 0.1% level would inject enough energy to eject significant amounts of gas from the halo, potentially preventing star formation within $\sim$1 kpc of the halo centre.

We present results of the coordinated observing campaign that made the first subarcsecond localization of a Fast Radio Burst, FRB 121102. During this campaign, we made the first simultaneous detection of an FRB burst by multiple telescopes: the VLA at 3 GHz and the Arecibo Observatory at 1.4 GHz. Of the nine bursts detected by the Very Large Array at 3 GHz, four had simultaneous observing coverage at other observatories. We use multi-observatory constraints and modeling of bursts seen only at 3 GHz to confirm earlier results showing that burst spectra are not well modeled by a power law. We find that burst spectra are characterized by a ~500 MHz envelope and apparent radio energy as high as $10^{40}$ erg. We measure significant changes in the apparent dispersion between bursts that can be attributed to frequency-dependent profiles or some other intrinsic burst structure that adds a systematic error to the estimate of DM by up to 1%. We use FRB 121102 as a prototype of the FRB class to estimate a volumetric birth rate of FRB sources $R_{FRB} \approx 5x10^{-5}/N_r$ Mpc$^{-3}$ yr$^{-1}$, where $N_r$ is the number of bursts per source over its lifetime. This rate is broadly consistent with models of FRBs from young pulsars or magnetars born in superluminous supernovae or long gamma-ray bursts, if the typical FRB repeats on the order of thousands of times during its lifetime.

In recent years, a revised set of solar abundances has led to a discrepancy in the sound-speed profile between helioseismology and theoretical solar models. Conventional solutions require additional mechanisms for energy transport within the Sun. Vincent et al. have recently suggested that dark matter with a momentum or velocity dependent cross section could provide a solution. In this work, we consider three models of dark matter with such cross sections and their effect on the stellar structure. In particular, the three models incorporate dark matter particles interacting through an electromagnetic dipole moment: an electric dipole, a magnetic dipole or an anapole. Each model is implemented in the \texttt{DarkStec} stellar evolution program, which incorporates the effects of dark matter capture and heat transport within the solar interior. We show that dark matter with an anapole moment of $\sim1\mathrm{GeV}^{-2}$ or magnetic dipole moment of $\sim10^{-3}\mu_p$ can improve the sound-speed profile, small frequency separations and convective zone radius with respect to the Standard Solar Model. However, the required dipole moments are strongly excluded by direct detection experiments.

The connection between accretion and jet production in accreting white dwarf binary systems, especially dwarf novae, is not well understood. Radio wavelengths provide key insights into the mechanisms responsible for accelerating electrons, including jets and outflows. Here we present densely-sampled radio coverage, obtained with the Arcminute MicroKelvin Imager Large Array, of the dwarf nova SS Cyg during its February 2016 anomalous outburst. The outburst displayed a slower rise (3 days per mag) in the optical than typical ones, and lasted for more than 3 weeks. Rapid radio flaring on timescales <1 hour was seen throughout the outburst. The most intriguing behavior in the radio was towards the end of the outburst where a fast, luminous ("giant"), flare peaking at ~20 mJy and lasting for 15 minutes was observed. This is the first time that such a flare has been observed in SS Cyg, and insufficient coverage could explain its non-detection in previous outbursts. These data, together with past radio observations, are consistent with synchrotron emission from plasma ejection events as being the origin of the radio flares. However, the production of the giant flare during the declining accretion rate phase remains unexplained within the standard accretion-jet framework and appears to be markedly different to similar patterns of behavior in X-ray binaries.

Interpreting samples from likelihood or posterior probability density functions is rarely as straightforward as it seems it should be. Producing publication-quality graphics of these distributions is often similarly painful. In this short note I describe pippi, a simple, publicly-available package for parsing and post-processing such samples, as well as generating high-quality PDF graphics of the results. Pippi is easily and extensively configurable and customisable, both in its options for parsing and post-processing samples, and in the visual aspects of the figures it produces. I illustrate some of these using an existing supersymmetric global fit, performed in the context of a gamma-ray search for dark matter. Pippi can be downloaded and followed at http://github.com/patscott/pippi .

After more than 26 years in quiescence, the black hole transient V404 Cyg went into a luminous outburst in June 2015, and additional activity was detected in late December of the same year. Here, we present an optical spectroscopic follow-up of the December mini-outburst, together with X-ray, optical and radio monitoring that spanned more than a month. Strong flares with gradually increasing intensity are detected in the three spectral ranges during the ~10 days following the Swift trigger. Our optical spectra reveal the presence of a fast outflowing wind, as implied by the detection of a P-Cyg profile (He I - 5876 A) with a terminal velocity of ~2500 km/s. Nebular-like spectra -- with an H_alpha equivalent width of ~500 A -- are also observed. All these features are similar to those seen during the main June 2015 outburst. Thus, the fast optical wind simultaneous with the radio jet is most likely present in every V404 Cyg outburst. Finally, we report on the detection of a strong radio flare in late January 2016, when X-ray and optical monitoring had stopped due to Sun constraints.

We present updated constraints on dark matter models with momentum-dependent or velocity-dependent interactions with nuclei, based on direct detection and solar physics. We improve our previous treatment of energy transport in the solar interior by dark matter scattering, leading to significant changes in fits to many observables. Based on solar physics alone, DM with a spin-independent $q^{4}$ coupling provides the best fit to data, and a statistically satisfactory solution to the solar abundance problem. Once direct detection limits are accounted for however, the best solution is spin-dependent $v^2$ scattering with a reference cross-section of 10$^{-35}$ cm$^2$ (at a reference velocity of $v_0=220$ km s$^{-1}$), and a dark matter mass of about 5 GeV.

Cosmological inflation generates primordial density perturbations on all scales, including those far too small to contribute to the cosmic microwave background. At these scales, isolated ultracompact minihalos of dark matter can form well before standard structure formation, if the perturbations have sufficient amplitude. Minihalos affect pulsar timing data and are potentially bright sources of gamma rays. The resulting constraints significantly extend the observable window of inflation in the presence of cold dark matter, coupling two of the key problems in modern cosmology.

We estimate the sensitivity of the upcoming CTA gamma-ray telescope to DM annihilation at the Galactic centre, improving on previous analyses in a number of significant ways. First, we perform a detailed analyses of all backgrounds, including diffuse astrophysical emission for the first time in a study of this type. Second, we present a statistical framework for including systematic errors and estimate the consequent degradation in sensitivity. These errors may come from e.g. event reconstruction, Monte Carlo determination of the effective area or uncertainty in atmospheric conditions. Third, we show that performing the analysis on a set of suitably optimised regions of interest makes it possible to partially compensate for the degradation in sensitivity caused by systematics and diffuse emission. To probe dark matter with the canonical thermal annihilation cross-section, CTA systematics like non-uniform variations in acceptance over a single field of view must be kept below the 0.3% level, unless the dark matter density rises more steeply in the centre of the Galaxy than predicted by a typical Navarro-Frenk-White or Einasto profile. For a contracted $r^{-1.3}$ profile, and systematics at the 1% level, CTA can probe annihilation to $b\bar{b}$ at the canonical thermal level for dark matter masses between 100 GeV and 10 TeV.

We critically examine recent claims of a high solar metallicity by von Steiger \& Zurbuchen (2016) based on \textit{in situ} measurements of the solar wind, rather than the standard spectroscopically-inferred abundances (Asplund et al. 2009). We test the claim by Vagnozzi et al. (2016) that a composition based on the solar wind enables one to construct a standard solar model in agreement with helioseismological observations and thus solve the decades-old solar modelling problem. We show that, although some helioseismological observables are improved compared to models computed with spectroscopic abundances, most are in fact worse. The high abundance of refractory elements leads to an overproduction of neutrinos, with a predicted $^8$B flux that is nearly twice its observed value, and $^7$Be and CNO fluxes that are experimentally ruled out at high confidence. A combined likelihood analysis shows that models using the vSZ16 abundances fare worse than AGSS09 despite a higher metallicity. We also present astrophysical and spectroscopic arguments showing the vSZ16 composition to be an implausible representation of the solar interior, identifying the first ionisation potential effect in the outer solar atmosphere and wind as the likely culprit.

We combine and extend the analyses of effective scalar, vector, Majorana and Dirac fermion Higgs portal models of dark matter (DM), in which DM couples to the Standard Model (SM) Higgs boson via an operator of the form $\mathcal{O}_{\textrm{DM}}\, H^\dagger H$. For the fermion models, we take an admixture of scalar $\overline{\psi} \psi$ and pseudoscalar $\overline{\psi} i\gamma_5 \psi$ interaction terms. For each model, we apply constraints on the parameter space based on the Planck measured DM relic density and the LHC limits on the Higgs invisible branching ratio. For the first time, we perform a consistent study of the indirect detection prospects for these models based on the WMAP7/Planck observations of the cosmic microwave background, a combined analysis of 15 dwarf spheroidal galaxies by Fermi-LAT and the upcoming Cherenkov Telescope Array (CTA). We also perform a correct treatment of the momentum-dependent direct search cross section that arises from the pseudoscalar interaction term in the fermionic DM theories. We find, in line with previous studies, that current and future direct search experiments such as LUX and XENON1T can exclude much of the parameter space, and we demonstrate that a joint observation in both indirect and direct searches is possible for high mass weakly interacting massive particles. In the case of a pure pseudoscalar interaction of a fermionic DM candidate, future gamma-ray searches are the only class of experiment capable of probing the high mass range of the theory.

We study the effects of energy transport in the Sun by asymmetric dark matter with momentum and velocity-dependent interactions, with an eye to solving the decade-old Solar Abundance Problem. We study effective theories where the dark matter-nucleon scattering cross-section goes as $v_{\rm rel}^{2n}$ and $q^{2n}$ with $n = -1, 0, 1 $ or $2$, where $v_{\rm rel}$ is the dark matter-nucleon relative velocity and $q$ is the momentum exchanged in the collision. Such cross-sections can arise generically as leading terms from the most basic nonstandard DM-quark operators. We employ a high-precision solar simulation code to study the impact on solar neutrino rates, the sound speed profile, convective zone depth, surface helium abundance and small frequency separations. We find that the majority of models that improve agreement with the observed sound speed profile and depth of the convection zone also reduce neutrino fluxes beyond the level that can be reasonably accommodated by measurement and theory errors. However, a few specific points in parameter space yield a significant overall improvement. A 3-5 GeV DM particle with $\sigma_{SI} \propto q^2$ is particularly appealing, yielding more than a $6\sigma$ improvement with respect to standard solar models, while being allowed by direct detection and collider limits. We provide full analytical capture expressions for $q$- and $v_{\rm rel}$-dependent scattering, as well as complete likelihood tables for all models.

We present an improved event-level likelihood formalism for including neutrino telescope data in global fits to new physics. We derive limits on spin-dependent dark matter-proton scattering by employing the new formalism in a re-analysis of data from the 79-string IceCube search for dark matter annihilation in the Sun, including explicit energy information for each event. The new analysis excludes a number of models in the weak-scale minimal supersymmetric standard model (MSSM) for the first time. This work is accompanied by the public release of the 79-string IceCube data, as well as an associated computer code for applying the new likelihood to arbitrary dark matter models.

One of the simplest models of dark matter is that where a scalar singlet field S comprises some or all of the dark matter, and interacts with the standard model through an HHSS coupling to the Higgs boson. We update the present limits on the model from LHC searches for invisible Higgs decays, the thermal relic density of S, and dark matter searches via indirect and direct detection. We point out that the currently allowed parameter space is on the verge of being significantly reduced with the next generation of experiments. We discuss the impact of such constraints on possible applications of scalar singlet dark matter, including a strong electroweak phase transition, and the question of vacuum stability of the Higgs potential at high scales.

Cosmic strings are expected to form loops. These can act as seeds for accretion of dark matter, leading to the formation of ultracompact minihalos (UCMHs). We perform a detailed study of the accretion of dark matter onto cosmic string loops and compute the resulting mass distribution of UCMHs. We then apply observational limits on the present-day abundance of UCMHs to derive corresponding limits on the cosmic string tension $G\mu$. The bounds are strongly dependent upon the assumed distribution of loop velocities and their impacts on UCMH formation. Under the assumption that a loop can move up to a thousand times its own radius and still form a UCMH, we find a limit of $G\mu\le 5\times10^{-8}$. We show, in opposition to previous results, that strong limits on the cosmic string tension are not obtainable from UCMHs when more stringent (and realistic) requirements are placed on loop velocities.

We present observations and analysis of a sample of 123 galaxy clusters from the 2013 Planck catalogue of Sunyaev-Zel'dovich sources with the Arcminute Microkelvin Imager (AMI), a ground-based radio interferometer. AMI provides an independent measurement with higher angular resolution, 3 arcmin compared to the Planck beams of 5-10 arcmin. The AMI observations thus provide validation of the cluster detections, improved positional estimates, and a consistency check on the fitted 'size' ($\theta_{s}$) and 'flux' ($Y_{\rm tot}$) parameters in the Generalised Navarro, Frenk and White (GNFW) model. We detect 99 of the clusters. We use the AMI positional estimates to check the positional estimates and error-bars produced by the Planck algorithms PowellSnakes and MMF3. We find that $Y_{\rm tot}$ values as measured by AMI are biased downwards with respect to the Planck constraints, especially for high Planck-SNR clusters. We perform simulations to show that this can be explained by deviation from the 'universal' pressure profile shape used to model the clusters. We show that AMI data can constrain the $\alpha$ and $\beta$ parameters describing the shape of the profile in the GNFW model for individual clusters provided careful attention is paid to the degeneracies between parameters, but one requires information on a wider range of angular scales than are present in AMI data alone to correctly constrain all parameters simultaneously.

QUIJOTE (Q-U-I JOint TEnerife) is a new polarimeter aimed to characterize the polarization of the Cosmic Microwave Background and other Galactic and extragalactic signals at medium and large angular scales in the frequency range 10-40 GHz. The multi-frequency (10-20~GHz) instrument, mounted on the first QUIJOTE telescope, saw first light on November 2012 from the Teide Observatory (2400~m a.s.l). During 2014 the second telescope has been installed at this observatory. A second instrument at 30~GHz will be ready for commissioning at this telescope during summer 2015, and a third additional instrument at 40~GHz is now being developed. These instruments will have nominal sensitivities to detect the B-mode polarization due to the primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r=0.05.

Broad disagreement persists between helioseismological observables and predictions of solar models computed with the latest surface abundances. Here we show that most of these problems can be solved by the presence of asymmetric dark matter coupling to nucleons as the square of the momentum $q$ exchanged in the collision. We compute neutrino fluxes, small frequency separations, surface helium abundances, sound speed profiles and convective zone depths for a number of models, showing more than a $6\sigma$ preference for $q^2$ models over others, and over the Standard Solar Model. The preferred mass (3\,GeV) and reference dark matter-nucleon cross-section ($10^{-37}$\,cm$^2$ at $q_0 = 40$\,MeV) are within the region of parameter space allowed by both direct detection and collider searches.

Anomalous microwave emission (AME) has been observed in numerous sky regions, in the frequency range ~10-60 GHz. One of the most scrutinized regions is G159.6-18.5, located within the Perseus molecular complex. In this paper we present further observations of this region (194 hours in total over ~250 deg^2), both in intensity and in polarization. They span four frequency channels between 10 and 20 GHz, and were gathered with QUIJOTE, a new CMB experiment with the goal of measuring the polarization of the CMB and Galactic foregrounds. When combined with other publicly-available intensity data, we achieve the most precise spectrum of the AME measured to date, with 13 independent data points being dominated by this emission. The four QUIJOTE data points provide the first independent confirmation of the downturn of the AME spectrum at low frequencies, initially unveiled by the COSMOSOMAS experiment in this region. We accomplish an accurate fit of these data using models based on electric dipole emission from spinning dust grains, and also fit some of the parameters on which these models depend. We also present polarization maps with an angular resolution of ~1 deg and a sensitivity of ~25 muK/beam. From these maps, which are consistent with zero polarization, we obtain upper limits of Pi<6.3% and <2.8% (95% C.L.) respectively at 12 and 18 GHz, a frequency range where no AME polarization observations have been reported to date. These constraints are compatible with theoretical predictions of the polarization fraction from electric dipole emission originating from spinning dust grains. At the same time, they rule out several models based on magnetic dipole emission from dust grains ordered in a single magnetic domain, which predict higher polarization levels. Future QUIJOTE data in this region may allow more stringent constraints on the polarization level of the AME.

(abridged) The Galactic Center is one of the most promising targets for indirect detection of dark matter with gamma rays. We investigate the sensitivity of the upcoming Cherenkov Telescope Array (CTA) to dark matter annihilation and decay in the Galactic Center. As the inner density profile of the Milky Way's dark matter halo is uncertain, we study the impact of the slope of the Galactic density profile, inwards of the Sun, on the prospects for detecting a dark matter signal with CTA. We find that the sensitivity achieved by CTA to annihilation signals is strongly dependent on the inner profile slope, whereas the dependence is more mild in the case of dark matter decay. Surprisingly, we find that the optimal choice of signal and background regions is virtually independent of the assumed density profile. For the fiducial case of a Navarro-Frenk-White profile, we find that CTA will be able to probe annihilation cross sections well below the canonical thermal relic value for dark matter masses from a few tens of GeV up to $\sim 5$ TeV for annihilation to $\tau^{+}\tau^{-}$, and will achieve only a slightly weaker sensitivity for annihilation to $b\bar{b}$ or $\mu^{+}\mu^{-}$. CTA will improve significantly on current sensitivity to annihilation signals for dark matter masses above $\sim 100$ GeV, covering parameter space that is complementary to that probed by searches with the Fermi Large Area Telescope. The interpretation of apparent excesses in the measured cosmic-ray electron and positron spectra as signals of dark matter decay will also be testable with CTA observations of the Galactic Center. We demonstrate that both for annihilation and for decay, including spectral information for hard channels (such as $\mu^{+}\mu^{-}$ and $\tau^{+}\tau^{-}$) leads to enhanced sensitivity for dark matter masses above $m_{\rm DM}\sim 200$ GeV.

We re-evaluate the abundances of the elements in the Sun from copper ($Z=29$) to thorium ($Z=90$). Our results are mostly based on neutral and singly-ionised lines in the solar spectrum. We use the latest 3D hydrodynamic solar model atmosphere, and in a few cases also correct for departures from local thermodynamic equilibrium (LTE) using non-LTE (NLTE) calculations performed in 1D. In order to minimise statistical and systematic uncertainties, we make stringent line selections, employ the highest-quality observational data and carefully assess oscillator strengths, hyperfine constants and isotopic separations available in the literature, for every line included in our analysis. Our results are typically in good agreement with the abundances in the most pristine meteorites, but there are some interesting exceptions. This analysis constitutes both a full exposition and a slight update of the relevant parts of the preliminary results we presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), including full line lists and details of all input data that we have employed.

The composition of the Sun is an essential piece of reference data for astronomy, cosmology, astroparticle, space and geo-physics. This article, dealing with the intermediate-mass elements Na to Ca, is the first in a series describing the comprehensive re-determination of the solar composition. In this series we severely scrutinise all ingredients of the analysis across all elements, to obtain the most accurate, homogeneous and reliable results possible. We employ a highly realistic 3D hydrodynamic solar photospheric model, which has successfully passed an arsenal of observational diagnostics. To quantify systematic errors, we repeat the analysis with three 1D hydrostatic model atmospheres (MARCS, MISS and Holweger & M\"{u}ller 1974) and a horizontally and temporally-averaged version of the 3D model ($\langle$3D$\rangle$). We account for departures from LTE wherever possible. We have scoured the literature for the best transition probabilities, partition functions, hyperfine and other data, and stringently checked all observed profiles for blends. Our final 3D+NLTE abundances are: $\log\epsilon_{\mathrm{Na}}=6.21\pm0.04$, $\log\epsilon_{\mathrm{Mg}}=7.59\pm0.04$, $\log\epsilon_{\mathrm{Al}}=6.43\pm0.04$, $\log\epsilon_{\mathrm{Si}}=7.51\pm0.03$, $\log\epsilon_{\mathrm{P}}=5.41\pm0.03$, $\log \epsilon_{\mathrm{S}}=7.13\pm0.03$, $\log\epsilon_{\mathrm{K}}=5.04\pm0.05$ and $\log\epsilon_{\mathrm{Ca}}=6.32\pm0.03$. The uncertainties include both statistical and systematic errors. Our results are systematically smaller than most previous ones with the 1D semi-empirical Holweger & M\"uller model. The $\langle$3D$\rangle$ model returns abundances very similar to the full 3D calculations. This analysis provides a complete description and a slight update of the Na to Ca results presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), with full details of all lines and input data.

We redetermine the abundances of all iron group nuclei in the Sun, based on neutral and singly-ionised lines of Sc, Ti, V, Mn, Fe, Co and Ni in the solar spectrum. We employ a realistic 3D hydrodynamic model solar atmosphere, corrections for departures from local thermodynamic equilibrium (NLTE), stringent line selection procedures and high quality observational data. We have scoured the literature for the best quality oscillator strengths, hyperfine constants and isotopic separations available for our chosen lines. We find $\log \epsilon_\mathrm{Sc}=3.16\pm0.04$, $\log \epsilon_\mathrm{Ti}=4.93\pm0.04$, $\log \epsilon_\mathrm{V}=3.89\pm0.08$, $\log \epsilon_\mathrm{Cr}=5.62\pm0.04$, $\log \epsilon_\mathrm{Mn}=5.42\pm0.04$, $\log \epsilon_\mathrm{Fe}=7.47\pm0.04$, $\log \epsilon_\mathrm{Co}=4.93\pm0.05$ and $\log \epsilon_\mathrm{Ni}=6.20\pm0.04$. Our uncertainties factor in both statistical and systematic errors (the latter estimated for possible errors in the model atmospheres and NLTE line formation). The new abundances are generally in good agreement with the CI meteoritic abundances but with some notable exceptions. This analysis constitutes both a full exposition and a slight update of the preliminary results we presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), including full line lists and details of all input data we employed.

We use the formalism of Gould and Raffelt to compute the dimensionless thermal conduction coefficients for scattering of dark matter particles with standard model nucleons via cross-sections that depend on the relative velocity or momentum exchanged between particles. Motivated by models invoked to reconcile various recent results in direct detection, we explicitly compute the conduction coefficients $\alpha$ and $\kappa$ for cross-sections that go as $v_{\rm rel}^2$, $v_{\rm rel}^4$, $v_{\rm rel}^{-2}$, $q^2$, $q^4$ and $q^{-2}$, where $v_{\rm rel}$ is the relative DM-nucleus velocity and $q$ is the momentum transferred in the collision. We find that a $v_{\rm rel}^{-2}$ dependence can significantly enhance energy transport from the inner solar core to the outer core. The same can true for any $q$-dependent coupling, if the dark matter mass lies within some specific range for each coupling. This effect can complement direct searches for dark matter; combining these results with state-of-the-art Solar simulations should greatly increase sensitivity to certain DM models. It also seems possible that the so-called Solar Abundance Problem could be resolved by enhanced energy transport in the solar core due to such velocity- or momentum-dependent scatterings.

The QUIJOTE (Q-U-I JOint Tenerife) CMB Experiment is designed to observe the polarization of the Cosmic Microwave Background and other Galactic and extragalactic signals at medium and large angular scales in the frequency range of 10-40 GHz. The first of the two QUIJOTE telescopes and the multi-frequency (10-20 GHz) instrument have been in operation since November 2012. In 2014 a second telescope and a new instrument at 30GHz will be ready for commissioning, and an additional instrument at 40 GHz is in its final design stages. After three years of effective observations, the data obtained by these telescopes and instruments will have the required sensitivity to detect a primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r = 0.05. At the moment, we have completed half of the wide Galactic survey with the multi-frequency instrument covering 18 000 square degrees of the Northern hemisphere. When we finish this survey in early 2014, we shall have reached approximately 14{\mu}K per one degree beam at 11, 13, 17 and 19 GHz, in both Q and U.

We describe a general procedure for using number counts of any object to constrain the probability distribution of the primordial fluctuations, allowing for generic weak non-Gaussianity. We apply this procedure to use limits on the abundance of primordial black holes and dark matter ultracompact minihalos (UCMHs) to characterize the allowed statistics of primordial fluctuations on very small scales. We present constraints on the power spectrum and the amplitude of the skewness for two different families of non-Gaussian distributions, distinguished by the relative importance of higher moments. Although primordial black holes probe the smallest scales, ultracompact minihalos provide significantly stronger constraints on the power spectrum and so are more likely to eventually provide small-scale constraints on non-Gaussianity.

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 AMI-detected clusters, we calculate large-scale 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 XCS-based core-temperature 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 uv-data 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 uv-data 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.

The cosmic microwave background provides constraints on the annihilation and decay of light dark matter at redshifts between 100 and 1000, the strength of which depends upon the fraction of energy ending up in the form of electrons and photons. The resulting constraints are usually presented for a limited selection of annihilation and decay channels. Here we provide constraints on the annihilation cross section and decay rate, at discrete values of the dark matter mass, for all the annihilation and decay channels whose secondary spectra have been computed using PYTHIA in arXiv:1012.4515 ("PPPC4DMID: A Poor Particle Physicist Cookbook for Dark Matter Indirect Detection"), namely e, mu, tau, V -> e, V -> mu, V -> tau, u, d, s, c, b, t, gamma, g, W, Z and h. By interpolating in mass, these can be used to find the CMB constraints and likelihood functions from WMAP7 and Planck for a wide range of dark matter models, including those with annihilation or decay into a linear combination of different channels.

We analyse the sensitivity of IceCube-DeepCore to annihilation of neutralino dark matter in the solar core, generated within a 25 parameter version of the minimally supersymmetric standard model (MSSM-25). We explore the 25-dimensional parameter space using scanning methods based on importance sampling and using DarkSUSY 5.0.6 to calculate observables. Our scans produced a database of 6.02 million parameter space points with neutralino dark matter consistent with the relic density implied by WMAP 7-year data, as well as with accelerator searches. We performed a model exclusion analysis upon these points using the expected capabilities of the IceCube-DeepCore Neutrino Telescope. We show that IceCube-DeepCore will be sensitive to a number of models that are not accessible to direct detection experiments such as SIMPLE, COUPP and XENON100, indirect detection using Fermi-LAT observations of dwarf spheroidal galaxies, nor to current LHC searches.

Dark halo substructure may reveal itself through secondary, small-scale gravitational lensing effects on light sources that are macrolensed by a foreground galaxy. Here, we explore the prospects of using Very Long Baseline Interferometry (VLBI) observations of multiply-imaged quasar jets to search for submilliarcsecond-scale image distortions produced by various forms of dark substructures in the 1e3-1e8 Msolar mass range. We present lensing simulations relevant for the angular resolutions attainable with the existing European VLBI Network (EVN), the global VLBI array, and an upcoming observing mode in which the Atacama Large Millimeter Array (ALMA) is connected to the global VLBI array. While observations of this type would not be sensitive to standard cold dark matter subhalos, they can be used to detect more compact forms of halo substructure predicted in alternative structure formation scenarios. By mapping ~5 strongly lensed systems, it should be possible to detect or robustly rule out primordial black holes in the 1e3-1e6 Msolar mass range if they constitute >1% percent of the dark matter in these lenses. Ultracompact minihalos are harder to detect using this technique, but 1e6-1e8 Msolar ultracompact minihalos could in principle be detected if they constitute >10% of the dark matter.

For a Gaussian spectrum of primordial density fluctuations, ultracompact minihalos (UCMHs) of dark matter are expected to be produced in much greater abundance than, e.g., primordial black holes. Forming shortly after matter-radiation equality, these objects would develop very dense and spiky dark matter profiles. In the standard scenario where dark matter consists of thermally-produced, weakly-interacting massive particles, UCMHs could thus appear as highly luminous gamma-ray sources, or leave an imprint in the cosmic microwave background by changing the reionisation history of the Universe. We derive corresponding limits on the cosmic abundance of UCMHs at different epochs, and translate them into constraints on the primordial power spectrum. We find the resulting constraints to be quite severe, especially at length scales much smaller than what can be directly probed by the cosmic microwave background or large-scale structure observations. We use our results to provide an updated compilation of the best available constraints on the power of density fluctuations on all scales, ranging from the present-day horizon to scales more than 20 orders of magnitude smaller.

We present 16-GHz Sunyaev-Zel'dovich observations using the Arcminute Microkelvin Imager (AMI) and subsequent Bayesian analysis of six galaxy clusters at redshift ($z \approx 1$) chosen from an X-ray and Infrared selected sample from Culverhouse et al. (2010). In the subsequent analysis we use two cluster models, an isothermal \beta-model and a Dark Matter GNFW (DM-GNFW) 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 low-radius X-ray 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 small-radius estimates introduce to large-radius SZ predictions. In addition, we follow-up the two detections with deep, single-pointed LA observations. We find a 3 sigma tentative decrement toward CLJ1415+3612 at high-resolution and a 5 sigma high-resolution decrement towards XMJ0830+5241.

It is well known that current spectroscopic determinations of the chemical composition of the Sun are starkly at odds with the metallicity implied by helioseismology. We investigate whether the discrepancy may be due to conversion of photons to a new light boson in the solar photosphere. We examine the impact of particles with axion-like interactions with the photon on the inferred photospheric abundances, showing that resonant axion-photon conversion is not possible in the region of the solar atmosphere in which line-formation occurs. Although non-resonant conversion in the line-forming regions can in principle impact derived abundances, constraints from axion-photon conversion experiments rule out the couplings necessary for these effects to be detectable. We show that this extends to hidden photons and chameleons (which would exhibit similar phenomenological behaviour), ruling out known theories of new light bosons as photospheric solutions to the solar abundance problem.

The first population III stars are predicted to form in minihalos at a redshift of approximately 10-30. The James Webb Space Telescope (JWST), tentatively scheduled for launch in 2018, will probably be able to detect some of the first galaxies, but whether it will also be able to detect the first stars remains more doubtful. Here, we explore the prospects of detecting an isolated population III star or a small cluster of population III stars down to redshift 2 in either lensed or unlensed fields. Our calculations are based on realistic stellar atmospheres and take into account the potential flux contribution from the surrounding HII region. We find that unlensed population III stars are beyond the reach of JWST, and that even lensed population III stars will be extremely difficult to detect. However, the main problem with the latter approach is not necessarily that the lensed stars are too faint, but that their surface number densities are too low. To detect even one population III star of 60 solar masses when pointing JWST through the galaxy cluster MACS J0717.5+3745, the lensing cluster with the largest Einstein radius detected so far, the cosmic star formation rate of population III stars would need to be approximately an order of magnitude higher than predicted by the most optimistic current models.

The AMI Galactic Plane Survey (AMIGPS) is a large area survey of the outer Galactic plane to provide arcminute resolution images at milli-Jansky sensitivity in the centimetre-wave 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 centimetre-wave frequencies greater than 1.4 GHz.

We present a fast likelihood method for including event-level neutrino telescope data in parameter explorations of theories for new physics, and announce its public release as part of DarkSUSY 5.0.6. Our construction includes both angular and spectral information about neutrino events, as well as their total number. We also present a corresponding measure for simple model exclusion, which can be used for single models without reference to the rest of a parameter space. We perform a number of supersymmetric parameter scans with IceCube data to illustrate the utility of the method: example global fits and a signal recovery in the constrained minimal supersymmetric standard model (CMSSM), and a model exclusion exercise in a 7-parameter phenomenological version of the MSSM. The final IceCube detector configuration will probe almost the entire focus-point region of the CMSSM, as well as a number of MSSM-7 models that will not otherwise be accessible to e.g. direct detection. Our method accurately recovers the mock signal, and provides tight constraints on model parameters and derived quantities. We show that the inclusion of spectral information significantly improves the accuracy of the recovery, providing motivation for its use in future IceCube analyses.

The Arcminute Microkelvin Imager (AMI) is a telescope specifically designed for high sensitivity measurements of low-surface-brightness features at cm-wavelength and has unique, important capabilities. It consists of two interferometer arrays operating over 13.5-18 GHz that image structures on scales of 0.5-10 arcmin with very low systematics. The Small Array (AMI-SA; ten 3.7-m antennas) couples very well to Sunyaev-Zel'dovich features from galaxy clusters and to many Galactic features. The Large Array (AMI-LA; eight 13-m antennas) has a collecting area ten times that of the AMI-SA and longer baselines, crucially allowing the removal of the effects of confusing radio point sources from regions of low surface-brightness, extended emission. Moreover AMI provides fast, deep object surveying and allows monitoring of large numbers of objects. In this White Paper we review the new science - both Galactic and extragalactic - already achieved with AMI and outline the prospects for much more.

We present 16-GHz 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 5-23sigma peak surface brightnesses. Cluster parameters are obtained using a Bayesian cluster analysis. We fit isothermal beta-models 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 X-ray 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.0-2.5xM_T(r500).(iii) Where we have found M_T X-ray (X) and weak-lensing (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 X-ray 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 X-ray masses. We stress that large-r 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 X-ray signal (and T_X) and, even at high r, mergers affect n^2-weighted X-ray data more than n-weighted SZ, implying significant shocking or clumping or both occur even in the outer parts of mergers.

We present 16-GHz 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 X-ray data. We use a fast, Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources 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 beta-model to our data and assumes the cluster follows the theoretical mass-temperature relation. Large-scale cluster parameters internal to r_{500} are derived under the assumption of hydrostatic equilibrium. Posterior distributions for the large-scale parameters of 8 of our clusters are given; SZ effects towards Abell 1704 and Zw0857.9+2107 were not detected and our spherical beta-profile was found to be an inadequate fit to the decrement on our map for Abell 2409.

Events like inflation or phase transitions can produce large density perturbations on very small scales in the early Universe. Probes of small scales are therefore useful for e.g. discriminating between inflationary models. Until recently, the only such constraint came from non-observation of primordial black holes (PBHs), associated with the largest perturbations. Moderate-amplitude perturbations can collapse shortly after matter-radiation equality to form ultracompact minihalos (UCMHs) of dark matter, in far greater abundance than PBHs. If dark matter self-annihilates, UCMHs become excellent targets for indirect detection. Here we discuss the gamma-ray fluxes expected from UCMHs, the prospects of observing them with gamma-ray telescopes, and limits upon the primordial power spectrum derived from their non-observation by the Fermi Large Area Space Telescope.

We present an interesting Sunyaev-Zel'dovich (SZ) detection in the first of the Arcminute Microkelvin Imager (AMI) 'blind', degree-square fields to have been observed down to our target sensitivity of 100{\mu}Jy/beam. In follow-up deep pointed observations the SZ effect is detected with a maximum peak decrement greater than 8 \times the thermal noise. No corresponding emission is visible in the ROSAT all-sky X-ray survey and no cluster is evident in the Palomar all-sky optical survey. Compared with existing SZ images of distant clusters, the extent is large (\approx 10') and complex; our analysis favours a model containing two clusters rather than a single cluster. Our Bayesian analysis is currently limited to modelling each cluster with an ellipsoidal or spherical beta-model, which do not do justice to this decrement. Fitting an ellipsoid to the deeper candidate we find the following. (a) Assuming that the Evrard et al. (2002) approximation to Press & Schechter (1974) correctly gives the number density of clusters as a function of mass and redshift, then, in the search area, the formal Bayesian probability ratio of the AMI detection of this cluster is 7.9 \times 10^4:1; alternatively assuming Jenkins et al. (2001) as the true prior, the formal Bayesian probability ratio of detection is 2.1 \times 10^5:1. (b) The cluster mass is MT,200 = 5.5+1.2\times 10^14h-1M\odot. (c) Abandoning a physical model with num- -1.3 70 ber density prior and instead simply modelling the SZ decrement using a phenomenological {\beta}-model of temperature decrement as a function of angular distance, we find a central SZ temperature decrement of -295+36 {\mu}K - this allows for CMB primary anisotropies, receiver -15 noise and radio sources. We are unsure if the cluster system we observe is a merging system or two separate clusters.

We present 16 GHz (1.9 cm) deep radio continuum observations made with the Arcminute Microkelvin Imager (AMI) of a sample of low-mass 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 sub-mm 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 free-free emission. We find an average spectral index in both Td scenarios consistent with free-free 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.

We discuss irreducible statistical limitations of future ton-scale dark matter direct detection experiments. We focus in particular on the coverage of confidence intervals, which quantifies the reliability of the statistical method used to reconstruct the dark matter parameters, and the bias of the reconstructed parameters. We study 36 benchmark dark matter models within the reach of upcoming ton-scale experiments. We find that approximate confidence intervals from a profile-likelihood analysis exactly cover or over-cover the true values of the WIMP parameters, and are hence conservative. We evaluate the probability that unavoidable statistical fluctuations in the data might lead to a biased reconstruction of the dark matter parameters, or large uncertainties on the reconstructed parameter values. We show that this probability can be surprisingly large, even for benchmark models leading to a large event rate of order a hundred counts. We find that combining data sets from two different targets leads to improved coverage properties, as well as a substantial reduction of statistical bias and uncertainty on the dark matter parameters.

We present an analysis of observations made with the Arcminute Microkelvin Imager (AMI) and the Canada-France-Hawaii Telescope (CFHT) of six galaxy clusters in a redshift range of 0.16--0.41. The cluster gas is modelled using the Sunyaev--Zel'dovich (SZ) data provided by AMI, while the total mass is modelled using the lensing data from the CFHT. In this paper, we: i) find very good agreement between SZ measurements (assuming large-scale virialisation and a gas-fraction prior) and lensing measurements of the total cluster masses out to r_200; ii) perform the first multiple-component weak-lensing analysis of A115; iii) confirm the unusual separation between the gas and mass components in A1914; iv) jointly analyse the SZ and lensing data for the relaxed cluster A611, confirming our use of a simulation-derived mass-temperature relation for parameterizing measurements of the SZ effect.

Most Sunyaev--Zel'dovich (SZ) and X-ray 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 beta-model 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 sub-mm 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 long-wavelength 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 sub-mm data to separate the contributions from thermal dust and radio emission at 1.8 cm and derive disk masses directly from the cm-wave 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 sub-mm 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.

The identity of dark matter is one of the key outstanding problems in both particle and astrophysics. In this thesis, I describe a number of complementary searches for particle dark matter. I discuss how the impact of dark matter on stars can constrain its interaction with nuclei, focussing on main sequence stars close to the Galactic Centre, and on the first stars as seen through the upcoming James Webb Space Telescope. The mass and annihilation cross-section of dark matter particles can be probed with searches for gamma rays produced in astronomical targets. Dwarf galaxies and ultracompact, primordially-produced dark matter minihalos turn out to be especially promising in this respect. I illustrate how the results of these searches can be combined with constraints from accelerators and cosmology to produce a single global fit to all available data. Global fits in supersymmetry turn out to be quite technically demanding, even with the simplest predictive models and the addition of complementary data from a bevy of astronomical and terrestrial experiments; I show how genetic algorithms can help in overcoming these challenges.

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 (AMI-LA). 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 inter-epoch 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 perform a detailed and systematic investigation of the possible impacts of dark stars upon the reionization history of the Universe, and its signatures in the cosmic microwave background (CMB). We compute hydrogen reionization histories, CMB optical depths and anisotropy power spectra for a range of stellar populations including dark stars. If dark stars capture large amounts of dark matter via nuclear scattering, reionization can be substantially delayed, leading to decreases in the integrated optical depth to last scattering and large-scale power in the EE polarization power spectrum. Using the integrated optical depth observed by WMAP7, in our canonical reionization model we rule out the section of parameter space where dark stars with high scattering-induced capture rates tie up more than ~90% of all the first star-forming baryons, and live for over ~250 Myr. When nuclear scattering delivers only moderate amounts of dark matter, reionization can instead be sped up slightly, modestly increasing the CMB optical depth. If dark stars do not obtain any dark matter via nuclear scattering, effects upon reionization and the CMB are negligible. The effects of dark stars upon reionization and its CMB markers can be largely mimicked or compensated for by changes in the existing parameters of reionization models, making dark stars difficult to disentangle from astrophysical uncertainties, but also widening the range of standard parameters in reionization models that can be made consistent with observations.

In order to place limits on dark matter (DM) properties using $\gamma$-ray observations, previous analyses have often assumed a very simple parametrisation of the $\gamma$-ray annihilation yield; typically, it has been assumed that annihilation proceeds through a single channel only. In realistic DM models, annihilation may occur into many different final states, making this quite a rough ansatz. With additional processes like virtual internal bremsstrahlung and final state radiation, this ansatz becomes even more incorrect, and the need for scans of explicit model parameter spaces becomes clear. Here we present scans of the parameter space of the Constrained Minimal Supersymmetric Standard Model (CMSSM), considering $\gamma$-ray spectra from three dwarf galaxies, the Galactic Centre region and the broader Galactic halo, as obtained with the High-Energy Stereoscopic System (H.E.S.S.). We present a series of likelihood scans combining the H.E.S.S. data with other experimental results. We show that observations of the Sagittarius, Carina and Sculptor dwarf galaxies disfavour the coannihilation region of the CMSSM and models with large annihilation cross-sections. This is true even under reasonable assumptions about the DM density profiles, and constitutes the strongest constraint to date on coannihilation models within the CMSSM. The Galactic halo has a similar, but weaker, effect. The Galactic Centre search is complicated by a strong (unknown) $\gamma$-ray source, and we see no effect on the CMSSM parameter space when assuming a realistic Galactic Centre DM density profile.

Models of weak-scale supersymmetry offer viable dark matter (DM) candidates. Their parameter spaces are however rather large and complex, such that pinning down the actual parameter values from experimental data can depend strongly on the employed statistical framework and scanning algorithm. In frequentist parameter estimation, a central requirement for properly constructed confidence intervals is that they cover true parameter values, preferably at exactly the stated confidence level when experiments are repeated infinitely many times. Since most widely-used scanning techniques are optimised for Bayesian statistics, one needs to assess their abilities in providing correct confidence intervals in terms of the statistical coverage. Here we investigate this for the Constrained Minimal Supersymmetric Standard Model (CMSSM) when only constrained by data from direct searches for dark matter. We construct confidence intervals from one-dimensional profile likelihoods and study the coverage by generating several pseudo-experiments for a few benchmark sets of pseudo-true parameters. We use nested sampling to scan the parameter space and evaluate the coverage for the benchmarks when either flat or logarithmic priors are imposed on gaugino and scalar mass parameters. The sampling algorithm has been used in the configuration usually adopted for exploration of the Bayesian posterior. We observe both under- and over-coverage, which in some cases vary quite dramatically when benchmarks or priors are modified. We show how most of the variation can be explained as the impact of explicit priors as well as sampling effects, where the latter are indirectly imposed by physicality conditions. For comparison, we also evaluate the coverage for Bayesian credible intervals, and observe significant under-coverage in those cases.

Weakly interacting massive particles (WIMPs) are amongst the most interesting dark matter (DM) candidates. Many DM candidates naturally arise in theories beyond the standard model (SM) of particle physics, like weak-scale supersymmetry (SUSY). Experiments aim to detect WIMPs by scattering, annihilation or direct production, and thereby determine the underlying theory to which they belong, along with its parameters. Here we examine the prospects for further constraining the Constrained Minimal Supersymmetric Standard Model (CMSSM) with future ton-scale direct detection experiments. We consider ton-scale extrapolations of three current experiments: CDMS, XENON and COUPP, with 1000 kg-years of raw exposure each. We assume energy resolutions, energy ranges and efficiencies similar to the current versions of the experiments, and include backgrounds at target levels. Our analysis is based on full likelihood constructions for the experiments. We also take into account present uncertainties on hadronic matrix elements for neutralino-quark couplings, and on halo model parameters. We generate synthetic data based on four benchmark points and scan over the CMSSM parameter space using nested sampling. We construct both Bayesian posterior PDFs and frequentist profile likelihoods for the model parameters, as well as the mass and various cross-sections of the lightest neutralino. Future ton-scale experiments will help substantially in constraining supersymmetry, especially when results of experiments primarily targeting spin-dependent nuclear scattering are combined with those directed more toward spin-independent interactions.

We present follow-up observations of two galaxy clusters detected blindly via the Sunyaev-Zel'dovich (SZ) effect and released in the Planck Early Release Compact Source Catalogue. We use the Arcminute Microkelvin Imager, a dual-array 14-18 GHz radio interferometer. After radio source subtraction, we find a SZ decrement of integrated flux density -1.08+/-0.10 mJy toward PLCKESZ G121.11+57.01, and improve the position measurement of the cluster, finding the centre to be RA 12 59 36.4, Dec +60 04 46.8, to an accuracy of 20 arcseconds. The region of PLCKESZ G115.71+17.52 contains strong extended emission, so we are unable to confirm the presence of this cluster via the SZ effect.

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 (AMI-LA). 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 have observed an area of approximatley 27 deg^2 to an rms noise level of less than 0.2 mJy at 15.7 GHz, using the Arcminute Microkelvin Imager Large Array. These observations constitute the most sensitive radio-source survey of any extent (greater than approximately 0.2 deg^2) above 1.4 GHz. This paper presents the techniques employed for observing, mapping and source extraction. We have used a systematic procedure for extracting information and producing source catalogues, from maps with varying noise and uv-coverage. We have performed simulations to test our mapping and source-extraction procedures, and developed methods for identifying extended, overlapping and spurious sources in noisy images. In an accompanying paper, AMI Consortium: Davies et al. 2010, the first results from the 10C survey, including the deep 15.7-GHz source count, are presented.

The first results from the Tenth Cambridge (10C) Survey of Radio Sources, carried out using the AMI Large Array (LA) at an observing frequency of 15.7 GHz, are presented. The survey fields cover an area of approximately 27 sq. degrees to a flux-density completeness of 1 mJy. Results for some deeper areas, covering approximately 12 sq. degrees, wholly contained within the total areas and complete to 0.5 mJy, are also presented. The completeness for both areas is estimated to be at least 93 per cent. The source catalogue contains 1897 entries and is available at www.mrao.cam.ac.uk/surveys/10C. It has been combined with that of the 9C Survey to calculate the 15.7-GHz source counts. A broken power law is found to provide a good parameterisation of the differential count between 0.5 mJy and 1 Jy. The measured count has been compared to that predicted by de Zotti et al. (2005). The model displays good agreement with the data at the highest flux densities but under-predicts the integrated count between 0.5 mJy and 1 Jy by about 30 per cent. Entries from the source catalogue have been matched to those contained in the catalogues of NVSS and FIRST (both of which have observing frequencies of 1.4 GHz). This matching provides evidence for a shift in the typical 1.4-to-15.7-GHz spectral index of the 15.7-GHz-selected source population with decreasing flux density towards sub-mJy levels - the spectra tend to become less steep. Automated methods for detecting extended sources have been applied to the data; approximately 5 per cent of the sources are found to be extended relative to the LA synthesised beam of approximately 30 arcsec. Investigations using higher-resolution data showed that most of the genuinely extended sources at 16 GHz are classical doubles, although some nearby galaxies and twin-jet sources were also identified.

We present 13.9-18.2 GHz observations of the Sunyaev-Zel'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 X-ray images suggests that they both have extended regions which lie approximately perpendicular to one another, with their emission peaks significantly displaced. These features indicate non-uniformities in the distributions of the gas temperature and pressure, indicative of a cluster merger. We use a Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources 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 beta-model, despite the inadequacy of this model for a merging system like Abell 2146, and assume the cluster follows the mass-temperature relation of a virialized, singular, isothermal sphere. We note that this model avoids inferring large-scale 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.

Isolated population III stars are postulated to exist at approximately z=10-30 and may attain masses up to a few hundred solar masses. The James Webb Space telescope (JWST) is the next large space based infrared telescope and is scheduled for launch in 2014. Using a 6.5 meter primary mirror, it will probably be able to detect some of the first galaxies forming in the early Universe. A natural question is whether it will also be able to see any isolated population III stars. Here, we calculate the apparent broadband AB-magnitudes for 300 solar masses population III stars in JWST filters at z=10-20. Our calculations are based on realistic stellar atmospheres and take into account the potential flux contribution from the surrounding HII region. The gravitational magnification boost achieved when pointing JWST through a foreground galaxy cluster is also considered. Using this machinery, we derive the conditions required for JWST to be able to detect population III stars in isolation. We find that a detection of individual population III stars with JWST is unlikely at these redshifts. However, the main problem is not necessarily that these stars are too faint, once gravitational lensing is taken into account, but that their surface number densities are too low.

The most compelling and popular models for dark matter predict that it should congregate and annihilate in stellar cores. Stars where annihilation contributes substantially to the total energy budget look very different to those with which we are familiar. Here I explain the general features of stars modified by dark matter annihilation with the help of a series of grids of 'dark' stellar evolutionary models, and describe the public code with which they were computed. I go on to discuss possible impacts of dark stars on the high-redshift Universe, including the history of reionisation. The preliminary reionisation calculations reproduced here are based on dedicated models for dark star atmospheres, and for the stellar populations to which dark stars would belong.

Some of the first stars could be cooler and more massive than standard stellar models would suggest, due to the effects of dark matter annihilation in their cores. It has recently been argued that such objects may attain masses in the 10^4--10^7 solar mass range, and that such supermassive dark stars should be within reach of the upcoming James Webb Space Telescope. Notwithstanding theoretical difficulties with this proposal, we argue here that some of these objects should also be readily detectable with both the Hubble Space Telescope and ground-based 8--10 m class telescopes. Existing survey data already place strong constraints on 10^7 solar mass dark stars at z~10. We show that such objects must be exceedingly rare or short-lived to have avoided detection.

We present observations from the Small Array of the Arcminute Microkelvin Imager (AMI) of eight high X-ray luminosity galaxy cluster systems selected from the Local Cluster Substructure Survey (LoCuSS) sample.We detect the Sunyaev-Zel'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 M-T 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'\iguez-Gonz\'alvez et al. (2011) and using large-radius X-ray 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 cm-wave 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_rad-L_IR correlation and suggest that radio flux density may be used as a proxy for the internal luminosity of low luminosity protostars.

Ultracompact minihalos have recently been proposed as a new class of dark matter structure. These minihalos would be produced by phase transitions in the early Universe or features in the inflaton potential, and constitute non-baryonic massive compact halo objects (MACHOs) today. We examine the prospect of detecting ultracompact minihalos in gamma-rays if dark matter consists of self-annihilating particles. We compute present-day fluxes from minihalos produced in the electron-positron annihilation epoch, and the QCD and electroweak phase transitions in the early Universe. Even at a distance of 100 pc, minihalos produced during the electron-positron annihilation epoch should be eminently detectable today, either by the Fermi satellite, current Air Cherenkov telescopes, or even in archival EGRET data. Within ~1 pc, minihalos formed in the QCD phase transition would have similar predicted fluxes to the dwarf spheroidal galaxies targeted by current indirect dark matter searches, so might also be detectable by present or upcoming experiments.

We present observations using the Small Array of the Arcminute Microkelvin Imager (AMI; 14-18 GHz) of four Abell and three MACS clusters spanning 0.171-0.686 in redshift. We detect Sunyaev-Zel'dovich (SZ) signals in five of these without any attempt at source subtraction, although strong source contamination is present. With radio-source measurements from high-resolution 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 more-distant 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.

The first stars in the history of the Universe are likely to form in the dense central regions of 10^5-10^6 Msolar cold dark matter halos at z=10-50. The annihilation of dark matter particles in these environments may lead to the formation of so-called dark stars, which are predicted to be cooler, larger, more massive and potentially more long-lived than conventional population III stars. Here, we investigate the prospects of detecting high-redshift dark stars with the upcoming James Webb Space Telescope (JWST). We find that dark stars at z>6 are intrinsically too faint to be detected by JWST. However, by exploiting foreground galaxy clusters as gravitational telescopes, certain varieties of cool (Teff < 30000 K) dark stars should be within reach at redshifts up to z=10. If the lifetimes of dark stars are sufficiently long, many such objects may also congregate inside the first galaxies. We demonstrate that this could give rise to peculiar features in the integrated spectra of galaxies at high redshifts, provided that dark stars make up at least 1 percent of the total stellar mass in such objects.

The Very Small Array (VSA) has been used to survey the l = 27 to 46 deg, |b|<4 deg region of the Galactic plane at a resolution of 13 arcmin. The survey consists of 44 pointings of the VSA, each with a r.m.s. sensitivity of ~90 mJy/beam. These data are combined in a mosaic to produce a map of the area. The majority of the sources within the map are HII regions. We investigated anomalous radio emission from the warm dust in 9 HII regions of the survey by making spectra extending from GHz frequencies to the FIR IRAS frequencies. Acillary radio data at 1.4, 2.7, 4.85, 8.35, 10.55, 14.35 and 94 GHz in addition to the 100, 60, 25 and 12 micron IRAS bands were used to construct the spectra. From each spectrum the free-free, thermal dust and anomalous dust emission were determined for each HII region. The mean ratio of 33 GHz anomalous flux density to FIR 100 micron flux density for the 9 selected HII regions was 1.10 +/-0.21x10^(-4). When combined with 6 HII regions previously observed with the VSA and the CBI, the anomalous emission from warm dust in HII regions is detected with a 33 GHz emissivity of 4.65 +/- 0.4 micro K/ (MJy/sr) at 11.5{\sigma}. The anomalous radio emission in HII regions is on average 41+/-10 per cent of the radio continuum at 33 GHz.

We report new cm-wave measurements at five frequencies between 15 and 18GHz of the continuum emission from the reportedly anomalous "region 4" of the nearby galaxy NGC6946. We find that the emission in this frequency range is significantly in excess of that measured at 8.5GHz, but has a spectrum from 15-18GHz consistent with optically thin free-free emission from a compact HII region. In combination with previously published data we fit four emission models containing different continuum components using the Bayesian spectrum analysis package radiospec. These fits show that, in combination with data at other frequencies, a model with a spinning dust component is slightly preferred to those that possess better-established emission mechanisms.

The Constrained Minimal Supersymmetric Standard Model (CMSSM) is one of the simplest and most widely-studied supersymmetric extensions to the standard model of particle physics. Nevertheless, current data do not sufficiently constrain the model parameters in a way completely independent of priors, statistical measures and scanning techniques. We present a new technique for scanning supersymmetric parameter spaces, optimised for frequentist profile likelihood analyses and based on Genetic Algorithms. We apply this technique to the CMSSM, taking into account existing collider and cosmological data in our global fit. We compare our method to the MultiNest algorithm, an efficient Bayesian technique, paying particular attention to the best-fit points and implications for particle masses at the LHC and dark matter searches. Our global best-fit point lies in the focus point region. We find many high-likelihood points in both the stau co-annihilation and focus point regions, including a previously neglected section of the co-annihilation region at large m_0. We show that there are many high-likelihood points in the CMSSM parameter space commonly missed by existing scanning techniques, especially at high masses. This has a significant influence on the derived confidence regions for parameters and observables, and can dramatically change the entire statistical inference of such scans.

We present 25 arcsecond resolution radio images of five Lynds Dark Nebulae (L675, L944, L1103, L1111 & L1246) at 16 GHz made with the Arcminute Microkelvin Imager (AMI) Large Array. These objects were previously observed with the AMI Small Array to have an excess of emission at microwave frequencies relative to lower frequency radio data. In L675 we find a flat spectrum compact radio counterpart to the 850 micron emission seen with SCUBA and suggest that it is cm-wave emission from a previously unknown deeply embedded young protostar. In the case of L1246 the cm-wave emission is spatially correlated with 8 micron emission seen with Spitzer. Since the MIR emission is present only in Spitzer band 4 we suggest that it arises from a population of PAH molecules, which also give rise to the cm-wave emission through spinning dust emission.

The dwarf galaxy Segue 1 is one of the most promising targets for the indirect detection of dark matter. Here we examine what constraints 9 months of Fermi-LAT gamma-ray observations of Segue 1 place upon the Constrained Minimal Supersymmetric Standard Model (CMSSM), with the lightest neutralino as the dark matter particle. We use nested sampling to explore the CMSSM parameter space, simultaneously fitting other relevant constraints from accelerator bounds, the relic density, electroweak precision observables, the anomalous magnetic moment of the muon and B-physics. We include spectral and spatial fits to the Fermi observations, a full treatment of the instrumental response and its related uncertainty, and detailed background models. We also perform an extrapolation to 5 years of observations, assuming no signal is observed from Segue 1 in that time. Results marginally disfavour models with low neutralino masses and high annihilation cross-sections. Virtually all of these models are however already disfavoured by existing experimental or relic density constraints.

The dust feature G159.6--18.5 in the Perseus region has previously been observed with the COSMOSOMAS experiment \citep{Watson:05} on angular scales of $\approx$ 1$^{\circ}$, and was found to exhibit anomalous microwave emission. We present new observations of this dust feature, performed with the Very Small Array (VSA) at 33 GHz, to help increase the understanding of the nature of this anomalous emission. On the angular scales observed with the VSA ($\approx$ 10 -- 40$^{\prime}$), G159.6--18.5 consists of five distinct components, each of which have been individually analysed. All five of these components are found to exhibit an excess of emission at 33 GHz, and are found to be highly correlated with far-infrared emission. We provide evidence that each of these compact components have anomalous emission that is consistent with electric dipole emission from very small, rapidly rotating dust grains. These components contribute $\approx$ 10 % to the flux density of the diffuse extended emission detected by COSMOSOMAS, and are found to have a similar radio emissivity.

In regions of very high dark matter density such as the Galactic centre, the capture and annihilation of WIMP dark matter by stars has the potential to significantly alter their evolution. We describe the dark stellar evolution code DarkStars, and present a series of grids of WIMP-influenced stellar models for main sequence stars. We describe changes in which occur as a function of the rate of energy injection by WIMPs, for stars of 0.3-2.0 solar masses and metallicities Z = 0.0003-0.02. We show what rates of energy injection can be obtained using realistic orbital parameters for stars at the Galactic centre, including detailed consideration of the velocity and density profiles of dark matter. Capture and annihilation rates are strongly boosted when stars follow elliptical rather than circular orbits. If there is a spike of dark matter induced by the supermassive black hole at the Galactic centre, single solar-mass stars following orbits with periods as long as 50 years and eccentricities as low as 0.9 could be significantly affected. Binary systems with similar periods about the Galactic centre could be affected on even less eccentric orbits. The most striking evidence of this scenario would be the existence of a binary consisting of a low-mass protostar and a higher-mass evolved star. The observation of low-mass stars and/or binaries on such orbits would either provide a detection of WIMP dark matter, or place stringent limits on the combination of the WIMP mass, spin-dependent nuclear-scattering cross-section, halo density and velocity distribution near the Galactic centre. In some cases, the limits on the WIMP mass and spin-dependent nuclear-scattering cross-section would be of comparable sensitivity to current direct-detection experiments.

The solar chemical composition is an important ingredient in our understanding of the formation, structure and evolution of both the Sun and our solar system. Furthermore, it is an essential reference standard against which the elemental contents of other astronomical objects are compared. In this review we evaluate the current understanding of the solar photospheric composition. In particular, we present a re-determination of the abundances of nearly all available elements, using a realistic new 3-dimensional (3D), time-dependent hydrodynamical model of the solar atmosphere. We have carefully considered the atomic input data and selection of spectral lines, and accounted for departures from LTE whenever possible. The end result is a comprehensive and homogeneous compilation of the solar elemental abundances. Particularly noteworthy findings are significantly lower abundances of carbon, nitrogen, oxygen and neon compared with the widely-used values of a decade ago. The new solar chemical composition is supported by a high degree of internal consistency between available abundance indicators, and by agreement with values obtained in the solar neighborhood and from the most pristine meteorites. There is, however, a stark conflict with standard models of the solar interior according to helioseismology, a discrepancy that has yet to find a satisfactory resolution.

Observations at 14.2 to 17.9 GHz made with the AMI Small Array towards fourteen Lynds Dark Nebulae with a resolution of 2' are reported. These sources are selected from the SCUBA observations of Visser et al. (2001) as small angular diameter clouds well matched to the synthesized beam of the AMI Small Array. Comparison of the AMI observations with radio observations at lower frequencies with matched uv-plane coverage is made, in order to search for any anomalous excess emission which can be attributed to spinning dust. Possible emission from spinning dust is identified as a source within a 2' radius of the Scuba position of the Lynds dark nebula, exhibiting an excess with respect to lower frequency radio emission. We find five sources which show a possible spinning dust component in their spectra. These sources have rising spectral indices in the frequency range 14.2--17.9 GHz. Of these five one has already been reported, L1111, we report one new definite detection, L675, and three new probable detections (L944, L1103 and L1246). The relative certainty of these detections is assessed on the basis of three criteria: the extent of the emission, the coincidence of the emission with the Scuba position and the likelihood of alternative explanations for the excess. Extended microwave emission makes the likelihood of the anomalous emission arising as a consequence of a radio counterpart to a protostar or a proto-planetary disk unlikely. We use a 2' radius in order to be consistent with the IRAS identifications of dark nebulae (Parker 1988), and our third criterion is used in the case of L1103 where a high flux density at 850 microns relative to the FIR data suggests a more complicated emission spectrum.

In our first paper we described three regions of the 9C survey of radio sources with the Ryle telescope at 15.2 GHz, constituting a total area of 520 deg^2 to a completeness limit of approximately 25 mJy. Here we report on a series of deeper regions, amounting to an area of 115 deg^2 complete to approximately 10 mJy and of 29 deg^2 complete to approximately 5.5 mJy. We have investigated the source counts and the distributions of the 1.4 to 15.2 GHz spectral index (alpha) for these deeper samples. The whole catalogue of 643 sources is available online. Down to our lower limit of 5.5 mJy we detect no evidence for any change in the differential source count from the earlier fitted count above 25 mJy. We have matched both our new and earlier catalogues with the NRAO VLA Sky Survey (NVSS) catalogue at 1.4 GHz. For samples of sources selected at 15.2 GHz, in three flux density ranges, we detect a significant shift in the median value of alpha; samples with higher flux densities have higher proportions of sources with flat and rising spectra. We suggest that this observed shift is consistent with a model containing two distinct source populations having differently sloped source counts. Samples selected at 1.4 GHz contain significantly smaller proportions of sources with flat and rising spectra. Also, in our area complete to approximately 10 mJy, we find 5 sources between 10 to 15 mJy, amounting to 4.3 per cent of sources in this range, with no counterpart in the NVSS catalogue. These results illustrate the problems inherent in using a low frequency catalogue to characterise the source population at a much higher frequency and emphasise the value of our blind 15.2 GHz survey.

Using the Arcminute Microkelvin Imager (AMI) at 16 GHz and the Very Small Array (VSA) at 33 GHz to make follow-up observations of sources in the New Extragalactic WMAP Point Source Catalogue, we have investigated the flux-density variability in a complete sample of 97 sources over timescales of a few months to approximately 1.5 years. We find that 53 per cent of the 93 sources, for which we have multiple observations, are variable, at the 99 per cent confidence level, above the flux density calibration uncertainties of approximately 4 per cent at 16 GHz; the fraction of sources having varied by more than 20 per cent is 15 per cent at 16 GHz and 20 per cent at 33 GHz. Not only is this common occurrence of variability at high frequency of interest for source physics, but strategies for coping with source contamination in CMB work must take this variability into account. There is no strong evidence of a correlation between variability and flux density for the sample as a whole. Using a maximum-likelihood method, we calculate the correlation in the variability at the two frequencies in a subset of sources classified as variable from both the AMI and VSA data and find the Pearson product-moment correlation coefficient to be very high (0.955 +/- 0.034). We also find the degree of variability at 16 GHz (0.202 +/- 0.028) to be very similar to that at 33 GHz (0.224 +/- 0.039). Finally, we have investigated the relationship between variability and spectral index, alpha_13.9^33.75 (where S proportional to nu^-alpha), and find a significant difference in the spectral indices of the variable sources (-0.06 +/- 0.05) and non-variable sources (0.13 +/- 0.04).

We present follow-up observations of 97 point sources from the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data, contained within the New Extragalactic WMAP Point Source (NEWPS) catalogue between declinations of -4 and +60 degrees; the sources form a flux-density-limited sample complete to 1.1 Jy (approximately 5 sigma) at 33 GHz. Our observations were made at 16 GHz using the Arcminute Microkelvin Imager (AMI) and at 33 GHz with the Very Small Array (VSA). 94 of the sources have reliable, simultaneous -- typically a few minutes apart -- observations with both telescopes. The spectra between 13.9 and 33.75 GHz are very different from those of bright sources at low frequency: 44 per cent have rising spectra (alpha < 0.0), where flux density is proportional to frequency^-alpha, and 93 per cent have spectra with alpha < 0.5; the median spectral index is 0.04. For the brighter sources, the agreement between VSA and WMAP 33-GHz flux densities averaged over sources is very good. However, for the fainter sources, the VSA tends to measure lower values for the flux densities than WMAP. We suggest that the main cause of this effect is Eddington bias arising from variability.

We announce the public release of the 'dark' stellar evolution code DarkStars. The code simultaneously solves the equations of WIMP capture and annihilation in a star with those of stellar evolution assuming approximate hydrostatic equilibrium. DarkStars includes the most extensive WIMP microphysics of any dark evolution code to date. The code employs detailed treatments of the capture process from a range of WIMP velocity distributions, as well as composite WIMP distribution and conductive energy transport schemes based on the WIMP mean-free path in the star. We give a brief description of the input physics and practical usage of the code, as well as examples of its application to dark stars at the Galactic centre.

We present observations between 1.4 and 18 GHz confirming that G64.5+0.9 is new Galactic shell supernova remnant, using the Very Large Array and the Arcminute Microkelvin Imager. The remnant is a shell ~8 arcmin in diameter with a spectral index of alpha = 0.47 +/- 0.03 (with alpha defined such that flux density S varies with frequency nu as S proportional to nu to the power of -alpha). There is also emission near the centre of the shell, ~1 arcmin in extent, with a spectral index of alpha = 0.81 +/- 0.02. We do not find any evidence for spectral breaks for either source within our frequency range. The nature of the central object is unclear and requires further investigation, but we argue that is most unlikely to be extragalactic. It is difficult to avoid the conclusion that it is associated with the shell, although its spectrum is very unlike that of known pulsar wind nebulae.

We present observations between 14.2 and 17.9 GHz of 12 reported supernova remnants (SNRs) made with the Arcminute Microkelvin Imager Small Array (AMI SA). In conjunction with data from the literature at lower radio frequencies, we determine spectra of these objects. For well-studied SNRs (Cas A, Tycho's SNR, 3C58 and the Crab Nebula), the results are in good agreement with spectra based on previous results. For the less well-studied remnants the AMI SA observations provide higher-frequency radio observations than previously available, and better constrain their radio spectra. The AMI SA results confirm a spectral turnover at ~11 GHz for the filled-centre remnant G74.9+1.2. We also see a possible steepening of the spectrum of the filled-centre remnant G54.1+0.3 within the AMI SA frequency band compared with lower frequencies. We confirm that G84.9+0.5, which had previously been identified as a SNR, is rather an HII region and has a flat radio spectrum.

Determinations of the solar oxygen content relying on the neutral forbidden transition at 630 nm depend upon the nickel abundance, due to a Ni I blend. Here we rederive the solar nickel abundance, using the same ab initio 3D hydrodynamic model of the solar photosphere employed in the recent revision of the abundances of C, N, O and other elements. Using 17 weak, unblended lines of Ni I together with the most accurate atomic and observational data available we find log epsilon_Ni = 6.17 +/- 0.02 (statistical) +/- 0.05 (systematic), a downwards shift of 0.06 to 0.08 dex relative to previous 1D-based abundances. We investigate the implications of the new nickel abundance for studies of the solar oxygen abundance based on the [O I] 630 nm line in the quiet Sun. Furthermore, we demonstrate that the oxygen abundance implied by the recent sunspot spectropolarimetric study of Centeno & Socas-Navarro needs to be revised downwards from log epsilon_O = 8.86 +/- 0.07 to 8.71 +/- 0.10. This revision is based on the new nickel abundance, application of the best available gf-value for the 630 nm forbidden oxygen line, and a more transparent treatment of CO formation. Determinations of the solar oxygen content relying on forbidden lines now appear to converge around log epsilon_O = 8.7.

We present observations of the Lynds' dark nebula LDN 1111 made at microwave frequencies between 14.6 and 17.2 GHz with the Arcminute Microkelvin Imager (AMI). We find emission in this frequency band in excess of a thermal free--free spectrum extrapolated from data at 1.4 GHz with matched uv-coverage. This excess is > 15 sigma above the predicted emission. We fit the measured spectrum using the spinning dust model of Drain & Lazarian (1998a) and find the best fitting model parameters agree well with those derived from Scuba data for this object by Visser et al. (2001).

The presence of large amounts of WIMP dark matter in stellar cores has been shown to have significant effects upon models of stellar evolution. We present a series of detailed grids of WIMP-influenced stellar models for main sequence stars, computed using the DarkStars code. We describe the changes in stellar structure and main sequence evolution which occur for masses ranging from 0.3 to 2.0 solar masses and metallicities from Z = 0.0003-0.02, as a function of the rate of energy injection by WIMPs. We then go on to show what rates of energy injection can be obtained using realistic orbital parameters for stars near supermassive black holes, including detailed considerations of dark matter halo velocity and density profiles. Capture and annihilation rates are strongly boosted when stars follow elliptical rather than circular orbits, causing WIMP annihilation to provide up to 100 times the energy of hydrogen fusion in stars at the Galactic centre.

We present the current status of the QUIJOTE (Q-U-I JOint TEnerife) CMB Experiment, a new instrument which will start operations early 2009 at Teide Observatory, with the aim of characterizing the polarization of the CMB and other processes of galactic and extragalactic emission in the frequency range 10-30 GHz and at large angular scales. QUIJOTE will be a valuable complement at low frequencies for the PLANCK mission, and will have the required sensitivity to detect a primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r=0.05.

The Arcminute Microkelvin Imager is a pair of interferometer arrays operating with six frequency channels spanning 13.9-18.2 GHz, with very high sensitivity to angular scales 30''-10'. The telescope is aimed principally at Sunyaev-Zel'dovich imaging of clusters of galaxies. We discuss the design of the telescope and describe and explain its electronic and mechanical systems.

Accurate calibration of data is essential for the current generation of CMB experiments. Using data from the Very Small Array (VSA), we describe procedures which will lead to an accuracy of 1 percent or better for experiments such as the VSA and CBI. Particular attention is paid to the stability of the receiver systems, the quality of the site and frequent observations of reference sources. At 30 GHz the careful correction for atmospheric emission and absorption is shown to be essential for achieving 1 percent precision. The sources for which a 1 percent relative flux density calibration was achieved included Cas A, Cyg A, Tau A and NGC7027 and the planets Venus, Jupiter and Saturn. A flux density, or brightness temperature in the case of the planets, was derived at 33 GHz relative to Jupiter which was adopted as the fundamental calibrator. A spectral index at ~30 GHz is given for each. Cas A,Tau A, NGC7027 and Venus were examined for variability. Cas A was found to be decreasing at $0.394 \pm 0.019$ percent per year over the period March 2001 to August 2004. In the same period Tau A was decreasing at $0.22\pm 0.07$ percent per year. A survey of the published data showed that the planetary nebula NGC7027 decreased at $0.16\pm 0.04$ percent per year over the period 1967 to 2003. Venus showed an insignificant ($1.5 \pm 1.3$ percent) variation with Venusian illumination. The integrated polarization of Tau A at 33 GHz was found to be $7.8\pm 0.6$ percent at pa $ = 148^\circ \pm 3^\circ$.}

We present interferometric imaging at 33 GHz, with the new superextended configuration of the Very Small Array (VSA), of a very deep decrement in the cosmic microwave background (CMB) temperature. This decrement is located in the direction of the Corona Borealis supercluster, at a position with no known galaxy clusters, and was discovered by a previous VSA survey (Genova-Santos et al.). A total area of 3 sq.deg. has now been imaged, with an angular resolution of 7 arcmin and a flux sensitivity of 5 mJy/beam. These observations confirm the presence of this strong and resolved negative spot at -41+/-5 mJy/beam (-258+/-29 muK), with a signal to noise level of 8. This structure is also present in the WMAP 5-year data. The temperature of the W-band (94 GHz) data at the position of the decrement agrees within 1.2-sigma with that observed by the VSA at 33 GHz, and within 0.2-sigma with the Sunyaev-Zel'dovich (SZ) spectrum. Our analyses show that it is a non-Gaussian feature in the CMB at a level of 4.8-sigma. The probability of finding such a deviation or larger in CMB Gaussian simulations is only 0.19 per cent. Therefore, an explanation other than primordial CMB is required. We have considered the possibility of an SZ effect generated in a diffuse, extended warm/hot gas distribution. This hypothesis is especially relevant, as the presence of such structures, if confirmed, could provide the location for a significant fraction of the missing baryons in the Local Universe.

We modify a stellar structure code to estimate the effect upon the main sequence of the accretion of weakly interacting dark matter onto stars and its subsequent annihilation. The effect upon the stars depends upon whether the energy generation rate from dark matter annihilation is large enough to shut off the nuclear burning in the star. Main sequence WIMP burners look much like protostars moving on the Hayashi track, although they are in principle completely stable. We make some brief comments about where such stars could be found, how they might be observed and more detailed simulations which are currently in progress. Finally we comment on whether or not it is possible to link the paradoxically young OB stars found at the galactic centre with WIMP burners.

We present an analysis of the spectral properties of the extragalactic radio sources in the nearly-complete VSA sample at 33GHz. Data from different surveys are used to study source spectra between 1.4 and 33GHz. We find that, in general, spectra can not be well described by a single power law in the range of frequencies considered. In particular, most of the VSA sources that are steep between 1.4 and 5GHz, show a spectral flattening at \nu>5GHz. We identify 20 objects (19% of the sample) clearly characterized by an upturn spectrum, i.e., a spectrum falling at low frequencies (\nu<5GHz) and inverted at higher frequencies. Spectra with high-frequency flattening or upturn shape are supposed to occur when the emission from the AGN compact core begins to dominate over the component from extended lobes. This picture fits well with the AGN unified scheme, for objects observed at intermediate viewing angles of the AGN jet. Finally, we discuss implications that this class of sources can have on future CMB observations at high resolution.

We present observations between 14.2 and 17.9 GHz of sixteen Galactic HII regions made with the Arcminute Microkelvin Imager (AMI). In conjunction with data from the literature at lower radio frequencies we investigate the possibility of a spinning dust component in the spectra of these objects. We conclude that there is no significant evidence for spinning dust towards these sources and measure an average spectral index of 0.15+/-0.07 between 1.4 and 17.9 GHz for the sample.

Recent work has indicated that WIMP annihilation in stellar cores has the potential to contribute significantly to a star's total energy production. We report on progress in simulating the effects of WIMP capture and annihilation upon stellar structure and evolution near supermassive black holes, using the new DarkStars code. Preliminary results indicate that low-mass stars are the most influenced by WIMP annihilation, which could have consequences for upcoming observational programs.

We present results from observations made at 33 GHz with the Very Small Array (VSA) telescope towards potential candidates in the Galactic plane for spinning dust emission. In the cases of the diffuse HII regions LPH96 and NRAO591 we find no evidence for anomalous emission and, in combination with Effelsberg data at 1.4 and 2.7 GHz, confirm that their spectra are consistent with optically thin free--free emission. In the case of the infra-red bright SNR 3C396 we find emission inconsistent with a purely non-thermal spectrum and discuss the possibility of this excess arising from either a spinning dust component or a shallow spectrum PWN, although we conclude that the second case is unlikely given the strong constraints available from lower frequency radio images.

Aims: We describe a 6-12 GHz analogue correlator that has been developed for use in radio interferometers. Methods: We use a lag-correlator technique to synthesis eight complex spectral channels. Two schemes were considered for sampling the cross-correlation function, using either real or complex correlations, and we developed prototypes for both of them. We opted for the ``add and square'' detection scheme using Schottky diodes over the more commonly used active multipliers because the stability of the device is less critical. Results: We encountered an unexpected problem, in that there were errors in the lag spacings of up to ten percent of the unit spacing. To overcome this, we developed a calibration method using astronomical sources which corrects the effects of the non-uniform sampling as well as gain error and dispersion in the correlator.

We have observed the Corona Borealis Supercluster with the Millimeter and Infrared Testa grigia Observatory (MITO), located in the Italian Alps, at 143, 214, 272, and 353 GHz. We present a description of the measurements, data analysis, and results of the observations together with a comparison with observations performed at 33 GHz with the Very Small Array (VSA) interferometer situated at the Teide Observatory (Tenerife-Spain). Observations have been made in the direction of the supercluster towards one Cosmic Microwave Background (CMB) cold spot previously detected in a VSA temperature map. Observational strategy and data analysis are described in detail, explaining the procedures used to disentangle primary and secondary anisotropies in the resulting maps. From a first level of data analysis we find evidence in MITO data of primary anisotropy, however still with room for the presence of secondary anisotropy, especially when VSA results are included. With a second level of data analysis using map-making and the maximum entropy method we claim a weak detection of a faint signal compatible with a SZ effect characterized at most by a Comptonization parameter $y=(7.8^{+5.3}_{-4.4})\times10^{-6}$ 68% CL. The low level of confidence in the presence of a SZ signal invite us to study this sky region with higher sensitivity and angular resolution experiments like the already planned upgraded versions of VSA and MITO.

CO spectral line formation in the Sun has long been a source of consternation for solar physicists, as have the elemental abundances it seems to imply. We modelled solar CO line formation using a realistic, ab initio, time-dependent 3D radiative-hydrodynamic model atmosphere. Results were compared with observations from the space-based ATMOS experiment. We employed weak 12C16O, 13C16O and 12C18O lines from the fundamental and first overtone bands to determine the solar carbon abundance, as well as the 12C/13C and 16O/18O isotopic ratios. A weighted carbon abundance of log epsilonC = 8.39 +-0.05 was found. We note with satisfaction that the derived abundance is identical to our recent 3D determination based on CI, [C I], C2 and CH lines. Identical calculations were carried out using 1D models, but only the 3D model was able to produce abundance agreement between different CO lines and the other atomic and molecular diagnostics. Solar 12C/13C and 16O/18O ratios were measured as 86.8+3.9-3.7 (delta13C = 30+46-44) and 479+29-28 (delta18O = 41+67-59), respectively. These values may require current theories of solar system formation to be revised. Excellent agreement was seen between observed and predicted weak CO line shapes, without invoking micro- or macroturbulence. Agreement breaks down for the strongest CO lines however, which are formed in very high atmospheric layers. The simplest explanation is that temperatures are overestimated in the highest layers of the 3D simulation. Thus, our analysis supports the presence of a COmosphere above the traditional photospheric temperature minimum, with an average temperature of less than 4000K. The shortcoming of the model atmosphere is not surprising, given that it was never intended to properly describe such high layers.

(Abridged) We have used the Rayner & Best (1989) smooth tests of goodness-of-fit to study the Gaussianity of the Very Small Array (VSA) data. Out of the 41 published VSA individual pointings dedicated to cosmological observations, 37 are found to be consistent with Gaussianity, whereas four pointings show deviations from Gaussianity. In two of them, these deviations can be explained as residual systematic effects of a few visibility points which, when corrected, have a negligible impact on the angular power spectrum. The non-Gaussianity found in the other two (adjacent) pointings seems to be associated to a local deviation of the power spectrum of these fields with respect to the common power spectrum of the complete data set, at angular scales of the third acoustic peak (l = 700-900). No evidence of residual systematics is found in this case, and unsubstracted point sources are not a plausible explanation either. If those visibilities are removed, a cosmological analysis based on this new VSA power spectrum alone shows no differences in the parameter constraints with respect to our published results, except for the physical baryon density, which decreases by 10 percent. Finally, the method has been also used to analyse the VSA observations in the Corona Borealis supercluster region (Genova-Santos et al. 2005), which show a strong decrement which cannot be explained as primordial CMB. Our method finds a clear deviation (99.82%) with respect to Gaussianity in the second-order moment of the distribution, and which can not be explained as systematic effects. A detailed study shows that the non-Gaussianity is produced in scales of l~500, and that this deviation is intrinsic to the data (in the sense that can not be explained in terms of a Gaussian field with a different power spectrum).

We present coincident observations of the Cosmic Microwave Background (CMB) from the Very Small Array (VSA) and Cosmic Background Imager (CBI) telescopes. The consistency of the full datasets is tested in the map plane and the Fourier plane, prior to the usual compression of CMB data into flat bandpowers. Of the three mosaics observed by each group, two are found to be in excellent agreement. In the third mosaic, there is a 2 sigma discrepancy between the correlation of the data and the level expected from Monte Carlo simulations. This is shown to be consistent with increased phase calibration errors on VSA data during summer observations. We also consider the parameter estimation method of each group. The key difference is the use of the variance window function in place of the bandpower window function, an approximation used by the VSA group. A re-evaluation of the VSA parameter estimates, using bandpower windows, shows that the two methods yield consistent results.

We report the first detection of a Sunyaev-Zel'dovich (S-Z) decrement with the Arcminute Microkelvin Imager (AMI). We have made commissioning observations towards the cluster A1914 and have measured an integrated flux density of -8.61 mJy in a uv-tapered map with noise level 0.19 mJy/beam. We find that the spectrum of the decrement, measured in the six channels between 13.5-18GHz, is consistent with that expected for a S-Z effect. The sensitivity of the telescope is consistent with the figures used in our simulations of cluster surveys with AMI.

We present interferometric imaging at 33 GHz of the Corona Borealis supercluster, using the extended configuration of the Very Small Array. A total area of 24 deg^2 has been imaged, with an angular resolution of 11 arcmin and a sensitivity of 12 mJy/beam. The aim of these observations is to search for Sunyaev-Zel'dovich (SZ) detections from known clusters of galaxies in this supercluster and for a possible extended SZ decrement due to diffuse warm/hot gas in the intercluster medium. We measure negative flux values in the positions of the ten richest clusters in the region. Collectively, this implies a 3.0-sigma detection of the SZ effect. In the clusters A2061 and A2065 we find decrements of approximately 2-sigma. Our main result is the detection of two strong and resolved negative features at -70+-12 mJy/beam (-157+-27 microK) and -103+-10 mJy/beam (-230+-23 microK), respectively, located in a region with no known clusters, near the centre of the supercluster. We discuss their possible origins in terms of primordial CMB anisotropies and/or SZ signals related to either unknown clusters or to a diffuse extended warm/hot gas distribution. Our analyses have revealed that a primordial CMB fluctuation is a plausible explanation for the weaker feature (probability of 37.82%). For the stronger one, neither primordial CMB (probability of 0.33%) nor SZ can account alone for its size and total intensity. The most reasonable explanation, then, is a combination of both primordial CMB and SZ signal. Finally, we explore what characteristics would be required for a filamentary structure consisting of warm/hot diffuse gas in order to produce a significant contribution to such a spot taking into account the constraints set by X-ray data.

We describe the source subtraction strategy and observations for the extended Very Small Array, a CMB interferometer operating at 33 GHz. A total of 453 sources were monitored at 33 GHz using a dedicated source subtraction baseline. 131 sources brighter than 20 mJy were directly subtracted from the VSA visibility data. Some characteristics of the subtracted sources, such as spectra and variability, are discussed. The 33-GHz source counts are estimated from a sample selected at 15 GHz. The selection of VSA fields in order to avoid bright sources introduces a bias into the observed counts. This bias is corrected and the resulting source count is estimated to be complete in the flux-density range 20-114 mJy. The 33-GHz source counts are used to calculate a correction to the VSA power spectrum for sources below the subtraction limit.

We present deep Ka-band ($\nu \approx 33$ GHz) observations of the CMB made with the extended Very Small Array (VSA). This configuration produces a naturally weighted synthesized FWHM beamwidth of $\sim 11$ arcmin which covers an $\ell$-range of 300 to 1500. On these scales, foreground extragalactic sources can be a major source of contamination to the CMB anisotropy. This problem has been alleviated by identifying sources at 15 GHz with the Ryle Telescope and then monitoring these sources at 33 GHz using a single baseline interferometer co-located with the VSA. Sources with flux densities $\gtsim 20$ mJy at 33 GHz are subtracted from the data. In addition, we calculate a statistical correction for the small residual contribution from weaker sources that are below the detection limit of the survey. The CMB power spectrum corrected for Galactic foregrounds and extragalactic point sources is presented. A total $\ell$-range of 150-1500 is achieved by combining the complete extended array data with earlier VSA data in a compact configuration. Our resolution of $\Delta \ell \approx 60$ allows the first 3 acoustic peaks to be clearly delineated. The is achieved by using mosaiced observations in 7 regions covering a total area of 82 sq. degrees. There is good agreement with WMAP data up to $\ell=700$ where WMAP data run out of resolution. For higher $\ell$-values out to $\ell = 1500$, the agreement in power spectrum amplitudes with other experiments is also very good despite differences in frequency and observing technique.

We estimate the bispectrum of the Very Small Array data from the compact and extended configuration observations released in December 2002, and compare our results to those obtained from Gaussian simulations. There is a slight excess of large bispectrum values for two individual fields, but this does not appear when the fields are combined. Given our expected level of residual point sources, we do not expect these to be the source of the discrepancy. Using the compact configuration data, we put an upper limit of 5400 on the value of f_NL, the non-linear coupling parameter, at 95 per cent confidence. We test our bispectrum estimator using non-Gaussian simulations with a known bispectrum, and recover the input values.

We estimate cosmological parameters using data obtained by the Very Small Array (VSA) in its extended configuration, in conjunction with a variety of other CMB data and external priors. Within the flat $\Lambda$CDM model, we find that the inclusion of high resolution data from the VSA modifies the limits on the cosmological parameters as compared to those suggested by WMAP alone, while still remaining compatible with their estimates. We find that $\Omega_{\rm b}h^2=0.0234^{+0.0012}_{-0.0014}$, $\Omega_{\rm dm}h^2=0.111^{+0.014}_{-0.016}$, $h=0.73^{+0.09}_{-0.05}$, $n_{\rm S}=0.97^{+0.06}_{-0.03}$, $10^{10}A_{\rm S}=23^{+7}_{-3}$ and $\tau=0.14^{+0.14}_{-0.07}$ for WMAP and VSA when no external prior is included.On extending the model to include a running spectral index of density fluctuations, we find that the inclusion of VSA data leads to a negative running at a level of more than 95% confidence ($n_{\rm run}=-0.069\pm 0.032$), something which is not significantly changed by the inclusion of a stringent prior on the Hubble constant. Inclusion of prior information from the 2dF galaxy redshift survey reduces the significance of the result by constraining the value of $\Omega_{\rm m}$. We discuss the veracity of this result in the context of various systematic effects and also a broken spectral index model. We also constrain the fraction of neutrinos and find that $f_{\nu}< 0.087$ at 95% confidence which corresponds to $m_\nu<0.32{\rm eV}$ when all neutrino masses are the equal. Finally, we consider the global best fit within a general cosmological model with 12 parameters and find consistency with other analyses available in the literature. The evidence for $n_{\rm run}<0$ is only marginal within this model.

We have tested Very Small Array (VSA) observations of three regions of sky for the presence of non-Gaussianity, using high-order cumulants, Minkowski functionals, a wavelet-based test and a Bayesian joint power spectrum/non-Gaussianity analysis. We find the data from two regions to be consistent with Gaussianity. In the third region, we obtain a 96.7% detection of non-Gaussianity using the wavelet test. We perform simulations to characterise the tests, and conclude that this is consistent with expected residual point source contamination. There is therefore no evidence that this detection is of cosmological origin. Our simulations show that the tests would be sensitive to any residual point sources above the data's source subtraction level of 20 mJy. The tests are also sensitive to cosmic string networks at an rms fluctuation level of $105 \mu K$ (i.e. equivalent to the best-fit observed value). They are not sensitive to string-induced fluctuations if an equal rms of Gaussian CDM fluctuations is added, thereby reducing the fluctuations due to the strings network to $74 \mu K$ rms . We especially highlight the usefulness of non-Gaussianity testing in eliminating systematic effects from our data.

We investigate the constraints on basic cosmological parameters set by the first compact-configuration observations of the Very Small Array (VSA), and other cosmological data sets, in the standard inflationary LambdaCDM model. Using a weak prior 40 < H_0 < 90 km/s/Mpc and 0 < tau < 0.5 we find that the VSA and COBE_DMR data alone produce the constraints Omega_tot = 1.03^{+0.12}_{-0.12}, Omega_bh^2 = 0.029^{+0.009}_{-0.009}, Omega_cdm h^2 = 0.13^{+0.08}_{-0.05} and n_s = 1.04^{+0.11}_{-0.08} at the 68 per cent confidence level. Adding in the type Ia supernovae constraints, we additionally find Omega_m = 0.32^{+0.09}_{-0.06} and Omega_Lambda = 0.71^{+0.07}_{-0.07}. These constraints are consistent with those found by the BOOMERanG, DASI and MAXIMA experiments. We also find that, by combining all the recent CMB experiments and assuming the HST key project limits for H_0 (for which the X-ray plus Sunyaev--Zel'dovich route gives a similar result), we obtain the tight constraints Omega_m=0.28^{+0.14}_{-0.07} and Omega_Lambda= 0.72^{+0.07}_{-0.13}, which are consistent with, but independent of, those obtained using the supernovae data.

The fields chosen for the first observations of the cosmic microwave background with the Very Small Array have been surveyed with the Ryle Telescope at 15 GHz. We have covered three regions around RA 00h20m Dec +30deg, RA 09h40m Dec +32deg and RA 15h40m Dec +43deg (J2000), an area of 520 deg^2. There are 465 sources above the current completeness limit of approximately 25 mJy, although a total of ~760 sources have been detected, some as faint as 10 mJy. This paper describes our techniques for observation and data analysis; it also includes source counts and some discussion of spectra and variability. Preliminary source lists are presented.

We have observed the cosmic microwave background (CMB) in three regions of sky using the Very Small Array (VSA) in an extended configuration with antennas of beamwidth 2 degrees at 34 GHz. Combined with data from previous VSA observations using a more compact array with larger beamwidth, we measure the power spectrum of the primordial CMB anisotropies between angular multipoles l = 160 - 1400. Such measurements at high l are vital for breaking degeneracies in parameter estimation from the CMB power spectrum and other cosmological data. The power spectrum clearly resolves the first three acoustic peaks, shows the expected fall off in power at high l and starts to constrain the position and height of a fourth peak.

The Very Small Array (VSA) is a synthesis telescope designed to image faint structures in the cosmic microwave background on degree and sub-degree angular scales. The VSA has key differences from other CMB interferometers with the result that different systematic errors are expected. We have tested the operation of the VSA with a variety of blank-field and calibrator observations and cross-checked its calibration scale against independent measurements. We find that systematic effects can be suppressed below the thermal noise level in long observations; the overall calibration accuracy of the flux density scale is 3.5 percent and is limited by the external absolute calibration scale.

We present the power spectrum of the fluctuations in the cosmic microwave background detected by the Very Small Array (VSA) in its first season of observations in its compact configuration. We find clear detections of first and second acoustic peaks at l~200 and l~550, plus detection of power on scales up to l=800. The VSA power spectrum is in very good agreement with the results of the Boomerang, Dasi and Maxima telescopes despite the differing potential systematic errors.

We constrain the basic comological parameters using the first observations by the Very Small Array (VSA) in its extended configuration, together with existing cosmic microwave background data and other cosmological observations. We estimate cosmological parameters for four different models of increasing complexity. In each case, careful consideration is given to implied priors and the Bayesian evidence is calculated in order to perform model selection. We find that the data are most convincingly explained by a simple flat Lambda-CDM cosmology without tensor modes. In this case, combining just the VSA and COBE data sets yields the 68 per cent confidence intervals Omega_b h^2=0.034 (+0.007, -0.007), Omega_dm h^2 = 0.18 (+0.06, -0.04), h=0.72 (+0.15,-0.13), n_s=1.07 (+0.06,-0.06) and sigma_8=1.17 (+0.25, -0.20). The most general model considered includes spatial curvature, tensor modes, massive neutrinos and a parameterised equation of state for the dark energy. In this case, by combining all recent cosmological data, we find, in particular, 95 percent limit on the tensor-to-scalar ratio R < 0.63 and on the fraction of massive neutrinos f_nu < 0.11; we also obtain the 68 per cent confidence interval w=-1.06 (+0.20, -0.25) on the equation of state of dark energy.

We have observed the cosmic microwave background temperature fluctuations in eight fields covering three separated areas of sky with the Very Small Array at 34 GHz. A total area of 101 square degrees has been imaged, with sensitivity on angular scales 3.6 - 0.4 degrees (equivalent to angular multipoles l=150-900). We describe the field selection and observing strategy for these observations. In the full-resolution images (with synthesised beam of FWHM ~ 17 arcmin) the thermal noise is typically 45 microK and the CMB signal typically 55 microK. The noise levels in each field agree well with the expected thermal noise level of the telescope, and there is no evidence of any residual systematic features. The same CMB features are detected in separate, overlapping observations. Discrete radio sources have been detected using a separate 15 GHz survey and their effects removed using pointed follow-up observations at 34 GHz. We estimate that the residual confusion noise due to unsubtracted radio sources is less than 14 mJy/beam (15 microK in the full-resolution images), which added in quadrature to the thermal noise increases the noise level by 6 %. We estimate that the rms contribution to the images from diffuse Galactic emission is less than 6 microK. We also present images which are convolved to maximise the signal-to-noise of the CMB features and are co-added in overlapping areas, in which the signal-to-noise of some individual CMB features exceeds 8.

We describe observations at frequencies near 15 GHz of the second 2x2 degree field imaged with the Cambridge Cosmic Anisotropy Telescope (CAT). After the removal of discrete radio sources, structure is detected in the images on characteristic scales of about half a degree, corresponding to spherical harmonic multipoles in the approximate range l= 330--680. A Bayesian analysis confirms that the signal arises predominantly from the cosmic microwave background (CMB) radiation for multipoles in the lower half of this range; the average broad-band power in a bin with centroid l=422 (theta = 51') is estimated to be Delta_T/T=2.1^{+0.4}_{-0.5} x 10^{-5}. For multipoles centred on l=615 (theta =35'), we find contamination from Galactic emission is significant, and constrain the CMB contribution to the measured power in this bin to be Delta_T/T <2.0 x 10^{-5} (1-sigma upper limit). These new results are consistent with the first detection made by CAT in a completely different area of sky. Together with data from other experiments, this new CAT detection adds weight to earlier evidence from CAT for a downturn in the CMB power spectrum on scales smaller than 1 degree. Improved limits on the values of H_0 and Omega are determined using the new CAT data.

We are contructing an interferometric telescope, the Very Small Array, to study the cosmic microwave background on angular scales 0.2--4.5 degrees. The physical layout and electronic design of the telescope are optimised to give maximum protection from systematic effects, while still providing sufficient sensitivity to make high signal-to-noise images. A prototype single baseline is currently being tested, with scientific results expected during 2000.