We present optical and ultraviolet spectra of SN 1885 (S And), visible in absorption against the bulge of the Andromeda galaxy (M31), using the Hubble Space Telescope's STIS spectrograph to probe the three dimensional arrangement of the supernova debris. Spectra covering 2900 to 5700 A taken using six 0.2" slit positions in two orientations show broad Ca II absorption with blue and red radial velocities out to at least 11,500 km/s, consistent with HST Ca II images of S And. Enhanced Ca II absorption is seen between 2000 to 6000 km/s suggestive of a Ca II-rich shell. The spectra also show strong, asymmetric Ca I 4227 A absorption extending out to +12,400 km/s, along with weak Fe I 3720 A absorption in a shell with velocities between 2000 and 9000 km/s. Ultraviolet spectra obtained revealed weak broad absorption shortward of 3000 A consistent with model predictions. The STIS spectra, together with previous HST images, show a layered structure with a well defined Ca-rich outer edge indicative of a delayed detonation phase. The remnant's clumpy inner Ca-rich shell plus only a handful of Fe-rich plumes is unlike morphologies expected from dynamical or violent merger scenarios. The small number of Fe-rich plumes suggest their formation during a deflagration phase by Rayleigh-Taylor instabilities but less well developed without extended mixing as expected from hydrodynamic calculations. The suppression of strong Rayleigh-Taylor instabilities is possibly the result of strong magnetic fields. We propose SN 1885 was an off-center, delayed detonation and slightly subluminous SN Ia similar to SN 1986g.

We propose a physically based, analytic model for intergalactic filaments during the first gigayear of the universe. The structure of a filament is based upon a gravitationally bound, isothermal cylinder of gas. The model successfully predicts for a cosmological simulation the total mass per unit length of a filament (dark matter plus gas) based solely upon the sound speed of the gas component, contrary to the expectation for collisionless dark matter aggregation. It argues that the gas, through its hydrodynamic properties, plays a key role in filament structure rather than being a passive passenger in a preformed dark matter potential. The dark matter of a galaxy follows the classic equation of collapse of a spherically symmetric overdensity in an expanding universe. In contrast, the gas usually collapses more slowly. The relative rates of collapse of these two components for individual galaxies can explain the varying baryon deficits of the galaxies under the assumption that matter moves along a single filament passing through the galaxy centre, rather than by spherical accretion. The difference in behaviour of the dark matter and gas can be simply and plausibly related to the model. The range of galaxies studied includes that of the so-called "too big to fail" galaxies, which are thought to be problematic for the standard Lambda-CDM model of the universe. The isothermal-cylinder model suggests a simple explanation for why these galaxies are, unaccountably, missing from the night sky.

We present Hubble Space Telescope (HST) ultraviolet Fe I and Fe II images of the remnant of Supernova 1885 (S And) which is observed in absorption against the bulge of the Andromeda galaxy, M31. We compare these Fe I and Fe II absorption line images to previous HST absorption images of S And, of which the highest quality and theoretically cleanest is Ca II H & K. Because the remnant is still in free expansion, these images provide a 2D look at the distribution of iron synthesized in this probable Type Ia explosion, thus providing insights and constraints for theoretical SN Ia models. The Fe I images show extended absorption offset to the east from the remnant's center as defined by Ca II images and is likely an ionization effect due to self-shielding. More significant is the remnant's apparent Fe II distribution which consists of four streams or plumes of Fe-rich material seen in absorption that extend from remnant center out to about 10,000 km/s. This is in contrast to the remnant's Ca II absorption, which is concentrated in a clumpy, roughly spherical shell at 1000 to 5000 km/s but which extends out to 12,500 km/s. The observed distributions of Ca and Fe rich ejecta in the SN 1885 remnant are consistent with delayed detonation white dwarf models. The largely spherical symmetry of the Ca-rich layer argues against a highly anisotropic explosion as might result from a violent merger of two white dwarfs.

This paper presents a concise derivation of a new set of solutions for the interior structure of accreting, rotating black holes. The solutions are conformally stationary, axisymmetric, and conformally separable. Hyper-relativistic counter-streaming between freely-falling collisionless ingoing and outgoing streams leads to mass inflation at the inner horizon, followed by collapse. The solutions fail at an exponentially tiny radius, where the rotational motion of the streams becomes comparable to their radial motion. The papers provide a fully nonlinear, dynamical solution for the interior structure of a rotating black hole from just above the inner horizon inward, down to a tiny scale.

(Abridged) A solution is obtained for the interior structure of an uncharged rotating black hole that accretes a collisionless fluid. The solution is conformally stationary, axisymmetric, and conformally separable, possessing a conformal Killing tensor. Hyper-relativistic counter-streaming between collisionless ingoing and outgoing streams drives inflation at (just above) the inner horizon, followed by collapse. As ingoing and outgoing streams approach the inner horizon, they focus into twin narrow beams directed along the ingoing and outgoing principal null directions, regardless of the initial angular motions of the streams. The radial energy-momentum of the counter-streaming beams gravitationally accelerates the streams even faster along the principal directions, leading to exponential growth in the streaming density and pressure, and in the Weyl curvature and mass function. At exponentially large density and curvature, inflation stalls, and the spacetime collapses. As the spacetime collapses, the angular motions of the freely-falling streams grow. When the angular motion has become comparable to the radial motion, which happens when the conformal factor has shrunk to an exponentially tiny scale, conformal separability breaks down, and the solution fails. The condition of conformal separability prescribes the form of the ingoing and outgoing accretion flows incident on the inner horizon. The dominant radial part of the solution holds provided that the densities of ingoing and outgoing streams incident on the inner horizon are uniform, independent of latitude; that is, the accretion flow is "monopole." The sub-dominant angular part of the solution requires a special non-radial pattern of angular motion of streams incident on the inner horizon. The prescribed angular pattern cannot be achieved if the collisionless streams fall freely from outside the horizon.

This paper extends to the case of charged rotating black holes the conformally stationary, axisymmetric, conformally separable solutions presented for uncharged rotating black holes in a companion paper. In the present paper, the collisionless fluid accreted by the black hole may be charged. The charge of the black hole is determined self-consistently by the charge accretion rate. As in the uncharged case, hyper-relativistc counter-streaming between ingoing and outgoing streams drives inflation at (just above) the inner horizon, followed by collapse. If both ingoing and outgoing streams are charged, then conformal separability holds during early inflation, but fails as inflation develops. If conformal separability is imposed throughout inflation and collapse, then only one of the ingoing and outgoing streams can be charged: the other must be neutral. Conformal separability prescribes a hierarchy of boundary conditions on the ingoing and outgoing streams incident on the inner horizon. The dominant radial boundary conditions require that the incident ingoing and outgoing number densities be uniform with latitude, but the charge per particle must vary with latitude such that the incident charge densities vary in proportion to the radial electric field. The sub-dominant angular boundary conditions require specific forms of the incident number- and charge-weighted angular motions. If the streams fall freely from outside the horizon, then the prescribed angular conditions can be achieved by the charged stream, but not by the neutral stream. Thus, as in the case of an uncharged black hole, the neutral stream must be considered to be delivered ad hoc to just above the inner horizon.

The Schweizer-Middleditch star, located behind the SN 1006 remnant and near its center in projection, provides the opportunity to study cold, expanding ejecta within the SN 1006 shell through UV absorption. Especially notable is an extremely sharp red edge to the Si II 1260 Angstrom feature, which stems from the fastest moving ejecta on the far side of the SN 1006 shell--material that is just encountering the reverse shock. Comparing HST far-UV spectra obtained with COS in 2010 and with STIS in 1999, we have measured the change in this feature over the intervening 10.5-year baseline. We find that the sharp red edge of the Si II feature has shifted blueward by 0.19 +/- 0.05 Angstroms, which means that the material hitting the reverse shock in 2010 was moving slower by 44 +/- 11 km/s than the material that was hitting it in 1999, a change corresponding to - 4.2 +/- 1.0 km/s/yr. This is the first observational confirmation of a long-predicted dynamic effect for a reverse shock: that the shock will work its way inward through expanding supernova ejecta and encounter ever slower material as it proceeds. We also find that the column density of shocked Si II (material that has passed through the reverse shock) has decreased by 7 +/- 2% over the ten-year period. The decrease could indicate that in this direction the reverse shock has been ploughing through a dense clump of Si,leading to pressure and density transients.

