If the dark matter is made of ultra-light axions, stable solitonic cores form at the centers of virialized halos. In some range for the mass $m$ of the axion particle, these cores are sufficiently compact and can mimic supermassive black holes (SMBH) residing at galactic nuclei. We use the solitonic core--halo mass relation, validated in numerical simulations, to constrain a new range of allowed axion mass from measurements of the SMBH mass in (pseudo)bulge and bulgeless galaxies. These limits are based on observations of galactic nuclei on scales smaller than 10 pc. Our analysis suggests that $m < 10^{-18}$ eV is ruled out by the data. We briefly discuss whether an attractive self-interaction among axions could alleviate this constraint.

A scenario for achieving a low velocity dispersion for the galaxy [KKS 2000]04 (aka NGC 1052-DF2) and similar galaxies is presented. A progenitor halo corresponding to a $z=0$ halo of mass $\sim 5\times 10^{10}\; \textrm{M}_\odot$ and a low concentration parameter (but consistent with cosmological simulations) infalls onto a Milky Way-size host at early times. {Substantial removal of cold gas} from the inner regions by supernova feedback and ram pressure, assisted by tidal stripping of the dark matter in the outer regions, leads to a substantial reduction of the velocity dispersion of stars within one effective radius. In this framework, the observed stellar content of [KKS 2000]04 is associated with a progenitor mass close to that inferred from the global stellar-to-halo-mass ratio. As far as the implications of kinematics are concerned, even if at a $\sim 20 $ Mpc distance, it is argued that [KKS 2000]04 is no more peculiar than numerous early type galaxies with seemingly little total dark-matter content.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Recently, \cite{vanDokkum2018} have presented an important discovery of an ultra diffuse galaxy, NGC1052-DF2, with a dark matter content significantly less than predicted from its stellar mass alone. The analysis relies on measured radial velocities of 10 Globular Clusters (GCs), of estimated individual masses of a few $ \times 10^6 M_\odot$. This is about $1\%$ of the inferred mass of NGC1052-DF2 of $2\times 10^8 M_\odot$ within a half-light radius, $R_\mathrm{e}=2.2\, \mathrm{kpc}$. The large relative mass and the old age of these objects imply that they might be susceptible to orbital decay by dynamical friction. Using analytic estimates and N-body simulations of an isolated system matching the inferred mass profile of NGC1052-DF2, we show that orbits of the most massive GCs should already have decayed on a time scale of a few Gyrs. These findings should help in constraining mass profile and formation scenarios of NGC1052-DF2.

It is demonstrated that the kinematics of the 10 star clusters in NGC2052-DF2 is compatible with a high dynamical mass close to those implied by the standard stellar-to-halo-mass ratio (SHMR). The analysis relies on a convenient form for the distribution function (DF) of projected phase space data, capturing non-gaussian features in the spread of true velocities of the mass tracers. A key ingredient is tidal stripping by the gravity of the apparently nearby larger galaxy, NGC 1052. Tidal stripping decreases the range of velocities of mass tracers, while only mildly lowering the total mass inside the trimming radius $r_{tr}$. The analysis is performed assuming halo profiles consistent with simulations of the $\Lambda$CDM model. For the fiducial value $r_{tr}=10$ kpc, we find that the virial mass of the pre-trimmed halo is $M<1.6\times 10^{10}M_\odot$ at $2\sigma $ ($95\%$) and $M<8.6\times 10^{9}M_\odot$ at $1.64\sigma$ ($90\%$). For the mass within 10 kpc we obtain, $M_\mathrm{10kpc}<3.9\times 10^{9}M_\odot$ and $<2.9\times 10^{9}M_\odot$ at $2\sigma$ and $1.64\sigma$, respectively. The $2\sigma$ upper limit on the virial mass is roughly a factor of 3-5 below the mean SHMR relation.Taking $r_{tr}=20$ kpc, lowers the $2\sigma$ virial mass limits by a factor of $\sim 4 $, bringing our results closer to those of Wasserman et al. (2018) without their SHMR prior.

Using cosmological N-body simulations, we study the abundance of local maxima (peaks) and minima (dips) identified in the smoothed distribution of halos and dark matter (DM) on scales of $10-100$s Mpcs. The simulations include Gaussian and local-type $f_{\rm NL}$ non-Gaussian initial conditions. The expression derived in the literature for the abundance (irrespective of height) of peaks for Gaussian fields is surprisingly accurate for the evolved halo and DM density fields for all initial conditions considered. Furthermore, the height distribution is very well fitted by a log-normal on quasi-linear scales. The abundance as a function of scale depends on the cosmological parameters ($H_0$ and background matter densities) through the shape of the power spectrum, but it is insensitive to the clustering amplitude. Further, the abundance in the smoothed halo distribution is substantially different in the non-Gaussian from the Gaussian simulations. The interpretation of this effect is straightforward in terms of the scale dependence of halo bias in non-Gaussian models. The abundance of extrema extracted from three-dimensional large galaxy redshift surveys could be a competitive probe of the cosmological parameters and initial non-Gaussianity. It breaks the degeneracy between $f_{\rm NL}$ and the clustering amplitude, making it complementary to counts of galaxy clusters and peaks in weak-lensing maps.

The peculiar velocity of a mass tracer is on average aligned with the dipole modulation of the surrounding mass density field. We present a first measurement of the correlation between radial peculiar velocities of objects in the cosmicflows-3 catalog and the dipole moment of the 2MRS galaxy distribution in concentric spherical shells centered on these objects. Limiting the analysis to cosmicflows-3 objects with distances of $100 \rm Mpc h^{-1}$, the correlation function is detected at a confidence level $> 4\sigma$. The measurement is found consistent with the standard $\Lambda$CDM model at $< 1.7\sigma$ level. We formally derive the constraints $0.32<\Omega^{0.55}\sigma_8<0.48$ ($68\% $ confidence level) or equivalently $0.34<\Omega^{0.55}/b<0.52$, where $b$ is the galaxy bias factor. Deeper and improved peculiar velocity catalogs will substantially reduce the uncertainties, allowing tighter constraints from this type of correlations.

We develop the framework for testing Lorentz invariance in the dark matter sector using galactic dynamics. We consider a Lorentz violating (LV) vector field acting on the dark matter component of a satellite galaxy orbiting in a host halo. We introduce a numerical model for the dynamics of satellites in a galactic halo and for a galaxy in a rich cluster to explore observational consequences of such an LV field. The orbital motion of a satellite excites a time dependent LV force which greatly affects its internal dynamics. Our analysis points out key observational signatures which serve as probes of LV forces. These include modifications to the line of sight velocity dispersion, mass profiles and shapes of satellites. With future data and a more detailed modeling these signatures can be exploited to constrain a new region of the parameter space describing the LV in the dark matter sector.

Observed galaxy luminosities (derived from redshifts) hold information on the large-scale peculiar velocity field in the form of spatially correlated scatter, which allows for bounds on bulk flows and the growth rate of matter density perturbations using large galaxy redshift surveys. We apply this luminosity approach to galaxies from the recent SDSS Data Release 13. Our goal is twofold. First, we take advantage of the recalibrated photometry to identify possible systematic errors relevant to our previous analysis of earlier data. Second, we seek improved constraints on the bulk flow and the normalized growth rate f$\sigma_{8}$ at z ~ 0.1. Our results confirm the robustness of our method. Bulk flow amplitudes, estimated in two redshift bins with 0.02 < z$_{1}$ < 0.07 < z$_{2}$ < 0.22, are generally smaller than in previous measurements, consistent with both the updated photometry and expectations for the $\Lambda$CDM model. The obtained growth rate, f$\sigma_{8}$ = 0.48 +/- 0.16, is larger than, but still compatible with, its previous estimate, and closer to the reference value of Planck. Rather than precision, the importance of these results is due to the fact that they follow from an independent method that relies on accurate photometry, which is a top requirement for next-generation photometric catalogs.

We assess the effect of the local large scale structure on the estimation of two-point statistics of the observed radial peculiar velocities of galaxies. A large N-body simulation is used to examine these statistics from the perspective of random observers as well as "Local Group (LG)-like" observers conditioned to reside in an environment resembling the observed universe within 20 Mpc. The local environment systematically distorts the shape and amplitude of velocity statistics with respect to ensemble-averaged measurements made by a Copernican (random) observer. The Virgo cluster has the most significant impact, introducing large systematic deviations in all the statistics. For a simple "top-hat" selection function, an idealized survey extending to $\sim 160h^{-1}\,{\rm Mpc}$ or deeper is needed to completely mitigate the effects of the local environment. Using shallower catalogues leads to systematic deviations of the order of $50$ to $200\%$ depending on the scale considered. For a flat redshift distribution similar to the one of the CosmicFlows-3 survey, the deviations are even more prominent in both the shape and amplitude at all separations considered $({\stackrel{<}{{}_\sim}} 100h^{-1}\,{\rm Mpc})$. Conclusions based on statistics calculated without taking into account the impact of the local environment should be revisited.

NoAM for "No Action Method" is a framework for reconstructing the past orbits of observed tracers of the large scale mass density field. It seeks exact solutions of the equations of motion (EoM), satisfying initial homogeneity and the final observed particle (tracer) positions. The solutions are found iteratively reaching a specified tolerance defined as the RMS of the distance between reconstructed and observed positions. Starting from a guess for the initial conditions, NoAM advances particles using standard N-body techniques for solving the EoM. Alternatively, the EoM can be replaced by any approximation such as Zel'dovich and second order perturbation theory (2LPT). NoAM is suitable for billions of particles and can easily handle non-regular volumes, redshift space, and other constraints. We implement NoAM to systematically compare Zel'dovich, 2LPT, and N-body dynamics over diverse configurations ranging from idealized high-res periodic simulation box to realistic galaxy mocks. Our findings are (i) Non-linear reconstructions with Zel'dovich, 2LPT, and full dynamics perform better than linear theory only for idealized catalogs in real space. For realistic catalogs, linear theory is the optimal choice for reconstructing velocity fields smoothed on scales > 5 Mpc/h. (ii) all non-linear back-in-time reconstructions tested here, produce comparable enhancement of the baryonic oscillation signal in the correlation function.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The goodness-of-fit (GoF) of the Cosmicflows-2 (CF2) database of peculiar velocities with the LCDM standard model of cosmology is presented. Standard application of the Chi^2 statistics of the full database, of its 4,838 data points, is hampered by the small scale nonlinear dynamics which is not accounted for by the (linear regime) velocity power spectrum. The bulk velocity constitutes a highly compressed representation of the data which filters out the small scales non-linear modes. Hence the statistics of the bulk flow provides an efficient tool for assessing the GoF of the data given a model. The particular approach introduced here is to use the (spherical top-hat window) bulk velocity extracted from the Wiener filter reconstruction of the 3D velocity field as a linear low pass filtered highly compressed representation of the CF2 data. An ensemble 2250 random linear realizations of the WMAP/LCDM model has been used to calculate the bulk velocity auto-covariance matrix. We find that the CF2 data is consistent with the WMAP/LCDM model to better than the 2 sigma confidence limits. This provides a further validation that the CF2 database is consistent with the standard model of cosmology.

We present a realistic modeling of the dipole component of the projected sky distribution of NVSS radio galaxies. The modeling relies on mock catalogs generated within the context of $\Lambda$CDM cosmology, in the linear regime of structure formation. After removing the contribution from the solar motion, the mocks show that the remaining observed signal is mostly (70\%) due to structures within $z\lesssim0.1$. The amplitude of the model signal depends on the bias factor $b$ of the NVSS mock galaxies. For sources with flux density, $ S> 15 \; \rm mJy$, the bias recipe inferred from higher order moments is consistent with the observed dipole signal at $2.12\sigma$. Flux thresholds above $20 \; \rm mJy$ yield a disagreement close to the $3\sigma $ level. A constant high bias, $b=3$ is needed to mitigate the tension to the $\sim 2.3\sigma$ level.

An interesting test of Einstein's equivalence principle (EEP) relies on the observed lag in arrival times of photons emitted from extragalactic transient sources. Attributing the lag between photons of different energies to the gravitational potential of the Milky Way (MW), several authors derive new constraints on deviations from EEP. It is shown here that potential fluctuations from the large scale structure are at least two orders of magnitude larger than the gravitational potential of the MW. Combined with the larger distances, for sources at redshift $z\gtsim 0.5$ the {\it rms} of the contribution from these fluctuations exceeds the MW by more than 4 orders of magnitude. We provide actual constraints for several objects based on a statistical calculation of the large scale fluctuations in the standard $\Lambda$CDM cosmological model.

