The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$\alpha$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$\alpha$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

The Nearby Supernova Factory has made spectrophotometric observations of Type Ia supernovae since $2004$. This work presents an interim version of the data produced, including $210$ supernovae observed between $2004$ and $2013$.

With a booming number of Type Ia supernovae (SNe Ia) discovered within a few days of their explosions, a fraction of SNe Ia that show luminosity excess in the early phase (early-excess SNe Ia) have been confirmed. In this article, we report early-phase observations of seven photometrically normal SNe Ia (six early detections and one deep non-detection limit) at the COSMOS field through a half-year transient survey as a part of the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP). In particular, a blue light-curve excess was discovered for HSC17bmhk, a normal SN Ia with rise time longer than 18.8 days, during the first four days after the discovery. The blue early excess in optical wavelength can be explained not only by interactions with a non-degenerate companion or surrounding dense circumstellar matter but also radiation powered by radioactive decays of $^{56}$Ni at the surface of the SN ejecta. Given the growing evidence of the early-excess discoveries in normal SNe Ia that have longer rise times than the average and a similarity in the nature of the blue excess to a luminous SN Ia subclass, we infer that early excess discovered in HSC17bmhk and other normal SNe Ia are most likely attributed to radioactive $^{56}$Ni decay at the surface of the SN ejecta. In order to successfully identify normal SNe Ia with early excess similar to that of HSC17bmhk, early UV photometries or high-cadence blue-band surveys are necessary.

Rapidly evolving transients form a new class of transients which show shorter timescales of the light curves than those of typical core-collapse and thermonuclear supernovae. We performed a systematic search for rapidly evolving transients using the deep data taken with the Hyper Suprime-Cam Subaru Strategic Program Transient Survey. By measuring the timescales of the light curves of 1824 transients, we identified 5 rapidly evolving transients. Our samples are found in a wide range of redshifts (0.3 $\le$ z $\le$ 1.5) and peak absolute magnitudes ($-$17 $\ge$ $M_i$ $\ge$ $-$20). The properties of the light curves are similar to those of the previously discovered rapidly evolving transients. They show a relatively blue spectral energy distribution, with the best-fit blackbody of 8,000 - 18,000 K. We show that some of the transients require power sources other than the radioactive decays of $^{56}$Ni because of their high peak luminosities and short timescales. The host galaxies of all the samples are star-forming galaxies, suggesting a massive star origin for the rapidly evolving transients. The event rate is roughly estimated to be $\sim$4,000 events yr$^{-1}$ Gpc$^{-3}$, which is about 1 $\%$ of core-collapse supernovae.

Type Ia Supernovae (SNe Ia) are widely used to measure the expansion of the Universe. Improving distance measurements of SNe Ia is one technique to better constrain the acceleration of expansion and determine its physical nature. This document develops a new SNe Ia spectral energy distribution (SED) model, called the SUpernova Generator And Reconstructor (SUGAR), which improves the spectral description of SNe Ia, and consequently could improve the distance measurements. This model is constructed from SNe Ia spectral properties and spectrophotometric data from The Nearby Supernova Factory collaboration. In a first step, a PCA-like method is used on spectral features measured at maximum light, which allows us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties are used to extract the average extinction curve. Third, an interpolation using Gaussian Processes facilitates using data taken at different epochs during the lifetime of a SN Ia and then projecting the data on a fixed time grid. Finally, the three steps are combined to build the SED model as a function of time and wavelength. This is the SUGAR model. The main advancement in SUGAR is the addition of two additional parameters to characterize SNe Ia variability. The first is tied to the properties of SNe Ia ejecta velocity, the second is correlated with their calcium lines. The addition of these parameters, as well as the high quality the Nearby Supernova Factory data, makes SUGAR an accurate and efficient model for describing the spectra of normal SNe Ia as they brighten and fade. The performance of this model makes it an excellent SED model for experiments like ZTF, LSST or WFIRST.

We report our observations of HSC16aayt (SN 2016jiu), which was discovered by the Subaru/Hyper Suprime-Cam (HSC) transient survey conducted as part of Subaru Strategic Program (SSP). It shows very slow photometric evolution and its rise time is more than 100 days. The optical magnitude change in 400 days remains within 0.6 mag. Spectra of HSC16aayt show a strong narrow emission line and we classify it as a Type IIn supernova. The redshift of HSC16aayt is 0.6814 +/- 0.0002 from the spectra. Its host galaxy center is at 5 kpc from the supernova location and HSC16aayt might be another example of isolated Type IIn supernovae, although the possible existence of underlying star forming activity of the host galaxy at the supernova location is not excluded.

We present an overview of a deep transient survey of the COSMOS field with the Subaru Hyper Suprime-Cam (HSC). The survey was performed for the 1.77 deg$^2$ ultra-deep layer and 5.78 deg$^2$ deep layer in the Subaru Strategic Program over 6- and 4-month periods from 2016 to 2017, respectively. The ultra-deep layer shows a median depth per epoch of 26.4, 26.3, 26.0, 25.6, and 24.6 mag in $g$, $r$, $i$, $z$, and $y$ bands, respectively; the deep layer is $\sim0.6$ mag shallower. In total, 1,824 supernova candidates were identified. Based on light curve fitting and derived light curve shape parameter, we classified 433 objects as Type Ia supernovae (SNe); among these candidates, 129 objects have spectroscopic or COSMOS2015 photometric redshifts and 58 objects are located at $z > 1$. Our unique dataset doubles the number of Type Ia SNe at $z > 1$ and enables various time-domain analyses of Type II SNe, high redshift superluminous SNe, variable stars, and active galactic nuclei.

We present Keck spectroscopic observations of three probable high redshift superluminous supernovae (SLSNe) from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), confirming redshifts of 1.851, 1.965 and 2.399. The host galaxies were selected for transient monitoring from multi-band photometric redshifts. The supernovae are detected during their rise, and the classically scheduled spectra are collected near maximum light. The restframe far-ultraviolet (FUV;$\sim$1000A-2500A) spectra include a significant host galaxy flux contribution and we compare our host galaxy subtracted spectra to UV-luminous SNe from the literature. While the signal-to-noise ratios of the spectra presented here are sufficient for redshift confirmation, supernova spectroscopic type confirmation remains inconclusive. The success of the first SHIZUCA Keck spectroscopic follow-up program demonstrates that campaigns such as SHIZUCA are capable of identifying high redshift SLSNe with sufficient accuracy, speed and depth for rapid, well-cadenced and informative follow-up.

We report our first discoveries of high-redshift supernovae from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), the transient survey using Subaru/Hyper Suprime-Cam. We report the discovery of three supernovae at the spectroscopically-confirmed redshifts of 2.399 (HSC16adga), 1.965 (HSC17auzg), and 1.851 (HSC17dbpf), and two supernova candidates with the host-galaxy photometric redshifts of 3.2 (HSC16apuo) and 4.2 (HSC17dsid), respectively. In this paper, we present their photometric properties and the spectroscopic properties of the confirmed high-redshift supernovae are presented in the accompanying paper (Curtin et al. 2019). The supernovae with the confirmed redshifts of z ~ 2 have the rest ultraviolet peak magnitudes close to -21 mag and they are likely superluminous supernovae. The discovery of three supernovae at z ~ 2 roughly corresponds to the approximate event rate of ~ 900 +/- 520 Gpc-3 yr-1 with Poisson error, which is consistent with the total superluminous supernova rate estimated by extrapolating the local rate based on the cosmic star-formation history. Adding unconfirmed superluminous supernova candidates would increase the event rate. Our superluminous supernova candidates at the redshifts of around 3 and 4 indicate the approximate superluminous supernova rates of ~ 400 +/- 400 Gpc-3 yr-1 (z ~ 3) and ~ 500 +/- 500 Gpc-3 yr-1 (z ~ 4) with Poisson errors. Our initial results demonstrate the outstanding capability of Hyper Suprime-Cam to discover high-redshift supernovae.

