We investigate the possibility of dark matter being a pure geometrical effect, rather than a particle or a compact object, by exploring a specific modified gravity model: mimetic dark matter. We present an alternative formulation of the theory, closer to the standard cosmological perturbation theory framework. We make manifest the presence of arbitrary parameters and extra functions, both at background level and at first order in perturbation theory. We present the full set of independent equations of motion for this model, and we discuss the amount of tuning needed to match predictions of the theory to actual data. By using the matter power spectrum and cosmic microwave background angular power spectra as benchmark observables, we explicitly show that since there is no natural mechanism to generate adiabatic initial conditions in this specific model, extra fine-tuning is required. We modify the publicly available Boltzmann code \texttt{CLASS} to make accurate predictions for the observables in mimetic dark matter. Our modified version of \texttt{CLASS} is available on GitHub. We have used mimetic dark matter as an illustration of how much one is allowed to change the initial conditions before contradicting observations when modifying the laws of gravity as described by General Relativity but we point out that modifying gravity without providing a natural mechanism to generate adiabatic initial conditions will always lead to highly fine-tuned models.

In the light of the discovery of numerous (almost) dark galaxies from the ALFALAFA and LITTLE THINGS surveys, we revisit the predictions of Jimenez et al. 1997, based on the Toomre stability of rapidly-spinning gas disks. We have updated the predictions for $\Lambda$CDM with parameters given by Planck18, computing the expected number densities of dark objects, and their spin parameter and mass distributions. Comparing with the data is more challenging, but where the spins are more reliably determined, the spins are close to the threshold for disks to be stable according to the Toomre criterion, where the expected number density is highest, and reinforces the concept that there is a bias in the formation of luminous galaxies based on the spin of their parent halo.

Using the quantum information picture to describe the early universe as a time dependent quantum density matrix, with time playing the role of a stochastic variable, we compute the non-gaussian features in the distribution of primordial fluctuations. We use a quasi de Sitter model to compute the corresponding quantum Fisher information function as the second derivative of the relative entanglement entropy for the density matrix at two different times. We define the curvature fluctuations in terms of the time quantum estimator. Using standard quantum estimation theory we compute the non-gaussian features in the statistical distribution of primordial fluctuations. Our approach is model independent and only relies on the existence of a quasi de Sitter phase. We show that there are primordial non-gaussianities, both in the form of squeezed and equilateral shapes. The squeezed limit gives a value of $f_{\rm NL} \sim n_s-1$. In the equilateral limit we find that $f_{\rm NL} \sim 0.03$. The equilateral non-gaussianity is due to the non-linearity of Einstein's equation. On the other hand, the squeezed one is due to the quantum nature of clock synchronization and thus real and cannot be gauged away as a global curvature. We identify a new effect: {\it clock bias} which is a pure quantum effect and introduces a bias in the spectral tilt and running of the power spectrum of order $\sim 10^{-4}$, which could be potentially measurable and yield precious information on the quantum nature of the early Universe.

During the red giant phase, stars loose mass at the highest rate since birth. The mass-loss rate is not fixed, but varies from star-to-star by up to 5%, resulting in variations of the star's luminosity at the tip of the red giant branch (TRGB). Also, most stars, during this phase, engulf part of their planetary system, including their gas giant planets. Gas giant planet masses range between 0.1 to 2% of the host star mass. The engulfing of their gas giants planets can modify their luminosity at the TRGB, i.e. the point at which the He-core degeneracy is removed. We show that the increase in mass of the star by the engulfing of the gas giant planets only modifies the luminosity of a star at the TRGB by less than 0.1%, while metallicity can modify the luminosity of a star at the TRGB by up to 0.5%. However, the increase in turbulence of the convective envelope of the star, i.e., modification of the mixing length, has a more dramatic effect, on the star's luminosity, which we estimate could be as large as 5%. The effect is always in the direction to increase the turbulence and thus the mixing length which turns into a systematic decrease of the luminosity of the star at the TRGB. We find that the star-to-star variation of the mass-loss rate will dominate the variations in the luminosity of the TRGB with a contribution at the 5% level. If the star-to-star variation is driven by environmental effects --as it is reasonable to assume--, the same effects can potentially create an environmentally-driven mean effect on the luminosity of the tip of the red giant branch of a galaxy. Finally, we touch upon how to infer the frequency, and identify the engulfment, of exoplanets in low-metallicity RGB stars through high resolution spectroscopy as well as how to quantify mass loss rate distributions from the morphology of the horizontal branch.

We suggest a general relation between the position of the cosmic microwave background temperature power spectrum peaks and the inflationary slow roll parameter $\epsilon$. This relation is based on interpreting the variable setting the position of the peaks as the quantum distance between the end of inflation and recombination. This distance is determined by the primordial cosmological Fisher information introduced in arXiv:2002.04294. The observational constraints set by cosmic microwave background temperature data lead to a very stringent prediction for the value of the tensor-to-scalar ratio: $r=0.01 \pm 0.002$. Future polarization data of the cosmic microwave background should be able to measure this signal and corroborate or discard our model.

The evolution of differential ages of passive galaxies at different redshifts (cosmic chronometers) has been proved to be a method potentially able to constrain the Hubble parameter in a cosmology-independent way, but the systematic uncertainties must be carefully evaluated. In this paper, we compute the contribution to the full covariance matrix of systematic uncertainties due to the choice of initial mass function, stellar library, and metallicity, exploring a variety of stellar population synthesis models; the effect of a residual young component and star formation history was discussed elsewhere. Through simulations in the redshift range 0<z<1.5 we find that the choice of the stellar population synthesis model dominates the total error budget on H(z), with contributions at ~4.5% level, discarding the most discordant model. The contribution due to the choice of initial mass function is <0.5%, while that due to the stellar library is ~6.6% on average. In parallel, we also assess the impact of an uncertainty in the determination of the stellar metallicity, finding that an error of ~10% (5%) on the stellar metallicity propagates to a 9% (4%) error on H(z). These results are used to provide the combined contribution of these systematic effects on the error budget. For current H(z) measurements, where the uncertainties due to metallicity and star formation history were already included, we show that, using the more modern stellar libraries, the additional systematic uncertainty is between 5.4% (at z=0.2) and 2.3% (at z=1.5). To reach the goal of keeping the systematic error budget below the 1% level we discuss that efforts needed to obtain higher resolution and signal-to-noise spectra and improvements in the modeling of stellar population synthesis.

We describe inflation in terms of a time dependent quantum density matrix with time playing the role of a stochastic variable. Using a quasi-de Sitter model we compute the corresponding quantum Fisher information function as the second derivative of the relative entanglement entropy for the density matrix at two different times. Employing standard quantum estimation theory we evaluate the minimal variance of quantum scalar fluctuations that reproduces the power spectrum and the corresponding tilt in the slow roll limit. The Jeffreys prior associated with such Fisher information can be used to define the probabilities on the set of initial conditions defined by the slow roll parameter $\epsilon$ and the initial Shannon information.

We derive the dynamics at galaxy scales (Virial theorem) for modified gravity models in the general Hordensky class. Our aim is to develop tests that can distinguish whether dark matter is part of the theory of gravity or a particle. The recent discovery that a small number of dwarf galaxies lack any dark matter provides us with a stringent test for models of modified gravity. If there are only a few galaxies without dark matter, it is interesting to envision how to explain this if dark matter is a modification of gravity without adding epicycles to the model. We quantify this, and our main finding is that there will always be an extra contribution to the Virial theorem coming from the modification of gravity, even if a certain galaxy shows very small, if not negligible, trace of dark matter, as has been reported recently [1-3]. Thus, if these and more galaxies are confirmed as devoid (or negligible) of dark matter, while other similar galaxies have abundant dark matter, it seems interesting to find modifications of gravity to describe dark matter. Our result can be used by future astronomical surveys to put constraints on the parameters of modified gravity models at astrophysical scales where dark matter is described as such.

We propose a method to probe the homogeneity of a general universe, without assuming symmetry. We show that isotropy can be tested at remote locations on the past lightcone by comparing the line-of-sight and transverse expansion rates, using the time dependence of the polarization of Cosmic Microwave Background photons that have been inverse-Compton scattered by the hot gas in massive clusters of galaxies. This probes a combination of remote transverse and parallel components of the expansion rate of the metric, and we may use radial baryon acoustic oscillations or cosmic clocks to measure the parallel expansion rate. Thus we can test remote isotropy, which is a key requirement of a homogeneous universe. We provide explicit formulas that connect observables and properties of the metric.

The ages of the oldest stellar objects in our galaxy provide an independent test of the current cosmological model as they give a lower limit to the age of the Universe. Recent accurate parallaxes by the Gaia space mission, accurate measurements of the metallicity of stars, via individual elemental abundances, and advances in the modelling of stellar evolution, provide new, higher-precision age estimates of the oldest stellar populations in the galaxy: globular clusters and very-low-metallicity stars. The constraints on the age of the Universe, $t_U$, so obtained are determined from the local Universe and at late time. It is well known that local and early-Universe determinations of another cosmological parameter closely related to the age of the Universe, the Hubble constant $H_0$, show a $\gtrsim 3 \sigma$ tension. In the standard cosmological model, $\Lambda$CDM, $t_U$ and $H_0$ are related by the matter density parameter $\Omega_{m,0}$. We propose to combine local $t_U$ constraints with late-time $\Omega_{m,0}$ estimates in a $\Lambda$CDM framework, to obtain a low-redshift $H_0$ determination that does not rely on early Universe physics. A proof-of-principle of this approach with current data gives $H_0=71\pm2.8$ ($H_0= 69.3 \pm 2.7$) km s$^{-1}$ Mpc$^{-1}$ from globular clusters (very-low-metallicity stars) with excellent prospects for improved constraints in the near future.

The expansion history of the Universe can be constrained in a cosmology-independent way by measuring the differential age evolution of cosmic chronometers. This yields a measurement of the Hubble parameter $H(z)$ as a function of redshift. The most reliable cosmic chronometers known so far are extremely massive and passively evolving galaxies. Age-dating these galaxies is, however, a difficult task, and even a small contribution of an underlying young stellar population could, in principle, affect the age estimate and its cosmological interpretation. We present several spectral indicators to detect, quantify and constrain such contamination in old galaxies, and study how their combination can be used to maximize the purity of cosmic chronometers selection. In particular, we analyze the CaII H/K ratio, the presence (or absence) of H$\alpha$ and [OII] emission lines, higher order Balmer absorption lines, and UV flux; each indicator is especially sensitive to a particular age range, allowing us to detect young components ranging between 10 Myr and 1 Gyr. The combination of these indicators minimizes the contamination to a level below 1% in the case of ideal data. More importantly, it offers a way to control the systematic error on $H(z)$ as a function of the contamination by young stellar populations. We show that for our previous measurements of the Hubble parameter, the possible bias induced by the presence of a younger component is well below the current errors. We envision that these indicators will be instrumental in paving the road for a robust and reliable dating of the old population and its cosmological interpretation.

The determination of the inflationary energy scale represents one of the first step towards the understanding of the early Universe physics. The (very mild) non-Gaussian signals that arise from any inflation model carry information about the energy scale of inflation and may leave an imprint in some cosmological observables, for instance on the clustering of high-redshift, rare and massive collapsed structures. In particular, the graviton exchange contribution due to interactions between scalar and tensor fluctuations leaves a specific signature in the four-point function of curvature perturbations, thus on clustering properties of collapsed structures. We compute the contribution of graviton exchange on two- and three-point function of halos, showing that at large scales $k\sim 10^{-3}\ \mathrm{Mpc}^{-1}$ its magnitude is comparable or larger to that of other primordial non-Gaussian signals discussed in the literature. This provides a potential route to probe the existence of tensor fluctuations which is alternative and highly complementary to B-mode polarisation measurements of the cosmic microwave background radiation.

We compute in detail how deviations from Einstein gravity at the inflation energy scale could appear as non-Gaussian features in the sky. To illustrate this we use multi-field $\alpha-$attractor models in the framework of supergravity to realise inflation. We find no obvious obstacle for having choices of model's parameters that generate non-Gaussian features of the equilateral and local type at the $\cal O$(1) level in the $f_{\rm NL}$ non-Gaussianity parameter, thus being potentially detectable in future cosmological surveys. This non-Gaussianity has its origin in either the non-canonical kinetic term (which, in turn, is an immediate consequence of assuming an hyperbolic geometry of the moduli space), the interactions of the fields in the potential or the $\alpha-$parameter, or a combination of these three. This opens up the exciting possibility of constraining the law of gravity at energy scales close to the Planck one.

The number of observed dwarf galaxies, with dark matter mass $\lesssim 10^{11}$ M$_{\odot}$ in the Milky Way or the Andromeda galaxy does not agree with predictions from the successful $\Lambda$CDM paradigm. To alleviate this problem a suppression of dark matter clustering power on very small scales has been conjectured. However, the abundance of dark matter halos outside our immediate neighbourhood (the Local Group) seem to agree with the $\Lambda$CDM--expected abundance. Here we connect these problems to observations of weak lensing cosmic shear, pointing out that cosmic shear can make significant statements about the missing satellites problem in a statistical way. As an example and pedagogical application we use recent constraints on small-scales power suppression from measurements of the CFHTLenS data. We find that, on average, in a region of $\sim $Gpc$^3$ there is no significant small-scale power suppression. This implies that suppression of small-scale power is not a viable solution to the `missing satellites problem' or, alternatively, that on average in this volume there is no `missing satellites problem' for dark matter masses $\gtrsim 5 \times 10^9$ M$_{\odot}$. Further analysis of current and future weak lensing surveys will probe much smaller scales, $k > 10h$ Mpc$^{-1}$ corresponding roughly to masses $M < 10^9 M_{\odot}$.

An interesting test on the nature of the Universe is to measure the global spatial curvature of the metric in a model independent way, at a level of $|\Omega_k|<10^{-4}$, or, if possible, at the cosmic variance level of the amplitude of the CMB fluctuations $|\Omega_k|\approx10^{-5}$. A limit of $|\Omega_k|<10^{-4}$ would yield stringent tests on several models of inflation. Further, improving the constraint by an order of magnitude would help in reducing "model confusion" in standard parameter estimation. Moreover, if the curvature is measured to be at the value of the amplitude of the CMB fluctuations, it would offer a powerful test on the inflationary paradigm and would indicate that our Universe must be significantly larger than the current horizon. On the contrary, in the context of standard inflation, measuring a value above CMB fluctuations will lead us to conclude that the Universe is not much larger than the current observed horizon, this can also be interpreted as the presence of large fluctuations outside the horizon. However, it has proven difficult, so far, to find observables that can achieve such level of accuracy, and, most of all, be model-independent. Here we propose a method that can in principle achieve that, this is done by making minimal assumptions and using distance probes that are cosmology-independent: gravitational waves, redshift drift and cosmic chronometers. We discuss what kind of observations are needed in principle to achieve the desired accuracy.

The configuration of the three neutrino masses can take two forms, known as the normal and inverted hierarchies. We compute the Bayesian evidence associated with these two hierarchies. Previous studies found a mild preference for the normal hierarchy, and this was driven by the asymmetric manner in which cosmological data has confined the available parameter space. Here we identify the presence of a second asymmetry, which is imposed by data from neutrino oscillations. By combining constraints on the squared-mass splittings with the limit on the sum of neutrino masses of $\Sigma m_\nu < 0.13$ eV, and using a minimally informative prior on the masses, we infer odds of 42:1 in favour of the normal hierarchy, which is classified as "strong" in the Jeffreys' scale. We explore how these odds may evolve in light of higher precision cosmological data, and discuss the implications of this finding with regards to the nature of neutrinos. Finally the individual masses are inferred to be $m_1 = 3.80^{+26.2}_{-3.73} \, \text{meV}, m_2 = 8.8^{+18}_{-1.2} \, \text{meV}, m_3 = 50.4^{+5.8}_{-1.2} \, \text{meV}$ ($95\%$ credible intervals).

In this paper we discuss the commonly-used limiting cases, or approximations, for two-point cosmic shear statistics. We discuss the most prominent assumptions in this statistic: the flat-sky (small angle limit), the Limber (Bessel-to-delta function limit) and the Hankel transform (large l-mode limit) approximations; that the vast majority of cosmic shear results to date have used simultaneously. We find that the combined effect of these approximations can suppress power by >1% on scales of l<40. A fully non-approximated cosmic shear study should use a spherical-sky, non-Limber-approximated power spectrum analysis; and a transform involving Wigner small-d matrices in place of the Hankel transform. These effects, unaccounted for, would constitute at least 11% of the total budget for systematic effects for a power spectrum analysis of a Euclid-like experiment; but they are unnecessary.

Cosmological observables show a dependence with the neutrino mass, which is partially degenerate with parameters of extended models of gravity. We study and explore this degeneracy in Horndeski generalized scalar-tensor theories of gravity. Using forecasted cosmic microwave background and galaxy power spectrum datasets, we find that a single parameter in the linear regime of the effective theory dominates the correlation with the total neutrino mass. For any given mass, a particular value of this parameter approximately cancels the power suppression due to the neutrino mass at a given redshift. The extent of the cancellation of this degeneracy depends on the cosmological large-scale structure data used at different redshifts. We constrain the parameters and functions of the effective gravity theory and determine the influence of gravity on the determination of the neutrino mass from present and future surveys.

We investigate the spin evolution of dark matter haloes and their dependence on the number of connected filaments from the cosmic web at high redshift (spin-filament relation hereafter). To this purpose, we have simulated $5000$ haloes in the mass range $5\times10^{9}h^{-1}M_{\odot}$ to $5\times10^{11}h^{-1}M_{\odot}$ at $z=3$ in cosmological N-body simulations. We confirm the relation found by Prieto et al. 2015 where haloes with fewer filaments have larger spin. We also found that this relation is more significant for higher halo masses, and for haloes with a passive (no major mergers) assembly history. Another finding is that haloes with larger spin or with fewer filaments have their filaments more perpendicularly aligned with the spin vector. Our results point to a picture in which the initial spin of haloes is well described by tidal torque theory and then gets subsequently modified in a predictable way because of the topology of the cosmic web, which in turn is given by the currently favoured LCDM model. Our spin-filament relation is a prediction from LCDM that could be tested with observations.

We investigate our knowledge of early universe cosmology by exploring how much additional energy density can be placed in different components beyond those in the $\Lambda$CDM model. To do this we use a method to separate early- and late-universe information enclosed in observational data, thus markedly reducing the model-dependency of the conclusions. We find that the 95\% credibility regions for extra energy components of the early universe at recombination are: non-accelerating additional fluid density parameter $\Omega_{\rm MR} < 0.006$ and extra radiation parameterised as extra effective neutrino species $2.3 < N_{\rm eff} < 3.2$ when imposing flatness. Our constraints thus show that even when analyzing the data in this largely model-independent way, the possibility of hiding extra energy components beyond $\Lambda$CDM in the early universe is seriously constrained by current observations. We also find that the standard ruler, the sound horizon at radiation drag, can be well determined in a way that does not depend on late-time Universe assumptions, but depends strongly on early-time physics and in particular on additional components that behave like radiation. We find that the standard ruler length determined in this way is $r_{\rm s} = 147.4 \pm 0.7$ Mpc if the radiation and neutrino components are standard, but the uncertainty increases by an order of magnitude when non-standard dark radiation components are allowed, to $r_{\rm s} = 150 \pm 5$ Mpc.

Assuming the existence of standard rulers, standard candles and standard clocks, requiring only the cosmological principle, a metric theory of gravity, a smooth expansion history, and using state-of-the-art observations, we determine the length of the "low-redshift standard ruler". The data we use are a compilation of recent Baryon acoustic oscillation data (relying on the standard ruler), Type 1A supernov\ae\ (as standard candles), ages of early type galaxies (as standard clocks) and local determinations of the Hubble constant (as a local anchor of the cosmic distance scale). In a standard $\Lambda$CDM cosmology the "low-redshift standard ruler" coincides with the sound horizon at radiation drag, which can also be determined --in a model dependent way-- from CMB observations. However, in general, the two quantities need not coincide. We obtain constraints on the length of the low-redshift standard ruler: $r^h_{\rm s}=101.0 \pm 2.3 h^{-1}$ Mpc, when using only Type 1A supernov\ae\ and Baryon acoustic oscillations, and $r_{\rm s}=150.0\pm 4.7 $ Mpc when using clocks to set the Hubble normalisation, while $r_{\rm s}=141.0\pm 5.5 $ Mpc when using the local Hubble constant determination (using both yields $r_{\rm s}=143.9\pm 3.1 $ Mpc). The low-redshift determination of the standard ruler has an error which is competitive with the model-dependent determination from cosmic microwave background measurements made with the {\em Planck} satellite, which assumes it is the sound horizon at the end of baryon drag.

We demonstrate that a scalar field is unable to reverse its direction of motion while continuously exchanging energy with another fluid. If the rate of transfer is modulated by the scalar's acceleration, the field can undergo a rapid process of freezing, despite being displaced from the local minimum of its potential. This enables dark energy to form from any potential, regardless of its shape. The field's equation of state mimicks that of a cosmological constant. We present a physically motivated realisation in the form of a derivative neutrino-majoron coupling. Coherent motions, which form only once the neutrinos become non-relativistic, could be responsible for instigating the freezing process. This would provide a natural resolution to the dark energy coincidence problem, while avoiding the dynamical instabilities associated with mass-varying neutrino models. Finally we discuss possible means by which this model could be experimentally verified.

Recent anomalies found in cosmological datasets such as the low multipoles of the Cosmic Microwave Background or the low redshift amplitude and growth of clustering measured by e.g., abundance of galaxy clusters and redshift space distortions in galaxy surveys, have motivated explorations of models beyond standard $\Lambda$CDM. Of particular interest are models where general relativity (GR) is modified on large cosmological scales. Here we consider deviations from $\Lambda$CDM+GR within the context of Horndeski gravity, which is the most general theory of gravity with second derivatives in the equations of motion. We adopt a parametrization in which the four additional Horndeski functions of time $\alpha_i(t)$ are proportional to the cosmological density of dark energy $\Omega_{DE}(t)$. Constraints on this extended parameter space using a suite of state-of-the art cosmological observations are presented for the first time. Although the theory is able to accommodate the low multipoles of the Cosmic Microwave Background and the low amplitude of fluctuations from redshift space distortions, we find no significant tension with $\Lambda$CDM+GR when performing a global fit to recent cosmological data and thus there is no evidence against $\Lambda$CDM+GR from an analysis of the value of the Bayesian evidence ratio of the modified gravity models with respect to $\Lambda$CDM, despite introducing extra parameters. The posterior distribution of these extra parameters that we derive return strong constraints on any possible deviations from $\Lambda$CDM+GR in the context of Horndeski gravity. We illustrate how our results can be applied to a more general frameworks of modified gravity models.

