We present results of a deep spectroscopic survey designed to quantify the statistics of the escape of ionizing photons from star-forming galaxies at z~3. We measure the ratio of ionizing to non-ionizing UV flux density <f900/f1500>_obs, where f900 is the mean flux density evaluated over the range [880,910] A. We quantify the emergent ratio of ionizing to non-ionizing UV flux density by analyzing high-S/N composite spectra formed from sub-samples with common observed properties and numbers sufficient to reduce the statistical uncertainty in the modeled IGM+CGM correction to obtain precise values of <f900/f1500>_out, including a full-sample average <f900/f1500>_out=$0.057\pm0.006$. We further show that <f900/f1500>_out increases monotonically with Ly$\alpha$ rest equivalent width, inducing an inverse correlation with UV luminosity as a by-product. We fit the composite spectra using stellar spectral synthesis together with models of the ISM in which a fraction f_c of the stellar continuum is covered by gas with column density N(HI). We show that the composite spectra simultaneously constrain the intrinsic properties of the stars (L900/L1500)_int along with f_c, N(HI), E(B-V), and $f_{esc,abs}$, the absolute escape fraction of ionizing photons. We find a sample-averaged $f_{esc,abs} =0.09\pm0.01$, and that subsamples fall along a linear relation $\langle f_{esc,abs}\rangle \sim 0.75[W(Ly\alpha)/110 A]$. We use the FUV luminosity function, the distribution function $n[W(Ly\alpha)]$, and the relationship between $W(Ly\alpha)$ and <f900/f1500>_out to estimate the total ionizing emissivity of $z\sim3$ star-forming galaxies with Muv < -19.5: $\epsilon_{LyC}\sim 6\times10^{24}$ ergs/s/Hz/Mpc$^3$, exceeding the contribution of QSOs by a factor of $\sim 3$, and accounting for $\sim50$% of the total $\epsilon_{LyC}$ at $z\sim3$ estimated using indirect methods.

We use a large sample of galaxies at z~3 to establish a relationship between reddening, neutral gas covering fraction (fcov(HI)), and the escape of ionizing photons at high redshift. Our sample includes 933 galaxies at z~3, 121 of which have very deep spectroscopic observations (>7 hrs) in the rest-UV (lambda=850-1300 A) with Keck/LRIS. Based on the high covering fraction of outflowing optically-thick HI indicated by the composite spectra of these galaxies, we conclude that photoelectric absorption, rather than dust attenuation, dominates the depletion of ionizing photons. By modeling the composite spectra as the combination of an unattenuated stellar spectrum including nebular continuum emission with one that is absorbed by HI and reddened by a line-of-sight extinction, we derive an empirical relationship between E(B-V) and fcov(HI). Galaxies with redder UV continua have larger covering fractions of HI characterized by higher line-of-sight extinctions. Our results are consistent with the escape of Lya through gas-free lines-of-sight. Covering fractions based on low-ionization interstellar absorption lines systematically underpredict those deduced from the HI lines, suggesting that much of the outflowing gas may be metal-poor. We develop a model which connects the ionizing escape fraction with E(B-V), and which may be used to estimate the escape fraction for an ensemble of high-redshift galaxies. Alternatively, direct measurements of the escape fraction for our data allow us to constrain the intrinsic 900-to-1500 A flux density ratio to be >0.20, a value that favors stellar population models that include weaker stellar winds, a flatter initial mass function, and/or binary evolution. Lastly, we demonstrate how the framework discussed here may be used to assess the pathways by which ionizing radiation escapes from high-redshift galaxies. [Abridged]

We present the first measurements of the shape of the far-ultraviolet (far-UV; lambda=950-1500 A) dust attenuation curve at high redshift (z~3). Our analysis employs rest-frame UV spectra of 933 galaxies at z~3, 121 of which have very deep spectroscopic observations (>7 hrs) at lambda=850-1300 A, with the Low Resolution Imaging Spectrograph on the Keck Telescope. By using an iterative approach in which we calculate the ratios of composite spectra in different bins of continuum color excess, E(B-V), we derive a dust curve that implies a lower attenuation in the far-UV for a given E(B-V) than those obtained with standard attenuation curves. We demonstrate that the UV composite spectra of z~3 galaxies can be modeled well by assuming our new attenuation curve, a high covering fraction of HI, and absorption from the Lyman-Werner bands of H2 with a small (<20%) covering fraction. The low covering fraction of H2 relative to that of the HI and dust suggests that most of the dust in the ISM of typical galaxies at z~3 is unrelated to the catalysis of H2, and is associated with other phases of the ISM (i.e., the ionized and neutral gas). The far-UV dust curve implies a factor of ~2 lower dust attenuation of Lyman continuum (ionizing) photons relative to those inferred from the most commonly assumed attenuation curves for L* galaxies at z~3. Our results may be utilized to assess the degree to which ionizing photons are attenuated in HII regions or, more generally, in the ionized or low column density (N(HI)<10^17.2 cm^-2) neutral ISM of high-redshift galaxies.

