The Baryon Oscillation Spectroscopic Survey of SDSS-III has provided an unprecedentedly large sample of Damped \lya systems (DLAs), the largest repositories of neutral hydrogen in the Universe. This DLA sample has been used to determine the DLA bias factor from their cross-correlation with the \lya forest absorption in \cite{FontRibera2012,Perez2018}, showing that DLAs are associated with relatively massive halos. However, the low resolution and signal-to-noise of BOSS spectra do not allow precise measurements of the DLA metal lines. We define a metal strength parameter, $S$, based on combining equivalent widths of 17 metal lines, that can be measured with an optimal signal-to-noise ratio for individual DLAs in BOSS spectra, allowing for the classification of the DLA population into subgroups of different $S$. We present the distribution of this DLA metal strength and the dependence of its mean value on $N_{\rm HI}$ and redshift. We search for systematic effects and variations in the catalogue purity by examining the dependence of the $S$ distribution on the spectral signal-to-noise and the estimated error on $S$. A catalogue of DLAs with measured equivalent widths for the selected 17 metal lines and the value of $S$ is made publicly available, which will be used to measure the dependence of the DLA bias factor on the $S$ parameter. The relation of the metal strength on the gas metal abundances and velocity dispersion can be constrained by studying the stacked metal absorption spectra of DLAs as a function of $S$, allowing for future determinations of the dependence of the bias factor on the metallicity and velocity dispersion of DLAs.

In this work, we investigate the abundance and distribution of metals in the intergalactic medium (IGM) at $\langle z \rangle \simeq 2.8$ through the analysis of an ultra-high signal-to-noise ratio UVES spectrum of the quasar HE0940-1050. In the CIV forest, our deep spectrum is sensitive at $3\,\sigma$ to lines with column density down to $\log N_{\rm CIV} \simeq 11.4$ and in 60 per cent of the considered redshift range down to $\simeq11.1$. In our sample, all HI lines with $\log N_{\rm HI} \ge 14.8$ show an associated CIV absorption. In the range $14.0 \le \log N_{\rm HI} <14.8$, 43 per cent of HI lines has an associated CIV absorption. At $\log N_{\rm HI} < 14.0$, the detection rates drop to $<10$ per cent, possibly due to our sensitivity limits and not to an actual variation of the gas abundance properties. In the range $\log N_{\rm HI} \ge 14$, we observe a fraction of HI lines with detected CIV a factor of 2 larger than the fraction of HI lines lying in the circum-galactic medium (CGM) of relatively bright Lyman-break galaxies hosted by dark matter haloes with $\langle M\rangle \sim10^{12}$ M$_{\odot}$. The comparison of our results with the output of a grid of photoionization models and of two cosmological simulations implies that the volume filling factor of the IGM gas enriched to a metallicity $\log Z/Z_{\odot} \ge -3$ should be of the order of $\sim 10-13$ percent. In conclusion, our results favour a scenario in which metals are found also outside the CGM of bright star-forming galaxies, possibly due to pollution by lower mass objects and/or to an early enrichment by the first sources.

We measure the sum of the neutrino particle masses using the three-dimensional galaxy power spectrum of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) CMASS galaxy sample. Combined with the cosmic microwave background (CMB), supernova (SN) and additional baryonic acoustic oscillation (BAO) data, we find upper 95 percent confidence limits of the neutrino mass $\Sigma m_{\nu}<0.340$ eV within a flat $\Lambda$CDM background, and $\Sigma m_{\nu}<0.821$ eV, assuming a more general background cosmological model. The number of neutrino species is measured to be $N_{\rm eff}=4.308\pm0.794$ and $N_{\rm eff}=4.032^{+0.870}_{-0.894}$ for these two cases respectively. We study and quantify the effect of several factors on the neutrino measurements, including the galaxy power spectrum bias model, the effect of redshift-space distortion, the cutoff scale of the power spectrum, and the choice of additional data. The impact of neutrinos with unknown masses on other cosmological parameter measurements is investigated. The fractional matter density and the Hubble parameter are measured to be $\Omega_M=0.2796\pm0.0097$, $H_0=69.72^{+0.90}_{-0.91}$ km/s/Mpc (flat $\Lambda$CDM) and $\Omega_M=0.2798^{+0.0132}_{-0.0136}$, $H_0=73.78^{+3.16}_{-3.17}$ km/s/Mpc (more general background model). Based on a Chevallier-Polarski-Linder (CPL) parametrisation of the equation-of-state $w$ of dark energy, we find that $w=-1$ is consistent with observations, even allowing for neutrinos. Similarly, the curvature \Omega_K and the running of the spectral index $\alpha_s$ are both consistent with zero. The tensor-to-scaler ratio is constrained down to $r<0.198$ (95 percent CL, flat $\Lambda$ CDM) and $r<0.440$ (95 percent CL, more general background model).

The 511 keV electron-positron annihilation line, most recently characterized by the INTEGRAL/SPI experiment, is highly concentrated towards the Galactic centre. Its origin remains unknown despite decades of scrutiny. We propose a novel scenario in which known extragalactic positron sources such as radio jets of active galactic nuclei (AGN) fill the intergalactic medium with MeV e+e- pairs, which are then accreted into the Milky Way. We show that interpreting the diffuse cosmic radio background (CRB) as arising from radio sources with characteristics similar to the observed cores and radio lobes in powerful AGN jets suggests that the intergalactic positron-to-electron ratio could be as high as 10^{-5}, although this can be decreased if the CRB is not all produced by pairs and if not all positrons escape to the intergalactic medium. Assuming an accretion rate of one solar mass per year of matter into the Milky Way, a positron-to-electron ratio of only 10^{-7} is already enough to account for much of the 511 keV emission of the Galaxy. A simple spherical accretion model predicts an emission profile highly peaked in the central bulge, consistent with INTEGRAL observations. However, a realistic model of accretion with angular momentum would likely imply a more extended emission over the disk, with uncertainties depending on the magnetic field structure and turbulence in the galactic halo.

We present the BOSS Lyman-alpha (Lya) Forest Sample from SDSS Data Release 9, comprising 54,468 quasar spectra with zqso > 2.15 suitable for Lya forest analysis. This data set probes the intergalactic medium with absorption redshifts 2.0 < z_alpha < 5.7 over an area of 3275 square degrees, and encompasses an approximate comoving volume of 20 h^-3 Gpc^3. With each spectrum, we have included several products designed to aid in Lya forest analysis: improved sky masks that flag pixels where data may be unreliable, corrections for known biases in the pipeline estimated noise, masks for the cores of damped Lya systems and corrections for their wings, and estimates of the unabsorbed continua so that the observed flux can be converted to a fractional transmission. The continua are derived using a principal component fit to the quasar spectrum redwards of restframe Lya (lambda > 1216 Ang), extrapolated into the forest region and normalized by a linear function to fit the expected evolution of the Lya forest mean-flux. The estimated continuum errors are ~5% rms. We also discuss possible systematics arising from uncertain spectrophotometry and artifacts in the flux calibration; global corrections for the latter are provided. Our sample provides a convenient starting point for users to analyze clustering in BOSS Lya forest data, and it provides a fiducial data set that can be used to compare results from different analyses of baryon acoustic oscillations in the Lya forest. The full data set is available from the SDSS-III DR9 web site.

We propose that the cloud moving on a highly eccentric orbit near the central black hole in our Galaxy, reported by Gillessen et al., is formed by a photoevaporation wind originating in a disk around a star that is tidally perturbed and shocked at every peribothron passage. The disk is proposed to have formed when a stellar black hole flew by the star, tidally disrupted its envelope, and placed the star on its present orbit with some of the tidal debris forming a disk. A disrupting encounter at the location of the observed cloud is most likely to be caused by a stellar black hole because of the expected dynamical mass segregation; the rate of these disk-forming encounters may be as high as $\sim 10^{-6}$ per year. The star should also be spun up by the encounter, so the disk may subsequently expand by absorbing angular momentum from the star. Once the disk expands up to the tidal truncation radius, the tidal perturbation of the outer disk edge at every peribothron may place gas streams on larger orbits which can give rise to a photoevaporation wind that forms the cloud at every orbit. This model predicts that, after the cloud is disrupted at the next peribothron passage in 2013, a smaller unresolved cloud will gradually grow around the star on the same present orbit. An increased infrared luminosity from the disk may also be detectable when the peribothron is reached. We also note that this model revives the encounter theory for planet formation.

We measure the quasar two-point correlation function over the redshift range 2.2<z<2.8 using data from the Baryon Oscillation Spectroscopic Survey. We use a homogeneous subset of the data consisting of 27,129 quasars with spectroscopic redshifts---by far the largest such sample used for clustering measurements at these redshifts to date. The sample covers 3,600 square degrees, corresponding to a comoving volume of 9.7(Gpc/h)^3 assuming a fiducial LambdaCDM cosmology, and it has a median absolute i-band magnitude of -26, k-corrected to z=2. After accounting for redshift errors we find that the redshift space correlation function is fit well by a power-law of slope -2 and amplitude s_0=(9.7\pm 0.5)Mpc/h over the range 3<s<25Mpc/h. The projected correlation function, which integrates out the effects of peculiar velocities and redshift errors, is fit well by a power-law of slope -1 and r_0=(8.4\pm 0.6)Mpc/h over the range 4<R<16Mpc/h. There is no evidence for strong luminosity or redshift dependence to the clustering amplitude, in part because of the limited dynamic range in our sample. Our results are consistent with, but more precise than, previous measurements at similar redshifts. Our measurement of the quasar clustering amplitude implies a bias factor of b~3.5 for our quasar sample. We compare the data to models to constrain the manner in which quasars occupy dark matter halos at z~2.4 and infer that such quasars inhabit halos with a characteristic mass of <M>~10^{12}Msun/h with a duty cycle for the quasar activity of 1 per cent.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

Lyman-alpha (Lya) photons that escape the interstellar medium of star-forming galaxies may be resonantly scattered by neutral hydrogen atoms in the circumgalactic and intergalactic media, thereby increasing the angular extent of the galaxy's Lya emission. We present predictions of this extended, low surface brightness Lya emission based on radiative transfer modeling in a cosmological reionization simulation. The extended emission can be detected from stacked narrowband images of Lya emitters (LAEs) or of Lyman break galaxies (LBGs). Its average surface brightness profile has a central cusp, then flattens to an approximate plateau beginning at an inner characteristic scale below ~0.2 Mpc (comoving), then steepens again beyond an outer characteristic scale of ~1 Mpc. The inner scale marks the transition from scattered light of the central source to emission from clustered sources, while the outer scale marks the spatial extent of scattered emission from these clustered sources. Both scales tend to increase with halo mass, UV luminosity, and observed Lya luminosity. The extended emission predicted by our simulation is already within reach of deep narrowband photometry using large ground-based telescopes. Such observations would test radiative transfer models of emission from LAEs and LBGs, and they would open a new window on the circumgalactic environment of high-redshift star-forming galaxies.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

The SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), a five-year spectroscopic survey of 10,000 deg^2, achieved first light in late 2009. One of the key goals of BOSS is to measure the signature of baryon acoustic oscillations in the distribution of Ly-alpha absorption from the spectra of a sample of ~150,000 z>2.2 quasars. Along with measuring the angular diameter distance at z\approx2.5, BOSS will provide the first direct measurement of the expansion rate of the Universe at z > 2. One of the biggest challenges in achieving this goal is an efficient target selection algorithm for quasars over 2.2 < z < 3.5, where their colors overlap those of stars. During the first year of the BOSS survey, quasar target selection methods were developed and tested to meet the requirement of delivering at least 15 quasars deg^-2 in this redshift range, out of 40 targets deg^-2. To achieve these surface densities, the magnitude limit of the quasar targets was set at g <= 22.0 or r<=21.85. While detection of the BAO signature in the Ly-alpha absorption in quasar spectra does not require a uniform target selection, many other astrophysical studies do. We therefore defined a uniformly-selected subsample of 20 targets deg^-2, for which the selection efficiency is just over 50%. This "CORE" subsample will be fixed for Years Two through Five of the survey. In this paper we describe the evolution and implementation of the BOSS quasar target selection algorithms during the first two years of BOSS operations. We analyze the spectra obtained during the first year. 11,263 new z>2.2 quasars were spectroscopically confirmed by BOSS. Our current algorithms select an average of 15 z > 2.2 quasars deg^-2 from 40 targets deg^-2 using single-epoch SDSS imaging. Multi-epoch optical data and data at other wavelengths can further improve the efficiency and completeness of BOSS quasar target selection. [Abridged]

We study the clustering properties of z~5.7 Lyman-alpha emitters (LAEs) in a cosmological reionization simulation with a full Lya radiative transfer calculation. Lya radiative transfer substantially modifies the intrinsic Lya emission properties, compared to observed ones, depending on the density and velocity structure environment around the Lya emitting galaxy. This environment-dependent Lya selection introduces new features in LAE clustering, suppressing (enhancing) the line-of-sight (transverse) density fluctuations and giving rise to scale-dependent galaxy bias. In real space, the contours of the three-dimensional two-point correlation function of LAEs appear to be prominently elongated along the line of sight on large scales, an effect that is opposite to and much stronger than the linear redshift-space distortion effect. The projected two-point correlation function is greatly enhanced in amplitude by a factor of up to a few, compared to the case without the environment dependent selection effect. The new features in LAE clustering can be understood with a simple, physically motivated model, where Lya selection depends on matter density, velocity, and their gradients. We discuss the implications and consequences of the effects on galaxy clustering from Lya selection in interpreting clustering measurements and in constraining cosmology and reionization from LAEs.

We combine a cosmological reionization simulation with box size of 100Mpc/h on a side and a Monte Carlo Lyman-alpha (Lya) radiative transfer code to model Lyman Alpha Emitters (LAEs) at z~5.7. The model introduces Lya radiative transfer as the single factor for transforming the intrinsic Lya emission properties into the observed ones. Spatial diffusion of Lya photons from radiative transfer results in extended Lya emission and only the central part with high surface brightness can be observed. Because of radiative transfer, the appearance of LAEs depends on density and velocity structures in circumgalactic and intergalactic media as well as the viewing angle, which leads to a broad distribution of apparent (observed) Lya luminosity for a given intrinsic Lya luminosity. Radiative transfer also causes frequency diffusion of Lya photons. The resultant Lya line is asymmetric with a red tail. The peak of the Lya line shifts towards longer wavelength and the shift is anti-correlated with the apparent to intrinsic Lya luminosity ratio. The simple radiative transfer model provides a new framework for studying LAEs. It is able to explain an array of observed properties of z~5.7 LAEs in Ouchi et al. (2008), producing Lya spectra, morphology, and apparent Lya luminosity function (LF) similar to those seen in observation. The broad distribution of apparent Lya luminosity at fixed UV luminosity provides a natural explanation for the observed UV LF, especially the turnover towards the low luminosity end. The model also reproduces the observed distribution of Lya equivalent width (EW) and explains the deficit of UV bright, high EW sources. Because of the broad distribution of the apparent to intrinsic Lya luminosity ratio, the model predicts effective duty cycles and Lya escape fractions for LAEs.

There is a remarkable synergy between requirements for Dark Energy probes by cosmic shear measurements and planet hunting by microlensing. Employing weak and strong gravitational lensing to trace and detect the distribution of matter on cosmic and Galactic scales, but as well as to the very small scales of exoplanets is a unique meeting point from cosmology to exoplanets. It will use gravity as the tool to explore the full range of masses not accessible by any other means. EUCLID is a 1.2m telescope with optical and IR wide field imagers and slitless spectroscopy, proposed to ESA Cosmic Vision to probe for Dark Energy, Baryonic acoustic oscillation, galaxy evolution, and an exoplanet hunt via microlensing. A 3 months microlensing program will already efficiently probe for planets down to the mass of Mars at the snow line, for free floating terrestrial or gaseous planets and habitable super Earth. A 12+ months survey would give a census on habitable Earth planets around solar like stars. This is the perfect complement to the statistics that will be provided by the KEPLER satellite, and these missions combined will provide a full census of extrasolar planets from hot, warm, habitable, frozen to free floating.

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.

Gaztanaga et al. have recently claimed to measure the Baryon Acoustic Oscillation (BAO) scale in the radial direction from the publicly available SDSS DR6 data. They focus on the correlation function of Luminous Red Galaxies (LRG) close to the line-of-sight direction to find a feature that they identify as the BAO peak, arguing that a magnification bias effect from gravitational lensing increases the amplitude of the BAO peak, facilitating its detection. In this Comment, we clarify that lensing has a negligible impact on the measurement of the BAO peak, and that the interpretation by Gaztanaga et al. is incorrect. The feature they identify in the LRG correlation function near the line-of-sight cannot be explained by any known physical effect and is in fact consistent with noise.

We predict the 21-cm signature of the first metal-free stars. The soft X-rays emitted by these stars penetrate the atomic medium around their host halos, generating Lyman alpha photons that couple the spin and kinetic temperatures. These creates a region we call the Lyman alpha sphere, visible in 21-cm against the CMB, which is much larger than the HII region produced by the same star. The spin and kinetic temperatures are strongly coupled before the X-rays can substantially heat the medium, implying that a strong 21-cm absorption signal from the adiabatically cooled gas in Hubble expansion around the star is expected when the medium has not been heated previously. A central region of emission from the gas heated by the soft X-rays is also present although with a weaker signal than the absorption. The Lyman alpha sphere is a universal signature that should be observed around any first star illuminating its vicinity for the first time. The 21-cm radial profile of the Lyman alpha sphere can be calculated as a function of the luminosity, spectrum and age of the star. For a star of a few hundred solar masses and zero metallicity (as expected for the first stars), the physical radius of the Lyman alpha sphere can reach tens of kiloparsecs. The first metal-free stars should be strongly clustered because of high cosmic biasing; this implies that the regions producing a 21-cm absorption signal may contain more than one star and will generally be irregular and not spherical, because of the complex distribution of the gas. We discuss the feasiblity of detecting these Lyman alpha spheres, which would be present at redshifts $z\sim 30$ in the Cold Dark Matter model. Their observation would represent a direct proof of the detection of a first star.

Recent observations support the idea that nuclear black holes grew by gas accretion while shining as luminous quasars at high redshift, and they establish a relation of the black hole mass with the host galaxy's spheroidal stellar system. We develop an analytic model to calculate the expected impact of mergers on the masses of black holes in massive clusters of galaxies. We use the extended Press-Schechter formalism to generate Monte Carlo merger histories of halos with a mass 10^{15} h^{-1} Msun. We assume that the black hole mass function at z=2 is similar to that inferred from observations at z=0 (since quasar activity declines markedly at z<2), and we assign black holes to the progenitor halos assuming a monotonic relation between halo mass and black hole mass. We follow the dynamical evolution of subhalos within larger halos, allowing for tidal stripping, the loss of orbital energy by dynamical friction, and random orbital perturbations in gravitational encounters with subhalos, and we assume that mergers of subhalos are followed by mergers of their central black holes. Our analytic model reproduces numerical estimates of the subhalo mass function. We find that the most massive black holes in massive clusters typically grow by a factor ~ 2 by mergers after gas accretion has stopped. In our ten realizations of 10^{15} h^{-1} Msun clusters, the highest initial (z=2) black hole masses are 5-7 x 10^9 Msun, but four of the clusters contain black holes in the range 1-1.5 x 10^{10} Msun at z=0. Satellite galaxies may host black holes whose mass is comparable to, or even greater than, that of the central galaxy. Thus, black hole mergers can significantly extend the very high end of the black hole mass function.

Two brightest cluster galaxies (BCGs), namely Abell 1836-BCG and Abell 3565-BCG, were observed with the Advanced Camera for Surveys (ACS) and the Space Telescope Imaging Spectrograph (STIS) on board the Hubble Space Telescope. By modeling the available photometric and kinematic data, it resulted that the mass of Abell 1836-BCG and Abell 3565-BCG are M_bh=4.8(+0.8,-0.7)x10^9 M_sun and M_bh=1.3(+0.3,-0.4)x10^9 M_sun at 1 sigma confidence level, respectively.

Astronomical observations have shown that the expansion of the universe is at present accelerating, in a way consistent with the presence of a positive cosmological constant. This is a major puzzle, because we do not understand: why the cosmological constant is so small; why, being so small, it is not exactly zero; and why it has precisely the value it must have to make the expansion start accelerating just at the epoch when we are observing the universe. We present a new model of cosmology, which we call the stryngbohtyk model, that solves all these problems and predicts exactly the value that the cosmological constant must have. The predicted value agrees with the observed one within the measurement error. We show that in the stryngbohtyk model, the fact the cosmological constant starts being important at the present epoch is not a coincidence at all, but a necessity implied by our origin in a planet orbiting a star that formed when the age of the universe was of the same order as the lifetime of the star.