Type Ia supernova cosmology depends on the ability to fit and standardize observations of supernova magnitudes with an empirical model. We present here a series of new models of Type Ia Supernova spectral time series that capture a greater amount of supernova diversity than possible with the models that are currently customary. These are entitled SuperNova Empirical MOdels (\textsc{SNEMO}\footnote{https://snfactory.lbl.gov/snemo}). The models are constructed using spectrophotometric time series from $172$ individual supernovae from the Nearby Supernova Factory, comprising more than $2000$ spectra. Using the available observations, Gaussian Processes are used to predict a full spectral time series for each supernova. A matrix is constructed from the spectral time series of all the supernovae, and Expectation Maximization Factor Analysis is used to calculate the principal components of the data. K-fold cross-validation then determines the selection of model parameters and accounts for color variation in the data. Based on this process, the final models are trained on supernovae that have been dereddened using the Fitzpatrick and Massa extinction relation. Three final models are presented here: \textsc{SNEMO2}, a two-component model for comparison with current Type~Ia models; \textsc{SNEMO7}, a seven component model chosen for standardizing supernova magnitudes which results in a total dispersion of $0.100$~mag for a validation set of supernovae, of which $0.087$~mag is unexplained (a total dispersion of $0.113$~mag with unexplained dispersion of $0.097$~mag is found for the total set of training and validation supernovae); and \textsc{SNEMO15}, a comprehensive $15$ component model that maximizes the amount of spectral time series behavior captured.

Hydrogen-poor superluminous supernovae (SLSN-I) are a class of rare and energetic explosions discovered in untargeted transient surveys in the past decade. The progenitor stars and the physical mechanism behind their large radiated energies ($\sim10^{51}$ erg) are both debated, with one class of models primarily requiring a large rotational energy, while the other requires very massive progenitors to either convert kinetic energy into radiation via interaction with circumstellar material (CSM), or engender a pair-instability explosion. Observing the structure of the CSM around SLSN-I offers a powerful test of some scenarios, though direct observations are scarce. Here, we present a series of spectroscopic observations of the SLSN-I iPTF16eh, which reveal both absorption and time- and frequency-variable emission in the Mg II resonance doublet. We show that these observations are naturally explained as a resonance scattering light echo from a circumstellar shell. Modeling the evolution of the emission, we find a shell radius of 0.1 pc and velocity of 3300 km s$^{-1}$, implying the shell was ejected three decades prior to the supernova explosion. These properties match theoretical predictions of pulsational pair-instability shell ejections, and imply the progenitor had a He core mass of $\sim 50-55~{\rm M}_{\odot}$, corresponding to an initial mass of $\sim 115~{\rm M}_{\odot}$.

Faint star-forming galaxies at $z\sim 2-3$ can be used as alternative background sources to probe the Lyman-$\alpha$ forest in addition to quasars, yielding high sightline densities that enable 3D tomographic reconstruction of the foreground absorption field. Here, we present the first data release from the COSMOS Lyman-Alpha Mapping And Mapping Observations (CLAMATO) Survey, which was conducted with the LRIS spectrograph on the Keck-I telescope. Over an observational footprint of 0.157$\mathrm{deg}^2$ within the COSMOS field, we used 240 galaxies and quasars at $2.17<z<3.00$, with a mean comoving transverse separation of $2.37\,h^{-1}\,\mathrm{Mpc^3}$, as background sources probing the foreground Lyman-$\alpha$ forest absorption at $2.05<z<2.55$. The Lyman-$\alpha$ forest data was then used to create a Wiener-filtered tomographic reconstruction over a comoving volume of $3.15\,\times 10^5\,h^{-3}\,\mathrm{Mpc^3}$ with an effective smoothing scale of $2.5\,h^{-1}\,\mathrm{Mpc}$. In addition to traditional figures, this map is also presented as a virtual-reality YouTube360 video visualization and manipulable interactive figure. We see large overdensities and underdensities that visually agree with the distribution of coeval galaxies from spectroscopic redshift surveys in the same field, including overdensities associated with several recently-discovered galaxy protoclusters in the volume. This data release includes the redshift catalog, reduced spectra, extracted Lyman-$\alpha$ forest pixel data, and tomographic map of the absorption.

As part of an on-going effort to identify, understand and correct for astrophysics biases in the standardization of Type Ia supernovae (SNIa) for cosmology, we have statistically classified a large sample of nearby SNeIa into those located in predominantly younger or older environments. This classification is based on the specific star formation rate measured within a projected distance of 1kpc from each SN location (LsSFR). This is an important refinement compared to using the local star formation rate directly as it provides a normalization for relative numbers of available SN progenitors and is more robust against extinction by dust. We find that the SNeIa in predominantly younger environments are DY=0.163\pm0.029 mag (5.7 sigma) fainter than those in predominantly older environments after conventional light-curve standardization. This is the strongest standardized SN Ia brightness systematic connected to host-galaxy environment measured to date. The well-established step in standardized brightnesses between SNeIa in hosts with lower or higher total stellar masses is smaller at DM=0.119\pm0.032 mag (4.5 sigma), for the same set of SNeIa. When fit simultaneously, the environment age offset remains very significant, with DY=0.129\pm0.032 mag (4.0 sigma), while the global stellar mass step is reduced to DM=0.064\pm0.029 mag (2.2 sigma). Thus, approximately 70% of the variance from the stellar mass step is due to an underlying dependence on environment-based progenitor age. Standardization using only the SNeIa in younger environments reduces the total dispersion from 0.142\pm0.008 mag to 0.120\pm0.010 mag. We show that as environment ages evolve with redshift a strong bias on measurement of the dark energy equation of state parameters can develop. Fortunately, data to measure and correct for this effect is likely to be available for many next-generation experiments. [abstract shorten]

We present the discovery and measurements of a gravitationally lensed supernova (SN) behind the galaxy cluster MOO J1014+0038. Based on multi-band Hubble Space Telescope and Very Large Telescope (VLT) photometry of the supernova, and VLT spectroscopy of the host galaxy, we find a 97.5% probability that this SN is a SN Ia, and a 2.5% chance of a CC SN. Our typing algorithm combines the shape and color of the light curve with the expected rates of each SN type in the host galaxy. With a redshift of 2.2216, this is the highest redshift SN Ia discovered with a spectroscopic host-galaxy redshift. A further distinguishing feature is that the lensing cluster, at redshift 1.23, is the most distant to date to have an amplified SN. The SN lies in the middle of the color and light-curve shape distributions found at lower redshift, disfavoring strong evolution to z = 2.22. We estimate an amplification due to gravitational lensing of 2.8+0.6-0.5 (1.10 +- 0.23 mag)---compatible with the value estimated from the weak-lensing-derived mass and the mass-concentration relation from LambdaCDM simulations---making it the most amplified SN Ia discovered behind a galaxy cluster.

Type Ia Supernovae (SNe Ia) are widely used to measure the expansion of the Universe. To perform such measurements the luminosity and cosmological redshift ($z$) of the SNe Ia have to be determined. The uncertainty on $z$ includes an unknown peculiar velocity, which can be very large for SNe Ia in the virialized cores of massive clusters. We determine which SNe Ia exploded in galaxy clusters. We then study how the correction for peculiar velocities of host galaxies inside the clusters improves the Hubble residuals. Using 145 SNe Ia from the Nearby Supernova Factory we found 11 candidates for membership in clusters. To estimate the redshift of a cluster we applied the bi-weight technique. Then, we use the galaxy cluster redshift instead of the host galaxy redshift to construct the Hubble diagram. For SNe Ia inside galaxy clusters the dispersion around the Hubble diagram when peculiar velocities are taken into account is smaller in comparison with a case without peculiar velocity correction, with a $wRMS=0.130\pm0.038$ mag instead of $wRMS=0.137\pm0.036$ mag. The significance of this improvement is 3.58 $\sigma$. If we remove the very nearby Virgo cluster member SN2006X ($z<0.01$) from the analysis, the significance decreases to 1.34 $\sigma$. The peculiar velocity correction is found to be highest for the SNe Ia hosted by blue spiral galaxies, with high local specific star formation rate and smaller stellar mass, seemingly counter to what might be expected given the heavy concentration of old, massive elliptical galaxies in clusters. As expected, the Hubble residuals of SNe Ia associated with massive galaxy clusters improve when the cluster redshift is taken as the cosmological redshift of the SN. This fact has to be taken into account in future cosmological analyses in order to achieve higher accuracy for cosmological redshift measurements. Here we provide an approach to do so.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

Context. Observations of Type Ia supernovae (SNe Ia) can be used to derive accurate cosmological distances through empirical standardization techniques. Despite this success neither the progenitors of SNe Ia nor the explosion process are fully understood. The U-band region has been less well observed for nearby SNe, due to technical challenges, but is the most readily accessible band for high-redshift SNe. Aims. Using spectrophotometry from the Nearby Supernova Factory, we study the origin and extent of U-band spectroscopic variations in SNe Ia and explore consequences for their standardization and the potential for providing new insights into the explosion process. Methods. We divide the U-band spectrum into four wavelength regions {\lambda}(uNi), {\lambda}(uTi), {\lambda}(uSi) and {\lambda}(uCa). Two of these span the Ca H&K {\lambda}{\lambda} 3934, 3969 complex. We employ spectral synthesis using SYNAPPS to associate the two bluer regions with Ni/Co and Ti. Results. (1) The flux of the uTi feature is an extremely sensitive temperature/luminosity indicator, standardizing the SN peak luminosity to 0.116 $\pm$ 0.011 mag RMS. A traditional SALT2.4 fit on the same sample yields a 0.135 mag RMS. Standardization using uTi also reduces the difference in corrected magnitude between SNe originating from different host galaxy environments. (2) Early U-band spectra can be used to probe the Ni+Co distribution in the ejecta, thus offering a rare window into the source of lightcurve power. (3) The uCa flux further improves standardization, yielding a 0.086 $\pm$ 0.010 mag RMS without the need to include an additional intrinsic dispersion to reach {\chi}$^2$/dof $\sim$ 1. This reduction in RMS is partially driven by an improved standardization of Shallow Silicon and 91T-like SNe.

We report the discovery of stars that show spectra very close to the black-body radiation. We found 17 such stars out of 798,593 stars in the Sloan Digital Sky Survey (SDSS) spectroscopic data archives. We discuss the value of these stars for the calibration of photometry, whatever is the physical nature of these stars. This gives us a chance to examine the accuracy of the zero point of SDSS photometry across various passbands: we conclude that the zero point of SDSS photometric system is internally consistent across its five passbands to the level below 0.01 mag. We may also examine the consistency of the zero-points between UV photometry of Galaxy Evolution Explorer and SDSS, and IR photometry of Wide-field Infrared Survey Explorer against SDSS. These stars can be used as not only photometric but spetrophotometric standard stars. We suggest that these stars showing the featureless black-body like spectrum of the effective temperature of $10000\pm1500$K are consistent with DB white dwarfs with the temperature too low to develop helium absorption features.

Type Ia supernovae (SNe Ia) arise from the thermonuclear explosion of carbon-oxygen white dwarfs. Though the uniformity of their light curves makes them powerful cosmological distance indicators, long-standing issues remain regarding their progenitors and explosion mechanisms. Recent detection of the early ultraviolet pulse of a peculiar subluminous SN Ia has been claimed as new evidence for the companion-ejecta interaction through the single-degenerate channel. Here, we report the discovery of a prominent but red optical flash at $\sim$ 0.5 days after the explosion of a SN Ia which shows hybrid features of different SN Ia sub-classes: a light curve typical of normal-brightness SNe Ia, but with strong titanium absorptions, commonly seen in the spectra of subluminous ones. We argue that the early flash of such a hybrid SN Ia is different from predictions of previously suggested scenarios such as the companion-ejecta interaction. Instead it can be naturally explained by a SN explosion triggered by a detonation of a thin helium shell either on a near-Chandrasekhar-mass white dwarf ($\gtrsim$ 1.3 M$_{\odot}$) with low-yield $^{56}$Ni or on a sub-Chandrasekhar-mass white dwarf ($\sim$ 1.0 M$_{\odot}$) merging with a less massive white dwarf. This finding provides compelling evidence that one branch of the previously proposed explosion models, the helium-ignition scenario, does exist in nature, and such a scenario may account for explosions of white dwarfs in a wider mass range in contrast to what was previously supposed.

The scientific yield of current and future optical surveys is increasingly limited by systematic uncertainties in the flux calibration. This is the case for Type Ia supernova (SN Ia) cosmology programs, where an improved calibration directly translates into improved cosmological constraints. Current methodology rests on models of stars. Here we aim to obtain flux calibration that is traceable to state-of-the-art detector-based calibration. We present the SNIFS Calibration Apparatus (SCALA), a color (relative) flux calibration system developed for the SuperNova Integral Field Spectrograph (SNIFS), operating at the University of Hawaii 2.2 m (UH 88) telescope. By comparing the color trend of the illumination generated by SCALA during two commissioning runs, and to previous laboratory measurements, we show that we can determine the light emitted by SCALA with a long-term repeatability better than 1%. We describe the calibration procedure necessary to control for system aging. We present measurements of the SNIFS throughput as estimated by SCALA observations. The SCALA calibration unit is now fully deployed at the UH\,88 telescope, and with it color-calibration between 4000 {\AA} and 9000 {\AA} is stable at the percent level over a one-year baseline.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

Correction of Type Ia Supernova brightnesses for extinction by dust has proven to be a vexing problem. Here we study the dust foreground to the highly reddened SN 2012cu, which is projected onto a dust lane in the galaxy NGC 4772. The analysis is based on multi-epoch, spectrophotometric observations spanning 3,300 - 9,200 {\AA}, obtained by the Nearby Supernova Factory. Phase-matched comparison of the spectroscopically twinned SN 2012cu and SN 2011fe across 10 epochs results in the best-fit color excess of (E(B-V), RMS) = (1.00, 0.03) and total-to-selective extinction ratio of (RV , RMS) = (2.95, 0.08) toward SN 2012cu within its host galaxy. We further identify several diffuse interstellar bands, and compare the 5780 {\AA} band with the dust-to-band ratio for the Milky Way. Overall, we find the foreground dust-extinction properties for SN 2012cu to be consistent with those of the Milky Way. Furthermore we find no evidence for significant time variation in any of these extinction tracers. We also compare the dust extinction curve models of Cardelli et al. (1989), O'Donnell (1994), and Fitzpatrick (1999), and find the predictions of Fitzpatrick (1999) fit SN 2012cu the best. Finally, the distance to NGC4772, the host of SN 2012cu, at a redshift of z = 0.0035, often assigned to the Virgo Southern Extension, is determined to be 16.6$\pm$1.1 Mpc. We compare this result with distance measurements in the literature.

We present spatial correlations of galaxies and IGM HI in the COSMOS/UltraVISTA 1.62 deg$^2$ field. Our data consist of 13,415 photo-$z$ galaxies at $z\sim2-3$ with $K_s<23.4$ and the Ly$\alpha$ forest absorptions in the background quasar spectra selected from SDSS data with no signature of damped Ly$\alpha$ system contamination. We estimate a galaxy overdensity $\delta_{gal}$ in an impact parameter of 2.5 pMpc, and calculate the Ly$\alpha$ forest fluctuations $\delta_{\langle F\rangle}$ whose negative values correspond to the strong Ly$\alpha$ forest absorptions. We identify weak evidence of an anti-correlation between $\delta_{gal}$ and $\delta_{\langle F\rangle}$ with a Spearman's rank correlation coefficient of $-0.39$ suggesting that the galaxy overdensities and the Ly$\alpha$ forest absorptions positively correlate in space at the $\sim90\%$ confidence level. This positive correlation indicates that high-$z$ galaxies exist around an excess of HI gas in the Ly$\alpha$ forest. We find four cosmic volumes, dubbed $A_{obs}$-$D_{obs}$, that have extremely large (small) values of $\delta_{gal} \simeq0.8$ ($-1$) and $\delta_{\langle F\rangle}$ $\simeq0.1$ ($-0.4$), three out of which, $B_{obs}$-$D_{obs}$, significantly depart from the correlation, and weaken the correlation signal. We perform cosmological hydrodynamical simulations, and compare with our observational results. Our simulations reproduce the correlation, agreeing with the observational results. Moreover, our simulations have model counterparts of $A_{obs}$-$D_{obs}$, and suggest that the observations pinpoint, by chance, a galaxy overdensity like a proto-cluster, gas filaments lying on the sightline, a large void, and orthogonal low-density filaments. Our simulations indicate that the significant departures of $B_{obs}$-$D_{obs}$ are produced by the filamentary large-scale structures and the observation sightline effects.

We report on the diversity in quasar spectra from the Baryon Oscillation Spectroscopic Survey. After filtering the spectra to mitigate selection effects and Malmquist bias associated with a nearly flux-limited sample, we create high signal-to-noise ratio composite spectra from 58,656 quasars (2.1 \le z \le 3.5), binned by luminosity, spectral index, and redshift. With these composite spectra, we confirm the traditional Baldwin effect (BE, i.e., the anticorrelation of C IV equivalent width (EW) and luminosity) that follows the relation W_\lambda \propto L^{\beta_w} with slope \beta_w = -0.35 \pm 0.004, -0.35 \pm 0.005, and -0.41 \pm 0.005 for z = 2.25, 2.46, and 2.84, respectively. In addition to the redshift evolution in the slope of the BE, we find redshift evolution in average quasar spectral features at fixed luminosity. The spectroscopic signature of the redshift evolution is correlated at 98% with the signature of varying luminosity, indicating that they arise from the same physical mechanism. At a fixed luminosity, the average C IV FWHM decreases with increasing redshift and is anti-correlated with C IV EW. The spectroscopic signature associated with C IV FWHM suggests that the trends in luminosity and redshift are likely caused by a superposition of effects that are related to black hole mass and Eddington ratio. The redshift evolution is the consequence of a changing balance between these two quantities as quasars evolve toward a population with lower typical accretion rates at a given black hole mass.

We compare the isophote shape parameter $a_{4}$ of early-type galaxies (ETGs) between $z\sim1$ and 0 as a proxy for dynamics to investigate the epoch at which the dynamical properties of ETGs are established, using cluster ETG samples with stellar masses of $\log(M_{*}/M_{\odot})\geq10.5$ which have spectroscopic redshifts. We have 130 ETGs from the Hubble Space Telescope Cluster Supernova Survey for $z\sim1$ and 355 ETGs from the Sloan Digital Sky Survey for $z\sim0$. We have developed an isophote shape analysis method which can be used for high-redshift galaxies and has been carefully compared with published results. We have applied the same method for both the $z\sim1$ and $0$ samples. We find similar dependence of the $a_{4}$ parameter on the mass and size at $z\sim1$ and 0; the main population of ETGs changes from disky to boxy at a critical stellar mass of $\log(M_{*}/M_{\odot})\sim11.5$ with the massive end dominated by boxy. The disky ETG fraction decreases with increasing stellar mass both at $z\sim1$ and $0$, and is consistent between these redshifts in all stellar mass bins when the Eddington bias is taken into account. Although uncertainties are large, the results suggest that the isophote shapes and probably dynamical properties of ETGs in massive clusters are already in place at $z>1$ and do not significantly evolve in $z<1$, despite significant size evolution in the same galaxy population. The constant disky fraction favors less violent processes than mergers as a main cause of the size and morphological evolution of intermediate mass ETGs in $z<1$.

We present an application of machine-learning (ML) techniques to source selection in the optical transient survey data with Hyper Suprime-Cam (HSC) on the Subaru telescope. Our goal is to select real transient events accurately and in a timely manner out of a large number of false candidates, obtained with the standard difference-imaging method. We have developed the transient selector which is based on majority voting of three ML machines of AUC Boosting, Random Forest, and Deep Neural Network. We applied it to our observing runs of Subaru-HSC in 2015 May and August, and proved it to be efficient in selecting optical transients. The false positive rate was 1.0% at the true positive rate of 90% in the magnitude range of 22.0--25.0 mag for the former data. For the latter run, we successfully detected and reported ten candidates of supernovae within the same day as the observation. From these runs, we learned the following lessons: (1) the training using artificial objects is effective in filtering out false candidates, especially for faint objects, and (2) combination of ML by majority voting is advantageous.

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

The Prime Focus Spectrograph (PFS) is a wide-field, multi-object spectrograph accommodating 2394 fibers to observe the sky at the prime focus of the Subaru telescope. The software system to operate a spectroscopic survey is structured by the four packages: Instrument control software, exposure targeting software, data reduction pipeline, and survey planning and tracking software. In addition, we operate a database system where various information such as properties of target objects, instrument configurations, and observation conditions is stored and is organized via a standardized data model for future references to update survey plans and to scientific researches. In this article, we present an overview of the software system and describe the workflows that need to be performed in the PFS operation, with some highlights on the database that organizes various information from sub-processes in the survey operation, and on the process of fiber configuration from the software perspectives.

The Baryon Oscillation Spectroscopic Survey (BOSS) has collected more than 150,000 $2.1 \leq z \leq 3.5$ quasar spectra since 2009. Using this unprecedented sample, we create a composite spectrum in the rest-frame of 102,150 quasar spectra from 800 \AA\ to 3300 \AA\ at a signal-to-noise ratio close to 1000 per pixel ($\Delta v$ of 69 km~s$^{-1}$). Included in this analysis is a correction to account for flux calibration residuals in the BOSS spectrophotometry. We determine the spectral index as a function of redshift of the full sample, warp the composite spectrum to match the median spectral index, and compare the resulting spectrum to SDSS photometry used in target selection. The quasar composite matches the color of the quasar population to within 0.02 magnitudes in $g-r$, 0.03 magnitudes in $r-i$, and 0.01 magnitudes in $i-z$ over the redshift range $2.2<z<2.6$. The composite spectrum deviates from the imaging photometry by 0.05 magnitudes around $z = 2.7$, likely due to differences in target selection as the quasar colors become similar to the stellar locus at this redshift. Finally, we characterize the line features in the high signal-to-noise composite and identify nine faint lines not found in the previous composite spectrum from SDSS.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

We introduce a method for identifying "twin" Type Ia supernovae, and using them to improve distance measurements. This novel approach to Type Ia supernova standardization is made possible by spectrophotometric time series observations from the Nearby Supernova Factory (SNfactory). We begin with a well-measured set of supernovae, find pairs whose spectra match well across the entire optical window, and then test whether this leads to a smaller dispersion in their absolute brightnesses. This analysis is completed in a blinded fashion, ensuring that decisions made in implementing the method do not inadvertently bias the result. We find that pairs of supernovae with more closely matched spectra indeed have reduced brightness dispersion. We are able to standardize this initial set of SNfactory supernovae to 0.083 +/- 0.012 magnitudes, implying a dispersion of 0.072 +/- 0.010 magnitudes in the absence of peculiar velocities. We estimate that with larger numbers of comparison SNe, e.g, using the final SNfactory spectrophotometric dataset as a reference, this method will be capable of standardizing high-redshift supernovae to within 0.06-0.07 magnitudes. These results imply that at least 3/4 of the variance in Hubble residuals in current supernova cosmology analyses is due to previously unaccounted-for astrophysical differences among the supernovae

Using moderate-resolution optical spectra from 58 background Lyman-break galaxies and quasars at $z\sim 2.3-3$ within a $11.5'\times13.5'$ area of the COSMOS field ($\sim 1200\,\mathrm{deg}^2$ projected area density or $\sim 2.4\,h^{-1}\,\mathrm{Mpc}$ mean transverse separation), we reconstruct a 3D tomographic map of the foreground Ly$\alpha$ forest absorption at $2.2<z<2.5$ with an effective smoothing scale of $\sigma_{3d}\approx3.5\,h^{-1}\,\mathrm{Mpc}$ comoving. Comparing with 61 coeval galaxies with spectroscopic redshifts in the same volume, we find that the galaxy positions are clearly biased towards regions with enhanced IGM absorption in the tomographic map. We find an extended IGM overdensity with deep absorption troughs at $z=2.45$ associated with a recently-discovered galaxy protocluster at the same redshift. Based on simulations matched to our data, we estimate the enclosed dark matter mass within this IGM overdensity to be $M_{\rm dm} (z=2.45) = (9\pm4)\times 10^{13}\,h^{-1}\,\mathrm{M_\odot}$, and argue based on this mass and absorption strength that it will form at least one $z\sim0$ galaxy cluster with $M(z=0) = (3\pm 2) \times 10^{14}\,h^{-1}\mathrm{M_\odot}$, although its elongated nature suggests that it will likely collapse into two separate clusters. We also point out a compact overdensity of six MOSDEF galaxies at $z=2.30$ within a $r\sim 1\,h^{-1}\,\mathrm{Mpc}$ radius and $\Delta z\sim 0.006$, which does not appear to have a large associated IGM overdensity. These results demonstrate the potential of Ly$\alpha$ forest tomography on larger volumes to study galaxy properties as a function of environment, as well as revealing the large-scale IGM overdensities associated with protoclusters and other features of large-scale structure.

By applying fractional calculus to the equation proposed by M. Planck in 1900, we obtain a new blackbody radiation law described by a Mittag-Leffler (ML) function. We have analyzed NASA COBE data by means of a non-extensive formula with a parameter $(q-1)$, a formula proposed by Ertik et al. with a fractional parameter $(\alpha-1)$, and our new formula including a parameter $(p-1)$, as well as the Bose-Einstein distribution with a dimensionless chemical potential $\mu$. It can be said that one role of the fractional parameter $(p-1)$ is almost the same as that of chemical potential $(\mu)$ as well as that of the parameter $(q-1)$ in the non-extensive approach.

NICMOS 2 observations are crucial for constraining distances to most of the existing sample of z > 1 SNe Ia. Unlike the conventional calibration programs, these observations involve long exposure times and low count rates. Reciprocity failure is known to exist in HgCdTe devices and a correction for this effect has already been implemented for high and medium count-rates. However observations at faint count-rates rely on extrapolations. Here instead, we provide a new zeropoint calibration directly applicable to faint sources. This is obtained via inter-calibration of NIC2 F110W/F160W with WFC3 in the low count-rate regime using z ~ 1 elliptical galaxies as tertiary calibrators. These objects have relatively simple near-IR SEDs, uniform colors, and their extended nature gives superior signal-to-noise at the same count rate than would stars. The use of extended objects also allows greater tolerances on PSF profiles. We find ST magnitude zeropoints (after the installation of the NICMOS cooling system, NCS) of 25.296 +- 0.022 for F110W and 25.803 +- 0.023 for F160W, both in agreement with the calibration extrapolated from count-rates 1,000 times larger (25.262 and 25.799). Before the installation of the NCS, we find 24.843 +- 0.025 for F110W and 25.498 +- 0.021 for F160W, also in agreement with the high-count-rate calibration (24.815 and 25.470). We also check the standard bandpasses of WFC3 and NICMOS 2 using a range of stars and galaxies at different colors and find mild tension for WFC3, limiting the accuracy of the zeropoints. To avoid human bias, our cross-calibration was "blinded" in that the fitted zeropoint differences were hidden until the analysis was finalized.

Previously we used the Nearby Supernova Factory sample to show that SNe~Ia having locally star-forming environments are dimmer than SNe~Ia having locally passive environments.Here we use the \constitution\ sample together with host galaxy data from \GALEX\ to independently confirm that result. The effect is seen using both the SALT2 and MLCS2k2 lightcurve fitting and standardization methods, with brightness differences of $0.094 \pm 0.037\ \mathrm{mag}$ for SALT2 and $0.155 \pm 0.041\ \mathrm{mag}$ for MLCS2k2 with $R_V=2.5$. When combined with our previous measurement the effect is $0.094 \pm 0.025\ \mathrm{mag}$ for SALT2. If the ratio of these local SN~Ia environments changes with redshift or sample selection, this can lead to a bias in cosmological measurements. We explore this issue further, using as an example the direct measurement of $H_0$. \GALEX{} observations show that the SNe~Ia having standardized absolute magnitudes calibrated via the Cepheid period--luminosity relation using {\textit{HST}} originate in predominately star-forming environments, whereas only ~50% of the Hubble-flow comparison sample have locally star-forming environments. As a consequence, the $H_0$ measurement using SNe~Ia is currently overestimated. Correcting for this bias, we find a value of $H_0^{corr}=70.6\pm 2.6\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$ when using the LMC distance, Milky Way parallaxes and the NGC~4258 megamaser as the Cepheid zeropoint, and $68.8\pm 3.3\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$ when only using NGC~4258. Our correction brings the direct measurement of $H_0$ within $\sim 1\,\sigma$ of recent indirect measurements based on the CMB power spectrum.

The Ly$\alpha$ forest transmission probability distribution function (PDF) is an established probe of the intergalactic medium (IGM) astrophysics, especially the temperature-density relationship of the IGM. We measure the transmission PDF from 3393 Baryon Oscillations Spectroscopic Survey (BOSS) quasars from SDSS Data Release 9, and compare with mock spectra that include careful modeling of the noise, continuum, and astrophysical uncertainties. The BOSS transmission PDFs, measured at $\langle z \rangle = [2.3,2.6,3.0]$, are compared with PDFs created from mock spectra drawn from a suite of hydrodynamical simulations that sample the IGM temperature-density relationship, $\gamma$, and temperature at mean-density, $T_0$, where $T(\Delta) = T_0 \Delta^{\gamma-1}$. We find that a significant population of partial Lyman-limit systems with a column-density distribution slope of $\beta_\mathrm{pLLS} \sim -2$ are required to explain the data at the low-transmission end of transmission PDF, while uncertainties in the mean Ly$\alpha$ forest transmission affect the high-transmission end. After modelling the LLSs and marginalizing over mean-transmission uncertainties, we find that $\gamma=1.6$ best describes the data over our entire redshift range, although constraints on $T_0$ are affected by systematic uncertainties. Within our model framework, isothermal or inverted temperature-density relationships ($\gamma \leq 1$) are disfavored at a significance of over 4$\sigma$, although this could be somewhat weakened by cosmological and astrophysical uncertainties that we did not model.

We develop a new framework for use in exploring Type Ia Supernova (SN Ia) spectra. Combining Principal Component Analysis (PCA) and Partial Least Square analysis (PLS) we are able to establish correlations between the Principal Components (PCs) and spectroscopic/photometric SNe Ia features. The technique was applied to ~120 supernova and ~800 spectra from the Nearby Supernova Factory. The ability of PCA to group together SNe Ia with similar spectral features, already explored in previous studies, is greatly enhanced by two important modifications: (1) the initial data matrix is built using derivatives of spectra over the wavelength, which increases the weight of weak lines and discards extinction, and (2) we extract time evolution information through the use of entire spectral sequences concatenated in each line of the input data matrix. These allow us to define a stable PC parameter space which can be used to characterize synthetic SN Ia spectra by means of real SN features. Using PLS, we demonstrate that the information from important previously known spectral indicators (namely the pseudo-equivalent width (pEW) of Si II 5972 / Si II 6355 and the line velocity of S II 5640 / Si II 6355) at a given epoch, is contained within the PC space and can be determined through a linear combination of the most important PCs. We also show that the PC space encompasses photometric features like B or V magnitudes, B-V color and SALT2 parameters c and x1. The observed colors and magnitudes, that are heavily affected by extinction, cannot be reconstructed using this technique alone. All the above mentioned applications allowed us to construct a metric space for comparing synthetic SN Ia spectra with observations.

We presented optical and near-infrared multi-band linear polarimetry of the highly reddened Type Ia SN~2014J appeared in M82. SN~2014J exhibits large polarization at shorter wavelengths, e.g., $4.8$\% in $B$ band, and the polarization decreases rapidly at longer wavelengths, with the position angle of the polarization remaining at approximately $40^{\circ}$ over the observed wavelength range. These polarimetric properties suggest that the observed polarization is likely to be caused predominantly by the interstellar dust within M82. Further analysis shows that the polarization peaks at a wavelengths much shorter than those obtained for the Galactic dust. The wavelength dependence of the polarization can be better described by an inverse power law rather than by Serkowski law for Galactic interstellar polarization. These suggests that the nature of the dust in M82 may be different from that in our Galaxy, with polarizing dust grains having a mean radius of $<0.1\ \mu$m.

We present the first observations of foreground Lyman-$\alpha$ forest absorption from high-redshift galaxies, targeting 24 star-forming galaxies (SFGs) with $z\sim 2.3-2.8$ within a $5' \times 15'$ region of the COSMOS field. The transverse sightline separation is $\sim 2\,h^{-1}\mathrm{Mpc}$ comoving, allowing us to create a tomographic reconstruction of the 3D Ly$\alpha$ forest absorption field over the redshift range $2.20\leq z\leq 2.45$. The resulting map covers $6\,h^{-1}\mathrm{Mpc} \times 14\,h^{-1}\mathrm{Mpc}$ in the transverse plane and $230\,h^{-1}\mathrm{Mpc}$ along the line-of-sight with a spatial resolution of $\approx 3.5\,h^{-1}\mathrm{Mpc}$, and is the first high-fidelity map of large-scale structure on $\sim\mathrm{Mpc}$ scales at $z>2$. Our map reveals significant structures with $\gtrsim 10\,h^{-1}\mathrm{Mpc}$ extent, including several spanning the entire transverse breadth, providing qualitative evidence for the filamentary structures predicted to exist in the high-redshift cosmic web. Simulated reconstructions with the same sightline sampling, spectral resolution, and signal-to-noise ratio recover the salient structures present in the underlying 3D absorption fields. Using data from other surveys, we identified 18 galaxies with known redshifts coeval with our map volume enabling a direct comparison to our tomographic map. This shows that galaxies preferentially occupy high-density regions, in qualitative agreement with the same comparison applied to simulations. Our results establishes the feasibility of the CLAMATO survey, which aims to obtain Ly$\alpha$ forest spectra for $\sim 1000$ SFGs over $\sim 1 \,\mathrm{deg}^2$ of the COSMOS field, in order to map out IGM large-scale structure at $\langle z \rangle \sim 2.3$ over a large volume $(100\,h^{-1}\mathrm{Mpc})^3$.

KA1858+4850 is a narrow-line Seyfert 1 galaxy at redshift 0.078 and is among the brightest active galaxies monitored by the Kepler mission. We have carried out a reverberation mapping campaign designed to measure the broad-line region size and estimate the mass of the black hole in this galaxy. We obtained 74 epochs of spectroscopic data using the Kast Spectrograph at the Lick 3-m telescope from February to November of 2012, and obtained complementary V-band images from five other ground-based telescopes. We measured the H-beta light curve lag with respect to the V-band continuum light curve using both cross-correlation techniques (CCF) and continuum light curve variability modeling with the JAVELIN method, and found rest-frame lags of lag_CCF = 13.53 (+2.03, -2.32) days and lag_JAVELIN = 13.15 (+1.08, -1.00) days. The H-beta root-mean-square line profile has a width of sigma_line = 770 +/- 49 km/s. Combining these two results and assuming a virial scale factor of f = 5.13, we obtained a virial estimate of M_BH = 8.06 (+1.59, -1.72) x 10^6 M_sun for the mass of the central black hole and an Eddington ratio of L/L_Edd ~ 0.2. We also obtained consistent but slightly shorter emission-line lags with respect to the Kepler light curve. Thanks to the Kepler mission, the light curve of KA1858+4850 has among the highest cadences and signal-to-noise ratios ever measured for an active galactic nucleus; thus, our black hole mass measurement will serve as a reference point for relations between black hole mass and continuum variability characteristics in active galactic nuclei.

This document presents the results from the Distances subgroup of the Cosmic Frontier Community Planning Study (Snowmass 2013). We summarize the current state of the field as well as future prospects and challenges. In addition to the established probes using Type IA supernovae and baryon acoustic oscillations, we also consider prospective methods based on clusters, active galactic nuclei, gravitational wave sirens and strong lensing time delays.

Using three magnified Type Ia supernovae (SNe Ia) detected behind CLASH clusters, we perform a first pilot study to see whether standardizable candles can be used to calibrate cluster mass maps created from strong lensing observations. Such calibrations will be crucial when next generation HST cluster surveys (e.g. FRONTIER) provide magnification maps that will, in turn, form the basis for the exploration of the high redshift Universe. We classify SNe using combined photometric and spectroscopic observations, finding two of the three to be clearly of type SN Ia and the third probable. The SNe exhibit significant amplification, up to a factor of 1.7 at $\sim5\sigma$ significance (SN-L2). We conducted this as a blind study to avoid fine tuning of parameters, finding a mean amplification difference between SNe and the cluster lensing models of $0.09 \pm 0.09^{stat} \pm 0.05^{sys}$ mag. This impressive agreement suggests no tension between cluster mass models and high redshift standardized SNe Ia. However, the measured statistical dispersion of $\sigma_{\mu}=0.21$ mag appeared large compared to the dispersion expected based on statistical uncertainties ($0.14$). Further work with the supernova and cluster lensing models, post unblinding, reduced the measured dispersion to $\sigma_{\mu}=0.12$. An explicit choice should thus be made as to whether SNe are used unblinded to improve the model, or blinded to test the model. As the lensed SN samples grow larger, this technique will allow improved constraints on assumptions regarding e.g. the structure of the dark matter halo.

We present measurements of the rates of high-redshift Type Ia supernovae derived from the Subaru/XMM-Newton Deep Survey (SXDS). We carried out repeat deep imaging observations with Suprime-Cam on the Subaru Telescope, and detected 1040 variable objects over 0.918 deg$^2$ in the Subaru/XMM-Newton Deep Field. From the imaging observations, light curves in the observed $i'$-band are constructed for all objects, and we fit the observed light curves with template light curves. Out of the 1040 variable objects detected by the SXDS, 39 objects over the redshift range $0.2 < z < 1.4$ are classified as Type Ia supernovae using the light curves. These are among the most distant SN Ia rate measurements to date. We find that the Type Ia supernova rate increase up to $z \sim 0.8$ and may then flatten at higher redshift. The rates can be fitted by a simple power law, $r_V(z)=r_0(1+z)^\alpha$ with $r_0=0.20^{+0.52}_{-0.16}$(stat.)$^{+0.26}_{-0.07}$(syst.)$\times 10^{-4} {\rm yr}^{-1}{\rm Mpc}^{-3}$, and $\alpha=2.04^{+1.84}_{-1.96}$(stat.)$^{+2.11}_{-0.86}$(syst.).

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

We report the discovery of a redshift 1.71 supernova in the GOODS North field. The Hubble Space Telescope (HST) ACS spectrum has almost negligible contamination from the host or neighboring galaxies. Although the rest frame sampled range is too blue to include any Si ii line, a principal component analysis allows us to confirm it as a Type Ia supernova with 92% confidence. A recent serendipitous archival HST WFC3 grism spectrum contributed a key element of the confirmation by giving a host-galaxy redshift of 1.713 +/- 0.007. In addition to being the most distant SN Ia with spectroscopic confirmation, this is the most distant Ia with a precision color measurement. We present the ACS WFC and NICMOS 2 photometry and ACS and WFC3 spectroscopy. Our derived supernova distance is in agreement with the prediction of LambdaCDM.

We present the BOSS Lyman-alpha (Lya) Forest Sample from SDSS Data Release 9, comprising 54,468 quasar spectra with zqso > 2.15 suitable for Lya forest analysis. This data set probes the intergalactic medium with absorption redshifts 2.0 < z_alpha < 5.7 over an area of 3275 square degrees, and encompasses an approximate comoving volume of 20 h^-3 Gpc^3. With each spectrum, we have included several products designed to aid in Lya forest analysis: improved sky masks that flag pixels where data may be unreliable, corrections for known biases in the pipeline estimated noise, masks for the cores of damped Lya systems and corrections for their wings, and estimates of the unabsorbed continua so that the observed flux can be converted to a fractional transmission. The continua are derived using a principal component fit to the quasar spectrum redwards of restframe Lya (lambda > 1216 Ang), extrapolated into the forest region and normalized by a linear function to fit the expected evolution of the Lya forest mean-flux. The estimated continuum errors are ~5% rms. We also discuss possible systematics arising from uncertain spectrophotometry and artifacts in the flux calibration; global corrections for the latter are provided. Our sample provides a convenient starting point for users to analyze clustering in BOSS Lya forest data, and it provides a fiducial data set that can be used to compare results from different analyses of baryon acoustic oscillations in the Lya forest. The full data set is available from the SDSS-III DR9 web site.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $\Lambda$CDM cosmology with H0 = 70 km/s/Mpc, $\Omega_{\rm M}$ = 0.3, and $\Omega_{\Lambda}$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

We present photometric recalibration of the Subaru Deep Field (SDF) and Subaru/XMM-Newton Deep Survey (SXDS). Recently, Yamanoi et al. (2012) suggested the existence of a discrepancy between the SDF and SXDS catalogs. We have used the Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) catalog and compared stars in common between SDF/SXDS and SDSS. We confirmed that there exists a 0.12 mag offset in B-band between the SDF and SXDS catalogs. Moreover, we found that significant zero point offsets in i-band (~ 0.10 mag) and z-band (~ 0.14 mag) need to be introduced to the SDF/SXDS catalogs to make it consistent with the SDSS catalog. We report the measured zero point offsets of five filter bands of SDF/SXDS catalogs. We studied the potential cause of these offsets, but the origins are yet to be understood.

There is a consensus that Type-Ia supernovae (SNe Ia) arise from the thermonuclear explosion of white dwarf stars that accrete matter from a binary companion. However, direct observation of SN Ia progenitors is lacking, and the precise nature of the binary companion remains uncertain. A temporal series of high-resolution optical spectra of the SN Ia PTF 11kx reveals a complex circumstellar environment that provides an unprecedentedly detailed view of the progenitor system. Multiple shells of circumsteller are detected and the SN ejecta are seen to interact with circumstellar material (CSM) starting 59 days after the explosion. These features are best described by a symbiotic nova progenitor, similar to RS Ophiuchi.

The nearby supernova SN 2011fe can be observed in unprecedented detail. Therefore, it is an important test case for Type Ia supernova (SN Ia) models, which may bring us closer to understanding the physical nature of these objects. Here, we explore how available and expected future observations of SN 2011fe can be used to constrain SN Ia explosion scenarios. We base our discussion on three-dimensional simulations of a delayed detonation in a Chandrasekhar-mass white dwarf and of a violent merger of two white dwarfs-realizations of explosion models appropriate for two of the most widely-discussed progenitor channels that may give rise to SNe Ia. Although both models have their shortcomings in reproducing details of the early and near-maximum spectra of SN 2011fe obtained by the Nearby Supernova Factory (SNfactory), the overall match with the observations is reasonable. The level of agreement is slightly better for the merger, in particular around maximum, but a clear preference for one model over the other is still not justified. Observations at late epochs, however, hold promise for discriminating the explosion scenarios in a straightforward way, as a nucleosynthesis effect leads to differences in the 55Co production. SN 2011fe is close enough to be followed sufficiently long to study this effect.

Using the sample of Type Ia supernovae (SNe Ia) discovered by the Hubble Space Telescope (HST) Cluster Supernova Survey and augmented with HST-observed SNe Ia in the GOODS fields, we search for correlations between the properties of SNe and their host galaxies at high redshift. We use galaxy color and quantitative morphology to determine the red sequence in 25 clusters and develop a model to distinguish passively evolving early-type galaxies from star-forming galaxies in both clusters and the field. With this approach, we identify six SN Ia hosts that are early-type cluster members and eleven SN Ia hosts that are early-type field galaxies. We confirm for the first time at z>0.9 that SNe Ia hosted by early-type galaxies brighten and fade more quickly than SNe Ia hosted by late-type galaxies. We also show that the two samples of hosts produce SNe Ia with similar color distributions. The relatively simple spectral energy distributions (SEDs) expected for passive galaxies enable us to measure stellar masses of early-type SN hosts. In combination with stellar mass estimates of late-type GOODS SN hosts from Thomson & Chary (2011), we investigate the correlation of host mass with Hubble residual observed at lower redshifts. Although the sample is small and the uncertainties are large, a hint of this relation is found at z>0.9. By simultaneously fitting the average cluster galaxy formation history and dust content to the red-sequence scatters, we show that the reddening of early-type cluster SN hosts is likely E(B-V) <~ 0.06. The similarity of the field and cluster early-type host samples suggests that field early-type galaxies that lie on the red sequence may also be minimally affected by dust. Hence, the early-type hosted SNe Ia studied here occupy a more favorable environment to use as well-characterized high-redshift standard candles than other SNe Ia.

We report a measurement of the Type Ia supernova (SN Ia) rate in galaxy clusters at 0.9 < z < 1.45 from the Hubble Space Telescope (HST) Cluster Supernova Survey. This is the first cluster SN Ia rate measurement with detected z > 0.9 SNe. Finding 8 +/- 1 cluster SNe Ia, we determine a SN Ia rate of 0.50 +0.23-0.19 (stat) +0.10-0.09 (sys) SNuB (SNuB = 10^-12 SNe L_{sun,B}^-1 yr^-1). In units of stellar mass, this translates to 0.36 +0.16-0.13 (stat) +0.07-0.06 (sys) SNuM (SNuM = 10^-12 SNe M_sun^-1 yr^-1). This represents a factor of approximately 5 +/- 2 increase over measurements of the cluster rate at z < 0.2. We parameterize the late-time SN Ia delay time distribution with a power law (proportional to t^s). Under the assumption of a cluster formation redshift of z_f = 3, our rate measurement in combination with lower-redshift cluster SN Ia rates constrains s = -1.41 +0.47/-0.40, consistent with measurements of the delay time distribution in the field. This measurement is generally consistent with expectations for the "double degenerate" scenario and inconsistent with some models for the "single degenerate" scenario predicting a steeper delay time distribution at large delay times. We check for environmental dependence and the influence of younger stellar populations by calculating the rate specifically in cluster red-sequence galaxies and in morphologically early-type galaxies, finding results similar to the full cluster rate. Finally, the upper limit of one host-less cluster SN Ia detected in the survey implies that the fraction of stars in the intra-cluster medium is less than 0.47 (95% confidence), consistent with measurements at lower redshifts.

The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of 4 in the U-band continuum and a factor of 2 in the H-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f=5.25, the virial estimate of the black hole mass is (3.2+-0.5)*10^7 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics.

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate out to z ~ 1.6 from the Hubble Space Telescope Cluster Supernova Survey. In observations spanning 189 orbits with the Advanced Camera for Surveys we discovered 29 SNe, of which approximately 20 are SNe Ia. Twelve of these SNe Ia are located in the foregrounds and backgrounds of the clusters targeted in the survey. Using these new data, we derive the volumetric SN Ia rate in four broad redshift bins, finding results consistent with previous measurements at z > 1 and strengthening the case for a SN Ia rate that is equal to or greater than ~0.6 x 10^-4/yr/Mpc^3 at z ~ 1 and flattening out at higher redshift. We provide SN candidates and efficiency calculations in a form that makes it easy to rebin and combine these results with other measurements for increased statistics. Finally, we compare the assumptions about host-galaxy dust extinction used in different high-redshift rate measurements, finding that different assumptions may induce significant systematic differences between measurements.

Type Ia supernovae (SNe Ia) have been used empirically as standardized candles to reveal the accelerating universe even though fundamental details, such as the nature of the progenitor system and how the star explodes, remained a mystery. There is consensus that a white dwarf star explodes after accreting matter in a binary system, but the secondary could be anything from a main sequence star to a red giant, or even another white dwarf. The uncertainty stems from the fact that no recent SN Ia has been discovered close enough to detect the stars before explosion. Here we report early observations of SN 2011fe (PTF11kly) in M101 at a distance of 6.4 Mpc, the closest SN Ia in the past 25 years. We find that the exploding star was likely a carbon-oxygen white dwarf, and from the lack of an early shock we conclude that the companion was most likely a main sequence star. Early spectroscopy shows high-velocity oxygen that varies on a time scale of hours and extensive mixing of newly synthesized intermediate mass elements in the outermost layers of the supernova. A companion paper uses pre-explosion images to rule out luminous red giants and most helium stars as companions.

Continuum fitting is an important aspect of Lyman-alpha forest science, since errors in the estimated optical depths scale with the fractional continuum error. However, traditional methods of estimating continua in noisy and moderate-resolution spectra (S/N < 10 pixel^-1 and R ~ 2000, respectively, such as SDSS) using power-law extrapolation or the mean spectrum, achieve no better than ~ 10-15% RMS accuracy. To improve on this, we introduce mean-flux regulated/principal component analysis (MF-PCA) continuum fitting. In this technique, PCA fitting is carried out redwards of the quasar Lyman-alpha line in order to provide a prediction for the shape of the Lyman-alpha forest continuum. The slope and amplitude of this continuum prediction is then corrected using external constraints for the Lyman-alpha forest mean-flux. From tests on mock spectra, we find that MF-PCA reduces the errors to 8% RMS in S/N ~ 2 spectra, and < 5% RMS in spectra with S/N > 5. The residual Fourier power in the continuum is decreased by a factor of a few in comparison with dividing by the mean continuum, enabling Lyman-alpha flux power spectrum measurements to be extended to ~2x larger scales. Using this new technique, we make available continuum fits for 12,069 z>2.3 Lyman-alpha forest spectra from SDSS DR7 for use by the community. This technique is also applicable to future releases of the ongoing BOSS survey, which is obtaining spectra for ~ 150,000 Lyman-alpha forest spectra at low signal-to-noise (S/N ~ 2).

PTF10ops is a Type Ia supernova (SN Ia), whose lightcurve and spectral properties place it outside the current SN Ia subtype classifications. Its spectra display the characteristic lines of subluminous SNe Ia, but it has a normal-width lightcurve with a long rise-time, typical of normal luminosity SNe Ia. The early-time optical spectra of PTF10ops were modelled using a spectral fitting code and found to have all the lines typically seen in subluminous SNe Ia, without the need to invoke more uncommon elements. The host galaxy environment of PTF10ops is also unusual with no galaxy detected at the position of the SN down to an absolute limiting magnitude of r \geq -12.0 mag, but a very massive galaxy is present at a separation of ~148 kpc and at the same redshift as suggested by the SN spectral features. The progenitor of PTF10ops is most likely a very old star, possibly in a low metallicity environment, which affects its explosion mechanism and observational characteristics. PTF10ops does not easily fit into any of the current models of either subluminous or normal SN Ia progenitor channels.

We present observational constraints on the nature of dark energy using the Supernova Legacy Survey three year sample (SNLS3) of Guy et al. (2010) and Conley et al. (2011). We use the 472 SNe Ia in this sample, accounting for recently discovered correlations between SN Ia luminosity and host galaxy properties, and include the effects of all identified systematic uncertainties directly in the cosmological fits. Combining the SNLS3 data with the full WMAP7 power spectrum, the Sloan Digital Sky Survey luminous red galaxy power spectrum, and a prior on the Hubble constant H0 from SHOES, in a flat universe we find omega_m=0.269+/-0.015 and w=-1.061+0.069-0.068 -- a 6.5% measure of the dark energy equation-of-state parameter w. The statistical and systematic uncertainties are approximately equal, with the systematic uncertainties dominated by the photometric calibration of the SN Ia fluxes -- without these calibration effects, systematics contribute only a ~2% error in w. When relaxing the assumption of flatness, we find omega_m=0.271+/-0.015, omega_k=-0.002+/-0.006, and w=-1.069+0.091-0.092. Parameterizing the time evolution of w as w(a)=w_0+w_a(1-a), gives w_0=-0.905+/-0.196, w_a=-0.984+1.094-1.097 in a flat universe. All of our results are consistent with a flat, w=-1 universe. The size of the SNLS3 sample allows various tests to be performed with the SNe segregated according to their light curve and host galaxy properties. We find that the cosmological constraints derived from these different sub-samples are consistent. There is evidence that the coefficient, beta, relating SN Ia luminosity and color, varies with host parameters at >4sigma significance (in addition to the known SN luminosity--host relation); however this has only a small effect on the cosmological results and is currently a sub-dominant systematic.

We present weak gravitational lensing analysis of 22 high-redshift (z >~1) clusters based on Hubble Space Telescope images. Most clusters in our sample provide significant lensing signals and are well detected in their reconstructed two-dimensional mass maps. Combining the current results and our previous weak-lensing studies of five other high-z clusters, we compare gravitational lensing masses of these clusters with other observables. We revisit the question whether the presence of the most massive clusters in our sample is in tension with the current LambdaCDM structure formation paradigm. We find that the lensing masses are tightly correlated with the gas temperatures and establish, for the first time, the lensing mass-temperature relation at z >~ 1. For the power law slope of the M-TX relation (M propto T^{\alpha}), we obtain \alpha=1.54 +/- 0.23. This is consistent with the theoretical self-similar prediction \alpha=3/2 and with the results previously reported in the literature for much lower redshift samples. However, our normalization is lower than the previous results by 20-30%, indicating that the normalization in the M-TX relation might evolve. After correcting for Eddington bias and updating the discovery area with a more conservative choice, we find that the existence of the most massive clusters in our sample still provides a tension with the current Lambda CDM model. The combined probability of finding the four most massive clusters in this sample after marginalization over current cosmological parameters is less than 1%.

We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

The SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), a five-year spectroscopic survey of 10,000 deg^2, achieved first light in late 2009. One of the key goals of BOSS is to measure the signature of baryon acoustic oscillations in the distribution of Ly-alpha absorption from the spectra of a sample of ~150,000 z>2.2 quasars. Along with measuring the angular diameter distance at z\approx2.5, BOSS will provide the first direct measurement of the expansion rate of the Universe at z > 2. One of the biggest challenges in achieving this goal is an efficient target selection algorithm for quasars over 2.2 < z < 3.5, where their colors overlap those of stars. During the first year of the BOSS survey, quasar target selection methods were developed and tested to meet the requirement of delivering at least 15 quasars deg^-2 in this redshift range, out of 40 targets deg^-2. To achieve these surface densities, the magnitude limit of the quasar targets was set at g <= 22.0 or r<=21.85. While detection of the BAO signature in the Ly-alpha absorption in quasar spectra does not require a uniform target selection, many other astrophysical studies do. We therefore defined a uniformly-selected subsample of 20 targets deg^-2, for which the selection efficiency is just over 50%. This "CORE" subsample will be fixed for Years Two through Five of the survey. In this paper we describe the evolution and implementation of the BOSS quasar target selection algorithms during the first two years of BOSS operations. We analyze the spectra obtained during the first year. 11,263 new z>2.2 quasars were spectroscopically confirmed by BOSS. Our current algorithms select an average of 15 z > 2.2 quasars deg^-2 from 40 targets deg^-2 using single-epoch SDSS imaging. Multi-epoch optical data and data at other wavelengths can further improve the efficiency and completeness of BOSS quasar target selection. [Abridged]