The purpose of this paper is to present a number of proposals about the interior structure of a rotating black hole that is accreting slowly, but in an arbitrary time- and space-dependent fashion. The proposals could potentially be tested with numerical simulations. Outgoing and ingoing particles free-falling in the parent Kerr geometry become highly focused along the principal outgoing and ingoing null directions as they approach the inner horizon, triggering the mass inflation instability. The original arguments of Barrabes, Israel & Poisson (1990) regarding inflation in rotating black holes are reviewed, and shown to be based on Raychauduri's equation applied along the outgoing and ingoing null directions. It is argued that gravitational waves should behave in the geometric optics limit, and consequently that the spacetime should be almost shear-free. A full set of shear-free equations is derived. A specific line-element is proposed, which is argued should provide a satisfactory approximation during early inflation. Finally, it is argued that super-Planckian collisions between outgoing and ingoing particles will lead to entropy production, bringing inflation to an end, and precipitating collapse.

Using a cosmological simulation at redshift 5, we find that the baryon-rich cores of intergalactic filaments extending outward from galaxies commonly form isothermal gas cylinders in regions favorable to their formation. The central gas density is typically about 500 times the cosmic mean total density, and the temperature is typically 1-2 times 10^4 K, just above the Lyman alpha cooling floor. These findings argue that the hydrodynamic properties of the gas are more important than the dark matter in determining this structure. It is noteworthy that the temperature and ionization state of the gas completely determine a finite total mass per unit length of an isothermal cylinder. Our findings may have implications for understanding the "cold mode" mechanism of gas transport into galaxies.

Particles dropped into a rotating black hole can collide near the inner horizon with enormous energies. The entropy produced by these collisions can be several times larger than the increase in the horizon entropy due to the addition of the particles. In this paper entropy is produced by releasing large numbers of neutrons near the outer horizon of a rotating black hole such that they collide near the inner horizon at energies similar to those achieved at the Relativistic Heavy Ion Collider. The increase in horizon entropy is approximately 80 per dropped neutron pair, while the entropy produced in the collisions is 160 per neutron pair. The collision entropy is produced inside the horizon, so this excess entropy production does not violate Bousso's bound limiting the entropy that can go through the black hole's horizon. The generalized laws of black hole thermodynamics are obeyed. No individual observer inside the black hole sees a violation of the second law of thermodynamics

