Simulations of thermally broadened HI Lya absorption arising in the warm-hot intergalactic medium

Recent far-ultraviolet (FUV) absorption line measurements of low-redshift quasars have unveiled a population of intervening broad HI Lya absorbers (BLAs) with large Doppler parameters (b> 40 km/s). If the large width of these lines is dominated by thermal line broadening, the BLAs may trace highly-ionized gas in the warm-hot intergalactic medium (WHIM) in the temperature range T ~ 10^5-10^6 K, a gas phase that is expected to contain a large fraction of the baryons at low redshift. In this paper we use a hydrodynamical simulation to study frequency, distribution, physical conditions, and baryon content of the BLAs at z=0. From our simulated spectra we derive a number of BLAs per unit redshift of (dN/dz)_BLA ~ 38 for HI absorbers with log (N(cm^-2)/b(km/s))>10.7, b>40 km/s, and log N(HII)<20.5. The baryon content of these systems is Omega_b(BLA)=0.0121/h_65, which represents ~25 percent of the total baryon budget in our simulation. Our results thus support the idea that BLAs represent a significant baryon reservoir at low redshift. BLAs predominantly trace shock-heated collisionally ionized WHIM gas at temperatures log T~4.4-6.2. About 27 percent of the BLAs in our simulation originate in the photoionized Lya forest (log T<4.3) and their large line widths are determined by non-thermal broadening effects such as unresolved velocity structure and macroscopic turbulence. Our simulation implies that for a large-enough sample of BLAs in FUV spectra it is possible to obtain a reasonable approximation of the baryon content of these systems solely from the measured HI column densities and b values.Comment: 11 pages, 8 figures; minor modifications; accepted for publication in A&

The Nature of Lyman Break Galaxies in Cosmological Hydrodynamic Simulations

What type of objects are being detected as $z\sim 3$ "Lyman break galaxies"? Are they predominantly the most massive galaxies at that epoch, or are many of them smaller galaxies undergoing a short-lived burst of merger-induced star formation? We attempt to address this question using high-resolution cosmological hydrodynamic simulations including star formation and feedback. Our $\Lambda$CDM simulation, together with Bruzual-Charlot population synthesis models, reproduces the observed number density and luminosity function of Lyman break galaxies when dust is incorporated. The inclusion of dust is crucial for this agreement. In our simulation, these galaxies are predominantly the most massive objects at this epoch, and have a significant population of older stars. Nevertheless, it is possible that our simulations lack the resolution and requisite physics to produce starbursts, despite having a physical resolution of \la 700 pc at z=3. Thus we cannot rule out merger-induced starburst galaxies also contributing to the observed population of high-redshift objects.Comment: 5 pages, contribution to the Proceedings of Rencontres Internationales de l'IGRAP, Clustering at High Redshift, Marseille 199

Tracing baryons in the warm-hot intergalactic medium with broad Ly alpha absorption

We discuss physical properties and baryonic content of broad Ly alpha absorbers (BLAs) at low redshift. These absorption systems, recently discovered in high-resolution, high-signal to noise quasar absorption line spectra, possibly trace the warm-hot intergalactic medium (WHIM) in the temperature range between 10^5 and 10^6 K. To extend previous BLA measurements we have analyzed STIS data of the two quasars H 1821+643 and PG 0953+415 and have identified 13 BLA candidates along a total (unblocked) redshift path of dz=0.440. Combining our measurements with previous results for the lines of sight toward PG 1259+593 and PG 1116+215, the resulting new BLA sample consists of 20 reliably detected systems as well as 29 additional tentative cases, implying a BLA number density of dN/dz=22-53. We estimate that the contribution of BLAs to the baryon density at z=0 is Omega_b(BLA)>0.0027 h_70^-1 for absorbers with log (N/b)>11.3. This number indicates that WHIM broad Ly alpha absorbers contain a substantial fraction of the baryons in the local Universe. (Abridged abstract)Comment: 17 pages, 7 figures; Accepted for publication in A&