We combine high redshift Type Ia supernovae from the first 3 years of the Supernova Legacy Survey (SNLS) with other supernova (SN) samples, primarily at lower redshifts, to form a high-quality joint sample of 472 SNe (123 low-$z$, 93 SDSS, 242 SNLS, and 14 {\it Hubble Space Telescope}). SN data alone require cosmic acceleration at >99.9% confidence, including systematic effects. For the dark energy equation of state parameter (assumed constant out to at least $z=1.4$) in a flat universe, we find $w = -0.91^{+0.16}_{-0.20}(\mathrm{stat}) ^{+0.07}_{-0.14} (\mathrm{sys})$ from SNe only, consistent with a cosmological constant. Our fits include a correction for the recently discovered relationship between host-galaxy mass and SN absolute brightness. We pay particular attention to systematic uncertainties, characterizing them using a systematics covariance matrix that incorporates the redshift dependence of these effects, as well as the shape-luminosity and color-luminosity relationships. Unlike previous work, we include the effects of systematic terms on the empirical light-curve models. The total systematic uncertainty is dominated by calibration terms. We describe how the systematic uncertainties can be reduced with soon to be available improved nearby and intermediate-redshift samples, particularly those calibrated onto USNO/SDSS-like systems.

The Seyfert 1 galaxy Zw 229-015 is among the brightest active galaxies being monitored by the Kepler mission. In order to determine the black hole mass in Zw 229-015 from H-beta reverberation mapping, we have carried out nightly observations with the Kast Spectrograph at the Lick 3m telescope during the dark runs from June through December 2010, obtaining 54 spectroscopic observations in total. We have also obtained nightly V-band imaging with the Katzman Automatic Imaging Telescope at Lick Observatory and with the 0.9m telescope at the Brigham Young University West Mountain Observatory over the same period. We detect strong variability in the source, which exhibited more than a factor of 2 change in broad H-beta flux. From cross-correlation measurements, we find that the H-beta light curve has a rest-frame lag of 3.86(+0.69,-0.90) days with respect to the V-band continuum variations. We also measure reverberation lags for H-alpha and H-gamma and find an upper limit to the H-delta lag. Combining the H-beta lag measurement with a broad H-beta width of sigma = 1590+/-47 km/s measured from the root-mean-square variability spectrum, we obtain a virial estimate of M_BH = 1.00(-0.24,+0.19)*10^7 solar masses for the black hole in Zw 229-015. As a Kepler target, Zw 229-015 will eventually have one of the highest-quality optical light curves ever measured for any active galaxy, and the black hole mass determined from reverberation mapping will serve as a benchmark for testing relationships between black hole mass and continuum variability characteristics in active galactic nuclei.