Deriving the expansion history of the Universe is a major goal of modern cosmology. To date, the most accurate measurements have been obtained with Type Ia Supernovae and Baryon Acoustic Oscillations, providing evidence for the existence of a transition epoch at which the expansion rate changes from decelerated to accelerated. However, these results have been obtained within the framework of specific cosmological models that must be implicitly or explicitly assumed in the measurement. It is therefore crucial to obtain measurements of the accelerated expansion of the Universe independently of assumptions on cosmological models. Here we exploit the unprecedented statistics provided by the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 to provide new constraints on the Hubble parameter $H(z)$ using the em cosmic chronometers approach. We extract a sample of more than 130000 of the most massive and passively evolving galaxies, obtaining five new cosmology-independent $H(z)$ measurements in the redshift range $0.3<z<0.5$, with an accuracy of $\sim$11-16\% incorporating both statistical and systematic errors. Once combined, these measurements yield a 6\% accuracy constraint of $H(z=0.4293)=91.8\pm5.3$ km/s/Mpc. The new data are crucial to provide the first cosmology-independent determination of the transition redshift at high statistical significance, measuring $z_{t}=0.4\pm0.1$, and to significantly disfavor the null hypothesis of no transition between decelerated and accelerated expansion at 99.9\% confidence level. This analysis highlights the wide potential of the cosmic chronometers approach: it permits to derive constraints on the expansion history of the Universe with results competitive with standard probes, and most importantly, being the estimates independent of the cosmological model, it can constrain cosmologies beyond -and including- the $\Lambda$CDM model.

We use the latest compilation of observational H(z) measurements obtained with cosmic chronometers in the redshift range $0<z<2$ to place constraints on cosmological parameters. We consider the sample alone and in combination with other state-of-the art cosmological probes: CMB data from the latest Planck 2015 release, the most recent estimate of the Hubble constant $H_{0}$, a compilation of recent BAO data, and the latest SNe sample. Since cosmic chronometers are independent of the assumed cosmological model, we are able to provide constraints on the parameters that govern the expansion history of the Universe in a way that can be used to test cosmological models. We show that the H(z) measurements obtained with cosmic chronometer from the BOSS survey provide enough constraining power in combination with CMB data to constrain the time evolution of dark energy, yielding constraints competitive with those obtained using SNe and/or BAO. From late-Universe probes alone we find that $w_0=-0.9\pm0.18$ and $w_a=-0.5\pm1.7$, and when combining also CMB data we obtain $w_0=-0.98\pm0.11$and $w_a=-0.30\pm0.4$. These new constraints imply that nearly all quintessence models are disfavoured, only phantom models or a pure cosmological constant being allowed. For the curvature we find $\Omega_k=0.003\pm0.003$, including CMB data. Cosmic chronometers data are important also to constrain neutrino properties by breaking or reducing degeneracies with other parameters. We find that $N_{eff}=3.17\pm0.15$, thus excluding the possibility of an extra (sterile) neutrino at more than $5\sigma$, and put competitive limits on the sum of neutrino masses, $\Sigma m_{\nu}< 0.27$ eV at 95% confidence level. Finally, we constrain the redshift evolution of dark energy, and find w(z) consistent with the $\Lambda$CDM model at the 40% level over the entire redshift range $0<z<2$. [abridged]

There is currently a discrepancy in the measured value of the amplitude of matter clustering, parameterised using $\sigma_8$, inferred from galaxy weak lensing, and cosmic microwave background data, which could be an indication of new physics, such as massive neutrinos or a modification to the gravity law, or baryon feedback. In this paper we make the assumption that the cosmological parameters are well determined by Planck, and use weak lensing data to investigate the implications for baryon feedback and massive neutrinos, as well as possible contributions from intrinsic alignments and biases in photometric redshifts. We apply a non-parametric approach to model the baryonic feedback on the dark matter clustering, which is flexible enough to reproduce the OWLS and Illustris simulation results. The statistic we use, 3D cosmic shear, is a method that extracts cosmological information from weak lensing data using a spherical-Bessel function power spectrum approach. We analyse the CFHTLenS weak lensing data and, assuming best fit cosmological parameters from the Planck CMB experiment, find that there is no evidence for baryonic feedback on the dark matter power spectrum, but there is evidence for a bias in the photometric redshifts in the CFHTLenS data, consistent with a completely independent analysis by Choi et al. (2015), based on spectroscopic redshifts; and that these conclusions are robust to assumptions about the intrinsic alignment systematic. We also find an upper limit on the sum of neutrino masses conditional on other $\Lambda$CDM parameters being fixed, of $< 0.28$ eV ($1\sigma$).

Recent constrains on the sum of neutrino masses inferred by analyzing cosmological data, show that detecting a non-zero neutrino mass is within reach of forthcoming cosmological surveys, implying a direct determination of the absolute neutrino mass scale. The measurement relies on constraining the shape of the matter power spectrum below the neutrino free streaming scale: massive neutrinos erase power at these scales. Detection of a lack of small-scale power, however, could also be due to a host of other effects. It is therefore of paramount importance to validate neutrinos as the source of power suppression at small scales. We show that, independent on hierarchy, neutrinos always show a footprint on large, linear scales; the exact location and properties can be related to the measured power suppression (an astrophysical measurement) and atmospheric neutrinos mass splitting (a neutrino oscillation experiment measurement). This feature can not be easily mimicked by systematic uncertainties or modifications in the cosmological model. The measurement of such a feature, up to 1% relative change in the power spectrum, is a smoking gun for confirming the determination of the absolute neutrino mass scale from cosmological observations. It also demonstrates the synergy of astrophysics and particle physics experiments.

Galactic Gamma-Ray Bursts (GRBs) are copious sources of gamma-rays that can pose a threat to complex life. Using recent determinations of their rate and the probability of GRBs causing massive extinction, we explore what type of universes are most likely to harbour advanced forms of life. For this purpose we use cosmological N-body simulations to determine at what time and for what value of the cosmological constant ($\Lambda$) the chances of life being unaffected by cosmic explosions are maximised. We find that $\Lambda-$dominated universes favour the survival of life against GRBs. Within a $\Lambda$CDM cosmology, the parameters that govern the likelihood of life survival to GRBs are dictated by the value of $\Lambda$ and the age of the Universe. We find that we seem to live in a favorable point in this parameter phase space which minimises the exposure to cosmic explosions, yet maximises the number of main sequence (hydrogen-burning) stars around which advanced life forms can exist.

Applying a transformation to a non-Gaussian field can enhance the information content of the resulting power spectrum, by reducing the correlations between Fourier modes. In the context of weak gravitational lensing, it has been shown that this gain in information content is significantly compromised by the presence of shape noise. We apply clipping to mock convergence fields, a technique which is known to be robust in the presence of noise and has been successfully applied to galaxy number density fields. When analysed in isolation the resulting convergence power spectrum returns degraded constraints on cosmological parameters. However substantial gains can be achieved by performing a combined analysis of the power spectra derived from both the original and transformed fields. Even in the presence of realistic levels of shape noise, we demonstrate that this approach is capable of reducing the area of likelihood contours within the $\Omega_m - \sigma_8$ plane by more than a factor of three.

Old Globular Clusters (GCs) in the Milky Way have ages of about 13 Gyr, placing their formation time in the reionization epoch. We propose a novel scenario for the formation of these systems based on the merger of two or more atomic cooling halos at high-redshift (z>6). First generation stars are formed as an intense burst in the center of a minihalo that grows above the threshold for hydrogen cooling (halo mass M_h~10^8 Msun) by undergoing a major merger within its cooling timescale (~150 Myr). Subsequent minor mergers and sustained gas infall bring new supply of pristine gas at the halo center, creating conditions that can trigger new episodes of star formation. The dark-matter halo around the GC is then stripped during assembly of the host galaxy halo. Minihalo merging is efficient only in a short redshift window, set by the LCDM parameters, allowing us to make a strong prediction on the age distribution for old GCs. From cosmological simulations we derive an average merging redshift <z>=9 and narrow distribution Dz=2, implying average GC age <t_age>=13.0+/-0.2 Gyr including ~0.2 Gyr of star formation delay. Qualitatively, our scenario reproduces other general old GC properties (characteristic masses and number of objects, metallicity versus galactocentric radius anticorrelation, radial distribution), but unlike age, these generally depend on details of baryonic physics. In addition to improved age measurements, direct validation of the model at z~10 may be within reach of ultradeep gravitationally lensed observations with the James Webb Space Telescope.

We present a detailed study of second-order matter perturbations for the general Horn- deski class of models. Being the most general scalar-tensor theory having second-order equations of motion, it includes many known gravity and dark energy theories and General Relativity with a cosmological constant as a specific case. This enables us to estimate the leading order dark matter bispectrum generated at late-times by gravitational instability. We parametrize the evolution of the first and second-order equations of motion as proposed by Bellini and Sawicki (2014), where the free functions of the theory are assumed to be proportional to the dark energy density. We show that it is unnatural to have large 10% ( 1%) deviations of the bispectrum introducing even larger ~ 30% (~ 5%) deviations in the linear growth rate. Considering that measurements of the linear growth rate have much higher signal-to-noise than bispectrum measurements, this indicates that for Horndeski models which reproduce the expansion history and the linear growth rate as predicted by GR the dark matter bispectrum kernel can be effectively modelled as the standard GR one. On the other hand, an observation of a large bispectrum deviation that can not be explained in terms of bias would imply either that the evolution of perturbations is strongly different than the evolution predicted by GR or that the theory of gravity is exotic (e.g., breaks the weak equivalence principle) and/or fine-tuned.

We exploit cosmological-model independent measurements of the expansion history of the Universe to provide a cosmic distance ladder. These are supernovae type Ia used as standard candles (at redshift between 0.01 and 1.3) and baryon acoustic oscillations (at redshifts between 0.1 and 0.8) as standard rulers. We calibrate (anchor) the ladder in two ways: first using the local $H_0$ value as an anchor at $z$ = 0 (effectively calibrating the standard candles) and secondly using the cosmic microwave background-inferred sound-horizon scale as an anchor (giving the standard ruler length) as an inverse distance ladder. Both methods are consistent, but the uncertainty in the expansion history $H(z)$ is smaller if the sound horizon scale is used. We present inferred values for the sound horizon at radiation drag $r_d$ which do not rely on assumptions about the early expansion history nor on cosmic microwave background measurements but on the cosmic distance ladder and baryon acoustic oscillations measurements. We also present derived values of $H_0$ from the inverse distance ladder and we show that they are in very good agreement with the extrapolated value in a $\Lambda$CDM model from Planck cosmic microwave background data.

Long-Duration Gamma-Ray Bursts (GRBs) are powerful probes of the Universe star formation history, but correlation between the two depends on the highly debated presence/strength of a metallicity bias. To investigate this correlation, we use a phenomenological model that successfully describes star formation rates, luminosities and stellar masses of star forming galaxies, applying it to GRB production. We predict luminosities, stellar masses, and metallicities of host galaxies depending on the metallicity bias. Our best-fitting model includes a moderate metallicity bias, broadly consistent with the large majority of long-duration GRBs in metal-poor environments originating from collapsars (probability ~83%), but with a secondary contribution (~17%) from metal-independent production channels, such as binary evolution. Because of the mass-metallicity relation of galaxies, the maximum likelihood model predicts that the metal-independent channel becomes dominant at z<2, where hosts have higher metallicities and collapsars are suppressed. This possibly explains why some studies find no clear evidence of a metal-bias based on low-z samples. However, while metallicity predictions match observations well at high redshift, there is tension with low redshift observations, since a significant fraction of GRB hosts are predicted to have (near-)solar metallicity. This is in contrast to observations, unless obscured, metal-rich hosts are preferentially missed in current datasets, and suggests that lower efficiencies of the metal-independent GRB channel might be preferred following a comprehensive fit from complete samples. Overall, we are able to establish the presence of a metallicity bias for GRB production, but continued characterization of GRB host galaxies is needed to quantify its strength.

We measure the length of the Baryon Acoustic Oscillation (BAO) feature, and the expansion rate of the recent Universe, from low-redshift data only, almost model-independently. We make only the following minimal assumptions: homogeneity and isotropy; a metric theory of gravity; a smooth expansion history, and the existence of standard candles (supernov\ae) and a standard BAO ruler. The rest is determined by the data, which are compilations of recent BAO and Type IA supernova results. Making only these assumptions, we find for the first time that the standard ruler has length $103.9 \pm 2.3\, h^{-1}$ Mpc. The value is a measurement, in contrast to the model-dependent theoretical prediction determined with model parameters set by Planck data ($99.3 \pm 2.1 \, h^{-1}$ Mpc). The latter assumes $\Lambda$CDM, and that the ruler is the sound horizon at radiation drag. Adding passive galaxies as standard clocks or a local Hubble constant measurement allows the absolute BAO scale to be determined ($142.8\pm 3.7$ Mpc), and in the former case the additional information makes the BAO length determination more precise ($101.9\pm 1.9 \, h^{-1}\,$Mpc). The inverse curvature radius of the Universe is weakly constrained and consistent with zero, independently of the gravity model, provided it is metric. We find the effective number of relativistic species to be $N_{\rm eff} = 3.53\pm 0.32$, independent of late-time dark energy or gravity physics.

As a copious source of gamma-rays, a nearby Galactic Gamma-Ray Burst (GRB) can be a threat to life. Using recent determinations of the rate of GRBs, their luminosity function and properties of their host galaxies, we estimate the probability that a life-threatening (lethal) GRB would take place. Amongst the different kinds of GRBs, long ones are most dangerous. There is a very good chance (but no certainty) that at least one lethal GRB took place during the past 5 Gyr close enough to Earth as to significantly damage life. There is a 50% chance that such a lethal GRB took place during the last 500 Myr causing one of the major mass extinction events. Assuming that a similar level of radiation would be lethal to life on other exoplanets hosting life, we explore the potential effects of GRBs to life elsewhere in the Galaxy and the Universe. We find that the probability of a lethal GRB is much larger in the inner Milky Way (95% within a radius of 4 kpc from the galactic center), making it inhospitable to life. Only at the outskirts of the Milky Way, at more than 10 kpc from the galactic center, this probability drops below 50%. When considering the Universe as a whole, the safest environments for life (similar to the one on Earth) are the lowest density regions in the outskirts of large galaxies and life can exist in only ~ 10% of galaxies. Remarkably, a cosmological constant is essential for such systems to exist. Furthermore, because of both the higher GRB rate and galaxies being smaller, life as it exists on Earth could not take place at $z > 0.5$. Early life forms must have been much more resilient to radiation.

In order to elucidate the origin of spin in both dark matter and baryons in galaxies, we have performed hydrodynamical simulations from cosmological initial conditions. We study atomic cooling haloes in the redshift range $100 > z > 9$ with masses of order $10^9{\rm M_{\odot}}$ at redshift $z=10$. We assume that the gas has primordial composition and that ${\rm H_2}$-cooling and prior star-formation in the haloes have been suppressed. We present a comprehensive analysis of the gas and dark matter properties of four halos with very low ($\lambda \approx 0.01$), low ($\lambda \approx 0.04$), high ($\lambda \approx 0.06$) and very high ($\lambda \approx 0.1$) spin parameter. Our main conclusion is that the spin orientation and magnitude is initially well described by tidal torque linear theory, but later on is determined by the merging and accretion history of each halo. We provide evidence that the topology of the merging region, i.e. the number of colliding filaments, gives an accurate prediction for the spin of dark matter and gas: halos at the center of knots will have low spin while those in the center of filaments will have high spin. The spin of a halo is given by $\lambda \approx 0.05 \times \left(\frac{7.6}{\rm number\,\,\, of \,\,\, filaments}\right)^{5.1}$.

We show that modifications of Einstein gravity during inflation could leave potentially measurable imprints on cosmological observables in the form of non-Gaussian perturbations. This is due to the fact that these modifications appear in the form of an extra field that could have non-trivial interactions with the inflaton. We show it explicitly for the case $R+\alpha R^2$, where nearly scale-invariant non-Gaussianity at the level of $f_{\rm NL} \approx -(1$ to $30)$ can be obtained, in a quasi-local configuration.

We give an analytical form for the weighted correlation function of peaks in a Gaussian random field. In a cosmological context, this approach strictly describes the formation bias and is the main result here. Nevertheless, we show its validity and applicability to the evolved cosmological density field and halo field, using Gaussian random field realisations and dark matter N-body numerical simulations. Using this result from peak theory we compute the bias of peaks (and dark matter halos) and show that it reproduces results from the simulations at the ${\mathcal O}(10\%)$ level. Our analytical formula for the bias predicts a scale-dependent bias with two characteristics: a broad band shape which, however, is most affected by the choice of weighting scheme and evolution bias, and a more robust, narrow feature localised at the BAO scale, an effect that is confirmed in simulations. This scale-dependent bias smooths the BAO feature but, conveniently, does not move it. We provide a simple analytic formula to describe this effect. We envision that our analytic solution will be of use for galaxy surveys that exploit galaxy clustering.

We use cosmology-independent measurements of the expansion history in the redshift range 0.1 < z <1.2 and compare them with the Cosmic Microwave Background-derived expansion history predictions. The motivation is to investigate if the tension between the local (cosmology independent) Hubble constant H0 value and the Planck-derived H0 is also present at other redshifts. We conclude that there is no tension between Planck and cosmology independent-measurements of the Hubble parameter H(z) at 0.1 < z < 1.2 for the LCDM model (odds of tension are only 1:15, statistically not significant). Considering extensions of the LCDM model does not improve these odds (actually makes them worse), thus favouring the simpler model over its extensions. On the other hand the H(z) data are also not in tension with the local H0 measurements but the combination of all three data-sets shows a highly significant tension (odds ~ 1:400). Thus the new data deepen the mystery of the mismatch between Planck and local H0 measurements, and cannot univocally determine wether it is an effect localised at a particular redshift. Having said this, we find that assuming the NGC4258 maser distance as the correct anchor for H0, brings the odds to comfortable values. Further, using only the expansion history measurements we constrain, within the LCDM model, H0 = 68.5 +- 3.5 and Omega_m = 0.32 +- 0.05 without relying on any CMB prior. We also address the question of how smooth the expansion history of the universe is given the cosmology independent data and conclude that there is no evidence for deviations from smoothness on the expansion history, neither variations with time in the value of the equation of state of dark energy.

We use recent measurements of the expansion history of the universe to place constraints on the parameter space of cubic Galileon models, in particular we concentrate on those models which contain the simplest Galileon term plus a linear potential. This gives strong constraints on the Lagrangian of these models. Most dynamical terms in the Galileon Lagrangian are constraint to be small and the acceleration is effectively provided by a constant term in the scalar potential, thus reducing, effectively, to a LCDM model for current acceleration. The effective equation of state is indistinguishable from that of a cosmological constant w = -1 and the data constraint it to have no temporal variations of more than at the few % level. The energy density of the Galileon can contribute only to about 10% of the acceleration energy density, being the other 90% a cosmological constant term. This demonstrates how useful direct measurements of the expansion history of the universe are at constraining the dynamical nature of dark energy.

There has been significant recent progress in observational cosmology. This, in turn, has provided an unprecedented picture of the early universe and its evolution. In this review I will present a (biased) view of how one can use these observational results to constraint fundamental physics and in particular physics beyond the standard model.

It has been shown by Shchekinov & Vasiliev2006 (SV06) that HD molecules can be an important cooling agent in high redshift z >10 haloes if they undergo mergers under specific conditions so suitable shocks are created. Here we build upon Prieto et al. (2012) who studied in detail the merger-generated shocks, and show that the conditions for HD cooling can be studied by combining these results with a suite of dark-matter only simulations. We have performed a number of dark matter only simulations from cosmological initial conditions inside boxes with sizes from 1 to 4 Mpc. We look for haloes with at least two progenitors of which at least one has mass M > M_cr (z), where M_cr (z) is the SV06 critical mass for HD over-cooling. We find that the fraction of over-cooled haloes with mass between M_cr (z) and 10^{0.2} M_cr (z), roughly below the atomic cooling limit, can be as high as ~ 0.6 at z ~ 10 depending on the merger mass ratio. This fraction decreases at higher redshift reaching a value ~0.2 at z ~ 15. For higher masses, i.e. above 10^{0.2} M_cr (z) up to 10^{0.6} M_cr (z), above the atomic cooling limit, this fraction rises to values ~ 0.8 until z ~ 12.5. As a consequence, a non negligible fraction of high redshift z > 10 mini-haloes can drop their gas temperature to the Cosmic Microwave Background temperature limit allowing the formation of low mass stars in primordial environments.

We show how general initial conditions for small field inflation can be obtained in multi-field models. This is provided by non-linear angular friction terms in the inflaton that provide a phase of non-slow-roll inflation before the slow-roll inflation phase. This in turn provides a natural mechanism to star small-field slow-roll at nearly zero velocity for arbitrary initial conditions. We also show that there is a relation between the scale of SUSY breaking sqrt (f) and the amount of non-gaussian fluctuations generated by the inflaton. In particular, we show that in the local non-gaussian shape there exists the relation sqrt (f) = 10^{13} GeV sqrt (f_NL). With current observational limits from Planck, and adopting the minimum amount of non-gaussian fluctuations allowed by single-field inflation, this provides a very tight constraint for the SUSY breaking energy scale sqrt (f) = 3-7 x 10^{13} GeV at 95% confidence. Further limits, or detection, from next year's Planck polarisation data will further tighten this constraint by a factor of two. We highlight that the key to our approach is to identify the inflaton with the scalar component of the goldstino superfield. This superfield is universal and implements the dynamics of SUSY breaking as well as superconformal breaking.

We use the latest Planck constraints, and in particular constraints on the derived parameters (Hubble constant and age of the Universe) for the local universe and compare them with local measurements of the same quantities. We propose a way to quantify whether cosmological parameters constraints from two different experiments are in tension or not. Our statistic, T, is an evidence ratio and therefore can be interpreted with the widely used Jeffrey's scale. We find that in the framework of the LCDM model, the Planck inferred two dimensional, joint, posterior distribution for the Hubble constant and age of the Universe is in "strong" tension with the local measurements; the odds being ~ 1:50. We explore several possibilities for explaining this tension and examine the consequences both in terms of unknown errors and deviations from the LCDM model. In some one-parameter LCDM model extensions, tension is reduced whereas in other extensions, tension is instead increased. In particular, small total neutrino masses are favored and a total neutrino mass above 0.15 eV makes the tension "highly significant" (odds ~ 1:150). A consequence of accepting this interpretation of the tension is that the degenerate neutrino hierarchy is highly disfavoured by cosmological data and the direct hierarchy is slightly favored over the inverse.

While there are numerous indications that GRBs arise from the death of massive stars, the GRB rate does not follow the global cosmic star formation rate and, within their hosts, GRBs are more concentrated in regions of very high star formation. We explain both puzzles here. Using the publicly available VESPA database of SDSS Data Release 7 spectra, we explore a multi-parameter space in galaxy properties, like stellar mass, metallicity, dust etc. to find the sub-set of galaxies that reproduce the recently obtained GRB rate by Wanderman & Piran (2010). We find that only galaxies with present stellar masses below < 10^{10} Msun and low metallicity reproduce the observed GRB rate. This is consistent with direct observations of GRB hosts and provides an independent confirmation of the nature of GRB hosts. Because of the significantly larger sample of SDSS galaxies, we compute their correlation function and show that they are anti-biased with respect to the dark matter: they are in filaments and voids. Using recent observations of massive stars in local dwarfs we show how the fact that GRB hosts galaxies are dwarfs can explain the observation that GRBs are more concentrated in regions of high star formation than SNe. Finally we explain these results using new theoretical advances in the field of star formation.