We explore similarities and differences between several estimators of the cosmological bulk flow, $\bf B$, from the observed radial peculiar velocities of galaxies. A distinction is made between two theoretical definitions of $\bf B$ as a dipole moment of the velocity field weighted by a radial window function. One definition involves the three dimensional (3D) peculiar velocity, while the other is based on its radial component alone. Different methods attempt at inferring $\bf B$ for either of these definitions which coincide only for a constant velocity field. We focus on the Wiener Filtering (WF, Hoffman et al. 2015) and the Constrained Minimum Variance (CMV,Feldman et al. 2010) methodologies. Both methodologies require a prior expressed in terms of the radial velocity correlation function. Hoffman et al. compute $\bf B$ in Top-Hat windows from a WF realization of the 3D peculiar velocity field. Feldman et al. infer $\bf B$ directly from the observed velocities for the second definition of $\bf B$. The WF methodology could easily be adapted to the second definition, in which case it will be equivalent to the CMV with the exception of the imposed constraint. For a prior with vanishing correlations or very noisy data, CMV reproduces the standard Maximum Likelihood (ML, Kaiser 1988) estimation for $\bf B$ of the entire sample independent of the radial weighting function. Therefore, this estimator is likely more susceptible to observational biases that could be present in measurements of distant galaxies. Finally, two additional estimators are proposed.

Spatial variations in the distribution of galaxy luminosities, estimated from redshifts as distance proxies, are correlated with the peculiar velocity field. Comparing these variations with the peculiar velocities inferred from galaxy redshift surveys is a powerful test of gravity and dark energy theories on cosmological scales. Using ~ 2 $\times$ 10$^{5}$ galaxies from the SDSS Data Release 7, we perform this test in the framework of gravitational instability to estimate the normalized growth rate of density perturbations f$\sigma_{8}$ = 0.37 +/- 0.13 at z ~ 0.1, which is in agreement with the $\Lambda$CDM scenario. This unique measurement is complementary to those obtained with more traditional methods, including clustering analysis. The estimated accuracy at z ~ 0.1 is competitive with other methods when applied to similar datasets.

We study the angular clustering of $\sim 6\times 10^5$ NVSS sources on scales $\gtrsim 50 h^{-1}$ Mpc in the context of the $\Lambda$CDM scenario. The analysis partially relies on the redshift distribution of 131 radio galaxies, inferred from the Hercules and CENSORS survey, and an empirical fit to the stellar to halo mass (SHM) relation. For redshifts $z\lesssim 0.7$, the fraction of radio activity versus stellar mass evolves as $f_{_{\rm RL}}\sim M_*^{\alpha_0+\alpha_1 z}$ where $\alpha_0=2.529{\pm0.184}$ and $\alpha_1=1.854^{+0.708}_{-0.761}$. The estimate on $\alpha_0$ is largely driven by the results of Best et al. (2005), while the constraint on $\alpha_1$ is new. We derive a biasing factor $b(z=0.5)=2.093^{+0.164}_{-0.109}$ between radio galaxies and the underlying mass.The function $b(z)=0.33z^2+0.85 z +1.6$ fits well the redshift dependence. We also provide convenient parametric forms for the redshift dependent radio luminosity function, which are consistent with the redshift distribution and the NVSS source count versus flux.

Comparison of galaxy flows with those predicted from the local galaxy distribution ended as an active field after two analyses came to vastly different conclusions 25 years ago, but that was due to faulty data. All the old results are therefore suspect. With new data collected in the last several years, the problem deserves another look. For this we analyze the gravity field inferred from the enormous data set derived from the 2MASS collection of galaxies, and compare it to the velocity field derived from the well calibrated SFI++ Tully-Fisher catalog. Using the "Inverse Method" to minimize Malmquist biases, within 10,000 km/s the gravity field is seen to predict the velocity field to remarkable consistency. This is a beautiful demonstration of linear perturbation theory and is fully consistent with standard values of the cosmological variables.

The phase space distribution of matter out to $\sim 100 \rm Mpc$ is probed by two types of observational data: galaxy redshift surveys and peculiar motions of galaxies. Important information on the process of structure formation and deviations from standard gravity have been extracted from the accumulating data. The remarkably simple Zel'dovich approximation is the basis for much of our insight into the dynamics of structure formation and the development of data analyses methods. Progress in the methodology and some recent results is reviewed.

Peculiar motion introduces systematic variations in the observed luminosity distribution of galaxies. This allows one to constrain the cosmic peculiar velocity field from large galaxy redshift surveys. Using around half a million galaxies from the SDSS Data Release 7 at z ~ 0.1, we demonstrate the applicability of this approach to large datasets and obtain bounds on peculiar velocity moments and $\sigma_{8}$, the amplitude of the linear matter power spectrum. Our results are in good agreement with the $\Lambda$CDM model and consistent with the previously reported ~ 1% zero-point tilt in the SDSS photometry. Finally, we discuss the prospects of constraining the growth rate of density perturbations by reconstructing the full linear velocity field from the observed galaxy clustering in redshift space.

Spatial modulations in the distribution of observed luminosities (computed using redshifts) of ~ 5 $\times$ 10$^5$ galaxies from the SDSS Data Release 7, probe the cosmic peculiar velocity field out to z ~ 0.1. Allowing for luminosity evolution, the r-band luminosity function, determined via a spline-based estimator, is well represented by a Schechter form with M$^{\star}$(z) - 5log$_{10}$h = -20.52 - 1.6(z - 0.1) $\pm$ 0.05 and $\alpha^{\star}$ = -1.1 $\pm$ 0.03. Bulk flows and higher velocity moments in two redshift bins, 0.02 < z < 0.07 and 0.07 < z < 0.22, agree with the predictions of the $\Lambda$CDM model, as obtained from mock galaxy catalogs designed to match the observations. Assuming a $\Lambda$CDM model, we estimate $\sigma_{8}$ $\approx$ 1.1 $\pm$ 0.4 for the amplitude of the linear matter power spectrum, where the low accuracy is due to the limited number of galaxies. While the low-z bin is robust against coherent photometric uncertainties, the bias of results from the second bin is consistent with the ~ 1% magnitude tilt reported by the SDSS collaboration. The systematics are expected to have a significantly lower impact in future datasets with larger sky coverage and better photometric calibration.

Maximum Likelihood estimation of the bulk flow from radial peculiar motions of galaxies, generally assumes a constant velocity field inside the survey volume. The assumption is inconsistent with the definition of the bulk flow as the average of the peculiar velocity field over the relevant volume. This follows from a straightforward mathematical relation between the bulk flow of a sphere and the velocity potential on its surface. The inconsistency exists also for ideal data with exact radial velocities and full spatial coverage. Based on the same relation we propose a simple modification to correct for this inconsistency.

There is a $\sim 150 km s^{-1}$ discrepancy between the measured motion of the Local Group of galaxies (LG) with respect to the CMB and the linear theory prediction based on the gravitational force field of the large scale structure in full-sky redshift surveys. We perform a variety of tests which show that the LG motion cannot be recovered to better than $150-200 km s^{-1}$ in amplitude and within a $\approx10^\circ$ in direction. The tests rely on catalogs of mock galaxies identified in the Millennium simulation using semi-analytic galaxy formation models. We compare these results to the $K_s=11.75$ Two-Mass Galaxy Redshift Survey, which provides the deepest, widest and most complete spatial distribution of galaxies available so far. In our analysis we use a new, concise relation for deriving the LG motion and bulk flow from the true distribution of galaxies in redshift space. Our results show that the main source of uncertainty is the small effective depth of surveys like the 2MRS that prevents a proper sampling of the large scale structure beyond $\sim100 h^{-1} Mpc$. Deeper redshift surveys are needed to reach the "convergence scale" of $\approx 250 h^{-1}Mpc $ in a $\Lambda$CDM universe. Deeper survey would also mitigate the impact of the "Kaiser rocket" which, in a survey like 2MRS, remains a significant source of uncertainty. Thanks to the quiet and moderate density environment of the LG, purely dynamical uncertainties of the linear predictions are subdominant at the level of $\sim 90 km s^{-1}$. Finally, we show that deviations from linear galaxy biasing and shot noise errors provide a minor contribution to the total error budget.

The observed density field in redshift space is directly affected by the radial motions generated from mass fluctuations outside the volume occupied by a given galaxy redshift survey. These motions introduce redshift space anisotropies which are more pronounced at larger distances from the survey's center, thus offering clues to the nature of mass fluctuations on super-survey scales. Furthermore, we note that all estimates of the growth factor derived from redshift space distortions are based on relations which explicitly assume that the velocity field is generated by mass fluctuations inside the survey volume. This may cause uncertainties in these estimates which are on the order of a few percent.

We constrain the most likely range of masses for the Milky Way and M31 using an application of the Numerical Action Method (NAM) that optimizes the fit to observed parameters over a large ensemble of NAM-generated solutions. Our 95% confidence level mass ranges, 1.5-4.5 x 10^12 m_sun for MW and 1.5-5.5 x 10^12 m_sun for M31, are consistent with the upper range of estimates from other methods and suggests that a larger proportion of the total mass becomes detectable when the peculiar motions of many nearby satellites are taken into account in the dynamical analysis. We test the method against simulated Local Group catalogs extracted from the Millennium Run to confirm that mass predictions are consistent with actual galaxy halo masses.

The cosmological peculiar velocity field (deviations from the pure Hubble flow) of matter carries significant information on dark energy, dark matter and the underlying theory of gravity on large scales. Peculiar motions of galaxies introduce systematic deviations between the observed galaxy redshifts z and the corresponding cosmological redshifts z_cos. A novel method for estimating the angular power spectrum of the peculiar velocity field based on observations of galaxy redshifts and apparent magnitudes m (or equivalently fluxes) is presented. This method exploits the fact that a mean relation between z_cos and m of galaxies can be derived from all galaxies in a redshift-magnitude survey. Given a galaxy magnitude, it is shown that the z_cos(m) relation yields its cosmological redshift with a 1-sigma error of sigma_z~0.3 for a survey like Euclid (~10^9 galaxies at z<~2), and can be used to constrain the angular power spectrum of z-z_cos(m) with a high signal-to-noise ratio. At large angular separations corresponding to l<~15, we obtain significant constraints on the power spectrum of the peculiar velocity field. At 15<~l<~60, magnitude shifts in the z_cos(m) relation caused by gravitational lensing magnification dominate, allowing us to probe the line-of-sight integral of the gravitational potential. Effects related to the environmental dependence in the luminosity function can easily be computed and their contamination removed from the estimated power spectra. The amplitude of the combined velocity and lensing power spectra at z~1 can be measured with <~5% accuracy.

Using redshifts as a proxy for galaxy distances, estimates of the 2D transverse peculiar velocities of distant galaxies could be obtained from future measurements of proper motions. We provide the mathematical framework for analyzing 2D transverse motions and show that they offer several advantages over traditional probes of large scale motions. They are completely independent of any intrinsic relations between galaxy properties, hence they are essentially free of selection biases. They are free from homogeneous and inhomogeneous Malmquist biases that typically plague distance indicator catalogs. They provide additional information to traditional probes which yield line-of-sight peculiar velocities only. Further, because of their 2D nature, fundamental questions regarding vorticity of large scale flows can be addressed. Gaia for example is expected to provide proper motions of at least bright galaxies with high central surface brightness, making proper motions a likely contender traditional probes based on current and future distance indicator measurements.

It is widely accepted that within the framework of LCDM a significant fraction of giant-disk galaxies has recently experienced a violent galactic merger. We present numerical simulations of such major mergers of gas-rich pure disk galaxies, and focus on the innermost stellar component (bulge) of the disk remnants. The simulations have high spatial and mass resolutions, and resolve regions deep enough to allow bulge classification according to standard kinematical and structural characteristics. In agreement with recent studies we find that these bulges are dominated by stars formed in the final coalescence process. In contrast to the common interpretation of such components as classical bulges (i.e. similar to intermediate luminosity ellipticals), we find they are supported by highly coherent rotations and have Sersic indices n<2, a result leading to their classification as pseudo-bulges. Pseudo-bulge formation by gas rich major mergers of pure disks is a novel mode of pseudo-bulge formation; It complements pseudo-bulge growth by secular evolution, and it could help explain the high fractions of classically bulge-less giant disk galaxies, and pseudo-bulges found in giant Sc galaxies.

