The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

The true nature of type Ia supernovae (SNe Ia) remains largely unknown, and as survey statistics increase, the question of astrophysical systematic uncertainties rises, notably that of the SN Ia population evolution. In this paper, we study the dependence with redshift of the SN Ia SALT2.4 lightcurve stretch, a purely intrinsic SN property, to probe its potential redshift drift. The SN stretch has been shown to strongly correlate with the SN environment, notably with stellar age tracers. We model the underlying stretch distribution as a function of redshift, using the evolution of fraction of young and old SNe Ia as predicted by Rigault et al. (2018), and assuming constant underlying stretch distribution for each age population made of Gaussian mixtures. We test our prediction against published samples chosen to have negligible magnitude selection effects, so that any observed change is indeed of astrophysical and not observational origin. We clearly demonstrate that the underlying SN Ia stretch distribution is evolving as a function of redshift, and that the young/old drifting model is a much better description of the data than any time-constant model, including the sample-based asymmetric distributions usually used to correct Malmquist bias. The favored underlying stretch model is the bimodal one derived from Rigault et al. (2018): a high-stretch mode shared by both young and old environments, and a low-stretch mode exclusive to old environments. The precise impact of the redshift evolution of the SN Ia population intrinsic properties on cosmology remains to be studied. Yet,the astrophysical drift of the SN stretch distribution does affect current Malmquist bias corrections and thereby distances derived from SN affected by selection effects. We highlight that such a bias will increase with surveys covering increasingly larger redshift ranges, which is particularly important for LSST.

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$.

Type Ia Supernovae have yet again the opportunity to revolutionize the field of cosmology as the new generation of surveys are acquiring thousands of nearby SNeIa opening a new era in cosmology: the direct measurement of the growth of structure parametrized by $fD$. This method is based on the SNeIa peculiar velocities derived from the residual to the Hubble law as direct tracers of the full gravitational potential caused by large scale structure. With this technique, we could probe not only the properties of dark energy, but also the laws of gravity. In this paper we present the analytical framework and forecasts. We show that ZTF and LSST will be able to reach 5\% precision on $fD$ by 2027. Our analysis is not significantly sensitive to photo-typing, but known selection functions and spectroscopic redshifts are mandatory. We finally introduce an idea of a dedicated spectrograph that would get all the required information in addition to boost the efficiency to each SNeIa so that we could reach the 5\% precision within the first two years of LSST operation and the few percent level by the end of the survey.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

Slitless spectroscopy has long been considered as a complicated and confused technique. Nonetheless, with the advent of Hubble Space Telescope (HST) instruments characterized by a low sky background level and a high spatial resolution (most notably WFC3), slitless spectroscopy has become an adopted survey tool to study galaxy evolution from space. We investigate its application to single object studies to measure not only redshift and integrated spectral features, but also spatially resolved quantities such as galaxy kinematics. We build a complete forward model to be quantitatively compared to actual slitless observations. This model depends on a simplified thin cold disk galaxy description -- including flux distribution, intrinsic spectrum and kinematic parameters -- and on the instrumental signature. It is used to improve redshifts and constrain basic rotation curve parameters, i.e. plateau velocity $v_{0}$ (in $\rm km.s^{-1}$) and central velocity gradient $w_{0}$ (in $\rm km.s^{-1}.arcsec^{-1}$). The model is tested on selected observations from 3D-HST and GLASS surveys, to estimate redshift and kinematic parameters on several galaxies measured with one or more roll angles. Our forward approach allows to mitigate the self-contamination effect, a primary drawback of slitless spectroscopy, and therefore has the potential to increase precision on redshifts. In a limited sample of well-resolved spiral galaxies from HST surveys, it is possible to significantly constrain galaxy rotation curve parameters. This proof-of-concept work is promising for future large slitless spectroscopic surveys such as EUCLID and WFIRST.

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.

As a candidate 'super-Chandrasekhar' or 09dc-like Type Ia supernova (SN Ia), SN 2012dn shares many characteristics with other members of this remarkable class of objects but lacks their extraordinary luminosity. Here, we present and discuss the most comprehensive optical data set of this SN to date, comprised of a densely sampled series of early-time spectra obtained within the Nearby Supernova Factory project, plus photometry and spectroscopy obtained at the VLT about 1 yr after the explosion. The light curves, colour curves, spectral time series and ejecta velocities of SN 2012dn are compared with those of other 09dc-like and normal SNe Ia, the overall variety within the class of 09dc-like SNe Ia is discussed, and new criteria for 09dc-likeness are proposed. Particular attention is directed to additional insight that the late-phase data provide. The nebular spectra show forbidden lines of oxygen and calcium, elements that are usually not seen in late-time spectra of SNe Ia, while the ionisation state of the emitting iron plasma is low, pointing to low ejecta temperatures and high densities. The optical light curves are characterised by an enhanced fading starting ~60 d after maximum and very low luminosities in the nebular phase, which is most readily explained by unusually early formation of clumpy dust in the ejecta. Taken together, these effects suggest a strongly perturbed ejecta density profile, which might lend support to the idea that 09dc-like characteristics arise from a brief episode of interaction with a hydrogen-deficient envelope during the first hours or days after the explosion.

This letter presents a new approach using the cosmic peculiar velocity field to characterize the morphology and size of large scale structures in the local Universe. The algorithm developed uses the three-dimensional peculiar velocity field to compute flow lines, or streamlines. The local Universe is then partitioned into volumes corresponding to gravitational basins, also called watersheds, among the different end-points of the velocity flow lines. This new methodology is first tested on numerical cosmological simulations, used as benchmark for the method, and then applied to the {\it Cosmic-Flows} project observational data in order to to pay particular attention to the nearby superclusters including ours. More extensive tests on both simulated and observational data will be discussed in an accompanying paper.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

Current time domain facilities are discovering hundreds of new galactic and extra-galactic transients every week. Classifying the ever-increasing number of transients is challenging, yet crucial to further our understanding of their nature, discover new classes, or ensuring sample purity, for instance, for Supernova Ia cosmology. The Zwicky Transient Facility is one example of such a survey. In addition, it has a dedicated very-low resolution spectrograph, the SEDMachine, operating on the Palomar 60-inch telescope. This spectrograph's primary aim is object classification. In practice most, if not all, transients of interest brighter than ~19 mag are typed. This corresponds to approximately 10 to 15 targets a night. In this paper, we present a fully automated pipeline for the SEDMachine. This pipeline has been designed to be fast, robust, stable and extremely flexible. pysedm enables the fully automated spectral extraction of a targeted point source object in less than 5 minutes after the end of the exposure. The spectral color calibration is accurate at the few percent level. In the 19 weeks since pysedm entered production in early August of 2018, we have classified, among other objects, about 400 Type Ia supernovae and 140 Type II supernovae. We conclude that low resolution, fully automated spectrographs such as the `SEDMachine with pysedm' installed on 2-m class telescopes within the southern hemisphere could allow us to automatically and simultaneously type and obtain a redshift for most (if not all) bright transients detected by LSST within z<0.2, notably potentially all Type Ia Supernovae. In comparison to the current SEDM design, this would require higher spectral resolution (R~1000) and slightly improved throughput. With this perspective in mind, pysedm has been designed to easily be adaptable to any IFU-like spectrograph (see https://github.com/MickaelRigault/pysedm).

A hierarchical Bayesian model is applied to the Cosmicflows-3 catalog of galaxy distances in order to derive the peculiar velocity field and distribution of matter within $z \sim 0.054$. The model assumes the $\Lambda$CDM model within the linear regime and includes the fit of the galaxy distances together with the underlying density field. By forward modeling the data, the method is able to mitigate biases inherent to peculiar velocity analyses, such as the Homogeneous Malmquist bias or the log-normal distribution of peculiar velocities. The statistical uncertainty on the recovered velocity field is about 150 km/s depending on the location, and we study systematics coming from the selection function and calibration of distance indicators. The resulting velocity field and related density fields recover the cosmography of the Local Universe which is presented in an unprecedented volume of universe 10 times larger than previously reached. This methodology open the doors to reconstruction of initial conditions for larger and more accurate constrained cosmological simulations. This work is also preparatory to larger peculiar velocity datasets coming from Wallaby, TAIPAN or LSST.

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.

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]

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.

