The vector-tensor (VT) theory of gravitation revisited in this article was studied in previous papers, where it was proved that VT works and deserves attention. New observational data and numerical codes have motivated further development which is presented here. New research has been planed with the essential aim of proving that current cosmological observations, including Planck data, baryon acoustic oscillations (BAO), and so on, may be explained with VT, a theory which accounts for a kind of dark energy which has the same equation of state as vacuum. New versions of the codes CAMB and COSMOMC have been designed for applications to VT, and the resulting versions have been used to get the cosmological parameters of the VT model at suitable confidence levels. The parameters to be estimated are the same as ingeneral relativity (GR), plus a new parameter $D$. For $D=0$, VT linear cosmological perturbations reduces to those of GR, but the VT background may explain dark energy. The fits between observations and VT predictions lead to non vanishing $|D|$ upper limits at the $1\sigma $ confidence level. The value $D=0$ is admissible at this level, but this value is not that of the best fit in any case. Results strongly suggest that VT may explain current observations, at least, as well as GR; with the advantage that, as it is proved in this paper, VT has an additional parameter which facilitates adjustments to current observational data.}

We study relevant cosmological topics in the framework of a certain vector-tensor theory of gravitation (hereafter VT). This theory is first compared with the so-called extended electromagnetism (EE). These theories have a notable resemblance and both explain the existence of a cosmological constant. It is shown that, in EE, a positive dark energy density requires a Lagrangian leading to quantum ghosts, whereas VT is free from these ghosts. On account of this fact, the remainder of the paper is devoted to study cosmology in the framework of VT. Initial conditions, at high redshift, are used to solve the evolution equations of all the VT scalar modes. In particular, a certain scalar mode characteristic of VT -which does not appear in general relativity (GR)- is chosen in such a way that it evolves separately. In other words, the scalar modes of the standard model based on GR do not affect the evolution of the VT characteristic mode; however, this scalar mode influences the evolution of the standard GR ones. Some well known suitable codes (CMBFAST and COSMOMC) have been modified to include our VT initial conditions and evolution equations, which are fully general. One of the resulting codes -based on standard statistical methods- has been used to fit VT predictions and observational evidences about both Ia supernovae and cosmic microwave background anisotropy. Seven free parameters are used in this fit. Six of them are often used in GR cosmology and the seventh one is characteristic of VT. From the statistical analysis it follows that VT seems to be advantageous against GR in order to explain cosmological observational evidences.

A certain vector-tensor (VT) theory is revisited. It was proposed and analyzed as a theory of electromagnetism without the standard gauge invariance. Our attention is first focused on a detailed variational formulation of the theory, which leads to both a modified Lorentz force and the true energy momentum tensor of the vector field. The theory is then applied to cosmology. A complete gauge invariant treatment of the scalar perturbations is presented. For appropriate gauge invariant variables describing the scalar modes of the vector field (A-modes), it is proved that the evolution equations of these modes do not involve the scalar modes appearing in General Relativity (GR-modes), which are associated to the metric and the energy momentum tensor of the cosmological fluids. However, the A-modes modify the standard gauge invariant equations describing the GR-modes. By using the new formalism, the evolution equations of the A-perturbations are derived and separately solved and, then, the correction terms --due to the A-perturbations-- appearing in the evolution equations of the GR-modes are estimated. The evolution of these correction terms is studied for an appropriate scale. The relevance of these terms depends on both the spectra and the values of the normalization constants involved in extended electromagnetism. Further applications of the new formalism will be presented elsewhere.

We estimate the impact of weak lensing by strongly nonlinear cosmological structures on the cosmic microwave background. Accurate calculation of large $\ell$ multipoles requires N-body simulations and ray-tracing schemes with both high spatial and temporal resolution. To this end we have developed a new code that combines a gravitational Adaptive Particle-Particle, Particle-Mesh (AP3M) solver with a weak lensing evaluation routine. The lensing deviations are evaluated while structure evolves during the simulation so that all evolution steps--rather than just a few outputs--are used in the lensing computations. The new code also includes a ray-tracing procedure that avoids periodicity effects in a universe that is modeled as a 3-D torus in the standard way. Results from our new simulations are compared with previous ones based on Particle-Mesh simulations. We also systematically investigate the impact of box volume, resolution, and ray-tracing directions on the variance of the computed power spectra. We find that a box size of $512 h^{-1}$ Mpc is sufficient to provide a robust estimate of the weak lensing angular power spectrum in the $\ell$-interval (2,000--7,000). For a reaslistic cosmological model the power $[\ell(\ell+1)C_{\ell}/2\pi]^{1/2}$ takes on values of a few $\mu K$ in this interval, which suggests that a future detection is feasible and may explain the excess power at high $\ell$ in the BIMA and CBI observations.

