We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

We present constrained radiative transfer calculations of Lyman Alpha (Lya) photons propagating through clumpy, dusty, large scale outflows, and explore whether we can quantitatively explain the Lya halos that have been observed around Lyman Break Galaxies. We construct phenomenological models of large-scale outflows which consist of cold clumps that are in pressure equilibrium with a constant-velocity hot wind. First we consider models in which the cold clumps are distributed symmetrically around the galaxy, and in which the clumps undergo a continuous acceleration in its 'circumgalactic' medium (CGM). We constrain the properties of the cold clumps (radius, velocity, HI column density, & number density) by matching the observed Lya absorption strength of the CGM in the spectra of background galaxies. We then insert a Lya source in the center of this clumpy outflow, which consists of 1e5-1e6 clumps, and compute observable properties of the scattered Lya photons. In these models, the scattered radiation forms halos that are significantly more concentrated than is observed. In order to simultaneously reproduce the observed Lya absorption line strengths and the Lya halos, we require - preferably bipolar - outflows in which the clumps decelerate after their initial acceleration. This deceleration is predicted naturally in 'momentum-driven' wind models of clumpy outflows. In models that simultaneously fit the absorption and emission line data, the predicted linear polarization is ~30-40% at a surface brightness contour of S=1e-18 erg/s/cm^2/arcsec^2. Our work illustrates clearly that Lya emission line halos around starforming galaxies provide valuable constraints on the cold gas distribution & kinematics in their circumgalactic medium, and that these constraints complement those obtained from absorption line studies alone.

In the near future, measurements of metal absorption features in the intergalactic medium (IGM) will become an important constraint on models of the formation and evolution of the earliest galaxies, the properties of the first stars, and the reionization and enrichment of the IGM. The first measurement of a metal abundance in the IGM at a redshift approaching the epoch of reionization already offers intriguing hints. Between z=5.8 and 4.7 (a 0.3 Gyr interval only 1 Gyr after the big bang), the measured density of CIV absorbers in the IGM increased by a factor of ~ 3.5 (Ryan-Weber et al. 2009; Becker, Rauch & Sargent 2009). If these values prove to be accurate, they pose two puzzles: (1) The total amount of CIV at z=5.8 implies too little star formation to reionize the IGM by z=6 or to match the WMAP electron scattering optical depth (tau). (2) The rapid growth from z = 6-5 is faster than the buildup of stellar mass or the increase in the star formation rate density over the same interval. We show that a delay of ~ 0.4-0.7 Gyr between the instantaneous production of ionizing photons and the later production of metal absorption features (added to the delay due to stellar lifetimes) can provide the full explanation for both puzzles. We calculate the delay in metal production due to finite stellar lifetimes alone and find that it is too short to explain the rapid CIV density increase. The additional delay could naturally be explained as the result of ~ 200 km/s outflows carrying carbon to distances of ~ 100 kpc, the typical distance between galaxies and CIV absorbers in enrichment simulations, and the typical outflow or absorption region scale observed at z ~ 2-3.

Arguably the best hope of understanding the tail end of the reionization of the intergalactic medium (IGM) at redshift z > 6 is through the detection and characterization of the Gunn-Peterson (GP) damping wing absorption of the IGM in bright quasar spectra. However, the use of quasar spectra to measure the IGM damping wing requires a model of the quasar's intrinsic Lyman-alpha emission line. Here we quantify the uncertainties in the intrinsic line shapes, and how those uncertainties affect the determination of the IGM neutral fraction. We have assembled a catalog of high-resolution HST spectra of the emission lines of unobscured low-redshift quasars, and have characterized the variance in the shapes of their lines. We then add simulated absorption from the high-redshift IGM to these quasar spectra in order to determine the corresponding uncertainties in reionization constraints using current and future samples of z > 6 quasar spectra. We find that, if the redshift of the Lyman-alpha emission line is presumed to coincide with the systemic redshift determined from metal lines, the inferred IGM neutral fraction is systematically biased to low values due to a systematic blueshift of the Lyman-alpha line relative to the metal lines. If a similar blueshift persists in quasars at z > 6, this bias strengthens previous claims of a significant neutral hydrogen fraction at z ~ 6. This technique is capable of making a robust distinction between a highly ionized (x_IGM ~ 10^-3) and a neutral (x_IGM = 1) IGM with even a few bright quasars.

