Studying the Large Scale Structure and interstellar Medium of Galaxies During the Epochs of Peak Cosmic Star formation and Reionization With Infrared Fine Structure Lines

Infrared (IR) fine-structure (FS) lines from trace metals in the interstellar medium (ISM) of galaxies are valuable diagnostics of the physical conditions in a broad range of astrophysical environments, such gas irradiated by stellar far-ultraviolet (FUV) photons or X-rays from accreting supermassive black holes, called active galactic nuclei (AGN). The transparency of these lines to dust and their high escape fractions into the intergalactic medium (IGM) render them as useful probes to study the epochs of peak cosmic star formation (SF) and Reionization. Chapter 1 of this thesis is a study of the ISM of the Cloverleaf quasar. Observations of IR FS lines from singly ionized carbon and neutral oxygen have allowed us to assess the physical conditions—parametrized by their gas density and the impingent FUV flux—prevalent in atomic gas heated by stellar FUV photons. We find that UV heating from local SF is not sufficient to explain the measured FS and molecular luminosities, and suggest that X-ray heating from the AGN is required to simultaneously explain both sets of data. The general picture of the Cloverleaf ISM that emerges from our composite model is one where the [CII] and [OI]63 line emission is produced primarily within PDRs and HII regions of a 1.3-kpc wide starburst, which is embedded in a denser XDR component that is the dominant source of heating for the CO gas. The fact that the star-forming PDR and HII region gas is co-spatial with the XDR—and within ∼ 650 pc of the accreting black hole—provides strong evidence that SF is ongoing while immersed in a strong X-ray radiation field provided by the nearby AGN. This finding has implications for the co-evolution of supermassive black holes and their host galaxies. The work in this chapter will be submitted for first-author publication imminently. In Chapter 2, we explore the possibility of studying the redshifted far-IR fine-structure line emission using the three-dimensional (3-D) power spectra obtained with an imaging spectrometer. The intensity mapping approach measures the spatio-spectral fluctuations due to line emission from all galaxies, including those below the individual detection threshold. The technique provides 3-D measurements of galaxy clustering and moments of the galaxy luminosity function. Furthermore, the linear portion of the power spectrum can be used to measure the total line emission intensity including all sources through cosmic time with redshift information naturally encoded. Total line emission, when compared to the total star formation activity and/or other line intensities reveals evolution of the interstellar conditions of galaxies in aggregate. As a case study, we consider measurement of [CII] autocorrelation in the 0.5 \u3c z \u3c 1.5 epoch, where interloper lines are minimized, using far-IR/submm balloon-borne and future space-borne instruments with moderate and high sensitivity, respectively. In this context, we compare the intensity mapping approach to blind galaxy surveys based on individual detections. We find that intensity mapping is nearly always the best way to obtain the total line emission because blind, wide-field galaxy surveys lack sufficient depth and deep pencil beams do not observe enough galaxies in the requisite luminosity and redshift bins. Also, intensity mapping is often the most efficient way to measure the power spectrum shape, depending on the details of the luminosity function and the telescope aperture. The work in this chapter has been published in Uzgil et al. (2014). In the final Chapter, we consider the extension of intensity mapping experiments targeting IR FS lines to the late stages of the Epoch of Reionization (EoR), at z ∼ 7. Intensity mapping experiments of emission lines from the ISM of galaxies are highly complementary to experiments that are aiming to detect the 21 cm power spectrum during the same epoch, as the former is a direct probe of the sources of Reionization, and the latter is a probe of the effect of those sources on the surrounding IGM. Since current and planned observations are limited by cosmic variance at the bright end of the galaxy luminosity function, and will not be able to detect the faintest galaxies responsible for a significant fraction of the ionizing photon supply during EoR, intensity mapping is an appealing approach to study the nature and evolution of galaxies during this stage in the history of the Universe. Again, the utility of FS lines as ISM diagnostics, combined with the ability of intensity mapping to measure redshift-evolution in mean intensity of individual lines or the evolution of line ratios (constructed from multiple cross-power spectra), presents a unique and tantalizing opportunity to directly observe changes in properties of interstellar medium (such as hardness of the ionizing spectrum in galaxies and metallicity) that are important to galaxy evolution studies

Statistics of Cosmological Black Hole Jet Sources: Blazar Predictions for GLAST

A study of the statistics of cosmological black-hole jet sources is applied to EGRET blazar data, and predictions are made for GLAST. Black-hole jet sources are modeled as collimated relativistic plasma outflows with radiation beamed along the jet axis due to strong Doppler boosting. The comoving rate density of blazar flares is assumed to follow a blazar formation rate (BFR), modeled by analytic functions based on astronomical observations and fits to EGRET data. The redshift and size distributions of gamma-ray blazars observed with EGRET, separated into BL Lac object (BL) and flat spectrum radio quasar (FSRQ) distributions, are fit with monoparametric functions for the distributions of the jet Lorentz factor \Gamma, comoving directional power l'_e, and spectral slope. A BFR factor ~10 x greater at z ~ 1 than at present is found to fit the FSRQ data. A smaller comoving rate density and greater luminosity of BL flares at early times compared to the present epoch fits the BL data. Based on the EGRET observations, ~1000 blazars consisting of ~800 FSRQs and FR2 radio galaxies and ~200 BL Lacs and FR1 radio galaxies will be detected with GLAST during the first year of the mission. Additional AGN classes, such as hard-spectrum BL Lacs that were mostly missed with EGRET, could add more GLAST sources. The FSRQ and BL contributions to the EGRET gamma-ray background at 1 GeV are estimated at the level of ~10 - 15% and ~2 - 4%, respectively. EGRET and GLAST sensitivities to blazar flares are considered in the optimal case, and a GLAST analysis method for blazar detection is outlined.Comment: 17 pages, 9 figures, ApJ, in press, v.660, May 1, 2007 (minor changes from previous version