We present photometric properties and distance measurements of 252 high redshift Type Ia supernovae (0.15 < z < 1.1) discovered during the first three years of the Supernova Legacy Survey (SNLS). These events were detected and their multi-colour light curves measured using the MegaPrime/MegaCam instrument at the Canada-France-Hawaii Telescope (CFHT), by repeatedly imaging four one-square degree fields in four bands. Follow-up spectroscopy was performed at the VLT, Gemini and Keck telescopes to confirm the nature of the supernovae and to measure their redshifts. Systematic uncertainties arising from light curve modeling are studied, making use of two techniques to derive the peak magnitude, shape and colour of the supernovae, and taking advantage of a precise calibration of the SNLS fields. A flat LambdaCDM cosmological fit to 231 SNLS high redshift Type Ia supernovae alone gives Omega_M = 0.211 +/- 0.034(stat) +/- 0.069(sys). The dominant systematic uncertainty comes from uncertainties in the photometric calibration. Systematic uncertainties from light curve fitters come next with a total contribution of +/- 0.026 on Omega_M. No clear evidence is found for a possible evolution of the slope (beta) of the colour-luminosity relation with redshift.

We report on work to increase the number of well-measured Type Ia supernovae (SNe Ia) at high redshifts. Light curves, including high signal-to-noise HST data, and spectra of six SNe Ia that were discovered during 2001 are presented. Additionally, for the two SNe with z>1, we present ground-based J-band photometry from Gemini and the VLT. These are among the most distant SNe Ia for which ground based near-IR observations have been obtained. We add these six SNe Ia together with other data sets that have recently become available in the literature to the Union compilation (Kowalski et al. 2008). We have made a number of refinements to the Union analysis chain, the most important ones being the refitting of all light curves with the SALT2 fitter and an improved handling of systematic errors. We call this new compilation, consisting of 557 supernovae, the Union2 compilation. The flat concordance LambdaCDM model remains an excellent fit to the Union2 data with the best fit constant equation of state parameter w=-0.997^{+0.050}_{-0.054} (stat) ^{+0.077}_{-0.082} (stat+sys\ together) for a flat universe, or w=-1.035^{+0.055}_{-0.059} (stat)^{+0.093}_{-0.097} (stat+sys together) with curvature. We also present improved constraints on w(z). While no significant change in w with redshift is detected, there is still considerable room for evolution in w. The strength of the constraints depend strongly on redshift. In particular, at z > 1, the existence and nature of dark energy are only weakly constrained by the data.