Using the nearly full sky Ks=11.75 2MASS Redshift Survey [2MRS]of ~45,000 galaxies we reconstruct the underlying peculiar velocity field and constrain the cosmological bulk flow within ~100. These results are obtained by maximizing the probability to estimate the absolute magnitude of a galaxy given its observed apparent magnitude and redshift. At a depth of ~60 Mpc/h we find a bulk flow Vb=(90\pm65,-230\pm65,50\pm65) km/s in agreement with the theoretical predictions of the LCDM model. The reconstructed peculiar velocity field that maximizes the likelihood is characterized by the parameter beta=0.323 +/- 0.08. Both results are in agreement with those obtained previously using the ~23,000 galaxies of the shallower Ks=11.25 2MRS survey. In our analysis we find that the luminosity function of 2MRS galaxies is poorly fitted by the Schechter form and that luminosity evolves such that objects become fainter with increasing redshift according to L(z)=L(z=0)(1+z)^(+2.7 +/-0.15).

Given a redshift survey of galaxies with measurements of apparent magnitudes, we present a novel method for measuring the growth rate $f(\Omega)$ of cosmological linear perturbations. We use the galaxy distribution within the survey to solve for the peculiar velocity field which depends in linear perturbation theory on $\beta=f(\Omega)/b$, where $b$ is the bias factor of the galaxy distribution. The recovered line-of-sight peculiar velocities are subtracted from the redshifts to derive the distances, which thus allows an estimate of the absolute magnitude of each galaxy. A constraint on $\beta$ is then found by minimizing the spread of the estimated magnitudes from their distribution function. We apply the method to the all sky $K = 11.25$ Two-MASS Redhsift Survey (2MRS) and derive $\beta=0.35\pm 0.1$ at $z\sim 0$, remarkably consistent with our previous estimate from the velocity-velocity comparison. The method could easily be applied to subvolumes extracted from the SDSS survey to derive the growth rate at $z \sim 0.1$. Further, it should also be applicable to ongoing and future spectroscopic redshift surveys to trace the evolution of $f(\Omega)$ to $z\sim1$. Constraints obtained from this method are entirely independent from those obtained from the two-dimensional distortion of $\xi(s)$ and provide an important check on $f(\Omega)$, as alternative gravity models predict observable differences.

We present a simple method for measuring cosmological bulk flows from large redshift surveys, based on the apparent dimming or brightening of galaxies due to their peculiar motion. It is aimed at estimating bulk flows of cosmological volumes containing large numbers of galaxies. Constraints on the bulk flow are obtained by minimizing systematic variations in galaxy luminosities with respect to a reference luminosity function measured from the whole survey. This method offers two advantages over more popular bulk flow estimators: it is independent of error-prone distance indicators and of the poorly-known galaxy bias. We apply the method to the 2MASS redshift survey (2MRS) to measure the local bulk flows of spherical shells centered on the Milky Way (MW). The result is consistent with that obtained by Nusser and Davis (2011) using the SFI++ catalogue of Tully-Fisher distance indicators. We also make an assessment of the ability of the method to constrain bulk flows at larger redshifts ($z=0.1-0.5$) from next generation datasets. As a case study we consider the planned EUCLID survey. Using this method we will be able to measure a bulk motion of $ \sim 200 \kms$ of $10^6$ galaxies with photometric redshifts, at the $3\sigma$ level for both $z\sim 0.15$ and $z\sim 0.5$. Thus the method will allow us to put strong constraints on dark energy models as well as alternative theories for structure formation.

We derive estimates for the cosmological bulk flow from the SFI++ Tully-Fisher (TF) catalog. For a sphere of radius $40 \hmpc$ centered on the MW, we derive a bulk flow of $333 \pm 38\kms $ towards Galactic $ (l,b)=(276^\circ,14^\circ)$ within a $3^\circ$ $1\sigma$ error. Within a $ 100\hmpc$ we get $ 257\pm 44\kms$ towards $(l,b)=(279^\circ, 10^\circ)$ within a $6^\circ$ error. These directions are at a $40^\circ$ with the Supergalactic plane, close to the apex of the motion of the Local Group of galaxies after the Virgocentric infall correction. Our findings are consistent with the $\Lambda$CDM model with the latest WMAP best fit cosmological parameters. But the bulk flow allows independent constraints. For WMAP inferred Hubble parameter $h=0.71$ and baryonic mean density parameter $\Omega_b=0.0449$, the constraint from the bulk flow on the matter density $\Omega_m$, the normalization of the density fluctuations, $\sigma_8$, and the growth index, $\gamma$, can be expressed as $\sigma_8\Omega_m^{\gamma-0.55}(\Omega_m/0.266)^{0.28}=0.86\pm 0.11$ (for $\Omega_m\approx 0.266$). Fixing $\sigma_8=0.8$ and $\Omega_m=0.266$ as favored by WMAP, we get $\gamma=0.495\pm 0.096$. The constraint derived here rules out popular DGP models at more than the 99% confidence level. Our results are based on a method termed \ace\ (All Space Constrained Estimate) which reconstructs the bulk flow from an all space three dimensional peculiar velocity field constrained to match the TF measurements. At large distances \ace\ generates a robust bulk flow from the SFI++ that is insensitive to the assumed prior. For comparison, a standard straightforward maximum likelihood estimate leads to very similar results.

(abridged) We perform a reconstruction of the cosmological large scale flows in the nearby Universe using two complementary observational sets. The first, the SFI++ sample of Tully-Fisher (TF) measurements of galaxies, provides a direct probe of the flows. The second, the whole sky distribution of galaxies in the 2MASS redshift survey (2MRS), yields a prediction of the flows given the cosmological density parameter, $\Omega$, and a biasing relation between mass and galaxies. We aim at an unbiased comparison between the peculiar velocity fields extracted from the two data sets and its implication on the cosmological parameters and the biasing relation. We expand the fields in a set of orthonormal basis functions, each representing a plausible realization of a cosmological velocity field. Our analysis completely avoids the strong error covariance in the smoothed TF velocities by the use of orthonormal basis functions and employs elaborate realistic mock data sets to extensively calibrate the errors in 2MRS predicted velocities. We relate the 2MRS galaxy distribution to the mass density field by a linear bias factor, $b$, and include a luminosity dependent, $\propto L^\alpha$, galaxy weighting. We assess the agreement between the fields as a function of $\alpha$ and $\beta=f(\Omega)/b$, where $f$ is the growth factor of linear perturbations. The agreement is excellent with a reasonable $\chi^2$ per degree of freedom. For $\alpha=0$, we derive $0.28<\beta<0.37$ and $0.24<\beta<0.43$, respectively, at the 68.3% and 95.4% confidence levels (CLs). For $\beta=0.33$, we get $\alpha<0.25$ and $\alpha<0.5$, respectively, at the 68.3% and 95.4% CLs. We set a constraint on the fluctuation normalization, finding $\sigma_8 = 0.73 \pm 0.1$, in very good agreement with the latest WMAP results.

Force balance considerations put a limit on the rate of AGN radiation momentum output, $L/c$, capable of driving galactic superwinds and reproducing the observed $\mbh -\sigma $ relation between black hole mass and spheroid velocity dispersion. We show that black holes cannot supply enough momentum in radiation to drive the gas out by pressure alone. Energy-driven winds give a $\mbh -\sigma $ scaling favoured by a recent analysis but also fall short energetically once cooling is incorporated. We propose that outflow-triggering of star formation by enhancing the intercloud medium turbulent pressure and squeezing clouds can supply the necessary boost, and suggest possible tests of this hypothesis. Our hypothesis simultaneously can account for the observed halo baryon fraction.

The great advances in the network of cosmological tests show that the relativistic Big Bang theory is a good description of our expanding universe. But the properties of nearby galaxies that can be observed in greatest detail suggest a still better theory would more rapidly gather matter into galaxies and groups of galaxies. This happens in theoretical ideas now under discussion.

A long-range force acting only between nonbaryonic particles would be associated with a large violation of the weak equivalence principle. We explore cosmological consequences of this idea, which we label ReBEL (daRk Breaking Equivalence principLe). A high resolution hydrodynamical simulation of the distributions of baryons and dark matter confirms our previous findings that a ReBEL force of comparable strength to gravity on comoving scales of about 1 Mpc/h causes voids between the concentrations of large galaxies to be more nearly empty, suppresses accretion of intergalactic matter onto galaxies at low redshift, and produces an early generation of dense dark matter halos. A preliminary analysis indicates the ReBEL scenario is consistent with the one-dimensional power spectrum of the Lyman-Alpha forest and the three-dimensional galaxy auto-correlation function. Segregation of baryons and DM in galaxies and systems of galaxies is a strong prediction of ReBEL. ReBEL naturally correlates the baryon mass fraction in groups and clusters of galaxies with the system mass, in agreement with recent measurements.

A method is proposed for constraining the Galactic gravitational potential from high precision observations of the phase space coordinates of a system of relaxed tracers. The method relies on an "ergodic" assumption that the observations are representative of the state of the system at any other time. The observed coordinates serve as initial conditions for moving the tracers forward in time in an assumed model for the gravitational field. The validity of the model is assessed by the statistical equivalence between the observations and the distribution of tracers at randomly selected times. The applicability of this ergodic method is not restricted by any assumption on the form or symmetry of the potential. However, it requires high recision observations as those that will be obtained from missions like SIM and GAIA.

Numerical simulations of the effect of a long-range scalar interaction (LRSI) acting only on nonbaryonic dark matter, with strength comparable to gravity, show patterns of disruption of satellites that can agree with what is seen in the Milky Way. This includes the symmetric Sagittarius stellar stream. The exception presented here to the Kesden and Kamionkowski demonstration that an LRSI tends to produce distinctly asymmetric streams follows if the LRSI is strong enough to separate the stars from the dark matter before tidal disruption of the stellar component, and if stars dominate the mass in the luminous part of the satellite. It requires that the Sgr galaxy now contains little dark matter, which may be consistent with the Sgr stellar velocity dispersion, for in the simulation the dispersion at pericenter exceeds virial. We present other examples of simulations in which a strong LRSI produces satellites with large mass-to-light ratio, as in Draco, or free streams of stars, which might be compared to "orphan" streams.

We present a method for extracting the expected cosmological 21-cm signal from the epoch of reionization, taking into account contaminating radiations and random instrumental noise. The method is based on the maximum a-posteriori probability (MAP) formalism and employs the coherence of the contaminating radiation along the line-of-sight and the three-dimensional correlations of the cosmological signal. We test the method using a detailed and comprehensive modeling of the cosmological 21-cm signal and the contaminating radiation. The signal is obtained using a high resolution N-body simulation where the gas is assumed to trace the dark matter and is reionized by stellar radiation computed from semi-analytic galaxy formation recipes. We model contaminations to the cosmological signal from synchrotron and free-free galactic foregrounds and extragalactic sources including active galactic nuclei, radio haloes and relics, synchrotron and free-free emission from star forming galaxies, and free-free emission from dark matter haloes and the intergalactic medium. We provide tests of the reconstruction method for several rms values of instrumental noise from $\sigma_{N}=1$ to 250 mK. For low instrumental noise, the recovered signal, along individual lines-of-sight, fits the true cosmological signal with a mean rms difference of $d_{rms}\approx 1.7\pm 0.6$ for $\sigma_{N}=1$ mK, and $d_{rms}\approx 4.2\pm 0.4$ for $\sigma_{N}=5$ mK. The one-dimensional power spectrum is nicely reconstructed for all values of $\sigma_{N}$ considered here, while the reconstruction of the two-dimensional power spectrum and the Minkowski functionals is good only for noise levels of the order of few mK.