If you fall into a real astronomical black hole (choosing a supermassive black hole, to make sure that the tidal forces don't get you first), then you will probably meet your fate not at a central singularity, but rather in the exponentially growing, relativistic counter-streaming instability at the inner horizon first pointed out by Poisson & Israel (1990), who called it mass inflation. The chief purpose of this paper is to present a clear exposition of the physical cause and consequence of inflation in spherical, charged black holes. Inflation acts like a particle accelerator in that it accelerates cold ingoing and outgoing streams through each other to prodigiously high energies. Inflation feeds on itself: the acceleration is powered by the gravity produced by the streaming energy. The paper: (1) uses physical arguments to develop simple approximations that follow the evolution of inflation from ignition, through inflation itself, to collapse; (2) confirms that the simple approximations capture accurately the results of fully nonlinear one- and two-fluid self-similar models; (3) demonstrates that, counter-intuitively, the smaller the accretion rate, the more rapidly inflation exponentiates; (4) shows that in single perfect-fluid models, inflation occurs only if the sound speed equals the speed of light, supporting the physical idea that inflation in single fluids is driven by relativistic counter-streaming of waves; (5) shows that what happens during inflation up to the Planck curvature depends not on the distant past or future, but rather on events happening only a few hundred black hole crossing times into the past or future; (6) shows that, if quantum gravity does not intervene, then the generic end result of inflation is not a general relativistic null singularity, but rather a spacelike singularity at zero radius.

Using cosmological simulations with a dynamic range in excess of 10 million, we study the transport of gas mass and angular momentum through the circumnuclear region of a disk galaxy containing a supermassive black hole (SMBH). The simulations follow fueling over relatively quiescent phases of the galaxy's evolution (no mergers) and without feedback from active galactic nuclei (AGNs), as part of the first stage of using state-of-the-art, high-resolution cosmological simulations to model galaxy and black hole co-evolution. We present results from simulations at different redshifts (z=6, 4, and 3) and three different black hole masses (30 million, 90 million, and 300 million solar masses; at z=4), as well as a simulation including a prescription that approximates optically thick cooling in the densest regions. The interior gas mass throughout the circumnuclear disk shows transient and chaotic behavior as a function of time. The Fourier transform of the interior gas mass follows a power law with slope -1 throughout the region, indicating that, in the absence of the effects of galaxy mergers and AGN feedback, mass fluctuations are stochastic with no preferred timescale for accretion over the duration of each simulation (~ 1-2 Myr). The angular momentum of the gas disk changes direction relative to the disk on kiloparsec scales over timescales less than 1 Myr, reflecting the chaotic and transient gas dynamics of the circumnuclear region. Infalling clumps of gas, which are driven inward as a result of the dynamical state of the circumnuclear disk, may play an important role in determining the spin evolution of an SMBH, as has been suggested in stochastic accretion scenarios.

A detailed non-linear analysis of the internal structure of spherical, charged black holes that are accreting scalar matter is performed in the framework of the Brans-Dicke theory of gravity. We choose the lowest value of the Brans-Dicke parameter that is compatible with observational constraints. First, the homogeneous approximation is used. It indicates that mass inflation occurs and that the variations of the Brans-Dicke scalar inside the black hole, which could in principle be large in the absence of mass inflation, become small when mass inflation does occur. Then, a full non-linear numerical study of the black hole interior perturbed by a self-gravitating massless uncharged scalar-field is performed. We use an algorithm with adaptive mesh refinement capabilities. In this way, the changes in the internal structure of the black hole caused by mass inflation are determined, as well as the induced variations of the Brans-Dicke scalar, confirming, qualitatively, the indications given by the homogeneous approximation.

We illustrate and discuss the view seen by an observer inside the horizon of a Schwarzschild black hole. The view as the observer approaches the central singularity is of particular interest because, according to ideas arising from "observer complementarity," points in opposite directions along the observer's past lightcone are at "the edge of locality," where standard flat-space quantum-field-theory commutation rules may be at the brink of failure. Near the singularity, the observer's view is aberrated by the diverging tidal force into a horizontal plane. The view in the horizontal plane is highly blueshifted, but all directions other than horizontal appear highly redshifted. We argue that the affine distance provides a canonical measure of distance along a light ray from emitter to observer. Since the affine distance is not directly measurable by the observer, we also consider perceptual distances, and argue that the trinocular distance (binocular distance is inadequate) provides an estimate of affine distance that would allow tree-leaping apes to survive in highly curved spacetime.

Heating processes inside large black holes can produce tremendous amounts of entropy. Locality requires that this entropy adds on space-like surfaces, but the resulting entropy (10^10 times the Bekenstein-Hawking entropy in an example presented in the companion paper) exceeds the maximum entropy that can be accommodated by the black hole's degrees of freedom. Observer complementarity, which proposes a proliferation of non-local identifications inside the black hole, allows the entropy to be accommodated as long as individual observers inside the black hole see less than the Bekenstein-Hawking entropy. In the specific model considered with huge entropy production, we show that individual observers do see less than the Bekenstein-Hawking entropy, offering strong support for observer complementarity.

We show that the entropy created by Ohmic dissipation inside an accreting charged black hole may exceed the Bekenstein-Hawking entropy by a large factor. If the black hole subsequently evaporates, radiating only the Bekenstein-Hawking entropy, then the black hole appears to destroy entropy, violating the second law of thermodynamics. A companion paper discusses the implications of this startling result. Bousso's covariant entropy bound is not violated.

The internal structure of a slowly rotating, charged black hole that is undergoing mass inflation at its inner horizon is derived. The equations governing the angular behavior decouple from the radial behavior, so all conclusions regarding inflation in a spherical charged black hole carry through unchanged for a slowly-rotating black hole. Quantities inflate only in the radial direction, not in the angular direction. Exact self-similar solutions are obtained. For sufficiently small accretion rates, the instantaneous angular motion of the accretion flow has negligible effect on the angular spacetime structure of the black hole, even if the instantaneous angular momentum of the accretion flow is large and arbitrarily oriented.

As galaxy surveys become larger and more complex, keeping track of the completeness, magnitude limit, and other survey parameters as a function of direction on the sky becomes an increasingly challenging computational task. For example, typical angular masks of the Sloan Digital Sky Survey contain about N=300,000 distinct spherical polygons. Managing masks with such large numbers of polygons becomes intractably slow, particularly for tasks that run in time O(N^2) with a naive algorithm, such as finding which polygons overlap each other. Here we present a "divide-and-conquer" solution to this challenge: we first split the angular mask into predefined regions called "pixels," such that each polygon is in only one pixel, and then perform further computations, such as checking for overlap, on the polygons within each pixel separately. This reduces O(N^2) tasks to O(N), and also reduces the important task of determining in which polygon(s) a point on the sky lies from O(N) to O(1), resulting in significant computational speedup. Additionally, we present a method to efficiently convert any angular mask to and from the popular HEALPix format. This method can be generically applied to convert to and from any desired spherical pixelization. We have implemented these techniques in a new version of the mangle software package, which is freely available at http://space.mit.edu/home/tegmark/mangle/, along with complete documentation and example applications. These new methods should prove quite useful to the astronomical community, and since mangle is a generic tool for managing angular masks on a sphere, it has the potential to benefit terrestrial mapmaking applications as well.

The semi-classical nature of braneworld cosmological models does not account for any quantum gravitational effects. In this letter we use the gauge/gravity correspondence to argue that quantum string corrections cannot be ignored in any study of braneworld stability. As an example, we find, by analysing the quantum gravitational backreaction, that a closed universe is unstable to radiation into the bulk.

Using a high resolution cosmological simulation of reionization we have examined the differing structures formed by gas and dark matter at a redshift of 5.1. Baryon-rich regions form a small number of filaments, which connect the largest galaxies in the simulation. More detailed examination of the ten largest galaxies reveals long, slender gaseous filaments about 5 proper kpc in width radiating from the galaxy centers. Extending out from each filament are a few smooth, thin, nearly planar gaseous sheets. By contrast, the dark matter concentrates into quasi-spherical bodies. The results have implications for our understanding of structure formation in the early universe and of the Lyman alpha forest.

Using a hydrodynamic adaptive mesh refinement code, we simulate the growth and evolution of a galaxy, which could potentially host a supermassive black hole, within a cosmological volume. Reaching a dynamical range in excess of 10 million, the simulation follows the evolution of the gas structure from super-galactic scales all the way down to the outer edge of the accretion disk. Here, we focus on global instabilities in the self-gravitating, cold, turbulence-supported, molecular gas disk at the center of the model galaxy, which provide a natural mechanism for angular momentum transport down to sub-pc scales. The gas density profile follows a power-law scaling as r^-8/3, consistent with an analytic description of turbulence in a quasi-stationary circumnuclear disk. We analyze the properties of the disk which contribute to the instabilities, and investigate the significance of instability for the galaxy's evolution and the growth of a supermassive black hole at the center.

The theory of perturbations of spherically symmetric self-similar black holes is presented, in the Newman-Penrose formalism. It is shown that the wave equations for gravitational, electromagnetic, and scalar waves are separable, though not decoupled. A generalization of the Teukolsky equation is given. Monopole and dipole modes are treated. The Newman-Penrose wave equations governing polar and axial spin-0 perturbations are explored.

We present HST images of the remnant of SN 1885 seen in absorption against M31's bulge via resonance lines of Ca I, Ca II, Fe I, and Fe II. Viewed in CaII H & K line absorption, the remnant appears as a nearly black circular spot with an outermost angular radius of 0.40" +/- 0.025" implying r = 1.52 pc and a 120 yr average expansion velocity of 12400 +/-1400 km/s. The strongest Ca II absorption is organized in a broken ring structure with a radius of 0.20" (=6000 km/s) with several apparent absorption `clumps' of an angular size near the pixel scale of 0.05" (= 1500 km/s). The detection of Ca II clumps is the first direct evidence for some instabilities and the existence of a deflagration phase in SNe Ia or, alternatively, mixing induced by radioactive decay of 56^Ni over time scales of seconds or days. However, the degree of mixing allowed by the observed images is much smaller than current 3D calculations for Rayleigh-Taylor dominated deflagration fronts. The images also require a central region of no or little Ca but iron group elements indicative of burning under sufficiently high densities for electron capture taking place, i.e., burning prior to a significant pre-expansion of the WD.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

This paper presents an under-appreciated way to conceptualize stationary black holes, which we call the river model. The river model is mathematically sound, yet simple enough that the basic picture can be understood by non-experts. %that can by understood by non-experts. In the river model, space itself flows like a river through a flat background, while objects move through the river according to the rules of special relativity. In a spherical black hole, the river of space falls into the black hole at the Newtonian escape velocity, hitting the speed of light at the horizon. Inside the horizon, the river flows inward faster than light, carrying everything with it. We show that the river model works also for rotating (Kerr-Newman) black holes, though with a surprising twist. As in the spherical case, the river of space can be regarded as moving through a flat background. However, the river does not spiral inward, as one might have anticipated, but rather falls inward with no azimuthal swirl at all. Instead, the river has at each point not only a velocity but also a rotation, or twist. That is, the river has a Lorentz structure, characterized by six numbers (velocity and rotation), not just three (velocity). As an object moves through the river, it changes its velocity and rotation in response to tidal changes in the velocity and twist of the river along its path. An explicit expression is given for the river field, a six-component bivector field that encodes the velocity and twist of the river at each point, and that encapsulates all the properties of a stationary rotating black hole.

We show how to estimate the covariance of the power spectrum of a statistically homogeneous and isotropic density field from a single periodic simulation, by applying a set of weightings to the density field, and by measuring the scatter in power spectra between different weightings. We recommend a specific set of 52 weightings containing only combinations of fundamental modes, constructed to yield a minimum variance estimate of the covariance of power. Numerical tests reveal that at nonlinear scales the variance of power estimated by the weightings method substantially exceeds that estimated from a simple ensemble method. We argue that the discrepancy is caused by beat-coupling, in which products of closely spaced Fourier modes couple by nonlinear gravitational growth to the beat mode between them. Beat-coupling appears whenever nonlinear power is measured from Fourier modes with a finite spread of wavevector, and is therefore present in the weightings method but not the ensemble method. Beat-coupling inevitably affects real galaxy surveys, whose Fourier modes have finite width. Surprisingly, the beat-coupling contribution dominates the covariance of power at nonlinear scales, so that, counter-intuitively, it is expected that the covariance of nonlinear power in galaxy surveys is dominated not by small scale structure, but rather by beat-coupling to the largest scales of the survey.

We measure the covariance of the non-linear matter power spectrum from N-body simulations using two methods. In the first case, the covariance of power is estimated from the scatter over many random realizations of the density field. In the second, we use a novel technique to measure the covariance matrix from each simulation individually by re-weighting the density field with a carefully chosen set of functions. The two methods agree at linear scales, but unexpectedly they disagree substantially at increasingly non-linear scales. Moreover, the covariance of non-linear power measured using the re-weightings method changes with box size. The numerical results are consistent with an explanation given in a companion paper, which argues that the cause of the discrepancy is beat-coupling, in which products of Fourier modes separated by a small wavevector couple by gravitational growth to the large-scale beat mode between them. We calculate the information content of the non-linear power spectrum (about the amplitude of the initial, linear power spectrum) using both methods and confirm the result of a previous paper, that at translinear scales the power spectrum contains little information over and above that in the linear power spectrum, but that there is a marked increase in information at non-linear scales. We suggest that, in real galaxy surveys, the covariance of power at non-linear scales is likely to be dominated by beat-coupling to the largest scales of the survey and that, as a result, only part of the information potentially available at non-linear scales is actually measurable from real galaxy surveys.

We report a high resolution, far-ultraviolet, STIS E140M spectrum of the strong, broad Si II, III, and IV features produced by the ejecta of SN1006 seen in absorption against the background Schweizer-Middleditch star. The spectrum confirms the extreme sharpness of the red edge of the redshifted Si II 1260 A feature, supporting the idea that this edge represents the location of the reverse shock moving into the freely expanding ejecta. The expansion velocity of ejecta at the reverse shock is measured to be 7026 +-3(relative) +-10(absolute) km/s. If the shock model is correct, then the expansion velocity should be decreasing at the observable rate of 2.7 +-0.1 km/s per year. The pre-shock velocity, post-shock velocity, and post-shock velocity dispersion are all measured from the Si II 1260 A feature, and consistency of these velocities with the shock jump conditions implies that there is little or no electron heating in this fast (2680 km/s) Si-rich shock.

Absorption-line spectroscopy is an effective probe for cold ejecta within an SNR, provided that suitable background UV sources can be identified. For the SN 1006 remnant we have identified four such sources, in addition to the much-studied Schweitzer-Middleditch (SM) star. We have used STIS on HST to obtain UV spectra of all four sources, to study "core samples" of the SN 1006 interior. The line of sight closest to the center of the SNR shell, passing only 2.0 arcmin away, is to a V = 19.5 QSO at z = 1.026. Its spectrum shows broad Fe II absorption lines, asymmetric with red wings broader than blue. The similarity of these profiles to those seen in the SM star, which is 2.8 arcmin from the center in the opposite direction, confirms the existence of a bulge on the far side of SN 1006. The Fe II equivalent widths in the QSO spectrum are ~ 50% greater than in the SM star, suggesting that somewhat more iron may be present within SN 1006 than studies of the SM star alone have indicated, but this is still far short of what most SNIa models require. The absorption spectrum against a brighter z = 0.337 QSO seen at 57% of the shell radius shows broad silicon absorption lines but no iron other than narrow, probably interstellar lines. The cold iron expanding in this direction must be confined within v <~ 5200 km/s, also consistent with a high-velocity bulge on the far side only. The broad silicon lines indicate that the silicon layer has expanded beyond this point, and that it has probably been heated by a reverse shock. Finally, the spectra of two ~ A0V stars near the southern shell rim show no broad or unusually strong absorption lines, suggesting that the low-ionization ejecta are confined within 83% of the shell radius, at least at the azimuths of these background sources.

We use cosmological simulations to explore the large-scale effects of reionization. Since reionization is a process that involves a large dynamic range - from galaxies to rare bright quasars - we need to be able to cover a significant volume of the universe in our simulation without losing hte important small scale effects from galaxies. Here we have taken an approach that uses clumping factors derived from small scale simulations to approximate the radiative transfer on the sub-cell scales. Using this technique, we can cover a simulation size up to $1280 h^{-1} Mpc$ with $10 h^{-1} Mpc$ cells. This allows us to construct synthetic spectra of quasars similar to observed spectra of SDSS quasars at high reshifts and compare them to the observational data. These spectra can then be analyzed for HII region sizes, the presence of the Gunn-Peterson trough and the Lyman-$\alpha$ forest.

We present a new, galaxy-halo model of large-scale structure, in which the galaxies entering a given sample are the fundamental objects. Haloes attach to galaxies, in contrast to the standard halo model, in which galaxies attach to haloes. The galaxy-halo model pertains mainly to the relationships between the power spectra of galaxies and mass, and their cross-power spectrum. With surprisingly little input, an intuition-aiding approximation to the galaxy-matter cross-correlation coefficient R(k) emerges, in terms of the halo mass dispersion. This approximation seems valid to mildly non-linear scales (k < ~3 h/Mpc), allowing measurement of the bias and the matter power spectrum from measurements of the galaxy and galaxy-matter power spectra (or correlation functions). This is especially relevant given the recent advances in precision in measurements of the galaxy-matter correlation function from weak gravitational lensing. The galaxy-halo model also addresses the issue of interpreting the galaxy-matter correlation function as an average halo density profile, and provides a simple description of galaxy bias as a function of scale.

(Abridged) This is the second of two companion papers on the interior structure of self-similar accreting charged black holes. In the first paper, the black hole was allowed to accrete only a single fluid of charged baryons. In this second paper, the black hole is allowed to accrete in addition a neutral fluid of almost non-interacting dark matter. Relativistic streaming between outgoing baryons and ingoing dark matter leads to mass inflation near the inner horizon. When enough dark matter has been accreted that the center of mass frame near the inner horizon is ingoing, then mass inflation ceases and the fluid collapses to a central singularity. A null singularity does not form on the Cauchy horizon. Although the simultaneous presence of ingoing and outgoing fluids near the inner horizon is essential to mass inflation, reducing one or other of the ingoing dark matter or outgoing baryonic streams to a trace relative to the other stream makes mass inflation more extreme, not the other way round as one might naively have expected. Consequently, if the dark matter has a finite cross-section for being absorbed into the baryonic fluid, then the reduction of the amount of ingoing dark matter merely makes inflation more extreme, the interior mass exponentiating more rapidly and to a larger value before mass inflation ceases. However, if the dark matter absorption cross-section is effectively infinite at high collision energy, so that the ingoing dark matter stream disappears completely, then the outgoing baryonic fluid can drop through the Cauchy horizon. In all cases, the solutions do not depend on what happens in the infinite past or future. We discuss in some detail the physical mechanism that drives mass inflation.

An extensive investigation is made of the interior structure of self-similar accreting charged black holes. In this, the first of two papers, the black hole is assumed to accrete a charged, electrically conducting, relativistic baryonic fluid. The mass and charge of the black hole are generated self-consistently by the accreted material. The accreted baryonic fluid undergoes one of two possible fates: either it plunges directly to the spacelike singularity at zero radius, or else it drops through the Cauchy horizon. The baryons fall directly to the singularity if the conductivity either exceeds a certain continuum threshold kappa_oo, or else equals one of an infinite spectrum kappa_n of discrete values. Between the discrete values kappa_n, the solution is characterized by the number of times that the baryonic fluid cycles between ingoing and outgoing. If the conductivity is at the continuum threshold kappa_oo, then the solution cycles repeatedly between ingoing and outgoing, displaying a discrete self-similarity reminiscent of that observed in critical collapse. Below the continuum threshold kappa_oo, and except at the discrete values kappa_n, the baryonic fluid drops through the Cauchy horizon, and in this case undergoes a shock, downstream of which the solution terminates at an irregular sonic point where the proper acceleration diverges, and there is no consistent self-similar continuation to zero radius. As far as the solution can be followed inside the Cauchy horizon, the radial direction is timelike. If the radial direction remains timelike to zero radius (which cannot be confirmed because the self-similar solutions terminate), then there is presumably a spacelike singularity at zero radius inside the Cauchy horizon, which is distinctly different from the vacuum solution for a charged black hole.

We use an ensemble of N-body simulations of the currently favoured (concordance) cosmological model to measure the amount of information contained in the non-linear matter power spectrum about the amplitude of the initial power spectrum. Two surprising results emerge from this study: (i) that there is very little independent information in the power spectrum in the translinear regime (k ~ 0.2-0.8 Mpc/h at the present day) over and above the information at linear scales and (ii) that the cumulative information begins to rise sharply again with increasing wavenumber in the non-linear regime. In the fully non-linear regime, the simulations are consistent with no loss of information during translinear and non-linear evolution. If this is indeed the case then the results suggest a picture in which translinear collapse is very rapid, and is followed by a bounce prior to virialization, impelling a wholesale revision of the HKLM-PD formalism.

This paper and its companion form an extended version of notes provided to participants in the Valencia September 2004 summer school on Data Analysis in Cosmology. The papers offer a pedagogical introduction to the problem of estimating the power spectrum from galaxy surveys. The intention is to focus on concepts rather than on technical detail, but enough mathematics is provided to point the student in the right direction. This first paper presents background material. It collects some essential definitions, discusses traditional methods for measuring power, notably the Feldman-Kaiser-Peacock (1994) method, and introduces Bayesian analysis, Fisher matrices, and maximum likelihood. For pedagogy and brevity, several derivations are set as exercises for the reader. At the summer school, multiple choice questions, included herein, were used to convey some didactic ideas, and provoked a little lively debate.

This second paper of two companion papers on the estimation of power spectra specializes to the topic of estimating galaxy power spectra at large, linear scales using maximum likelihood methods. As in the first paper, the aims are pedagogical, and the emphasis is on concepts rather than technical detail. The paper covers most of the salient issues, including selection functions, likelihood functions, Karhunen-Loeve compression, pair-integral bias, Local Group flow, angular or radial systematics (arising for example from extinction), redshift distortions, quadratic compression, decorrelation, and disentanglement. The procedures are illustrated with results from the IRAS PSCz survey. Most of the PSCz graphics included in this paper have not been published elsewhere.

We have developed an algorithm to find haloes in an N-body dark matter simulation, called VOBOZ (VOronoi BOund Zones), which has as little dependence on free parameters as we can manage. By using the Voronoi diagram, we achieve nonparametric, `natural' measurements of each particle's density and set of neighbors. We then eliminate much of the ambiguity in merging sets of particles together by identifying every possible density peak, and measuring the probability that each does not arise from Poisson noise. This probability could be interpreted as a measure of a halo's sub-haloness, and might be of use to researchers of the halo model. The first parameter in VOBOZ controls the subtlety of particle unbinding, and may be eliminated if one is cavalier with processor time; even if one is not, the results saturate to the parameter-free answer when the parameter is sufficiently small. The only parameter which remains, an outer density cut-off, does not influence whether or not haloes are identified, nor does it have any effect on subhaloes; it only affects the masses returned for supercluster haloes.

We show that the amount of particle production in an arbitrary cosmological background can be determined using only the late-time positive-frequency modes. We don't refer to modes at early times, so there is no need for a Bogolubov transformation. We also show that particle production can be extracted from the Feynman propagator in an auxiliary spacetime. This provides a first-quantized formalism for computing particle production which, unlike conventional Bogolubov transformations, may be amenable to a string-theoretic generalization.

By comparing the PSCz galaxy power spectrum with the results of nested pure dark matter N-body simulations, we try to understand how infrared-selected galaxies populate dark-matter haloes, paying special attention to the method of halo identification in the simulations. We thus test the hypothesis that baryonic physics negligibly affects the distribution of galaxies down to the smallest scales yet observed. We are successful in reproducing the PSCz power spectrum on scales < ~40 h/Mpc, near our resolution limit, by imposing a central density cut-off on simulated haloes, which gives a rough minimum mass and circular velocity of haloes in which PSCz galaxies formed.

We measure the large-scale real-space power spectrum P(k) using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 square degrees with mean redshift z~0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.02 h/Mpc < k < 0.3h/Mpc. We pay particular attention to modeling, quantifying and correcting for potential systematic errors, nonlinear redshift distortions and the artificial red-tilt caused by luminosity-dependent bias. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k<0.1h/Mpc, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the WMAP satellite. As a simple characterization of the data, our measurements are well fit by a flat scale-invariant adiabatic cosmological model with h Omega_m =0.201+/- 0.017 and L* galaxy sigma_8=0.89 +/- 0.02 when fixing the baryon fraction Omega_b/Omega_m=0.17 and the Hubble parameter h=0.72; cosmological interpretation is given in a companion paper.

We investigate anomalies reported in the Cosmic Microwave Background maps from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite on very large angular scales and discuss possible interpretations. Three independent anomalies involve the quadrupole and octopole: 1. The cosmic quadrupole on its own is anomalous at the 1-in-20 level by being low (the cut-sky quadrupole measured by the WMAP team is more strikingly low, apparently due to a coincidence in the orientation of our Galaxy of no cosmological significance); 2. The cosmic octopole on its own is anomalous at the 1-in-20 level by being very planar; 3. The alignment between the quadrupole and octopole is anomalous at the 1-in-60 level. Although the a priori chance of all three occurring is 1 in 24000, the multitude of alternative anomalies one could have looked for dilutes the significance of such a posteriori statistics. The simplest small universe model where the universe has toroidal topology with one small dimension of order half the horizon scale, in the direction towards Virgo, could explain the three items above. However, we rule this model out using two topological tests: the S-statistic and the matched circle test.

We perform an independent foreground analysis of the WMAP maps to produce a cleaned CMB map (available online) useful for cross-correlation with, e.g., galaxy and X-ray maps. We use a variant of the Tegmark & Efstathiou (1996) technique that is completely blind, making no assumptions about the CMB power spectrum, the foregrounds, WMAP detector noise or external templates. Compared with the foreground-cleaned internal linear combination map produced by the WMAP team, our map has the advantage of containing less non-CMB power (from foregrounds and detector noise) outside the Galactic plane. The difference is most important on the the angular scale of the first acoustic peak and below, since our cleaned map is at the highest (13') rather than lowest (49') WMAP resolution. We also produce a Wiener filtered CMB map, representing our best guess as to what the CMB sky actually looks like, as well as CMB-free maps at the five WMAP frequencies useful for foreground studies. We argue that our CMB map is clean enough that the lowest multipoles can be measured without any galaxy cut, and obtain a quadrupole value that is slightly less low than that from the cut-sky WMAP team analysis. This can be understood from a map of the CMB quadrupole, which shows much of its power falling within the Galaxy cut region, seemingly coincidentally. Intriguingly, both the quadrupole and the octopole are seen to have power suppressed along a particular spatial axis, which lines up between the two, roughly towards (l,b) \~ (-110,60) in Virgo.

This paper presents a scheme to deal accurately and efficiently with complex angular masks, such as occur typically in galaxy surveys. An angular mask is taken to be an arbitrary union of arbitrarily weighted angular regions bounded by arbitrary numbers of edges. The restrictions on the mask are (i) that each edge must be part of some circle on the sphere (but not necessarily a great circle), and (ii) that the weight within each subregion of the mask must be constant. The scheme works by resolving a mask into disjoint polygons, convex angular regions bounded by arbitrary numbers of edges. The polygons may be regarded as the `pixels' of a mask, with the feature that the pixels are allowed to take a rather general shape, rather than following some predefined regular pattern. Among other things, the scheme includes facilities to compute the spherical harmonics of the angular mask, and Data-Random <DR> and Random-Random <RR> angular integrals. A software package mangle which implements this scheme, along with complete software documentation, is available at http://casa.colorado.edu/~ajsh/mangle/ .

By comparing the PSCz galaxy power spectrum with haloes from nested and phased N-body simulations, we try to understand how IRAS infrared-selected galaxies populate dark-matter haloes. We pay special attention to the way we identify haloes in the simulations.

We report on the linear matter power spectrum reconstructed from the Lyman alpha forest, on baryonic wiggles in the galaxy power spectrum, and on how to measure real space (as opposed to redshift space) power.

We compute the real-space power spectrum and the redshift-space distortions of galaxies in the 2dF 100k galaxy redshift survey using pseudo-Karhunen-Loeve eigenmodes and the stochastic bias formalism. Our results agree well with those published by the 2dFGRS team, and have the added advantage of producing easy-to-interpret uncorrelated minimum-variance measurements of the galaxy-galaxy, galaxy-velocity and velocity-velocity power spectra in 27 k-bands, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.8h/Mpc. We find no significant detection of baryonic wiggles, although our results are consistent with a standard flat Omega_Lambda=0.7 ``concordance'' model and previous tantalizing hints of baryonic oscillations. We measure the galaxy-matter correlation coefficient r > 0.4 and the redshift-distortion parameter beta=0.49+/-0.16 for r=1 (beta=0.47+/- 0.16 without finger-of-god compression). Since this is an apparent-magnitude limited sample, luminosity-dependent bias may cause a slight red-tilt in the power spectum. A battery of systematic error tests indicate that the survey is not only impressive in size, but also unusually clean, free of systematic errors at the level to which our tests are sensitive. Our measurements and window functions are available at http://www.hep.upenn.edu/~max/2df.html together with the survey mask, radial selection function and uniform subsample of the survey that we have constructed.

We investigate possible systematic errors in the recent measurement of the matter power spectrum from the Lyman-alpha forest by Croft et al. (2001). We find that for a large set of prior cosmological models the Croft et al. result holds quite well, with systematic errors being comparable to random ones, when a dependence of the recovered matter power spectrum on the cosmological parameters at z~3 is taken into account. We find that peculiar velocities cause the flux power spectrum to be smoothed over about 100-300 km/s, dependng on scale. Consequently, the recovered matter power spectrum is a smoothed version of the underlying true power spectrum. Uncertainties in the recovered power spectrum are thus correlated over about 100-300 km/s. As a side effect, we find that residual fluctuations in the ionizing background, while having almost no effect on the recovered matter power spectrum, significantly bias estimates of the baryon density from the Lyman-alpha forest data. We therefore conclude that the Croft et al. result provides a powerful new constraint on cosmological parameters and models of structure formation.

We report a measurement of the real space (not redshift space) power spectrum of galaxies over four and a half decades of wavenumber, 0.01 to 300 h/Mpc, from the IRAS Point Source Catalog Redshift Survey (PSCz). Since estimates of power are highly correlated in the nonlinear regime, we also report results for the prewhitened power spectrum, which is less correlated. The inferred bias between optically-selected APM and IRAS-selected PSCz galaxies is about 1.15 at linear scales < 0.3 h/Mpc, increasing to about 1.4 at nonlinear scales > 1 h/Mpc. The nonlinear power spectrum of PSCz shows a near power-law behaviour to the smallest scales measured, with mild upward curvature in the broad vicinity of 1.5 h/Mpc. Contrary to the prediction of unbiased Dark Matter models, there is no prominent inflection at the linear-nonlinear transition scale, and no turnover at the transition to the stable clustering regime. The nonlinear power spectrum of PSCz requires scale-dependent bias: all Dark Matter models without scale-dependent bias are ruled out with high confidence.

We constrain the spectrum of primordial curvature perturbations by using recent Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) data. Specifically, we consider CMB data from the COBE, Boomerang and Maxima experiments, the real space galaxy power spectrum from the IRAS PSCz survey, and the linear matter power spectrum inferred from Ly-alpha forest spectra. We study the case of single field slow roll inflationary models, and we extract bounds on the scalar spectral index, n, the tensor to scalar ratio, r, and the running of the scalar spectral index, dlnk, for various combinations of the observational data. We find that CMB data, when combined with data from Lyman-alpha forest, place strong constraints on the inflationary parameters. Specifically, we obtain n \approx 0.9, r < 0.3 and dlnk \approx 0, indicating that single field hybrid models are ruled out.