We have discovered an anomalous behavior of CCD readout electronics that affects their use in many astronomical applications. An offset in the digitization of the CCD output voltage that depends on the binary encoding of one pixel is added to pixels that are read out one, two and/or three pixels later. One result of this effect is the introduction of a differential offset in the background when comparing regions with and without flux from science targets. Conventional data reduction methods do not correct for this offset. We find this effect in 16 of 22 instruments investigated, covering a variety of telescopes and many different front-end electronics systems. The affected instruments include LRIS and DEIMOS on the Keck telescopes, WFC3-UVIS and STIS on HST, MegaCam on CFHT, SNIFS on the UH88 telescope, GMOS on the Gemini telescopes, HSC on Subaru, and FORS on VLT. The amplitude of the introduced offset is up to 4.5 ADU per pixel, and it is not directly proportional to the measured ADU level. We have developed a model that can be used to detect this "binary offset effect" in data and correct for it. Understanding how data are affected and applying a correction for the effect is essential for precise astronomical measurements.

We present a new approach to automate the spectroscopic redshift reliability assessment based on machine learning (ML) and characteristics of the redshift probability density function (PDF). We propose to rephrase the spectroscopic redshift estimation into a Bayesian framework, in order to incorporate all sources of information and uncertainties related to the redshift estimation process, and produce a redshift posterior PDF that will be the starting-point for ML algorithms to provide an automated assessment of a redshift reliability. As a use case, public data from the VIMOS VLT Deep Survey is exploited to present and test this new methodology. We first tried to reproduce the existing reliability flags using supervised classification to describe different types of redshift PDFs, but due to the subjective definition of these flags, soon opted for a new homogeneous partitioning of the data into distinct clusters via unsupervised classification. After assessing the accuracy of the new clusters via resubstitution and test predictions, unlabelled data from preliminary mock simulations for the Euclid space mission are projected into this mapping to predict their redshift reliability labels.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

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.

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.

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 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

Context. Our Local Group of galaxies appears to be moving relative to the cosmic microwave background with the source of the peculiar motion still uncertain. While in the past this has been studied mostly using galaxies as distance indicators, the weight of type Ia supernovae (SNe Ia) has increased recently with the continuously improving statistics of available low-redshift supernovae. Aims. We measured the bulk flow in the nearby universe ($0.015 < z < 0.1$) using 117 SNe Ia observed by the Nearby Supernova Factory, as well as the Union2 compilation of SN Ia data already in the literature. Methods. The bulk flow velocity was determined from SN data binned in redshift shells by including a coherent motion (dipole) in a cosmological fit. Additionally, a method of spatially smoothing the Hubble residuals was used to verify the results of the dipole fit. To constrain the location and mass of a potential mass concentration (e.g., the Shapley supercluster) responsible for the peculiar motion, we fit a Hubble law modified by adding an additional mass concentration. Results. The analysis shows a bulk flow that is consistent with the direction of the CMB dipole up to $z \sim 0.06$, thereby doubling the volume over which conventional distance measures are sensitive to a bulk flow. We see no significant turnover behind the center of the Shapley supercluster. A simple attractor model in the proximity of the Shapley supercluster is only marginally consistent with our data, suggesting the need for another, more distant source. In the redshift shell $0.06 < z < 0.1$, we constrain the bulk flow velocity to $< 240~\textrm{km s}^{-1}$ (68% confidence level) for the direction of the CMB dipole, in contradiction to recent claims of the existence of a large-amplitude dark flow.

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.

We estimate systematic errors due to K-corrections in standard photometric analyses of high redshift Type Ia supernovae. Errors due to K-correction occur when the spectral template model underlying the lightcurve fitter poorly represents the actual supernova spectral energy distribution, meaning that the distance modulus cannot be recovered accurately. In order to quantify this effect, synthetic photometry is performed on artificially redshifted spectrophotometric data from 119 low-redshift supernovae from the Nearby Supernova Factory, and the resulting lightcurves are fit with a conventional lightcurve fitter. We measure the variation in the standardized magnitude that would be fit for a given supernova if located at a range of redshifts and observed with various filter sets corresponding to current and future supernova surveys. We find significant variation in the measurements of the same supernovae placed at different redshifts regardless of filters used, which causes dispersion greater than $\sim0.05$ mag for measurements of photometry using the Sloan-like filters and a bias that corresponds to a $0.03$ shift in $w$ when applied to an outside data set. To test the result of a shift in supernova population or environment at higher redshifts, we repeat our calculations with the addition of a reweighting of the supernovae as a function of redshift and find that this strongly affects the results and would have repercussions for cosmology. We discuss possible methods to reduce the contribution of the K-correction bias and uncertainty.

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 present a sample of normal type Ia supernovae from the Nearby Supernova Factory dataset with spectrophotometry at sufficiently late phases to estimate the ejected mass using the bolometric light curve. We measure $^{56}$Ni masses from the peak bolometric luminosity, then compare the luminosity in the $^{56}$Co-decay tail to the expected rate of radioactive energy re- lease from ejecta of a given mass. We infer the ejected mass in a Bayesian context using a semi-analytic model of the ejecta, incorporating constraints from contemporary numerical models as priors on the density structure and distribution of $^{56}$Ni throughout the ejecta. We find a strong correlation between ejected mass and light curve decline rate, and consequently $^{56}$Ni mass, with ejected masses in our data ranging from 0.9-1.4 $M_\odot$. Most fast-declining (SALT2 $x_1 < -1$) normal SNe Ia have significantly sub-Chandrasekhar ejected masses in our fiducial analysis.