The time scale for cooling of the gravitationally bound gaseous intracluster medium (ICM) is not determined by radiative processes alone. If the ICM is in quasi-hydrostatic equilibrium in the fixed gravitational field of the dark matter halo then energy losses incurred by the gravitational potential energy of the gas should also be taken into account. This "gravitational heating" has been known for a while using explicit solutions to the equations of motion. Here, we re-visit this effect by applying the virial theorem to gas in quasi-hydrostatic equilibrium in an external gravitational field, neglecting the gravity of the gas. For a standard NFW form of halo profiles and for a finite gas density, the response of the gas temperature to changes in the total energy is significantly delayed. The effective cooling time could be prolonged by more than an order of magnitude inside the scale radius ($\rs$) of the halo. Gas lying at a distance twice the scale radius, has negative heat capacity so that the temperature increases as a result of energy losses. Although external heating (e.g. by AGN activity) is still required to explain the lack of cool ICM near the center, the analysis here may circumvent the need for heating in farther out regions where the effective cooling time could be prolonged to become larger than the cluster age and also explains the increase of temperature with radius in these regions.

The presence of dilute hot cavities in the intracluster medium (ICM) at the cores of clusters of galaxies changes the relation between gas temperature and its X-ray emission properties. Using the hydrostatic equations of a porous medium we solve for the ICM density for a given temperature as a function of the filling factor of dilute bubbles. We find that at a given temperature, the core X-ray luminosity increases with the filling factor. If the frequency of AGNs in clusters were higher in the past, then the filling factor could correspondingly be significant, with implications for the cluster scaling relations at high redshifts. This is especially important for the core properties, including the temperature-luminosity ($L_X-T$) relation and estimates of the core gas mass. The results imply an epoch-dependent sensitivity of the $L_X-T$ relation in the core to the porosity of the ICM. Detection of such an effect would give new insights into AGN feedback.

Weakly interacting massive particles (WIMPs) are a viable candidate for the relic abundance of dark matter (DM) produced in the early universe. So far WIMPs have eluded direct detection through interactions with baryonic matter. Neutrino emission from accumulated WIMP annihilations in the solar core has been proposed as a signature of DM, but has not yet been detected. These null results may be due to small scale DM density fluctuations in the halo with the density of our local region being lower than the average (around 0.3 GeV/cm^3 ). However, the accumulated neutrino signal from WIMP annihilations in the Galactic stellar disk would be insensitive to local density variations. Inside the disk, dark matter can be captured by stars causing an enhanced annihilation rate and therefore a potentially higher neutrino flux than what would be observed from elsewhere in the halo. We estimate a neutrino flux from the WIMP annihilations in the stellar disk to be enhanced by more than an order of magnitude compared to the neutrino fluxes from the halo. We offer a conservative estimate for this enhanced flux, based on the WIMP-nucleon cross-sections obtained from direct-detection experiments by assuming a density of around 0.3 GeV/cm^3 for the local DM. We also compare the detectability of these fluxes with a signal of diffuse high energy neutrinos produced in the Milky Way by the interaction of cosmic rays (CRs) with the interstellar medium (ISM). These comparative signals should be observable by large neutrino detectors.

Internal flows inside gravitationally stable astrophysical objects, such as the Sun, stars and compact stars are compressed and extremely subsonic. Such low Mach number flows are usually encountered when studying for example dynamo action in stars, planets, the hydro-thermodynamics of X-ray bursts on neutron stars and dwarf novae. Treating such flows is numerically complicated and challenging task. We aim to present a robust numerical tool that enables modeling the time-evolution or quasi-stationary of stratified low Mach number flows under astrophysical conditions. It is argued that astrophysical low Mach number flows cannot be considered as an asymptotic limit of incompressible flows, but rather as highly compressed flows with extremely stiff pressure terms. Unlike the pseudo-pressure in incompressible fluids, a Possion-like treatment for the pressure would smooth unnecessarily the physically induced acoustic perturbations, thereby violating the conservation character of the compressible equations. Moreover, classical dimensional splitting techniques, such as ADI or Line-Gauss-Seidel methods are found to be unsuited for modeling compressible flows with low Mach numbers. In this paper we present a nonlinear Newton-type solver that is based on the defect-correction iteration procedure and in which the Approximate Factorization Method (AFM) is used as a preconditioner. This solver is found to be sufficiently robust and is capable of capturing stationary solutions for viscous rotating flows with Mach number as small as $\mcal{M} \approx 10^{-3},$ i.e., near the incompressibility limit.

The dynamics of cosmological gravitating system is governed by the Euler and the Poisson equations. Tiny fluctuations near the big bang singularity are amplified by gravitational instability into the observed structure today. Given the current distribution of galaxies and assuming initial homogeneity, dynamical reconstruction methods have been developed to derive the cosmic density and velocity fields back in time. The reconstruction method described here is based a least action principle formulation of the dynamics of collisionless particle (representing galaxies). Two observational data sets will be considered. The first is the distribution of galaxies which is assumed to be an honest tracer of the mass density field of the dark matter. The second set is measurements of the peculiar velocities (deviations from pure Hubble flow) of galaxies. Given the first data set, the reconstruction method recovers the associated velocity field which can then be compared with the second data set. This comparison constrains the nature of the dark matter and the relation between mass and light in the Universe.

Assuming the dark matter is made entirely from neutralinos, we re-visit the role of their annihilation on the temperature of diffuse gas in the high redshift universe. We consider neutralinos of particle mass 36 GeV and 100 GeV, respectively. The former is able to produce ~7 electron/positron particles per annihilation through the fremionic channel, and the latter ~53 particles assuming a purely bosonic channel. High energy electron/positron particles up-scatter the Cosmic Microwave Background (CMB) photons into higher energies via the inverse-Compton scattering. The process produces a power-law electron/positron energy spectrum of index -1 in the energy range of interest, independent of the initial energy distribution. The corresponding energy spectrum of the up-scattered photons is a power-law of index -1/2, if absorption by the gas is not included. The scattered photons photo-heat the gas by releasing electrons which deposit a fraction (14%)of their energy as heat into the ambient medium. For uniformly distributed neutralinos the heating is insignificant. The effect is greatly enhanced by the clumping of neutralinos into dense haloes. We use a time-dependent clumping model which takes into account the damping of density fluctuations on mass scales smaller than ~10^{-6}M_sol. With this clumping model, the heating mechanism boosts the gas temperature above that of the CMB after a redshift of z 30. By z\approx 10 the gas temperature is nearly 100 times its temperature when no heating is invoked. Similar increase is obtained for the two neutralino masses considered.

paper removed: double submission

We present a simplified dynamical model of the ``Bullet'' system of two colliding clusters. The model constrains the masses of the system by requiring that the orbits of the main and sub components satisfy the cosmological initial conditions of vanishing physical separation a Hubble time ago. This is also known as the timing argument. The model considers a system embedded in an over-dense region. We argue that a relative speed of $4500 \rm km/s$ between the two components is consistent with cosmological conditions if the system is of a total mass of $2.8\times 10^{15}h^{-1} M_\odot$ is embedded in a region of a (mild) over-density of 10 times the cosmological background density. Combining this with the lensing measurements of the projected mass, the model yields a ratio of 3:1 for the mass of the main relative to that of the subcomponent. The effect of the background weakens as the relative speed between the two components is decreased. For relative speeds lower than $\sim 3700\rm km/s$, the timing argument yields masses which are too low to be consistent with lensing.

We address the question of whether or not assembly bias arises in the absence of highly non-linear effects such as tidal stripping of halos near larger mass concentrations. Therefore, we use a simplified dynamical scheme where these effects are not modeled. We choose the punctuated Zel'dovich (PZ) approximation, which prevents orbit mixing by coalescing particles coming within a critical distance of each other. A numerical implementation of this approximation is fast, allowing us to run a large number of simulations to study assembly bias. We measure an assembly bias from 60 PZ simulations, each with 512^3 cold particles in a 128h^-1 Mpc cubic box. The assembly bias estimated from the correlation functions at separations < 5h^-1 Mpc for objects (halos) at z=0 is comparable to that obtained in full N-body simulations. For masses 4x10^11 h^-1 Mo the "oldest" 10% haloes are 3-5 times more correlated than the "youngest" 10%. The bias weakens with increasing mass, also in agreement with full N-body simulations. We find that that halo ages are correlated with the dimensionality of the surrounding linear structures as measured by the parameter (\lambda_1+\lambda_2+\lambda_3)/ (\lambda_1^2+\lambda_2^2+\lambda_3^2)^{1/2} where \lambda_i are proportional to the eigenvalues of the velocity deformation tensor. Our results suggest that assembly bias may already be encoded in the early stages of the evolution.

(Abridged) Existing models invoking AGN activty to resolve the cooling flow conundrum in galaxy clusters focus exclusively on the role of the central galaxy. Such models require fine-tuning of highly uncertain microscopic transport properties to distribute the thermal thermal over the entire cluster cooling core. We propose that the ICM is instead heated by multiple, spatially distributed AGNs. There is mounting observational evidence for multiple AGNs in cluster environments. Active AGNs drive bubbles into the ICM. We identify three distinct interactions between the bubble and the ICM: (1) Upon injection, the bubbles expand rapidly in situ to reach pressure equilibrium with their surroundings, generating shocks and waves whose dissipation is the principal source of ICM heating. (2) Once inflated, the bubbles rise buoyantly at rate determined by a balance with the viscous drag force, which itself results in some additional heating. (3) Rising bubbles expand and compress their surroundings. This process is adiabatic and does not contribute to any additional heating; rather, the increased ICM density due to compression enhances cooling. Our model sidesteps the ``transport'' issue by relying on the spatially distributed galaxies to heat the cluster core. We include self regulation in our model by linking AGN activity in a galaxy to cooling characteristics of the surrounding ICM. We use a spherically symmetric one-dimensional hydrodynamical code to carry out a preliminary study illustrating the efficacy of the model. Our self-regulating scenario predicts that there should be enhanced AGN activity of galaxies inside the cooling regions compared to galaxies in the outer parts of the cluster. This prediction remains to be confirmed or refuted by observations.

We present general relativistic solutions for self-similar spherical perturbations in an expanding cosmological background of cold pressure-less gas. We focus on solutions having shock discontinuities propagating in the surrounding cold gas. The pressure, $p$, and energy-density, $\mu$, in the shock-heated matter are assumed to obey $p=w\mu$, where $w$ is a positive constant. Consistent solutions are found for shocks propagating from the symmetry center of a region of a positive density excess over the background. In these solutions, shocks exist outside the radius marking the event horizon of the black hole which would be present in a shock-less collapse. For large jumps in the energy-density at the shock, a black hole is avoided altogether and the solutions are regular at the center. The shock-heated gas does not contain any sonic points, provided the motion of the cold gas ahead of the shock deviates significantly from the Hubble flow. For shocks propagating in the uniform background, sonic points always appear for small jumps in the energy-density. We also discuss self-similar solutions without shocks in fluids with $w<-1/3$.

We present a measurement of the probability distribution function (PDF) of the transmitted flux in the Lya forest from a sample of 3492 quasars included in the SDSS DR3 data release. Our intention is to investigate the sensitivity of the Lya flux PDF as measured from low resolution and low signal-to-noise data to a number of systematic errors such as uncertainties in the mean flux, continuum and noise estimate. The quasar continuum is described by the superposition of a power law and emission lines. We perform a power law continuum fitting on a spectrum-by-spectrum basis, and obtain an average continuum slope of 0.59 +/- 0.36 in the redshift range 2.5<z<3.5. Taking into account the variation in the continuum indices increases the mean flux by 3 and 7 per cent at z=3 and 2.4, respectively, as compared to the values inferred with a single (mean) continuum slope. We compare our measurements to the PDF obtained with mock lognormal spectra, whose statistical properties have been constrained to match the observed Lya flux PDF and power spectrum of high resolution data. Using our power law continuum fitting and the SDSS pipeline noise estimate yields a poor agreement between the observed and mock PDFs. Allowing for a break in the continuum slope and, more importantly, for residual scatter in the continuum level substantially improves the agreement. A decrease of 10-15 per cent in the mean quasar continuum with a typical rms variance at the 20 per cent level can account for the data, provided that the noise excess correction is no larger than 10 per cent.