We examine the distributions of eccentricity and host star metallicity of exoplanets as a function of their mass. Planets with M sin i >~ 4 M_J have an eccentricity distribution consistent with that of binary stars, while planets with M sin i <~ 4 M_J are less eccentric than binary stars and more massive planets. In addition, host star metallicities decrease with planet mass. The statistical significance of both of these trends is only marginal with the present sample of exoplanets. To account for these trends, we hypothesize that there are two populations of gaseous planets: the low-mass population forms by gas accretion onto a rock-ice core in a circumstellar disk and is more abundant at high metalliticities, and the high-mass population forms directly by fragmentation of a pre-stellar cloud. Planets of the first population form in initially circular orbits and grow their eccentricities later, and may have a mass upper limit from the total mass of the disk that can be accreted by the core. The second population may have a mass lower limit resulting from opacity-limited fragmentation. This would roughly divide the two populations in mass, although they would likely overlap over some mass range. If most objects in the second population form before the pre-stellar cloud becomes highly opaque, they would have to be initially located in orbits larger than ~30 AU, and would need to migrate to the much smaller orbits in which they are observed. The higher mean orbital eccentricity of the second population might be caused by the larger required intervals of radial migration, and the brown dwarf desert might be due to the inability of high-mass brown dwarfs to migrate inwards sufficiently in radius.

We observed a sample of three Brightest Cluster Galaxies (BCGs), Abell 1836-BCG, Abell 2052-BCG, and Abell 3565-BCG, with the Advanced Camera for Surveys (ACS) and the Imaging Spectrograph (STIS) on board the Space Telescope. For each target galaxy we obtained high-resolution spectroscopy of the Halpha and [NII]lambda6583 emission lines at three slit positions, to measure the central ionized-gas kinematics. ACS images in three different filters (F435W, F625W, and FR656N) have been used to determine the optical depth of the dust, stellar mass distribution near the nucleus, and intensity map. We present supermassive black hole (SBH) mass estimates for two galaxies which show regular rotation curves and strong central velocity gradients, and an upper limit on the SBH mass of the third one. For the SBHs of Abell 1836-BCG and Abell 3565-BCG, we derived M_bh=4.8(-0.7,+0.8)10**9 M_sun and M_bh=1.3(-0.4,+0.3)10**9 M_sun at 1 sigma confidence level, respectively. For the SBH of Abell 2052-BCG, we found M_bh < 7.3 10**9 M_sun.

Precise proper motion measurements (sigma_mu ~ 10 mkas/yr) of the recently discovered hyper-velocity star (HVS) SDSS J090745.0+024507 would yield significant constraints on the axis ratios and orientation of a triaxial model for the Galactic halo. Triaxiality of dark matter halos is predicted by Cold Dark Matter models of galaxy formation and may be used to probe the nature of dark matter. However, unless the distance to this star is determined to better than 10%, these constraints suffer from one-dimensional degeneracies, which we quantify. We show how proper motion measurements of several HVSs could simultaneously resolve the distance degeneracies of all such stars and produce a detailed picture of the triaxial halo. Additional HVSs may be found from radial velocity surveys or from parallax/proper-motion data derived from GAIA. High-precision proper-motion measurements of these stars using the Space Interferometry Mission (SIM PlanetQuest) would substantially tighten the constraints they yield on the Galactic potential.

We estimate the maximum contribution to reionization from the first generation of massive stars, with zero metallicity, under the assumption that one of these stars forms with a fixed mass in every collapsed halo in which metal-free gas is able to cool. We assume that any halo that has already had stars previously formed in one of their halo progenitors will form only stars with metals, which are assigned an emissivity of ionizing radiation equal to that determined at z=4 from the measured intensity of the ionizing background. We examine the impact of molecular hydrogen photodissociation (which tends to reduce cooling when a photodissociating background is produced by the first stars) and X-Ray photoheating (which heats the atomic medium, raising the entropy of the gas before it collapses into halos). We find that in the CDM$\Lambda$ model supported by present observations, and even assuming no negative feedbacks for the formation of metal-free stars, a reionized mass fraction of 50% is not reached until redshift $z\sim 12$. The combination of ordinary metal-enriched stars and early metal-free stars can yield a CMB optical depth to electron scattering of at most 0.13. The contribution of metal-free stars to the present Cosmic Infrarred Background is negligibly small.

We present predictions of Ly-alpha forest-galaxy correlations at z=3 from Eulerian simulations that include the effects of galactic winds, driven primarily by supernova explosions. Galactic winds produce expanding bubbles of shock-heated gas within ~0.5 comoving Mpc/h of luminous galaxies in the simulation, which have space density similar to that of observed Lyman break galaxies (LBGs). However, most of the low-density intergalactic gas that determines the observed properties of the Ly-alpha forest is unaffected by winds. The impact of winds on the Ly-alpha optical depth near galaxies is less dramatic than their impact on gas temperature because winds heat only a small fraction of the gas present in the turnaround regions surrounding galaxies. Hence, Ly-alpha absorption from gas outside the wind bubbles is spread out over the same velocity range occupied by the wind-heated gas. In general, Ly-alpha absorption is expected to be stronger than average near galaxies because of the high gas density. Winds result in a modest reduction of this expected increase of Ly-alpha absorption. Our predictions can be compared to future observations to detect the wind effects and infer their strength, although with the caveat that the results are still dependent on the correspondence of simulated galaxies and observed LBGs. We find that wind effects in our simulations are not strong enough to reproduce the high Ly-alpha transmission within 0.5 Mpc/h comoving of galaxies that has been suggested by recent observations; powerful galactic explosions or ejecta with hyper-escape velocities would be required, but these are unlikely to be produced by ordinary star formation and supernovae alone

We use cosmological simulations of reionization to predict the effect of large HII regions around individual high-redshift quasars on the possible signal from the redshifted 21cm line of neutral hydrogen in the pre-reionization era. We show that these HII regions appear as ``spectral dips'' in frequency space with equivalent widths in excess of 3 mK MHz or depths in excess of about 1.5 mK, and that they are by far the most prominent cosmological signals in the redshifted HI distribution. These spectral dips are expected to be present in almost every line of sight. If the spectral dips of a large enough sample of HII regions are well resolved in frequency space, the distribution of line depth and equivalent width in frequency with a known observing beamsize can be used to infer the HII region size distribution and the mean difference in neutral hydrogen density between the \HII regions (which may contain self-shielded neutral gas clumps) and the surrounding medium, providing a powerful test for models of reionization.

An upper limit to the importance of local sources of radiation compared to the cosmic background in cosmological absorption systems is derived, as a simple consequence of the conservation of surface brightness. The limit depends only on the rate of incidence of the absorbers and the mean free path of the radiation. It is found that, on average, the ionizing radiation intensity from local sources in Lyman limit systems at z>2 must be less than half of the intensity of the cosmic background. In absorbers with column densities much lower than Lyman limit systems, the local source contribution must be negligible. The limit on the ratio of local source to background intensities is then applied to the class of damped Lya absorption systems with detectable excited CII lines. A cooling rate of the gas in these systems has been measured by Wolfe et al., who assumed that the balancing heating source is photoelectric heating on dust by light at ~ 1500 A . The intensity from local star formation at this wavelength in this class of damped Lya systems is found to be at most ~ 3 times the background intensity. If the heating source is indeed photoelectric heating of dust, the background created by sources associated with damped Lya systems can then be estimated from the average cooling rates measured in the absorbers. Current results yield a background intensity higher than previous estimates based on observed galaxy and quasar luminosity functions, although with a large uncertainty. The possibility of other sources of heating, such as shock-heating in a turbulent medium, should be explored.

A new theory of quasars is presented in which the matter of thin accretion disks around black holes is supplied by stars that plunge through the disk. Stars in the central part of the host galaxy are randomly perturbed to highly radial orbits, and as they repeatedly cross the disk they lose orbital energy by drag, eventually merging into the disk. Requiring the rate of stellar mass capture to equal the mass accretion rate into the black hole, a relation between the black hole mass and the stellar velocity dispersion is predicted of the form M_{BH} \propto sigma_*^{30/7}. The normalization depends on various uncertain parameters such as the disk viscosity, but is consistent with observation for reasonable assumptions. We show that a seed central black hole in a newly formed stellar system can grow at the Eddington rate up to this predicted mass via stellar captures by the accretion disk. Once this mass is reached, star captures are insufficient to maintain an Eddington accretion rate, and the quasar may naturally turn off as the accretion switches to a low-efficiency advection mode. The model provides a mechanism to deliver mass to the accretion disk at small radius, probably solving the problem of gravitational instability to star formation in the disk at large radius. We note that the matter from stars that is incorporated to the disk has an average specific angular momentum that is very small or opposite to that of the disk, and discuss how a rotating disk may be maintained as it captures this matter if a small fraction of the accreted mass comes from stellar winds that form a disk extending to larger radius. We propose several observational tests and consequences of this theory.

The recent discovery of luminous quasars up to a redshift z=6.43 has renewed interest in the formation of black holes massive enough to power quasars. If black holes grow by Eddington-limited gas accretion with a radiative efficiency of at least 10%, the time required to grow from a stellar black hole to ~10^9 msun is ~10^9 years, close to the age of the universe at z=6. Black hole mergers may accelerate the rate of mass growth, but can also completely eject black holes from halo centers owing to the gravitational wave recoil effect. Recently, Haiman concluded that black hole ejections likely do not allow black holes to grow to ~10^9 msun by z=6.43. We reexamine this problem and show that, by using a different halo escape velocity, accounting for the dependence of the recoil velocity on the black hole binary mass ratio and spins, and allowing seed black holes to form in all halos down to virial temperatures of 2000 K, black hole masses may reach ~10^9 msun as early as z=9 starting from stellar seeds, without super-Eddington accretion. In this particular case, we find that these massive black holes form from the merger of ~10^4 stellar black holes formed in low-mass halos at z~20, which must all grow close to the maximum Eddington rate over most of the time available from their birth to z~6. The alternative is that black holes can grow more rapidly by super-Eddington accretion.