The combined CfA2/SSRS redshift catalog has recently been improved and made public. We compute the redshift-space power spectrum of this catalog using the regions at -2.5 < delta < 50 north (8h < alpha < 17h) and south (20h < alpha < 3h, b < -13) of the Galactic plane, where it is 98% complete down to Zwicky magnitude 15.5 and contains 13,681 galaxies. Our analysis uses Heavens-Taylor mode expansion, Karhunen-Loeve data compression and the Fisher matrix technique to compute quadratic band power estimates. This allows an exact calculation of window functions, including the integral constraint, in addition to the production of a power spectrum with uncorrelated error bars. Our results with this larger data set agree well with previous studies. We analyze 101Mpc/h and 130Mpc/h volume-limited subsets in addition to the full magnitude-limited sample, and our results are well fit by, e.g., Lambda-CDM models with bias b=1.2, 1.4 and 1.4, respectively. We estimate the effect of extinction using the Schlegel, Finkbeiner & Davis dust map. Our results are exclusively for the redshift space galaxy power spectrum, so they can only be compared with theoretical predictions if appropriate corrections are made for biasing and redshift space distortions.

Formulae are presented for the linear growth factor D/a and its logarithmic derivative dlnD/dlna in expanding Friedmann-Robertson-Walker Universes with arbitrary matter and vacuum densities. The formulae permit rapid and stable numerical evaluation. A fortran program is available at http://casa.colorado.edu/~ajsh/growl/ .