The morphology of the total gas and the neutral hydrogen (H I) distributions during the cosmological epoch of reionization can be quantified with Minkowski Functionals (MFs) of isodensity surfaces. We compute the MFs from the output of a high-resolution numerical simulation which includes explicit treatment of transfer of UV ionizing radiation. "Galaxies" identified in the simulation using semi-analytic models of galaxy formation are assumed to be the sole sources of UV photons. The MFs of the total gas distribution are well described by the analytic expressions derived for lognormal random fields. The statistical properties of the diffuse H I depend on the gas distribution and on the way ionized regions propagate in the inter-galactic medium (IGM). The deviations of the MFs of the H I distribution from those of a lognormal random field are, therefore, caused by reionization. We use the MFs to discriminate between the various stages of reionization in the simulation. We suggest a simple model of reionization which reproduces the MFs derived from this simulation. Using random realizations of lognormal density fields, we also assess the ability of MFs to distinguish between different reionization scenarios. Our results are relevant to the analysis of future cosmological twenty-one centimeter maps.

We measure the matter power spectrum from the Lyman-alpha forest of 31 quasar spectra spanning the redshift range of 1.6-3.6. The optical depth, $\tau$, for absorption of the intergalactic medium is obtained from the flux using the inversion method of Nusser & Haehnelt (1999). The optical depth is converted to density by using a simple power law relation, $\tau \propto (1+\delta)^\alpha$. The non-linear 1D power spectrum of the gas density is then inferred with a method that makes simultaneous use of the 1 and 2 point statistics of the flux and compared against theoretical models with a likelihood analysis. A Cold Dark Matter (CDM) model with standard cosmological parameters fits the data well. The power spectrum amplitude is measured to be (assuming a flat Universe), $\sigma_8 = (0.9 \pm 0.09)\times (\Omega_{\rm m}/0.3)^{-0.3}$, with $\alpha$ varying in the range of $1.56-1.8$ with redshift. Enforcing the same cosmological parameters in all four redshift bins, the likelihood analysis suggests some evolution in the density-temperature relation and the thermal smoothing length of the gas. The inferred evolution is consistent with that expected if reionization of \hep occurred at $z\sim 3.2$. A joint analysis with the WMAP results together with a prior on the Hubble constant as suggested by the HST key project data, yields values of $\Omega_{\rm m}$ and $\sigma_8$ that are consistent with the cosmological concordance model. We also perform a further inversion to obtain the linear 3D power spectrum of the matter density fluctuations.

We have modelled the process of reionization of the IGM by photoionization by galaxies, in order to predict the epoch of reionization. We use a sophisticated semi-analytic model to track the formation of galaxies. Our study represents a much more complete and physically consistent modelling of reionization than has been conducted in the past, containing significant improvements in the modelling of the collapse of baryons into dark matter halos, and in the model for the cooling and condensation of gas within halos (including photoheating from a self-consistently computed ionizing background and cooling due to H2). We find that reionization can be achieved by z~10-20 in a LCDM cosmological model with sigma8~0.9. However, a cosmological model with a running spectral index is only able to achieve reionization before z~9 if extreme assumptions are made about the physics of feedback at high redshifts. We also consider the galaxy formation model discussed by Baugh et al., and find that it is able to reionize the Universe by z~12. The previous results assume that all of the ionizing photons produced by stars in galaxies are able to escape and ionize the IGM. If this is not the case, then the redshift of reionization could be substantially reduced. We find that extended periods of partial reionization and double reionizations can occur in models in which the first stars formed via cooling by H2 molecules, are very massive, and in which the escape fraction of ionizing photons ~10-30%. Such models do not fully reionize until z~6-7, but predict an electron scattering optical depth as large as 0.15. Models with lower sigma8=0.7-0.8 as suggested by the recent WMAP three year data have reduced redshifts of reionization, but can be consistent with the lower optical depth also suggested by the WMAP three year data.

We test the ability of the numerical action method (NAM) to recover the individual orbit histories of mass tracers in an expanding universe in a region of radius 26Mpc/h, given the masses and redshift-space coordinates at the present epoch. The mass tracers are represented by dark matter haloes identified in a high resolution N-body simulation of the standard LCDM cosmology. Since previous tests of NAM at this scale have traced the underlying distribution of dark matter particles rather than extended haloes, our study offers an assessment of the accuracy of NAM in a scenario which more closely approximates the complex dynamics of actual galaxy haloes. We show that NAM can recover present-day halo distances with typical errors of less than 3 per cent, compared to 5 per cent errors assuming Hubble flow distances. The total halo mass and the linear bias were both found to be constained at the 50 per cent level. The accuracy of individual orbit reconstructions was limited by the inability of NAM, in some instances, to correctly model the positions of haloes at early times solely on the basis of the redshifts, angular positions, and masses of the haloes at the present epoch. Improvements in the quality of NAM reconstructions may be possible using the present-day three-dimensional halo velocities and distances to further constrain the dynamics. This velocity data is expected to become available for nearby galaxies in the coming generations of observations by SIM and GAIA.

An excess of sight-lines close to Lyman-break galaxies (LBGs) with little or no absorption in QSO absorption spectra has been reported and has been interpreted as the effect of galactic winds on the Intergalactic Medium. We use here numerical simulations to investigate the flux probability function close to plausible sites of LBGs. We show that the flux distribution near our LBGs in the simulation depends strongly on redshift, and is very sensitive to the averaging procedure. We show that a model without galactic winds and a model with a wind bubble size of 0.5Mpc/h (comoving) are equally consistent with the new determination of the conditional flux distribution by Adelberger et al. (2005). Models with the larger bubble sizes suggested by the previous observations of Adelberger et al. (2003) based on a much smaller sample at higher redshift are not consistent with the new data. We therefore argue that the volume filling factor of galactic winds driven by LBGs may be much smaller than previously thought, and that most of the metals responsible for the metal absorption associated with the low column density Lya forest are unlikely to have been ejected by LBGs.

Protons gain energy in short range collisions with heavier dark matter particles (DMPs) of comparable velocity dispersion. We examine the conditions under which the heating of baryons by scattering off DMPs can offset radiative cooling in the cores of galaxy clusters. Collisions with a constant cross section independent of the relative velocity of the colliding particles, cannot produce stable thermal balance. In this case, avoiding an unrealistic increase of the central temperatures yields the upper bound on the cross-section, \sigma_xp<10^-25 cm^2 (m_x/m_p), where m_x and m_p are the DMP and proton mass, respectively. A stable balance, however, can be achieved for a power law dependence on the relative velocity, V, of the form \sigma_xp \propto V^a with a<-3. An advantage of this heating mechanism is that it preserves the metal gradients observed in clusters.

We examine the possibility of constraining the cosmological mean mass and dark energy densities by an application of the Alcock-Paczynski test on redshifted 21-cm maps of the epoch of reionization. The 21-cm data will be provided as a function of frequency and angular positions on the sky. The ratio of the frequency to angular distance scales can be determined by inspecting the anisotropy pattern of the correlation function of brightness temperature. We assess the sensitivity of the distance ratio to the cosmological parameters and present a technique for disentangling geometrical distortions from redshift space distortions caused by peculiar motions.

We present a Monte-Carlo model of He II reionization by QSOs and its effect on the thermal state of the clumpy intergalactic medium (IGM). The model assumes that patchy reionization develops as a result of the discrete distribution of QSOs. It includes various recipes for the propagation of the ionizing photons, and treats photo-heating self-consistently. The model provides the fraction of He III, the mean temperature in the IGM, and the He II mean optical depth -- all as a function of redshift. It also predicts the evolution of the local temperature versus density relation during reionization. Our findings are as follows: The fraction of He III increases gradually until it becomes close to unity at $z\sim 2.8-3.0$. The He II mean optical depth decreases from $\tau\sim 10$ at $z\geq 3.5$ to $\tau\leq 0.5$ at $z\leq 2.5$. The mean temperature rises gradually between $z\sim 4$ and $z\sim 3$ and declines slowly at lower redshifts. The model predicts a flattening of the temperature-density relation with significant increase in the scatter during reionization at $z\sim 3$. Towards the end of reionization the scatter is reduced and a tight relation is re-established. This scatter should be incorporated in the analysis of the Ly$\alpha$ forest at $z\leq 3$. Comparison with observational results of the optical depth and the mean temperature at moderate redshifts constrains several key physical parameters.

We presents results of joint modeling of the probability distribution function (PDF) and the one-dimensional power spectrum (PS) of the Lya forest flux decrement. The sensitivity of these statistical measures to the shape and amplitude of the linear matter power spectrum is investigated using two variants of the LCDM cosmology. In the first model, the linear power spectrum has a scale-invariant spectral index, whereas in the second, it has a negative running index (RSI). We generate mock catalogs of QSO spectra, and compare their statistical properties to those of the observations at z=3. We perform a joint fit of the power spectrum and PDF. A scale-invariant model with \sigma_8=0.9 matches well the data if the mean IGM temperature is T\leq 15000K. For higher temperature, it tends to overestimate the flux power spectrum over scales k < 0.01s/km. The discrepancy is less severe when the PS alone is fitted. However, models matching the PS alone do not yield a good fit to the PDF. A joint analysis of the flux PS and PDF tightens the constraints on the model parameters and reduces systematic biases. The RSI model is consistent with the observed PS and PDF only if the temperature is T\geq 20000K. The best fit models reproduce the slope and normalisation of the column density distribution, irrespective of the shape and amplitude of the linear power spectrum. They are also consistent with the observed line-width distribution given the large uncertainties. Our joint analysis suggests that \sigma_8 is likely to be in the range 0.7 - 0.9 for a temperature 10000K<T<20000K and a reasonable reionization history.

We discuss the stability of galactic disks in which the energy of interstellar clouds is gained in encounters with expanding supernova remnants and lost in inelastic collisions. Energy gain and loss processes introduce a phase difference between the pressure and density perturbations, making disks unstable on small scales for several recipes of star formation. This is in contrast to the standard stability analysis in which small scale perturbations are stabilized by pressure. In the limit of small scales the dispersion relation for the growth rate reduces to that of thermal instabilities in a fluid without gravity. If instabilities lead to star formation, then our results imply a secondary mode of star formation which operates on small scales and feeds on the existence of a primary mode on intermediate scales. This may be interpreted as a positive feedback. Further, the standard stability criterion on intermediate scales is significantly modified.

We use a one dimensional hydrodynamical code to study the evolution of spherically symmetric perturbations in the framework of Modified Newtonian Dynamics (MOND). The code evolves spherical gaseous shells in an expanding Universe by employing a MOND-type relationship between the fluctuations in the density field and the gravitational force, $g$. We focus on the evolution of initial density perturbations of the form $\delta_{i}\sim r_{i}^{-s}$ for $0<s<3$. A shell is initially cold and remains so until it encounters the shock formed by the earlier collapse of shells nearer to the centre. During the early epochs $g$ is sufficiently large and shells move according to Newtonian gravity. As the physical size of the perturbation increases with time, $g$ gets smaller and the evolution eventually becomes MOND dominated. However, the density in the inner collapsed regions is large enough that they re-enter the Newtonian regime. The evolved gas temperature and density profiles tend to a universal form that is independent of the the slope, $s$, and of the initial amplitude. An analytic explanation of this intriguing result is offered. Over a wide range of scales, the temperature, density and entropy profiles in the simulations, depend on radius roughly like $r^{0.5}$, $r^{-1.5}$ and $r^{1.5}$, respectively. We compare our results with XMM-{\it Newton} and {\it Chandra} observations of clusters. The temperature profiles of 16 observed clusters are either flat or show a mild decrease at $R\gtsim 200\rm kpc$. MOND profiles show a significant increase that cannot reconciled with the data. Our simulated MOND clusters are substantially denser than the observed clusters. It remains to be seen whether these difficulties persist in three-dimensional hydrodynamical simulations with generic initial conditions.

We re-examine the role of collisions in decoupling the HI 21-cm spin temperature from the cosmic microwave background (CMB). The cross section for de-exciting the 21-cm trasition in collisions with free electrons is more than 10 times larger than it is in collisions with other atoms. If the fraction of free electrons in the diffuse cosmic gas is between 10 and 30 per cent then collisions alone can decouple the spin temperature from the cosmic microwave background (CMB), even in moderately under-dense regions at z >15. This decoupling is especially important during the very early stages of re-ionization when a Ly-alpha continuum background had yet to be established. As a detailed example, we develop a semi-analytic model to quantify 21-cm emission signatures from a diffuse gas which is partially ionized at z>10 by an X-ray background. We find 21-cm differential brightness temperature fluctuations with a mean of ~8 mK and a rms value as large as 5 mK, for a frequency resolution bandwidth of 100 KHz and a beamsize of 3 arcmin. Another example where free electron-atom collisions are important is during the recombination of bubbles ionized by short-lived UV sources. When the ionized fraction in these bubbles drops to 10-20 per cent their differential temperature can be as high as 10 mK. We also consider the correlation of the 21-cm temperature as a function of frequency. By isotropy, a comparison with the angular projected correlation in the plane of the sky should yield an estimate for the matter density parameter, Omega_m, and the dark energy component. This variant of the Alcock-Paczyski test is applicable for any reionization scenario.