We study the Lya profiles of 36 spectroscopically-detected Lya-emitters (LAEs) at z2-3, using Keck MOSFIRE to measure systemic redshifts and velocity dispersions from rest-frame optical nebular emission lines. The sample has a median optical magnitude R=26.0, and ranges from R23 to R>27, corresponding to rest-frame UV absolute magnitudes M_UV-22 to M_UV>-18.2. Dynamical masses range from M_dyn<1.3e8 Msun to M_dyn=6.8e9 Msun, with a median value of M_dyn=6.3e8 Msun. Thirty of the 36 Lya emission lines are redshifted with respect to the systemic velocity with at least 1 sigma significance, and the velocity offset with respect to systemic dv_Lya is correlated with R-band magnitude, M_UV, and the velocity dispersion measured from nebular emission lines with >3 sigma significance: brighter galaxies with larger velocity dispersions tend to have larger values of dv_Lya. We also make use of a comparison sample of 122 UV-color-selected R<25.5 galaxies at z2, all with Lya emission and systemic redshifts measured from nebular emission lines. Using the combined LAE and comparison samples for a total of 158 individual galaxies, we find that dv_Lya is anti-correlated with the Lya equivalent width with 7 sigma significance. Our results are consistent with a scenario in which the Lya profile is determined primarily by the properties of the gas near the systemic redshift; in such a scenario, the opacity to Lya photons in lower mass galaxies may be reduced if large gaseous disks have not yet developed and if the gas is ionized by the harder spectrum of young, low metallicity stars.

It has been claimed that the nova-like cataclysmic variable (CV) AE Aquarii (AE Aqr) is a very-high-energy (VHE, $E>$100 GeV) source both on observational and theoretical grounds. We aim to search for VHE gamma-ray emission from AE Aqr during different states of the source at several wavelengths to confirm or rule out previous claims of detection of gamma-ray emission from this object. We report on observations of AE Aqr performed by MAGIC. The source was observed during 12 hours as part of a multiwavelength campaign carried out between May and June 2012 covering the optical, X-ray, and gamma-ray ranges. Besides MAGIC, the other facilities involved were the KVA, Skinakas, and Vidojevica telescopes in the optical and Swift in X-rays. We calculated integral upper limits coincident with different states of the source in the optical. We computed upper limits to the pulsed emission limiting the signal region to 30% of the phaseogram and we also searched for pulsed emission at different frequencies applying the Rayleigh test. AE Aqr was not detected at VHE energies during the multiwavelength campaign. We establish integral upper limits at the 95\% confidence level for the steady emission assuming the differential flux proportional to a power-law function d\phi/dE \propto E^{-Gamma}, with a Crab-like photon spectral index of Gamma=2.6. The upper limit above 200 GeV is 6.4\times10^{-12} cm^{-2}s^{-1} and above 1 TeV is 7.4\times10^{-13} cm^{-2}s^{-1}. We obtained an upper limit for the pulsed emission of 2.6\times10^{-12} cm^{-2}s^{-1} for energies above 200 GeV. Applying the Rayleigh test for pulsed emission at different frequencies we did not find any significant signal. Our results indicate that AE Aqr is not a VHE gamma-ray emitter at the level of emission previously claimed. We have established the most constraining upper limits for the VHE gamma-ray emission of AE Aqr.