The transverse proximity effect is the expected decrease in the strength of the Lya forest absorption in a QSO spectrum when another QSO lying close to the line of sight enhances the photoionization rate above that due to the average cosmic ionizing background. We select three QSOs from the Early Data Release of the Sloan Digital Sky Survey that have nearby foreground QSOs, with proper line of sight tangential separations of 0.50, 0.82, and 1.10 h^{-1} Mpc. We estimate that the ionizing flux from the foreground QSO should increase the photoionization rate by a factor (94, 13, 13) in these three cases, which would be clearly detectable in the first QSO and marginally so in the other two. We do not detect the transverse proximity effect. Three possible explanations are provided: an increase of the gas density in the vicinity of QSOs, time variability, and anisotropy of the QSO emission. We find that the increase of gas density near QSOs can be important if they are located in the most massive halos present at high redshift, but is not enough to fully explain the absence of the transverse proximity effect. Anisotropy requires an unrealistically small opening angle of the QSO emission. Variability demands that the luminosity of the QSO with the largest predicted effect was much lower 10^6 years ago, whereas the transverse proximity effect observed in the HeII Lya absorption in QSO 0302-003 by Jakobsen et al. (2003) implies a lifetime longer than 10^7 years. A combination of all three effects may better explain the lack of Lya absorption reduction. A larger sample of QSO pairs may be used to diagnose the environment, anisotropy and lifetime distribution of QSOs.

We investigate the interaction of Lyman alpha photons produced by the first stars in the universe with the intergalactic medium (IGM) prior to reionization. The background Lyman alpha spectral profile is obtained by solving a Fokker-Planck equation. Accurate values of the heating and scattering rates, and the spin-kinetic temperature coupling coefficient, are presented. We show that the heating rate induced by the Lyman alpha scatterings is much lower than found previously, and is basically negligible. The dominant heating source is most likely the X-rays from the first ionizing sources, which are able to penetrate into the atomic medium. The scattering of Lyman alpha photons couples the hydrogen spin temperature to the kinetic temperature. If the first ionizing sources in the universe did not emit significant X-rays, the spin temperature would be rapidly brought down to the very low gas kinetic temperature, and a 21cm absorption signal against the CMB larger than 100 mK would be predicted. However, we argue that sufficient X-rays are likely to have been emitted by the first stellar population, implying that the gas kinetic temperature should rapidly increase, turning a reduced and brief absorption signal into emission, with a smaller amplitude of about 10 mK. The detection of the 21cm absorption and emission feature would be a hallmark in unravelling the history of the ``dark age'' before reionization.

The Dark Age is the period between the time when the cosmic microwave background was emitted and the time when the evolution of structure in the universe led to the gravitational collapse of objects in which the first stars were formed. The period of reionization started with the ionizing light from the first stars, and it ended when all the atoms in the intergalactic medium had been reionized. The most distant sources of light known at present are galaxies and quasars at redshift z ~ 6, and their spectra indicate that the end of reionization was occurring just at that time. The Cold Dark Matter theory for structure formation predicts that the first sources formed much earlier.

The average rate of emission of ionizing radiation per unit volume (or emissivity) in the universe can be calculated as the ratio of the intensity of the ionizing background to the mean free path of ionizing photons. The intensity of the background is measured from the mean transmitted flux of the Lya forest, and the mean free path is measured from the abundance of Lyman limit systems, which has been observed so far up to z=4. This yields an emissivity that is not larger than 7 ionizing photons per Hubble time for each atom in the universe at z=4, which may reasonably arise from QSOs and star-forming galaxies. In order for the reionization to end by z=6, and assuming that the clumping factor of ionized gas during the reionization epoch is close to unity, this ionizing emissivity cannot decline from z=4 up to z ~ 9 by more than a factor 1.5. If the clumping factor were much larger than unity, then the emissivity would need to rapidly increase with redshift. Unless the ionizing emissivity increases substantially from z=4 to z ~ 10 - 20, the Thomson optical depth to the CMB must be in the range 0.045 < tau < 0.09.

We calculate the reionization history in Cold Dark Matter (CDM) models. The epoch of the end of reionization and the Thomson scattering optical depth to the cosmic microwave background depend on the power spectrum amplitude on small scales and on the ionizing photon emissivity per unit mass in collapsed halos. We calibrate the emissivity to reproduce the measured ionizing background intensity at z=4. Models in which all CDM halos have either a constant emissivity or a constant energy emitted per Hubble time, per unit mass, predict that reionization ends near z~6 and the optical depth is in the range 0.05 < tau_e < 0.09, consistent with WMAP results at the 1 to 2 sigma level. If the optical depth is as high as 0.17 (as suggested by WMAP), halos of velocity dispersion ~ 3-30 km/s at z>15 must have ionizing emissivities per unit mass larger by a factor >~ 50 compared to the more massive halos that produce the ionizing emissivity at z=4. This factor increases to 100 if the CDM power spectrum amplitude is required to agree with the Croft et al. (2002) measurement from the Lyman alpha forest. If tau_e >~ 0.17 were confirmed, a higher ionizing emissivity at z>15 compared to z=4 might arise from an enhanced star formation rate or quasar abundance per unit mass and an increased escape fraction for ionizing photons; the end of reionization could have been delayed to z~6 because of the suppression of gas accretion and star formation in low-mass halos as the medium was reionized.

A Monte Carlo code to solve the transfer of Lyman alpha (Lya) photons is developed, which can predict the Lya image and two-dimensional Lya spectra of a hydrogen cloud with any given geometry, Lya emissivity, neutral hydrogen density distribution, and bulk velocity field. We apply the code to several simple cases of a uniform cloud to show how the Lya image and emitted line spectrum are affected by the column density, internal velocity gradients, and emissivity distribution. We then apply the code to two models for damped Lya absorption systems: a spherical, static, isothermal cloud, and a flattened, axially symmetric, rotating cloud. If the emission is due to fluorescence of the external background radiation, the Lya image should have a core corresponding to the region where hydrogen is self-shielded. The emission line profile has the characteristic double peak with a deep central trough. We show how rotation of the cloud causes the two peaks to shift in wavelength as the slit is perpendicular to the rotation axis, and how the relative amplitude of the two peaks is changed. In reality, damped Lya systems are likely to have a clumpy gas distribution with turbulent velocity fields, which should smooth the line emission profile, but should still leave the rotation signature of the wavelength shift across the system.

We calculate the column density distribution of damped Lyman alpha systems, modeled as spherical isothermal gaseous halos ionized by the external cosmic background. The effects of self-shielding introduce a hump in this distribution, at a column density N_{HI} \sim 1.6x10^{17} X^{-1} cm^{-2}, where X is the neutral fraction at the radius where self-shielding starts being important. The most recent compilation of the column density distribution by Storrie-Lombardi & Wolfe shows marginal evidence for the detection of this feature due to self-shielding, suggesting a value X \sim 10^{-3}. Assuming a photoionization rate \Gamma \sim 10^{-12} s^{-1} from the external ionizing background, the radius where self-shielding occurs is inferred to be about 3.8kpc. If damped Lyman alpha systems consist of a clumpy medium, this should be interpreted as the typical size of the gas clumps in the region where they become self-shielding. Clumps of this size with typical column densities N_H \sim 3x10^{20} cm^{-2} would be in hydrostatic equilibrium at the characteristic photoionization temperature \sim 10^4 K if they do not contain dark matter. Since this size is similar to the overall radius of damped \lya systems in Cold Dark Matter models, where all halos are assumed to contain similar gas clouds producing damped absorbers, this suggests that the gas in damped absorbers is in fact not highly clumped.

We show that at any given time, the Galactocentric black hole Sgr A* is expected to be microlensing N(lens)=1.7 bulge stars if the threshold of detectability of the fainter image is K(thr)=21, and about N(lens)=8 sources if K(thr)=23. The lensed images then provide a unique way to detect stellar-mass black holes predicted to cluster around Sgr A*. If a black hole passes close to a microlensed image, it will give rise to a short (weeks long) microlensing event. We show that the mass and projected velocity of the stellar-mass black hole can both be measured unambiguously from such an event, provided that either a caustic crossing is observed or the astrometric displacement is measured. For K(thr)=23 and moderate magnifications by Sgr A*, the microlensing event rate from a cluster of 20000 black holes within a radius of 0.7 pc is only 0.06 events per year; however, if highly magnified images of a star were found, the rate of events by the stellar black holes would be much higher. In addition, the N(lens) sources lensed by Sgr A* provide a unique probe of extinction behind the Galactic center along 2N(lens) lines of sight.

The recent discovery of a planetary transit in the star HD 209458, and the subsequent highly precise observation of the transit lightcurve with Hubble Space Telescope, is encouraging to search for any phenomena that might induce small changes in the light curve. Here we consider the effect of the quadrupole moment of the parent star and of a possible second planet perturbing the orbit of the transiting planet. Both of these cause a precession of the orbital plane and of the periastron of the planet, which result in a long-term variation of the duration and the period of the transits. For a transiting planet at 0.05 AU, either a quadrupole moment similar to that of the Sun or the gravitational tug from an Earth-like planet on an orbit of semimajor axis ~ 0.2 AU and a relative inclination near the optimal 45 degrees would cause a transit duration time derivative of ~ 1 second per year.

Spergel & Steinhardt proposed the possibility that the dark matter particles are self-interacting, as a solution to two discrepancies between the predictions of cold dark matter models and the observations: first, the observed dark matter distribution in some dwarf galaxies has large, constant-density cores, as opposed to the predicted central cusps; and second, small satellites of normal galaxies are much less abundant than predicted. The dark matter self-interaction would produce isothermal cores in halos, and would also expel the dark matter particles from dwarfs orbiting within large halos. However, another inevitable consequence of the model is that halos should become spherical once most particles have interacted. Here, I rule out this model by the fact that the innermost regions of dark matter halos in massive clusters of galaxies are elliptical, as shown by gravitational lensing and other observations. The absence of collisions in the lensing cores of massive clusters implies that any dark matter self-interaction is too weak to have affected the observed density profiles in the dark-matter dominated dwarf galaxies, or to have eased the destruction of dwarf satellites in galactic halos. If $s_x$ is the cross section and $m_x$ the mass of the dark matter particle, then $s_x/m_x < 10^{-25.5} \cm^2/\gev$.