We discuss whether some perturbed Friedmann-Robertson-Walker (FRW) universes could be creatable, i. e., could have vanishing energy, linear momentum and angular momentum, as it could be expectable if the Universe arose as a quantum fluctuation. On account of previous results, the background is assumed to be either closed (with very small curvature) or flat. In the first case, fully arbitrary linear perturbations are considered; whereas in the flat case, it is assumed the existence of: (i) inflationary scalar perturbations, that is to say, Gaussian adiabatic scalar perturbations having an spectrum close to the Harrison-Zel'dovich one, and (ii) arbitrary tensor perturbations. We conclude that, any closed perturbed universe is creatable, and also that, irrespective of the spectrum and properties of the inflationary gravitational waves, perturbed flat FRW universes with standard inflation are not creatable. Some considerations on pre-inflationary scalar perturbations are also presented. The creatable character of perturbed FRW universes is studied, for the first time, in this paper.

Recent observations have pointed out various anomalies in some multipoles (small $\ell $) of the cosmic microwave background (CMB). In this paper, it is proved that some of these anomalies could be explained in the framework of a modified concordance model, in which, there is an appropriate distribution of vector perturbations with very large spatial scales. Vector modes are associated with divergenceless (vortical) velocity fields. Here, the generation of these modes is not studied in detail (it can be done "a posteriori"); on the contrary, we directly look for the distributions of these vector modes which lead to both alignments of the second and third multipoles and a planar octopole. A general three-dimensional (3D) superimposition of vector perturbations does not produce any alignment, but we have found rather general 2D superimpositions leading to anomalies similar to the observed ones; in these 2D cases, the angular velocity has the same direction at any point of an extended region and, moreover, this velocity has the same distribution in all the planes orthogonal to it. Differential rotations can be seen as particular cases, in which, the angular velocity only depends on the distance to a rotation axis. Our results strongly suggest that appropriate mixtures of scalar and vector modes with very large spatial scales could explain the observed CMB anomalies.

We discuss the emission of gravitational radiation produced in encounters of dark matter galactic halos. To this aim we perform a number of numerical simulations of typical galaxy mergers, computing the associated gravitational radiation waveforms as well as the energy released in the processes. Our simulations yield dimensionless gravitational wave amplitudes of the order of $10^{-13}$ and gravitational wave frequencies of the order of $10^{-16}$ Hz, when the galaxies are located at a distance of 10 Mpc. These values are of the same order as those arising in the gravitational radiation originated by strong variations of the gravitational field in the early Universe, and therefore, such gravitational waves cannot be directly observed by ground-based detectors. We discuss the feasibility of an indirect detection by means of the B-mode polarization of the Cosmic Microwave Background (CMB) induced by such waves. Our results show that the gravitational waves from encounters of dark matter galactic halos leave much too small an imprint on the CMB polarization to be actually observed with ongoing and future missions.

Various effects produced by vector perturbations (vortical peculiar velocity fields) of a flat Friedmann-Robertson-Walker background are considered. In the presence of this type of perturbations, the polarization vector rotates. A formula giving the rotation angle is obtained and, then, it is used to prove that this angle depends on both the observation direction and the emission redshift. Hence, rotations are different for distinct quasars and also for the Cosmic Microwave Background (CMB) radiation coming along different directions (from distinct points of the last scattering surface). As a result of these rotations, some correlations could appear in an initially random field of quasar polarization orientations. Furthermore, the polarization correlations of the CMB could undergo alterations. Quasars and CMB maps are both considered in this paper. In the case of linear vector modes with very large spatial scales, the maximum rotation angles appear to be of a few degrees for quasars (located at redshifts z < 2.6) and a few tenths of degree for the CMB. These last rotations produce contributions to the B-mode of the CMB polarization which are too small to be observed with PLANCK (in the near future); however, these contributions are large enough to be observed with the next generation of satellites, which are being designed to detect the small B-mode produced by primordial gravitational waves.