Numerical simulations show that a long-range scalar interaction in a single species of massive dark matter particles causes voids between the concentrations of large galaxies to be more nearly empty, suppresses accretion of intergalactic matter onto galaxies at low redshift, and produces an early generation of dense dark matter halos. These three effects, in moderation, seem to be improvements over the Lambda CDM model predictions for cosmic structure formation. Because the scalar interaction in this model has negligible effect on laboratory physics and the classical cosmological tests, it offers an observationally attractive example of cosmology with complicated physics in the dark sector, notably a large violation of the weak equivalence principle.

It is shown that a conversion of neurons into strange quark matter in neutron stars is possible by means of a combustion process with a well-defined front. Conditions for the realization of a specific combustion mode, whether deflagration, detonation, or fast combustion, are discussed for several forms of the equations of state of neutron and strange matters.

We describe the feasibility of detecting the gravitational deflection of light emitted by stars moving under the influence of the massive object at the Galactic center. Light emitted by a star orbiting behind the central mass has a smaller impact parameter than the star itself, and suffers the effect of gravitational lensing, providing a closer probe of the central mass distribution and hence a stricter test of the black hole hypothesis. A mass of $4.3\times 10^{6} M_{\odot}$ causes a $0.1-2\rm mas$ deviation in the apparent position of orbiting stars projected within $10^{\circ}$ of the line of sight to the galactic center. In addtion, we may uniquely constrain the distance to the center of the galaxy because lensing deflections constrain the ratio $\rg/R_{0}$ of the Schwarzschild radius to the distance to the black hole, $R_{o}$, whereas the ratio $\rg/R_{o}^{3}$ is obtained by fitting the orbit.

We present a model for one-dimensional (1D) matter power spectra in redshift space as estimated from data provided along individual lines of sight. We derive analytic expressions for these power spectra in the linear and nonlinear regimes, focusing on redshift distortions arising from peculiar velocities. In the linear regime, redshift distortions enhance the 1D power spectra only on small scales, and do not affect the power on large scales. This is in contrast to the effect of distortions on three-dimensional (3D) power spectra estimated from data in 3D space, where the enhancement is independent of scale. For CDM cosmologies, the 1D power spectra in redshift and real space are similar for wavenumbers $q<0.1h/Mpc$ where both have a spectral index close to unity, independent of the details of the 3D power spectrum. Nonlinear corrections drive the 1D power spectrum in redshift space into a nearly universal shape over scale $q<10h/Mpc$, and suppress the power on small scales as a result of the strong velocity shear and random motions. The redshift space, 1D power spectrum is mostly sensitive to the amplitude of the initial density perturbations. Our results are useful in particular for power spectra computed from the SDSS quasars sample.

We study the conversion of neutron matter into strange matter as a detonation wave. The detonation is assumed to originate from a central region in a spherically symmetric background of neutrons with a varying radial density distribution. We present self-similar solutions for the propagation of detonation in static and collapsing backgrounds of neutron matter. The solutions are obtained in the framework of general relativistic hydrodynamics, and are relevant for the possible transition of neutron into strange stars. Conditions for the formation of either bare or crusted strange stars are discussed.

We study the effect of galactic outflows on the Lya forest. The winds are modelled as fully ionised spherical bubbles centered around galactic haloes. The observed flux probability distribution and flux power spectrum limit the volume filling factor of bubbles to be less than 10%. We have compared the mean flux as a function of distance from haloes with the Adelberger etal (ASSP) measurement. A model with bubbles of constant radius of 1.5Mpc/h surrounding the most massive haloes matches the high transmissivity at separations <0.5Mpc/h, but the increase of the transmissivity at small scales is more gradual than observed. At separations <0.5Mpc/h, our estimate of the cosmic variance error is $\Delta F\sim 0.3$, 30% higher than that of ASSP. The difficulty in matching the rise in the transmissivity at small separations is caused by residual absorption of neutral hydrogen lying physically outside the bubbles but having a redshift position similar to the haloes. The flux level is thus sensitive to the amplitude of the coherent velocity shear near haloes which increases with halo mass in the simulation. A model where LBGs are starbursts in small mass haloes matches the observations with smaller bubble size than one where more massive haloes host the LBGs. A starburst model with a bubble radius of 1Mpc/h and a volume filling factor of 2% is consistent with the ASSP measurements at the 1-1.5$\sigma$ level. If this model is correct the sharp rise of the transmissivity at small separation in the ASSP sample is due to cosmic variance and is expected to become more moderate for a larger sample.

We revisit the process of gravitational sedimentation of helium and heavy elements in the intra-cluster medium. We find that helium applies an inward drag force on heavy elements, boosting their sedimentation speed to nearly half its own. This speed is almost independent of the mass and the electric charge of heavy elements. In the absence of small-scale magnetic fields, helium sedimentation can increase the He/H abundance ratio in the cores of hot clusters by three orders of magnitude. It also steepens the baryonic density profile yielding a higher X-ray luminosity, which offers an explanation of the observed luminosity-temperature relation. If the primordial He/H ratio is assumed, then the gas density inferred from the observed X-ray emissivity might be underestimated by 30% in the cores of clusters and overestimated by 7% in the outer regions. The dark matter density on the other hand might be overestimated by a factor of 8/3 in the cores and underestimated by 18% in outer regions.

We derive similarity solutions for the expansion of negative initial density perturbations \delta M/M \propto r^-s (s>0) in an Einstein-de Sitter universe filled with collisional baryonic gas and collisionless matter. The solutions are obtained for planar, cylindrical, and spherical perturbations. For steep perturbations, central cavities surrounded by over-dense shells expanding as t^{2(s+1)/3s} develop in both matter components. We also consider the case when baryonic shells are driven by internal pressure. Before redshift \sim 8-10 baryonic shells cool by inverse Compton scattering with the cosmic background radiation and can form bound objects of mass \sim 10^11\Omega_b^{-1/2} M_{\odot} (1+z)^{-3/2}. Hot shocks of radius R cause a Compton y-distortion of order y=10^{-5}\Omega_b (R/10 Mpc)^2 f_F, where f_F is the volume filling factor of these shells.

We present predictions for two statistical measures of the hydrogen reionization process at high redshift. The first statistic is the number of neutral segments identified in spectra of high redshift QSOs as a function of their length. The second is the cross-correlation of neutral regions with possible sources of ionizing radiation. These independent probes are sensitive to the topology of the ionized regions. If reionization proceeded from high to low density regions then the cross-correlation will be negative, while if voids were ionized first then we expect a positive correlation and a relatively small number of long neutral segments. We test the sensitivity of these statistics for reionization by stars in high redshift galaxies. The flux of ionizing radiation emitted from stars is estimated by identifying galaxies in an N-body simulation using a semi-analytic galaxy formation model. The spatial distribution of ionized gas is traced in various models for the propagation of the ionization fronts. A model with ionization proceeding from high to low density regions is consistent with the observations of Becker et al. (2001), while models in which ionization begins in the lowest density regions appear to be inconsistent with the present data.

We use a simple analytic model to compute the angular correlation function of clusters identified in upcoming thermal SZ effect surveys. We then compute the expected fraction of close pairs of clusters on the sky that are also close along the line of sight. We show how the expected number of cluster pairs as a function of redshift is sensitive to the assumed biasing relation between the cluster and the mass distribution. We find that, in a LambdaCDM model, the fraction of physically associated pairs is 70% for angular separations smaller than 20 arcmin and clusters with specific flux difference larger than 200 mJy at 143 GHz. The agreement of our analytic results with the Hubble volume N-body simulations is satisfactory. These results quantify the feasibility of using SZ surveys to compile catalogues of superclusters at any redshifts.

We present extensive tests of the Fast Action Method (FAM) for recovering the past orbits of mass tracers in an expanding universe from their redshift-space coordinates at the present epoch. The tests focus on the reconstruction of present-day peculiar velocities using mock catalogs extracted from high resolution $N$-body simulations. The method allows for a self-consistent treatment of redshift-space distortions by direct minimization of a modified action for a cosmological gravitating system. When applied to ideal, volume limited catalogs, FAM recovers unbiased peculiar velocities with a 1-D, 1\sigma error of ~220 km/s, if velocities are smoothed on a scale of 5 Mpc/h. Alternatively, when no smoothing is applied, FAM predicts nearly unbiased velocities for objects residing outside the highest density regions. In this second case the 1\sigma$error decreases to a level of ~150 km/s. The correlation properties of the peculiar velocity fields are also correctly recovered on scales larger than 5 Mpc/h. Similar results are obtained when FAM is applied to flux limited catalogs mimicking the IRAS PSCz survey. In this case FAM reconstructs peculiar velocities with similar intrinsic random errors, while velocity-velocity correlation properties are well reproduced beyond scales of ~8 Mpc/h. We also show that FAM provides better velocity predictions than other, competing methods based on linear theory or Zel'dovich approximation. These results indicate that FAM can be successfully applied to presently available galaxy redshift surveys such as IRAS PSCz.

We examine the implications of Modified Newtonian Dynamics (MOND) on the large scale structure in a Friedmann-Robertson-Walker universe. We employ a ``Jeans swindle'' to write a MOND-type relationship between the fluctuations in the density and the gravitational force, $\vg$. In linear Newtonian theory, $|\vg|$ decreases with time and eventually becomes $<g_0$, the threshold below which MOND is dominant. If the Newtonian initial density field has a power-law power-spectrum of index $n<-1$, then MOND domination proceeds from small to large scale. At early times MOND tends to drive the density power-spectrum towards $k^{-1}$, independent of its shape in the Newtonian regime. We use N-body simulations to solve the MOND equations of motion starting from initial conditions with a CDM power-spectrum. MOND with the standard value $g_0=10^{-8} cm s^{-2}$, yields a high clustering amplitude that can match the observed galaxy distribution only with strong (anti-) biasing. A value of $g_0 \approx 10^{-9}cm s^{-2}$, however, gives results similar to Newtonian dynamics and can be consistent with the observed large scale structure.

We study the signal in the Cosmic Microwave Background (CMB) polarization anisotropy resulting from patchy reionization. It is well known that the primordial polarization of the CMB is very sensitive to the details of reionization. Combining a semi-analytic model of galaxy formation, in which the optical depth to the reionization epoch is in the range 0.014 to 0.048, with a high resolution N-body simulation we find that reionization generates a peak with amplitude 0.05~0.15 \mu K at large angular scales. The position of this peak reveals the size of the horizon at reionization, whilst its amplitude is a measure of the optical depth to reionization. On small scales, ionized patches prior to full reionization create a second order polarization signal due to the coupling of the free electron density fluctuation with the quadrupole moment of the temperature anisotropy. Careful study reveals that the coupling generates the same power spectra for electric and magnetic modes, whose amplitudes of polarization anisotropies from this process are predicted to be ~10nK. The amplitude depends strongly on the total baryon density $\Omega_b$ and on the spatial correlations of the free electron density in the ionized regions, and weakly on the fraction of ionizing photons able to escape their source galaxy, $f_{esc}$. The first- and second-order signals are therefore sensitive to the details of how the reionization occurred. Detection of these signals will place important constraints on the reionization history of the Universe.

We present predictions for the reionization of the intergalactic medium by stars in high-redshift galaxies. We calculate ionizing luminosities of galaxies, including the effects of absorption by interstellar gas and dust on the escape fraction, in order to calculate the filling factor of ionized hydrogen in the IGM. For a $\Lambda{\rm CDM}$ cosmology and a physical calculation of the escape fraction, we find that the hydrogen in the IGM will be reionized at redshift z=6.1 if the IGM has uniform density, but only by z=4.5 if the IGM is clumped. We combine our model with an N-body simulation of the distribution of dark matter in the universe in order to calculate the secondary temperature anisotropies induced in the cosmic microwave background (CMB) by scattering off free electrons. The models predict a spectrum of secondary anisotropies covering a broad range of angular scales, with fractional temperature fluctuations $\sim 10^{-7}-10^{-6}$ on arcminute scales. The amplitude depends somewhat on the geometry of reionization, with models in which the regions of highest gas density are reionized first giving larger CMB fluctuations than the case where galaxies ionize surrounding spherical regions, and models where low density regions reionize first giving the smallest fluctuations. Measurement of these anisotropies can therefore put important constraints on the reionization process, in particular, the redshift evolution of the filling factor, and should be a primary objective of a next generation submillimeter telescope such as the Atacama Large Millimeter Array. (abrdiged)