We calculate the neutral hydrogen column density of self-absorption in QSOs predicted in a model where the QSOs are located in the same halos that contain the gas in damped Lya absorption systems. The model is parameterized by the probability P_0 that any halo has an active QSO. We assume that the QSOs ionize the gas, but do not expel or heat it. The derived HI column densities produce negligible Lyman limit absorption, even in the lowest luminosity QSOs, with an optical depth of only ~ 10% for luminosity L=0.01 L_*, when P_0=10^{-2}. We also compute the HeII Lyman limit self-absorption, which is slightly higher but still negligible. The self-absorption can be higher if the gas is highly clumped; only in this case the overall emissivity from QSOs could be significantly reduced due to absorption by the known damped Lya systems, to affect the predicted intensity of the ionizing background or the epoch of HeII reionization. The presence of the gas associated with damped absorption systems around QSOs could also be detected from the narrow Lya emission line, which should have and angular extent of 0.1 to 1'' in typical high-redshift QSOs.

We use the redshift evolution of the lyman-alpha forest mean transmitted flux at z>2.5 to infer the evolution of the intensity of the ionizing background, using theoretical predictions for the density distribution of the intergalactic gas giving rise to the lyman-alpha forest. The proper background intensity declines gradually with redshift, decreasing by a factor \sim 3 between z=3 and z=5. The gravitational structure formation theory of the lyman-alpha forest implies that, given the observed mean transmitted flux of \bar{F}=0.1 in the highest redshift quasar (with a mean lyman-alpha forest redshift z=5.2), only about 3% of the lyman-alpha spectrum should have a flux higher than 50% of the continuum, and less than 0.5% of pixels should have flux above 75% of the continuum, assuming a uniform ionizing background intensity. We show that this is consistent with the spectrum shown in Fan et al.

A large number of microlensing events have been observed in the direction of the Galactic bulge, with a measured optical depth in the range 2 - 3 x 10^{-6}. It has been shown that most of these events are due to bulge stars being lensed by other bulge stars or by foreground disk stars. Among the stars observed in the bulge fields, there should also be disk stars located behind the bulge; here, we consider their effect on the microlensing rates. The optical depth of background disk stars is much higher than that of typical bulge stars, reaching 10^{-5} at 6 Kpc behind the bulge. Thus, although background disk stars are a very small fraction of the stars in Baade's window, we find that ~ 5 to 10 % of the optical depth should be due to disk stars more than 3 kpc behind the bulge. This fraction is sensitive to the luminosity function of disk stars at large scale-height, to the magnitude cutoff of the survey, and to the amplification bias effect causing large numbers of ``blended'' events. We consider also the effect of a warp and flare in the disk at large distances behind the bulge; this could increase the optical depth from the background disk to ~ 20 % of the total. Events on background disk stars should on average be longer than other events and could be distinguished also by measuring the proper motion or distance of the stars that have been microlensed. The number of these events could be an interesting probe to the structure and stellar population of the far-side of the Galactic disk.

The evolution of the temperature in the intergalactic medium is related to the reionization of hydrogen and helium, and has important consequences for our understanding of the Lya forest and of galaxy formation in gravitational models of large-scale structure. We measure the temperature-density relation of intergalactic gas from Lya forest observations of eight quasar spectra with high resolution and signal-to-noise ratio, using a new line fitting technique to obtain a lower cutoff of the distribution of line widths from which the temperature is derived. We carefully test the accuracy of this technique to recover the gas temperature with a hydrodynamic simulation. The temperature at redshift z=(3.9, 3.0, 2.4) is best determined at densities slightly above the mean: T_star=(20200\pm2700, 20200\pm1300, 22600\pm1900)K (statistical error bars) for gas density (in units of the mean density) Delta_star=(1.42\pm0.08, 1.37\pm0.11, 1.66\pm0.11). The power-law index of the temperature-density relation, defined by T=T_star(Delta/Delta_star)^{gamma-1}, is gamma-1= (0.43\pm0.45, 0.29\pm0.30, 0.52\pm0.14) for the same three redshifts. The temperature at the fixed over-density Delta=1.4 is T_1.4=(20100\pm2800, 20300\pm1400, 20700\pm1900)K. These temperatures are higher than expected for photoionized gas in ionization equilibrium with a cosmic background, and can be explained by a gradual additional heating due to on-going HeII reionization. The measurement of the temperature reduces one source of uncertainty in the lower limit to the baryon density implied by the observed mean flux decrement. We find that the temperature cannot be reliably measured for under-dense gas, because the velocities due to expansion always dominate the widths of the corresponding weak lines.

