4MOST, the 4m Multi Object Spectroscopic Telescope, is the optical, fibre-fed, MOS facility for the VISTA telescope at ESO's Paranal Observatory in Chile. Its main science drivers are in the fields of galactic archeology, high-energy physics, galaxy evolution and cosmology. The 4MOST consortium consists of several institutes in Europe and Australia under leadership of the Leibniz-Institut fur Astrophysik Potsdam (AIP). This paper focuses on the successful testing, installation and technical commissioning of the Low Resolution Spectrographs (LRS-A and B) for the 4MOST instrument at ESO's Paranal Observatory, Chile. This work was completed on October 18. 2025. Details on the assembly, integration, and performance of both 4MOST spectrographs from their arrival in the integration hall through to the telescope installation are provided. Attention is given to the optimization of procedures implemented to enhance performance and meet the expected top-level requirements. The 4MOST LRS features 2436 fibres split into two low-resolution spectrographs LRS-A and LRS-B (1624 fibres, three arms, 370-950 nm, R > 4000) and one high-resolution spectrograph (812 fibres, three arms, 44-69 nm coverage each, R > 18000). The fibre positioner covers a hexagonal field of view of 4.1 deg2. The fibers are 85 μm core with an output beam at f/3. The Centre de Recherche Astrophysique de Lyon (CRAL) had the full responsibility for the two low-resolution spectrographs. Each of them is composed of an off-axis Schmidt collimator that produces a 200 mm beam, which is split into three spectral arms by dichroics and directed to F/1.73 cameras with standard 6k x 6k 15 microns pixel CCD detectors.

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.

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 autoencoder that is interpreted probabilistically after training using a normalizing flow. 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 multistage 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 the 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 ± 0.010 mag, which corresponds to 0.074 ± 0.010 mag if peculiar velocity contributions are removed. Trained models and codes are released at https://github.com/georgestein/suPAErnova.

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.