Kim et al. (2013) [K13] introduced a new methodology for determining peak-brightness absolute magnitudes of type Ia supernovae from multi-band light curves. We examine the relation between their parameterization of light curves and Hubble residuals, based on photometry synthesized from the Nearby Supernova Factory spectrophotometric time series, with global host-galaxy properties. The K13 Hubble residual step with host mass is $0.013\pm 0.031$ mag for a supernova subsample with data coverage corresponding to the K13 training; at $\ll 1\sigma$, the step is not significant and lower than previous measurements. Relaxing the data coverage requirement the Hubble residual step with host mass is $0.045\pm 0.026$ mag for the larger sample; a calculation using the modes of the distributions, less sensitive to outliers, yields a step of 0.019 mag. The analysis of this article uses K13 inferred luminosities, as distinguished from previous works that use magnitude corrections as a function of SALT2 color and stretch parameters: Steps at $>2\sigma $ significance are found in SALT2 Hubble residuals in samples split by the values of their K13 $x(1)$ and $x(2)$ light-curve parameters. $x(1)$ affects the light-curve width and color around peak (similar to the $\Delta m_{15}$ and stretch parameters), and $x(2)$ affects colors, the near-UV light-curve width, and the light-curve decline 20 to 30 days after peak brightness. The novel light-curve analysis, increased parameter set, and magnitude corrections of K13 may be capturing features of SN~Ia diversity arising from progenitor stellar evolution.

(Abridged) We study the host galaxy regions in close proximity to Type Ia supernovae (SNe Ia) to analyze relations between the properties of SN Ia events and environments most similar to where their progenitors formed. We focus on local H\alpha\ emission as an indicator of young environments. The Nearby Supernova Factory has obtained flux-calibrated spectral timeseries for SNe Ia using integral field spectroscopy, allowing the simultaneous measurement of the SN and its immediate vicinity. For 89 SNe Ia we measure H\alpha\ emission tracing ongoing star formation within a 1 kpc radius around each SN. This constitutes the first direct study of the local environment for a large sample of SNe Ia also having accurate luminosity, color and stretch measurements. We find that SNe Ia with local H\alpha\ emission are redder by 0.036+/-0.017 mag, and that the previously-noted correlation between stretch and host mass is entirely driven by the SNe Ia coming from passive regions. Most importantly, the mean standardized brightness for SNe Ia with local H\alpha\ emission is 0.094+/-0.031 mag fainter than for those without. This offset arises from a bimodal structure in the Hubble residuals, that also explains the previously-known host-mass bias. We combine this bimodality with the cosmic star-formation rate to predict changes with redshift in the mean SN Ia brightness and the host-mass bias. This change is confirmed using high-redshift SNe Ia from the literature. These environmental dependences point to remaining systematic errors in SNe Ia standardization. The observed brightness offset is predicted to cause a significant bias in measurements of the dark energy equation of state. Recognition of these effects offers new opportunities to improve SNe Ia as cosmological probes - e.g. SNe Ia having local H\alpha\ emission are more homogeneous, having a brightness dispersion of 0.105+/-0.012 mag.

We present photometric and spectroscopic observations of galaxies hosting Type Ia supernovae (SNe Ia) observed by the Nearby Supernova Factory (SNfactory). Combining GALEX UV data with optical and near infrared photometry, we employ stellar population synthesis techniques to measure SN Ia host galaxy stellar masses, star-formation rates (SFRs), and reddening due to dust. We reinforce the key role of GALEX UV data in deriving accurate estimates of galaxy SFRs and dust extinction. Optical spectra of SN Ia host galaxies are fitted simultaneously for their stellar continua and emission lines fluxes, from which we derive high precision redshifts, gas-phase metallicities, and Halpha-based SFRs. With these data we show that SN Ia host galaxies present tight agreement with the fiducial galaxy mass-metallicity relation from SDSS for stellar masses log(M_*/M_Sun)>8.5 where the relation is well-defined. The star-formation activity of SN Ia host galaxies is consistent with a sample of comparable SDSS field galaxies, though this comparison is limited by systematic uncertainties in SFR measurements. Our analysis indicates that SN Ia host galaxies are, on average, typical representatives of normal field galaxies.

We examine the relationship between Type Ia Supernova (SN Ia) Hubble residuals and the properties of their host galaxies using a sample of 115 SNe Ia from the Nearby Supernova Factory (SNfactory). We use host galaxy stellar masses and specific star-formation rates fitted from photometry for all hosts, as well as gas-phase metallicities for a subset of 69 star-forming (non-AGN) hosts, to show that the SN Ia Hubble residuals correlate with each of these host properties. With these data we find new evidence for a correlation between SN Ia intrinsic color and host metallicity. When we combine our data with those of other published SN Ia surveys, we find the difference between mean SN Ia brightnesses in low and high mass hosts is 0.077 +- 0.014 mag. When viewed in narrow (0.2 dex) bins of host stellar mass, the data reveal apparent plateaus of Hubble residuals at high and low host masses with a rapid transition over a short mass range (9.8 <= log(M_*/M_Sun) <= 10.4). Although metallicity has been a favored interpretation for the origin of the Hubble residual trend with host mass, we illustrate how dust in star-forming galaxies and mean SN Ia progenitor age both evolve along the galaxy mass sequence, thereby presenting equally viable explanations for some or all of the observed SN Ia host bias.

We present a novel class of models for Type Ia supernova time-evolving spectral energy distributions (SED) and absolute magnitudes: they are each modeled as stochastic functions described by Gaussian processes. The values of the SED and absolute magnitudes are defined through well-defined regression prescriptions, so that data directly inform the models. As a proof of concept, we implement a model for synthetic photometry built from the spectrophotometric time series from the Nearby Supernova Factory. Absolute magnitudes at peak $B$ brightness are calibrated to 0.13 mag in the $g$-band and to as low as 0.09 mag in the $z=0.25$ blueshifted $i$-band, where the dispersion includes contributions from measurement uncertainties and peculiar velocities. The methodology can be applied to spectrophotometric time series of supernovae that span a range of redshifts to simultaneously standardize supernovae together with fitting cosmological parameters.

We present 32 epochs of optical (3300-9700 \AA) spectrophotometric observations of the nearby quintessential "normal" type Ia supernova (SN Ia) SN 2011fe in the galaxy M101, extending from -15 to +97 d with respect to B-band maximum, obtained by the Nearby Supernova Factory collaboration. SN 2011fe is the closest (\mu = 29.04) and brightest (Bmax = 9.94 mag) SN Ia observed since the advent of modern large scale programs for the intensive periodic followup of supernovae. Both synthetic light curve measurements and spectral feature analysis attest to the normality of SN 2011fe. There is very little evidence for reddening in its host galaxy. The homogeneous calibration, intensive time sampling, and high signal-to-noise ratio of the data set make it unique. Thus it is ideal for studying the physics of SN Ia explosions in detail, and for furthering the use of SNe Ia as standardizable candles for cosmology. Several such applications are shown, from the creation of a bolometric light curve and measurement of the 56Ni mass, to the simulation of detection thresholds for unburned carbon, direct comparisons with other SNe Ia, and existing spectral templates.