We present a method for calculating large numbers of power spectra C_l and P(k) that accelerates CMBfast by a factor around 10^3 without appreciable loss of accuracy, then apply it to constrain 11 cosmological parameters from current Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) data. While the CMB alone still suffers from several degeneracies, allowing, e.g., closed models with strong tilt and tensor contributions, the shape of the real space power spectrum of galaxies from the IRAS Point Source Catalogue Redshift (PSCz) survey breaks these degeneracies and helps place strong constraints on most parameters. At 95% confidence, the combined CMB and LSS data imply a baryon density 0.020 < omega_b < 0.037, dark matter density 0.10 < omega_d < 0.32 with a neutrino fraction f_nu < 38%, vacuum density Omega_Lambda < 0.76, curvature -0.19 < Omega_k < 0.10, scalar tilt 0.86 < n_s < 1.16 and reionization optical depth tau < 0.44. These joint constraints are quite robust, changing little when we impose priors on the Hubble parameter, tilt, flatness, gravity waves or reionization. Adding nucleosynthesis and neutrino priors on the other hand tightens parameters considerably, requiring Omega_Lambda > 0.49 and a red-tilt. The analysis allows a number of consistency tests to be made, all of which pass. At the 95% level, the flat scalar ``concordance model'' with Omega_Lambda=0.62, omega_d=0.13, omega_b=0.02, f_nu~0, n_s=0.9, tau=0.1, h=0.63 is consistent with external constraints. The inferred PSCz bias b~1.2 agrees with the value estimated independently from redshift space distortions. The inferred cosmological constant value agrees with the one independently derived from SN 1a studies. Cosmology seems to be on the right track.

