The Nonlinear Redshift Space Power Spectrum: Omega from Redshift Surveys

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.Comment: uuencoded compressed postscript. The preprint is also available at http://www.ast.cam.ac.uk/preprint/PrePrint.htm

The Velocity Field from Type Ia Supernovae Matches the Gravity Field from Galaxy Surveys

We compare the peculiar velocities of nearby SNe Ia with those predicted by the gravity fields of full sky galaxy catalogs. The method provides a powerful test of the gravitational instability paradigm and strong constraints on the density parameter beta = Omega^0.6/b. For 24 SNe Ia within 10,000 km/s we find the observed SNe Ia peculiar velocities are well modeled by the predictions derived from the 1.2 Jy IRAS survey and the Optical Redshift Survey (ORS). Our best $\beta$ is 0.4 from IRAS, and 0.3 from the ORS, with beta>0.7 and beta<0.15 ruled out at 95% confidence levels from the IRAS comparison. Bootstrap resampling tests show these results to be robust in the mean and in its error. The precision of this technique will improve as additional nearby SNe Ia are discovered and monitored.Comment: 16 pages (LaTex), 3 postscript figure

Relativistic Model of Detonation Transition from Neutron to Strange Matter

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.Comment: 16 pages, 4 figures. Submitted to IJMP

POTENT Reconstruction from Mark III Velocities