High-redshift quasars (z >~ 6) drive ionization fronts into the intergalactic medium (IGM). If the thickness of the front can be measured, it can provide a novel constraint on the ionizing spectral energy distribution (SED). Here we follow the propagation of an I-front into a uniform IGM, and compute its thickness for a range of possible quasar spectra and ages. We also explore the effects of uniform and non-uniform ionizing backgrounds. We find that even for hard spectra, the fronts are initially thin, with a thickness much smaller than the mean free path of ionizing photons, but the thickness increases as the front approaches equilibrium in 10^8 - 10^9 years, and can eventually significantly exceed simple estimates based on the mean free path. With a high intrinsic hydrogen column density obscuring the source (log(N_H/cm^-2) >~ 19.2) or a hard power-law spectrum combined with some obscuration (e.g. dlog(F_\nu)/dlog(\nu) >~ -1.2 at log(N_H/cm^-2) >~ 18.0), the thickness of the front exceeds ~1 physical Mpc and may be measurable from the morphology of its redshifted 21cm signal. We find that the highly ionized inner part of the front, which may be probed by Lyman line absorption spectra, remains sharp for bright quasars unless a large obscuring column (log(N_H/cm^-2) >~ 19.2) removes most of their ionizing photons up to ~40 eV. For obscured sources with log(N_H/cm^-2) >~ 19.8, embedded in a significantly neutral IGM, the black Lyman-alpha trough (where the neutral fraction is ~10^-3) underestimates the size of the HII region by a factor of >~4.

The reionization history of the intergalactic medium (IGM) at high redshift (z > 6) was likely strongly shaped by several global feedback processes. Because the earliest ionizing sources formed at the locations of the rare density peaks, their spatial distribution was strongly clustered. Here we demonstrate that this clustering significantly boosts the impact of feedback processes operating at high redshift. We build a semi-analytical model to include feedback and clustering simultaneously, and apply this model to the suppression of star-formation in minihalos due to photoionization. The model is built on the excursion-set-based formalism of Furlanetto, Zaldarriaga and Hernquist (2004), which incorporates the clustering of ionizing sources, and which we here extend to include suppression of star formation in minihalos. We find that clustering increases the mean HII bubble size by a factor of several, and it dramatically increases the fraction of minihalos that are suppressed, by a factor of up to 60 relative to a randomly distributed population. This enhanced suppression can significantly reduce the electron scattering optical depth, as required by the three-year data from the Wilkinson Microwave Anisotropy Probe (WMAP). We argue that source clustering is likely to similarly boost the importance of a variety of other feedback mechanisms.

X-ray line profile analysis has proved to be the most direct diagnostic of the kinematics and spatial distribution of the very hot plasma around O stars. In this paper we apply several analysis techniques to the emission lines in the Chandra HETGS spectrum of the late-O supergiant zeta Ori (O9.7 Ib), including the fitting of a simple line-profile model. We show that there is distinct evidence for blue shifts and profile asymmetry, as well as broadening in the X-ray emission lines of zeta Ori. These are the observational hallmarks of a wind-shock X-ray source, and the results for zeta Ori are very similar to those for the earlier O star, zeta Pup, which we have previously shown to be well-fit by the same wind-shock line-profile model. The more subtle effects on the line-profile morphologies in zeta Ori, as compared to zeta Pup, are consistent with the somewhat lower density wind in this later O supergiant. In both stars, the wind optical depths required to explain the mildly asymmetric X-ray line profiles imply reductions in the effective opacity of nearly an order of magnitude, which may be explained by some combination of mass-loss rate reduction and large-scale clumping, with its associated porosity-based effects on radiation transfer. In the context of the recent reanalysis of the helium-like line intensity ratios in both zeta Ori and zeta Pup, and also in light of recent work questioning the published mass-loss rates in OB stars, these new results indicate that the X-ray emission from zeta Ori can be understood within the framework of the standard wind-shock scenario for hot stars.

X-ray emission line profiles provide the most direct insight into the dynamics and spatial distribution of the hot, X-ray-emitting plasma above the surfaces of OB stars. The O supergiant zeta Puppis shows broad, blueshifted, and asymmetric line profiles, generally consistent with the wind-shock picture. We model the profiles of eight lines in the Chandra HETGS spectrum of this prototypical hot star. The fitted lines indicate that the plasma is distributed throughout the wind starting close to the photosphere, that there is significantly less attenuation of the X-rays by the overlying wind than is generally supposed, and that there is not a strong trend in wind absorption with wavelength.

The launch of high-spectral-resolution x-ray telescopes (Chandra, XMM) has provided a host of new spectral line diagnostics for the astrophysics community. In this paper we discuss Doppler-broadened emission line profiles from highly supersonic outflows of massive stars. These outflows, or winds, are driven by radiation pressure and carry a tremendous amount of kinetic energy, which can be converted to x rays by shock-heating even a small fraction of the wind plasma. The unshocked, cold wind is a source of continuum opacity to the x rays generated in the shock-heated portion of the wind. Thus the emergent line profiles are affected by transport through a two-component, moving, optically thick medium. While complicated, the interactions among these physical effects can provide quantitative information about the spatial distribution and velocity of the x-ray-emitting and absorbing plasma in stellar winds. We present quantitative models of both a spherically-symmetric wind and a wind with hot plasma confined in an equatorial disk by a dipole magnetic field.