We provide an exact expression for the multi-variate joint probability distribution function of non-Gaussian fields primordially arising from local transformations of a Gaussian field. This kind of non-Gaussianity is generated in many models of inflation. We apply our expression to the non- Gaussianity estimation from Cosmic Microwave Background maps and the halo mass function where we obtain analytical expressions. We also provide analytic approximations and their range of validity. For the Cosmic Microwave Background we give a fast way to compute the PDF which is valid up to 7{\sigma} for fNL values (both true and sampled) not ruled out by current observations, which consists of expressing the PDF as a combination of bispectrum and trispectrum of the temperature maps. The resulting expression is valid for any kind of non-Gaussianity and is not limited to the local type. The above results may serve as the basis for a fully Bayesian analysis of the non-Gaussianity parameter.

We explore how local, cosmology-independent measurements of the Hubble constant and the age of the Universe help to provide a powerful consistency check of the currently favored cosmological model (flat LambdaCDM) and model-independent constraints on cosmology. We use cosmic microwave background (CMB) data to define the model-dependent cosmological parameters, and add local measurements to assess consistency and determine whether extensions to the model are justified. At current precision, there is no significant tension between the locally measured Hubble constant and age of the Universe (with errors of 3% and 5% respectively) and the corresponding parameters derived from the CMB. However, if errors on the local measurements could be decreased by a factor of two, one could decisively conclude if there is tension or not. We also compare the local and CMB data assuming simple extensions of the flat, $\Lambda$CDM model (including curvature, dark energy with a constant equation of state parameter not equal to -1, non-zero neutrino masses and a non-standard number of neutrino families) and find no need for these extra parameters; in particular, we constrain the effective number of neutrino species to be Neff < 4 at 95% confidence. We show that local measurements provide constraints on the curvature and equation of state of dark energy nearly orthogonal to those of the CMB. We argue that cosmology-independent measurements of local quantities at the percent level would be very useful to explore cosmology in a model-independent way.

We have assembled a compilation of observational Hubble parameter measurements estimated with the differential evolution of cosmic chronometers, in the redshift range 0<z<1.75. This sample has been used, in combination with CMB data and with the most recent estimate of the Hubble constant H_0, to derive new constraints on several cosmological parameters. The new Hubble parameter data are very useful to break some of the parameter degeneracies present in CMB-only analysis, and to constrain possible deviations from the standard (minimal) flat \Lambda CDM model. The H(z) data are especially valuable in constraining \Omega_k and \Omega_DE in models that allow a variation of those parameters, yielding constraints that are competitive with those obtained using Supernovae and/or baryon acoustic oscillations. We also find that our H(z) data are important to constrain parameters that do no affect directly the expansion history, by breaking or reducing degeneracies with other parameters. We find that Nrel=3.45\pm0.33 using WMAP 7-years data in combination with South Pole Telescope data and our H(z) determinations (Nrel=3.71\pm0.45 using Atacama Cosmology Telescope data instead of South Pole Telescope). We exclude Nrel>4 at 95% CL (74% CL) using the same datasets combinations. We also put competitive limits on the sum of neutrino masses, \Sigma m_\nu<0.24 eV at 68% confidence level. These results have been proven to be extremely robust to many possible systematic effects, such as the initial choice of stellar population synthesis model adopted to estimate H(z) and the progenitor-bias.

We present new improved constraints on the Hubble parameter H(z) in the redshift range 0.15 < z < 1.1, obtained from the differential spectroscopic evolution of early-type galaxies as a function of redshift. We extract a large sample of early-type galaxies (\sim11000) from several spectroscopic surveys, spanning almost 8 billion years of cosmic lookback time (0.15 < z < 1.42). We select the most massive, red elliptical galaxies, passively evolving and without signature of ongoing star formation. Those galaxies can be used as standard cosmic chronometers, as firstly proposed by Jimenez & Loeb (2002), whose differential age evolution as a function of cosmic time directly probes H(z). We analyze the 4000 {\AA} break (D4000) as a function of redshift, use stellar population synthesis models to theoretically calibrate the dependence of the differential age evolution on the differential D4000, and estimate the Hubble parameter taking into account both statistical and systematical errors. We provide 8 new measurements of H(z) (see Tab. 4), and determine its change in H(z) to a precision of 5-12% mapping homogeneously the redshift range up to z \sim 1.1; for the first time, we place a constraint on H(z) at z \neq 0 with a precision comparable with the one achieved for the Hubble constant (about 5-6% at z \sim 0.2), and covered a redshift range (0.5 < z < 0.8) which is crucial to distinguish many different quintessence cosmologies. These measurements have been tested to best match a \Lambda CDM model, clearly providing a statistically robust indication that the Universe is undergoing an accelerated expansion. This method shows the potentiality to open a new avenue in constrain a variety of alternative cosmologies, especially when future surveys (e.g. Euclid) will open the possibility to extend it up to z \sim 2.

We have performed cosmo-hydro simulations using the RAMSES code to study atomic cooling (ACHs) haloes at z=10 with masses 5E7Msun<~M<~2E9Msun. We assume primordial gas and H2-cooling and prior star-formation have been suppressed. We analysed 19 haloes (gas and DM) at a resolution of ~10 (proper) pc, selected from a total volume of ~2E3 (comoving) Mpc3. This is the largest statistical hydro-sim. study of ACHs at z>10 to date. We examine the morphology, angular momentum (AM), thermodynamic, and turbulence of these haloes, in order to assess the prevalence of disks and supermassive black holes (SMBHs). We find no correlation between either the magnitude or the direction of the AM of the gas and its parent DM halo. Only 3 haloes form rotationally supported cores. Two of the most massive haloes form massive, compact overdense blobs. These blobs have an accretion rate ~0.5 Msun/yr (at a distance of 100 pc), and are possible sites of SMBH formation. Our results suggest that the degree of rotational support and the fate of the gas in a halo is determined by its large-scale environment and merger history. In particular, the two haloes forming blobs are located at knots of the cosmic web, cooled early on, and experienced many mergers. The gas in these haloes is lumpy and highly turbulent, with Mach N. >~ 5. In contrast, the haloes forming rotationally supported cores are relatively more isolated, located midway along filaments, cooled more recently, and underwent fewer mergers. Thus, the gas in these haloes is less lumpy and less turbulent (Mach <~ 4), and could retain most of its AM. The remaining 14 haloes have intermediate properties. If verified in a larger sample of haloes and with additional physics, our results will have implications for observations of the highest-redshift galaxies and quasars with JWST.

We investigate the accuracy of Eulerian perturbation theory for describing the matter and galaxy power spectra in real and redshift space in light of future observational probes for precision cosmology. Comparing the analytical results with a large suite of N-body simulations (160 independent boxes of 13.8 (Gpc/h)^3 volume each, which are publicly available), we find that re-summing terms in the standard perturbative approach predicts the real-space matter power spectrum with an accuracy of < 2% for k < 0.20 h/Mpc at redshifts z < 1.5. This is obtained following the widespread technique of writing the resummed propagator in terms of 1-loop contributions. We show that the accuracy of this scheme increases by considering higher-order terms in the resummed propagator. By combining resummed perturbation theories with several models for the mappings from real to redshift space discussed in the literature, the multipoles of the dark-matter power spectrum can be described with sub-percent deviations from N-body results for k < 0.15h/Mpc at z < 1. As a consequence, the logarithmic growth rate, f, can be recovered with sub-percent accuracy on these scales. Extending the models to massive dark-matter haloes in redshift space, our results describe the monopole term from N-body data within 2% accuracy for scales k < 0.15 h/Mpc at z < 0.5; here f can be recovered within < 5% when the halo bias is known. We conclude that these techniques are suitable to extract cosmological information from future galaxy surveys.

Observationally confirming spatial homogeneity on sufficiently large cosmological scales is of importance to test one of the underpinning assumptions of cosmology, and is also imperative for correctly interpreting dark energy. A challenging aspect of this is that homogeneity must be probed inside our past lightcone, while observations take place on the lightcone. The star formation history (SFH) in the galaxy fossil record provides a novel way to do this. We calculate the SFH of stacked Luminous Red Galaxy (LRG) spectra obtained from the Sloan Digital Sky Survey. We divide the LRG sample into 12 equal area contiguous sky patches and 10 redshift slices (0.2 < z < 0.5), which correspond to 120 blocks of volume 0.04Gpc3. Using the SFH in a time period which samples the history of the Universe between look-back times 11.5 to 13.4 Gyrs as a proxy for homogeneity, we calculate the posterior distribution for the excess large-scale variance due to inhomogeneity, and find that the most likely solution is no extra variance at all. At 95% credibility, there is no evidence of deviations larger than 5.8%.

We work out an effective theory of accelerated expansion to describe general phenomena of inflation and acceleration (dark energy) in the Universe. Our aim is to determine from theoretical grounds, in a physically-motivated and model independent way, which and how many (free) parameters are needed to broadly capture the physics of a theory describing cosmic acceleration. Our goal is to make as much as possible transparent the physical interpretation of the parameters describing the expansion. We show that, at leading order, there are five independent parameters, of which one can be constrained via general relativity tests. The other four parameters need to be determined by observing and measuring the cosmic expansion rate only, H(z). Therefore we suggest that future cosmology surveys focus on obtaining an accurate as possible measurement of $H(z)$ to constrain the nature of accelerated expansion (dark energy and/or inflation).

We present new observational constraints on inhomogeneous models based on observables independent of the CMB and large-scale structure. Using Bayesian evidence we find very strong evidence for the homogeneous LCDM model, thus disfavouring inhomogeneous models. Our new constraints are based on quantities independent of the growth of perturbations and rely on cosmic clocks based on atomic physics and on the local density of matter.

We compare the stellar population properties in the central regions of visually classified non-starforming spiral and elliptical galaxies from Galaxy Zoo and SDSS DR7. The galaxies lie in the redshift range $0.04<z<0.1$ and have stellar masses larger than $logM_*=10.4$. We select only face-on spiral galaxies in order to avoid contamination by light from the disk in the SDSS fiber and enabling the robust visual identification of spiral structure. Overall, we find that galaxies with larger central stellar velocity dispersions, regardless of morphological type, have older ages, higher metallicities, and an increased overabundance of alpha-elements. Age and alpha-enhancement, at fixed velocity dispersion, do not depend on morphological type. The only parameter that, at a given velocity dispersion, correlates with morphological type is metallicity, where the metallicity of the bulges of spiral galaxies is 0.07 dex higher than that of the ellipticals. However, for galaxies with a given total stellar mass, this dependence on morphology disappears. Under the assumption that, for our sample, the velocity dispersion traces the mass of the bulge alone, as opposed to the total mass (bulge+disk) of the galaxy, our results imply that the formation epoch of galaxy and the duration of its star-forming period are linked to the mass of the bulge. The extent to which metals are retained within the galaxy, and not removed as a result of outflows, is determined by the total mass of the galaxy.

We examine the fraction of early-type (and spiral) galaxies found in groups and clusters of galaxies as a function of dark matter halo mass. We use morphological classifications from the Galaxy Zoo project matched to halo masses from both the C4 cluster catalogue and the Yang et al (2007) group catalogue. We find that the fraction of early-type (or spiral) galaxies remains constant (changing by less than 10%) over three orders of magnitude in halo mass (13<log MH/Msol/h<15.8). This result is insensitive to our choice of halo mass measure, from velocity dispersions or summed optical luminosity. Furthermore, we consider the morphology-halo mass relations in bins of galaxy stellar mass M*, and find that while the trend of constant fraction remains unchanged, the early-type fraction amongst the most massive galaxies (11<log M*/Msol/h <12) is a factor of three greater than lower mass galaxies (10<logM*/Msol/h<10.7). We compare our observational results with those of simulations presented in De Lucia et al (2011), as well as previous observational analyses, and semi-analytic bulge (or disc) dominated galaxies from the Millennium Simulation. We find the simulations recover similar trends as observed, but may over-predict the abundances of the most massive bulge dominated (early-type) galaxies. Our results suggest that most morphological transformation is happening on the group scale before groups merge into massive clusters. However, we show that within each halo a morphology-density relation remains: it is summing the total fraction to a self-similar scaled radius which results in a flat morphology-halo mass relationship.

We identify a novel physical mechanism that may be responsible for energy release in $\gamma$-ray bursts. Radial perturbations in the neutron core, induced by its collision with collapsing outer layers during the early stages of supernova explosions, can trigger a gravitational shock, which can readily eject a small but significant fraction of the collapsing material at ultra-relativistic speeds. The development of such shocks is a strong-field effect arising in near-critical collapse in General Relativity and has been observed in numerical simulations in various contexts, including in particular radially perturbed neutron star collapse, albeit for a tiny range of initial conditions. Therefore, this effect can be easily missed in numerical simulations if the relevant parameter space is not exhaustively investigated. In the proposed picture, the observed rarity of $\gamma$-ray bursts would be explained if the relevant conditions for this mechanism appear in only about one in every $10^4-10^5$ core collapse supernovae. We also mention the possibility that near-critical collapse could play a role in powering the central engines of Active Galactic Nuclei.

We have performed cosmological N-body simulations which include the effect of the masses of the individual neutrino species. The simulations were aimed at studying the effect of different neutrino hierarchies on the matter power spectrum. Compared to the linear theory predictions, we find that non-linearities enhance the effect of hierarchy on the matter power spectrum at mildly non-linear scales. The difference between the different hierarchies is about 0.5% for a sum of neutrino masses of 0.1eV. Albeit this is a small effect, it is potentially measurable from upcoming surveys. In combination with neutrinoless double-beta decay experiments, this opens up the possibility of using the sky to determine if neutrinos are Majorana or Dirac fermions.

In this paper we present an improved fitting formula for the dark matter bispectrum motivated by the previous phenomenological approach of Scoccimarro & Couchman (2001). We use a set of LCDM simulations to calibrate the fitting parameters in the k-range of 0.03 h/Mpc<k<0.4 h/Mpc and in the redshift range of 0<z<1.5. This new proposed fit describes well the BAO-features although it was not designed to. The deviation between the simulations output and our analytic prediction is typically less than 5% and in the worst case is never above 10%. We envision that this new analytic fitting formula will be very useful in providing reliable predictions for the non-linear dark matter bispectrum for LCDM models.

We discuss the phenomenology of the minimal inflation scenario. We concentrate on two aspects: inflationary trajectories and particle production. Our findings can be summarized in two main results: first, that inflationary trayectories that are very flat and provide enough number of e-foldings are natural in the scenario without fine tunning. We present a general formalism to identify attractors in multi-field inflation regardless of trajectories fulfilling the slow-roll conditions. We then explore particle production in the model and show how the inflaton naturally transmutes into a dark matter particle. One interesting feature of our model is that it provides a novel mechanism to generate particles and entropy in the universe: the filling of the Fermi sphere up to a given momentum p_F due to the sea of goldstinos that are an important part of the matter generated after inflation. With this mechanism in hand we predict that the gravitino should have a mass of > 100-1000 TeV. Another interesting feature of our model is that the predicted level of gravity waves is r=0.1 - 0.001, which is in the range of detectability from Planck and upcoming CMB polarization experiments.

Conversion of photons into axions under the presence of a strong magnetic field can dim the radiation from magnetized astrophysical objects. Here we perform a detailed calculation aimed at quantifying the signatures of photon-axion conversion in the spectra, light curves, and polarization of neutron stars (NSs). We take into account the energy and angle-dependence of the conversion probability and the surface thermal emission from NSs. The latter is computed from magnetized atmosphere models that include the effect of photon polarization mode conversion due to vacuum polarization. The resulting spectral models, inclusive of the general-relativistic effects of gravitational redshift and light deflection, allow us to make realistic predictions for the effects of photon to axion conversion on observed NS spectra, light curves, and polarization signals. We identify unique signatures of the conversion, such as an increase of the effective area of a hot spot as it rotates away from the observer line of sight. For a star emitting from the entire surface, the conversion produces apparent radii that are either larger or smaller (depending on axion mass and coupling strength) than the limits set by NS equations of state. For an emission region that is observed phase-on, photon-axion conversion results in an inversion of the plane of polarization with respect to the no-conversion case. While the quantitative details of the features that we identify depend on NS properties (magnetic field strength, temperature) and axion parameters, the spectral and polarization signatures induced by photon-axion conversion are distinctive enough to make NSs very interesting and promising probes of axion physics.

Using observational data on the expansion rate of the universe (H(z)) we constrain the effective Lagrangian of the current accelerated expansion. Our results show that the effective potential is consistent with being flat i.e., a cosmological constant; it is also consistent with the field moving along an almost flat potential like a pseudo-Goldstone boson. We show that the potential of dark energy does not deviate from a constant at more than 6% over the redshift range 0 < z < 1. The data can be described by just a constant term in the Lagrangian and do not require any extra parameters; therefore there is no evidence for augmenting the number of parameters of the LCDM paradigm. We also find that the data justify the effective theory approach to describe accelerated expansion and that the allowed parameters range satisfy the expected hierarchy. Future data, both from cosmic chronometers and baryonic acoustic oscillations, that can measure H(z) at the % level, could greatly improve constraints on the flatness of the potential or shed some light on possible mechanisms driving the accelerated expansion. Besides the above result, it is shown that the effective Lagrangian of accelerated expansion can be constrained from cosmological observations in a model-independent way and that direct measurements of the expansion rate H(z) are most useful to do so.

We critically investigate current statistical tests applied to high redshift clusters of galaxies in order to test the standard cosmological model and describe their range of validity. We carefully compare a sample of high-redshift, massive, galaxy clusters with realistic Poisson sample simulations of the theoretical mass function, which include the effect of Eddington bias. We compare the observations and simulations using the following statistical tests: the distributions of ensemble and individual existence probabilities (in the >M,>z sense), the redshift distributions, and the 2d Kolmogorov-Smirnov test. Using seemingly rare clusters from Hoyle et al. (2011), and Jee et al. (2011) and assuming the same survey geometry as in Jee et al. (2011, which is less conservative than Hoyle et al. 2011), we find that the (>M,>z) existence probabilities of all clusters are fully consistent with LCDM. However assuming the same survey geometry, we use the 2d K-S test probability to show that the observed clusters are not consistent with being the least probable clusters from simulations at >95% confidence, and are also not consistent with being a random selection of clusters, which may be caused by the non-trivial selection function and survey geometry. Tension can be removed if we examine only a X-ray selected sub sample, with simulations performed assuming a modified survey geometry.

In this paper we analyze a suite of cosmological simulations of modified gravitational action f(R) models, where cosmic acceleration is induced by a scalar field that acts as a fifth force on all forms of matter. In particular, we focus on the bispectrum of the dark matter density field on mildly non-linear scales. For models with the same initial power spectrum, the dark matter bispectrum shows significant differences for cases where the final dark matter power spectrum also differs. Given the different dependence on bias of the galaxy power spectrum and bispectrum, bispectrum measurements can close the loophole of galaxy bias hiding differences in the power spectrum. Alternatively, changes in the initial power spectrum can also hide differences. By constructing LCDM models with very similar final non-linear power spectra, we show that the differences in the bispectrum are reduced (<4%) and are comparable with differences in the imperfectly matched power spectra. These results indicate that the bispectrum depends mainly on the power spectrum and less sensitively on the gravitational signatures of the f(R) model. This weak dependence of the matter bispectrum on gravity makes it useful for breaking degeneracies associated with galaxy bias, even for models beyond general relativity.

We elaborate on a minimal inflation scenario based entirely on the general properties of supersymmetry breaking in supergravity models. We identify the inflaton as the scalar component of the Goldstino superfield. We write plausible candidates for the effective action describing this chiral superfield. In particular the theory depends (apart from parameters of O(1)) on a single free parameter: the scale of supersymmetry breaking. This can be fixed using the amplitude of CMB cosmological perturbations and we therefore obtain the scale of supersymmetry breaking to be 10^{12-14} GeV. The model also incorporates explicit R-symmetry breaking in order to satisfy the slow roll conditions. In our model the eta-problem is solved without extra fine-tuning. We try to obtain as much information as possible in a model independent way using general symmetry properties of the theory's effective action, this leads to a new proposal on how to exit the inflationary phase and reheat the Universe.

We have performed a cosmological numerical simulation of primordial baryonic gas collapsing onto a $3\times10^7$M$_{\odot}$ dark matter (DM) halo. We show that the large scale baryonic accretion process and the merger of few $\sim10^6$ M$_{\odot}$ DM halos, triggered by the gravitational potential of the biggest halo, is enough to create super sonic ($\mathcal{M}>10$) shocks and develop a turbulent environment. In this scenario the post shocked regions are able to produced both H$_2$ and HD molecules very efficiently reaching maximum abundances of $n_\mathrm{H_2}\sim10^{-2}n_\mathrm{H}$ and $n_\mathrm{HD}\sim \mathrm{few}\times10^{-6}n_\mathrm{H}$, enough to cool the gas below 100K in some regions. The kinetic energy spectrum of the turbulent primordial gas is close to a Burgers spectrum, $\hat{E}_k\propto k^{-2}$, which could favor the formation of low mass primordial stars. The solenoidal to total kinetic energy ratio is $0.65\la R_k\la0.7$ for a wide range of wave numbers; this value is close to $R_k\approx 2/3$ natural equipartition energy value of a random turbulent flow. In this way turbulence and molecular cooling seem to work together in order to produce potential star formation regions of cold and dense gas in primordial environments. We conclude that both the mergers and the collapse process onto the main DM halo provide enough energy to develop super sonic turbulence which favor the molecular coolants formation: this mechanism, which could be universal and the main route toward formation of the first galaxies, is able to create potential star forming regions at high redshift.

The standard Friedmann model of cosmology is based on the Copernican Principle, i.e. the assumption of a homogeneous background on which structure forms via perturbations. Homogeneity underpins both general relativistic and modified gravity models and is central to the way in which we interpret observations of the CMB and the galaxy distribution. It is therefore important to probe homogeneity via observations. We describe a test based on the fossil record of distant galaxies: if we can reconstruct key intrinsic properties of galaxies as functions of proper time along their worldlines, we can compare such properties at the same proper time for our galaxy and others. We achieve this by computing the lookback time using radial Baryon Acoustic Oscillations, and the time along galaxy world line using stellar physics, allowing us to probe homogeneity, in principle anywhere inside the past light cone. Agreement in the results would be an important consistency test -- although it would not in itself prove homogeneity. Any significant deviation in the results however would signal a breakdown of homogeneity.

We present a method to minimize, or even cancel out, the nuisance parameters affecting a measurement. Our approach is general and can be applied to any experiment or observation. We compare it with the bayesian technique used to deal with nuisance parameters: marginalization, and show how the method compares and improves by avoiding biases. We illustrate the method with several examples taken from the astrophysics and cosmology world: baryonic acoustic oscillations, cosmic clocks, Supernova Type Ia luminosity distance, neutrino oscillations and dark matter detection. By applying the method we recover some known results but also find some interesting new ones. For baryonic acoustic oscillation (BAO) experiments we show how to combine radial and angular BAO measurements in order to completely eliminate the dependence on the sound horizon at radiation drag. In the case of exploiting SN1a as standard candles we show how the uncertainty in the luminosity distance by a second parameter modeled as a metallicity dependence can be eliminated or greatly reduced. When using cosmic clocks to measure the expansion rate of the universe, we demonstrate how a particular combination of observables nearly removes the metallicity dependence of the galaxy on determining differential ages, thus removing the age-metallicity degeneracy in stellar populations. We hope that these findings will be useful in future surveys to obtain robust constraints on the dark energy equation of state.