We present a state-of-the-art linear redshift distortion analysis of the recently published IRAS Point Source Catalog Redshift Survey (PSCz). The procedure involves linear compression into 4096 Karhunen-Loeve modes culled from a potential pool of about 3 x 10^5 modes, followed by quadratic compression into three separate power spectra, the galaxy-galaxy, galaxy-velocity, and velocity-velocity power spectra. Least squares fitting to the decorrelated power spectra yields a linear redshift distortion parameter beta = Omega_m^0.6/b = 0.41(+0.13,-0.12).

As experimentalists step up their pursuit of cold dark matter particles and neutrino masses, cosmological constraints are tightening. We compute the joint constraints on 11 cosmological parameters from cosmic microwave background and large scale structure data, and find that at 95% confidence, the total (cold+hot) dark matter density is h^2 Omega_dark = 0.20 +.12 -.10 with at most 38% of it being hot (due to neutrinos). A few assumptions, including negligible neutrinos, tighten this measurement to h^2 Omega_dark = 0.13 +0.04 -0.02, i.e, 2.4e-27 kg/m^3 give or take 20%.

Nonlinear evolution causes the galaxy power spectrum to become broadly correlated over different wavenumbers. It is shown that prewhitening the power spectrum - transforming the power spectrum in such a way that the noise covariance becomes proportional to the unit matrix - greatly narrows the covariance of power. The eigenfunctions of the covariance of the prewhitened nonlinear power spectrum provide a set of almost uncorrelated nonlinear modes somewhat analogous to the Fourier modes of the power spectrum itself in the linear, Gaussian regime. These almost uncorrelated modes make it possible to construct a near minimum variance estimator and Fisher matrix of the prewhitened nonlinear power spectrum analogous to the Feldman-Kaiser-Peacock estimator of the linear power spectrum. The paper concludes with summary recipes, in gourmet, fine, and fastfood versions, of how to measure the prewhitened nonlinear power spectrum from a galaxy survey in the FKP approximation. An Appendix presents FFTLog, a code for taking the fast Fourier or Hankel transform of a periodic sequence of logarithmically spaced points, which proves useful in some of the manipulations.

