The two-point correlation function of randomly distributed Lyman alpha clouds

It is often assumed that Ly-alpha forest clouds are randomly distributed, intergalactic objects that are highly ionized by the UV background produced by quasars. If these assumptions are true, fluctuations in the UV background should produce a nonzero two point correlation function in the Ly-alpha forest. This effect, which is really just a generalization of the proximity effect, is more significant at high redshift (z is approximately 3-4) because the mean free path for UV photons is smaller there, and the fluctuations correspondingly larger. This effect was studied using both the semi-analytic techniques of Zuo's recent papers and Monte Carlo simulations. The correlation function is expected to have a small yet potentially measurable amplitude that is consistent with current upper limits. Furthermore, the signature of this effect is distinctive because the nonzero correlation function extends over the photon mean free path, which is larger than the expected scale of large-scale structure. Observations or upper limits on this effect could provide information about the source of the ionizing background at high redshifts and the nature of the Ly-alpha forest clouds

Sub-mm Galaxies in Cosmological Simulations

We study the predicted sub-mm emission from massive galaxies in a Lambda-CDM universe, using hydrodynamic cosmological simulations. Assuming that most of the emission from newly formed stars is absorbed and reradiated in the rest-frame far-IR, we calculate the number of galaxies that would be detected in sub-mm surveys conducted with SCUBA. The predicted number counts are strongly dependent on the assumed dust temperature and emissivity law. With plausible choices for SED parameters (e.g., T=35 K, beta=1.0), the simulation predictions reproduce the observed number counts above ~ 1 mJy. The sources have a broad redshift distribution with median z ~ 2, in reasonable agreement with observational constraints. However, the predicted count distribution may be too steep at the faint end, and the fraction of low redshift objects may be larger than observed. In this physical model of the sub-mm galaxy population, the objects detected in existing surveys consist mainly of massive galaxies (several M_*) forming stars fairly steadily over timescales ~ 10^8-10^9 years, at moderate rates ~100 Msun/yr. The typical descendants of these sub-mm sources are even more massive galaxies, with old stellar populations, found primarily in dense environments. While the resolution of our simulations is not sufficient to determine galaxy morphologies, these properties support the proposed identification of sub-mm sources with massive ellipticals in the process of formation. The most robust and distinctive prediction of this model, stemming directly from the long timescale and correspondingly moderate rate of star formation, is that the far-IR SEDs of SCUBA sources have a relative high 850 micron luminosity for a given bolometric luminosity. [Abridged]Comment: Submitted to ApJ. 34 pages including 8 PS figure

Galactic outflows and the kinematics of damped Lyman alpha absorbers

The kinematics of damped Lyman alpha absorbers (DLAs) are difficult to reproduce in hierarchical galaxy formation models, particularly the preponderance of wide systems. We investigate DLA kinematics at z=3 using high-resolution cosmological hydrodynamical simulations that include a heuristic model for galactic outflows. Without outflows, our simulations fail to yield enough wide DLAs, as in previous studies. With outflows, predicted DLA kinematics are in much better agreement with observations. Comparing two outflow models, we find that a model based on momentum-driven wind scalings provides the best match to the observed DLA kinematic statistics of Prochaska & Wolfe. In this model, DLAs typically arise a few kpc away from galaxies that would be identified in emission. Narrow DLAs can arise from any halo and galaxy mass, but wide ones only arise in halos with mass >10^11 Mo, from either large central or small satellite galaxies. This implies that the success of this outflow model originates from being most efficient at pushing gas out from small satellite galaxies living in larger halos. This increases the cross-section for large halos relative to smaller ones, thereby yielding wider kinematics. Our simulations do not include radiative transfer effects or detailed metal tracking, and outflows are modeled heuristically, but they strongly suggest that galactic outflows are central to understanding DLA kinematics. An interesting consequence is that DLA kinematics may place constraints on the nature and efficiency of gas ejection from high-z galaxies.Comment: submitted to MNRA

Generation of mock tidal streams