We present a new atmospheric extinction curve for Mauna Kea spanning 3200--9700 \AA. It is the most comprehensive to date, being based on some 4285 standard star spectra obtained on 478 nights spread over a period of 7 years obtained by the Nearby SuperNova Factory using the SuperNova Integral Field Spectrograph. This mean curve and its dispersion can be used as an aid in calibrating spectroscopic or imaging data from Mauna Kea, and in estimating the calibration uncertainty associated with the use of a mean extinction curve. Our method for decomposing the extinction curve into physical components, and the ability to determine the chromatic portion of the extinction even on cloudy nights, is described and verified over the wide range of conditions sampled by our large dataset. We demonstrate good agreement with atmospheric science data obtain at nearby Mauna Loa Observatory, and with previously published measurements of the extinction above Mauna Kea.

We present optical photometry and spectroscopy of five type Ia supernovae discovered by the Nearby Supernova Factory selected to be spectroscopic analogues of the candidate super-Chandrasekhar-mass events SN 2003fg and SN 2007if. Their spectra are characterized by hot, highly ionized photospheres near maximum light, for which SN 1991T supplies the best phase coverage among available close spectral templates. Like SN 2007if, these supernovae are overluminous (-19.5 < M_V < -20) and the velocity of the Si II 6355 absorption minimum is consistent with being constant in time from phases as early as a week before, and up to two weeks after, $B$-band maximum light. We interpret the velocity plateaus as evidence for a reverse-shock shell in the ejecta formed by interaction at early times with a compact envelope of surrounding material, as might be expected for SNe resulting from the mergers of two white dwarfs. We use the bolometric light curves and line velocity evolution of these SNe to estimate important parameters of the progenitor systems, including nickel-56 mass, total progenitor mass, and masses of shells and surrounding carbon/oxygen envelopes. We find that the reconstructed total progenitor mass distribution of the events (including SN 2007if) is bounded from below by the Chandrasekhar mass, with SN 2007if being the most massive. We discuss the relationship of these events to the emerging class of super-Chandrasekhar-mass SNe Ia, estimate the relative rates, compare the mass distribution to that expected for double-degenerate SN Ia progenitors from population synthesis, and consider implications for future cosmological Hubble diagrams.

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.

We present convincing evidence of unburned carbon at photospheric velocities in new observations of 5 Type Ia supernovae (SNe Ia) obtained by the Nearby Supernova Factory. These SNe are identified by examining 346 spectra from 124 SNe obtained before +2.5 d relative to maximum. Detections are based on the presence of relatively strong C II 6580 absorption "notches" in multiple spectra of each SN, aided by automated fitting with the SYNAPPS code. Four of the 5 SNe in question are otherwise spectroscopically unremarkable, with ions and ejection velocities typical of SNe Ia, but spectra of the fifth exhibits high-velocity (v > 20,000 km/s) Si II and Ca II features. On the other hand, the light curve properties are preferentially grouped, strongly suggesting a connection between carbon-positivity and broad band light curve/color behavior: Three of the 5 have relatively narrow light curves but also blue colors, and a fourth may be a dust-reddened member of this family. Accounting for signal-to-noise and phase, we estimate that 22 +10/-6% of SNe Ia exhibit spectroscopic C II signatures as late as -5 d with respect to maximum. We place these new objects in the context of previously recognized carbon-positive SNe Ia, and consider reasonable scenarios seeking to explain a physical connection between light curve properties and the presence of photospheric carbon. We also examine the detailed evolution of the detected carbon signatures and the surrounding wavelength regions to shed light on the distribution of carbon in the ejecta. Our ability to reconstruct the C II 6580 feature in detail under the assumption of purely spherical symmetry casts doubt on a "carbon blobs" hypothesis, but does not rule out all asymmetric models. A low volume filling factor for carbon, combined with line-of-sight effects, seems unlikely to explain the scarcity of detected carbon in SNe Ia by itself.

We employ 76 type Ia supernovae with optical spectrophotometry within 2.5 days of B-band maximum light obtained by the Nearby Supernova Factory to derive the impact of Si and Ca features on supernovae intrinsic luminosity and determine a dust reddening law. We use the equivalent width of Si II {\lambda}4131 in place of light curve stretch to account for first-order intrinsic luminosity variability. The resultant empirical spectral reddening law exhibits strong features associated with Ca II and Si II {\lambda}6355. After applying a correction based on the Ca II H&K equivalent width we find a reddening law consistent with a Cardelli extinction law. Using the same input data, we compare this result to synthetic rest-frame UBVRI-like photometry in order to mimic literature observations. After corrections for signatures correlated with Si II {\lambda}4131 and Ca II H&K equivalent widths, and introducing an empirical correlation between colors, we determine the dust component in each band. We find a value of the total-to-selective extinction ratio, RV = 2.8 \pm 0.3. This agrees with the Milky Way value, in contrast to the low RV values found in most previous analyses. This result suggests that the long-standing controversy in interpreting SN Ia colors and their compatibility with a classical extinction law, critical to their use as cosmological probes, can be explained by the treatment of the dispersion in colors, and by the variability of features apparent in SN Ia spectra.

We present Keck LRIS spectroscopy and $g$-band photometry of the metal-poor, low-luminosity host galaxy of the super-Chandrasekhar mass Type Ia supernova SN 2007if. Deep imaging of the host reveals its apparent magnitude to be $m_g=23.15\pm0.06$, which at the spectroscopically-measured redshift of $z_{helio}=0.07450\pm0.00015$ corresponds to an absolute magnitude of $M_g=-14.45\pm0.06$. Galaxy $g-r$ color constrains the mass-to-light ratio, giving a host stellar mass estimate of $\log(M_*/M_\odot)=7.32\pm0.17$. Balmer absorption in the stellar continuum, along with the strength of the 4000\AA\ break, constrain the age of the dominant starburst in the galaxy to be $t_\mathrm{burst}=123^{+165}_{-77}$ Myr, corresponding to a main-sequence turn-off mass of $M/M_\odot=4.6^{+2.6}_{-1.4}$. Using the R$_{23}$ method of calculating metallicity from the fluxes of strong emission lines, we determine the host oxygen abundance to be $12+\log(O/H)_\mathrm{KK04}=8.01\pm0.09$, significantly lower than any previously reported spectroscopically-measured Type Ia supernova host galaxy metallicity. Our data show that SN 2007if is very likely to have originated from a young, metal-poor progenitor.

We present photometric and spectroscopic observations of SN 2007if, an overluminous (M_V = -20.4), red (B-V = 0.16 at B-band maximum), slow-rising (t_rise = 24 days) type Ia supernova in a very faint (M_g = -14.10) host galaxy. A spectrum at 5 days past B-band maximum light is a direct match to the super-Chandrasekhar-mass candidate SN Ia 2003fg, showing Si II and C II at ~9000 km/s. A high signal-to-noise co-addition of the SN spectral time series reveals no Na I D absorption, suggesting negligible reddening in the host galaxy, and the late-time color evolution has the same slope as the Lira relation for normal SNe Ia. The ejecta appear to be well mixed, with no strong maximum in I-band and a diversity of iron-peak lines appearing in near-maximum-light spectra. SN2007 if also displays a plateau in the Si II velocity extending as late as +10 days, which we interpret as evidence for an overdense shell in the SN ejecta. We calculate the bolometric light curve of the SN and use it and the \ion{Si}{2} velocity evolution to constrain the mass of the shell and the underlying SN ejecta, and demonstrate that SN2007 if is strongly inconsistent with a Chandrasekhar-mass scenario. Within the context of a "tamped detonation" model appropriate for double-degenerate mergers, and assuming no host extinction, we estimate the total mass of the system to be 2.4 +/- 0.2 solar masses, with 1.6 +/- 0.1 solar masses of nickel-56 and with 0.3-0.5 solar masses in the form of an envelope of unburned carbon/oxygen. Our modeling demonstrates that the kinematics of shell entrainment provide a more efficient mechanism than incomplete nuclear burning for producing the low velocities typical of super-Chandrasekhar-mass SNeIa.

