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

Wavelet functions allow the sparse and efficient representation of a signal at different scales. Recently the application of wavelets to the denoising of maps of cosmic microwave background (CMB) fluctuations has been proposed. The maximum-entropy method (MEM) is also often used for enhancing astronomical images and has been applied to CMB data. In this paper, we give a systematic discussion of combining these two approaches by the use of the MEM in wavelet bases for the denoising and deconvolution of CMB maps and more general images. Certain types of wavelet transforms, such as the a trous transform, can be viewed as a multi-channel intrinsic correlation function (ICF). We find that the wavelet MEM has lower reconstruction residuals than conventional pixel-basis MEM in the case when the signal-to-noise ratio is low and the point spread function narrow. Furthermore, the Bayesian evidence for the wavelet MEM reconstructions is generally higher for a wide range of images. From a Bayesian point of view, the wavelet basis thus provides a better model of the image.

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.

Interferometers require accurate determination of the array configuration in order to produce reliable observations. A method is presented for finding the maximum-likelihood estimate of the telescope geometry, and of other instrumental parameters, astrometrically from the visibility timelines obtained from observations of celestial calibrator sources. The method copes systematically with complicated and unconventional antenna and array geometries, with electronic bandpasses that are different for each antenna radiometer, and with low signal-to-noise ratios for the calibrators. The technique automatically focusses on the geometry errors that are most significant for astronomical observation. We apply this method to observations made with the Very Small Array and constrain some 450 telescope parameters, such as the antenna positions, effective observing frequencies and correlator amplitudes and phase shifts; this requires only ~ 1 h of CPU time on a typical workstation.

A maximum-likelihood method is presented for estimating the power spectrum of anisotropies in the cosmic microwave background (CMB) from interferometer observations. The method calculates flat band-power estimates in separate bins in l-space, together with confidence intervals on the power in each bin. For multifrequency data, the power spectrum of the other foreground components may also be recovered. Advantage is taken of several characteristic properties of interferometer data, which together allow the fast calculation of the likelihood function. The method may be applied to single-field or mosaiced observations, and proper account can be taken of non-coplanar baselines. The method is illustrated by application to simulated data from the Very Small Array.