SIGAME simulations of the [CII], [OI] and [OIII] line emission from star forming galaxies at z ~ 6

Of the almost 40 star forming galaxies at z>~5 (not counting QSOs) observed in [CII] to date, nearly half are either very faint in [CII], or not detected at all, and fall well below expectations based on locally derived relations between star formation rate (SFR) and [CII] luminosity. Combining cosmological zoom simulations of galaxies with SIGAME (SImulator of GAlaxy Millimeter/submillimeter Emission) we have modeled the multi-phased interstellar medium (ISM) and its emission in [CII], [OI] and [OIII], from 30 main sequence galaxies at z~6 with star formation rates ~3-23Msun/yr, stellar masses ~(0.7-8)x10^9Msun, and metallicities ~(0.1-0.4)xZsun. The simulations are able to reproduce the aforementioned [CII]-faintness at z>5, match two of the three existing z>~5 detections of [OIII], and are furthermore roughly consistent with the [OI] and [OIII] luminosity relations with SFR observed for local starburst galaxies. We find that the [CII] emission is dominated by the diffuse ionized gas phase and molecular clouds, which on average contribute ~66% and ~27%, respectively. The molecular gas, which constitutes only ~10% of the total gas mass is thus a more efficient emitter of [CII] than the ionized gas making up ~85% of the total gas mass. A principal component analysis shows that the [CII] luminosity correlates with the star formation activity as well as average metallicity. The low metallicities of our simulations together with their low molecular gas mass fractions can account for their [CII]-faintness, and we suggest these factors may also be responsible for the [CII]-faint normal galaxies observed at these early epochs.Comment: 24 pages, 14 figures. Accepted for publication in the Astrophysical Journa

Molecular Star Formation Rate Indicators in Galaxies

We derive a physical model for the observed relations between star formation rate (SFR) and molecular line (CO and HCN) emission in galaxies, and show how these observed relations are reflective of the underlying star formation law. We do this by combining 3D non-LTE radiative transfer calculations with hydrodynamic simulations of isolated disk galaxies and galaxy mergers. We demonstrate that the observed SFR-molecular line relations are driven by the relationship between molecular line emission and gas density, and anchored by the index of the underlying Schmidt law controlling the SFR in the galaxy. Lines with low critical densities (e.g. CO J=1-0) are typically thermalized and trace the gas density faithfully. In these cases, the SFR will be related to line luminosity with an index similar to the Schmidt law index. Lines with high critical densities greater than the mean density of most of the emitting clouds in a galaxy (e.g. CO J=3-2, HCN J=1-0) will have only a small amount of thermalized gas, and consequently a superlinear relationship between molecular line luminosity and mean gas density. This results in a SFR-line luminosity index less than the Schmidt index for high critical density tracers. One observational consequence of this is a significant redistribution of light from the small pockets of dense, thermalized gas to diffuse gas along the line of sight, and prodigious emission from subthermally excited gas. At the highest star formation rates, the SFR-Lmol slope tends to the Schmidt index, regardless of the molecular transition. The fundamental relation is the Kennicutt-Schmidt law, rather than the relation between SFR and molecular line luminosity. We use these results to make imminently testable predictions for the SFR-molecular line relations of unobserved transitions.Comment: ApJ Accepted - Results remain same as previous version. Content clarified with Referee's comment

The Flux Auto- and Cross-Correlation of the Lyman-alpha Forest. II. Modelling Anisotropies with Cosmological Hydrodynamic Simulations