We have combined the large SN Ia database of the Canada-France-Hawaii Telescope Supernova Legacy Survey and catalogs of galaxies with photometric redshifts, VLA 1.4 GHz radio sources, and Spitzer infrared sources. We present eight SNe Ia in early-type host galaxies which have counterparts in the radio and infrared source catalogs. We find the SN Ia rate in subsets of radio and infrared early-type galaxies is ~1-5 times the rate in all early-type galaxies, and that any enhancement is always <~ 2 sigma. Rates in these subsets are consistent with predictions of the two component "A+B" SN Ia rate model. Since infrared properties of radio SN Ia hosts indicate dust obscured star formation, we incorporate infrared star formation rates into the "A+B" model. We also show the properties of SNe Ia in radio and infrared galaxies suggest the hosts contain dust and support a continuum of delay time distributions for SNe Ia, although other delay time distributions cannot be ruled out based on our data.

We present spectra of high-redshift supernovae (SNe) that were taken with the Subaru low resolution optical spectrograph, FOCAS. These SNe were found in SN surveys with Suprime-Cam on Subaru, the CFH12k camera on the Canada-France-Hawaii Telescope (CFHT), and the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). These SN surveys specifically targeted z>1 Type Ia supernovae (SNe Ia). From the spectra of 39 candidates, we obtain redshifts for 32 candidates and spectroscopically identify 7 active candidates as probable SNe Ia, including one at z=1.35, which is the most distant SN Ia to be spectroscopically confirmed with a ground-based telescope. An additional 4 candidates are identified as likely SNe Ia from the spectrophotometric properties of their host galaxies. Seven candidates are not SNe Ia, either being SNe of another type or active galactic nuclei. When SNe Ia are observed within a week of maximum light, we find that we can spectroscopically identify most of them up to z=1.1. Beyond this redshift, very few candidates were spectroscopically identified as SNe Ia. The current generation of super red-sensitive, fringe-free CCDs will push this redshift limit higher.