Small maps of the Rees-Sciama (RS) effect are simulated by using an appropriate N-body code and a certain ray-tracing procedure. A method designed for the statistical analysis of cosmic microwave background (CMB) maps is applied to study the resulting simulations. These techniques, recently proposed --by our team-- to consider lens deformations of the CMB, are adapted to deal with the RS effect. This effect and the deviations from Gaussianity associated to it seem to be too small to be detected in the near future. This conclusion follows from our estimation of both the RS angular power spectrum and the RS reduced n-direction correlation functions for n<7.

A method for ray-tracing through n-body simulations has been recently proposed. It is based on a periodic universe covered by simulation boxes. Photons move along appropriate directions to avoid periodicity effects. Here, an improved version of this method is applied to simulate lensed CMB maps and maps of lens deformations. Particle mesh n-body simulations with appropriate boxes and resolutions are used to evolve the nonlinear inhomogeneities until present time. The resulting maps are statistically analyzed to look for deviations from Gaussianity. These deviations are measured --for the first time-- using correlations for configurations of n directions (n being smaller than 7). A wide range of angular scales are considered. Some interesting prospects are outlined.

We work in the framework of an inflationary cold dark matter universe with cosmological constant, in which the cosmological inhomogeneities are considered as gravitational lenses for the CMB photons. This lensing deforms the angular distribution of the CMB maps in such a way that the induced deformations are not Gaussian. Our main goal is the estimation of the deviations with respect to Gaussianity appeared in the distribution of deformations. In the new approach used in this paper, matter is evolved with a particle-mesh N-body code and, then, an useful ray-tracing technique designed to calculate the correlations of the lens deformations induced by nonlinear structures is applied. Our approach is described in detail and tested. Various correlations are estimated at an appropriate angular scale. The resulting values point out both deviations with respect to Gaussian statistics and a low level of correlation in the lens deformations.

In this paper, we focus our attention on the following question: How well can we recover the power spectrum of the cosmic microwave background from the maps of a given experiment?. Each experiment is described by a a pixelization scale, a beam size, a noise level and a sky coverage. We use accurate numerical simulations of the microwave sky and a cold dark matter model for structure formation in the universe. Angular scales smaller than those of previous simulations are included. The spectrum obtained from the simulated maps is appropriately compared with the theoretical one. Relative deviations between these spectra are estimated. Various contributions to these deviations are analyzed. The method used for spectra comparisons is discussed.

We describe and compare two types of microwave sky simulations which are good for small angular scales. The first type uses expansions in spherical harmonics, and the second one is based on plane waves and the Fast Fourier Transform. The angular power spectrum is extracted from maps corresponding to both types of simulations, and the resulting spectra are appropriately compared. In this way, the features and usefulness of Fourier simulations are pointed out. For $\ell \geq 100$, all the simulations lead to similar accuracies; however, the CPU cost of Fourier simulations is $\sim 10$ times smaller than that for spherical harmonic simulations. For $\ell \leq 100$, the simulations based on spherical harmonics seem to be preferable.

The Tolman-Bondi solution of the Einstein equations is used in order to model the time evolution of the void observed in Bo\"otes. The present density contrast of the central region ($\sim -0.75$) and its radius ($\sim 30h^{-1} \ Mpc$) are fixed, while the density parameter of the Universe, the amplitude of the density contrast inside the void wall, the width of this wall and the distance from the void centre to the Local Group are appropriately varied. The microwave background anisotropy produced by Bo\"otes-like voids is estimated for a significant set of locations. All the voids are placed far from the last scattering surface. It is shown that the anisotropy generated by these voids strongly depends on the density parameter, the wall structure and the void location. The Doppler dipole and quadrupole are subtracted and the residual anisotropy is calculated. In the case of some isolated Bo\"otes-like voids placed at redshifts between 1 and 10 in an open universe with density parameter $\Omega_{0}=0.2$, the residual anisotropy appears to be a few times $10^{-6}$ on scales of a few degrees. This anisotropy is about one order of magnitude greater than previous estimates corresponding to other cases. The anisotropy produced by a distribution of voids is qualitatively studied in the light of this result. Comparisons with previous estimates are discussed.