To date, 14 high-redshift (z>1.0) galaxy clusters with mass measurements have been observed, spectroscopically confirmed and are reported in the literature. These objects should be exceedingly rare in the standard LCDM model. We conservatively approximate the selection functions of these clusters' parent surveys, and quantify the tension between the abundances of massive clusters as predicted by the standard LCDM model and the observed ones. We alleviate the tension considering non-Gaussian primordial perturbations of the local type, characterized by the parameter fnl and derive constraints on fnl arising from the mere existence of these clusters. At the 95% confidence level, fnl>467 with cosmological parameters fixed to their most likely WMAP5 values, or fnl > 123 (at 95% confidence) if we marginalize over WMAP5 parameters priors. In combination with fnl constraints from Cosmic Microwave Background and halo bias, this determination implies a scale-dependence of fnl at approx. 3 sigma. Given the assumptions made in the analysis, we expect any future improvements to the modeling of the non-Gaussian mass function, survey volumes, or selection functions to increase the significance of fnl>0 found here. In order to reconcile these massive, high-z clusters with an fnl=0, their masses would need to be systematically lowered by 1.5 sigma or the sigma8 parameter should be approx. 3 sigma higher than CMB (and large-scale structure) constraints. The existence of these objects is a puzzle: it either represents a challenge to the LCDM paradigme or it is an indication that the mass estimates of clusters is dramatically more uncertain than we think.

We present a new methodology to determine the expansion history of the Universe analyzing the spectral properties of early type galaxies (ETG). We found that for these galaxies the 4000\AA break is a spectral feature that correlates with the relative ages of ETGs. In this paper we describe the method, explore its robustness using theoretical synthetic stellar population models, and apply it using a SDSS sample of $\sim$14 000 ETGs. Our motivation to look for a new technique has been to minimise the dependence of the cosmic chronometer method on systematic errors. In particular, as a test of our method, we derive the value of the Hubble constant $H_0 = 72.6 \pm 2.8$ (stat) $\pm2.3$ (syst) (68% confidence), which is not only fully compatible with the value derived from the Hubble key project, but also with a comparable error budget. Using the SDSS, we also derive, assuming w=constant, a value for the dark energy equation of state parameter $w = -1 \pm 0.2$ (stat) $\pm0.3$ (syst). Given the fact that the SDSS ETG sample only reaches $z \sim 0.3$, this result shows the potential of the method. In future papers we will present results using the high-redshift universe, to yield a determination of H(z) up to $z \sim 1$.

We explore how the Peccei-Quinn (PQ) axion parameter space can be constrained by the frequency-dependence dimming of radiation from astrophysical objects. To do so we perform accurate calculations of photon-axion conversion in the presence of a variable magnetic field. We propose several tests where the PQ axion parameter space can be explored with current and future astronomical surveys: the observed spectra of isolated neutron stars, occultations of background objects by white dwarfs and neutron stars, the light curves of eclipsing binaries containing a white dwarf. We find that the lack of dimming of the light curve of a detached eclipsing white dwarf binary recently observed, leads to relevant constraints on the photon-axion conversion. Current surveys designed for Earth-like planet searches are well matched to strengthen and improve the constraints on the PQ axion using astrophysical objects radiation dimming.

We present a modification of the standard halo model with the goal of providing an improved description of galaxy clustering. Recent surveys, like the Sloan Digital Sky Survey (SDSS) and the Anglo-Australian Two-degree survey (2dF), have shown that there seems to be a correlation between the clustering of galaxies and their properties such as metallicity and star formation rate, which are believed to be environment-dependent. This environmental dependence is not included in the standard halo model where the host halo mass is the only variable specifying galaxy properties. In our approach, the halo properties i.e., the concentration, and the Halo Occupation Distribution --HOD-- prescription, will not only depend on the halo mass (like in the standard halo model) but also on the halo environment. We examine how different environmental dependence of halo concentration and HOD prescription affect the correlation function. We see that at the level of dark matter, the concentration of haloes affects moderately the dark matter correlation function only at small scales. However the galaxy correlation function is extremely sensitive to the HOD details, even when only the HOD of a small fraction of haloes is modified.

Using the publicly available VESPA database of SDSS Data Release 7 spectra, we calculate the stellar Mass Weighted Age (hereafter MWA) as a function of local galaxy density and dark matter halo mass. We compare our results with semi-analytic models from the public Millennium Simulation. We find that the stellar MWA has a large scatter which is inherent in the data and consistent with that seen in semi-analytic models. The stellar MWA is consistent with being independent (to first order) with local galaxy density, which is also seen in semi-analytic models. As a function of increasing dark matter halo mass (using the SDSS New York Value Added Group catalogues), we find that the average stellar MWA for member galaxies increases, which is again found in semi-analytic models. Furthermore we use public dark matter Mass Accretion History (MAH) code calibrated on simulations, to calculate the dark matter Mass Weighted Age as a function of dark matter halo mass. In agreement with earlier analyses, we find that the stellar MWA and the dark matter MWA are anti correlated for large mass halos, i.e, dark matter accretion does not seem to be the primary factor in determining when stellar mass was compiled. This effect can be described by down-sizing.

We report results from a cosmological simulation with non-equilibrium chemistry of 21 species, including H2, HD, and LiH molecular cooling. Starting from cosmological initial conditions, we focus on the evolution of the central 1.8 Kpc region of a 3 x 10^7 Msun halo. The crossing of a few 10^6 Msun halos and the gas accretion through larger scale filaments generate a turbulent environment within this region. Due to the short cooling time caused by the non-equilibrium formation of H2, the supersonic turbulence results in a very fragmented mass distribution, where dense, gravitationally unstable clumps emerge from a complex network of dense filaments. At z=10.87, we find approximately 25 well defined, gravitationally unstable clumps, with masses of 4 x 10^3-9 x 10^5 Msun, temperatures of approximately 300K, and cooling times much shorter than the free-fall time. Only the initial phase of the collapse of individual clumps is spatially resolved in the simulation. Depending on the density reached in the collapse, the estimated average Bonnor-Ebert masses are in the range 200-800 Msun. We speculate that each clump may further fragment into a cluster of stars with a characteristic mass in the neighborhood of 50 Msun. This process at z ~ 11 may represent the dominant mode of Pop. III star formation, causing a rapid chemical enrichment of the protogalactic environment.

We illustrate how recently improved low-redshift cosmological measurements can tighten constraints on neutrino properties. In particular we examine the impact of the assumed cosmological model on the constraints. We first consider the new HST H0 = 74.2 +/- 3.6 measurement by Riess et al. (2009) and the sigma8*(Omegam/0.25)^0.41 = 0.832 +/- 0.033 constraint from Rozo et al. (2009) derived from the SDSS maxBCG Cluster Catalog. In a Lambda CDM model and when combined with WMAP5 constraints, these low-redshift measurements constrain sum mnu<0.4 eV at the 95% confidence level. This bound does not relax when allowing for the running of the spectral index or for primordial tensor perturbations. When adding also Supernovae and BAO constraints, we obtain a 95% upper limit of sum mnu<0.3 eV. We test the sensitivity of the neutrino mass constraint to the assumed expansion history by both allowing a dark energy equation of state parameter w to vary, and by studying a model with coupling between dark energy and dark matter, which allows for variation in w, Omegak, and dark coupling strength xi. When combining CMB, H0, and the SDSS LRG halo power spectrum from Reid et al. 2009, we find that in this very general model, sum mnu < 0.51 eV with 95% confidence. If we allow the number of relativistic species Nrel to vary in a Lambda CDM model with sum mnu = 0, we find Nrel = 3.76^{+0.63}_{-0.68} (^{+1.38}_{-1.21}) for the 68% and 95% confidence intervals. We also report prior-independent constraints, which are in excellent agreement with the Bayesian constraints.

Cosmological analysis based on currently available observations are unable to rule out a sizeable coupling among the dark energy and dark matter fluids. We explore a variety of coupled dark matter-dark energy models, which satisfy cosmic microwave background constraints, in light of low redshift and near universe observations. We illustrate the phenomenology of different classes of dark coupling models, paying particular attention in distinguishing between effects that appear only on the expansion history and those that appear in the growth of structure. We find that while a broad class of dark coupling models are effectively models where general relativity (GR) is modified --and thus can be probed by a combination of tests for the expansion history and the growth of structure--, there is a class of dark coupling models where gravity is still GR, but the growth of perturbations is, in principle modified. While this effect is small in the specific models we have considered, one should bear in mind that an inconsistency between reconstructed expansion history and growth may not uniquely indicate deviations from GR. Our low redshift constraints arise from cosmic velocities, redshift space distortions and dark matter abundance in galaxy voids. We find that current data constrain the dimensionless coupling to be |xi|<0.2, but prospects from forthcoming data are for a significant improvement. Future, precise measurements of the Hubble constant, combined with high-precision constraints on the growth of structure, could provide the key to rule out dark coupling models which survive other tests. We shall exploit as well weak equivalence principle violation arguments, which have the potential to highly disfavour a broad family of coupled models.

Comparing clustering of differently biased tracers of the dark matter distribution offers the opportunity to reduce the cosmic variance error in the measurement of certain cosmological parameters. We develop a formalism that includes bias non-linearities and stochasticity. Our formalism is general enough that can be used to optimise survey design and tracers selection and optimally split (or combine) tracers to minimise the error on the cosmologically interesting quantities. Our approach generalises the one presented by McDonald & Seljak (2009) of circumventing sample variance in the measurement of $f\equiv d \ln D/d\ln a$. We analyse how the bias, the noise, the non-linearity and stochasticity affect the measurements of $Df$ and explore in which signal-to-noise regime it is significantly advantageous to split a galaxy sample in two differently-biased tracers. We use N-body simulations to find realistic values for the parameters describing the bias properties of dark matter haloes of different masses and their number density. We find that, even if dark matter haloes could be used as tracers and selected in an idealised way, for realistic haloes, the sample variance limit can be reduced only by up to a factor $\sigma_{2tr}/\sigma_{1tr}\simeq 0.6$. This would still correspond to the gain from a three times larger survey volume if the two tracers were not to be split. Before any practical application one should bear in mind that these findings apply to dark matter haloes as tracers, while realistic surveys would select galaxies: the galaxy-host halo relation is likely to introduce extra stochasticity, which may reduce the gain further.

The Atacama Cosmology Telescope (ACT) is currently observing the cosmic microwave background with arcminute resolution at 148 GHz, 218 GHz, and 277 GHz. In this paper, we present ACT's first results. Data have been analyzed using a maximum-likelihood map-making method which uses B-splines to model and remove the atmospheric signal. It has been used to make high-precision beam maps from which we determine the experiment's window functions. This beam information directly impacts all subsequent analyses of the data. We also used the method to map a sample of galaxy clusters via the Sunyaev-Zel'dovich (SZ) effect, and show five clusters previously detected with X-ray or SZ observations. We provide integrated Compton-y measurements for each cluster. Of particular interest is our detection of the z = 0.44 component of A3128 and our current non-detection of the low-redshift part, providing strong evidence that the further cluster is more massive as suggested by X-ray measurements. This is a compelling example of the redshift-independent mass selection of the SZ effect.

We present a series of colour evolution models for Luminous Red Galaxies (LRGs) in the 7th spectroscopic data release of the Sloan Digital Sky Survey (SDSS), computed using the full-spectrum fitting code VESPA on high signal-to-noise stacked spectra. The colour-evolution models are computed as a function of colour, luminosity and redshift, and we do not a-priori assume that LRGs constitute a uniform population of galaxies in terms of stellar evolution. By computing star-formation histories from the fossil record, the measured stellar evolution of the galaxies is decoupled from the survey's selection function, which also evolves with redshift. We present these evolutionary models computed using three different sets of Stellar Population Synthesis (SPS) codes. We show that the traditional fiducial model of purely passive stellar evolution of LRGs is broadly correct, but it is not sufficient to explain the full spectral signature. We also find that higher-order corrections to this model are dependent on the SPS used, particularly when calculating the amount of recent star formation. The amount of young stars can be non-negligible in some cases, and has important implications for the interpretation of the number density of LRGs within the selection box as a function of redshift. Dust extinction, however, is more robust to the SPS modelling: extinction increases with decreasing luminosity, increasing redshift, and increasing r-i colour. We are making the colour evolution tracks publicly available at http://www.icg.port.ac.uk/~tojeiror/lrg_evolution/.

We present a catalog of 105 rich and massive ($M>3\times10^{14}M_{\sun}$) optically-selected clusters of galaxies extracted from 70 square-degrees of public archival griz imaging from the Blanco 4-m telescope acquired over 45 nights between 2005 and 2007. We use the clusters' optically-derived properties to estimate photometric redshifts, optical luminosities, richness, and masses. We complement the optical measurements with archival XMM-Newton and ROSAT X-ray data which provide additional luminosity and mass constraints on a modest fraction of the cluster sample. Two of our clusters show clear evidence for central lensing arcs; one of these has a spectacular large-diameter, nearly-complete Einstein Ring surrounding the brightest cluster galaxy. A strong motivation for this study is to identify the massive clusters that are expected to display prominent signals from the Sunyaev-Zeldovich Effect (SZE) and therefore be detected in the wide-area mm-band surveys being conducted by both the Atacama Cosmology Telescope and the South Pole Telescope. The optical sample presented here will be useful for verifying new SZE cluster candidates from these surveys, for testing the cluster selection function, and for stacking analyzes of the SZE data.

We present a measurement of the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz. The measurement uses maps with 1.4' angular resolution made with data from the Atacama Cosmology Telescope (ACT). The observations cover 228 square degrees of the southern sky, in a 4.2-degree-wide strip centered on declination 53 degrees South. The CMB at arcminute angular scales is particularly sensitive to the Silk damping scale, to the Sunyaev-Zel'dovich (SZ) effect from galaxy clusters, and to emission by radio sources and dusty galaxies. After masking the 108 brightest point sources in our maps, we estimate the power spectrum between 600 < \ell < 8000 using the adaptive multi-taper method to minimize spectral leakage and maximize use of the full data set. Our absolute calibration is based on observations of Uranus. To verify the calibration and test the fidelity of our map at large angular scales, we cross-correlate the ACT map to the WMAP map and recover the WMAP power spectrum from 250 < ell < 1150. The power beyond the Silk damping tail of the CMB is consistent with models of the emission from point sources. We quantify the contribution of SZ clusters to the power spectrum by fitting to a model normalized at sigma8 = 0.8. We constrain the model's amplitude ASZ < 1.63 (95% CL). If interpreted as a measurement of sigma8, this implies sigma8^SZ < 0.86 (95% CL) given our SZ model. A fit of ACT and WMAP five-year data jointly to a 6-parameter LCDM model plus terms for point sources and the SZ effect is consistent with these results.

Using light curves and host galaxy spectra of 101 Type Ia supernovae (SNe Ia) with redshift $z \lesssim 0.3$ from the SDSS Supernova Survey (SDSS-SN), we derive the SN Ia rate as a function of progenitor age (the delay time distribution, or DTD). We use the VESPA stellar population synthesis algorithm to analyze the SDSS spectra of all galaxies in the field searched by SDSS-SN, giving us a reference sample of 77,000 galaxies for our SN Ia hosts. Our method does not assume any a priori shape for the DTD and therefore is minimally parametric. We present the DTD in physical units for high stretch (luminous, slow declining) and low stretch (subluminous, fast declining) supernovae in three progenitor age bins. We find strong evidence of two progenitor channels: one that produces high stretch SNe Ia $\lesssim 400$ Myr after the birth of the progenitor system, and one that produces low stretch SNe Ia with a delay $\gtrsim 2.4$ Gyr. We find that each channel contributes roughly half of the Type Ia rate in our reference sample. We also construct the average spectra of high stretch and low stretch SN Ia host galaxies, and find that the difference of these spectra looks like a main sequence B star with nebular emission lines indicative of star formation. This supports our finding that there are two populations of SNe Ia, and indicates that the progenitors of high stretch SNe are at the least associated with very recent star formation in the last few tens of Myr. Our results provide valuable constraints for models of Type Ia progenitors and may help improve the calibration of SNe Ia as standard candles.

We update constraints on cosmic opacity by combining recent SN Type Ia data compilation with the latest measurements of the Hubble expansion at redshifts between 0 and 2. The new constraint on the parameter $\epsilon$ parametrising deviations from the luminosity-angular diameter distance relation ($d_L=d_A(1+z)^{2+\epsilon}$), is $\epsilon=-0.04_{-0.07}^{+0.08}$ (2-$\sigma$). For the redshift range between 0.2 and 0.35 this corresponds to an opacity $\Delta\tau<0.012$ (95% C.L.), a factor of 2 stronger than the previous constraint. Various models of beyond the standard model physics that predict violation of photon number conservation contribute to the opacity and can be equally constrained. In this paper we put new limits on axion-like particles, including chameleons, and mini-charged particles.

Cosmological probes are steadily reducing the total neutrino mass window, resulting in constraints on the neutrino-mass degeneracy as the most significant outcome. In this work we explore the discovery potential of cosmological probes to constrain the neutrino hierarchy, and point out some subtleties that could yield spurious claims of detection. This has an important implication for next generation of double beta decay experiments, that will be able to achieve a positive signal in the case of degenerate or inverted hierarchy of Majorana neutrinos. We find that cosmological experiments that nearly cover the whole sky could in principle distinguish the neutrino hierarchy by yielding 'substantial' evidence for one scenario over the another, via precise measurements of the shape of the matter power spectrum from large scale structure and weak gravitational lensing.

We report on the finding of a correlation between galaxies' past star formation activity and the degree to which neighbouring galaxies rotation axes are aligned. This is obtained by cross-correlating star formation histories, derived with MOPED, and spin direction (chirality), as determined by the Galaxy Zoo project, for a sample of SDSS galaxies. Our findings suggest that spiral galaxies which formed the majority of their stars early (z > 2) tend to display coherent rotation over scales of ~10 Mpc/h. The correlation is weaker for galaxies with significant recent star formation. We find evidence for this alignment at more than the 5-sigma level, but no correlation with other galaxy stellar properties. This finding can be explained within the context of hierarchical tidal-torque theory if the SDSS galaxies harboring the majority of the old stellar population where formed in the past, in the same filament and at about the same time. Galaxies with significant recent star formation instead are in the field, thus influenced by the general tidal field that will align them in random directions or had a recent merger which would promote star formation, but deviate the spin direction.

Using the universal X superfield that measures in the UV the violation of conformal invariance we build up a model of multifield inflation. The underlying dynamics is the one controlling the natural flow of this field in the IR to the Goldstino superfield once SUSY is broken. We show that flat directions satisfying the slow roll conditions exist only if R-symmetry is broken. Naturalness of our model leads to scales of SUSY breaking of the order of 10^{11-13} Gev, a nearly scale-invariant spectrum of the initial perturbations and negligible gravitational waves. We obtain that the inflaton field is lighter than the gravitino by an amount determined by the slow roll parameter eta. The existence of slow-roll conditions is directly linked to the values of supersymmetry and R-symmetry breaking scales. We make cosmological predictions of our model and compare them to current data.

We applied the VESPA algorithm to the Sloan Digital Sky Survey final data release of the Main Galaxies and Luminous Red Galaxies samples. The result is a catalogue of stellar masses, detailed star formation and metallicity histories and dust content of nearly 800,000 galaxies. We make the catalogue public via a T-SQL database, which is described in detail in this paper. We present the results using a range of stellar population and dust models, and will continue to update the catalogue as new and improved models are made public. The data and documentation are currently online, and can be found at http://www-wfau.roe.ac.uk/vespa/. We also present a brief exploration of the catalogue, and show that the quantities derived are robust: luminous red galaxies can be described by one to three populations, whereas a main galaxy sample galaxy needs on average two to five; red galaxies are older and less dusty; the dust values we recover are well correlated with measured Balmer decrements and star formation rates are also in agreement with previous measurements.

The recent weak lensing measurement of the dark matter mass of the high-redshift galaxy cluster XMMUJ2235.3-2557 of (8.5 +- 1.7) x 10^{14} Msun at z=1.4, indicates that, if the cluster is assumed to be the result of the collapse of dark matter in a primordial gaussian field in the standard LCDM model, then its abundance should be < 0.002 clusters in the observed area. Here we investigate how to boost the probability of XMMUJ2235.3-2557 in particular resorting to deviations from Gaussian initial conditions. We show that this abundance can be boosted by factors > 3-10 if the non-Gaussianity parameter f^local_NL is in the range 150-200. This value is comparable to the limit for f_NL obtained by current constraints from the CMB. We conclude that mass determination of high-redshift, massive clusters can offer a complementary probe of primordial non-gaussianity.

We present the first lensing mass measurements of Sunyaev-Zel'dovich (SZ) selected clusters. Using optical imaging from the Southern Cosmology Survey (SCS), we present weak lensing masses for three clusters selected by their SZ emission in the South Pole Telescope survey (SPT). We confirm that the SZ selection procedure is successful in detecting mass concentrations. We also study the weak lensing signals from 38 optically-selected clusters in ~8 square degrees of the SCS survey. We fit Navarro, Frenk and White (NFW) profiles and find that the SZ clusters have amongst the largest masses, as high as 5x10^14 Msun. Using the best fit masses for all the clusters, we analytically calculate the expected SZ integrated Y parameter, which we find to be consistent with the SPT observations.

We present new determinations of the cosmic expansion history from red-envelope galaxies. We have obtained for this purpose high-quality spectra with the Keck-LRIS spectrograph of red-envelope galaxies in 24 galaxy clusters in the redshift range 0.2 < z < 1.0. We complement these Keck spectra with high-quality, publicly available archival spectra from the SPICES and VVDS surveys. We improve over our previous expansion history measurements in Simon et al. (2005) by providing two new determinations of the expansion history: H(z) = 97 +- 62 km/sec/Mpc at z = 0.5 and H(z) = 90 +- 40 km/sec/Mpc at z = 0.8. We discuss the uncertainty in the expansion history determination that arises from uncertainties in the synthetic stellar-population models. We then use these new measurements in concert with cosmic-microwave-background (CMB) measurements to constrain cosmological parameters, with a special emphasis on dark-energy parameters and constraints to the curvature. In particular, we demonstrate the usefulness of direct H(z) measurements by constraining the dark- energy equation of state parameterized by w0 and wa and allowing for arbitrary curvature. Further, we also constrain, using only CMB and H(z) data, the number of relativistic degrees of freedom to be 4 +- 0.5 and their total mass to be < 0.2 eV, both at 1-sigma.

We present a spectroscopic catalog of (mostly) red galaxies in 24 galaxy clusters in the redshift range 0.17 < z < 0.92 obtained with the LRIS instrument on the Keck I telescope. Here we describe the observations and the galaxy spectra, including the discovery of three cD galaxies with LINER emission spectra, and the spectroscopic discovery of four new galaxy-galaxy lenses in cluster environments.

We explore consistency among different distance measures, including Supernovae Type Ia data, measurements of the Hubble parameter, and determination of the Baryon acoustic oscillation scale. We present new constraints on the cosmic transparency combining $H(z)$ data together with the latest Supernova Type Ia data compilation. This combination, in the context of a flat $\Lambda$CDM model, improves current constraints by nearly an order of magnitude. We re-examine the recently reported tension between the Baryon acoustic oscillation scale and Supernovae data in light of possible deviations from transparency, concluding that the source of the discrepancy may most likely be found among systematic effects of the modelling of the low redshift data or a simple $\sim 2-\sigma$ statistical fluke, rather than in exotic physics. Finally, we attempt to draw model-independent conclusions about the recent accelerated expansion, determining the acceleration redshift to be $z_{acc}=0.35^{+0.20}_{-0.13}$ (1-$\sigma$).