We present a new method to standardize Type Ia supernova (SN Ia) luminosities to ~<0.13 magnitudes using flux ratios from a single flux-calibrated spectrum per SN. Using Nearby Supernova Factory spectrophotomery of 58 SNe Ia, we performed an unbiased search for flux ratios which correlate with SN Ia luminosity. After developing the method and selecting the best ratios from a training sample, we verified the results on a separate validation sample and with data from the literature. We identified multiple flux ratios whose correlations with luminosity are stronger than those of light curve shape and color, previously identified spectral feature ratios, or equivalent width measurements. In particular, the flux ratio R(642/443) = F(642 nm) / F(443 nm) has a correlation of 0.95 with SN Ia absolute magnitudes. Using this single ratio as a correction factor produces a Hubble diagram with a residual scatter standard deviation of 0.125 +- 0.011 mag, compared with 0.161 +- 0.015 mag when fit with the SALT2 light curve shape and color parameters x1 and c. The ratio R(642/443) is an effective correction factor for both extrinsic dust reddening and instrinsic variations such as those of SN 1991T-like and SN 1999aa-like SNe. When combined with broad-band color measurements, spectral flux ratios can standardize SN Ia magnitudes to ~0.12 mag. These are the first spectral metrics that improve over the standard normalization methods based upon light curve shape and color and they provide among the lowest scatter Hubble diagrams ever published.

The use of Type Ia supernovae as distance indicators led to the discovery of the accelerating expansion of the universe a decade ago. Now that large second generation surveys have significantly increased the size and quality of the high-redshift sample, the cosmological constraints are limited by the currently available sample of ~50 cosmologically useful nearby supernovae. The Nearby Supernova Factory addresses this problem by discovering nearby supernovae and observing their spectrophotometric time development. Our data sample includes over 2400 spectra from spectral timeseries of 185 supernovae. This talk presents results from a portion of this sample including a Hubble diagram (relative distance vs. redshift) and a description of some analyses using this rich dataset.

We present four spectra of the Type Ia supernova (SN Ia) 2006D extending from -7 to +13 days with respect to B-band maximum. The spectra include the strongest signature of unburned material at photospheric velocities observed in a SN Ia to date. The earliest spectrum exhibits C II absorption features below 14,000 km/s, including a distinctive C II \lambda 6580 absorption feature. The carbon signatures dissipate as the SN approaches peak brightness. In addition to discussing implications of photospheric-velocity carbon for white dwarf explosion models, we outline some factors that may influence the frequency of its detection before and around peak brightness. Two effects are explored in this regard, including depopulation of the C II optical levels by non-LTE effects, and line-of-sight effects resulting from a clumpy distribution of unburned material with low volume-filling factor.

We have obtained optical spectrophotometry of the evolution of comet 9P/Tempel 1 after the impact of the Deep Impact probe, using the Supernova Integral Field Spectrograph (SNIFS) at the UH 2.2m telescope, as well as simultaneous optical and infrared spectra using the Lick Visible-to-Near-Infrared Imaging Spectrograph (VNIRIS) spectrograph. The spatial distribution and temporal evolution of the "violet band" CN (0-0) emission and of the 630 nm [OI] emission was studied. We found that CN emission centered on the nucleus increased in the two hours after impact, but that this CN emission was delayed compared to the light curve of dust-scattered sunlight. The CN emission also expanded faster than the cloud of scattering dust. The emission of [OI] at 630 nm rose similarly to the scattered light, but then remained nearly constant for several hours after impact. On the day following the impact, both CN and [OI] emission concentrated on the comet nucleus had returned nearly to pre-impact levels. We have also searched for differences in the scattering properties of the dust ejected by the impact compared to the dust released under normal conditions. Compared to the pre-impact state of the comet, we find evidence that the color of the comet was slightly bluer during the post-impact rise in brightness. Long after the impact, in the following nights, the comet colors returned to their pre-impact values. This can be explained by postulating a change to a smaller particle size distribution in the ejecta cloud, in agreement with the findings from mid-infrared observatons, or by postulating a large fraction of clean ice particles, or by a combination of these two.

We present a study of the type IIn supernova (SN) 2005gl, in the relatively nearby (d~66 Mpc) galaxy NGC 266. Photometry and spectroscopy of the SN indicate it is a typical member of its class. Pre-explosion Hubble Space Telescope (HST) imaging of the location of the SN, along with a precise localization of this event using the Laser-Guide-Star assisted Adaptive Optics (LGS-AO) system at Keck Observatory, are combined to identify a luminous (M_V=-10.3) point source as the possible progenitor of SN 2005gl. If the source is indeed a single star, it was likely a member of the class of luminous blue variable stars (LBVs). This finding leads us to consider the possible general association of SNe IIn with LBV progenitors. We find this is indeed supported by observations of other SNe, and the known properties of LBV stars. For example, we argue that should the prototypical Galactic LBV eta Carina explode in a phase similar to its current state, it will likely produce a type IIn SN. We discuss our findings in the context of current ideas about the evolution of massive stars, and review the census of SNe with identified progenitors. We introduce the concept of the progenitor-SN map as a convenient means to discuss the present status and future prospects of direct searches for SN progenitors. We conclude that this field has matured considerably in recent years, and the transition from anecdotal information about rare single events to robust associations of progenitor classes with specific SN types has already begun.

We report Nearby Supernova Factory observations of SN 2005gj, the second confirmed case of a "hybrid" Type Ia/IIn supernova. Our early-phase photometry of SN 2005gj shows that the interaction is much stronger than for the prototype, SN 2002ic. Our first spectrum shows a hot continuum with broad and narrow H-alpha emission. Later spectra, spanning over 4 months from outburst, show clear Type Ia features combined with broad and narrow H-gamma, H-beta, H-alpha and HeI 5876,7065 in emission. At higher resolution, P Cygni profiles are apparent. Surprisingly, we also observe an inverted P Cygni profile for [OIII] 5007. We find that the lightcurve and measured velocity of the unshocked circumstellar material imply mass loss as recently as 8 years ago. The early lightcurve is well-described by a flat radial density profile for the circumstellar material. However, our decomposition of the spectra into Type Ia and shock emission components allows for little obscuration of the supernova, suggesting an aspherical or clumpy distribution for the circumstellar material. We suggest that the emission line velocity profiles arise from electron scattering rather than the kinematics of the shock. This is supported by the inferred high densities, and the lack of evidence for evolution in the line widths. Ground- and space-based photometry, and Keck spectroscopy, of the host galaxy are used to ascertain that the host galaxy has low metallicity Z/Zsun < 0.3; (95% confidence) and that this galaxy is undergoing a significant star formation event that began roughly 200+/-70 Myr ago. We discuss the implications of these observations for progenitor models and cosmology using Type Ia supernovae.

