Search

In the next decade, many optical surveys will aim to tackle the question of dark energy nature, measuring its equation of state parameter at the permil level. This requires trusting the photometric calibration of the survey with a precision never reached so far, controlling many sources of systematic uncertainties. The measurement of the on-site atmospheric transmission for each exposure, or on average for each season or for the full survey, can help reach the permil precision for magnitudes. This work aims at proving the ability to use slitless spectroscopy for standard star spectrophotometry and its use to monitor on-site atmospheric transmission as needed, for example, by the Vera C. Rubin Observatory Legacy Survey of Space and Time supernova cosmology program. We fully deal with the case of a disperser in the filter wheel, which is the configuration chosen in the Rubin Auxiliary Telescope. The theoretical basis of slitless spectrophotometry is at the heart of our forward model approach to extract spectroscopic information from slitless data. We developed a publicly available software called Spectractor (https://github.com/LSSTDESC/Spectractor) that implements each ingredient of the model and finally performs a fit of a spectrogram model directly on image data to get the spectrum. We show on simulations that our model allows us to understand the structure of spectrophotometric exposures. We also demonstrate its use on real data, solving specific issues and illustrating how our procedure allows the improvement of the model describing the data. Finally, we discuss how this approach can be used to directly extract atmospheric transmission parameters from data and thus provide the base for on-site atmosphere monitoring. We show the efficiency of the procedure on simulations and test it on the limited data set available., Comment: 30 pages, 36 figures, published version in A&A

We apply the color-magnitude intercept calibration method (CMAGIC) to the Nearby Supernova Factory SNe Ia spectrophotometric dataset. The currently existing CMAGIC parameters are the slope and intercept of a straight line fit to the first linear region in the color-magnitude diagram, which occurs over a span of approximately 30 days after maximum brightness. We define a new parameter, $\omega_{XY}$, the size of the ``bump'' feature near maximum brightness for arbitrary filters $X$ and $Y$. We find a significant correlation between the slope of the first linear region, $\beta_{XY, 1}$, in the CMAGIC diagram and $\omega_{XY}$. These results may be used to our advantage, as they are less affected by extinction than parameters defined as a function of time. Additionally, $\omega_{XY}$ is computed independently of templates. We find that current empirical templates are successful at reproducing the features described in this work, particularly SALT3, which correctly exhibits the negative correlation between slope and bump size seen in our data. In 1-D simulations, we show that the correlation between the size of the bump feature and $\beta_{XY, 1}$ can be understood as a result of chemical mixing due to large-scale Rayleigh-Taylor instabilities., Comment: 19 pages, 9 figures

We construct a physically-parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an Auto-Encoder (AE) which is interpreted probabilistically after training using a Normalizing Flow (NF). We demonstrate that the PAE learns a low-dimensional latent space that captures the nonlinear range of features that exists within the population, and can accurately model the spectral evolution of SNe Ia across the full range of wavelength and observation times directly from the data. By introducing a correlation penalty term and multi-stage training setup alongside our physically-parameterized network we show that intrinsic and extrinsic modes of variability can be separated during training, removing the need for the additional models to perform magnitude standardization. We then use our PAE in a number of downstream tasks on SNe Ia for increasingly precise cosmological analyses, including automatic detection of SN outliers, the generation of samples consistent with the data distribution, and solving the inverse problem in the presence of noisy and incomplete data to constrain cosmological distance measurements. We find that the optimal number of intrinsic model parameters appears to be three, in line with previous studies, and show that we can standardize our test sample of SNe Ia with an RMS of $0.091 \pm 0.010$ mag, which corresponds to $0.074 \pm 0.010$ mag if peculiar velocity contributions are removed. Trained models and codes are released at \href{https://github.com/georgestein/suPAErnova}{github.com/georgestein/suPAErnova}, Comment: 23 pages, 8 Figures, 1 Table. Accepted to ApJ

We calibrate spectrophotometric optical spectra of 32 stars commonly used as standard stars, referenced to 14 stars already on the HST-based CALSPEC flux system. Observations of CALSPEC and non-CALSPEC stars were obtained with the SuperNova Integral Field Spectrograph over the wavelength range 3300 A to 9400 A as calibration for the Nearby Supernova Factory cosmology experiment. In total, this analysis used 4289 standard-star spectra taken on photometric nights. As a modern cosmology analysis, all pre-submission methodological decisions were made with the flux scale and external comparison results blinded. The large number of spectra per star allows us to treat the wavelength-by-wavelength calibration for all nights simultaneously with a Bayesian hierarchical model, thereby enabling a consistent treatment of the Type Ia supernova cosmology analysis and the calibration on which it critically relies. We determine the typical per-observation repeatability (median 14 mmag for exposures >~ 5 s), the Maunakea atmospheric transmission distribution (median dispersion of 7 mmag with uncertainty 1 mmag), and the scatter internal to our CALSPEC reference stars (median of 8 mmag). We also check our standards against literature filter photometry, finding generally good agreement over the full 12-magnitude range. Overall, the mean of our system is calibrated to the mean of CALSPEC at the level of ~ 3 mmag. With our large number of observations, careful crosschecks, and 14 reference stars, our results are the best calibration yet achieved with an integral-field spectrograph, and among the best calibrated surveys., Comment: Accepted for publication in ApJS

We study the spectral diversity of Type Ia supernovae (SNe Ia) at maximum light using high signal-to-noise spectrophotometry of 173 SNe Ia from the Nearby Supernova Factory. We decompose the diversity of these spectra into different extrinsic and intrinsic components, and we construct a nonlinear parameterization of the intrinsic diversity of SNe Ia that preserves pairings of "twin" SNe Ia. We call this parameterization the "Twins Embedding". Our methodology naturally handles highly nonlinear variability in spectra, such as changes in the photosphere expansion velocity, and uses the full spectrum rather than being limited to specific spectral line strengths, ratios or velocities. We find that the time evolution of SNe Ia near maximum light is remarkably similar, with 84.6% of the variance in common to all SNe Ia. After correcting for brightness and color, the intrinsic variability of SNe Ia is mostly restricted to specific spectral lines, and we find intrinsic dispersions as low as ~0.02 mag between 6600 and 7200 A. With a nonlinear three-dimensional model plus one dimension for color, we can explain 89.2% of the intrinsic diversity in our sample of SNe Ia, which includes several different kinds of "peculiar" SNe Ia. A linear model requires seven dimensions to explain a comparable fraction of the intrinsic diversity. We show how a wide range of previously-established indicators of diversity in SNe Ia can be recovered from the Twins Embedding. In a companion article, we discuss how these results an be applied to standardization of SNe Ia for cosmology., Comment: Accepted to ApJ

We show how spectra of Type Ia supernovae (SNe Ia) at maximum light can be used to improve cosmological distance estimates. In a companion article, we used manifold learning to build a three-dimensional parameterization of the intrinsic diversity of SNe Ia at maximum light that we call the "Twins Embedding". In this article, we discuss how the Twins Embedding can be used to improve the standardization of SNe Ia. With a single spectrophotometrically-calibrated spectrum near maximum light, we can standardize our sample of SNe Ia with an RMS of $0.101 \pm 0.007$ mag, which corresponds to $0.084 \pm 0.009$ mag if peculiar velocity contributions are removed and $0.073 \pm 0.008$ mag if a larger reference sample were obtained. Our techniques can standardize the full range of SNe Ia, including those typically labeled as peculiar and often rejected from other analyses. We find that traditional light curve width + color standardization such as SALT2 is not sufficient. The Twins Embedding identifies a subset of SNe Ia including but not limited to 91T-like SNe Ia whose SALT2 distance estimates are biased by $0.229 \pm 0.045$ mag. Standardization using the Twins Embedding also significantly decreases host-galaxy correlations. We recover a host mass step of $0.040 \pm 0.020$ mag compared to $0.092 \pm 0.024$ mag for SALT2 standardization on the same sample of SNe Ia. These biases in traditional standardization methods could significantly impact future cosmology analyses if not properly taken into account., Comment: Accepted to ApJ

We study astrometric residuals from a simultaneous fit of Hyper Suprime-Cam images. We aim to characterize these residuals and study the extent to which they are dominated by atmospheric contributions for bright sources. We use Gaussian process interpolation, with a correlation function (kernel), measured from the data, to smooth and correct the observed astrometric residual field. We find that Gaussian process interpolation with a von K\'arm\'an kernel allows us to reduce the covariances of astrometric residuals for nearby sources by about one order of magnitude, from 30 mas$^2$ to 3 mas$^2$ at angular scales of ~1 arcmin, and to halve the r.m.s. residuals. Those reductions using Gaussian process interpolation are similar to recent result published with the Dark Energy Survey dataset. We are then able to detect the small static astrometric residuals due to the Hyper Suprime-Cam sensors effects. We discuss how the Gaussian process interpolation of astrometric residuals impacts galaxy shape measurements, in particular in the context of cosmic shear analyses at the Rubin Observatory Legacy Survey of Space and Time.

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$., Comment: 5 pages

Tearing patterns affecting flat field frames in CCDs are a visually striking obstacle to performing Pixel Response Non-Uniformity corrections. These patterns can be explained by lateral field distortions, caused by the non-uniform distribution of holes in the channel stops between sensor columns. Over the course of LSST camera development, a number of practical fixes have been suggested to get rid of tearing. But applying these fixes to our 16-channel Teledyne-e2v sensors leaves at best a distortion pattern at the vertical edges of every segment. Our working hypothesis is that the origin of the tearing is the parallel clocking itself, which moves the holes that are present in the channel stops regions. The efficiency of these transfers depends strongly on the details of the clocking operations, resulting in the observed variety of distortion patterns. Removal of most of the distortion patterns can therefore be achieved by executing a purge operation, which flattens back the hole distribution in the channel stops, immediately before acquiring a frame. A more effective solution is to switch all clocking operations to use a bipolar voltage set., Comment: 11 pages, proceeding of ISPA2018

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., Comment: White paper CNRS/IN2P3 Prospectives 2020

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., Comment: 25 pages, 27 figures

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., Comment: 16 pages, 11 figures, MNRAS in press, missing line added in Table 1

The Photon Transfer Curve (PTC) of a CCD depicts the variance of uniform images as a function of their average. It is now well established that the variance is not proportional to the average, as Poisson statistics would indicate, but rather flattens out at high flux. This "variance deficit", related to the brighter-fatter effect, feeds correlations between nearby pixels, that increase with flux, and decay with distance. We propose an analytical expression for the PTC shape, and for the dependence of correlations with intensity, and relate both to some more basic quantities related to the electrostatics of the sensor, that are commonly used to correct science images for the brighter-fatter effect. We derive electrostatic constraints from a large set of flat field images acquired with a CCD e2v 250, and eventually question the generally-admitted assumption that boundaries of CCD pixels shift by amounts proportional to the source charges. Our results show that the departure of flat field statistics from Poisson law is entirely compatible with charge redistribution during the drift in the sensor., Comment: 16 pages. Accepted for publication in A&A

In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $\gamma$, which captures the effect of gravity on the linear growth of structure through the relation $f=\Omega_M^\gamma$. We present the first projections for the precision in $\gamma$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $\gamma$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity., Comment: Decadal Survey White Paper submission

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

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., Comment: 51 page, 19 figures, accepted in ApJ

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], Comment: Published in Astronomy and Astrophysics Rigault et al. 2020, A&A, 644, A176 | arxiv Rigault et al. 2018

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., Comment: 13 pages, 11 figures, 2 tables, accepted for publication in A&A

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., Comment: 19 pages, 17 figures, accepted for publication in A&A

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible., Comment: 312 pages, 90 figures. Browse the current version at https://github.com/LSSTScienceCollaborations/ObservingStrategy, new contributions welcome!

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., Comment: 17 pages, 15 figures, accepted for publication by A&A in 8 of August 2017

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., Comment: 48 pages, 13 figures. Accepted for publication in The Astrophysical Journal. The spectral time series data presented in this article can be found at http://snfactory.lbl.gov/snf/data/

To take advantage of the astrophysical potential of Gamma-Ray Bursts (GRBs), Chinese and French astrophysicists have engaged the SVOM mission (Space-based multi-band astronomical Variable Objects Monitor). Major advances in GRB studies resulting from the synergy between space and ground observations, the SVOM mission implements space and ground instrumentation. The scientific objectives of the mission put a special emphasis on two categories of GRBs: very distant GRBs at z$>$5 which constitute exceptional cosmological probes, and faint/soft nearby GRBs which allow probing the nature of the progenitors and the physics at work in the explosion. These goals have a major impact on the design of the mission: the on-board hard X-ray imager is sensitive down to 4 keV and computes on line image and rate triggers, and the follow-up telescopes on the ground are sensitive in the NIR. At the beginning of the next decade, SVOM will be the main provider of GRB positions and spectral parameters on very short time scale. The SVOM instruments will operate simultaneously with a wide range of powerful astronomical devices. This rare instrumental conjunction, combined with the relevance of the scientific topics connected with GRB studies, warrants a remarkable scientific return for SVOM. In addition, the SVOM instrumentation, primarily designed for GRB studies, composes a unique multi-wavelength observatory with rapid slew capability that will find multiple applications for the whole astronomy community beyond the specific objectives linked to GRBs. This report lists the scientific themes that will benefit from observations made with SVOM, whether they are specific GRB topics, or more generally all the issues that can take advantage of the multi-wavelength capabilities of SVOM., Comment: 67 pages, 1 frontispiece, 21 figures. Report on the Scientific prospects of the SVOM mission. Proceedings of the Workshop held from 11th to 15th April 2016 at Les Houches School of Physics, France

We employ electrostatic conversion drift calculations to match CCD pixel signal covariances observed in flat field exposures acquired using candidate sensor devices for the LSST Camera. We thus constrain pixel geometry distortions present at the end of integration, based on signal images recorded. We use available data from several operational voltage parameter settings to validate our understanding. Our primary goal is to optimize flux point-spread function (FPSF) estimation quantitatively, and thereby minimize sensor-induced errors which may limit performance in precision astronomy applications. We consider alternative compensation scenarios that will take maximum advantage of our understanding of this underlying mechanism in data processing pipelines currently under development. To quantitatively capture the pixel response in high-contrast/high dynamic range operational extrema, we propose herein some straightforward laboratory tests that involve altering the time order of source illumination on sensors, within individual test exposures. Hence the word {\it hysteretic} in the title of this paper., Comment: 19 pages, 8 figures, 4 appendices. Presented at SPIE Astronomical Telescopes + Instrumentation (paper 9915-45) 26 June-1 July 2016, Edinburgh, Scotland. Fixed an error in right hand expression of Eq.5 and resubmitted to arXiv 160912

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, Comment: 37 pages, 9 figures, 5 tables. Accepted for publication in ApJ. Fixed typo in arXiv abstract

We describe the design, operation, and first results of a photometric calibration project, called DICE (Direct Illumination Calibration Experiment), aiming at achieving precise instrumental calibration of optical telescopes. The heart of DICE is an illumination device composed of 24 narrow-spectrum, high-intensity, light-emitting diodes (LED) chosen to cover the ultraviolet-to-near-infrared spectral range. It implements a point-like source placed at a finite distance from the telescope entrance pupil, yielding a flat field illumination that covers the entire field of view of the imager. The purpose of this system is to perform a lightweight routine monitoring of the imager passbands with a precision better than 5 per-mil on the relative passband normalisations and about 3{\AA} on the filter cutoff positions. The light source is calibrated on a spectrophotometric bench. As our fundamental metrology standard, we use a photodiode calibrated at NIST. The radiant intensity of each beam is mapped, and spectra are measured for each LED. All measurements are conducted at temperatures ranging from 0{\deg}C to 25{\deg}C in order to study the temperature dependence of the system. The photometric and spectroscopic measurements are combined into a model that predicts the spectral intensity of the source as a function of temperature. We find that the calibration beams are stable at the $10^{-4}$ level -- after taking the slight temperature dependence of the LED emission properties into account. We show that the spectral intensity of the source can be characterised with a precision of 3{\AA} in wavelength. In flux, we reach an accuracy of about 0.2-0.5% depending on how we understand the off-diagonal terms of the error budget affecting the calibration of the NIST photodiode. With a routine 60-mn calibration program, the apparatus is able to constrain the passbands at the targeted precision levels., Comment: 25 pages, 27 figures, accepted for publication in A&A

Charge-coupled devices (CCDs) are widely used in astronomy to carry out a variety of measurements, such as for flux or shape of astrophysical objects. The data reduction procedures almost always assume that ther esponse of a given pixel to illumination is independent of the content of the neighboring pixels. We show evidence that this simple picture is not exact for several CCD sensors. Namely, we provide evidence that localized distributions of charges (resulting from star illumination or laboratory luminous spots) tend to broaden linearly with increasing brightness by up to a few percent over the whole dynamic range. We propose a physical explanation for this "brighter-fatter" effect, which implies that flatfields do not exactly follow Poisson statistics: the variance of flatfields grows less rapidly than their average, and neighboring pixels show covariances, which increase similarly to the square of the flatfield average. These covariances decay rapidly with pixel separation. We observe the expected departure from Poisson statistics of flatfields on CCD devices and show that the observed effects are compatible with Coulomb forces induced by stored charges that deflect forthcoming charges. We extract the strength of the deflections from the correlations of flatfield images and derive the evolution of star shapes with increasing flux. We show for three types of sensors that within statistical uncertainties,our proposed method properly bridges statistical properties of flatfields and the brighter-fatter effect.

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., Comment: 3 Figures ; Submitted to ApJ: Oct.~30, 2014 -- Accepted: Dec.~17, 2014

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., Comment: 22 pages, 26 figures, 3 tables, accepted for publication in MNRAS

Near-future astronomical survey experiments, such as LSST, possess system requirements of unprecedented fidelity that span photometry, astrometry and shape transfer. Some of these requirements flow directly to the array of science imaging sensors at the focal plane. Availability of high quality characterization data acquired in the course of our sensor development program has given us an opportunity to develop and test a framework for simulation and modeling that is based on a limited set of physical and geometric effects. In this paper we describe those models, provide quantitative comparisons between data and modeled response, and extrapolate the response model to predict imaging array response to astronomical exposure. The emergent picture departs from the notion of a fixed, rectilinear grid that maps photo-conversions to the potential well of the channel. In place of that, we have a situation where structures from device fabrication, local silicon bulk resistivity variations and photo-converted carrier patterns still accumulating at the channel, together influence and distort positions within the photosensitive volume that map to pixel boundaries. Strategies for efficient extraction of modeling parameters from routinely acquired characterization data are described. Methods for high fidelity illumination/image distribution parameter retrieval, in the presence of such distortions, are also discussed., Comment: 13 pages, 8 figures

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., Comment: 20 pages, 10 figures, accepted to MNRAS

We present evidence that spots imaged using astronomical CCDs do not exactly scale with flux: bright spots tend to be broader than faint ones, using the same illumination pattern. We measure that the linear size of spots or stars, of typical size 3 to 4 pixels FWHM, increase linearly with their flux by up to $2\%$ over the full CCD dynamic range. This brighter-fatter effect affects both deep-depleted and thinned CCD sensors. We propose that this effect is a direct consequence of the distortions of the drift electric field sourced by charges accumulated within the CCD during the exposure and experienced by forthcoming light-induced charges in the same exposure. The pixel boundaries then become slightly dynamical: overfilled pixels become increasingly smaller than their neighbors, so that bright star sizes, measured in number of pixels, appear larger than those of faint stars. This interpretation of the brighter-fatter effect implies that pixels in flat-fields should exhibit statistical correlations, sourced by Poisson fluctuations, that we indeed directly detect. We propose to use the measured correlations in flat-fields to derive how pixel boundaries shift under the influence of a given charge pattern, which allows us to quantitatively predict how star shapes evolve with flux. This physical model of the brighter-fatter effect also explains the commonly observed phenomenon that the spatial variance of CCD flat-fields increases less rapidly than their average., Comment: Proceedings of the workshop "Precision Astronomy with Fully Depleted CCDs". Accepted for publication in JINST. 15 pages

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., Comment: 17 pages, 6 figures. Accepted by Astrophysical Journal