We propose a novel mechanism for powering the central engines of Active Galactic Nuclei through super-critical (type II) black hole collapse. In this picture, ~$10^3 M_\odot$ of material collapsing at relativistic speeds can trigger a gravitational shock, which can eject a large percentage of the collapsing matter at relativistic speeds, leaving behind a "light" black hole. In the presence of a poloidal magnetic field, the plasma collimates along two jets, and the associated electron synchrotron radiation can easily account for the observed radio luminosities, sizes and durations of AGN jets. For Lorentz factors of order 100 and magnetic fields of a few hundred $\mu G$, synchrotron electrons can shine for $10^6$ yrs, producing jets of sizes of order 100 kpc. This mechanism may also be relevant for Gamma Ray Bursts and, in the absence of magnetic field, supernova explosions.

We present first results from the Southern Cosmology Survey, a new multiwavelength survey of the southern sky coordinated with the Atacama Cosmology Telescope (ACT), a recently commissioned ground-based mm-band Cosmic Microwave Background experiment. This article presents a full analysis of archival optical multi-band imaging data covering an 8 square degree region near right ascension 23 hours and declination -55 degrees, obtained by the Blanco 4-m telescope and Mosaic-II camera in late 2005. We describe the pipeline we have developed to process this large data volume, obtain accurate photometric redshifts, and detect optical clusters. Our cluster finding process uses the combination of a matched spatial filter, photometric redshift probability distributions and richness estimation. We present photometric redshifts, richness estimates, luminosities, and masses for 8 new optically-selected clusters with mass greater than $3\times10^{14}M_{\sun}$ at redshifts out to 0.7. We also present estimates for the expected Sunyaev-Zel'dovich effect (SZE) signal from these clusters as specific predictions for upcoming observations by ACT, the South Pole Telescope and Atacama Pathfinder Experiment.

We use a non-equilibrium chemical network to revisit and study the effect of H_{2}, HD and LiH molecular cooling on a primordial element of gas. We solve both the thermal and chemical equations for a gas element with an initial temperature T\approx 1000K and a gas number density in the range n_{tot}=1-10^{4} cm^{-3}. At low densities, n_{tot}<10^{2} cm^{-3}, the gas reaches temperatures \sim 100K and the main coolant is H_{2}, but at higher densities, n_{tot}>10^{2} cm^{-3}, the HD molecule dominates the gas temperature evolution. The effect of LiH is negligible in all cases. We studied the effect of D abundance on the gas cooling. The D abundance was set initially to be in the range n_{D}/n_{H}=10^{-7}-10^{-4.5}, with n_{HD}/n_{H}={D^{+}}/n_{H}=10^{-10}. The simulations show that at n_{tot}>10^{2} cm^{-3} the HD cooling dominates the temperature evolution for D abundances greater than 10^{-5}n_{H}. This number decrease at higher densities. Furthermore, we studied the effect of electrons and ionized particules on the gas temperature. We followed the gas temperature evolution with n_{H_{+}}/n_{H}=10^{-4}-10^{-1} and n_{D^{+}}/n_{H^{+}}=10^{-5}. The gas temperature reached lower values at high ionization degree because electrons, H^{+} and D^{+} are catalizers in the formation paths of the H_{2} and HD molecules, which are the main coolers at low temperatures. Finaly, we studied the effect of an OB star, with T_{eff}=4\times 10^{4}K, would have on gas cooling. It is very difficult for a gas with n_{tot} in the range between 1-100 cm^{-3} to drop its temperature if the star is at a distance less than 100 pc.

Baryon Acoustic Oscillations (BAO) in the radial direction offer a method to directly measure the Universe expansion history, and to set limits to space curvature when combined to the angular BAO signal. In addition to spectroscopic surveys, radial BAO might be measured from accurate enough photometric redshifts obtained with narrow-band filters. We explore the requirements for a photometric survey using Luminous Red Galaxies (LRG) to competitively measure the radial BAO signal and discuss the possible systematic errors of this approach. If LRG were a highly homogeneous population, we show that the photo-z accuracy would not substantially improve by increasing the number of filters beyond $\sim 10$, except for a small fraction of the sources detected at high signal-to-noise, and broad-band filters would suffice to achieve the target $\sigma_z = 0.003 (1+z)$ for measuring radial BAO. Using the LRG spectra obtained from SDSS, we find that the spectral variability of LRG substantially worsens the achievable photometric redshift errors, and that the optimal system consists of $\sim$ 30 filters of width $\Delta \lambda / \lambda \sim 0.02$. A $S/N > 20$ is generally necessary at the filters on the red side of the $H\alpha$ break to reach the target photometric accuracy. We estimate that a 5-year survey in a dedicated telescope with etendue in excess of 60 ${\rm m}^2 {\rm deg}^2$ would be necessary to obtain a high enough density of galaxies to measure radial BAO with sufficiently low shot noise up to $z= 0.85$. We conclude that spectroscopic surveys have a superior performance than photometric ones for measuring BAO in the radial direction.

We present catalogs of QSO candidates selected using photometry from GALEX combined with SDSS in the Stripe 82 region and Blanco Cosmology Survey (BCS) near declination -55 degrees. The SDSS region contains ~700 objects with magnitude i < 20 and ~3600 objects with i < 21.5 in a ~60 square degree sky region, while the BCS region contains ~280 objects with magnitude i < 20 and ~2000 objects with i < 21.5 for a 11 square degree sky region that is being observed by three current microwave Sunyaev-Zeldovich surveys. Our QSO catalog is the first one in the BCS region. Deep GALEX exposures (~2000 seconds in FUV and NUV, except in three fields) provide high signal-to-noise photometry in the GALEX bands (FUV, NUV < 24.5 mag). From this data, we select QSO candidates using only GALEX and optical r-band photometry, using the method given by Atlee and Gould (2008). In the Stripe 82 field, 60% (30%) of the GALEX selected QSO's with optical magnitude i<20 (i<21.5) also appear in the Richards et al. (2008) QSO catalog constructed using 5-band optical SDSS photometry. Comparison with the same catalog by Richards et al. shows that the completeness of the sample is approximately 40%(25%). However, for regions of the sky with very low dust extinction, like the BCS 23hr field and the Stripe 82 between 0 and 10 degrees in RA, our completeness is close to 95%, demonstrating that deep GALEX observations are almost as efficient as multi-wavelength observations at finding QSO's. GALEX observations thus provide a viable alternate route to QSO catalogs in sky regions where u-band optical photometry is not available. The full catalog is available at http://www.ice.csic.es/personal/jimenez/PHOTOZ

We discuss the cosmological degeneracy between the Hubble parameter H(z), the age of the universe and cosmological parameters describing simple variations from the minimal LCDM model. We show that independent determinations of the Hubble parameter H(z) such as those recently provided by Simon et al., Phys. Rev. D 71 (2005) 123001, combined with Cosmic Microwave Background data, can provide stringent constraints on possible deviations from the LCDM model. In particular we find that this data combination constrains at the 68% (95%) confidence level the following parameters: sum of the neutrino masses < 0.5 (1.0) eV, number of relativistic neutrino species N_rel = 4.1^{+0.4}_{-0.9} (^{+1.1}_{-1.5}), dark energy equation of state parameter w = -0.95 +- 0.17 (+- 0.32), and curvature \Omega_k = 0.002 +- 0.006 (+- 0.014), in excellent agreement with dataset combinations involving Cosmic Microwave Background, Supernovae and Baryon Acoustic Oscillations. This offers a valuable consistency check for systematic errors.

The availability of high-quality spectra for a large number of galaxies in the SDSS survey allows for a more sophisticated extraction of information about their stellar populations than, e.g., the luminosity weighted age. Indeed, sophisticated and robust techniques to fully analyze galaxy spectra have now reached enough maturity as to trust their results and findings. By reconstructing the star formation and metallicity history of galaxies from the SDSS fossil record and analyzing how it relates to its environment, we have learned how to classify galaxies: to first order the evolution of a galaxy is determined by its present stellar mass, which in turn seems to be governed by the merger rate of dark halos.

We use the VESPA algorithm and spectra from the Sloan Digital Sky Survey to investigate the star formation history of the host galaxies of 257 Type Ia supernovae. We find 5$\sigma$ evidence for a short-lived population of progenitors with lifetimes of less than 180 Myr, indicating a Type Ia supernova channel arising from stars in the mass range $\sim$3.5-8 $M_\odot$. As standardizeable candles, Type Ia supernovae play an important role in determining the expansion history of the Universe, but to be useful for future cosmological surveys, the peak luminosity needs to be free of uncorrected systematic effects at the level of 1-2%. If the different progenitor routes lead to supernovae with even moderately small differences in properties, then these need to be corrected for separately, or they could lead to a systematic bias in future supernovae surveys, as the prompt route is likely to increase in importance at high redshift. VESPA analysis of hosts could be a valuable tool in this, by identifying which progenitor route is most likely.

Using the SDSS spectroscopic sample, we estimate the dark matter halo spin parameter lambda for ~53,000 disk galaxies for which MOPED star formation histories are available. We investigate the relationship between spin and total stellar mass, star formation history, and environment. First, we find a clear anti-correlation between stellar mass and spin, with low mass galaxies generally having high dark matter spins. Second, galaxies which have formed more than ~5% of their stars in the last 0.2 Gyr have more broadly distributed and typically higher spins (including a significant fraction with lambda > 0.1) than galaxies which formed a large fraction of their stars more than 10 Gyr ago. Finally, we find little or no correlation between the value of spin of the dark halo and environment as determined both by proximity to a new cluster catalog and a marked correlation study. This agrees well with the predictions from linear hierarchical torquing theory and numerical simulations.

We use a Press-Schechter-like calculation to study how the abundance of voids changes in models with non-Gaussian initial conditions. While a positive skewness increases the cluster abundance, a negative skewness does the same for the void abundance. We determine the dependence of the void abundance on the non-Gaussianity parameter fnl for the local-model bispectrum-which approximates the bispectrum in some multi-field inflation models-and for the equilateral bispectrum, which approximates the bispectrum in e.g. string-inspired DBI models of inflation. We show that the void abundance in large-scale-structure surveys currently being considered should probe values as small as fnl < 10 and fnl^eq < 30, over distance scales ~10 Mpc.

We compare the results of the mark correlation analysis of galaxies in a sample from the Sloan Digital Sky Survey and from two galaxy catalogs obtained by semi-analytical galaxy formation models implemented on the Millennium Simulation. We use the MOPED method to retrieve the star formation history of observed galaxies and use star formation parameters as weights to the mark correlations. We find an excellent match between models and observations when the mark correlations use stellar mass and luminosity as weights. The most remarkable result is related to the mark correlations associated to the evolution of mass assembly through star formation in galaxies, where we find that semi-analytical models are able to reproduce the main trends seen in the observational data. In addition, we find a good agreement between the redshift evolution of the mean total mass formed by star formation predicted by the models and that measured by MOPED. Our results show that close galaxy pairs today formed more stellar mass ~10 Gyr ago than the average, while more recently this trend is the opposite, with close pairs showing low levels of star formation activity. We also show a strong correlation in simulations between the shape and time evolution of the star formation marks and the number of major mergers experienced by galaxies, which drive the environmental dependence in galaxy formation by regulating the star formation process.

We present the time evolution of the stellar metallicity for SDSS galaxies, a sample that spans five orders of magnitude in stellar mass (10^7 - 10^{12} Msun). Assuming the BC03 stellar population models, we find that more massive galaxies are more metal-rich than less massive ones at all redshifts; the mass-metallicity relation is imprinted in galaxies from the epoch of formation. For galaxies with present stellar masses > 10^{10} Msun, the time evolution of stellar metallicity is very weak, with at most 0.2-0.3 dex over a Hubble time- for this reason the mass-metallicity relation evolves little with redshift. However, for galaxies with present stellar masses < 10^{10} Msun, the evolution is significant, with metallicity increasing by more than a decade from redshift 3 to the present. By being able to recover the metallicity history, we have managed to identify the origin of a recent discrepancy between the metallicity recovered from nebular lines and absorption lines. As expected, we show that the young population dominates the former while the old population the latter. We have investigated the dependence on the stellar models used and find that older stellar population synthesis codes do not produce a clear result. Finally, we have explored the relationship between cluster environment and metallicity, and find a strong correlation in the sense that galaxies in high density regions have high metallicity.

Four large-area Sunyaev-Zeldovich (SZ) experiments -- APEX-SZ, SPT, ACT, and Planck -- promise to detect clusters of galaxies through the distortion of Cosmic Microwave Background photons by hot (> 10^6 K) cluster gas (the SZ effect) over thousands of square degrees. A large observational follow-up effort to obtain redshifts for these SZ-detected clusters is under way. Given the large area covered by these surveys, most of the redshifts will be obtained via the photometric redshift (photo-z) technique. Here we demonstrate, in an application using ~3000 SDSS stripe 82 galaxies with r<20, how the addition of GALEX photometry (FUV, NUV) greatly improves the photometric redshifts of galaxies obtained with optical griz or ugriz photometry. In the case where large spectroscopic training sets are available, empirical neural-network-based techniques (e.g., ANNz) can yield a photo-z scatter of $\sigma_z = 0.018 (1+z)$. If large spectroscopic training sets are not available, the addition of GALEX data makes possible the use simple maximum likelihood techniques, without resorting to Bayesian priors, and obtains $\sigma_z=0.04(1+z)$, accuracy that approaches the accuracy obtained using spectroscopic training of neural networks on ugriz observations. This improvement is especially notable for blue galaxies. To achieve these results, we have developed a new set of high resolution spectral templates based on physical information about the star formation history of galaxies. We envision these templates to be useful for the next generation of photo-z applications. We make our spectral templates and new photo-z catalogs available to the community at http://www.ice.csic.es/personal/jimenez/PHOTOZ .

Using the star formation rates from the SDSS galaxy sample, extracted using the MOPED algorithm, and the empirical Kennicutt law relating star formation rate to gas density, we calculate the time evolution of the gas fraction as a function of the present stellar mass. We show how the gas-to-stars ratio varies with stellar mass, finding good agreement with previous results for smaller samples at the present epoch. For the first time we show clear evidence for progressive gas loss with cosmic epoch, especially in low-mass systems. We find that galaxies with small stellar masses have lost almost all of their cold baryons over time, whereas the most massive galaxies have lost little. Our results also show that the most massive galaxies have evolved faster and turned most of their gas into stars at an early time, thus strongly supporting a downsizing scenario for galaxy evolution.

We discuss the cosmological degeneracy between the age of the Universe, the Hubble parameter and the effective number of relativistic particles N_eff. We show that independent determinations of the Hubble parameter H(z) as those recently provided by Simon,Verde, Jimenez (2006), combined with other cosmological data sets can provide the most stringent constraint on N_eff, yielding N_eff=3.7 (-1.2) (+1.1) at 95% confidence level. A neutrino background is detected with high significance: N_eff >1.8 at better than 99% confidence level. Constraints on the age of the universe in the framework of an extra background of relativistic particles are improved by a factor 3.

We show that cosmological observables can constrain the topology of the compact additional dimensions predicted by string theory. To do this, we develop a general strategy for relating cosmological observables to the microscopic parameters of the potentials and field-dependent kinetic terms of the multiple scalar fields that arise in the low-energy limit of string theory. We apply this formalism to the Large Volume Scenarios in Type IIB flux compactifications where analytical calculations are possible. Our methods generalize to other settings.

We use nearby K dwarf stars to measure the helium-to-metal enrichment ratio, a diagnostic of the chemical history of the Solar Neighbourhood. Our sample of K dwarfs has homogeneously determined effective temperatures, bolometric luminosities and metallicities, allowing us to fit each star to the appropriate stellar isochrone and determine its helium content indirectly. We use a newly computed set of Padova isochrones which cover a wide range of helium and metal content. Our theoretical isochrones have been checked against a congruous set of main sequence binaries with accurately measured masses, to discuss and validate their range of applicability. We find that the stellar masses deduced from the isochrones are usually in excellent agreement with empirical measurements. Good agreement is also found with empirical mass-luminosity relations. Despite fitting the masses of the stars very well, we find that anomalously low helium content (lower than primordial helium) is required to fit the luminosities and temperatures of the metal poor K dwarfs, while more conventional values of the helium content are derived for the stars around solar metallicity. We have investigated the effect of diffusion in stellar models and LTE assumption in deriving metallicities. Neither of these is able to resolve the low helium problem alone and only marginally if the cumulated effects are included, unless we assume a mixing-length which is strongly decreasing with metallicity. Further work in stellar models is urgently needed. The helium-to-metal enrichment ratio is found to be Delta Y / Delta Z = 2.1 +/- 0.9 around and above solar metallicity, consistent with previous studies, whereas open problems still remain at the lowest metallicities. Finally, we determine the helium content for a set of planetary host stars.

At the epoch of reionization, when the high-redshift inter-galactic medium (IGM) is being enriched with metals, the 63.2 micron fine structure line of OI is pumped by the ~ 1300 AA soft UV background and introduces a spectral distortion in the Cosmic Microwave Background (CMB). Here we use a toy model for the spatial distribution of neutral oxygen, assuming metal bubbles surround dark matter halos, and compute the fluctuations of this distortion, and the angular power spectrum it imprints on the CMB. We discuss the dependence of the power spectrum on the velocity of the winds polluting the IGM with metals, the minimum mass of the halos producing these winds, and on the cosmic epoch when the OI pumping occurs. We find that, although the clustering signal of the CMB distortion is weak \delta y_{rms} ~ 10^{-7} (roughly corresponding to a temperature anisotropy of few nK), it may be reachable in deep integrations with high-sensitivity infrared detectors. Even without a detection, these instruments should be able to useful constraints on the heavy element enrichment history of the IGM.

We introduce VErsatile SPectral Analysis (VESPA): a new method which aims to recover robust star formation and metallicity histories from galactic spectra. VESPA uses the full spectral range to construct a galaxy history from synthetic models. We investigate the use of an adaptative parametrization grid to recover reliable star formation histories on a galaxy-by-galaxy basis. Our goal is robustness as opposed to high resolution histories, and the method is designed to return high time resolution only where the data demand it. In this paper we detail the method and we present our findings when we apply VESPA to synthetic and real Sloan Digital Sky Survey (SDSS) spectroscopic data. We show that the number of parameters that can be recovered from a spectrum depends strongly on the signal-to-noise, wavelength coverage and presence or absence of a young population. For a typical SDSS sample of galaxies, we can normally recover between 2 to 5 stellar populations. We find very good agreement between VESPA and our previous analysis of the SDSS sample with MOPED.

We compute the effect of local electrons on the CMB temperature anisotropies. The number density and distribution of free electrons in our Galaxy has been accurately measured from pulsar dispersion measurements. Because of their distribution, the dynamics of our Galaxy and the Galaxy peculiar velocity with respect to the Hubble flow, these free electrons leave a frequency-independent imprint on the cosmic microwave background (CMB). In particular, the coherent motion of the free electrons respect to us and to the CMB rest frame produce a kinetic Sunyaev-Zeldovich signal. We compute this effect and we note that the large-scale antisymmetry of the signal gives it an angular power spectrum with a sawtooth pattern where even multipoles are suppressed with respect to the odd ones. We find the signal to be small (~ 2 uK) and sub-dominant compared to the primary CMB and other foreground signals. However, since there are no free parameters in the modeling of this signal, it can be taken into account if more precise measurements of the primordial signal are required.

We present the results of a MOPED analysis of ~3 x 10^5 galaxy spectra from the Sloan Digital Sky Survey Data Release Three (SDSS DR3), with a number of improvements in data, modelling and analysis compared with our previous analysis of DR1. The improvements include: modelling the galaxies with theoretical models at a higher spectral resolution of 3\AA; better calibrated data; an extended list of excluded emission lines, and a wider range of dust models. We present new estimates of the cosmic star formation rate, the evolution of stellar mass density and the stellar mass function from the fossil record. In contrast to our earlier work the results show no conclusive peak in the star formation rate out to a redshift around 2 but continue to show conclusive evidence for `downsizing' in the SDSS fossil record. The star formation history is now in good agreement with more traditional instantaneous measures. The galaxy stellar mass function is determined over five decades of mass, and an updated estimate of the current stellar mass density is presented. We also investigate the systematic effects of changes in the stellar population modelling, the spectral resolution, dust modelling, sky lines, spectral resolution and the change of data set. We find that the main changes in the results are due to the improvements in the calibration of the SDSS data, changes in the initial mass function and the theoretical models used.

We consider the pumping of the $63.2 \mu$m fine structure line of neutral OI in the high--redshift intergalactic medium (IGM), in analogy with the Wouthuysen--Field effect for the 21cm line of cosmic HI. We show that the soft UV background at $\sim 1300$\AA can affect the population levels, and if a significant fraction of the IGM volume is filled with ``fossil HII regions'' containing neutral OI, then this can produce a non--negligible spectral distortion in the cosmic microwave background (CMB). OI from redshift $z$ is seen in emission at $(1+z)63.2\mu$m, and between $7<z<10$ produces a mean spectral distortion of the CMB with a $y$--parameter of $y=(10^{-9} - 3\times10^{-8}) (Z/10^{-3}{\rm Z_{\odot}}) (I_{UV})$, where $Z$ is the mean metallicity of the IGM and $I_{UV}$ is the UV background at 1300\AA in units of $10^{-20}$ erg/s/Hz/cm$^2$/sr. Because O is in charge exchange equilibrium with H, a measurement of this signature can trace the metallicity at the end of the dark ages, prior to the completion of cosmic reionization and is complementary to cosmological 21cm studies. While future CMB experiments, such as Planck could constrain the metallicity to the $10^{-2} Z_{\odot}$ level, specifically designed experiments could potentially achieve a detection. Fluctuations of the distortion on small angular scale may also be detectable.

We present the first results from a study designed to test whether, given high-quality spectrophotometry spanning the mid-UV--optical wavelength regime, it is possible to distinguish the metal content (Z) and star-formation history (sfh) of individual elliptical galaxies with sufficient accuracy to establish whether their formation history is linked to their detailed morphology and position on the Fundamental Plane. From a detailed analysis of UV-optical spectrophotometry of the `cuspy' elliptical galaxy NGC 3605 and the giant elliptical NGC 5018 we find that: 1) optical spectra with l > 3500 A may not contain sufficient data to robustly uncover all the stellar populations present in individual galaxies, even in such relatively passive objects as ellipticals, 2) the addition of the UV data approaching l = 2500 A holds the key to establishing well-constrained sfhs, from which we can infer a formation and evolution history which is consistent with their photometric properties, 3) despite the superficial similarity of their spectra, the two galaxies have very different `recent' sfhs -- the smaller, cuspy elliptical NGC 3605 contains a high-Z population of age ~= 1 Gyr, and has a position on the fundamental plane typical of the product of a low-z gas-rich merger (most likely at z ~ 0.08), while the giant elliptical NGC 5018, with a sub-solar secondary population, appears to have gained its more recent stars via mass transfer / accretion of gas from its spiral companion, 4) despite these differences in detailed history, more than 85% of the stellar mass of both galaxies is associated with an old (9-12 Gyr) stellar population of near-solar Z. This pilot study provides strong motivation for the construction and analysis of high-quality UV-optical spectra for a substantial sample of ellipticals spanning the Fundamental Plane.