4MOST, the 4m Multi Object Spectroscopic Telescope, is an upcoming optical, fiber-fed, MOS facility for the VISTA telescope at ESO's Paranal Observatory in Chile. Its main science drivers are in the fields of galactic archeology, high-energy physics, galaxy evolution and cosmology. The 4MOST consortium consists of several institutes in Europe and Australia under the leadership of the Leibniz-Institut für Astrophysik Potsdam (AIP). 4MOST is currently in its Assembly, Integration and Tests Phase with an expected start of science operations in 2023. The design of 4MOST features 2436 fibers split into two low-resolution spectrographs (1624 fibers, 370-950 nm, R < 4000) and a high-resolution spectrograph (812 fibers, ~44-69 nm coverage, R < 18000). A fiber positioner covers a hexagonal field of view of ~4.1 deg². CRAL has the full responsibility of the Low Resolution Spectrographs. Each of them is composed of an off-axis collimator collecting the f/3 beam coming from the fibers and providing a 200 mm collimated beam, which is then split into three arms thanks to two dichroics. Each arm is composed of a Schmidt corrector, a VPHG, an f/1.73 camera and a standard 6k by 6k 15µm pixel detector. Manufacturing, Assembly, Integration and Test (MAIT) phase of the LRS started in May 2018. The first LRS has now been fully integrated, aligned and tested and its Local Acceptance Review on site at CRAL passed successfully in December 2021. This paper describes the procedures developed to test the spectrograph and demonstrate its compliance with the requirements. Performances achieved are presented throughout the paper, especially in terms of spectral resolving power and sampling, spectral purity, cross-talk, throughput, straylight, light tightness and stability.

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 ± 0.007 mag, which corresponds to 0.084 ± 0.009 mag if peculiar velocity contributions are removed and to 0.073 ± 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 ± 0.045 mag. Standardization using the Twins Embedding also significantly decreases host-galaxy correlations. We recover a host mass step of 0.040 ± 0.020 mag compared to 0.092 ± 0.026 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.

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 Å. 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 can be applied to the standardization of SNe Ia for cosmology.

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 (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 part of an on-going effort to identify, understand and correct for astrophysics biases in the standardization of Type Ia supernovae (SN Ia) for cosmology, we have statistically classified a large sample of nearby SNe Ia into those that are located in predominantly younger or older environments. This classification is based on the specific star formation rate measured within a projected distance of 1 kpc 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 SNe Ia in predominantly younger environments are ΔY = 0.163 ± 0.029 mag (5.7σ) fainter than those in predominantly older environments after conventional light-curve standardization. This is the strongest standardized SN Ia brightness systematic connected to the host-galaxy environment measured to date. The well-established step in standardized brightnesses between SNe Ia in hosts with lower or higher total stellar masses is smaller, at ΔM = 0.119 ± 0.032 mag (4.5σ), for the same set of SNe Ia. When fit simultaneously, the environment-age offset remains very significant, with ΔY = 0.129 ± 0.032 mag (4.0σ), while the global stellar mass step is reduced to ΔM = 0.064 ± 0.029 mag (2.2σ). Thus, approximately 70% of the variance from the stellar mass step is due to an underlying dependence on environment-based progenitor age. Also, we verify that using the local star formation rate alone is not as powerful as LsSFR at sorting SNe Ia into brighter and fainter subsets. Standardization that only uses the SNe Ia in younger environments reduces the total dispersion from 0.142 ± 0.008 mag to 0.120 ± 0.010 mag. We show that as environment-ages evolve with redshift, a strong bias, especially on the measurement of the derivative of the dark energy equation of state, can develop. Fortunately, data that measure and correct for this effect using our local specific star formation rate indicator, are likely to be available for many next-generation SN Ia cosmology experiments.

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.

Author(s): Rigault, M; Copin, Y; Aldering, G; Antilogus, P; Aragon, C; Bailey, S; Baltay, C; Bongard, S; Buton, C; Canto, A; Cellier-Holzem, F; Childress, M; Chotard, N; Fakhouri, HK; Feindt, U; Fleury, M; Gangler, E; Greskovic, P; Guy, J; Kim, AG; Kowalski, M; Lombardo, S; Nordin, J; Nugent, P; Pain, R; Pecontal, E; Pereira, R; Perlmutter, S; Rabinowitz, D; Runge, K; Saunders, C; Scalzo, R; Smadja, G; Tao, C; Thomas, RC; Weaver, BA; Factory, Nearby Supernova

Author(s): Rigault, M; Copin, Y; Aldering, G; Antilogus, P; Aragon, C; Bailey, S; Baltay, C; Bongard, S; Buton, C; Canto, A; Cellier-Holzem, F; Childress, M; Chotard, N; Fakhouri, HK; Feindt, U; Fleury, M; Gangler, E; Greskovic, P; Guy, J; Kim, AG; Kowalski, M; Lombardo, S; Nordin, J; Nugent, P; Pain, R; Pecontal, E; Pereira, R; Perlmutter, S; Rabinowitz, D; Runge, K; Saunders, C; Scalzo, R; Smadja, G; Tao, C; Thomas, RC; Weaver, BA; Factory, Nearby Supernova

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 introduce a method for identifying “twin” Type Ia supernovae, and using them to improve distance measurements. This novel approach to Type Ia supernova stan1 Physics Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA, 94720 2 Department of Physics, University of California Berkeley, 366 LeConte Hall MC 7300, Berkeley, CA, 94720-7300 3 Laboratoire de Physique Nucléaire et des Hautes Énergies, Université Pierre et Marie Curie Paris 6, Université Paris Diderot Paris 7, CNRS-IN2P3, 4 place Jussieu, 75252 Paris Cedex 05, France 4 Department of Physics, Yale University, New Haven, CT, 06250-8121 5 Tsinghua Center for Astrophysics, Tsinghua University, Beijing 100084, China 6 Université de Lyon, F-69622, Lyon, France ; Université de Lyon 1, Villeurbanne ; CNRS/IN2P3, Institut de Physique Nucléaire de Lyon. 7 Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia 8 Institut fur Physik, Humboldt-Universitat zu Berlin, Newtonstr. 15, 12489 Berlin 9 Centre de Physique des Particules de Marseille, Aix-Marseille Université , CNRS/IN2P3, 163 avenue de Luminy Case 902 13288 Marseille Cedex 09, France 10 Clermont Université, Université Blaise Pascal, CNRS/IN2P3, Laboratoire de Physique Corpusculaire, BP 10448, F-63000 Clermont-Ferrand, France 11 DESY, D-15735 Zeuthen, Germany 12 Centre de Recherche Astronomique de Lyon, Université Lyon 1, 9 Avenue Charles André, 69561 Saint Genis Laval Cedex, France 13 Department of Physics, Florida State University, 315 Keen Building, Tallahassee, FL 32306-4350 14 Computational Cosmology Center, Computational Research Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road MS 50B-4206, Berkeley, CA, 94720 15 Center for Cosmology and Particle Physics, New York University, 4 Washington Place, New York, NY 10003, USA

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 λ(uNi), λ(uTi), λ(uSi) and λ(uCa). Two of these span the Ca H&K λλ 3934, 3969 complex. We employ spectral synthesis using SYNAPPS to associate the two bluer regions with Ni/Co and Ti. Results. The flux of the uTi feature is an extremely sensitive temperature/luminosity indicator, standardizing the SN peak luminosity to 0.116 ± 0.011 mag root mean square (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. 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 light curve power. The uCa flux further improves standardization, yielding a 0.086 ± 0.010 mag RMS without the need to include an additional intrinsic dispersion to reach χ2∕dof ~ 1. This reduction in RMS is partially driven by an improved standardization of Shallow Silicon and 91T-like SNe.