The isotropy of the Lyman-alpha forest in real-space uniquely provides a measurement of cosmic geometry at z > 2. The angular diameter distance for which the correlation function along the line of sight and in the transverse direction agree corresponds to the correct cosmological model. However, the Lyman-alpha forest is observed in redshift-space where distortions due to Hubble expansion, bulk flows, and thermal broadening introduce anisotropy. Similarly, a spectrograph's line spread function affects the autocorrelation and cross-correlation differently. In this the second paper of a series on using the Lyman-alpha forest observed in pairs of QSOs for a new application of the Alcock-Paczynski (AP) test, these anisotropies and related sources of potential systematic error are investigated with cosmological hydrodynamic simulations. Three prescriptions for galactic outflow were compared and found to have only a marginal effect on the Lyman-alpha flux correlation (which changed by at most 7% with use of the currently favored variable-momentum wind model vs. no winds at all). An approximate solution for obtaining the zero-lag cross-correlation corresponding to arbitrary spectral resolution directly from the zero-lag cross-correlation computed at full-resolution (good to within 2% at the scales of interest) is presented. Uncertainty in the observationally determined mean flux decrement of the Lyman-alpha forest was found to be the dominant source of systematic error; however, this is reduced significantly when considering correlation ratios. We describe a simple scheme for implementing our results, while mitigating systematic errors, in the context of a future application of the AP test.Comment: 20 page

Tracing the Warm Hot Intergalactic Medium in the local Universe

We present a simple method for tracing the spatial distribution and predicting the physical properties of the Warm-Hot Intergalactic Medium (WHIM), from the map of galaxy light in the local universe. Under the assumption that biasing is local and monotonic we map the ~ 2 Mpc/h smoothed density field of galaxy light into the mass density field from which we infer the spatial distribution of the WHIM in the local supercluster. Taking into account the scatter in the WHIM density-temperature and density-metallicity relation, extracted from the z=0 outputs of high-resolution and large box size hydro-dynamical cosmological simulations, we are able to quantify the probability of detecting WHIM signatures in the form of absorption features in the X-ray spectra, along arbitrary directions in the sky. To illustrate the usefulness of this semi-analytical method we focus on the WHIM properties in the Virgo Cluster region.Comment: 16 pages 11 Figures. Discussion clarified, alternative methods proposed. Results unchanged. MNRAS in pres

Line versus Flux Statistics -- Considerations for the Low Redshift Lyman-alpha Forest

The flux/transmission power spectrum has become a popular statistical tool in studies of the high redshift ($z > 2$) Lyman-alpha forest. At low redshifts, where the forest has thinned out into a series of well-isolated absorption lines, the motivation for flux statistics is less obvious. Here, we study the relative merits of flux versus line correlations, and derive a simple condition under which one is favored over the other on purely statistical grounds. Systematic errors probably play an important role in this discussion, and they are outlined as well.Comment: 6 pages, to appear in "The IGM/Galaxy Connection: The Distribution of Baryons at z=0", eds. J. L. Rosenberg and M. E. Putma

The Role of Galactic Winds on Molecular Gas Emission from Galaxy Mergers

We assess the impact of starburst and AGN feedback-driven winds on the CO emission from galaxy mergers, and, in particular, search for signatures of these winds in the simulated CO morphologies and emission line profiles. We do so by combining a 3D non-LTE molecular line radiative transfer code with smoothed particle hydrodynamics (SPH) simulations of galaxy mergers that include prescriptions for star formation, black hole growth, a multiphase interstellar medium (ISM), and the winds associated with star formation and black hole growth. Our main results are: (1) Galactic winds can drive outflows of masses ~10^8-10^9 Msun which may be imaged via CO emission line mapping. (2) AGN feedback-driven winds are able to drive imageable CO outflows for longer periods of time than starburst-driven winds owing to the greater amount of energy imparted to the ISM by AGN feedback compared to star formation. (3) Galactic winds can control the spatial extent of the CO emission in post-merger galaxies, and may serve as a physical motivation for the sub-kiloparsec scale CO emission radii observed in local advanced mergers. (4) Secondary emission peaks at velocities greater than the circular velocity are seen in the CO emission lines in all models. In models with winds, these high velocity peaks are seen to preferentially correspond to outflowing gas entrained in winds, which is not the case in the model without winds. The high velocity peaks seen in models without winds are typically confined to velocity offsets (from the systemic) < 1.7 times the circular velocity, whereas the models with AGN feedback-driven winds can drive high velocity peaks to ~2.5 times the circular velocity.Comment: Accepted by ApJ; Minor revisions; Resolution tests include