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.

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.

This paper is concerned with small angular scale experiments for the observation of cosmic microwave background anisotropies. In the absence of beam, the effects of partial coverage and pixelisation are disentangled and analyzed (using simulations). Then, appropriate maps involving the CMB signal plus the synchrotron and dust emissions from the Milky Way are simulated, and an asymmetric beam --which turns following different strategies-- is used to smooth the simulated maps. An associated circular beam is defined to estimate the deviations in the angular power spectrum produced by beam asymmetry without rotation and, afterwards, the deviations due to beam rotation are calculated. For a certain large coverage, the deviations due to pure asymmetry and asymmetry plus rotation appear to be very systematic (very similar in each simulation). Possible applications of the main results of this paper to data analysis in large coverage experiments --as PLANCK-- are outlined.

Non-linear evolution of cosmological energy density fluctuations triggers deviations from Gaussianity in the temperature distribution of the cosmic microwave background. A method to estimate these deviations is proposed. N-body simulations -- in a $\Lambda$CDM cosmology -- are used to simulate the strongly non-linear evolution of cosmological structures. It is proved that these simulations can be combined with the potential approximation to calculate the statistical moments of the CMB anisotropies produced by non-linear gravity. Some of these moments are computed and the resulting values are different from those corresponding to Gaussianity.

This paper is concerned with experiments which measure CMB anisotropies on small angular scales. A certain coverage, a beam structure and a level of uncorrelated noise define each experiment. We focus our atention on the reversion of the beam average. In each experiment, we look for the best pixelization for reversion, namely, for the pixelization that -after reversion- leads to good maps containing right spectra for the most wide range of angular scales. Squared pixels having different sizes "smaller" than the beam radius are considered. For a given size, the following question arises: How well can we assign a temperature to each pixel? Various mathematical methods are used to show that, in practice, this assignation -beam reversion or deconvolution- only leads to right spectra for pixel sizes greater than a certain lower limit close to the half of the beam radius. This limit is estimated for negligible and relevant levels of noise and also for spherically symmetric and asymmetric beams. After this general study, we focus our attention on two feasible detectors (which have been proposed to be on board of PLANCK satellite). For each of them, we estimate the size of the most appropriate pixelization compatible with beam reversion, difraction, observational strategy et cetera and, then, we answer the following question: Which is the part of the angular power spectrum which can be extracted from appropriately pixelized maps after deconvolution?

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.

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.