In this paper we discuss a method for the generation of mock tidal streams. Using an ensemble of simulations in an isochrone potential where the actions and frequencies are known, we derive an empirical recipe for the evolving satellite mass and the corresponding mass loss rate, and the ejection conditions of the stream material. The resulting stream can then be quickly generated either with direct orbital integration, or by using the action-angle formalism. The model naturally produces streaky features within the stream. These are formed due to the radial oscillation of the progenitor and the bursts of stars emitted near pericenter, rather than clumping at particular oscillation phases as sometimes suggested. When detectable, these streaky features are a reliable diagnostic for the stream's direction of motion and encode other information on the progenitor and its orbit. We show several tests of the recipe in alternate potentials, including a case with a chaotic progenitor orbit which displays a marked effect on the width of the stream. Although the specific ejection recipe may need adjusting when elements such as the orbit or satellite density profile are changed significantly, our examples suggest that model tidal streams can be quickly and accurately generated by models of this general type for use in Bayesian sampling.Comment: 20 pages, 11 figures, 1 table; submitted to MNRA

Intergalactic Helium Absorption toward High-Redshift Quasars

The recent Hubble Space Telescope (HST) observations of the z(q) = 3.286 quasar Q0302-003 (Jakobsen et at. 1994) and the z(q) = 3.185 quasar Q1935-67 by Tytler (1995) show absorption edges at the redshifted wavelength of He II 304 A. A key goal is to distinguish between contributions from discrete Ly-alpha forest clouds and a smoothly distributed intergalactic medium (IGM). We model the contributions from each of these sources of He II absorption, including the distribution of line Doppler widths and column densities, the 'He II proximity effect' from the quasar, and a self-consistent derivation of the He II opacity of the universe as a function of the spectrum of ionizing sources, with the assumption that both the clouds and the IGM are photoionized. The He II edge can be fully accounted for by He II line blanketing for reasonable distributions of line widths and column densities in the Ly-alpha forest, provided that the ionizing sources have spectral index alpha(s) greater than 1.5, and any He II proximity effect is neglected. Even with some contribution from a diffuse IGM, it is difficult to account for the edge observed by Jakobsen et al. (1994) with a 'hard' source spectrum (alpha(s) less than 1.3). The proximity effect modifies the relative contributions of the clouds and IGM to tau(He II) near the quasar (z approx. less than z(q)) and markedly increases the amount of He II absorption required. This implies, for example, that to account for the He II edge with line blanketing alone, the minimum spectral index alpha(s) must be increased from 1.5 to 1.9. We demonstrate the need for higher resolution observations that characterize the change in transmission as z approaches z(q) and resolve line-free gaps in the continuum. We set limits on the density of the diffuse IGM and suggest that the IGM and Ly-alpha clouds are likely to be a significant repository for dark baryons

Parametrising Star Formation Histories

We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.Comment: 28 pages, 18 figure

Connecting the Milky Way potential profile to the orbital timescales and spatial structure of the Sagittarius Stream

Recent maps of the halo using RR Lyrae from Pan-STARRS1 have clearly depicted the spatial structure of the Sagittarius stream. These maps show the leading and trailing stream apocenters differ in galactocentric radius by a factor of two, and also resolve substructure in the stream at these apocenters. Here we present dynamical models that reproduce these features of the stream in simple Galactic potentials. We find that debris at the apocenters must be dynamically young, in the sense of being stripped off in the last two pericentric passages, while the Sagittarius dwarf is currently experiencing a third passage. The ratio of apocenters is sensitive to both dynamical friction and the outer slope of the Galactic rotation curve. These dependences can be understood with simple regularities connecting the apocentric radii, circular velocities, and orbital period of the progenitor. The effect of dynamical friction on the stream can be constrained using substructure within the leading apocenter. Our ensembles of models are not intended as statistically proper fits to the stream. Nevertheless, out of the range of models we consider, we consistently find the mass within 100 kpc to be $\sim 7 \times 10^{11} \, M_{\odot}$, with a nearly flat rotation curve between 50 and 100 kpc. This points to a more extended Galactic halo than assumed in some current models. As in previous work, we find prolate or triaxial halos ease agreement with the track of the leading stream. We display the behavior of our models in various observational spaces and characterize the substructure expected within the stream. In particular, the young trailing stream visible near trailing apocenter should exhibit a tight trend of velocity with distance separate from the older debris, and we suggest that this will serve as an especially useful probe of the outer Galactic potential.Comment: Submitted to MNRA