It is shown how to decorrelate the (prewhitened) power spectrum measured from a galaxy survey into a set of high resolution uncorrelated band-powers. The treatment includes nonlinearity, but not redshift distortions. Amongst the infinitely many possible decorrelation matrices, the square root of the Fisher matrix, or a scaled version thereof, offers a particularly good choice, in the sense that the band-power windows are narrow, approximately symmetric, and well-behaved in the presence of noise. We use this method to compute band-power windows for, and the information content of, the Sloan Digital Sky Survey, the Las Campanas Redshift Survey, and the IRAS 1.2 Jy Survey.

Ultraviolet imaging of the remnant of Supernova 1885 in M31 with the Hubble Space Telescope using the F255W filter on the WFPC2 reveals a dark spot of Fe II absorption at the remnant's known position in the bulge of M31. The diameter of the absorbing spot is 0"55 +- 0"15, slightly smaller than, but consistent with, the 0"70 +- 0"05 diameter measured in the higher quality WFPC2 Ca II absorption image previously reported by us. The measured ratio of flux inside to outside SNR 1885 in the Fe II image is 0.24 +- 0.17, consistent with the ratio 0.33 +- 0.04 expected on the basis of a model fit to the previously obtained near-UV FOS spectrum. The observed depth of Fe II absorption suggests that Fe II is fully saturated, with an iron mass in the range M_Fe = 0.1-1.0 Msun. Besides Fe, ion species Mg I, Mg II, and Mn I probably make some contribution to the absorption from the SN 1885 remnant in the F255W image.

Near UV HST images of the remnant of SN 1885 (S And) in M31 show a 0"70 +- 0"05 diameter absorption disk silhouetted against M31's central bulge, at SN 1885's historically reported position. The disk's size corresponds to a linear diameter of 2.5 +- 0.4 pc at a distance of 725 +- 70 kpc, implying an average expansion velocity of 11000 +- 2000 km/s over 110 years. Low-dispersion FOS spectra over 3200-4800 A; reveal that the absorption arises principally from Ca II H & K (equivalent width ~215 A;) with weaker absorption features of Ca I 4227 A; and Fe I 3720 A;. The flux at Ca II line center indicates a foreground starlight fraction of 0.21, which places SNR 1885 some 64 pc to the near side of the midpoint of the M31 bulge, comparable to its projected 55 pc distance from the nucleus. The absorption line profiles suggest an approximately spherically symmetric, bell-shaped density distribution of supernova ejecta freely expanding at up to 13100 +- 1500 km/s. We estimate Ca I, Ca II, and Fe I masses of 2.9(+2.4,-0.6) x 10^-4 M_o, 0.005(+0.016,-0.002) M_o, and 0.013(+0.010,-0.005) M_o respectively. If the ionization state of iron is similar to the observed ionization state of calcium, M_CaII/M_CaI = 16(+42,-5), then the mass of Fe II is 0.21(+0.74,-0.08) M_o, consistent with that expected for either normal or subluminous SN Ia.

Precision measurements of the galaxy power spectrum P(k) require a data analysis pipeline that is both fast enough to be computationally feasible and accurate enough to take full advantage of high-quality data. We present a rigorous discussion of different methods of power spectrum estimation, with emphasis on the traditional Fourier method, the linear (Karhunen-Loeve; KL), and quadratic data compression schemes, showing in what approximations they give the same result. To improve speed, we show how many of the advantages of KL data compression and power spectrum estimation may be achieved with a computationally faster quadratic method. To improve accuracy, we derive analytic expressions for handling the integral constraint, since it is crucial that finite volume effects are accurately corrected for on scales comparable to the depth of the survey. We also show that for the KL and quadratic techniques, multiple constraints can be included via simple matrix operations, thereby rendering the results less sensitive to galactic extinction and mis-estimates of the radial selection function. We present a data analysis pipeline that we argue does justice to the increases in both quality and quantity of data that upcoming redshift surveys will provide. It uses three analysis techniques in conjunction: a traditional Fourier approach on small scales, a pixelized quadratic matrix method on large scales and a pixelized KL eigenmode analysis to probe anisotropic effects such as redshift-space distortions.

Peculiar velocities induce apparent line of sight displacements of galaxies in redshift space, distorting the pattern of clustering in the radial versus transverse directions. On large scales, the amplitude of the distortion yields a measure of the dimensionless linear growth rate $\beta \approx \Omega^{0.6}/b$, where $\Omega$ is the cosmological density and $b$ the linear bias factor. To make the maximum statistical use of the data in a wide angle redshift survey, and for the greatest accuracy, the spherical character of the distortion needs to be treated properly, rather than in the simpler plane parallel approximation. In the linear regime, the redshift space correlation function is described by a spherical distortion operator acting on the true correlation function. It is pointed out here that there exists an operator, which is essentially the logarithmic derivative with respect to pair separation, which both commutes with the spherical distortion operator, and at the same time defines a characteristic scale of separation. The correlation function can be expanded in eigenfunctions of this operator, and these eigenfunctions are eigenfunctions of the distortion operator. Ratios of the observed amplitudes of the eigenfunctions yield measures of the linear growth rate $\beta$ in a manner independent of the shape of the correlation function. More generally, the logarithmic derivative $\partial/\partial\ln r$ with respect to depth $r$, along with the square $L^2$ and component $L_z$ of the angular momentum operator, form a complete set of commuting operators for the spherical distortion operator acting on the density. The eigenfunctions of this complete set of operators are spherical waves about the observer, with radial part lying in logarithmic real or Fourier space.

Redshift space distortions on large scales can be used to measure the linear growth rate parameter $\ff \approx \Omega^{0.6}/b$. I report here measurements of such distortions in the IRAS 2 Jy, 1.2 Jy, and QDOT redshift surveys, finding $\ff = 0.69^{+ .21}_{- .19}$ from a merged QDOT plus 1.2 Jy catalogue. Unfortunately, confidence in this result is undermined by a marked ($4\sigma$) change in the pattern of clustering in QDOT beyond about $80 h^{-1} Mpc$. A similar effect may be present at a mild level in the 1.2 Jy survey. The effect may be caused by systematic variation in the effective flux limit of the IRAS PSC over the sky, with a dispersion of $\sim 0.1$ Jy on scales $\sim 7^{\circ}$. If so, then the value of $\ff$ inferred from redshift distortions in IRAS surveys may be systematically underestimated.