We derive the asymptotic mass profile near the collapse center of an initial spherical density perturbation, $\delta \propto M^{-\epsilon}$, of collision-less particles with non-radial motions. We show that angular momenta introduced at the initial time do not affect the mass profile. Alternatively, we consider a scheme in which a particle moves on a radial orbit until it reaches its turnaround radius, r_*. At turnaround the particle acquires an angular momentum $L={\cal L} \sqrt{GM_* r_*}$ per unit mass, where M_* is the mass interior to r_*. In this scheme, the mass profile is $M\propto r^{3/(1+3\epsilon)}$ for all $\epsilon >0$, in the region $r/r_t\ll {\cal L}$, where r_t is the current turnaround radius. If ${\cal L} \ll 1$ then the profile in the region ${\cal L} \ll r/r_t \ll$ is $M\propto r$ for $\epsilon <2/3$. The derivation relies on a general property of non-radial orbits which is that ratio of the pericenter to apocenter is constant in a force field k(t) r^{n} with k(t) varying adiabatically.

Similarity solutions are found for the adiabatic collapse of density perturbations $\delta M/M \propto r^{-s}$ $(s>0)$ in a flat universe containing collisional gas only. The solutions are obtained for planar, cylindrical, and spherical perturbations with zero initial pressure. For adiabatic index $\gamma\ge 4/3$, a shock develops at a fixed fraction of the current turnaround distance. Near the center of a spherical perturbations with $\gamma>4/3$ and $s>1/2$, the gas is in quasi-hydrostatic equilibrium (pressure supported) and has an asymptotic power law density profile, $\rho\sim r^{-3s/(s+1)}$, independent of $\gamma$. For $s\le 1/2$, the profile depends on $\gamma$, the pressure is finite, the temperature decreases inward, and gravity dominates pressure causing a continuous inward flow. Although for $1/2<s<2$ the temperature decreases at the center, the gas is pressure supported. The pressure is finite in cylindrical perturbations for $s\le 2(\gamma-1)/(3\gamma-4)$, and in planar perturbations for any $s>0$. We also derive the asymptotic behaviour of the gas variables near the center in a universe dominated by collisionless matter. In such a universe, the gas in a spherical perturbation with $s<2$ cannot be pressure supported and the temperature approaches a constant near the center. The solutions and the asymptotic behaviour are relevant for modelling the gas distribution in galaxy clusters and pancake-like superclusters, and determining the structure of haloes of self-interacting dark matter with large interaction cross section.

We present a comparison between the peculiar velocity field measured from the ENEAR all-sky $D_n-\sigma$ catalog and that derived from the galaxy distribution of the IRAS PSCz redshift survey. The analysis is based on a modal expansion of these data in redshift space by means of spherical harmonics and Bessel functions. The effective smoothing scale of the expansion is almost linear with redshift reaching 1500km/s at 3000km/s. The general flow patterns in the filtered ENEAR and PSCz velocity fields agree well within 6000km/s, assuming a linear biasing relation between the mass and the PSCz galaxies. The comparison allows us to determine the parameter $\beta=\Omega^{0.6}/b$, where $\Omega$ is the cosmological density parameter and $b$ is the linear biasing factor. A likelihood analysis of the ENEAR and PSCz modes yields $\beta=0.5 +- 0.1$, in good agreement with values obtained from Tully-Fisher surveys.

We identify high density regions of supercluster size in high resolution N-body simulations of a representative volume of three Cold Dark Matter Universes. By assuming that (1) the density and peculiar velocities of baryons trace those of the dark matter, and (2) the temperature of plasma is proportional to the velocity dispersion of the dark matter particles in regions where the crossing times is smaller than the supercluster free-fall time, we investigate how thermal motions of electrons in the intra-cluster medium and peculiar velocity of clusters can affect the secondary anisotropies in the cosmic microwave background (CMB). We show that the thermal effect dominates the kinematic effect and that the largest thermal decrements are associated with the most massive clusters in superclusters. Thus, searching for the presence of two or more close large CMB decrements represents a viable strategy for identifying superclusters at cosmological distances. Moreover, maps of the kinematic effect in superclusters are characterized by neighboring large peaks of opposite signs. These peaks can be as high as ~ 10 microK at the arcminute angular resolution. Simultaneous pointed observations of superclusters in the millimeter and submillimeter bands with upcoming sensitive CMB experiments can separate between the thermal and kinematic effect contributions and constrain the evolution of the velocity field in large overdense regions.

Analytic solutions for the evolution of cosmological linear density perturbations in the baryonic gas and collisionless dark matter are derived. The solutions are expressed in a closed form in terms of elementary functions, for arbitrary baryonic mass fraction. They are obtained assuming $\Omega=1$ and a time independent comoving Jeans wavenumber, $k_J$. By working with a time variable $\tau\equiv \ln(t^{2/3})$, the evolution of the perturbations is described by linear differential equations with constant coefficients. The new equations are then solved by means of Laplace transformation assuming that the gas and dark matter trace the same density field before a sudden heating epoch. In a dark matter dominated Universe, the ratio of baryonic to dark matter density perturbation decays with time roughly like $\exp(-5\tau/4)\propto t^{-5/6}$ to the limiting value $1/[1+(k/k_J)^2]$. For wavenumbers $k>k_J/\sqrt{24}$, the decay is accompanied with oscillations of a period $ 8\pi/\sqrt{24 (k/k_J)^2 -1}$ in $\tau$. In comparison, as $\tau $ increases in a baryonic matter dominated Universe, the ratio approaches $1-(k/k_J)^2$ for $k\le k_J$, and zero otherwise.