We present an improved POTENT method for reconstructing the velocity and mass density fields from radial peculiar velocities, test it with mock catalogs, and apply it to the Mark III Catalog. Method improvments: (a) inhomogeneous Malmquist bias is reduced by grouping and corrected in forward or inverse analyses of inferred distances, (b) the smoothing into a radial velocity field is optimized to reduce window and sampling biases, (c) the density is derived from the velocity using an improved nonlinear approximation, and (d) the computational errors are made negligible. The method is tested and optimized using mock catalogs based on an N-body simulation that mimics our cosmological neighborhood, and the remaining errors are evaluated quantitatively. The Mark III catalog, with ~3300 grouped galaxies, allows a reliable reconstruction with fixed Gaussian smoothing of 10-12 Mpc/h out to ~60 Mpc/h. We present maps of the 3D velocity and mass-density fields and the corresponding errors. The typical systematic and random errors in the density fluctuations inside 40 Mpc/h are \pm 0.13 and \pm 0.18. The recovered mass distribution resembles in its gross features the galaxy distribution in redshift surveys and the mass distribution in a similar POTENT analysis of a complementary velocity catalog (SFI), including the Great Attractor, Perseus-Pisces, and the void in between. The reconstruction inside ~40 Mpc/h is not affected much by a revised calibration of the distance indicators (VM2, tailored to match the velocities from the IRAS 1.2Jy redshift survey). The bulk velocity within the sphere of radius 50 Mpc/h about the Local Group is V_50=370 \pm 110 km/s (including systematic errors), and is shown to be mostly generated by external mass fluctuations. With the VM2 calibration, V_50 is reduced to 305 \pm 110 km/s.Comment: 60 pages, LaTeX, 3 tables and 27 figures incorporated (may print the most crucial figures only, by commenting out one line in the LaTex source

Maximum-Likelihood Comparisons of Tully-Fisher and Redshift Data: Constraints on Omega and Biasing

We compare Tully-Fisher (TF) data for 838 galaxies within cz=3000 km/sec from the Mark III catalog to the peculiar velocity and density fields predicted from the 1.2 Jy IRAS redshift survey. Our goal is to test the relation between the galaxy density and velocity fields predicted by gravitational instability theory and linear biasing, and thereby to estimate $\beta_I = \Omega^{0.6}/b_I,$ where $b_I$ is the linear bias parameter for IRAS galaxies. Adopting the IRAS velocity and density fields as a prior model, we maximize the likelihood of the raw TF observables, taking into account the full range of selection effects and properly treating triple-valued zones in the redshift-distance relation. Extensive tests with realistic simulated galaxy catalogs demonstrate that the method produces unbiased estimates of $\beta_I$ and its error. When we apply the method to the real data, we model the presence of a small but significant velocity quadrupole residual (~3.3% of Hubble flow), which we argue is due to density fluctuations incompletely sampled by IRAS. The method then yields a maximum likelihood estimate $\beta_I=0.49\pm 0.07$ (1-sigma error). We discuss the constraints on $\Omega$ and biasing that follow if we assume a COBE-normalized CDM power spectrum. Our model also yields the 1-D noise noise in the velocity field, including IRAS prediction errors, which we find to be be 125 +/- 20 km/sec.Comment: 53 pages, 20 encapsulated figures, two tables. Submitted to the Astrophysical Journal. Also available at http://astro.stanford.edu/jeff

The accretion mechanism in low-power radio galaxies

We study a sample of 44 low-luminosity radio-loud AGN, which represent a range of nuclear radio-power spanning 5 orders of magnitude, to unveil the accretion mechanism in these galaxies. We estimate the accretion rate of gas associated with their hot coronae by analyzing archival Chandra data, to derive the deprojected density and temperature profiles in a spherical approximation. Measuring the jet power from the nuclear radio-luminosity, we find that the accretion power correlates linearly with the jet power, with an efficiency of conversion from rest mass into jet power of ~0.012. These results strengthen and extend the validity of the results obtained by Allen and collaborators for 9 radio galaxies, indicating that hot gas accretion is the dominant process in FR I radio galaxies across their full range of radio-luminosity. We find that the different levels of nuclear activity are driven by global differences in the structure of the galactic hot coronae. A linear relation links the jet power with the host X-ray surface brightness. This implies that a substantial change in the jet power must be accompanied by a global change in its ISM properties, driven for example by a major merger. This correlation provides a simple widely applicable method to estimate the jet-power of a given object by observing the intensity of its host X-ray emission. To maintain the mass flow in the jet, the fraction of gas that crosses the Bondi radius reaching the accretion disk must be > 0.002. This implies that the radiative efficiency of the disk must be < 0.005, an indication that accretion in these objects occurs not only at a lower rate, but also at lower efficiency than in standard accretion disks.Comment: Accepted for publication in A&

Reconstruction Analysis of Galaxy Redshift Surveys: A Hybrid Reconstruction Method

In reconstruction analysis of galaxy redshift surveys, one works backwards from the observed galaxy distribution to the primordial density field in the same region, then evolves the primordial fluctuations forward in time with an N-body code. This incorporates assumptions about the cosmological parameters, the properties of primordial fluctuations, and the biasing relation between galaxies and mass. These can be tested by comparing the reconstruction to the observed galaxy distribution, and to peculiar velocity data. This paper presents a hybrid reconstruction method that combines the `Gaussianization'' technique of Weinberg(1992) with the dynamical schemes of Nusser & Dekel(1992) and Gramann(1993). We test the method on N-body simulations and on N-body mock catalogs that mimic the depth and geometry of the Point Source Catalog Redshift Survey and the Optical Redshift Survey. This method is more accurate than Gaussianization or dynamical reconstruction alone. Matching the observed morphology of clustering can limit the bias factor b, independent of Omega. Matching the cluster velocity dispersions and z-space distortions of the correlation function xi(s,mu) constrains the parameter beta=Omega^{0.6}/b. Relative to linear or quasi-linear approximations, a fully non-linear reconstruction makes more accurate predictions of xi(s,mu) for a given beta, thus reducing the systematic biases of beta measurements and offering further scope for breaking the degeneracy between Omega and b. It also circumvents the cosmic variance noise that limits conventional analyses of xi(s,mu). It can also improve the determination of Omega and b from joint analyses of redshift & peculiar velocity surveys as it predicts the fully non-linear peculiar velocity distribution at each point in z-space.Comment: 72 pages including 33 figures, submitted to Ap

Cold collapse and the core catastrophe

We show that a universe dominated by cold dark matter fails to reproduce the rotation curves of dark matter dominated galaxies, one of the key problems that it was designed to resolve. We perform numerical simulations of the formation of dark matter halos, each containing \gsim 10^6 particles and resolved to 0.003 times the virial radius, allowing an accurate comparison with rotation curve data. A good fit to both galactic and cluster sized halos can be achieved using the density profile rho(r) \propto [(r/r_s)^1.5(1+(r/r_s)^1.5)]^-1, where r_s is a scale radius. This profile has a steeper asymptotic slope, rho(r) \propto r^-1.5, and a sharper turnover than found by lower resolution studies. The central structure of relaxed halos that form within a hierarchical universe has a remarkably small scatter (unrelaxed halos would not host disks). We compare the results with a sample of dark matter dominated, low surface brightness (LSB) galaxies with circular velocities in the range 100-300 km/s. The rotation curves of disks within cold dark matter halos rise too steeply to match these data which require a constant mass density in the central regions. The same conclusion is reached if we compare the scale free shape of observed rotation curves with the simulation data. It is important to confirm these results using stellar rather than HI rotation curves for LSB galaxies. We test the effects of introducing a cut-off in the power spectrum that may occur in a universe dominated by warm dark matter. In this case halos form by a monolithic collapse but the final density profile hardly changes, demonstrating that the merger history does not play a role in determining the halo structure.Comment: Latex 13 pages, 4 figures. Submitted to MNRAS. High resolution colour version of figure 4 and other N-body images here: http://star-www.dur.ac.uk:80/~moore/images

Metabolic and hormonal studies of Type 1 (insulin-dependent) diabetic patients after successful pancreas and kidney transplantation

Long-term normalization of glucose metabolism is necessary to prevent or ameliorate diabetic complications. Although pancreatic grafting is able to restore normal blood glucose and glycated haemoglobin, the degree of normalization of the deranged diabetic metabolism after pancreas transplantation is still questionable. Consequently glucose, insulin, C-peptide, glucagon, and pancreatic polypeptide responses to oral glucose and i.v. arginine were measured in 36 Type 1 (insulin-dependent) diabetic recipients of pancreas and kidney allografts and compared to ten healthy control subjects. Despite normal HbA1 (7.2±0.2%; normal <8%) glucose disposal was normal only in 44% and impaired in 56% of the graft recipients. Normalization of glucose tolerance was achieved at the expense of hyperinsulinaemia in 52% of the subjects. C-peptide and glucagon were normal, while pancreatic polypeptide was significantly higher in the graft recipients. Intravenous glucose tolerance (n=21) was normal in 67% and borderline in 23%. Biphasic insulin release was seen in patients with normal glucose tolerance. Glucose tolerance did not deteriorate up to 7 years post-transplant. In addition, stress hormone release (cortisol, growth hormone, prolactin, glucagon, catecholamines) to insulin-induced hypoglycaemia was examined in 20 graft recipients and compared to eight healthy subjects. Reduced blood glucose decline indicates insulin resistance, but glucose recovery was normal, despite markedly reduced catecholamine and glucagon release. These data demonstrate the effectiveness of pancreatic grafting in normalizing glucose metabolism, although hyperinsulinaemia and deranged counterregulatory hormone response are observed frequently

Fast Large Volume Simulations of the 21 cm Signal from the Reionization and pre-Reionization Epochs

While limited to low spatial resolution, the next generation low-frequency radio interferometers that target 21 cm observations during the era of reionization and prior will have instantaneous fields-of-view that are many tens of square degrees on the sky. Predictions related to various statistical measurements of the 21 cm brightness temperature must then be pursued with numerical simulations of reionization with correspondingly large volume box sizes, of order 1000 Mpc on one side. We pursue a semi-numerical scheme to simulate the 21 cm signal during and prior to Reionization by extending a hybrid approach where simulations are performed by first laying down the linear dark matter density field, accounting for the non-linear evolution of the density field based on second-order linear perturbation theory as specified by the Zel'dovich approximation, and then specifying the location and mass of collapsed dark matter halos using the excursion-set formalism. The location of ionizing sources and the time evolving distribution of ionization field is also specified using an excursion-set algorithm. We account for the brightness temperature evolution through the coupling between spin and gas temperature due to collisions, radiative coupling in the presence of Lyman-alpha photons and heating of the intergalactic medium, such as due to a background of X-ray photons. The hybrid simulation method we present is capable of producing the required large volume simulations with adequate resolution in a reasonable time so a large number of realizations can be obtained with variations in assumptions related to astrophysics and background cosmology that govern the 21 cm signal.Comment: 14 pages and 15 figures. New version to match accepted version for MNRAS. Code available in: http://www.SimFast21.or