We use the spectra of ~ 22,000 early-type galaxies, selected from the Sloan Digital Sky Survey, to infer the ages, metallicities and star formation histories of these galaxies. We find clear evidence of "downsizing", i.e. galaxies with larger velocity dispersion have older stellar populations. In particular, most early-type galaxies with velocity dispersion exceeding 200 km s-1 formed more than 90% of their current stellar mass at redshift z > 2.5. Therefore, star formation was suppressed around this redshift. We also show that chemical enrichment was rapid, lasting 1-2 Gyr and find evidence that [Fe/H] is sub-solar. We study the robustness of these results by comparing three different approaches: using (i) Lick absorption line indices; (ii) fitting a single-burst stellar population model to the whole spectrum (lines+continuum); and (iii) reconstructing the star formation and metallicity histories in multiple age-bins, providing a method to measure mass-weighted ages and metallicities. We find good agreement between the luminosity-weighted ages and metallicities computed with these three methods.

We use the spectra of 22,000 nearby early-type galaxies from the Sloan Digital Sky Survey (SDSS) to determine the age distribution of these galaxies as a function of their velocity dispersion sigma_v in the range 100 km/s < sigma_v < 280 km/s. We then combine the inferred age-distributions with the local abundance of spheroids, including early-type galaxies and late-type bulges, to predict the evolution of the quasar luminosity function (LF) in the redshift range 0<z<6. We make the following simple assumptions: (i) the formation of stars in each galaxy, at the epoch identified with the mean mass-weighted stellar age, is accompanied by the prompt assembly of the nuclear supermassive black hole (SMBH); (ii) the mass of the SMBH obeys the M_bh-sigma_v correlation observed in nearby galaxies; (iii) the SMBH radiates at a fraction f_Edd of the Eddington limit for a fixed duration t_Q, and is identified as a luminous quasar during this epoch, (iv) the intrinsic dispersions in the Eddington ratio and the M_bh-sigma_v relation produce a combined scatter of Delta(log L_Q) around the mean logarithmic quasar luminosity <log L_Q> at fixed sigma_v. These assumptions require that the SMBH remnants of quasars with bolometric luminosity below L_bol=10^{12.5} f_Edd L_sun reside predominantly in bulges of late type galaxies. We find that evolution of the observed quasar LF can be fit over the entire redshift range in this simple model, 0<z<6 with the choices of Delta(log L_Q)=0.6-0.9, t_Q= (6-8)x10^7 yr, and <f_Edd>=0.3-0.5. We find no evidence that any of the model parameters evolves with redshift, supporting the strong connection between the formation of stars and nuclear SMBHs in spheroids.

Using cosmological hydrodynamical simulations, we investigate the prospects of the thermal Sunyaev-Zel'dovich (tSZ) effect to detect the missing baryons in the local universe. We find that at least 80% of the tSZ luminosity is generated in collapsed structures, and that $\sim$ 70% of the remaining diffuse tSZ luminosity (i.e., $\sim 15$% of the total) comes from overdense regions with $\delta_{gas}>$10, such as filaments and superclusters. The gas present in slightly overdense and underdense regions with $\delta_{gas} < 10$, despite making up 50% of the total baryon budget, leaves very little tSZ signature: it gives rise to only $\sim$ 5% of the total tSZ luminosity. Thus, future Cosmic Microwave Background (CMB) observations will be sensitive to, at best, one half of the missing baryons, improving the current observational status, but still leaving one half unobserved. Since most of the tSZ is generated in haloes, we find a tight correlation between gas pressure and galaxy number density. This allows us to predict the CMB Comptonization from existing galaxy surveys and to forecast the tSZ effect from the local structures probed by the Two Micron All Sky Survey (2MASS) galaxy catalog.

A four-dimensional universe, arising from a flux compactification of Type IIB string theory, contains scalar fields with a potential determined by topological and geometric parameters of the internal -hidden- dimensions. We show that inflation can be realized via rolling towards the large internal volume minima that are generic in these scenarios, and we give explicit formulae relating the microscopic parameters (e.g., the Euler number of the internal space) to the cosmological observables (e.g., the spectral index). We find that the tensor-to-scalar ratio, the running of the spectral index, and the potential energy density at the minimum are related by consistency relations and are exponentially small in the number of e-foldings. Further, requiring that these models arise as low-energy limits of string theory eliminates most of them, even if they are phenomenologically valid. In this context, this approach provides a strategy for systematically falsifying stringy inflation models.

We present a mark correlation analysis of the galaxies in the Sloan Digital Sky Survey using weights provided by MOPED. The large size of the sample permits statistically significant statements about how galaxies with different metallicities and star formation histories are spatially correlated. Massive objects formed a larger fraction of their stars at higher redshifts and over shorter timescales than did less massive objects (sometimes called down-sizing). We find that those galaxies which dominated the cosmic star formation at z~3 are predominantly in clusters today, whereas galaxies which dominate the star formation at z~0 inhabit substantially lower mass objects in less dense regions today. Hence, our results indicate that star formation and chemical enrichment occured first in the denser regions of the Universe, and moved to less dense regions at later times.

We study the constraints that future multifrequency Cosmic Microwave Background (CMB) experiments will be able to set on the metal enrichment history of the Inter Galactic Medium at the epoch of reionisation. We forecast the signal to noise ratio for the detection of the signal introduced in the CMB by resonant scattering off metals at the end of the Dark Ages. We take into account systematics associated to inter-channel calibration, PSF reconstruction errors and innacurate foreground removal. We develop an algorithm to optimally extract the signal generated by metals during reionisation and to remove accurately the contamination due to the thermal Sunyaev-Zel'dovich effect. Although demanding levels of foreground characterisation and control of systematics are required, they are very distinct from those encountered in HI-21cm studies and CMB polarization, and this fact encourages the study of resonant scattering off metals as an alternative way of conducting tomography of the reionisation epoch. An ACT-like experiment with optimistic assumtions on systematic effects, and looking at clean regions of the sky, can detect changes of 3%-12% (95% c.l.) of the OIII abundance (with respect its solar value) in the redshift range $z\in$ [12,22], for reionization redshift $z_{\rm re}>10$. However, for $z_{\rm re} <10$, it can only set upper limits on NII abundance increments of $\sim$ 60% its solar value in the redshift range $z\in$ [5.5,9], (95% c.l.). These constraints assume that inter-channel calibration is accurate down to one part in $10^{4}$, which constitutes the most critical technical requirement of this method, but still achievable with current technology.

Four recent observational results have challenged our understanding of high--redshift galaxies, as they require the presence of far more ultraviolet photons than should be emitted by normal stellar populations. First, there is significant ultraviolet emission from Lyman Break Galaxies (LBGs) at wavelenghts shorter than 912\AA. Second, there is strong Lyman alpha emission from extended ``blobs'' with little or no associated apparent ionizing continuum. Third, there is a population of galaxies with unusually strong Lyman-alpha emission lines. And fourth, there is a strong HeII (1640 \AA) emission line in a composite of LBGs. The proposed explanations for the first three observations are internally inconsistent, and the fourth puzzle has remained hitherto unexplained. Here we show that all four problems are resolved simultaneously if 10-30 percent of the stars in many galaxies at z ~ 3-4 are mainly primordial - unenriched by elements heavier than helium ('metals'). Most models of hierarchical galaxy formation assume efficient intra--galactic metal mixing, and therefore do not predict metal-free star formation at redshifts significantly below z ~5. Our results imply that micro-mixing of metals within galaxies is inefficient on a ~ Gyr time-scale, a conclusion that can be verified with higher resolution simulations, and future observations of the HeII emission line.

We quantify the limiting factors in optimizing current-technology cosmic microwave background (CMB) polarization experiments in order to learn about inflationary physics. We consider space-based, balloon-borne and ground-based experiments. We find that foreground contamination and residuals from foreground subtraction are ultimately the limiting factors in detecting a primordial gravity wave signal. For full-sky space-based experiments, these factors hinder the detection of tensor-to-scalar ratios of r < 10^{-3} on large scales, while for ground-based experiments these factors impede the ability to apply delensing techniques. We consider ground-based/balloon-borne experiments of currently planned or proposed designs and find that it is possible for a value of r=0.01 to be measured at ~ 3-sigma level. A small space-based CMB polarization experiment, with current detector technology and full sky coverage, can detect r ~ 10^{-3} at the ~ 3-sigma level, but a markedly improved knowledge of polarized foregrounds is needed. We advocate using as wide a frequency coverage as possible in order to carry out foreground subtraction at the percent level, which is necessary to measure such a small primordial tensor amplitude. To produce a clearly detectable (>3-sigma) tensor component in a realistic CMB experiment, inflation must either involve large-field variations, \Delta\phi >~ 1 or multi-field/hybrid models. Hybrid models can be easily distinguished from large-field models due to their blue scalar spectral index. Therefore, an observation of a tensor/scalar ratio and n < 1 in future experiments with the characteristics considered here may be an indication that inflation is being driven by some physics in which the inflaton cannot be described as a fundamental field.

Irradiation of a stellar atmosphere by an external source (e.g. an AGN) changes its structure and therefore its spectrum. Using a state-of-the-art stellar atmosphere code, we calculate the infrared spectra of such irradiated and transformed stars. We show that the original spectrum of the star, which is dominated by molecular bands, changes dramatically when irradiated even by a low-luminosity AGN ($L_{\rm X} = 10^{33}$ erg s$^{-1}$), becoming dominated by atomic lines in absorption. We study the changes in the spectrum of low-mass carbon- and oxygen-rich giant stars as they are irradiated by a modest AGN, similar to the one at the Galactic center (GC). The resulting spectra are similar to those of the faintest S-cluster stars observed in the GC. The spectrum of a star irradiated by a much brighter AGN, like that powered by a tidally disrupted star, is very different from that of any star currently observed near the GC. For the first time we have discovered that the structure of the atmosphere of an irradiated giant changes dramatically and induces a double inversion layer. We show that irradiation at the current level can explain the observed trend of CO band intensities decreasing as a function of increasing proximity to Sg $A^{*}$. This may indicate that (contrary to previous claims) there is no paucity of old giants in the GC, which coexist simultaneously with young massive stars.

We show how to place constraints on cluster physics by stacking the weak lensing signals from multiple clusters found through the Sunyaev-Zeldovich (SZ) effect. For a survey that covers about 200 sq. deg. both in SZ and weak lensing observations, the slope and amplitude of the mass vs. SZ luminosity relation can be measured with few percent error for clusters at z~0.5. This can be used to constrain cluster physics, such as the nature of feedback. For example, we can distinguish a pre-heated model from a model with a decreased accretion rate at more than 5sigma. The power to discriminate among different non-gravitational processes in the ICM becomes even stronger if we use the central Compton parameter y_0, which could allow one to distinguish between models with pre-heating, SN feedback and AGN feedback, for example, at more than 5sigma. Measurement of these scaling relations as a function of redshift makes it possible to directly observe e.g., the evolution of the hot gas in clusters. With this approach the mass-L_SZ relation can be calibrated and its uncertainties can be quantified, leading to a more robust determination of cosmological parameters from clusters surveys. The mass-L_SZ relation calibrated in this way from a small area of the sky can be used to determine masses of SZ clusters from very large SZ-only surveys and is nicely complementary to other techniques proposed in the literature.

In the context of a cosmological study of the bulk flows in the Universe, we present a detailed study of the statistical properties of the kinematic Sunyaev-Zel'dovich (kSZ) effect. We first compute analytically the correlation function and the power spectrum of the projected peculiar velocities of galaxy clusters. By taking into account the spatial clustering properties of these sources, we perform a line-of-sight computation of the {\em all-sky} kSZ power spectrum and find that at large angular scales ($l<10$), the local bulk flow should leave a visible signature above the Poisson-like fluctuations dominant at smaller scales, while the coupling of density and velocity fluctuations should give much smaller contribution. We conduct an analysis of the prospects of future high resolution CMB experiments (such as ACT and SPT) to detect the kSZ signal and to extract cosmological information and dark energy constraints from it. We present two complementary methods, one suitable for ``deep and narrow'' surveys such as ACT and one suitable for ``wide and shallow'' surveys such as SPT. Both methods can constraint the equation of state of dark energy $w$ to about 5-10% when applied to forthcoming and future surveys, and probe $w$ in complementary redshift ranges, which could shed some light on its time evolution. These determinations of $w$ do not rely on the knowledge of cluster masses, although they make minimal assumptions on cluster physics.

We reconstruct the shape of the primordial power spectrum (PPS) using a smoothing spline. Our adapted smoothing spline technique provides a complementary method to existing efforts to search for smooth features in the PPS, such as a running spectral index. With this technique we find no significant indication with WMAP first-year data that the PPS deviates from Harrison-Zeldovich and no evidence for loss of power on large scales. We also examine the effect on the cosmological parameters of the additional PPS freedom. Smooth variations in the PPS are not significantly degenerate with other cosmological parameters, but the spline reconstruction greatly increases the errors on the optical depth and baryon fraction.

We present an attempt to constrain deviations from the gravitational inverse-square law on large-scale structure scales. A perturbed law modifies the Poisson equation, which implies a scale-dependent growth of overdensities in the linear regime and thus modifies the power spectrum shape. We use two large-scale structure surveys (the Sloan Digital Sky survey and the Anglo-Australian Two-degree field galaxy redshift survey) to constrain the parameters of two simple modifications of the inverse-square law. We find no evidence for deviations from normal gravity on the scales probed by these surveys (~ 10^(23) m.)

We present an algorithm that allows fast and efficient detection of transits, including planetary transits, from light-curves. The method is based on building an ensemble of fiducial models and compressing the data using the MOPED algorithm. We describe the method and demonstrate its efficiency by finding planet-like transits in simulated Pan-STARRS light-curves. We show that that our method is independent of the size of the search space of transit parameters. In large sets of light-curves, we achieve speed up factors of order of $10^{8}$ times over the full $\chi2$ search. We discuss how the algorithm can be used in forthcoming large surveys like Pan-STARRS and LSST and how it may be optimized for future space missions like Kepler and COROT where most of the processing must be done on board.

We develop a formalism to characterize the redshift evolution of the dark energy potential. Our formalism makes use of quantities similar to the Horizon-flow parameters in inflation and is general enough that can deal with multiscalar quintessence scenarios, exotic matter components, and higher order curvature corrections to General Relativity. We show how the shape of the dark energy potential can be recovered non parametrically using this formalism and we present approximations analogous to the ones relevant to slow-roll inflation. Since presently available data do not allow a non-parametric and exact reconstruction of the potential, we consider a general parametric description. This reconstruction can also be used in other approaches followed in the literature (e.g., the reconstruction of the redshift evolution of the dark energy equation of state w(z)). Using observations of passively evolving galaxies and supernova data we derive constraints on the dark energy potential shape in the redshift range 0.1 < z < 1.8. Our findings show that at the 1sigma level the potential is consistent with being constant, although at the same level of confidence variations cannot be excluded with current data. We forecast constraints achievable with future data from the Atacama Cosmology Telescope.

Using the spectroscopic sample of the SDSS DR1 we measure how gas was transformed into stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a clear correlation between early star formation activity and present-day stellar mass: the more massive galaxies have formed about 80% of their stars at $z>1$, while for the less massive ones the value is only about 20%. By comparing the BCT to the dark matter merger tree, we find indications that star formation efficiency at $z>1$ had to be about a factor of two higher than today ($\sim 10%$) in galaxies with present-day stellar mass larger than $2 \times 10^{11}M_\odot$, if this early star formation occurred in the main progenitor. Therefore, the LCDM paradigm can accommodate a large number of red objects. On the other hand, in galaxies with present-day stellar mass less than $10^{11}$ M$_{\odot}$, efficient star formation seems to have been triggered at $z \sim 0.2$. We show that there is a characteristic mass (M$_* \sim 10^{10}$ M$_{\odot}$) for feedback efficiency (or lack of star formation). For galaxies with masses lower than this, feedback (or star formation suppression) is very efficient while for higher masses it is not. The BCT, determined here for the first time, should be an important observable with which to confront theoretical

Using the MOPED algorithm we determine non-parametrically the Stellar Mass Function of 96,545 galaxies from the Sloan Digital Sky Survey data release one. By using the reconstructed spectrum due to starlight we can eliminate contamination from either emission lines or AGN components. Our results give excellent agreement with previous works, but extend their range by more than two decades in mass to 10^{7.5} < M_s/h^{-2}M_\odot < 10^{12}. We present both a standard Schechter fit and a fit modified to include an extra, high-mass contribution, possibly from cluster cD galaxies. The Schechter fit parameters are $\phi^*=(7.7\pm 0.8)\times 10^{-3} h^3 Mpc^{-3}$, $M^*=(7.53 \pm 0.04) \times 10^{10} h^{-2}M_\odot$ and $\alpha=-1.167\pm 0.004$. Our sample also yields an estimate for the contribution from baryons in stars to the critical density of $\Omega_{b*}h=(2.40 \pm 0.04) \times 10^{-3}$, in good agreement with other indicators. Error bars are statistical and a Salpeter IMF is assumed throughout. We find no evolution of the mass function in the redshift range $0.05 < z < 0.34$ indicating that almost all present day stars were already formed at $z \sim 0.34$ with little or no star formation activity since then and that the evolution seen in the luminosity function must be largely due to stellar fading.

Supersonic turbulent flows of magnetized gas are believed to play an important role in the dynamics of star-forming clouds in galaxies. Understanding statistical properties of such flows is crucial for developing a theory of star formation. In this letter we propose a unified approach for obtaining the velocity scaling in compressible and super--Alfv\'{e}nic turbulence, valid for arbitrary sonic Mach number, \ms. We demonstrate with numerical simulations that the scaling can be described with the She--L\'{e}v\^{e}que formalism, where only one parameter, interpreted as the Hausdorff dimension of the most intense dissipative structures, needs to be varied as a function of \ms. Our results thus provide a method for obtaining the velocity scaling in interstellar clouds once their Mach numbers have been inferred from observations.

We improve the algorithm of Kosowsky, Milosavljevic, and Jimenez (2002) for computing power spectra of the cosmic microwave background. The present algorithm computes not only the temperature power spectrum but also the E-mode polarization and the temperature-polarization cross power spectra, providing the accuracy required for current cosmological parameter estimation. We refine the optimum set of cosmological parameters for computing the power spectra as perturbations around a fiducial model, leading to an accuracy better than 0.5% for the temperature power spectrum throughout the region of parameter space within WMAP's first-year 3sigma confidence region. This accuracy is comparable to the difference between the widely-used CMBFAST code (Seljak and Zaldarriaga 1996) and Boltzmann codes. Our algorithm (CMBwarp) makes possible a full exploration of the likelihood region for eight cosmological parameters in about one hour on a laptop computer. We provide the code to compute power spectra as well as the Markov Chain Monte Carlo algorithm for cosmological parameters estimation at http://www.physics.upenn.edu/~raulj/CMBwarp

The determination of the star-formation history of the Universe is a key goal of modern cosmology, as it is crucial to our understanding of how structure in the Universe forms and evolves. A picture has built up over recent years, piece-by-piece, by observing young stars in distant galaxies at different times in the past. These studies indicated that the stellar birthrate peaked some 8 billion years ago, and then declined by a factor of around ten to its present value. Here we report on a new study which obtains the complete star formation history by analysing the fossil record of the stellar populations of 96545 nearby galaxies. Broadly, our results support those derived from high-redshift galaxies elsewhere in the Universe. We find, however, that the peak of star formation was more recent - around 5 billion years ago. Our study also shows that the bigger the stellar mass of the galaxy, the earlier the stars were formed. This striking result indicates a very different formation history for high- and low-mass formation.

We present a new set of stellar interior and synthesis models for predicting the integrated emission from stellar populations in star clusters and galaxies of arbitrary age and metallicity. This work differs from existing spectral synthesis codes in a number of important ways, namely (1) the incorporation of new stellar evolutionary tracks, with sufficient resolution in mass to sample rapid stages of stellar evolution; (2) a physically consistent treatment of evolution in the HR diagram, including the approach to the main sequence and the effects of mass loss on the giant and horizontal-branch phases. Unlike several existing models, ours yield consistent ages when used to date a coeval stellar population from a wide range of spectral features and colour indexes. We rigorously discuss degeneracies in the age-metallicity plane and show that inclusion of spectral features blueward of 4500 AA, suffices to break any remaining degeneracy and that with moderate S/N spectra (10 per 20AA, resolution element) age and metallicity are not degenerate. We also study sources of systematic errors in deriving the age of a single stellar population and conclude that they are not larger than 10-15%. We illustrate the use of single stellar populations by predicting the colors of primordial proto-galaxies and show that one can first find them and then deduce the form of the IMF for the early generation of stars in the universe. Finally, we provide accurate analytic fitting formulas for ultra fast computation of colors of single stellar populations. The models can be found at http://www.physics.upenn.edu/~raulj

The magnetic field strength in molecular clouds is a fundamental quantity for theories of star formation. It is estimated by Zeeman splitting measurements in a few dense molecular cores, but its volume--averaged value within large molecular clouds (over several parsecs) is still uncertain. In this work we provide a new method to constrain the average magnetic field strength in molecular clouds. We compare the power spectrum of gas density of molecular clouds with that of two $350^3$ numerical simulations of supersonic MHD turbulence. The numerical simulation with approximate equipartition of kinetic and magnetic energies (model A) yields the column density power spectrum $P(k)\propto k^{-2.25\pm 0.01}$, the super--Alfv\'{e}nic simulation (model B) $P(k)\propto k^{-2.71\pm 0.01}$. The column density power spectrum of the Perseus, Taurus and Rosetta molecular cloud complexes is found to be well approximated by a power law, $P_{\rm o}(k)\propto k^{-a}$, with $a=2.74\pm 0.07$, $2.74\pm 0.08$ and $2.76\pm 0.08$ respectively. We conclude that the observations are consistent with the presence of super--Alfv\'{e}nic turbulence in molecular clouds (model B) while model A is inconsistent (more than 99% confidence) with the observations.