We report on observations of the nova-like cataclysmic variable AE Aqr performed by MAGIC. The observations were part of a quasi-simultaneous multi-wavelength campaign carried out between 2012 May and June covering the optical, UV, X-ray and gamma-ray ranges. MAGIC conducted the campaign and observed the source during 12 hours. The other instruments involved were KVA, Skinakas, and Vidojevica in the optic and Swift in the X-ray. We also used optical data from the AAVSO. The goals were to: monitor the variability of the source at different wavelengths, perform gamma-ray studies coincident with the highest states of the source at the other wavelengths, and confirm or rule out previous claims of detection of very-high-energy emission from this object. We report on a search for steady TeV emission during the whole observation, for variable TeV emission coincident with the highest optical and X-ray states and periodic TeV emission at the 33.08 s rotation period (30.23 mHz rotation frequency) of the white dwarf and its first harmonic (60.46 mHz rotation frequency). These are the first observations under good weather conditions performed by the present generation of IACTs for this object.

The timing and duration of the reionization epoch is crucial to the emergence and evolution of structure in the universe. The relative roles that star-forming galaxies, active galactic nuclei and quasars play in contributing to the metagalactic ionizing background across cosmic time remains uncertain. Deep quasar counts provide insights into their role, but the potentially crucial contribution from star-formation is highly uncertain due to our poor understanding of the processes that allow ionizing radiation to escape into the intergalactic medium (IGM). The fraction of ionizing photons that escape from star-forming galaxies is a fundamental free parameter used in models to "fine-tune" the timing and duration of the reionization epoch that occurred somewhere between 13.4 and 12.7 Gyrs ago (redshifts between 12 > z > 6). However, direct observation of Lyman continuum (LyC) photons emitted below the rest frame \ion{H}{1} ionization edge at 912 \AA\ is increasingly improbable at redshifts z > 3, due to the steady increase of intervening Lyman limit systems towards high z. Thus UV and U-band optical bandpasses provide the only hope for direct, up close and in depth, observations of the types of environment that favor LyC escape. By quantifying the evolution over the past 11 billion years (z < 3) of the relationships between LyC escape and local and global parameters ..., we can provide definitive information on the LyC escape fraction that is so crucial to answering the question of, how did the universe come to be ionized? Here we provide estimates of the ionizing continuum flux emitted by "characteristic" (L_{uv}^*) star-forming galaxies as a function of look back time and escape fraction, finding that at z = 1 (7.6 Gyrs ago) L_{uv}^* galaxies with an escape fraction of 1% have a flux of 10^{-19} ergs cm^{-2} s^{-1} \AA^{-1}.

We present results from the Keck Baryonic Structure Survey (KBSS), a unique spectroscopic survey designed to explore the connection between galaxies and intergalactic baryons. The KBSS is optimized for the redshift range z ~ 2-3, combining S/N ~ 100 Keck/HIRES spectra of 15 hyperluminous QSOs with densely sampled galaxy redshift surveys surrounding each QSO sightline. We perform Voigt profile decomposition of all 6000 HI absorbers within the full Lya forest in the QSO spectra. Here we present the distribution, column density, kinematics, and absorber line widths of HI surrounding 886 star-forming galaxies with 2.0 < z < 2.8 and within 3 Mpc of a QSO sightline. We find that N_HI and the multiplicity of HI components increase rapidly near galaxies. The strongest HI absorbers within ~ 100 physical kpc of galaxies have N_HI ~ 3 dex higher than those near random locations in the IGM. The circumgalactic zone of most enhanced HI absorption (CGM) is found within 300 kpc and 300 km/s of galaxies. Nearly half of absorbers with log(N_HI) > 15.5 are found within the CGM of galaxies meeting our photometric selection, while their CGM occupy only 1.5% of the cosmic volume. The spatial covering fraction, multiplicity of absorption components, and characteristic N_HI remain elevated to transverse distances of 2 physical Mpc. Absorbers with log(N_HI) > 14.5 are tightly correlated with the positions of galaxies, while absorbers with lower N_HI are correlated only on Mpc scales. Redshift anisotropies on Mpc scales indicate coherent infall toward galaxies, while on scales of ~100 physical kpc peculiar velocities of 260 km/s are indicated. The median Doppler widths of absorbers within 1-3 virial radii of galaxies are ~50% larger than randomly chosen absorbers of the same N_HI, suggesting higher gas temperatures and/or increased turbulence likely caused by accretion shocks and/or galactic winds.