We report the discovery of a multiply lensed Ly Alpha (Lya) emitter at z = 3.90 behind the massive galaxy cluster WARPS J1415.1+3612 at z = 1.026. Images taken by the Hubble Space Telescope(HST) using ACS reveal a complex lensing system that produces a prominent, highly magnified arc and a triplet of smaller arcs grouped tightly around a spectroscopically confirmed cluster member. Spectroscopic observations using FOCAS on Subaru confirm strong Lya emission in the source galaxy and provide redshifts for more than 21 cluster members, from which we obtain a velocity dispersion of 807+/-185 km/s. Assuming a singular isothermal sphere profile, the mass within the Einstein ring (7.13+/-0.38") corresponds to a central velocity dispersion of 686+15-19 km/s for the cluster, consistent with the value estimated from cluster member redshifts. Our mass profile estimate from combining strong lensing and dynamical analyses is in good agreement with both X-ray and weak lensing results.

We present data from the Supernova Cosmology Project for five high redshift Type Ia supernovae (SNe Ia) that were obtained using the NICMOS infrared camera on the Hubble Space Telescope. We add two SNe from this sample to a rest-frame I-band Hubble diagram, doubling the number of high redshift supernovae on this diagram. This I-band Hubble diagram is consistent with a flat universe (Omega_Matter, Omega_Lambda= 0.29, 0.71). A homogeneous distribution of large grain dust in the intergalactic medium (replenishing dust) is incompatible with the data and is excluded at the 5 sigma confidence level, if the SN host galaxy reddening is corrected assuming R_V=1.75. We use both optical and infrared observations to compare photometric properties of distant SNe Ia with those of nearby objects. We find generally good agreement with the expected color evolution for all SNe except the highest redshift SN in our sample (SN 1997ek at z=0.863) which shows a peculiar color behavior. We also present spectra obtained from ground based telescopes for type identification and determination of redshift.

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

Using deep two-band imaging from the Hubble Space Telescope, we measure the color-magnitude relations (CMR) of E/S0 galaxies in a set of 9 optically-selected clusters principally from the Red-Sequence Cluster Survey (RCS) at 0.9 < z < 1.23. We find that the mean scatter in the CMR in the observed frame of this set of clusters is 0.049 +/- 0.008, as compared to 0.031 +/- 0.007 in a similarly imaged and identically analyzed X-ray sample at similar redshifts. Single-burst stellar population models of the CMR scatter suggest that the E/S0 population in these RCS clusters truncated their star-formation at z~1.6, some 0.9 Gyrs later than their X-ray E/S0 counterparts which were truncated at z~2.1. The notion that this is a manifestation of the differing evolutionary states of the two populations of cluster galaxies is supported by comparison of the fraction of bulge-dominated galaxies found in the two samples which shows that optically-selected clusters contain a smaller fraction of E/S0 galaxies at the their cores

We present a study of the morphological fractions and color-magnitude relation in the most distant X-ray selected galaxy cluster currently known, XMMXCS J2215.9-1738 at z=1.46, using a combination of optical imaging data obtained with the Hubble Space Telescope Advanced Camera for Surveys, and infrared data from the Multi-Object Infrared Camera and Spectrograph, mounted on the 8.2m Subaru telescope. We find that the morphological mix of the cluster galaxy population is similar to clusters at z~1: approximately ~62% of the galaxies identified as likely cluster members are ellipticals or S0s; and ~38% are spirals or irregulars. We measure the color-magnitude relations for the early type galaxies, finding that the slope in the z_850-J relation is consistent with that measured in the Coma cluster, some ~9 Gyr earlier, although the uncertainty is large. In contrast, the measured intrinsic scatter about the color-magnitude relation is more than three times the value measured in Coma, after conversion to rest frame U-V. From comparison with stellar population synthesis models, the intrinsic scatter measurements imply mean luminosity weighted ages for the early type galaxies in J2215.9-1738 of ~3 Gyr, corresponding to the major epoch of star formation coming to an end at z_f = 3-5. We find that the cluster exhibits evidence of the `downsizing' phenomenon: the fraction of faint cluster members on the red sequence expressed using the Dwarf-to-Giant Ratio (DGR) is 0.32+/-0.18. This is consistent with extrapolation of the redshift evolution of the DGR seen in cluster samples at z < 1. In contrast to observations of some other z > 1 clusters, we find a lack of very bright galaxies within the cluster.

Timmes, Brown & Truran found that metallicity variations could theoretically account for a 25% variation in the mass of 56Ni synthesized in Type Ia supernovae (SNe Ia), and thus account for a large fraction of the scatter in observed SN Ia luminosities. Higher-metallicity progenitors are more neutron-rich, producing more stable burning products relative to radioactive 56Ni. We develop a new method for estimating bolometric luminosity and 56Ni yield in SNe Ia and use it to test the theory with data from the Supernova Legacy Survey. We find that the average 56Ni yield does drop in SNe Ia from high metallicity environments, but the theory can only account for 7%--10% of the dispersion in SN Ia 56Ni mass, and thus luminosity. This is because the effect is dominant at metallicities significantly above solar, whereas we find that SN hosts have predominantly subsolar or only moderately above-solar metallicities. We also show that allowing for changes in O/Fe with the metallicity [Fe/H] does not have a major effect on the theoretical prediction of Timmes, Brown & Truran, so long as one is using the O/H as the independent variable. Age may have a greater effect than metallicity -- we find that the luminosity weighted age of the host galaxy is correlated with 56Ni yield, and thus more massive progenitors give rise to more luminous explosions. This is hard to understand if most SNe Ia explode when the primaries reach the Chandrasekhar mass. Finally, we test the findings of Gallagher et al., that the residuals of SNe Ia from the Hubble diagram are correlated with host galaxy metallicity, and we find no such correlation.

We present observations of SCP 06F6, an unusual optical transient discovered during the Hubble Space Telescope Cluster Supernova Survey. The transient brightened over a period of ~100 days, reached a peak magnitude of ~21.0 in both i_775 and z_850, and then declined over a similar timescale. There is no host galaxy or progenitor star detected at the location of the transient to a 3 sigma upper limit of i_775 = 26.4 and z_850 = 26.1, giving a corresponding lower limit on the flux increase of a factor of ~120. Multiple spectra show five broad absorption bands between 4100 AA and 6500 AA and a mostly featureless continuum longward of 6500 AA. The shape of the lightcurve is inconsistent with microlensing. The transient's spectrum, in addition to being inconsistent with all known supernova types, is not matched to any spectrum in the Sloan Digital Sky Survey (SDSS) database. We suggest that the transient may be one of a new class.

The recent robust and homogeneous analysis of the world's supernova distance-redshift data, together with cosmic microwave background and baryon acoustic oscillation data, provides a powerful tool for constraining cosmological models. Here we examine particular classes of scalar field, modified gravity, and phenomenological models to assess whether they are consistent with observations even when their behavior deviates from the cosmological constant Lambda. Some models have tension with the data, while others survive only by approaching the cosmological constant, and a couple are statistically favored over LCDM. Dark energy described by two equation of state parameters has considerable phase space to avoid Lambda and next generation data will be required to constrain such physics.

We present a new compilation of Type Ia supernovae (SNe Ia), a new dataset of low-redshift nearby-Hubble-flow SNe and new analysis procedures to work with these heterogeneous compilations. This ``Union'' compilation of 414 SN Ia, which reduces to 307 SNe after selection cuts, includes the recent large samples of SNe Ia from the Supernova Legacy Survey and ESSENCE Survey, the older datasets, as well as the recently extended dataset of distant supernovae observed with HST. A single, consistent and blind analysis procedure is used for all the various SN Ia subsamples, and a new procedure is implemented that consistently weights the heterogeneous data sets and rejects outliers. We present the latest results from this Union compilation and discuss the cosmological constraints from this new compilation and its combination with other cosmological measurements (CMB and BAO). The constraint we obtain from supernovae on the dark energy density is $\Omega_\Lambda= 0.713^{+0.027}_{-0.029} (stat)}^{+0.036}_{-0.039} (sys)}$, for a flat, LCDM Universe. Assuming a constant equation of state parameter, $w$, the combined constraints from SNe, BAO and CMB give $w=-0.969^{+0.059}_{-0.063}(stat)^{+0.063}_{-0.066} (sys)$. While our results are consistent with a cosmological constant, we obtain only relatively weak constraints on a $w$ that varies with redshift. In particular, the current SN data do not yet significantly constrain $w$ at $z>1$. With the addition of our new nearby Hubble-flow SNe Ia, these resulting cosmological constraints are currently the tightest available.

The Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS) has created a large homogeneous database of intermediate redshift (0.2 < z < 1.0) type Ia supernovae (SNe Ia). The SNLS team has shown that correlations exist between SN Ia rates, properties, and host galaxy star formation rates. The SNLS SN Ia database has now been combined with a photometric redshift galaxy catalog and an optical galaxy cluster catalog to investigate the possible influence of galaxy clustering on the SN Ia rate, over and above the expected effect due to the dependence of SFR on clustering through the morphology-density relation. We identify three cluster SNe Ia, plus three additional possible cluster SNe Ia, and find the SN Ia rate per unit mass in clusters at intermediate redshifts is consistent with the rate per unit mass in field early-type galaxies and the SN Ia cluster rate from low redshift cluster targeted surveys. We also find the number of SNe Ia in cluster environments to be within a factor of two of expectations from the two component SNIa rate model.

We present a new measurement of the volumetric rate of Type Ia supernova up to a redshift of 1.7, using the Hubble Space Telescope (HST) GOODS data combined with an additional HST dataset covering the North GOODS field collected in 2004. We employ a novel technique that does not require spectroscopic data for identifying Type Ia supernovae (although spectroscopic measurements of redshifts are used for over half the sample); instead we employ a Bayesian approach using only photometric data to calculate the probability that an object is a Type Ia supernova. This Bayesian technique can easily be modified to incorporate improved priors on supernova properties, and it is well-suited for future high-statistics supernovae searches in which spectroscopic follow up of all candidates will be impractical. Here, the method is validated on both ground- and space-based supernova data having some spectroscopic follow up. We combine our volumetric rate measurements with low redshift supernova data, and fit to a number of possible models for the evolution of the Type Ia supernova rate as a function of redshift. The data do not distinguish between a flat rate at redshift > 0.5 and a previously proposed model, in which the Type Ia rate peaks at redshift >1 due to a significant delay from star-formation to the supernova explosion. Except for the highest redshifts, where the signal to noise ratio is generally too low to apply this technique, this approach yields smaller or comparable uncertainties than previous work.

We list and analyze H I absorption lines at redshifts 2 < z < 4 with column density (12 < log(N_HI) < 19) in 40 high-resolutional (FWHM = 8.0 km/s) quasar spectra obtained with the Keck+HIRES. We de-blend and fit all H I lines within 1,000 km/s of 86 strong H I lines whose column densities are log(N_HI/[cm^-2]) > 15. Unlike most prior studies, we use not only Lya but also all visible higher Lyman series lines to improve the fitting accuracy. This reveals components near to higher column density systems that can not be seen in Lya. We list the Voigt profile fits to the 1339 H I components that we found. We examined physical properties of H I lines after separating them into several sub-samples according to their velocity separation from the quasars, their redshift, column density and the S/N ratio of the spectrum. We found two interesting trends for lines with 12 < log(N_HI) < 15 which are within 200-1000 km/s of systems with log(N_HI) > 15. First, their column density distribution becomes steeper, meaning relatively fewer high column density lines, at z < 2.9. Second, their column density distribution also becomes steeper and their line width becomes broader by about 2-3 km/s when they are within 5,000 km/s of their quasar.

We use Keck/HIRES spectra of 37 optically bright quasars at z=2-4 to study narrow absorption lines that are intrinsic to the quasars (intrinsic NALs, produced in gas that is physically associated with the quasar central engine). We identify 150 NAL systems, that contain 124 C IV, 12 N V, and 50 Si IV doublets, of which 18 are associated systems (within 5,000 km/s of the quasar redshift). We use partial coverage analysis to separate intrinsic NALs from NALs produced in cosmologically intervening structures. We find 39 candidate intrinsic systems, (28 reliable determinations and 11 that are possibly intrinsic). We estimate that 10-17% of C IV systems at blueshifts of 5,000-70,000 km/s relative to quasars are intrinsic. At least 32% of quasars contain one or more intrinsic C IV NALs. Considering N V and Si IV doublets showing partial coverage as well, at least 50% of quasars host intrinsic NALs. This result constrains the solid angle subtended by the absorbers to the background source(s). We identify two families of intrinsic NAL systems, those with strong N V absorption, and those with negligible absorption in N V, but with partial coverage in the C IV doublet. We discuss the idea that these two families represent different regions or conditions in accretion disk winds. Of the 26 intrinsic C IV NAL systems, 13 have detectable low-ionization absorption lines at similar velocities, suggesting that these are two-phase structures in the wind rather than absorbers in the host galaxy. We also compare possible models for quasar outflows, including radiatively accelerated disk-driven winds, magnetocentrifugally accelerated winds, and pressure-driven winds, and we discuss ways of distinguishing between these models observationally.

We present statistical and monitoring results of narrow absorption lines that are physically related to quasars (i.e., intrinsic NALs). We use Keck/HIRES spectra of 37 optically bright quasars at z=2-4, and identify 150 NAL systems that contain 124 C IV, 12 N V, and 50 Si IV doublets. Among them, 39 are classified as intrinsic systems based on partial coverage analysis. At least 50% of quasars host intrinsic NALs. We identify two families of intrinsic systems based on their ionization state. Some intrinsic systems have detectable low-ionization NALs at similar velocities as higher-ionization NALs, although such low-ionization lines are rare in broad absorption line (BAL) systems. We also have observed an optically bright quasar, HS1603+3820, eight times with Subaru/HDS and HET/MRS over an interval of 4.2 years (1.2 years in the quasar rest frame), for the purpose of monitoring a variable C IV mini-BAL system. We find that all the troughs of the system vary in concert. However, no other correlations are seen between the variations of different profile parameters. We propose that the observed variations are either (i) a result of rapid continuum fluctuations, caused by a clumpy screen of variable optical depth located between the continuum source and the mini-BAL gas, or (ii) a result of variable scattering of continuum photons around the absorber.

Powerful outflows from star-forming galaxies are expected to push away the neutral intergalactic medium (IGM) around those galaxies, and produce absorption-free gaps in the Ly-alpha forest. We analyze the abundance of gaps of various sizes in three high resolution spectra of quasars at z ~ 3 - 3.5. The gap statistics agrees well with a model in which galactic halos above a minimum mass scale of M_min ~ 10^10 M_sun produce bubbles with a characteristic radius of R_b ~ 0.48 Mpc/h. Both numbers are consistent with naive theoretical expectations, where the minimum galaxy mass reflects the threshold for infall of gas out of a photo-ionized IGM. The observed gaps are typically bounded by deep absorption features as expected from the accumulation of swept-up gas on the bubble walls.

We use high quality echelle spectra of 24 QSOs to provide a calibrated measurement of the total amount of Lyman-alpha forest absorption (DA) over the redshift range 2.2 < z < 3.2. Our measurement of DA excludes absorption from metal lines or the Lyman-alpha lines of Lyman limit systems and damped Lyman-alpha systems. We use artificial spectra with realistic flux calibration errors to show that we are able to place continuum levels that are accurate to better than 1%. When we combine our results with our previous results between 1.6 < z < 2.2, we find that the redshift evolution of DA is well described over 1.6 < z < 3.2 as A*(1+z)^gamma, where A = 0.0062 and gamma = 2.75. We detect no significant deviations from a smooth power law evolution over the redshift range studied. We find less H I absorption than expected at z=3, implying that the UV background is about 40% higher than expected. Our data appears to be consistent with an H I ionization rate of Gamma ~ 1.4E-12 1/s.

We investigate the variety in quasar UV spectra (1020-1600A) with emphasis on the weak emission lines in the Ly alpha forest region using principal component analysis (PCA). This paper is a continuation of Suzuki et al. (2005, Paper I), but with a different approach. We use 50 smooth continuum fitted quasar spectra (0.14 < z < 1.04) taken by the Hubble Space Telescope (HST) Faint Object Spectrograph. Three weak emission lines in the Ly alpha forest are identified as FeII 1070.95A, FeII+FeIII 1123.17A, and CIII* 1175.88A. Using first two standardized PCS coefficients, we introduce five classifications: Class Zero and Class I-IV. These classifications will guide us in finding the continuum level in the Ly alpha forest, and we discuss the characteristics of each class. We show weak emission lines in the Ly alpha forest become eminent for Class III and IV. By actively using PCS, we can generate artificial quasar spectra which are useful to test the detection of quasars, DLAs, and the continuum calibration. We provide 10,000 artificially generated spectra. We show that the power-law extrapolated continuum is inadequate to perform precise measurements of the mean flux in the Ly alpha forest because of the weak emission lines and the extended tails of Ly alpha and Ly beta/OVI emission lines. We introduce a correction factor dF such that the true mean flux <F> can be related to <F_{Power-Law}> that is measured by using power-law continuum extrapolation: <F> = <F_{Power-Law}> dF. The correction factor dF ranges from 0.84 to 1.05 with a mean of 0.947 and a standard deviation of 0.031 for our 50 quasars. This result means that we miss 5.3% of flux on average and we show that there are cases where we would miss 16% of flux using a power-law extrapolation.

We present 40 fully hydrodynamical numerical simulations of the intergalactic gas that gives rise to the Ly-alpha forest. We make artificial spectra from each and measure three output parameters: the mean flux, a measure of the most common Ly-alpha line width (b-value), and the 1D power spectrum of the flux on scales from 0.01 - 0.1 s/km. We investigate how these output parameters are correlated with the input parameters that describe the simulations. These input parameters include the main cosmological parameters; and two astrophysical parameters gammah and gammahe. The parameter gammah controls the rate of ionization of H I, He I and He II while gammahe controls the rate of heating. We compare the output from the simulations to data from HIRES and UVES spectra and we give a new measurement of the flux power from these spectra at z=1.95. We find that simulations with a wide variety of sigma_8 values, from at least 0.8 - 1.1, can fit the small scale flux power and b-values when we adjust gammahe to compensate for the sigma_8 change. We can also use gammah to adjust the H I ionization rate to simultaneously match the mean flux. We derive scaling relations that give the output parameter values expected for a variety of input parameters with errors of 2% for optical depth and 1.4% for b-balues. We can break the degeneracy between gammahe and sigma_8 using large scale power or other data to fix sigma_8.

Using H I absorption alone, we attempt to separate H I absorption lines in quasar spectra into two categories; HDLs (Higher Density Lines) and LDLs (Lower Density Lines), and we discuss the difference in their physical properties. We deblend and fit all H I lines with Voigt profiles, and make an unbiased sample of H I lines covering a wide column density range (12 < log N_HI < 19 cm^-2). To reduce the influence of line blending, we simultaneously fit several Lyman series lines. As a result of a two-point correlation analysis, we found that higher column density H I lines are clustering at $\Delta$ v < 200 km/s, while lower ones at $\Delta$ v < 100 km/s. We define HDLs as H I lines with 15 < log N_HI < 19 cm^-2 and all H I lines within +/- 200 km/s of a line with log N_HI > 15 cm^-2, and LDLs as others with 12 < log N_HI < 15 cm^-2. We found that the HDLs have smaller minimum b-values for a given column density than the LDLs. This difference is successfully reproduced by our Hydrodynamic simulation. The LDLs seem to be cool or shock-heated diffuse IGM gas, while the HDLs are likely to be cooler dense gas near to galaxies.

We present a moderate resolution (~1.15 Angstroms/pixel) survey of 79 quasars obtained using the Kast spectrograph on the Shane 3m telescope at Lick observatory. The spectra span the wavelength range of 3175-5880 Angstroms, and have typical signal to noise of 6-20 in the regions of the spectra showing Lyman alpha forest absorption. The quasars have a mean emission redshift of z=2.17, and nearly all cover the entire Lyman alpha forest between Lyman alpha and Lyman beta. Although the quasars were selected to avoid BAL, two quasars in the survey are BAL, one of which is a new discovery. We list the HI and metal ions observed in a total of 140 absorption systems. We also identify 526 emission lines, and list their observed wavelengths, along with new redshifts of the quasars. We determine the rest wavelengths of 3 emission lines or line blends in the forest to be 1070.95 +/- 1.00, 1123.13 +/- 0.51, and 1175.88 +/- 0.30 Angstroms.

We identify a concordant model for the intergalactic medium (IGM) at redshift z=1.9 that uses popular values for cosmological and astrophysical parameters and accounts for all baryons with an uncertainty of 6%. We have measured the amount of absorption, DA, in the Ly-alpha forest at redshift 1.9 in spectra of 77 QSO from the Kast spectrograph. We calibrated the continuum fits with realistic artificial spectra, and we found that averaged over all 77 QSOs the mean continuum level is within 1-2% of the correct value. Absorption from all lines in the Ly-alpha forest at z=1.9 removes DA=15.1 +/- 0.7% of the flux between 1070 and 1170 (rest) Angstroms. This is the first measurement using many QSOs at this z, and the first calibrated measurement at any redshift. Metal lines absorb 2.3 +/- 0.5%, and LLS absorb 1.0 +/- 0.4% leaving 11.8 +/- 1.0% from the lower density bulk of the IGM. Averaging over Delta z=0.1 or 154 Mpc, the dispersion is 6.1 +/- 0.3% including LLS and metal lines, or 3.9 (+0.5, -0.7)% for the lower density IGM alone, consistent with the usual description of large scale structure. LLS and metal lines are major contributors to the variation in the mean flux, and they make the flux field significantly non-Gaussian. We find that a hydrodynamic simulation on a 1024 cubed grid in a 75.7 Mpc box reproduces the observed DA from the low density IGM with parameters values H_o=71 km/s/Mpc, Omega_Lambda=0.73, Omega_m=0.27, Omega_b=0.044, sigma_8=0.9 and a UV background that has an ionization rate that is 1.08 +/- 0.08 times the prediction by Madau, Haardt & Rees (1999).

We present a method to make predictions with sets of correlated data values, in this case QSO flux spectra. We predict the continuum in the Lyman-Alpha forest of a QSO, from 1020 -- 1216 A, using the spectrum of that QSO from 1216 -- 1600 A . We find correlations between the unabsorbed flux in these two wavelengths regions in the HST spectra of 50 QSOs. We use principal component analysis (PCA) to summarize the variety of these spectra and we relate the weights of the principal components for 1020 -- 1600 A to the weights for 1216 -- 1600 A, and we apply this relation to make predictions. We test the method on the HST spectra, and we find an average absolute flux error of 9%, with a range 3 -- 30%, where individual predictions are systematically too low or too high. We mention several ways in which the predictions might be improved.

We describe a new method to calibrate the relative flux levels in spectra from the HIRES echelle spectrograph on the Keck-I telescope. Standard data reduction techniques that transfer the instrument response between HIRES integrations leave errors in the flux of 5 - 10%, because the effective response varies. The flux errors are most severe near the ends of each spectral order, where there can be discontinuous jumps. The source of these errors is uncertain, but may include changes in the vignetting connected to the optical alignment. Our new flux calibration method uses a calibrated reference spectrum of each target to calibrate individual HIRES integrations. We determine the instrument response independently for each integration, and hence we avoid the need to transfer the instrument response between HIRES integrations. The procedure can be applied to any HIRES spectrum, or any other spectrum. While the accuracy of the method depends upon many factors, we have been able to flux calibrate a HIRES spectrum to 1% over scales of 200 A that include order joins. We illustrate the method with spectra of Q1243+3047 towards which we have measured the deuterium to hydrogen abundance ratio.

We report the detection of Deuterium absorption at redshift 2.525659 towards Q1243+3047. We describe improved methods to estimate the Deuterium to Hydrogen abundance ratio (D/H) in absorption systems, including improved modeling of the continuum level, the Ly-alpha forest and the velocity structure of the absorption. Together with improved relative flux calibration, these methods give D/H = 2.42^+0.35_-0.25 x 10^-5 cm^-2 from our Keck-I HIRES spectra of Q1243+3047, where the error is from the uncertainty in the shape of the continuum level and the amount of D absorption in a minor second component. The measured D/H is likely the primordial value because the [O/H] = -2.79 +/- 0.05. This absorption system has a neutral Hydrogen column density of 19.73 +/- 0.04 cm^-2, it shows five D lines and is mostly ionized. The best estimate of the primordial D/H is 2.78^+0.44_-0.38 x 10^-5, from the log D/H values towards five QSOs. The dispersion in the five values is larger than we expect from their individual measurement errors and we suspect this is because some of these errors were underestimated. We observe a trend in D/H with neutral H column density that we also suspect is spurious. The D/H corresponds to a baryon-to-photon ratio ETA = 5.9 +/- 0.5 x 10^-10 and a cosmological baryon density Omega_b h^2 = 0.0214 +/- 0.0020 (9.3%) that agrees with values from the pre-MAP measurements of the anisotropy of the Cosmic Microwave Background.

We present 14 N^0 measurements from our HIRES/Keck database of damped Lya abundances. These data are combined with measurements from the recent and past literature to build an homogeneous, uniform set of observations. We examine photoionization diagnostics like Fe^++ and Ar^0 in the majority of the complete sample and assess the impact of ionization corrections on N/alpha and alpha/H values derived from observed ionic column densities of N^0, Si^+, H^0, and S^+. Our final sample of 19 N/alpha, alpha/H pairs appears bimodal; the majority of systems show N/alpha values consistent with metal-poor emission regions in the local universe but a small sub-sample exhibit significantly lower N/alpha ratios. Contrary to previous studies of N/alpha in the damped systems, our sample shows little scatter within each sub-sample. We consider various scenarios to explain the presence of the low N/alpha sightlines and account for the apparent bimodality. We favor a model where at least some galaxies undergo an initial burst of star formation with suppressed formation of intermediate-mass stars. We found a power-law IMF with slope 0.10 or a mass cut of ~5-8 Msolar would successfully reproduce the observed LN-DLA values. If the bimodal distribution is confirmed by a larger sample of measurements, this may present the first observational evidence for a top heavy initial mass function in some early stellar populations.

We investigate the ionization state of the damped Lya system at z=2.62 toward GB1759+7539 through an analysis of ionic ratios sensitive to photoionization: ArI/SII, FeIII/FeII, NII/NI, AlIII/AlII. Approximately half of the metals arise in a mostly neutral velocity component with HI/H > 0.9, based on FeIII/FeII < 0.013. In contrast, the remaining half exhibits FeIII/FeII~0.3 indicative of a partially ionized medium with HI/H~0.5. These conclusions are supported by the observed NII/NI, AlIII/AlII, and ArI/SII ratios. We assess ionization corrections for the observed column densities through photoionization models derived from the CLOUDY software package. In the neutral gas, the ionization corrections are negligible except for ArI. However for the partially ionized gas, element abundance ratios differ from the ionic ratios by 0.1-0.3 dex for (SiII, SII, NiII, AlII)/FeII ratios and more for (NI, ArI)/FeII. Independent of the shape of the photoionizing spectrum and assumptions on the number of ionization phases, these ionization corrections have minimal impact (<0.1dex) on the total metallicity inferred for this damped Lya system. Measurements on the relative elemental abundances of the partially ionized gas, however, have a greater than ~0.15 dex uncertainty which hides the effects of nucleosynthesis and dust depletion. We caution the reader that this damped system is unusual for a number of reasons (e.g. a very low ArI/SII ratio) and we believe its ionization properties are special but not unique. Nevertheless, it clearly shows the value of examining photoionization diagnostics like FeIII/FeII in a larger sample of systems.

The z ~ 0.701 absorption system towards the QSO PG1718+4807 is the only example of a QSO absorption system which might have a deuterium/hydrogen ratio approximately ten times the value found towards other QSOs. We have obtained new STIS spectra from the Hubble Space Telescope of the Lyman alpha and Lyman limit regions of the system. These spectra give the redshift and velocity dispersion of the neutral hydrogen which produces most of the observed absorption. The Lyman alpha line is too narrow to account for all of the observed absorption. It was previously known that extra absorption is needed on the blue side of the main H I near to the expected position of deuterium. The current data suggests with a 98% confidence level that the extra absorption is not deuterium. Some uncertainty persists because we have a low signal to noise ratio and the extra absorption - be it deuterium or hydrogen - is heavily blended with the Lyman alpha absorption from the main hydrogen component.

We report the measurement of the primordial D/H abundance ratio towards QSO \object. The column density of the hydrogen in the $z \simeq 2.536$ Lyman limit system is high, \lnhi $= 19.422 \pm 0.009$ \cmm, allowing for the deuterium to be seen in 5 Lyman series transitions. The measured value of the D/H ratio towards QSO \object is found to be D/H$ = 2.54 \pm 0.23 \times 10^{-5}$. The metallicity of the system showing D/H is found to be $\simeq 0.01$ solar, indicating that the measured D/H is the primordial D/H within the measurement errors. The gas which shows D/H is neutral, unlike previous D/H systems which were more highly ionized. Thus, the determination of the D/H ratio becomes more secure since we are measuring it in different astrophysical environments, but the error is larger because we now see more dispersion between measurements. Combined with prior measurements of D/H, the best D/H ratio is now D/H$ = 3.0 \pm 0.4 \times 10^{-5}$, which is 10% lower than the previous value. The new values for the baryon to photon ratio, and baryonic matter density derived from D/H are $\eta = 5.6 \pm 0.5 \times 10^{-10} $ and \ob $=0.0205 \pm 0.0018$ respectively.

Big Bang Nucleosynthesis (BBN) is the synthesis of the light nuclei, Deuterium, He3, He4 and Li7, during the first few minutes of the universe. This review concentrates on recent improvements in the measurement of the primordial (after BBN, and prior to modification) abundances of these nuclei. We mention improvement in the standard theory, and the non-standard extensions which are limited by the data. (abridged)