Using MOPED we determine non-parametrically the star-formation and metallicity history of over 37,000 high-quality galaxy spectra from the Sloan Digital Sky Survey (SDSS) early data release. We use the entire spectral range, rather than concentrating on specific features, and we estimate the complete star formation history without prior assumptions about its form (by constructing so-called `population boxes'). The main results of this initial study are that the star formation rate in SDSS galaxies has been in decline for ~6 Gyr; a metallicity distribution for star-forming gas which is peaked ~3 Gyr ago at about solar metallicity, inconsistent with closed-box models, but consistent with infall models. We also determine the infall rate of gas in SDSS and show that it has been significant for the last 3 Gyr. We investigate errors using a Monte-Carlo Markov Chain algorithm. Further, we demonstrate that recovering star formation and metallicity histories for such a large sample becomes intractable without data compression methods, particularly the exploration of the likelihood surface. By exploring the whole likelihood surface we show that age-metallicity degeneracies are not as severe as by using only a few spectral features. We find that 65% of galaxies contain a significant old population (with an age of at least 8 Gyr), including recent starburst galaxies, and that over 97% have some stars older than 2 Gyr. It is the first time that the past star formation history has been determined from the fossil record of the present-day spectra of galaxies.

From recent CMB and Large Scale Structure observations the value of the equation of state of dark energy, assuming it to be constant in time, is constrained to be -1.3<w<-0.7 at the 95% confidence level: consistent with dark energy being a classical vacuum term. Here we describe two novel and independent methods, sensitive to different systematics, that give the same value for w and similar confidence regions. This suggests that systematics are not an issue in current determinations of w. The first method yields a measurement of $w$ that relies on the minimum number of model-dependent parameters; the second method is a non-parametric measurement of the time dependence of w(z). We also present a method to statistically determine the edge of a distribution.

The photoelectric heating is believed to be the main heating mechanism in cool HI clouds. The heating rate can be estimated through observations of the [CII] line emission, since this is the main coolant in regions where the photoelectric effect dominates the heating. Comparison of the [CII] emission with the far-infrared (FIR) emission allows to constrain the efficiency of the photoelectric heating, using model calculations that take into account the strength of the radiation field. Recent [CII] observations carried out with the ISO satellite have made this study possible. In this work we study the correlation between FUV absorption and FIR emission using three-dimensional models. The density distributions are obtained with numerical simulations of compressible magneto-hydrodynamic turbulence, with rms sonic Mach numbers 0.6<M<10. The FIR intensities are solved with detailed radiative transfer calculations. The [CII] line radiation is estimated assuming the [CII] line cooling equals the FUV absorption multiplied by the efficiency of the photoelectric heating, epsilon. The average ratio between the predicted [CII] and FIR emissions is found to be remarkably constant between different models, implying that the derived values of epsilon should not depend on the rms Mach number. The comparison with empirical data from translucent, high latitude clouds yields an estimate of the photoelectric heating efficiency of 2.9 10^-2. This value confirms previous theoretical predictions. Our models show that most of the scatter in the observed [CII] and FIR intensities can be understood as a result of the highly fragmented density field in turbulent HI clouds. The scatter can be reproduced with models with supersonic turbulence, while subsonic turbulence fails to generate the observed scatter.

We estimate the cosmic production rate of helium relative to metals ($\Delta Y/\Delta Z$) using K dwarf stars in the Hipparcos catalog with accurate spectroscopic metallicities. The best fitting value is $\Delta Y/\Delta Z=2.1 \pm 0.4$ at the 68% confidence level. Our derived value agrees with determinations from HII regions and with theoretical predictions from stellar yields with standard assumptions for the initial mass function. The amount of helium in stars determines how long they live and therefore how fast they will enrich the insterstellar medium with fresh material.

We place tight constraints on the redshift-averaged, effective value of the equation of state of dark energy, w, using only the absolute ages of Galactic stars and the observed position of the first peak in the angular power spectrum of the CMB. We find w<-0.8 at the 68% confidence level. If we further consider that w > -1, this finding suggests that within our uncertainties, dark energy is indistinguishable from a classical vacuum energy term. We detect a correlation between the ages of the oldest galaxies and their redshift. This opens up the possibility of measuring w(z) by computing the relative ages of the oldest galaxies in the universe as a function of redshift, dz/dt. We show that this is a realistic possibility by computing dz/dt at z~0 from SDSS galaxies and obtain an independent estimate for the Hubble constant, H_0 = 69 \pm 12 km s-1 Mpc-1. The small number of galaxies considered at z>0.2 does not yield, currently, a precise determination of w(z), but shows that the age--redshift relation is consistent with a Standard LCDM universe with $w=-1$.

We use a simple model of spheroid formation to explore the relationship between the creation of stars and dust in a massive proto galaxy and the growth of its central black hole. This model predicts that sub-mm luminosity peaks after only ~ 0.2 Gyr. However, without a very massive seed black hole, Eddington-limited growth means that a black hole mass of 10^9 solar masses, and hence very luminous AGN activity, cannot be produced until > 0.5 Gyr after the formation of the first massive stars in the halo. The model thus predicts a time-lag between the peak of sub-mm luminosity and AGN luminosity in a massive proto-elliptical of a few times 10^8 years. For a formation redshift z = 5, this means that powerful AGN activity is delayed until z = 3.5, by which time star formation in the host is 90% complete, and sub-mm luminosity has declined to ~ 25% of its peak value. This provides a natural explanation for why successful sub-mm detections of luminous radio galaxies are largely confined to z > 2.5. Conversely the model also predicts that while all high-redshift luminous sub-mm-selected sources should contain an active (and growing) black hole, the typical luminosity of the AGN in such objects is ~1000 times smaller than that of the most powerful AGN. This is consistent with the almost complete failure to detect sub-mm selected galaxies with existing X-ray surveys. Finally the model yields a black-hole:spheroid mass ratio which evolves rapidly in the first Gyr, but asymptotes to ~ 0.001-0.003 in agreement with results at low redshift. This ratio arises not because the AGN terminates star formation, but because fueling of the massive black hole is linked to the total mass of gas available for star formation in the host.

We explore whether the rest-frame near-UV spectral region, observable in high-redshift galaxies via optical spectroscopy, contains sufficient information to allow the degeneracy between age and metallicity to be lifted. We do this by testing the ability of evolutionary synthesis models to reclaim the correct metallicity when fitted to the near-UV spectra of F stars of known (sub-solar and super-solar) metallicity. F stars are of particular interest because the rest-frame near-UV spectra of the oldest known elliptical galaxies at z > 1 appear to be dominated by F stars near to the main-sequence turnoff. We find that, in the case of the F stars, where the HST ultraviolet spectra have high signal:noise, model-fitting with metallicity allowed to vary as a free parameter is rather successful at deriving the correct metallicity. As a result, the estimated turnoff ages of these stars yielded by the model fitting are well constrained. Encouraged by this we have fitted these same variable- metallicity models to the deep, optical spectra of the z \simeq 1.5 mJy radio galaxies 53W091 and 53W069 obtained with the Keck telescope. While the age-metallicity degeneracy is not so easily lifted for these galaxies, we find that even when metallicity is allowed as a free parameter, the best estimates of their ages are still \geq 3 Gyr, with ages younger than 2 Gyr now strongly excluded. Furthermore, we find that a search of the entire parameter space of metallicity and star formation history using MOPED (Heavens et al., 2000) leads to the same conclusion. Our results therefore continue to argue strongly against an Einstein-de Sitter universe, and favour a lambda-dominated universe in which star formation in at least these particular elliptical galaxies was completed somewhere in the redshift range z = 3 - 5.

Deep spectroscopy of the two mJy radio galaxies LBDS 53W069 and 53W091 has previously shown them to have old (>3 Gyr) stellar populations at z=1.5. Here we present the results of HST observations of these galaxies. 53W069 has an r^1/4 profile in both the F814W & F110W data with an effective radius of 0.30 arcsec (2.7 kpc). The restframe U-B colour gradient requires a stellar population of super-solar (3Z_sun) metallicity that formed on a very short timescale at high redshift (z>5). 53W091 has a regular r^1/4 profile in F110W with an effective radius of 0.32 arcsec (2.9 kpc). The F814W profile is more extended and is consistent with a blue exponential disk that contributes 20% of the flux within r_e. The restframe U-B colour gradient is significantly larger than that observed in field ellipticals at z<1, implying a stellar population of mixed metallicity (1-3Z_sun) that formed in a high-redshift rapid burst. We have compared these two LBDS radio galaxies with the Kormendy relation of ten 3CR radio galaxies at z=0.8. The LBDS galaxies follow the same relation as the more radio-luminous 3CR galaxies, assuming passive evolution of their stellar populations, although they are smaller than the 3CR galaxies whose mean effective radius is 12 kpc. Their sizes and radio luminosities are consistent with scaling relations applied to the 3CR galaxies, in which both radio power and effective radius scale with galaxy mass. Compared with a sample of z=0.4 cluster ellipticals, 53W069 & 53W091 lie well within the scatter of the Kormendy relation. We conclude that the hosts of these mJy radio sources at z=1.5 are passively-evolving elliptical galaxies that will evolve into ordinary L* ellipticals by the present day. (Abridged)

We revisit the issue of cosmological parameter estimation in light of current and upcoming high-precision measurements of the cosmic microwave background power spectrum. Physical quantities which determine the power spectrum are reviewed, and their connection to familiar cosmological parameters is explicated. We present a set of physical parameters, analytic functions of the usual cosmological parameters, upon which the microwave background power spectrum depends linearly (or with some other simple dependence) over a wide range of parameter values. With such a set of parameters, microwave background power spectra can be estimated with high accuracy and negligible computational effort, vastly increasing the efficiency of cosmological parameter error determination. The techniques presented here allow calculation of microwave background power spectra $10^5$ times faster than comparably accurate direct codes (after precomputing a handful of power spectra). We discuss various issues of parameter estimation, including parameter degeneracies, numerical precision, mapping between physical and cosmological parameters, and systematic errors, and illustrate these considerations with an idealized model of the MAP experiment.

The annihilation rate of weakly interacting cold dark matter particles at the galactic center could be greatly enhanced by the growth of a density spike around the central supermassive black hole (SBH). Here we discuss the effects of hierarchical mergers on the central spike. Mergers between halos containing SBHs lead to the formation of SBH binaries which transfer energy to the dark matter particles, lowering their density. The predicted flux of annihiliation photons from the galactic center is several orders of magnitude smaller than in models that ignore the effects of SBHs and mergers. Measurement of the annihilation radiation could in principle be used to constrain the merger history of the galaxy.

We present a comparison of the predictions of stellar models with the luminosity of the lower main sequence ($5.5 < M_V < 7.3$) using K dwarfs in the Hipparcos catalog. The parallaxes of our comparison stars are known to better than 15% and metallicities have been recently determined from photometry. A major advantage of our comparison is that distances in our sample are known with good accuracy, while tests that involve open and globular clusters are constrained by potentially inaccurate distances. We show that the luminosity of the lower main sequence relative to a fiducial (solar metallicity) isochrone is a simple function of metallicity: $\Delta M_V = 0.84375 \times {\rm [ Fe/H]} - 0.04577$. We compare the data with a range of isochrones from the literature. None of the models fit all the data, although some models do clearly better than others. In particular, metal rich isochrones seem to be difficult to construct. The relationship between luminosity, colour and metallicity for K dwarfs is found to be very tight. We are thus able to derive metallicities for K dwarfs based on their position in the Hipparcos colour-magnitude diagram with accuracies better than 0.1 dex. The metallicity-luminosity relation for K dwarfs leads to a new distance indicator with a wide range of possible applications.

We show that gas in a large fraction of low mass dark matter halos may form Toomre stable disks, if angular momentum is conserved when the gas contracts. Such halos would be stable to star formation and therefore remain dark. This may potentially explain the discrepancy between the predicted and observed number of dwarf satellites in the Local Group, as well as the deviation between the predicted and the observed faint end slope of the luminosity function. The above mechanism does not require a strong variation of the baryon fraction with the virial mass of the dark halo. We show that model fits to rotation curves are also consistent with this hypothesis: none of the observed galaxies lie in the region of parameter space forbidden by the Toomre stability criterion.

We study a large set of high spatial resolution optical rotation curves of galaxies with the goal of determining the model parameters for a disk embedded within a CDM halo that we model either with a NFW profile or pseudo-isothermal profile. We show that parameter degeneracies present in lower resolution data are lifted at these higher resolutions. 34% of the galaxies do not have a meaningful fit when using the NFW profile and 32% when using the pseudoisothermal profile, however only 14% do not have a meaningful fit in either model. In both models we find correlations between the disk baryon fraction $f_d$ and the spin parameter of the halo $\lambda'$, between $f_d$ and the dark halo mass $M_{200}$, and between $M_{200}$ and the concentration parameter $c$. We show that the distribution of the concentration parameters $c$, for a NFW halo, is in good agreement with CDM predictions; no significant galaxy population is found with very low values of $c$. The overall distribution of $\lambda'$ is in good agreement with theoretical predictions from hierarchical tidal torque theory. The whole sample is also well fitted by a pseudo-isothermal dark halo with a core, but the size of the core is rather small (for 70% of the sample the core size is less than 2 kpc). Thus we conclude that the profile of dark matter is steep ($r^{-1}$ or steeper) down to this radius; large dark matter cores (and therefore very low dark matter central densities) seem to be excluded. LSBs tend to have higher values of $\lambda'$ for a given $f_d$ and lower values of $c$ for a given mass than HSBs. In an appendix we give some useful formula for pseudo-isothermal profile halos and discuss in detail the issue of parameter degeneracies.

I present a biased review of when the epoch of formation of galaxies (both disks and ellipticals) maybe took place. I base my arguments in simple (mostly) analytic models that have been recently developed to reproduce most of the observed photometric, chemical and dynamical properties of galaxies both at low and high redshift.

In this letter we address the problem of the origin of blue cores and inverse color gradients in early-type galaxies reported in the Hubble Deep Field North and South (HDFs) by Menanteau, Abraham & Ellis 2001. We use a multi-zone single collapse model. This model accounts for the observed blue cores by adopting a broad spread in formation redshifts for ellipticals, allowing some of these galaxies begin forming no more than ~1 Gyr before the redshift of observation. The single-zone collapse model then produces cores that are bluer than the outer regions because of the increase of the local potential well toward the center which makes star-formation more extended in the central region of the galaxy. We compare the predicted V-I(r) color gradients with the observed ones using the redshift of formation ($z_F$) of the elliptical as the only free parameter. We find that the model can account with relatively good agreement for the blue cores and inverse color gradients found in many spheroidals and at the same time for the red and smooth colors profiles reported. Based on the model our analysis suggests two populations of field ellipticals, one formed recently, within $\lesssim1$Gyr and another much older formed $\gtrsim4$Gyr since the redshift of observation.

We propose to use relative galaxy ages as a means of constraining cosmological parameters. By measuring the age difference between two ensembles of old galaxies at somewhat different redshifts, one could determine the derivative of redshift with respect to cosmic time, dz/dt. At high redshifts, z=1-2, this measurement would constrain the equation-of-state of the dark energy, while at low redshifts, z< 0.2, it would determine the Hubble constant, H_0. The selected galaxies need to be passively-evolving on a time much longer than their age difference.

We determine the distribution of total energy emitted by gamma-ray bursts for bursts with fluences and distance information. Our core sample consists of eight bursts with BATSE spectra and spectroscopic redshifts. We extend this sample by adding four bursts with BATSE spectra and host galaxy R magnitudes. From these R magnitudes we calculate a redshift probability distribution; this method requires a model of the host galaxy population. From a sample of ten bursts with both spectroscopic redshifts and host galaxy R magnitudes (some do not have BATSE spectra) we find that the burst rate is proportional to the galaxy luminosity at the epoch of the burst. Assuming that the total energy emitted has a log-normal distribution, we find that the average emitted energy (assumed to be radiated isotropically) is $<E_{\gamma iso} > = 1.3^{+1.2}_{-1.0} \times 10^{53}$ ergs (for H$_0$ = 65 km s$^{-1}$ Mpc$^{-1}$, $\Omega_m=0.3$ and $\Omega_\Lambda=0.7$); the distribution has a logarithmic width of $\sigma_\gamma=1.7^{+0.7}_{-0.3}$. The corresponding distribution of X-ray afterglow energy (for seven bursts) has $<E_{X iso} > = 4.0^{+1.6}_{-1.8} \times 10^{51}$ergs and $\sigma_X = 1.3^{+0.4}_{-0.3}$. For completeness, we also provide spectral fits for all bursts with BATSE spectra for which there were afterglow searches.

We investigate the relative sensitivities of several tests for deviations from Gaussianity in the primordial distribution of density perturbations. We consider models for non-Gaussianity that mimic that which comes from inflation as well as that which comes from topological defects. The tests we consider involve the cosmic microwave background (CMB), large-scale structure (LSS), high-redshift galaxies, and the abundances and properties of clusters. We find that the CMB is superior at finding non-Gaussianity in the primordial gravitational potential (as inflation would produce), while observations of high-redshift galaxies are much better suited to find non-Gaussianity that resembles that expected from topological defects. We derive a simple expression that relates the abundance of high-redshift objects in non-Gaussian models to the primordial skewness.

We predict the Tully-Fisher (TF) and surface-brightness--magnitude relation for disk galaxies at z=3 and discuss the origin of these scaling relations and their scatter. We show that the variation of the TF relation with redshift can be a potentially powerful discriminator of galaxy-formation models. In particular, the TF relation at high redshift might be used to break parameter degeneracies among galactosynthesis models at z=0, as well as to constrain the redshift distribution of collapsing dark-matter halos, the star-formation history and baryon fraction in the disk, and the distribution of halo spins.

We derive physical parameters of galaxies from their observed spectrum, using MOPED, the optimized data compression algorithm of Heavens, Jimenez & Lahav 2000. Here we concentrate on parametrising galaxy properties, and apply the method to the NGC galaxies in Kennicutt's spectral atlas. We focus on deriving the star formation history, metallicity and dust content of galaxies. The method is very fast, taking a few seconds of CPU time to estimate 17 parameters, and so specially suited to study of large data sets, such as the Anglo-Australian 2 degree field galaxy survey and the Sloan Digital Sky Survey. Without the power of MOPED, the recovery of star formation histories in these surveys would be impractical. In the Kennicutt atlas, we find that for the spheroidals a small recent burst of star formation is required to provide the best fit to the spectrum. There is clearly a need for theoretical stellar atmospheric models with spectral resolution better than 1\AA if we are to extract all the rich information that large redshift surveys contain in their galaxy spectra.

I review the statistical techniques needed to extract information about physical parameters of galaxies from their observed spectra. This is important given the sheer size of the next generation of large galaxy redshift surveys. Going to the opposite extreme I review what we can learn about the nature of the primordial density field from observations of high--redshift objects.

Given sufficiently deep optical spectroscopy, the age estimation of high-redshif t ($z > 1$) galaxies has been claimed to be a relatively robust process (e.g. Dunlop et al. 1996) due to the fact that, for ages $< 5$Gyr, the near-ultraviolet light of a stellar population is expected to be dominated by `well-understood' main-sequence (MS) stars. Recently, however, the reliability of this process has been called into question by Yi et al (2000), who claim to have developed models in which the spectrum produced by the main sequence reddens much more rapidly than in the models of Jimenez et al (2000a), leading to much younger age estimates for the reddest known high-redshift ellipticals. In support of their revised age estimates, Yi et al cite the fact that their models can reproduce the spectrum of the Sun at an age of 5 Gyr, whereas the solar spectrum is not reproduced by the Jimenez et al models until $\simeq 10$ Gyr. Here we confirm this discrepancy, but point out that this is in fact a {\it strength} of the Jimenez et al models and indicative of some flaw in the models of Yi et al (which, in effect, imply that the Sun will turn into a red giant any minute now). We have also explored the models of Worthey (1994) (which are known to differ greatly from those of Jimenez et al in the treatment of post-MS evolution) and find that the main-sequence component of Worthey's models also cannot reproduce the solar spectrum until an age of 9-10 Gyr. We conclude that either the models of Yi et al are not as main-sequence dominated at 4-5 Gyr as claimed, or that the stellar evolutionary timescale in these models is in error by a factor possibly as high as two. (abridged)

We present the results of fitting deep off-nuclear optical spectroscopy of radio-quiet quasars, radio-loud quasars and radio galaxies at z ~ 0.2 with evolutionary synthesis models of galaxy evolution. Our aim was to determine the age of the dynamically dominant stellar populations in the hos t galaxies of these three classes of powerful AGN. Some of our spectra display residual nuclear contamination at the shortest wavelengths, but the detailed quality of the fits longward of the 4000A break provide unequivocal proof, if further proof were needed, that quasars lie in massive galaxies with (at least at z ~ 0.2) evolved stellar populations. By fitting a two-component model we have separated the very blue (starburst and/or AGN contamination) from the redder underlying spectral energy distribution, and find that the hosts of all three classes of AGN are dominated by old stars of age 8 - 14 Gyr. If the blue component is attributed to young stars, we find that, at most, 1% of the baryonic mass of these galaxies is involved in star-formation activity at the epoch of observation. These results strongly support the conclusion reached by McLure et al. (1999) that the host galaxies of luminous quasars are massive ellipticals which formed prior to the peak epoch of quasar activity at z ~ 2.5.

We present images and spectroscopy of objects close to the sightline of Q2302+029 in order to search for galaxies responsible for the remarkable z=0.7 high-ionization absorption line system found by Jannuzi et al (1996). This system shows `normal' narrow OVI, NV, and CIV lines superimposed on broader (3,000-5,000 km/s wide), unsaturated absorption troughs some 56,000 km/s away from the QSO emission redshift (z=1.052). Despite reaching sensitivities sufficient to detect 1/10 L* galaxies in the optical and 1/20 L* in the infra-red, we are unable to detect any obvious bright galaxies which might be responsible for the absorption beyond ~ 6 h-1 kpc of the sightline. This supports the hypotheses that the absorption is either intrinsic to the QSO, or arises in intracluster gas. Adopting either explanation is problematic: in the first case, `associated' absorption at such high ejection velocities is hard to understand, and challenges the conventional discrimination between intrinsic and intervening absorbers; in the second case, the gas must reside in a ~ 40 h-1 Mpc long filament aligned along the line of sight in order to reproduce the broad absorption. Spectroscopy of objects beyond the immediate vicinity of the QSO sightline reveals a galaxy cluster at z=0.59, which coincides with strong Lya and more narrow high ionization lines in the quasar spectrum. Here too, the lack of galaxies at distances comparable to those found for, e.g., Lya-absorbing galaxies, suggests that the absorption may arise from intracluster gas unassociated with any individual galaxies.

We investigate the effects of a variety of ingredients that must enter into a realistic model for disk-galaxy formation, focusing primarily on the Tully-Fisher (TF) relation and its scatter in several wavebands. Our main findings are: (a) the slope, normalization, and scatter of the TF relation across various wavebands is determined {\em both} by halo properties and star formation in the disk; (b) TF scatter owes primarily to the spread in formation redshifts. The scatter can be measurably reduced by chemical evolution, and also in some cases by the weak anti-correlation between peak height and spin; (c) multi-wavelength constraints can be important in distinguishing between models which appear to fit the TF relation in I or K; (d) successful models seem to require that the bulk of disk formation cannot occur too early (z>2) or too late (z<0.5), and are inconsistent with high values of $\Omega_0$; (e) a realistic model with the above ingredients can reasonably reproduce the observed z=0 TF relation in {\em all} bands (B, R, I, and K). It can also account for the z=1 B-band TF relation and yield rough agreement with the local B and K luminosity functions and B-band surface-brightness--magnitude relation. The remarkable agreement with observations suggests that the amount of gas that is expelled or poured into a disk galaxy must be small, and that the specific angular momentum of the baryons must roughly equal that of the halo; there is little room for angular momentum transfer. In an appendix we present analytic fits to stellar-population synthesis models.

We present a method for radical linear compression of datasets where the data are dependent on some number $M$ of parameters. We show that, if the noise in the data is independent of the parameters, we can form $M$ linear combinations of the data which contain as much information about all the parameters as the entire dataset, in the sense that the Fisher information matrices are identical; i.e. the method is lossless. We explore how these compressed numbers fare when the noise is dependent on the parameters, and show that the method, although not precisely lossless, increases errors by a very modest factor. The method is general, but we illustrate it with a problem for which it is well-suited: galaxy spectra, whose data typically consist of $\sim 10^3$ fluxes, and whose properties are set by a handful of parameters such as age, brightness and a parametrised star formation history. The spectra are reduced to a small number of data, which are connected to the physical processes entering the problem. This data compression offers the possibility of a large increase in the speed of determining physical parameters. This is an important consideration as datasets of galaxy spectra reach $10^6$ in size, and the complexity of model spectra increases. In addition to this practical advantage, the compressed data may offer a classification scheme for galaxy spectra which is based rather directly on physical processes.

The observed abundance of high-redshift galaxies and clusters contains precious information about the properties of the initial perturbations. We present a method to compute analytically the number density of objects as a function of mass and redshift for a range of physically motivated non-Gaussian models. In these models the non-Gaussianity can be dialed from zero and is assumed to be small. We compute the probability density function for the smoothed dark matter density field and we extend the Press and Schechter approach to mildly non-Gaussian density fields. The abundance of high-redshift objects can be directly related to the non-Gaussianity parameter and thus to the physical processes that generated deviations from the Gaussian behaviour. Even a skewness parameter of order 0.1 implies a dramatic change in the predicted abundance of $z\gap 1$ objects. Observations from NGST and X-ray satellites (XMM) can be used to accurately measure the amount of non-Gaussianity in the primordial density field.

The first stars to form in protogalaxies must have primordial chemical composition. We refer to a protogalaxy that is forming stars of primordial composition, or very low metallicity ($Z\le 0.01 Z_{\odot}$), as ``Primordial Proto-Galaxy'' (PPG). PPGs contain little or no dust, and therefore their spectral energy distribution can be modelled from the rest-frame ultraviolet to the infrared without accounting for dust extinction and emission. We present the results of computing the photometric properties of high redshift PPGs at near-d-to-mid infrared wavelengths, that will soon be available with the new generation of infrared space telescopes, such as SIRTF and the NGST. We show that: i) PPGs at very high redshift (5<z<10) should be easily selected from deep near/mid IR surveys with a colour-colour criterion; ii) PPGs at redshift 5<z<10 can be detected at 8 $\mu$m with the NGST, if they have a constant star formation rate of at least 100 M$_{\odot}/$yr; iii) once the redshift of a PPG photometric candidate is determined, its near-to-mid infrared colors should provide strong constraints on the stellar IMF at zero or very low metallicity.