We present a generalisation of surface photometry to the higher-order moments of the line-of-sight velocity distribution of galaxies observed with integral-field spectrographs. The generalisation follows the approach of surface photometry by determining the best fitting ellipses along which the profiles of the moments can be extracted and analysed by means of harmonic expansion. The assumption for the odd moments (e.g. mean velocity) is that the profile along an ellipse satisfies a simple cosine law. The assumption for the even moments (e.g velocity dispersion) is that the profile is constant, as it is used in surface photometry. We find that velocity profiles extracted along ellipses of early-type galaxies are well represented by the simple cosine law (with 2% accuracy), while possible deviations are carried in the fifth harmonic term which is sensitive to the existence of multiple kinematic components, and has some analogy to the shape parameter of photometry. We compare the properties of the kinematic and photometric ellipses and find that they are often very similar. Finally, we offer a characterisation of the main velocity structures based only on the kinemetric parameters which can be used to quantify the features in velocity maps (abridged).

We present the stellar kinematics of 48 representative elliptical and lenticular galaxies obtained with our custom-built integral-field spectrograph SAURON operating on the William Herschel Telescope. The data were homogeneously processed through a dedicated reduction and analysis pipeline. All resulting SAURON datacubes were spatially binned to a constant minimum signal-to-noise. We have measured the stellar kinematics with an optimized (penalized pixel-fitting) routine which fits the spectra in pixel space, via the use of optimal templates, and prevents the presence of emission lines to affect the measurements. We have thus generated maps of the mean stellar velocity, the velocity dispersion, and the Gauss-Hermite moments h3 and h4 of the line-of-sight velocity distributions. The maps extend to approximately one effective radius. Many objects display kinematic twists, kinematically decoupled components, central stellar disks, and other peculiarities, the nature of which will be discussed in future papers of this series.

We present SAURON integral-field observations of the S0 galaxy NGC7332. Existing broad-band ground-based and HST photometry reveals a double disk structure and a boxy bulge interpreted as a bar viewed close to edge-on. The SAURON two-dimensional stellar kinematic maps confirm the existence of the bar and inner disk but also uncover the presence of a cold counter-rotating stellar component within the central 250 pc. The Hbeta and [OIII] emission line maps show that the ionised gas has a complex morphology and kinematics, including both a component counter-rotating with respect to the stars and a fainter co-rotating one. Analysis of the absorption line-strength maps show that NGC7332 is young everywhere. The presence of a large-scale bar can explain most of those properties, but the fact that we see a significant amount of unsettled gas, together with a few peculiar features in the maps, suggest that NGC7332 is still evolving. Interactions as well as bar-driven processes must thus have played an important role in the formation and evolution of NGC7332, and presumably of S0 galaxies in general.

The Nearby Supernova Factory (SNfactory) is an ambitious project to find and study in detail approximately 300 nearby Type Ia supernovae (SNe~Ia) at redshifts 0.03<z<0.08. This program will provide an exceptional data set of well-studied SNe in the nearby smooth Hubble flow that can be used as calibration for the current and future programs designed to use SNe to measure the cosmological parameters. The first key ingredient for this program is a reliable supply of Hubble-flow SNe systematically discovered in unprecedented numbers using the same techniques as those used in distant SNe searches. In 2002, 35 SNe were found using our test-bed pipeline for automated SN search and discovery. The pipeline uses images from the asteroid search conducted by the Near Earth Asteroid Tracking group at JPL. Improvements in our subtraction techniques and analysis have allowed us to increase our effective SN discovery rate to ~12 SNe/month in 2003.

We present a unique set of nested stellar kinematical maps of NGC 3377 obtained with the integral-field spectrographs OASIS and SAURON. We then construct general axisymmetric dynamical models for this galaxy, based on the Schwarzschild numerical orbit superposition technique applied to these complementary measurements. We show how these two datasets constrain the mass of the central massive object and the overall mass-to-light ratio of the galaxy by probing the inner and outer regions respectively. The simultaneous use of both datasets leads us to confirm the presence of a massive black hole with a mass of $M_{BH} = 7_{-5}^{+4} 10^{7} M_\sun$ (99.7% confidence level), with a best-fit stellar mass-to-light ratio $\Upsilon_I = 2.1 \pm 0.2$ (for an assumed edge-on inclination).

We present a summary of high-spatial resolution follow-up observations of the elliptical (E) and lenticular (S0) galaxies in the SAURON survey using the OASIS integral field spectrograph. The OASIS observations explore the central 8x10" regions of these galaxies using a spatial sampling four times higher than SAURON, often revealing previously undiscovered features. Around 75% (31/48) of the SAURON E/S0s with central velocity dispersion >= 120 km/s were observed with OASIS, covering well the original SAURON representative sample. We present here an overview of this follow-up survey, and some preliminary results on individual objects, including a previously unreported counter-rotating core in NGC 4382; the decoupled stellar and gas velocity fields of NGC 2768; and the strong age gradient towards the centre of NGC 3489.

We describe our program for the dynamical modeling of early-type galaxies observed with the panoramic integral-field spectrograph SAURON. We are using Schwarzschild's numerical orbit superposition method to reproduce in detail all kinematical and photometric observables, and recover the intrinsic orbital structure of the galaxies. Since catastrophes are the most prominent features in the orbital observables, two-dimensional kinematical coverage is essential to constrain the dynamical models.

We present new techniques to perform adaptive spatial binning of Integral-Field Spectroscopic (IFS) data to reach a chosen constant signal-to-noise ratio per bin. These methods are required for the proper analysis of IFS observations, but can also be used for standard photometric imagery or any other two-dimensional data. Various schemes are tested and compared by binning and extracting the stellar kinematics of the Sa galaxy NGC2273 from spectra obtained with the panoramic IFS SAURON.

We have measured the two-dimensional kinematics and line-strength distributions of 72 representative nearby early-type galaxies, out to approximately one effective radius, with our panoramic integral-field spectrograph SAURON. The resulting maps reveal a rich variety in kinematical structures and linestrength distributions, indicating that early-type galaxies are more complex systems than often assumed. We are building detailed dynamical models for these galaxies, to derive their intrinsic shape and dynamical structure, and to determine the mass of the supermassive central black hole. Here we focus on two examples, the compact elliptical M32 and the E3 galaxy NGC4365. These objects represent two extreme cases: M32 has very regular kinematics which can be represented accurately by an axisymmetric model in which all stars rotate around the short axis, while NGC4365 is a triaxial galaxy with a prominent kinematically decoupled core, with an inner core that rotates about an axis that is nearly perpendicular to the rotation axis of the main body of the galaxy. Our dynamical models for these objects demonstrate that two-dimensional observations are essential for deriving the intrinsic orbital structure and dark matter content of galaxies.