Redshift maps of galaxies in the Universe are distorted by the peculiar velocities of galaxies along the line of sight. The amplitude of the distortions on large, linear scales yields a measurement of the linear redshift distortion parameter, which is $\beta \approx \Omega_0^{0.6}/b$ in standard cosmology with cosmological density $\Omega_0$ and light-to-mass bias $b$. All measurements of $\beta$ from linear redshift distortions published up to mid 1997 are reviewed. The average and standard deviation of the reported values is $\beta_{optical} = 0.52 \pm 0.26$ for optically selected galaxies, and $\beta_{IRAS} = 0.77 \pm 0.22$ for IRAS selected galaxies. The implied relative bias is $b_{optical}/b_{IRAS} \approx 1.5$. If optical galaxies are unbiased, then $\Omega_0 = 0.33^{+0.32}_{-0.22}$, while if IRAS galaxies are unbiased, then $\Omega_0 = 0.63^{+0.35}_{-0.27}$.

We describe how to compute estimates of the power spectrum C_l from Cosmic Microwave Background (CMB) maps that not only retain all the cosmological information, but also have uncorrelated error bars and well-behaved window functions. We apply this technique to the 4-year COBE/DMR data, and tabulate eight uncorrelated band-powers.

We have obtained a far-ultraviolet spectrum (1150 - 1600 Ang.) of a hot subdwarf star behind the remnant of SN 1006 with the Faint Object Spectrograph (FOS) on the Hubble Space Telescope. The high-quality spectrum is used to test previous identifications of the strong absorption features discovered with the International Ultraviolet Explorer. These features have FWHM = 4000 (+/- 300) km/sec and are not at the rest wavelengths of known interstellar lines, as opposed to the broader (8000 km/sec FWHM) Fe II lines from the remnant centered at zero km/sec in near-UV FOS spectra. We confirm that the broad absorption features are principally due to redshifted Si II, Si III, and Si IV lines, which are centered at a radial velocity of 5100 (+/- 200) km/sec. The Si II 1260.4 profile is asymmetric, with a nearly flat core and sharp red wing, unlike the Si II 1526.7 and Si IV 1393.8, 1402.8 profiles. One possible explanation is additional absorption from another species. Previous work has suggested that S II 1250.6, 1253.8, 1259.5 at a radial velocity of 6000 km/sec is responsible, but this would require a sulfur to silicon abundance ratio that is at least a factor of ten higher than expected. Another possible explanation is that the Si II and Si IV profiles are intrinsically different, but this does not explain the symmetric (albeit weaker) Si II 1526.7 profile.

This is the second of two papers which address the problem of measuring the unredshifted power spectrum of fluctuations from a galaxy survey in optimal fashion. A key quantity is the Fisher matrix, which is the inverse of the covariance matrix of minimum variance estimators of the power spectrum of the survey. It is shown that bases of kernels which give rise to complete sets of statistically orthogonal windowed power spectra are obtained in general from the eigenfunctions of the Fisher matrix scaled by some arbitrary positive definite scaling matrix. Among the many possible bases of kernels, there is a basis, obtained by applying an infinitely steep scaling function, which leads to kernels which are positive and compact in Fourier space. This basis of kernels, along with the associated minimum variance pair weighting derived in the previous paper, would appear to offer a solution to the problem of how to measure the unredshifted power spectrum optimally. Illustrative kernels are presented for the case of the PSCz survey.

This is the first of a pair of papers which address the problem of measuring the unredshifted power spectrum in optimal fashion from a survey of galaxies, with arbitrary geometry, for Gaussian or non-Gaussian fluctuations, in real or redshift space. In this first paper, that pair weighting is derived which formally minimizes the expected variance of the unredshifted power spectrum windowed over some arbitrary kernel. The inverse of the covariance matrix of minimum variance estimators of windowed power spectra is the Fisher information matrix, which plays a central role in establishing optimal estimators. Actually computing the minimum variance pair window and the Fisher matrix in a real survey still presents a formidable numerical problem, so here a perturbation series solution is developed. The properties of the Fisher matrix evaluated according to the approximate method suggested here are investigated in more detail in the second paper.

We present a theoretical interpretation of the broad silicon and iron UV absorption features observed with the Hubble Space Telescope in the spectrum of the Schweizer-Middleditch star behind the remnant of Supernova 1006. These features are caused by supernova ejecta in SN1006. We propose that the redshifted SiII2 1260 A feature consists of both unshocked and shocked SiII. The sharp red edge of the line at 7070 km/s indicates the position of the reverse shock, while its Gaussian blue edge reveals shocked Si with a mean velocity of 5050 km/s and a dispersion of 1240 km/s, implying a reverse shock velocity of 2860 km/s. The measured velocities satisfy the energy jump condition for a strong shock, provided that all the shock energy goes into ions, with little or no collisionless heating of electrons. The line profiles of the SiIII and SiIV absorption features indicate that they arise mostly from shocked Si. The total mass of shocked and unshocked Si inferred from the SiII, SiIII and SiIV profiles is M_Si = 0.25 \pm 0.01 Msun on the assumption of spherical symmetry. Unshocked Si extends upwards from 5600 km/s. Although there appears to be some Fe mixed with the Si at lower velocities < 7070 km/s, the absence of FeII absorption with the same profile as the shocked SiII suggests little Fe mixed with Si at higher (before being shocked) velocities. The column density of shocked SiII is close to that expected for SiII undergoing steady state collisional ionization behind the reverse shock, provided that the electron to SiII ratio is low, from which we infer that most of the shocked Si is likely to be of a fairly high degree of purity, unmixed with other elements. We propose that the ambient interstellar density on the far side of SN1006 is anomalously low compared to the density around the rest of the remnant. This

We present an expression for the nonlinear evolution of the cosmological power spectrum based on following Lagrangian trajectories. This is simplified using the Zel'dovich approximation to trace particle displacements, assuming Gaussian initial conditions. The model is found to exhibit the transfer of power from large to small scales expected in self- gravitating fields. We have extended this analysis into redshift-space and found a solution for the nonlinear, anisotropic redshift-space power spectrum in the limit of plane--parallel redshift distortions. The quadrupole-to- monopole ratio is calculated for the case of power-law initial spectra. We find that the shape of this ratio depends on the shape of the initial spectrum, but when scaled to linear theory depends only weakly on the redshift-space distortion parameter, $\beta$. The point of zero-crossing of the quadrupole, $k_0$, is found to obey a scaling relation. This model is found to be in agreement with $N$-body simulations on scales down to the zero-crossing of the quadrupole, although the wavenumber at zero-crossing is underestimated. These results are applied to the quadrupole--monopole ratio found in the merged QDOT+1.2 Jy IRAS redshift survey. We have estimated that the distortion parameter is constrained to be $\beta>0.5$ at the $95 \%$ level. The local primordial spectral slope is not well constrained, but analysis suggests $n \approx -2$ in the translinear regime. The zero-crossing scale of the quadrupole is $k_0=0.5 \pm 0.1 h/Mpc$ and from this we infer the amplitude of clustering is $\sigma_8=0.7 \pm 0.05$. We suggest that the success of this model is due to nonlinear redshift--space effects arising from infall onto caustics and is not dominated by virialised cluster cores.

We derive the evolution equations for the electric and magnetic parts of the Weyl tensor for cold dust from both general relativity and Newtonian gravity. In a locally inertial frame at rest in the fluid frame, the Newtonian equations agree with those of general relativity. We give explicit expressions for the electric and magnetic parts of the Weyl tensor in the Newtonian limit. In general, the magnetic part does not vanish, implying that the Lagrangian evolution of the fluid is not purely local.