We have computed the evolving ultraviolet-millimeter spectral energy distribution (SED) produced by proto-galactic starbursts at high redshift, incorporating the chemical evolution of the interstellar medium in a consistent manner. Dust extinction is calculated in a novel way, that is not based on empirical calibrations of extinction curves, but rather on the lifetime of molecular clouds which delays the emergence of each successive generation of stars at ultraviolet wavelengths by typically 15 Myr. The predicted rest-frame far-infrared to millimeter-wave emission includes the calculation of molecular emission-line luminosities ($^{12}$CO and O$_{2}$ among other molecules) consistent with the evolving chemical abundances. Here we present details of this new model, along with the results of comparing its predictions with several high-redshift observables, namely: the ultraviolet SEDs of Lyman-limit galaxies, the high-redshift radiogalaxies 4C41.17 and 8C1435, the SCUBA sub-millimeter survey of the Hubble Deep Field (HDF), and the SEDs of intermediate-redshift elliptical galaxies. With our new reddening method, we are able to fit the spectrum of the Lyman-limit galaxy 1512-cB58, and we find an extinction of about 1.9 magnitudes at 1600 {\AA}. The model also predicts that most Lyman-limit galaxies should have a value of $\alpha$ inside that range, as it is observed. The 850 $\mu$m flux density of a typical Lyman-limit galaxy is expected to be only $\simeq 0.5$ mJy, and therefore the optical counterparts of the most luminous sub-mm sources in the HDF (or any other currently feasible sub-mm survey) are unlikely to be Lyman-break galaxies.

In this invited talk I show how the reddest and rarest galaxies at high redshift ($z \simeq 1.5$) can be used to set the stellar evolution clock. I argue that one can confidently compute the collapse redshift of these objects. This yields to a high collapse redshift ($z>6$) and therefore their age is well constrained (in all cosmologies) between 3 and 4 Gyr. I also show that this is, indeed, the age derived using a variety of synthetic stellar population models when proper statistical tools are used to analyse their observed spectral energy distribution. This allows me to conclude that all stellar population models yield to the same consistent age for these galaxies, i.e. about 3.5 Gyr and that the stellar clock is properly set. Low ages are therefore excluded with high confidence.

We show that the discrepancy between the Tully-Fisher relation and the luminosity function predicted by most phenomenological galaxy formation models is mainly due to overmerging of galaxy haloes. We have circumvented this overmerging problem, which is inherent in both the Press-Schechter formalism and dissipationless N-body simulations, by including a specific galaxy halo formation recipe into an otherwise standard N-body code. This numerical technique provides the merger trees which, together with simplified gas dynamics and star formation physics, constitute our implementation of a phenomenological galaxy formation model. Resolving the overmerging problem provides us with the means to match both the I-band Tully-Fisher relation and the B and K band luminosity functions within an EdS sCDM structure formation scenario. It also allows us to include models for chemical evolution and starbursts, which improves the match to observational data and renders the modelling more realistic. We show that the inclusion of chemical evolution into the modelling requires a significant fraction of stars to be formed in short bursts triggered by merging events.

We investigate the influence of the star formation rate on the Tully-Fisher relation. We find that a simple model which combines the empirically-determined star-formation rate with the expected properties of galaxy halos provides a remarkably good fit to the absolute magnitude-rotation speed correlation. We find that the power-law nature, its slope, normalisation and scatter, are all readily accounted for if the Universe has a low density parameter, with or without a cosmological constant and disks are assembled at $z \sim 1 - 1.5$. Moreover, this agreement is found simultaneously in 4 wavebands. An Einstein-de Sitter Universe produces disks which are too faint unless the disks are assembled at $z\sim 0.5$. The scatter in the relation is due to a combination of the expected range of spin parameters of the halos and the range of formation redshifts. The source of the scatter opens up possibilities of a better galaxy distance indicator, if spectroscopic observations of globular clusters can be used to determine the halo rotation.

It has recently been argued that single-collapse high-redshift models for elliptical galaxy formation can be rejected because they predict large numbers of very red galaxies at intermediate redshifts which are not seen in deep optical-infrared surveys. We argue, however, that this conclusion is premature since, while much effort has been invested in refining the predictions of hierarchical CDM models, only very simplistic models have been used to study the evolution of galaxies in other cosmogonies (e.g. isocurvature models). We demonstrate that the use of a more realistic multi-zone chemo-dynamical single-collapse model, yields colours at intermediate redshifts which are much bluer than inferred from the one-zone model, and indeed are comparable to those predicted by hierarchical merging despite still allowing $> 90%$ of the final stellar mass of elliptical galaxies to be formed in the first Gyr of their evolution. We, therefore, conclude that the one-zone model should be avoided to predict the colours of high-redshift galaxies and that the use of realistic multi-zone models allows the existence of ellipticals at high redshift, being their dismissal premature.

We use a model of galaxy disk formation to explore the metallicities, dust content, and neutral-gas mass density of damped Lyman-$\alpha$ (D\lya) absorbers. We find that the [Zn/H] abundance measurements of D\lya systems now available can be reproduced either by a population of low surface brightness (LSB) galaxies forming at redshifts $z > 3$, whose chemical contents evolve slowly with time and whose star formation rates are described by continuous bursts, or by high surface brightness (HSB) galaxies which form continuously over an interval of $z \sim 0.5-3$ (and no higher). Although, in reality, a mixture of galaxy types may be responsible for low-z D\lya systems, our models predict that HSB galaxies form more dust, more rapidly, than LSB galaxies, and that HSB galaxies may therefore obscure background QSOs and not give rise to D\lya lines, as suggested by other researchers. Significantly, we find that the rate at which HSB disks consume neutral gas is too fast to explain the observed evolution in the neutral gas mass density with redshift, and that the consumption of hydrogen by LSB galaxies better fits the data. This further suggests that LSB disks may dominate the D\lya population at high-redshift.

In the past 40 years a considerable effort has been focused in determining the age of the Universe at zero redshift using several stellar clocks. In this review I will describe the best theoretical methods to determine the age of the oldest Galactic Globular Clusters (GC). I will also argue that a more accurate age determination may come from passively evolving high-redshift ellipticals. In particular, I will review two new methods to determine the age of GC. These two methods are more accurate than the classical isochrone fitting technique. The first method is based on the morphology of the horizontal branch and is independent of the distance modulus of the globular cluster. The second method uses a careful binning of the stellar luminosity function which determines simultaneously the distance and age of the GC. It is found that the oldest GCs have an age of $13.5 \pm 2$ Gyr. The absolute minimum age for the oldest GCs is 10.5 Gyr and the maximum is 16.0 Gyr (with 99% confidence). Therefore, an Einstein-De Sitter Universe ($\Omega=1$) is not totally ruled out if the Hubble constant is about $65 \pm 10$ km s$^{-1}$ Mpc$^{-1}$. On the other hand, the newly discovered red elliptical 53W069 ($z=1.43$) provides an stronger constraint since its minimun age is 3.2 Gyr, thus ruling out an Einstein-De Sitter Universe unless the Hubble constant is $\leq 45 \pm 10$ km s$^{-1}$. Using 53W069 we find an age at $z=0$ of $13 \pm 2$ Gyr, in excellent agreement with the GC determination.

Accurate synthetic models of stellar populations are constructed and used in evolutionary models of stellar populations in forming galaxies. Following their formation, the late type galaxies are assumed to follow the Schmidt law for star formation, while early type galaxies are normalized to the present-day fundamental plane relations assumed to mimic the metallicity variations along their luminosity sequence. We then compute predictions of these models for the observational data at early epochs for various cosmological parameters $\Omega, \Omega_\Lambda$ and $H_0$. We find good match to the metallicity data from the damped $L_\alpha$ systems and the evolution of the luminosity density out to $z\simeq 1$. Likewise, our models provide good fits for low values of $\Omega$ to the deep number counts of galaxies in all bands where data is available; this is done without assuming existence of extra populations of galaxies at high $z$. Our models also match the data on the redshift distribution of galaxy counts in $B$ and $K$ bands. We compute the predicted mean levels and angular distribution of the cosmic infrared background produced from the early evolution of galaxies. The predicted fluxes and fluctuations are still below the current observational limits, but not by a large factor. Finally, we find that the recent detection of the diffuse extragalactic light in the visible bands requires for our models high redshift of galaxy formation, $z_f \geq$(3-4); otherwise the produced flux of the extragalactic light at optical bands exceeds the current observational limits.

We investigate in detail the hypothesis that low surface brightness galaxies (LSB) differ from ordinary galaxies simply because they form in halos with large spin parameters. We compute star formation rates using the Schmidt law, assuming the same gas infall dependence on surface density as used in models of the Milky Way. We build stellar population models, predicting colours, spectra, and chemical abundances. We compare our predictions with observed values of metallicity and colours for LSB galaxies and find excellent agreement with all observables. In particular, integrated colours, colour gradients, surface brightness and metallicity match very well to the observed values of LSBs for models with ages larger than 7 Gyr and high values ($\lambda > 0.05$) for the spin parameter of the halos. We also compute the global star formation rate (SFR) in the Universe due to LSBs and show that it has a flatter evolution with redshift than the corresponding SFR for normal discs. We furthermore compare the evolution in redshift of $[Zn/H]$ for our models to those observed in Damped Lyman $\alpha$ systems by \scite{Pettini+97} and show that Damped Lyman $\alpha$ systems abundances are consistent with the predicted abundances at different radii for LSBs. Finally, we show how the required late redshift of collapse of the halo may constrain the power spectrum of fluctuations.

We have discovered a population of extremely red galaxies at $z\simeq 1.5$ which have apparent stellar ages of $\gs 3$ Gyr, based on detailed spectroscopy in the rest-frame ultraviolet. In order for galaxies to have existed at the high collapse redshifts indicated by these ages, there must be a minimum level of power in the density fluctuation spectrum on galaxy scales. This paper compares the required power with that inferred from other high-redshift populations. If the collapse redshifts for the old red galaxies are in the range $z_c\simeq 6$ -- 8, there is general agreement between the various tracers on the required inhomogeneity on 1-Mpc scales. This level of small-scale power requires the Lyman-limit galaxies to be approximately $\nu\simeq 3.0$ fluctuations, implying a very large bias parameter $b\simeq 6$. The high collapse redshifts of the red galaxies as deduced from gravitational collapse provides independent support for the ages estimated from their stellar populations. Such early-forming galaxies are rare, and their contribution to the cosmological stellar density is consistent with an extrapolation to higher redshifts of the star-formation rate measured at $z<5$; there is no evidence for a general era of spheroid formation at extreme redshifts.

We have used the Hipparcos satellite colour magnitude diagram to determine the age of the Galactic disc. We first measure the metallicities of Hipparcos giants using DDO photometry (H{\o}g & Flynn 1997). With accurate metallicities and distances we compute the minimum age of the Galactic disc as $11\pm1$ Gyr using the colour of the reddest clump stars as a function of metallicity. In conjunction with the new ages derived for Globular Clusters using the same method (Jimenez et al 1996) we show that any delay between the formation of the halo and the disc was only 2-3 Gyr. Since disc galaxies may have been detected to redshifts as high as $z=3.15$ (Lu, Sargent and Barlow 1997) an age for our own Galactic disc of 11 Gyr favours open or non-zero $\Lambda$ cosmological models. We show that a Reimers mass-loss law is sufficient to explain the morphology of the red clump.

Gravitational lensing studies suggest that the Universe may contain a population of dark galaxies; we investigate this intriguing possibility and propose a mechanism to explain their nature. In this mechanism a dark galaxy is formed with a low density disk in a dark halo of high spin parameter; such galaxies can have surface densities below the critical Toomre value for instabilities to develop, and following Kennicutt's work we expect these galaxies to have low star formation rates. The only stellar component of the galaxies is a halo system, formed during the collapse of the proto-galactic cloud. We compute synthetic stellar population models and show that, at a redshift $z=0.5$, such galaxies have apparent magnitudes $B \simeq 28, R \simeq 26$ and $I \simeq 25$, and could be unveiled by deep searches with the Hubble Space Telescope. Dark galaxies have an initial short blue phase and then become essentially invisible, therefore they may account for the blue population of galaxies at high redshift. We find a strong mass-dependence in the fraction of dark galaxies, and predict that spiral galaxies will not be found in halos with masses less than about $10^9 M_\odot$, if $\Omega=1$. Above about $10^{12}M_\odot$, all halos can produce luminous disks. The mass-dependence of the galaxy-formation efficiency introduces the possibility of `spin bias' -- luminous galaxies being associated preferentially with strongly-clustered high-mass halos. A further prediction is that the slope of the faint-end luminosity function for galaxies will be flatter than the associated halo mass function.

The weak radio source LBDS 53W091 is associated with a very faint (R approx. 24.5) red (R-K approx. 5.8) galaxy. Long spectroscopic integrations with the W. M. Keck telescope have provided an absorption-line redshift, z=1.552 +/- 0.002. The galaxy has a rest frame ultraviolet spectrum very similar to that of an F6V star, and a single-burst old stellar population that matches the IR colors, the optical energy distribution and the spectral discontinuities has a minimum age of 3.5 Gyr. We present detailed population synthesis analyses of the observed spectrum in order to estimate the time since the last major epoch of star formation. We discuss the discrepancies in these estimates resulting from using different models, subjecting the UV spectrum of M32 to the same tests as a measure of robustness of these techniques. The models most consistent with the data tend to yield ages at z=1.55 of >3.5Gyr, similar to that inferred for the intermediate-age population in M32. Depending upon the assumed Hubble constant and the value of Omega_0, only certain cosmological expansion times are consistent with the age of LBDS 53W091; in particular, for Omega_0=1, only models with H_0 < 45 km/s/Mpc are permitted. For H_0=50 km/s/Mpc and Omega_0=0.2, we derive a formation redshift, z_f >= 5.

In this article I review the main methods for determining the age of the Universe. I describe how to determine the age of the oldest known systems at z=0, the system of galactic globular clusters, using different techniques. I also describe how to date the Universe using the decay of radioactive elements (Cosmochronology). Finally, I focus on how to determine the age of the Universe at different redshifts and specially the age of radio-quiet galaxies at high redshift. I finish by arguing that the most probable age for the Universe is $14 \pm 2$ Gyr.

We present new age and distance determinations for the Galactic Globular Clusters M55 and M5, using the luminosity function method (Jimenez & Padoan 1996, Padoan & Jimenez 1997). We find an age of $11.8 \pm 1.5$ Gyr for M55 and $11.1 \pm 0.7$ Gyr for M5. This confirms previous results (Jimenez et al. 1996, Sandquist et al. 1996) and allows to conclude that the oldest stars in the Universe are not older than 14 Gyr. We also find $m-M=14.13 \pm 0.11$ for M55, and $m-M=14.49 \pm 0.06$ for M5. These values agree with the ones obtained using the tip of the red giant branch (Jimenez et al. 1996) and the sub-dwarf fitting method (Sandquist et al. 1996).

The recently discovered high redshift galaxy, 53W091, with accurate age measurements (Dunlop et al 1996) provides a measure of the small-scale power of the primordial density field and, as we show, a crucial test of the inflation-inspired models. It allows $\Omega$=1 cosmologies only for low values of $H_0$, but then pushes formation of that galaxy to redshift much greater than allowed for by the cold-dark-matter density field. Cold-dark-matter (CDM) models with cosmological constant ($\Lambda$) and low $\Omega$ would decrease the redshift at which this galaxy has collapsed. However, in CDM models decreasing $\Omega$ suppresses the small scale power in the density field and this effect turns out to be dominant. We estimate the mass of the galaxy and show that it represents a very rare and unlikely event in the density field of such models. Similar problems would occur in other modifications of the CDM cosmogonies.

The epoch of galaxy formation provides an important additional test of cosmological theories. Cold-dark-matter (CDM) models with cosmological constant ($\Lambda$) are designed to account for the observed excess power in galaxy distribution, but at the same time suppress the power on small scales pushing galaxy formation to recent epochs. We point out that the recently discovered high redshift galaxy, 53W091, with accurate age measurements (Dunlop et al 1996) provides a particularly important test of these models. In the flat $\Lambda$-dominated Universe, the redshift of formation of 53W091 decreases with decreasing $\Omega$. However, in the modified CDM models decreasing $\Omega$ suppresses the small scale power in the density field and this effect turns out to be dominant. We estimate the mass of the galaxy and show that it represents a very rare and unlikely event in the density field of such models. Similar problems would occur in other modifications of the CDM cosmogonies.

We reconsider the problem of the age of the stellar discs of late-type Low Surface Brightness (LSB) galaxies by making use of a new IMF recently derived from numerical fluid dynamical simulations (Padoan, Nordlund and Jones, 1997). While a Miller-Scalo IMF cannot adequately describe the photometric properties of LSBs, when we apply the new Padoan et al. (1997) IMF to a simple exponential disc model with parameters appropriate to LSBs, we get excellent fits of the colors and color gradients. We then conclude that: a) The star formation history of LSB disc galaxies can be described by an initial burst of a few times $10^7$ yr followed by a quiescent period with only sporadic star formation; b) LSBs' discs are not young. The age of the LSB disc galaxies for which colors have been measured are all larger than about 10 Gyr$.

Using a new physical model for star formation (Padoan 1995) we have tested the possibility that globular clusters (GCs) are formed from primordial mass fluctuations, whose mass scale ($10^8$ - $10^9$ M$_{\odot}$) is selected out of a CDM spectrum by the mechanism of non-equilibrium formation of $H_2$. We show that such clouds are able to convert about 0.003 of their total mass into a bound system (GC) and about 0.02 into halo stars. The metal enriched gas is dispersed away from the GC by supernova explosions and forms the galactic disk. These mass ratios between GCs, halo and disk depend on the predicted IMF which is a consequence of the universal statistics of fluid turbulence. They also depend on the ratio of baryonic over non-baryonic mass ,$X_b$, and are comparable with the values observed in typical spiral galaxies for $X_b \approx 0.1-0.2$. The computed mass and radius for a GC ( $5\times 10^5$ M$_{\odot}$ and 30 pc) are in good agreement with the average values in the Galaxy. The model predicts an exponential cut off in the stellar IMF below 0.1 M$_{\odot}$ in GCs and 0.6 M$_{\odot}$ in the halo. The quite massive star formation in primordial clouds leads to a large number of supernovae and to a high blue luminosity during the first two Gyr of the life of every galaxy.

We extend our previous method to determine globular cluster ages using the luminosity function (Jimenez \& Padoan 1996). We show that the luminosity function depends on both age and distance modulus and that it is possible to distinguish between the two. This method provides at the same time independent determinations of distance and age of a GC by simply counting the number of stars found inside specified luminosity bins. The main uncertainties in other traditional methods for determining GCs ages are absent (e.g. mixing length, color-$T_{\rm eff}$ calibration, morphology of the color-magnitude diagram ). The distance modulus is the biggest uncertainty in determining the age of GCs. Here we show that the age can be determined with small uncertainty for any value of distance modulus using the LF and that the LF allows a determination of the distance modulus itself. This is explained by the fact that the luminosity function is affected by a change in distance-modulus in a way that is different from its time evolution. If GC stellar counts with statistical errors not larger than $3\%$ are available, the age can be determined with an uncertainty of about 0.4 Gyr (independent of distance modulus, mixing length and color calibration) and the distance modulus with an uncertainty of about 0.04 mag.

We have applied a new method to analyze the horizontal branch (HB) morphology in relation to the distribution of stars near the red giant branch (RGB) tip for the globular clusters M22, M5, M68, M107, M72, M92, M3 and 47 Tuc. This new method permits determination of the cluster ages to greater accuracy than conventional isochrone fitting. Using the method in conjunction with our new high-quality photometric data for RGB and HB stars in the first five of these clusters, we discuss the origins of the spread in color on the HB and its relation to the `second parameter' problem. The oldest clusters in our sample are found to have relatively low ages ($13.5 \pm 2$ Gyr). A 1$\sigma$ uncertainty in each of the parameters of mass and helium content combined with the effects of helium diffusion gives a lower limit for the age of the oldest clusters of 9.7 Gyr.

We present a new method to compute stellar ages in Globular Clusters (GC) that is ten times more precise than the traditional isochrone fitting procedure. The method relies on accurate stellar evolutionary tracks and on photometry for GCs complete down to the main sequence, and it is based on counting number of stars in two different regions of the CMD: the red giant branch and the main-sequence. We have applied this method to the globular cluster M68 and found an age of 16.4$\pm$0.2 Gyr for $(m-M)_V=15.3$. This new method reduces the error associated to the uncertainty in the distance modulus by a factor of two, the error due to the choice of the value for the mixing length parameter to almost zero and the error due to the colour-$T_{\rm eff}$ transformation to zero.

We present here results of stellar evolution calculations that include the latest advances in radiative opacities and neutrino cooling, and discuss on the basis of these models how the internal stellar structure responds to mass-loss from the stellar surface. This problem has particular importance for the development of semi-analytical algorithms for effi cient calculation of synthetic stellar populations with realistic (and hence complex) mass-loss scenarios. We therefore compare our numerical results with test calculations based on a semi-analytical stellar evolution method developed by us. Although small, but important, differences between results from the two methods are revealed, the evolutionary tracks in the HR-diagram predicted by the two approaches are almost identical.