Extended nebulae of Ly-alpha emission ("Ly-alpha blobs") are known to be associated with overdense regions at high redshift. Here we present six large Ly-alpha blobs in a previously known protocluster with galaxy overdensity delta ~ 7 at z = 2.3; this is the richest field of giant Ly-alpha blobs detected to date. The blobs have linear sizes >~100 kpc and Ly-alpha luminosities of ~10^43 erg/s. The positions of the blobs define two linear filaments with an extent of at least 12 comoving Mpc; these filaments intersect at the center of one of the blobs. Measurement of the position angles of the blobs indicates that five of the six are aligned with these filaments to within ~10 degrees, suggesting a connection between the physical processes powering extended Ly-alpha emission and those driving structure on larger scales.

We present the results of a photometric and spectroscopic survey of 321 Lyman break galaxies (LBGs) at z ~ 3 to investigate systematically the relationship between Lya emission and stellar populations. Lya equivalent widths (EW) were calculated from rest-frame UV spectroscopy and optical/near-infrared/Spitzer photometry was used in population synthesis modeling to derive the key properties of age, dust extinction, star formation rate (SFR), and stellar mass. We directly compare the stellar populations of LBGs with and without strong Lya emission, where we designate the former group (EW > 20 AA) as Lya emitters (LAEs) and the latter group (EW < 20 AA) as non-LAEs. This controlled method of comparing objects from the same UV luminosity distribution represents an improvement over previous studies in which the stellar populations of LBGs and narrowband-selected LAEs were contrasted, where the latter were often intrinsically fainter in broadband filters by an order of magnitude simply due to different selection criteria. Using a variety of statistical tests, we find that Lya equivalent width and age, SFR, and dust extinction, respectively, are significantly correlated in the sense that objects with strong Lya emission also tend to be older, lower in star formation rate, and less dusty than objects with weak Lya emission, or the line in absorption. We accordingly conclude that, within the LBG sample, objects with strong Lya emission represent a later stage of galaxy evolution in which supernovae-induced outflows have reduced the dust covering fraction. We also examined the hypothesis that the attenuation of Lya photons is lower than that of the continuum, as proposed by some, but found no evidence to support this picture.

We have conducted a spectroscopic survey to find faint quasars (-26.0 < M_{1450} < -22.0) at redshifts z=3.8-5.2 in order to measure the faint end of the quasar luminosity function at these early times. Using available optical imaging data from portions of the NOAO Deep Wide-Field Survey and the Deep Lens Survey, we have color-selected quasar candidates in a total area of 3.76 deg^2. Thirty candidates have R <= 23 mags. We conducted spectroscopic followup for 28 of our candidates and found 23 QSOs, 21 of which are reported here for the first time, in the 3.74 < z <5.06 redshift range. We estimate our survey completeness through detailed Monte Carlo simulations and derive the first measurement of the density of quasars in this magnitude and redshift interval. We find that the binned luminosity function is somewhat affected by the K-correction used to compute the rest-frame absolute magnitude at 1450A. Considering only our R <= 23 sample, the best-fit single power-law (Phi \propto L^beta) gives a faint-end slope beta = -1.6+/-0.2. If we consider our larger, but highly incomplete sample going one magnitude fainter, we measure a steeper faint-end slope -2 < beta < -2.5. In all cases, we consistently find faint-end slopes that are steeper than expected based on measurements at z ~ 3. We combine our sample with bright quasars from the Sloan Digital Sky Survey to derive parameters for a double-power-law luminosity function. Our best fit finds a bright-end slope, alpha = -2.4+/-0.2, and faint-end slope, beta = -2.3+/-0.2, without a well-constrained break luminosity. This is effectively a single power-law, with beta = -2.7+/-0.1. We use these results to place limits on the amount of ultraviolet radiation produced by quasars and find that quasars are able to ionize the intergalactic medium at these redshifts.

We report the discovery of two low-luminosity quasars at z~4, both of which show prominent N IV] 1486A emission. This line is extremely rare in quasar spectra at any redshift; detecting it in two out of a sample of 23 objects (i.e., ~ 9% of the sample) is intriguing and is likely due to the low-luminosity, high-redshift quasar sample we are studying. This is still a poorly explored regime, where contributions from associated, early starbursts may be significant. One interpretation of this line posits photoionization by very massive young stars. Seeing N IV] 1486A emission in a high-redshift quasar may thus be understood in the context of co-formation and early co-evolution of galaxies and their supermassive black holes. Alternatively, we may be seeing a phenomenon related to the early evolution of quasar broad emission line regions. The non-detection (and possibly even broad absorption) of N V 1240A line in the spectrum of one of these quasars may support that interpretation. These two objects may signal a new faint quasar population or an early AGN evolutionary stage at high redshifts.

