We discuss in the framework of the excursion set formalism a recent discovery from N-body simulations that the clustering of haloes of given mass depends on their formation history. We review why the standard implementation of this formalism is unable to explain such dependencies, and we show that this can, in principle, be rectified by implementing in full an ellipsoidal collapse model where collapse depends not only on the overdensity but also on the shape of the initial density field. We also present an alternative remedy for this deficiency, namely the inclusion of collapse barriers for pancakes and filaments, together with the assumption that formation history depends on when these barriers are crossed. We implement both these extensions in a generalised excursion set method, and run large Monte Carlo realisations to quantify the effects. Our results suggest that effects as large as those found in simulations can only arise in the excursion set formalism if the formation history of a halo does indeed depend on the size of its progenitor filaments and pancakes. We also present conditional distributions of progenitor pancakes and filaments for low-mass haloes identified at present epoch, and discuss a recent claim by Mo et.al. that most low-mass haloes were embedded in massive pancakes at $z\sim 2$.

We describe a method to estimate the mass distribution of a gravitational lens and the position of the sources from combined strong and weak lensing data. The algorithm combines weak and strong lensing data in a unified way producing a solution which is valid in both the weak and strong lensing regimes. We study how the result depends on the relative weighting of the weak and strong lensing data and on choice of basis to represent the mass distribution. We find that combining weak and strong lensing information has two major advantages: it eliminates the need for priors and/or regularization schemes for the intrinsic size of the background galaxies (this assumption was needed in previous strong lensing algorithms) and it corrects for biases in the recovered mass in the outer regions where the strong lensing data is less sensitive. The code is implemented into a software package called WSLAP (Weak & Strong Lensing Analysis Package) which is publicly available at http://darwin.cfa.harvard.edu/SLAP/

We revisit the issue of non-parametric gravitational lens reconstruction and present a new method to obtain the cluster mass distribution using strong lensing data without using any prior information on the underlying mass. The method relies on the decomposition of the lens plane into individual cells. We show how the problem in this approximation can be expressed as a system of linear equations for which a solution can be found. Moreover, we propose to include information about the null space. That is, make use of the pixels where we know there are no arcs above the sky noise. The only prior information is an estimation of the physical size of the sources. No priors on the luminosity of the cluster or shape of the halos are needed thus making the results very robust. In order to test the accuracy and bias of the method we make use of simulated strong lensing data. We find that the method reproduces accurately both the lens mass and source positions and provide error estimates.

We present the mass distribution in the central area of the cluster A1689 by fitting over 100 multiply lensed images with the non-parametric Strong Lensing Analysis Package (SLAP, Diego et al. 2004). The surface mass distribution is obtained in a robust way finding a total mass of 0.25E15 M_sun/h within a 70'' circle radius from the central peak. Our reconstructed density profile fits well an NFW profile with small perturbations due to substructure and is compatible with the more model dependent analysis of Broadhurst et al. (2004a) based on the same data. Our estimated mass does not rely on any prior information about the distribution of dark matter in the cluster. The peak of the mass distribution falls very close to the central cD and there is substructure near the center suggesting that the cluster is not fully relaxed. We also examine the effect on the recovered mass when we include the uncertainties in the redshift of the sources and in the original shape of the sources. Using simulations designed to mimic the data, we identify some biases in our reconstructed mass distribution. We find that the recovered mass is biased toward lower masses beyond 1 arcmin (150 kpc) from the central cD and that in the very center we may be affected by degeneracy problems. On the other hand, we confirm that the reconstructed mass between 25'' and 70'' is a robust, unbiased estimate of the true mass distribution and is compatible with an NFW profile.

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

We present a method for ultra-fast confrontation of the WMAP cosmic microwave background observations with theoretical models, implemented as a publicly available software package called CMBfit, useful for anyone wishing to measure cosmological parameters by combining WMAP with other observations. The method takes advantage of the underlying physics by transforming into a set of parameters where the WMAP likelihood surface is accurately fit by the exponential of a quartic or sextic polynomial. Building on previous physics based approximations by Hu et.al., Kosowsky et.al. and Chu et.al., it combines their speed with precision cosmology grade accuracy. A Fortran code for computing the WMAP likelihood for a given set of parameters is provided, pre-calibrated against CMBfast, accurate to Delta lnL ~ 0.05 over the entire 2sigma region of the parameter space for 6 parameter ``vanilla'' Lambda CDM models. We also provide 7-parameter fits including spatial curvature, gravitational waves and a running spectral index.

Despite the interest in dark matter and dark energy, it has never been shown that they are in fact two separate substances. We provide the first strong evidence that they are separate by ruling out a broad class of so-called unified dark matter models that have attracted much recent interest. We find that they produce oscillations or exponential blowup of the matter power spectrum inconsistent with observation. For the particular case of generalized Chaplygin gas models, 99.999% of the previously allowed parameter space is excluded, leaving essentially only the standard Lambda-CDM limit allowed.

We study inhomogeneous cosmological variations in the fine structure 'constant', $\alpha ,$ in Friedmann universes. Inhomogeneous motions of the scalar field driving changes in $\alpha $ display spatial oscillations that decrease in amplitude with increasing time. The inhomogeneous evolution quickly approaches that found for exact Friedmann universes. We prove a theorem to show that oscillations of $\alpha $ in time (or redshift) cannot occur in Friedmann universes in the BSBM theories considered here.