If the stellar population of the bulge contains black holes formed in the final core collapse of ordinary stars with M \ga 30 M_{\odot}, then about 25,000 stellar mass black holes should have migrated by dynamical friction into the central parsec of the Milky Way, forming a black hole cluster around the central supermassive black hole. These black holes can be captured by the central black hole when they randomly reach a highly eccentric orbit due to relaxation, either by direct capture (when their Newtonian peribothron is less than 4 Schwarzschild radii), or after losing orbital energy through gravitational waves. The overall depletion timescale is ~ 30 Gyr, so most of the 25,000 black holes remain in the central cluster today. The presence of this black hole cluster would have several observable consequences. First, the low-mass, old stellar population should have been expelled from the region occupied by the black hole cluster due to relaxation, implying a core in the profile of solar-mass red giants with a radius of ~ 2 pc (i.e., 1'). The observed central density cusp (which has a core radius of only a few arc seconds) should be composed primarily of young (\la 1 Gyr) stars. Second, flares from stars being captured by supermassive black holes in other galaxies should be rarer than usually expected because the older stars will have been expelled from the central regions by the black hole clusters of those galaxies. Third, the young (\la 2 Gyr) stars found at distances ~ 3 - 10 pc from the Galactic center should be preferentially on highly eccentric orbits. Fourth, if future high-resolution $K$-band images reveal sources microlensed by the Milky Way's central black hole, then the cluster black holes could give rise to secondary (``planet-like'') perturbations on the main event.

Weak gravitational lensing observations probe the spectrum and evolution of density fluctuations and the cosmological parameters which govern them but are currently limited to small fields and subject to selection biases. We show how the expected signal from large-scale structure arises from the contributions from and correlations between individual halos. We determine the convergence power spectrum as a function of the maximum halo mass and so provide the means to interpret results from surveys that lack high mass halos either through selection criteria or small fields. Since shot noise from rare massive halos is mainly responsible for the sample variance below 10', our method should aid our ability to extract cosmological information from small fields.

Several emerging links between high-redshift observational cosmology and the Galactic fossil evidence found in the kinematics, metallicities and ages of Milky Way stars are discussed. In a flat Cold Dark Matter model with $\Omega\simeq 0.3$ that agrees with present large-scale structure observations, the oldest stars in the Milky Way should have formed in the first halos where gas was able to cool, at $z\simeq 20$. These earliest, weakly bound dwarf galaxies probably turned only a small fraction of their gas to stars, which should be metal-poor. However, the merging rate in the early universe was much faster than the present one, so massive halos with more efficient star formation and metallicities up to the highest values present today in the bulge could have formed less than $10^9$ years after the oldest stars. The mean metallicity produced in the universe by a given redshift is related to the mean surface brightness of star-forming galaxies above this redshift, and also to the reionization epoch if galaxies were the dominant sources of ionizing radiation. The biased distribution of the early dwarf galaxies where the first stars formed should result in an age gradient with radius of the low-metallicity stars in the Milky Way, with the oldest ones concentrated in the bulge and the youngest in the outer halo.

A sample of eight quasars observed at high resolution and signal-to-noise is used to determine the probability distribution function (PDF), the power spectrum, and the correlation function of the transmitted flux in the \lya forest, in three redshift bins centered at z=2.41, 3.00, and 3.89. All the results are presented in tabular form, with full error covariance matrices to allow for comparisons with any numerical simulations and with other data sets. The observations are compared with a numerical simulation of the \lya forest of a Lambda-CDM model with Omega=0.4, known to agree with other large-scale structure observational constraints. There is excellent agreement for the PDF, if the mean transmitted flux is adjusted to match the observations. A small difference between the observed and predicted PDF is found at high fluxes and low redshift, which may be due to the uncertain effects of fitting the spectral continuum. Using the numerical simulation, we show how the flux power spectrum can be used to recover the initial power spectrum of density fluctuations. From our sample of eight quasars, we measure the amplitude of the mass power spectrum to correspond to a linear variance per unit ln(k) of $\Delta^2_\rho(k)=0.72\pm0.09$ at k=0.04(km/s)^{-1} and z=3, and the slope of the power spectrum near the same k to be $n_p=-2.55\pm0.10$ (statistical error bars). The results are statistically consistent with Croft et. al. (1999), although our value for the rms fluctuation is lower by a factor 0.75. For the Lambda-CDM model we use, the implied primordial slope is $n=0.93\pm0.10$, and the normalization is $\sigma_8=0.68+1.16(0.95-n)\pm0.04$.

The Lyman alpha absorption from intergalactic, once-ionized helium (HeII) has been measured with HST in four quasars over the last few years, over the redshift range 2.4 < z < 3.2. These observations have indicated that the HeII reionization may not have been completed until z\simeq 2.8, and that large fluctuations in the intensity of the HeII-ionizing background were present before this epoch. The detailed history of HeII reionization at higher redshifts is, however, model-dependent and difficult to determine from these observations, because the IGM can be completely optically thick to Lya photons when only a small fraction of the helium remains as HeII. In addition, finding quasars in which the HeII Lya absorption can be observed becomes increasingly difficult at higher redshift, owing to the large abundance of hydrogen Lyman limit systems. It is pointed out here that HeII in the IGM should also cause detectable continuum absorption in the soft X-rays. The spectrum of a high-redshift source seen behind the IGM when most of the helium was HeII should recover from the HeII Lyman continuum absorption at an observed energy \sim 0.1 keV. Galactic absorption will generally be stronger, but not by a large factor; the intergalactic HeII absorption can be detected as an excess over the expected Galactic absorption from the 21cm HI column density. In principle, this method allows a direct determination of the fraction of helium that was singly ionized as a function of redshift, if the measurement is done on a large sample of high-redshift sources over a range of redshift.

This review presents a brief discussion of the theory of Cold Dark Matter of structure formation in the universe, describing the main processes determining the power spectrum, the non-linear gravitational collapse, the formation of galaxies and the evolution of the intergalactic medium. Recent advances in the observations of high redshift galaxies, and their interpretation in the context of this theory, are then summarized. Este art\'\i culo de revisi\'on presenta una breve exposici\'on de la teor\'\i a de la Materia Invisible Fr\'\i a de la formaci\'on de estructura en el universo, describiendo los procesos m\'as importantes que determinan el espectro de potencia, el colapso gravitatorio no lineal, la formaci\'on de galaxias y la evoluci\'on del medio intergal\'actico. Resumimos tambi\'en los avances recientes en las observaciones de galaxias a alto corrimiento al rojo, y la interpretaci\'on en el contexto de esa teor\'\i a.

A model of the density distribution in the intergalactic medium, motivated by that found in numerical simulations, is used to demonstrate the effect of a clumpy IGM and discrete sources on the reionization of the universe. In an inhomogeneous universe reionization occurs outside-in, starting in voids and gradually penetrating into overdense regions. Reionization should not be sudden but gradual, with a continuous rise of the photon mean free path over a fair fraction of the Hubble time as the emissivity increases. We show that a hydrogen Gunn-Peterson trough should be present at $z\simeq 6$ unless the emissivity increases with redshift at $z>4$. However, the epoch of overlap of cosmological \hii regions could have occurred at a higher redshift if sources of low luminosity reionized the IGM; the Gunn-Peterson trough at $z\sim 6$ would then appear because even the most underdense voids have a large enough neutral fraction in ionization equilibrium to be optically thick to \lya photons. Cosmological \hii regions near the epoch of overlap can produce gaps of transmitted flux only if luminous quasars contributed to the reionization. Despite the clumpiness of the matter distribution, recombinations are not very important during the reionization of hydrogen because the high density gas is not ionized until a late time. We show that the \heii reionization was most likely delayed relative to the hydrogen reionization, but should be completed by $z\sim 3$, the redshift where observations are available. The reported large optical depth fluctuations of \heii are probably not due to an incomplete \heii reionization, but arise from a combination of density fluctuations and the variations in the intensity of the ionizing background due to luminous QSO's.

A model of damped lyman alpha systems is presented based on randomly moving clouds in spherical halos. We use the Press-Schechter model for the abundance of halos, and assume that each halo has a similar population of clouds, with total mass and spatial distribution constrained to fit observations of the column density distribution. We show that the kinematics of the multiple absorbing components revealed in absorption profiles of the low-ionization lines, presented by Prochaska & Wolfe, are consistent with our spherical halo model. The presence of multiple absorbing components with a large covering factor, combined with the small impact parameters of the systems predicted in our analytical model and in numerical simulations, implies a high rate of energy dissipation in cloud collisions. We calculate the rate of energy dissipation in our model, and show that it is far greater than the rate at which energy can be supplied by gravitational mergers of halos. This poses a problem for the model of merging protogalactic clumps of Haehnelt et al., based on numerical simulations. We also present new constraints on the amplitude of the power spectrum in hierarchical theories required to account for the observed velocity dispersion in the absorbers. We find that the linearly extrapolated rms fluctuation on spheres of radius HR = 100 km/s at z=4 must be greater than 0.75. Although this limit is obtained only for our specific model of the absorbing components, it should not be highly model-dependent because the velocity dispersion of the absorbers is essentially determined by the velocity dispersion of the halos where the gas is moving.

Gamma-ray bursts are now known to be a cosmological population of objects, which are often accompanied by X-ray and optical afterglows. The total energy emitted in the afterglow can be similar to the energy radiated in the gamma-ray burst itself. If a galaxy containing a large column density of dust is near the line of sight to a gamma-ray burst, small-angle scattering of the X-rays due to diffraction by the dust grains will give rise to an X-ray echo of the afterglow. A measurement of the angular size of the echo at a certain time after the afterglow is observed yields a combination of the angular diameter distances to the scattering galaxy and the gamma-ray burst that can be used to constrain cosmological models in the same way as a time delay in a gravitational lens. The scattering galaxy will generally cause gravitational lensing as well, and this should modify the shape of the X-ray echo from a circular ring. The main difficulty in detecting this phenomenon is the very low flux expected for the echo. The flux can be increased when the gamma-ray burst is highly magnified by gravitational lensing, or when the deflecting galaxy is at low redshift. X-ray echos of continuous (but variable) sources, such as quasars, may also be detectable with high-resolution instruments and would allow similar measurements.

We discuss the feasibility of measuring the cosmological metric using the redshift space correlation function of the Lya forest in multiple lines of sight, as a function of angular and velocity separation. The geometric parameter that is measured is f(z) = H(z) D(z)/c, where H(z) is the Hubble constant and D(z) the angular diameter distance at redshift z. The correlation function is computed in linear theory. We describe a method to measure it from observations with the Gaussianization procedure of Croft et al (1998) to map the Lya forest transmitted flux to an approximation of the linear density field. The effect of peculiar velocities on the shape of the recovered power spectrum is pointed out. We estimate the error in recovering the f(z) factor from observations due to the variance in the Lya absorbers. We show that ~ 20 pairs of quasars (separations < 3') are needed to distinguish a flat \Omega_0=1 universe from a universe with \Omega_0=0.2, \Omega_\Lambda=0.8. A second parameter that is obtained from the correlation function of the Lya forest is \beta \simeq \Omega(z)^{0.6}/b (affecting the magnitude of the peculiar velocities), where b is a linear theory bias of the Lya forest. The statistical error of f(z) is reduced if b can be determined independently from numerical simulations, reducing the number of quasar pairs needed for constraining cosmology to approximately six. On small scales, where the correlation function is higher, f(z) should be measurable with fewer quasars, but non-linear effects must then be taken into account. The anisotropy of the non-linear redshift space correlation function as a function of scale should also provide a precise quantitative test of the gravitational instability theory of the Lya forest.

We use numerical simulations of structure formation in a Cold Dark Matter model to predict the absorption lines in the soft X-rays produced by heavy elements in the shock-heated intergalactic medium at low redshift. The simulation incorporates a model for metal production in galaxies and the subsequent dispersion of the metals to the intergalactic medium. We analyze in particular absorption lines produced by oxygen, and calculate the ionization stage taking into account the observed X-ray background at the present time. We find that oxygen is fully ionized by the X-ray background in low-density voids, and is mostly in the form of OVII and OVIII in the sheets and filamentary regions. Strong absorption lines of OVII and OVIII with equivalent widths W ~ 100 km/s are produced in filamentary regions of overdensities ~ 100 and temperatures ~ 10^6 K, located in the outskirts of groups and clusters of galaxies. The OVII line at E = 574 eV is generally the strongest one in these systems. Our model predicts that any X-ray source (such as a quasar) should typically show about one OVII absorption line with W > 100 km/s in the interval from z=0 to z=0.3. These lines could be detected with the upcoming generation of X-ray telescopes, and their origin in intervening systems could be confirmed by the association with groups of galaxies and X-ray emitting halos near the line-of-sight at the same redshift. The hot intergalactic medium may be one of the main reservoirs of baryons in the present universe, and the heavy element X-ray absorption lines offer a promising possibility of detecting this new component in the near future.

If the universe was reionized by O and B stars in an early population of galaxies, the associated supernovae should have enriched the universe to a mean metallicity $\bar Z = 10^{-5} (1+ n_{rec})$, where $n_{rec}$ is the mean number of times that each baryon recombined during the reionization era. This is consistent with recent observations of the metallicity in the Lya forest at $z\simeq 3$. The mean surface brightness observable at present from the galaxies that produced these heavy elements, in the rest-frame wavelengths $1216\AA < \lambda \lesssim 2500\AA$, should be $\sim 10^{-6} (\bar Z/10^{-4}) (\Omega_b h^2/0.02) ~ photons \cm^{-2} s^{-1} arcsec^{-2}$. Most of this radiation was emitted at $z>5$, before reionization was complete. These high-redshift galaxies may be detectable in near-infrared photometric surveys, identifying them via the Gunn-Peterson trough (analogous to the use of the Lyman limit cutoff to search for galaxies at $z\sim 3$, where the Lya forest blanketing is smaller). Their spectrum may also be characterized by a strong Lya emission line. However, the spectra of galaxies seen behind intervening gas that is still neutral should show the red damping wing of the Gunn-Peterson trough, with a predictable profile that obstructs part of the Lya emission. The low-mass galaxies formed before reionization might constitute a distinctive population; we discuss the signature that this population could have in the faint number counts. Although most of these galaxies should have merged into larger ones, those that survived to the present could be dwarf spheroidals.

We propose a new method of analysis for the \lya forest, namely to measure the 1-point and 2-point joint probability distribution of the transmitted flux. The results for a sample of seven observed quasars and from two simulations of structure formation are shown and compared. Statistically significant differences in the 2-point function between the results of the numerical simulations and the observations are easily found. The analysis we suggest is very simple to apply to observed data sets, and we discuss its superiority over the traditional Voigt-profile fitting algorithms for accurate comparison to the predictions of theoretical models.

Observations of high-redshift quasars show that the IGM must have been reionized at some redshift $z>5$. If a source of radiation could be observed at the rest-frame Lya wavelength, at a sufficiently high redshift where some of the IGM in the line-of-sight was not yet reionized, the Gunn-Peterson trough should be present. Longward of the Lya wavelength, a damping wing should be observed caused by the neutral IGM whose absorption profile can be predicted. Measuring the shape of this damping wing would provide irrefutable evidence of the observation of the IGM before reionization, and a determination of the density of the neutral IGM. This measurement might be hindered by the possible presence of a dense absorption system associated with the source. Shortward of the \lya wavelength, absorption should be seen from the patchy structure of the IGM in the process of reionization. We show that a complete Gunn-Peterson trough is most likely to continue to be observed through the epoch where the IGM is partially ionized. The damping wings of the neutral patches should overlap if the proper pathlength through an ionized region is less than 1 h^{-1} Mpc; even in larger ionized regions, the characteristic background intensity should be low enough to yield a very high optical depth due to the residual neutral fraction, although occasionally some flux may be transmitted through large, underdense voids within an ionized region. The case of the HeII reionization is also discussed, and we argue that helium was already doubly ionized by z=3 throughout the IGM. The recently discovered afterglows of gamma-ray bursts might soon be observed at very high redshifts. Their featureless continuum spectrum and high luminosities make them ideal sources for studying absorption by the IGM.

Recent evidence on the metal content of the high-redshift Lyman alpha forest seen in quasar spectra suggests that an early generation of galaxies enriched the intergalactic medium (IGM) at z ~ 5. We calculate the number of supernovae that need to have taken place to produce the observed metallicity. The progenitor stars of the supernovae should have emitted $\sim 20$ ionizing photons for each baryon in the universe, i.e., more than enough to ionize the IGM. We calculate that the rate of these supernovae is such that about one of them should be observable at any time per square arc minute. Their fluxes are, of course, extremely faint: at z=5, the peak magnitude should be K=27 with a duration of ~ 1 year. However, these supernovae should still be the brightest objects in the universe beyond some redshift, because the earliest galaxies should form before quasars and they should have very low mass, so their luminosities should be much lower than that of a supernova. We also show that, under the assumption of a standard initial mass function, a significant fraction of the stars in the Galactic halo should have formed in the early galaxies that reionized and enriched the IGM, and which later must have merged with our Galaxy. These stars should have a more extended radial distribution than the observed halo stars.

We derive analytic lower bounds on the cosmic baryon density by requiring that the high-redshift IGM contain enough neutral hydrogen to produce the observed \lya absorption in quasar spectra. The key theoretical assumption is that absorbing structures are no more extended in redshift space than in real space. This assumption might not hold if \lya clouds are highly overdense and thermally broadened, but it is likely to hold in the gravitational instability picture for the \lya forest suggested by cosmological simulations, independently of the details of the cosmological model. The other ingredients in these bounds are an estimate of the UV background from quasars, a temperature $T~10^4K$ for the "warm" photoionized IGM that produces most of the \lya absorption, a value of the Hubble constant, and observational estimates of the mean \lya flux decrement $\Dbar$ or, for a more restrictive bound, the distribution function $P(\tau)$ of \lya optical depths. With plausible parameter values, the mean decrement bound implies a baryon density parameter $\Omb \ga 0.0125/h^2$. With conservative values, the bound weakens to $\Omb \ga 0.005/h^2$, but the required clustering of the IGM is then incompatible with other properties of quasar spectra. A recent observational determination of $P(\tau)$ implies $\Omb \ga 0.0125/h^2$ even for a conservative estimate of the UV background, and $\Omb \ga 0.018/h^2$ for a more reasonable estimate. These bounds are consistent with recent low estimates of the primordial deuterium-to-hydrogen ratio $\dtoh$, which imply $\Omb \approx 0.025/h^2$, but the $P(\tau)$ bound can only be reconciled with high $\dtoh$ estimates by abandoning standard big bang nucleosynthesis or the gravitational instability picture for the origin of the \lya forest. (Shortened abstract.)

We have measured the distribution function of the flux decrement caused by Lyman alpha forest absorption in a new sample of high resolution QSO spectra. The observations are compared to the results from two simulations of the Lya forest: an Eulerian Lambda-CDM model, and an SPH standard CDM model. Good agreement between the shapes of simulated and observed distributions is achieved by globally scaling the optical depth to match the mean flux decrements. This procedure amounts to a measurement of the parameter Omega_b^2 h^3 / Gamma (where Omega_b is the baryonic matter density and Gamma is the HI ionization rate). Estimating a lower limit Gamma > 7 10^{-13} s^{-1} from the abundance of known QSOs, we derive a lower limit to the baryon density, Omega_b h^2>0.021(0.017) for the Lambda-CDM (SCDM) model. In both cases the large values are inconsistent with some recent D/H determinations (Rugers & Hogan 1996a,b), favoring a low deuterium abundance as reported by Tytler, Fan & Burles (1996). Adopting a fixed Omega_b, we can determine the evolution of the ionizing radiation field. Our models predict the intensity to be approximately constant with redshift, consistent with the assumption that the ionizing background is produced by known quasars for z < 3. However, additional sources of ionizing photons are required at higher redshift.

We propose that relative variability on short time-scales of the multiple images of a lensed quasar, after removal of the time delay, may be caused by hot spots or other moving structures in the accretion disk crossing microlens caustics caused by stellar mass objects in the lensing galaxy. Such variability has been reported in the two images of 0957+561. The short durations would be due to the high rotation speed of the disk (v/c ~ 0.1), rather than planetary mass objects in the slowly moving (v/c ~ 0.001) lens. This interpretation could be confirmed by finding periodicity, or correlations of the spectral and flux variations due to the Doppler effect in the disk. We also propose another signature of stationary accretion disks (with no intrinsic variability): the gradient of the magnification over the accretion disk should cause a relative color change between the images whose sign and amplitude are correlated with the time derivative of the flux difference between the images. Other color terms induced by the radial variation of disk colors are of second order in the magnification gradient. The methods proposed here can be used first to verify that accretion disks near supermassive black holes are the source of the continuum radiation from quasars, and then to study them.

It has recently been shown that the highest energy cosmic rays (CRs) may originate in the same cosmological objects producing $\gamma$-ray bursts. This model requires the presence of intergalactic magnetic fields (IGMF) to delay the arrival times of $\sim 10^{20}$ eV CRs by 50 years or longer relative to the $\gamma$-rays, of an amplitude that is consistent with other observational constraints. Sources of CRs coming from individual bursts should be resolved with the planned ``Auger'' experiment, with as many as hundreds of CRs for the brightest sources. We analyze here the apparent angular and energy distribution of CRs from bright sources below the pion production threshold (in the energy range $10^{19}{\rm eV} < E < 4\times10^{19}{\rm eV}$) expected in this model. This observable distribution depends on the structure of the IGMF: the apparent spectral width $\Delta E$ is small, $\Delta E/E\lesssim1\%$, if the intergalactic field correlation length $\lambda$ is much larger than $1{\rm Mpc}$, and large, $\Delta E/E=0.3$, in the opposite limit $\lambda\ll 1{\rm Mpc}$. The apparent angular size is also larger for smaller $\lambda$. If the sources of CRs we predict are found, they will corroborate the bursting model and they will provide us with a technique to investigate the structure of the IGMF.

Microlensing events are now regularly being detected by monitoring the flux of a large number of potential sources and measuring the combined magnification of the images. This phenomenon could also be detected directly from the gravitational deflection, by means of high precision astrometry using interferometry. Relative astrometry at the level of $10\muas$ may become possible in the near future. The gravitational deflection can be measured by astrometric monitoring of a bright star having a background star within a small angular separation. This type of monitoring program will be carried out for the independent reasons of discovering planets from the angular motion they induce on the nearby star around which they are orbiting, and for measuring parallaxes, proper motions and orbits of binary stars. We discuss three applications of the measurement of gravitational deflections by astrometric monitoring: measuring the mass of the bright stars that are monitored, measuring the mass of brown dwarfs or giant planets around the bright stars, and detecting microlensing events by unrelated objects near the line of sight to the two stars. We discuss the number of stars whose mass could be measured by this procedure. We also give expressions for the number of expected microlensing events by unrelated objects, which could be stars, brown dwarfs, or other compact objects accounting for dark matter in the halo or in the disk.

A model of Lyman limit QSO absorption systems is investigated where they are produced in gaseous galactic halos with a two-phase structure: a hot phase at the halo virial temperature in approximate pressure equilibrium, and a cold, photoionized phase in the form of clouds confined by the pressure of the hot medium, and falling through it to the halo center. We show that the masses of these clouds must be over a relatively narrow range so that they are stable against gravitational collapse and evaporation. These masses, near $10^6 \msun$, also lead to a covering factor near unity. The hot phase is required to have a core radius such that its cooling time in the core is equal to the age of the halo, and the mass in the cold phase is determined by the rate at which the hot gas cools. We calculate the number of Lyman limit systems arising in halos of different masses in a CDM model and their impact parameters implied by these assumptions, and find them to be in reasonable agreement with observations. The observed low ionization systems such as MgII are well reproduced as arising from the photoionized phase, and are mostly in halos around massive galaxies, while CIV selected systems are predicted to be more commonly associated with lower mass galaxies and larger impact parameters. The hot phase may give rise to detectable absorption lines in OVI, while the column densities predicted for other highly ionized species are low and difficult to observe.

We derive observational consequences of the hypothesis that cosmic rays (CR's) of energy $>10^{19}eV$ originate in the same cosmological objects producing gamma-ray bursts (GRB's). Inter-galactic magnetic fields $\gtrsim 10^{-12} G$ are required in this model to allow CR's to be observed continuously in time by producing energy dependent delays in the CR arrival times. This results in individual CR sources having very narrow observed spectra, since at any given time only those CR's having a fixed time delay are observed. Thus, the brightest CR sources should be different at different energies. The average number of sources contributing to the total CR flux decreases with energy much more rapidly than in a model of steady CR sources, dropping to one at $E_{crit} \simeq2\times10^{20}$~eV with very weak sensitivity to the inter-galactic magnetic field strength. Below $E_{crit}$, a very large number of sources is expected, consistent with observations. Above $E_{crit}$, a source may be observed with a flux considerably higher than the time-averaged CR flux from all sources, if a nearby GRB occurred recently. If such a source is present, its narrow spectrum may produce a ``gap'' in the total spectrum. These signatures should be detectable by the planned ``Auger'' CR experiment.

Cosmological simulations with gas dynamics suggest that the Lyman-alpha forest is produced mainly by "small scale structure" --- filaments and sheets that are the high redshift analog of today's galaxy superclusters. There is no sharp distinction between Lyman-alpha clouds and "Gunn-Peterson" absorption produced by the fluctuating IGM -- the Lyman-alpha forest {\it is} the Gunn-Peterson effect. Lyman limit and damped Lyman-alpha absorption arises in the radiatively cooled gas of forming galaxies. At $z~2-3$, most of the gas is in the photoionized, diffuse medium associated with the Lyman-alpha forest, but most of the {\it neutral} gas is in damped Lyman-alpha systems. We discuss generic evolution of cosmic gas in a hierarchical scenario of structure formation, with particular attention to the prospects for detecting 21cm emission from high redshift HI. A scaling argument based on the present-day cluster mass function suggests that objects with $M_{HI} >~ 5e11 h^{-1} \msun$ should be extremely rare at $z~3$, so detections with existing instruments will be difficult. An instrument like the proposed Square Kilometer Array could detect individual damped Lyman-alpha systems at high redshift, making it possible to map structure in the high redshift universe in much the same way that today's galaxy redshift surveys map the local large scale structure.

The X-ray emission by hot gas at the centers of clusters of galaxies is commonly modeled assuming the existence of steady-state, inhomogeneous cooling flows. We derive the metallicity profiles of the intracluster medium expected from such models. The inflowing gas is chemically enriched by type Ia supernovae and stellar mass loss in the outer parts of the central galaxy, which may give rise to a substantial metallicity gradient. The amplitude of the expected metallicity enhancement towards the cluster center is proportional to the ratio of the central galaxy luminosity to the mass inflow rate. The metallicity of the hotter phases is expected to be higher than that of the colder, denser phases. The metallicity profile expected for the Centaurus cluster is in good agreement with the metallicity gradient recently inferred from ASCA measurements (Fukazawa et al. 1994). However, current data do not rule out alternative models where cooling is balanced by some heat source. The metallicity gradient does not need to be present in all clusters, depending on the recent merging history of the gas around the central cluster galaxy, and on the ratio of the stellar mass in the central galaxy to the gas mass in the cooling flow.

We model the $\lya$ forest in a spatially flat, CDM model with $\Omega=0.4$, with an Eulerian hydrodynamic simulation, and find that the intergalactic, photoionized gas collapses into sheet-like and filamentary structures which produce absorption lines like those observed in the $\lya$ forest. A typical filament is $\sim 1 h^{-1} \mpc$ long with thickness $\sim 50-100 h^{-1} \kpc$ (in proper units), and baryonic mass $\sim 10^{10} h^{-1} \msun$. The gas temperature increases with time as structures with larger velocities collapse gravitationally. The predicted distributions of column densities and b-parameters of the absorption lines agree with the observed ones, and their evolution is consistent with our model if the ionizing background has an approximately constant intensity between $z=2$ and $z=4$. We also predict the distribution of transmitted flux and its correlation along a spectrum and on parallel spectra, and the \heii flux decrement as a function of redshift. We predict a correlation length of $\sim 80 h^{-1}\kpc$ perpendicular to the line of sight for features in the Lyman alpha forest. In order to reproduce the observed number of lines and average flux transmission, the baryon content of the clouds may need to be significantly higher than in previous models because of the low densities we predict. If the background intensity $\jhi$ is at least that predicted from the observed quasars, $\Omega_b$ needs to be $\sim 0.025 h^{-2}$, higher than expected by light element nucleosynthesis; the model also predicts that most of the baryons at $z>2$ are in $\lya$ clouds, and that the rate at which the baryons move to more overdense regions is slow.

Cosmological simulations with gas provide a detailed description of the intergalactic medium, making possible predictions of neutral hydrogen absorption in the spectra of background QSOs. We present results from a high-resolution calculation of an $\Omega=1$ cold dark matter model. Our simulation reproduces many of the observed properties of the \lya forest surprisingly well. The distribution of HI column densities agrees with existing data to within a factor of $\sim$ two over most of the range from $10^{14}\cm^{-2}$ to $10^{22}\cm^{-2}$; i.e., from unsaturated \lya forest lines to damped \lya systems. The equivalent width distribution matches the observed exponential form with a characteristic width $W_* \approx 0.3$ angstroms. The distribution of b-parameters appears consistent with that derived from QSO spectra. Most of the low column density absorption arises in large, flattened structures of moderate or even relatively low overdensity, so there is no sharp distinction between the \lya forest and the ``Gunn-Peterson'' absorption produced by the smooth intergalactic medium. Our results demonstrate that a \lya forest like that observed develops naturally in a hierarchical clustering scenario with a photoionizing background. Comparison between simulations and high-resolution QSO spectra should open a new regime for testing theories of cosmic structure formation.

We study the formation of damped \lya and Lyman limit absorbers in a hierarchical clustering scenario using a gas dynamical simulation of an $\Omega = 1$, cold dark matter universe. In the simulation, these high column density systems are associated with forming galaxies. Damped \lya absorption, $N_{HI} \simgt 10^{20.2}\cm^{-2}$, arises along lines of sight that pass near the centers of relatively massive, dense protogalaxies. Lyman limit absorption, $10^{17}\cm^{-2} \simlt N_{HI} \simlt 10^{20.2}\cm^{-2}$, develops on lines of sight that pass through the outer parts of such objects or near the centers of smaller protogalaxies. The number of Lyman limit systems is less than observed, while the number of damped \lya systems is quite close to the observed abundance. Damped absorbers are typically $\sim 10$ kpc in radius, but the population has a large total cross section because the systems are much more numerous than present day $L_*$ galaxies. Our results demonstrate that high column density systems like those observed arise naturally in a hierarchical theory of galaxy formation and that it is now possible to study these absorbers directly from numerical simulations.

This paper is intended as an introduction to the theory of weak lensing. A review of the inversion formula introduced by Kaiser and Squires is presented. We then prove the formula of the aperture densitometry method in a simple way that allows a clear understanding of where the various terms come from. This is particularly useful to quantitatively measure masses in any region of a lens. We then summarize what has been learned from observations of strong lensing about the dark matter distribution; weak lensing should provide similar information on larger scales in clusters of galaxies.

A mass model of the rich cluster A2218 is presented based on optical and new near-infrared images of arc(let)s and the ROSAT/HRI X-ray map. A lensing model is proposed with two dark matter clumps, centered on the two bright galaxies around which the various arcs are observed. Two of them, close to the most massive clump around the brightest galaxy, are required to be fold images, given their position angle with respect to the major axis. A third image for each of the two arcs is then predicted to be on the opposite side of the major axis. Both images are indeed found at the expected positions, with consistent magnitudes, colors and shapes. One of the arcs has a known redshift of $0.7$, and the model predicts that the other arc should be at higher redshift, $z \approx 3$. Photometric data in UBrzJK' indicates that a redshift close to 3 is indeed probable for this arc. The dark matter halo is found to be sharply peaked on the giant central galaxy. The ROSAT/HRI X-ray map shows a strongly clumped image, poorly correlated with the cD light and the mass model on scales below $\sim$ 200$h_{50}^{-1}$kpc, but similar to the distribution of galaxies and the mass model when smoothed on this scale. The large degree of central substructure implied helps to account for the discrepancy previously found (in the two massive clusters A2218 \& A1689) between the mass derived from the gravitational lenses and that derived from the X-ray profile and temperature distributions, in a model of hydrostatic equilibrium. Most of the discrepancy can probably be explained by bulk motions and acceleration of the gas, and projection effects enhancing the lensing power of the cluster.

We use the Zel'dovich approximation and another analytical approximation to calculate the evolution under gravitational instability of the underdense regions of a photoionized intergalactic medium (IGM). We find that over most of the spectrum of a quasar, the optical depth to $\lya$ scattering originates from gas in underdense regions, or voids. This causes the ratio of the median Gunn-Peterson (GP) absorption to the value for a uniform medium containing all the baryons in the universe to be very small, and to decrease as gravitational collapse proceeds. We apply our calculations to the observations of the intensity distribution in a $z=4.11$ quasar by Webb and coworkers. We show that if $\lya$ clouds arise from gravitational collapse, their observations must be interpreted as the first detection of the (fluctuating) GP effect, with a median value $\tau_{GP}\simeq 0.06$ at $z=4$. For typical low-density ($\Omega\sim 0.4$) cosmological models, this is consistent with the predicted baryon density from primordial nucleosynthesis, and the intensity of the ionizing background derived from the proximity effect. From the numerical simulations of Cen \etal, such models also predict correctly the number of $\lya$ absorption lines observed.

If gravitational clustering is a hierarchical process, the present large-scale structure of the galaxy distribution implies that structures on smaller scales must have formed at high redshift. We simulate the formation of small-scale structure (average cell mass: $\Delta \bar m_b=10^{4.2}$M$_\odot$) and the evolution of photoionized gas, in the specific case of a CDM model with a cosmological constant. The photoionized gas has a natural minimal scale of collapse, the Jeans scale ($m_{b,J}\simeq 10^{9}$M$_\odot$). We find that low column density ($\nhi \le 10^{14}\cm^{-2}$) lines originate in regions resembling Zel'dovich pancakes, where gas with overdensities in the range $3 - 30$ is enclosed by two shocks but is typically re-expanding at approximately the Hubble velocity. However, higher column density ($\nhi \ge 10^{15}\cm^{-2}$) lines stem from more overdense regions where the shocked gas is cooling. We show that this model can probably account for the observed number of lines, their distribution in column density and b-parameters, as well as the cloud physical sizes as observed in gravitationally lensed quasars. We find a redshift evolution that is

Long arcs in clusters of galaxies, produced by gravitational lensing, can be used to estimate the mass interior to the arcs and therefore, constrain the cluster mass distribution. The radial density distribution of the intracluster gas (ICM) can be extracted from the X-ray surface brightness observations. If the gas temperature is also known, it is then possible to probe the dynamical state of the gas and test whether the ICM is in hydro- static equilibrium within the gravitational potential of the cluster as a result of thermal pressure support. We analyze three clusters that exhibit large arcs, whose X-ray surface brightness profiles have been observed, and whose gas temperatures have been determined. In two of the clusters, A2218 and A1689, the central mass implied by lensing is a factor of $2$--$2.5$ too large for the gas at the observed temperature to be in hydrostatic equilibrium solely due to thermal pressure support. In other words, if we accept the mass estimate derived from the lensing analysis and demand that the X-ray surface brightness profile be consistent with the observations, the required gas temperature is a factor of $2$--$2.5$ higher than observed. The results for the third cluster, A2163 (the most luminous and the hottest cluster known), are more ambiguous. The discrepancy between the X-ray and the lensing mass estimates arise because the presence of arcs imply a highly concentrated cluster mass distribution whereas the observed X-ray profiles imply a more extended mass distribution. The large X-ray core radii are not the result of the limited resolution of the X-ray detectors. We consider various possibilities that could account for the discrepancy.

We examine the constraints on theories of galaxy formation that are obtained from observations of damped $\lya$ (DL) systems, assuming they are gaseous protodisks in dark matter halos. Using the Press-Schechter formalism, we find that the mixed dark matter model, with $\ohdm = 0.3$, $\ocdm = 0.65$, $\obaryon = 0.05$, and $h=0.5$, is ruled out because the number of galactic halos at $z\simeq 3$ is too small to account for the total gaseous mass in DL systems, even under the assumption that all the gas in collapsed halos has settled into disks of neutral gas. The standard CDM model can account for the gas in DL systems if the bias is $b\lsim 2$; the same is true for the CDM model with a cosmological constant, if $b\lsim 1.5$ for $\Lambda = 0.8$. However, one still needs to assume that a fraction $\gsim 0.4$ of the baryons in collapsed halos at $z\simeq 3$ is in the form of neutral gas in disks. We also calculate the column density distribution $f(\nhi)$ of the DL systems, in terms of the surface density profiles of disks and the distribution of their central column densities. It is shown that the form of $f(\nhi)$ at the high end of column density is a diagnostic for the nature of DL systems.