The QSO evolution derived from the HBQS and other complete QSO surveys

An ESO Key programme dedicated to an Homogeneous Bright QSO Survey (HBQS) has been completed. 327 QSOs (Mb<-23, 0.3<z<2.2) have been selected over 555 deg^2 with 15<B<18.75. For B<16.4 the QSO surface density turns out to be a factor 2.2 higher than what measured by the PG survey, corresponding to a surface density of 0.013+/-.006 deg^{-2}. If the Edinburgh QSO Survey is included, an overdensity of a factor 2.5 is observed, corresponding to a density of 0.016+/-0.005 deg^{-2}. In order to derive the QSO optical luminosity function (LF) we used Monte Carlo simulations that take into account of the selection criteria, photometric errors and QSO spectral slope distribution. The LF can be represented with a Pure Luminosity Evolution (L(z)\propto(1+z)^k) of a two power law both for q_0=0.5 and q_0=0.1. For q_0=0.5 k=3.26, slower than the previous Boyle's (1992) estimations of k=3.45. A flatter slope beta=-3.72 of the bright part of the LF is also required. The observed overdensity of bright QSOs is concentrated at z<0.6. It results that in the range 0.3<z<0.6 the luminosity function is flatter than observed at higher redshifts. In this redshift range, for Mb<-25, 32 QSOs are observed instead of 19 expected from our best-fit PLE model. This feature requires a luminosity dependent luminosity evolution in order to satisfactorily represent the data in the whole 0.3<z<2.2 interval.Comment: Invited talk in "Wide Field Spectroscopy" (20-24 May 1996, Athens), eds. M. Kontizas et al. 6 pages and 3 eps figures, LaTex file, uses epfs.sty and crckapb.sty (included

On the evolutionary behaviour of BL Lac objects

We present a new well defined sample of BL Lac objects selected from the ROSAT All-Sky Survey (RASS). The sample consists of 39 objects with 35 forming a flux limited sample down to f_X = 8 x 10^{-13} cgs, redshifts are known for 33 objects (and 31 of the complete sample). X-ray spectral properties were determined for each object individually with the RASS data. The luminosity function of RASS selected BL Lac objects is compatible with results provided by objects selected with the Einstein observatory, but the RASS selected sample contains objects with luminosities at least tenfold higher. Our analysis confirms the negative evolution for X-ray selected BL Lac objects found in a sample by the Einstein observatory, the parameterization provides similar results. A subdivision of the sample into halves according to the X-ray to optical flux ratio yielded unexpected results. The extremely X-ray dominated objects have higher redshifts and X-ray luminosities and only this subgroup shows clear signs of strong negative evolution. The evolutionary behaviour of objects with an intermediate spectral energy distribution between X-ray and radio dominated is compatible with no evolution at all. Consequences for unified schemes of X-ray and radio selected BL Lac objects are discussed.We suggest that the intermediate BL Lac objects are the basic BL Lac population. The distinction between the two subgroups can be explained if extreme X-ray dominated BL Lac objects are observed in a state of enhanced X-ray activity.Comment: 14 pages incl. 8 figures, accepted by A&

Follow-Up Chandra Observations of Three Candidate Tidal Disruption Events

Large-amplitude, high-luminosity soft X-ray flares were detected by the ROSAT All-Sky Survey in several galaxies with no evidence of Seyfert activity in their ground-based optical spectra. These flares had the properties predicted for a tidal disruption of a star by a central supermassive black hole. We report Chandra observations of three of these galaxies taken a decade after their flares that reveal weak nuclear X-ray sources that are from 240 to 6000 times fainter than their luminosities at peak, supporting the theory that these were special events and not ongoing active galactic nucleus (AGN) variability. The decline of RX J1624.9+7554 by a factor of 6000 is consistent with the (t-t_D)^(-5/3) decay predicted for the fall-back phase of a tidal disruption event, but only if ROSAT was lucky enough to catch the event exactly at its peak in 1990 October. RX J1242.6-1119A has declined by a factor of 240, also consistent with (t-t_D)^(-5/3). In the H II galaxy NGC 5905 we find only resolved, soft X-ray emission that is undoubtedly associated with starburst activity. When accounting for the starburst component, the ROSAT observations of NGC 5905, as well as the Chandra upper limit on its nuclear flux, are consistent with a (t-t_D)^(-5/3) decay by at least a factor of 1000. Although we found weak Seyfert~2 emission lines in Hubble Space Telescope spectra of NGC 5905, indicating that a low-luminosity AGN was present prior to the X-ray flare, we favor a tidal disruption explanation for the flare itself.Comment: 17 pages, 4 figures, to appear in ApJ April 1 200