We present adaptive optics assisted OASIS integral field spectrography of the S0 galaxy NGC 4621. Two-dimensional stellar kinematical maps (mean velocity and dispersion) reveal the presence of a 60 pc diameter counter-rotating core (CRC), the smallest observed to date. The OASIS data also suggests that the kinematic center of the CRC is slightly offset from the center of the outer isophotes. This seems to be confirmed by archival HST/STIS data. We also present the HST/WFPC2 V-I colour map, which exhibits a central elongated red structure, also slightly off-centered in the same direction as the kinematic centre. We then construct an axisymmetric model of NGC 4621: the two-integral distribution function is derived using the Multi-Gaussian Expansion and the Hunter & Qian (1993) formalisms. Although the stellar velocities are reasonably fitted, including the region of the counter-rotating core, significant discrepancies between the model and the observations demonstrate the need for a more general model (e.g. a three-integral model).

We present results from a new and unique integral-field spectrograph, SAURON. It has a large field of view and high throughput and is primarily built for the study of stellar & gaseous kinematics and stellar populations in galaxies. Its aim is to carry out a systematic survey of the velocity fields, velocity dispersions, and line-strength distributions of nearby ellipticals, lenticular galaxies and spiral bulges. Its wide field is especially useful for the study of complicated velocity structures. Together with other spectroscopic data, images, and dynamical modelling, SAURON will help to constrain the intrinsic shapes, mass-to-light ratios, and stellar populations of early-type galaxies and spiral bulges.

We present dynamical models of the nearby compact elliptical galaxy M32, using high quality kinematical measurements, obtained with the integral-field spectrograph SAURON mounted on the William Herschel Telescope on La Palma. We also include STIS data obtained by Joseph et al. We find a best-fit black hole mass of (2.5 +/- 0.5) million solar masses and a stellar I-band mass-to-light ratio of (1.85 +/- 0.15) in solar units. For the first time, we are also able to constrain the inclination along which M32 is observed to (70 +/- 5) degrees. Combined with an averaged observed flattening of 0.73, this corresponds to an intrinsic flattening of approximately 0.68 +/- 0.03. These tight constraints are mainly caused by the use of integral-field data. We show this quantitatively by comparing with models that are constrained by multiple slits only. We show the phase-space distribution and intrinsic velocity structure of the best-fit model and investigate the effect of regularisation on the orbit distribution.

We present the SAURON project, which is aimed at studying the morphology, two-dimensional kinematics and stellar populations of a representative sample of elliptical galaxies and spiral bulges. SAURON, a dedicated integral-field spectrograph that is optimized for wide-field observations and has high throughput, was built in Lyon and is now operated at the WHT 4.2m telescope. At present, we have observed approximately two thirds of the seventy-two sample galaxies with SAURON. A comparison with published long-slit measurements demonstrates that the SAURON-data is of equal or better quality, and provides full two-dimensional coverage. The velocity and velocity dispersion fields exhibit a large variety of morphologies: from simple rotating systems to cylindrical, disky and triaxial velocity fields, bars and decoupled cores. Most of these kinematical signatures do not have counterparts in the light distribution. While some galaxies are consistent with axisymmetry, most are more complex systems than assumed previously. This suggests that the kinematical properties of nearby E/S0 galaxies do not agree with the often assumed simplistic two-family model, in which the giant non-rotating triaxial ellipticals are opposed to the fast-rotating axisymmetric faint ellipticals and S0s.

The panoramic integral-field spectrograph SAURON is currently being used to map the stellar kinematics, gaseous kinematics, and stellar populations of a large number of early-type galaxies and bulges. Here, we describe SAURON observations of cold stellar disks embedded in spheroids (NGC3384, NGC4459, NGC4526), we illustrate the kinematics and ionization state of large-scale gaseous disks (NGC4278, NGC7742), and we show preliminary comparisons of SAURON data with barred galaxy N-body simulations (NGC3623).

The SAURON project will deliver two-dimensional spectroscopic data of a sample of nearby early-type galaxies with unprecedented quality. In this paper, we focus on the mapping of their stellar populations using the SAURON data, and present some preliminary results on a few prototypical cases.

We present new techniques to perform adaptive spatial binning of two-dimensional (2D) data to reach a chosen constant signal-to-noise ratio per bin. These methods are required particularly for the proper analysis of Integral Field Spectrograph (IFS) observations, but can also be used for standard photometric imagery. Various schemes are tested and compared using data obtained with the panoramic IFS SAURON.

Early results are reported from the SAURON survey of the kinematics and stellar populations of a representative sample of nearby E, S0 and Sa galaxies. The survey is aimed at determining the intrinsic shape of the galaxies, their orbital structure, the mass-to-light ratio as a function of radius, the age and metallicity of the stellar populations, and the frequency of kinematically decoupled cores and nuclear black holes. The construction of the representative sample is described, and its properties are illustrated. A comparison with long-slit spectroscopic data establishes that the SAURON measurements are comparable to, or better than, the highest-quality determinations. Comparisons are presented for NGC 3384 and NGC 4365 where stellar velocities and velocity dispersions are determined to a precision of 6 km/s, and the h3 and h4 parameters of the line-of-sight velocity distribution to a precision of better than 0.02. Extraction of accurate gas emission-line intensities, velocities and line widths from the datacubes is illustrated for NGC 5813. Comparisons with published line-strengths for NGC 3384 and NGC 5813 reveal uncertainties of < 0.1 A on the measurements of the Hbeta, Mgb and Fe5270 indices. Integral-field mapping uniquely connects measurements of the kinematics and stellar populations to the galaxy morphology. The maps presented here illustrate the rich stellar kinematics, gaseous kinematics, and line-strength distributions of early-type galaxies. The results include the discovery of a thin, edge-on, disk in NGC 3623, confirm the axisymmetric shape of the central region of M32, illustrate the LINER nucleus and surrounding counter-rotating star-forming ring in NGC 7742, and suggest a uniform stellar population in the decoupled core galaxy NGC 5813.

A summary of the SAURON project and its current status is presented. SAURON is a panoramic integral-field spectrograph designed to study the stellar kinematics, gaseous kinematics, and stellar populations of spheroids. Here, the sample of galaxies and its properties are described. The instrument is detailed and its capabilities illustrated through observational examples. These includes results on the structure of central stellar disks, the kinematics and ionization state of gaseous disks, and the stellar populations of galaxies with decoupled cores.

We describe a new technique, kinemetry, to quantify kinematic maps of early-type galaxies in an efficient way. We present the first applications to velocity fields obtained with the integral-field spectrograph SAURON.

Early results are reported from the SAURON survey of the kinematics and stellar populations of nearby cluster and field E, S0, and Sa galaxies. We present maps of the Hbeta and [OIII] emission-line distribution and kinematics for NGC 5813 and NGC 7742.

A new integral-field spectrograph, SAURON, is described. It is based on the Tiger principle, and uses a lenslet array. SAURON has a large field of view and high throughput, and allows simultaneous sky subtraction. Its design is optimized for studies of the stellar kinematics, gas kinematics, and line-strength distributions of nearby early-type galaxies. The instrument design and specifications are described, as well as the extensive analysis software which was developed to obtain fully calibrated spectra, and the associated kinematic and line-strength measurements. A companion paper reports on the first results obtained with SAURON on the William Herschel Telescope.