We present observations of an optically-faint quasar, RD J114816.2+525339, discovered from deep multi-color observations of the field around the z = 6.42 quasar SDSS J1148+5251. The two quasars have a projected separation of 109 arcsec and both are outliers in r-z versus z-J color-color space. Keck spectroscopy reveals RD J114816.2+525339 to be a broad-absorption line quasar at z = 5.70. With z_AB = 23.0, RD J114816.2+525339 is 3.3 mag fainter than SDSS J1148+5251, making it the faintest quasar known at z>5.5. This object was identified in a survey of ~2.5 square degrees. The implied surface density of quasars at these redshifts and luminosities is broadly consistent with previous extrapolations of the faint end of the quasar luminosity function and supports the idea that active galaxies provide only a minor component of the reionizing ultraviolet flux at these redshifts.

Observations of QSOs at z ~ 5.7 - 6.4 show the appearance of Gunn-Peterson troughs around z ~ 6, and a change in the slope of the IGM optical depth tau(z) near z ~ 5.5. These results are interpreted as a signature of the end of the reionization era, which probably started at considerably higher redshifts. However, there also appears to be a substantial cosmic variance in the transmission of the IGM, both along some lines of sight, and among different lines of sight, in this intriguing redshift regime. We suggest that this is indicative of a spatially uneven reionization, possibly caused by the bias-driven primordial clustering of the reionization sources. There is also some independent evidence for a strong clustering of QSOs at z ~ 4 - 5 and galaxies around them, supporting the idea of the strong biasing of the first luminous sources at these redshifts. Larger samples of high-z QSOs are needed in order to provide improved, statistically significant constraints for the models of these phenomena. We expect that the Palomar-Quest (PQ) survey will soon provide a new set of QSOs to be used as cosmological probes in this redshift regime.

One of the new frontiers of astronomical research is the exploration of time variability on the sky at different wavelengths and flux levels. We have carried out a pilot project using DPOSS data to study strong variables and transients, and are now extending it to the new Palomar-QUEST synoptic sky survey. We report on our early findings and outline the methodology to be implemented in preparation for a real-time transient detection pipeline. In addition to large numbers of known types of highly variable sources (e.g., SNe, CVs, OVV QSOs, etc.), we expect to find numerous transients whose nature may be established by a rapid follow-up. Whereas we will make all detected variables publicly available through the web, we anticipate that email alerts would be issued in the real time for a subset of events deemed to be the most interesting. This real-time process entails many challenges, in an effort to maintain a high completeness while keeping the contamination low. We will utilize distributed Grid services developed by the GRIST project, and implement a variety of advanced statistical and machine learning techniques.

Exploration of the time variability on the sky over a broad range of flux levels and wavelengths is rapidly becoming a new frontier of astronomical research. We describe here briefly the Palomar-QUEST survey being carried out from the Samuel Oschin 48-inch Schmidt telescope at Palomar. The following features make the survey an attractive candidate for studying time variability: anticipated survey area of 12,000 - 15,000 sq. degrees in the drift scan mode, point source depth of 21st mag. in I under good conditions, near simultaneous observations in four filters, and at least four passes per year at each location covered. The survey will yield a large number of transients and highly variable sources in the near future and in that sense is a prototype of LSST and Pan-STARRS. We briefly outline our strategy for searching such objects and the proposed pipeline for detecting transients in real-time.

The challenge in searching for non-radio-pulsing isolated neutron stars (INSs) is in excluding association with objects in the very large error boxes (~13", 1 sigma radius) typical of sources from the largest X-ray all-sky survey, the ROSAT All-Sky-Survey/Bright Source Catalog (RASS/BSC). We search for candidate INSs using statistical analysis of optical (USNO-A2), infrared (IRAS), and radio (NVSS) sources near the ROSAT X-ray localization, and show that this selection would find 20% of the INSs in the RASS/BSC. This selection finds 32 candidates at declinations greater than -39 deg, among which are two previously known INSs, seventeen sources which we show are not INSs, and thirteen the classification of which are as yet undetermined. These results require a limit of <67 INSs (90% confidence, full sky, assuming isotropy) in the RASS/BSC. This limit modestly constrains a naive and optimistic model for cooling NSs in the galaxy.