The Ly-alpha forest provides important constraints on the smoothness of the universe on large scales. We calculate the flux distribution along line-of-sight to quasars in a universe made of randomly distributed clumps, each of them with a Rayleigh-L`evi fractal structure. We consider the probability distribution function of the normalised flux in the line-of-sight to quasars. We show that the truncated clustering hierarchy model shows far too many voids along the line-of-sight to quasars compared with the observed flux distribution and the distribution in a Cold Dark Matter model. This supports the common view that on large scales the universe is homogeneous, rather than fractal-like.

Given an irrotational (vorticity free) velocity field in real space, we prove that, in the distant observer limit and in the absence of multi-valued zones, the associated velocity field in redshift space is also irrotational. The proof does not rely on any approximation to gravitational dynamics. The result can be particularly useful for the analysis of redshift distortions and for reconstruction methods of cosmological velocity fields from galaxy redshift surveys, in the nonlinear regime. Although the proof is restricted to the distant observer limit, we show that the POTENT method can be modified to derive the full real space velocity field as a function of the redshift space coordinate, thus avoiding spatial Malmquist biases.

Given the present distribution of mass tracing objects in an expanding universe, we develop and test a fast method for recovering their past orbits using the least action principle. In this method, termed FAM for Fast Action Minimization, the orbits are expanded in a set of orthogonal time-base functions satisfying the appropriate boundary conditions at the initial and final times. The conjugate gradient method is applied to locate the extremum of the action in the space of the expansion coefficients of the orbits. The TREECODE gravity solver routine is used for computing the gravitational fields appearing in the action and its gradients. The time integration of the Lagrangian is done using Gaussian quadratures. FAM allows us to increase the number of galaxies used in previous numerical action principle implementations by more than one order of magnitude. For example, orbits for the $\sim 15,000$ IRAS PSC$z$ galaxies can be recovered in $\sim 12,000$ CPU seconds on a 400MHz DEC-Alpha machine. FAM can recover the present peculiar velocities of particles and the initial fluctuations field. It successfully recovers the flow field down to clusters scales, where deviations of the flow from the Zel'dovich solution are significant. We also show how to recover orbits from the present distribution of objects as in redshift space by direct minimization of a modified action, without iterating the solution between real and redshift spaces.

The real-space optical depth distribution along the line of sight to the QSO Q1422+231 is recovered from two HIRES spectra. The first two moments of the truncated optical depth distribution are used to constrain the density fluctuation amplitude of the intergalactic medium (IGM). The $rms$ of the IGM density at $z\approx 3.25$ estimated from the first spectrum is $\sigma = \sqrt{\exp{[(1.8 \pm 0.27)^2/\alpha^2]}-1}$, with $1.56 <\alpha <2$ for plausible reionization histories. This corresponds to $0.9 \la \sigma \la 2.1$ with $\sigma(\alpha =1.7)= 1.44\pm 0.3 $. The values obtained from the second spectrum are higher by $\approx 20 %$. If the IGM density traces the dark matter (DM) as suggested by numerical simulations we have measured the fluctuation amplitude of the DM density at an effective Jeans scale of about a hundred to two hundred (comoving) kpc. For CDM-like power spectra the amplitude of dark matter fluctuations on these small scales depends on the cosmological density parameter $\Omega$. For power spectra normalized to reproduce the space density of present-day clusters and with a slope parameter of $\Gamma=0.21$ consistent with the observed galaxy power spectrum, the inferred $\Omega$ can be expressed as: $\Omega= 0.61 (\alpha/1.7)^{1.3}(x_{_{\rm J}} /0.62)^{-0.6}$ for a flat universe, and $\Omega= 0.91(\alpha/1.7)^{1.3} (x_{_{\rm J}}/0.62)^{-0.7}$ for a $\lambda=0$ universe. $x_{_{\rm J}}$ is the effective Jeans scale in (comoving) $\mpc$. Based on a suit of detailed mock spectra the 1-$\sigma$ error is $\approx 25 %$. The estimates increase with increasing $\Gamma$. For the second spectrum we obtain 15% lower values.

Using third-order perturbation theory, we derive a relation between the mean divergence of the peculiar velocity given density and the density itself. Our calculations assume Gaussian initial conditions and are valid for Gaussian filtering of the evolved density and velocity fields. The mean velocity divergence turns out to be a third-order polynomial in the density contrast. We test the power spectrum dependence of the coefficients of the polynomial for scale-free and standard CDM spectra and find it rather weak. Over scales larger than about 5 megaparsecs, the scatter in the relation is small compared to that introduced by random errors in the observed density and velocity fields. The relation can be useful for recovering the peculiar velocity from the associated density field, and also for non-linear analyses of the anisotropies of structure in redshift surveys.

A method for the recovery of the real space line-of-sight mass density field from Lyman absorption in QSO spectra is presented. The method makes use of a Lucy-type algorithm for the recovery of the HI density. The matter density is inferred from the HI density assuming that the absorption is due to a photoionized intergalactic medium which traces the mass distribution as suggested by recent numerical simulations. Redshift distortions are corrected iteratively from a simultaneous estimate of the peculiar velocity. The method is tested with mock spectra obtained from N-body simulations. The density field is recovered reasonably well up to densities where the absorption features become strongly saturated. The method is an excellent tool to study the density probability distribution and clustering properties of the mass density in the (mildly) non-linear regime. Combined with redshift surveys along QSO sightlines the method will make it possible to relate the clustering of high-redshift galaxies to the clustering of the underlying mass density. We further show that accurate estimates for \Omega_{bar}h^2)^2 J^{-1} H(z)^{-1} and higher order moments of the density probability function can be obtained despite the missing high density tail of the density distribution if a parametric form for the probability distribution of the mass density is assumed.

We develop a model for the growth of dark matter halos and use it to study their evolved density profiles. In this model, halos are spherical and form by quiescent accretion of matter in clumps, called satellites. The halo mass as a function of redshift is given by the mass of the most massive progenitor, and is determined from Monte-Carlo realizations of the merger-history tree. Inside the halo, satellites move under the action of the gravitational force of the halo and a dynamical friction drag force. The associated equation of motion is solved numerically. The energy lost to dynamical friction is transferred to the halo in the form of kinetic energy. As they sink into the halo, satellites continually lose matter as a result of tidal stripping. The stripped matter moves inside the halo free of dynamical friction. The evolved density profiles are steeper than those obtained by assuming that, once they have been accreted onto the parent halo, satellites remain at a fixed distance from the halo center. We find that the final density profile depends mainly on the rate of infall of matter onto the halo. This, in turn, depends on the initial fluctuation field as well as on cosmology. For mass scales where the effective spectral index of the initial density field is less than -1, the model predicts a profile which can only approximately be matched by the one parameter family of curves suggested by Navarro, Frenk and White (1997). For scale-free power-spectra with initial slope $n$, the density profile within about 1% of the virial radius is $\rho\propto r^{-\beta}$, with $3(3+n)/(5+n)\le\beta\le 3(3+n)/(4+n)$.

We present a comparison between the peculiar velocity fields measured from the SFI all-sky Sbc-Sc Tully-Fisher catalog and that derived from the \iras 1.2 Jy redshift survey galaxy distribution. The analysis is based on the expansion of these data in redshift space using smooth orthonormal functions and is performed using low and high resolution expansions, with an effective smoothing scale which increases almost linearly with redshift. The effective smoothing scales at 3000 \kms are 1500\kms and 1000\kms for the low and high resolution filters. The agreement between the high and low resolution SFI velocity maps is excellent. The general features in the filtered SFI and \iras velocity fields agree remarkably well within 6000\kms. This good agreement between the fields allows us to determine the parameter $\beta=\Omega^{0.6}/b$, where $\Omega $ is the cosmological density parameter and $b$ is the linear biasing factor. From a likelihood analysis on the SFI and \iras modes we find that $\beta=0.6\pm 0.1$ independently of the resolution of the modal expansion. For this value of $\beta$, the residual fields for the two filters show no systematic variations within 6000 \kms. Most remarkable is the lack of any coherent, redshift dependent dipole flow in the residual field.

We show that the equations of motion governing the evolution of a collisionless gravitating system of particles in an expanding universe can be cast in a form which is almost independent of the cosmological density parameter, $\Omega$, and the cosmological constant, $\Lambda$. The new equations are expressed in terms of a time variable $\tau\equiv \ln D$, where $D$ is the linear rate of growth of density fluctuations. The weak dependence on the density parameter is proportional to $\epsilon=\Omega^{-0.2}-1$ times the difference between the peculiar velocity (with respect to $\tau$) of particles and the gravity field. In the general case, the effect of this weak $\Omega$ dependence is to enhance the rate of evolution of density perturbations in dense regions. In a flat universe with $\Lambda\ne 0$, this enhancement is less pronounced than in an open universe with $\Lambda=0$ and the same $\Omega$. Using the spherical collapse model, we find that the increase of the $rms$ density fluctuations in a low $\Omega$ universe relative to that in a flat universe with the same linear normalization is $\sim 0.01 \epsilon(\Omega) < \delta^3 >$, where $\delta$ is the density field in the flat universe. The equations predict that the smooth average velocity field scales like $\Omega^{0.6}$ while the local velocity dispersion (rms value) scales, approximately, like $\Omega^{0.5}$. High resolution N-body simulations confirm these results and show that density fields, when smoothed on scales slightly larger than clusters, are insensitive to the cosmological model. Halos in an open model simulation are more concentrated than halos of the same $M/\Omega$ in a flat model simulation.

We consider a measure of the peculiar velocity field derived from the Mark III compilation of 2900 spiral galaxies (Willick \etal 1996b), using an analysis that is substantially free of bias (Nusser and Davis 1995). We expand the velocity field in a set of orthogonal, smooth modes, reducing the data to a set of 56 coefficients fitted to a maximum redshift of 6000 km/s, and maximum spherical harmonic of l=3. Equivalent mode coefficients can be computed for the gravity field derived from any whole-sky redshift catalog of galaxies, such as the \iras 1.2 Jy survey (Fisher \etal 1995). Detailed intercomparison shows the two independent fields to be remarkably aligned in general. There are, however, systematic discrepancies in the fields that lead to coherence in the residuals between them. These residuals take the form of a dipole field in the LG frame that grows with distance; it is not consistent with a bulk flow residual. A chi-squared analysis of the mode-mode comparison designed to determine which value of $\beta\equiv \Omega^{0.6}/b$ for the \iras gravity field best fits the Mark III velocity field shows instead that no value of $\beta$ is acceptable. This is in contrast to the simulated catalogs, and suggests that there must be some systematic discrepancy between the two fields.

We outline a simple approach to understanding the physical origin of bias in the distribution of galaxies relative to that of dark matter. The first step is to specify how collapsed, virialized halos of dark matter trace the overall matter distribution. The next step is to make a connection between halos and the luminous galaxies we observe. We appeal to the results of semi-analytic models of galaxy formation that are tuned to fit the observed luminosity functions of local groups and clusters. We have also used a high-resolution N-body simulation of a cold dark matter (CDM) universe to study the bias relation in more detail. The differences between the galaxy and dark matter distributions are quantified using a number of different clustering statistics. We arrive at the following general conclusions: 1) A comparison of the galaxy and dark matter density fields shows that linear biasing is a good description on large scales. 2) The bias factor b depends on the shape and normalization of the power spectrum. The lower the normalization, the larger the bias. More bias is obtained for spectra with more power on large scales. For "realistic" models, b ranges from 1 to 2.5. 3) Galaxies of different luminosity or morphology have different bias factors. 4) The scale dependence of the bias factor is weak.

We examine the anisotropies in the power spectrum by the mapping of real to redshift space. Using the Zel'dovich approximation, we obtain an analytic expression for the nonlinear redshift space power spectrum in the distant observer limit. For a given unbiased galaxy distribution in redshift space, the anisotropies in the power spectrum depend on the parameter $f(\Omega)\approx \Omega^{0.6}$, where $\Omega$ is the density parameter. We quantify these anisotropies by the ratio, $R$, of the quadrupole to monopole angular moments of the power spectrum. In contrast to linear theory, the Zel'dovich approximation predicts a decline in $R$ with decreasing scale. This departure from linear theory is due to nonlinear dynamics and not a result of incoherent random velocities. The rate of decline depends strongly on $\Omega$ and the initial power spectrum. However, we find a {\it universal} relation between the quantity $R/R_{lin}$ (where $R_{lin}$ the linear theory value of $R$) and the dimensionless variable $k/k_{nl}$, where $k_{nl}$ is a wavenumber determined by the scale of nonlinear structures. The universal relation is in good agreement with a large N-body simulation. This universal relation greatly extends the scales over which redshift distortions can be used as a probe of $\Omega$. A preliminary application to the 1.2 Jy IRAS yields $\Omega\sim 0.4$ if IRAS galaxies are unbiased.

Recent progress in the measurement of relative distances to galaxies has been quite substantial, and catalogs of 3000 galaxies with distances are soon to become available. The peculiar {\it velocity} field (deviations from Hubble flow) derivable from these catalogs, when compared to the peculiar {\it gravity} field derived from all sky redshift surveys of galaxies such as the 1.2Jy IRAS survey, leads to a unique and extremely powerful test of the density parameter $\beta \equiv \Omega^{0.6}/b_I$, where $b_I$ is the possible linear bias of the IRAS selected galaxies relative to the mass fluctuations. We review the status of these large scale flow measurements and present a new methodology to describe the two fields by means of an expansion in a set of orthogonalized functions describing a general potential flow to any chosen resolution. The parameters of the flow can be estimated by minimization of the $\chi^2$ describing the scatter of observed versus predicted linewidths from an inverse Tully-Fisher relation. By this method one can intercompare the gravity and velocity fields coefficient by coefficient, deriving a precise fit for the density parameter and an assessment of the degree of coherence between the fields. The present situation is transitory-- different analyses of the same data are not yielding consistent results. Until this embarassment is untangled, estimates of $\beta$ should be taken with a large grain of salt.

We present a method for deriving a smoothed estimate of the peculiar velocity field of a set of galaxies with measured circular velocities $\eta\equiv {\rm log} \Delta v$ and apparent magnitudes $m$. The method is based on minimizing the scatter of a linear inverse Tully-Fisher relation $\eta= \eta(M)$ where the absolute magnitude of each galaxy is inferred from its redshift $z$, corrected by a peculiar velocity field, $M \propto m - 5\log(z-u)$. We describe the radial peculiar velocity field $u({\bf z})$ in terms of a set of orthogonal functions which can be derived from any convenient basis set; as an example we take them to be linear combinations of low order spherical harmonic and spherical Bessel functions. The model parameters are then found by maximizing the likelihood function for measuring a set of observed $\eta$. The predicted peculiar velocities are free of Malmquist bias in the absence of multi-streaming, provided no selection criteria are imposed on the measurement of circular velocities. This procedure can be considered as a generalized smoothing algorithm of the peculiar velocity field, and is particularly useful for comparison to the smoothed gravity field derived from full-sky galaxy redshift catalogs such as the IRAS surveys. We demonstrate the technique using a catalog of ``galaxies" derived from an N-body simulation. Increasing the resolution of the velocity smoothing beyond a certain level degrades the correlation of fitted velocities against the velocities calculated from linear theory methods, which have finite resolution,

We study the quasilinear evolution of the one-point probability density functions (PDFs) of the smoothed density and velocity fields in a cosmological gravitating system beginning with Gaussian initial fluctuations. Our analytic results are based on the Zel'dovich approximation and laminar flow. A numerical analysis extends the results into the multistreaming regime using the smoothed fields of a CDM N-body simulation. We find that the PDF of velocity, both Lagrangian and Eulerian, remains Gaussian under the laminar Zel'dovich approximation, and it is almost indistinguishable from Gaussian in the simulations. The PDF of mass density deviates from a normal distribution early in the quasilinear regime and it develops a shape remarkably similar to a lognormal distribution with one parameter, the \rms density fluctuation $\sigma$. Applying these results to currently available data we find that the PDFs of the velocity and density fields, as recovered by the \pot\ procedure from observed velocities assuming $\Omega=1$, or as deduced from a redshift survey of \iras\ galaxies assuming that galaxies trace mass, are consistent with Gaussian initial fluctuations.

We present a new method for recovering the underlying velocity field from an observed distribution of galaxies in redshift space. The method is based on a kinematic Zel'dovich relation between the velocity and density fields in redshift space. This relation is expressed in a differential equation slightly modified from the usual Poisson equation, and which depends non-trivially on $\beta \equiv \Omega^{0.6}/b$. The linear equation can be readily solved by standard techniques of separation of variables by means of spherical harmonics. One can also include a term describing the ``rocket effect" discussed by Kaiser (1987). From this redshift space information alone, one can generate a prediction of the peculiar velocity field for each harmonic ($l,m$) as a function of distance. We note that for the quadrupole and higher order moments, the equation is a boundary value problem with solutions dependent on both the interior and exterior mass distribution. However, for a shell at distance $r$, the dipole, as well as the monopole, of the velocity field in the Local Group frame is fully determined by the interior mass distribution. This implies that the shear of the measured velocity field, when fit to a dipole distortion, should be aligned and consistent with the gravity field inferred from the well determined local galaxy distribution. As a preliminary application we compute the velocity dipole of distant shells as predicted from the 1.2Jy IRAS survey compared to the measured velocity dipole on shells, as inferred from a recent POTENT analysis. The coherence between the two fields is good, yielding a best estimate of $\beta= 0.6 \pm 0.2$.

We recover the one-point probability distribution function of the {\it initial} density fluctuations (\ipdf) from the quasi-linear galaxy density field of the 1.2 Jy \iras\ redshift survey smoothed by $10\hmpc$. The recovery, using the laminar, Eulerian, Zel'dovich approximation of Nusser and Dekel, is independent of $\Omega$ and is relatively insensitive to linear galaxy biasing in the range $0.5 \leq b \leq 2$. Errors due to discrete sampling, the limited survey volume, and the method of recovering the \ipdf\ are evaluated by comparing the \ipdf\ determined from the \iras\ data with Monte-Carlo \iras-like catalogs ``observed" from $N$-body simulations of CDM models with Gaussian initial conditions. Eight sensitive statistical tests used in this comparison find the \iras\ \ipdf\ to be consistent with Gaussian. We provide observational constraints on possible deviations from Gaussianity, which should be obeyed by any theoretical model.