We present observations with the adaptive optics assisted integral field spectrograph OASIS of the M 31 double nucleus at a spatial resolution better than 0.5 arcsec FWHM. These data are used to derive the two-dimensional stellar kinematics within the central 2 arcsec. Archival WFPC2/HST images are revisited to perform a photometric decomposition of the nuclear region. We also present STIS/HST kinematics obtained from the archive. The luminosity distribution of the central region is well separated into the respective contributions of the bulge, the nucleus including P1 and P2, and the so-called UV peak. We then show that the axis joining P1 and P2, the two local surface brightness maxima, does not coincide with the kinematic major-axis, which is also the major-axis of the nuclear isophotes (excluding P1). We also confirm that the velocity dispersion peak is offset by ~ 0.2 arcsec from the UV peak, assumed to mark the location of the supermassive black hole. The newly reduced STIS/HST velocity and dispersion profiles are then compared to OASIS and other published kinematics. We find significant offsets with previously published data. Simple parametric models are then built to successfully reconcile all the available kinematics. We finally interpret the observations using new N-body simulations. The nearly keplerian nuclear disk of M31 is subject to a natural m=1 mode, with a very slow pattern speed (3 km/s/pc for M_BH = 7 10^7~\Msun), that can be maintained during more than a thousand dynamical times. The resulting morphology and kinematics of the mode can reproduce the M~31 nuclear-disk photometry and mean stellar velocity, including the observed asymmetries. It requires a central mass concentration and a cold disk system representing between 20 and 40% of its mass. Abridged..

We report the first wide-field mapping of the kinematics and stellar populations in the E3 galaxy NGC 4365. The velocity maps extend previous long-slit work. They show two independent kinematic subsystems: the central 300 x 700 pc rotates about the projected minor axis, and the main body of the galaxy, 3 x 4 kpc, rotates almost at right angles to this. The line-strength maps show that the metallicity of the stellar population decreases from a central value greater than solar, to one-half solar at a radius of 2 kpc. The decoupled core and main body of the galaxy have the same luminosity-weighted age, of ~14 Gyr, and the same elevated magnesium-to-iron ratio. The two kinematically distinct components have thus shared a common star formation history. We infer that the galaxy underwent a sequence of mergers associated with dissipative star formation that ended >12 Gyr ago. The misalignment between the photometric and kinematic axes of the main body is unambiguous evidence of triaxiality. The similarity of the stellar populations in the two components suggests that the observed kinematic structure has not changed substantially in 12 Gyr.

We briefly describe the SAURON project, aimed at determining the intrinsic shape and internal dynamics of spheroids. We focus here on the ability of SAURON to identify gaseous and stellar disks and to measure their morphology and kinematics. We illustrate some of our results with complete maps of NGC3377, NGC3623, and NGC4365.

We have extended the spectral dynamics formalism introduced by Binney & Spergel, and have implemented a semi-analytic method to represent regular orbits in any potential, making full use of their regularity. We use the spectral analysis code of Carpintero & Aguilar to determine the nature of an orbit (irregular, regular, resonant, periodic) from a short-time numerical integration. If the orbit is regular, we approximate it by a truncated Fourier time series of a few tens of terms per coordinate. Switching to a description in action-angle variables, this corresponds to a reconstruction of the underlying invariant torus. We then relate the uniform distribution of a regular orbit on its torus to the non-uniform distribution in the space of observables by a simple Jacobian transformation between the two sets of coordinates. This allows us to compute, in a cell-independent way, all the physical quantities needed in the study of the orbit, including the density and in the line-of-sight velocity distribution, with much increased accuracy. The resulting flexibility in the determination of the orbital properties, and the drastic reduction of storage space for the orbit library, provide a significant improvement in the practical application of Schwarzschild's orbit superposition method for constructing galaxy models. We test and apply our method to two-dimensional orbits in elongated discs, and to the meridional motion in axisymmetric potentials, and show that for a given accuracy, the spectral dynamics formalism requires an order of magnitude fewer computations than the more traditional approaches.

In this paper, historical optical(UBVRI) data and newly observed data from the Yunnan Observatory of China(about100 years) are presented for BL Lacertae. Maximum variations in UBVRI: 5.12, 5.31, 4.73, 2.59, and 2.54 and color indices of U-B = -0.11 +/- 0.20, B-V= 1.0 +/- 0.11, V-R= 0.73 +/- 0.19, V-I= 1.42 +/- 0.25, R-I= 0.82 +/- 0.11, and B-I= 2.44 +/- 0.29 have been obtained from the literature; The Jurkevich method is used to investigate the existence of periods in the B band light curve, and a long-term period of 14 years is found. The 0.6 and 0.88 year periods reported by Webb et al.(1988) are confirmed. In addition, a close relation between B-I and B is found, suggesting that the spectra flattens when the source brightens.

In this paper, we will discuss the long-term variations in the optical and the infrared bands of some AGNs. It is interesting to note that the reported periods of AGNs are of the similar values (about 10 years, see Fan et al. 1998a, ApJ, 507 in press, and references therein). DCF method shows the optical and infrared bands are strongly correlated suggesting that the emission mechanisms in the two bands have a common origin.

In this paper, the correlation between the gamma-ray and the radio bands is investigated. The results show that there is a closer correlation between the gamma-ray emission and the high frequency (1.3mm, 230GHz) radio emission for maximum data than between the gamma-ray and the lower frequency (5GHz) radio emissions, which means that the gamma-ray is associated with the radio emission from the jet.

In this paper, the long-term historical optical (UBVRI) and near-infrared (JHK) data are presented with some new observations in the optical (February 1994-January 1995) and near-infrared (November 1995) bands included for BL Lac object OJ287. The new optical data in V-band are in agreement with the results reported by other authors (Sillanpaa et al. 1996a; Arimoto et al. 1997), a close correlation between the color index of B-V and the magnitude V has been obtained from our new observations; The new infrared observations presented here indicate that the source was at a high level in the infrared band during the observation period; From the available literature, we have got that the largest variations for UBVRIJHK bands are respectively: 4.72mag, 5.93mag, 5.18mag, 4.45mag, 4.07mag, 3.87mag, 3.78mag, and 3.54mag. A strong correlation is found between the optical and near-infrared bands when the DCF method is used, which suggests that these two bands have the same emission mechanism.

We have obtained 2640 CCD spectra with resolution ~4 Angstrom in the region 7250-9000 Angstroms for 976 stars lying near the red giant branches in color-magnitude diagrams of 52 Galactic globular clusters. Radial velocities of ~16 km/second accuracy per star determined from the spectra are combined with other criteria to assess quantitative membership probabilities. Measurements of the equivalent widths of the infrared calcium triplet lines yield a relative metal-abundance ranking with a precision that compares favorably to other techniques. Regressions between our system and those of others are derived. Our reduction procedures are discussed in detail, and the resultant catalog of derived velocities and equivalent widths is presented. The metal abundances derived from these data will be the subject of a future paper.