(Abridged) We present the analysis of the first 18 months of data obtained with the COSMOSOMAS experiment at the Teide Observatory (Tenerife). Three maps have been obtained at 12.7, 14.7 and 16.3 GHz covering 9000 square degrees each with a resolution of ~1 degree and with sensitivities 49, 59 and 115 muK per beam respectively. These data in conjuction with the WMAP first year maps have revealed that the Cosmic Microwave Background (CMB) is the dominant astronomical signal at high galatic latitude in the three COSMOSOMAS channels with an average amplitude of 29.7+/- 1.0 \muK (68% c.l. not including calibration errors). This value is in agreement with the predicted CMB signal in the COSMOSOMAS maps using the best fit Lambda-CDM model to the WMAP power spectrum. Cross-correlation of COSMOSOMAS data with the DIRBE map at 100 \mu m shows the existence of a common signal with amplitude 7.4+/- 1.1, 7.5+/- 1.1, and 6.5+/-2.3 muK in the 12.7, 14.7 and 16.3 GHz COSMOSOMAS maps at |b|>30^\deg. Using the WMAP data we find this DIRBE correlated signal rises from high to low frequencies flattening below ~20 GHz. At higher galactic latitudes the average amplitude of the correlated signal with the DIRBE maps decreases slightly. The frequency behaviour of the COSMOSOMAS/WMAP correlated signal with DIRBE is not compatible with the expected tendency for thermal dust. A study of the H-alpha emission maps do not support free-free as a major contributor to that signal. Our results provide evidence of a new galactic foreground with properties compatible with those predicted by the spinning dust models.

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.

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.

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.

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