We formulate a simple extension of general relativity which incorporates space-time variations in the Newtonian gravitation 'constant', $G$, and the fine structure 'constant', $\alpha $, which generalises Brans-Dicke theory and our theory of varying $\alpha .$ We determine the behaviour of Friedmann universes in this theory. In the radiation and dust-dominated eras $\alpha G$ approaches a constant value and the rate of variation of $\alpha $ is equal to the magnitude of the rate of variation in $G$. The expansion dynamics of the universe are dominated by the variation of $G$ but the variation of $G$ has significant effects upon the time variation of $\alpha .$ Time variations in $\alpha $ are extinguished by the domination of the expansion by spatial curvature or quintessence fields, as in the case with no $G$ variation.

In some cosmological theories with varying constants there are anthropic reasons why the expansion of the universe must not be too {\it close} to flatness or the cosmological constant too close to zero. Using exact theories which incorporate time-variations in $\alpha $ and in $G$ we show how the presence of negative spatial curvature and a positive cosmological constant play an essential role in bringing to an end variations in the scalar fields driving time change in these 'constants' during any dust-dominated era of a universe's expansion. In spatially flat universes with $\Lambda =0$ the fine structure constant grows to a value which makes the existence of atoms impossible.

Is the recent evidence for a time-varying fine structure 'constant' $\alpha$ to be interpreted as a varying $e$, $c$, $\hbar$, or a combination thereof? We consider the simplest varying electric charge and varying speed of light theories (VSL) and prove that for the same type of dark matter they predict opposite senses of variation in $\alpha$ over cosmological times. We also show that unlike varying $e$ theories, VSL theories do not predict violations of the weak equivalence principle (WEP). Varying $e$ theories which explain astronomical inferences of varying $\alpha$ predict WEP violations only an order of magnitude smaller than existing E\"otv\"os experiment limits but could be decisively tested by STEP. We finally exhibit a set of atomic-clock and related experiments for which {\it all} (hyperbolic) varying $\alpha$ theories predict non-null results. They provide independent tests of the recent astronomical evidence.

We investigate the cosmological consequences of a simple theory in which the electric charge $e$ is allowed to vary. The theory is locally gauge and Lorentz invariant, and satisfies general covariance. We find that in this theory the fine structure 'constant', $\alpha $, remains almost constant in the radiation era, undergoes small increase in the matter era, but approaches a constant value when the universe starts accelerating because of the presence of a positive cosmological constant. This model satisfies geonuclear, nucleosynthesis, and CMB constraints on time-variation in $\alpha $, while fitting simultaneously the observed accelerating universe and the recent high-redshift evidence for small $\alpha $ variations in quasar spectra. The model also places specific restrictions on the nature of the dark matter. Further tests, involving stellar spectra and the E\"otv\"os experiment, are proposed.

We determine the behaviour of a time-varying fine structure 'constant' $\alpha (t)$ during the early and late phases of universes dominated by the kinetic energy of changing $\alpha (t)$, radiation, dust, curvature, and lambda, respectively. We show that after leaving an initial vacuum-dominated phase during which $\alpha$ increases, $\alpha$ remains constant in universes like our own during the radiation era, and then increases slowly, proportional to a logarithm of cosmic time, during the dust era. If the universe becomes dominated by negative curvature or a positive cosmological constant then $\alpha$ tends rapidly to a constant value. The effect of an early period of de Sitter or power-law inflation is to drive $\alpha$ to a constant value. Various cosmological consequences of these results are discussed with reference to recent observational studies of the value of $\alpha$ from quasar absorption spectra and to the existence of life in expanding universes.

We consider soliton solutions to Bekenstein's theory, for which the fine structure constant $\alpha=e^2/(4\pi\hbar c)$ is allowed to vary due to the presence of a dielectric field pervading the vacuum. More specifically we investigate the effects of a varying $\alpha$ upon a complex scalar field with a U(1) electromagnetic gauge symmetry subject to spontaneous symmetry breaking. We find vortex solutions to this theory, similar to the Nielsen-Olesen vortex. Near the vortex core the electric charge is typically much larger than far away from the string, lending these strings a superconducting flavour. In general the dielectric field coats the usual local string with a global string envelope. We discuss the cosmological implications of networks of such strings, with particular emphasis on their ability to generate inhomogeneous recombination scenarios. We also consider the possibility of the dielectric being a charged free field. Even though the vacuum of such a field is trivial, we find that the dielectric arranges itself in the shape of a local string, with a quantized magnetic flux at the core -- presumably borrowing these topological features from the underlying Nielsen-Olesen vortex.

We extend a previous bispectrum analysis of the COBE-DMR 4 year maps, allowing for the presence of correlations between all possible angular scales. We find that the non-Gaussian signal found earlier for bispectrum components including adjacent modes does not extend to triplets of modes with larger separations. Indeed for all separations $\Delta\ell >1$ we find that the COBE-DMR data is very Gaussian. The implication seems to be that the previously detected non-Gaussian scale-scale correlation falls off very quickly with mode separation.

We solve numerically the Boltzmann equation describing the evolution of a cosmic string network which contains only loops. In Minkowski space time the equilibrium solution predicted by statistical mechanics is recovered, and we prove that this solution is stable to non-linear perturbations provided that their energy does not exceed the critical energy for the Hagedorn transition. In expanding Einstein - de Sitter Universes we probe the distribution of loops with length much smaller than the horizon. For these loops we discover stable scaling solutions both in the radiation and matter dominated epochs. The shape of these solutions is very different in the two eras, with much higher energy density in the radiation epoch, and a larger average loop length in the matter epoch. These results suggest that if the conditions for formation of loop networks are indeed satisfied, these could in principle be good candidates for structure formation.